Abstract

The sources of soil gases in the Tangshan fault zone, North China, were discussed, based on the soil gas compositions and isotopic ratios obtained by measurement in the field and sample analysis in the laboratory. Soil gas compositions and isotopic ratios indicate that air (A) end-member, limestone (L) end-member, and sediment (S) end-member are the major end-member components contributing to the soil gas in our study area, with fractional contributions in the range of 2-15 vol.%, 23-36 vol.%, and 62-65 vol.%, respectively, to CO2 from the gas wells. According to the relationship among the 3He/4He, average CO2 concentration, and He concentration of soil gas, the deepest depth the fault cut downward and the most developed fractures in the segment where the Heibeiligong (HBLG) well located were inferred, and the shallower depth the fault cut downward and the more developed fractures in the fault segments where the Weifengshan (WFS), Siwangzhuang (SWZ), Tianjingyice (TJYC), and Douhetai (DHT) wells located were inferred. Significant variations in CO2 concentration were observed in soil gases sampled in DHT, HBLG, and WFS soil gas wells in concomitance with a local seismic sequence by 2018 confirming for the first time a possible source of carbon dioxide generated in underlying limestones.

1. Introduction

Earthquakes are caused by tectonic evolution accompanied by matter and energy exchange deep inside the earth [13], which is mainly transmitted by fluid being released through active faults and fractures at different depths [4], owing to that fault and fracture are preferential migration pathways for gases (CO2, Rn, He, etc.) in the deep crust to migrate upward to the surface, due to their enhanced permeability and porosity relative to the surrounding rocks [57]. General overviews of the geochemical, structural, and seismic features in tectonically active areas have shown some evidence of correlation between soil gas geochemistry anomalies and tectonic activities, and soil gas discharge through fault and fracture in the active fault zones can be enhanced by fault and earthquake activity (Wakita et al., 1980; Toutain et al., 1992; [811]). As such, analyzing the origins and geochemical variations of soil gas in the seismically active areas could be a potential way to study faults and seismic activity [1214]. At present, soil gas concentration and flux surveys along active fault zones have been widely undertaken for earthquake research and prediction (Caracausi et al., 2003; [11, 15]).

Recent studies of soil gas geochemistry within fault zones and its relationship with earthquake activity further indicate that it can potentially be used to monitor earthquakes [10, 11, 16, 17]. Soil gas geochemical surveys in seismically active areas have been carried out across the world [5, 6, 18, 19], and some anomalies appearing before earthquakes have been identified [11, 20]. In the Arax Basin (Armenia), radon volume activities have been observed to vary before and after earthquakes [21]. Similarly, anomalously high radon volume activity fluctuations were observed several hours to a few days before an earthquake () that occurred in northern Taiwan [7].

Soil gas anomalies are usually complex as they are subject to multiple influences (atmospheric, biogenic, organic, and from the deep crust and mantle; [11, 2224]) and occasionally their relationships with seismic events may be ambiguous [4, 7, 25, 26]. Accordingly, if the observation of soil gas in fault zones is to become a useful method of fault activity assessment and seismic forecasting, further analysis of their relationships is required and the source identification for soil gas in the seismically active areas is a prerequisite.

The Tangshan fault is the main shock fault in which occurred the 1976 Tangshan 7.8 earthquake, one of the most destructive earthquakes ever recorded and caused 240000 deaths and property losses of over $ 1.4 billion [27]. In this region, earthquakes have occurred with increasing frequency (more than 20 earthquakes of per year since 2016) and high concentrations of soil gas components have been observed [5, 28].

In the present study, we aim to (1) investigate the origins of soil gas emitted from the Tangshan fault zone based on gas compositions and isotope ratio values, (2) analyze correlations between variations of soil gas component concentrations and earthquakes in the Tangshan fault zone, and (3) identify methods and components that are potential for earthquake monitoring in the Tangshan area of Northern China.

2. Geological Setting

The Tangshan fault zone is located in a fold depression region where the Northern China and Yanshan uplifts converge. It is considered to be a strain accumulation structure comprising three active faults, which are as follows, from north to south: Douhe fault (F1), Weishan-Changshan fault (F2), and Tangshan-Guye fault (F3). All three are normal faults with right lateral strike-slip and have NE trending (25–35°) and NW dipping with high dip angles (70–80°; Liu et al., 2011). The 1976 Tangshan MS 7.8 earthquake occurred in the middle segment of the Tangshan fault zone and was followed by 54 aftershocks of (Figure 1). This was interpreted as a complex tectonic process of strike-slip, thrust, and normal-fault events caused by detachment, fault propagation in the middle crust, and stress accumulation in the base of the upper crust [14].

The basement of the Tangshan area is composed of pre-Sinian metamorphic rocks with varying metamorphic grades whose outcrops rarely occur in the study area covered by the thick sediments of the Holocene (Q4) and upper Pleistocene (Q3), but Ordovician, Cambrian, and Proterozoic strata consisting mainly of limestone are widely exposed in the northern segment of the Tangshan fault zone, with δ13 C (PDB) and δ18 O (PDB) in the range of -0.65–1.40‰ and -12.34–7.47‰, respectively [29]; Mesozoic and Tertiary sediments are missing in the study area, and the thickness of Pleistocene and Holocene sediments increases from about 100 m in the north to 800 m in the south [14, 30, 31]. The total carbon contents of soil in the northern and southern segments of the Tangshan fault zone were 1.63% and 0.85%, respectively, and the Ra contents were 27.2 Bq kg-1 and 20.4 Bq kg-1, respectively [32]. Geological, seismic, and electrical surveys show that the active faults have dislocated the strata of the Cenozoic Era [14, 33]. The target area for the soil gas measurement in the field was about 2100 km, with the length of about 70 km and width of about 30 km, covered the three active faults and seven soil gas wells (Figure 1).

3. Methods

3.1. Soil Gas Measurements

Soil gas measurements were performed in the field at 388 measurement points covering the Tangshan fault zone (Figure 1). The CO2 concentration and Rn (radon) activity were measured at approximately 1–1.5 km intervals between every two measurement points. In order to avoid possible meteorological effects on the soil gas concentrations (Hinkle, 1994), the soil gas survey was carried out during a period of stable meteorological conditions from 8 April to 7 May 2010. At this time (spring), the air temperature () and wind velocity () remain relatively stable and the mean rainfall over the study period (approximately 1 month) is 21.3 mm, which is relatively low given that the mean monthly rainfall for the Tangshan area is 50 mm (http://data.cma.cn). The CO2 concentrations and Rn volume activity were measured by inserting a stainless steel sampling tube with a 3 cm diameter into the ground to a depth of 0.8 m [34]. The sampling tube was connected by rubber tubes to a radon detector (RAD 7, for Rn volume activity) and a portable infrared CO2 monitor (GXH-3010-E, for CO2 concentration; Figure 2). Rn activity was measured 15 min after sampling (the time necessary for Po and Rn nuclei to reach equilibrium, which is about five times the half-life of 218Po). The detection limit and measurement error of the RAD 7 radon detector were 3.7 Bq m-3 and ±5%, respectively. An inlet filter and molecular sieve were used to protect the detector from dust and soil moisture. The detection limit and measurement error of the GXH 3010-E CO2 monitor were 0.01% and ±2%, respectively. An inlet filter was used to protect the detector from dust.

3.2. Soil Gas Measurement in Selected Wells

Seven sites were selected for soil gas measurements and equipped with shallow depth wells in November 2017 (Figure 1). These sites were selected because they are characterized by (1) relatively high CO2 concentration and Rn activity, (2) from the fault, and (3) a lack of possible artificial sources of contamination. An inverted polytetrafluoroethylene (PTFE) circular accumulation hemispherical chamber with a volume of and a radius of 0.2 m was fixed 5.0 m below the ground. A 6-channel deconcentrator was installed on the inner wall of the chamber to reinject the circulating gas in order to ensure the immediate and homogeneous mixture of gas in the chamber. Two exhaust tubes (PTFE, 5.0 m in length, 10 cm in diameter) were used, with their lower ends connected to the top of the chamber and their other ends connected to the inlet and outlet of the detector. These devices form a closed circuit for the analysis of gases present in the first five meters of sediments (Figure 3). Gas measurements were performed repeatedly in the field at the end of each month from June to November 2018. Gas samples were also collected from the outlet of the apparatus following the methods described by Chen et al. [5] and Zhang et al. [35] in September 2018.

3.3. Laboratory Analysis

Gas compositions (CO2, CH4, Ar, H2S, N2, and O2) and isotopic ratios (δ13, 3He/4He (R/Ra), CO2/3He, and 4He/20Ne) were all measured in the Laboratory of Gas Geochemistry (Lanzhou, China), Institute of Geology and Geophysics, Chinese Academy of Sciences. Gas compositions were determined by a mass spectrometer (MAT 271) and characterized by relative standard deviations of <5%. The 3He/4He (reported as R/Ra, ) and 4He/20Ne values were determined by a different mass spectrometer (VG 5400), while δ13 C values were analyzed by a gas chromatography-pyrolysis-isotope ratio mass spectrometer (HP 6890-Delta Plus XL) with uncertainties of ±0.3‰. Analysis of all the samples was completed from December 5 to December 7, 2018, within 10 days of sampling. Air on the Gaolan Hill in the south of Lanzhou was used to calibrate the instrument [35].

3.4. Population Identification and Spatial Distribution Analysis

The graphical statistical method (GSA) of Sinclair [36] was used as described by Chiodini et al. [22] to distinguish different populations of soil gas from the Tangshan fault zone. This method involves plotting CO2 concentrations or Rn volume activity values against a cumulative probability and identifying inflection points within each dataset. The inflection points signify partitions among different statistical populations, which may represent discrete sources of soil gas [37]. The data obtained at the mobile measurement points were processed using the Kriging interpolation method in the SURFER 8.0 package; then, the contour maps were built with the data set to analyze the spatial distributions of CO2 concentrations and Rn volume activity. Kriging is a general approach for stochastic spatial interpolation in which the continuous regionalized variable of interest , sometimes referred to as the primary variable, is predicted at any unsampled location of the study area using the values of measured at different locations; , the predicted value at location , , is calculated as an affine linear combination of the observed values in such a way that it is a best linear unbiased prediction predictor [38].

4. Results

The CO2 concentrations and Rn activities obtained at the 388 measurement points were in the ranges of 514.0–10318.1 ppm and 882.4–17283.5 Bq m-3, respectively, with means of 2438.6 ppm and 4730.5 Bq m-3, respectively (Table 1).

The geochemical parameters of soil gases sampled in the seven shallow depth wells (DHT, TJYC, HBLG, LHD (Luhuadai), HET (Haiertun), SWZ, and WFS) in the Tangshan fault zone are listed in Table 2. The concentrations of CO2, CH4, Ar, H2S, N2, and O2 in soil gas sampled in the seven constructed shallow depth wells were in the ranges of 0.09–6.23%, 31–113 ppm, 0.92–1.15%, 0–631 ppm, 78.64–97.81%, and 0.45–19.90%, respectively. The Rn volume activities were in the ranges of 4.59–33.00 k Bq m-3. The atmospheric concentrations of CO2 and He in the study area were 0.04% and 5.35 ppm, respectively.

The 3He/4He (R/Ra), CO2/3He, and 4He/20Ne ratios of soil gas sampled in the seven shallow depth wells varied from 0.99 to 1.13, to , and 0.28 to 0.31, respectively. The δ13 C values of CO2 were in the range of −19.9 to −18.9‰. The 3He/4He (R/Ra), δ13, CO2/3He, and 4He/20Ne ratios of air in the study area were 0.99, −8.3‰, , and 0.29, respectively (Table 2).

5. Discussion

5.1. Sources of Soil Gas Sampled in the Tangshan Fault Zone

The probability plots of both CO2 concentration and Rn activity of soil gas collected at the 388 sampling points displayed three-segment distributions, characterized by the partial overlap of a low population (A), a medium population (B), and a high population (C; Figures 4 and 5). There were two inflection points at CO2 concentrations of 3692 ppm and 8269 ppm, indicated by arrows in the figures, which delineate the fractions of 80% for Population A, 16% for Population B, and 4% for Population C (Figure 4). There were also two inflection points at Rn volume activities of 7588.2 Bq m-3 and 13749.8 Bq m-3, delineating the fractions of 83% for Population A, 12% for Population B, and 5% for Population C (Figure 5). The means and 95% confidence intervals of CO2 concentration and Rn volume activity and the fractions of the three populations listed in Table 3 were determined using the GSA method [39].

Contour maps of soil gas CO2 concentrations and of Rn activities found in the 388 mobile measurement points were mapped in Figures 4 and 5. A consistent spatial coincidence of CO2 concentrations and of Rn activities was found (Figures 6 and 7), which may indicate common sources of CO2 and Rn, as CO2 can act as a carrier gas for Rn [4, 40]. In addition, almost all the CO2 and Rn values from Population C are concentrated near the fault zone; in contrast, almost all CO2 and Rn values from Population A are distributed outside the fault zone. Moreover, most of the CO2 and Rn values from Population B are located near the fault zone but with some outside it (Figures 6 and 7). Gases migrating from deep geological formations towards the surface through faults are usually characterized by enrichment in one or more components compared to atmospheric composition and are expected in fault zones. Outside the fault zones, gases are diffused in the superficial soil and show lower concentrations of nonatmospheric components [4143]. Hence, it could be inferred that the CO2 and Rn of Population A come from the superficial soil, while those of Population C uprise from deeper locations through faults. Meanwhile, the CO2 and Rn of Population B could have a complex source due to biogenic, thermogenic, mechanochemical, or rock weathering origins [9, 4446]. Accordingly, further investigations into eventual geochemical variations over time should be carried out chiefly in soil gases from Population C since they are more representative of the deep originating component.

Therefore, attention has been focused on identifying the origin of soil gases of Population C and seven soil gas wells with depths of 5.0 m were constructed along faults in November 2017 for soil gas measurement and sampling (Figures 1 and 3).

As shown in Tables 1 and 2, the CO2 concentrations (0.12–6.23%) and Rn activities (4.59–33.00 kBq m-3) of soil gas from the seven soil gas wells were much higher than those from the 388 mobile measurement points (514.0–10318.1 ppm for CO2 concentration and 882.4–17283.5 Bq m-3 for Rn volume activity), suggesting that the new wells were relatively less exposed to dilution processes of gases by the superficial soil gas and by the atmosphere.

The results of the total gas component analysis showed that the major chemical species of the soil gas samples were N2, O2, CO2, and Ar. The He concentrations in the soil gas samples from the well were slightly higher than those of the atmosphere except for that of DHT (Table 2), according to the results of 3He/4He (R/Ra) and 4He/20Ne analyses (Figure 8), which could indicate the minority contribution of gases upwelling from the deep earth to the soil gases of the Tangshan fault, although the soil gases could mainly derive from the atmosphere.

However, the CO2 concentrations of gas samples from the seven new soil gas wells were much higher than the atmospheric concentration (Table 2), suggesting that there were other sources of CO2. As reported in previous studies, underground CO2 is a crucial link in the interactions among gas, water, rocks, and biological activities [46]. The two main origins of underground CO2 are the microbial degradation of organic matter in sediment and limestone weathering, including water-limestone interactions and limestone decomposition [44, 45, 47, 48]. Thus, an approach analogous to that used for the identification and quantification of CO2 in volcanic and geothermal regions was adopted [35, 49, 50]. The mantle (M) end-member was replaced with an air (A) end-member (Figure 9); then, three major end-member components involving an air (A) end-member, limestone (L) end-member, and sediment (S) end-member were selected to identify and quantify the various sources contributing to CO2 from the seven new wells. The δ13 C value of −8.3‰ and the CO2/3He ratio of for the (A) end-member were obtained by measuring air samples (Table 1). For the (L) and (S) end-members, δ13 C values of 0‰ and −30‰, respectively, with a corresponding CO2/3He ratio of , were suggested by Sano and Marty [49]. The fractions contributed by the (A), (L), and (S) end-members to the CO2 at the seven new wells were estimated using the following mass balance equations: where  = observed, and is the fraction contributed by the (A), (L), and (S) end-members to the CO2 at the seven soil gas wells.

The results showed that 62–65% of the CO2 from the seven wells was derived from the (S) end-member, 23–36% was derived from the (L) end-member, and the remaining 2–15% came from the (A) end-member (Figure 9 and Table 4).

This suggests that the major sources of CO2 in the Tangshan fault zone are the microbial degradation of organic matter and the decomposition of carbonatic rocks in sediments. Furthermore, the higher N2/O2 ratios, together with the strong correlation between the depleted O2 concentrations (compared to the atmosphere) and the CO2 concentrations in the wells (correlation coefficient of 0.90; Table 2 and Figure 10), indicate that CO2 is mainly contributed by the microbial oxidation of organic matter in the sediments of the Tangshan fault zone. In particular, enriched N2 (93.65% and 97.81%, respectively) and depleted O2 (0.45% and 4.55%, respectively) were observed in soil gas wells LHD and HET (Table 2), which are located in the farmland. This should be attributed to the acceleration of nitrifying and denitrifying bacterial activity caused by artificial fertilization [51, 52].

The Tangshan fault zone is located in the well-known Bohai Bay sedimentary basin, where limestones are widely distributed in Ordovician, Cambrian, and Proterozoic strata in the deep earth [30, 53]. These limestones could release CO2 at high concentrations during its decomposition and reaction with groundwater and CO2, as per Eq. (2) [5457].

Above all, it should be noted that the decomposition of carbonates in limestone strata could be another important potential source for the CO2 in the wells at the Tangshan fault zone.

In addition, it was shown in the relationship between 3He/4He and the average values of CO2 concentration and He concentration of soil gas from the seven soil gas wells that there were three distinguishing populations (A, B, and C) for the seven soil gas wells (Figure 11). Population A included the soil gas well HBLG with relatively high He and CO2 concentrations, population B included the soil gas wells LHD and HET with relatively low CO2 concentration and medium He concentration, and population C included the soil gas wells TJYC, WFS, DHT, and SWZ with medium CO2 concentration and relatively low He concentration (Figure 11); higher fractional contributions of L-derived CO2 is also present in TJYC, DHT, SWZ, and WFS wells, compared with LHD and HET wells, and with the highest one present in HBLG well (Table 4). Faults are preferred, and high-speed passage ways for uprising of gases from the deep earth [58], through which the escape of gases included in rocks of the deeper crust to the surface soil becomes relative easily [4, 59], and the gases migrating from the deeper earth usually have higher gas concentrations [15, 60]. This indicated the deepest depth the fault cut downward and the most developed fractures in the segment where the HBLG well located, and the shallower depth the fault cut downward and the more developed fractures in the fault segments where the WFS, SWZ, TJYC, and DHT wells located [8, 61, 62].

5.2. Relationship between Earthquakes and Soil Gas Chemical Variations

In the period August 5–19, 2018, a seismic sequence (four earthquakes of magnitude 2.0–4.0) struck the Guye region in the northern section of the Tangshan fault zone. The distances between the epicenter of the main shock (MS 3.7, ) and the soil gas wells ranged between 6 and 45 km (Figure 12).

Abnormal soil gas signals related to earthquakes have been recognized in previous studies [7, 17, 63]. In this study, considerable variations in CO2 concentration were unexpectedly observed in the soil gas of the DHT, HBLG, TJYC, and WFS wells from July 28–September 28, 2018, during which time the seismic sequence occurred (Figure 13). Taking the mean values of CO2 concentration observed in June, October, and December 2018 as background values, the maximum increases in CO2 concentration were calculated to be 175%, 94%, 47%, and 33% at the DHT, HBLG, TJYC, and WFS wells, respectively. These variations of the CO2 concentration were beyond the measurement error of the GXH 3010-E CO2 monitor, which is ±2%.

The time series data of CO2 concentrations in soil gas from the TJYC soil gas well, atmospheric pressure, and local rainfall from June 15, 2018, to December 15, 2018 are shown in Figure 14. It is indicated that the atmospheric pressure and rainfall have little impact on CO2 concentrations in soil gas from the TJYC soil gas well, which may be attributed to the significant depth at which the accumulation chamber had been fixed (5.0 m below the ground).

Therefore, the abnormal variations in soil gas CO2 concentration at wells DHT, HBLG, TJYC, and WFS could be responses to the seismic sequence. The peaks in CO2 concentration could be attributed to increased contributions of gas from underlying limestones as observed in different geological contexts [24, 37].

6. Conclusions

The sources of soil gases in the Tangshan fault zone, North China, have been analyzed based on the soil gas chemical compositions and isotopic ratios obtained by measurement in the field and samples analysis in the laboratory.

The soil gases from the Tangshan fault zone were distinguished by three populations, based on the statistical analysis of CO2 concentration and Rn volume activity of soil gas collected at 388 sampling points using graphical statistical methods. Carbon dioxide of group (A) originates in the superficial soil while carbon dioxide of group (C) comes from deeper geological formations through faults. Carbon dioxide from group (B) could have a complex source, such as biogenic, thermogenic, mechanochemical, or be due to rock weathering.

Further information about the origin of soil gases was obtained according to the geochemical investigations of soil gases sampled in seven soil gas wells, drilled at the sites where group (C) soil gases had been observed. Soil gas compositions and isotopic ratios indicated that air (A) end-member, limestone (L) end-member, and sediment (S) end-member were identified as three major end-member components contributing to soil gas in the Tangshan fault zone; gas from the sediment (S) end-member was the primary source for the soil gases, and limestone (L) end-member was the secondary one.

The relationship between 3He/4He and the average values of CO2 concentration and He concentration of soil gas indicated the deepest depth the fault cut downward and the most developed fractures in fault segment were where the HBLG well was located and the shallower depth the fault cut downward and the more developed fractures in fault segments were where the WFS, SWZ, TJYC, and DHT wells were located.

Variations in CO2 concentration in the DHT, HBLG, and WFS wells observed in concomitance with a local seismic sequence confirm for the first time a possible contribution of carbon dioxide probably generated in underlying limestones.

Data Availability

The data for this paper are available in the text.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors are grateful to Prof. Liwu Li for analyzing the gas compositions and isotopic ratios (Institute of Geology and Geophysics, Chinese Academy of Sciences). This study is jointly supported by the Natural Science Foundation of China (Nos. 41573121, 41402298, and 41503114) and the Basic Science Research Plan of the Institute of Earthquake Science, China Earthquake Administration (Nos. 2016IES0102 and 2019IEF0703).