华北平原新生界盖层地温梯度图及其简要说明

本文报道新编比例尺为1:1500000的华北平原新生界盖层地温梯度图。该图以近4000口钻井的温度资料和对地温场控制因素的分析为基础,并结合地温场数学模拟计算结果编制而成。圈定全区地温梯度G>4℃/100m及大地热流q>62mw/m²的局部地热异常区44片,总面积为25000km²,为地热能勘探、开发远景规划提供了重要的科学依据。     

Thermal structure of lithosphere in Central Asian and Pacific belts and their adjacent cratons,from data of geoscience transects

We investigate the 2D thermal structure of lithosphere in the Central Asian and Pacific tectonic belts and adjacent cratonic areas of Siberia and North China using a synthesis of geothermal data from six geoscience transects covered by seismic, resistivity, and gravity surveys. The patterns of rock density, radiogenic heat production derived from U and Th abundances, thermal conductivity, temperatures, and respective heat flows reveal a layered structure. The model with layers distinguished according to density and thermal parameters includes well pronounced dome-shaped features in the crust which correlate with upwarps of the asthenospheric top. The domes are marked by high heat flows of 60–90 mW/m² with a mantle component higher than the crustal one (30–60 mW/m² against 20–30 mW/m²) and temperatures as high as 800–1100 °C at the Moho. Many of these features correspond to known and potential petroleum basins.     

Heat flow in the rifted continental margin of the South China Sea near Taiwan and its tectonic implications

Temperature measurements carried out on 9 hydrocarbon exploration boreholes together with Bottom Simulating Reflectors (BSRs) from reflection seismic images are used in this study to derive geothermal gradients and heat flows in the northern margin of the South China Sea near Taiwan. The method of Horner plot is applied to obtain true formation temperatures from measured borehole temperatures, which are disturbed by drilling processes. Sub-seafloor depths of BSRs are used to calculate sub-bottom temperatures using theoretical pressure/temperature phase boundary that marks the base of gas hydrate stability zone. Our results show that the geothermal gradients and heat flows in the study area range from 28 to 128 °C/km and 40 to 159 mW/m², respectively. There is a marked difference in geothermal gradients and heat flow beneath the shelf and slope regions. It is cooler beneath the shelf with an average geothermal gradient of 34.5 °C/km, and 62.7 mW/m² heat flow. The continental slope shows a higher average geothermal gradient of 56.4 °C/km, and 70.9 mW/m² heat flow. Lower heat flow on the shelf is most likely caused by thicker sediments that have accumulated there compared to the sediment thickness beneath the slope. In addition, the continental crust is highly extended beneath the continental slope, yielding higher heat flow in this region. A half graben exists beneath the continental slope with a north-dipping graben-bounding fault. A high heat-flow anomaly coincides at the location of this graben-bounding fault at the Jiulong Ridge, indicating vigorous vertical fluid convection which may take place along this fault.     

Numerical simulations to assess thermal potential at Tauranga low-temperature geothermal system,New Zealand

Tauranga low-temperature geothermal system (New Zealand) has been used for the last 40 years for direct uses including space heating, bathing and greenhouses. Warm-water springs in the area are between 22 and 39 °C, with well temperatures up to 67 °C at 750 m depth. A heat and fluid flow model of the system is used to determine reservoir properties and assess thermal potential. The model covers 130 km by 70 km to 2 km depth, and was calibrated against temperatures measured in 17 wells. Modelling shows that to maintain the observed primarily conductive heat flow regime, bulk permeability is ≤2.5?×?10^?14 m² in sedimentary cover and ≤1?×?10^?16 m² in the underlying volcanic rocks. The preferred model (R ²?=?0.9) corresponds to thermal conductivities of 1.25 and 1.8 W/m² for sedimentary and volcanic rocks, respectively, and maximum heat flux of 350 mW/m². The total surface heat flow is 258 MW over 2,200 km². Heat flux is highest under Tauranga City, which may be related to inferred geology. Model simulations give insights into rock properties and the dynamics of heat flow in this low-temperature geothermal system, and provide a basis to estimate the effects of extracting hot fluid.     

Spatial characteristics of groundwater temperature in the Ishikari Lowland, Hokkaido, northern Japan: analytical and numerical applications

In porous sediments of the Ishikari Lowland, there is a gradual increase in the background geothermal gradient from the Ishikari River (3–4 °C 100 m^–1) to the southwest highland area (10 °C 100 m^–1). However, the geothermal gradient at shallow depths differs in detail from the background distribution. In spite of convective heat-flow loss generally associated with groundwater flow, heat flow remains high (100 mW m^–2) in the recharge area in the southwestern part of the Ishikari basin, which is part of an active geothermal field. In the northeastern part of the lowland, heat flow locally reaches 140 mW m^–2, probably due to upward water flow from the deep geothermal field. Between the two areas the heat flow is much lower. To examine the role of hydraulic flow in the distortion of the isotherms in this area, thermal gradient vs. temperature analyses were made, and they helped to define the major components of the groundwater-flow system of the region. Two-dimensional simulation modeling aided in understanding not only the cause of horizontal heat-flow variations in this field but also the contrast between thermal properties of shallow and deep groundwater reservoirs. Electronic Publication     

Application of spectral analysis technique on ground magnetic data to calculate the Curie depth point of the eastern shore of the Gulf of Suez,Egypt

The bottom of the magnetized crust determined from the spectral analysis of magnetic anomaly is interpreted as a level of the Curie point isotherm. A spectral analysis technique was used to estimate the depth of the magnetic anomalies sources (Curie point depth analysis) of the eastern shore of the Gulf of Suez, Egypt. The depth to the tops and centers of the magnetic anomalies are calculated by azimuthally averaged power spectrum method for the whole area. The results obtained suggests from this study showed that the average depth to the top of the crustal block ranges between 1.15 and 1.9 km, whereas the average depth to the center of the deepest crustal block ranges between 9.1 and 12.7 km. Curie point depths in the study area range between 14.5 km in the northwestern part of the study area and 26 km in the southeastern part of the study area. The results imply a high geothermal gradient (34.7 °C/km) and corresponding high heat flow value (72.87 mW/m²) in the northwestern part of the study area. The southeastern part of the study area displays a low geothermal gradient (24.26 °C/km) and low heat flow value (50.9 mW/m²). These results are consistent with the existence of the possible promising geothermal reservoir in the eastern shore of the Gulf of Suez especially at Hammam Faraun area.     

Shallow geothermal regime in the Perth Metropolitan Area

Exploration of Perth's geothermal potential has been performed by the Western Australian Geothermal Centre of Excellence (WAGCoE). Detailed vertical temperature and gamma ray logging of 17 Western Australia Department of Water's (DoW) Artesian Monitoring (AM) wells was completed throughout the Perth Metropolitan Area (PMA). In addition, temperature logs from 53 DoW AM wells measured in the 1980s were digitised into LAS format. The logged data are available in the WAGCoE Data Catalogue.

Analysis of the gamma ray logs yielded the first estimates of radiogenic heat production in Perth Basin formations. Values by formation ranged between 0.24 and 1.065 μW m^?3. The temperature logs provide a picture of true formation temperatures within shallow sediments in the Perth Basin. A three-dimensional model of the temperature distribution was used to produce maps of temperature at depth and on the top of the Yarragadee aquifer.

The temperature data were interpreted with a one-dimensional conductive heat model. Significant differences between the model and the observations was indicative of heat moving via non-conductive mechanisms, such as advection or convection. Evidence of non-conductive or advective heat flow is demonstrated in most formations in the region, with significant effects in the aquifers. Average conductive geothermal gradients range from 13°C km^?1 to 39°C km^?1, with sandstone formations exhibiting average gradients of approximately 25°C km^?1, while insulating silt/shale formations show higher average gradients of over 30°C km^?1.

To produce preliminary heat flow estimates, temperature gradients were combined with thermal conductivities measured elsewhere. The geometric mean heat flow estimates range between 64 mW m^?2 to 91 mW m^?2, with the standard deviation of the arithmetic mean heat flow ranging between 15 and 23 mW m^?2.

The study characterises the shallow temperature regime in the Perth Metropolitan Area, which is of direct relevance towards developing commercial geothermal projects.     

Absence of a regional surface thermal high in the Baikal rift; new insights from detailed contouring of heat flow anomalies

Heat flow in active tectonic zones as the Baikal rift is a crucial parameter for evaluating deep anomalous structures and lithosphere evolution. Based on the interpretation of the existing datasets, the Baikal rift has been characterized in the past by either high heat flow, or moderately elevated heat flow, or even lacking a surface heat flow anomaly. We made an attempt to better constrain the geothermal picture by a detailed offshore contouring survey of known anomalies, and to estimate the importance of observed heat flow anomalies within the regional surface heat output. A total of about 200 new and close-spaced heat flow measurements were obtained in several selected study areas in the North Baikal Basin. With an outrigged and a violin-bow designed thermoprobe of 2–3-m length, both the sediment temperature and thermal conductivity were measured. The new data show at all investigated sites that the large heat flow highs are limited to local heat flow anomalies. The maximum measured heat flow reaches values of 300–35000 mW/m², but the extent of the anomalies is not larger than 2 to 4 km in diameter. Aside of these local anomalies, heat flow variations are restricted to near background values of 50–70 mW/m², except in the uplifted Academician zone. The extent of the local anomalies excludes a conductive source, and therefore heat transport by fluids must be considered. In a conceptual model where all bottom floor heat flow anomalies are the result of upflowing fluids along a conduit, an extra heat output of 20 MW (including advection) is estimated for all known anomalies in the North Baikal Basin. Relative to a basal heat flow of 55–65 mW/m², these estimations suggest an extra heat output in the northern Lake Baikal of only 5%, corresponding to a regional heat flow increase of 3 mW/m². The source of this heat can be fully attributed to a regional heat redistribution by topographically driven ground water flow. Thus, the surface heat flow is not expected to bear a signal of deeper lithospheric thermal anomalies that can be separated from heat flow typical for orogenically altered crust (40–70 mW/m²). The new insights on the geothermal signature in the Baikal rift once more show that continental rifting is not by default characterized by high heat flow.     

中国地下热水分布之特点及属性

本文依据我国近年来地下热水资源调研和勘探的新进展及其与有关的研究结果，简述我国地下热水形成和赋存的地热地质背景，分析和归纳构造隆起区中的温泉和中、新生代沉积盆地中的地下热水分布之基本特点，讨论主要水热带的地热学属性和沉积盆地热水的形成机制，为我国地热资源的勘探、开发和今后地球科学有关问题的深化研究提供参考。     

渤海地温场特点的初步研究

渤海位于华北断块区的东部,是华北古生代地台之上发育起来的中、新代断陷盆地,是陆上埕宁隆起和辽河、黄骅及济阳三大拗陷的自然延伸地带。渤海海域的地热特点,一直为人们所关注。在研究渤海的形成和演化时,有人对渤中凹陷的地热状况作了推论;有人根据少量的地热资料认为渤海及其邻近地区是一个高热流异常区。我们根据渤海海域的地温数据,和近年对华北平原地温场的研究结果进行对比,着重对本区地壳浅部的地热状况进行分析和讨论,提出一些初步的看法。     

Research on the characteristics and influencing factors of terrestrial heat flow in Guizhou Province

Terrestrial heat flow is an important physical parameter in the study of heat transfer and thermal structure of the earth and it has great significance in the genesis and development and utilization potential of regional geothermal resources. Although several breakthroughs in geothermal exploration have been made in Guizhou Province. The terrestrial heat flow in this area has not been properly measured, restricting the development of geothermal resources in the province. For this reason, the terrestrial heat flow in Guizhou was measured in this study, during which the characteristics of heat flow were determined using borehole thermometry, geothermal monitoring and thermal property testing. Moreover, the influencing factors of the terrestrial heat flow were analyzed. The results show that the thermal conductivity of rocks ranges from 2.0 W/(m·K) to 5.0 W/(m·K), with an average of 3.399 W/(m·K); the heat flow varies from 30.27 mW/m² to 157.55 mW/m², with an average of 65.26 ± 20.93 mW/m², which is slightly higher than that of the average heat flow in entire land area in China. The heat flow in Guizhou generally follows a dumbbell-shaped distribution, with high values present in the east and west and low values occurring in the north and south. The terrestrial heat flow is related to the burial depths of the Moho and Curie surface. The basaltic eruptions in the Emeishan led to a thinner lithosphere, thicker crust and lateral emplacement, which dominated the basic pattern of heat flow distribution in Guizhou. In addition, the dichotomous structure of regional active faults and concealed deep faults jointly control the heat transfer channels and thus influence the terrestrial heat flow.     

Present Geothermal Fields of the Dongpu Sag in the Bohai Bay Basin

The Dongpu sag is located in the south of the Bohai Bay basin, China, and has abundant oil and gas reserves. To date, there has been no systematic documentation of its geothermal fields. This study measured the rock thermal conductivity of 324 cores from 47 wells, and calculated rock thermal conductivity for different formations. The geothermal gradient and terrestrial heat flow were calculated for 192 wells on basis of 892 formation-testing data from 523 wells. The results show that the Dongpu sag is characterized by a medium-temperature geothermal field between stable and active tectonic areas, with an average geothermal gradient of 32.0°C/km and terrestrial heat flow of 65.6 mW/m2. The geothermal fields in the Dongpu sag is significantly controlled by the Changyuan, Yellow River, and Lanliao basement faults. They developed in the Paleogene and the Dongying movement occurred at the Dongying Formation depositional period. The geothermal fields distribution has a similar characteristic to the tectonic framework of the Dongpu sag, namely two subsags, one uplift, one steep slope and one gentle slope. The oil and gas distribution is closely associated with the present geothermal fields. The work may provide constraints for reconstructing the thermal history and modeling source rock maturation evolution in the Dongpu sag.     

Geothermal gradients and heat flow variations in parts of the eastern Niger delta,Nigeria

Geothermal gradients and present day heat flow values were evaluated for about seventy one wells in parts of the eastern Niger delta, using reservoir and corrected bottom–hole temperatures data and other data collected from the wells. The results showed that the geothermal gradients in the shallow/continental sections in the Niger delta vary between 10 - 18° C/km onshore, increasing to about 24° C/km seawards, southwards and eastwards. In the deeper (marine/paralic) section, geothermal gradients vary between 18 - 45° C/km. Heat flow values computed using Petromod 1–D modeling software and calibrated against corrected BHT and reservoir temperatures suggests that heat flow variations in this part of the Niger delta range from 29–55 mW/m² (0.69–1.31 HFU) with an average value of 42.5 mW/m² (1.00 HFU). Heat flow variations in the eastern Niger delta correspond closely to variations in geothermal gradients. Geothermal gradients increase eastwards, northwards and seawards from the coastal swamp. Vertically, thermal gradients in the Niger delta show a continuous and non-linear relationship with depth, increasing with diminishing sand percentages. As sand percentages decrease eastwards and seawards, thermal gradient increases. Lower heat flow values (< 40 mW/m²) occur in the western and north central parts of the study area. Higher heat flow values (40 - 55 mW/m²) occur in the eastern and northwestern parts of the study area. A significant regional trend of eastward increase in heat flow is observed in the area. Other regional heat flow trends includes; an eastwards and westwards increase in heat flow from the central parts of the central swamp and an increase in heat flow from the western parts of the coastal swamp to the shallow offshore. Vertical and lateral variations in thermal gradients and heat flow values in parts of the eastern Niger delta are influenced by certain mechanisms and geological factors which include lithological variations, variations in basement heat flow, temporal changes in thermal gradients and heat flow, related to thicker sedmentary sequence, prior to erosion and evidenced by unconformities, fluid redistribution by migration of fluids and different scales of fluid migration in the sub-surface and overpressures.     

Estimates of heat flow and heat production and a thermal model of the São Francisco craton

An updated analysis of geothermal data from the highland area of eastern Brazil has been carried out and the characteristics of regional variations in geothermal gradients and heat flow examined. The database employed includes results of geothermal measurements at 45 localities. The results indicate that the Salvador craton and the adjacent metamorphic fold belts northeastern parts of the study area are characterized by geothermal gradients in the range of 6–17°C/km. The estimated heat flow values fall in the range of 28–53 mW/m², with low values in the cratonic area relative to the fold belts. On the other hand, the São Francisco craton and the intracratonic São Francisco sedimentary basin in the southwestern parts are characterized by relatively higher gradient values, in the range of 14–42°C/km, with the corresponding heat flow values falling in the range of 36–89 mW/m². Maps of regional variations indicate that high heat flow anomaly in the São Francisco craton is limited to areas of sedimentary cover, to the west of the Espinhaço mountain belt. Crustal thermal models have been developed to examine the implications of the observed intracratonic variations in heat flow. The thermal models take into consideration variation of thermal conductivity with temperature as well as change of radiogenic heat generation with depth. Vertical distributions of seismic velocities were used in obtaining estimates of radiogenic heat production in crustal layers. Crustal temperatures are calculated based on a procedure that makes simultaneous use of the Kirchoff and Generalized Integral Transforms, providing thereby analytical solutions in 2D and 3D geometry. The results point to temperature variations of up to 300°C at the Moho depth, between the northern Salvador and southern São Francisco cratons. There are indications that differences in rheological properties, related to thermal field, are responsible for the contrasting styles of deformation patterns in the adjacent metamorphic fold belts.     

Present Terrestrial Heat Flow Measurements of the Geothermal Fields in the Chagan Sag of the YingenEjinaqi Basin,Inner Mongolia,China

Owing to the lack o f terrestrial heat flow data, studying lithospheric thermal structure and geodynamics of the Yingen-Ejinaqi Basin in Inner Mongolia is limited. In this paper, the terrestrial heat flow o f the Chagan sag in the YingenEjinaqi Basin were calculated by 193 system steady-state temperature measurements of 4 wells, and newly measuring 62 rock thermal conductivity and 20 heat production rate data on basis o f the original 107 rock thermal conductivity and 70 heat production data. The results show that the average thermal conductivity and heat production rate are 2.11 ±0.28 W/(m.K) and2.42±0.25 nW/m~3 in the Lower Cretaceous o f the Chagan sag. The average geothermal gradient from the Lower Suhongtu 2 Formation to the Suhongtu 1 Fonnation is 37.6 °C/km, and that o f the Bayingebi 2 Formation is 27.4 °C/km. Meanwhile, the average terrestrial heat flow in the Chagan sag is 70.6 mW/m~2. On the above results, it is clear that there is an obvious negative correlation between the thermal conductivity o f the stratum and its geothermal gradient. Moreover, it reveals that there is a geothermal state between tectonically stable and active areas. This work may provide geothermal parameters for further research o f lithospheric thermal structure and geodynamics in the Chagan sag.     

Velocity structure and geothermics of the Earth's crust along the European Geotraverse

Crustal isovelocity lines are constructed along the European Geotraverse for the seismic velocities 6.0, 6.4, 7.1 and 7.8 km/s. Using this velocity structure and a correlation between heat generation and seismic velocity for crustal rocks, the contribution of the crust to the surface heat flow density value is calculated. The heat flow density at the Moho varies from 5 to 40 mW/m² from Paleo-Europe in the north to Neo-Europe in the south, while the mantle heat flow density is close to zero beneath the Alps; the temperatures calculated for the Moho are 260°–390°C for Paleo- to Meso-Europe, 420°–520°C for Neo-Europe and 700°C for the mountain-root beneath the Alps.     

2-D Crustal thermal structure along Thuadara-Sindad DSS profile across Narmada-Son lineament,central India

Central India is traversed by a WSW-ENE trending Narmada-Son lineament (NSL) which is characterized by the presence of numerous hot springs, feeder dykes for Deccan Traps and seismicity all along its length. It is divided in two parts by the Barwani-Sukta Fault (BSF). To the west of this fault a graben exists, whereas to the east the basement is uplifted between Narmada North Fault (NNF) and Narmada South Fault (NSF). The present work deals with the 2-D thermal modeling to delineate the crustal thermal structure of the western part of NSL region along the Thuadara-Sindad Deep Seismic Sounding (DSS) profile which runs almost in the N-S direction across the NSL. Numerical results of the model reveal that the conductive surface heat flow value in the region under consideration varies between 45 and 47mW/m². Out of which 23mW/m² is the contribution from the mantle heat flow and the remaining from within the crust. The Curie depth is found to vary between 46 and 47 km and is in close agreement with the earlier reported Curie depth estimated from the analysis of MAGSAT data. The Moho temperature varies between 470 and 500°C. This study suggests that this western part of central Indian region is characterized by low mantle heat flow which in turn makes the lower crust brittle and amenable to the occurrence of deep focused earthquakes such as Satpura (1938) earthquake.     

Secular cooling and crystallization of partially molten Archaean continental crust over 1 Ga

The protracted tectonic and magmatic record of cratons over the Archaean Eon has been classically interpreted in terms of long-lived shallow-dipping subduction or repeated mantle plumes. In this paper, we use the 1D conductive heat equation to model the evolution of the geotherm of a generic felsic-dominated Archaean cratonic nuclei solely considering the secular decay of radioactive isotopes (²³⁸U, ²³⁵U, ²³²Th, and ⁴⁰K), responsible for heat production in the crust. Using a range of plausible parameters for crustal thickness, lithospheric thickness, and surface heat flux, this modelling shows that Archaean crust was characterized by an initially high geothermal gradient at 3.5 Ga, with a Moho temperature close to 900 °C, and that it might have remained partially molten for about one billion years. The existence of a partially molten crust for an extended period of time offers an alternative option to shallow-dipping subduction or repeated mantle plumes for the understanding of the peculiar tectonic evolution of Archaean cratons marked by (i) protracted high-temperature metamorphism and magmatism associated with crustal differentiation, and (ii) widespread deformation characterized by structural domes attributed to the development of crustal-scale gravitational instabilities.     

Vanished diamondiferous cratonic root beneath the Southern Superior province: evidence from diamond inclusions in the Wawa metaconglomerate

We studied diamonds from a 2.697–2.700 Ga Wawa metaconglomerate (Southern Superior craton) and identified mineral inclusions of high-Cr, low-Ca pyrope garnet, low-Ti Mg-chromite, olivine (Fo₉₃), and orthopyroxene (En₉₄). The diamonds have δ¹³C of ?2.5 to ?4.0 ‰ and derive from the spinel-garnet and garnet facia of harzburgite. Geothermobarometry on non-touching, coexisting garnet-olivine and garnet-orthopyroxene pairs constrains the maximum geothermal gradient of 41 mW/m² for the Neoarchean and a minimum lithosphere thickness of 190 km. The depleted harzburgitic paragenesis equilibrated at a relatively cold geotherm suggests the presence of a pre-2.7 Ga diamondiferous cratonic root beneath the northern Wawa terrane or the Opatica terrane of the Southern Superior craton, i.e., beneath terranes identified as sources for the metaconglomerate diamonds. Geophysical surveys, geothermal data, and petrology of mantle xenoliths emplaced in the Proterozoic-Mesozoic trace evolution of the mantle thermal regime and composition from the Archean to present. The root was thinned down to 150 km by the Jurassic, when the mantle was heated to 41–42 mW/m². The diamondiferous root destruction was accompanied by more significant heating and was complete by 1.1 Ga in areas adjacent to the Midcontinent Rift. The geometry of the current high-velocity root and spatial correlations with boundaries of crustal terranes that docked to the nuclei of the Superior protocraton in the Neoarchean suggest that the root destruction in the Southern Superior may have been associated with tectonic erosion, craton amalgamation, and ensuing ingress of asthenospheric fluids.     

Characteristics of the geothermal water in Changbai Mountain volcanic region,northeast of China

Geothermal water is plentiful in Changbai Mountain region, northeastern China, due to the volcanic activities and widespread faults. For the exploration of geothermal resources, this study uses quartz and cation geothermometer to estimate the temperatures of the geothermal reservoir and uses the tubular models to evaluate the thermal gradient. The hydrogeochemical characteristics of the geothermal resources were also evaluated by hydrogeochemical analysis. The results showed that the geothermal reservoir temperatures of the four major thermal springs in Changbai Mountain region range from 72 to 169 °C. The average geothermal reservoir temperatures of Jinjiang hot springs, Changbai hot springs I, Xianrenqiao hot springs, and Changbai hot springs II are 129.25, 169, 89, and 73.67 °C, respectively. The geothermal gradient values of the four major thermal springs have different characteristics. The geothermal gradient values of Jinjiang hot springs and Changbai hot springs I are 4.6 and 3.1 °C/100 m, respectively. The geothermal gradient values of Xianrenqiao thermal springs and Changbai thermal springs II are both lower than 1.5 °C/100 m, with the values of 1.1 and 1.4 °C/100 m. And the geothermal gradients are influenced by Changbai Mountain Tianchi volcano. In addition, the water chemical analyses showed that the geothermal water types are HCO₃-Na with higher concentrations of Na⁺, Cl^?, SO₄ ^2?, TDS, and HCO₃ ^? than the non-thermal waters, which suggested a deep and long water cycle of the thermal water, and therefore a sufficient water-rock interaction.