Spatial–temporal evolution of soil gas Rn before two Ms ≥ 5.0 earthquakes in the mid-eastern of the Qilian fault zone (QLF)

The mid-eastern segment of the Qilianshan fault zone (QLF) on the northeastern margin of the Qinghai–Tibet Plateau is considered one of the key seismic hazard areas. The Zhangye Ms5.0 earthquake and Menyuan Ms6.9 earthquake are the two Ms ≥ 5.0 earthquakes in recent years. The spatio-temporal evolution of Rn across the fault before the two Ms ≥ 5.0 earthquakes were explored by combining a solid seismogenic model and numerical simulation results in this study. The results demonstrates the spatial distribution of Rn concentration intensity varies over time, indicating the evolving characteristics of fracture zone activity. The time-series variation characteristics are closely related the Zhangye Ms5.0 earthquake and Menyuan Ms6.9 earthquake. Overall, in the seismic source area and surrounding medium area of Zhangye Ms5.0 earthquake, the soil gas Rn anomaly across faults characterized by a turning upward trend after continuous decline. The closer to the source area, the more obvious the upward trend. For Menyuan Ms6.9 earthquake, the survey line (HT1) located in the main fracture zone of the earthquake and the survey line (HT7,30km from the epicenter) closer to the epicenter also showed a similar trend, while the other measurement lines in far-field exhibit declining trend before the Menyuan Ms6.9 earthquake. Therefore, the continuous decline trend of soil gas may be crucial information for medium-term earthquake preparation in the seismogenic zone, and the trend of turning upward after continuous decline is a significant signal of short-term seismogenic event in far-field. This research could improve the understanding of the anomalous features of soil gas precursors and tracking the active sections of the fault. According to the model, the earthquake area canseismic source area, the surrounding medium area be divided into three sections: the seismic source area, the surrounding medium area, and the fracture fragmentation area.

Subsurface fluids constitute essential components of the Earth system, playing a critical role in Earth's evolution as an effective medium for interconnecting the planet's various layers [1][2][3][4][5][6] .Soil gas geochemical properties of deep gases discharged in geologically active areas has been extensively employed for precursors of tectonic (e.g., earthquakes and fault activity) [7][8][9][10][11][12][13][14][15] .The preparation and occurrence of earthquakes encompass a complex physical and chemical process, involving the energy transfer of deep materials, changes in medium conditions, and interactions between underground fluids, such as water and gas, and rocks under stress [13][14][15][16][17] .This process necessitates the migration, differentiation, and evolution of fluids [16][17][18] .Active fault zones serve as channels for underground fluid escape and are crucial locations for earthquake preparation and occurrence [16][17][18][19][20][21][22][23][24] .The gas geochemical properties of active fracture zones are often linked to the region's seismic activity [25][26][27] .In fact, these properties have been successfully utilized to elucidate key earthquake mechanisms along major active fault zones in Taiwan, Japan, and China [28][29][30][31][32][33] .Soil gases (CO 2 , CH 4 , Rn, Hg, etc.) in fracture zones frequently exhibit anomalies before and after earthquakes 21,34 , with the magnitude of these anomalies often fluctuating alongside seismic activity 35 .Research on the relationship between released CO 2 fluxes and seismic activity in the L' Aquila region and the Apennines in Italy has demonstrated that CO 2 gas released through rupture zones plays a crucial role in earthquake nucleation, occurrence, and aftershock activity 36 .Continuous monitoring of soil gas Rn concentrations on the Muzaffarabad fault in Pakistan abnormal Rn concentrations approximately 30 days prior to seven earthquakes, with magnitudes ranging from 0.8 to 4.9, impacting the rift 37 .Subsequent to the Wenchuan Ms8.0

Geological setting
The central and eastern sections of the Qilian Mountains are predominantly located along the northeastern edge of the Tibetan Plateau, which is subject to the ongoing extrusion and collisional effects of the Indian and Eurasian plates from a distance 41,42 .This region is characterized by tectonic deformation and intersected by numerous active faults, rendering it one of the most seismically active and intense areas in mainland China, with substantial crustal movements [43][44][45][46] .The primary active fracture zones in this study area include the Sunan-Qilian fault, the Yumushan fault, the Minle-Damaying fault, the Huangcheng-Shuangta fault, and the Lenglongling fault, all of which are crucial components of the Qilian Mountains' northern margin (Fig. 1).These faults exhibit a NWW orientation and consist of multiple compression-torsional fractures arranged parallel to one another and in oblique rows 47 .The fault zone demonstrates a high degree of thrust characteristics, with the thrust fault cutting through strata from various periods, leading to the overthrusting of pre-Cenozoic strata, metamorphic rocks, and Paleozoic granites on Neogene mudstone and early-middle Pleistocene unconsolidated gravel.The zone also displays significant recent tectonic activity 48 .The Lenglongling fault, part of the North Qilian active fault zone, is situated at the northeastern boundary of the Tibetan Plateau uplift zone.Geotectonically, the fault is positioned within the North Qilian Fold Belt, north of the corridor transition zone and south of the Middle Qilian Uplift Zone 49 .The fault's eastern end connects to the Gulang fault and the Maomaoshan fault, while the western end links to the northern edge of the Tolaishan fault 50 .Approximately 120 km in length, the Lenglongling fault exhibits distinct linear characteristics [51][52][53] and generally strikes at 110° N ~ 115° E. The fault's late Quaternary features include left-lateral strike-slip movement with a minor amount of dip-slip motion [54][55][56][57] .Although the Lenglongling fault's Quaternary slip rate remains contentious, it is still considered to have the highest slip rate within the Qilian-Haiyuan fault zone and plays a significant role in plateau deformation 58 .

Survey line layout and measuring instruments
Soil gas radon (Rn) measurements were conducted at nine survey sites (i.e., CGS, YLX, GJZ, SMC, QSZ, BDK, HC, MJZ, and GJT sites, numbered HT1 to HT9) situated on various NWW-oriented faults running parallel to each other and in oblique rows, covering the mid-eastern region of the Qilian fault (Table 1).The CGS (HT1) is located on the Sunan-Qilian fault, YLX (HT2) and GJT (HT3) are located on the Yumushan fault, SMC (HT4), QSZ (HT5), and BDK (HT6) are situated on the Minle-Damaying fault, HC (HT7) and MJZ (HT8) are on the Huangcheng-Shuangta fault, and GJT (HT9) is located on the Lenglongling fault.The survey line is oriented perpendicular to the fault strike, with 10-m intervals and 10-12 survey points arranged on each line to fully encompass the fault (Fig. 2).
To minimize the impact of meteorology on soil gas concentrations, surveys were conducted under stable meteorological conditions in June of each year from 2016 to 2021.Rn concentrations, temperature, humidity, and atmospheric pressure were simultaneously measured using an AlphaGUARD detector (model PQ2000, Saphymo GmbH, Germany).According to the detector's requirements, 15 to 20 values were measured during each assessment.All samples were collected using a stainless-steel probe inserted into the ground at a depth of 80-100 cm, depending on soil consistency and thickness, in order to minimize the effects of meteorological variables [59][60][61] .No extreme climatic variations were observed during the monitoring period, and the meteorological conditions at each measuring point remained relatively stable.

Data analysis method
In this paper,the concentration strength of each survey site was calculated.Concentration intensity refers to the measurement of the change degree of the of fault gas concentration near the fault plane relative to the background concentration.It can intuitively reflect several characteristics of fault media such as porosity, fragmentation, and fracture development, as well as the mechanics and motion state of fault activity It is the main characteristic parameter for analyzing the spatial characteristics of fault gas.The calculation formula is as follows: Concentration strength (S) = Anomaly threshold/background value.The background value is obtained by calculating the mean value of the remaining data after removing jump values, reflecting the normal accumulation level of specific elements or components in a given fault zone.The anomaly threshold can be expressed as the background value plus 'n' standard deviations (n = 1, 2, 3,…); in this  Calculations and evaluations led to the decision to use the maximum value method for determining the concentration strengthof soil gas Rn(S_max).This approach utilizes the maximum abnormal concentration values to calculate concentration intensity, better explaining the peak value of the concentration intensity for each profile.Additionally, this method minimizes differences in site-effect factors that influence gas release concentrations, such as soil and rock properties 62,63 .Details of the nine survey lines are listed in Table 1.

Soil temperature, humidity and atmospheric pressure
The average, maximum, and minimum values of soil temperature, humidity, and atmospheric pressure for each survey line are presented in Table 2. From 2016 to 2021, the soil temperature variation range for all measurement lines spanned from 11.7 °C to 35 °C.Soil temperatures were lower in the eastern section of the fracture zone  www.nature.com/scientificreports/measurement lines compared to the western section; however, the soil temperature variation remained more stable across six periods at each site.During this period, soil moisture variation from HT1 to HT9 ranged between 3.6 and 58.8.The humidity for all survey lines in 2018 was low, which is associated with that year's precipitation.The humidity differences in other years were minimal.As observed in Fig. 1s, the atmospheric pressure at each measurement line remained highly stable.Except for HT2, where the atmospheric pressure was concentrated around 820 mbar, the atmospheric pressure in the remaining measurement lines was approximately 800 mbar.This indicates that the measurement environment was relatively stable throughout the monitoring period.

Spatial distribution of Rn concentration intensity
The spatial distribution of Rn concentration intensity is depicted in Fig. 3.In 2016, the highest concentration intensity was observed at HT5 (5.42) in the middle section, followed by HT1(4.22) in the western section, and then HT9(3.73) and HT6 (3.28).The concentration intensities for HT2, HT3, HT4, and HT8 were 2.47, 2.2, 2.25, and 3.03, respectively.From 2017 to 2020, the spatial distribution characteristics of concentration intensity were stronger in the east than in the west, with the highest values observed at HT9 and HT8 in the eastern section.The HT9 values from 2017 to 2019 were 3.25, 3.18, and 3.31, respectively, while the HT8 value in 2020 was 2.99.Additionally, the western part of the QLF (Sunan section), including the Yumushan fault zone and the central part of the Minle Damaying fault zone, exhibited non-segmented Rn concentration intensity distributions between 1.74 and 2.63.In 2021, the spatial distribution was less distinct, with values across the entire fault zone uniformly ranging between 1.48 and 2.17.In conclusion, the spatial distribution of Rn concentration intensity varies over time, indicating the evolving characteristics of fracture zone activity.The degree of fault gas release is influenced by subsurface medium conditions, local stress states, and both current and previous seismic activity [64][65][66] .The fault's permeability is crucial for controlling the intensity of soil gas concentration changes.Consequently, soil gas concentration strength is relatively higher in more open sections.Our previous research findings on the Haiyuan fault zone, the northern margin of the West Qinling fault, and the Liupanshan fault zone all support this perspective.

Coupling relationship between fault soil Rn and fault activity in QLF
The spatial distribution of seismic activity is illustrated in Fig. 4. We divided the entire research area into two sections, east (A) and west (B).According to the time series statistics of seismicity frequency in the eastern and western sections, the seismicity in this region is very intensive, and there is a trend change of small amplitude enhancement in the time series.It can be seen that the seismicity of QLF has obvious segmentation that it is   www.nature.com/scientificreports/more frequent seismicity in the east section than that in the west section.The Huangcheng-Shuangta fault, the Lenglongling fault, and their intersection exhibit the strongest seismic activities.Notable events such as the Ms6.9 earthquake in 2022, the Menyuan Ms6.4 in 1986, and the Ms6. 4 earthquake in 2016 all occurred in this area.Additionally, earthquakes with Ms ≥ 5.0 and smaller seismic activities are more concentrated in this region.Moreover, using remote sensing image interpretation, field geological survey, fault displacement measurement, and geomorphic surface age determination, the recoil sliding rate of typical dislocation points of the Yumushan fault was found to be (0.55 ± 0.15) mm/a, and the left-sliding rate was (0.95 ± 0.11) mm/a 25 .The linear fitting results of the vertical sliding rate of the Minle-Damaying fault were (0.91 ± 0.09) mm/a 41 .The Huangcheng-Shuangta fault experienced significant activity in the Holocene, most notably the Gulang Ms8.0 earthquake in 1927 43 .The Quaternary slip rate of the Lenglongling fault was determined to be (4.3 ± 0.7) mm/a by radiocarbon dating methods, while the slip rate since the late Holocene was (4.3 ± 0.36) mm/a 49 .The Huangcheng-Shuangta and Lenglongling fracture zones exhibit more activity than other fracture zones, indicating a relatively open fault.This corresponds well with the spatial distribution pattern of soil gas Rn in the QLF.
Trend evolutionary characteristics of Rn related to Ms ≥ 5.0 earthquake in QLF Since the initiation of cross-fault soil gas observations in the eastern QLF, the Menyuan Ms6.9 earthquake on January 8, 2022, and the Zhangye Ms5.0 earthquake on September 16, 2019 have occurred.The seismogenic fault for the Menyuan Ms6.9 earthquake is the western segment of the Lenglongling fault, which ruptured to the Tolashan fault on the west side of its northwest end [67][68][69] , while the Zhangye Ms5.0 seismogenic fault is the northern margin fault of the Qilian Mountain (Sunan Qilian section).The Rn trend changes before and after 2019 have captured our particular interest, as they may be related to the Zhangye Ms5.0 earthquake(Table 4).The Rn concentration intensity of all lines exhibited an apparent "V" pattern trend, with a significant decrease followed by an upward turn before 2019 (Fig. 5).However, differences in the range of change were observed.Specifically, HT1, located in the Sunan-Qilian fault and 80 km from the epicenter of the Zhangye Ms5.0 earthquake, displayed a significant decrease in Rn concentration intensity from 2016 to 2018, followed by an ascent before the earthquake.HT2 and HT3, situated in the Yumushan fault and 2-27 km from the epicenter of the Zhangye Ms5.0 earthquake, showed different patterns.The Rn concentration intensity of HT2 decreased significantly from 2016 to 2017 and increased significantly from 2017 to 2019, while the Rn concentration intensity of HT3 declined from 2016 to 2018 and increased significantly from 2018 to 2019.HT4-HT6, located in the Minle-Damaying fault with a distance of approximately 16-67 km from the epicenter, exhibited varying trends.Although HT4 was closest to the epicenter, the trend of continuous decline was not apparent, but an upward turn was observed before the Zhangye Ms5.0 earthquake.The trends of HT5 and HT6 showed a significant continuous decline followed by a slight upward turn before the earthquake.HT7 and HT8 are located in the Huangcheng-Shuangta fault.The Rn concentration intensity of HT8 increased after a continuous decline.However, the trend of HT9, which is 250 km away, was evidently weaker than that of other survey lines before the Zhangye Ms5.0 earthquake.In addition, we observed that the Rn concentration intensity exhibited an ascending trend in 2020, following a declining background in HT1, HT4, and HT7 (Fig. 5 and Table 4).Specifically, these lines showed a decline from 2019 to 2020 and then an ascent from 2020 to 2021, forming a noticeable "V" pattern, with an upward turn following the decline.HT1 is situated in the primary rupture zone of the northern edge of the Qilian Mountains, where several Ms ≥ 3.0 events occurred after the Menyuan Ms6.9 earthquake.HT7 is located within the Huangcheng-Shuangta fault, approximately 30 km from the Menyuan Ms6.9 earthquake's epicenter.HT4 is positioned in the Minle-Damaying fault to the west of the Menyuan Ms6.9 earthquake, with an epicenter roughly 100 km away.Furthermore, with the exception of these three lines, all other measurement lines continued to exhibit a declining trend before the Menyuan Ms6.9 earthquake.We assumed that these trend changes may be related to the Menyuan Ms6.9 earthquake that occurred on January 8, 2022.
In summary, we found that Rn concentrations near the epicenter areas displayed an upward trend following a continuous decline before both Ms ≥ 5.0 earthquakes in QLF (Table 4 and Fig. 5).The persistent decline trend may be a crucial signal for moderate to strong earthquakes in the medium to long term.We have observed such trend characteristics before several moderately strong earthquakes on the northeastern margin of the Qinghai Tibet Plateau.

The mechanism of Rn spatial-temporal evolution before Ms ≥ 5.0 earthquakes
Earthquake precursors are phenomena that manifest at specific stages of earthquake preparation.It is essential to connect the spatial-temporal evolution of precursors with the earthquake process and mechanism.Since the 1970s, the Fracture Collusion model (IPE) and Dilatancy Hypothesis model (DD) have been proposed.However, these models mainly address the issue of epicenter precursors, and they cannot explain the complex spatiotemporal evolution process of precursors 70 .Through extensive research on large earthquakes (Ms ≥ 6.0), it has been observed that strong earthquakes typically do not occur in the high seismicity zones of fault systems, but rather in seismic gaps 71,72 .Subsequently, the concept of seismic gaps was introduced, and a solid seismogenic model was established to elucidate the conditions, processes, and associated precursors of strong earthquakes, as well as their spatial-temporal evolution 73,74 .This model is supported by a wealth of evidence from deep structures, mechanical analyses, rock fracture experiments, numerical simulation results, and observed facts 75 .According to the model, the earthquake area can be divided into three sections: the seismic source area, the surrounding medium area, and the fracture fragmentation area (Fig. 6a).The source area may comprise a block with a relatively intact macroscopic rupture or a locked segment within a fault zone.The earthquake preparation process consists of several stages, including the elastic accumulation stage, the nonlinear deformation and local hardening stage, the local hardening stage and expansion, the large-scale expansion stage, and the rupture stage 76,77 .
The seismogenic medium is a saturated two-phase medium composed of solid lithosphere and pore fluid.Under the influence of regional tectonic stress, the medium undergoes various stages of deformation, causing changes in the physical parameters of the rock 78,79 .Among these parameters, porosity is a significant quantity closely related to radon emission.As the rock's porosity increases, so does its surface area, radon ejection coefficient, and radon content in the rock's pore fluid.According to the solid seismogenic model, during the elastic accumulation stage and nonlinear deformation and local hardening stage in the seismic source area, porosity decreases due to compression 73,75,76 .Consequently, the soil gas radon anomaly exhibits a decreasing trend at the beginning and middle period of earthquake inception.However, during the local hardening stage, expansion, and large-scale expansion stage, porosity increases.Therefore, in the short-to medium-term of earthquake preparation, especially during the impending stage, the soil gas radon anomaly primarily trends upward.Furthermore, underground fluid anomalies first appear in the source area and surrounding medium area due to the earliest expansion occurring in the source region.During the medium to short-term seismogenic stage, the anomalous area expands, and numerous precursors emerge in the fracture fragmentation area, while the development of anomalies in the source area remains relatively stable.This phenomenon was confirmed by numerical simulation of the spatiotemporal evolution of underground fluid anomalies under the strong solid seismogenic model (Fig. 6b).Overall, in the seismic source area and surrounding medium area, the soil gas Rn anomaly across faults may exhibit a "V" shaped change trend, characterized by a turning upward trend after continuous decline 75 .The closer to the seismic source area, the earlier the turning upward time.Therefore, the continuous decline trend of soil gas may be crucial information for medium to long-term earthquake preparation, and the trend of turning upward after continuous decline is a significant signal of a medium and short-term seismogenic event.These important signals require a comprehensive understanding of the tectonic background and activity of fault zones.This model effectively explains the relationship between two Ms ≥ 5.0 earthquakes in the QLF.For the Zhangye Ms5.0 earthquake, the Rn concentration intensity of HT2 (27km away from the epicenter) decreased significantly from 2016 to 2017 and increased significantly from 2017 to 2019.In HT3 (2 km away from the epicenter), the Rn concentration intensity decreased from 2016 to 2018 and increased significantly from 2018 to 2019.Other survey lines also demonstrated some trend anomalies, but the short-term anomalies (turning upward trend after continuous decline) were not more pronounced than in HT2 and HT3.For the Menyuan Ms6.9 earthquake, HT1, HT4, and HT7 also displayed turning upward trends after continuous decline before the earthquake.HT1 is located on the the northern edge of the Qilian Mountains, the southern section of Sunan, 200 km from the epicenter.Despite this distance, it is situated along the main fault at the northern margin of the Qilian fault.The Lenglongling fault zone, where the Menyuan Ms6.9 earthquake occurred, is also an integral component of the northern Qilian fold belt in terms of geotectonics.Several aftershocks with a magnitude of Ms ≥ 3.0 were recorded in the region following the Menyuan Ms6.9 earthquake.Furthermore, except for three lines, all other measurement lines exhibited a continuing decline trend before the Menyuan Ms6.9 earthquake.The Menyuan Ms6.9 seismic activity and tectonic deformation directly expose the adjustment and transmission of tectonic deformation throughout the Qilian Mountain fault zone.As the survey lines did not directly cross the seismic fault of the Menyuan Ms6.9 earthquake, the abnormal features were not as pronounced as those observed in the Zhangye Ms5.0 earthquake.The response of soil gas Rn to seismic events is not only related to the magnitude and distance from the epicenter but also controlled by the tectonic stress of deep major faults.Thus, the present observation results and the analysis of seismic cases support the body seismogenic model and numerical simulation computation results.This model elucidates the relationship between tectonic geochemical characteristics and rupture blockage, offering theoretical support for determining the time, space, and intensity of seismogenesis.The robust body seismogenic model can clearly explain the relationship between tectonic geochemical features and fracture locking, providing theoretical support for estimating the time, space, and intensity of earthquake nucleation.

Conclusion
In this study, we explored the spatial-temporal evolution characteristics of soil gas Rn from 2016 to 2021 in the mid-eastern of the Qilian fault zone (QLF) before Zhangye Ms5.0 earthquake and Menyuan Ms6.9 earthquake by combining a solid seismogenic model and numerical simulation results.By analyzing the spatial-temporal evolution characteristics of soil gas Rn before two Ms ≥ 5.0 earthquakes in the Qilian fault zone, the following conclusions can be drawn: 1. From a spatial distribution perspective, in 2016, the highest concentration intensity was observed in the middle (HT5) and west section (HT1).From 2017 to 2020, the spatial distribution characteristics of concentration intensity were stronger in the east than in the west.In 2021, the spatial distribution was less distinct, with values across the entire fault zone uniformly.The spatial distribution of Rn concentration intensity varies over time, indicating the evolving characteristics of fracture zone activity.2. From a temporal evolution trend before earthquake,the soil gas Rn concentration intensity across mideastern of the Qilian fault zone(QLF) exhibited a significant response to the Zhangye Ms5.0 and Menyuan Ms6.9 earthquakes.For the Zhangye Ms5.0 earthquake, the Rn concentration intensity in major survey lines (except HT8 and HT9) showed a noticeable "V" temporal evolution trend, with an continue decline and then turn upward.For the Menyuan Ms6.9 earthquake, HT1 (located on the the northern edge of the Qilian Mountains) and HT4, and HT7() also displayed turning upward trends after continuous decline before the earthquake.While the other measurement lines exhibit a continued declining trend before the Menyuan Ms6.9 earthquake.3.These observational facts were supported by the sturdy body seismogenic model and numerical simulation results.Consequently, we hypothesized that the continuous decline trend of fault gas may serve as a reliable indicator of the interlocked section of fault tectonic activity, and the trend of turning upward after continuous decline may represent a significant signal of a medium and short-term seismogenic event in the source area.Furthermore, the sturdy body seismogenic model could potentially provide theoretical support for determining the time, space, and intensity of seismogenesis based on the relationship between the tectonic geochemical characteristics and rupture locking.
In summary, seismic prediction remains a difficult global scientific problem, and fluid geochemistry is one of the potential tools for earthquake prediction.The extraction of earthquake precursor information from soil gas must be based on a specific physical mechanism of seismogenesis.The results of long-term research have shown that, by studying the physical mechanisms of forecasting, the sensitivity and convenience of fault gas can be used to trace locked segments directly on dangerous faults for the identification and determination of precursor information, which can compensate for the limitations of springs and well-exposure locations used by subsurface fluid.

Figure 1 .
Figure 1.Tectonic settings of the middle-eastern section of the Qilian Mountains fault zone and the distribution of the soil gas measurement lines (Geology data set is provided by Geospatial Data Cloud site, Computer Network Information Center, Chinese Academy of Sciences (http:// www.gsclo ud.cn); Active faults data are provided by sharing Infrastructure of National Earthquake Data Center (http:// data.earth quake.cn)).

Figure 2 .
Figure 2. Schematic diagram of field sampling methods.(a) Laying of measurement lines; (b) Measurement Method.

Figure 5 .
Figure 5. Radon concentration intensity time-series trend related to two Ms≥5.0 earthquakes in the mideastern of the Qilian fault zone (QLF) (red arrows indicate the strend related to Zhangye Ms5.0 earthquake (red dashed line) and blue arrows related to Menyuan Ms6.9 earthquake).

Figure 6 .
Figure 6.Solid seismogenic model and schematic diagram of numerical model ((a) is schematic of numerical simulation calculation model, (b) is the variation of porosity with time at various representative points in different regions) t represents different stages of earthquake preparation 75 .

Table 1 .
List of soil gas Rn survey line in the mid-eastern section of QLF.

No Longitude (°E) Latitude (°N) Number of measurement points Length of survey line (m) Direction of line measuring (°) Fault zone
The purpose of this is to eliminate the measured value instability caused by environment.

Table 3
presents the soil gas Rn measurement results for 9 profile from 2016 to 2021 in the eastern section of the QLF.In 2016, Rn concentrations were relatively high at HT7, HT8, and HT9, with values of 36 Bq/L, 34 Bq/L, and 39 Bq/L, respectively.Concentrations were lower at HT2, HT3, and HT4, with values of 15 Bq/L, 9 Bq/L, and 13 Bq/L, respectively.Rn concentrations in the other lines ranged between 21 Bq/L and 23 Bq/L.In 2017, Rn concentrations at HT7, HT8, and HT9 were also higher than those at other locations, with values of 40 Bq/L, 50 Bq/L, and 60 Bq/L, respectively.The Rn concentration at HT1 was the lowest,with values of 18 Bq/L.In 2018, the easternmost sites HT8 and HT9 had the highest Rn concentrations, with values of 55 Bq/L and 72 Bq/L, respectively, while the westernmost sites HT1-HT3 had the lowest , with values of 17 Bq/L, 19 Bq/L, and 16 Bq/L, respectively.In 2019, the easternmost site HT9 had the highest Rn concentration of 72 Bq/L, the westernmost site HT1 had the lowest concentration at 12 Bq/L, and HT3 also exhibited a relatively low concentration of 15 Bq/L.The concentration range for other measurement lines was 24-37 Bq/L.In 2020, the Rn concentrations at the easternmost sites HT8 and HT9 were relatively high, at 41 Bq/L and 50 Bq/L, while the concentrations for other measurement lines ranged from 20 Bq/L to 27 Bq/L.In 2021, Rn concentrations at HT8 and HT9 reached 37 Bq/L and 48 Bq/L, respectively, surpassing those at other measurement lines.Rn concentration varied over time but also demonstrated relative stability, indicating the reliability of the measurement data.

Table 3 .
Results of the soil gas Rn concentrations (D_max) and strength (S_max) of the mid-eastern section of QLF from 2016 to 2021.

Table 4 .
Trend characteristics of Rn concentration intensity across the mid-eastern section of QLF from 2016 to 2021.