ALBERT

Description: An edited version of this paper was published by AGU. Copyright (2010) American Geophysical Union.

Description: We study how heterogeneous rupture propagation affects the coherence of shear– and Rayleigh–Mach wave fronts radiated by supershear earthquakes. We address this question using numerical simulations of ruptures on a planar, vertical strike–slip fault embedded in a three–dimensional, homogeneous, linear elastic half–space. Ruptures propagate spontaneously in accordance with a linear slip–weakening friction law through both homogeneous and heterogeneous initial shear stress fields. In the 3–D homogeneous case, rupture fronts are curved due to interactions with the free surface and the finite fault width; however, this curvature does not greatly diminish the coherence of Mach fronts relative to cases in which the rupture front is constrained to be straight, as studied by Dunham and Bhat (2008). Introducing heterogeneity in the initial shear stress distribution causes ruptures to propagate at speeds that locally fluctuate above and below the shear–wave speed. Calculations of the Fourier amplitude spectra (FAS) of ground velocity time histories corroborate the kinematic results of Bizzarri and Spudich (2008): 1) The ground motion of a supershear rupture is richer in high frequency with respect to a subshear one. 2) When a Mach pulse is present, its high frequency content overwhelms that arising from stress heterogeneity. Present numerical experiments indicate that a Mach pulse causes approximately an –1.7 high frequency falloff in the FAS of ground displacement. Moreover, within the context of the employed representation of heterogeneities and over the range of parameter space that is accessible with current computational resources, our simulations suggest that while heterogeneities reduce peak ground velocity and diminish the coherence of the Mach fronts, ground motion at stations experiencing Mach pulses should be richer in high frequencies compared to stations without Mach pulses. In contrast to the foregoing theoretical results, we find no average elevation of 5%–damped absolute response spectral accelerations (SA) in the period band 0.05–0.4 s observed at stations that presumably experienced Mach pulses during the 1979 Imperial Valley, 1999 Kocaeli, and 2002 Denali Fault earthquakes compared to SA observed at non–Mach pulse stations in the same earthquakes. A 20% amplification of short period SA is seen only at a few of the Imperial Valley stations closest to the fault. This lack of elevated SA suggests that either Mach pulses in real earthquakes are even more incoherent that in our simulations, or that Mach pulses are vulnerable to attenuation through nonlinear soil response. In any case, this result might imply that current engineering models of high frequency earthquake ground motions do not need to be modified by more than 20% close to the fault to account for Mach pulses, provided that the existing data are adequately representative of ground motions from supershear earthquakes.

Description: Published

Description: B08301

Description: 3.1. Fisica dei terremoti

Description: JCR Journal

Description: reserved

Description: Tsunami generation by offshore earthquakes is a problem of scientific interest and practical relevance, and one that requires numerical modelling for data interpretation and hazard assessment. Most numerical models utilize two-step methods with one-way coupling between separate earthquake and tsunami models, based on approximations that might limit the applicability and accuracy of the resulting solution. In particular, standard methods focus exclusively on tsunami wave modelling, neglecting larger amplitude ocean acoustic and seismic waves that are superimposed on tsunami waves in the source region. In this study, we compare four earthquake-tsunami modelling methods. We identify dimensionless parameters to quantitatively approximate dominant wave modes in the earthquake-tsunami source region, highlighting how the method assumptions affect the results and discuss which methods are appropriate for various applications such as interpretation of data from offshore instruments in the source region. Most methods couple a 3-D solid earth model, which provides the seismic wavefield or at least the static elastic displacements, with a 2-D depth-averaged shallow water tsunami model. Assuming the ocean is incompressible and tsunami propagation is negligible over the earthquake duration leads to the instantaneous source method, which equates the static earthquake seafloor uplift with the initial tsunami sea surface height. For longer duration earthquakes, it is appropriate to follow the time-dependent source method, which uses time-dependent earthquake seafloor velocity as a forcing term in the tsunami mass balance. Neither method captures ocean acoustic or seismic waves, motivating more advanced methods that capture the full wavefield. The superposition method of Saito et al. solves the 3-D elastic and acoustic equations to model the seismic wavefield and response of a compressible ocean without gravity. Then, changes in sea surface height from the zero-gravity solution are used as a forcing term in a separate tsunami simulation, typically run with a shallow water solver. A superposition of the earthquake and tsunami solutions provides an approximation to the complete wavefield. This method is algorithmically a two-step method. The complete wavefield is captured in the fully coupled method, which utilizes a coupled solid Earth and compressible ocean model with gravity. The fully coupled method, recently incorporated into the 3-D open-source code SeisSol, simultaneously solves earthquake rupture, seismic waves and ocean response (including gravity). We show that the superposition method emerges as an approximation to the fully coupled method subject to often well-justified assumptions. Furthermore, using the fully coupled method, we examine how the source spectrum and ocean depth influence the expression of oceanic Rayleigh waves. Understanding the range of validity of each method, as well as its computational expense, facilitates the selection of modelling methods for the accurate assessment of earthquake and tsunami hazards and the interpretation of data from offshore instruments.