4.1 Hydrochemical characteristics
Several geothermal wells in the study area were drilled for high temperature geothermal energy, among which the wellhead temperature of the ZK02 geothermal well in Yulin Village of Kangding reached 198 ℃ and the bottom hole temperature reached 208 ℃ at 300 m deep (Tang et al., 2017). However, owing to the influence of the surface water and quaternary sediments, the temperature of the hot spring exposed on the surface was characterized by low temperatures. Table 1 shows that the hot spring water temperature in the Daofu geothermal area was approximately 40.12–74.46 ℃. These temperatures were the natural outfall, as was the case with the data presented below. The average water temperature was 53.07 ℃, the pH value was 6.3 ~ 8.0, which was slightly alkaline. The river water temperature was 6.48 ℃~15.59 ℃. The hot spring water temperature in the Kangding geothermal area ranged from 15.86 ℃ to 74.81 ℃. The average water temperature was 50.72 ℃, the pH was 6.61 ~ 8.97 with fresh or brackish water, and the river water temperature was 8.52 ℃ to 19.88 ℃. The hot spring water temperature in the Moxi hydrothermal area ranged from 40.52 ℃ to 70.02 ℃, with an average temperature of 50.72 ℃ and a pH value of 6.4 ~ 9.0, which was slightly alkaline water.
The hot spring water in the study area contained a variety of elements, and the main cations were Na+ and Ca2+. However, their contents were different. The Na+ and Ca2+ content of the hot springs water in the Daofu hydrothermal zone ranged from 104.51 mg/L to 536.19 mg/L with an average of 234.66 mg/L, and 49.4 mg/L to 130.4 mg/L with an average of 70.50 mg/L, respectively. The Na+ and Ca2+ contents of the hot spring water of the Kangding hydrothermal zone ranged from 155.27mg/L to 371.70mg/L with an average of 266.27mg/L to 5.6mg/L to 238.9mg/L with an average of 79.91mg/L. The Na+ and Ca2+content of the hot spring water of the Moxi hydrothermal zone ranged from 38.90 mg/L to 295.60 mg/L with an average of 173.12 mg/L, and 1.23 mg/L to 117.30 mg/L with an average of 65.10 mg/L, respectively.
According to Shukalev's classification, where anions greater than 20% of milligram equivalents are invxaysolved in naming, the hot springs of XFZ can be divided into six chemical types, namely Na-HCO3, Na·Ca-HCO3, Na-HCO3·Cl, Na·Ca-HCO3·Cl, Na·Ca-HCO3·SO4, and Na-HCO3·SO4·Cl. A Piper diagram (Fig. 3) was drawn using Origin software to analyze the ionic components of the hot spring water in the XFZ. The results have shown that there are substantial differences in the chemical types of the hot spring water in the three hydrothermal areas. The hydrochemical types of the Daofu hydrothermal active area are Na-HCO3 and Na·Ca-HCO3, while the hydrochemical types of Kangding hydrothermal area are mainly Na-HCO3, Na·Ca-HCO3, and Na-HCO3·Cl. The hydrochemical types of the Moxi hydrothermal zone are Na-HCO3, Na·Ca-HCO3, Na·Ca-HCO3·SO4,Na·Ca-HCO3· Cl, and Na-HCO3·SO4·Cl.
4.2 H and O isotope characteristics
The H and O isotopic composition of the hot spring water can strongly reflect the source, recharge, water–rock interaction, and intensity of water alternation of hot spring water. It is an effective indicator for understanding the circulation path of hot spring water. Craig (1961) summarized the numerical variation relationship of δD and δ18O in meteoric precipitation in different regions worldwide and proposed the Global Meteoric Water Line equation δD = δ18O + 10 (Global Meteoric Water Line, GMWL). Since then, researchers have proposed a series of atmospheric precipitation lines for various regions. In this study, the atmospheric precipitation line equation of the eastern margin of the Qinghai–Tibet Plateau proposed by Wang et al. (2012) was selected, that is δD = 8.41δ18O + 16.72 (Local, LMWL) and the global Meteoric Water Line equation: δD = δ18O + 10. Based on the test results, a straight-line diagram of the δD-δ18O relationship of the hot spring water samples in the XFZ was established (Fig. 4).
Figure 4 shows that the H and O isotope values of the hot spring water samples in XFZ are close to those of the GMWL and LMWL. This indicates that the source of the hot spring recharge in this fault zone is meteoric precipitation.
4.3 Hot water recharge elevation
The XFZ is located on the eastern margin of the Qinghai–Tibet Plateau, and the elevation difference changed substantially. The δ18O and δD of the hot springs were negatively correlated with the exposed height (Zhang et al., 2021). The H and O isotope values were affected by the elevation effect and gradually reduced in value with increasing precipitation height. Therefore, this effect can be used to estimate the recharge height of the hot spring (Clark and Fritz, 1997). Given the slight oxygen drift in the hot spring water sample, which leads to the fractionation of oxygen isotopes, it is more effective to use δD as a calculation index in estimating the elevation of the hot spring recharge.
H = (δS − δP)/K + h (1)
where H is the replenishment elevation of the hot spring in units of m, δS is the δD value of the sampling point of the hot spring water (‰), δP is the δD value of precipitation (‰), and K is the δD isotope gradient value of precipitation in the study area (‰/km). H is the altitude of the hot spring water sampling point (m). The δD isotope gradient in The Sichuan–Tibetan–Guizhou region of China is − 26‰/km (Yu et al., 1984), and the δD value of precipitation in the Western Sichuan Plateau is approximately − 90‰. The relationship between the δD value of precipitation and altitude H in Sichuan, Guizhou, and Tibet in China is as follows (Zhu et al., 2008; Zhou et al., 2010):
δD = − 0.026H − 30.2 (2)
The elevation effect formula for the δD of precipitation in China (Zhou et al., 2010) is as follows:
δD = − 0.03HALT − 27 (3)
In the formula, HALT is the recharge elevation for the hot springs.
According to formulas (1), (2), and (3), the recharge elevation of the hot spring on the XFZ was calculated (Table 2).
4.4 Reservoir temperature and circulation depth
The commonly used geothermal temperature scales can be divided into two main categories, namely the cationic temperature scale and SiO2 temperature scale. The cationic geothermal temperature scale is established using the relationship between the ratio of cations in the chemical components of hot spring water and the temperature. The cations are less affected by dilution and boiling effects in the absence of external ion mixing. Giggenbach (1988) proposed the Na–K–Mg triangle diagram method to evaluate the water–rock equilibrium state and temperature. The hot spring water is divided into three zones, namely fully equilibrium water, partially equilibrium water and immature water, and whether it is suitable for cationic geothermal temperature scale is determined. The water–rock reaction diagram of the study area in Fig. 5 shows that most of the water sample points in the study area were located in the immature water area in the lower right corner of the triangle diagram. Only three water sample points were located in the partial equilibrium water area. This indicates that most of the water sample points did not reach the water–rock equilibrium state. This is likely because the hot water mixed with the shallow cold water during the rising process or the groundwater did not fully react with the surrounding rock in thermal storage. Therefore, the cationic geothermal temperature scale was not suitable for estimating the thermal storage temperature in the study area. The hot water temperature is below 110°C, and it is generally chalcedony that controls the silica content in hot spring water (Fournier and Rowe, 1966), calculated according to the method by (Arnórsson, 1983)
T = 1032/[4.969 − lg(CS )] − 273.15 (4)
In Eq. (4), CS represents the mass concentration of SiO2 in water, and the thermal storage temperature of the Daofu section in the study area of this study was calculated to be 103–166°C, the thermal storage temperature of the Kangding section was 70–120°C, and the thermal storage temperature of the Moxi section was 60–120°C. Given that some points in the hot spring water samples from the Kangding section were not sampled, the thermal storage temperature of the Kangding section was supplemented with other studies to become 70–300°C (Bian et al., 2018; Jiang et al., 2022).
The temperature of hot springs is often controlled by the geothermal heat in which they circulate. The depth of hot water circulation can be estimated using Eq. (5):
Z = Tgrad (T − T0 ) + Z0 (5)
Where Z0 is the circulation depth (m), Tgrad is the geothermal warming rate, which is taken as 20.4 m/°C (Zhang et al., 2021; Zhang et al., 2021), T is the thermal storage temperature (°C), T0 is the average annual temperature of the recharge area, which is taken as 7°C, and Z0 is the depth of the normal temperature zone, which is taken as 30 m. The circulation depth of the study area is calculated to be 1.1 km–6.0 km, which is consistent with that of previous studies (Zhang et al., 2017; Li et al., 2020).
4.5 Structural interpretation of XFZ
The XFZ shows a highly striking linear topography shape in the satellite remote sensing images (Fig. 6a). This is generally distributed in an arc direction toward the NE microridges. The remote sensing interpretation marks of this fault zone are clear, with wide valleys and depressions being formed along the fault zone (Fig. 6b) and the shore of the lake being controlled by the fault and having a straight line (Fig. 6c). The ridge and drainage system adjacent to the fault are twisted and deformed.
The XFZ generally shows strong compressive properties, and the shear properties are strengthened from the southeast to the northwest. Therefore, the mechanical characteristics of the XFZ were segmented and divided into three sections. The Daofu section west of Bamei town of Daofu County is extremely straight-line in shape (Fig. 7a) with a long and narrow fault valley. The fault strike is N55°W, the ridge and drainage system on both sides of the main fault show regular “elbow” twist, with sinistral shear properties. This section of the XFZ is composed of an pronounced major fault and lateral dispersion of the small faults. The Kangding section from Bamei town to Kangding city comprised of three diagonally dependent faults (Fig. 7b, c). From SW to NE there are the Zheduotang fault, Selaha fault, and Yala River faults. The fault surface is in the shape of a slow wave, the shores of seismic lakes controlled by these faults are rectilinear, and drainage dislocations have occurred. The overall characteristics of the fault are torsional compression properties, which are primarily torsional. The length of the Moxi section (Kangding–Tianwan) is approximately 40 km, and the fault has a smooth wavy shape along the strike, showing compression and torsion properties and predominantly compressional properties (Fig. 7d). In general, the torsion property of the XFZ increased from southeast to northwest, while the pressure property gradually decreased. This is because the strike of the southeast segment of the fault zone is deflected from NW to NNW, and the increase in lateral pressure from the sinestral movement has led to the obstruction of the east-trending horizontal slip of the southwest wall of the southeast segment of the fault and the conversion to a vertical motion component.