Relationship between gravity change and Yangbi M_(S)6.4 earthquake

Based on the relative and absolute gravity measurements in the southern South-North Seismic Belt since 2015,we analyzed the dynamic change of the regional gravity field and its relationship with the Yangbi M_(S)6.4 earthquake that occurred on May 21,2021.The results show that:(1)The regional gravity field changes are closely related to the Weixi-Qiaohou fault,which reflects the surface gravity field changes caused by the fault activity from 2015 to 2021;(2)The gravity field change related to the preparation of Yangbi earthquake has experienced the evolution process of"steady state-regional gravity anomaly-local gravity anomaly-four-quadrant distribution-large area positive anomaly-earthquake occurring in the reverse change process";(3)The cumulative and differential change images of the gravity field show that there were significant gravity changes in the two years preceding the Yangbi earthquake,and the earthquake occurred in the high-gradient belt of gravity variation,the center of the four-quadrant,and close to the zero contour turn;(4)The dynamic evolution image of the gravity field can well reflect the precursory phenomena during the preparation for the Yangbi MS6.4 earthquake.Based on the anomaly change of mobile gravity,a certain degree of medium-term prediction was made before the Yangbi M_(S)6.4 earthquake,especially the determination of strong earthquake location.

Recent gravity changes in China Mainland

Based on results of the mobile gravity measurements of the Crustal Movement Observation Network of China and Digital Earthquake Observation Network of China, this paper shows the pattern of temporal gravity changes in China mainland on a time scale of 2 - 3 years since 1998, and gives an analysis of the patterns. The result shows that the temporal gravity changes basically reflect the current mass movement and occurrence of strong earthquakes.

Detection of a half-microgal coseismic gravity change after the Ms7. 0 Lushan earthquake

Because only a small near-field coseismie gravity change signal remains after removal of noise from the accuracy of observations and the time and spatial resolution of the earth＇s surface gravity observation system, it is difficult to verify simulations of dislocation theory. In this study, it is shown that the GS15 gravimeter, located 99.5 km from the epicenter of the Ms7.0 Lushan earthquake on April 20, 2013 at 08 ： 04 UTC ＋ 8, showed the influence of the earthquake from 2013-04-16 to 2013-04-26 after a time calibration, tide correc- tions, drift correction, period correction and relaxation correction were applied to its data. The post-seismic relaxation process of the spring in the gravimeter took approximately 430 minutes and showed a 2. 5 ×10^-8 ms^-2 gravity change. After correcting for the relaxation process, it is shown that a coseismic gravity change of approximately ＋0.59 ＋-0. 4 ~ 10-Sms-2 was observed by the GS15 gravimeter; this agrees with the simulated gravity change of approximately 0.31 ~ 10 -8 ms-2. The rate of the coseismie gravity change and the coseismic vertical displacement, as measured by one-second and one-day sampling interval GPS units, is also consistent with the theoretical rate of change. Therefore, the GS15 gravimeter at the Pixian Station observed a coseismic gravity change after the Ms7.0 Lushan earthquake. This and similar measurements could be applied to test and confirm the theory used for these simulations.

Far-field coseismic gravity changes related to the 2015 MW7.8 Nepal(Gorkha)earthquake observed by superconducting gravimeters in Chinese mainland

Using data from five SGs at four stations in Chinese mainland,obvious permanent gravity changes caused by the 2015 MW7.8 Nepal(Gorkha)earthquake were detected.We analyzed the gravity effects from ground vertical deformation(VD)using co-site continuous GPS(cGPS)data collocated at the Lijiang and the Wuhan station,and hydrological effects using GLDAS models and groundwater level records.After removing these effects,SG observations before and after the earthquake revealed obvious permanent gravity changes:−3.0μGal,7.3μGal and 8.0μGal at Lhasa,Lijiang and Wuhan station,respectively.We found that the gravity changes cannot be explained by the results of dislocation theory.

Simulation of co-seismic gravity change and deformation of Wenchuan Ms8.0 earthquake

Surface co-seismic gravity changes and displacements caused by the Wenchuan Ms8.0 earthquake are calculated on the basis of the half-space dislocation theory and two fault models inversed, respectively, by Institute of Geophysics, CEA and USGS. The results show that 1 ） the dislocation consists of dip slip and rightlateral strike slip ;2 ）the co-seismic gravity change shows a four-quadrant pattern ,which is greatly controlled by the distribution of the vertical displacements, especially in the near-filed ; 3 ） the gravity change is generally less than 10 × 10^-8 ms^-2 in the far-field,but as high as several 100 × 10^-8 ms^-2 in the near-filed. These results basically agree with observational results.

GRACE detection of the medium-to far-field coseismic gravity changes caused by the 2004 Mw9.3 Sumatra-Andaman earthquake

Large earthquakes cause observable changes in the Earth’s gravity field, which have been detected by the Gravity Recovery and Climate Experiment (GRACE). Since most previous studies focus on the detection of near-field gravity effects, this study provides the results from the medium- to far-field gravity changes caused by the 2004 Sumatra-Andaman earthquake that are recorded within GRACE monthly solutions. Utilizing a spherical-earth dislocation model we documented that large-scale signals predominate in the global field of the coseismic gravity changes caused by the earthquake. After removing the near-field effects, the coseismic gravity changes show a negative anomaly feature with an average magnitude of -0.18×10-8 m·s-2 in the region ranging ～40° from the epicenter, which is considered as the 'medium ffield' in this study. From the GRACE data released by Center for Space Research from August 2002 to December 2008, we retrieved the large-scale gravity changes smoothed with 3 000 km Gaussian ffilter. The results show that the coseismic gravity changes detected by GRACE in the medium field have an average of (-0.20±0.06)×10-8 m·s-2, which agrees with the model prediction. The detection confirms that GRACE is sensitive to large-scale medium-field coseismic gravitational effects of mega earthquakes, and also validates the spherical-earth dislocation model in the medium field from the perspective of satellite gravimetry.

Gravity change observed in a local gravity network and its implication to seasonal precipitation in Dali county, Yunnan province, China

This study investigates data-processing methods and examines the precipitation effect on gravity measurements at the Dali gravity network, established in 2005. High-quality gravity data were collected during four measurement campaigns. To use the gravity data validly, some geophysical corrections must be considered carefully. We first discuss data-processing methods using weighted least- squares adjustment with the constraint of the absolute gravity datum. Results indicate that the gravity precision can be improved if all absolute gravity data are used as constraints and if calibration functions of relative gravi- meters are modeled within the observation function. Using this data-processing scheme, the mean point gravity pre- cision is better than 12 μgal. After determining the best data-processing scheme, we then process the gravity data obtained in the four measurement campaigns, and obtain gravity changes in three time periods. Results show that the gravity has a remarkable change of more than 50 pgal in the first time period from Apr-May of 2005 to Aug-Sept of 2007. To interpret the large gravity change, a mean water mass change （0.6 m in height） is assumed in the ETOPO1 topographic model. Calculations of the precipitation effect on gravity show that it can reach the same order of the observed gravity change. It is regarded as a main source of the remarkable gravity change in the Dali gravity network, suggesting that the precipitation effect on gravity mea- surements must be considered carefully.

Influence of fault asymmetric dislocation on the gravity changes

A fault is a planar fracture or discontinuity in a volume of rock, across which there has been significant displacement along the fractures as a result of earth movement. Large faults within the Earth＇s crust result from the action of plate tectonic forces, with the largest forming the boundaries between the plates, energy release associated with rapid movement on active faults is the cause of most earthquakes. The relationship between unevenness dislocation and gravity changes was studied on the theoretical thought of differential fault. Simulated observation values were adopted to deduce the gravity changes with the model of asymmetric fault and the model of Okada, respectively. The characteristic of unevenness fault momentum distribution is from two end points to middle by 0 according to a certain continuous functional increase. However, the fault momentum distribution in the fault length range is a constant when the Okada model is adopted. Numerical simulation experiments for the activities of the strike-slip fault, dip-slip fault and extension fault were carried out, respectively, to find that both the gravity contours and the gravity variation values are consistent when either of the two models is adopted. The apparent difference lies in that the values at the end points are 17.97% for the strike-slip fault, 25.58% for the dip-slip fault, and 24.73% for the extension fault.

Trends in gravity changes from 2009 to 2013 derived from ground-based gravimetry and GRACE data in North China

North China is a key region for studying geophysical progress. In this study, ground-based and Gravity Recovery and Climate Experiment（GRACE） gravity data from 2009 to 2013 are used to calculate the gravity change rate（GCR） using the polynomial fitting method. In general, the study area was divided into the Shanxi rift, Jing-Jin-Ji（Beijing-Tianjin-Hebei Province）, and Bohai Bay Basin（BBB） regions. Results of the distribution of the GCR determined from ground-based gravimetry show that the GCR appears to be ＂negativepositive-negative＂ from west to east, which indicates that different geophysical mechanisms are involved in the tectonic activities of these regions. However, GRACE solutions are conducted over a larger spatial scale and are able to show a difference between southern and northern areas and a mass redistribution of land water storage.

Coulomb stress and gravity changes associated with the 2016 Mw 7.8 Kaikoura Earthquake, New Zealand: Application for aftershock triggering and fault interaction process analysis

The Kaikoura earthquake on November 14,2016 is one of the largest and most complex earthquakes in New Zealand since 1947.Despite the fact that it has ruptured about 12 separate faults,triggered 2132 aftershocks within one week of the mainshock and induced considerable stress changes,few studies have been conducted to comprensively investigate the characteristics.The current study examines the horizontal and vertical displacements as well as the stress and gravity changes,aftershock distributions and also find out whether these changes affect the surrounding regions along the complex fault systems.The study covers the entire area affected by the Kaikoura event,which includes the northern part of the South Island and the southern part of the North Island.The dislocation theory was employed to evaluate the coseismic slip model on the multiple faults.The displacement results revealed that the maximum horizontal displacement is about 6 m and the vertical about 2 m,which are reasonably consistent with earlier study findings.Besides,the stress and gravity changes are quite complicated and inhomogeneous as evidenced by our coseismic model,demonstrating the complexity of the Kaikoura earthquake as well.Almost all the aftershocks are distributed in places where the stress and gravity change are found to be significant.In order to investigate the stability of our stress change models,we applied different friction coefficients and receiver fault parameters.The results justify the friction coefficient(μ=0.4)and the receiver fault parameters(230°,70°,150°)are suitable to define good stress change estimates.According to the stress change results at 15 km depth,the northern parts of the mainshock region,Hundalee fault,Humps fault and Jordan thrust areas together with the Wellington area are closer to failure and situated in a seismic risk zone.The multidimensional analysis adopted in this paper is helpful for making decisions and applications of stress and gravity change models in assessing seismic hazards.

Gravity change and its mechanism after the first water impoundment in Three Gorges Project

In this paper we have analyzed precise gravity survey and gravity effects resulted from water loading, crustal deformation, ground water level change and precipitation before and after the water impoundment in the Three Gorges Reservoir. We found that： ① In dam area of the reservoir, gravity effect resulted from water load increase is the most significant, maximum gravity change is 200×10^-8 m/s^2, but this effect is limited in amplitude and range. Gravity change can be observed about 5 km offshore. The gravity change caused by ground water level change is regional; and the impact of precipitation on it should not be neglected. ② At head area of the reservoir, the maximum gravity change is near Xiangxi. Monitoring the variation of gravity field and further study should continue in the future.

Gravity changes and surface deformations due to faults with different geometry

Based on the formulae of the gravity changes and surface deformations raised by the dislocation of a point source,the gravity changes and deformations caused by the dislocations of fault with arbitrary geometry are computed by using numerical method. The results show that both of the dislocation and the geometry of the fault plane are the basic elements that determine the gravity and deformation effects. Gravity changes, vertical deformations and apparent vertical deformations induced by the dislocation are alike in their characteristic patterns. The similarities and differences of these patterns provide us a probability in acquiring the gravity and deformation anomalies due to faulting from the observed data. Thus the geometric and kinematic features of the earthquake-generating faults can be appropriately distinguished and evaluated.

Study of the Gravity Changes before and after the Jingtai M_S5.9 Earthquake

By systematically analyzing the data of gravity reiteration in the Hexi region and taking a dynamic viewpoint, we have studied the evolution characters of gravity field during the preparation-occurrence of the Jingtai M S5.9 earthquake of June 6, 2000. The patterns of dynamic change of the gravity field clearly reflected how the gravity field evolved from the quasi-homogeneous state to non-homogeneous state for earthquake preparation and then the earthquake occurred. Besides, we have also studied the relation between the characteristic gravity change and strong earthquake activity.

Estimating Field Source Parameters of Gravity Change in North China Based on the Euler Deconvolution Method

Based on the absolute and relative gravity observations in North China from 2009 to 2014,spatial dynamic variations of the regional gravity field are obtained. We employed the Euler deconvolution method and the theoretical model to get the best estimates of parameters. Gravity field change caused by the depth and distribution in North China is calculated by back analysis. The results show the structural index that equals 1 is suitable for inversion of the gravity variation data. The inversion results indicate that the depths of anomaly field sources are spread over the Hetao fault. The research method of this paper can be used in the quantitative study on the field source and may shed new light on the interpretations of gravity change, and also provide quantitative basis for earthquake prediction index criterions based on the gravity change.

Gravity Change Caused by Heavy Rainfall Detected by A gPhone Gravimeter in Zhengzhou,China

The short-term effect of heavy rainfall on gPhone gravimeter observation at Zhengzhou Seismic Station is investigated.According to the observation data during Jul.17–20,2021,the corrected gravity residual reflects the gravimetric response caused by heavy rainfall.The observed gravity change is dominated by the local effect considering topographic effect on gravity.The deduced water depth near the observation station is about 300 mm.

Viscous lithospheric structure beneath Sumatra inferred from post-seismic gravity changes detected by GRACE

Satellite gravity data of the regional rheological structure of the lithosphere in the vicinity of Sumatra is used as evidence to enable a better understanding of the regional geodynamic environment.The data is interpreted using the theory of post-seismic viscoelastic relaxation.Coand post-seismic changes in the gravity field resulting from the 2004 M w 9.3 Sumatra earthquake were calculated from Gravity Recovery and Climate Experiment (GRACE) satellite data.A spatial Gaussian filter,500 km wide,was used in the calculation.The results indicate that there were significant co-seismic jumps in both uplifted and subducted regions.The magnitude of the jump in the subducted zone was ～9 10 8 m/s 2,more significant than the ～2 10 8 m/s 2 jump observed in the uplifted zone.However,a positive gravity change occurred in the uplifted zone very soon after the earthquake.The rheological structure of the lithosphere has a great effect on deformation and its determination is a fundamental part of developing reliable numerical simulations in geodynamics.Based on the temporally-variable gravity field observed by GRACE,the viscous lithospheric structure of the Sumatra area is investigated with a self-gravitating,half space,viscoelastic earth model.The estimated viscosity is of the order of 1.0 10 18 Pa·s and there are differences in the rheological parameters on the two sides of the fault.The factors that affect the viscosity are discussed in connection with the tectonic structure of the Sumatra area.

Coseismic gravity and displacement changes of Japan Tohoku earthquake(Mw 9.0)

The greatest earthquake in the modern history of Japan and probably the fourth greatest in the last 100 years in the world occurred on March 11, 2011 off the Pacific coast of Tohoku.Large tsunami and ground motions caused severe damage in wide areas, particularly many towns along the Pacific coast. So far, gravity change caused by such a great earthquake has been reported for the 1964 Alaska and the 2010 Maule events. However, the spatial-temporal resolution of the gravity data for these cases is insufficient to depict a co-seismic gravity field variation in a spatial scale of a plate subduction zone. Here, we report an unequivocal co-seismic gravity change over the Japanese Island, obtained from a hybrid gravity observation（combined absolute and relative gravity measurements）. The time interval of the observation before and after the earthquake is within 1 year at almost all the observed sites, including 13 absolute and 16 relative measurement sites, which deduced tectonic and environmental contributions to the gravity change. The observed gravity agrees well with the result calculated by a dislocation theory based on a self-gravitating and layered spherical earth model. In this computation, a co-seismic slip distribution is determined by an inversion of Global Positioning System（GPS） data. Of particular interest is that the observed gravity change in some area is negative where a remarkable subsidence is observed by GPS, which can not be explained by simple vertical movement of the crust. This indicated that the mass redistribution in the underground affects the gravity change. This result supports the result that Gravity Recovery and Climate Experiment（GRACE） satellites detected a crustal dilatation due to the 2004 Sumatra earthquake by the terrestrial observation with a higher spatial and temporal resolution.

Surface deformation and gravity changes caused by dilatancy in a layered elastic-vis-coelastic half space

Surface deformation and gravity changes caused by dilatancy with magmatic intrusion (Mogi model) in two elastic layers overlying a Maxwell viscoelastic half space were systematically studied. We focused on the numerical calculation technique and the discussion of the effects of rheological properties of Maxwell half spqce on vertical displacements and gravity changes. The results show that the rheological properties of materials has an apparent effect on the results, when the intrusion was placed under Moho discontinuity, and has smaller effect if the intrusion was placed inside the crust. Therefore the elastic model is valid for considering the gravity changes or deformation caused by dilatancy with magmatic intrusion within the crust. The model and numerical calculation technique in this study can be used to model and interpret the long term deformation and gravity observations in volcanic, seismic and geothermal regions.

Time-varying gravity field model of Sichuan-Yunnan region based on the equivalent mass source model

High-precision time-varying gravity field is an effective way to study the internal mass movement and understanding the spatio-temporal evolution process of the geodynamic system.Compared to the satellite gravity measurement,the repeated terrestrial gravity observation can provide a more high-order signal related to the shallow crust and subsurface.However,the suitable and unified method for gravity model estimation is a key problem for further applications.In this study,we introduce the spherical hexahedron element to simulate the field source mass and forward model the change of gravity field located at the Sichuan-Yunnan region(99—104°E,23—29°N)in the four epochs from 2015 to 2017.Compared to the experimental results based on Slepian or spherical harmonics frequency domain method,this alternative approach is suitable for constructing the equivalent mass source model of regional-scale gravity data,by introducing the first-order smooth prior condition of gravity time-varying signal to suppress the high-frequency component of the signal.The results can provide a higher spatial resolution reference for regional gravity field modeling in the Sichuan-Yunnan region.

Water storage changes in North America retrieved from GRACE gravity and GPS data

As global warming continues,the monitoring of changes in terrestrial water storage becomes increasingly important since it plays a critical role in understanding global change and water resource management.In North America as elsewhere in the world,changes in water resources strongly impact agriculture and animal husbandry.From a combination of Gravity Recovery and Climate Experiment（GRACE） gravity and Global Positioning System（GPS） data,it is recently found that water storage from August,2002 to March,2011 recovered after the extreme Canadian Prairies drought between 1999 and 2005.In this paper,we use GRACE monthly gravity data of Release 5 to track the water storage change from August,2002 to June,2014.In Canadian Prairies and the Great Lakes areas,the total water storage is found to have increased during the last decade by a rate of 73.8 ± 14.5 Gt/a,which is larger than that found in the previous study due to the longer time span of GRACE observations used and the reduction of the leakage error.We also find a long term decrease of water storage at a rate of-12.0 ± 4.2 Gt/a in Ungava Peninsula,possibly due to permafrost degradation and less snow accumulation during the winter in the region.In addition,the effect of total mass gain in the surveyed area,on present-day sea level,amounts to-0.18 mm/a,and thus should be taken into account in studies of global sea level change.