2000 Uglegorsk earthquake

The paper presents the results of waveform inversion of the Mw 6.8 August 4 (5), 2000 Uglegorsk earthquake (Sakhalin Island, Russia). The detailed rupture process of the 2000 Uglegorsk earthquake was simulated using the waveform inversion method. The average parameters were calculated for both nodal planes. Waveform inversion was carried out on the basis of Global Seismographic Network (GSN) data. Only P-waves from BHZ channels of all stations from the GSN were used. The simulated source parameters included a double-couple source, the scalar seismic moment, the source time function, and the slip directions. The performed studies made it possible to investigate the features of the rupture development and the amplitude of displacements along the east and west-dipping nodal planes of the August 4 (5), 2000 Uglegorsk earthquake. The obtained P-slip model for the 2000 Uglegorsk earthquake source area is in good agreement with the surface manifestations of the rupture according to the field geology data and the results of geodetic inversion.

Combined Geodetic and Seismological Study of the December 2020 Mw = 4.6 Thiva (Central Greece) Shallow Earthquake

On 2 December 2020, a moderate and shallow Mw = 4.6 earthquake occurred in Boeotia (Central Greece) near the city of Thiva. Despite its magnitude, the co-seismic ground deformation field was detectable and measurable by Sentinel-1, ascending and descending, synthetic aperture interferometry radar (InSAR) acquisitions. The closest available GNSS station to the epicenter, located 11 km west, measured no deformation, as expected. We proceeded to the inversion of the deformation source. Moreover, we reassessed seismological data to identify the activated zone, associated with the mainshock and the aftershock sequence. Additionally, we used the rupture plane information from InSAR to better determine the focal mechanism and the centroid location of the mainshock. We observed that the mainshock occurred at a shallower depth and the rupture then expanded downdip, as revealed by the aftershock distribution. Our geodetic inversion modelling indicated the activation of a normal fault with a small left-lateral component, length of 2.0 km, width of 1.7 km, average slip of 0.2 m, a low dip angle of 33°, and a SW dip-direction. The inferred fault top was buried at a depth of ~0.5 km, rooted at a depth of ~1.4 km, with its geodetic centroid buried at 1.0 km. It was aligned with the Kallithea fault. In addition, the dip-up projection of the modeled fault to the surface was located very close (~0.4 km SW) to the mapped (by existing geological observations) trace of the Kallithea fault. The ruptured area was settled in a transition zone. We suggest the installation of at least one GNSS and seismological station near Kallithea; as the activated zone (inferred by the aftershock sequence and InSAR results) could yield events with M≥5.0, according to empirical laws relating to rupture zone dimensions and earthquake magnitude.

Modeling surface deformation of geothermal environments with high-fidelity finite element models

<p>The geothermal environment is an assembly of heterogeneous geological settings and complex interactions among different phases of rock and fluid medium. The artificial activities of energy production could further interplay with the on-going natural processes within volcanos, hotspots, and other geothermal areas to result in spatiotemporal signatures of displacements near the surface. Here, we study the temporal ground deformation near a geothermal site through processing Interferometric Synthetic Aperture Radar (InSAR) time-series data obtained over the past decades, as reconciled with the nearby GPS station. To interpret these signals and potentially reveal the reservoir’s temporal activity, we employ state-of-art finite element models (FEMs) to simulate a more realistic crustal domain near the energy-production zone with irregular reservoir geometry, distributed rock materials, and surface topography. Linear Bayesian geodetic inversion and Green’s function library area were adopted to quantify the cause of surface subsidence, as compared to the documented production history. Our study demonstrates an unprecedented approach to precisely simulate the elastic deformation caused by geothermal energy extraction and pumping, providing an important platform to further explore the in-depth evolving stress state and its relation to surrounding induced and natural seismicity.</p>

Joint Inversion of Geodetic Observations and Relative Weighting—The 1999 Mw 7.6 Chi-Chi Earthquake Revisited

The Mw 7.6 Chi-Chi earthquake struck central western Taiwan in 1999. The rupture was complex with several dislocations along the 100-km long Chelungpu thrust fault. Revisiting this earthquake is a challenge, as the precision and coverage of the data sets available are quite poor. Furthermore, the topographic and vegetation coverage complexity of the area prevents coherent radar images. In this paper, radar and optical images, and terrestrial geodetic measurements, were utilised to study the fault. The Particle Swarm Optimization and Okada Inversion (PSOKINV) geodetic inversion package was used with the generalized Akaike’s Bayesian Information Criterion (gABIC) to precisely determine the slip distribution and relative weighting of datasets. Differences in results using the data sets jointly or separately (e.g., under-estimation due to InSAR, inconsistencies in SPOT offsets, smoother slip distribution with gABIC weighting) are observable. Most of the energy was released in the northern part of the fault, where the strike veers toward the east, and mainly at depths less than 4 km. The PSOKINV-gABIC approach is viable for the study of complicated cases such as the Chi-Chi earthquake and can significantly benefit the weight determination and physical realism of the fault geometry.

3-D joint geodetic and strong-motion finite fault inversion of the 2008 May 12, Wenchuan, China Earthquake

SUMMARY We present a source inversion of the 2008 Wenchuan, China earthquake, using strong-motion waveforms and geodetic offsets together with 3-D synthetic ground motions. We applied the linear multiple time window technique considering geodetic and dynamic Green's functions computed with the finite-element method and the reciprocity and Strain Green's Tensor formalism. All ground motion estimates, valid up to 1 Hz, accounted for 3-D effects, including the topography and the geometry of the Beichuan and Pengguan faults. Our joint inversion has a higher moment (M0) than a purely geodetic inversion and the slip distribution presents differences when compared to 1-D model source inversions. The moment is estimated to be M0 = 1.2 × 1021 N·m, slightly larger than other works. Our results show that considering a complex 3-D structure reduces the size of large areas of 10 m slip or greater by distributing it in wider zones, with reduced slips, in the central portion of the Beichuan and the Pengguan faults. Finally, we compare our source with a relocated aftershock catalogue and conclude that the 4–5 m slip contours approximately bound the absence or presence of aftershocks.

Geodetic inversion of complex fault geometries for strong earthquakes

<p>The fault geometry closely controls earthquake rupture process. Previous seismic inversion of the fault geometry is to derive the multiple-point moment tensor solutions. Because of the trade-off between the moment tensor and rupture velocity, the inversion has high instabilities. In contrast, geodetic inversion has less unknowns, since there is no need to solve for rupture velocity. But from the elastic dislocation theory, the relations between the surface deformation and sub-fault parameters (i.e. strike, dip and rake) are nonlinear. In this study, we develop a linear technique to invert geodetic data for sub-fault moment tensors. From the sub-fault moment tensor solutions, the strike, dip, rake, and their spatial variations can be estimated, which provide valuable information for assessing the complexities in fault geometry. We applied this technique to several significant earthquakes, i.e., the 2008 Mw7.9 Wenchuan earthquake, the 2015 Mw7.8 Gorkha earthquake, and the 2017 Mw6.5 Jiuzhaigou earthquake. The results of the 2008 Wenchuan earthquake suggest that the strike, dip and rake are all variable from southwest to northeast, which are well consistent with the aftershock distributions and mechanisms. The dip variations of the 2015 Gorkha earthquake suggest the earthquake has ruptured a listric fault (dep decreases with depth). Particularly, a dip anomaly appears in the northeast corner of the rupture area, indicating a geometric barrier accounting for the slip gap between the mainshock and largest Mw7.3 aftershock. For the 2017 Jiuzhaigou earthquake, two right-stepping and left-lateral strike-slip segments were distinguished. Accordingly, a compressional step-over was identified between the two segments.</p>

Multi-Sensor Geodetic Observations and Modeling of the 2017 Mw 6.3 Jinghe Earthquake

The Mw 6.3 Jinghe earthquake struck Xingjiang Province, China, on 8 August 2017 (05:15:04 UTC); the epicenter was near the Kusongmuxieke Piedmont Fault (KPF) of the northern Tian Shan Mountains. We used multi-source and multi-track satellite Synthetic Aperture Radar (SAR) imagery and Interferometric SAR (InSAR) techniques to reconstruct the coseismic displacement field from different line-of-sight geometries. To reduce the phase artifacts, we employed multi-temporal scenes acquired by Sentinel-1, and reconstructed the coseismic deformation through a temporal averaging strategy. Together with a single interferometric pair obtained using the Phased Array type L-band Synthetic Aperture Radar 2 (PALSAR2) sensor aboard the Advanced Land Observing Satellite 2 (ALOS2), we obtained five displacement maps with slightly different viewing geometries; all of which were used to constrain a geodetic inversion to retrieve the fault geometry parameters and slip distribution. Based on the focal mechanism and regional geology, we constructed multiple fault models that differ in dip direction (south and north dipping), and various striking angles. Both models fit the InSAR displacement maps, but have slip distributions of different depths. The slip depth of the south dipping model, with a dip of ~42°, is the most consistent with the relocated earthquake sequence and regional geological structure. Through the geodetic inversion, the maximum slip (0.25 m) occurred at 14.05 km and the associated rake was 89.56°. The result implies that the seismogenic fault is a blind thrust fault north of KPF (towards the foreland). Considering the relative locations of the suggested blind fault, the KPF, and the continuing north to south (N–S) shortening of the Tian Shan Mountains, this fault could be formed by the northward propagation of the regional fold-thrust belt.