Geology topics

G. C. Beroza

Publications and source records attributed to G. C. Beroza.

8 recordsLinked to original sources

USGS2003

Strong ground-motion prediction from Stochastic-dynamic source models

In the absence of sufficient data in the very near source, predictions of the intensity and variability of ground motions from future large earthquakes depend strongly on our ability to develop realistic models of the earthquake source. In this article we simulate near-fault strong ground motion using dynamic source models. We use a boundary integral method to simulate dynamic rupture of earthquakes by specifying dynamic source parameters (fracture energy and stress drop) as spatial random fields. We choose these quantities such that they are consistent with the statistical properties of slip heterogeneity found in finite-source models of past earthquakes. From these rupture models we compute theoretical strong-motion seismograms up to a frequency of 2 Hz for several realizations of a scenario strike-slip Mw 7.0 earthquake and compare empirical response spectra, spectra obtained from our dynamic models, and spectra determined from corresponding kinematic simulations. We find that spatial and temporal variations in slip, slip rise time, and rupture propagation consistent with dynamic rupture models exert a strong influence on near-source ground motion. Our results lead to a feasible approach to specify the variability in the rupture time distribution in kinematic models through a generalization of Andrews' (1976) result relating rupture speed to apparent fracture energy, stress drop, and crack length to 3D dynamic models. This suggests that a simplified representation of dynamic rupture may be obtained to approximate the effects of dynamic rupture without having to do full dynamic simulations.

Bulletin of the Seismological Society of America↗

USGS2003

Apparent break in earthquake scaling due to path and site effects on deep borehole recordings

We reexamine the scaling of stress drop and apparent stress, rigidity times the ratio between seismically radiated energy to seismic moment, with earthquake size for a set of microearthquakes recorded in a deep borehole in Long Valley, California. In the first set of calculations, we assume a constant Q and solve for the corner frequency and seismic moment. In the second set of calculations, we model the spectral ratio of nearby events to determine the same quantities. We find that the spectral ratio technique, which can account for path and site effects or nonconstant Q, yields higher stress drops, particularly for the smaller events in the data set. The measurements determined from spectral ratios indicate no departure from constant stress drop scaling down to the smallest events in our data set (Mw 0.8). Our results indicate that propagation effects can contaminate measurements of source parameters even in the relatively clean recording environment of a deep borehole, just as they do at the Earth's surface. The scaling of source properties of microearthquakes made from deep borehole recordings may need to be reevaluated.

Journal of Geophysical Research B: Solid Earth↗

USGS2002

High-resolution image of Calaveras Fault seismicity

By measuring relative earthquake arrival times using waveform cross correlation and locating earthquakes using the double difference technique, we are able to reduce hypocentral errors by 1 to 2 orders of magnitude over routine locations for nearly 8000 events along a 35-km section of the Calaveras Fault. This represents ∼92% of all seismicity since 1984 and includes the rupture zone of the M 6.2 1984 Morgan Hill, California, earthquake. The relocated seismicity forms highly organized structures that were previously obscured by location errors. There are abundant repeating earthquake sequences as well as linear clusters of earthquakes. Large voids in seismicity appear with dimensions of kilometers that have been aseismic over the 30-year time interval, suggesting that these portions of the fault are either locked or creeping. The area of greatest slip in the Morgan Hill main shock coincides with the most prominent of these voids, suggesting that this part of the fault may be locked between large earthquakes. We find that the Calaveras Fault at depth is extremely thin, with an average upper bound on fault zone width of 75 m. Given the location error, however, this width is not resolvably different from zero. The relocations reveal active secondary faults, which we use to solve for the stress field in the immediate vicinity of the Calaveras Fault. We find that the maximum compressive stress is at a high angle, only 13° from the fault normal, supporting previous interpretations that this fault is weak.

California↗

USGS2001

Inferring rate and state friction parameters from a rupture model of the 1995 Hyogo-ken Nanbu (Kobe) Japan earthquake

We consider the applicability of laboratory-derived rate-and state-variable friction laws to the dynamic rupture of the 1995 Kobe earthquake. We analyze the shear stress and slip evolution of Ide and Takeo's [1997] dislocation model, fitting the inferred stress change time histories by calculating the dynamic load and the instantaneous friction at a series of points within the rupture area. For points exhibiting a fast-weakening behavior, the Dieterich-Ruina friction law, with values of d c = 0.01–0.05 m for critical slip, fits the stress change time series well. This range of d c is 10–20 times smaller than the slip distance over which the stress is released, D c , which previous studies have equated with the slip-weakening distance. The limited resolution and low-pass character of the strong motion inversion degrades the resolution of the frictional parameters and suggests that the actual d c is less than this value. Stress time series at points characterized by a slow-weakening behavior are well fitted by the Dieterich-Ruina friction law with values of d c ≥0.01–0.05 m. The apparent fracture energy G c can be estimated from waveform inversions more stably than the other friction parameters. We obtain a G c ≈ 1.5×l0 6 J m −2 for the 1995 Kobe earthquake, in agreement with estimates for previous earthquakes. From this estimate and a plausible upper bound for the local rock strength we infer a lower bound for D c of about 0.008 m.

Kobe, Osaka↗

USGS1998

Observation of the seismic nucleation phase in the Ridgecrest, California, earthquake sequence

Near-source observations of five M 3.8–5.2 earthquakes near Ridgecrest, California are consistent with the presence of a seismic nucleation phase. These earthquakes start abruptly, but then slow or stop before rapidly growing again toward their maximum rate of moment release. Deconvolution of instrument and path effects by empirical Green's functions demonstrates that the initial complexity at the start of the earthquake is a source effect. The rapid growth of the P -wave arrival at the start of the seismic nucleation phase supports the conclusion of Mori and Kanamori [1996] that these earthquakes begin without a magnitude-scaled slow initial phase of the type observed by Iio [1992, 1995].

Geophysical Research Letters↗

USGS1996

Detailed observations of California foreshock sequences: Implications for the earthquake initiation process

We find that foreshocks provide clear evidence for an extended nucleation process before some earthquakes. In this study, we examine in detail the evolution of six California foreshock sequences, the 1986 Mount Lewis ( M L = 5.5), the 1986 Chalfant ( M L = 6.4), the 1986 Stone Canyon ( M L = 4.7), the 1990 Upland ( M L = 5.2), the 1992 Joshua Tree ( M W = 6.1), and the 1992 Landers ( M W = 7.3) sequence. Typically, uncertainties in hypocentral parameters are too large to establish the geometry of foreshock sequences and hence to understand their evolution. However, the similarity of location and focal mechanisms for the events in these sequences leads to similar foreshock waveforms that we cross correlate to obtain extremely accurate relative locations. We use these results to identify small-scale fault zone structures that could influence nucleation and to determine the stress evolution leading up to the mainshock. In general, these foreshock sequences are not compatible with a cascading failure nucleation model in which the foreshocks all occur on a single fault plane and trigger the mainshock by static stress transfer. Instead, the foreshocks seem to concentrate near structural discontinuities in the fault and may themselves be a product of an aseismic nucleation process. Fault zone heterogeneity may also be important in controlling the number of foreshocks, i.e., the stronger the heterogeneity, the greater the number of foreshocks. The size of the nucleation region, as measured by the extent of the foreshock sequence, appears to scale with mainshock moment in the same manner as determined independently by measurements of the seismic nucleation phase. We also find evidence for slip localization as predicted by some models of earthquake nucleation.

Journal of Geophysical Research B: Solid Earth↗

USGS1996

Properties of the seismic nucleation phase

Near-source observations show that earthquakes begin abruptly at the P-wave arrival, but that this beginning is weak, with a low moment rate relative to the rest of the main shock. We term this initial phase of low moment rate the seismic nucleation phase. We have observed the seismic nucleation phase for a set of 48 earthquakes ranging in magnitude from 1.1–8.1 The size and duration of the seismic nucleation phase scale with the total seismic moment of the earthquake, suggesting that the process responsible for the seismic nucleation phase carries information about the eventual size of the earthquake. The seismic nucleation phase is characteristically followed by quadratic growth in the moment rate, consistent with self-similar rupture at constant stress drop. In this paper we quantify the properties of the seismic nucleation phase and offer several possible explanations for it.

California↗

USGS1995

Seismic evidence for an earthquake nucleation phase

Near-source observations show that earthquakes initiate with a distinctive seismic nucleation phase that is characterized by a low rate of moment release relative to the rest of the event. This phase was observed for the 30 earthquakes having moment magnitudes 2.6 to 8.1, and the size and duration of this phase scale with the eventual size of the earthquake. During the nucleation phase, moment release was irregular and appears to have been confined to a limited region of the fault. It was characteristically followed by quadratic growth in the moment rate as rupture began to propagate away from the nucleation zone. These observations suggest that the nucleation process exerts a strong influence on the size of the eventual earthquake.

Science↗