Eric M. Dunham | Search USGS

USGS2019

Combining dynamic rupture simulations with ground motion data to characterize seismic hazard from Mw 3-5.8 earthquakes in Oklahoma and Kansas

Many seismically active areas suffer from a lack of near‐source ground‐motion recordings, making ground‐motion prediction difficult at distances within ∼ 40    km "> ∼ 40 km ∼40 km from an earthquake. We aim to aid the development of near‐source ground‐motion prediction equations (GMPEs) by generating synthetic ground‐motion data via simulation. Building on previous work using point‐source moment tensor sources to simulate small ( ⁠ M w "> M w Mw 3–4) earthquakes for a target region encompassing north central Oklahoma and south central Kansas, we perform dynamic rupture simulations of earthquakes up to M w "> M w Mw 5.8. We introduce complexity into the rupture process by adding stochastically generated heterogeneity to initial stress conditions. Our simulated ground‐motion data are added to a catalog of recorded ground‐motion data to construct a composite recorded‐simulated ground‐motion catalog that we use to develop a GMPE for the target region. This procedure can be generalized and used for GMPE development in other regions with near‐source ground‐motion data scarcity, which could directly benefit critical applications such as the National Seismic Hazard Maps produced by the U.S. Geological Survey (USGS).

Kansas, Oklahoma↗

USGS2013

Frictional-faulting model for harmonic tremor before Redoubt Volcano eruptions

Seismic unrest, indicative of subsurface magma transport and pressure changes within fluid-filled cracks and conduits, often precedes volcanic eruptions. An intriguing form of volcano seismicity is harmonic tremor, that is, sustained vibrations in the range of 0.5–5 Hz. Many source processes can generate harmonic tremor. Harmonic tremor in the 2009 eruption of Redoubt Volcano, Alaska, has been linked to repeating earthquakes of magnitudes around 0.5–1.5 that occur a few kilometres beneath the vent. Before many explosions in that eruption, these small earthquakes occurred in such rapid succession—up to 30 events per second—that distinct seismic wave arrivals blurred into continuous, high-frequency tremor. Tremor abruptly ceased about 30 s before the explosions. Here we introduce a frictional-faulting model to evaluate the credibility and implications of this tremor mechanism. We find that the fault stressing rates rise to values ten orders of magnitude higher than in typical tectonic settings. At that point, inertial effects stabilize fault sliding and the earthquakes cease. Our model of the Redoubt Volcano observations implies that the onset of volcanic explosions is preceded by active deformation and extreme stressing within a localized region of the volcano conduit, at a depth of several kilometres.

Alaska↗

USGS2010

Coherence of Mach fronts during heterogeneous supershear earthquake rupture propagation: Simulations and comparison with observations

We study how heterogeneous rupture propagation affects the coherence of shear and Rayleigh Mach wavefronts 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 owing 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 (2008a). 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 (2008a): (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.

Journal of Geophysical Research B: Solid Earth↗