An SPH Framework for Earthquake-Induced Liquefaction Hazard Assessment of Geotechnical Structures

Earthquake-induced soil liquefaction poses significant risks to the stability of geotechnical structures worldwide.An understanding of the liquefaction triggering,and the post-failure large deformation behaviour is essential for designing resilient infrastructure.The present study develops a Smoothed Particle Hydrodynamics(SPH)framework for earthquake-induced liquefaction hazard assessment of geotechnical structures.The coupled flowdeformation behaviour of soils subjected to cyclic loading is described using the PM4Sand model implemented in a three-phase,single-layer SPH framework.A staggered discretisation scheme based on the stress particle SPH approach is adopted to minimise numerical inaccuracies caused by zero-energy modes and tensile instability.Further,non-reflecting boundary conditions for seismic analysis of semi-infinite soil domains using the SPH method are proposed.The numerical framework is employed for the analysis of cyclic direct simple shear test,seismic analysis of a level ground site,and liquefaction-induced failure of the Lower San Fernando Dam.Satisfactory agreement for liquefaction triggering and post-failure behaviour demonstrates that the SPH framework can be utilised to assess the effect of seismic loading on field-scale geotechnical structures.The present study also serves as the basis for future advancements of the SPH method for applications related to earthquake geotechnical engineering.

Frequency characteristics and far-field effect of gravity perturbation before earthquake

We used high-pass filtering and the Fourier transform to analyze tidal gravity data prior to five earthquakes from four superconducting gravity stations around the world. A stable gravitational perturbation signal is received within a few days before the earthquakes. The gravitational perturbation signal before the Wenchuan earthquake on May 12, 2008 has main frequency of 0.1–0.3 Hz, and the other four have frequency bands of 0.12-0.17 Hz and 0.06-0.085 Hz. For earthquakes in continental and oceanic plate fault zones, gravity anomalies often appear on the superconducting gravimeters away from the epicenter, whereas the stations near the epicenter record small or no anomalies. The results suggest that this kind of gravitational perturbation signals correlate with earthquake occurrence, making them potentially useful earthquake predictors. The far-field effect of the gravitational perturbation signals may reveal the interaction mechanisms of the Earth’s tectonic plates. However, owing to the uneven distribution of gravity tide stations, the results need to be further confirmed in the future.

Seismic Analysis of Reinforced Concrete Silos under Far-Field and Near-Fault Earthquakes

Silos are strategical structures used to stockpile various types of granular materials.They are highly vulnerable to earthquake excitation and have been frequently reported to fail at a higher rate than any other industrial structure.The seismic response of silos within the near-fault region will suffer a complex combination of loadings due to the unique characteristics of the near-fault ground motions;which are usually associated with a large amplitude pulse at the beginning of either the velocity or the displacement time histories.This study aims to numerically evaluate the seismic response of reinforced concrete cylindrical silos under near-fault ground motions(NFGM)and far-field ground motions(FFGM).The assessment investigates the impact of the slenderness ratio and the para-meters’influence on the seismic behavior of reinforced concrete silos.The validity of the Eurocode provisions in the structural safety of silos will also be inspected.The nonlinear time history analysis is carried out through the finite element approach by examining four silos with different slenderness ratios.The concrete damage plas-ticity model is assigned to the silo wall to simulate the nonlinear behavior of concrete in the plastic zone;while,the behavior of the stored material is represented by the Drucker-Prager plastic model.The wall-granular material interaction is considered and defined by coulomb’s friction theory.The results of the near-fault records reveal a growth up to 72.8%in the hoop stress and 160.4%in the vertical stress compared to the far-field earthquakes.Consequently,the seismic response of reinforced concrete silos is highly sensitive to the type of ground motion,and slender silos tend to impose greater structural demand under the NFGM.Additionally,The Eurocode-8 seismic provisions were adequate in the conventional far-field ground motions and less effective in the near-fault zone.

The MW5.5 earthquake on August 6,2023,in Pingyuan,Shandong,China:A rupture on a buried fault

On August 6,2023,a magnitude MW5.5 earthquake struck Pingyuan County,Dezhou City,Shandong Province,China.This event was significant as no large earthquakes had been recorded in the region for over a century,and no active fault had been previously identified.This study collects 1309 P-wave arrival times and 866 S-wave arrival times from 74 seismic stations less than 200 km to the epicenter to constrain the spatial distribution of the mainshock and its 125 early aftershocks by the double difference earthquake relocation method,and selects 864 P-waveforms from 288 stations located within 800 km of the epicenter to constrain the focal mechanism solution of the mainshock through centroid moment tensor inversion.The relocation and the inversion indicate,the Pingyuan MW5.5 earthquake was caused by a rupture on a buried fault,likely an extensive segment of the Gaotang fault.This buried fault exhibited a dip of approximately 75°to the northwest,with a strike of 222°,similar to the Gaotang fault.The rupture initiated at the depth of 18.6 km and propagated upward and northeastward.However,the ground surface was not broken.The total duration of the rupture was~6.0 s,releasing the scalar moment of 2.5895×1017 N·m,equivalent to MW5.54.The moment rate reached the maximum only 1.4 seconds after the rupture initiation,and the 90%scalar moment was released in the first 4.6 s.In the first 1.4 seconds of the rupture process,the rupture velocity was estimated to be 2.6 km/s,slower than the local S-wave velocity.As the rupture neared its end,the rupture velocity decreased significantly.This study provides valuable insights into the seismic characteristics of the Pingyuan MW5.5 earthquake,shedding light on the previously unidentified buried fault responsible for the seismic activity in the region.Understanding the behavior of such faults is crucial for assessing seismic hazards and enhancing earthquake preparedness in the future.

Rapid report of source parameters of 2023 M6.2 Jishishan,Gansu earthquake sequence

The M6.2 earthquake in Jishishan,Gansu Province,on December 18,2023,caused extraordinary earthquake disasters.It was located in the northern part of the north−south seismic zone,which is a key area for earthquake monitoring in China.The newly built dense strong motion stations in this area provide unprecedented conditions for high-precision earthquake relocation,especially the earthquake focal depth.This paper uses the newly built strong motion and traditional broadband seismic networks to relocate the source locations of the M3.0 and above aftershocks and to invert their focal mechanisms.The horizontal error of earthquake location is estimated to be 0.5−1 km,and the vertical error is 1−2 km.The focal depth range of aftershocks is 9.6−14.6 km,distributed in a 12-km-long strip with SSE direction.Aftershocks in the south are more concentrated horizontally and vertically,while aftershocks in the north are more scattered.The focal mechanisms of the main shock and aftershocks are relatively consistent,and the P-axis orientation is consistent with the regional strain direction.There is a seismic blank area of M3.0 and above,about 3−5 km between the main shock and aftershocks.It is suggested that the energy released by the main shock rupture is concentrated in this area.Based on the earthquake location and focal mechanism of the main shock,it is inferred that the Northern Lajishan fault zone is the seismogenic structure of the main shock,and the main shock did not occur on the main fault,but on a secondary fault.The initial rupture depth and centroid depth of the main shock were 12.8 and 14.0 km,respectively.The source rupture depth may not be the main reason for the severe earthquake disaster.

Characterization and spatial analysis of coseismic landslides triggered by the Luding Ms 6.8 earthquake in the Xianshuihe fault zone, Southwest China

On September 5, 2022, a magnitude Ms 6.8 earthquake occurred along the Moxi fault in the southern part of the Xianshuihe fault zone located in the southeastern margin of the Tibetan Plateau,resulting in severe damage and substantial economic loss. In this study, we established a coseismic landslide database triggered by Luding Ms 6.8 earthquake, which includes 4794 landslides with a total area of 46.79 km^(2). The coseismic landslides primarily consisted of medium and small-sized landslides, characterized by shallow surface sliding. Some exhibited characteristics of high-position initiation resulted in the obstruction or partial obstruction of rivers, leading to the formation of dammed lakes. Our research found that the coseismic landslides were predominantly observed on slopes ranging from 30° to 50°, occurring at between 1000 m and 2500 m, with slope aspects varying from 90° to 180°. Landslides were also highly developed in granitic bodies that had experienced structural fracturing and strong-tomoderate weathering. Coseismic landslides concentrated within a 6 km range on both sides of the Xianshuihe and Daduhe fault zones. The area and number of coseismic landslides exhibited a negative correlation with the distance to fault lines, road networks, and river systems, as they were influenced by fault activity, road excavation, and river erosion. The coseismic landslides were mainly distributed in the southeastern region of the epicenter, exhibiting relatively concentrated patterns within the IX-degree zones such as Moxi Town, Wandong River basin, Detuo Town to Wanggangping Township. Our research findings provide important data on the coseismic landslides triggered by the Luding Ms 6.8 earthquake and reveal the spatial distribution patterns of these landslides. These findings can serve as important references for risk mitigation, reconstruction planning, and regional earthquake disaster research in the earthquake-affected area.

Rapid report of the December 18,2023 M_(S)6.2 Jishishan earthquake,Gansu,China

On December 18,2023,the Jishishan area in Gansu Province was jolted by a M_(S) 6.2 earthquake,which is the most powerful seismic event that occurred throughout the year in China.The earthquake occurred along the NWtrending Lajishan fault(LJSF),a large tectonic transformation zone.After this event,China Earthquake Networks Center(CENC)has timely published several reports about seismic sources for emergency responses.The earthquake early warning system issued the first alert 4.9 s after the earthquake occurrence,providing prompt notification that effectively mitigated panics,injuries,and deaths of residents.The near real-time focal mechanism solution indicates that this earthquake is associated with a thrust fault.The distribution of aftershocks,the rupture process,and the recorded amplitudes from seismic monitoring and GNSS stations,all suggest that the mainshock rupture predominately propagates to the northwest direction.The duration of the rupture process is~12 s,and the largest slip is located at approximately 6.3 km to the NNW from the epicenter,with a peak slip of 0.12 m at~8 km depth.Seismic station N0028 recorded the highest instrumental intensity,which is 9.4 on the Mercalli scale.The estimated intensity map shows a seismic intensity reaching up to IX near the rupture area,consistent with field survey results.The aftershocks(up to December 22,2023)are mostly distributed in the northwest direction within~20 km of the epicenter.This earthquake caused serious casualties and house collapses,which requires further investigations into the impact of this earthquake.

Rapid rupture characterization for the 2023 M_(S)6.2 Jishishan earthquake

On December 18, 2023, the M_(S)6.2 Jishishan earthquake occurred in the northeastern region of the QinghaiXizang Plateau, causing heavy casualties and property damage in Gansu and Qinghai Provinces. In this study,we integrate space imaging geodesy, finite fault inversion, and back-projection methods to decipher its rupture property, including fault geometry, coseismic slip distribution, rupture direction, and propagation speed. The results reveal that the seismogenic fault dips to the southwest at an angle of 29°. The major slip asperity is dominated by reverse slip and is concentrated within a depth range of 7–16 km, which explains the significant uplift near the epicenter observed by both the Sentinel-1 ascending and descending In SAR data. Moreover, the teleseismic array waveforms indicate a northwest propagating rupture with an overall slow rupture velocity of~1.91 km/s(AK array) or 1.01 km/s(AU array).

Far-field crustal movements before and after the 2011 Ms9.0 Japan earthquake from GPS observations

The paper gives an analysis of the displacement time series before and after the March 11,2011 Ms9.0 east Japan earthquake and co-seismic displacements observed at continuous GPS stations in and around China. The results showed a broad-scaled related elastic-rebound process and some premonitory horizontal crustal movements to this earthquake over this vast area.

Microearthquake reveals the lithospheric structure at midocean ridges and oceanic transform faults

Mid-ocean ridge and oceanic transforms are among the most prominent features on the seafloor surface and are crucial for understanding seafloor spreading and plate tectonic dynamics,but the deep structure of the oceanic lithosphere remains poorly understood.The large number of microearthquakes occurring along ridges and transforms provide valuable information for gaining an indepth view of the underlying detailed seismic structures,contributing to understanding geodynamic processes within the oceanic lithosphere.Previous studies have indicated that the maximum depth of microseismicity is controlled by the 600-℃isotherm.However,this perspective is being challenged due to increasing observations of deep earthquakes that far exceed this suggested isotherm along mid-ocean ridges and oceanic transform faults.Several mechanisms have been proposed to explain these deep events,and we suggest that local geodynamic processes(e.g.,magma supply,mylonite shear zone,longlived faults,hydrothermal vents,etc.)likely play a more important role than previously thought.

Seismic response of a mid-story isolated structure considering SSI in mountainous areas under long-period earthquakes

At present,there is not much research on mid-story isolated structures in mountainous areas.In this study,a model of a mid-story isolated structure considering soil-structure interaction(SSI)in mountainous areas is established along with a model that does not consider SSI.Eight long-period earthquake waves and two ordinary earthquake waves are selected as inputs for the dynamic time history analysis of the structure.The results show that the seismic response of a mid-story isolated structure considering SSI in mountainous areas can be amplified when compared with a structure that does not consider SSI.The structure response under long-period earthquakes is larger than that of ordinary earthquakes.The structure response under far-field harmonic-like earthquakes is larger than that of near-fault pulse-type earthquakes.The structure response under near-fault pulse-type earthquakes is larger than that of far-field non-harmonic earthquakes.When subjected to long-period earthquakes,the displacement of the isolated bearings exceeded the limit value,which led to instability and overturning of the structure.The structure with dampers in the isolated story could adequately control the nonlinear response of the structure,effectively reduce the displacement of the isolated bearings,and provide a convenient,efficient and economic method not only for new construction but also to retrofit existing structures.

Strong ground motion characteristics observed in the February 6,2023 M_(W)7.7 Türkiye earthquake

Türkiye is located in a seismically active region,where the Anatolian,African,and Arabian tectonic plates converge.High seismic hazards cause the region to be struck repeatedly by major earthquakes.On February 06,2023,a devastating M_(W)7.7 earthquake struck Türkiye at 01:17 am local time(01:17 UTC).In this regard,near and far-field ground motion data within the distance of 120 km are compiled and later characterized to identify the key ground motion intensity measures.Additionally,the vertical components of ground motions were examined to capture the complete three-dimensional nature of the seismic event.Moreover,the effect of Pulse-Like(PL)and Non-Pulse-Like(NPL)ground motion on a representative RC frame structure built as per the Türkiye code was investigated.The results indicate that PL behavior was observed in both horizontal and vertical components of ground motions and PL behavior were noted both near the epicenter and at higher distances from the epicenter.Moreover,the ratio of the peak vertical acceleration to peak horizontal acceleration at certain stations was found to be close to 1.Finally,the non-linear time history analysis of the representative reinforced concrete frame structure for ground motions recorded at stations located equidistant from the epicenter,indicated that PL ground motions led to more significant damage compared to NPL ground motions.

The September 16, 2015 Mw 8.3 Illapel, Chile Earthquake： characteristics of tsunami wave from near-field to far-field

On September 16, 2015, an earthquake with magnitude ofMw 8.3 occurred 46 km offshore from Illapel, Chile, generating a 4.4-m local tsunami measured at Coquimbo. In this study, the characteristics of tsunami are presented by a combination of analysis of observations and numerical simulation based on sources of USGS and NOAA. The records of 16 DART buoys in deep water, ten tidal gauges along coasts of near-field, and ten coastal gauges in the far-field are studied by applying Fourier analyses. The numerical simulation based on nonlinear shallow water equations and nested grids is carried out to provide overall tsunami propagation scenarios, and the results match well with the observations in deep water and but not well in coasts closed to the epicenter. Due to the short distance to the epicenter and the shelf resonance of southern Peru and Chile, the maximum amplitude ranged from 0.1 m to 2 m, except for Coquimbo. In deep water, the maximum amplitude of buoys decayed from 9.8 cm to 0.8 cm, suggesting a centimeter-scale Pacific-wide tsunami, while the governing period was 13-17 min and 32 min. Whereas in the far-field coastal region, the tsunami wave amplified to be around 0.2 m to 0.8 m, mostly as a result of run-up effect and resonance from coast reflection. Although the tsunami was relatively moderate in deep water, it still produced non-negligible tsunami hazards in local region and the coasts of farfield.

Coseismic deformation and fault slip distribution of the 2023 M_(W)7.8 and M_(W)7.6 earthquakes in Türkiye

On February 6,2023,a devastating earthquake with a moment magnitude of M_(W)7.8 struck the town of Pazarcik in south-central Türkiye,followed by another powerful earthquake with a moment magnitude of M_(W)7.6 that struck the nearby city of Elbistan 9 h later.To study the characteristics of surface deformation caused by this event and the influence of fault rupture,this study calculated the static coseismic deformation of 56 stations and dynamic displacement waveforms of 15 stations using data from the Turkish national fixed global navigation satellite system(GNSS)network.A maximum static coseismic displacement of 0.38 m for the M_(W)7.8 Kahramanmaras earthquake was observed at station ANTE,36 km from the epicenter,and a maximum dynamic coseismic displacement of 4.4 m for the M_(W)7.6 Elbistan earthquake was observed at station EKZ1,5 km from the epicenter.The rupture-slip distributions of the two earthquakes were inverted using GNSS coseismic deformation as a constraint.The results showed that the Kahramanmaras earthquake rupture segment was distinct and exposed on the ground,resulting in significant rupture slip along the Amanos and Pazarcik fault segments of the East Anatolian Fault.The maximum slip in the Pazarcik fault segment was 10.7 m,and rupture occurred at depths of 0–15 km.In the Cardak fault region,the Elbistan earthquake caused significant ruptures at depths of 0–12 km,with the largest amount of slip reaching 11.6 m.The Coulomb stress change caused by the Kahramanmaras earthquake rupture along the Cardak fault segment was approximately 2 bars,and the area of increased Coulomb stress corresponded to the subsequent rupture region of the M_(W)7.6 earthquake.Thus,it is likely that the M_(W)7.8 earthquake triggered or promoted the M_(W)7.6 earthquake.Based on the cumulative stress impact of the M_(W)7.8 and M_(W)7.6 events,the southwestern segment of the East Anatolian Fault,specifically the Amanos fault segment,experienced a Coulomb rupture stress change exceeding 2 bars,warranting further attention to assess its future seismic hazard risk.

Variations of shear-wave splitting parameters in the source region of the 2023 Türkiye doublet earthquakes

In this study,the shear-wave splitting parameters of local seismic events from the source regions of the 2023 Türkiye MW7.7 and MW7.6 doublet earthquakes(event 1 and event 2,respectively)were measured from June 1,2022,to April 25,2023,and their spatiotemporal characteristics were analyzed.The results revealed clear spatial and temporal differences.Spatially,the dominant fast-wave polarization direction at each station shows a strong correlation with the direction of the maximum horizontal principal compressive stress,as characterized by focal mechanism solutions of seismic events(MW≥3.5)near the station.The dominant fast-wave polarization direction and the regional stress field also showed a strong correlation with the intermovement of the Arabian Plate,African Plate,and Anatolian Block.Along the Nurdagi-Pazarcik fault zone,the seismic fault of event 1,stations closer to the middle of the fault where the mainshock occurred exhibited notably greater delay times than stations located towards the ends of the fault and far from the mainshock.In addition,the stations located to the east of the Nurdagi-Pazarcik fault and to the north of the Sürgüfault also exhibited large delay times.The spatial distribution of shear-wave splitting parameters obtained from each station indicates that the upper-crust anisotropy in the source area is mainly controlled by the regional stress field,which is closely related to the state of the block motion.During the seismogenic process of the MW7.7 earthquake,more stress accumulated in the middle of the Nurdagi-Pazarcik fault than at either end of the fault.Under the influence of the MW7.7 and MW7.6 events,the stress that accumulated during the seismogenic process of the earthquake doublet may have migrated towards some areas outside the aftershock intensive area after the earthquakes,and the crustal stress and its adjustment range near the outer stations increased significantly.With the exception of two stations with few effective events,all stations showed a consistent change in shear-wave splitting parameters over time.In particular,each station showed a decreasing trend in delay times after the doublet earthquakes,reflecting the obvious intensification of crustal stress adjustment in the seismogenic zone after the doublet earthquakes.With the occurrence of the earthquake doublet and a large number of aftershocks,the stress accumulated during the seismogenic process of the doublet earthquakes is gradually released,and then the adjustment range of crustal stress is also gradually reduced.

Analysis of faulting destruction and water supply pipeline damage from the first mainshock of the February 6,2023 Türkiye earthquake doublet

In 2023,two consecutive earthquakes exceeding a magnitude of 7 occurred in Türkiye,causing severe casualties and economic losses.The damage to critical urban infrastructure and building structures,including highways,railroads,and water supply pipelines,was particularly severe in areas where these structures intersected the seismogenic fault.Critical infrastructure projects that traverse active faults are susceptible to the influence of fault movement,pulse velocity,and ground motions.In this study,we used a unique approach to analyze the acceleration records obtained from the seismic station array(9 strong ground motion stations)located along the East Anatolian Fault(the seismogenic fault of the MW7.8 mainshock of the 2023 Türkiye earthquake doublet).The acceleration records were filtered and integrated to obtain the velocity and displacement time histories.We used the results of an on-site investigation,jointly conducted by China Earthquake Administration and Türkiye’s AFAD,to analyze the distribution of PGA,PGV,and PGD recorded by the strong motion array of the East Anatolian Fault.We found that the maximum horizontal PGA in this earthquake was 3.0 g,and the maximum co-seismic surface displacement caused by the East Anatolian Fault rupture was 6.50 m.As the fault rupture propagated southwest,the velocity pulse caused by the directional effect of the rupture increased gradually,with the maximum PGA reaching 162.3 cm/s.We also discussed the seismic safety of critical infrastructure projects traversing active faults,using two case studies of water supply pipelines in Türkiye that were damaged by earthquakes.We used a three-dimensional finite element model of the PE(polyethylene)water pipeline at the Islahiye State Hospital and fault displacement observations obtained through on-site investigation to analyze pipeline failure mechanisms.We further investigated the effect of the fault-crossing angle on seismic safety of a pipeline,based on our analysis and the failure performance of the large-diameter Thames Water pipeline during the 1999 Kocaeli earthquake.The seismic method of buried pipelines crossing the fault was summarized.

Determining the surface fault-rupture hazard zone for the Pazarcık segment of the East Anatolian fault zone through comprehensive analysis of surface rupture from the February 6,2023,Earthquake(Mw 7.7)

Following surface rupture observations in populated areas affected by the KahramanmaraşEarthquake(Mw 7.7)on February 6th,2023,along the Pazarcık segment of the East Anatolian Fault Zone(EAFZ),this study presents novel insights into physical criteria for delineating surface fault-rupture hazard zones(SRHZs)along ruptured strike-slip faults.To achieve this objective,three trench studies across the surface rupture were conducted on the Pazarcık segment of the EAFZ to collect field data,and earthquake recurrence intervals were interpreted using Bayesian statistics from previously conducted paleoseismological trenchings.The results of the proposed model indicate that the Pazarcık segment produced five significant surface-rupturing earthquakes in the last∼11 kyr:E1:11.13±1.74 kyr,E2:7.62±1.20 kyr,E3:5.34±1.05 kyr,E4:1.82±0.93 kyr,and E5:0.35±0.11 kyr.In addition,the recurrence intervals of destructive earthquakes on the subject in question range from 0.6 kyr to 4.8 kyr.Considering that the last significant earthquake occurred in 1513,the longest time since the most recent surface fault rupturing earthquake on this particular segment was 511 years.These results indicate that,in terms of the theoretical recurrence interval of earthquakes that can create surface ruptures on the Pazarcık segment,the period in which the February 6,2023,earthquake occurred was within the end of the expected return period.As a result,the potential for a devastating earthquake in the near future is not foreseen on the same fault.Finally,the SRHZ proposed for the Pazarcık section of Gölbaşıvillage was calculated as a 61-meter-wide offset on the fault lineament to reduce the negativities that may occur in the ruptured area in the future.It is recommended to take into account this width in the settlement of this area and nearby areas.

Cascading multi-segment rupture process of the 2023 Turkish earthquake doublet on a complex fault system revealed by teleseismic P wave back projection method

In this study,the vertical components of broadband teleseismic P wave data recorded by China Earthquake Network are used to image the rupture processes of the February 6th,2023 Turkish earthquake doublet via back projection analysis.Data in two frequency bands(0.5-2 Hz and 1-3 Hz)are used in the imaging processes.The results show that the rupture of the first event extends about 200 km to the northeast and about 150 km to the southwest,lasting~90 s in total.The southwestern rupture is triggered by the northeastern rupture,demonstrating a sequential bidirectional unilateral rupture pattern.The rupture of the second event extends approximately 80 km in both northeast and west directions,lasting~35 s in total and demonstrates a typical bilateral rupture feature.The cascading ruptures on both sides also reflect the occurrence of selective rupture behaviors on bifurcated faults.In addition,we observe super-shear ruptures on certain fault sections with relatively straight fault structures and sparse aftershocks.

Co-seismic displacements of 2011 Japan Mw9.0 earthquake recorded by far-field GPS stations

Co-seismic displacements of the 2011 Mw9.0 Japan earthquake recorded by GPS stations in China and surrounding areas showed a movement toward the epicenter. The horizontal displacements were up to 1 - 3 cm in northeastern China, 3 -8 mm in the North China, and 2 cm in the Korean peninsula. The vertical movements in China were small uplifts.

Far-field deformations caused by the 2004 Sumatra earthquake

Theoretical horizontal displacements caused by the 2004 Sumatra earthquake in the Sichuan-Yunnan area have been calculated according to a spherical dislocation theory and an earthquake-fault model. The results show that the theoretical displacements are basically consistent with the GPS observations in situ. On this basis,we have calculated the co-seismic displacements, strains, changes of gravity and geoid of the whole Earth, including China mainland and vicinity, caused by this earthquake. Key wards：