Geochemistry of Geothermal Fluids in the Three Rivers Lateral Collision Zone in Northwest Yunnan, China: Relevance for Tectonic Structure and Seismic Activity

: The Three Rivers Lateral Collision Zone (TRLCZ), situated at the southeastern margin of the Tibetan Plateau, is a crucial frontier where materials from the plateau ﬂow southeastward. This study extensively investigated the hydrochemical characteristics and origin of helium and carbon isotopes in 73 thermal springs within the TRLCZ. The analysis revealed dominant processes, including carbonate and silicate interactions, resulting in elevated concentrations of HCO 3 − and Na + . The impact of Ca/Mg-rich minerals, particularly dolomite, inﬂuenced the cation composition. Additionally, gypsum dissolution, notably in the Lancangjiang Fault and Weixi–Qiaohou Fault, was highlighted through Ca/SO 4 ratios. The positive correlation between SO 42 − and Cl − indicated dilution by shallow cold water, explaining the lower SO 42 − content in the Jingshajiang–Zhongdian Fault and Nujiang Fault compared to the Weixi–Qiaohou Fault and Lancangjiang Fault. The circulation depth of thermal spring water varied, with the northern Weixi–Qiaohou Fault exhibiting the shallowest circulation depth (~3 km), while the Jingshajiang–Zhongdian Fault and southern segments of the Nujiang Fault displayed deeper depths—ranging from 4 to 7 km. A positive correlation between the circulation depth and fault activity was also observed. The Rc/Ra ratios of free gas samples, predominantly indicating crustal origin, varied from 0.01 Ra to 0.53 Ra. Elevated Rc/Ra ratios in the research area suggested potential minor additions of mantle helium through faults and fractures. Crustal limestone was identiﬁed as the primary source of CO 2 -rich samples, with δ 13 C CO2 values ranging from − 1.6‰ to − 7.2‰, while trace amounts of mantle CO 2 were found. The spatial distribution of the H 2 concentration, CO 2 concentration, He concentration, and mantle He proportions in gases indicated that higher values of He concentration and mantle He% always occur near sampling points with deeper circulation depths. However, no similar correlation was observed for H 2 and CO 2 . Most earthquakes of magnitude 5 or greater occurred near the regions with high values of mantle source He release, highlighting the critical role of mantle ﬂuids in the occurrence of earthquakes in the region. In this study, a ﬂuid circulation model was developed to describe the process of ﬂuid (water and gas) circulation migration and earthquake generation in the TRLCZ.


Introduction
Geothermal fluids include hot water and steam that carry dissolved solutes and gases.These fluids traverse fractured conduits, ultimately reaching the surface; they provide important insights into the geological, chemical, and hydrological characteristics of geothermal systems [1][2][3][4].Some research suggests a substantial release of fluids to the surface, including 3 He, H 2 , N 2 , CO 2 , and CH 4 , through major strike-slip fault zones at considerable cutting depths.This phenomenon is particularly notable in regions undergoing extension Minerals 2024, 14, 274 2 of 25 and experiencing significant rates of active crustal deformation [5,6].Notable examples include the Xianshuihe Fault in China [7,8], the North Anatolian Fault in Turkey [9], and the San Andreas Fault in the United States [10].Additionally, deep-source fluid activity has a significant influence on earthquake occurrence [11,12].When these fluids infiltrate a fault zone, they can increase pore pressure, reduce frictional resistance within the fault zone, and promote fault movement, ultimately leading to earthquake generation [13][14][15][16][17][18].Geophysical evidence also indicates that high-pressure fluids from the upper mantle are a direct factor contributing to seismic activity along the San Andreas Fault near Parkfield and Cholame [19].Deep fluids within the fault system are closely linked to shallow fluids, and pressure pulses stemming from the release of deep fluids before earthquakes generate corresponding responses on the Earth's surface [20].Anomalies in mantle-derived helium have been documented to be associated with seismic events, such as the 2011 Tohoku-Oki M 9.0 earthquake [21], the 2016 Kumamoto M 7.3 earthquake in Japan [22], and the Wenchuan M 8.0 earthquake in China [23].Therefore, geothermal fluids are an effective and promising research direction for tectonic activity and deep-seismogenic environment detection.
The convergence and collision process between the Indian Plate and the Eurasian Plate has caused a large amount of material to escape southeastward [24], forming the lateral collision zone of three rivers (namely Nujiang River, Lancang River, and Jinsha River) [25].It is located in the southeastern margin of the Tibetan Plateau, mainly in western Yunnan Province, China (Figure 1a,b).As the forefront of the eastern collision between the Indian Plate and the Eurasian Plate, new tectonic movement in this area is mainly characterised by strong uplift and significant extrusion deformation of the crust and is closely linked to seismic activity [26][27][28].The slip rate of active fractures is a crucial parameter for assessing the potential seismic hazard [29].However, variations and uncertainties persist in the estimated slip rates of these faults in the region since the Quaternary period.The transport of deep-source volatiles ( 3 He, CO 2 ) through faults and fractures whose permeability is controlled by fault activity.These discrepancies stem from the use of different research methods on the same fault zones [30][31][32][33][34][35][36][37][38].The Three Rivers Lateral Collision Zone (TRLCZ) serves as a crucial pathway for the eastward transport of materials from the Tibetan Plateau.Recent seismic tomography results show that high-conductivity and low-velocity zones (LVZ) are widely distributed in the Sichuan Yunnan region [39][40][41][42][43]. Areas exhibiting high conductivity and low velocity within the fault zone are recognised as regions characterised by fault disruptions resulting from impactful sliding or shearing along the fault [44].These areas are typically filled with deep sources of magma and/or metamorphic fluids ( 3 He, CO 2 ).These zones and their associated fluid activity could be crucial for understanding a range of geological phenomena, including tectonic structures and earthquake occurrence [45,46].Previous studies in the TRLCZ mainly focused on the origin and evolution of one or several springs in this geothermal field, which can be used for the mineral exploration and utilisation of geothermal resources [47,48].As an inaccessible area with one of the most complex terrains and landforms in the world, 70% of the hot springs in Yunnan Province are located in the TRLCZ.Therefore, a comprehensive understanding of the groundwater evolution and genesis of geothermal waters is of critical importance.However, few investigations into geothermal fluids over the entire region have been carried out to date.Particularly, a comprehensive dataset comprising chemical and isotopic compositions is lacking, as well as an interpretation of the tectonic structural and seismic processes that underpin the distribution and release of geothermal fluids in the TRLCZ.
The goal of this study is to define the geochemical signatures, evolution, and genesis of geothermal fluids in the TRLCZ, and to analyse the relationships among the deep-source fluid release, tectonic structure, and seismic activity.Our study could act as a reference to better understand the spatial fault-fluid interactions in active geothermal systems, the deep-seismogenic environment in the TRLCZ, the context of deep-source fluid release, and geochemical parameters for use in seismic precursor monitoring.The goal of this study is to define the geochemical signatures, evolution, and genesis of geothermal fluids in the TRLCZ, and to analyse the relationships among the deepsource fluid release, tectonic structure, and seismic activity.Our study could act as a reference to better understand the spatial fault-fluid interactions in active geothermal systems, the deep-seismogenic environment in the TRLCZ, the context of deep-source fluid release, and geochemical parameters for use in seismic precursor monitoring.

Geological Setting
Owing to the lateral collision of the Eurasian Plate and the Indian Plate in the study area, the TRLCZ is characterised by complex geological structures, intense crustal deformation, and frequent seismicity [49].Within the TRLCZ, three rivers (Nujiang River, Lancang River, and Jinsha River) run roughly parallel to each other within a narrow band, with multiple active faults passing through.These faults include the Nujiang Fault (F1), Lancangjiang Fault (F2), Jinshajiang Fault (F3), Zhongdian Fault (F4), and Weixi-Qiaohou Fault (F5) (as shown in Figure 1b).
The Nujiang Fault Zone (NJFZ) is an east-oriented tensile fault that spreads N-S along the Gaoligong Mountain, where associated drag fold and tectonic breccia can be found along the west of the Nujiang River valley [50].The estimated slip rate of the NJFZ ranges from 2 to 3 mm/a based on studies by Tang et al. [32] and Song et al. [33], while inversion results from GPS observations suggest a higher slip rate of 4 to 6 mm/a [34].The

Geological Setting
Owing to the lateral collision of the Eurasian Plate and the Indian Plate in the study area, the TRLCZ is characterised by complex geological structures, intense crustal deformation, and frequent seismicity [49].Within the TRLCZ, three rivers (Nujiang River, Lancang River, and Jinsha River) run roughly parallel to each other within a narrow band, with multiple active faults passing through.These faults include the Nujiang Fault (F1), Lancangjiang Fault (F2), Jinshajiang Fault (F3), Zhongdian Fault (F4), and Weixi-Qiaohou Fault (F5) (as shown in Figure 1b).
The Nujiang Fault Zone (NJFZ) is an east-oriented tensile fault that spreads N-S along the Gaoligong Mountain, where associated drag fold and tectonic breccia can be found along the west of the Nujiang River valley [50].The estimated slip rate of the NJFZ ranges from 2 to 3 mm/a based on studies by Tang et al. [32] and Song et al. [33], while inversion results from GPS observations suggest a higher slip rate of 4 to 6 mm/a [34].The Lancangjiang Fault Zone (LCJFZ), which is the middle branch of the three rivers fault system, exhibits a near-NS trend in its northern part and an 'S'-shaped spreading pattern in the middle section.The southern part of the fault exhibits an NNW trend.The slip rate of the LCJFZ, obtained from the inversion of GPS data, is estimated to be 2.4 ± 1.2 mm/a [31].This fault system is a thrusting zone coupled with regional metamorphism and magma intrusion [51].The Jinshajiang Fault and Zhongdian Fault (referred to as the JZFZ), located at the western edge of the Sichuan-Yunnan rhombic block, follow an SN-trend and rightlateral strike.Magmatic and metamorphic rocks, such as monzonitic granite, alkali basalt, and mylonite, are typically distributed along this zone [52].The slip rate of the JZFZ is approximately 4 mm/a, according to Tang et al. [30] and Wang et al. [35], while Xu et al. [53] suggested a higher right-lateral slip rate of about 5-7 mm/a.The Weixi-Qiaohou Fault Zone (WQFZ), situated in the Lanping Basin, exhibits an NW-trending fault characterised by dextral strike-slip features.Its average movement rate is approximately 1.8-2.4mm/a horizontally and 0.3-0.35mm/a vertically since the late Pleistocene, as reported by Chang et al. [36,37].Li et al. [38] suggested a maximum right-slip rate of 0.89 mm/a based on the inversion of GPS data.Within the Lanping basin, 80% of the folded Precambrian crystalline basement is covered by late-Paleozoic and Mesozoic limestone, sandstone, mudstone, and conglomerates, where the Jurassic and Cretaceous sedimentary strata are most widely outcropped [54].
The northwestern Yunnan region has experienced numerous earthquakes with magnitudes of 5 or higher in the 20th century, indicating a higher probability of strong earthquakes compared to other regions in China.Some notable earthquakes in recent years include the 2013 Zhongdian earthquake (M 5.9), the 2013 Eryuan earthquake (M 5.5), the 2016 Yunlong earthquake (M 5.0), the 2017 Yangbi earthquake (M 5.1), and the 2021 Yangbi earthquake (M 6.4).It is worth noting that there is relatively less moderate to strong earthquake activity in the western LCJFZ and NJFZ compared to the eastern region.
The western Yunnan province is the southern segment of the Yunnan-Tibet geothermal belt, which is a component of the larger Mediterranean-Himalaya geothermal belt.Natural springs are distributed along deep active fractures in the study area, as shown in Figure 1b.The TRLCZ has a landscape of longitudinal range-gorge zones, where the height difference between the mountain summit and the river floor can reach 6000 m.The elevation of the study area is notably higher in the northeast and lower in the southeast, with the northern areas being permanently covered with snow and ice.The recharge sources for spring water predominantly originate from high mountain areas flanking both sides of river valleys.Meteoric water and snowmelt water significantly contribute to these sources [55].The chemical composition of a hot spring is a consequence of its source waters, with major components of groundwater largely determined by the mineralogy of aquifers.Carbonate precipitation is an important geochemical process as travertine is observed in springs [56].Given its distinctive features, the TRLCZ serves as an exceptional resource for studying structural and seismic processes, through the use of insights derived from geochemistry.

Gas Samples
For collecting the bubbling gas, a gas-collection device using the gas drainage collection method was employed.The sampling procedures, outlined in detail by Zhao et al. [59], were followed closely.Gas samples designated for chemical composition analysis were gathered in 500 mL aluminium foil gas sample bags.Simultaneously, gas samples intended for helium and carbon isotope composition analysis were collected in 125 mL sodium glass bottles.Each bottle was filled with 25 mL of thermal spring water and subsequently sealed with crimped rubber stoppers before being stored in an inverted position.Analysis of the 39 gas samples was completed at the Key Laboratory of Petroleum Resources Research, Chinese Academy of Sciences, Lanzhou, within one month of sampling.Gas compositions were determined by gas chromatography (Agilent 7890A, produced by Agilent Technologies, Santa Clara, CA, USA).The analysis accuracies (v/v) were He and H 2 = 5 ppm; CH 4 , Ar, and O 2 = 10 ppm; N 2 = 0.05%; and CO 2 = 0.5%.The 3 He/ 4 He and 4 He/ 20 Ne ratios were measured using a mass spectrometer (VG5400, produced by VG Instruments, Winsford, UK) with relative standard deviations of <5%.Atmospheric samples from Gaolan Hill in southern Lanzhou were used as standards for the He and Ne isotope analyses [60].The results are reported as R/Ra values (Ra being the isotopic ratio of He in the atmosphere, equal to 1.4 × 10 −6 ) [61] and were corrected for atmospheric He contamination using the air-normalized 4 He/ 20 Ne ratio (X-factor) multiplied by the ratio of Bunsen coefficients [62].Detailed methods are provided in the Supplementary Materials [63].The δ 13 C V-PDB (CO 2 ) values were analysed by gas chromatography (Agilent 6890)-pyrolysis isotope ratio mass spectrometry (GC-IRMS; Thermo-Fisher Scientific Delta Plus XP, produced by Thermo Fisher Scientific, Waltham, MA, USA) [64,65].The results are reported in per mil (‰) relative to the PDB scale, with an uncertainty of ±0.3‰.

Results
The physicochemical parameters and the geochemical data of thermal springs in the TRLCZ are presented in Tables S1 and S2, respectively.The 73 cold and thermal springs had temperatures ranging from 11 • C to 70 • C. The pH values and electrical conductivity (EC) were 5.68-8.99 and 135-5070 µs/cm, respectively.The calculated total dissolved solids (TDS) of water varied from 94 mg/L to 35,188 mg/L.The saline spring samples (No.L41) had relatively higher TDS values (35,188 mg/L) than others.
The chemical compositions and isotope compositions of the 39 spring gases sampled from the TRLCZ (Figure 1b) are listed in Table S3.The CO 2 concentration was in the range of 0.01-89.93vol%, and more than half of the samples had concentrations larger than 70 vol%.N 2 had concentrations ranging from 0.97 vol% to 93.13 vol%.The argon concentration ranged between 0.04% and 2.50%.The concentrations of He varied from 4 ppm to 1276 ppm, and samples J32 and L51 had much higher He concentrations of >1000 ppm.The H 2 concentrations ranged from 1 ppm to 482 ppm.The 4 He/ 20 Ne ratios varied between 1.03 and 65.36, far exceeding that of air (0.318).The measured 3 He/ 4 He ratios (reported as R/Ra, Ra = 1.4 × 10 −6 ) ranged from 0.04 to 0.64 Ra.The air-corrected 3 He/ 4 He ratios (reported as Rc/Ra, Ra = 1.4 × 10 −6 ) ranged from 0.01 to 0.53 Ra (Table 1).The measured δ 13 C CO2 levels in the gas samples ranged from 2.34 to −22.09%.

Characteristics of Geothermal Water
A box plot and Piper diagram were constructed to visually represent the physicochemical parameters and concentrations of major ions (K + , Na + , Ca 2+ , Mg 2+ , Cl − , SO 4 2− , and HCO 3 − ) in the springs, as shown in Figures 2 and 3.In the Piper diagram, the main ions in the JZFZ and NJFZ are Na + , Ca 2+ , and HCO 3 − .The main cations of spring water samples in the WQFZ and LCJFZ were Na + , while the main anions were HCO 3 − and SO 4 2− .According to the Shchukarev classification, the spring water samples collected along the JZFZ and NJFZ predominantly fall into the Na-HCO 3 and Ca-Mg-HCO 3 types.Conversely, within the WQFZ and LCJFZ, the water demonstrates primarily Ca-Na-HCO 3 -SO 4 , Ca-Mg-SO 4 , and Na-Cl-SO 4 types.
To elucidate the origin and evolution of major ions in thermal groundwater, an analysis of the relationships between key ions in the spring water was conducted, including Na + vs. Cl − ; HCO 3 − vs. Na + ; (Ca 2+ + Mg 2+ ) vs. (CO 3 2− + HCO 3 − ); SO 4 2− vs. Ca 2+ ; Cl − vs. SO 4 2− ; and Cl − vs. HCO 3 − (Figure 4).From Figure 4a, it can be seen that only sample (No.L41) aligns with the 1:1 Na/Cl line, indicating the influence of halite dissolution.Sample L41 is located on the LCJFZ, known for its abundant brine production.Along this fault zone, thermal spring water, dissolving Triassic salt-bearing formations, continuously emerges as brackish water rich in salts.Ratios such as γNa/γCl and γNa/γHCO 3 (Figure 4a,b) reveal that water samples with elevated HCO 3 − and Na + concentrations are predominantly distributed in the JZFZ and NJFZ.This implies the significant involvement of carbonate and silicate processes, considering reactions such as NaAlSi 3 O 8 + 8H 2 O → Na + + Al (OH) 4 − + 3H 4 SiO 4 [67].The study area, characterised by widely distributed carbonatite (Triassic marine limestone) and silicate rocks, facilitates abundant HCO 3 -through robust water-rock interactions.Cation exchange between Na-containing minerals and dissolved Ca + further contributes to the high Na + content in the water.Additionally, Figure 4f indicates that γHCO 3 /γCl ratios in spring water along the JZFZ and NJFZ exceed those in the WQFZ and LCJFZ, suggesting that elevated concentrations of HCO 3 -result not only from the dissolution of carbonate minerals but may also come from the introduction of deep-derived CO 2 through conduit-like fault fractures.Some δ 13 C CO2 values suggest the possible addition of mantle-derived CO 2 to the sampled geothermal fluids (see Section 5.3).From Figure 4c, it can be seen that the ratio of γCa + Mg/γCO 3 + HCO 3 in all water from the TRLCF falls below the 1:1 line, indicating the influence of Ca/Mg-rich minerals such as dolomite (Ca The ratios of γCa/γSO 4 can be roughly divided into two groups (Figure 4d), with one group from LCJFZ and WQFZ falling in and below the 1:1 line, suggesting gypsum dissolution as a significant source of Ca 2+ and SO 4 ).The dissolution of gypsum and anhydrite occurs commonly in the clastic sedimentary sandy conglomerates of the Lanping basin of the study area [68].SO 4 2− in spring water samples is related to the gypsum salt layer.The second group, from JZFZ and NJFZ, mostly falls above the 1:1 ratio line, indicating that gypsum dissolution is not the primary source of Ca 2+ and SO 4 2− .Granite, schist, sandstone, and limestone are well developed in the study area, and some springs may percolate into deeper Carboniferous carbonates.The interactions of groundwater with granite and carbonate formations likely contribute to the observed chemical characteristics of the spring water (N63-N69).Finally, the positive correlation between SO 4 2− and Cl − , as indicated in Figure 4e, suggests the dilution of thermal water by shallow cold water.This correlation aligns with the lower SO 4 2− content in spring water along the JZFZ and NJFZ, contrasting with the WQFZ and LCJFZ, as is consistently shown in Figure 3.   − near the WQFZ had slightly higher values than those of the JZFZ, LCJFZ, and NJFZ.This observation suggests that the circulation depth of springs in the WQFZ region is relatively shallow, making them more susceptible to mixing with shallow groundwater.Li + is a highly active element that is easily enriched in thermal spring water and is the signature element of deep liquid upwelling [69].Li + is relatively enriched in samples J14, J21, L41, L42, L43, N64, N71, N72, and N73, where volcanic rocks were found in certain amounts.Rb + and Cs + are typical postmagmatic components and indicate the involvement of magmatic sources.The Rb + and Cs + concentrations of samples J12, J14, J15, J16, J18, L41, N68, and N71 are significantly higher than those of other samples and are closely related to the distribution of volcanic rocks.The study area features the development of carbonate and clastic components, with the majority of springs exhibiting weak alkalinity.Due to the propensity for F − and Sr 2+ to accumulate in weak alkaline water, the concentrations of F − and Sr 2+ in the water bodies of the study region are relatively high [70,71].To elucidate the origin and evolution of major ions in thermal ground ysis of the relationships between key ions in the spring water was condu Na + vs. Cl − ; HCO3 − vs. Na + ; (Ca 2+ + Mg 2+ ) vs. (CO3 2− + HCO3 − ); SO4 2− vs. Ca 2 and Cl − vs. HCO3 − (Figure 4).From Figure 4a, it can be seen that only sam aligns with the 1:1 Na/Cl line, indicating the influence of halite dissolution located on the LCJFZ, known for its abundant brine production.Along thermal spring water, dissolving Triassic salt-bearing formations, continu as brackish water rich in salts.Ratios such as γNa/γCl and γNa/γHCO3 (F veal that water samples with elevated HCO3 − and Na + concentrations are distributed in the JZFZ and NJFZ.This implies the significant involveme

Reservoir Temperature and Circulation Depth
Traditionally, cation ratio geothermometers-such as Na-K, Na-K-Ca, and K-Mg geothermometers-and SiO 2 geothermometers are employed to estimate geothermal reservoir temperatures.The water-rock equilibrium state of thermal water is often visualised using the Na-K-Mg triangular diagram, categorizing spring waters into three types: fully equilibrated waters, partially equilibrated waters, and immature waters [72].Only fully equilibrated waters can be reliably used with cation geothermometers to estimate geothermal reservoir temperatures.In Figure 6, most data points are situated in the area representing immature water, which is evident in the thermal springs plot.Sample No. L41 is positioned in the partially equilibrated area of the triangular diagram, indicating that the wall rock has not achieved equilibrium, likely due to the influx of a significant quantity of cold water (shallow groundwater).This suggests that cation geothermometers may not be suitable for calculating reservoir temperatures in the study area.
Minerals 2024, 14, x FOR PEER REVIEW

Reservoir Temperature and Circulation Depth
Traditionally, cation ratio geothermometers-such as Na-K, Na-K-Ca, an othermometers-and SiO2 geothermometers are employed to estimate geothe voir temperatures.The water-rock equilibrium state of thermal water is often using the Na-K-Mg triangular diagram, categorizing spring waters into three equilibrated waters, partially equilibrated waters, and immature waters [72] equilibrated waters can be reliably used with cation geothermometers to estim mal reservoir temperatures.In Figure 6, most data points are situated in the senting immature water, which is evident in the thermal springs plot.Sample positioned in the partially equilibrated area of the triangular diagram, indicat wall rock has not achieved equilibrium, likely due to the influx of a significant cold water (shallow groundwater).This suggests that cation geothermometers suitable for calculating reservoir temperatures in the study area.When analysing the geothermal spring water samples, mineral saturation indices (SI) were calculated using the PHREEQC software; the results are visually represented in Figure 7.The SI variations for chrysotiles and halites in each spring water sample are notably diverse.The majority of spring water samples predominantly exhibit an unsaturated state (SI < 0) concerning SiO 2 (a), suggesting that amorphous silica does not primarily control the SiO 2 content in the solution.In contrast, chalcedony and quartz approach equilibrium (SI ≈ 0), with chalcedony exhibiting closer proximity to the saturation line than quartz.Fluctuations above and below the nearby equilibrium line for dolomite, aragonite, and calcite indicate variations in the characteristics of the surrounding rock.This phenomenon may signify higher mineral content and extended residence time in the aquifer system [73], as indicated by saturation.Conversely, an unsaturated state suggests that minerals are relatively soluble or have insufficient reaction time with thermal water.These nuanced variations highlight the intricate and complex nature of mineral-water interactions within the studied geothermal springs.
SiO 2 , due to its higher stability in comparison to other minerals, serves as an indicator of thermal reservoir temperature for immature water.Considering this, SiO 2 geothermometers [74] were selected based on Equations ( 1) and ( 2).However, it is important to note that, when the temperature is below 110 • C, chalcedony typically controls the SiO 2 content in the solution rather than quartz [75].Therefore, when the reservoir temperature calculated using Equation ( 1) is below 110 • C, it is recalculated using Equation (2).Additionally, springs with water temperatures below 20 • C are not included in the calculation.
θ chalcedony (  SiO2, due to its higher stability in comparison to other minerals, serves as an indicator of thermal reservoir temperature for immature water.Considering this, SiO2 geothermometers [74] were selected based on Equations ( 1) and ( 2).However, it is important to note that, when the temperature is below 110 °C, chalcedony typically controls the SiO2 content in the solution rather than quartz [75].Therefore, when the reservoir temperature calculated using Equation ( 1) is below 110 °C, it is recalculated using Equation (2).Additionally, springs with water temperatures below 20 °C are not included in the calculation.

𝜃
(℃) = .− 273.15 In Equations ( 1) and ( 2), SiO2 represents the SiO2 concentration in the thermal spring samples, and the corresponding results are presented in Table 1.The temperatures calculated using SiO2 geothermometers fall within the range of 25.6 °C to 157.4 °C.The circulation depth is estimated by utilising the local geothermal gradient and the reservoir temperature of the thermal spring water, as described in Equation (3).
where H is the circulation depth (km), h0 is the depth of the thermostatic zone (km), T is the reservoir temperature (℃), T0 is the temperature of the constant temperature zone (℃), In Equations ( 1) and ( 2), SiO 2 represents the SiO 2 concentration in the thermal spring samples, and the corresponding results are presented in Table 1.The temperatures calculated using SiO 2 geothermometers fall within the range of 25.6 • C to 157.4 • C. The circulation depth is estimated by utilising the local geothermal gradient and the reservoir temperature of the thermal spring water, as described in Equation (3).
where H is the circulation depth (km), h0 is the depth of the thermostatic zone (km), T is the reservoir temperature ( • C), T0 is the temperature of the constant temperature zone ( • C), and q is the geothermal gradient ( • C/km).Considering previous studies on groundwater in some areas of Yunnan Province, a geothermal gradient q of 20 • C/km, an annual average temperature T0 of 15.8 • C, and a constant temperature zone depth h0 of 30 m were assumed [76].Therefore, the final circulation depth of the 73 thermal springs ranged from 0.5 to 7.1 km, while the average circulation depth was 3.8 km (Table 1).

Origin of He
Helium has two stable isotopes with mass numbers 3 and 4. 3 He mainly escapes from the mantle.A large portion of 3 He is primordial He, which was trapped in the solid earth at the time of accretion [77].A small portion of 3 He is produced by the nuclear reaction of 6 Li (n, α) 3 H (−1 β) 3 He.Decaying crustal U and Th predominantly produce 4 He [78].An R/Ra versus 4 He/ 20 Ne diagram was used to identify the origin of the He in the TRLCZ (Figure 8).The He isotopic compositions of the air, MORB, and crust are already well established [79].By assuming R/Ra and 4 He/ 20 Ne ratios of 0.02 and 1000, 1 and 0.318, and 8 and 1000 for the crust, air, and MORB-type mantle, respectively, the mixing proportions were estimated using Sano and Wakita's [79] equations.The boundary lines of this model are calculated using the following mass balance equations:  The detection of mantle-derived helium on the Earth's surface requires both a mechanism to extract helium from the mantle and a pathway for its ascent [88].Geophysical investigations, such as seismic tomography, have provided insights into the origin and migration of mantle volatiles.Studies by Crossey et al. [89], Horiguchi et al. [90], and Zhang et al. [91] revealed that regions with higher emissions of mantle-derived helium correspond to negative anomaly bodies, indicating that mantle volatiles originate from the upwelling asthenosphere.Previous studies have shown that the upper mantle of the southeastern Tibetan Plateau has a series of low shear-wave velocity anomalies beneath the JZFZ, as well as the southern segments of the LCJFZ and NJFZ [40,92].Within the TRLCF, major faults can serve as conduits for the ascent of mantle volatiles.Furthermore, the clockwise movement of the Sichuan-Yunnan rhombic block around the eastern Himalayan syntaxis (EHS) suggests a coupling of deformations between the upper and lower crusts [93].Geophysical observations (e.g., ambient noise tomography and magnetotelluric imaging) and geodynamic modelling suggest that the existence of middle-to-lower crustal flow could contribute to material escape toward the southeast [39][40][41][42].The TRLCZ is located in the frontier zone where materials from the Tibetan Plateau flow southeastward.Recent seismic investigations have identified a low-velocity anomaly indicative of The He isotopic compositions of atmospheric, crustal, and mantle have their own typical characteristics.The 3 He/ 4 He ratio of air (1.4 × 10 −6 ) is defined as 1 Ra.The MORB mantle source has a 3 He/ 4 He ratio of ~8 Ra [80], whereas the production of radiogenic 4 He due to the decay of radioactive U and Th in the crust leads to lower 3 He/ 4 He ratios of 0.005 Ra to 0.02 Ra [81].To conduct a comprehensive comparison, additional data on helium (He) and carbon (C) isotopes from thermal springs near the TRLCZ were needed.To address this, we incorporated Zhou's results (Table S3) into our analysis [82].As illustrated in Figure 8, all data points from the samples fell on the binary mixing lines between the atmospheric and crustal components.However, 12 gas samples were excluded from consideration due to relatively high air contamination during sampling, exceeding 1%.It is worth noting that all the isotopic ratios of helium (Rc/Ra) were found to be less than 1 Ra, indicating that, along the TRLCZ, the helium isotope ratios have experienced significant dilution from crustal radiogenic 4 He.The Rc/Ra ratios of all free gas samples ranged from 0.01 Ra to 0.52 Ra (Table S3), while the Rc/Ra ratios of quite a few samples were higher at >0.1 Ra.The Rc/Ra ratio (>0.1 Ra) can be taken as evidence for the presence of mantle-derived fluids entrained by the hydrologic system [83,84].The contributions of mantle-derived He were calculated using Equations ( 4)- (6).Except for a higher addition of mantle-derived He (6.3%) from the JZFZ, for most of the other samples, the mantle-derived He proportions were lower than 5%.Multiple mechanisms can generate 3 He, including the decay of 6 Li or tritium 3 H, as well as contributions from the mantle [85].However, the decay of 6 Li alone cannot account for the relatively high 3 He/ 4 He ratios observed.This is because the crustal production ratio of 3 He/ 4 He is estimated to be 10 −8 (approximately 0.005-0.02Ra) [86].Furthermore, the contribution of tritiogenic helium to the total 3 He in thermal spring gases is unlikely to be significant due to the low natural tritium background levels in precipitation (usually ≤5 TU) [87].Additionally, there are no active Cenozoic volcanoes in the study area.Therefore, the relatively high 3 He/ 4 He values observed, along with the approximately 2%-5% contribution of mantle-derived helium in the research area, can be attributed to the transport of mantle helium to the surface through faults and fractures.
The detection of mantle-derived helium on the Earth's surface requires both a mechanism to extract helium from the mantle and a pathway for its ascent [88].Geophysical investigations, such as seismic tomography, have provided insights into the origin and migration of mantle volatiles.Studies by Crossey et al. [89], Horiguchi et al. [90], and Zhang et al. [91] revealed that regions with higher emissions of mantle-derived helium correspond to negative anomaly bodies, indicating that mantle volatiles originate from the upwelling asthenosphere.Previous studies have shown that the upper mantle of the southeastern Tibetan Plateau has a series of low shear-wave velocity anomalies beneath the JZFZ, as well as the southern segments of the LCJFZ and NJFZ [40,92].Within the TRLCF, major faults can serve as conduits for the ascent of mantle volatiles.Furthermore, the clockwise movement of the Sichuan-Yunnan rhombic block around the eastern Himalayan syntaxis (EHS) suggests a coupling of deformations between the upper and lower crusts [93].Geophysical observations (e.g., ambient noise tomography and magnetotelluric imaging) and geodynamic modelling suggest that the existence of middle-to-lower crustal flow could contribute to material escape toward the southeast [39][40][41][42].The TRLCZ is located in the frontier zone where materials from the Tibetan Plateau flow southeastward.Recent seismic investigations have identified a low-velocity anomaly indicative of partial melt and flow in the research area at depths ranging from 20 to 30 km [43,94].The presence of partial melting facilitates the release of 3 He, since solid-state processes are insufficiently rapid for the extraction of noble gases [83].

Origin of CO 2
The mantle-derived δ 13 C CO2 concentration ranges from −8‰ to −4‰ [95], and that from limestone thermometamorphism is ±2‰ [96].Mantle-derived CO 2 can exhibit strong variation in its isotopic composition (δ 13 C CO2 = 0‰ to −8‰), mainly due to the variable mixing of mantle-derived CO 2 and crustal contributions [97,98].Organic CO 2 is typically characterised by a strongly negative carbon isotopic signature from −70 to −25‰ [99].A CO 2 / 3 He-δ 13 C mixing model based on the methodology proposed by Sano and Marty [100], and a plot of δ 13 C (‰) vs. CO 2 , were employed to improve the identification of mixing processes.Most samples generally fell on a binary mixing line between the mantle and limestone endmembers (Figure 9a), suggesting the possible addition of mantle-derived CO 2 to the sampled geothermal fluids.As shown in Figure 9b, the increase in CO 2 content was accompanied by 13 C isotope enrichment.Further, carbonate-rich rocks are widespread in this area.Crustal limestone could serve as the primary source for these CO 2 -rich samples, exhibiting high δ 13 C CO2 values, ranging from −1.6‰ to −7.2‰, at the top right of Figure 9b.Considering our previous discussion suggesting that the higher 3 He/ 4 He values in the research area might be the result of minor additions of mantle helium, it is reasonable to consider that trace amounts of mantle CO 2 could also contribute to the geothermal volatiles in the TRLCZ.
in the research area might be the result of minor additions of mantle helium, it is reasonable to consider that trace amounts of mantle CO2 could also contribute to the geothermal volatiles in the TRLCZ.

Spatial Distribution of Hydrogeochemical Characteristics in TRLCF
The kriging spatial interpolation method in Surfer13.0(Surface Mapping System, Golden Software, Inc, Golden, CO, USA) was used to draw contour maps of the gas contents (H2, CO2, He, and mantle He %) of sampling sites in order to expound the spatial discrepancies in thermal waters from the study area (Figure 10).According to Figure 10, the geothermal system in the TRLCF is significantly influenced by the underlying structural characteristics.This influence is evident as the majority of thermal springs are located along four active NS-trending normal and strike-slip crustal faults.The presence of faults and fractures provides pathways for pressured fluids and gases to migrate towards the surface.Figure 11 demonstrates a positive correlation between the depth of circulation and fault activity.In our research area, the JZFZ exhibits a considerably higher average slip rate (~4 mm/a or 5-7 mm/a) compared to the WQFZ (0.89 mm/a or 1.8-2.4mm/a) [30,[35][36][37][38]53].The average circulation depth of thermal spring water along the JZFZ is 4 km, which is relatively deep.Conversely, the northern segment of the WQFZ has the shallowest circulation depth of thermal spring water (ranging from 1 to 3 km, excluding a maximum value of 6 km).The fracture junctions of the LCJFZ with the JZFZ exhibit higher circulation depths of thermal spring water, ranging from 4 to 7 km.The circulation depth of thermal spring water along the southern segments of the NJFZ surpasses that of the northern segments, with depths ranging from 5 to 7 km and 2 to 3 km, respectively.The strike-slip rates of the LCJFZ and NJFZ are approximately ~3 mm/a [31][32][33][34].Additionally, based on earthquake observation since 2009, it is evident that the frequency of Ms ≥ 4 earthquakes is notably higher in active fault zones (JZFZ), where thermal springs have deeper groundwater circulation (as shown in Figure 11).In the structural junction areas, such as the fracture junctions of the LCJFZ with the JZFZ, where the fluid circulation depth is relatively deep, there is a higher propensity for seismic activity, as illustrated in Figures 10d and 11.This observation suggests a close relationship between underground fluids and the seismogenic process in the study area.Deeper circulation depths correspond to more intense water-rock reactions.The seepage of groundwater into the deeper layers increases pore fluid pressure, intensifies the interactions between water and rock, weakens fault zones, and enhances earthquake susceptibility [101][102][103].Therefore, it is crucial to further investigate seismic activity in areas with deep geothermal water circulation.

Spatial Distribution of Hydrogeochemical Characteristics in TRLCF
The kriging spatial interpolation method in Surfer13.0(Surface Mapping System, Golden Software, Inc, Golden, CO, USA) was used to draw contour maps of the gas contents (H 2 , CO 2 , He, and mantle He %) of sampling sites in order to expound the spatial discrepancies in thermal waters from the study area (Figure 10).According to Figure 10, the geothermal system in the TRLCF is significantly influenced by the underlying structural characteristics.This influence is evident as the majority of thermal springs are located along four active NS-trending normal and strike-slip crustal faults.The presence of faults and fractures provides pathways for pressured fluids and gases to migrate towards the surface.Figure 11 demonstrates a positive correlation between the depth of circulation and fault activity.In our research area, the JZFZ exhibits a considerably higher average slip rate (~4 mm/a or 5-7 mm/a) compared to the WQFZ (0.89 mm/a or 1.8-2.4mm/a) [30,[35][36][37][38]53].The average circulation depth of thermal spring water along the JZFZ is 4 km, which is relatively deep.Conversely, the northern segment of the WQFZ has the shallowest circulation depth of thermal spring water (ranging from 1 to 3 km, excluding a maximum value of 6 km).The fracture junctions of the LCJFZ with the JZFZ exhibit higher circulation depths of thermal spring water, ranging from 4 to 7 km.The circulation depth of thermal spring water along the southern segments of the NJFZ surpasses that of the northern segments, with depths ranging from 5 to 7 km and 2 to 3 km, respectively.The strike-slip rates of the LCJFZ and NJFZ are approximately ~3 mm/a [31][32][33][34].Additionally, based on earthquake observation since 2009, it is evident that the frequency of Ms ≥ 4 earthquakes is notably higher in active fault zones (JZFZ), where thermal springs have deeper groundwater circulation (as shown in Figure 11).In the structural junction areas, such as the fracture junctions of the LCJFZ with the JZFZ, where the fluid circulation depth is relatively deep, there is a higher propensity for seismic activity, as illustrated in Figures 10d and 11.This observation suggests a close relationship between underground fluids and the seismogenic process in the study area.Deeper circulation depths correspond to more intense water-rock reactions.The seepage of groundwater into the deeper layers increases pore fluid pressure, intensifies the interactions between water and rock, weakens fault zones, and enhances earthquake susceptibility [101][102][103].Therefore, it is crucial to further investigate seismic activity in areas with deep geothermal water circulation.
Figure 10 shows the spatial distribution of the H 2 concentration, the CO 2 concentration, the He concentration, and mantle He proportions in gases emitted from the four faults of the TRLCF.A comparison of the spatial distribution of hydrogeochemical characteristics (H 2 , CO 2 , He, and mantle He %) and the circulation depth of thermal springs in the different faults shows that there is a positive correlation between the emission of He that has traversed the crust and the circulation depth (higher values of concentrations of He and mantle He% always occur near the sampling points with deeper circulation depths; Figure 10c,d).However, no such correlation existed between the emission of H 2 and CO 2 .A high concentration region of H 2 in the southern segment of the LCJFZ roughly corresponds to a higher circulation depth (Figure 10a).However, no such consistent pattern is observed for the JZFZ, WQFZ, and NJFZ.The elevated concentration of H 2 in the fracture zones is consistent with the idea that H 2 is generated by radical reactions between H 2 O and Si-O or Si-radicals, which appear on the fresh surfaces of fractured rocks [104].The more active the fracture zone is, the more developed the fresh rupture surface of silicate rock in the fracture zone, and the more H 2 is produced [105].Interestingly, the recent 2023 Baoshan M 5.2 earthquake occurred in the range of an anomalous area with high concentrations of H 2 release (Figure 10a).This implies that, in the monitoring of earthquake precursors, the measurement of hydrogen gas content could potentially serve as a favourable geochemical surface observation method.Figure 10 shows the spatial distribution of the H2 concentration, the CO2 concentration, the He concentration, and mantle He proportions in gases emitted from the four faults of the TRLCF.A comparison of the spatial distribution of hydrogeochemical characteristics (H2, CO2, He, and mantle He %) and the circulation depth of thermal springs in the different faults shows that there is a positive correlation between the emission of He that has traversed the crust and the circulation depth (higher values of concentrations of He and mantle He% always occur near the sampling points with deeper circulation depths; Figure 10c,d).However, no such correlation existed between the emission of H2 and CO2.A high concentration region of H2 in the southern segment of the LCJFZ roughly corresponds to a higher circulation depth (Figure 10a).However, no such consistent pattern is observed for the JZFZ, WQFZ, and NJFZ.The elevated concentration of H2 in the fracture zones is consistent with the idea that H2 is generated by radical reactions between H2O and Si-O or Si-radicals, which appear on the fresh surfaces of fractured rocks [104].The more active the fracture zone is, the more developed the fresh rupture surface of silicate rock in the fracture zone, and the more H2 is produced [105].Interestingly, the recent 2023 Baoshan M 5.2 earthquake occurred in the range of an anomalous area with high concentrations of H2 release (Figure 10a).This implies that, in the monitoring of earthquake precursors, the measurement of hydrogen gas content could potentially serve as a favourable geochemical surface observation method.
The concentration of carbon dioxide (CO2) is found to be high in the WQFZ, which is characterised by a shallow circulation depth and low fault activity (Figure 10b).This observation suggests that the elevated CO2 levels in the WQFZ originate from the surface, as there are no available higher-permeability pathways for fluids to migrate from deeper regions.One potential explanation for the heightened concentration of CO2 in the WQFZ is the prevalence of carbonate rocks, distinguishing it from other fault zones.Another possible explanation for the high CO2 concentration in low-activity fault zones is the occurrence of chemical reactions, as documented in previous studies, such as that by Garrels and Christ [106] (7-10).These reactions play a role in diminishing CO2 levels within areas characterised by high fault activity.Notably, the chemical reactions proceed at a faster pace in zones with heightened fault activity as compared to those with lower fault activity.This disparity in reaction rates can be attributed to the larger total surface area of minerals present in high fault activity zones, providing more opportunities for these reactions to take place.The concentration of carbon dioxide (CO 2 ) is found to be high in the WQFZ, which is characterised by a shallow circulation depth and low fault activity (Figure 10b).This observation suggests that the elevated CO 2 levels in the WQFZ originate from the surface, as there are no available higher-permeability pathways for fluids to migrate from deeper regions.One potential explanation for the heightened concentration of CO 2 in the WQFZ is the prevalence of carbonate rocks, distinguishing it from other fault zones.Another possible explanation for the high CO 2 concentration in low-activity fault zones is the occurrence of chemical reactions, as documented in previous studies, such as that by Garrels and Christ [106] (7)(8)(9)(10).These reactions play a role in diminishing CO 2 levels within areas characterised by high fault activity.Notably, the chemical reactions proceed at a faster pace in zones with heightened fault activity as compared to those with lower fault activity.This disparity in reaction rates can be attributed to the larger total surface area of minerals present in high fault activity zones, providing more opportunities for these reactions to take place.
Although helium (He) was a minor component of the gases, it increased in quantity at fracture junctions of the JZFZ (Figure 10c).The contributions of mantle helium (He) from the JZFZ, the southern segment of the NJFZ, and the LCJFZ were estimated to be in the range of 2-5% (Figure 10d).Several geophysical investigations (seismic tomography) revealed that higher emissions of mantle-derived helium are found above low-velocity bodies, implying that mantle volatiles originate from the upwelling asthenosphere [89][90][91].Previous studies have noted the presence of a series of low shear-wave velocity anomalies in the upper mantle beneath the JZFZ, as well as the southern segments of the LCJFZ and NJFZ (situated between 24 • and 26 • N) in the southeastern Tibetan Plateau (see Figure 12) [40,92].Within the JZFZ, there are large-scale low-resistivity zones in the middle and upper crust, which may be related to partial melting and heat flow [107,108].The spatial distribution of mantle He contribution (>2%) along the southern segments of the LCJFZ and NJFZ (positioned between 24 • and 26 • N) is approximately 70-80 km away from the Tengchong volcanic field (TC; refer to Figure 12).The Tengchong volcanic field (TC) is representative of a bottom boundary zone with upwelling in the asthenosphere.Its 3 He/ 4 He ratio consistently surpasses 1 Ra, reaching a maximum of 6 Ra, and the percentage of mantle-derived helium release exceeds 60% [109].This observation implies a discernible influence of the Tengchong volcanic field on the release of mantle-derived helium along the southern segments of the LCJFZ and NJFZ.Conversely, the blue domain in Figure 12 represents HVZs (high-velocity zones) beneath the WQFZ, implying that the tectonic activity in this region may be blocked by stiffer rocks.Additionally, Figure 12 provides visual evidence that regions characterised by higher 3 He/ 4 He ratios (Rc/Ra) in surface waters align with the areas where the lowvelocity anomaly zone (LVZ) is observed.This correlation strongly suggests that these regions are closely linked to the upwelling of hot material within the upper mantle.This conclusion is further supported by a correlation (R 2 = 0.9) between Ra and the mantle velocity in the TRLCF (Figure 12, inset).Areas characterised by high conductivity and low velocity within a fault zone are typically interpreted as zones of fault disruption resulting from impact sliding or shearing along the fault.These zones are often associated with the presence of deep-seated magma and/or metamorphic fluids.  He values, reported relative to air, from thermal springs and cold sprin study region are overlaid on S-wave mantle tomography in the upper mantle (modifie et al. [40]; Fu et al. [92]).Inset shows the correlation (R 2 = 0.9) of 3 He/ 4 He in each veloci to the mantle S-wave velocity (Wang et al. [40]).
The transport of He volatiles through faults and fractures whose permea trolled by fault activity, and the correlation of 3 He/ 4 He ratios with time-averag rates, suggest that higher slip rates are associated with more permeable ve search conducted by Wang et al. [110] has indicated that there is a positive between the emission of He that has traversed the crust and the fault slip rate In our research area, the 3 He/ 4 He ratios in thermal spring gases along the JZ LCJFZ, and NJFZ were 0.10-0.52Ra, 0.01-0.35Ra, 0.04-0.30Ra, and 0.10-0.33tively.The 3 He/ 4 He values (Rc/Ra) along the JZFZ exhibit the highest ratios four faults studied, while the WQFZ has the lowest 3 He/ 4 He values (Rc/Ra) c the other three faults.By utilising the linear trend observed in the graph, it i  4 He values, reported relative to air, from thermal springs and cold springs within the study region are overlaid on S-wave mantle tomography in the upper mantle (modified from Wang et al. [40]; Fu et al. [92]).Inset shows the correlation (R 2 = 0.9) of 3 He/ 4 He in each velocity bin relative to the mantle S-wave velocity (Wang et al. [40]).
The transport of He volatiles through faults and fractures whose permeability is controlled by fault activity, and the correlation of 3 He/ 4 He ratios with time-averaged fault slip rates, suggest that higher slip rates are associated with more permeable vents [91].Research conducted by Wang et al. [110] has indicated that there is a positive correlation between the emission of He that has traversed the crust and the fault slip rate (Figure 13).In our research area, the 3 He/ 4 He ratios in thermal spring gases along the JZFZ, WQFZ, LCJFZ, and NJFZ were 0.10-0.52Ra, 0.01-0.35Ra, 0.04-0.30Ra, and 0.10-0.33Ra, respectively.The 3 He/ 4 He values (Rc/Ra) along the JZFZ exhibit the highest ratios among the four faults studied, while the WQFZ has the lowest 3 He/ 4 He values (Rc/Ra) compared to the other three faults.By utilising the linear trend observed in the graph, it is possible to estimate the strike-slip rates based on the 3 He/ 4 He ratios (Rc/Ra).The calculated average slip rate for the JZFZ is notably higher at approximately 4 mm/a, whereas the WQFZ has a lower average slip rate of around 1 mm/a.The strike-slip rates of the LCJFZ and NJFZ fall within the range of 2-3 mm/a.These calculated strike-slip rates are consistent with some prior geological findings (see Section 2) [30][31][32][33][34][35][36][37][38], further supporting the notion that the release of mantle-derived helium in gas bubbling from thermal springs along active faults serves as an indicator of fault activity.The elevated helium isotope values of thermal spring gas in this region suggest that there has been a significant upwelling of deep fluids [84], which could have the effect of weakening the fracture [111,112].Moreover, it increases the pore pressure within the fracture, which is a key factor that encourages the incidence and breeding of earthquakes [113].By studying the distribution characteristics of seismic activity (M ≥ 5.0) in the TRLCZ (Figure 10d), it can be seen that the regions with seismic activity are often accompanied by mantle-derived He release, indicating that mantle fluids play a critical role in the formation and occurrence of earthquakes in this region.Therefore, thermal springs with mantle-derived characteristics can be used as observation objects for earthquake monitoring and prediction, and mantle He contributions, gas component content, and circulation depth can be employed as parameters of geochemical observations.

Fluid Circulation Model
Figure 14 illustrates the fluid circulation model governing groundwater migration and deep mass transfer in the TRLCZ.Meteoric water infiltrates fractured metamorphic and clastic rocks within high mountainous regions to be heated as thermal water.The prevalence of a thick red mudstone layer, serving as the cover for the geothermal system, contributes to poor permeability, fostering heat accumulation during groundwater movement [114].The existence of prevalent NW-and NE-trending fault systems, coupled with intense seismic activities, fractures the outcropped rocks in these regions.Consequently, meteoric water easily percolates in the depths through structurally weak zones.As discussed in Section 5.1, the chemo-diversity of the studied geothermal water is constrained by the characteristics of host rocks and the extent of water-rock interactions.Their evolution appears to be controlled by specific tectonic units.For instance, spring water samples collected along the JZFZ have Na-HCO3 and Ca-Mg-HCO3 compositions, respectively.In contrast, within the WQFZ (located in the Lanping Basin), the water changes to the Ca-Na-HCO3-SO4 and Na-SO4 types.This characteristic shift is indicative of increased interaction with carbonate-contained aquifers during the water's flow from northern mountain regions to southern basins.This characteristic change is also associated with the diverse lithology of aquifers, with the JZFZ predominantly composed of resistant high-grade metamorphic rocks and the WQFZ containing limestone, red-coloured terrestrial siliciclastic

Fluid Circulation Model
Figure 14 illustrates the fluid circulation model governing groundwater migration and deep mass transfer in the TRLCZ.Meteoric water infiltrates fractured metamorphic and clastic rocks within high mountainous regions to be heated as thermal water.The prevalence of a thick red mudstone layer, serving as the cover for the geothermal system, contributes to poor permeability, fostering heat accumulation during groundwater movement [114].The existence of prevalent NW-and NE-trending fault systems, coupled with intense seismic activities, fractures the outcropped rocks in these regions.Consequently, meteoric water easily percolates in the depths through structurally weak zones.As discussed in Section 5.1, the chemo-diversity of the studied geothermal water is constrained by the characteristics of host rocks and the extent of water-rock interactions.Their evolution appears to be controlled by specific tectonic units.For instance, spring water samples collected along the JZFZ have Na-HCO 3 and Ca-Mg-HCO 3 compositions, respectively.In contrast, within the WQFZ (located in the Lanping Basin), the water changes to the Ca-Na-HCO 3 -SO 4 and Na-SO 4 types.This characteristic shift is indicative of increased interaction with carbonatecontained aquifers during the water's flow from northern mountain regions to southern basins.This characteristic change is also associated with the diverse lithology of aquifers, with the JZFZ predominantly composed of resistant high-grade metamorphic rocks and the WQFZ containing limestone, red-coloured terrestrial siliciclastic rocks, gypsum, and halite.
water.This process leads to the outcrops of warm springs.On another note, fluids derived from the mantle, enriched with 3 He/ 4 He isotopes and CO2, ascend through the ductil crust due to the upwelling of thermal material in the upper mantle.The distribution o thermal springs in the TRLCF is significantly influenced by active fractures, acting as con duits for the ascent of mantle volatiles.The activity of these fractures also impacts ga emission activity [115].Fluid migration towards the focal region of seismic incubation fa cilitates the lubrication of faults, activating existing structures in the shallower crust and contributing to the occurrence of strong earthquakes [116].Consequently, in charac terised by intense seismic activity, deep faults exhibiting high permeability, and the up ward movement of mantle fluids from the depths, elevated releases of mantle-derived he lium frequently manifest at the Earth's surface.Quartz geothermometers indicate some higher reservoir temperatures (>130 • C) located in the JZFZ and the southern segments of the NJFZ.The JZFZ and NJFZ are active large-scale fault zone faults, and low-velocity zones within these zones [43,94], linked to partial melt in the upper-middle crust, can serve as deep water circulation and heat sources.Consequently, elevated temperatures enhance water-rock interaction, as observed in some springs in these areas.Thermal springs in other regions exhibit lower reservoir temperatures, mainly below 90 • C, indicating shallower depths, particularly notable in the WQFZ situated in the Lanping Basin.
The flow direction of thermal water in the TRLCZ is from the mountain region towards the interior of the basin.The thermal water obtains salinity from the diverse lithology of aquifers, subsequently rising and eventually mixing with shallow underground water.This process leads to the outcrops of warm springs.On another note, fluids derived from the mantle, enriched with 3 He/ 4 He isotopes and CO 2 , ascend through the ductile crust due to the upwelling of thermal material in the upper mantle.The distribution of thermal springs in the TRLCF is significantly influenced by active fractures, acting as conduits for the ascent of mantle volatiles.The activity of these fractures also impacts gas emission activity [115].Fluid migration towards the focal region of seismic incubation facilitates the lubrication of faults, activating existing structures in the shallower crust and contributing to the occurrence of strong earthquakes [116].Consequently, in areas characterised by intense seismic activity, deep faults exhibiting high permeability, and the upward movement of mantle fluids from the depths, releases of mantle-derived helium frequently manifest at the Earth's surface.

Conclusions
In this paper, we investigated the hydrochemical characteristics and origins of He and C isotopes in 73 thermal springs located in the TRLCZ.The analysis of the relationships between major ions in thermal groundwater provides a comprehensive understanding of their origin and evolution.Dominant processes include carbonate and silicate interactions, with elevated HCO 3 − and Na + concentrations observed in the JZFZ and NJFZ regions.The influence of Ca/Mg-rich minerals, such as dolomite, is apparent, impacting the cation composition of water.Additionally, the ratios of Ca/SO 4 indicate the significant role of gypsum dissolution, particularly in the LCJFZ and WQFZ.The positive correlation between SO 4 2− and Cl − implies dilution by shallow cold water.This study also investigated the depth of circulation in the fault zones, revealing variations among the different faults.The JZFZ faults exhibit a relatively deeper circulation depth, ranging from 4 to 7 km.On the other hand, the northern segment of the WQFZ has the shallowest circulation depth of thermal spring water (1-3 km).At fracture connections with the JZFZ, the LCJFZ revealed greater depths (4-7 km), and the circulation depth of thermal spring water along the southern segments of the NJFZ surpasses that of the northern segments, with depths ranging from 5 to 7 km and 2 to 3 km, respectively.Notably, a positive correlation exists between the depth of circulation and fault activity, suggesting that areas with deeper groundwater circulation are more prone to earthquakes.This finding emphasizes the importance of considering seismic activity in regions with deep geothermal water circulation for earthquake monitoring and prediction.
We analysed the Rc/Ra ratios of free gas samples, which ranged from 0.01 Ra to 0.52 Ra.Certain samples exhibited higher Rc/Ra ratios along with a greater contribution of mantlederived helium, potentially indicating minor additions of mantle helium through faults and fractures to the surface.Crustal limestone was identified as the primary source of CO 2 -rich samples, exhibiting high δ 13 C CO2 values ranging from −1.6‰ to −7.2‰, while the presence of a trace amount of mantle CO 2 cannot be ruled out as a contributing factor.The spatial distribution of the H 2 concentration, CO 2 concentration, He concentration, and mantle He proportions in gases showed that higher values of He and mantle He% concentrations always occur near the sampling points with deeper circulation depths.However, no similar correlation was observed in the emission of H 2 and CO 2 .Additionally, most earthquakes of magnitude 5 or greater occurred near the region of high values of mantle source He release, highlighting the critical role of mantle fluids in the formation and occurrence of earthquakes in the region.Hence, the inclusion of mantle-derived He contributions and the variations in relevant gas components should be considered as parameters in geochemical observations.Furthermore, a fluid circulation model was developed to describe the process of fluid (water and gas) circulation migration and earthquake generation in the TRLCZ.
Overall, this study provides valuable insights into the hydrochemical characteristics, origins of He and C isotopes, and their relationships with tectonic structures and seismic activity in the TRLCZ.The findings contribute to a better understanding of the geology and geodynamics of the region, which can aid in earthquake monitoring and prediction efforts.

Figure 2 .
Figure 2. Physicochemical parameters and the concentrations of major ions of thermal springs in the TRLCZ.

Figure 2 .
Figure 2. Physicochemical parameters and the concentrations of major ions of thermal springs in the TRLCZ.To analyse the differences in the contents of Li + , NH 4 + , Rb + , Cs + , Sr 2+ , F − , NO 2 − , Br − , and NO 3 − ions in thermal springs along JZFZ, WQFZ, LCJFZ, and NJFZ, a box plot was created, as shown in Figure 5. Ammonium is an important parameter of geothermal fluids and high NH 4 + value indicates that there may be input from anthropogenic sources.The presence of NH 4 + , NO 3 − , and NO 2 − can indicate the degree of groundwater contamination by surface water.Typically, deep groundwater has little to no NO 3 − present, whereas the presence of NH 4 + and NO 2 − indicates that the groundwater may have recently been contaminated.The content of NO 3− , NH 4 + , and NO 2 − near the WQFZ had slightly higher values than those of the JZFZ, LCJFZ, and NJFZ.This observation suggests that the circulation depth of springs in the WQFZ region is relatively shallow, making them more susceptible to mixing with shallow groundwater.Li + is a highly active element that is easily enriched in thermal spring water and is the signature element of deep liquid

Figure 2 .
Figure 2. Physicochemical parameters and the concentrations of major ions of th the TRLCZ.

Figure 5 .
Figure 5. Box plot of geothermal waters from the study area.

Figure 5 .
Figure 5. Box plot of geothermal waters from the study area.

Figure 6 .Figure 6 .
Figure 6.The Na-K-Mg trilinear equilibrium diagram of hydrothermal water of the TR

Figure 7 .
Figure 7. Saturation indices values of water samples with respect to minerals.

Figure 7 .
Figure 7. Saturation indices values of water samples with respect to minerals.

Figure 10 .
Figure 10.Spatial distribution of hydrogeochemical characteristics (H 2 , CO 2 , He, and mantle He %) and the circulation depth of thermal springs in the different faults of the TRLCF.(a) H 2 concentration; (b) CO 2 concentration; (c) He concentration; (d) mantle He proportions.

Figure 11 .
Figure 11.Plot of the depth of circulation versus the rate of slip (the size of bubbles reflects the frequency of earthquakes with a magnitude of four or higher occurring along the fault).

Figure 11 .
Figure 11.Plot of the depth of circulation versus the rate of slip (the size of bubbles reflects the frequency of earthquakes with a magnitude of four or higher occurring along the fault).

Figure 12 .
Figure 12. 3 He/4 He values, reported relative to air, from thermal springs and cold sprin study region are overlaid on S-wave mantle tomography in the upper mantle (modifie et al.[40]; Fu et al.[92]).Inset shows the correlation (R 2 = 0.9) of 3 He/4 He in each veloci to the mantle S-wave velocity (Wang et al.[40]).

Figure 12 .
Figure 12. 3 He/4 He values, reported relative to air, from thermal springs and cold springs within the study region are overlaid on S-wave mantle tomography in the upper mantle (modified from Wang et al.[40]; Fu et al.[92]).Inset shows the correlation (R 2 = 0.9) of 3 He/4 He in each velocity bin relative to the mantle S-wave velocity (Wang et al.[40]).

Figure 14 .
Figure 14.Fluid circulation model in the TRLCZ.Figure 14. Fluid circulation model in the TRLCZ.

Figure 14 .
Figure 14.Fluid circulation model in the TRLCZ.Figure 14. Fluid circulation model in the TRLCZ.

Table 1 .
The estimated reservoir temperatures of thermal springs in the TRLCZ.