Heat generation effects from shear friction along Xianshui river strike-slip fault in western Sichuan, China
Introduction
The high-temperature geothermal belt in western Sichuan is one of the most active geothermal regions in China. It is located in the Tethys tectonic region, which lies in the eastern part of the Himalayan Geothermal Belt (HGB) (Fig. 1a) and is an area with intense tectonic activity. Investigation has indicated that the average value of the heat flow is higher than 80 mW/m2. Recent results of hydrogeochemistry research and drilling work (Cao et al., 2006; Guo et al., 2017; Jiao et al., 2018; Luo et al., 2017) have proven the existence of numerous high-temperature geothermal systems in the west Sichuan high-temperature geothermal belt, which can be subdivided into the Batang, Litang, and Kangding geothermal belts, according to the distribution of hot springs from the west to the east (Tang et al., 2017). All the types of high-temperature surface geothermal manifestations that occur within China can be found in this region, including boiling springs, boiling fountains, steam vents, steaming ground, and a hydrothermal exploration zone (Liao, 1990a; Liao, 1990b), which are distributed with banding structures.
High-temperature geothermal systems are closely related to volcanic and seismic activities and can act as superficial supports and constraints for corresponding geodynamic processes. Therefore, defining the origin of high-temperature geothermal systems is not only important for conducting geothermal and hydrogeological groundwork, but it is also beneficial for resolving geodynamic scientific issues. In addition, modeling and analyzing the generation process of geothermal systems is important for the evaluation and development of geothermal resources. Furthermore, the modern generation of hydrothermal mineralization has an evident relation to geothermal systems. Thus, the results of associated studies can be used as a reference in order to predict the presence of metal and mineral resources.
Although the distribution of the springs in the studied area has been thoroughly explored and detailed, the generation of the geothermal system process has not been completely apprehended. Thus, exploring the heat source is the key to analyzing the problems associated with the generation of geothermal systems. The high crustal heat beneath the HGB has been explained and discussed (Bird, 1978; England et al., 1992; Hochstein and Regenauer-Lieb, 1998). The exposed volcanic rocks, generated more than 30 Ma ago, having anomalously low 3He/4He ratios of geothermal gases that were collected from several high-temperature systems within the HGB show that the high heat flow in the HGB (including the western Sichuan geothermal systems) has an intra-crustal origin, and is not related to the intrusions of sub-crustal melts that have been occurring during the collisional period (Giggenbach et al., 1983). Research (Regenauer-Lieb, 1992) has shown that plastic deformation of the ductile crust may also have played a role in forming the HGB and the enriched radioactive elements in crust due to the crustal thickening during collision may have increased the radioactive heat flow. The uplift of the thickening crust related to the isostatic compensation has led to the elevation of the isothermal surface and increased the heat flow background value. It has also been proposed that shear heating occurs in the crustal thrust interface between the Indian and Asian plates beneath the HGB(Bird, 1978; England et al., 1992; Shi and Zhu, 1993). Furthermore, heat can be stored within magma — in the form of latent heat when partial melting is induced by fault shear frictional heat — and can affect the heat flow to some extent. Heat from the late-magmatic intrusion (<10 Ma) and released melting latent heat generated from magma solidification can also be considered. Overall, the heat source is mostly located within the crust in the western Sichuan high-temperature geothermal systems.
The focus of this study is the Xianshui river high-temperature geothermal systems, controlled by the Xianshui river strike-slip fault, where the Yulingong in Kangding, Yuemachaqu, and Qimeikechaqu in Daofu hot springs are located along the fault. The well was drilled within Kangding City during 2013, which has shown that temperatures at the top and bottom of the well reach 140 ℃ and 208 ℃ at 267 m, respectively; this classifies the springs as high-temperature systems and implies that the Kangding geothermal systems are possibly significant latent heat resources. Studies on the hydrogeochemical and isotopic characteristics of hot springs and gases in the Kangding geothermal area have also been conducted, and it has been found that the reservoir temperature reaches approximately 260 °C (Guo et al., 2017; Luo et al., 2017). Furthermore, the heat-storage value of the Kangding geothermal systems has also been estimated (Liao, 2018; Wang et al., 2019).
Nevertheless, it is currently unclear how a geothermal system such as Kangding is controlled by the Xianshui river strike-slip fault; therefore, this study conducts quantitative simulations of temperature evolution beneath the Xianshui river strike-slip fault using COMSOL Multiphysics software to elucidate this. The contributions of mantle heat flow, radioactive heat flow, shear frictional heat flow, block uplift heat flow, and melting latent heat flow to surface heat flow, are also respectively calculated. The simulations consider several important developments: first, fault geometry inferred from analysis of historical seismic events since 1870 is used to constrain the model structure; second, different kinematic histories of the Xianshui river strike-slip fault are tested; third, high-temperature and pressure experimental petrology physics data are used to assist in understanding the important feedback processes between temperature, shear heat generation, and the thermal structure. As a typical representation of the western Sichuan geothermal belt, the modeling results for the Kangding geothermal system may provide an insight into the deep dynamic mechanism and thermal controlling structures of the system and other analogous high-temperature geothermal systems in the Tethys tectonic region.
Section snippets
large-scale setting
The western Sichuan plateau is located at the junction of the Qiangtang, Songpan-ganzi, and Yangtze active blocks (Fig. 1a) (Copley et al., 2010; Harrowfield and Wilson, 2005; Zhang, 2008). It is part of the eastern Tibet Plateau that formed during the collision and continuous compression of the Indian and Eurasian plates (beginning approximately 55 Ma ago) (Clark and Royden, 2000; Tapponnier et al., 2001; Xu et al., 2006), which resulted in the formation of numerous deep faults extending down
Model geometric constraints
In this study, we simplified the investigation by focusing on a well-constrained, two-dimensional, vertical, cross-sectional model of the Kangding geothermal system which has a length and depth of 100 km and 30 km, respectively. The position of the cross-section is represented by the black solid line in Fig. 1. For an initial approximation, we simulated the Xianshui river fault friction heating process using a kinematic model, which is broadly consistent with geological and geophysical
Results
According to the key parameters discussed in Section 3, such as geometric constraints, kinematic parameters, and rock rheology, we conducted a series of frictional heat simulations for the Xianshui river fault zone.
Fig. 3a shows the simulated results of surface heat flow values along the section, where the horizontal coordinate of 0 m corresponds to the position of the Xianshui river fault zone, and the different colors represent the surface heat flow values calculated under different
Discussion
Unlike many previous studies, this model achieves a coupling between thermal conductivity, thermal capacity, and the temperature of rock in the crust, which share complex feedback processes between them. In this respect, the fault slip generates shear frictional heat on the contact surface of the fault, and the temperature of the surrounding rock increases. Thermal conductivity then sharply decreases while the heat capacity increases, and the thermal diffusion coefficient further decreases with
Conclusions
Based on geometric constraints, kinematic parameters, rheology, and other parameters, a 2D model of the Xianshui river fault zone is established and a series of thermal simulations are conducted based on the use of different parameter values. Combined with geological, hydrological, geothermal, and geophysical data, the heat generation effects of shear friction along the Xianshui river strike-slip fault on the high temperature geothermal activity in the studied region is also discussed. In
Funding
This work was supported by the Strategic Priority Research Program of the Chinese Academy of Sciences, Grant No. XDB42020104 and the National Natural Science Foundation of China [No. 41430319, 41574074].
CRediT authorship contribution statement
Yifei Ai: Conceptualization, Formal analysis, Methodology, Software, Writing - original draft. Jian Zhang: Conceptualization, Writing - review & editing. Miao Dong: Validation. Beiyu Wang: Investigation. Gui Fang: Investigation.
Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
References (83)
The fracture strength and frictional strength of Weber Sandstone
International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts
(1975)On the magnitude and the fault length of the great Wenchuan earthquake
Science and Technology Review. (in Chinese)
(2008)- et al.
Isotopic and chemical composition of parbati valley geothermal discharges, North-West Himalaya, India
Geothermics
(1983) - et al.
Fluid geochemistry and geothermometry applications of the Kangding high-temperature geothermal system in Eastern Himalayas
Applied Geochemistry
(2017) - et al.
Indosinian deformation of the Songpan Garze Fold Belt, northeast Tibetan Plateau
Journal of Structural Geology
(2005) - et al.
Heat generation associated with collision of two plates: the Himalayan geothermal belt
Journal of Volcanology and Geothermal Research
(1998) Permian tectonic framework and palaeogeography of SE Asia
Journal of Asian Earth Sciences
(2002)- et al.
Miocene enplacement and deformation of the Konga Shan granite(Xianshuihe fault zone, west Sichuan, China): geodynamic implications
Earth and Planetary Science Letters
(1995) - et al.
Some thermatectonic aspects of the Tibetan plateau
Tectonophysics
(1993) - et al.
A revised method for heat flux measurement with applications to the fracture-controlled Kangding geothermal system in the Eastern Himalayan Syntaxis
Geothermics
(2019)
Cooling history of the Gongga batholith: Implications for the Xianshuihe Fault and Miocene kinematics of SE Tibet
Earth and Planetary Science Letters
On the computations of the effect of the attraction of the mountain masses as disturbing the apparent astronomical latitude of stations in geodetic surveys
Philos.Trans. R. Soc. London Ser. A
Field study of a highly active fault zone: The Xianshuihe fault of southwestern China
Geological Society of America Bulletin
Crustal deformation of the eastern Tibetan plateau revealed by magnetotelluric imaging
Nature Geoscience.
Crustal deformation of the eastern Tibetan plateau revealed by magnetotelluric imaging
Nature Geoscience
Himalayan tectonics explained by extrusion of a low-viscosity crustal channel coupled to focused surface denudation
Nature
Initiation of intracontinental subduction in Himalaya
Journal of Geophysical Research
Frictional slip of granite at hydrothermal conditions
Journal of Geophysical Research Solid Earth
Comparision of geochemical features pf warm springs between Chongqing and Kangding.
Surface rupture zones of the 2001 earthquake Ms8.1 west of Kunlun pass, northern Qinghai-Xizang plateau
Quaternary Sciences. (in Chinese)
The Ms7.1 Yushu earthquake surface ruptures and historical earthquakes
Chinese Sci Bull (in Chinese)
Late Cenozoic episodic uplifting in southeastern part of the Tibetan plateau-evidence from Ar-Ar thermochronology
Acta Petrologica Sinica (in Chinese)
Topographic ooze: Building the eastern margin of Tibet by lower crustal flow: Geology, 28
no.
India-Asia collision and the Cenozoic slowdown of the Indian plate: Implications for the forces driving plate motions
Journal of Geophysical Research-Solid Earth
Seismic activities and earthquake potential in the Tibetain Plateau
Chinese J. Geophys. (in Chinese)
Heat-sources for Tertiary metamorphism and anatexis in the Annapurna-Manaslu region central Nepal
Journal of Geophysical Research-Solid Earth
The study of gravity-magnetic anomaly and tectonic background in Sichuan west region
Chinese J. Geophys. (in Chinese)
Deep conductivity characteristics of the southern Xianshuihe fault zone Chinese
J. Geophys. (in Chinese)
The water-undersaturated and dry Qz-Ab-Or system revisited
Experimental results at very low water activities and geological implications: Contributions to Mineralogy and Petrology
Pn-wave velocity structure and anisotropy in the upper mantle of sichuan, yunnan and adjacent regions
Geochemistry of geothermal fluids with implications on the sources of water and heat recharge to the Rekeng high-temperature geothermal system in the Eastern Himalayan Syntax
Geothermics
Distribution of active rock uplift along the eastern margin of the Tibetan Plateau: inferences from bedrock channel longitudinal profiles
Journal of Geophysical Research
Late Cenozoic evolution of the eastern margin of the Tibetan Plateau: Inferences from Ar-40/Ar-39 and (U-Th)/He thermochronology
Tectonics
Co-seisimic surface rupture and dextral-slip oblique Thrusting of the Ms8.0 Wenchuan earthquake
Acta Geologica Sinica. (in Chinese)
Yunnan-Tibet geothermal belt
Paper read at The 6th Annual Academic Meeting of the Chinese Geophysical Society, at Wuhan, Hubei, China.
Yunnan-Tibet geothermal belt
Paper read at The 6th Annual Academic Meeting of the Chinese Geophysical Society, at Wuhan, Hubei, China.
Thermal Springs and Geothermal Energy in the Qinghai-Tibetan Plateau and the Surroundings.
Geothermal potential evaluation and development prioritization based on geochemistry of geothermal waters from Kangding area, western Sichuan
China: Environmental Earth Sciences
Inquisition of the distribution and cause of the hot springs in western Sichuan
Journal of Chongqing Teachers College (Natural Science Edition) (in Chinese)
Cited by (20)
Lithospheric and asthenospheric properties of the saharan platform inferred from potential field, geoid and heat flow data
2024, Journal of African Earth SciencesHydrogeochemical characterization and precursor anomalies of hot springs in the North Tianshan orogen
2023, Applied GeochemistryTerrestrial heat flow measurement and numerical simulation of lithospheric thermal structure in western Sichuan, China
2023, Geoenergy Science and EngineeringModeling and interpretation of geomagnetic data related to geothermal sources, Northwest of Delijan
2022, Renewable EnergyCitation Excerpt :Geothermal systems are closely related to seismic events and volcanic activity, which reflect zones of complex geodynamic processes [1].
Heat flux measurements and thermal potential of the Garze geothermal area in the eastern Himalayan Syntaxis
2021, Journal of Volcanology and Geothermal ResearchCitation Excerpt :Moreover, the geothermal area should have received additional heat source as its heat flux to exceed that of some geothermal areas on the edge with magmatic activity. Such a replenishment heat source is probably tectonic friction heat (Wang et al., 2019; Ai et al., 2021). It is clear that the heat flux obtained by using eddy covariance tower has a large variation of about 66% on average, such as the Solfatara geothermal area (Fig. 11b).