Intrinsic and scattering attenuations of the Sichuan-Yunnan region in China from S coda waves

Seismic attenuation is a fundamental property of the Earth's media.Attenuation structure for the complicated geological structures with strong seismicity in the Sichuan-Yunnan region is poorly studied.In this study,we collected 108,399 waveforms of 11,517 local small earthquakes with magnitudes between 1.5 and 3.5 from January 2014 to September 2021 in the Sichuan-Yunnan region and its adjacent areas.We employed an envelope inversion technique for separating the intrinsic and scattering attenuations of the S coda wave,and obtained the intrinsic and scattering attenuation structures for frequencies between 0.25 and 8.00 Hz.The attenuation structures correlate well with the geological units,and some major faults mark the attenuation variations where historic large earthquakes have occurred.The regional average attenuation shows a negative frequency dependence.The average scattering attenuation has a faster descending rate than the average intrinsic attenuation,and is dominant at low frequencies,while at high frequencies the average intrinsic attenuation is stronger.The lateral variation in the intrinsic attenuation is consistent with the variation in heat flow,the scattering attenuation may be related to the scatter distribution and size.The total attenuation is consistent with the previous studies in this region,and the separate intrinsic and scattering attenuation may be useful in understanding regional tectonics and important in earthquake prevention and disaster reduction.

Magnitude and distance distribution of strong aftershocks in Sichuan-Yunnan region

研究了川滇地区1966年以来12次主震震级Ms≥6.5的地震序列中，主震与强余震（本文定义为所有Ms≥5的余震）震级差分布特征和强余震与主震距离的分布特征。结果表明，强余震与主震震级差服从截断的指数分布，据此推导出了强余震与主震震级差的概率密度函数；强余震距离分布的优势范围是距主震10-39km，且强余震与主震震中的距离服从正态分布。

Determination of the degree 120 time-variable gravity field in the Sichuan-Yunnan region using Slepian functions and terrestrial measurements

The terrestrial time-variable gravity measurements are characterized by a high signal-to-noise ratio and sensitivity to the sources of mass change in the Earth's crust.These gravity data have many applications,such as surface deformation,groundwater storage changes,and mass migration before and after earthquakes.Based on repeated terrestrial gravity measurements at 198 gravity stations in the Sichuan-Yunnan region(SYR)from 2015 to 2017,we determine a time series of degree 120 gravity fields using the localized spherical harmonic(Slepian)basis functions.Our results show that adopting the first 6 Slepian basis functions is sufficient for effective localized Slepian modeling in the SYR.The differences between two gravity campaigns at the same time of year show an obvious correlation with tectonic features.The degree 120 timevariable gravity models presented in this paper will benefit the study of the regional mass migration inside the crust of the SYR and supplement the existing geophysical models for the China Seismic Experimental Site.

Influence of 2008 Wenchuan earthquake on earthquake occurrence trend of active faults in Sichuan-Yunnan region

The Wenchuan earthquake coseismic deformation field is inferred from the coseismic dislocation data based on a 3-D geometric model of the active faults in Sichuan-Yunnan region. Then the potential dislocation displacement is inverted from the deformation field in the 3-D geometric model. While the faults＇ slip velocities are inverted from GPS and leveling data, which can be used as the long-term slip vector. After the potential dislocation displacements are projected to long-term slip direction, we have got the influence of Wenchuan earthquake on active faults in Sichuan-Yunnan region. The results show that the northwestern segment of Longmenshan fault, the southern segments of Xianshuihe fault, Anninghe fault, Zemuhe fault, northern and southern segments of Daliangshan fault, Mabian fault got earthquake risks advanced of 305, 19, 12, 9.1 and 18, 51 years respectively in the eastern part of Sichuan and Yunnan. The Lijiang-Xiaojinhe fault, Nujiang fault, Longling-Lancang fault, Nantinghe fault and Zhongdian fault also got earthquake risks advanced in the western part of Sichuan-Yunnan region. Whereas the northwestern segment of Xianshuihe fault and Xiaojiang fault got earthquake risks reduced after the Wenchuan earthquake.

Wavelet analysis and interpretation of gravity data in Sichuan-Yunnan region, China

The Bouguer gravity anomaly data of Sichuan-Yunnan region and its vicinity were analyzed with wavelet trans- formation method. In the process, complete orthogonal wavelet function system with good symmetry and higher vanishing moment was selected to decompose the gravity anomaly into two parts. With the power spectral analysis on the decomposed anomalies, we interpreted that the two parts of anomalies represent the density variation in upper and middle crust, and in deep crust and uppermost mantle, respectively. The two parts of anomalies indicate the difference between shallow and deep tectonics. The results of shallow-layer apparent density mapping reveal that: a) the crustal density in Sichuan basin is higher than that in Songpan-Garze orogenic zone; b) the density of Kangdian rhombic block is heterogeneous; c) the boundary faults of Kangdian block are of different density fea- tures, suggesting different tectonic signification. The results of deep-layer apparent density mapping show a similar, but not the same, density distribution pattern as the shallow results, and indicate that the tectonics of shallow and deep crust are different, they may be in a status of incomplete coupling. Our results also show that the earthquakes in this area are controlled not only by the fracture zones but also by the deep density distribution.

Determinations of directions of the mean stress field in Sichuan-Yunnan region from a number of focal mechanism solutions

Based on the spatial orientation and slip direction of the fault plane solutions, we present the expression of corresponding mechanical axis tensor in geographic coordinate system, and then put forward a method for calculating average mechanical axis tensor and its eigenvalues, which involves solving the corresponding eigenequation. The method for deducing mean stress field from T, B, and P axes parameters of a number of focal mechanism solutions has been verified by inverting data of mean stress fields in Fuyun region and in Tangshan region with fitting method of slip direction, and both results are consistent. To study regional average stress field, we need to choose a population of focal mechanism solutions of earthquakes in the massifs where there are significant tectonic structures. According to the focal mechanism solutions of 256 moderate-strong earthquakes occurred in 13 seismic zones of Sichuan-Yunnan region, the quantitative analysis results of stress tensor in each seismic zone have been given. The algorithm of such method is simple and convenient, which makes the method for analyzing tectonic stress field with large amount of focal mechanism solution data become quantified.

Numerical simulation of the influence of lower-crustal flow on the deformation of the Sichuan-Yunnan Region

On the basis of distribution of active fault and regional rheological structure, a three-dimensional finite element model of Sichuan-Yunnan region, China, is constructed to simulate contemporary crustal motion and stress distribution and discuss the dynamic mechanism of crustal motion and deformation in the Sichuan-Yunnan region. Linear Maxwell visco-elastic model is applied, which includes the active fault zones, the elastic upper crust and viscous lower crust and upper mantle. Four different models with different boundary conditions and deep structure are calculated. Some conclusions are drawn through comparison. Firstly, the crustal rotation about the eastern syntaxis of the Himalaya in the Sicuan-Yunnan region may be controlled by the special dynamic boundary condition. The drag force of the lower-crust on the upper crust is not negligible. At the same time, the main active fault zones play an important role in the contemporary crustal motion and deformation in Sichuan-Yunnan region.

Study on Postseismic Impact of Wenchuan Earthquake on the Sichuan-Yunnan Region Based on Three-Dimensional Viscoelastic Finite Element Method

Based on the lateral segmentation and vertical stratification characteristics of the crustal medium in Sichuan-Yunnan region,and the asymmetry of the static dislocations on the coseismic fault plane of the Wenchuan M_S8.0 earthquake,we built a three-dimensional viscoelastic finite element model of the crust in the Sichuan-Yunnan region.The postseismic impact of the Wenchuan M_S8.0 earthquake on the Sichuan-Yunnan region was studied.The results show that:①The far-field horizontal deformation caused by the viscoelastic relaxation of the medium in the 10 years after the earthquake is about 0-20 mm within the Sichuan-Yunnan diamond-shaped block,which has a greater influence on north side and smaller on south side.②In the 10 years after the earthquake,the far-field vertical deformation caused by the viscoelastic relaxation effect of the medium is small,and it shows an increase of about 0-4 mm in most areas of the Sichuan-Yunnan diamond-shaped block.③The Xianshuihe fault and the eastern segment of the East Kunlun fault,which are close to the seismogenic fault,show a high gradient on deformation fields after the earthquake.④In order to compare with the strong earthquake activity in the SichuanYunnan region after the Wenchuan earthquake,the horizontal stress state and the Coulomb failure stress change of the active block boundary are also calculated.From the spatial distribution of the coseismic and postseismic displacement field,the fault activity characteristics reflected by the stress state and the stress loading of the fault layer reflected by the Coulomb failure stress change,there is a certain correlation with the spatial distribution of strong earthquake activity in this region.

Characteristics of isostatic gravity anomaly in Sichuan-Yunnan region,China

Sichuan-Yunnan region in China, a tectonic transition belt where earthquakes occurred frequently and intensely, has a distinct variation characteristic of gradient zone of Bouguer gravity anomaly （BGA）. Many deep faults and epicenters of severe earthquake scatter along the BGA gradient zones. Here we apply two forward models （Airy model and Vening Meinesz model） of isostatic gravity mechanisms （local versus regional） in this region to calculated the isostatic gravity anomaly （IGA）. Afterwards, the relationship between IGA and distribution of faults as well as seismicity is also illustrated. The IGA results show that the two models are similar and most parts of the study area are in an isostatic state. Most featured faults are distributed along the steep anomaly gradient zones; earthquakes tend to occur in the non-isostatic area and steep gradient belt of IGA. The distribution of root thickness based on regional mechanism can be associated with the main trend of BGA variation. The regional mechanism is more plausible and closer to the reality because of its relatively further consideration of the horizontal forces derived from adjacent particles in the crust. Then we analyze the effect of isostasy on the tectonic movements and find that the isostatic adjustment is not the main cause of the continuous uplift process of Longmenshan Mountain fault zone. which is due to the Indian-Eurasian continental collision.

Analysis of Ground Motion Attenuation Characterization for Moderate Earthquakes in the Sichuan-Yunnan Region

The ground motion attenuation models for PGA,PGV and response spectrum at rock sites and soils sites are derived separately from the digital strong motion records of moderate earthquakes in the Sichuan-Yunnan region after 2007. By comparison with Yu Yanxiangs attenuation model for rock sites in western China,reliability of the model is verified for moderate earthquake. According to the distribution of strong motion data against magnitude and epicentral distance,the applicability and reliability of the results in this paper are discussed.

Strong Earthquake Activity and Its Relation to Regional Neotectonic Movement in Sichuan-Yunnan Region

Based on analyzing space inhomogeneous image of strong earthquake activity, the image of source rupture and the mechanical property of the source fault in Sichuan-Yunnan region, the relations among the strong earthquake activity, active fault, modern movement status of active blocks and structural background of the deep media have been discussed, and the characteristics of strong earthquake activity and possible mechanism have been also discussed.

Study on Integrated Recurrence Behaviors of Strong Earthquakes Along Entire Active Fault Zones in the Sichuan-Yunnan Region, China

Based on historical earthquake data, we use statistical methods to study integrated recurrence behaviors of strong earthquakes along 7 selected active fault zones in the Sichuan-Yunnan region. The results show that recurrences of strong earthquakes in the 7 fault zones display near-random, random and clustering behaviors. The recurrence processes are never quasi-periodic, and are neither strength-time nor time-strength dependent. The more independent segments for strong earthquake rupturing a fault zone has, the more complicated the corresponding recurrence process is. And relatively active periods and quiescent periods for earthquake activity occur alternatively. Within the active periods, the distribution of recurrence time intervals between earthquakes has relatively large discretion, and can be modelled well by a Weibull distribution. The time distribution of the quiescent periods has relatively small discretion, and can be approximately described by some distributions as the normal. Both the durations of the active periods and the numbers of strong earthquakes within the active periods vary obviously cycle by cycle, leading to the relatively active periods having never repeated quasi-periodically. Therefore, the probabilistic assessment for middle- and long-term seismic hazard for entireties of active fault zones based on data of historical strong earthquakes on the fault zones still faces difficulty.

Recent tectonic stress field zoning in Sichuan-Yunnan region and its dynamic interest

In this paper, we have carefully determined the stress zones in the Sichuan-Yunnan region with reference to the in-situ stress data of hydraulic fracturing and the inverted fault slip data by using the step-by-step convergence method for stress zoning based on focal mechanism solutions. The results indicate that the tectonic stress field in the Sichuan-Yunnan region is divided into 3 stress zones by 2 approximately parallel NNW-trending stress transition belts. The area between the 2 belts is the Sichuan-Yunnan stress zone where the maximum principal stress σ1 is just in the NNW direction. The eastern boundary of Sichuan-Yunnan stress zone （the eastern stress transition belt） is basically consistent with the eastern boundary of Sichuan-Yunnan rhombic block. The western boundary of Sichuan-Yunnan stress zone （the western stress transition belt） is not totally consistent with the western boundary of Sichuan-Yunnan rhombic block. The northern segment of the western stress transition belt extends basically along the Jinshajiang fault and accords with the western boundary of Sichuan-Yunnan rhombic block, while its southern segment does not extend along the southwestern boundary of the rhombic block, i.e., Honghe fault and converge with the eastern stress transition belt, but stretches continuously in the NNW direction and accords with the Yingpanshan fault. We therefore consider that under the combined influence from the northward motion of India Plate, the southeastward shift of east Qinghai-Xizang Plateau and the strong obstruction of South China block, the tectonic stress field in the Sichuan-Yunnan region might not be totally controlled by the previous tectonic frame and new stress transition belt may have possibly formed.

Tidal Response Characteristics of the Well Water Level in the Sichuan-Yunnan Region

In this paper,the long time series data of the well water-level data of 12 wells in the Sichuan and Yunnan area is analyzed by the Baytap-G tidal analysis software,and well water level tidal response characteristic parameters( amplitude ratio and phase change)are extracted. We analyzed the features of the shape and stage change,and characteristic parameters of the tidal response of well water level before and after the earthquakes,which can provide a new method and approach to analyzing the response relationships between well water level and earth tide and barometric pressure. The results show that Luguhu Well and 9 other wells are affected by earth tides,and their well water level amplitude ratios and phases are relatively stable; the Nanxi Well and Dayao Well water level changes are affected by the barometric pressure combined with tide force,and their well water level amplitude ratios and phases are more discrete. The water level amplitude ratios and phases of Jiangyou Well,Luguhu Well and Dongchuan Well are significant to large earthquakes,and the relationship between seismic energy density and water level amplitude ratios and phases of M_2 wave of the three wells are presented.

The Two-dimensional Time Coordinate System and Time Prediction Research of M≥6.7 Strong Earthquakes in the Sichuan-Yunnan Region

Since the 20 thcentury,the time intervals of M ≥6.7 strong earthquakes in the SichuanYunnan region show obvious regularity.Using the years of the strong events,a twodimensional time coordinate system is generated,based on which,the time prediction model is constructed for strong earthquakes in the Sichuan-Yunnan region.Prediction analysis shows that there is risk of generating four earthquakes with M ≥ 6.7 in the Sichuan-Yunnan region in the future 16 years,and there are strong signals for M ≥6.7earthquakes for periods 2012-2021 and 2025-2029.The strong earthquakes may occur around 2014-2015,2019 and 2027.

The Evolutionary Characteristics of Seismic Activity in the Sichuan-Yunnan Region,China

Based on the earthquake activity characteristics of the diamond block in the Sichuan- Yunnan region and by using the method of the meso-scope damage dynamics and damage evolution, we studied the damage evolution process for moderately strong earthquakes along two seismic belts. The original combination patterns of all the units which illuminate the changes from stable state to destroyed state are given. All these patterns can direct the earthquake prediction practice in this region.

Study on the Relationship Between Recent Crustal Deformation and the Temporal-Spatial Distribution of Strong Earthquakes (M_s≥6.0) in the Sichuan-Yunnan Region

Analysis of deformation data measured across the faults, regional vertical deformation data and GPS measurements in the Sichuan-Yunnan region made since the 1980s permitted us to conclude that the crustal deformation in the region during this period of time was relatively weak and caused the occurrence of earthquakes (M S≥6.0), which were not distributed along the major boundary active faults in the region after the 1981 Dawu M S 6.9 earthquake and that the seismic activity is characterized by quasi-clockwise migration. Thus, it follows that earthquake prediction research should be focused on the central part of the Sichuan-Yunnan region in the coming years. Finally, a concept of temporal division of the region into active blocks is suggested and the preliminary result of the division is given in the paper.

Construction and Application of a Regional Kilometer-Scale Carbon Source and Sink Assimilation Inversion System(CCMVS-R)

CO_(2)is one of the most important greenhouse gases(GHGs)in the earth’s atmosphere.Since the industrial era,anthropogenic activities have emitted excessive quantities of GHGs into the atmosphere,resulting in climate warming since the 1950s and leading to an increased frequency of extreme weather and climate events.In 2020,China committed to striving for carbon neutrality by 2060.This commitment and China’s consequent actions will result in significant changes in global and regional anthropogenic carbon emissions and therefore require timely,comprehensive,and objective monitoring and verification support(MVS)systems.The MVS approach relies on the top-down assimilation and inversion of atmospheric CO_(2)concentrations,as recommended by the Intergovernmental Panel on Climate Change(IPCC)Inventory Guidelines in 2019.However,the regional high-resolution assimilation and inversion method is still in its initial stage of development.Here,we have constructed an inverse system for carbon sources and sinks at the kilometer level by coupling proper orthogonal decomposition(POD)with four-dimensional variational(4DVar)data assimilation based on the weather research and forecasting-greenhouse gas(WRF-GHG)model.Our China Carbon Monito ring and Verification Support at the Regional level(CCMVS-R)system can continuously assimilate information on atmospheric CO_(2)and other related information and realize the inversion of regional and local anthropogenic carbon emissions and natural terrestrial ecosystem carbon exchange.Atmospheric CO_(2)data were collected from six ground-based monito ring sites in Shanxi Province,China to verify the inversion effect of regio nal anthropogenic carbon emissions by setting ideal and real experiments using a two-layer nesting method(at 27 and 9 km).The uncertainty of the simulated atmospheric CO_(2)decreased significantly,with a root-mean-square error of CO_(2)concentration values between the ideal value and the simulated after assimilation was close to 0.The total anthropogenic carbon emissions in Shanxi Province in 2019 from the assimilated inversions were approximately 28.6%(17%-38%)higher than the mean of five emission inventories using the bottomup method,showing that the top-down CCMVS-R system can obtain more comprehensive information on anthropogenic carbon emissions.

Investigation of Groundwater Quality with Borehole Depth in the Basin Granitoids of the Ashanti Region of Ghana

The dependence of groundwater quality on borehole depth is usually debatable in groundwater studies, especially in complex geological formations where aquifer characteristics vary spatially with depth. This study therefore seeks to investigate the relationship between borehole depth and groundwater quality across the granitoid aquifers within the Birimian Supergroup in the Ashanti Region. Physicochemical analysis records of groundwater quality data were collected from 23 boreholes of public and private institutions in the Ashanti Region of Ghana, and the parametric values of iron, fluoride, total hardness, pH, nitrate, and nitrite were used to study the groundwater quality-depth relationship. The results showed that the depth-to-groundwater quality indicated a marginal increase in water quality in the range of 30 to 50 m, which is mathematically represented by the low-value correlation coefficient (r2 = 0.026). A relatively significant increase occurs in the depth range of 50 to 80 m, which is given by a correlation coefficient of r2 = 0.298. The mean percent parameter compatibility was 74%, 82%, 89%, and 97% at 50, 60, 70, and 80 m depths, respectively. The variations in groundwater quality per depth ratio ranged from 1.48, 1.37, 1.27, and 1.21 for 50, 60, 70, and 80 m depth, respectively. The recommended minimum borehole depth for excellent groundwater quality is suggested with a compatibility per meter depth ratio of 1.37. This results in a range between 50 and 70 m as the most desirable drilling depth for excellent groundwater quality within the granitoids of the Birimian Supergroup of the Ashanti Region in Ghana.

Uncertainties of ENSO-related Regional Hadley Circulation Anomalies within Eight Reanalysis Datasets

El Nino-Southern Oscillation(ENSO),the leading mode of global interannual variability,usually intensifies the Hadley Circulation(HC),and meanwhile constrains its meridional extension,leading to an equatorward movement of the jet system.Previous studies have investigated the response of HC to ENSO events using different reanalysis datasets and evaluated their capability in capturing the main features of ENSO-associated HC anomalies.However,these studies mainly focused on the global HC,represented by a zonal-mean mass stream function(MSF).Comparatively fewer studies have evaluated HC responses from a regional perspective,partly due to the prerequisite of the Stokes MSF,which prevents us from integrating a regional HC.In this study,we adopt a recently developed technique to construct the three-dimensional structure of HC and evaluate the capability of eight state-of-the-art reanalyses in reproducing the regional HC response to ENSO events.Results show that all eight reanalyses reproduce the spatial structure of HC responses well,with an intensified HC around the central-eastern Pacific but weakened circulations around the Indo-Pacific warm pool and tropical Atlantic.The spatial correlation coefficient of the three-dimensional HC anomalies among the different datasets is always larger than 0.93.However,these datasets may not capture the amplitudes of the HC responses well.This uncertainty is especially large for ENSO-associated equatorially asymmetric HC anomalies,with the maximum amplitude in Climate Forecast System Reanalysis(CFSR)being about 2.7 times the minimum value in the Twentieth Century Reanalysis(20CR).One should be careful when using reanalysis data to evaluate the intensity of ENSO-associated HC anomalies.