Stress-Induced Apparent Resistivity Variations at the Kalpin Observatory and the Correlation with the 2020 Mw 6.0 Jiashi Earthquake

: Stress may induce apparent resistivity changes. Clarifying the deformation process of the source media is critical for determining the correlations between resistivity variations and earthquake occurrence. In this study, the stress state of a medium was analyzed by integrating GPS measurements, the spatiotemporal evolution of the load/unload response ratio (LURR), geochemical monitoring, and synchronous apparent resistivity changes preceding the 2020 Mw 6.0 Jiashi earthquake. The medium hosting the Kalpin Observatory underwent elastic deformation before 2019, and the synchronous decreases in the E–W and N–S apparent resistivities from 2015 can be attributed to N–S-dominated compressive stress. The microdamage stage occurred in 2019, with subsequent E–W apparent resistivity variation amplitudes that were ~0.4 Ω m higher than those in previous years. This difference is a result of microdamage to the medium owing to tensile stress during the seismogenic process. The spatiotemporal evolution of the LURR and gas seepage monitoring data also indicate that the medium was damaged prior to the earthquake. Variations in the apparent resistivity measured at the Kalpin Observatory indicate that the medium underwent elastic deformation, followed by microdamage, until stress triggered the earthquake.


Introduction
Resistivity is an important physical property of geotechnical materials that varies with stress [1][2][3][4]. According to data from laboratory experiments, mechanical stress can narrow or widen cracks in rocks and alter the distribution of conductive solutions in cracks, thereby affecting the rock resistivity [5][6][7][8][9]. Earthquakes result from deformation and damage to a medium via continuous stress loading. The seismogenic process is accompanied by changes in the physical parameters of the underground media [10]. A 50% decrease in the resistivity occurred along deep portions of the San Andreas Fault from 1986 [11]. Four years prior to the 1976 Tangshan earthquake, a gradual decrease in the apparent resistivity of the shallow crust was observed at multiple sites [12,13]. A gradual apparent resistivity decrease of up to 6.7% several years prior to the 2008 Wenchuan earthquake was recorded at Pixian station [14]. Researchers have determined that some variations in the apparent resistivity were associated with earthquakes in the Garm region [15]. Temporal changes in resistivity before and after the Izmit (Mw 7.6) earthquake have also been reported [16]. Hautot [17] reported resistivity changes in the Sur-Frêtes ridge, which correlated with mechanical stress caused by changes in the lake level. Based on the fluid distribution underneath the hypocentral region, which is well known as a fluid valve model of Sibson [18], some researchers have claimed that temporal resistivity changes observed using magnetotellurics can reflect fluid movement in the hypocentral region [19][20][21].
In China, shallow apparent resistivity detection began in 1966. There are more than 80 continuous apparent resistivity stations currently in operation. Prior to some large earthquakes, significant resistivity changes have been observed and reported later [12,14,[22][23][24]. Thus, studying changes in shallow apparent resistivity related to earthquakes is still of great significance for capturing seismic signals.
Although many studies have linked resistivity changes to stress [13,25,26], capturing the precursory changes that occur prior to an earthquake remains challenging. This is mainly because the relationship between resistivity changes and stress is complex, with varying resistivity characteristics present in different stress states. For example, although the resistivity of a water-bearing rock decreases significantly under high stress, both decreases and increases in resistivity have been reported at low stress levels [1,2,5,8,[27][28][29][30][31]. The complexity of resistivity variations is likely associated with the difficulty of obtaining the source properties, mechanical state, and deformation state of the medium during the seismogenic process [32]. Therefore, determining the deformation state of a medium is an important step for establishing correlations between changes in apparent resistivity and earthquakes.
Based on the deformation states of rocks during stress loading, Ma et al. [33] partitioned the seismogenic process into the following stages: elastic deformation, yield, meta-instability, and instability (failure). According to Scholz [34], a medium experiences elastic strain buildup, dilatancy-dominated, and water influx-dominated stages. Considering the constitutive law for rocks [35], during the initial stress loading stage (when the rock is in the elastic stage), the loading and unloading response rates (deformation modulus) are approximately equal. However, when the stress accumulation increases, the rock enters the yield stage and its loading response rate exceeds that of unloading. This decreases the strength of the rock or damages it, eventually causing failure. This mechanism has been exploited to propose a method for quantitatively describing the stage of damage to a medium, known as the load/unload response ratio (LURR) [35][36][37]. Abnormal LURRs have been reported for months to years prior to large earthquakes [35,36,38]. Therefore, the evolution of the LURR reflects the deformation state and the evolution of damage within a medium.
On 19 January 2020, the Jiashi Mw 6.0 earthquake occurred in the Atushi-Kalpin Thrust Zone, southern Tien Shan Mountain at the depth of 14.9 km. This area ( Figure 1) hosts the Kalpin Apparent Resistivity Observatory, a baseline geodetic survey observatory, a geochemical monitoring campaign, and several GPS monitoring sites. To determine the correlation between apparent resistivity changes and the Jiashi Mw 6.0 earthquake, the spatiotemporal evolution of the LURR, GPS strain monitoring, and synchronous gas seepage data were utilized in this study to investigate the deformation state of the medium in the seismogenic region. The variations in the apparent resistivity under different stress states were also analyzed.

Tectonic Setting
The Tien Shan orogenic belt undergoes strong active deformation that is characterized by nearly N-S crustal shortening and uplift. Several reverse faults and fold belts are present in the northern and southern Tien Shan Mountains [39]. The latest crustal defor mation in this region has been attributed to remote strain caused by the collision of the Indian and Eurasian tectonic plates. The Atushi-Kalpin Thrust Zone is located in the western South Tien Shan Seismic Belt and the Kuche Depression, which includes the Kashi Depression and the Bachu Uplift [40]. The Kalpin thrust system is divided into two parts by the transform-like strike-slip Piqiang Fault. The western part of the system comprises the Kalpintag (KPT), Aozigertawu (AZT), Toksanatanengbayler (TST), Kekebuke Sanshan (KST), and Aoyibulake (AYT) thrust faults. The eastern part of the system includes the KPT, Tataiertage (TAT), Yimugantawu (YMT), TST, KST, and AYT thrust faults from south to north [39] (Figure 1). According to Yang et al. [41], the crust in the Kalpin Thrust Zone has shortened by 40-45 km since the Cenozoic Era, indicating a shortening rate of 33-37%.
The KPT is the foremost (southern) thrust fault in the Atushi-Kalpin Thrust Zone, comprising an overturned fold-and-front thrust that extends over a lateral distance of 220 km [42][43][44]. Based on an integrated analysis of the surface geology, topography, focal mechanisms, seismic reflection profiles, and InSAR co-seismic deformation, Yao et al. [45] determined that the 19 January 2020 Jiashi Mw 6.0 earthquake occurred on the KPT Thrust Fault.

Kalpin Apparent Resistivity Observatory
The Kalpin Apparent Resistivity Observatory is located in the Atushi-Kalpin Thrust Zone. A ZD8M resistivity monitoring system, with two quasi-Schlumberger arrays (transmitting dipole AB and receiving dipole MN) aligned N-S and E-W, is installed at the observatory. The distances between transmitting dipoles A and B in the N-S and E-W directions are 1004 m and 1007 m, respectively ( Figure 2). This autonomous system records hourly resistivity data. The monitoring field is on the edge of the Gobi Desert, and an abandoned hydropower station, a sinking pond, an irrigation well, and a transformer is located around the station. Considering that these sites are located outside the range of the two quasi-Schlumberger arrays, their interference with the apparent resistivity observations is negligible. Owing to the high-porosity medium present in the observation field, the apparent resistivity is sensitive to stress.
The KPT is the foremost (southern) thrust fault in the Atushi-Kalpin Thrust Zone, comprising an overturned fold-and-front thrust that extends over a lateral distance of 220 km [42][43][44]. Based on an integrated analysis of the surface geology, topography, focal mechanisms, seismic reflection profiles, and InSAR co-seismic deformation, Yao et al. [45] determined that the 19 January 2020 Jiashi Mw 6.0 earthquake occurred on the KPT Thrust Fault.

Kalpin Apparent Resistivity Observatory
The Kalpin Apparent Resistivity Observatory is located in the Atushi-Kalpin Thrust Zone. A ZD8M resistivity monitoring system, with two quasi-Schlumberger arrays (transmitting dipole AB and receiving dipole MN) aligned N-S and E-W, is installed at the observatory. The distances between transmitting dipoles A and B in the N-S and E-W directions are 1004 m and 1007 m, respectively ( Figure 2). This autonomous system records hourly resistivity data. The monitoring field is on the edge of the Gobi Desert, and an abandoned hydropower station, a sinking pond, an irrigation well, and a transformer is located around the station. Considering that these sites are located outside the range of the two quasi-Schlumberger arrays, their interference with the apparent resistivity observations is negligible. Owing to the high-porosity medium present in the observation field, the apparent resistivity is sensitive to stress. The resistivities of rocks in the upper crust mainly depend on the mobility of ions and the distribution of fluids. Conductivity is generated in sediments via current flowing through fluid-filled pores or along the surfaces of clay minerals [25]. Therefore, precipitation and other nontectonic processes affect the resistivity of sediment. Considering the climatic characteristics of mainland China, elevated summer precipitation often decreases the resistivity. Therefore, resistivity data from most stations decreased during the rainy season (July-September), followed by an increase during the dry season (December-March). However, data from some stations (e.g., Pixian Station in Sichuan Province) exhibited an inverse annual pattern, in which the apparent resistivity increased during the rainy season and decreased during the dry season [46]. The apparent resistivity data from the Kalpin Observatory were also characterized by an abnormal annual pattern, i.e., the apparent resistivity was high when precipitation was high and low when the precipitation was low. were extracted using a polynomial fitting method. We fitted the data using the first-order polynomials to obtain the decreasing trend of apparent resistivity. Then the decreasing trend was subtracted from the data. We finally plotted remaining data with the same time scale of one year, as shown in Figure 3b As precipitation typically affects the apparent resistivity measurements, we assessed the precipitation data collected at the Kalpin Meteorological Observatory, which is located less than 5 km from the apparent resistivity observatory, between 1 January 2015 and 30 April 2020 to exclude the potential interference ( Figure 3e). We found that there was a good correlation between precipitation and apparent resistivity, i.e., when the precipitation was high, the resistivity was also high, and when precipitation was low, the apparent resistivity was low. The data yielded an average annual precipitation of 160.32 mm from 2015 to 2019, with annual averages of 148.9 and 171.8 mm for 2018 and 2019, respectively, which are comparable to the amounts recorded in previous years. Atmosphere 2021, 12, x FOR PEER REVIEW 6 of 16 As precipitation typically affects the apparent resistivity measurements, we assessed the precipitation data collected at the Kalpin Meteorological Observatory, which is located less than 5 km from the apparent resistivity observatory, between 1 January 2015 and 30 April 2020 to exclude the potential interference (Figure 3e). We found that there was a good correlation between precipitation and apparent resistivity, i.e., when the

Local Stress Environment
In order to describe the local crustal movement in more detail, we calculated the GPS strain rates based on the distance weighting method [39,40]. The study area was divided into uniform grids and the horizontal velocity gradient of each grid point was calculated one by one.
where, U i denotes the velocity of the GPS sites, t i denotes the fitted velocity of the grid point, x j denotes the distance between the grid point and the corresponding GPS sites, and e ij is the velocity gradient tensor. After adjustment of the calculation, the maximum and minimum principal strain rates were calculated by the following formula.
where ε ij is the strain rate tensor, and ε 1,2 are the maximum and minimum principal strain rates.
We calculated the strain rates for the Tien Shan orogenic belt using GPS data collected from 2015 to 2019. The results indicated that N-S compression occurred in the Tien Shan Mountains. The strain rates to the SW and NE of the Atushi-Kalpin Thrust Zone were high, with maximum principal strain rates that exceeded −5 × 10 −8 /y. In the central region, the maximum principal strain rate decreased to approximately −1 × 10 −8 /y (Figure 4a). The crust in the Tarim Basin, which is located southeast of Kalpin, was characterized by tension. The Kalpin Apparent Resistivity Observatory is located in a transition zone between the compressional and tensional regimes, which explains its relatively low strain rate of -1.89 × 10 −8 /yr. The maximum principal strain at Kalpin also exhibited a N-S trend, with strain rates of -1.80 × 10 −8 and -0.57 × 10 −8 /yr for the N-S and E-W components, respectively (Figure 4b).
Atmosphere 2021, 12, x FOR PEER REVIEW 8 of 16 Wang and Shen [47] processed GPS data collected during the past 25 years (until 2016) for continental China and found that the dilatation rates indicate areal contraction of up to 40 nanostrains/y across the Tien Shan region. This orogenic belt has been characterized by continuous compression [47] and the GPS-derived data used herein also indicated that the Atushi-Kalpin Thrust Zone was under compressive deformation, similar to the period before 2016.
Geodetic baselines in the Atushi-Kalpin Thrust Zone (Figure 1) include the NNE trending Wuqia-Yecheng and Bachu-Wushi, and the E-W trending Bachu-Wuqia. Li et al. [48] analyzed data from these baselines, which indicated that N-S compression has remained the dominant stress component since 2017 in the Atushi-Kalpin Thrust Zone.
Three recent earthquakes of Ms ≥ 5.0, excluding the Jiashi Mw 6.0 earthquake, occurred in the Atushi-Kalpin Thrust Zone (Figure 1). Based on the focal mechanism solutions issued by the U.S. Geological Survey and Harvard University (see Data and Resources), the principal stress direction of the Jiashi Ms 5.0 earthquake on 4 September 2018, was roughly N-S, whereas those of the Atushi Ms 5.1 on 4 November 2018, and the Wushi Ms 5.0 earthquake on 27 October 2019 were NW. GPS strain measurements, geodetic baseline monitoring, and focal mechanism solutions indicate that the Atushi-Kalpin Thrust Zone is dominated by N-S compressive stress, with a subordinate E-W component. In the N-S direction, which has a greater compressive stress, the apparent resistivity decreased more rapidly. Wang and Shen [47] processed GPS data collected during the past 25 years (until 2016) for continental China and found that the dilatation rates indicate areal contraction of up to 40 nanostrains/y across the Tien Shan region. This orogenic belt has been characterized by continuous compression [47] and the GPS-derived data used herein also indicated that the Atushi-Kalpin Thrust Zone was under compressive deformation, similar to the period before 2016.
Geodetic baselines in the Atushi-Kalpin Thrust Zone (Figure 1) include the NNE trending Wuqia-Yecheng and Bachu-Wushi, and the E-W trending Bachu-Wuqia. Li et al. [48] analyzed data from these baselines, which indicated that N-S compression has is dominated by N-S compressive stress, with a subordinate E-W component. In the N-S direction, which has a greater compressive stress, the apparent resistivity decreased more rapidly.

Relationship between Apparent Resistivity and Stress
Zhao et al. [49] investigated the response of apparent resistivity to tensile and compressive stresses under different conditions (e.g., in situ experiments and mining) and found that the apparent resistivity usually decreased under compression (Figure 5a) and increased under tension (Figure 5b). As discussed above, the Atushi-Kalpin Thrust Zone is dominated by N-S compression, with a minor E-W component. Therefore, decreasing apparent resistivity values are expected in the Kalpin Thrust Zone, at a higher rate in the N-S direction relative to the E-W direction. Our observations were consistent with this expectation (Figure 3a,c). According to the focal mechanism solution obtained from the USGS, the Jiashi Mw 6.0 earthquake on 19 January 2020 was a thrust rupture event. The detailed source model of the main shock was strike = 221°, dip = 20°, rake = 72°, and depth = 14.9 km, respectively ( Figure 5c). We projected the location of the Kalpin Apparent Resistivity Observatory onto the focal sphere and found that both observatories were located in the quadrant where the polarity of the first P wave arrival was positive. In other words, the area hosting these stations experienced tension during the process of earthquake preparation (Figure 5c). According to the response of the apparent resistivity to the tensile stress (Figure 5b), the apparent resistivity should increase during the seismogenic period. The distance between the observation station and the epicenter is about 170 km, which is far greater than the earthquake depth, so it is difficult to establish a definite relationship between the earthquake depth and the change in Kalpin Apparent Resistivity. After September 2019, the annual amplitude of apparent resistivity in the E-W direction at Kalpin Observatory increased by ~0.4 Ωm relative to the previous value (Figure 3b).
Based on the stress recorded in the Kalpin Observatory area and the focal mechanism solution, the monotonic decrease in the apparent resistivity recorded in both directions can be attributed to compressive stress. The resistivity increase recorded in the E-W direction preceding the 2020 Mw 6.0 Jiashi Earthquake was assigned to the tensional stress that occurred in the region during the short-term seismogenic stage.

Spatiotemporal Evolution of LURR in the Seismogenic Region
The state of earthquake preparation during a large earthquake (M > 6.0) is determined by changes in the constitutive relations of the medium in the seismogenic region.
As shown in Figure 6, when a nonlinear system is in the elastic phase, nearly the same responses (e.g., the strain, energy release, and so on) to the loading and unloading are observed, and LURR ≈1.0; whereas when the system is stressed beyond the elastic limit, the system enters the inelastic stage, microcracks occur in the media, the responses become significantly different, and LURR > 1.0. Thus, the LURR value might serve as a According to the focal mechanism solution obtained from the USGS, the Jiashi Mw 6.0 earthquake on 19 January 2020 was a thrust rupture event. The detailed source model of the main shock was strike = 221 • , dip = 20 • , rake = 72 • , and depth = 14.9 km, respectively ( Figure 5c). We projected the location of the Kalpin Apparent Resistivity Observatory onto the focal sphere and found that both observatories were located in the quadrant where the polarity of the first P wave arrival was positive. In other words, the area hosting these stations experienced tension during the process of earthquake preparation (Figure 5c). According to the response of the apparent resistivity to the tensile stress (Figure 5b), the apparent resistivity should increase during the seismogenic period. The distance between the observation station and the epicenter is about 170 km, which is far greater than the earthquake depth, so it is difficult to establish a definite relationship between the earthquake depth and the change in Kalpin Apparent Resistivity. After September 2019, the annual amplitude of apparent resistivity in the E-W direction at Kalpin Observatory increased by~0.4 Ωm relative to the previous value (Figure 3b).
Based on the stress recorded in the Kalpin Observatory area and the focal mechanism solution, the monotonic decrease in the apparent resistivity recorded in both directions can be attributed to compressive stress. The resistivity increase recorded in the E-W direction preceding the 2020 Mw 6.0 Jiashi Earthquake was assigned to the tensional stress that occurred in the region during the short-term seismogenic stage.

Spatiotemporal Evolution of LURR in the Seismogenic Region
The state of earthquake preparation during a large earthquake (M > 6.0) is determined by changes in the constitutive relations of the medium in the seismogenic region.
As shown in Figure 6, when a nonlinear system is in the elastic phase, nearly the same responses (e.g., the strain, energy release, and so on) to the loading and unloading are observed, and LURR ≈1.0; whereas when the system is stressed beyond the elastic limit, the system enters the inelastic stage, microcracks occur in the media, the responses become significantly different, and LURR > 1.0. Thus, the LURR value might serve as a useful index to assess the damaged state of the media.
Atmosphere 2021, 12, x FOR PEER REVIEW 10 of 16 Figure 6. Schematic diagram of the constitutive relation of rock. P and R represent, respectively, the load and response. At stage I, the responses (∆R) to the small changes of ∆P in the loading and unloading are almost the same, while at the stage II, the response to the loading is significantly greater than the unloading, and the stage Ⅲ is characterized with macroscopic damage..
The energy (or the Benioff strain) of the small earthquakes was used for LURR calculation through the following formula: where, denotes the Benioff strain released by the i-th event, denotes the energy of the event, and '+' and '−' represent, respectively, the loading and unloading periods derived from the change of Coulomb failure stress (CFS) induced by earth tide in tectonically preferred slip direction [35,36,[49][50][51][52][53].
We calculated LURR by using the small earthquakes of M < 4.0 in the Southern Tianshan Mountains (38° N~44° N, 72° E~84° E) from 1 January 2016 to 1 February 2020. Figure 7e shows the time series of LURR produced with small earthquakes within a circular region of 150 km radius center at the Kalpin Observatory. In 2016, a high LURR occurred near the epicenter of the Jiashi Mw 6.0 earthquake, whereas no similar phenomenon was present in the Kalpin Apparent Resistivity Observatory region (Figure 7a). In the next three years (until the end of 2018), the high LURR adjacent to the epicenter progressively decreased, whereas only a slight increase was observed around Kalpin. However, in 2019, the LURR around the Kalpin Observatory increased rapidly ( Figure  7d). The apparent resistivity values in the E-W direction also increased and the Jiashi Mw 6.0 earthquake occurred during this period. Combining this with Figure 3a, we note that the value of LURR started to increase in early 2019, while the apparent anomaly occurred in September 2019. The difference is possibly caused by sensitivity of different methods. The occurrence rate of small earthquakes during the loading and unloading processes can induce the change of LURR. However, it takes a certain amount of time for the liquid to make these tiny pores well connected and thus significantly change the apparent resistivity. Figure 6. Schematic diagram of the constitutive relation of rock. P and R represent, respectively, the load and response. At stage I, the responses (∆R) to the small changes of ∆P in the loading and unloading are almost the same, while at the stage II, the response to the loading is significantly greater than the unloading, and the stage III is characterized with macroscopic damage.
The energy (or the Benioff strain) of the small earthquakes was used for LURR calculation through the following formula: where, B i denotes the Benioff strain released by the i-th event, E i denotes the energy of the event, and '+' and '−' represent, respectively, the loading and unloading periods derived from the change of Coulomb failure stress (CFS) induced by earth tide in tectonically preferred slip direction [35,36,[49][50][51][52][53]. We calculated LURR by using the small earthquakes of M < 4.0 in the Southern Tianshan Mountains (38 • N~44 • N, 72 • E~84 • E) from 1 January 2016 to 1 February 2020. Figure 7e shows the time series of LURR produced with small earthquakes within a circular region of 150 km radius center at the Kalpin Observatory. In 2016, a high LURR occurred near the epicenter of the Jiashi Mw 6.0 earthquake, whereas no similar phenomenon was present in the Kalpin Apparent Resistivity Observatory region (Figure 7a). In the next three years (until the end of 2018), the high LURR adjacent to the epicenter progressively decreased, whereas only a slight increase was observed around Kalpin. However, in 2019, the LURR around the Kalpin Observatory increased rapidly (Figure 7d). The apparent resistivity values in the E-W direction also increased and the Jiashi Mw 6.0 earthquake occurred during this period. Combining this with Figure 3a, we note that the value of LURR started to increase in early 2019, while the apparent anomaly occurred in September 2019. The difference is possibly caused by sensitivity of different methods. The occurrence rate of small earthquakes during the loading and unloading processes can induce the change of LURR. However, it takes a certain amount of time for the liquid to make these tiny pores well connected and thus significantly change the apparent resistivity. Therefore, prior to 2019 (Figure 7a-c), continuous stress caused damage to the medium around the epicenter. However, in the area around the Kalpin Apparent Resistivity Observatory located 169 km from the epicenter, the crustal deformation was mainly characterized by elastic deformation. Conversely, 2019 was characterized by inelastic deformation and microdamage that caused the increase in the apparent resistivity observed in the Kalpin Observatory region until the Jiashi Mw 6.0 earthquake occurred. Figure 8 shows a survey line of a geochemical monitoring campaign that was conducted 23 km south of the Kalpin Apparent Resistivity Observatory by the Earthquake Administration Bureau of the Xinjiang Uygur Autonomous Region. This survey line, which traversed the KPT in an almost N-S direction, included 12 monitoring sites at intervals of 10-15 m. Radon and hydrogen monitoring were performed quarterly. Figure 8a shows the concentrations of Rn and H2, respectively, at each point since 2016. As expected, higher concentrations were generally measured at sites p5 to p7, which are closer to the KPT. Using p5 as an example, we analyzed the changes in the seepage gas near the KPT. Figure 8c,e display the concentrations of Rn and H2, respectively, measured at site p5 since 2016. The highest H2 concentration recorded before the Jiashi Mw 6.0 earthquake occurred in September 2019, while the corresponding Rn concentration was subordinate Therefore, prior to 2019 (Figure 7a-c), continuous stress caused damage to the medium around the epicenter. However, in the area around the Kalpin Apparent Resistivity Observatory located 169 km from the epicenter, the crustal deformation was mainly characterized by elastic deformation. Conversely, 2019 was characterized by inelastic deformation and microdamage that caused the increase in the apparent resistivity observed in the Kalpin Observatory region until the Jiashi Mw 6.0 earthquake occurred. Figure 8 shows a survey line of a geochemical monitoring campaign that was conducted 23 km south of the Kalpin Apparent Resistivity Observatory by the Earthquake Administration Bureau of the Xinjiang Uygur Autonomous Region. This survey line, which traversed the KPT in an almost N-S direction, included 12 monitoring sites at intervals of 10-15 m. Radon and hydrogen monitoring were performed quarterly. Figure 8a shows the concentrations of Rn and H 2 , respectively, at each point since 2016. As expected, higher concentrations were generally measured at sites p5 to p7, which are closer to the KPT. Using p5 as an example, we analyzed the changes in the seepage gas near the KPT. Figure 8c,e display the concentrations of Rn and H 2 , respectively, measured at site p5 since 2016. The highest H 2 concentration recorded before the Jiashi Mw 6.0 earthquake occurred in September 2019, while the corresponding Rn concentration was subordinate only to that preceding the Jiashi Mw 5.5 earthquake on 4 September 2018. To highlight this anomaly, we employed the seepage gas intensity, which is the ratio of the maximum to the average concentration measured along the survey line. The seepage H 2 intensity in September 2019 was 20% higher than the average value (Figure 8f), whereas the Rn intensity was 12.5% higher than the average value (Figure 8d).

Seepage Gas Monitoring
Atmosphere 2021, 12, x FOR PEER REVIEW 12 of 16 only to that preceding the Jiashi Mw 5.5 earthquake on 4 September 2018. To highlight this anomaly, we employed the seepage gas intensity, which is the ratio of the maximum to the average concentration measured along the survey line. The seepage H2 intensity in September 2019 was 20% higher than the average value (Figure 8f), whereas the Rn intensity was 12.5% higher than the average value ( Figure 8d). Seepage gases such as Rn are useful geochemical indicators [54]. Their concentrations can be elevated mainly through two processes: (1) increased surface area, owing to the formation of new cracks in the rock and (2) pore water flow. Both of these processes can be explained by microscopic damage to rocks [35,55]. Therefore, the enhanced Rn and H2 concentrations measured along the September 2019 survey line in the Kalpin region further support the presence of microdamage in the surrounding medium[56]. Seepage gases such as Rn are useful geochemical indicators [54]. Their concentrations can be elevated mainly through two processes: (1) increased surface area, owing to the formation of new cracks in the rock and (2) pore water flow. Both of these processes can be explained by microscopic damage to rocks [35,55]. Therefore, the enhanced Rn and H 2 concentrations measured along the September 2019 survey line in the Kalpin region further support the presence of microdamage in the surrounding medium [56].
Stress will attenuate drastically as it transfers from the epicenter to the Kalpin Apparent Resistivity Observatory. In theory, if the observatory was located close to the epicenter, the results would be more convincing [57]. However, this region is mainly the sparsely populated Gobi Desert, and additional observatories are not available to verify this hypothesis. In addition, the environmental influence on apparent resistivity would be better described if shallow water level data were available. Despite this, a reasonable explanation can be obtained from existing data for the apparent resistivity variations observed at Kalpin, which indicates that the apparent resistivity of Kalpin is very sensitive to earthquakes that occur on the KPT. Further quantitative research is required to verify the sensitivity of the shallow apparent resistivity to seismic stress, which can help capture earthquake precursor information.

Conclusions
The Kalpin Apparent Resistivity Observatory is located in the transitional region between compressional and tensional zones. This complex stress environment caused the value of the apparent resistivity to synchronously and monotonically decrease in both the E-W and N-S directions since 2015, and the annual amplitude in the E-W direction has increased since 2019.
N-S compressive stress dominates the Atushi-Kalpin Thrust Zone. This accounted for the decline in the apparent resistivities in the N-S and E-W directions in the Kalpin region since 2015. The rate of decrease in the N-S direction was approximately twice that in the E-W direction. Prior to 2019, no significant LURR values were measured in the Kalpin Apparent Resistivity Observatory region, and the medium was characterized by elastic deformation.
In early 2019, the value of LURR started to rise, indicating that the medium entered the inelastic stage and microcracks appeared inside. The apparent resistivity anomaly and gas seepage that occurred in September 2019 may be due to the continuous crack growth or the tiny pores being connected. Considering that the Kalpin region was under tensional stress during the seismogenic period of the Jiashi Mw 6.0 earthquake, the annual amplitude of the apparent resistivity in the E-W direction since September 2019 surpassed the previous average values.
Variations in the apparent resistivity recorded at the Kalpin Observatory before the Jiashi Mw 6.0 earthquake on 19 January 2020 reflect physical processes in the medium throughout the seismogenic region. The medium underwent a transformation from elastic deformation to a microdamaged state, until the stress triggered the earthquake. Further quantitative research is required to confirm the relationship between shallow apparent resistivity and seismic stress.