SUMMARY
We introduce a scheme for probabilistic hypocentre inversion with Stein variational inference. Our approach uses a differentiable forward model in the form of a physics informed neural ...network, which we train to solve the Eikonal equation. This allows for rapid approximation of the posterior by iteratively optimizing a collection of particles against a kernelized Stein discrepancy. We show that the method is well-equipped to handle highly multimodal posterior distributions, which are common in hypocentral inverse problems. A suite of experiments is performed to examine the influence of the various hyperparameters. Once trained, the method is valid for any seismic network geometry within the study area without the need to build traveltime tables. We show that the computational demands scale efficiently with the number of differential times, making it ideal for large-N sensing technologies like Distributed Acoustic Sensing. The techniques outlined in this manuscript have considerable implications beyond just ray tracing procedures, with the work flow applicable to other fields with computationally expensive inversion procedures such as full waveform inversion.
Our understanding of subsurface processes suffers from a profound observation bias: seismometers are sparse and clustered on continents. A new seismic recording approach, distributed acoustic sensing ...(DAS), transforms telecommunication fiber‐optic cables into sensor arrays enabling meter‐scale recording over tens of kilometers of linear fiber length. We analyze cataloged earthquake observations from three DAS arrays with different horizontal geometries to demonstrate some possibilities using this technology. In Fairbanks, Alaska, we find that stacking ground motion records along 20 m of fiber yield a waveform that shows a high degree of correlation in amplitude and phase with a colocated inertial seismometer record at 0.8–1.6 Hz. Using an L‐shaped DAS array in Northern California, we record the nearly vertically incident arrival of an earthquake from The Geysers Geothermal Field and estimate its backazimuth and slowness via beamforming for different phases of the seismic wavefield. Lastly, we install a fiber in existing telecommunications conduits below Stanford University and show that little cable‐to‐soil coupling is required for teleseismic P and S phase arrival detection.
Plain Language Summary
A new seismic recording technology, called distributed acoustic sensing (DAS), turns common telecommunications fiber‐optic cables into a very long (tens of kilometers) array of single‐component seismometers. We catalog and analyze earthquake observations from three different DAS experiments and demonstrate some of the possibilities of using this technology for earthquake seismology. We compare DAS records with a conventional seismometer and then show how the array style response of DAS enables recording of additional information such as the direction of the seismic energy. We also install a fiber in a typical telecommunications conduit below Stanford University, as opposed to directly burying the fiber in soil, and show that main seismic waves are still detected despite the hypothesized decrease in coupling between the fiber and the ground.
Key Points
Distributed acoustic sensing (DAS) transforms fiber optics into dense seismic arrays (1 sensor/meter) using laser scattering
Earthquake observations indicate that DAS and seismometer sensitivity is similar at 0.8‐1.6 Hz
Fiber installed in telecommunication conduit shows sensitivity to propagating seismic waves
SUMMARY Meteorite impacts have proved to be a significant source of seismic signal on the Moon, and have now been recorded on Mars by InSight seismometers. Understanding how impacts produce seismic ...signal is key to the interpretation of this unique data, and to improve their identification in continuous seismic records. Here, we use the seismic Representation Theorem, and particularly the stress glut theory, to model the seismic motion resulting from impact cratering. The source is described by equivalent forces, some resulting from the impactor momentum transfer, and others from the stress glut, which represents the mechanical effect of plasticity and non linear processes in the source region. We condense these equivalent forces into a point-source with a time-varying single force and nine-component moment tensor. This analytical representation bridges the gap between the complex dynamics of crater formation, and the linear point-source representation classically used in seismology. Using the multiphysics modelling software HOSS, we develop a method to compute the stress glut of an impact, and the associated point-source from hypervelocity impact simulations. For a vertical and an oblique impact at 1000 m s−1, we show that the moment tensor presents a significant deviatoric component. Hence, the source is not an ideal isotropic explosion contrary to previous assumptions, and draws closer to a double couple for the oblique impact. The contribution of the point force to the seismic signal appears negligible. We verify this model by comparing two signals: (1) HOSS is coupled to SPECFEM3D to propagate the near-source signal elastically to remote seismic stations; (2) the point-source model derived from the stress-glut theory is used to generate displacements at the same distance. The comparison shows that the point-source model is accurately simulating the low-frequency impact seismic waveform, and its seismic moment is in trend with Lunar and Martian impact data. High-frequencies discrepancies exist, which are partly related to finite-source effects, but might be further explained by the difference in mathematical framework between classical seismology and HOSS’ numerical modelling.
Recent seismological studies challenge the traditional view that the interface between the core and mantle is a straightforward discontinuity. As seismology is pushed to its observational limits, a ...complex - potentially compositionally layered - region between the core and mantle is emerging.
On 6 February 2023, the Msub.W 7.8 Pazarcik and the Msub.W 7.5 Elbistan earthquakes occurred in southeastern Turkey, close to the Syrian border, causing many deaths and a great deal of property ...destruction. The Pazarcik earthquake mainly damaged the East Anatolian Fault Zone (EAFZ). The Elbistan earthquake mainly damaged the Cardak fault (CF) and the Doğanşehir fault (DF). In this study, Sentinel-1A ascending (ASC) and descending (DES) orbit image data and pixel offset tracking (POT) were used to derive surface deformation fields in the range and azimuth directions induced by the Pazarcik and Elbistan earthquakes (hereinafter referred to as the Turkey double earthquakes). Utilizing GPS coordinate sequence data, we computed the three-dimensional surface deformation resulting from the Turkey double earthquakes. The surface deformation InSAR and GPS results were combined to invert the coseismic slip distribution of the EAFZ, CF, and DF using a layered earth model. The results show that the coseismic ruptures of the Turkey double earthquakes were dominated by left-lateral strike-slips. The maximum slip was 7.76 m on the EAFZ and about 8.2 m on the CF. Both the earthquakes ruptured the surface. The Coulomb failure stress (CFS) was computed based on the fault slip distribution and the geometric parameters of all the active faults within 300 km of the Msub.W 7.8 Pazarcik earthquake’s epicenter. The CFS change resulting from the Pazarcik earthquake suggests that the subsequent Elbistan earthquake was triggered by the Pazarcik earthquake. The Antakya fault experienced an increase in CFS of 8.4 bars during this double-earthquake event. Therefore, the Msub.W 6.3 Uzunbağ earthquake on 20 February 2023 was jointly influenced by the Turkey double earthquakes. Through stress analysis of all the active faults within 300 km of the Msub.W 7.8 Pazarcik earthquake’s epicenter, the Ecemis segment, Camliyayla fault, Aadag fault, Ayvali fault, and Pula segment were all found to be under stress loading. Particularly, the Ayvali fault and Pula segment exhibited conspicuous stress loading, signaling a higher risk of future seismic activity.
SUMMARY
We present results on radiated seismic energy during simulations of dynamic ruptures in a continuum damage-breakage rheological model incorporating evolution of damage within the seismic ...source region. The simulations vary in their initial damage zone width and rate of damage diffusion with parameter values constrained by observational data. The radiated energy recorded at various positions around the source is used to calculate seismic potency and moment. We also calculate the normalized radiated energy from the source, in a way that allows comparing between results of different simulations and highlighting aspects related to the dilatational motion during rupture. The results show that at high-frequencies, beyond the dominant frequency of the source ($( {f > 3{f}_d} )$, the damage process produces an additional burst of energy mainly in the Pwaves. This eccess of high-frequency energy is observed by comparing the radiated energy to a standard Brune's model with a decay slope of the radiated energy of n = 2. While the Swaves show good agreement with the n = 2 slope, the Pwaves have a milder slope of n = 1.75 or less depending on the damage evolution at the source. In the used damage-breakage rheology, the rate of damage diffusivity governs the damage evolution perpendicular to the rupture direction and dynamic changes of the damage zone width. For increasing values of damage diffusivity, dilatational energy becomes more prominent during rupture, producing a high-frequency dilatational signature within the radiation pattern. The high-frequency radiation pattern of the Pwaves includes two main lobes perpendicular to the rupture direction, reflecting high-rate local tensile cracking during the overall shear rupture process. Analysing the possible existence and properties of such high-frequency radiation pattern in observed Pwaves could provide important information on earthquake source processes.