Geology topics

Ryan D. Gold

Publications and source records attributed to Ryan D. Gold.

At least 37 records · Page 2Linked to original sources

USGS2022

Portable optically stimulated luminescence age map of a paleoseismic exposure

The quality and quantity of geochronologic data used to constrain the history of major earthquakes in a region exerts a first-order control on the accuracy of seismic hazard assessments that affect millions of people. However, evaluations of geochronological data are limited by uncertainties related to inherently complex depositional processes that may vary spatially and temporally. To improve confidence in models of earthquake timing, we use a high-density suite of radiocarbon and optically stimulated luminescence (OSL) ages with a grid of 342 portable OSL samples to explore spatiotemporal trends in geochronological data across an exemplary normal fault colluvial wedge exposure. The data reveal a two-dimensional age map of the paleoseismic exposure and demonstrate how vertical and horizontal trends in age relate to dominant sedimentary facies and soil characteristics at the site. Portable OSL data provide critical context for the interpretation of 14 C and OSL ages, show that geochronologic age boundaries between pre- and post-earthquake deposits do not match stratigraphic contacts, and provide the basis for selecting alternate Bayesian models of earthquake timing. Our results demonstrate the potential to use emergent, portable OSL methods to dramatically improve paleoseismic constraints on earthquake timing.

Utah↗

USGS2022

Luminescence sediment tracing reveals the complex dynamics of colluvial wedge formation

Paleoearthquake studies that inform seismic hazard rely on assumptions of sediment transport that remain largely untested. Here, we test a widespread conceptual model and a new numerical model on the formation of colluvial wedges, a key deposit used to constrain the timing of paleoearthquakes. We perform this test by applying luminescence, a sunlight-sensitive sediment tracer, at a field site displaying classic colluvial wedge morphostratigraphy. The model and data comparison reveals complex sediment transport processes beyond the predictions of either conceptual or numerical models, including periods of simultaneous debris and wash facies forming processes, erosion, and reworking. These processes could lead to preservation bias, such as incomplete or overinterpretable paleoearthquake records, given the right environmental conditions. Attention to the site-specific mechanics of fault zone depositional systems, such as via sediment tracing, may buffer against the possible effects of preservation bias on paleoseismic study.

Science↗

USGS2021

Geophysical constraints on the crustal architecture of the transtensional Warm Springs Valley fault zone, northern Walker Lane, western Nevada, USA

The Walker Lane is a zone of distributed transtension where normal faults are overprinted by strike-slip motion. We use two newly-acquired high-resolution seismic reflection profiles and a reprocessed Consortium for Continental Reflection Profiling (COCORP) deep crustal reflection profile to assess the subsurface geometry of the Holocene-active, transtensional Warm Springs Valley fault zone (WSVFZ) near Reno, Nevada, USA. Our multi-scale observations extend to 12 km depth and suggest that the WSVFZ is more complex in the subsurface than implied by late Pleistocene surface fault traces. Two ~4-km-long high-resolution profiles image to a depth of ~2 km and reveal moderately dipping reflections and truncations, some of which project to mapped scarps formed in late Pleistocene surfaces. The shallow lines are co-located with COCORP profile NV 08 along ~40° N latitude. Re-analysis of the COCORP data reveals previously unidentified coherent reflections to a depth of ~12 km and a previously mapped ~30 west-dipping fault at 8-12 km. From these seismic profiles, the WSVFZ is not a simple, sub-vertical fault zone extending through the entire seismogenic crust. Instead, the reflections are consistent with a zone of steep- and moderately-dipping faults that simplify and steepen with depth before intersecting a mid-crustal, low angle (~25-30°) fault. The complex fault geometry of the WSVFZ implies that crustal shear is accommodated by a mix of dipping and subvertical faults in the transtensional northern Walker Lane. If so, transtensional fault zones may present challenges to paleoseismic and geodetic studies and require careful treatment when included in seismic hazard analyses.

Nevada↗

USGS2021

STEPS: Slip time earthquake path simulations applied to the San Andreas and Toe Jam Hill Faults to redefine geologic slip rate uncertainty

Geologic slip rates are a time-averaged measurement of fault displacement calculated over hundreds to million-year time scales and are a primary input for probabilistic seismic hazard analyses, which forecast expected ground shaking in future earthquakes. Despite their utility for seismic hazard calculations, longer-term geologic slip rates represent a time-averaged measure of the tempo of strain release and do not measure variability across earthquake cycles. We have developed a numerical approach called STEPS (Slip Time Earthquake Path Simulations), which is built upon field-based observations and explicitly incorporates realistic variations in displacement per event and variability in the recurrence interval between earthquakes. The STEPS approach, which simulates strain release through time, relies on representing earthquake cycles as stairsteps, rather than straight-line paths, connecting per earthquake time-displacement coordinates. We simulate earthquake histories based on these input constraints using two examples: the Carrizo section of the San Andreas fault and the Toe Jam Hill fault of the Seattle fault zone. We find that modeled slip rate distributions agree with slip rates reported for the sites of interest by the original investigators, while providing a slip rate distribution that reflects the variability of earthquake frequency and displacement. The STEPS approach provides an estimate of fault slip rate uncertainty based on a synthetic suite of plausible time-displacement paths resulting from individual earthquakes, rather than measurement uncertainties associated with offset features. When considering this simulated earthquake behavior between measurements, the uncertainty associated with earthquake paths is greater than that calculated by the long-term rate.

California, Washington↗

USGS2021

Diffuse deformation and surface faulting distribution from sub-metric image correlation along the 2019 Ridgecrest ruptures (California, USA)

The 2019 M w "> M w 6.4 and 7.1 Ridgecrest, California, earthquake sequence (July 2019) ruptured consecutively a system of high‐angle strike‐slip cross faults (northeast‐ and northwest‐trending) within 34 hr. The complex rupture mechanism was illuminated by seismological and geodetic data, bringing forward the issue of the interdependency of the two fault systems both at depth and at the surface, and of its effect on the final surface displacement pattern. Here, we use high‐resolution (WorldView and Pleiades) optical satellite image correlation to measure the near‐fault horizontal and vertical surface displacement fields at 0.5 m ground resolution for the two earthquakes. We point out significant differences with previous geodetic‐ and geologic‐based measurements, and document the essential role of distributed faulting and diffuse deformation in producing the observed surface displacement patterns. We derive strain fields from the horizontal displacement maps, and highlight the predominant role of rotation and shear strain in the surface rupture process. We discuss the segmentation of the rupture based on the fault geometry and along‐strike slip variations. We also image several northeast‐trending faults with similar orientation to the deeply embedded shear fabric identified in aftershock studies, and show that these cross faults are present all along the rupture, including at a scale <100 m. Finally, we compare our results to kinematic slip inversions, and show that the surface diffuse deformation is primarily associated with areas of shallow slip deficit; however, this diffuse deformation cannot be explained using elastic modeling. We conclude that inelastic processes play an important role in contributing to the total surface deformation associated with the 2019 Ridgecrest sequence.

California↗

USGS2021

Quaternary reelfoot fault deformation in the Obion River Valley, Tennessee, USA

Blind reverse faults are challenging to detect, and earthquake records can be elusive because deep fault slip does not break the surface along readily recognized scarps. The blind Reelfoot fault in the New Madrid seismic zone in the central United States has been the subject of extensive prior investigation; however, the extent of slip at the southern portion of the fault remains unconstrained. In this study, we use lidar to map terraces and lacustrine landforms in the Obion River valley and investigate apparent broad folding resulting from slip on the buried Reelfoot fault. We compare remote surface mapping results with three auger boreholes in the ∼24 ka Finley terrace and interpret apparent warping as due to tectonic folding and not stratigraphic thickening. We combine our results with historical records of coseismic lake formation that indicate surface deformation dammed the Obion River in the 1812 CE earthquake. Older terraces (deposited at least 35–55 ka) record progressive fold scarps ≥1, ≥2, and ≥8 m high indicating a long record of earthquakes predating the existing paleoseismic record. Broad, distributed folding above the Reelfoot fault into the Obion River valley is consistent with a deep active fault tip along the southern reaches of the fault. Our analyses indicate the entire length of the fault (≥70 km) is capable of rupture and is more consistent with longer rupture scenarios.

Tennessee↗

USGS2021

Quick and dirty (and accurate) 3-D paleoseismic trench models using coded scale bars

Structure‐from‐motion (SfM) modeling has dramatically increased the speed of generating geometrically accurate orthophoto mosaics of paleoseismic trenches, but some aspects of this technique remain time and labor intensive. Model accuracy relies on control points to establish scale, reduce distortion, and orient 3D models. Traditional SfM methods use total station or Global Navigation Satellite System (GNSS) surveys to constrain models, but collecting control points along a vertical trench wall is often inhibited by poor line of sight to the survey sensor or limited sky view and requires many hours in the field and office. We used physical scale bars printed with coded targets to constrain SfM models of a dusty, 46‐m‐long trench excavation across the Teton fault (Wyoming, U.S.A.). We present a workflow for generating quick and accurate 3D SfM models and orthophoto mosaics and compare the effectiveness of using scale bar, GNSS, and total‐station control in the models. Our results show that the scale bar model deviates from total station survey points by an average of 3.1 cm (maximum of 5.3 cm). In addition, the scale‐bar model only deviates an average of 1.7 cm (maximum 3.5 cm) when compared to the best model alternative, the SfM model controlled by the total station survey. Scale bars eliminate several hours needed to collect and incorporate control points from total station or GNSS surveys and significantly simplify the workflow, at the cost of slightly increased 3D model and orthophoto mosaic error. Our results further suggest that trench models can be constrained with at least four physical scale bars, but using five to six physical scale bars provides redundant control for field deployment and model optimization. The scale bar method for paleoseismic trenches proves to be portable and fast, minimizes the need for specialized survey equipment, and maintains model accuracy needed for mapping trench walls.

Seismological Research Letters↗

USGS2021

Coseismic surface displacement in the 2019 ridgecrest earthquakes: Comparison of field measurements and optical image correlation results

A fundamental topic in earthquake studies is understanding the extent to which fault rupture at the surface is localized on primary fault strands as opposed to distributed tens to hundreds of meters away from primary ruptures through off‐fault deformation (OFD) via a combination of discrete secondary faulting and bulk deformation. The 2019 Ridgecrest, CA M w 6.4 and M w 7.1 earthquakes provide an opportunity to explore this problem via comparison of published field‐based and mostly on‐fault offset measurements, with new lateral displacement measurements made over length scales of hundreds of meters across the primary rupture using WorldView satellite images collected before and after the earthquakes. A mean fit to the field observations underestimates net slip (pixel‐correlation results) by an average of 0.4 m along the M w 6.4 rupture (∼41% of net left‐lateral displacement) and 0.7 m along the M w 7.1 rupture (∼65% of net right‐lateral displacement). We attribute these differences to substantial OFD along the lengths of the M w 6.4 (∼59%) and M w 7.1 (∼35%) ruptures. A maximum fit to the field observations provides an improved match to the pixel correlation results (difference of 0.1–0.2 m or ∼76%–98% of net displacement). OFD may in part depend on displacement magnitude, where locations along the rupture with smaller displacements (<2 m) are characterized by a greater percentage of OFD. Furthermore, OFD increases with rupture zone width, especially at the scale of individual rupture strands (hundreds of meters). These findings contribute to a growing understanding of the important role of OFD in accommodating deformation near the surface.

California↗

USGS2021

Preface to the Focus Section on the 2020 Intermountain West earthquakes

The Intermountain West region of the United States extends from the eastern margin of the Sierra Nevada and Cascade Mountains in the west to the Rocky Mountains in the east. The region is characterized by dextral shear along the eastern margin of the Sierra Nevada and nearly east-west extension in the Basin and Range. This region experienced four significant earthquake sequences in the first half of 2020. The most significant mainshocks were the 18 March 2020 Mw 5.7 earthquake north of Magna, Utah (a suburb of Salt Lake City), the 31 March 2020 Mw 6.5 earthquake northwest of Stanley, Idaho, the 15 May 2020 Mw 6.5 earthquake in the Monte Cristo Range, northwest of Tonopah, Nevada, and the 24 June 2020 Mw 5.8 earthquake near Lone Pine, California. The 15 articles appearing in this focus section explore timely and important topics associated with these sequences, including kinematic rupture models, near-field ground motions, aftershock statistics, geologic observations, seismic hazard implications, and seismotectonics. It is noteworthy that the efforts to record and characterize these earthquake sequences took place during travel and work restrictions necessitated by the COVID-19 pandemic.

Seismological Research Letters↗

USGS2020

Holocene paleoseismology of the Steamboat Mountain Site: Evidence for full‐Llngth rupture of the Teton Fault, Wyoming

The 72‐km‐long Teton fault in northwestern Wyoming is an ideal candidate for reconstructing the lateral extent of surface‐rupturing earthquakes and testing models of normal‐fault segmentation. To explore the history of earthquakes on the northern Teton fault, we hand‐excavated two trenches at the Steamboat Mountain site, where the east‐dipping Teton fault has vertically displaced west‐sloping alluvial‐fan surfaces. The trenches exposed glaciofluvial, alluvial‐fan, and scarp‐derived colluvial sediments and stratigraphic and structural evidence of two surface‐rupturing earthquakes (SM1 and SM2). A Bayesian geochronologic model for the site includes three optically stimulated luminescence ages ( ⁠ ∼ 12 – 17 ka ⁠ ) for the glaciofluvial units and 16 radiocarbon ages ( ⁠ ∼ 1.2 – 8.6 ka ⁠ ) for the alluvial‐fan and colluvial units and constrains SM1 and SM2 to 5.5 ± 0.2 ka , 1 σ (5.2–5.9 ka, 95%) and 9.7 ± 0.9 ka , 1 σ (8.5–11.5 ka, 95%), respectively. Structural, stratigraphic, and geomorphic relations yield vertical displacements for SM1 ( ⁠ 2.0 ± 0.6 m , 1 σ ⁠ ) and SM2 ( ⁠ 2.0 ± 1.0 m , 1 σ ⁠ ). The Steamboat Mountain paleoseismic chronology overlaps temporally with earthquakes interpreted from previous terrestrial and lacustrine paleoseismic data along the fault. Integrating these data, we infer that the youngest Teton fault rupture occurred at ∼ 5.3 ka ⁠ , generated 1.7 ± 1.0 m , 1 σ of vertical displacement along 51–70 km of the fault, and had a moment magnitude ( ⁠ M w ⁠ ) of ∼ 7.0 – 7.2 ⁠ . This rupture was apparently unimpeded by structural complexities along the Teton fault. The integrated chronology permits a previous full‐length rupture at ∼ 10 ka and possible partial ruptures of the fault at ∼ 8 – 9 ka ⁠ . To reconcile conflicting terrestrial and lacustrine paleoseismic data, we propose a hypothesis of alternating full‐ and partial‐length ruptures of the Teton fault, including M w ∼ 6.5 – 7.2 earthquakes every ∼ 1.2 ky ⁠ . Additional paleoseismic data for the northern and central sections of the fault would serve to test this bimodal rupture hypothesis.

Wyoming↗

USGS2020

Seismic reflection imaging of the low-angle Panamint normal fault system, eastern California

Shallowly dipping (<30°) low‐angle normal faults (LANFs) have been documented globally; however, examples of active LANFs in continental settings are limited. The western margin of the Panamint Range in eastern California is defined by a LANF that dips west beneath Panamint Valley and has evidence of Quaternary motion. In addition, high‐angle dextral‐oblique normal faults displace middle to late Quaternary alluvial fans near the range front. To image shallow (<1 km depth), crosscutting relationships between the low‐ and high‐angle faults along the range front, we acquired two high‐resolution P wave seismic reflection profiles. The northern, 4.6‐km‐long profile crosses the 2‐km‐wide Wildrose graben and the southern, 0.8‐km‐long profile extends onto the Panamint Valley playa, ~7.5 km S of Ballarat, CA. The profile across the Wildrose graben reveals a robust, low‐angle reflector interpreted to represent the LANF separating Plio‐Pleistocene alluvial fanglomerate and Proterozoic metasedimentary deposits. High‐angle faults interpreted in the seismic profile correspond to fault scarps on Quaternary alluvial fan surfaces. Interpretation of the reflection data suggests that the high‐angle faults vertically displace the LANF up to 80 m within the Wildrose graben. Similarly, the profile south of Ballarat reveals a low‐angle reflector, which appears both rotated and displaced up to 260 m by high‐angle faults. These results suggest that near the Panamint range front, the high‐angle faults are the dominant active structures. We conclude that at least at shallow (<1 km) depths, the LANF we imaged is not active today.

California↗

USGS2020

Localized fault-zone dilatancy and surface inelasticity of the 2019 Ridgecrest earthquakes

Earthquakes produce a spectrum of elastic and inelastic deformation processes that are reflected across various length and time scales. While elasticity has long dominated research assumptions in active tectonics, increasing interest has focused on the inelastic characteristics of earthquakes, particularly those of the surface fault rupture zone itself, and how they relate to ground rupture hazard and the mechanics of damage zones. Here we present detailed co-seismic surface-strain analysis of the 2019 Ridgecrest, California, earthquakes. We derive three-dimensional high-resolution surface displacements from satellite optical imagery, which we then invert for the co-seismic surface-strain tensors. We show that fault-zone dilation is pervasive throughout these earthquakes and that inelastic failure is present but relatively localized (median width of 31 m). The width of the inelastic failure zone is not correlated to off-fault deformation, surface geology or displacement magnitude. Instead, the extent and kinematics of inelastic failure reflect active, mylonitic deformation of the fault damage zone that is influenced by rupture velocity and fault maturity. These results highlight how a single earthquake contributes to the long-term, permanent geologic record of faulting.

California↗

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

Evidence of previous faulting along the 2019 Ridgecrest, California earthquake ruptures

The July 2019 Ridgecrest earthquake sequence in southeastern California was characterized as surprising because only ~35% of the rupture occurred on previously mapped faults. Employing more detailed inspection of pre-event high-resolution topography and imagery in combination with field observations, we document evidence of active faulting in the landscape along the entire fault system. Scarps, deflected drainages, and lineaments and contrasts in topography, vegetation, and ground color demonstrate previous slip on a dense network of orthogonal faults, consistent with patterns of surface rupture observed in 2019. Not all of these newly mapped fault strands ruptured in 2019. Outcrop-scale field observations additionally reveal tufa lineaments and sheared Quaternary deposits. Neotectonic features are commonly short (<2 km), discontinuous, and display en echelon patterns along both the M 6.4 and M 7.1 ruptures. These features are generally more prominent and better preserved outside the late Pleistocene lake basins. Fault expression may also be related to deformation style: scarps and topographic lineaments are more prevalent in areas where substantial vertical motion occurred in 2019. Where strike-slip displacement dominated in 2019, the faults are mainly expressed by less prominent tonal and vegetation features. Both the NE- and NW-trending active fault systems are subparallel to regional bedrock fabrics that were established as early as ~150 Ma, and may be reactivating these older structures. Overall, we estimate that 50-70% (i.e., an additional 15-35%) of the 2019 surface ruptures could have been recognized as active faults with detailed inspection of pre-event data. Similar detailed mapping of potential neotectonic features could help improve seismic hazard analyses in other regions of eastern California and elsewhere that have distributed faulting or incompletely mapped faults. In areas where faults cannot be resolved as single thoroughgoing structures, a zone of potential faulting should be used as a hazard model input.

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↗

USGS2020

Evidence for late Quaternary deformation along Crowley's Ridge, New Madrid seismic zone

The New Madrid seismic zone has been the source of multiple major (M ~7.0–7.5) earthquakes in the past 2 ka, yet the surface expression of recent deformation remains ambiguous. Crowleys Ridge, a linear ridge trending north‐south for 300+ km through the Mississippi Embayment, has been interpreted as either a fault‐bounded uplift or a nontectonic erosional remnant. New and previously published seismic reflection and shallow resistivity data show discontinuities at the ridge margins in Plio‐Pleistocene strata, yet the timing of most recent faulting and the lateral extent of these faults remain unknown. To assess Pleistocene‐to‐recent tectonic activity of Crowleys Ridge, we perform landscape‐scale geomorphic analyses, such as relief, slope, hypsometry, and drainage basin shape, on a 10‐m digital elevation model (DEM). North‐to‐south variations in geomorphic indices indicate Pleistocene‐to‐recent tectonic uplift of the southern ridge. Moreover, mapping on a <1‐m lidar‐derived DEM reveals scarps on late Pleistocene geomorphic surfaces. The scarps are primarily located along the southern ridge, trend parallel to the ridge margin discontinuously for 0.1–1 km, and vertically offset <56 ka surfaces 0.4 m with up to 6 m of tilting. These landscape‐scale patterns and scarps, integrated with discontinuities in the seismic reflection and resistivity data, provide evidence of low‐rate (<0.2 mm/year) late Quaternary tectonic activity along the southern segment of Crowleys Ridge. The interpretations agree with recent tectonic models suggesting southern Crowleys Ridge is a compressional step over in a right‐lateral fault system within the Reelfoot Rift.

Arkansas, Illinois, Kentucky, Missouri, Tennessee↗

USGS2019

Offset channels may not accurately record strike-slip fault displacement: Evidence from landscape evolution models

Slip distribution, slip rate, and slip per event for strike‐slip faults are commonly determined by correlating offset stream channels—under the assumption that they record seismic slip—but offset channels are formed by the interplay of tectonic and geomorphic processes. To constrain offset channel development under known tectonic and geomorphic conditions, we use numerical landscape evolution simulations along a theoretical strike‐slip fault with uniform and steady uplift, erosion, and diffusion. We investigate the influence of four tectonic parameters (fault zone width, earthquake recurrence interval, variance of the recurrence interval, and total slip relative to channel spacing) on offset channel development through multiple earthquake cycles. Analysis of >3,000 automatically measured offsets from >135 simulations suggests ~30% variability in individual measurements, but modeled displacement is recovered by averaging multiple measurements. However, the average of multiple offset measurements systematically underestimates modeled slip except when the fault zone is less than ~5 m wide, total slip is less than channel spacing, and offsets are measured shortly after an earthquake. In these simulations, postearthquake landscape evolution widens the geomorphic expression of the fault zone and modifies apparent channel offsets. We distinguish this “geomorphic fault zone” from the tectonic fault zone (zone of coseismic distributed deformation). This study highlights the capability of landscape evolution models to explore a range of conditions not easily defined in natural examples and the importance of averaging multiple measurements. Our results verify that paleoseismic studies must consider how geomorphic change has modified offset markers and use caution interpreting slip histories with multiple earthquakes.

Journal of Geophysical Research Solid Earth↗

USGS2019

Holocene earthquake history and slip rate of the southern Teton fault, Wyoming, USA

The 72-km-long Teton normal fault bounds the eastern base of the Teton Range in northwestern Wyoming, USA. Although geomorphic surfaces along the fault record latest Pleistocene to Holocene fault movement, the postglacial earthquake history of the fault has remained enigmatic. We excavated a paleoseismic trench at the Buffalo Bowl site along the southernmost part of the fault to determine its Holocene rupture history and slip rate. At the site, ∼6.3 m of displacement postdates an early Holocene (ca. 10.5 ka) alluvial-fan surface. We document evidence of three surface-faulting earthquakes based on packages of scarp-derived colluvium that postdate the alluvial-fan units. Bayesian modeling of radiocarbon and luminescence ages yields earthquake times of ca. 9.9 ka, ca. 7.1 ka, and ca. 4.6 ka, forming the longest, most complete paleoseismic record of the Teton fault. We integrate these data with a displaced deglacial surface 4 km NE at Granite Canyon to calculate a postglacial to mid-Holocene (14.4−4.6 ka) slip rate of ∼1.1 mm/yr. Our analysis also suggests that the postglacial to early Holocene (14.4−9.9 ka) slip rate exceeds the Holocene (9.9−4.6 ka) rate by a factor of ∼2 (maximum of 3); however, a uniform rate for the fault is possible considering the 95% slip-rate errors. The ∼5 k.y. elapsed time since the last rupture of the southernmost Teton fault implies a current slip deficit of ∼4−5 m, which is possibly explained by spatially/temporally incomplete paleoseismic data, irregular earthquake recurrence, and/or variable per-event displacement. Our study emphasizes the importance of minimizing slip-rate uncertainties by integrating paleoseismic and geomorphic data sets and capturing multiple earthquake cycles.

Wyoming↗