Geology topics

James H. Dieterich

Publications and source records attributed to James H. Dieterich.

At least 19 recordsLinked to original sources

USGS2021

Toward physics-based nonergodic PSHA: A prototype fully-deterministic seismic hazard model for southern California

We present a nonergodic framework for probabilistic seismic‐hazard analysis (PSHA) that is constructed entirely of deterministic, physical models. The use of deterministic ground‐motion simulations in PSHA calculations is not new (e.g., CyberShake), but prior studies relied on kinematic rupture generators to extend empirical earthquake rupture forecasts. Fully dynamic models, which simulate rupture nucleation and propagation of static and dynamic stresses, are still computationally intractable for the large simulation domains and many seismic cycles required to perform PSHA. Instead, we employ the Rate‐State earthquake simulator (RSQSim) to efficiently simulate hundreds of thousands of years of M ≥ 6.5 earthquake sequences on the California fault system. RSQSim produces full slip‐time histories for each rupture, which, unlike kinematic models, emerge from frictional properties, fault geometry, and stress transfer; all intrinsic variability is deterministic. We use these slip‐time histories directly as input to a 3D wave‐propagation code within the CyberShake platform to obtain simulated F max = 0.5 Hz ground motions. The resulting 3 s spectral acceleration ground motions closely match empirical ground‐motion model (GMM) estimates of median and variability of shaking. When computed over a range of sources and sites, the variability is similar to that of ergodic GMMs. Variability is reduced for individual pairs of sources and sites that repeatedly sample a single path, which is expected for a nonergodic model. This results in increased exceedance probabilities for certain characteristic ground motions for a source–site pair, while decreasing probabilities at the extreme tails of the ergodic GMM predictions. We present these comparisons and preliminary fully deterministic physics‐based RSQSim–CyberShake hazard curves, as well as a new technique for estimating within‐ and between‐event variability through simulation.

California↗

USGS2018

A physics-based earthquake simulator replicates seismic hazard statistics across California

Seismic hazard models are important for society, feeding into building codes and hazard mitigation efforts. These models, however, rest on many uncertain assumptions and are difficult to test observationally because of the long recurrence times of large earthquakes. Physics-based earthquake simulators offer a potentially helpful tool, but they face a vast range of fundamental scientific uncertainties. We compare a physics-based earthquake simulator against the latest seismic hazard model for California. Using only uniform parameters in the simulator, we find strikingly good agreement of the long-term shaking hazard compared with the California model. This ability to replicate statistically based seismic hazard estimates by a physics-based model cross-validates standard methods and provides a new alternative approach needing fewer inputs and assumptions for estimating hazard.

California↗

USGS2017

Delayed seismicity rate changes controlled by static stress transfer

On 15 June 2010, a M w 5.7 earthquake occurred near Ocotillo, California, in the Yuha Desert. This event was the largest aftershock of the 4 April 2010 M w 7.2 El Mayor-Cucapah (EMC) earthquake in this region. The EMC mainshock and subsequent Ocotillo aftershock provide an opportunity to test the Coulomb failure hypothesis (CFS). We explore the spatiotemporal correlation between seismicity rate changes and regions of positive and negative CFS change imparted by the Ocotillo event. Based on simple CFS calculations we divide the Yuha Desert into three subregions, one triggering zone and two stress shadow zones. We find the nominal triggering zone displays immediate triggering, one stress shadowed region experiences immediate quiescence, and the other nominal stress shadow undergoes an immediate rate increase followed by a delayed shutdown. We quantitatively model the spatiotemporal variation of earthquake rates by combining calculations of CFS change with the rate-state earthquake rate formulation of Dieterich (1994), assuming that each subregion contains a mixture of nucleation sources that experienced a CFS change of differing signs. Our modeling reproduces the observations, including the observed delay in the stress shadow effect in the third region following the Ocotillo aftershock. The delayed shadow effect occurs because of intrinsic differences in the amplitude of the rate response to positive and negative stress changes and the time constants for return to background rates for the two populations. We find that rate-state models of time-dependent earthquake rates are in good agreement with the observed rates and thus explain the complex spatiotemporal patterns of seismicity.

California↗

USGS2012

Velocity dependent friction of granite over a wide range of conditions

Direct shear sliding experiments on bare ground surfaces of Westerly granite have been conducted over an exceptionally wide range of sliding rates (10 −4 µm/s to 10³ µm/s) at unconfined normal stresses (σ n ) of 5, 15, 30, 70, and 150 MPa. A new sample configuration was developed that permitted measurements at normal stresses of 70 and 150 MPa without immediate sample failure. Measurements of steady-state velocity dependence of friction at velocities between 10 −4 and 1 µm/s show similar velocity weakening behavior at all normal stresses, with more negative dependence at lower slip rates. However, at rates above 10 µm/s, velocity weakening is observed only at σ n =30, 70 and 150 MPa, while velocity neutral behavior is observed at σ n =15 MPa and velocity strengthening is observed at σ n =5 MPa. The greater velocity weakening observed at velocities below 10 −2 µm/s may suggest a transition in competing deformation mechanisms, or the influence of additional mechanisms. The transition to velocity strengthening at high velocity and low normal stress implies that rapid slip on shallow faults could be arrested before resulting in true stick-slip behavior. Stable fault creep and creep events observed at shallow levels on some natural faults may result from this transition in velocity dependence.

Geophysical Research Letters↗

USGS2012

Potential for geophysical experiments in large scale tests

Potential research applications for large-specimen geophysical experiments include measurements of scale dependence of physical parameters and examination of interactions with heterogeneities, especially flaws such as cracks. In addition, increased specimen size provides opportunities for improved recording resolution and greater control of experimental variables. Large-scale experiments using a special purpose low stress (<40 MPa) biaxial apparatus demonstrate that a minimum fault length is required to generate confined shear instabilities along pre-existing faults. Experimental analysis of source interactions for simulated earthquakes consisting of confined shear instabilities on a fault with gouge appears to require large specimens (˜1m) and high confining pressures (> 100 MPa).

Geophysical Research Letters↗

USGS2012

A comparison among observations and earthquake simulator results for the allcal2 California fault model

In order to understand earthquake hazards we would ideally have a statistical description of earthquakes for tens of thousands of years. Unfortunately the ∼100‐year instrumental, several 100‐year historical, and few 1000‐year paleoseismological records are woefully inadequate to provide a statistically significant record. Physics‐based earthquake simulators can generate arbitrarily long histories of earthquakes; thus they can provide a statistically meaningful history of simulated earthquakes. The question is, how realistic are these simulated histories? This purpose of this paper is to begin to answer that question. We compare the results between different simulators and with information that is known from the limited instrumental, historic, and paleoseismological data. As expected, the results from all the simulators show that the observational record is too short to properly represent the system behavior; therefore, although tests of the simulators against the limited observations are necessary, they are not a sufficient test of the simulators’ realism. The simulators appear to pass this necessary test. In addition, the physics‐based simulators show similar behavior even though there are large differences in the methodology. This suggests that they represent realistic behavior. Different assumptions concerning the constitutive properties of the faults do result in enhanced capabilities of some simulators. However, it appears that the similar behavior of the different simulators may result from the fault‐system geometry, slip rates, and assumed strength drops, along with the shared physics of stress transfer. This paper describes the results of running four earthquake simulators that are described elsewhere in this issue of Seismological Research Letters. The simulators ALLCAL (Ward, 2012), VIRTCAL (Sachs et al., 2012), RSQSim (Richards‐Dinger and Dieterich, 2012), and ViscoSim (Pollitz, 2012) were run on our most recent all‐California fault model, allcal2. With the exception of ViscoSim, which ran for 10,000 years, all the simulators ran for 30,000 years. Presentations containing content similar to this paper can be found at http://scec.usc.edu/research/eqsims/.

Seismological Research Letters↗

USGS2012

Modeling of rock friction 1. Experimental results and constitutive equations

Direct shear experiments on ground surfaces of a granodiorite from Raymond, California, at normal stresses of ??6 MPa demonstrate that competing time, displacement, and velocity effects control rock friction. It is proposed that the strength of the population of points of contacts between sliding surfaces determines frictional strength and that the population of contacts changes continuously with displacements. Previous experiments demonstrate that the strength of the contacts increases with the age of the contacts. The present experiments establish that a characteristic displacement, proportional to surface roughness, is required to change the population of contacts. Hence during slip the average age of the points of contact and therefore frictional strength decrease as slip velocity increases. Displacement weakening and consequently the potential for unstable slip occur whenever displacement reduces the average age of the contacts. In addition to this velocity dependency, which arises from displacement dependency and time dependency, the experiments also show a competing but transient increase in friction whenever slip velocity increases. Creep of the sliding surface at stresses below that for steady state slip is also observed. Constitutive relationships are developed that permit quantitative simulation of the friction versus displacement data as a function of surface roughness and for different time and velocity histories. Unstable slip in experiments is controlled by these constitutive effects and by the stiffness of the experimental system. It is argued that analogous properties control earthquake instability. Copyright ?? 1979 by the American Geophysical Union.

Conference Paper↗

USGS2000

Heightened odds of large earthquakes near Istanbul: an interaction-based probability calculation

We calculate the probability of strong shaking in Istanbul, an urban center of 10 million people, from the description of earthquakes on the North Anatolian fault system in the Marmara Sea during the past 500 years and test the resulting catalog against the frequency of damage in Istanbul during the preceding millennium, departing from current practice, we include the time-dependent effect of stress transferred by the 1999 moment magnitude M = 7.4 Izmit earthquake to faults nearer to Istanbul. We find a 62 ± 15% probability (one standard deviation) of strong shaking during the next 30 years and 32 ± 12% during the next decade.

Science↗

USGS2000

High magma storage rates before the 1983 eruption of Kilauea, Hawaii

After a magnitude 7.2 earthquake in 1975 and before the start of the ongoing eruption in 1983, deformation of Kilauea volcano was the most rapid ever recorded. Three-dimensional numerical modeling shows that this deformation is consistent with the dilation of a dike within Kilauea's rift zones coupled with creep over a narrow area of a low-angle fault beneath the south flank. Magma supply is estimated to be 0.18 cubic kilometers per year, twice that of previous estimates. The 1983 eruption may be a direct consequence of the high rates of magma storage within the rift zone that followed the 1975 earthquake.

Science↗

USGS1998

Stress transferred by the 1995 Mw = 6.9 Kobe, Japan, shock: Effect on aftershocks and future earthquake probabilities

The Kobe earthquake struck at the edge of the densely populated Osaka-Kyoto corridor in southwest Japan. We investigate how the earthquake transferred stress to nearby faults, altering their proximity to failure and thus changing earthquake probabilities. We find that relative to the pre-Kobe seismicity, Kobe aftershocks were concentrated in regions of calculated Coulomb stress increase and less common in regions of stress decrease. We quantify this relationship by forming the spatial correlation between the seismicity rate change and the Coulomb stress change. The correlation is significant for stress changes greater than 0.2–1.0 bars (0.02–0.1 MPa), and the nonlinear dependence of seismicity rate change on stress change is compatible with a state- and rate-dependent formulation for earthquake occurrence. We extend this analysis to future mainshocks by resolving the stress changes on major faults within 100 km of Kobe and calculating the change in probability caused by these stress changes. Transient effects of the stress changes are incorporated by the state-dependent constitutive relation, which amplifies the permanent stress changes during the aftershock period. Earthquake probability framed in this manner is highly time-dependent, much more so than is assumed in current practice. Because the probabilities depend on several poorly known parameters of the major faults, we estimate uncertainties of the probabilities by Monte Carlo simulation. This enables us to include uncertainties on the elapsed time since the last earthquake, the repeat time and its variability, and the period of aftershock decay. We estimate that a calculated 3-bar (0.3-MPa) stress increase on the eastern section of the Arima-Takatsuki Tectonic Line (ATTL) near Kyoto causes fivefold increase in the 30-year probability of a subsequent large earthquake near Kyoto; a 2-bar (0.2-MPa) stress decrease on the western section of the ATTL results in a reduction in probability by a factor of 140 to 2000. The probability of a M w = 6.9 earthquake within 50 km of Osaka during 1997–2007 is estimated to have risen from 5–6% before the Kobe earthquake to 7–11% afterward; during 1997–2027, it is estimated to have risen from 14–16% before Kobe to 16–22%.

Journal of Geophysical Research B: Solid Earth↗

USGS1996

Implications of fault constitutive properties for earthquake prediction

The rate- and state-dependent constitutive formulation for fault slip characterizes an exceptional variety of materials over a wide range of sliding conditions. This formulation provides a unified representation of diverse sliding phenomena including slip weakening over a characteristic sliding distance D(c), apparent fracture energy at a rupture front, time- dependent healing after rapid slip, and various other transient and slip rate effects. Laboratory observations and theoretical models both indicate that earthquake nucleation is accompanied by long intervals of accelerating slip. Strains from the nucleation process on buried faults generally could not be detected if laboratory values of D, apply to faults in nature. However, scaling of D(c) is presently an open question and the possibility exists that measurable premonitory creep may precede some earthquakes. Earthquake activity is modeled as a sequence of earthquake nucleation events. In this model, earthquake clustering arises from sensitivity of nucleation times to the stress changes induced by prior earthquakes. The model gives the characteristic Omori aftershock decay law and assigns physical interpretation to aftershock parameters. The seismicity formulation predicts large changes of earthquake probabilities result from stress changes. Two mechanisms for foreshocks are proposed that describe observed frequency of occurrence of foreshock-mainshock pairs by time and magnitude. With the first mechanism, foreshocks represent a manifestation of earthquake clustering in which the stress change at the time of the foreshock increases the probability of earthquakes at all magnitudes including the eventual mainshock. With the second model, accelerating fault slip on the mainshock nucleation zone triggers foreshocks.

Conference Paper↗

USGS1996

Imaging surface contacts: Power law contact distributions and contact stresses in quartz, calcite, glass and acrylic plastic

A procedure has been developed to obtain microscope images of regions of contact between roughened surfaces of transparent materials, while the surfaces are subjected to static loads or undergoing frictional slip. Static loading experiments with quartz, calcite, soda-lime glass and acrylic plastic at normal stresses to 30 MPa yield power law distributions of contact areas from the smallest contacts that can be resolved (3.5 μm 2 ) up to a limiting size that correlates with the grain size of the abrasive grit used to roughen the surfaces. In each material, increasing normal stress results in a roughly linear increase of the real area of contact. Mechanisms of contact area increase are by growth of existing contacts, coalescence of contacts and appearance of new contacts. Mean contacts stresses are consistent with the indentation strength of each material. Contact size distributions are insensitive to normal stress indicating that the increase of contact area is approximately self-similar. The contact images and contact distributions are modeled using simulations of surfaces with random fractal topographies. The contact process for model fractal surfaces is represented by the simple expedient of removing material at regions where surface irregularities overlap. Synthetic contact images created by this approach reproduce observed characteristics of the contacts and demonstrate that the exponent in the power law distributions depends on the scaling exponent used to generate the surface topography.

Tectonophysics↗

USGS1996

An unusual pattern of recurring seismic quiescence at Kalapana, Hawaii

An unusual pattern of recurring seismic quiescence is observed in the Kalapana, Hawaii region of Kilauea Volcano. Statistically significant intervals of quiescence preceded the Kalapana earthquakes of 1975 (M7.2) and 1989 (M6.1) and a third quiescence is presently underway. The sensitivity of the volcano flank to continuing magmatic activity in the nearby east rift zone complicates interpretation of these observations. The current quiescence episode may be caused by magmatic processes in the east rift zone or by changes within the flank of Kilauea. The latter possibility, if correct, may represent a precursor to another earthquake.

Geophysical Research Letters↗

USGS1995

Mechanical modeling of circumferential and radial dike intrusion on Galapagos volcanoes

A distinctive and unusual pattern of eruptive fissures is observed on the active volcanoes of the Galapagos islands, reflecting circumferential dike intrusion near the calderas and radial dike intrusion on the volcano flanks. Elastic finite-element models were used to investigate how a stress field could be produced and maintained to promote both circumferential and radial dike emplacement. Modeling results show that magma reservoirs of Galapagos volcanoes are probably diapiric, because this shape promotes both circumferential and radial intrusions, but magma pressure alone cannot create the observed pattern of dikes. Loading by volcano growth and magma reservoir pressure could produce a stress field of suitable orientation but insufficient magnitude. The intrusion of circumferential dikes could alter the stress field in a way that promotes future radial diking, and vice versa. Faulting or slumping within the calderas or on the volcano flanks in response to repeated intrusions could also create a stress field conducive to continued intrusion.

Galápagos Islands, Fernandina, Isabela↗

USGS1994

Direct observation of frictional contacts: New insights for state-dependent properties

Rocks and many other materials display a rather complicated, but characteristic, dependence of friction on sliding history. These effects are well-described by empirical rate- and state-dependent constitutive formulations which have been utilized for analysis of fault slip and earthquake processes. We present a procedure for direct quantitative microscopic observation of frictional contacts during slip. The observations reveal that frictional state dependence represents an increase of contact area with contact age. Transient changes of sliding resistance correlate with changes in contact area and arise from shifts of contact population age. Displacement-dependent replacement of contact populations is shown to cause the diagnostic evolution of friction over a characteristic sliding distance that occurs whenever slip begins or sliding conditions change. ?? 1994 Birkha??user Verlag.

Pure and Applied Geophysics PAGEOPH↗

USGS1992

Earthquake nucleation on faults with rate-and state-dependent strength

Faults with rate- and state-dependent constitutive properties reproduce a range of observed fault slip phenomena including spontaneous nucleation of slip instabilities at stresses above some critical stress level and recovery of strength following slip instability. Calculations with a plane-strain fault model with spatially varying properties demonstrate that accelerating slip precedes instability and becomes localized to a fault patch. The dimensions of the fault patch follow scaling relations for the minimum critical length for unstable fault slip. The critical length is a function of normal stress, loading conditions and constitutive parameters which include D c , the characteristic slip distance. If slip starts on a patch that exceeds the critical size, the length of the rapidly accelerating zone tends to shrink to the characteristic size as the time of instability approaches. Solutions have been obtained for a uniform, fixed-patch model that are in good agreement with results from the plane-strain model. Over a wide range of conditions, above the steady-state stress, the logarithm of the time to instability linearly decreases as the initial stress increases. Because nucleation patch length and premonitory displacement are proportional to D c , the moment of premonitory slip scales by D 3 c . The scaling of D c is currently an open question. Unless D c for earthquake faults is significantly greater than that observed on laboratory faults, premonitory strain arising from the nucleation process for earthquakes may by too small to detect using current observation methods. Excluding the possibility that D c in the nucleation zone controls the magnitude of the subsequent earthquake, then the source dimensions of the smallest earthquakes in a region provide an upper limit for the size of the nucleation patch.

Tectonophysics↗

USGS1989

The frictional properties of a simulated gouge having a fractal particle distribution

The frictional properties of a layer of simulated Westerly granite fault gouge sandwiched between sliding blocks of Westerly granite have been measured in a high-speed servo-controlled double-direct shear apparatus. Most gouge layers were prepared to have a self-similar particle distribution with a fractal dimension of 2.6. The upper fractal limit was varied between 45 and 710 μm. Some gouges were prepared with all particles in the range between 360 and 710 μm. In each experiment the sliding velocity was cyclically alternated between 1 and 10 μms −1 and the coefficient of friction μm and its transient parameters a, b and D c were measured as functions of displacement. In addition to the particle size distribution, the following experimental variables were also investigated: the layer thickness (1 and 3 mm), the roughness of the sliding surfaces (Nos 60 and 600 grit) and the normal stress (10 and 25 MPa). Some of the sample assemblies were epoxy impregnated following a run so the gouge structure could be microscopically examined in thin section. We observed that gouges which were initially non-fractal evolved to a fractal distribution with dimension 2.6. Gouges which had an initial fractal distribution remained fractal. When the sliding blocks had smooth surfaces, the coefficient of friction was relatively low and was independent of the particle distribution. In these cases, strong velocity weakening was observed throughout the experiment and the transient parameters a, b and D c , remained almost constant. When the sliding blocks had rough surfaces, the coefficient of friction was larger and more dependent on the particle distribution. Velocity strengthening was observed initially but evolved to velocity weakening with increased sliding displacement. All three transient parameters changed with increasing displacement. The a and b values were about three times as large for rough surfaces as for smooth. The characteristic displacement D c was not sensitive to surface roughness but was the only transient parameter which was sensitive to the normal stress. For the case of rough surfaces, the coefficient of friction of the 1 mm thick gouge was significantly larger than that for the 3 mm thick layers. Many of these observations can be explained by a micromechanical model in which the stress in the gouge layer is heterogeneous. The applied normal and shear stresses are supported by ‘grain bridges’ which span the layer and which are continually forming and failing. In this model, the frictional properties of the gouge are largely determined by the dominant failure mode of the bridging structures.

Journal of Structural Geology↗

USGS1987

Nucleation and triggering of earthquake slip: Effect of periodic stresses

Results of stability analyses for spring and slider systems, with state variable constitutive properties, are applied to slip on embedded fault patches. Unstable slip may nucleate only if the slipping patch exceeds some minimum size. Subsequent to the onset of instability the earthquake slip may propagate well beyond the patch. It is proposed that the seismicity of a volume of the earth's crust is determined by the distribution of initial conditions on the population of fault patches that nucleate earthquake slip, and the loading history acting upon the volume. Patches with constitutive properties inferred from laboratory experiments are characterized by an interval of self-driven accelerating slip prior to instability, if initial stress exceeds a minimum threshold. This delayed instability of the patches provides an explanation for the occurrence of aftershocks and foreshocks including decay of earthquake rates by time −1 . A population of patches subjected to loading with a periodic component results in periodic variation of the rate of occurrence of instabilities. The change of the rate of seismicity for a sinusoidal load is proportional to the amplitude of the periodic stress component and inversely proportional to both the normal stress acting on the fault patches and the constitutive parameter, A 1 , that controls the direct velocity dependence of fault slip. Values of A 1 representative of laboratory experiments indicate that in a homogeneous crust, correlation of earthquake rates with earth tides should not be detectable at normal stresses in excess of about 8 MPa. Correlation of earthquakes with tides at higher normal stresses can be explained if there exist inhomogeneities that locally amplify the magnitude of the tidal stresses. Such amplification might occur near magma chambers or other soft inclusions in the crust and possibly near the ends of creeping fault segments if the creep or afterslip rates vary in response to tides. Observations of seismicity rate variations associated with seasonal fluctuations of reservoir levels appear to be consistent with the model.

Tectonophysics↗