Search USGS⌕ Search

Geology topics

Joan S. Gomberg

Publications and source records attributed to Joan S. Gomberg.

At least 19 recordsLinked to original sources

Rupture continuity through intermittent pauses in Cascadia slow slip events

Cascadia slow slip events (SSEs) are often envisioned as smooth, continuous ruptures, progressively activating tremor asperities as they propagate. Macroscopically, geodetic inversions and spatiotemporal maps of tremor epicenters show steady, uniform migration. In detail tremor is more chaotic and discontinuous. Larger long-term SSEs observed in daily geodetic solutions are inferred to exhibit intermittent pauses that reflect temporary re-locking of the fault, but this temporal resolution limits tests for similar re-locking on shorter timescales. We use temporal measurements of the areal growth and radiated energy of tremor clusters to investigate SSE intermittence. We find that ruptures mirror tremor pauses. Areal growth rate, however, does not reset, and removing the pauses results in smoother and more similar growth measurements among all SSEs. The rupture similarity occurs regardless of size or location and hints at an underlying uniformity and lack of predeterminism in eventual SSE size. Epicentral uncertainty precludes quantifying early rupture stages, but for larger events areal growth follows a power-law and slows with increasing size. Temporal correlations in tremor energy with inferred SSE propagation velocities and tremor rates suggest its use as a proxy for slip velocity. We find that tremor energy is tidally modulated at daily and sub-daily frequencies, and this modulation is continuous through pauses, suggesting a memory of slip state is sustained through them. We argue these pauses reflect unsteady propagation of the slip front, marked by rapid re- and un-locking, and excluding them removes rupture complexity to reveal a diffusive-like slip process and underlying universality in growth.

California, Oregon, Washington↗

Submarine avalanche deposits hold clues to past earthquakes

Earthquakes and other natural events sometimes shake the seafloor near coastlines severely enough to cause underwater avalanches that rush down steep slopes, scouring the seabed and carrying sediment to greater depths. These fast-moving sediment-laden flows, called turbidity currents , have at times damaged underwater infrastructure like pipelines and communications cables, as they did, for example, in snapping transatlantic cables off the coast of Newfoundland after the 1929 Grand Banks earthquake.

EOS Transactions↗

Turbidite correlation for paleoseismology

Marine turbidite paleoseismology relies on the assumption of synchronous triggering of turbidity currents by earthquake shaking to infer rupture extent and recurrence. Such inference commonly depends on age dating and correlation of the physical stratigraphy of deposits carried by turbidity currents (i.e., turbidites) across great distances. Along the Cascadia subduction zone, which lies offshore the Pacific Northwest, USA, turbidite facies in core photographs, X-ray computed tomography images, and magnetic susceptibility (MS) data exhibit differences in character over relatively short distances, which implies that not all deposits can be correlated with confidence. Thus, subjective correlation based on expected similarity over great distances and weak age constraints does not independently support paleoseismic models. We present a new method for correlating turbidites along the Cascadia margin that can yield a more objective and repeatable stratigraphic framework to underpin earthquake recurrence. We use dynamic time warping to correlate MS logs and measure correlation coefficients of core pairs to evaluate correlation strength. We then compare these measures to a distribution of correlation coefficients of randomly generated turbidite sequences and find that only a small number of core pairs can be correlated more confidently than randomly stacked turbidites. This methodology promises a more robust correlation strategy for future stratigraphic studies.

Oregon, Washington↗

Laboratory hydrofractures as analogs to tectonic tremors

The fracture of Earth materials occurs over a wide range of time and length scales. Physical conditions, particularly the stress field and Earth material properties, may condition rupture in a specific fracture regime. In nature, fast and slow fractures occur concurrently: tectonic tremor events are fast enough to emit seismic waves and frequently accompany slow earthquakes, which are too slow to emit seismic waves and are referred to as aseismic slip events. In this study, we generate simultaneous seismic and aseismic processes in a laboratory setting by driving a penny-shaped crack in a transparent sample with pressurized fluid. We leverage synchronized high-speed imaging and high-frequency acoustic emission (AE) sensing to visualize and listen to the various sequences of propagation (breaks) and arrest (sticks) of a fracture undergoing stick-break instabilities. Slow radial crack propagation is facilitated by fast tangential fractures. Fluid viscosity and pressure regulate the fracture dynamics of slow and fast events, and control the inter-event time and the energy released during individual fast events. These AE signals share behaviors with observations of episodic tremors in Cascadia, United States; these include: (a) bursty or intermittent slow propagation, and (b) nearly linear scaling of radiated energy with area. Our laboratory experiments provide a plausible model of tectonic tremor as an indicative of hydraulic fracturing facilitating shear slip during slow earthquakes.

AGU Advances↗

Slow slip detectability in seafloor pressure records offshore Alaska

In subduction zones worldwide, seafloor pressure data are used to observe tectonic deformation, particularly from megathrust earthquakes and slow slip events (SSEs). However, such measurements are also sensitive to oceanographic circulation-generated pressures over a range of frequencies that conflate with tectonic signals of interest. Using seafloor pressure and temperature data from the Alaska Amphibious Community Seismic Experiment, and sea surface height data from satellite altimetry, we evaluate the efficacy of various seasonal and oceanographic pressure signal proxy corrections and conduct synthetic tests to determine their impact on the timing and amplitude prediction of ramp-like signals typical of SSEs. We find that subtracting out the first mode of the complex empirical orthogonal functions of the pressure records on either the shelf or slope yields signal root-mean-square error (RMS) reductions up to 73% or 80%, respectively. Additional correction with proxies that exploit the depth-dependent spatial coherence of pressure records provides cumulative variance reductions up to 83% and 93%, respectively. Our detectability tests show that the timing and amplitude of synthetic SSE-like ramps can be well constrained for ramp amplitudes ≥4 cm on the shelf and ≥2 cm on the slope, using a fully automated detector. The principal limits on detectability are residual abrupt changes in pressure that occur as part of the transition to and from summer to winter conditions but are not adequately characterized by our seasonal corrections, as well as the inability to properly account for instrumental drift, which is not readily separated from the seasonal signal.

Alaska↗

Working together to advance subduction zone science within SZ4D and the USGS

OnSunday, December 10, 2023, the U.S. Geological Survey (USGS) and Subduction Zones in Four Dimensions (SZ4D) co-hosted a workshop titled “Working Together to Advance Subduction Zone Science within SZ4D and the USGS”. The workshop attracted ~80 participants, with ~11 attendees from the USGS and ~17 participants from outside the United States Organizations represented included the Cascadia Region Earthquake Science Center (CRESCENT), the Cascadia Coastlines and Peoples (CoPes) Hub, and the Community Network for Volcanic Eruption Response (CONVERSE). Goals for this workshop included: ○ Discuss ideas and potential topics for SZ4D-USGS collaborations ○ Identify mechanisms and opportunities for potential collaboration ○ Identify barriers to collaboration (with possible solutions) ○ Compile topics addressed in the workshop and provide this information to SZ4D and USGS Leadership and the broader SZ4D and USGS Communities for their use.

Conference Paper↗

The finicky nature of earthquake shaking-triggered submarine sediment slope failures and sediment gravity flows

Since 2011, seafloor temperatures, pressures, and seismic ground motions have been measured by the seafloor cabled Dense Oceanfloor Network system for Earthquakes and Tsunamis (DONET) on the Nankai margin. These measurements, high-resolution bathymetry, and abundant contextual information make the DONET region seem ideally suited to provide constraints on seismic shaking-triggered sediment slope failures and gravity flows, particularly since numerous published studies have linked paleo- to modern earthquakes to failures and flows within the DONET. The occurrences of the local 2016 M6.0 Mie-ken and regional M7.0 Kumamoto earthquakes within and at regional distances, respectively, from the DONET data set provided an opportunity to explore this potential. We used DONET seismic recordings of the posited triggering shaking and to search for submarine slide signals and continuous temperature and pressure data to detect pulses of warm and densified water indicative of passing flows. We developed and applied a variety of analytical methods to eliminate signals generated by water column processes, while leaving slope failures and sediment gravity flow anomalies as residuals. Our explorations yielded no evidence that earthquake shaking initiated either phenomenon, which we suggest reflects the finicky nature both of the detection of and the physical processes that contribute to slope failures and flows (i.e., both require satisfying precise suites of conditions). Nonetheless, this negative result, our analyses, and the estimates of physical properties we derived for them, provide useful lessons and inputs for future studies.

Journal of Geophysical Research↗

The productivity of Cascadia aftershock sequences

This study addresses questions about the productivity of Cascadia mainshock–aftershock sequences using earthquake catalogs produced by the Geological Survey of Canada and the Pacific Northwest Seismic Network. Questions concern the likelihood that future moderate to large intermediate depth intraslab earthquakes in Cascadia would have as few detectable aftershocks as those documented since 1949. More broadly, for Cascadia, we consider if aftershock productivities vary spatially, if they are outliers among global subduction zones, and if they are consistent with a physical model in which aftershocks are clock‐advanced versions of tectonically driven background seismicity. A practical motivation for this study is to assess the likely accuracy of aftershock forecasts based on productivities derived from global data that are now being issued routinely by the U.S. Geological Survey. For this reason, we estimated productivity following the identical procedures used in those forecasts and described in Page et al. (2016) . Results indicate that in Cascadia we can say that the next intermediate depth intraslab earthquake will likely have just a few detectable aftershocks and that aftershock productivity appears to be an outlier among global subduction zones, with rates that on average are lower by more than half, except for mainshocks in the upper plate. Our results are consistent with a clock‐advance model; productivities may be related to the proximity of mainshocks to a population of seismogenic fault patches and correlate with background seismicity rates. The latter and a clear correlation between productivities with mainshock depth indicate that both factors may have predictive value for aftershock forecasting.

British Columbia, California, Oregon, Washington↗

Toward an integrative geological and geophysical view of Cascadia subduction zone earthquakes

The Cascadia subduction zone (CSZ) is an exceptional geologic environment for recording evidence of land level changes, tsunamis, and ground motion that reveals at least 19 great megathrust earthquakes over the past 10 kyr. Such earthquakes are among the most impactful natural hazards on Earth, transcend national boundaries, and can have global impact. Reducing the societal impacts of future events in the U.S. Pacific Northwest and coastal British Columbia, Canada requires improved scientific understanding of megathrust earthquake rupture, recurrence, and corresponding hazards. Despite substantial knowledge gained from decades of research, large uncertainties remain about the characteristics and frequencies of past CSZ earthquakes. In this review, we summarize geological, geophysical, and instrumental evidence relevant to understanding megathrust earthquakes along the CSZ and associated uncertainties. We discuss how the evidence constrains various models of great megathrust earthquake recurrence in Cascadia and identify potential paths forward for the earthquake science community.

British Columbia, Oregon, Washington↗

The ocean's impact on slow slip events

We test the hypothesis that ocean seafloor pressures impart stresses that alter the initiation or termination of transient slow slip events (SSEs) on shallow submarine and near-coastal faults, using simulated seafloor pressures and a new catalog of SSEs in the Hikurangi subduction zone. We show that seafloor pressures may be represented by an average time history over the ~100-km dimensions of the study area. We account for SSE uncertainties and the multiplicity of processes that affect SSEs statistically by estimating the probabilities of rejecting the null hypothesis that SSE initiation or termination pressures are those to be expected by chance sampling of known (modeled) seafloor pressures, with low probabilities indicating some causal connection. No impact of ocean pressure changes on SSE initiation is detectable, but a correlation with their terminations is suggested. SSE slip that weakens the fault and makes it more sensitive to small stress changes may explain results.

Pacific Ocean↗

Separating sea and slow slip signals on the seafloor

Seafloor pressure measurements hold promise for estimating vertical displacements from transient slow slip events on submarine faults. We assess the accuracy of pressure offset estimates that evolve over days to weeks and the confidence with which they may be attributed to tectonic deformation or to the ocean water column. One common approach to resolve this ambiguity assumes water column pressures vary insignificantly over the study region and are represented by stable reference site pressures. Assessing the validity of this assumption requires independent evidence. Correlations between pressures and colocated temperatures collected during the Hikurangi Ocean Bottom Investigation of Tremor and Slow Slip experiment suggest temperatures might provide a useful independent proxy for water column pressures. We compared offsets estimated using several methods, with temperature and other proxies. The use of a temperature proxy was unsuccessful, because seafloor temperatures did not track the seasonal signal that contributes significantly to seafloor pressure changes over the slow slip event period. Regardless of the estimation method, offsets varied within a few cm around some uncertain reference level. Commonly used statistical measures are shown not to be reliable indicators of offset accuracy since offsets contribute minimally to the total variance. Offsets estimated using identical methods but with seafloor pressures simulated using a regional ocean model were larger than those derived from the data but had a similar pattern. Since the model simulates only water column processes, this suggests a significant fraction of the estimated pressure offsets are due to seasonal water column signal and are not of tectonic origin.

Journal of Geophysical Research- Solid Earth↗

Earthquake arrival association with backprojection and graph theory

The association of seismic‐wave arrivals with causative earthquakes becomes progressively more challenging as arrival detection methods become more sensitive, and particularly when earthquake rates are high. For instance, seismic waves arriving across a monitoring network from several sources may overlap in time, false arrivals may be detected, and some arrivals may be of unknown phase (e.g., P or S waves). We propose an automated method to associate arrivals with earthquake sources and obtain source locations applicable to such situations. To do so, we use a pattern detection metric based on the principle of backprojection to reveal candidate sources followed by graph‐theory‐based clustering and an integer linear optimization routine to associate arrivals with the minimum number of sources necessary to explain the data. This method solves for all sources and phase assignments simultaneously, rather than in a sequential greedy procedure as is common in other association routines. We demonstrate our method on both synthetic and real data from the Integrated Plate Boundary Observatory Chile seismic network of northern Chile. For the synthetic tests, we report results for cases with varying complexity, including rates of 500 earthquakes/day and 500 false arrivals/station/day, for which we measure true positive detection accuracy of >95% . For the real data, we develop a new catalog between 1 January 2010 and 31 December 2017 containing 817,548 earthquakes, with detection rates on average 279 earthquakes/day and a magnitude‐of‐completion of M∼1.8 ⁠ . A subset of detections are identified as sources related to quarry and industrial site activity, and we also detect thousands of foreshocks and aftershocks of the 1 April 2014 M w 8.2 Iquique earthquake. During the highest rate of aftershock activity, >600 earthquakes/day are detected in the vicinity of the Iquique earthquake rupture zone.

Bulletin of the Seismological Society of America↗

Reducing risk where tectonic plates collide—U.S. Geological Survey subduction zone science plan

The U.S. Geological Survey (USGS) serves the Nation by providing reliable scientific information and tools to build resilience in communities exposed to subduction zone earthquakes, tsunamis, landslides, and volcanic eruptions. Improving the application of USGS science to successfully reduce risk from these events relies on whole community efforts, with continuing partnerships among scientists and stakeholders, including researchers from universities, other government labs and private industry, land-use planners, engineers, policy-makers, emergency managers and responders, business owners, insurance providers, the media, and the general public. Motivated by recent technological advances and increased awareness of our growing vulnerability to subduction-zone hazards, the USGS is uniquely positioned to take a major step forward in the science it conducts and products it provides, building on its tradition of using long-term monitoring and research to develop effective products for hazard mitigation. This science plan provides a blueprint both for prioritizing USGS science activities and for delineating USGS interests and potential participation in subduction zone science supported by its partners. The activities in this plan address many USGS stakeholder needs: High-fidelity tools and user-tailored information that facilitate increasingly more targeted, neighborhood-scale decisions to mitigate risks more cost-effectively and ensure post-event operability. Such tools may include maps, tables, and simulated earthquake ground-motion records conveying shaking intensity and frequency. These facilitate the prioritization of retrofitting of vulnerable infrastructure; Information to guide local land-use and response planning to minimize development in likely hazardous zones (for example, databases, maps, and scenario documents to guide evacuation route planning in communities near volcanoes, along coastlines vulnerable to tsunamis, and built on landslide-prone terrain); New tools to assess the potential for cascading hazards, such as landslides, tsunamis, coastal changes, and flooding caused by earthquakes or volcanic eruptions; Geospatial models of permanent, widespread land- and sea-level changes that may occur in the immediate aftermath of great ( M ≥8.0) subduction zone earthquakes; Strong partnerships between scientists and public safety providers for effective decision making during periods of elevated hazard and risk; Accurate forecasts of far-reaching hazards (for example, ash clouds, tsunamis) to avert catastrophes and unnecessary disruptions in air and sea transportation; Aftershock forecasts to guide decisions about when and where to re-enter, repair, or rebuild buildings and infrastructure, for all types of subduction zone earthquakes.

Circular↗

Reducing risk where tectonic plates collide

Most of the world’s earthquakes, tsunamis, landslides, and volcanic eruptions are caused by the continuous motions of the many tectonic plates that make up the Earth’s outer shell. The most powerful of these natural hazards occur in subduction zones, where two plates collide and one is thrust beneath another. The U.S. Geological Survey’s (USGS) “Reducing Risk Where Tectonic Plates Collide—A USGS Plan to Advance Subduction Zone Science” is a blueprint for building the crucial scientific foundation needed to inform the policies and practices that can make our Nation more resilient to subduction zone-related hazards.

Fact Sheet↗

Sediment gravity flows triggered by remotely generated earthquake waves

Recent great earthquakes and tsunamis around the world have heightened awareness of the inevitability of similar events occurring within the Cascadia Subduction Zone of the Pacific Northwest. We analyzed seafloor temperature, pressure, and seismic signals, and video stills of sediment-enveloped instruments recorded during the 2011–2015 Cascadia Initiative experiment, and seafloor morphology. Our results led us to suggest that thick accretionary prism sediments amplified and extended seismic wave durations from the 11 April 2012 M w 8.6 Indian Ocean earthquake, located more than 13,500 km away. These waves triggered a sequence of small slope failures on the Cascadia margin that led to sediment gravity flows culminating in turbidity currents. Previous studies have related the triggering of sediment-laden gravity flows and turbidite deposition to local earthquakes, but this is the first study in which the originating seismic event is extremely distant (> 10,000 km). The possibility of remotely triggered slope failures that generate sediment-laden gravity flows should be considered in inferences of recurrence intervals of past great Cascadia earthquakes from turbidite sequences. Future similar studies may provide new understanding of submarine slope failures and turbidity currents and the hazards they pose to seafloor infrastructure and tsunami generation in regions both with and without local earthquakes.

Cascadia subduction zone↗

Alternative source models of very low frequency events

We present alternative source models for very low frequency (VLF) events, previously inferred to be radiation from individual slow earthquakes that partly fill the period range between slow slip events lasting thousands of seconds and low-frequency earthquakes (LFE) with durations of tenths of a second. We show that VLF events may emerge from bandpass filtering a sum of clustered, shorter duration, LFE signals, believed to be the components of tectonic tremor. Most published studies show VLF events occurring concurrently with tremor bursts and LFE signals. Our analysis of continuous data from Costa Rica detected VLF events only when tremor was also occurring, which was only 7% of the total time examined. Using analytic and synthetic models, we show that a cluster of LFE signals produces the distinguishing characteristics of VLF events, which may be determined by the cluster envelope. The envelope may be diagnostic of a single, dynamic, slowly slipping event that propagates coherently over kilometers or represents a narrowly band-passed version of nearly simultaneous arrivals of radiation from slip on multiple higher stress drop and/or faster propagating slip patches with dimensions of tens of meters (i.e., LFE sources). Temporally clustered LFE sources may be triggered by single or multiple distinct aseismic slip events or represent the nearly simultaneous chance occurrence of background LFEs. Given the nonuniqueness in possible source durations, we suggest it is premature to draw conclusions about VLF event sources or how they scale.

Journal of Geophysical Research B: Solid Earth↗

Reconsidering earthquake scaling

The relationship (scaling) between scalar moment, M 0 , and duration, T , potentially provides key constraints on the physics governing fault slip. The prevailing interpretation of M 0 - T observations proposes different scaling for fast (earthquakes) and slow (mostly aseismic) slip populations and thus fundamentally different driving mechanisms. We show that a single model of slip events within bounded slip zones may explain nearly all fast and slow slip M 0 - T observations, and both slip populations have a change in scaling, where the slip area growth changes from 2-D when too small to sense the boundaries to 1-D when large enough to be bounded. We present new fast and slow slip M 0 - T observations that sample the change in scaling in each population, which are consistent with our interpretation. We suggest that a continuous but bimodal distribution of slip modes exists and M 0 - T observations alone may not imply a fundamental difference between fast and slow slip.

Geophysical Research Letters↗

Demonstration of the Cascadia G‐FAST geodetic earthquake early warning system for the Nisqually, Washington, earthquake

A prototype earthquake early warning (EEW) system is currently in development in the Pacific Northwest. We have taken a two‐stage approach to EEW: (1) detection and initial characterization using strong‐motion data with the Earthquake Alarm Systems (ElarmS) seismic early warning package and (2) the triggering of geodetic modeling modules using Global Navigation Satellite Systems data that help provide robust estimates of large‐magnitude earthquakes. In this article we demonstrate the performance of the latter, the Geodetic First Approximation of Size and Time (G‐FAST) geodetic early warning system, using simulated displacements for the 2001 M w 6.8 Nisqually earthquake. We test the timing and performance of the two G‐FAST source characterization modules, peak ground displacement scaling, and Centroid Moment Tensor‐driven finite‐fault‐slip modeling under ideal, latent, noisy, and incomplete data conditions. We show good agreement between source parameters computed by G‐FAST with previously published and postprocessed seismic and geodetic results for all test cases and modeling modules, and we discuss the challenges with integration into the U.S. Geological Survey’s ShakeAlert EEW system.

Washington↗