Geology topics

Stephen J. Angster

Publications and source records attributed to Stephen J. Angster.

15 recordsLinked to original sources

USGS2026

Latest Pleistocene to 19th-century earthquakes on bending-moment reverse faults of the Seattle fault zone, Washington

Fault-related folds and their associated secondary faults play a critical yet often underrecognized role in accommodating strain and generating earthquakes in active fold-and-thrust belts. In the Seattle fault zone (SFZ), Washington, USA, we present new paleoseismic, geomorphic, and geophysical evidence for late Pleistocene and Holocene earthquakes on shallow, south-dipping secondary faults—the Lytle Beach and Vasa Park faults—that lie within the hanging wall of the greater SFZ and are on trend with the primary, blind Blakely Harbor fault. Our data show that these structures have ruptured independently, producing localized uplift and deformation at the surface, with the most recent event (RH2) likely occurring in the early nineteenth century. While a temporal overlap between the late Pleistocene RH1 and VP1 earthquakes raises the possibility of a ≥35 km rupture along the Blakely Harbor fault, structural and temporal evidence instead supports independent rupture on individual faults related to folding. We interpret these faults as bending-moment reverse faults that formed within a synclinal hinge zone of the main fault, reflecting mechanical and kinematic influences of the broader fault system. Combined with prior studies, our findings indicate that faulting related to folding dominates the mode of strain release within the SFZ since the late Pleistocene with more frequent earthquake recurrence (∼350 yr) over the past ∼2500 yr.

Washington↗

USGS2024

Debris avalanches in the northern California Coast Range triggered by plate boundary earthquakes

Determining the timing and cause for ancient hillslope failures proves difficult in the western United States, yet critical as it ties directly into groundmotion estimates for hazardous events. This knowledge gap is important to confront as hillslope failures are candidates to be triggered by earthquakes along active plate boundaries. We identify two prehistoric, i.e., preinstrumental history, debris avalanches (3 10^6 and 6 10^6 m^3) in the Coast Ranges of northern California. These debris avalanches are well suited for studying the timing (to the exact year) and trigger as the densely forested landscape enables effective use of dendrochronology and high-resolution radiocarbon to pinpoint time of failure, and the legacy of landslide studies in the region provides context for evaluating a seismic trigger as the most likely triggering mechanism. Neither debris avalanche was triggered by the largest precipitation event of the instrumental record in 1964 CE, and the two debris avalanche sites are physiographically suited to accommodate topographic amplification of seismic shaking. Through a suite of geologic and dendrochronologic evidence, we establish the likely time of failure of the two debris avalanches as 1906 CE and 933 CE. The failure of the younger landslide coincides with the 1906 CE San Francisco earthquake on the Pacific/North American plate boundary, while the older landslide dates to 933 CE falling within the broad age range (850-966 CE) of a likely megathrust earthquake on the Cascadia subduction zone as recorded in coastal marshes in Humboldt Bay, California. The precise age on the 933 CE debris avalanche could trim the uncertainty on a Cascadia subduction zone event from 116 years to one year. Utilizing debris avalanche records from sites suitable to record seismic shaking improves understanding of plate-boundary earthquake timing and extent of shaking.

Bulletin of the Seismological Society of America↗

USGS2024

Field observations and logs from the Rose Hip trench exposure across a north-facing scarp within the Seattle Fault Zone, southern Bainbridge Island, Washington

The Seattle Fault Zone is an approximately 70-km-long, east-west-trending zone of south-dipping blind reverse faults within the Puget lowland region in Washington. Because of the proximity, the Seattle Fault Zone poses a significant earthquake hazard to the Puget sound and Seattle metropolitan regions. We present preliminary mapping and trench-site information from a paleoseismic investigation across a newly identified active fault scarp located within the hanging wall of the Seattle Fault Zone on southern Bainbridge Island, Washington. The trench exposed monoclinally folded Miocene bedrock, fractured and faulted glacial-related deposits, and laminated lacustrine deposits capped by slope-derived colluvium. The observations from this investigation record late Pleistocene to Holocene north-vergent folding and faulting along this new fault scarp.

Washington↗

USGS2023

Logs and data from the Starthistle trench across a scarp within the Wallula Fault Zone, southeastern Washington

Introduction The Wallula Fault Zone is composed of a series of northwest-trending faults and folds that coincide with a prominent magnetic anomaly that extends uninterrupted for approximately 120 kilometers within the Cascadia back arc of southeastern Washington and northeastern Oregon. It is part of the geologic structures associated with the topographic lineament known as the Olympic-Wallowa lineament and represents a relatively narrow zone of active faulting and seismicity that trend along the northern flank of the Horse Heaven Hills. Prior paleoseismic study of the Wallula Fault at an exposure at Finley Quarry indicates multiple Quaternary ruptures, including a Holocene liquefaction event, demonstrating that this fault zone poses a seismic hazard to the Tri-Cities region (Richland, Pasco, Kennewick) in southeastern Washington. Recent airborne light detection and ranging (lidar) data coverage east of the Columbia River reveals an approximately 0.5-meter-high south-facing scarp east of the Columbia River that extends, almost continuously, for about 25 km to the east along the base of the Horse Heaven Hills. As part of an effort to assess and characterize the seismic hazard posed by the Wallula Fault Zone, we excavated and studied a trench exposure across the scarp to understand its origin and the potential history of rupture along the fault zone. We present preliminary mapping and trench site information from a paleoseismic investigation. These field and laboratory data may support development of a history of the latest Pleistocene and Holocene surface rupture within the Wallula Fault Zone.

Washington↗

USGS2022

Simplifying complex fault data for systems-level analysis: Earthquake geology inputs for U.S. NSHM 2023

As part of the U.S. National Seismic Hazard Model (NSHM) update planned for 2023, two databases were prepared to more completely represent Quaternary-active faulting across the western United States: the NSHM23 fault sections database (FSD) and earthquake geology database (EQGeoDB). In prior iterations of NSHM, fault sections were included only if a field-measurement-derived slip rate was estimated along a given fault. By expanding this inclusion criteria, we were able to assess a larger set of faults for use in NSHM23. The USGS Quaternary Fault and Fold Database served as a guide for assessing possible additions to the NSHM23 FSD. Reevaluating available data from published sources yielded an increase of fault sections from ~650 faults in NSHM18 to ~1,000 faults proposed for use in NSHM23. EQGeoDB, a companion dataset linked to NSHM23 FSD, contains geologic slip rate estimates for fault sections included in FSD. Together, these databases serve as common input data used in deformation modeling, earthquake rupture forecasting, and additional downstream uses in NSHM development.

Scientific Data↗

USGS2021

Improved scaling relationships for seismic moment and average slip of strike-slip earthquakes incorporating fault slip rate, fault width and stress drop

We develop a self‐consistent scaling model relating magnitude M w to surface rupture length ( ⁠L E ⁠ ), surface displacement D E ⁠ , and rupture width W E ⁠ , for strike‐slip faults. Knowledge of the long‐term fault‐slip rate S F improves magnitude estimates. Data are collected for 55 ground‐rupturing strike‐slip earthquakes that have geological estimates of L E ⁠ , D E ⁠ , and S F⁠ , and geophysical estimates of W E ⁠ . We begin with the model of Anderson et al. (2017) , which uses a closed form equation for the seismic moment of a surface‐rupturing strike‐slip fault of arbitrary aspect ratio and given stress drop, Δτ C ⁠ . Using W E estimates does not improve M w estimates. However, measurements of D E plus the relationship between Δτ C and surface slip provide an alternate approach to study W E ⁠ . A grid of plausible stress drop and width pairs were used to predict displacement and earthquake magnitude. A likelihood function was computed from within the uncertainty ranges of the corresponding observed M w and D E values. After maximizing likelihoods over earthquakes in length bins, we found the most likely values of W E for constant stress drop; these depend on the rupture length. The best‐fitting model has the surprising form W E ∝logL E —a gentle increase in width with rupture length. Residuals from this model are convincingly correlated to the fault‐slip rate and also show a weak correlation with the crustal thickness. The resulting model thus supports a constant stress drop for ruptures of all lengths, consistent with teleseismic observation. The approach can be extended to test other observable factors that might improve the predictability of magnitude from a mapped fault for seismic hazard analyses.

Bulletin of the Seismological Society of America↗

USGS2020

Characterizing strain between rigid crustal blocks in the southern Cascadia forearc: Quaternary faults and folds of the northern Sacramento Valley, California

Topographic profiles across late Quaternary surfaces in the northern Sacramento Valley (California, USA) show offset and progressive folding on series of active east- and northeast—trending faults and folds. Optically stimulated luminescence ages on deposits draping a warped late Pleistocene river terrace yielded differential incision rates along the Sacramento River and indicate tectonic uplift equal to 0.2 ± 0.1 and 0.6 ± 0.2 mm/yr above the anticline of the Inks Creek fold system and Red Bluff fault, respectively. Uplift rates correspond to a total of 1.3 ± 0.4 mm/yr of north-directed crustal shortening, accounting for all of the geodetically observed contractional strain in the northern Sacramento Valley, but only part of the far-field contraction between the Sierra Nevada–Great Valley and Oregon Coast blocks. These structures define the southern limit of the transpressional transition between the two blocks.

California↗

USGS2020

Post-glacial Mw 7.0-7.5 earthquakes on the North Olympic fault zone, Washington

Holocene crustal faulting in the northern Olympic Peninsula of Washington State manifests in a zone of west‐northwest‐striking crustal faults herein named the North Olympic fault zone, which extends for ∼ 80    km "> ∼ 80 km ∼80 km along strike and includes the Lake Creek–Boundary Creek fault to the east and the Sadie Creek fault and newly discovered scarps to the west. This study focuses on the Sadie Creek fault, which extends for > 14    km "> > 14 km >14 km west‐northwest from Lake Crescent. Airborne light detection and ranging (lidar) imagery reveals the trace of the Sadie Creek fault and offset postglacial landforms showing a history of Holocene surface‐rupturing earthquakes dominated by dextral displacement along a steeply dipping fault zone. Paleoseismic trenches at two sites on the Sadie Creek fault reveal till and outwash overlain by progressively buried forest and wetland soils developed on scarp‐derived colluvial wedges. Trench exposures of complex faulting with subhorizontal slickenlines indicate dextral displacement with lesser dip slip. Correlation of broadly constrained time intervals for earthquakes at the Sadie Creek sites and those to the east along the Lake Creek–Boundary Creek fault is consistent with rupture of much of the length of the North Olympic fault zone three to four times: at about 11, 7, 3, and 1 ka, with a shorter rupture at about 8.5 ka. Dated ruptures from trenches only partially coincide with coseismic landslides and megaturbidites in Lake Crescent, indicating that some earthquakes did not trigger megaturbidites, and some turbidites were unrelated to local fault rupture. Landform mapping suggests single‐event dextral displacement of 4 ± 1    m "> 4 ± 1 m 4±1 m on the Sadie Creek fault. Inferred maximum rupture length and single‐event slip imply earthquake magnitudes M w "> M w Mw 7.0–7.5. Dextral slip rates of 1.3 – 2.3    mm / yr "> 1.3 – 2.3 mm / yr 1.3–2.3 mm/yr and the ∼ 11 , 000    yr "> ∼ 11 , 000 yr ∼11,000 yr slip history suggest that the North Olympic fault zone is a prominent contributor to permanent strain in the northern Cascadia fore‐arc.

Washington↗

USGS2020

Documentation of Surface Fault Rupture and Ground‐Deformation Features Produced by the 4 and 5 July 2019 Mw 6.4 and Mw 7.1 Ridgecrest Earthquake Sequence

The M w "> M w Mw 6.4 and M w "> M w Mw 7.1 Ridgecrest earthquake sequence occurred on 4 and 5 July 2019 within the eastern California shear zone of southern California. Both events produced extensive surface faulting and ground deformation within Indian Wells Valley and Searles Valley. In the weeks following the earthquakes, more than six dozen scientists from government, academia, and the private sector carefully documented the surface faulting and ground‐deformation features. As of December 2019, we have compiled a total of more than 6000 ground observations; approximately 1500 of these simply note the presence or absence of fault rupture or ground failure, but the remainder include detailed descriptions and other documentation, including tens of thousands of photographs. More than 1100 of these observations also include quantitative field measurements of displacement sense and magnitude. These field observations were supplemented by mapping of fault rupture and ground‐deformation features directly in the field as well as by interpreting the location and extent of surface faulting and ground deformation from optical imagery and geodetic image products. We identified greater than 68 km of fault rupture produced by both earthquakes as well as numerous sites of ground deformation resulting from liquefaction or slope failure. These observations comprise a dataset that is fundamental to understanding the processes that controlled this earthquake sequence and for improving earthquake hazard estimates in the region. This article documents the types of data collected during postearthquake field investigations, the compilation effort, and the digital data products resulting from these efforts.

California↗

USGS2020

Surface displacement distributions for the July 2019 Ridgecrest, California earthquake ruptures

Surface rupture in the 2019 Ridgecrest, California, earthquake sequence occurred along two orthogonal cross faults and includes dominantly left‐lateral and northeast‐striking rupture in the M w "> M w 6.4 foreshock and dominantly right‐lateral and northwest‐striking rupture in the M w "> M w 7.1 mainshock. We present > 650 "> > 650 field‐based, surface‐displacement observations for these ruptures and synthesize our results into cumulative along‐strike displacement distributions. Using these data, we calculate displacement gradients and compare our results with historical strike‐slip ruptures in the eastern California shear zone. For the M w "> M w 6.4 rupture, we report 96 displacements measured along 18 km of northeast‐striking rupture. Cumulative displacement curves for the rupture yield a mean left‐lateral displacement of 0.3–0.5 m and maximum of 0.7–1.6 m. Net mean vertical displacement based on the difference of down‐to‐the‐west (DTW) and down‐to‐the‐east (DTE) displacement curves is close to zero (0.02 m DTW). The M w "> M w 6.4 displacement distribution shows that the majority of displacement occurred southwest of the intersection with the M w "> M w 7.1 rupture. The M w "> M w 7.1 rupture is northwest‐striking and 50 km long based on 576 field measurements. Displacement curves indicate a mean right‐lateral displacement of 1.2–1.7 m and a maximum of 4.3–7.0 m. Net vertical displacement in the rupture averages 0.3 m DTW. The M w "> M w 7.1 displacement distributions demonstrate that maximum displacement occurred along a 12‐km‐long portion of the fault near the M w "> M w 7.1 epicenter, releasing 66% of the geologically based seismic moment along 24% of the total rupture length. Using our displacement distributions, we calculate kilometer‐scale displacement gradients for the M w "> M w 7.1 rupture. The steepest gradients ( ⁠ ∼ 1 – 3    m / km "> ∼ 1 – 3 m / km ⁠ ) flank the 12‐km‐long region of maximum displacement. In contrast, gradients for the 1992 M w "> M w 7.3 Landers and 1999 M w "> M w 7.1 Hector Mine earthquakes are < 0.6    m / km "> < 0.6 m / km . Our displacement distributions are important for understanding the influence of cross‐fault rupture on M w "> M w 6.4 and 7.1 rupture length and displacement and will facilitate comparisons with distributions generated remotely and at broader scales.

California↗

USGS2017

P- and S-wave velocity models incorporating the Cascadia subduction zone for 3D earthquake ground motion simulations, Version 1.6—Update for Open-File Report 2007–1348

In support of earthquake hazard studies and ground motion simulations in the Pacific Northwest, three-dimensional P- and S-wave velocity ( V P and V S , respectively) models incorporating the Cascadia subduction zone were previously developed for the region encompassed from about 40.2°N. to 50°N. latitude, and from about 122°W. to 129°W. longitude. This report describes updates to the Cascadia velocity property volumes of model version 1.3 (V1.3), herein called model version 1.6 (V1.6). As in model V1.3, the updated V1.6 model volume includes depths from 0 kilometers (mean sea level) to 60 kilometers, and it is intended to be a reference for researchers who have used, or are planning to use, this model in their Earth science investigations. To this end, it is intended that the V P and V S property volumes of model V1.6 will be considered a template for a community velocity model of the Cascadia region as additional results become available. With the recent and ongoing development of the National Crustal Model, we envision any future versions of this model will be directly integrated with that effort.

Oregon, Vancouver, Washington↗

USGS2017

3D ground‐motion simulations of Mw 7 earthquakes on the Salt Lake City segment of the Wasatch fault zone: Variability of long‐period (T≥1 s) ground motions and sensitivity to kinematic rupture parameters

We examine the variability of long‐period ( T ≥1 s) earthquake ground motions from 3D simulations of M w 7 earthquakes on the Salt Lake City segment of the Wasatch fault zone, Utah, from a set of 96 rupture models with varying slip distributions, rupture speeds, slip velocities, and hypocenter locations. Earthquake ruptures were prescribed on a 3D fault representation that satisfies geologic constraints and maintained distinct strands for the Warm Springs and for the East Bench and Cottonwood faults. Response spectral accelerations (SA; 1.5–10 s; 5% damping) were measured, and average distance scaling was well fit by a simple functional form that depends on the near‐source intensity level SA 0 ( T ) and a corner distance R c :SA( R , T )=SA 0 ( T )(1+( R / R c )) −1 . Period‐dependent hanging‐wall effects manifested and increased the ground motions by factors of about 2–3, though the effects appeared partially attributable to differences in shallow site response for sites on the hanging wall and footwall of the fault. Comparisons with modern ground‐motion prediction equations (GMPEs) found that the simulated ground motions were generally consistent, except within deep sedimentary basins, where simulated ground motions were greatly underpredicted. Ground‐motion variability exhibited strong lateral variations and, at some sites, exceeded the ground‐motion variability indicated by GMPEs. The effects on the ground motions of changing the values of the five kinematic rupture parameters can largely be explained by three predominant factors: distance to high‐slip subevents, dynamic stress drop, and changes in the contributions from directivity. These results emphasize the need for further characterization of the underlying distributions and covariances of the kinematic rupture parameters used in 3D ground‐motion simulations employed in probabilistic seismic‐hazard analyses.

Utah↗

USGS2017

A paleoseismic transect across the northwestern Basin and Range Province, northwestern Nevada and northeastern California, USA

We use new and existing data to compile a record of ∼18 latest Quaternary large-magnitude surface-rupturing earthquakes on 7 fault zones in the northwestern Basin and Range Province of northwestern Nevada and northeastern California. The most recent earthquake on all faults postdates the ca. 18–15 ka last glacial highstand of pluvial Lake Lahontan and other pluvial lakes in the region. These lacustrine data provide a window in which we calculate latest Quaternary vertical slip rates and compare them with rates of modern deformation in a global positioning system (GPS) transect spanning the region. Average vertical slip rates on these fault zones range from 0.1 to 0.8 mm/yr and total ∼2 mm/yr across a 265-km-wide transect from near Paradise Valley, Nevada, to the Warner Mountains in California. We converted vertical slip rates to horizontal extension rates using fault dips of 30°–60°, and then compared the extension rates to GPS-derived rates of modern (last 7–9 yr) deformation. Our preferred fault dip values (45°–55°) yield estimated long-term extension rates (1.3–1.9 mm/yr) that underestimate our modern rate (2.4 mm/yr) by ∼21%–46%. The most likely sources of this underestimate are geologically unrecognizable deformation from moderate-sized earthquakes and unaccounted-for coseismic off-fault deformation from large surface-rupturing earthquakes. However, fault dip values of ≤40° yield long-term rates comparable to or greater than modern rates, so an alternative explanation is that fault dips are closer to 40° than our preferred values. We speculate that the large component of right-lateral shear apparent in the GPS signal is partitioned on faults with primary strike-slip displacement, such as the Long Valley fault zone, and as not easily detected oblique slip on favorably oriented normal faults in the region.

California, Nevada↗

USGS2014

Uplift and subsidence reveal a nonpersistent megathrust rupture boundary (Sitkinak Island, Alaska)

We report stratigraphic evidence of land-level change and tsunami inundation along the Alaska-Aleutian megathrust during prehistoric and historical earthquakes west of Kodiak Island. On Sitkinak Island, cores and tidal outcrops fringing a lagoon reveal five sharp lithologic contacts that record coseismic land-level change. Radiocarbon dates, 137 Cs profiles, CT scans, and microfossil assemblages are consistent with rapid uplift ca. 290-0, 520-300, and 1050-790 cal yr BP, and subsidence in AD 1964 and ca. 640-510 cal yr BP. Radiocarbon, 137 Cs, and 210 Pb ages bracketing a sand bed traced 1.5 km inland and evidence for sudden uplift are consistent with Russian accounts of an earthquake and tsunami in AD 1788. The mixed uplift and subsidence record suggests that Sitkinak Island sits above a non-persistent boundary near the southwestern limit of the AD 1964 Mw 9.2 megathrust rupture.

Alaska↗

USGS2012

Modeling of wave driven circulation and water quality in nearshore environments

In order to investigate the effects of nearshore discharges of water quality degrading substances and bacteria in coastal environments, models capable of predicting nearshore circulation due to local wave and tide conditions are required. One of the larger challenges to nearshore coastal modeling is accurately reproducing nearshore circulation due to wave action. Local wave action not only drives circulation through processes such as longshore transport and rip currents, but also contributes significantly to the mixing of water quality constituents. In the present work, a wave model was used to calculate radiation shear stresses and dissipation due to wave action. The shear stresses and dissipation were incorporated into a hydrodynamic model to force circulation in the nearshore environment. The model was applied to a site in Santa Cruz, CA where site specific current data was available. The model reproduces the nearshore current structure observed in the region and was used to study the transport of dredge disposal plumes in the region which could have deleterious effects on local beaches. This presentation will outline the nearshore circulation model development and application.

Conference Paper↗