Search

An Ms 6.0 earthquake struck the Luxian County in southern Sichuan basin, China, on 15 September 2021, which was the largest local event thus far and caused heavy casualties and serious property losses. Using data of 107 temporary and six permanent stations from five months before to four and a half months after the time of the mainshock in the Rongchang-Luxian region, we utilize the machine learning-based workflow for phase picking, earthquake detection and relocation, and obtain precise hypocenter data for 22,806 events. We found that earthquakes increase significantly since July 2021, showing north-northeast-, northeast-, east-northeast- and east-southeast-directed seismic strips. Almost all events are located within the sedimentary layers at 1–7 km depth above the crystalline basement, and more than half of them occurred in Silurian and Ordovician formations. The Ms 6.0 Luxian earthquake sequence occurred in an area where no clear seismicity was observed before the mainshock. On map view, the aftershock zone shows a convex polygon with an overall strike of 112°. The focal depths of aftershocks southwest to the mainshock are distributed in a narrower range of 4–7 km and extend to a wider range of 2–7 km on the northeast side. It is speculated that the seismogenic fault is a blind thrust fault dipping to the southwest at 2–7 km depth. The focal mechanisms of 13 ML 3.0+ earthquakes show either pure reverse faulting or reverse faulting dominated mechanisms with a strike-slip component. The centroid depths fall in the range of 1.95–6.75 km. The maximum and middle principal stress axes are nearly horizontal; the azimuth of the maximum horizontal stress is about 120° and the stress shape ratio is 0.78. These results provide an important reference for reunderstanding the local seismic risk and for detailed studies on the seismological aspects about the unusual Luxian earthquake in the future.

The spatiotemporal heterogeneity of b values has great potential for helping in understanding the seismogenic process and assessing seismic hazard. However, there is still much controversy about whether it exists or not, and an important reason is that the choice of subjective parameters has eroded the foundations of much research. To overcome this problem, we used a recently developed non-parametric method based on a data-driven concept to calculate b values. The major steps of this method include (1) performing a large number of Voronoi tessellations and Bayesian information criterion (BIC) value calculation, selecting the optimal models for the study area, and (2) using the ensemble median (Q2) and median absolute deviation (MAD) value to represent the final b value and its uncertainty. We investigated spatiotemporal variations in b values before and after the 2019 Changning MS=6.0 earthquake in the Sichuan Basin, China. The results reveal a spatial volume with low pre-mainshock b values near the mainshock source region, and its size corresponds roughly with the rupture area of the mainshock. The anomalously high pre-mainshock b values distributed in the NW direction of the epicenter were interpreted to be related to fluid invasion. The decreases in b values during the aftershock sequence along with the occurrences of several strong aftershocks imply that b values could be an indicator of the stress state. In addition, we found that although the distribution characteristics of b values obtained from different methods of investigation are qualitatively consistent, they differ significantly in terms of their specific values, suggesting that the best way to study the heterogeneous pattern of b values is in the joint dimension of space-time rather than separately in time and space. Overall, our study emphasizes the importance of b-value studies in assessing earthquake hazards.

High-frequency (>1 Hz) ambient noise is usually closely related to anthropogenic activities. During the outbreak and spread of the COVID-19, as various anthropogenic activities are restricted, high-frequency ambient noise level has been observed to be reduced on a worldwide scale. The continuous waveform data at dense broadband seismic stations from ChinArray in eastern North China provides a good opportunity to study the temporal and spatial patterns of the ambient noise level in the region, and to further study the influencing factors, such as the topography and the population density. In this study, we calculated the average power spectral density of ambient noise at each station ±90 days around the Spring Festival in 2019 and in 2020, analyzed the noise level at different stations through normal times, Spring Festivals, epidemic control period, and recovery period, and studied the influencing factors of the noise level. We found that normally high-frequency (1–10 Hz) ambient noise correlates well with the surrounding sedimentary thickness: The noise level is higher when the surrounding sedimentary layer is thicker and vice versa. It correlates moderately with local population density and is time-varying due to anthropogenic activities. During the Spring Festival in 2019 and in 2020, and the epidemic control period after the Spring Festival in 2020, the reduction extent of the noise level correlates moderately with both the sedimentary thickness and population density; the ambient noise level reduces more significantly to the south of 40° N than to the north of it in the study region. Considering that the sedimentary thickness beneath each station is not time-varying, the variation in ambient noise level due to anthropogenic activities is clearly amplified by the sedimentary layer.

Anomalous changes in groundwater have been long reported prior to some sizable earthquakes. Here we examined the water temperature and water level data at the Lijiang well in Yunnan, China from November 1, 2007 to September 30, 2017, and obtained the correlation coefficients between them varying with time by using different time windows. We then use Molchan’s error diagram to obtain the area skill scores for different alarm sizes and distances, which are used for evaluating the prediction efficiency of local sizable earthquakes. We found that (i) when the time window is 90 days and the correlation coefficient is greater than 0.10, a sizable earthquake (M ≥ 5.0) is more likely to happen within 200 km during the following 45 days; (ii) when the time window is 90 days and the correlation coefficient is greater than 0.65, a sizable earthquake (M ≥ 5.0) is more likely to happen within 150 km during the following 45 days. Because there are two separate aquifers and the water temperatures become warmer with depths larger than 200 m, we speculate that normally the well water is recharged mainly by colder rainfall through the shallow aquifer, and increasing water level causes decreasing water temperature and vice versa. However, before a local sizeable earthquake (M ≥ 5.0) within 200 km, there may be deep hotter water flowing into and out of the well through the deep aquifer and breaking the negative correlation between them, and positive correlation coefficients are even observed. This study may provide some useful clues for short to medium-term forecast of sizable earthquakes in the Lijiang region and can provide methodological reference for similar studies in other regions.

The spatiotemporal heterogeneity of b values has great potential for understanding the seismogenic process and assessing the seismic hazard. However, there is still much controversy about whether it exists or not, and an important reason is that the choice of subjective parameters has eroded the foundations of many researches. To overcome this problem, we used a recent developed non-parametric method based on the data-driven concept to calculate b values. The major steps of this method include: 1) perform a large number of Voronoi tessellation, Bayesian information criterion (BIC) value calculation and selection of the optimal models for the study area, and 2) use the ensemble median (Q2) and median absolute deviation (MAD) value to represent the final b value and its uncertainty. We investigated spatiotemporal variations of b values before and after the 2019 Changning MS 6.0 earthquake in Sichuan Basin, China. The results reveal a spatial volume with low pre-mainshock b values near the mainshock source region, and its size corresponds roughly with the rupture area of the mainshock. The anomalously high pre-mainshock b values distributed in the NE direction of the epicenter was interpreted to be related with fluid invasion or increased pore pressure. The decreases of b values during the aftershock sequence along with the occurrences of several strong aftershocks imply that b values could be an indicator of stress state. In addition, we found that although the distribution characteristics of b values obtained from different way of investigating are qualitatively consistent, they differ significantly in terms of their specific values, suggesting that the best way to study the heterogeneous pattern of b values is in the joint dimension of space-time rather than alone in time and space. Overall, our study emphasizes the importance of b value studies on assessing the earthquake hazards.

We examined water level data at the Rongchang (RC) well in the Rongchang natural gas field in China's Sichuan Basin from November 2007 through 2019, where injection-induced seismicity with a magnitude of up to 5.2 has been reported. Coseismic water level changes caused by earthquakes from 3 km to 3500 km away, as well as changes in water temperature, provide a good opportunity to systematically study the mechanisms underlying the response of water level to earthquakes. Among 41 ML ≥ 3.0 local earthquakes, 20 events exceeded the magnitude-distance threshold for the water level to change. Eleven events caused water levels to rise, and only one earthquake led to a decrease in water level. For the 22 distant earthquakes with a seismic energy density greater than the threshold of 0.0017 J/m3, ten events caused rises in the water level and only the MW 7.9 Wenchuan Earthquake caused a temporary, but large, decrease in water level. In total, eight local events and 11 distant events, which were expected to cause coseismic water level changes, actually resulted in no changes. The findings imply that an earlier sizable earthquake may decrease the sensitivity of the well for a certain period of time (approximately two months). The dominant mechanisms of coseismic water level changes in response to near-, mid-, and far-field earthquakes are different. For near-field earthquakes, the drop in coseismic water level is governed by the combined effects of static strain changes and the release of trapped gas bubbles due to the strong vibration of seismic waves in the aquifer; however, only static strain changes are dominant after mid-field earthquakes. For far-field earthquakes, the dominant mechanism is aquifer softening due to growth of gas bubbles in the groundwater that is associated with weak and low frequency seismic shaking of passing surface waves.

Crustal thickness and Poisson's ratio are important parameters to characterize regional isostasy state and material composition or state. Using the teleseismic waveform data from 141 permanent stations and 785 temporary stations in southwest China, we obtain the crustal thickness and average Poisson's ratio by the H-κ stacking of receiver functions. In the west (the SE Tibetan plateau and the Yunnan-Burma-Thailand block) and southeast (the Cathaysian block and southern Yangtze craton) of the study region, there are high correlation coefficients for the crustal thicknesses between what we obtain from the receiver functions and what we calculate from the Airy isostasy model, indicating that a state of isostasy can be achieved at the crust-mantle boundary beneath these two regions. In the northeast (northern Yangtze craton), the correlation coefficient is lower, indicating that the effect of the lithosphere needs to be considered for a regional isostasy. Intermediate Poisson's ratios (0.26≤σ≤0.28) are found beneath the northern Panzhihua-Emeishan region. Combing the high velocity features from previous study, we speculate that it may be related to the Emeishan large igneous province. High Poisson's ratios (σ>0.28)- are found beneath the SE Tibetan plateau and the nearby strike-slip faults, such as the Anninghe-Zemuhe fault and the northern Xiaojiang fault. Combing the low velocity zones from previous study, we speculate that there may be partially melted and lower crustal flow.

The sedimentary and crustal thicknesses and Poisson's ratios of the NE Tibetan Plateau and its adjacent regions are estimated by the h – κ stacking and CCP image of receiver functions from the data of 1,317 stations. The horizontal resolution of the obtained results is as high as 0.5° × 0.5°, which can be used for further high resolution model construction in the region. The crustal thicknesses from Airy's equilibrium are smaller than our results in the Sichuan Basin, Qilian tectonic belt, northern Alxa block and Qaidam Basin, which is consistent with the high densities in the mantle lithosphere and may indicate that the high-density lithosphere drags crust down overall. High Poisson's ratios and low velocity zones are found in the mid- and lower crust beneath eastern Qilian tectonic belt and the boundary areas of the Ordos block, indicating that partial melting may exist in these regions. Low Poisson's ratios and low-velocity anomalies are observed in the crust in the NE Tibetan Plateau, implying that the mafic lower crust is thinning or missing and that the mid- and lower crust does not exhibit melting or partial melting in the NE Tibetan Plateau, and weak flow layers are not likely to exist in this region.

We examined the water level data at the LGH well in Sichuan, China, from December 2007 to July 2015 and their responses to multiple large earthquakes with seismic energy densities greater than 10− 4 J/m3. Co-seismic water level declines were observed in response to eleven earthquakes out of twelve in the farfield, and co-seismic water level increase was observed in one nearfield case. The water level declines in the farfield showed a linear relation with the common logarithm of the seismic energy densities, whereas the water level increase in the nearfield fell away from this relation, indicating that the farfield responses and the nearfield response were produced by distinct mechanisms. We used the phase shift of tidal responses as a proxy for permeability and found that permeability enhancements were observed both in the farfield and nearfield. The co-seismic water level declines in response to the distant earthquakes could be explained by permeability enhancements caused by the passage of seismic waves through the mobilization of colloidal particles; the co-seismic water level increase in response to the nearfield case could be caused both by the compression of the static stress and by the seismic waves.

Using the surface wave records of 504 teleseismic events at 50 temporary and 92 permanent seismic stations in southwest China, we obtain phase velocity maps at 10, 15, 25, 40, 60, and 75 s at a grid spacing of 0.5° × 0.5° from the interstation correlation method and surface wave tomography. We also obtain the S wave velocity structures beneath three profiles using the joint inversion of receiver functions and surface waves. At short periods (10 and 15 s), high-velocity zones (HVZs) are found in the Panzhihua-Emeishan region, the Sichuan basin and the Weixi-Lijiang region surrounding the low-velocity zones (LVZs) from Songpan-Ganzi block to the east of Lijiang, where the elevations are significantly higher. At long periods (40, 60, and 75 s), HVZs are found in the Weixi-Lijiang region, the Panzhihua-Chuxiong basin,and the Kunming-Tonghai region, which form a belt in the center part of the study area. The fast polarization directions on both sides of the belt defined by the shear wave splitting of teleseismic SKS waves vary significantly and indicate that the flow of material from the plateau is blocked in two different depth intervals and leads to different horizontal extents. The long-period maps and the structures along the three profiles show that LVZs are present in the upper mantle beneath rapidly slipping fault zones, such as the Anninghe-Zemuhe-Xiaojiang fault zone, the Red River fault zone, and the Xiaojinhe fault zone, implying that these faults penetrate deep into the mantle.

[1] The transfer functions of the well-aquifer systems response to atmospheric loading and Earth tide can be used to calculate the well-aquifer properties. Due to the low signal-to-noise ratio, the study on barometric response at frequencies higher than 8 cycles per day (cpd for short) is almost blank. Using the recorded water level and barometric pressure as well as the corresponding theoretical tidal volumetric strain at 17 well stations in China, we obtained continuous barometric and tidal responses of the well-aquifer systems by cross-spectra estimation. It shows that the barometric responses are stable at low frequencies (0.1–0.5 cpd), while fluctuate at tidal frequencies (0.5–8 cpd). In the high-frequency band (8–30 cpd), by stacking the transfer functions to suppress the noise, we obtained stable barometric responses for the first time. According to the low- and high-frequency barometric responses, we can better judge whether the aquifers are confined in the timescale that we focused and whether the wellbore storage effect can be ignored, and expect to determine the fluid flow properties of the aquifers more reliably. For the three aquifers whose water table drainage and wellbore storage effects are ignored, we used the barometric and tidal responses to estimate their formation properties, which are consistent with previous results. The tidal strain sensitivities are related to the aquifer lithology, which are mainly controlled by the compressibility of the porous matrixes with different porosities and different aspect ratio fractures.

Taking the M2 wave as calibration signals, we extract the phase shifts of the water level relative to the Earth tide in the Zhouzhi well by utilizing the cross-correlation function. And we further obtain the apparent permeability variation in the aquifer of the Zhouzhi well in 2008. Comparison with the commonly used tidal analysis software Baytap-G shows that phase shifts obtained by cross-correlation function are more stable. The resulting apparent permeability of the Zhouzhi well aquifer fluctuates with time, indicating it is a dynamically controlled parameter. The 2008 Wenchuan earthquake caused the apparent permeability increasing drastically, which is interpreted as the combination effects of effective stress changes and the barriers removing in the flow channel due to seismic wave pressure pulse. After the Wenchuan earthquake, the effective stress began to recover and the impurities deposited gradually, causing the apparent permeability to decrease a month later and almost recover to the pre-earthquake level in six months.

Based on isotropic linear poroelastic theory and under the undrained condition, we summarize three equations connecting the Skempton’s coefficient B with the groundwater level. After analysis, we propose a method to calculate the Skempton’s coefficient B according to the relationship between water level and tidal strain. With this method we can get the value of B without the earthquake occurrence, which can provide the high frequency waves for research. Besides, we can also get the in-suit Skempton’s coefficient B without the experiment of rock physics. In addition, we analyze the observed data of Changping station recorded in groundwater monitoring network (abv., GMN) before and after the Wenchuan MS8.0 with this method, and find out there’s a slight change of the value of B after the seismic waves passed by, which implies that the propagation of seismic waves may have brought some variations to the poroelastic medium of the well.