Regional wave propagation using the discontinuous Galerkin method

We present an application of the discontinuous Galerkin (DG) method to regional wave propagation. The method makes use of unstructured tetrahedral meshes, combined with a time integration scheme solving the arbitrary high-order derivative (ADER) Riemann problem. This ADER-DG method is high-order accurate in space and time, beneficial for reliable simulations of high-frequency wavefields over long propagation distances. Due to the ease with which tetrahedral grids can be adapted to complex geometries, undulating topography of the Earth's surface and interior interfaces can be readily implemented in the computational domain. The ADER-DG method is benchmarked for the accurate radiation of elastic waves excited by an explosive and a shear dislocation source. We compare real data measurements with synthetics of the 2009 L'Aquila event (central Italy). We take advantage of the geometrical flexibility of the approach to generate a European model composed of the 3-D <i>EPcrust</i> model, combined with the depth-dependent <i>ak135</i> velocity model in the upper mantle. The results confirm the applicability of the ADER-DG method for regional scale earthquake simulations, which provides an alternative to existing methodologies

Verification of an ADER-DG method for complex dynamic rupture problems

We present results of thorough benchmarking of an arbitrary high-order derivative discontinuous Galerkin (ADER-DG) method on unstructured meshes for advanced earthquake dynamic rupture problems. We verify the method by comparison to well-established numerical methods in a series of verification exercises, including dipping and branching fault geometries, heterogeneous initial conditions, bimaterial interfaces and several rate-and-state friction laws. We show that the combination of meshing flexibility and high-order accuracy of the ADER-DG method makes it a competitive tool to study earthquake dynamics in geometrically complicated setups

Tilt effects on moment tensor inversion in the near field of active volcanoes

Dynamic tilts (rotational motion around horizontal axes) change the projection of local gravity onto the horizontal components of seismometers. This causes sensitivity of these components to tilt, especially at low frequencies. We analyse the consequences of this effect onto moment tensor inversion for very long period (vlp) events in the near field of active volcanoes on the basis of synthetic examples using the station distribution of a real deployed seismic network and the topography of Mt. Merapi volcano (Java, Indonesia). The examples show that for periods in the vlp range of 10-30 s tilt can have a strong effect on the moment tensor inversion, although its effect on the horizontal seismograms is significant only for few stations. We show that tilts can be accurately computed using the spectral element method and include them in the Green's functions. The (simulated) tilts might be largely influenced by strain-tilt coupling (stc). However, due to the frequency dependence of the tilt contribution to the horizontal seismograms, only the largest tilt signals affect the source inversion in the vlp frequency range. As these are less sensitive to stc than the weaker signals, the effect of stc can likely be neglected in this application. In the converse argument, this is not necessarily true for longer periods, where the horizontal seismograms are dominated by the tilt signal and rotational sensors would be necessary to account for it. As these are not yet commercially available, this study underlines the necessity for the development of such instrument

On the initiation of sustained slip-weakening ruptures by localized stresses

Numerical simulations of dynamic earthquake rupture require an artificial initiation procedure, if they are not integrated in long-term earthquake cycle simulations. A widely applied procedure involves an ‘overstressed asperity’, a localized region stressed beyond the static frictional strength. The physical properties of the asperity (size, shape and overstress) may significantly impact rupture propagation. In particular, to induce a sustained rupture the asperity size needs to exceed a critical value. Although criteria for estimating the critical nucleation size under linear slip-weakening friction have been proposed for 2-D and 3-D problems based on simplifying assumptions, they do not provide general rules for designing 3-D numerical simulations. We conduct a parametric study to estimate parameters of the asperity that minimize numerical artefacts (e.g. changes of rupture shape and speed, artificial supershear transition, higher slip-rate amplitudes). We examine the critical size of square, circular and elliptical asperities as a function of asperity overstress and background (off-asperity) stress. For a given overstress, we find that asperity area controls rupture initiation while asperity shape is of lesser importance. The critical area obtained from our numerical results contrasts with published theoretical estimates when background stress is low. Therefore, we derive two new theoretical estimates of the critical size under low background stress while also accounting for overstress. Our numerical results suggest that setting the asperity overstress and area close to their critical values eliminates strong numerical artefacts even when the overstress is large. We also find that properly chosen asperity size or overstress may significantly shorten the duration of the initiation. Overall, our results provide guidelines for determining the size of the asperity and overstress to minimize the effects of the forced initiation on the subsequent spontaneous rupture propagation

Coupled, Physics-Based Modeling Reveals Earthquake Displacements are Critical to the 2018 Palu, Sulawesi Tsunami

The September 2018, Mw 7.5 Sulawesi earthquake occurring on the Palu-Koro strike-slip fault system was followed by an unexpected localized tsunami. We show that direct earthquake-induced uplift and subsidence could have sourced the observed tsunami within Palu Bay. To this end, we use a physics-based, coupled earthquake–tsunami modeling framework tightly constrained by observations. The model combines rupture dynamics, seismic wave propagation, tsunami propagation and inundation. The earthquake scenario, featuring sustained supershear rupture propagation, matches key observed earthquake characteristics, including the moment magnitude, rupture duration, fault plane solution, teleseismic waveforms and inferred horizontal ground displacements. The remote stress regime reflecting regional transtension applied in the model produces a combination of up to 6 m left-lateral slip and up to 2 m normal slip on the straight fault segment dipping 65∘ East beneath Palu Bay. The time-dependent, 3D seafloor displacements are translated into bathymetry perturbations with a mean vertical offset of 1.5 m across the submarine fault segment. This sources a tsunami with wave amplitudes and periods that match those measured at the Pantoloan wave gauge and inundation that reproduces observations from field surveys. We conclude that a source related to earthquake displacements is probable and that landsliding may not have been the primary source of the tsunami. These results have important implications for submarine strike-slip fault systems worldwide. Physics-based modeling offers rapid response specifically in tectonic settings that are currently underrepresented in operational tsunami hazard assessment