The coseismic displacements of the 2013 Lushan Mw6.6 earthquake determined using continuous global positioning system measurements

Based on Continuous GPS （CGPS） observation data of the Crustal Movement Observation Network of China （CMONOC） and the Sichuan Continuous Operational Reference System （ SCCORS）, we calculated the horizontal coseismic displacements of CGPS sites caused by the 2013 Lushan Mw 6.6 earthquake. The resuits indicate that the horizontal coseismic deformations of CGPS stations are consistent with thrust-compression rupture. Furthermore, the sites closest to the epicenter underwent significant coseismic displacements. Three network stations exhibited displacements greater than 9 mm （ the largest is 20.9 mm at SCTQ） , while the others were displaced approximately 1 -4 mm.

Coseismic displacements and inospheric changes of the 2013 Ms7. 0 Lushan earthquake from GPS measurements

By inverting GPS data recorded at stations of the Crustal Movement Observation Network of China （CMONOC） near the 2013 Lushan Ms7.0 earthquake, we found a horizontal displacement of 22 mm at a site about 32 kin SW of the epicenter and vertical displacements of as much as 12.4 mm at several sites. The vertical displacements were generally uplift on the west side of the nearby Longmenshan fault zone and subsidence on the east side. We also found coseismic ionospheric disturbances about 0.5 to 0.9 TECU in amplitude that lasted for about one hour.

Co and postseismic fault slip models of the 2022 M_(W)6.7 Menyuan earthquake reveal conjugated faulting tectonics at the central section of the Lenglongling fault

The 2022 M_(W)6.7 Menyuan earthquake ruptured the western end of the Tianzhu seismic gap,providing an opportunity to study the regional seismogenic characteristics and seismic hazards.Here we use interferometric synthetic aperture radar(InSAR)and seismic data to study the mainshock rupture,early afterslip and the second largest aftershock of the 2022 Menyuan earthquake sequences.Our modeling results show that the mainshock ruptured the Lenglongling fault and the Tuolaishan fault with a maximum slip of~3 m.Rapid postseismic transient deformation occurred at the center of the Lenglongling fault.Our afterslip modeling reveals that the majority of afterslip occurred in the deeper part of the Lenglongling fault.A high-angle conjugated faulting event is found at the middle section of the Lenglongling fault.We use the stress inversion to investigate the possible triggering mechanism of the conjugated rupture event.The results indicate the maximum principal stress direction is in~222°,forming a~22°angle between the conjugated fault of second largest aftershock and the mainshock.The calculated normal stress changes indicate the region is within a pull-apart stress field,which favors such a conjugated rupturing event.Our study will help understand the rupture behavior of such kind of conjugated fault in other regions.

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.

A method to retrieve coseismic displacement by integrating GPS and strong motion data

We introduce a method of integrating Global Positioning System(GPS)and accelerometer data for high-rate seismogeodesy.This method is based on the GPS variometric approach which can obtain seismic waves in real time using only the readily available broadcast products and a single receiver.Collocated 5 Hz GPS and 200 Hz accelerometer data from the 2010 M_(W) 7.2 EI Mayor-Cucapaph earthquake were analyzed to verify the effectiveness of this method.Results reveal that this method can provide broadband and more accurate displacements qualified to avoid the baseline drifting caused by strong motion.Moreover,this method can effectively avoid the aliasing problem present in the 5-Hz GPS waveforms.We therefore conclude that this method can be a powerful tool to capture seismic waves in real time,which is crucial to tsunami early warning and earthquake rapid response.

Recent advances of computing coseismic deformations in theory and applications

This paper reviews the recent advances in computing coseismic deformations,and their contributions to seismology and geodesy. At first,an overview on the history of the dislocation theory development is given in the introduction section. Then,emphasis are given on some new developments through few examples in the following sections,such as the new dislocation theory for a 3D Earth model,a new computing scheme on coseismic deflection change of vertical,the relation of dislocation Love number and the conventional Love numbers,the application of dislocation theory applied in satellite gravity observations,the coseismic deformations observed by GRACE,and a new method to determine dislocation Love numbers by GRACE. Furthermore,some advanced theoretical and cases studies are introduced to illustrate how dislocation theory is important in interpret geodetic data,or invert seismic slip for co- and post-seismic processes,using seismic and geodetic data. Final remarks are given in the last section,with discussions,conclusions,comments on existing problems,and expected methods to solve them.

Numerical Simulation on Coseismic Effect of the November 14,2001 Great Kunlun Earthquake,Northern Tibet,China

The November 14,2001 M_S8.1 Kunlun Mountains earthquake in northern Tibet is the largest earthquake occurring on the Chinese mainland since 1950.We apply a three-dimensional(3-D)finite element numerical procedure to model the coseismic displacement and stress fields of the earthquake based on field investigations.We then further investigate the stress interaction between the M_S8.1 earthquake and the intensive aftershocks.Our primary calculation shows that the coseismic displacement field is centralized around the east Kunlun fault zone.And the attenuation of coseismic displacements on the south side of Kunlun fault zone is larger than that on the north side.The calculated coseismic stress field also indicates that the calculated maximal shear stress field is centralized around the east Kunlun fault zone;the directions of the coseismic major principal stress are opposite to that of the background crustal stress field of the Qinghai-Xizang(Tibet)Plateau.It indicates that the earthquake relaxes the crustal stress state in the Qinghai-Xizang(Tibet)Plateau.Finally,we study the stress interaction between M_S8.1 earthquake and its intensive aftershocks.The calculated Coulomb stress changes of the M_S8.1 great earthquake are in favor of triggering 4 aftershocks.

Measuring ground deformations caused by 2015 Mw7.8 Nepal earthquake using high-rate GPS data

The April 25, 2015 Mw7.8 Nepal earthquake was successfully recorded by Crustal Movement Observation Network of China （CMONOC） and Nepal Geodetic Array （NGA）. We processed the high-rate GPS data （1 Hz and 5 Hz） by using relative kinematic positioning and derived dynamic ground motions caused by this large earthquake. The dynamic displacements time series clearly indicated the displacement amplitude of each station was related to the rupture directivity. The stations which located in the di- rection of rupture propagation had larger displacement amplitudes than others. Also dynamic ground displacement exceeding 5 cm was detected by the GPS station that was 2000 km away from the epicenter. Permanent coseismic displacements were resolved from the near-field high-rate GPS stations with wavelet decomposition-reconstruction method and P-wave arrivals were also detected with S transform method. The results of this study can be used for earthquake rupture process and Earthquake Early Warning studies.

Effects of Earth's curvature and radial heterogeneity in dislocation studies:Case studies of the 2008 Wenchuan earthquake and the 2004 Sumatra earthquake

Recently, effects of Earth＇s curvature and radial heterogeneity on coseismic deformations are often investi- gated based on the 2004 Sumatra earthquake. However, such effects are strongly related to earthquake types. As a low dip angle event, the 2004 Sumatra earthquake is not a good seismic case for such a topic since the effects for moderate dip angle events are much bigger. In this study, the half-space and spherical dislocation theories are used, respectively, to calculate co- seismic displacements caused by the 2008 Wenchuan earthquake and the 2004 Sumatra earthquake. Effects of Earth＇s curva- ture and stratification are investigated through the discrepancies of results calculated using the two dislocation theories. Re- sults show that the effects of Earth＇s curvature and stratification for the 2008 Wenchuan earthquake are much larger than those for the 2004 Sumatra earthquake. Ignoring the effects will cause errors up to 100%-200% in far field displacements for a moderate dip angle event like the 2008 Wenchuan earthquake. Such great effects are much bigger than those conclusions of previous studies. Besides, comparison with observations verifies that spherical dislocation theories yield better results than half-space ones in far fields.

Dynamic evolution of crustal horizontal deformation before the Ms6.4 Menyuan earthquake

An Ms6.4 earthquake occurred in the Menyuan county of Qinghai Province on Jan 21, 2016. In order to recognize the development of horizontal deformation and distinguish precursory deformation anomalies, we obtained coordinates time series, velocity and strain model around the seismic zones from processing of continuous observations from 2010 and 6 times of surveying Global Positioning System （GPS） data since 2009. The results show that, before the earthquake, the eastern segmentation of the Qilian tectonic zone where the Lenglongling Fault located is in strong crustal shortening and compressional strain state with dilatational rates of -15 to -25 ppb. The Lenglongling Fault has a strike-slip rate of 3.1 mm/a and a far-field differential orthogonal rate of 7 mm/a, while differential rate is only 1.2 mm/a near the fault, which reflects its locking feature with strain energy accumulation and high seismic risks. Dynamic evolution of deformation model shows that preevent dilatational rates around the seismic zones increases from 15 ppb/a to -20 ppb/a with its center moving to the source areas. Time series of N components of G337 station, which is 13.7 km away from the Lenglongling Fault, exhibit a 5 mm/a acceleration anomaly. Time series of base-station QHME （in Menyuan） displays a reverse acceleration from the end of Sep. to Dec., 2016 when it comes to a largest deviation, and the accumulative displacement is more than 4 mm and the value reverse till the earthquake. In our results, coseismic displacement of N, E, U components in QHME site are 3.0 mm, 3.0 mm, -5.4 mm, respectively. If we profile these values onto the Lenglongling Fault, we can achieve a 1.1 mm of strike slip and 4.1 mm updip slip relative to the hanging wall.

The Study of Paleoearthquakes on the Weihe Fault Zone

It is indicated by historical records and the exploratory trench on the Weihe fault that the Yaodian-Zhangjiawan segment of the Weihe fault zone has experienced a historical earthquake and 3 paleoearthquake events in the past 9110a. The historical earthquake, namely, event Ⅳ, occurred between 1487 and 1568 AD. The date of paleoseismic event Ⅰ is （9110 ＋ 90） a, and the ages of events Ⅱ and Ⅲ are unknown. The coseismic vertical displacement of events Ⅰ, Ⅱ and Ⅲ is 0.5m, 0.5m and 0.2m, respectively. The exploratory trench also indicates that the Yaodian-Zhangjiawan segment of the Weihe fault was active in the Holocene.

The relationship between three-dimensional coseismic displacement and distribution of coseismic landslides

In mountainous areas,landslides induced by destructive earthquakes are one of the main causes of human casualties,which is an important link in the chain of earthquake hazards.Earthquake-triggered landslides are mainly controlled by three factors,namely seismic property,topography,and geology.Many studies have been conducted on these controlling factors of earthquake-triggered landslides.However,little is known about the effect of coseismic displacement on the distribution of landslides under different slope aspects and slope angles,hindering our understanding of the mechanism of inducing landslides by the combination of surface displacement and slope geometry at the local scale and leading to controversial opinions about the abnormal number of earthquake-triggered landslides in several cases.Here,we took the 2008 Wenchuan M_(w) 7.9 earthquake in China,the 2015 Gorkha M_(w) 7.8 earthquake in Nepal,and the 2016 Kaikōura M_(w) 7.8 earthquake in New Zealand as examples to investigate the relationship between the distribution of large earthquake-triggered landslides and the three-dimensional (3D)coseismic displacement field.We divided the landslide-prone area around the epicenter into regular grids and calculated the 3D coseismic displacement in each grid according to the radar satellite images and slip distribution model.Then,the 3D coseismic displacement was projected to two coordinate systems related to the slope where the landslides were located for statistical analysis.We determined that the surface uplift perpendicular to the slope is more likely to induce landslides,particularly when combined with large slope angles.Meanwhile,the number of landslides will be significantly reduced where the subsidence occurs.Regardless of uplift or subsidence,landslides are more likely to occur when the direction of coseismic horizontal displacement is far from the slope.The larger the slope angles are,the greater the effects of horizontal displacement and slope aspect.A dominant slope aspect also exists for earthquake-triggered landslides,which is different from the mean slope aspect calculated from the background topography.This dominant aspect angle is related to the focal mechanism and striking angle of surface rupture.These results indicate that we can simulate the 3D coseismic displacement field from known fault location and earthquake mechanism and combine the topographic data for landslide risk assessment in earthquake-prone mountainous areas to minimize the damage caused by possible earthquake-triggered landslides.

Coseismic Slip Distribution of the 2021 Mw7.4 Maduo,Qinghai Earthquake Estimated from InSAR and GPS Measurements

On 22 May 2021,the Maduo Earthquake occurred on a branch fault of the East Kunlun fault in the Bayan Har Block,which provides opportunity to constrain fault geometry and strain accumulation and release for assessment of earthquake hazards.We processed the Sentinal-1A/B SAR images acquired before and after the earthquake,with which we constrained a finite fault model to best fit to the combined data set of downsampled InSAR image and GPS displacements.The inversion indicates that the Maduo event ruptured a 160 km long section striking 286.5°and a dipping 81.39°with rake angle of 4.62°.The model suggests three compact rupture areas with the slip amplitude exceeding 4 m on the main rupture section and the largest slip region is in the east of the epicenter with a slip of approximately 4.6 m below the surface,in a good agreement with the field geological survey.The total geodetic moment is 1.67×10^(20) N·m equivalent to Mw7.44,slightly larger than estimate of the USGS.