Geofluids
Volume 2019 (2019), Article ID 1945450, 12 pages
https://doi.org/10.1155/2019/1945450

Research Article

Evidence of Multiple Sources of Soil Gas in the Tangshan Fault Zone, North China

Zhi Chen,¹ Ying Li,¹ Zhaofei Liu,¹ Chang Lu,¹ Yuanxin Zhao,¹ and Jiang Wang²

¹CEA Key Laboratory of Earthquake Prediction (Institute of Earthquake Forecasting), China Earthquake Administration, 63 Fuxing Road, Haidian District, Beijing 100036, China
²Hebei Earthquake Agency, Shijiazhuang 050022, China

Correspondence should be addressed to Ying Li; subduction6@hotmail.com

Received 20 May 2019; Revised 11 July 2019; Accepted 5 August 2019; Published 12 September 2019

Academic Editor: Andrea Brogi

Copyright © 2019 Zhi Chen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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 CO₂ from the gas wells. According to the relationship among the ³He/⁴He, average CO₂ 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 CO₂ 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 [1–3], 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 (CO₂, 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 [5–7]. 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; [8–11]). 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 [12–14]. 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, 22–24]) 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 (F₁), Weishan-Changshan fault (F₂), and Tangshan-Guye fault (F₃). 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].

Figure 1: A map of the study area, with inset showing its location in NE China. F₁: Douhe fault; F₂: Weishan-Tangshan fault; F₃: Tangshan-Guye fault.

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 (Q₄) and upper Pleistocene (Q₃), but Ordovician, Cambrian, and Proterozoic strata consisting mainly of limestone are widely exposed in the northern segment of the Tangshan fault zone, with δ¹³ C (PDB) and δ¹⁸ 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⁻¹ and 20.4 Bq kg⁻¹, 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 CO₂ 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 CO₂ 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 CO₂ monitor (GXH-3010-E, for CO₂ 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 ²¹⁸Po). The detection limit and measurement error of the RAD 7 radon detector were 3.7 Bq m⁻³ 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 CO₂ monitor were 0.01% and ±2%, respectively. An inlet filter was used to protect the detector from dust.

Figure 2: Sketches of the methods used for soil gas measurement at the mobile measurement points.

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 CO₂ 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.

Figure 3: A sketch of the constructed soil gas well.

3.3. Laboratory Analysis

Gas compositions (CO₂, CH₄, Ar, H₂S, N₂, and O₂) and isotopic ratios (δ¹³, ³He/⁴He (R/Ra), CO₂/³He, and ⁴He/²⁰Ne) 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 ³He/⁴He (reported as R/Ra, ) and ⁴He/²⁰Ne values were determined by a different mass spectrometer (VG 5400), while δ¹³ 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 CO₂ 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 CO₂ 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 CO₂ 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⁻³, respectively, with means of 2438.6 ppm and 4730.5 Bq m⁻³, respectively (Table 1).

Table 1: CO₂ concentrations and Rn activities in soil sampling points.

Gas component
No.
Mean
Min
Max
Standard deviation

Rn(Bq m⁻³)
388
4730.5
882.4
17283.5
3465.3

CO₂(ppm)
388
2438.6
514.0
10318.1
2068.0

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 CO₂, CH₄, Ar, H₂S, N₂, and O₂ 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⁻³. The atmospheric concentrations of CO₂ and He in the study area were 0.04% and 5.35 ppm, respectively.

Table 2: Geochemical compositions of gases sampled in shallow depth wells.

No.
Site
CO₂ (%)
Rn (kBq m⁻³)
He (ppm)
N₂ (%)
O₂ (%)
H₂S (ppm)
Ar (%)
CO₂ (%)
CH₄ (ppm)
δ¹³ C (‰)
³He/⁴He (R/Ra)
⁴He/²⁰Ne
CO₂/³He

6.2018
7.2018
8.2018
9.2018
10.2018
11.2018
6.2018
7.2018
8.2018
9.2018
10.2018
11.2018
9.2018

1
DHT
1.91
3.60
5.17
4.54
1.72
2.01
18.40
16.10
19.80
16.16
19.00
25.80
5.04
79.32
18.18
0
0.92
1.57
43
-19.8
0.56
0.31
1.94E+09

2
TJYC
3.76
4.40
3.83
3.30
2.08
2.11
9.69
9.14
8.20
6.17
6.06
6.80
5.98
78.73
19.90
0
0.92
0.46
31
-19.9
0.54
0.30
5.95E+08

3
HBLG
3.59
4.01
4.32
6.23
2.55
3.50
23.20
22.30
22.40
13.80
33.00
27.20
8.11
79.66
17.14
0
0.94
2.25
38
-18.9
0.55
0.31
2.96E+09

4
LHD
0.62
0.59
0.74
0.85
0.89
0.89
6.72
6.47
6.7
5.88
4.59
7.15
7.33
93.65
4.55
56
1.09
0.70
87
-19.8
0.56
0.29
3.68E+08

5
HET
0.14
0.12
0.09
0.13
0.15
0.18
8.53
9.31
5.34
5.93
4.37
7.86
7.12
97.81
0.45
631
1.15
0.32
113
-19.7
0.55
0.28
8.13E+08

6
SWZ
2.04
1.62
2.50
2.82
3.66
3.48
8.71
8.08
9.33
10.20
11.80
12.90
6.17
80.96
15.61
0
0.94
2.49
38
-19.7
0.49
0.30
3.18E+09

7
WFS
1.15
2.45
2.78
3.28
3.26
2.99
31.00
26.10
26.10
25.70
31.80
31.60
6.22
78.64
18.48
0
0.92
1.96
32
-19.4
0.52
0.31
2.48E+09

The ³He/⁴He (R/Ra), CO₂/³He, and ⁴He/²⁰Ne 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 δ¹³ C values of CO₂ were in the range of −19.9 to −18.9‰. The ³He/⁴He (R/Ra), δ¹³, CO₂/³He, and ⁴He/²⁰Ne 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 CO₂ 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 CO₂ 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⁻³ and 13749.8 Bq m⁻³, 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 CO₂ concentration and Rn volume activity and the fractions of the three populations listed in Table 3 were determined using the GSA method [39].

Figure 4: Probability plots of soil gas CO₂ concentrations in the 388 sampling points.

Figure 5: Probability plots of soil gas Rn activities in the 388 sampling points.

Table 3: Statistical parameters of partitioned CO₂ concentrations and Rn volume activities at the mobile measurement points.

Population
Fraction
Mean concentration
95% confidence interval
Mean activity
95% confidence interval

CO₂
Rn
CO₂ (ppm)
Rn (Bq m⁻³)

Population A
0.80
0.83
1513
1491-1643
3894.1
3549.0-4240.2

Population B
0.16
0.12
5415
5396-6109
10229.0
9655.5-10803.0

Population C
0.04
0.05
8893
8679-9328
15471.0
14861.0-16082.0

Contour maps of soil gas CO₂ concentrations and of Rn activities found in the 388 mobile measurement points were mapped in Figures 4 and 5. A consistent spatial coincidence of CO₂ concentrations and of Rn activities was found (Figures 6 and 7), which may indicate common sources of CO₂ and Rn, as CO₂ can act as a carrier gas for Rn [4, 40]. In addition, almost all the CO₂ and Rn values from Population C are concentrated near the fault zone; in contrast, almost all CO₂ and Rn values from Population A are distributed outside the fault zone. Moreover, most of the CO₂ 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 [41–43]. Hence, it could be inferred that the CO₂ and Rn of Population A come from the superficial soil, while those of Population C uprise from deeper locations through faults. Meanwhile, the CO₂ and Rn of Population B could have a complex source due to biogenic, thermogenic, mechanochemical, or rock weathering origins [9, 44–46]. 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.

Figure 6: A contour map of soil gas CO₂ concentrations in the 388 sampling points.

Figure 7: A contour map of soil gas Rn volume activities in the 388 sampling points.

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 CO₂ concentrations (0.12–6.23%) and Rn activities (4.59–33.00 kBq m⁻³) 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 CO₂ concentration and 882.4–17283.5 Bq m⁻³ 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 N₂, O₂, CO₂, 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 ³He/⁴He (R/Ra) and ⁴He/²⁰Ne 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.

Figure 8: Plots showing ³He/⁴He (R/Ra) versus ⁴He/²⁰Ne values for gas samples from the seven soil gas wells.

However, the CO₂ 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 CO₂. As reported in previous studies, underground CO₂ is a crucial link in the interactions among gas, water, rocks, and biological activities [46]. The two main origins of underground CO₂ 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 CO₂ 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 CO₂ from the seven new wells. The δ¹³ C value of −8.3‰ and the CO₂/³He ratio of for the (A) end-member were obtained by measuring air samples (Table 1). For the (L) and (S) end-members, δ¹³ C values of 0‰ and −30‰, respectively, with a corresponding CO₂/³He ratio of , were suggested by Sano and Marty [49]. The fractions contributed by the (A), (L), and (S) end-members to the CO₂ 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 CO₂ at the seven soil gas wells.

Figure 9: Plot showing CO₂/³He versus δ¹³ C values for gas samples from the shallow depth wells. A: air end-member; L: limestone end-member; S: sediment end-member.

The results showed that 62–65% of the CO₂ 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).

Table 4: Fractional CO₂ contributions by different end-members to the soil gases.

Site
Proportions of A-derived CO₂
Proportions of L-derived CO₂
Proportions of S-derived CO₂

(%)

DHT
3
32
65

TJYC
5
31
64

HBLG
2
36
62

LHD
15
23
62

HET
8
28
64

SWZ
2
33
65

WFS
2
34
64

This suggests that the major sources of CO₂ in the Tangshan fault zone are the microbial degradation of organic matter and the decomposition of carbonatic rocks in sediments. Furthermore, the higher N₂/O₂ ratios, together with the strong correlation between the depleted O₂ concentrations (compared to the atmosphere) and the CO₂ concentrations in the wells (correlation coefficient of 0.90; Table 2 and Figure 10), indicate that CO₂ is mainly contributed by the microbial oxidation of organic matter in the sediments of the Tangshan fault zone. In particular, enriched N₂ (93.65% and 97.81%, respectively) and depleted O₂ (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].

Figure 10: Relationship between depleted O₂ concentrations (compared to atmosphere) and CO₂ concentrations from the wells.

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 CO₂ at high concentrations during its decomposition and reaction with groundwater and CO₂, as per Eq. (2) [54–57].

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

In addition, it was shown in the relationship between ³He/⁴He and the average values of CO₂ 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 CO₂ concentrations, population B included the soil gas wells LHD and HET with relatively low CO₂ concentration and medium He concentration, and population C included the soil gas wells TJYC, WFS, DHT, and SWZ with medium CO₂ concentration and relatively low He concentration (Figure 11); higher fractional contributions of L-derived CO₂ 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].

Figure 11: Relationship between ³He/⁴He and the average values of CO₂ concentration and He concentration of soil gas. Red dots: values for wells DHT, HBLG, and WFS. Blue dots: values for wells TJYC, LHD, HET, and SWZ.

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).

Figure 12: Locations of earthquakes in the Tangshan fault zone, August 5–19, 2018.

Abnormal soil gas signals related to earthquakes have been recognized in previous studies [7, 17, 63]. In this study, considerable variations in CO₂ 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 CO₂ concentration observed in June, October, and December 2018 as background values, the maximum increases in CO₂ concentration were calculated to be 175%, 94%, 47%, and 33% at the DHT, HBLG, TJYC, and WFS wells, respectively. These variations of the CO₂ concentration were beyond the measurement error of the GXH 3010-E CO₂ monitor, which is ±2%.

Figure 13: Relationship between earthquakes and soil gas chemical variations.

The time series data of CO₂ 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 CO₂ 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).

Figure 14: CO₂ concentrations in soil gas from the TJYC well, atmospheric pressure, and local rainfall.

Therefore, the abnormal variations in soil gas CO₂ concentration at wells DHT, HBLG, TJYC, and WFS could be responses to the seismic sequence. The peaks in CO₂ 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 CO₂ 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 ³He/⁴He and the average values of CO₂ 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 CO₂ 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).

References

Y. Y. Kagan and D. D. Jackson, “Earthquake forecasting in diverse tectonic zones of the globe,” Pure and Applied Geophysics, vol. 167, no. 6–7, pp. 709–719, 2010.

M. D. Zobak, “Seismology: Earthquake Forecasting and Warning,” Science, vol. 222, no. 4621, p. 319, 1983.

R. J. Zhou, R. C. Tang, Z. Z. Huang, Y. Lin, T. Yang, and W. Zhou, “Tectonic deformation and strong earthquake activities on the east border of Tibet Plateau,” Earth Science Frontiers, vol. z1, pp. 297-298, 2000.

G. Ciotoli, S. Lombardi, and A. Annunziatellis, “Geostatistical analysis of soil gas data in a high seismic intermontane basin: Fucino Plain, Central Italy,” Journal of Geophysical Research, vol. 112, no. B5, article B05407, 2007.

Z. Chen, Y. Li, Z. Liu, J. Wang, X. Zhou, and J. du, “Radon emission from soil gases in the active fault zones in the capital of China and its environmental effects,” Scientific Reports, vol. 8, no. 1, 2018.

T. Doğan, H. Sumino, K. Nagao, K. Notsu, M. K. Tuncer, and C. Çelik, “Adjacent releases of mantle helium and soil CO₂ from active faults: Observations from the Marmara region of the North Anatolian Fault zone, Turkey,” Geochemistry, Geophysics, Geosystems, vol. 10, no. 11, 2009.

C. C. Fu, T. F. Yang, C. H. Chen et al., “Spatial and temporal anomalies of soil gas in northern Taiwan and its tectonic and seismic implications,” Journal of Asian Earth Sciences, vol. 149, pp. 64–77, 2017.

G. Ciotoli, G. Etiope, M. Guerra, and S. Lombardi, “The detection of concealed faults in the Ofanto Basin using the correlation between soil–gas fracture surveys,” Tectonophysics, vol. 301, no. 3-4, pp. 321–332, 1999.

F. Italiano, G. Martinelli, and P. Plescia, “CO₂ degassing over seismic areas: the role of mechanochemical production at the study case of Central Apennines,” Pure and Applied Geophysics, vol. 165, no. 1, pp. 75–94, 2008.

C. Y. King, “Gas geochemistry applied to earthquake prediction: an overview,” Journal of Geophysical Research Solid Earth, vol. 91, no. B12, pp. 12269–12281, 1986.

A. Sciarra, B. Cantucci, and M. Coltorti, “Learning from soil gas change and isotopic signatures during 2012 Emilia seismic sequence,” Scientific Reports, vol. 7, no. 1, 2017.

G. Iovine, I. Guagliardi, C. Bruno et al., “Soil-gas radon anomalies in three study areas of Central-Northern Calabria (Southern Italy),” Natural Hazards, vol. 91, no. 1, pp. 193–219, 2018.

J. L. Lewicki, W. C. Evans, G. E. Hilley, M. L. Sorey, J. D. Rogie, and S. L. Brantley, “Shallow soil CO₂ flow along the San Andreas and Calaveras Faults, California,” Journal of Geophysical Research: Solid Earth, vol. 108, no. B4, p. 2187, 2003.

Y. Li, J. du, X. Wang, X. Zhou, C. Xie, and Y. Cui, “Spatial Variations of Soil Gas Geochemistry in the Tangshan Area of Northern China,” Terrestrial, Atmospheric and Oceanic Sciences, vol. 24, no. 3, pp. 323–332, 2013.

C. C. Fu, T. F. Yang, V. Walia et al., “Variations of soil–gas composition around the active Chihshang Fault in a plate suture zone, eastern Taiwan,” Radiation Measurements, vol. 44, no. 9-10, pp. 940–944, 2008.

J. P. Toutain and J. C. Baubron, “Gas geochemistry and seismotectonics: a review,” Tectonophysics, vol. 304, no. 1-2, pp. 1–27, 1999.

T. F. Yang, C. C. Fu, V. Walia et al., “Seismo-Geochemical Variations in SW Taiwan: Multi-Parameter Automatic Gas Monitoring Results,” Pure and Applied Geophysics, vol. 163, no. 4, pp. 693–709, 2006.

L. L. Chyi, T. J. Quick, T. F. Yang, and C. H. Chen, “The experimental investigation of soil gas radon migration mechanisms and its implication in earthquake forecast,” Geofluids, vol. 10, no. 4, 563 pages, 2010.

C. Y. King, “Radon emanation on San Andreas Fault,” Nature, vol. 271, no. 5645, pp. 516–519, 1978.

S. Džeroski, L. Todorovski, B. Zmazek, J. Vaupotic, and I. Kobal, “Modelling soil radon concentration for earthquake prediction,” Discovery Science, vol. 2843, pp. 87–99, 2003.

K. Kharatian, Y. Travi, and V. Igoumnov, “Activité du gaz radon dans l’air du sol et sismicité locale: exemple du bassin de l’Arax (Arménie),” Comptes Rendus Geoscience, vol. 334, no. 3, pp. 179–185, 2002.

G. Chiodini, S. Caliro, C. Cardellini, R. Avino, D. Granieri, and A. Schmidt, “Carbon isotopic composition of soil CO₂ efflux, a powerful method to discriminate different sources feeding soil CO₂ degassing in volcanic–hydrothermal areas,” Earth and Planetary Science Letters, vol. 274, no. 3-4, pp. 372–379, 2008.

M. C. Ray, D. R. Hilton, J. Muñoz, T. P. Fischer, and A. M. Shaw, “The effects of volatile recycling, degassing and crustal contamination on the helium and carbon geochemistry of hydrothermal fluids from the southern volcanic zone of Chile,” Chemical Geology, vol. 266, no. 1-2, pp. 38–49, 2009.

G. Yuce, C. C. Fu, W. D'Alessandro et al., “Geochemical characteristics of soil radon and carbon dioxide within the Dead Sea Fault and Karasu Fault in the Amik Basin (Hatay), Turkey,” Chemical Geology, vol. 469, pp. 129–146, 2017.

J. Heinicke, U. Koch, and G. Martinelli, “CO₂ and radon measurements in the Vogtland Area (Germany) ‐ A contribution to earthquake prediction research,” Geophysical Research Letters, vol. 22, no. 7, pp. 771–774, 1995.

V. S. Iakovleva and N. K. Ryzhakova, “Spatial and temporal variations of radon concentration in soil air,” Radiation Measurements, vol. 36, no. 1-6, pp. 385–388, 2003.

Q. LIU, J. WANG, J. CHEN, S. LI, and B. GUO, “Seismogenic tectonic environment of 1976 great Tangshan earthquake: results from dense seismic array observations,” Earth Science Frontiers, vol. 14, no. 6, pp. 205–212, 2007.

X. Wang, Y. Li, J. du, and X. Zhou, “Correlations between radon in soil gas and the activity of seismogenic faults in the Tangshan area, North China,” Radiation Measurements, vol. 60, no. 1, pp. 8–14, 2014.

W. Shen, Y. Fang, Q. Ni, Y. Lin, and Y. Liu, “Geochemical variation across the Cambrian-Ordovician boundary at Changshan section of Zhaogezhuang, Tangshan,” Geological Journal of China Universities, vol. 2, no. 3, pp. 324–331, 1996.

P. Feng, Tangshan engineering geological layer set in the key areas, Shijiazhuang University of Economics, Hebei, 2013, Master thesis.

H. Guo, W. Jiang, and X. Xie, “Late-Quaternary strong earthquakes on the seismogenic fault of the 1976 Ms7.8 Tangshan earthquake, Hebei, as revealed by drilling and trenching,” Science China Earth Sciences, vol. 54, no. 11, pp. 1696–1715, 2011.

J. Yang, Soil gas geochemistry characteristics in the capital area of China, Institute of Earthquake Forecasting, China Earthquake Administration, Beijing, 2018, Master thesis.

S. R. Mei, C. H. Hu, C. Z. Zhu, J. Ma, Z. C. Zhang, and M. Y. Yang, The 1976 Tangshan Earthquake, Seismological Press, Beijing, 1982.

N. Segovia, J. L. Seidel, and M. Monnin, “Variations of radon in soils induced by external factors,” Journal of Radioanalytical and Nuclear Chemistry Letters, vol. 119, no. 3, pp. 199–209, 1987.

W. Zhang, J. du, X. Zhou, and F. Wang, “Mantle volatiles in spring gases in the Basin and Range Province on the west of Beijing, China: constraints from helium and carbon isotopes,” Journal of Volcanology and Geothermal Research, vol. 309, no. 1, pp. 45–52, 2016.

A. J. Sinclair, “Selection of threshold values in geochemical data using probability graphs,” Journal of Geochemical Exploration, vol. 3, no. 2, pp. 129–149, 1974.

M. M. Parks, S. Caliro, G. Chiodini et al., “Distinguishing contributions to diffuse CO2 emissions in volcanic areas from magmatic degassing and thermal decarbonation using soil gas 222Rn–δ13C systematics: Application to Santorini volcano, Greece,” Earth and Planetary Science Letters, vol. 377-378, pp. 180–190, 2013.

R. Borgoni, V. Tritto, C. Bigliotto, and D. de Bartolo, “A geostatistical approach to assess the spatial association between indoor radon concentration, geological features and building characteristics: the case of Lombardy, northern Italy,” International Journal of Environmental Research and Public Health, vol. 8, no. 5, pp. 1420–1440, 2011.

G. Chiodini, R. Cioni, M. Guidi, B. Raco, and L. Marini, “Soil CO₂ flux measurements in volcanic and geothermal areas,” Applied Geochemistry, vol. 13, no. 5, pp. 543–552, 1998.

V. Walia, S. J. Lin, C. C. Fu et al., “Soil–gas monitoring: a tool for fault delineation studies along Hsinhua Fault (Tainan), Southern Taiwan,” Applied Geochemistry, vol. 25, no. 4, pp. 602–607, 2010.

G. F. Birchard and W. F. Libby, “Soil radon concentration changes preceding and following four magnitude 4.2–4.7 earthquakes on the San Jacinto Fault in southern California,” Journal of Geophysical Research, vol. 85, no. B6, pp. 3100–3106, 1980.

G. Etiope and S. Lombardi, “Evidence for radon transport by carrier gas through faulted clays in Italy,” Journal of Radioanalytical and Nuclear Chemistry Articles, vol. 193, no. 2, pp. 291–300, 1995.

Y. Yang, Y. Li, Z. Guan et al., “Correlations between the radon concentrations in soil gas and the activity of the Anninghe and the Zemuhe faults in Sichuan, southwestern of China,” Applied Geochemistry, vol. 89, pp. 23–33, 2018.

P. J. Hanson, N. T. Edwards, C. T. Garten, and J. A. Andrews, “Separating root and soil microbial contributions to soil respiration: a review of methods and observations,” Biogeochemistry, vol. 48, no. 1, pp. 115–146, 2000.

R. Ramnarine, C. Wagner-Riddle, K. E. Dunfield, and R. P. Voroney, “Contributions of carbonates to soil CO₂ emissions,” Canadian Journal of Soil Science, vol. 92, no. 4, pp. 599–607, 2012.

C. Song, G. Han, N. Zhuo, and C. Q. Liu, “The characteristics and origin of CO₂ in unsaturated zone at loess tableland of Northwestern China,” Quaternary Sciences, vol. 37, no. 6, pp. 1172–1181, 2017.

W. Gong, X. Yan, and J. Wang, “The effect of chemical fertilizer on soil organic carbon renewal and CO₂ emission–a pot experiment with maize,” Plant and Soil, vol. 353, no. 1-2, pp. 85–94, 2012.

G. Tamir, M. Shenker, H. Heller, P. R. Bloom, P. Fine, and A. Bar-Tal, “Dissolution and Re-crystallization Processes of Active Calcium Carbonate in Soil Developed on Tufa,” Soil Science Society of America Journal, vol. 76, no. 5, p. 1606, 2012.

Y. Sano and B. Marty, “Origin of carbon in fumarolic gas from island arcs,” Chemical Geology, vol. 119, no. 1–4, pp. 265–274, 1995.

M. TAO, “Characteristics of mantle degassing and deep–seated geological structures in different typical fault zones of China,” Science in China Series D, vol. 48, no. 7, pp. 1074–1088, 2005.

R. L. Mulvaney, S. A. Khan, and C. S. Mulvaney, “Nitrogen fertilizers promote denitrification,” Biology and Fertility of Soils, vol. 24, no. 2, pp. 211–220, 1997.

W. A. N. G. Jing, C. H. E. N. G. Yi, C. A. I. Zucong, and Z. H. A. N. G. Jinbo, “Effects of long–term fertilization on key processes of soil nitrogen cycling in agricultural soil: a review,” Acta Pedologica Sinica, vol. 53, no. 2, pp. 292–304, 2016.

X. M. Yu, C. G. Yang, G. M. Dong, and D. Zhan, “Drilling and completion technology in the first geothermal well construction in Jixian system in Tangshan,” Exploration Engineering Media, vol. 41, no. 2, pp. 41–44, 2014.

Z. Chen, Y. Li, C. G. Wang, and L. Yan, “Geochemical characteristics of water from mud volcano in the Wenquan County, Xinjiang,” Earthquake Research in Sichuan, vol. 2, pp. 12–15, 2015.

Z. Chen, Y. Li, Z. Liu et al., “CH₄ and CO₂ Emissions From Mud Volcanoes on the Southern Margin of the Junggar Basin, NW China: Origin, Output, and Relation to Regional Tectonics,” Journal of Geophysical Research: Solid Earth, vol. 124, no. 5, pp. 5030–5044, 2019.

P. Rovira and V. R. Vallejo, “Changes in δ¹³ C composition of soil carbonates driven by organic matter decomposition in a Mediterranean climate: a field incubation experiment,” Geoderma, vol. 144, no. 3-4, pp. 517–534, 2008.

G. Tamir, M. Shenker, H. Heller, P. R. Bloom, P. Fine, and A. Bar-Tal, “Can soil carbonate dissolution lead to overestimation of soil respiration?” Soil Science Society of America Journal, vol. 75, no. 4, p. 1414, 2011.

T. F. Yang, V. Walia, L. L. Chyi et al., “Variations of soil radon and thoron concentrations in a fault zone and prospective earthquakes in SW Taiwan,” Radiation Measurements, vol. 40, no. 2-6, pp. 496–502, 2005.

T. F. Yang, C. Y. Chou, C. H. Chen, L. L. Chyi, and J. H. Jiang, “Exhalation of radon and its carrier gases in SW Taiwan,” Radiation Measurements, vol. 36, no. 1-6, pp. 425–429, 2003.

M. H. Al-Tamimi and K. M. Abumurad, “Radon anomalies along faults in North of Jordan,” Radiation Measurements, vol. 34, no. 1-6, pp. 397–400, 2001.

C. E. Bond, Y. Kremer, G. Johnson et al., “The physical characteristics of a CO₂ seeping fault: the implications of fracture permeability for carbon capture and storage integrity,” International Journal of Greenhouse Gas Control, vol. 61, pp. 49–60, 2017.

Y. Sun, X. Zhou, G. Zheng et al., “Carbon monoxide degassing from seismic fault zones in the Basin and Range Province, west of Beijing, China,” Journal of Asian Earth Sciences, vol. 149, pp. 41–48, 2017.

V. Walia, T. F. Yang, S. J. Lin et al., “Temporal variation of soil gas compositions for earthquake surveillance in Taiwan,” Radiation Measurements, vol. 50, no. 2013, pp. 154–159, 2013.