ORIGINAL RESEARCH article

Front. Earth Sci., 08 June 2022

Sec. Solid Earth Geophysics

Volume 10 - 2022 | https://doi.org/10.3389/feart.2022.854232

Lithospheric Thermal Structure in Jinggangshan City: Implications for High Geothermal Background

  • 1. Institute of Hydrogeology and Environmental Geology, Chinese Academy of Geological Sciences, Shijiazhuang, China

  • 2. Technology Innovation Center for Geothermal and Hot Dry Rock Exploration and Development Ministry of Natural Resources, Shijiazhuang, China

Article metrics

View details

5

Citations

2,8k

Views

616

Downloads

Abstract

Jinggangshan City, which is located in the southwestern part of Jiangxi Province, is rich in hot springs. However, lack of geothermal studies has limited the exploration and utilization of the geothermal resources in this city. This study estimated the terrestrial heat flow in Jinggangshan City through well logging and analyses of rock thermal conductivity and heat production. Based on this and the surrounding geoscience transect and the Crust 1.0 model, this study constructed a one-dimensional lithospheric thermal structure by solving the steady-state heat conduction equation. Furthermore, the deep temperature distribution in Jinggangshan City was obtained to better understand the geodynamic condition of the geothermal resources in this city. The results show that the heat flow in this city is around 83.52 mW/m2, indicating the high heat background for the formation of geothermal resources. According to the lithospheric thermal structure of Jinggangshan City, the mantle contributes more to the terrestrial heat flow (qm/qc > 1) than the crust. The temperature of the Mohorovicic discontinuity (the Moho) is 671.7°C, which is consistent with that below the Tanlu deep fault (620–690°C). Moreover, the calculated depth of the Curie surface (585°C) is 27 km, which is consistent with the Curie isotherm depth estimated from aeromagnetic data. This consistency verifies the validity of the lithospheric thermal structure of Jinggangshan City constructed in this study. In summary, the high heat background plays an important role in the formation of geothermal resources in Jinggangshan City.

1 Introduction

Lithospheric thermal structure plays a significant role in the potential assessment and exploration of geothermal energy (Sclater et al., 1980; Davies and Davies, 2010; Liu et al., 2016; Li et al., 2017). It integrates information on crustal and mantle heat flow, ground temperature distribution, and thermal conductivity and heat productivity of rocks (Liu et al., 2017; Zhang et al., 2018). Terrestrial heat flow density (also referred to as heat flow) is a key parameter used to constrain the lithospheric thermal structure (Jaupart and Mareschal, 2015; Jiang et al., 2019; Prol-Ledesma and Moran-Zenteno, 2019). It is the heat flux per unit time per unit area flowing from the interior to the surface of the Earth (unit: mW/m2 (Furlong and Chapman, 2013; Wang et al., 2015) and comprises crustal and mantle heat flow (Jaupart and Mareschal, 2015; Liu et al., 2021).

Jiangxi Province, which lies in the geothermal district of Southeast China, is rich in hot springs ((Zhu and Zhu, 1992). The geothermal resources in this province are predominantly controlled by faults and have middle–low temperatures (Xiao et al., 2018; Li, 2019). Since heat flow distribution is closely related to geological structures, previous heat flow values in Jiangxi Province were discussed based on the tectonic setting (Hu et al., 1992; Li et al., 1992; Jiang et al., 2019). For example, the heat flow in deep structural variation zones has a high average of 71.9 (n = 14) mW/m2, whereas that in uplifts of the upper mantle has a similar average of 70.2 mW/m2 but a relatively wide range of 56.4–98.8 mW/m2, and that in depressions of the upper mantle (only that in Xiushui, northeastern Jiangxi Province has been measured) is 57.8 mW/m2 (Figure 1). Therefore, deep structural variation zones have the highest heat flow background, with an average of 71.9 mW/m2, which is higher than that of continental China (60.4 ± 12.3 mW/m2 (Jiang et al., 2019)). On the other hand, the geothermal resources in Jiangxi Province, which mostly consists of uplifts and mountains, are predominantly controlled by the distribution of the NE–NNE-trending faults (Zhang, 2017; Li, 2019). The southeastern part of Jinggangshan City is intersected by the Huang’ao Fault—a deep NE-trending fault (Figure 2). Although Jinggangshan City in southwestern Jiangxi Province is located in the depression of the upper mantle, its location in an NE deep fault indicates high potential of high heat flow background (Li et al., 1992). Thus, the southeastern part of Jinggangshan City has great potential for the formation of geothermal resources (Zhu, 2007; Yang et al., 2018). Furthermore, Huang’ao Fault in Jinggangshan City belongs to Tanlu (Tancheng-Lujiang) deep fault (>3,000 km), extending NE along eastern China (Lin, 1992). The Tanlu fault, which experienced multistage deformation events since the Mesozoic, is inferred to be a channel for the percolation of melt and fluid and controls the lithosphere evolution in eastern China (Deng et al., 2013). As exhibited in Figure 2, massive Yanshanian granite (Mesozoic) is distributed along the Huangao fault in the southern part of Jinggangshan City (Fan et al., 2014). The study of the lithospheric thermal structure in Jinggangshan City, therefore, shed light on the heat flow anomaly led by Tanlu fault and its impact on the formation of geothermal resource.

FIGURE 1

FIGURE 2

This study focuses on Huang’ao and Qutan villages in the southeastern part of Jinggangshan City. The heat flow in the study area was acquired through well temperature logging and analyses of rock thermal properties and heat generation rates. Furthermore, by combining the geophysical studies of surrounding areas, the crustal thermal structure of Jinggangshan City was analyzed. The results provide insights into the genesis of the geothermal resources and have important implications for geothermal exploration in Jiangxi Province.

2 Geological Setting

Jinggangshan City is located in the middle part of the Caledonian South China orogenic belt bordering the Neoproterozoic Jiangnan Orogenic Belt (Figure 1; Shen et al. (2009)). As shown in Figure 2, Huang’ao and Qutan villages, which are situated in the southeastern part of Jinggangshan City, are intersected by the regional deep Huang’ao fault (Fan et al., 2014; Liu et al., 2018). The Huang’ao deep fault is an NE-trending (∼45°C) compressional fault with a dip angle of 60°. Its width decreases from the southwest (Qutan village) to the northeast (Huang’ao village), as shown in Figure 2. Several types of alteration have extensively developed in this fault zone, including siliconization, schistosity, mylonitization, and brecciation. The Huang’ao deep fault is a geological boundary between the granites on its west side and the metamorphic rocks on its east side. The granites comprise Caledonian and Yanshanian granites (Figure 1). The Caledonian granites are grayish-white, medium-to-fine-grained equigranular biotite granites, containing quartz (∼35 vo1%), plagioclase (∼30 vol%), alkali feldspar (∼25 vo1%), and biotite (∼5 vol%). The Yanshanian granites are grayish-white, fine-grained, equigranular two-mica granites with a mineral composition similar to that of the Caledonian granites. The strata in the study area are dominated by Precambrian and Ordovician metamorphic rocks, which are subjected to strong folding (Zhu, 2007; Yang et al., 2018). The Precambrian strata mainly consist of the Niujiaohe (Є1n), Gaotan (Є2g), and Shuishi (Є3s) formations. The Niujiaohe formation comprises feldspar–quartz sandy slates, metamorphic sandstones, lamellate silicalites, and organic carbon-bearing phyllites. The Gaotan formation comprisesmedium-to-fine-grained sandy slates and carbonaceous phyllites. The Shuishi formation mostly comprises interbeds consisting of metamorphic feldspar–quartz slates and metamorphic sandstones, with some lenticular limestones at the top. Ordovician strata include the Jueshangou, Qixiling, and Shikou formations. The Jueshangou formation comprises metamorphosed (sandy-silt) slates, the Qixiling formation consists of sericite-slate and feldspar–quartz sandstones, and the Shikou formation includes metamorphic slates and sandstones.

3 Methods

3.1 Measurement of the Terrestrial Heat Flow

Instead of direct measurement, the surface terrestrial heat flow was calculated by multiplying the geothermal gradient by the thermal conductivity of corresponding rocks (Jaupart and Mareschal, 2015; Wang et al., 2015).

3.1.1 Temperature Measurement of Well

The temperature of well JGS-ZK02 was measured using the well logging system (PSJ-2) equipped with a well temperature probe PS2512, which has a measurement range of 0–65°C and an accuracy of 0.1°C. The temperature was measured according to the Chinese industrial standard DZ/T0080-2010 Coal Geophysical Logging Specification. To ensure that the well temperature was in a (quasi-) steady state, the well temperature was repeatedly logged at a time interval of 12 h after its deviation was less than 0.5°C following the completion of drilling.

The temperature measurement of wells in Qutan and Huang’ao villages was carried out using WD-016A-embedded sensors and TD-016C data acquisition instruments. The embedded sensors have a temperature measurement range of -50–100°C and an accuracy of 0.1°C at 0–80 and 0.5°C at 50–0°C and 81–125°C. These temperature sensors (diameter: 11 mm) were implanted into a PVC cable (diameter: 12 mm) at an interval of 5 m (distance 0–50 m) or 10 m (>50 m).

3.1.2 Thermal Conductivity

Thermal conductivity (unit: W/(mK)) is a measure of a given material’s ability to transfer heat, representing the transferred heat per unit area at per unit time when the temperature gradient is 1°C per unit length (Wang et al., 2015). The thermal conductivities of the rock samples were measured using thermal conductivity scanning (TCS) in the East China University of Technology, with a measurement range of 0.2–25 W/(mK) and an accuracy of 3%.

Anand et al. (1973) suggested that thermal conductivity would increase with a decrease in temperature. The temperature at the burial depth of rocks before sampling was probably different from room temperature (25°C), at which the thermal conductivity was measured. Therefore, the measured thermal conductivity at room temperature should be corrected according to the in situ temperature.

Regarding exposed rocks or rocks buried at a relatively shallow depth that could be obtained by drilling, the measured temperatures should be corrected according to their in situ temperature using the following empirical formulas (Sass et al., 1992):where T denotes the in situ temperature of rock samples and K (0), K (25), and K(T) denote the thermal conductivity values of rock samples at 0°C, 25°C, and the in situ temperature, respectively.

For the rocks buried at a depth that could not be acquired through drilling, their thermal conductivity was corrected using the following formula (Cermak et al., 1990; Zang et al., 2002; Xu et al., 2011):where k0 is the thermal conductivity at the top surface of each rock stratum; c is 0.0015; z is the depth of each rock stratum; and b is 0.0015, 0.0001, and 0.0001 for the upper, middle, and lower crust, respectively.

3.2 Lithospheric Thermal Structure

3.2.1 Lithospheric Thermal Structure

The crust is generally divided into four layers: a sedimentary layer, and an upper, middle, and lower crust. Its lower boundary—the Moho—is the contact surface between the lower crust and the mantle.

As the primary seismic wave (P wave, Vp) is distinct in each layer due to different geophysical characteristics, the distribution of Vp in the crust provides evidence of the contact surface between adjacent crust layers and crust structure. In the case where no geophysical analysis of an area has been conducted, the model of Crust 1.0 (Laske et al., 2013) is effective in establishing the crust structure in the area (Liu et al., 2021). Crust 1.0 is a global crustal model specified on a 1 × 1 degree grid and can be used for compiling global artificial seismic data. It provides the crustal structure and related parameters of each crust layer, such as Vp, density and boundary depths.

Since no geophysical study has been carried out in Jinggangshan City, Crust 1.0 is an effective tool for establishing the crust structure of the city. Nevertheless, the Crust 1.0 model has low resolution. To better constrain the crust structure, the geophysical information from a surrounding geoscience transect was collected. The geoscience transect from Menyuan in Qinghang Province to Ningde in Fujian Province (perpendicular distance: ∼340 km) is similar to that of Jinggangshan City (Editorial Committee of Geoscience Transect, 1994; Wang et al., 1995). Though this geoscience transect is a little far away from the study area, it provides more accurate geophysical information.

3.2.2 Vertical Distribution of Heat Flow

The surface heat flow (q0), which is calculated according to the thermal gradient and thermal conductivity (as mentioned in Section 3.1), consists of the crustal heat flow (qc) and mantle heat flow (qm).

The crustal heat flow is the sum of the radioactive heat generated by the heat production elements in the crust (U, Th, and K). The layer-stripping method is usually used to quantitatively estimate the radioactive heat generated in each crust layer and the whole crustal heat flow. The calculation formulas are as follows:where qi is the radioactive heat flow in each crust layer; Di is the thickness of each crust layer; Ai is the heat production rate of each crust layer; and qc is the crustal heat flow.

The heat production rate (A) is traditionally calculated based on the U, Th, and K concentrations of representative rock samples using the following formula (Rybach and Buntebarth, 1981; Artemieva et al., 2017):where A is the heat production rate (μW/m3); ρ is rock density (g/cm3); and CU (ppm), CK (ppm), and CTh (ppm) are the concentrations of U, Th, and K, respectively.

The U, Th, and K contents of rock samples were determined by the Changsha Uranium Geology Research Institute. The K content was measured using a flame atomic absorption spectrophotometer, with a lower detection limit of 0.17%, and the U and Th contents were analyzed using the LA-ICP-MS method, with lower detection limits of 0.003 ug/g and 0.8 ug/g, respectively.

However, this method is not suitable for rocks at depths that are difficult to acquire through drilling. The heat production rate of these rocks can be estimated based on the functional relationship between Vp and Ai (Rybach and Buntebarth, 1984):where A is the heat production rate (μW/m3) and Vp is the primary seismic wave (P wave). In this formula, the in situ seismic velocities should be corrected according to laboratory conditions (100 MPa, room temperature) using the method proposed by Rybach and Buntebarth (1984).

The surface heat flow (q0) minus the crustal heat flow (qi) is the mantle heat flow (qm). The mantle heat flow and the qi/qm ratio are important for understanding deep geodynamic processes (Wang et al., 2015; Zhang et al., 2018).

3.2.3 Deep Temperature Distribution

Well logging can be used to directly measure the temperature of shallow strata that can be unearthed by boreholes. By contrast, the temperature of deep strata should be calculated using the following one-dimensional steady-state heat conduction formula (Lachenbruch, 1970; Liu et al., 2021):where Tz (°C) is the temperature at depth z (m); T0 and q are the temperature and heat flow (mW/m2) at the top of each crust layer, respectively; and D, k (W/(mK)), and A are the thickness (km), thermal conductivity (µW·m−3), and heat production rate of each crust layer, respectively.

4 Results and Discussion

4.1 Terrestrial Heat Flow

4.1.1 Temperature Measurement and Geothermal Gradient

Temperature logging has been conducted for wells JGSDW-B006 and JGS-ZK02 in Huang’ao village and wells JGSDW-B002 and JGSDW-B004 in Qutan village (Table 1). The results are shown in Figure 2. Data fitting was completed using the linear least-squares regression method. The slope of the fitting line was numerically equal to the corresponding temperature gradient (Figure 3 and Table 2). Well JGSDW-B006 had the highest thermal gradient of 26°C·km−1, while well JGS-ZK02 had the lowest thermal gradient of 7.3°C·km−1. The geothermal gradients of the two wells in Qutan village were 9.3°C·km−1 and 13.1°C·km−1.

FIGURE 3

TABLE 1

JGSDW-B006JGS-ZK02JGSDW-B002JGSDW-B004
DTDTDTDTDT
018.9032.226033.8036.0041.7
1019.21032.227033.91036.21041.7
2019.32032.3280342036.32041.8
3019.63032.4290343036.33042.0
40204032.430034.14036.54042.2
5020.35032.531034.25036.55042.3
6020.86032.632034.36036.56042.3
7021.57032.633034.47036.77042.4
8020.88032.734034.58036.88042.4
9021.59032.735034.59037.19042.7
10021.610032.836034.610037.010042.7
1102211032.837034.711037.411042.9
12022.312032.938034.812037.612043.1
13022.613032.939034.913037.613043.0
14022.81403340034.914037.814043.0
15023.115033.14103515038.015043.2
16022.816033.142035.116038.216043.2
17023.417033.243035.217038.217043.4
18023.618033.244035.318038.318043.4
19024.119033.345035.419038.519043.5
20033.446035.520038.720043.5
21033.547035.621038.821043.6
22033.548035.722039.022043.8
23033.649035.823039.223043.8
24033.75003624039.124044.0
25033.725039.1
26039.4
27039.5
28039.6
29039.7
30039.8
31040.0
32040.1
33040.2
34040.3

Temperature logging of geothermal wells in the study area.

Note: D denotes the depth of each well (unit: m) and T denotes the measured temperature (unit °C).

TABLE 2

AreaWellD (m)GG (°C·km−1)DRT (m)
Huang’ao CountryJGSDW-B006451.5260–190
Qutan CountryJGS-ZK025007.30–500
JGSDW-B00244613.10–350
JGSDW-B0044509.30–250

Geothermal gradients and their depth ranges of temperature measurement.

Note: D denotes well depth; GG denotes the geothermal gradient of a well; and DRT denotes the depth range of temperature measurement.

4.1.2 Thermal Conductivity

Five rock samples were collected in total to analyze thermal conductivity. These samples consisted of granites, metamorphic sandstones, slates, limestones, and quartz sandstones (Figure 3).

As shown in Table 3, the metamorphic sandstones had the highest thermal conductivity of 2.75–3.53 W/(mK) (average: 3.79 W/(mK)). Compared with metamorphic sandstones, the thermal conductivity of the Caledonian granites had a similar average of 3.42 W/(mK) but a relatively wider range of 2.06–5.15 W/(mK). The Yanshanian granites had the lowest thermal conductivity, with an average of 2.68 W/(mK). The limestones, sandstones, and slates had approximate thermal conductivity, with an average of 2.77 W/(mK), 2.86 W/(mK), and 2.72 W/(mK), respectively.

TABLE 3

SampleEraLithologyAreaLongitudeLatitudeHeat ConductivityDensityAverageN
JGSDZ-B010OrdovicianMeta sandstoneHuang’ao114°1526°31′3.442.783.243
JGSDZ-B029Meta sandstoneLongshi Town114°0626°46′3.533.08
PCS-8Meta sandstoneNashan114°13′26°39′2.753.25
QT-1YanshanianGraniteHuang’ao114°1126°25′2.662.822.684
QT-2GraniteHuang’ao114°1126°25′2.623.5
YY4GraniteHuang’ao114°12′26°28′2.552.57
dz002GraniteNashan114°10′26°38′2.92.7
JGSDZ-B016CaledonianGraniteHuang’ao114°12′26°28′3.572.673.428
dz003-1GraniteLongshi Town114°00′26°39′3.562.66
dz003-2GraniteLongshi Town114°00′26°39′4.742.67
JGSDZ-B033GraniteLongshi Town114°05′26°443.552.7
JGSDZ-B036GraniteLongshi Town114°04′26°43′2.263.11
JGSDZ-B042GraniteLongshi Town114°03′26°45′2.452.62
JGSDZ-C017GraniteLongshi Town114°05′26°40′5.152.68
JGSDZ-C019GraniteLongshi Town114°08′26°41′2.062.97
PCS-4DevonianLimestoneNashan114°14′26°39′2.813.432.774
PCS-6LimestoneNashan114°14′26°39′2.812.77
PCS-7LimestoneNashan114°14′26°39′2.72.6
PCS-1LimestoneNashan114°14′26°39′2.783.28
PCS-3OrdovicianMetamorphic slateNashan114°14′26°39′2.812.942.722
R-5Metamorphic slateReshuizhou114°12′26°26′2.632.56
JGSDZ-B020DevonianSandstoneNashan114°17′26°42′2.482.42.864
JGSDZ-B021SandstoneNashan114°17′26°42′2.412.41
JGSDZ-B025SandstoneNashan114°10′26°36′2.992.67
JGSDZ-B026SandstoneNashan114°10′26°35′3.543.01

Heat conductivity of rocks in the study area.

Note: the units of heat conductivity and density are W/(mK) and g/cm3, respectively; average denotes the average of heat flux; and N denotes the number of heat flux values involved in the calculation of the average.

4.1.3 Surface Heat Flow

The heat flow was calculated based on the thermal gradient and thermal conductivity of the representative rocks (Table 4).

TABLE 4

AreaWellGG (°C·km−1)DRT (m)HC_C (W/(m·K))Heat flow (mW/m2)
QutanJGSDW-B00213.10–3502.6435.108
JGSDW-B0049.30–2502.6324.924
Shijiao CountryJGSDW-B00626.60–1903.1483.52
Guangming CountryJGS-ZK027.30–5003.5325.77

Calculated heat flow of each well in the study area.

Note: GG denotes the geothermal gradient of a well; DRT denotes the depth range of temperature measurement; and HC_C denotes the corrected heat conductivity.

As presented in Figure 3 and Table 4, the strata at valid depths of 0–190 m mainly consisted of metamorphic sandstones and slates, whose corrected thermal conductivity determined using formulas (1) and (2) was 3.8 W/(mK) and 2.73 W/(mK), respectively. By combining their thicknesses, their weighted thermal conductivity was estimated to be 3.57 W/(mK). The surface heat flow of well JGSDW-B006, which equaled the product of the geothermal gradient (26°C·km−1) and the weighted thermal conductivity, was calculated to be 83.52 mW/m2. This surface heat flow is much higher than the average heat flow (60.4 ± 12.3 mW/m2) of continental China. This finding indicates a high geothermal background in the study area, which is genetically associated to the formation of geothermal resources. In addition, the surface heat flow of wells JGSDW-B002, JGSDW-B004, and JGS-ZK02 was all low (24.9–35.1 mW/m2), which was attributable to the disturbance of vertical water movement (Fan et al., 2014; Wang et al., 2015).

The surface heat flow value of well JGSDW-B006 (83.52 mW/m2) is plotted in Figure 1, with the others in Jiangxi province. As shown in Figure 1, Jinggangshan City is located in the depression of the upper mantle where the heat flow background is expected to be low; however, the measured heat flow value (83.52 mW/m2) there ranks second in all the measured values in Jiangxi Province. It is suggested that Huang’ao deep fault (NE trending) plays a key role in the high heat flow background. This deep fault is probably a channel of heat transfer from the heat source at a depth.

4.2 Crustal and Mantle Heat Flow

4.2.1 Boundaries of the Lithospheric Thermal Structure

The geoscience transect from the Menyuan area in Qinghai Province to the Ningde area in Fujian Province is shown in Figure 4. Jinggangshan City is located to the east of the Tanlu (Tancheng-Lujiang) deep fault, where the boundaries of the upper and the middle crusts have blurred due to the strong thrusting process (Editorial Committee of Geoscience Transect C E a E, 1994; Wang et al., 1995). Therefore, the crust was divided into a sedimentary layer and an upper–middle and lower layer. Based on the distribution of the Vp as shown in Figure 4, the depth of the interface between the middle and lower crusts in the study area was determined to be 20.97 km (Table 5). Similarly, the depth of the Moho was 31.97 km, and the lower boundary of the sedimentary layer was 0.62 km.

FIGURE 4

TABLE 5

LayerH (km)D (km)
Sedimentary layer−0.620.62
Upper–middle crust−20.9720.35
Lower crust−31.9711

Geothermal structure in Jinggangshan City.

Note: H is the depth from the bottom of each layer to the ground surface in Jinggangshan City. The ground surface elevation is 330 m according to the logging of well JGSDW-B006. D is the thickness of each layer.

Based on the Crust 1.0 model, the depth of the interface between the middle and lower crusts was 22.06 km and the Moho depth was 32.9 km. By averaging the results calculated based on the Vp distribution and the Crust 1.0 model, the final depth of the interface between the middle and lower crusts was 21.52 km and the final depth of the Moho surface was 32.44 km (Table 7).

4.2.2 Layer-Stripping Calculation of Heat Flow

As mentioned in Section 3.2.2, the layer-stripping method can be used to calculate the crustal and mantle heat flow. The three layers of the crust (i.e., the sedimentary layer, the upper-middle crust, and the lower crust) were further divided into several portions according to the well logging (Figure 3) and the isolines of P-wave velocity (Vp; Figure 4). The thickness (Di) represents the vertical distance between the top and bottom isolines of Vp.

For the sedimentary layer (thickness: 620 m), the well logging of well JGSDW-B006 (depth: 500 m) showed that this layer mainly consisted of slates and metamorphic sandstones (Fan et al., 2014), with the former located at a depth of < 62 m and the latter distributed at a depth of >62 m. Therefore, the sedimentary layer was subdivided into a slate layer and a metamorphic sandstone layer. The heat production rates (Ai) of these two sublayers were obtained based on the U, Th, and K contents of rock samples. Eleven samples were collected using a hammer to analyze the thermal conductivity (Table 6). The average heat production of the metamorphic slates and metamorphic sandstones was 3.63 μW/m3 and 4.17 μW/m3, respectively.

TABLE 6

SampleEraLithologyAreaLongitudeLatitudeUThKHeat fluxAverageN
Dz001-1OrdovicianMeta sandstoneNashan114°07′10.80″26°37′19.10″13.75.61.484.174.171
PCS-3OrdovicianMetamorphic slateNashan114°13′31.04″26°39′21.25″9.1921.25.634.423.632
R-5Metamorphic slateReshuizhou114°12′00.83″26°26′15.98″7.110.92.052.84

Heat production of rocks in the study area.

Note: the units of U, Th, K and heat flux are ppm, ppm, %, and μW/m3, respectively; average denotes the average value of heat flux; and N denotes the number of heat flux values involved in the calculation of the average value.

Regarding the upper–middle crust and the lower crust, Vp-layer was the average velocity of the top and bottom isolines of Vp (Table 7). The heat flow on the top surface of each layer was calculated using formulas (4), (5), and (7). It is noteworthy that the Vp-layer of the middle and lower crust was corrected by multiplying the correction factor of 1.01 according to the laboratory conditions (Rybach and Buntebarth, 1984). The heat flow of the mantle and crust was 45.63 mW/m2 and 37.89 mW/m2, respectively. Therefore, the qc/qm ratio was less than 1 (0.83), suggesting a hot-mantle–cold-crust type of lithospheric thermal structure. In other words, the heat flow mainly originated from the heat contribution of the mantle. Based on Table 7, a conceptual model is constructed for the lithospheric thermal structure in the study area (Figure 5). As shown in this model, the values of heat flow and heat production rate in each layer could be easily identified and their changes between each other. The heat production and heat flow values in the upper–middle crust is obviously higher than those of the lower crust. Their values have a decreasing trend with increasing depth.

TABLE 7

LayerNZDKVp-bVp-layerVp-cAAiqTz
kmkmW/(m·K)m·s−1m·s−1m·s−1μW/m3mW/m2mW/m2
Sedimentary layerCTZ−0.0283.5215
Slate1−0.0620.0422.733.630.1583.3716.3
Meta-sandstone2−0.620.563.85.84.172.3381.0428.4
Upper–middle crust3−8.517.892.7765.95.962.1617.0264.02234.5
4−12.774.262.766.16.056.111.556.6257.40328.3
5−17.234.462.766.26.156.211.255.5651.84416.7
6−21.524.292.756.56.356.410.803.4548.39494.8
Lower crust7−23.251.732.846.86.656.720.420.7247.66524.0
8−25.642.392.866.96.856.920.270.6447.02563.6
9−30.524.882.9176.957.020.221.0545.97641.6
10−32.441.922.937.17.057.120.170.3345.63671.7

Geoscience transect around the study area.

Note: CTZ denotes the constant temperature zone; N denotes the sequence number of the layer; Z denotes the bottom depth of each layer; D denotes the thickness of each layer; K denotes the thermal conductivity of each layer; Vp-b denotes the P-wave velocity at the bottom of each layer; Vp-layer denotes the P-wave velocity of each layer, numerically equal to the average of vp values of the top and bottom of each layer; Vp-c denotes the corrected Vp (correction factor is 1.01) (Rybach and Buntebarth, 1984); A denotes the heat production rate of each layer; Ai denotes the heat production of each layer; q denotes the heat flow at the bottom of each layer; and Tz denotes the temperature at depth Z.

FIGURE 5

4.2.3 Deep Temperature Distribution

According to Fan et al. (2014), the constant temperature zone in Jinggangshan City (depth: 20 m) has a temperature of 15°C. The temperature below the constant temperature zone was thereby calculated using formulas (1), (2), (3), and (8). As shown in Table 7, the crust was divided into 10 layers to estimate the deep temperature distribution. The calculated temperature of the Moho was 671.7°C, which is consistent with the estimated temperature of the Moho below the Tanlu deep fault (620–690°C (He et al., 2006)).

The Curie temperature isotherm is the basal surface of the magnetite crust consisting of magnetitic minerals and corresponds to the temperature at which magnetic minerals become paramagnetic (Xiong et al., 2016; Yu-fei et al., 2021). Ilmenites and pyrrhotites lose their ferromagnetism at a temperature of 300–350°C. The demagnetized temperature of magnetite is 585°C (Ross et al., 2006), while that of Co–Ni–Fe–bearing minerals is up to 760–800°C. The Curie temperature isotherm is an important parameter of the lithospheric thermal structure since it can be used to estimate the temperature at a depth (Mayhew, 1985; Ross et al., 2006; Yang, 2015). It is a supplement to heat flow in constructing the lithospheric thermal structure. The Curie point isotherm in China (Xiong et al., 2016) shows that the Curie surface in the study area is at a depth of 19–30 km. Based on the deep temperature distribution of the lithospheric thermal structure listed in Table 7, the Curie surface depth corresponding to the demagnetized temperature of magnetite (585°C) was estimated to be 27 km using the interpolation method. The consistency of the Curie surface depth further dominates the validity of the lithospheric thermal structure and deep temperature distribution in Jinggangshan City.

5 Conclusion

This study is the first comprehensive investigation of the lithospheric thermal structure in Jinggangshan City, contributing to a better understanding of the genesis and exploitation potential of the geothermal resources therein. Based on the abovementioned analyses, the following conclusions can be drawn:

  • 1) The heat flow of Jinggangshan City is 83.52 mW/m2, which is much higher than the average heat flow (60.4 ± 12.3 mW/m2) of continental China. The high heat flow indicates a high geothermal background in the study area.

  • 2) As estimated from the lithospheric thermal structure of Jinggangshan City, qc and qm are 37.89 mW/m2 and 45.63 mW/m2, respectively. Therefore, the qc/qm ratio is less than 1, indicating that the terrestrial heat flow in Jinggangshan City mainly originates from the mantle.

  • 3) As shown in the deep temperature distribution calculated using the one-dimensional steady-state equation, the temperature of the Moho in Jinggangshan City is 671.7°C, which is consistent with the temperature of the Moho below the Tanlu deep fault. Furthermore, the calculated depth of the Curie temperature isotherm in Jinggangshan City is 27 km, which is consistent with the depth estimated based on aeromagnetic data. The consistency verifies the validity of the lithospheric thermal structure established in this study.

  • 4) The high heat flow background, which mainly originates from the heat supply of the mantle, contributes significantly to the formation of the geothermal resources in Jinggangshan City.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Author contributions

YL: ideas, field investigation, data analysis, picture drawing, creation of models, manuscript writing, and revision; YL: ideas, field investigation, review, and financial support; FL: field investigation, review, and financial support; SW: field investigation; and HD: picture drawing.

Funding

This work was financially supported by the China Geological Survey Project (DD20190128, 20190505, and DD20160190), the Basal Science Research Fund from the Institute of Hydrogeology and Environmental Geology (Grant Nos. SK202008 and SK202104), the S&T Program of Hebei China (No. 20374201D), and the Natural Science Foundation of Hebei Province China (No. D2019330003).

Acknowledgments

The authors would like to thank Xiaoxue Yan, Haonan Gan, Lingxia Liu, Ruoxi, Yuan, Liu Tao, and Ronghua Hu for their assistance in field work.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

  • 1

    AnandJ.SomertonW. H.GomaaE. (1973). Predicting Thermal Conductivities of Formations from Other Known Properties. Soc. Petroleum Eng. J.13, 267273. 10.2118/4171-PA

  • 2

    ArtemievaI. M.ThyboH.JakobsenK.SørensenN. K.NielsenL. S. K. (2017). Heat Production in Granitic Rocks: Global Analysis Based on a New Data Compilation Granite2017. Earth-Science Rev.172, 126. 10.1016/j.earscirev.2017.07.003

  • 3

    CermakV.BodriL.RybachL.BuntebarthG. (1990). Relationship between Seismic Velocity and Heat Production: Comparison of Two Sets of Data and Test of Validity. Earth Planet. Sci. Lett.99 (1), 4857. 10.1016/0012-821X(90)90069-A

  • 4

    DaviesJ. H.DaviesD. R. (2010). Earth's Surface Heat Flux. Solid earth.1 (1), 524. 10.5194/se-1-5-2010

  • 5

    DengY.FanW.ZhangZ.BadalJ. (2013). Geophysical Evidence on Segmentation of the Tancheng-Lujiang Fault and its Implications on the Lithosphere Evolution in East China. J. Asian. Earth. Sci.78 ,263-276. 10.1016/j.jseaes.2012.11.006

  • 6

    Editorial Committee of Geoscience Transect C E a E (1994). Geoscience Transect from Menyuan, Qinghai Province to Ningde, Fujian Province: Seismological press.

  • 7

    FanB.LiuX.YeH.XvL.YuS. (2014). Geothermal Water Survey Report in Qutan, Xiaqi Country, Jinggangshan City, Jiangxi Province. Nanchang: Jiangxi Institute of Geo-Environment Monitoring.

  • 8

    FurlongK. P.ChapmanD. S. (2013). Heat Flow, Heat Generation, and the Thermal State of the Lithosphere. Annu. Rev. Earth Planet. Sci.41 (1), 385410. 10.1146/annurev.earth.031208.100051

  • 9

    HeL.HuS.YangW.WangJ.YangS.YuanY.ChengZ. (2006). Heating Generation Model for the Continental Crust Based on the Main Hole of the Chinese Continental Scientific Drilling Project. Acta Petrol. Sin.22 (11), 7

  • 10

    HuS. B.XiongL. P.WangJ. Y.ShenJ. Y.LiJ. L (1992). The First Group of Heat Flow Data from Jiangxi Province. Chin. Sci. Bull.19, 1791.

  • 11

    JaupartC.MareschalJ. C. (2015). “6.05 - Heat Flow and Thermal Structure of the Lithosphere.” in Treatise on Geophysics. Editors SchubertG. (Amsterdam: Elsevier), 217251

  • 12

    JiangG.HuS.ShiY.ZhangC.WangZ.HuD. (2019). Terrestrial Heat Flow of Continental China: Updated Dataset and Tectonic Implications. Tectonophysics753, 3648. 10.1016/j.tecto.2019.01.006

  • 13

    LachenbruchA. H. (1970). Crustal Temperature and Heat Production: Implications of the Linear Heat-Flow Relation. J. Geophys. Res.75 (17), 32913300. 10.1029/JB075i017p03291

  • 14

    LaskeG.MastersG.MaZ.PasyanosM. E. (2013). Update on Crust1.0 - a 1-Degree Global Model of Earth's Crust. United States: Jolla

  • 15

    LiJ. (2019). A Discussion on the Symbolic Role of Secondary Fault Structure in Geothermal Prospecting in Southern Mountainous Areas - to Take the Longgui Mountain Geothermal Field in Ningdu County, Jiangxi Province as an Example. J. Geol. Hazards Environ. Preserv.030 (4), 85

  • 16

    LiX.ShiX.ZhouW.SunZ.ZhangW.NiuX.et al (1992). Terrestrial Heat Flow in Jiangxi Province. Sci. Geol. Sin.A12, 383

  • 17

    LiZ.ZuoY.QiuN.GaoJ. (2017). Meso-Cenozoic Lithospheric Thermal Structure in the Bohai Bay Basin, Eastern North China Craton. Geosci. Front.8 (5), 977987. 10.1016/j.gsf.2016.09.003

  • 18

    LinZ. (1992). Geological Transects from Menyuan in Qinghang to Ningde in Fujian. Beijing: Seismological Press.

  • 19

    LiuC.ChenL.WangZ. (2018). Geological Characteristics of Yangkeng Tungsten Ore Deposit in Jinggangshan City, Jiangxi Province. Min. Eng.16 (6), 9

  • 20

    LiuF.LangX.LuC.LinW.TongJ.WangG. (2017). Thermophysical Parameters and Lithospheric Thermal Structure in Guide Basin, Northeast Qinghai-Tibet Plateau. Environ. Earth Sci.76 (5), 199. 10.1007/s12665-017-6503-2

  • 21

    LiuF.WangA.WangG.ZhangW.LiaoY.TongJ. (2021). Characteristics of Terrestrial Heat Flow and Lithospheric Thermal Structure in Typical Intermountain Basin at the Eastern Foot of Yanshan Mountain, North China. Front. Earth Sci.9, 758605. 10.3389/feart.2021.758605

  • 22

    LiuQ.ZhangL.ZhangC.HeL. (2016). Lithospheric Thermal Structure of the North China Craton and its Geodynamic Implications. J. Geodyn.102, 139150. 10.1016/j.jog.2016.09.005

  • 23

    MayhewM. A. (1985). Curie Isotherm Surfaces Inferred from High-Altitude Magnetic Anomaly Data. J. Geophys. Res.90 (B3), 26472654. 10.1029/JB090iB03p02647

  • 24

    Prol-LedesmaR. M.Morán-ZentenoD. J. (2019). Heat Flow and Geothermal Provinces in Mexico. Geothermics78 (MAR.), 183200. 10.1016/j.geothermics.2018.12.009

  • 25

    RossH. E.BlakelyR. J.ZobackM. D. (2006). Testing the Use of Aeromagnetic Data for the Determination of Curie Depth in California. Geophysics71 (5), L51L59. 10.1190/1.2335572

  • 26

    RybachL.BuntebarthG. (1984). The Variation of Heat Generation, Density and Seismic Velocity with Rock Type in the Continental Lithosphere. Tectonophysics103 (1-4), 335344. 10.1016/0040-1951(84)90095-7

  • 27

    RybachL.BuntebarthG. (1981). Heat-Generating Radioelements in Granitic Magmas. J. Volcanol. Geotherm. Res.10 (4), 395404. 10.1016/0377-0273(81)90088-3

  • 28

    SassJ. H.LachenbruchA. H.MosesT. H.MorganP. (1992). Heat Flow from a Scientific Research Well at Cajon Pass, California. J. Geophys. Res.97, 50175030. 10.1029/91JB01504

  • 29

    SclaterJ. G.JaupartC.GalsonD. (1980). The Heat Flow through Oceanic and Continental Crust and the Heat Loss of the Earth. Rev. Geophys.18 (1), 269311. 10.1029/RG018i001p00269

  • 30

    ShenW.LingH.ShuL.ZhangF.XiangL. (2009). Sm-Nd Isotopic Compositions of Cambrian-Ordovician Strata at the Jinggangshan Area in Jiangxi Province: Tectonic Implications. Sci. Bull.54 (010), 17501758. 10.1007/s11434-009-0214-3

  • 31

    WangC. Y.LinZ. Y.ChenX. B. (1995). Comprehensive Study of Geophysics on Geoscience Transect from Menyuan, Qinghai Province, to Ningde, Fujian Province, China. Chin. J. Geophys.14, 58

  • 32

    WangJ.PangZ.HuS.HeL.HuangS.QiuN. (2015). Geothermics and its Application. Zealand: Supports open access

  • 33

    XiaoZ.WangJ.HouH. (2018). Analysis of Characteristics and Causes of Geothermal Water in the East of Southern Jiangxi Province. J. East China Univ. Technol.41 (3), 58

  • 34

    XiongS.YangH.DingY.LiZ. (2016). Characteristics of Chinese Continent Curie Point Isotherm. Chin. J. Geophys.59 (10), 14. 10.6038/cjg20161008

  • 35

    XuM.ZhuC.RaoS.HuS. (2011). Difference of Thermal Structure between Eastern Edge of Tibet Plateau and Weasten Sichuan Basin. Chin. J. Geol.46 (1), 203212. 10.1007/s11769-011-0451-7

  • 36

    YangH. (2015). Characteristics of Chinese Continental Curie Point Isotherm. Sichuan, Chengdu: Chengdu University of Technology

  • 37

    YangL.ZhouY.YuY.YangW. (2018). Analysis on Genesis and Hydrochemical Characteristics of Geothermal Water in Qutan Area,Jinggangshan,Jiangxi Province. Sci. Technol. Eng.18 (23), 1318. 10.1007/s12665-021-09577-8

  • 38

    Yu-Fei -X.Ya-Bo -Z.Yuen- D. (2021). Geothermal Structure Revealed by Curie Isothermal Surface under Guangdong Province, China. J. Groundw. Sci. Eng.9, 114. 10.19637/j.cnki.2305-7068.2021.02.003

  • 39

    ZangS.LiuY.NingJ. (2002). Thermal Structure of the Lithosphere in North China. Chin. J. Geophys.45 (1), 5666. 10.1002/cjg2.216

  • 40

    ZhangC.JiangG.ShiY.WangZ.WangY.LiS.Jiax.HuS., (2018). Terrestrial Heat Flow and Crustal Thermal Structure of the Gonghe-Guide Area, Northeastern Qinghai-Tibetan Plateau. Geothermics72, 182192. 10.1016/j.geothermics.2017.11.011

  • 41

    ZhangZ. (2017). Hydrogeochemistry of Nitric Hot Springs in Southern Jiangxi Province. Fuzhou, China: East China University of Technology

  • 42

    ZhuB.ZhuL. (1992). Geochemical Exploration of Geothermal Fields. Bath, UK: Geological Publishing House.10.1016/S0375-6505(97)00012-6

  • 43

    ZhuZ. (2007). Hot-Spring in the Reshuizhou Region of Jiangxi Province. Geol. Mineral Resour. South China03, 63

Summary

Keywords

Jinggangshan City, geothermal gradient, heat flow, lithospheric thermal structure, moho temperature

Citation

Liao Y, Liu Y, Liu F, Wei S and Duan H (2022) Lithospheric Thermal Structure in Jinggangshan City: Implications for High Geothermal Background. Front. Earth Sci. 10:854232. doi: 10.3389/feart.2022.854232

Received

13 January 2022

Accepted

13 April 2022

Published

08 June 2022

Volume

10 - 2022

Edited by

Dawei Hu, Institute of Rock and Soil Mechanics, (CAS), China

Reviewed by

Xianchun Tang, Chinese Academy of Geological Science, China

Funda Bilim, Cumhuriyet University, Turkey

Updates

Copyright

*Correspondence: Yanguang Liu,

This article was submitted to Solid Earth Geophysics, a section of the journal Frontiers in Earth Science

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics