the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Modelling transient thermal processes in the lithosphere: application to the NW Pannonian basin
Abstract. The reconstruction of thermal evolution in sedimentary basins is a key input for constraining geodynamic processes and geo-energy resource potential. We present a methodology to reproduce the most important transient thermal footprints accompanying basin formation: lithosphere extension and sedimentation. The forward model is extended with data assimilation to constrain models with temperature measurements. We apply the methodology to the NW part of Hungary. Realistic past- and present-day temperature predictions for the entire lithosphere are achieved, suggesting the relatively uniform, but strong attenuation of the mantle lithosphere through extension, and relatively small variations in the present-day thermal lithosphere thickness. The new temperature model allows an improved estimation of lithosphere rheology and the interpretation of mantle xenolith origins.
- Preprint
(1377 KB) - Metadata XML
- BibTeX
- EndNote
Status: final response (author comments only)
-
CC1: 'Comment on egusphere-2024-308', Giacomo Medici, 13 Feb 2024
General comments
It’s always good review original paper on large-scale hydro, and thermal models from Hungary! The research is also original and can be exported to many other areas of geothermal interests worldwide. Please, follow my comments to improve the manuscript.
Specific comments
Abstract
Line 10. “The forward model is extended”. Please, be more specific. The object of the sentence is unclear and the abstract is short with obvious chance for clarifications
Introduction
Lines 17-61. Did you consider adding a general statement to steady state and transient modelling in other fields of geo-science? Many large scale (deep and large in plant view) flow models have been developed in the Pannonia Basin. Your country has an original and well recognized academic tradition on this aspect of geo-science.
Lines 16-20. “Understanding...thermal evolution pattern”. Long statement without references. Please, insert recent review papers in the field of geothermal energy for characterization, production and modelling:
- Review of Discrete Fracture Network Characterization for Geothermal Energy Extraction. Frontiers in Earth Science, 11, 1328397
- Direct utilization of geothermal energy 2020 worldwide review. Geothermics, 90, 101915.
Line 44. Clearly state the other hot basins in Europe (e.g., Rhine Graben, Tyrrhenian Sea). They are not so many and you can avoid vague sentences in that way.
Line 61. Specify the 3 to 4 specific objectives of your research by using numbers (e.g., i, ii, and iii).
Data and methods
Line 127. “We calibrated the thermal model with subsurface temperature measurements from hydrocarbon and geothermal wells”. Please, specify the depth of the temperature data used for the calibration. 0.2 - 5.0 km based on geothermal and hydrocarbon observations?
Line 127. If we assume observations 0.2 - 5.0 km, did you discuss reliability/validity of the model much deeper? The model should not be very sensitive in the deeper part.
Line 127 – onwards. Do you need to add some detail on the sensitivity of your model with respect to the parameters?
Line 127. Link the depth range of temperature observations to Figure 3a
Line 181-222. The time steps of your transient model should be much more clear when you describe the methodology. They should be clear from the first lines. Do you need a link with the Table 2?
Discussion
Line 342. “It has already been”. Avoid to start a new sentence with “it”. Please, revise the language.
Line 347. “These factors”. Difficult to follow. Please, remind the specific factors to the reader.
Line 408. I suggest “considering this scenario”. Avoid to use the word “this” alone.
Conclusions
Line 451. Insert a connector such as “indeed” to link the last two sentences.
References
Lines 477-639. Please, integrate relevant literature as suggested above.
Figures and tables
Figure 3a. Please, increase the graphic resolution. Some details are difficult to read.
Figures 5 and 6. Make the figures larger.
Figure 8. Make the letters of the labels larger.
Citation: https://doi.org/10.5194/egusphere-2024-308-CC1 -
AC3: 'Reply on CC1', Eszter Békési, 07 Aug 2024
Dear Giacomo Medici,
Thank you for providing useful suggestions for the manuscript. We incorporated most of the suggestions, which improved the presentation of modelling procedure, description of results and discussion. Please find our detailed point-by-point responses to the comments in the attachment.
Kind regards,
Eszter Békési and co-authors
-
AC3: 'Reply on CC1', Eszter Békési, 07 Aug 2024
-
RC1: 'Comment on egusphere-2024-308', Anonymous Referee #1, 06 Mar 2024
The manuscript by Bekesi et al. entitled ‘Modelling transient thermal processes in the lithosphere: application to the NW Pannonian basin’ presents a simplified modeling study on the thermal evolution of the NW part of the extensional Pannonian basin considering distinct crustal and mantle thinning factors and sedimentation. The calculated new thermal field is then used to present a 2D yield stress section of the lithosphere. Finally, the manuscript contains a brief discussion on mantle xenoliths. Given the large number of major issues of the manuscript, I suggest substantial revision before considering it for publication.
- The title does not reflect the content of the manuscript. Reconstructing the thermal evolution of the lithosphere and particularly the deep lithospheric mantle is challenging, indeed, because of the large number of transient effects, i.e. partial melting, melt emplacement, phase changes, shear heating, non-uniform upper crustal, lower crustal and mantle thinning, paleo-surface temperature variations, basin inversion and related deformation, water circulation, etc. This manuscript uses the stretching factors approach of Royden and Keen (1980) to somehow consider crustal and mantle thinning in a simplified way, but none of the other transient effects are taken into account.
- The abstract and the manuscript claims that one of the main goals is to better quantify the thermal field in the entire lithosphere. There are too problems with this: (1) the model does not use any observational constraints from the deep basins, crust or lithosphere, and likely it is not sensitive to temperature variations at great depths; therefore, the goal cannot be reached with this method. (2) While the manuscript presents one possible model result, , a sensitivity analysis, assessing the role of different initial and boundary conditions and input parameters are missing, therefore, it is not an easy task to see how robust or reliable is the model. No model limitation section is included, despite the large number of assumptions the authors made.
- The used model parameters: Many parameters are not justified and seem to be far from reality. The source of other input data is not clarified and thus cannot be checked. A model needs to be reproduceable by the community and you need to make available the most important input data. 1. Initial crustal thickness: the authors assume a constant 35 km thickness. What is the source of this parameter? The study area includes metamorphic core complexes, their formation requires a thick and hot crust, which infers that your chosen initial values are lower than it should be. Previous reconstructions (e.g. van Hinsbergen et al. 2020) reported a much larger amounts of extension and a thicker initial crust. Geochemical studies based on xenoliths inferred a much larger initial crustal thickness (Torok et al.). Finally, the basement units of the region derived from the Alps, likely having a much thicker crust than 35 km in the Early Miocene. The initial lithospheric thickness of 120km: what is the constraint on this and how much role does it have? Lithologies: what is the source of this? For instance, the sand to shale ratio is proposed to be 1:9 for the ‘Lower Pannonian‘ (Upper Miocene). How is this constrained? After a brief google search, well logs published by Stano et al. 2016 shows a sand to shale ratio of at least 50%. This means that your applied thermal conductivities are wrong, and this is a major issue. The timing of extension: in the model a uniform timing for rifting is assumed between 18-10. Most structures are inferred to be active only until the Middle Miocene (e.g. Majcin et al. 2015), a few small-offset normal faults would not have influenced lithospheric thinning.
- The model result: this is already a mixture of discussion and describing some results. How is it possible that nearly 0 crustal thinning is calculated for the Rechnitz core complex area, that must have undergone substantial crustal thinning? This is a sign of the wrong model parameters. In Figure 6, it is not possible to read the values of the mismatch between the well and model data, but it still seems to be a significant error. About sediment blanketing: in the results of the shallow temperature field chapter you write: “Positive anomalies are the reflection of sediment blanketing, meaning the insulating effect of sediments with low thermal conductivity.” – The deposition of cold sediment would lead to decreased temperature values at shallow depth and higher thermal values in the basement because of the blanketing of low conductivity sediments.
- How did you consider the uplift of the basin margins linked to the ongoing inversion of the basin (e.g. Bada et al. 2007)? Likely it would have a major impact.
- Structure of the manuscript: The results and their discussion are not separated. You should make clear which parameters and which model outputs are well constrained and what is the sensitivity of others.
- Comparison with previous studies: this manuscript doesn’t even mention previous modelling efforts on the crustal and mantle thinning, surface heat flow and basin temperature evolution. In the detailed comments below, you find many useful papers that can be used to compare your results with previous inferences. Besides well data, vitrinite information is also widely available in the region that should be used to validate such models.
- Because of the many limitations listed above, the final sentences on the new stress envelope or comparison with xenoliths remain elusive and in general they don’t really connect with the manuscript. Instead, you should discuss the sensitivity and reliability of the thermal model and compare it with previous inferences and with other similar regions.
further detailed comments:
Title: it does not reflect the content of the manuscript
Abstract: reliable thermal evolution is not modelled for the entire lithosphere due to the limitations of the modelling approach and lack of constraints
ln 19-20: not all the sedimentary basins are extensional
ln 22: Royden and Keen 1980
ln 32: most thermo-mechanical models are constrained by observations, e.g. Lescoutre et al. 2019; Heckenbach et al. 2021; many others
ln 37: in the upper crust
ln 44: sometimes you include Late Miocene, in other places you write Early to Middle Miocene. Which is true?
ln 47: how is this inversion stage considered in the model?
ln 53: i.e. compositional changes through sedimentation: what does this mean?
ln 54-56: you should reflect on the large number of previous thermal modelling efforts in the region, including, but not limited to: Lankreijer et al. 1999; Majcin et al. 2015; Bartha et al. 2018; Balasz et al. 2021; Rybar and Kotulova 2023
ln 59: "high precision" - can you elaborate?
ln 59: for (not to)
ln 73: lower plate with respect to what? Out of context.
fig. 2: what is the sedimentary basin on the right side? Also indicate the orientation of the section.
ln 87: justification?
ln 91: delete -
ln 95: where is this thickness map presented, shown?
ln 100-101: justification?
ln 118-120: rephrase
table 1: what about paleogene rocks?
table 1: what is the source of information behind this data?
ln 128: is it available or the most important data now made available with this manuscript?
ln 135: meters?
ln 136: in fact I cannot see too many wells in the deep basins. Elaborate
Figure 3: scale of the basement depth map?
ln 161: what about different amounts of upper and lower crustal thinning, likely affecting the Rechnitz region?
ln173: sensitivity of this assumption? What if the initial lithopsheric thickness was lower or higher?
ln 174-175: In this model, when the lithosphere was thinned to ca. 60 km, you had a 60 km depth domain of constant temperature beneath? How reliable is this? Why dont you use a constant heat flow lower boundary condition?
Ln 180: so your model is only accurate until 5-10 km depth?
Ln 190: instead of this, it would be more useful to write about the thinning factors of the study area
Ln 196: grammar
Ln 198-200: what is the limitation of this?
Table 2: Lab: 120 meters?
Ln 204: There are many other studies calculating different crustal and mantle thinning, e.g.: Lankreijer et al. 1995; 1999; Majcin et al. 2015; Bartha et al. 2018; Balasz et al. 2021; Rybar and Kotulova 2023
Ln 205: Primary?
Ln 206: what does past-extension mean?
Ln 209: why 35 km?
Fig. 5: how would you discuss these patterns?
Ln 259-263: this is discussion, not results
Ln 275: I would respectfully challenge this statement. How can you be sure that the deposition of cold sediments would increase the temperature in such shallow depth? It would increase at larger depth. Of course, you have higher temperature values, where the crust is thinner and therefore the mantle is more elevated.
Figure 7: which wells are these, what is the source of information? Is it open-source? At least the used and presented well data should be better documented and shared with this manuscript. It is also a warning sign how the errors increase with depth which questions the reliability of the models.
Ln 294: ref
Fig. 10: on the well data the basin was much shallower, which is right? Furthermore, it is not likely that the crust would be laterally homogenous, therefore it is difficult to understand the value of this cross-section.
Citation: https://doi.org/10.5194/egusphere-2024-308-RC1 -
AC2: 'Reply on RC1', Eszter Békési, 07 Aug 2024
Dear Reviewer,
We would like to express our gratitude for the thorough review of our work. We received many insightful comments that have helped us to significantly improve the manuscript. We revised the model input parameters and performed a sensitivity analysis to a selection of model input parameters. The temperature predictions have slightly changed in the revised models, together with the estimated amount of lithosphere extension, due to the revised input parameters we applied. Please find our detailed point-by-point responses to the comments in the attachment.
Kind regards,
Eszter Békési and co-authors
-
EC1: 'Reply on RC1', Patrice Rey, 03 Sep 2024
Dear colleague,
Following your comments and advice, Eszter Békési and co-authors have now submitted their revised version.
It looks to me that they have done a reasonable job, nevertheless, could you please have a look at the revised version and let us know your expert opinion?
Kind regards
PatriceCitation: https://doi.org/10.5194/egusphere-2024-308-EC1
-
CC2: 'Comment on egusphere-2024-308', Nicolas Coltice, 12 May 2024
The manuscript present a thermal reconstruction study of the Hungarian area of the Pannonian basin. It describes a new inverse methodology in order to obtain tectonic information on lithosphere thinning in the area. First of all, I state here that I am more a specialist of modelling than on the tectonics and geothermics of this area.
My point of view is the qualities of the manuscript lie in:
- the new methodology employed to get information of the deep lithospheric structure from temperature measurements.
- pushing the result towards interpretations on rheology and xenolith depth origin.My opinion is that the shortcomings of this manuscript are:
- it is difficult to estimate if the method is able to improve the knowledge of the deep lithospheric structure, especially in a hot and thin crust area in which hydrothermalism and deformation/melting are present. Before inversion, the prior has already a very small misfit (1.33°C). I guess that the uncertainties on the depth of the different layers can introduce such misfit on its own (the thermal gradient is around 40°C/km). The misfit is improved through the process (0.43°C), but is it significant? Since we don’t have here an analysis of how varying the properties of rocks and depth of interface within uncertainties impact the mist, it is hard to know if the authors can resolve the deep lithospheric thermal structure. Given the low value of the misfit prior to inversion, I would say no.
- the method is not a data assimilation method. Data assimilation, which is mostly used for chaotic models with butterfly effect, means that there are new data than can be assimilated (correction of the model) in time. Here, the observations are present-day only. So it is a classical inversion problem with a new methodology. This is a detail but it is worth to use the proper terms.
- most of the figures/captions require additional information.Minor details along the text:
- Fig.1: explain more clearly why the country borders are used for the study (it is stated later in the text but it would be good to have it here)
- Fig.2: orientation is missing (NW - SE)
- line 100: how is the LAB defined here? The study is thermal, so it would be good to explain.
- Table 1: the table is not very informative. Is it possible to either a graph or more details on how the variations are produced?
- Figure 3: Why gray, green and black circles for the same information? What does the color mean?
- line 146: that would be nice to have more details for the errors. Citing the papers of the first authors does not seem enough to evaluate where they come from.
- line 159: remove statement on the inverse modeling. This is the forward model section and it is fundamental to distinguish the difference between the forward and inverse model.
- line 173: why 120km for the LAB?
- line 178: typo ‘preduction’
- line 204: what is unrealistic? More details are needed here to evaluate how the authors rule out a model.
- equation 3: misplaced parenthesis
- section 3.4: more theoritical details on the inversion method would be good. Why this one and not another one? Where does equation (4) come from and why is it adapted to the problem?
- line 248: explain what a variogram is? Provide a figure?
- Figure 6: large errors in the hottest spots. Explain please.
- Figure 7: the choice of colors make it difficult to read (black and blue lines especially)
- line 294: ref missing
- Figure 8: same for colors
- Figure 10: what are the units?Citation: https://doi.org/10.5194/egusphere-2024-308-CC2 -
AC4: 'Reply on CC2', Eszter Békési, 07 Aug 2024
Dear Nicolas Coltice,
Thank you for providing useful suggestions for the manuscript. We agree that using the correct terms is important, so we revised the text to use the term inversion instead of data assimilation. We agree that the reported uncertainties resulting from the inversion procedure are not sufficient to represent the overall uncertainties of the model. In the revised ms., we show multiple models to assess the effect changes in input parameters to modelling results and discuss further the influence of assumptions and simplifications we made on the temperature estimates. Please find our detailed point-by-point responses to the comments below.
Kind regards,
Eszter Békési and co-authors
-
AC4: 'Reply on CC2', Eszter Békési, 07 Aug 2024
-
RC2: 'Comment on egusphere-2024-308', László Lenkey, 15 Jun 2024
Dear Eszter and coauthors,
The manuscript about modeling transient thermal processes in the lithosphere in the NW part of Hungary presents a method to assess deep lithosphere temperatures. The transient conductive heat transport equation is solved, and the calculated temperatures are fitted to observed temperature data to constrain the subcrustal stretching factor. The transient model considers the two most relevant processes, which were active during the evolution of the area: lithospheric stretching and sedimentation. The modeled lithosphere temperatures are used to deduce rheological inferences and estimate the depth of origin of mantle xenoliths found in the region. It is a valuable manuscript, but I suggest modifications and clarifications before publication.
The past and present temperatures are calculated by solving the transient conductive heat transport equation (Eq. 2). The calculated temperature depends on the initial conditions, the thermal parameters and the vertical velocity vz. You fix these quantities except the vertical velocities related to stretching, thus the results are valid for this specific model. Other choice of the quantities probably would result different stretching factors and different lithosphere temperatures. In the following I will discuss the effects of thermal conductivity of sediments and the sedimentation rate.
The thermal conductivity (TC) of sediments in the model ranges from 1.2 W/mK to 2 W/mK (Table 1). These values are lower than the ones measured on clastic sedimentary rocks from Hungary (ranging from 1.5 W/mK to 5 W/mK in Dövényi and Horváth, 1988; Dövényi et al., 1983, shales: 2.3 ± 0.56 W/mK, sandstones: 3.75 ± 0.71 W/mK summarized in Mihályka et al. 2024). Based on the measured TC values Dövényi and Horváth (1988) established TC-depth trends for shales and sandstones, which result in higher TC than used in the study. The TC of carbonates is also lower than the measured ones. For conductivity of carbonates in Hungary see Dövényi et al. (1983). (limestones: 2.7 W/mK, dolomites: 4.4 W/mK).
As the model temperatures are fitted to the observed values, we would expect lower heat flow using the model TC’s compared to the observed heat flow. You provided the model results and the thermal parameters in an asset for reviewers, and I calculated the model heat flow in the depth interval 1000-1200 m (Fig R1). The modeled heat flow is uniform in the area: 70 ± 3 mW/m2. Except the Zala basin and part of the Danube basin, where the observed values are 90 mW/m2 and 85 mW/m2, respectively, the modeled heat flow is close the to the observed one (disregarding also the Transdanubian range, where groundwater flow occurs). Better fit to the heat flow could be achieved by varying the TC, heat production and subcrustal stretching factor, but it was not the purpose of the study as you mentioned in the manuscript.
In Eq. 2 constant sedimentation rate is applied. Over the center of the basins, where sedimentation took place in Quaternary, the constant sedimentation rate is a good approximation. In the peripheral parts of the basins erosion has been taken place since Pliocene due to basin inversion and uplift. (See e.g. the seismic sections published in Szafián et al. (1999), a paper you refer.) Erosion increases the subsurface temperature thus, in the peripheral parts less lithospheric stretching is required to obtain fit to the observed temperatures.
As it is demonstrated varying the TC of sediments and sedimentation/erosion rate would change the modeled T and heat flow. The question arises: how much is the uncertainty of the derived stretching factors and the calculated deep lithosphere temperature?
Please make an estimate of the uncertainty of the stretching factor and the lithosphere temperature.
My questions and notes related to the text are the following.
Lines 107-108The sediment bulk thermal conductivities were finally obtained using the geometric mean of the bulk matrix conductivities and the thermal conductivity of the pore fluid.
How much was the porosity?
Line 148 outflow temperatures were marked by uncertainties of ±5 ºC,
Do you mean outflow temperature as temperature at the well head? Well head T is unreliable, the error can be much more than ±5 C.
Line 184 What kind of numerical method did you use to integrate Eq. 2?
The last term in the equation is a partial derivation.
Line 211 Eq.3 is mistyped. It is correctly: (Zcrust init- Zbasement )/ZMoho present.
Line 259 temperatures up to 170 ºC, meaning a geothermal gradient of ~45 ºC/km
gradT=(170-12)/4=39.5 C/km
Lines 262-263 Negative anomalies can be attributed to outcropping/near-surface basement rocks (mostly carbonates) having significantly higher thermal conductivities
Near Sopron and Rechnitz the lower temperature is partly caused by lower lithospheric stretching relative to the basin areas.
Lines 264-265 The conductive assumption is although not fully valid for parts of the Transdanubian range built up by fractured and karstified carbonate rocks.
A larger area than the Transdanubian rage is affected by groundwater flow and convective heat transport, because groundwater flow also occurs in the carbonate rocks covered by sediments.
Line 270, Fig. 6. The grey color code to visualize the difference between model and observed temperatures is not suitable to quantify the difference.
Lines 322-324 Towards the Transdanubian Range (Balaton Highland), predicted model temperatures are slightly higher in the deeper part of the model compared to the NW part (Sopron Mts.). This might be explained by the shift in the timing of active rifting, that migrated from NW towards SE (e.g. Balázs et al., 2016).
The rifting time was fixed in the model, so the temperature difference has a different reason, e.g. different subcrustal stretching factors.
Line 329 we compared the overall misfit between modelled and observed temperatures
How did you calculate the difference between the modeled and observed temperatures as they belong to different depths?
Lines 349-351 It must be noted that the predicted subcrustal stretching might not be entirely correct due to changes in the timing of stretching throughout the study area but provide a realistic picture for the degree of lithosphere attenuation.
It is not only the timing of stretching, which influences the stretching factors, but all parameters used in Eq. 2.
Lines 366-368 These differences in shallow temperature predictions can partly be explained by the different calibration datasets used by Lenkey et al. (2017) and (Lenkey et al., 2021), excluding temperature measurements from (recent) geothermal wells documented in the OGRE database.
OGRe (2020) is a very useful database to get quick-look temperature data. However, no information is given about the conditions of the measurement. E.g. it is not known if the BHT value is corrected or not. In the Geothermal Database of Hungary (Dövényi, 1994, Lenkey et al., 2021) the observed data are corrected if possible, and every temperature data is quality checked, and depending on the type and conditions of the measurement they are ranked into a quality categories.
Lines 398-399 give references to Porkoláb et al.
Citation: https://doi.org/10.5194/egusphere-2024-308-RC2 -
AC1: 'Reply on RC2', Eszter Békési, 07 Aug 2024
Dear László,
Thank you for reviewing our manuscript, providing constructive comments as an expert both in thermal modelling and in the study area. We incorporated most of the suggestions, which significantly improved the models and the whole manuscript. Most importantly, the revised thermal conductivities we applied provide more realistic input for the models. We also tested the effect of selected input parameters on the model and provided an estimate on the uncertainties, including both a quantitative (resulting from the inversion) and qualitative (discussion on the effect of model assumptions and fixed input parameters) assessment. The temperature predictions have slightly changed in the revised models, together with the estimated amount of lithosphere extension, due to the revised input parameters we applied. Please find our detailed point-by-point responses to the comments in the attached file.
Kind regards,
Eszter and co-authors
-
AC1: 'Reply on RC2', Eszter Békési, 07 Aug 2024
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
395 | 109 | 39 | 543 | 12 | 13 |
- HTML: 395
- PDF: 109
- XML: 39
- Total: 543
- BibTeX: 12
- EndNote: 13
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1