Global S-wave tomography using receiver pairs: An alternative to get rid of earthquake mislocation

International audienceGlobal seismic tomography suffers from uncertainties in earthquake parameters routinely published in seismic catalogues. In particular, errors in earthquake location and origin-time may lead to strong biases in measured body wave delay-times and significantly pollute tomographic models. Common ways of dealing with this issue are to incorporate source parameters as additional unknowns into the linear tomographic equations, or to seek combinations of data to minimize the influence of source mislocations. We propose an alternative, physically-based method to desensitize direct S-wave delay-times to errors in earthquake location and origin-time. Our approach takes advantage of the fact that mislocation delay-time biases depend to first order on the earthquake-receiver azimuth, and to second order on the epicentral distance. Therefore, for every earthquake, we compute S-wave differential delay-times between optimized receiver pairs, such that a large part of their mislocation delay-time biases cancels out (for example origin-time fully subtracts out), while the difference of their sensitivity kernels remains sensitive to the model parameters of interest. Considering realistic, randomly distributed source mislocation vectors, as well as various levels of data noise and different synthetic Earths, we demonstrate that mislocation-related model errors are highly reduced when inverting for such differential delay-times, compared to absolute ones. The reduction is particularly rewarding for imaging the upper-mantle and transition zone. We conclude that using optimized receiver pairs is a suitable, low cost alternative to get rid of errors on earthquake location and origin-time for teleseismic direct S-wave traveltimes. Moreover, it can partly remove unilateral rupture propagation effects in cross-correlation delay-times, since they are similar to mislocation effects

The 2015 Gorkha earthquake: A large event illuminating the Main Himalayan Thrust fault

International audienceThe 2015 Gorkha earthquake sequence provides an outstanding opportunity to better characterize the geometry of the Main Himalayan Thrust (MHT). To overcome limitations due to unaccounted lateral heterogeneities, we perform Centroid Moment Tensor inversions in a 3-D Earth model for the main shock and largest aftershocks. In parallel, we recompute S-toP and P-to-S receiver functions from the Hi-CLIMB data set. Inverted centroid locations fall within a low-velocity zone at 10–15 km depth and corresponding to the subhorizontal portion of the MHT that ruptured during the Gorkha earthquake. North of the main shock hypocenter, receiver functions indicate a north dipping feature that likely corresponds to the midcrustal ramp connecting the flat portion to the deep part of the MHT. Our analysis of the main shock indicates that long-period energy emanated updip of high-frequency radiation sources previously inferred. This frequency-dependent rupture process might be explained by different factors such as fault geometry and the presence of fluids

Earthquake in a Maze: Compressional Rupture Branching During the 2012 M_w 8.6 Sumatra Earthquake

Seismological observations of the 2012 moment magnitude 8.6 Sumatra earthquake reveal unprecedented complexity of dynamic rupture. The surprisingly large magnitude results from the combination of deep extent, high stress drop, and rupture of multiple faults. Back-projection source imaging indicates that the rupture occurred on distinct planes in an orthogonal conjugate fault system, with relatively slow rupture speed. The east-southeast–west-northwest ruptures add a new dimension to the seismotectonics of the Wharton Basin, which was previously thought to be controlled by north-south strike-slip faulting. The rupture turned twice into the compressive quadrant, against the preferred branching direction predicted by dynamic Coulomb stress calculations. Orthogonal faulting and compressional branching indicate that rupture was controlled by a pressure-insensitive strength of the deep oceanic lithosphere

Aseismic slip and seismogenic coupling along the central San Andreas Fault

International audienceWe use high-resolution Synthetic Aperture Radar- and GPS-derived observations of surfacedisplacements to derive the first probabilistic estimates of fault coupling along the creeping section of theSan Andreas Fault, in between the terminations of the 1857 and 1906 magnitude 7.9 earthquakes. Usinga fully Bayesian approach enables unequaled resolution and allows us to infer a high probability ofsignificant fault locking along the creeping section. The inferred discreet locked asperities are consistentwith evidence for magnitude 6+ earthquakes over the past century in this area and may be associated withthe initiation phase of the 1857 earthquake. As creeping segments may be related to the initiation andtermination of seismic ruptures, such distribution of locked and creeping asperities highlights the centralrole of the creeping section on the occurrence of major earthquakes along the San Andreas Fault

Transient self-potential anomalies associated with recent lava flows at Piton de la Fournaise volcano (Réunion Island, Indian Ocean)

International audienceSelf-potential signals are sensitive to various phenomena including ground water flow (streaming potential), thermal gradients (thermoelectric potential), and potentially rapid fluid disruption associated with vaporization of water. We describe transient self-potential anomalies observed over recent (< 9 years) lava flows at Piton de la Fournaise volcano (Reunion Island, Indian Ocean). Repeated self-potential measurements are used to determine the decay of the self-potential signals with time since the emplacement of a set of lava flow. We performed a 9 km-long self-potential profile in February 2004 in the Grand Brûlé area. This profile was repeated in July–August 2006. The second repetition of this profile crossed eight lava flows emplaced between 1998 and 2005 during seven eruptions of Piton de la Fournaise volcano. The self-potential data show clear positive anomalies (up to 330 mV) and spatially correlated with the presence of recent lava flows. The amplitude of the self-potential anomalies decreases exponentially with the age of the lava flows with a relaxation time of not, vert, similar 44 months. We explain these anomalies by the shallow convection of meteoric water and the associated streaming potential distribution but we cannot exclude possible contributions from the thermoelectric effect and the rapid fluid disruption mechanism. This field case evidences for the first time transient self-potential signals associated with recent volcanic deposits. It can be also a shallow analogue to understand the variation of self-potential signals in active geothermal areas and transient self-potential signals associated with dike intrusion at larger depths. The empirical equation we proposed can also be used to diagnose the cooling of recent lava flow on shield volcanoes

The 2013 M_w 7.7 Balochistan Earthquake: Seismic Potential of an Accretionary Wedge

Great earthquakes rarely occur within active accretionary prisms, despite the intense long‐term deformation associated with the formation of these geologic structures. This paucity of earthquakes is often attributed to partitioning of deformation across multiple structures as well as aseismic deformation within and at the base of the prism (Davis et al., 1983). We use teleseismic data and satellite optical and radar imaging of the 2013 M_w 7.7 earthquake that occurred on the southeastern edge of the Makran plate boundary zone to study this unexpected earthquake. We first compute a multiple point‐source solution from W‐phase waveforms to estimate fault geometry and rupture duration and timing. We then derive the distribution of subsurface fault slip from geodetic coseismic offsets. We sample for the slip posterior probability density function using a Bayesian approach, including a full description of the data covariance and accounting for errors in the elastic structure of the crust. The rupture nucleated on a subvertical segment, branching out of the Chaman fault system, and grew into a major earthquake along a 50° north‐dipping thrust fault with significant along‐strike curvature. Fault slip propagated at an average speed of 3.0  km/s for about 180 km and is concentrated in the top 10 km with no displacement on the underlying décollement. This earthquake does not exhibit significant slip deficit near the surface, nor is there significant segmentation of the rupture. We propose that complex interaction between the subduction accommodating the Arabia–Eurasia convergence to the south and the Ornach Nal fault plate boundary between India and Eurasia resulted in the significant strain gradient observed prior to this earthquake. Convergence in this region is accommodated both along the subduction megathrust and as internal deformation of the accretionary wedge

Strain budget of the Ecuador–Colombia subduction zone: A stochastic view

The 2016 Pedernales earthquake (M_W=7.8) ruptured a portion of the Colombia–Ecuador subduction interface where several large historical earthquakes have been documented since the great 1906 earthquake (M=8.6). Considering all significant ruptures that occurred in the region, it has been suggested that the cumulative moment generated co-seismically along this part of the subduction over the last century exceeds the moment deficit accumulated inter-seismically since 1906. Such an excess challenges simple models with earthquakes resetting the elastic strain accumulated inter-seismically in locked asperities. These inferences are however associated with large uncertainties that are generally unknown. The impact of spatial smoothing constraints on co-seismic and inter-seismic models also prevents any robust assessment of the strain budget. We propose a Bayesian kinematic slip model of the 2016 Pedernales earthquake using the most comprehensive dataset to date including InSAR and GPS offsets, tsunami waveforms, and kinematic records from high-rate GPS and strong-motions. In addition, we use inter-seismic geodetic velocities to produce a probabilistic inter-seismic coupling model of the subduction interface. Our stochastic co-seismic and inter-seismic solutions include the ensemble of all plausible models consistent with our prior information and that fit the observations within uncertainties. The analysis of these model ensembles indicates that an excess of co-seismic moment during the 1906–2016 period is likely in Central Ecuador only if we assume that 1942 and 2016 earthquakes are colocated. If this assumption is relaxed, we show that this conclusion no longer holds given uncertainties in co- and inter-seismic processes. The comparison of 1942 and 2016 teleseismic records reveals large uncertainties in the location of the 1942 event, hampering our ability to draw strong conclusions on the unbalanced moment budget in the region. Our results also show a heterogeneous coupling of the subduction interface that coincides with two slip asperities in our co-seismic model for the 2016 Pedernales earthquake and with the location of historical ruptures in 1958, 1979 and 1998. The spatial variability in coupling and complexity in earthquake history suggest strong heterogeneities in frictional properties of the subduction megathrust

Citizen seismology helps decipher the 2021 Haiti earthquake

5 pages, 4 figures, supplementary materials https://doi.org/10.1126/science.abn1045.-- Data and materials availability: All data and code used in this study are openly available. RADAR data can be obtained through ESA (Sentinel) or JAXA (Alos-2). Aftershock data can be obtained from https://ayiti.unice.fr/ayiti-seismes/ (7). The codes used to process or model the data are published and public (8). The catalog of high-precision earthquake relocated with the NLL-SSST-coherence procedure (SM4) is available as supplementary dataOn 14 August 2021, the moment magnitude (Mw) 7.2 Nippes earthquake in Haiti occurred within the same fault zone as its devastating 2010 Mw 7.0 predecessor, but struck the country when field access was limited by insecurity and conventional seismometers from the national network were inoperative. A network of citizen seismometers installed in 2019 provided near-field data critical to rapidly understand the mechanism of the mainshock and monitor its aftershock sequence. Their real-time data defined two aftershock clusters that coincide with two areas of coseismic slip derived from inversions of conventional seismological and geodetic data. Machine learning applied to data from the citizen seismometer closest to the mainshock allows us to forecast aftershocks as accurately as with the network-derived catalog. This shows the utility of citizen science contributing to our understanding of a major earthquakeThis work was supported by the Centre National de la Recherche Scientifique (CNRS) and the Institut de Recherche pour le Développement (IRD) through their “Natural Hazard” program (E.C., S.S., T.M., B.D., F.C., J.P.A., J.C., A.D., D.B., S.P.); the FEDER European Community program within the Interreg Caraïbes “PREST” project (E.C., S.S., D.B.); Institut Universitaire de France (E.C., R.J.); Université Côte d’Azur and the French Embassy in Haiti (S.P.); the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (grant no. 758210, Geo4D project to R.J. and grant no. 805256 to Z.D.); the French National Research Agency (project ANR-21-CE03-0010 “OSMOSE” to E.C. and ANR-15-IDEX-01 “UCAJEDI Investments in the Future” to Q.B.); the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (grant no. 949221 to Q.B.); and HPC resources of IDRIS (under allocations 2020-AD011012142, 2021-AP011012536, and 2021-A0101012314 to Q.B.With the institutional support of the ‘Severo Ochoa Centre of Excellence’ accreditation (CEX2019-000928-S)Peer reviewe