Search USGS⌕ Search

Geology topics

J. Gomberg

Publications and source records attributed to J. Gomberg.

35 records · Page 2Linked to original sources

Lithology and shear-wave velocity in Memphis, Tennessee

We have derived a new three-dimensional model of the lithologic structure beneath the city of Memphis , Tennessee , and examined its correlation with measured shear - wave velocity profiles. The correlation is sufficiently high that the better-constrained lithologic model may be used as a proxy for shear - wave velocities, which are required to calculate site-amplification for new seismic hazard maps for Memphis . The lithologic model and its uncertainties are derived from over 1200 newly compiled well and boring logs, some sampling to 500 m depth, and a moving-least-squares algorithm. Seventy-six new shear - wave velocity profiles have been measured and used for this study, most sampling to 30 m depth or less. All log and velocity observations are publicly available via new web sites.

Tennessee↗

The Mw 7.7 Bhuj earthquake: Global lessons for earthquake hazard in intra-plate regions

The Mw 7.7 Bhuj earthquake occurred in the Kachchh District of the State of Gujarat, India on 26 January 2001, and was one of the most damaging intraplate earthquakes ever recorded. This earthquake is in many ways similar to the three great New Madrid earthquakes that occurred in the central United States in 1811-1812, An Indo-US team is studying the similarities and differences of these sequences in order to learn lessons for earthquake hazard in intraplate regions. Herein we present some preliminary conclusions from that study. Both the Kutch and New Madrid regions have rift type geotectonic setting. In both regions the strain rates are of the order of 10-9/yr and attenuation of seismic waves as inferred from observations of intensity and liquefaction are low. These strain rates predict recurrence intervals for Bhuj or New Madrid sized earthquakes of several thousand years or more. In contrast, intervals estimated from paleoseismic studies and from other independent data are significantly shorter, probably hundreds of years. All these observations together may suggest that earthquakes relax high ambient stresses that are locally concentrated by rheologic heterogeneities, rather than loading by plate-tectonic forces. The latter model generally underlies basic assumptions made in earthquake hazard assessment, that the long-term average rate of energy released by earthquakes is determined by the tectonic loading rate, which thus implies an inherent average periodicity of earthquake occurrence. Interpreting the observations in terms of the former model therefore may require re-examining the basic assumptions of hazard assessment.

Gujarat↗

Aftershock triggering by complete Coulomb stress changes

We examine the correlation between seismicity rate change following the 1992, M 7.3, Landers, California, earthquake and characteristics of the complete Coulomb failure stress (CFS) changes (ΔCFS( t )) that this earthquake generated. At close distances the time-varying “dynamic” portion of the stress change depends on how the rupture develops temporally and spatially and arises from radiated seismic waves and from permanent coseismic fault displacement. The permanent “static” portion (ΔCFS) depends only on the final coseismic displacement. ΔCFS diminishes much more rapidly with distance than the transient, dynamic stress changes. A common interpretation of the strong correlation between ΔCFS and aftershocks is that load changes can advance or delay failure. Stress changes may also promote failure by physically altering properties of the fault or its environs. Because it is transient, ΔCFS( t ) can alter the failure rate only by the latter means. We calculate both ΔCFS and the maximum positive value of ΔCFS( t ) (peak ΔCFS( t )) using a reflectivity program. Input parameters are constrained by modeling Landers displacement seismograms. We quantify the correlation between maps of seismicity rate changes and maps of modeled ΔCFS and peak ΔCFS( t ) and find agreement for both models. However, rupture directivity, which does not affect ΔCFS, creates larger peak ΔCFS( t ) values northwest of the main shock. This asymmetry is also observed in seismicity rate changes but not in ΔCFS. This result implies that dynamic stress changes are as effective as static stress changes in triggering aftershocks and may trigger earthquakes long after the waves have passed.

Journal of Geophysical Research B: Solid Earth↗

The failure of earthquake failure models

In this study I show that simple heuristic models and numerical calculations suggest that an entire class of commonly invoked models of earthquake failure processes cannot explain triggering of seismicity by transient or "dynamic" stress changes, such as stress changes associated with passing seismic waves. The models of this class have the common feature that the physical property characterizing failure increases at an accelerating rate when a fault is loaded (stressed) at a constant rate. Examples include models that invoke rate state friction or subcritical crack growth, in which the properties characterizing failure are slip or crack length, respectively. Failure occurs when the rate at which these grow accelerates to values exceeding some critical threshold. These accelerating failure models do not predict the finite durations of dynamically triggered earthquake sequences (e.g., at aftershock or remote distances). Some of the failure models belonging to this class have been used to explain static stress triggering of aftershocks. This may imply that the physical processes underlying dynamic triggering differs or that currently applied models of static triggering require modification. If the former is the case, we might appeal to physical mechanisms relying on oscillatory deformations such as compaction of saturated fault gouge leading to pore pressure increase, or cyclic fatigue. However, if dynamic and static triggering mechanisms differ, one still needs to ask why static triggering models that neglect these dynamic mechanisms appear to explain many observations. If the static and dynamic triggering mechanisms are the same, perhaps assumptions about accelerating failure and/or that triggering advances the failure times of a population of inevitable earthquakes are incorrect.

Journal of Geophysical Research B: Solid Earth↗

Earthquake triggering by seismic waves following the landers and hector mine earthquakes

The proximity and similarity of the 1992, magnitude 7.3 Landers and 1999, magnitude 7.1 Hector Mine earthquakes in California permit testing of earthquake triggering hypotheses not previously possible. The Hector Mine earthquake confirmed inferences that transient, oscillatory 'dynamic' deformations radiated as seismic waves can trigger seismicity rate increases, as proposed for the Landers earthquake1-6. Here we quantify the spatial and temporal patterns of the seismicity rate changes7. The seismicity rate increase was to the north for the Landers earthquake and primarily to the south for the Hector Mine earthquake. We suggest that rupture directivity results in elevated dynamic deformations north and south of the Landers and Hector Mine faults, respectively, as evident in the asymmetry of the recorded seismic velocity fields. Both dynamic and static stress changes seem important for triggering in the near field with dynamic stress changes dominating at greater distances. Peak seismic velocities recorded for each earthquake suggest the existence of, and place bounds on, dynamic triggering thresholds. These thresholds vary from a few tenths to a few MPa in most places, depend on local conditions, and exceed inferred static thresholds by more than an order of magnitude. At some sites, the onset of triggering was delayed until after the dynamic deformations subsided. Physical mechanisms consistent with all these observations may be similar to those that give rise to liquefaction or cyclic fatigue.

Nature↗

On rate-state and Coulomb failure models

We examine the predictions of Coulomb failure stress and rate-state frictional models. We study the change in failure time (clock advance) Δt due to stress step perturbations (i.e., coseismic static stress increases) added to "background" stressing at a constant rate (i.e., tectonic loading) at time t 0 . The predictability of Δt implies a predictable change in seismicity rate r(t)/r 0 , testable using earthquake catalogs, where r 0 is the constant rate resulting from tectonic stressing. Models of r(t)/r 0 , consistent with general properties of aftershock sequences, must predict an Omori law seismicity decay rate, a sequence duration that is less than a few percent of the mainshock cycle time and a return directly to the background rate. A Coulomb model requires that a fault remains locked during loading, that failure occur instantaneously, and that Δt is independent of t 0 . These characteristics imply an instantaneous infinite seismicity rate increase of zero duration. Numerical calculations of r(t)/r 0 for different state evolution laws show that aftershocks occur on faults extremely close to failure at the mainshock origin time, that these faults must be "Coulomb-like," and that the slip evolution law can be precluded. Real aftershock population characteristics also may constrain rate-state constitutive parameters; a may be lower than laboratory values, the stiffness may be high, and/or normal stress may be lower than lithostatic. We also compare Coulomb and rate-state models theoretically. Rate-state model fault behavior becomes more Coulomb-like as constitutive parameter a decreases relative to parameter b. This is because the slip initially decelerates, representing an initial healing of fault contacts. The deceleration is more pronounced for smaller a, more closely simulating a locked fault. Even when the rate-state Δt has Coulomb characteristics, its magnitude may differ by some constant dependent on b. In this case, a rate-state model behaves like a modified Coulomb failure model in which the failure stress threshold is lowered due to weakening, increasing the clock advance. The deviation from a non-Coulomb response also depends on the loading rate, elastic stiffness, initial conditions, and assumptions about how state evolves.

Journal of Geophysical Research B: Solid Earth↗

Triggering of earthquake aftershocks by dynamic stresses

It is thought that small 'static' stress changes due to permanent fault displacement can alter the likelihood of, or trigger, earthquakes on nearby faults. Many studies of triggering in the nearfield, particularly of aftershocks, rely on these static changes as the triggering agent and consider them only in terms of equivalent changes in the applied load on the fault. Here we report a comparison of the aftershock pattern of the moment magnitude MW = 7.3 Landers earthquake, not only with static stress changes but also with transient, oscillatory stress changes transmitted as seismic waves (that is, 'dynamic' stresses). Dynamic stresses do not permanently change the applied load and thus can trigger earthquakes only by altering the mechanical state or properties of the fault zone. These dynamically weakened faults may fail after the seismic waves have passed by, and might even cause earthquakes that would not otherwise have occurred. We find similar asymmetries in the aftershock and dynamic stress patterns, the latter being due to rupture propagation, whereas the static stress changes lack this asymmetry. Previous studies have shown that dynamic stresses can promote failure at remote distances, but here we show that they can also do so nearby.

Nature↗

The strain in the array is mainly in the plane (waves below ~1 Hz)

We compare geodetic and single-station methods of measuring dynamic deformations and characterize their causes in the frequency bands 0.5-1.0 Hz and 4.0-8.0 Hz. The geodetic approach utilizes data from small-aperture seismic arrays, applying techniques from geodesy. It requires relatively few assumptions and a priori information. The single-station method uses ground velocities recorded at isolated or single stations and assumes all the deformation is due to plane-wave propagation. It also requires knowledge of the azimuth and horizontal velocity of waves arriving at the recording station. Data employed come from a small-aperture, dense seismic array deployed in Geyokcha, Turkmenistan, and include seismograms recorded by broadband STS2 and short-period L28 sensors. Poor agreement between geodetic and single-station estimates in the 4.0-8.0 Hz passband indicates that the displacement field may vary nonlinearly with distance over distances of ∼50 m. STS2 geodetic estimates provide a robust standard in the 0.5-1.0 Hz passband because they appear to be computationally stable and require fewer assumptions than single-station estimates. The agreement between STS2 geodetic estimates and single-station L28 estimates is surprisingly good for the S -wave and early surface waves, suggesting that the single-station analysis should be useful with commonly available data. These results indicate that, in the 0.5 to 1.0 Hz passband, the primary source of dynamic deformation is plane-wave propagation along great-circle source-receiver paths. For later arriving energy, the effects of scattering become important. The local structure beneath the array exerts a strong control on the geometry of the dynamic deformation, implying that it may be difficult to infer source characteristics of modern or paleoearthquakes from indicators of dynamic deformations. However, strong site control also suggests that the dynamic deformations may be predictable, which would be useful for engineering seismically resistant structures.

Bulletin of the Seismological Society of America↗

The initial subevent of the 1994 Northridge, California, earthquake: Is earthquake size predictable?

We examine the initial subevent (ISE) of the M 6.7, 1994 Northridge, California, earthquake in order to discriminate between two end-member rupture initiation models: the ‘preslip’ and ‘cascade’ models. Final earthquake size may be predictable from an ISE's seismic signature in the preslip model but not in the cascade model. In the cascade model ISEs are simply small earthquakes that can be described as purely dynamic ruptures. In this model a large earthquake is triggered by smaller earthquakes; there is no size scaling between triggering and triggered events and a variety of stress transfer mechanisms are possible. Alternatively, in the preslip model, a large earthquake nucleates as an aseismically slipping patch in which the patch dimension grows and scales with the earthquake's ultimate size; the byproduct of this loading process is the ISE. In this model, the duration of the ISE signal scales with the ultimate size of the earthquake, suggesting that nucleation and earthquake size are determined by a more predictable, measurable, and organized process. To distinguish between these two end-member models we use short period seismograms recorded by the Southern California Seismic Network. We address questions regarding the similarity in hypocenter locations and focal mechanisms of the ISE and the mainshock. We also compare the ISE's waveform characteristics to those of small earthquakes and to the beginnings of earthquakes with a range of magnitudes. We find that the focal mechanisms of the ISE and mainshock are indistinguishable, and both events may have nucleated on and ruptured the same fault plane. These results satisfy the requirements for both models and thus do not discriminate between them. However, further tests show the ISE's waveform characteristics are similar to those of typical small earthquakes in the vicinity and more importantly, do not scale with the mainshock magnitude. These results are more consistent with the cascade model.

Journal of Seismology↗

Possible cause for an improbable earthquake: The 1997 MW 4.9 southern Alabama earthquake and hydrocarbon recovery

Circumstantial and physical evidence indicates that the 1997 MW 4.9 earthquake in southern Alabama may have been related to hydrocarbon recovery. Epicenters of this earthquake and its aftershocks were located within a few kilometers of active oil and gas extraction wells and two pressurized injection wells. Main shock and aftershock focal depths (2-6 km) are within a few kilometers of the injection and withdrawal depths. Strain accumulation at geologic rates sufficient to cause rupture at these shallow focal depths is not likely. A paucity of prior seismicity is difficult to reconcile with the occurrence of an earthquake of MW 4.9 and a magnitude-frequency relationship usually assumed for natural earthquakes. The normal-fault main-shock mechanism is consistent with reactivation of preexisting faults in the regional tectonic stress field. If the earthquake were purely tectonic, however, the question arises as to why it occurred on only the small fraction of a large, regional fault system coinciding with active hydrocarbon recovery. No obvious temporal correlation is apparent between the earthquakes and recovery activities. Although thus far little can be said quantitatively about the physical processes that may have caused the 1997 sequence, a plausible explanation involves the poroelastic response of the crust to extraction of hydrocarbons.

Geology↗

Earthquake triggering by transient and static deformations

Observational evidence for both static and transient near-field and far-field triggered seismicity are explained in terms of a frictional instability model, based on a single degree of freedom spring-slider system and rate- and state-dependent frictional constitutive equations. In this study a triggered earthquake is one whose failure time has been advanced by Δ t (clock advance) due to a stress perturbation. Triggering stress perturbations considered include square-wave transients and step functions, analogous to seismic waves and coseismic static stress changes, respectively. Perturbations are superimposed on a constant background stressing rate which represents the tectonic stressing rate. The normal stress is assumed to be constant. Approximate, closed-form solutions of the rate-and-state equations are derived for these triggering and background loads, building on the work of Dieterich [1992, 1994]. These solutions can be used to simulate the effects of static and transient stresses as a function of amplitude, onset time t 0 , and in the case of square waves, duration. The accuracies of the approximate closed-form solutions are also evaluated with respect to the full numerical solution and t 0 . The approximate solutions underpredict the full solutions, although the difference decreases as t 0 approaches the end of the earthquake cycle. The relationship between Δ t and t 0 differs for transient and static loads: a static stress step imposed late in the cycle causes less clock advance than an equal step imposed earlier, whereas a later applied transient causes greater clock advance than an equal one imposed earlier. For equal Δ t , transient amplitudes must be greater than static loads by factors of several tens to hundreds depending on t 0 . We show that the rate-and-state model requires that the total slip at failure is a constant, regardless of the loading history. Thus a static load applied early in the cycle, or a transient applied at any time, reduces the stress at the initiation of failure, whereas static loads that are applied sufficiently late raise it. Rate-and-state friction predictions differ markedly from those based on Coulomb failure stress changes (ΔCFS) in which Δ t equals the amplitude of the static stress change divided by the background stressing rate. The ΔCFS model assumes a stress failure threshold, while the rate-and-state equations require a slip failure threshold. The complete rate-and-state equations predict larger Δ t than the ΔCFS model does for static stress steps at small t 0 , and smaller Δ t than the ΔCFS model for stress steps at large t 0 . The ΔCFS model predicts nonzero Δ t only for transient loads that raise the stress to failure stress levels during the transient. In contrast, the rate-and-state model predicts nonzero Δ t for smaller loads, and triggered failure may occur well after the transient is finished. We consider heuristically the effects of triggering on a population of faults, as these effects might be evident in seismicity data. Triggering is manifest as an initial increase in seismicity rate that may be followed by a quiescence or by a return to the background rate. Available seismicity data are insufficient to discriminate whether triggered earthquakes are “new” or clock advanced. However, if triggering indeed results from advancing the failure time of inevitable earthquakes, then our modeling suggests that a quiescence always follows transient triggering and that the duration of increased seismicity also cannot exceed the duration of a triggering transient load. Quiescence follows static triggering only if the population of available faults is finite.

Journal of Geophysical Research B: Solid Earth↗

Dynamic deformations of shallow sediments in the Valley of Mexico, Part I: Three-dimensional strains and rotations recorded on a seismic array

We study the spatial variation in earthquake ground motions, or equivalently the dynamic displacement gradient field, using a novel analysis procedure borrowed from geodesy. Seismic data recorded in the Valley of Mexico by a microarray of three three-component surface accelerographs and two three-component accelerographs at depths of 30 m and 102 m constrain our estimates of the dynamic displacement gradient field (from which strains and rotations derive) for four moderate earthquakes at distances of 250 to 300 km. Our study focuses on the effects of low-velocity surface materials on the deformation. At the surface, the gradients corresponding to deformation across vertical planes dominate, and vertical-axis rotations are of similar magnitudes as strains. The greatest peak surface gradient we observed was 206 μ strain for the 14 September 1995 M W 7.5 earthquake at a distance of ∼300 km. However, much larger gradients occur across horizontal planes ( ∂u/∂z , where u is a horizontal displacement and z is depth) at some depth between 0 and 30 m. These values are about a factor 10 greater than the corresponding gradient components at the surface. ∂u/∂z for the 14 September earthquake equaled or exceeded 665 μ strain at depth. The dynamic deformations experienced in Mexico City undoubtedly have occurred before and will occur again in other densely populated areas. However, in many other regions, the sediment response will not remain linear and elastic, resulting instead in liquefaction and ground failure.

Valley of Mexico↗

Dynamic deformations of shallow sediments in the Valley of Mexico, Part II: Single-station estimates

We develop simple relations to estimate dynamic displacement gradients (and hence the strains and rotations) during earthquakes in the lake-bed zone of the Valley of Mexico, where the presence of low-velocity, high-water content clays in the uppermost layers cause dramatic amplification of seismic waves and large strains. The study uses results from a companion article (Bodin et al. , 1997) in which the data from an array at Roma, a lake-bed site, were analyzed to obtain displacement gradients. In this article, we find that the deformations at other lake-bed sites may differ from those at Roma by a factor of 2 to 3. More accurate estimates of the dominant components of the deformation at an individual instrumented lake-bed site may be obtained from the maximum horizontal velocity and displacement, ν max and u max , at the surface. The maximum surface strain ɛ max is related to ν max by ɛ max = ν max / C , with C ∼ 0.6 km/sec. From the analysis of data from sites equipped with surface and borehole sensors, we find that the vertical gradient of peak horizontal displacement (Δ u max /Δ z ) computed from sensors at 0 and 30 m equals ( u max ) z=0 /Δ z , Δ z = 30 m, within a factor of 1.5. This is the largest gradient component, and the latter simple relation permits its estimation from surface records alone. The observed profiles of u max versus depth suggest a larger gradient in some depth range of 10 to 20 m, in agreement with synthetic calculations presented in Bodin et al. (1997).

Valley of Mexico↗

Transient triggering of near and distant earthquakes

We demonstrate qualitatively that frictional instability theory provides a context for understanding how earthquakes may be triggered by transient loads associated with seismic waves from near and distance earthquakes. We assume that earthquake triggering is a stick-slip process and test two hypotheses about the effect of transients on the timing of instabilities using a simple spring-slider model and a rate- and state-dependent friction constitutive law. A critical triggering threshold is implicit in such a model formulation.

Bulletin of the Seismological Society of America↗

On the origin of diverse aftershock mechanisms following the 1989 Loma Prieta earthquake

We test the hypothesis that the origin at the diverse suite of aftershock mechanisms Following the 1989 M 7.1 Loma Prieta, California, earthquake is related to the post-main-shock static stress field We use a 3-D boundary-element algorithm to calculate static stresses combined with a Coulomb failure criterion to calculate conjugate failure planes at aftershock locations. The post-main-shock static stress field is taken as the of a sum of a pre-existing stress field and changes in stress due to the heterogeneous slip across the Loma Prieta rupture plane. The background stress field is assumed to be either a simple shear parallel to the regional trend of the San Andreas fault or approximately fault-normal compression A suite of synthetic aftershock mechanisms from the conjugate failure planes is generated and quantitatively compared (allowing for uncertainties in both mechanism parameters and earthquake locations) to well-constrained mechanisms reported in the US Geological Survey Northern California Seismic Network catalogue. We also compare calculated rakes with those observed by resolving the calculated stress tensor onto observed focal mechanism nodal planes, assuming either plane to be a likely rupture plane.

Geophysical Journal International↗

Stress/strain changes and triggered seismicity at The Geysers, California

The principal results of this study of remotely triggered seismicity in The Geysers geothermal field are the demonstration that triggering (initiation of earthquake failure) depends on a critical strain threshold and that the threshold level increases with decreasing frequency, or, equivalently, depends on strain rate. This threshold function derives from (1) analyses of dynamic strains associated with surface waves of the triggering earthquakes, (2) statistically measured aftershock zone dimensions, and (3) analytic functional representations of strains associated with power production and tides. The threshold is also consistent with triggering by static strain changes and implies that both static and dynamic strains may cause aftershocks. The observation that triggered seismicity probably occurs in addition to background activity also provides an important constraint on the triggering process. Assuming the physical processes underlying earthquake nucleation to be the same, Gomberg [this issue] discusses seismicity triggered by the M W 7.3 Landers earthquake, its constraints on the variability of triggering thresholds with site, and the implications of time delays between triggering and triggered earthquakes. Our results enable us to reject the hypothesis that dynamic strains simply nudge prestressed faults over a Coulomb failure threshold sooner than they would have otherwise. We interpret the rate-dependent triggering threshold as evidence of several competing processes with different time constants, the faster one(s) facilitating failure and the other(s) inhibiting it. Such competition is a common feature of theories of slip instability. All these results, not surprisingly, imply that to understand earthquake triggering one must consider not only simple failure criteria requiring exceedence of some constant threshold but also the requirements for generating instabilities.

Journal of Geophysical Research B: Solid Earth↗

Seismicity remotely triggered by the magnitude 7.3 Landers, California, earthquake

The magnitude 7.3 Landers earthquake of 28 June 1992 triggered a remarkably sudden and widespread increase in earthquake activity across much of the western United States. The triggered earthquakes, which occurred at distances up to 1250 kilometers (17 source dimensions) from the Landers mainshock, were confined to areas of persistent seismicity and strike-slip to normal faulting. Many of the triggered areas also are sites of geothermal and recent volcanic activity. Static stress changes calculated for elastic models of the earthquake appear to be too small to have caused the triggering. The most promising explanations involve nonlinear interactions between large dynamic strains accompanying seismic waves from the mainshock and crustal fluids (perhaps including crustal magma).

California↗