Geology topics

Don Helmberger

Publications and source records attributed to Don Helmberger.

6 recordsLinked to original sources

USGS2020

Shallow basin structure and attenuation are key to predicting long shaking duration in Los Angeles Basin

Ground motions in the Los Angeles Basin during large earthquakes are modulated by earthquake ruptures, path effects into the basin, basin effects, and local site response. We analyzed the direct effect of shallow basin structures on shaking duration at a period of 2–10 s in the Los Angeles region through modeling small magnitude, shallow, and deep earthquake pairs. The source depth modulates the basin response, particularly the shaking duration, and these features are a function of path effect and not site condition. Three-dimensional simulations using the CVM-S4.26.M01 velocity model show good fitting to the initial portion of the waveforms at periods of 5 s and longer but fail to predict the long shaking duration during shallow events, especially at periods less than 5 s. Simulations using CVM-H do not match the timing of the initial arrivals as well as CVM-S4.26.M01, and the strong late arrivals in the CVM-H simulation travel with an apparent velocity slower than observed. A higher-quality factor than traditionally assumed may produce synthetics with longer durations but is unable to accurately match the amplitude and phase. Beamforming analysis using dense array data further reveals the long duration surface waves have the same back azimuth as the direct arrivals and are generated at the basin edges, while the later coda waves are scattered from off-azimuth directions, potentially due to strong, sharp boundaries offshore. Improving the description of these shallow basin structures and attenuation model will enhance our capability to predict long-period ground motions in basins.

California↗

USGS2018

The 2015 Gorkha (Nepal) Earthquake sequence: I. Source modeling and deterministic 3D ground shaking

To better quantify the relatively long period (< 0.3 Hz) shaking experienced during the 2015 Gorkha (Nepal) earthquake sequence, we study the finite rupture processes and the associated 3D ground motion of the Mw7.8 mainshock and the Mw7.2 aftershock . The 3D synthetics are then used in the broadband ground shaking in Kathmandu with a hybrid approach, summarized in a companion paper (Chen and Wei, 2017, submitted together). We determined the coseismic rupture process of the mainshock by joint inversion of InSAR/SAR, GPS (static and high-rate), strong motion and teleseismic waveforms. Our inversion for the mainshock indicates unilateral rupture towards the ESE, with an average rupture speed of 3.0 km/s and a total duration of ~ 60 s. Additionally, we find that the beginning part of the rupture (5–18 s) has about 40% longer rise time than the rest of the rupture, as well as slower rupture velocity. Our model shows two strong asperities occurring ~ 24 s and ~ 36 s after the origin and located ~ 30 km to the northwest and northeast of the Kathmandu valley, respectively. In contrast, the Mw7.2 aftershock is more compact both in time and space, as revealed by joint inversion of teleseismic body waves and InSAR data. The different rupture features between the mainshock and the aftershock could be related to difference in fault zone structure. The mainshock and aftershock ground motions in the Kathmandu valley, recorded by both strong motion and high-rate GPS stations, exhibited strong amplification around 0.2 Hz. A simplified 3D basin model, calibrated by an Mw5.2 aftershock, can match the observed waveforms reasonably well at 0.3 Hz and lower frequency. The 3D simulations indicate that the basin structure trapped the wavefield and produced an extensive ground vibration. Our study suggests that the combination of rupture characteristics and propagational complexity are required to understand the ground shaking produced by hazardous earthquakes such as the Gorkha event.

Tectonophysics↗

USGS2017

Evidence for strong lateral seismic velocity variation in the lower crust and upper mantle beneath the California margin

Regional seismograms from earthquakes in Northern California show a systematic difference in arrival times across Southern California where long period (30–50 seconds) SH waves arrive up to 15 seconds earlier at stations near the coast compared with sites towards the east at similar epicentral distances. We attribute this time difference to heterogeneity of the velocity structure at the crust-mantle interface beneath the California margin. To model these observations, we propose a fast seismic layer, with thickness growing westward from the San Andreas along with a thicker and slower continental crust to the east. Synthetics generated from such a model are able to match the observed timing of SH waveforms better than existing 3D models. The presence of a strong upper mantle buttressed against a weaker crust has a major influence in how the boundary between the Pacific plate and North American plate deforms and may explain the observed asymmetric strain rate across the boundary.

California↗

USGS2015

The 2014 Mw6.1 South Napa Earthquake: A unilateral rupture with shallow asperity and rapid afterslip

The M w 6.1 South Napa earthquake occurred near Napa, California on August 24, 2014 (UTC), and was the largest inland earthquake in Northern California since the 1989 M w 6.9 Loma Prieta earthquake. The first report of the earthquake from the Northern California Earthquake Data Center (NCEDC) indicates a hypocentral depth of 11.0km with longitude and latitude of (122.3105°W, 38.217°N). Surface rupture was documented by field observations and Lidar imaging (Brooks et al. 2014; Hudnut et al. 2014; Brocher et al., 2015), with about 12 km of continuous rupture starting near the epicenter and extending to the northwest. The southern part of the rupture is relatively straight, but the strike changes by about 15° at the northern end over a 6-km segment. The peak dextral offset was observed near the Buhman residence with right-.‐lateral motion of 46 cm, near the location where the strike of fault begins to rotate clock-.‐wise (Hudnut et al., 2014). The earthquake was well recorded by the strong motion network operated by the NCEDC, the California Geological Survey and the U.S. Geological Survey (USGS). There are about 12 sites within an epicentral distance of 15km, with relatively good azimuthal coverage (Fig.1). The largest peak-ground-velocity (PGV) of nearly 100 cm/s was observed on station 1765, which is the closest station to the rupture and lies about 3 km east of the northern segment (Fig. 1). The ground deformation associated with the earthquake was also well recorded by the high-resolution COSMO-SkyMed satellite and Sentinel-1A satellite, providing independent static observations.

California↗

USGS2013

Rupture complexity of the Mw 8.3 sea of okhotsk earthquake: Rapid triggering of complementary earthquakes?

We derive a finite slip model for the 2013 M w 8.3 Sea of Okhotsk Earthquake ( Z = 610 km) by inverting calibrated teleseismic P waveforms. The inversion shows that the earthquake ruptured on a 10° dipping rectangular fault zone (140 km × 50 km) and evolved into a sequence of four large sub-events (E1–E4) with an average rupture speed of 4.0 km/s. The rupture process can be divided into two main stages. The first propagated south, rupturing sub-events E1, E2, and E4. The second stage (E3) originated near E2 with a delay of 12 s and ruptured northward, filling the slip gap between E1 and E2. This kinematic process produces an overall slip pattern similar to that observed in shallow swarms, except it occurs over a compressed time span of about 30 s and without many aftershocks, suggesting that sub-event triggering for deep events is significantly more efficient than for shallow events.

Geophysical Research Letters↗

USGS2012

Sources of shaking and flooding during the Tohoku-Oki earthquake: a mixture of rupture styles

Modeling strong ground motions from great subduction zone earthquakes is one of the great challenges of computational seismology. To separate the rupture characteristics from complexities caused by 3D sub-surface geology requires an extraordinary data set such as provided by the recent Mw9.0 Tohoku-Oki earthquake. Here we combine deterministic inversion and dynamically guided forward simulation methods to model over one thousand high-rate GPS and strong motion observations from 0 to 0.25 Hz across the entire Honshu Island. Our results display distinct styles of rupture with a deeper generic interplate event (~Mw8.5) transitioning to a shallow tsunamigenic earthquake (~Mw9.0) at about 25 km depth in a process driven by a strong dynamic weakening mechanism, possibly thermal pressurization. This source model predicts many important features of the broad set of seismic, geodetic and seafloor observations providing a major advance in our understanding of such great natural hazards.

Earth and Planetary Science Letters↗