Geology Reports⌕ Search

Geology topics

J. Gomberg

Publications and source records attributed to J. Gomberg.

At least 19 recordsLinked to original sources
USGS2012

Seismicity and shear strain in the southern Great Basin of Nevada and California

This study examines the relationship between the distribution of small earthquakes (M L ≤4.3) and mechanisms of strain accumulation and relaxation in an area with long repeat times between large events, the Southern Great Basin Seismic Network (SGBSN) region. The Great Basin is a unique continental extensional province characterized by normal and strike-slip faulting, high heat flow, crust of thin to normal thickness, and high elevations. The SGBSN is operated to provide data to address suitability issues pertaining to Yucca Mountain, Nevada which is being evaluated as a potential site for a national mined geologic nuclear waste repository. Suitability issues include estimation of the probability of occurrence of future damaging earthquakes, the characterization of the mechanisms that drive hydrologic flow, and the identification of fractures (faults) that might act as flow conduits or barriers. This study attempts to explain the distribution of small earthquakes in terms of spatial variations in the shear strain field; where strain concentrates there should be a greater number of small earthquakes. Strain field models are constructed under the assumption that long term fault behavior perturbs an otherwise uniform strain field. These strain field models are then interpreted with regard to the regional tectonics and site suitability issues. Modeling results provide one possible explanation of why earthquake clusters cover regions much larger than the surface projections of any of mapped major faults; clusters in a wide band along and extending northeast of the northern half of the Furnace Creek fault may correspond to elevated shear strains along the fault and a broad cluster in the Pahranagat Shear Zone may be associated with shear strain arising from a distribution of smaller localized faults. The relatively large number of small earthquakes in the southern and eastern portions of the Nevada Test Site is consistent with the strain field models. A minimum in shear strain at Yucca Mountain is predicted by all models consistent with an almost total lack of earthquakes observed there. The region to the west of the Death Valley/Furnace Creek fault system, the portion of the study area with the most active deformation but few small earthquakes, is an area of low shear strain. A possible reason for this is that the fault configuration in the area is optimal for accommodating regional deformation via large earthquakes or creep. While there is also a relative lack of earthquakes at Yucca Mountain, this may be indicative of a lack of accumulating strain energy and thus, a lower potential for a large earthquake.

Journal of Geophysical Research↗
USGS2012

Seismicity and detection/location threshold in the southern Great Basin seismic network

A spatially varying model of the detection/location capabilities of the Southern Great Basin seismic network (SGBSN) has been derived that is based on simple empirical relations and statistics. This permits use of almost all the catalog data gathered; instead of ignoring data that are below the threshold of completeness, a spatially varying threshold model is developed so that subregions having lower completeness levels than the network as a whole can be outlined and the completeness level of each sub-region determined. Such a model is required to unambiguously identify regions that are aseismic due to natural processes rather than to limited detection and/or location capabilities. Accounting for spatial variations in detection/location threshold is also important for studies in which magnitude-frequency distributions are interpreted in terms of source scaling properties. The characteristics of the spatial distribution of earthquakes, where earthquake clusters and aseismic regions locate, appear to be stable at all magnitude levels so that inferences about where strain is being accommodated will be the same whether numbers of earthquakes or strain estimated from seismic moments are examined. For the southern Great Basin region these principal characteristics include clusters at the northern end of the Furnace Creek fault and in the Pahranagat Shear Zone, and a relatively large number of earthquakes in the northern and southeastern portions of the Nevada Test Site. These clusters cover regions much larger than the surface projections of any of the mapped faults. The extent to which seismicity is induced by nuclear testing is unclear. The predominantly aseismic regions include the area west of the Death Valley/Furnace Creek fault system and an almost complete absence of events at Yucca Mountain. Finally, a considerable number of isolated events in the SGBSN catalog cannot be correlated with mapped faults.

Journal of Geophysical Research↗
USGS2012

Contemporary seismicity in and around the Yakima-Fold-and-Thrust Belt in eastern Washington

We examined characteristics of routinely cataloged seismicity from 1970 to the present in and around the Yakima fold‐and‐thrust belt (YFTB) in eastern Washington to determine if the characteristics of contemporary seismicity provide clues about regional‐scale active tectonics or about more localized, near‐surface processes. We employed new structural and hydrologic models of the Columbia River basalts (CRB) and found that one‐third to one‐half of the cataloged earthquakes occur within the CRB and that these CRB earthquakes exhibit significantly more clustered, and swarmlike, behavior than those outside. These results and inferences from published studies led us to hypothesize that clustered seismicity is likely associated with hydrologic changes in the CRB, which hosts the regional aquifer system. While some general features of the regional groundwater system support this hypothesis, seismicity patterns and mapped long‐term changes in groundwater levels and present‐day irrigation neither support nor refute it. Regional tectonic processes and crustal‐scale structures likely influence the distribution of earthquakes both outside and within the CRB as well. We based this inference on qualitatively assessed alignments between the dominant northwest trends in the geologic structure and the seismicity generally and between specific faults and characteristics of the 2009 Wooded Island swarm and aseismic slip, which is the only cluster studied in detail and the most vigorous since regional monitoring began.

Washington↗
USGS2011

InSAR observations of aseismic slip associated with an earthquake swarm in the Columbia River flood basalts

In 2009 a swarm of small shallow earthquakes occurred within the basalt flows of the Columbia River Basalt Group (CRBG). The swarm occurred within a dense seismic network in the U.S. Department of Energys Hanford Site. Data from the seismic network along with interferometric synthetic aperture radar (InSAR) data from the European Space Agencys (ESA) ENVISAT satellite provide insight into the nature of the swarm. By modeling the InSAR deformation data we constructed a model that consists of a shallow thrust fault and a near horizontal fault. We suggest that the near horizontal lying fault is a bedding-plane fault located between basalt flows. The geodetic moment of the modeled fault system is about eight times the cumulative seismic moment of the swarm. Precise location estimates of the swarm earthquakes indicate that the area of highest slip on the thrust fault, ???70mm of slip less than ???0.5km depth, was not located within the swarm cluster. Most of the slip on the faults appears to have progressed aseismically and we suggest that interbed sediments play a central role in the slip process. Copyright 2011 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2011

Seismic and geodetic signatures of fault slip at the Slumgullion Landslide Natural Laboratory

We tested the hypothesis that the Slumgullion landslide is a useful natural laboratory for observing fault slip, specifically that slip along its basal surface and side-bounding strike-slip faults occurs with comparable richness of aseismic and seismic modes as along crustal- and plate-scale boundaries. Our study provides new constraints on models governing landslide motion. We monitored landslide deformation with temporary deployments of a 29-element prism array surveyed by a robotic theodolite and an 88-station seismic network that complemented permanent extensometers and environmental instrumentation. Aseismic deformation observations show that large blocks of the landslide move steadily at approximately centimeters per day, possibly punctuated by variations of a few millimeters, while localized transient slip episodes of blocks less than a few tens of meters across occur frequently. We recorded a rich variety of seismic signals, nearly all of which originated outside the monitoring network boundaries or from the side-bounding strike-slip faults. The landslide basal surface beneath our seismic network likely slipped almost completely aseismically. Our results provide independent corroboration of previous inferences that dilatant strengthening along sections of the side-bounding strike-slip faults controls the overall landslide motion, acting as seismically radiating brakes that limit acceleration of the aseismically slipping basal surface. Dilatant strengthening has also been invoked in recent models of transient slip and tremor sources along crustal- and plate-scale faults suggesting that the landslide may indeed be a useful natural laboratory for testing predictions of specific mechanisms that control fault slip at all scales.

Colorado↗
USGS2011

Tiny intraplate earthquakes triggered by nearby episodic tremor and slip in Cascadia

Episodic tremor and slip (ETS) has been observed in many subduction zones, but its mechanical underpinnings as well as its potential for triggering damaging earthquakes have proven difficult to assess. Here we use a seismic array in Cascadia of unprecedented density to monitor seismicity around a moderate 16 day ETS episode. In the 4 months of data we examine, we observe five tiny earthquakes within the subducting slab during the episode and only one more in the same area, which was just before and nearby the next ETS burst. These earthquakes concentrate along the sides and updip edge of the ETS region, consistent with greater stress concentration there than near the middle and downdip edge of the tremor area. Most of the seismicity is below the megathrust, with a similar depth extent to the background intraslab seismicity. The pattern of earthquakes that we find suggests slow slip has a more continuous temporal and spatial pattern than the tremor loci, which notoriously appear in bursts, jumps, and streaks. Copyright 2011 by the American Geophysical Union.

Geochemistry, Geophysics, Geosystems↗
USGS2009

Inducing in situ, nonlinear soil response applying an active source

[1] It is well known that soil sites have a profound effect on ground motion during large earthquakes. The complex structure of soil deposits and the highly nonlinear constitutive behavior of soils largely control nonlinear site response at soil sites. Measurements of nonlinear soil response under natural conditions are critical to advancing our understanding of soil behavior during earthquakes. Many factors limit the use of earthquake observations to estimate nonlinear site response such that quantitative characterization of nonlinear behavior relies almost exclusively on laboratory experiments and modeling of wave propagation. Here we introduce a new method for in situ characterization of the nonlinear behavior of a natural soil formation using measurements obtained immediately adjacent to a large vibrator source. To our knowledge, we are the first group to propose and test such an approach. Employing a large, surface vibrator as a source, we measure the nonlinear behavior of the soil by incrementally increasing the source amplitude over a range of frequencies and monitoring changes in the output spectra. We apply a homodyne algorithm for measuring spectral amplitudes, which provides robust signal-to-noise ratios at the frequencies of interest. Spectral ratios are computed between the receivers and the source as well as receiver pairs located in an array adjacent to the source, providing the means to separate source and near-source nonlinearity from pervasive nonlinearity in the soil column. We find clear evidence of nonlinearity in significant decreases in the frequency of peak spectral ratios, corresponding to material softening with amplitude, observed across the array as the source amplitude is increased. The observed peak shifts are consistent with laboratory measurements of soil nonlinearity. Our results provide constraints for future numerical modeling studies of strong ground motion during earthquakes.

Journal of Geophysical Research B: Solid Earth↗
USGS2009

A search in strainmeter data for slow slip associated with triggered and ambient tremor near Parkfield, California

We test the hypothesis that, as in subduction zones, slow slip facilitates triggered and ambient tremor in the transform boundary setting of California. Our study builds on the study of Peng et al. (2009) of triggered and ambient tremor near Parkfield, California during time intervals surrounding 31, potentially triggering, M ≥ 7.5 teleseismic earthquakes; waves from 10 of these triggered tremor and 29 occurred in periods of ambient tremor activity. We look for transient slow slip during 3-month windows that include 11 of these triggering and nontriggering teleseisms, using continuous strain data recorded on two borehole Gladwin tensor strainmeters (GTSM) located within the distribution of tremor epicenters. We model the GTSM data assuming only tidal and “drift” signals are present and find no detectable slow slip, either ongoing when the teleseismic waves passed or triggered by them. We infer a conservative detection threshold of about 5 nanostrain for abrupt changes and about twice this for slowly evolving signals. This could be lowered slightly by adding analyses of other data types, modeled slow slip signals, and GTSM data calibration. Detection of slow slip also depends on the slipping fault's location and size, which we describe in terms of equivalent earthquake moment magnitude, M. In the best case of the GTSM above a very shallow slipping fault, detectable slip events must exceed M~2, and if the slow slip is beneath the seismogenic zone (below ~15 km depth), even M~5 events are likely to remain hidden.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

A model of earthquake triggering probabilities and application to dynamic deformations constrained by ground motion observations

We have used observations from Felzer and Brodsky (2006) of the variation of linear aftershock densities (i.e., aftershocks per unit length) with the magnitude of and distance from the main shock fault to derive constraints on how the probability of a main shock triggering a single aftershock at a point, P(r, D), varies as a function of distance, r, and main shock rupture dimension, D. We find that P(r, D) becomes independent of D as the triggering fault is approached. When r ??? D P(r, D) scales as Dm where m-2 and decays with distance approximately as r-n with n = 2, with a possible change to r-(n-1) at r > h, where h is the closest distance between the fault and the boundaries of the seismogenic zone. These constraints may be used to test hypotheses about the types of deformations and mechanisms that trigger aftershocks. We illustrate this using dynamic deformations (i.e., radiated seismic waves) and a posited proportionality with P(r, D). Deformation characteristics examined include peak displacements, peak accelerations and velocities (proportional to strain rates and strains, respectively), and two measures that account for cumulative deformations. Our model indicates that either peak strains alone or strain rates averaged over the duration of rupture may be responsible for aftershock triggering.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Effects of acoustic waves on stick-slip in granular media and implications for earthquakes

It remains unknown how the small strains induced by seismic waves can trigger earthquakes at large distances, in some cases thousands of kilometres from the triggering earthquake, with failure often occurring long after the waves have passed. Earthquake nucleation is usually observed to take place at depths of 10-20 km, and so static overburden should be large enough to inhibit triggering by seismic-wave stress perturbations. To understand the physics of dynamic triggering better, as well as the influence of dynamic stressing on earthquake recurrence, we have conducted laboratory studies of stick-slip in granular media with and without applied acoustic vibration. Glass beads were used to simulate granular fault zone material, sheared under constant normal stress, and subject to transient or continuous perturbation by acoustic waves. Here we show that small-magnitude failure events, corresponding to triggered aftershocks, occur when applied sound-wave amplitudes exceed several microstrain. These events are frequently delayed or occur as part of a cascade of small events. Vibrations also cause large slip events to be disrupted in time relative to those without wave perturbation. The effects are observed for many large-event cycles after vibrations cease, indicating a strain memory in the granular material. Dynamic stressing of tectonic faults may play a similar role in determining the complexity of earthquake recurrence. ??2007 Nature Publishing Group.

Nature↗
USGS2008

Induced dynamic nonlinear ground response at Gamer Valley, California

We present results from a prototype experiment in which we actively induce, observe, and quantify in situ nonlinear sediment response in the near surface. This experiment was part of a suite of experiments conducted during August 2004 in Garner Valley, California, using a large mobile shaker truck from the Network for Earthquake Engineering Simulation (NEES) facility. We deployed a dense accelerometer array within meters of the mobile shaker truck to replicate a controlled, laboratory-style soil dynamics experiment in order to observe wave-amplitude-dependent sediment properties. Ground motion exceeding 1g acceleration was produced near the shaker truck. The wave field was dominated by Rayleigh surface waves and ground motions were strong enough to produce observable nonlinear changes in wave velocity. We found that as the force load of the shaker increased, the Rayleigh-wave phase velocity decreased by as much as ???30% at the highest frequencies used (up to 30 Hz). Phase velocity dispersion curves were inverted for S-wave velocity as a function of depth using a simple isotropic elastic model to estimate the depth dependence of changes to the velocity structure. The greatest change in velocity occurred nearest the surface, within the upper 4 m. These estimated S-wave velocity values were used with estimates of surface strain to compare with laboratory-based shear modulus reduction measurements from the same site. Our results suggest that it may be possible to characterize nonlinear soil properties in situ using a noninvasive field technique.

Bulletin of the Seismological Society of America↗
USGS2007

Non-volcanic tremor driven by large transient shear stresses

Non-impulsive seismic radiation or 'tremor' has long been observed at volcanoes and more recently around subduction zones. Although the number of observations of non-volcanic tremor is steadily increasing, the causative mechanism remains unclear. Some have attributed non-volcanic tremor to the movement of fluids, while its coincidence with geodetically observed slow-slip events at regular intervals has led others to consider slip on the plate interface as its cause. Low-frequency earthquakes in Japan, which are believed to make up at least part of non-volcanic tremor, have focal mechanisms and locations that are consistent with tremor being generated by shear slip on the subduction interface. In Cascadia, however, tremor locations appear to be more distributed in depth than in Japan, making them harder to reconcile with a plate interface shear-slip model. Here we identify bursts of tremor that radiated from the Cascadia subduction zone near Vancouver Island, Canada, during the strongest shaking from the moment magnitude Mw = 7.8, 2002 Denali, Alaska, earthquake. Tremor occurs when the Love wave displacements are to the southwest (the direction of plate convergence of the overriding plate), implying that the Love waves trigger the tremor. We show that these displacements correspond to shear stresses of approximately 40 kPa on the plate interface, which suggests that the effective stress on the plate interface is very low. These observations indicate that tremor and possibly slow slip can be instantaneously induced by shear stress increases on the subduction interface - effectively a frictional failure response to the driving stress. ??2007 Nature Publishing Group.

Nature↗
USGS2005

Time-dependent earthquake probabilities

We have attempted to provide a careful examination of a class of approaches for estimating the conditional probability of failure of a single large earthquake, particularly approaches that account for static stress perturbations to tectonic loading as in the approaches of Stein et al. (1997) and Hardebeck (2004). We have loading as in the framework based on a simple, generalized rate change formulation and applied it to these two approaches to show how they relate to one another. We also have attempted to show the connection between models of seismicity rate changes applied to (1) populations of independent faults as in background and aftershock seismicity and (2) changes in estimates of the conditional probability of failures of different members of a the notion of failure rate corresponds to successive failures of different members of a population of faults. The latter application requires specification of some probability distribution (density function of PDF) that describes some population of potential recurrence times. This PDF may reflect our imperfect knowledge of when past earthquakes have occurred on a fault (epistemic uncertainty), the true natural variability in failure times, or some combination of both. We suggest two end-member conceptual single-fault models that may explain natural variability in recurrence times and suggest how they might be distinguished observationally. When viewed deterministically, these single-fault patch models differ significantly in their physical attributes, and when faults are immature, they differ in their responses to stress perturbations. Estimates of conditional failure probabilities effectively integrate over a range of possible deterministic fault models, usually with ranges that correspond to mature faults. Thus conditional failure probability estimates usually should not differ significantly for these models. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2005

A frictional population model of seismicity rate change

We study models of seismicity rate changes caused by the application of a static stress perturbation to a population of faults and discuss our results with respect to the model proposed by Dieterich (1994). These models assume distribution of nucleation sites (e.g., faults) obeying rate-state frictional relations that fail at constant rate under tectonic loading alone, and predicts a positive static stress step at time to will cause an immediate increased seismicity rate that decays according to Omori's law. We show one way in which the Dieterich model may be constructed from simple general idead, illustratted using numerically computed synthetic seismicity and mathematical formulation. We show that seismicity rate change predicted by these models (1) depend on the particular relationship between the clock-advanced failure and fault maturity, (2) are largest for the faults closest to failure at to, (3) depend strongly on which state evolution law faults obey, and (4) are insensitive to some types of population hetrogeneity. We also find that if individual faults fail repeatedly and populations are finite, at timescales much longer than typical aftershock durations, quiescence follows at seismicity rate increase regardless of the specific frictional relations. For the examined models the quiescence duration is comparable to the ratio of stress change to stressing rate ????/??,which occurs after a time comparable to the average recurrence interval of the individual faults in the population and repeats in the absence of any new load may pertubations; this simple model may partly explain observations of repeated clustering of earthquakes. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2004

Earthquake nucleation by transient deformations caused by the M = 7.9 Denali, Alaska, earthquake

The permanent and dynamic (transient) stress changes inferred to trigger earthquakes are usually orders of magnitude smaller than the stresses relaxed by the earthquakes themselves, implying that triggering occurs on critically stressed faults. Triggered seismicity rate increases may therefore be most likely to occur in areas where loading rates are highest and elevated pore pressures, perhaps facilitated by high-temperature fluids, reduce frictional stresses and promote failure. Here we show that the 2002 magnitude M = 7.9 Denali, Alaska, earthquake triggered wide-spread seismicity rate increases throughout British Columbia and into the western United States. Dynamic triggering by seismic waves should be enhanced in directions where rupture directivity focuses radiated energy, and we verify this using seismic and new high-sample GPS recordings of the Denali mainshock. These observations are comparable in scale only to the triggering caused by the 1992 M = 7.4 Landers, California, earthquake, and demonstrate that Landers triggering did not reflect some peculiarity of the region or the earthquake. However, the rate increases triggered by the Denali earthquake occurred in areas not obviously tectonically active, implying that even in areas of low ambient stressing rates, faults may still be critically stressed and that dynamic triggering may be ubiquitous and unpredictable.

Nature↗
USGS2003

Using 1-Hz GPS data to measure deformations caused by the denali fault earthquake

The 3 November 2002 moment magnitude 7.9 Denali fault earthquake generated large, permanent surface displacements in Alaska and large-amplitude surface waves throughout western North America. We find good agreement between strong ground-motion records integrated to displacement and 1-hertz Global Positioning System (GPS) position estimates collected ??? 140 kilometers from the earthquake epicenter. One-hertz GPS receivers also detected seismic surface waves 750 to 3800 kilometers from the epicenter, whereas these waves saturated many of the seismic instruments in the same region. High-frequency GPS increases the dynamic range and frequency bandwidth of ground-motion observations, providing another tool for studying earthquake processes.

Science↗