Geology topics

J. Gomberg

Publications and source records attributed to J. Gomberg.

27 records · Page 2Linked to original sources

USGS1999

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↗

USGS1999

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↗

USGS1998

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↗

USGS1997

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↗

USGS1997

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↗

USGS1997

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↗

USGS1997

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↗

USGS1996

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↗

USGS1993

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↗