ORIGINAL RESEARCH article

Front. Earth Sci., 09 May 2023

Sec. Geochemistry

Volume 11 - 2023 | https://doi.org/10.3389/feart.2023.1117862

Soil gas CO2 emissions from active faults: a case study from the Anninghe—Zemuhe fault, Southeastern Tibetan Plateau, China

  • 1. CEA Key Laboratory of Earthquake Prediction, Institute of Earthquake Forecasting, China Earthquake Administration, Beijing, China

  • 2. School of Water Resources and Environment, China University of Geosciences (Beijing), Beijing, China

  • 3. University of Sydney Business School, Sydney, NSW, Australia

Abstract

Introduction: Carbon dioxide emissions from non-volcanic areas are undervalued in the carbon cycle.

Methods: First estimates of diffuse CO2 flux from the Anninghe—Zemuhe fault (AZF), Southeastern Tibetan Plateau, China, which suggests this could equal 15% emissions from all volcanoes in China. Following the accumulation chamber method, CO2 flux was investigated at 1,483 points, and along 67 profiles crossing the AZF.

Results and discussion: Total CO2 emissions from the AZF were estimated 1.2 Mt yr-1. The relationship between soil gas CO2 fluxes, earthquakes, and fault activity was discussed. The intense fault activity in the southern part of the Zemuhe fault (ZMHF) and the northern part of the Anninghe fault (ANH) was inferred, which could have enhanced the porosity of the soil, and accelerated the water-rock interactions and soil gas emission within the fault zone. The chemical and isotopic data indicated that biogenic CO2 was the primary source of CO2 from the AZF. Produced by interactions between groundwaters and carbonates, soil gas CO2 could migrate to the near surface through cracks. Spatial variations of CO2 flux in soil gas indicate that seismic activity could be responsible for the jumpy variations of CO2 flux. The diffuse CO2 from deep faults may contribute considerably to the greenhouse gas cycles.

1 Introduction

Consistently, Global warming is a central focus of climate change and studies have suggested that the growing number of carbon emissions has been considered as an essential factor for global warming (; ; ). Nowadays, studying the distribution and magnitude of global carbon sources and sinks is an international mainstream for a better understanding of rates and mechanism of carbon cycling in and out of earth (; ; ; Xu et al., 2022). In addition to human activities, natural factors (e.g., volcanic eruptions, biological respiration, seismic activity) will also contribute to the increase of greenhouse gas concentration in the atmosphere (; ). Therefore, a great interest has been focused on studying of CO2 Earth degassing. There are many targets for such research, such as the relationships between the CO2 flux and the tectonic structures (; ), the researches of volcanic degassing (; ; ), and the quantification of deeply derived CO2 discharging into the atmosphere (; ; ; ; ).

However, the contribution of CO2 released from non-volcanic areas, particularly deep and large faults, is also worthy of attention (; ; Zhang et al., 2021). Based on a study of 3He and CO2 in springs and wells, estimated a mantle-derived CO2 flux of ∼0.02 gm-2 d−1 for the entire San Andreas Fault system. The relationship between active crustal stress and soil CO2 flux in southern Italy suggested that crustal stress associated with the seism genic process can effectively modulate the gas release in a seismically active area (). On the Chaozhou fault in China, significant anomalies in soil gas He and Rn were observed before the earthquake, which may reflect changes in the stress field before the earthquake (); The sudden increase in CO2 and Rn before the earthquake was observed on a fault in the earthquake cluster area in northwestern Bohemia, Czech Republic, which may indicate an increase in porosity and the opening of transport channels due to stress redistribution (). Weinlich and others concluded that underground gas components can objectively and sensitively reflect the stress and tectonic activity changes of the earth’s crust, and usually show various anomalies before or after the occurrence of earthquakes.

AZF is a seismic zone with vigorous seismic activity, which has high slip rates and developed fractures providing a channel for fluid migration deep and CO2 degassing. Therefore, the estimation of CO2 degassing from the AZF is of great significance to understand the impact of the Earth’s natural degassing on the environment. Yang et al. (2028) studied three characteristic points in AZF for soil gas CO2 concentration measurements. The results showed that the mean and maximum values of soil gas CO2 concentration in the northern section of the ZMHF were much higher than those in the southern section of the ANHF. The soil gas CO2 geochemistry produced different anomalous features, which, combined with the degree of occlusion in the AZF, reflect the gas-bearing characteristics of the AZF as well as the seismic hazard. However, there is no research on the characteristics of AZF soil gas CO2 flux and its causes.

Taken together, the relationship between soil gas CO2 fluxes and active fault tectonic features is considered a hot topic for discussion in the field of geochemistry. In this paper, measurements of flux of soil gas CO2 were performed at 1,483 sampling points along 67 profile lines oriented perpendicularly to the AZF. The geochemical characteristics of soil gas CO2 emissions in study areas and the association between the earthquakes and degassing process were investigated. This investigation aims to evaluate the total output of the CO2 from the AZF and elucidate the tectonic association between the flux of soil gas CO2 and the activity of earthquakes.

2 Geological setting

The AZF is the main active rift on the eastern boundary of the Sichuan-Yunnan block (; Zhang et al., 2003). Under the combined effect of deep dynamic processes, large shear displacement and deformation zones are formed near the Xianshui fault, the ANHF and the Xiaojiang fault, which control the tectonic pattern of the Sichuan-Yunnan region (Figure 1).

FIGURE 1

The ANHF starts from Shimian in the north, extending through Xichang to Huili (Figure 2), with an overall northsouth trend. The overall NS trend is dominated by left-slip motion, with a maximum left-slip. The maximum left-slip rate since the Cenozoic is about 6.2 years−1 (), and the maximum extrusion rate is about 1.4 years−1 (Ran et al. 2008; ). The Shimian-Xichang area in the southern section of the ANHF was the main rupture site of the earthquake with MS 71/2 in 1,536, but no more significant earthquakes have occurred in the ANHF since an Ms 6 earthquake occurred south of Shimian in 1952 (Wen, 2000).

FIGURE 2

The Zemu River rift zone is connected with the ANHF in the north and the Xiaojiang fault in the south, extending from Xichang to Qiaojia (Figure 2), Yunnan via Puge and Ningnan, with a general strike of 330°. The Zemu River rift zone is dominated by left-slip motion, with a Holocene left-slip rate of about (6.4 ± 0.6) year−1 (Wen, 2000). Historically, there were earthquakes of an Ms 7 in Xichang in 1814 and Ms 7 in Dajiaoliangzi in 1850, and the most recent moderately strong earthquake was the Xichang Ms 5.1 earthquake on 31 October 2018.

3 Materials and methods

3.1 Site description

Soil gas flux CO2 surveys were performed at ANHF and ZMHF (26.8°N-29.2°N, 102°E – 103°E). During the layout of the CO2 flux measuring line, a large number of field investigations were carried out on the fault zone in the field. First, we found the location of the active fault fracture zone, and then the suitable measuring line was selected. A total of 67 measuring lines were arranged in ANHF and ZMHF, recording the longitude and latitude of the starting point and end point of the measuring line. Finally, the measuring lines designed for the investigation were drawn on the map indoors, according to the line density and feasibility and the final line scheme was obtained. Overall, sampling for soil gas CO2 flux was conducted at 1,483 sampling points along the 67 profiles that crossed the target faults (Figure 3). The measuring line is perpendicular to the fracture zone of the active fault zone and each line is generally more than 20 measuring points according to the landform. The spacing of measuring points near the active fault fracture zone is generally 10 m. After leaving the fracture zone, the spacing of measuring points was expanded by tens of meters to hundreds of meters according to the terrain (Figure 4).

FIGURE 3

FIGURE 4

3.2 Measurement methods of soil CO2 flux

The details and principles of the closed-chamber method were described by Zhou et al. (2016). Analytical error associated with a single measurement was about ±5% and the reproducibility was about ±10% for the range of 100–10,000 gm−2 d−1.

Soil temperature at a depth of 10 cm was taken at every site, using a pocket digital thermometer (Henxin AZ8821), which has an accuracy of ±0.1°C. To minimize the inherent variability of gas flux due to soil water and humidity, sampling was performed during dry stable weather. Ambient air temperature and barometric pressure were recorded every 25 samples. Each survey was executed over three consecutive days, and under rain-free conditions, to avoid environmental effects. Both soil temperature and atmospheric pressure, recorded at every measurement, were used for calibrating soil flux values. At each observation site, the wind speed and barometric pressure was measured.

3.3 Measurement methods of carbon isotope

Four samples of CO2 in soil gas from ANHF and four samples from ZMHF were collected for measurement of values of δ13CCO2(V-PDB). Glass bottles (500 mL) were used to collect soil gas in order to improve the accuracy of carbon isotopes. The gas samples were collected as follows: first, the glass bottle was filled with saturated brine and then inverted into a bucket filled with saturated brine and kept the glass bottle from tipping over. Soil gas is pumped into the glass bottle through a rubber hose. Until the filled saturated brine was replaced by gas in one-half of the bottle, the glass bottles were sealed in water with rubber blocks to avoid air contamination. All gas samples were analyzed at the Key Laboratory of Petroleum Resources Research, Institute of Geology and Geophysics, Chinese Academy of Sciences. Carbon isotope analysis was performed using a Delta Plus XL mass spectrometer. It is manufactured by ThermoFinnigan Inc., consisting of the HP6890 gas chromatograph, combustion/conversion oven, interface and DeltaPlusXP mass spectrometer (). The stable carbon isotope composition is generally expressed as δ13C and has a precision of ±0.2‰.

4 Results

Overall, 734 points along 32 profiles across the ANHF and 728 points along 35 profiles across the ZMHF were sampled for the CO2 flux analysis (Figure 5). There are two main patterns in the distribution of soil gas CO2 fluxes on the AZF, one with anomalously high values near the fractures, which are mainly because the fault fractures are well developed at this location, creating many gas escape routes. The other is where the anomaly occurs on both sides of the fault, because the fault is mature and its fractured core is well developed, while the surrounding fragmentation is so high that a peak occurs above the lateral fault. This bimodal distribution is also observed in the San Andreas Fault profile ().

FIGURE 5

In ANHF, the CO2 flux in soil gas varied in the range of 0–1996.2 g m2d−1, and the arithmetic mean values were 73.7 g ±2.36 m2d−1 (Figure 6); In ZMHF, the CO2 flux in soil gas varied in the range of 0–442.3 gm2d−1, and the arithmetic mean values were 87.5 g ±4.21 m2d−1. The δ13CCO2(V-PDB) were in the range of −25.0‰ ∼ −19.7‰ (Table 1). The statistics data on the fluxes of CO2 in soil gas from the AZF are listed in Table 2. The mean fluxes of CO2 in each soil gas survey line were in the range of 7.23–272.22 gm−2 d−1 (Figure 7).

FIGURE 6

TABLE 1

NameLine numberConcentration (%)δ13C-CO2 (‰)
Soilgas01ANH152.7−22.3
Soilgas02ANH204.8−23.9
Soilgas03ANH261.2−20.8
Soilgas04ANH271.2−22
Soilgas05ZMH031.7−25.6
Soilgas06ZMH061.3−22.8
Soilgas07ZMH200.7−19.7
Soilgas08ZMH343.1−25

Values of Soil gas CO2 carbon isotope.

TABLE 2

Line numberLength of survey lines(m)Number of measuring pointsMeanMaximumMinimum
(gm−2d−1)(gm−2d−1)(gm−2d−1)
ANH015552893.84201.128.83
ANH0242529125.46688.520
ANH034752874.59160.540
ANH045393073.68224.6418.75
ANH051,1462779.31266.4119.09
ANH063672354.28131.930
ANH071251334.75100.8313.43
ANH0842839114.6213.4943.42
ANH091,1492932.0376.180
ANH1010629272.221996.210
ANH118852937.09187.280
ANH1223511121.74288.9860.84
ANH134411481.89231.6419.99
ANH147561940.4677.5120.15
ANH15255257.2342.840.66
ANH168783952.05112.914.76
ANH172,2724262.47143.6511.77
ANH184722035.8662.458.29
ANH191,6953256.29239.4916.9
ANH204392152.48173.1323.6
ANH211,3392873.23173.010
ANH221,6773668.98169.3820.5
ANH237182065.25154.7712.05
ANH2419842173.77179.8720.61
ANH251,1302433.5100.8113.67
ANH26123877.34329.9522.36
ANH276631866.63144.6626.03
ANH2892813112.51196.1945
ANH297602359.37109.227.75
ANH305351892.92212.6134.27
ANH312961426.4457.20
ANH321731466.73135.0122.57
ZMH011,2042353.79124.7413.87
ZMH023671655.42101.6114.11
ZMH0370727113.35241.0626.62
ZMH043541846.8777.9223.38
ZMH051,0562588.52232.3624.32
ZMH06193235124.17414.350
ZMH0718162482.33229.427.45
ZMH081,2572260.45136.7416.11
ZMH091,47926112.87378.990
ZMH101,1241987.98197.5633.71
ZMH118441448.32107.2522.46
ZMH125741652.84133.819.98
ZMH139112246.58108.378.01
ZMH141,3632369.56140.154.8
ZMH151,1962092.94293.4120.97
ZMH164781539.38110.416.33
ZMH172351366.13176.3120.68
ZMH186281968.9160.5311.14
ZMH192951445.5189.6125.61
ZMH205592081.38236.274.78
ZMH211,2752789.63264.7211.11
ZMH223581682.61121.7626.49
ZMH238941985.71124.814.8
ZMH245852093.19244.6655.09
ZMH2539216109.49223.5541.42
ZMH267872099.74186.8150.55
ZMH271,0932084.78140.239.49
ZMH2892524117.03442.2657.96
ZMH293141555.2190.223.47
ZMH3061620143.79229.7637.73
ZMH319972286.73271.5923.28
ZMH3247417111.38329.6625.21
ZMH332,11627114.28353.4837.43
ZMH341,3902492.56276.1814.06
ZMH359652099.29257.7824.72

The value of each measurement line on AZF.

FIGURE 7

5 Discussion

5.1 Sources of CO2 degassing

The CO2 in subsurface fluids can be classified into two categories based on the type of parent material formed, organic and inorganic. The organic mainly formed by organic matter decomposition and bacterial activity, and the inorganic mainly from mantle/magmatic activity, thermal decomposition of carbonate rocks and dissolution of carbonate rocks (Wycherley et al., 1999). Different types of CO2 exhibit different carbon isotopic signatures. δ13CCO2 (V-PDB) vs. 1/CO2 is regarded to be an indicator for gas sources region (; Yuce et al., 2017; ). The values of end members for δ13CCO2 of mantle and crust are −6.5‰ and 0‰, respectively, while for atmospheric and biogenic end members are −7‰ and −26‰, respectively (; ; ).

From Figure 8, it can be seen that the values of δ13CCO2 values of the hot spring gas and soil gas samples are located in the composition mixing range between the mixed end of the crustal mantle and the biological end. The δ13CCO2 values of the soil gas samples are located between the hot spring gas and the end of biological, which indicates that the CO2 in the surface AZF mainly originates from the biological and the mixed end of the crustal mantle. Compared to the shallow depth of soil air circulation, the circulation of hot spring gas can reach depths of tens of kilometers. The CO2 in soil gas is closer to the end of biological, which indicates that the oxidation of organic matter during aerobic microbial respiration plays the most important role in the production of CO2 in soil gas.

FIGURE 8

Based on the above results, an attempt is made to explain the conceptual model of CO2 sources and transport in the active fault zone of AZF. Fluids from the crust and mantle that accumulate in the lower crust and upper mantle rise through deep faults with a small amount of crust- and mantle-derived gas diffusing into near-surface soil gas. Although biogenic CO2 is the primary source of soil gas, crust-derived and mantle-derived CO2 may be a secondary source. In addition, air can intrude into faulted soil gas due to fluctuations in air pressure (; ).

5.2 Contribution to the atmosphere from CO2 degassing in the AZHF

Global warming, which has received widespread social attention, is thought to be closely related to the rapid increase in atmospheric CO2 concentrations over the past 100 years (Joos et al., 1999; solomon et al.,2009; ). However, many scientists currently question the warming caused by increasing atmospheric CO2 concentrations, and the main debate focuses on what are the main drivers of warming; how accurate predictions of future climate trends based on existing climate models; and what is the magnitude of the impact of climate change (; Iqbal et al.,2009), all of which are subject to considerable uncertainty. There is a great deal of uncertainty in these issues. Therefore, it is necessary not only to reduce greenhouse gas emissions caused by human activities, but also to study the contribution of natural factors to atmospheric greenhouse gases, so as to distinguish between natural and anthropogenic factors leading to changes in atmospheric greenhouse gas concentrations, and to correctly understand the impact of Earth degassing on the increase of atmospheric greenhouse gas concentrations. In addition, the degassing areas of seismically active fracture zones are also “natural analogue” sites for studying the leakage of CO2 geological storage, especially the risk of sudden leakage of stored gases to the biosphere ().

The annual contribution of CO2 degassing from the fault zone to the atmosphere is equal to the amount of the annual average released flux each section of the fault zone multiplied by the area of each section of the fault zone. Based on the AZF and degassing characteristics, the AZF was divided into 5 segments, and the average value of the length and degassing flux of each segment was calculated based on the measurements (Table 3). The width of the fracture zone of the fault is generally tens of meters to hundreds of meters, and the surface avoidance zone of the active fault is more than 50 m (Xu et al., 2002). After the Wenchuan Ms8.0 earthquake, the width of the surface fracture zone of Longmenshan fault is about 200 m (Zhou et al., 2017), Several earthquakes of magnitude 6 or greater have occurred in the history of AZF, with the strongest one occurring north of Xichang with a magnitude of 7.5. The flux data from all sampling points on the AZF were superimposed on a single graph, which could be seen that the outliers with more than double the variance are concentrated within 200 m (Figure 8). High values of soil gas are usually exhibited near the fracture zone (Yuce et al., 2017; ), So the width of AZF was supposed as 200 m in this study.

TABLE 3

NumberNumber of survey linesNamelength (km)Size (106m2)Mean flux (gm-2d-1)emission (t·y-1)
116ANHF Northern Section601280.95354,561
216ANHF Southern Section651363.99303,632.6
313ZMHF Northern Section489.673.29256,808.1
412ZMHF Middle Section26.55.377.04149,033.9
510ZMHF Southern Section22.54.5100.48165,038.4
Total6722244.4395.751,229,074

CO2 emissions by segment.

By calculating the average flux of each survey line of ANHF and ZMHF, the annual contribution of ANHF and ZMHF zone to the atmosphere is 1.2 ±0.4 Mt (Table 3). The total flux of greenhouse gases from typical Cenozoic volcanic areas in Chinese Mainland to the atmosphere is about 8.13 Mt · a−1 equivalent to about 6% of the total greenhouse gas emissions caused by global (102 ∼ 103 Mt·y−1) volcanic activities (). The annual emission of AZF is relatively small compared with that of typical volcanoes in China. However, there are only 8 active volcanoes in China (), while there are more than 200 active fault zones in the intraplate area of Chinese Mainland (), among which the Tanlu fault zone, Altun fault zone, Haiyuan fault zone and Xianshuihe fault zone are more than 1000 km long. The total emission of CO2 from all active fault zones would be huge, which may be higher than that from active volcanoes in China. Only a few of researches were aimed to calculate carbon emissions from faults in the world (Table 4), but it is far from enough. From a public health risk assessment perspective, detailed CO2 degassing maps should also be performed to evaluate the potential for CO2 release from fault.

TABLE 4

LocationCountryAverage CO2 fluxMax CO2 fluxReferences
(gm−2d−1)(gm−2d−1)
ANHFChina72.471996.21This study
ZMHFChina82.93442.26This study
Dead Sea FaultTurkey26.755.4Yuce et al. (2017)
San Andreas FaultUnited States19631.40
Calaveras FaultUnited States56428
Arbia FaultItaly78204.98
Chukuo FaultChina1729

Study on soil gas CO2 fluxes in fracture zones in the world.

5.3 Diffuse degassing in the ANHF and ZMHF

It is well known that understanding of crustal dynamics and degassing processes can be improved by studying soil gas in fracture zones (). Soil gas CO2 flux is an important indicator of crustal stress accumulation and fault activity (; Zhou et al., 2010; Yuce et al., 2017). Changes in tectonic stress leads to expansion or contraction of bedrock, altering the degree of water-rock response and thus affecting CO2 degassing. The strength of degassing along active faults is related to seismic activity, which can reflect the state of stress accumulation to some extent (; Manaker et al.,2008). Hence, there could be some degree of correlation between CO2 emissions and historical earthquakes and nowadays earthquakes.

According to the location and magnitude of historical strong earthquakes on the fault zone, the seismic gap on the fault zone, and the spatial distribution of current seismic activity, the soil gas flux on the AZF has obvious segmented features (Yu et al., 2014; Yu et al., 2018; Wen, 2000; Zhou et al., 2017; ; ).

In this paper, our results show that high soil gas CO2 fluxes are found in the northern section of ANHF, southern section of AMHF and the central section of AZF (Figure 9), which well correlates with the sites that have high seismic activities. Notably, the higher soil gas concentration values in the northern section of the ANHF relative to the southern section of the ANHF, combined with the high flux values at point A (Figure 9), may indicate that the higher degree of fracture fragmentation in the area leads to elevated bedrock permeability, which promotes CO2 migration from deeper to the surface. B point (Figure 9) is located at the intersection of ANHF and ZMHF. On the one hand, the area has a high stress background value under the interaction of multiple sets of ruptures (Wen, 2000); on the other hand, the study shows that the area is highly occluded and belongs to the seismic gap area for both historical and present-day earthquakes (). Therefore, the fault fragmentation in this area is low, the fracture is not fully developed, and the overall soil gas CO2 flux values are in the middle. C point (Figure 9) is characterized by the same high seismic activity, but the soil gas CO2 flux is weaker than point A. This may be due to the “Y” shaped intersection of the ZMHF, Daliangshan fault, and Xiaojiang fault in the southern section, which may not have the tectonic conditions for high stress concentration. Although the frequency of seismic activity is high and the magnitude is large, the stress drop value of the earthquake source is generally in the middle. Therefore, the degree of fault fragmentation in this area is weak, resulting in relatively low bedrock permeability and fewer soil gas CO2 migration channels.

FIGURE 9

The higher degree of active fault fragmentation and the increased permeability of bedrock and soil, provide an upward channel for gas transport (). However, under the action of tectonic stress, the degree of water-rock response varies at different stages of deformation evolution, and the permeability within the fault changes temporally and spatially (; King, 1986; ; ; Annunziatellis et al., 2008). It is well known that faults can be both an upward pathway for transporting gases and can also block deep gas leaks (Caine et al., 1996; ; ). Moreover, fault permeability is subject to self-sealing processes that might limit or inhibit fluid gas motion (). The pore pressure of faulted soils usually depends on the crustal stresses associated with tectonic activity, which promotes the development of fractures and the formation of microfractures, resulting in the migration and redistribution of gases within the pore space (; Zhang et al., 2003). CO2 transport to the surface along well-developed fractures driven by pressure, which explains the higher soil gas CO2 fluxes are often found near faults.

6 Conceptual model of soil gas in the AZF

The cross-cutting and interaction of fracture zones of different scales on the AZF, with some strongly extruded passages forming blockages, leads to stress concentration, earthquake nucleation and breeding (). The geometric structure and activity habits of different sections of the fault zone is also controlled by the general dynamical background of the east boundary of the Sichuan-Yunnan active block (; Toda et al., 2008). This is reflected in the uncoordinated deformation movement of the tectonic units. This incoordination is probably the main reason for the type and size of earthquakes and the local stress distribution pattern. In discussing the seismic activity and stress distribution in this multi-crossing fault zone, it is necessary to consider the characteristics of the cross-faults (Brune et al., 1970; ). The interaction of cross-faults on the local stress field and seismic hazard needs to be considered when discussing the seismic activity and stress distribution characteristics in this multi-cross-fault region.

The northern section of the ANHF meets the Xianshuihe fault and the Daliangshan fault at an acute angle, which has local conditions of high stress concentration, high seismic activity, large magnitude, and high source stress, so the soil gas CO2 flux is the highest (Figure 10) (Wang et al., 1998; Toda et al., 2008); the ANHF and ZMHF meet at an obtuse angle, which has weak seismic activity and low source stress in general, so the soil gas CO2 flux is the lowest; the southern section of the AZF, ZMHF meets the Xiaojiang fault and the Daliangshan fault at an obtuse angle, which may not have high source stress (Xu et al., 2022) (Figure 10). The southern section of AZF, ZMHF and Xiaojiang fault and Daliangshan fault meet at an obtuse angle, so they may not have the tectonic condition of high stress concentration, and although the seismic activity is frequent, the source stress is relatively not high (Wu, 2020). Therefore, the overall performance of soil gas CO2 flux is at medium level.

FIGURE 10

It has been suggested that degassing features similar to the AZF can be found in worldwide geologic extensional contexts (), and spatial variations in CO2 emissions from deep sources may be attributed to tectonic activity (Irwin and Barnes, 1980). In addition, it has also been demonstrated that areas with enhanced crustal permeability may show precursors of CO2 degassing activity. Therefore, the geochemical exploration method proposed in this paper for AZF would be applied to other CO2 degassing areas on Earth (). Meanwhile, tectonic activity in possible future earthquakes can be monitored by implementing soil gas CO2 flux and concentration measurements on active faults.

7 Conclusion

In this study, geochemical characteristics and causes of the fluids released from the AZF were systematically analyzed, and the seismic activity of various segments of the AZF was discussed. The main conclusions are as follows:

  • (1) The investigation of carbon isotopic of underground fluid in AZF highlighted that: there are multiple sources of underground fluids in AZF; Soil gas CO2 was demonstrated to be mainly biogenesis in AZF, while the spring gas is mainly of crustal and mantle origin.

  • (2) Soil gas CO2 emissions surveys in this study and previous studies reveal a substantial range of average soil CO2 fluxes (1996.2 gm−2 d−1) for the AZF, which may have close affnities with the influence of regional fault systems and the seismic activity. Total soil CO2 output of the AZF is about 1.2 Mt a−1, suggesting high CO2 outgassing flux of the eastern boundary of the Chuan–Dian rhombic block. There are many longer fracture zones than AZF in the worldwide, so more attention should be given to CO2 emissions from faults.

  • (3) The variation of CO2 emissions in the study area is consistent with the strength of regional seismicity, and CO2 fluxes are also relatively low in areas with weaker tectonic stresses. These results indicate that seismicity and regional stress are the major triggering mechanism for CO2 emissions from the AZF.

Statements

Data availability statement

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

Author contributions

XZ, YL, and JD conceived the study. FL performed the literature review with support from YY, JT, JL, MH, OY, and SZ. KL, BY, ZZ, and YW participated in data processing. All authors contributed to interpretation and writing of the manuscript.

Funding

The work was funded by the National Key Research and Development Project (2018YFE0109700 and 2019YFC1509203), Central Public-interest Scientific Institution Basal Research Fund (CEAIEF2022030205, CEAIEF20220213, CEAIEF2022030200, and CEAIEF20220507), and the National Natural Science Foundation of China (41673106, 42073063, 4193000170, and U2039207), IGCP Project 724.

Acknowledgments

Special thanks are addressed to Chang Lu and Zhaofei Liu for the support during the fieldwork.

Conflict of interest

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

Publisher’s note

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

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/feart.2023.1117862/full#supplementary-material

References

  • 1

    AnnunziatellisA.BeaubienS. E.BigiS.CiotoliG.ColtellaM.LombardiS. (2008). Gas migration along fault systems and through the vadose zone in the Latera caldera (central Italy): Implications for CO2 geological storage[J]. Inter. J. Greenhouse Gas Cont.2 (3), 353372.

  • 2

    ApollaroC.DotsikaE.MariniL.BarcaD.BloiseA.De RosaR.et al (2012). Chemical and isotopic characterization of the thermo mineral water of Terme Sibarite springs (Northern Calabria, Italy). Geochem. J.46, 117129. 10.2343/geochemj.1.0166

  • 3

    BekaertD. V.TurnerS. J.BroadleyM. W.BarnesJ.HalldórssonS.LabidiJ.et al (2021). Subduction-driven volatile recycling: A global mass balance. Annu. Rev. Earth Planet. Sci.49, 3770. 10.1146/annurev-earth-071620-055024

  • 4

    BruneJ. N. (1970). Tectonic stress and the spectra of seismic shear waves from earthquakes[J]. J. Geophys. Res.75 (6), 49975009.

  • 5

    CamardaM.De GregorioS.Di MartinoR. M. R.FavaraR. (2016). Temporal and spatial correlations between soil CO2 flux and crustal stress[J]. J. Geophys. Res. Solid Earth121 (10), 70717085. 10.1002/2016jb013297

  • 6

    CaineJ. SEvansJ. P.ForsterC. B. (1996). Fault zone architecture and permeability structure[J]. Geology24(11), 10251028.

  • 7

    ChenZ.LiY.LiuZ.ZhengG.XuW.YanW.et al (2019). CH4 and CO2 emissions from mud volcanoes on the southern margin of the Junggar basin, NW China: Origin, output, and relation to regional tectonics[J]. J. Geophys. Res. Solid Earth124 (5), 50305044. 10.1029/2018jb016822

  • 8

    ChenZ.LiY.MartinelliG.LiuZ.LuC.ZhaoY. (2020). Spatial and temporal variations of CO2 emissions from the active fault zones in the capital area of China[J]. Appl. Geochem.112, 104489. 10.1016/j.apgeochem.2019.104489

  • 9

    ChiodiniG.CioniR.GuidiM.RacoB.MariniL. (1998). Soil CO2 flux measurements in volcanic and geothermal areas[J]. Appl. Geochem.13 (5), 543552. 10.1016/s0883-2927(97)00076-0

  • 10

    CiotoliG.SciarraA.RuggieroL.AnnunziatellisA.BigiS. (2016). Soil gas geochemical behaviour across buried and exposed faults during the 24 August 2016 central Italy earthquake. Ann. Geophys.59 (5), 1. 10.4401/ag-7242

  • 11

    DaiJ.SongY.DaiC.YanS.Chun-senD.Da-ruiW.et al (1996). Geochemistry and accumulation of carbon dioxide gases in China[J]. AAPG Bull.80 (10), 16151625. 10.1306/64EDA0D2-1724-11D7-8645000102C1865D

  • 12

    DaiJ. X. (1995). Abiogenic gas in oil-gas bearing basins in China and its reservoirs. Nat. Gas. Ind.15 (3), 2227.

  • 13

    DengQ.ZhangP.RanY.XiaopingY.WeieM.QuanzhiC.et al (2003). Basic characteristics of active tectonics of China[J]. Sci. China Ser. D Earth Sci.46 (4), 356372. 10.1360/03yd9032

  • 14

    DoğanT.SuminoH.NagaoK.NotsuK.TuncerM. K.ÇelikC. (2009). Adjacent releases of mantle helium and soil CO2 from active faults: Observations from the Marmara region of the North Anatolian Fault zone, Turkey[J]. Geochem. Geophys. Geosystems10 (11). 10.1029/2009gc002745

  • 15

    EtiopeG.BeneduceP.CalcaraM.FavaliP.FrugoniF.SchiattarellaM.et al (1999). Structural pattern and CO2–CH4 degassing of Ustica Island, Southern Tyrrhenian basin[J]. J. Volcanol. Geotherm. Res.88 (4), 291304. 10.1016/s0377-0273(99)00010-4

  • 16

    FischerT. P.AiuppaA. (2020). AGU centennial grand challenge: volcanoes and deep carbon global CO2 emissions from subaerial volcanism—recent progress and future challenges[J]. Geochem. Geophys. Geosystems21 (3), e2019GC008690. 10.1029/2019GC008690

  • 17

    FischerT. P.ArellanoS.CarnS. (2019). The emissions of CO2 and other volatiles from the world’s subaerial volcanoes[J]. Sci. Rep.9 (1), 111.

  • 18

    FuC. C.YangT. F.DuJ.WaliaV.ChenY. G.LiuT. K.et al (2008). Variations of helium and radon concentrations in soil gases from an active fault zone in southern Taiwan. Radiat. Meas.43, S348S352. 10.1016/j.radmeas.2008.03.035

  • 19

    GuoZ.ZhangM.ChengZ.ZhangL.LiuJ. (2014). Fluxes and Genesis of greenhouse gases emissions from typical volcanic fields in China. Acta Petrol. Sin.30 (11), 3463480. ( in Chinese).

  • 20

    HeH.OguchiT. (2008). Late Quaternary activity of the Zemuhe and Xiaojiang faults in southwest China from geomorphological mapping. Geomorphology96 (1-2), 6285. 10.1016/j.geomorph.2007.07.009

  • 21

    HiettC. D.NewellD. L.JessupM. J.GramblingT. A.ScottB. E.UpinH. E. (2022). Deep CO2 and N2 emissions from Peruvian hot springs: Stable isotopic constraints on volatile cycling in a flat-slab subduction zone[J]. Chem. Geol.595, 120787. 10.1016/j.chemgeo.2022.120787

  • 22

    HolubR. F.BradyB. T. (1981). The effect of stress on radon emanation from rock. J. Geophys. Res.86, 17761784. 10.1029/jb086ib03p01776

  • 23

    HuoC.HameedJ.SharifA.OhoudA.NoufA.Noor ul-ainB.et al (2022). Recent scenario and nexus of globalization to CO2 emissions: Evidence from wavelet and Quantile on Quantile Regression approach. Environ. Res., 113067. 10.1016/j.envres.2022.113067

  • 24

    InguaggiatoS.MazotA.DilibertoI. S.InguaggiatoC.MadoniaP.RouwetD.et al (2012). Total CO2 output from Vulcano island (Aeolian Islands, Italy)[J]. Geochem. Geophys. Geosystems13 (2). 10.1029/2011gc003920

  • 25

    IqbalJ.LinS.HuR.FengM. (2009). Temporal variability of soil-atmospheric CO2 and CH4 fluxes from different land uses in mid-subtropical China[J]. Atmosph. Environ.43(37), 58655875.

  • 26

    IrwinW. P.BarnesI. (1980). Tectonic relations of carbon dioxide discharges and earthquakes. J. Geophys. Res.85, 31153121. 10.1029/JB085iB06p03115

  • 27

    IssonT. T.PlanavskyN. J.CooganL. A.et al (2020). Evolution of the global carbon cycle and climate regulation on Earth[J]. Glob. Biogeochem. Cycles34 (2), 128. 10.1029/2018GB006061

  • 28

    ItalianoF.BonfantiP.DittaM.PetriniR.SlejkoF., Helium and carbon isotopes in the dissolved gases of Friuli Region (NE Italy): Geochemical evidence of CO2 production and degassing over a seismically active area, Chem. Geol.,266,Issues 76–85, 2009, Pages 7685. ISSN 0009-2541, 10.1016/j.chemgeo.2009.05.022

  • 29

    ItalianoF.BonfantiP.PizzinoL.QuattrocchiF. (2010). Geochemistry of fluids discharged over the seismic area of the Southern Apennines (Calabria region, Southern Italy): Implications for Fluid-Fault relationshipsfluids discharged over the seismic area of the southern apennines (Calabria region,southern Italy): Implications for fluid-fault relationships. Appl. Geochem.25, 540554. 10.1016/j.apgeochem.2010.01.011

  • 30

    Jacome-PazM. P.Gonzalez-RomoI. A.Prol-LedesmaR. M.Torres VeraM.Pérez-ZárateD.Rodríguez-DíazA.et al (2020). Multivariate analysis of CO2, H2S and CH4 diffuse degassing and correlation with fault systems in Agua Caliente-Tzitzio, Michoacan, Mexico[J]. J. Volcanol. Geotherm. Res.394, 106808. 10.1016/j.jvolgeores.2020.106808

  • 31

    JiangG. Y.XuX. W.ChenY. H.LiuY. J.FukahataY.WangH.et al (2015). Geodetic imaging of potential seismogenic asperities on the Xianshuihe-Anninghe-Zemuhe fault system, southwest China, with a new 3-D viscoelastic interseismic coupling model. J. Geophys. Res.120, 18551873. 10.1002/2014JB011492

  • 32

    JingJ.TangZ.YangY.MaL. (2019). Impact of formation slope and fault on CO2 storage efficiency and containment at the Shenhua CO2 geological storage site in the Ordos Basin, China[J]. Int. J. Greenh. Gas Control88, 209225. 10.1016/j.ijggc.2019.06.013

  • 33

    JolieE.MoeckI.FauldsJ. E. (2015). Quantitative structural–geological exploration of fault-controlled geothermal systems—A case study from the Basin-and-Range Province, Nevada (USA). Geothermics54, 5467. 10.1016/j.geothermics.2014.10.003

  • 34

    JoosF.PlattnerG. K.StockerT. F.MarchalO.SchmittnerA. (1999). Global warming and marine carbon cycle feedbacks on future atmospheric CO2[J]. Sci.284 (5413), 464467.

  • 35

    KämpfH.BräuerK.SchumannJ.HahneK.StrauchG. (2013). CO2 discharge in an active, non-volcanic continental rift area (Czech Republic): characterisation (δ13C, 3He/4He) and quantification of diffuse and vent CO2 emissions[J]. Chem. Geol.339, 7183. 10.1016/j.chemgeo.2012.08.005

  • 36

    KangY. J.YunS. T.YuS.DoH. K.ChaeG. (2020). Quantitative assessment of deep-seated CO2 leakage around CO2-rich springs with low soil CO2 efflux using end-member mixing analysis and carbon isotopes[J]. J. Environ. Manag.276, 111333. 10.1016/j.jenvman.2020.111333

  • 37

    KennedyB. M.KharakaY. K.EvansW. C.EllwoodA.DePaoloD. J.ThordsenJ.et al (1997). Mantle Fluids in the San Andreas Fault System, California. Science278 (5341), 12781281. 10.1126/science.278.5341.1278

  • 38

    KingC. Y. (1986). Gas geochemistry applied to earthquake prediction: An overview[J]. J. Geophys. Research: Solid Earth91 (B12), 1226912281. 10.1038/ngeo2622

  • 39

    LeeH.MuirheadJ. D.FischerT. P.EbingerC. J.KattenhornS. A.SharpZ. D.et al (2016). Massive and prolonged deep carbon emissions associated with continental rifting. Nat. Geosci.9 (2), 145149. 10.1038/ngeo2622

  • 40

    LiY. C.ShanX. J.QuC. Y.ZhangY. F.SongX. G.JiangY.et al (2017). Elastic block and strain modeling of GPS data around the Haiyuan-Liupanshan fault, northeastern Tibetan Plateau. J. Asian Earth Sci.150, 8797. 10.1016/j.jseaes.2017.10.010

  • 41

    LiuJ. (1999). Volcanoes hazards and montoring[J]. Quadrangular Res. (05), 414422. (in Chinese abstract).

  • 42

    ManakerD. M.CalaisE.FreedA. M.AliS. T.PrzybylskiP.MattioliG.et al (2008). Interseismic plate coupling and strain partitioning in the northeastern Caribbean[J]. Chem. Geol.469, 176184. 10.1016/j.chemgeo.2017.01.006

  • 43

    MartinelliG.DadomoA. (2017). Factors constraining the geographic distribution of earthquake geochemical and fluid-related precursors. Chem. Geol.469, 176184. 10.1016/j.chemgeo.2017.01.006

  • 44

    MuirheadJ. D.FischerT. P.OlivaS. J.LaizerA.van WijkJ.CurrieC. A.et al (2020). Displaced cratonic mantle concentrates deep carbon during continental rifting. Nature582 (7810), 6772. 10.1038/s41586-020-2328-3

  • 45

    NeriG.BarberiG.OlivaG.OrecchioB.PrestiD. (2006). A Possible Seismic Gap within a Highly Seismogenic Belt Crossing Calabria and Eastern Sicily, Italy. Italy. Bull. Seismol. Soc. Am.96, 13211331. 10.1785/0120050170

  • 46

    PeiX. Y.WangX. M.ZhangC. G. (1998). Basic character of segmentation of the Quaternary movement on the Anninghe fault[J]. Earthq. Res. Sichuan4, 5261.

  • 47

    RahillyK. E.FischerT. P. (2021). Total diffuse CO2 flux from Yellowstone caldera incorporating high CO2 emissions from cold degassing sites[J]. J. Volcanol. Geotherm. Res.419, 107383. 10.1016/j.jvolgeores.2021.107383

  • 48

    RanY. K.ChengJ. W.GongH. L.ChenL. (2008). Late Quaternary geomorphic deformation and displacement rates of the Anninghe fault around Zimakua[J]. Seismology and Geology30 (1), 8698.

  • 49

    RenZ.LinA. (2010). Deformation characteristics of co-seismic surface ruptures produced by the 1850 M 7.5 Xichang earthquake on the eastern margin of the Tibetan Plateau. J. Asian Earth Sci.38 (1-2), 113. 10.1016/j.jseaes.2009.12.008

  • 50

    RoviraP.VallejoV. (2008). Changes in δ13C composition of soil carbonates driven by organic matter decomposition in a Mediterranean climate: a field incubation experiment. Geoderma144 (3), 517534. 10.1016/j.geoderma.2008.01.006

  • 51

    SanoY.MartyB. (1995). Origin of carbon in fumarolic gas from island arcs. Chem. Geol.119 (1–4), 265274. 10.1016/0009-2541(94)00097-r

  • 52

    ShenZ. K.J.WangM.BürgmannR. (2005). Contemporary crustal deformation around the southeast borderland of the Tibetan Plateau. J. Geophys. Res. Solid Earth110. 10.1029/2004JB003421

  • 53

    SizovaE.HauzenbergerC.FritzH.FaryadS. W.GeryaT. (2019). Late orogenic heating of (ultra)high pressure rocks: Slab rollback vs. slab breakoff. Geosciences9, 499. 10.3390/geosciences9120499

  • 54

    SolomonS.PlattnerG. K.KnuttiR.FriedlingsteinP. (2009). Irreversible climate change due to carbon dioxide emissions[J]. Proc. Natl. Acad. Sci.106 (6), 17041709.

  • 55

    TamirG.ShenkerM.HellerH.BloomP. R.FineP.BartalA. (2012). Dissolution and recrystallization processes of active calcium carbonate in soil developed on tufa. Soil. Sci. Soc. Am. J.76 (5), 16061613. 10.2136/sssaj2012.0041

  • 56

    TodaS.LinJ.MeghraouiM.SteinR. S. (2008). May 2008 M = 7.9 Wenchuan, China, earthquake calculated to increase failure stress and seismicity rate on three major fault systems[J]. Geophy. Res. Let.35 (17).

  • 57

    ViveirosF.MarcosM.FariaC.JoãoL. G.CatarinaS.TeresaF.et al (2017). Soil CO2 degassing path along volcano-tectonic structures in the Pico-Faial-São Jorge islands (Azores archipelago, Portugal)[J]. Front. Earth Sci.5, 50. 10.3389/feart.2017.00050

  • 58

    WaliaV.LinJ.FuC.YangF.HongW.WenL.et al2010. Soil-gas monitoring: a tool for fault delineation studies along Hsinhua Fault (Tainan),Southern Taiwan. Appl. Geochem.25 (4), 602607. 10.1016/j.apgeochem.2010.01.017

  • 59

    WangQ. L.WangW. P.LiangW. F.ChenB. (1998). Stress-dissipative heat geotemperature precursor mechanism of earthquakes[J]. Acta Seismologica Sinica11, 615621.

  • 60

    WeinlichF. H.FaberE.BouškováA.HorálekJ.TeschnerM.PoggenburgJ. (2006). Seismically induced variations in Mariánské Lázně fault gas composition in the NW Bohemian swarm quake region, Czech Republic — A continuous gas monitoring. Tectonophysics421 (1-2), 89110. 10.1016/j.tecto.2006.04.012

  • 61

    WenX. Z. (2000). Character of rupture segmentation of the Xianshuihe-Anninghe-Zemuhe fault zone, Western Sichuan. Seismol. Geol.22 (3), 239249. 10.3969/j.issn.0253-4967.2000.03.005

  • 62

    WycherleyH.FleetA.ShawH. (1999). Some observations on the origins of large volumes of carbon dioxide accumulations in sedimentary basins. Mar. Petroleum Geol.16 (6), 489494. 10.1016/s0264-8172(99)00047-1

  • 63

    XuS.GuanL.ZhangM.ZhongJ.LiuW.XieX.et al (2022). Degassing of deep-sourced CO2 from Xianshuihe-Anninghe fault zones in the eastern Tibetan Plateau[J]. Sci. China Earth Sci.65 (1), 139155. 10.1007/s11430-021-9810-x

  • 64

    XuX.YuG.MaW.YongkangR.YouH.-C.ZhangL.-F.et al (2002). Evidence and methods for determining the safety distance from the potential earthquake surface rupture on active fault[J]. Seismol. EGOLOGY24 (4), 470.

  • 65

    YangY.LiY.GuanZ.ChenZ.ZhangL.LvC. J.et al (2018). Correlations between the radon concentrations in soil gas and the activity of the Anninghe and the Zemuhe faults in Sichuan, southwestern of China. Appl. Geochem.89, 2333. 10.1016/j.apgeochem.2017.11.006

  • 66

    YuX.YangC.DongG.ZhanD. (2014). Drilling and completion technology in the first geothermal well construction in Jixian system in Tangshan. Explor. Eng. Media41 (2), 4144.

  • 67

    YuceG.FuC.DʼAlessandroW.GulbayA.WaliaV.BellomoS.et al (2017). Geochemical characteristics of soil radon and carbon dioxide within the Dead Sea fault and Karasu fault in the Amik basin (Hatay), Turkey. Chem. Geol.469, 129146. 10.1016/j.chemgeo.2017.01.003

  • 68

    ZhangM.XuS.ZhouX.CaracausiA.SanoY.GuoZ.et al (2021). Deciphering a mantle degassing transect related with India-Asia continental convergence from the perspective of volatile origin and outgassing. Geochimica Cosmochimica Acta310, 6178. 10.1016/j.gca.2021.07.010

  • 69

    ZhangP.DengQ.ZhangG. (2003). Active tectonic blocks and strong earthquakes in the continent of China. Sci. China Ser. D Earth Sci.46 (2), 1324. 10.1360/03dz0002

  • 70

    ZhouX.SunF.ChenZ. (2017). Degassing of CO2, CH4, Rn and Hg in the rupture zones produced by Wenchuan Ms 8.0 earthquake. Acta Petrol. Sin.33 (1), 291303. (in Chinese with English abstract). Web of Science Id WOS:000393257900023 .

  • 71

    ZhouX. C.DuJ. G.ChenZ.ChengJ. W.TangY.YangL. M.et al (2010). Geochemistry of soil gas in the seismic fault zone produced by the Wenchuan Ms 8.0 earthquake, southwestern China. Geochem. Trans.11, 5. 10.1186/1467-4866-11-5

  • 72

    ZhouX.ChenZ.CuiY. (2016). Environmental impact of CO2, Rn, Hg degassing from the rupture zones produced by Wenchuan Ms 8.0 earthquake in Western Sichuan,China. Environ. Geochem. Health38 (5), 10671082. 10.1007/s10653-015-9773-1

Summary

Keywords

soilgas, carbon dioxide, geochemistry, diffuse degassing, AZF

Citation

Liu F, Zhou X, Dong J, Yan Y, Tian J, Li J, Ouyang S, He M, Liu K, Yao B, Wang Y, Zeng Z and Zhang Y (2023) Soil gas CO2 emissions from active faults: a case study from the Anninghe—Zemuhe fault, Southeastern Tibetan Plateau, China. Front. Earth Sci. 11:1117862. doi: 10.3389/feart.2023.1117862

Received

06 December 2022

Accepted

11 April 2023

Published

09 May 2023

Volume

11 - 2023

Edited by

Ziming Yang, Oakland University, United States

Reviewed by

Nunzia Voltattorni, Istituto Nazionale di Geofisica e Vulcanologia (INGV), Italy

Kan Li, Woods Hole Oceanographic Institution, United States

Updates

Copyright

*Correspondence: Xiaocheng Zhou, ; Jinyuan Dong,

Disclaimer

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

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics