Earthquake source parameter estimation using synthetic waveform modelling

Fault plane solutions and focal depths for three crustal events occurring in the Himalayan collision zone have been obtained using synthetic waveform modelling. Two crustal events with their epicenters in the Tibetan plateau show large component of normal faulting with east-west trading T-axes. The third event with It’s epicenter north of Main Boundary Thrust (MBT) shows reverse faulting with the nodal planes paralleling the local structural trend. All the three crustal events studied have occurred at shallow focal depths of less than 15 km. The Inferred source parameters of these events are discussed In the light of active tectonics of the region.

Structural controls on earthquake rupture revealed by the 2020 Mw 6.0 Jiashi earthquake (Kepingtag belt, SW Tian Shan, China)

The Kepingtag (Kalpin) fold-and-thrust belt of the southern Chinese Tian Shan is characterized by active shortening and intense seismic activity. Geological cross-sections and seismic reflection profiles suggest thin-skinned, northward-dipping thrust sheets detached in an Upper Cambrian décollement. The January 19 2020 Mw 6.0 Jiashi earthquake provides an opportunity to investigate how coseismic deformation is accommodated in this structural setting. Coseismic surface deformation resolved with Sentinel-1 Interferometric Synthetic Aperture Radar (InSAR) is centered on the back limb of the frontal Kepingtag anticline. Elastic dislocation modelling suggests that the causative fault is located at ~7 km depth and dips ~7° northward, consistent with the inferred position of the décollement. The narrow slip pattern (length ~37 km but width only ~9 km) implies that there is a strong structural or lithological control on the rupture extent, with up-dip slip propagation possibly halted by an abrupt change in dip angle where the Kepingtag thrust is inferred to branch off the décollement. A depth discrepancy between mainshock slip constrained by InSAR and teleseismic waveform modelling (~7 km) and well-relocated aftershocks (~10-20 km) may imply that sediments above the décollement are velocity strengthening. We also relocate 148 regional events from 1977 to 2020 to characterize the broader distribution of seismicity across the Kepingtag belt. The calibrated hypocenters combined with previous teleseismic waveform models show that thrust and reverse faulting earthquakes cluster at relatively shallow depths of ~7-15 km but include abundant out-of-sequence events both north and south of the frontal Kepingtag fault.

Aftershock signature of the M7.5 Palu 2018 supershear rupture from a rapidly deployed nodal array

<p>Supershear earthquakes have significant implications for seismic hazard, in terms of  ground shaking and aftershock pattern. It has been suggested that supershear ruptures are associated with fewer aftershocks on the supershear rupture segment, however this needs to be tested using high resolution event locations. Current aftershock catalogues for the M7.5 Palu 2018 supershear rupture are not of sufficient resolution to identify any characteristic aftershock pattern. Additionally it is unclear whether the supershear rupture speed occurred from the time of earthquake initiation, or at a later time on a certain segment of the fault.</p><p>We deployed a nodal array to record aftershocks following the main event. The array comprised of twenty short-period nodes, which can be deployed rapidly, making them ideal for post-rupture investigations in areas of sparse coverage. We expand the earthquake catalogue by applying template matching to the nodal array data. We then relocate seismicity recorded by the array using a double difference method. We also relocate seismicity that occurred before the array was active, using a relative relocation method. To do this, we calibrate the more distant permanent stations using events well-located by the nodal array. We further derive moment tensors for the largest events by waveform modelling using short-period and broadband records.</p><p>Our results show that the aftershocks cluster at the northern and southern extents of rupture. There is a relative dearth of aftershocks in the middle part of the rupture, particularly in the Palu valley, where rupture terminated to the surface. The fault here is a long and straight distinctive geomorphic feature. Many secondary faults were triggered, particularly in the southern Sapu valley fault system. An earthquake swarm was triggered 1 month after the main event on a strike-slip fault 200km away.</p>

The 25 October, 2018 Zakynthos (Greece) earthquake: seismic activity at the transition between a transform fault and a subduction zone.

Summary The properties of the Mw = 6.7 earthquake that took place on 25 October 2018, 22:54:51 UTC, ∼50 km SW of the Zakynthos Island, Greece, are thoroughly examined. The main rupture occurred on a dextral strike-slip, low-angle, east-dipping fault at a depth of 12 km, as determined by teleseismic waveform modelling. Over 4000 aftershocks were manually analysed for a period of 158 days. The events were initially located with an optimal 1D velocity model and then relocated with the double-difference method to reveal details of their spatial distribution. The latter spreads in an area spanning 80 km NNW-SSE and ∼55 km WSW-ENE. Certain parts of the aftershock zone present strong spatial clustering, mainly to the north, close to Zakynthos Island, and at the southernmost edge of the sequence. Focal mechanisms were determined for 61 significant aftershocks using regional waveform modelling. The results revealed characteristics similar to the mainshock, with few aftershocks exhibiting strike-slip faulting at steeper dip angles, possibly related to splay faults on the accretionary prism. The slip vectors that correspond to the east-dipping planes are compatible with the long-term plate convergence and with the direction of coseismic displacement on the Zakynthos Island. Fault-plane solutions in the broader study area were inverted for the determination of the regional stress-field. The results revealed a nearly horizontal, SW-NE to E-W-trending S1 and a more variable S3 axis, favouring transpressional tectonics. Spatial clusters at the northern and southern ends of the aftershock zone coincide with the SW extension of sub-vertical along-dip faults of the segmented subducting slab. The mainshock occurred in an area where strike-slip tectonics, related to the Cephalonia Transform Fault and the NW Peloponnese region, gradually converts into reverse faulting at the western edge of the Hellenic subduction. Plausible scenarios for the 2018 Zakynthos earthquake sequence include a rupture on the subduction interface, provided the slab is tilted eastwards in that area, or the reactivation of an older east-dipping thrust as a low-angle strike-slip fault that contributes to strain partitioning.

Monitoring the seasonal changes of an englacial conduit network using repeated ground-penetrating radar measurements

Englacial conduits act as water pathways to feed surface meltwater into the subglacial drainage system. A change of meltwater into the subglacial drainage system can alter the glacier's dynamics. Between 2012 and 2019, repeated 25 MHz ground-penetrating radar (GPR) surveys were carried out over an active englacial conduit network within the ablation area of the temperate Rhonegletscher, Switzerland. In 2012, 2016, and 2017 GPR measurements were carried out only once a year, and an englacial conduit was detected in 2017. In 2018 and 2019 the repetition survey rate was increased to monitor seasonal variations in the detected englacial conduit. The resulting GPR data were processed using an impedance inversion workflow to compute GPR reflection coefficients and layer impedances, which are indicative of the conduit's infill material. The spatial and temporal evolution of the reflection coefficients also provided insights into the morphology of the Rhonegletscher's englacial conduit network. During the summer melt seasons, we observed an active, water-filled, sediment-transporting englacial conduit network that yielded large negative GPR reflection coefficients (<-0.2). The GPR surveys conducted during the summer provided evidence that the englacial conduit was 15–20 m±6 m wide, ∼0.4m±0.35m thick, ∼250m±6m long with a shallow inclination (2∘), and having a sinusoidal shape from the GPR data. We speculate that extensional hydraulic fracturing is responsible for the formation of the conduit as a result of the conduit network geometry observed and from borehole observations. Synthetic GPR waveform modelling using a thin water-filled conduit showed that a conduit thickness larger than 0.4 m (0.3× minimum wavelength) thick can be correctly identified using 25 MHz GPR data. During the winter periods, the englacial conduit no longer transports water and either physically closed or became very thin (<0.1 m), thereby producing small negative reflection coefficients that are caused by either sediments lying within the closed conduit or water within the very thin conduit. Furthermore, the englacial conduit reactivated during the following melt season at an identical position as in the previous year.

Variational full-waveform inversion

SUMMARY Seismic full-waveform inversion (FWI) can produce high-resolution images of the Earth’s subsurface. Since full-waveform modelling is significantly nonlinear with respect to velocities, Monte Carlo methods have been used to assess image uncertainties. However, because of the high computational cost of Monte Carlo sampling methods, uncertainty assessment remains intractable for larger data sets and 3-D applications. In this study, we propose a new method called variational FWI, which uses Stein variational gradient descent to solve FWI problems. We apply the method to a 2-D synthetic example and demonstrate that the method produces accurate approximations to those obtained by Hamiltonian Monte Carlo. Since variational inference solves the problem using optimization, the method can be applied to larger data sets and 3-D applications by using stochastic optimization and distributed optimization.

Probing mantle plumes using seismic arrays

&lt;p&gt;Elucidating the role of deep mantle plumes in mantle convection is challenging because their influence on seismic waveforms &amp;#8211; which could be used to map their location &amp;#8211; is subtle. Previous seismic studies have mainly focused on waveform modelling and inversion (i.e. tomography). In this study we instead consider the potential visibility of mantle plumes using array methods. We investigate, in particular, how plumes deviate seismic energy from the great-circle path. This requires a multidisciplinary approach: first, we perform geodynamic modelling to generate thermochemical plumes, and convert them to &amp;#8220;seismic&amp;#8221; plumes via thermodynamic modelling of mineral physics data. Next, spectral element methods are used to model the interaction of seismic waves with the plumes and generate synthetic seismograms. These seismograms are divided into arrays and we generate slowness-backazimuth plots for each array. With recent advances in computational methods and resources, we investigate wave behaviour at previously unattainable frequencies.&amp;#160; We find that plumes do indeed cause seismic waves to change direction, although the exact behaviour may be frequency-dependent, and at low frequencies we observe waves apparently bending around the plume conduit.&amp;#160; We consider how and where these results may be applied to real seismic arrays, to provide new constraints on the location and structure of mantle plumes.&lt;/p&gt;