Geology topics

Ryan D. Gold

Publications and source records attributed to Ryan D. Gold.

At least 55 records · Page 3Linked to original sources

USGS2019

Relaxing segmentation on the Wasatch Fault Zone: Impact on seismic hazard

The multisegment Wasatch fault zone is a well-studied normal fault in the western United States that has paleoseismic evidence of recurrent Holocene surface-faulting earthquakes. Along the 270-km-long central part of the fault, four primary structural complexities provide possible along-strike limits to these ruptures and form the basis for models of fault segmentation. Here, we assess the impact that the Wasatch fault segmentation model has on seismic hazard by evaluating the time-independent long-term rate of ruptures on the fault that satisfy fault slip rates and paleoseismic event rates, adapting standard inverse theory used in the Uniform California Earthquake Rupture Forecast 3 (UCERF3), and implementing a segmentation constraint where ruptures across primary structural complexities are penalized. We define three models with varying degrees of rupture penalization: (1) segmented (ruptures confined to individual segments), (2) penalized (multi-segment ruptures allowed, but penalized), and (3) unsegmented (all ruptures allowed). Seismic-hazard results show that on average, hazard is highest for the segmented model, where seismic moment is accommodated by frequent moderate (moment magnitude, M w 6.2–6.8) earthquakes. The unsegmented model yields the lowest average seismic hazard because part of the seismic moment is accommodated by large (M w 6.9–7.9), but infrequent ruptures. We compare these results to model differences derived from other inputs such as slip rate and magnitude scaling relationships and conclude that segmentation exerts a primary control on seismic hazard. This study demonstrates the need for additional geologic constraints on rupture extent and methods by which these observations can be included in hazard-modeling efforts.

Utah↗

USGS2019

Variable normal-fault rupture behavior, northern Lost River fault zone, Idaho, USA

The 1983 Mw 6.9 Borah Peak earthquake generated ∼36 km of surface rupture along the Thousand Springs and Warm Springs sections of the Lost River fault zone (LRFZ, Idaho, USA). Although the rupture is a well-studied example of multisegment surface faulting, ambiguity remains regarding the degree to which a bedrock ridge and branch fault at the Willow Creek Hills influenced rupture progress. To explore the 1983 rupture in the context of the structural complexity, we reconstruct the spatial distribution of surface displacements for the northern 16 km of the 1983 rupture and prehistoric ruptures in the same reach of the LRFZ using 252 vertical-separation measurements made from high-resolution (5–10-cm-pixel) digital surface models. Our results suggest the 1983 Warm Springs rupture had an average vertical displacement of ∼0.3–0.4 m and released ∼6% of the seismic moment estimated for the Borah Peak earthquake and <12% of the moment accumulated on the Warm Springs section since its last prehistoric earthquake. The 1983 Warm Springs rupture is best described as the moderate-displacement continuation of primary rupture from the Thousand Springs section into and through a zone of structural complexity. Historical and prehistoric displacements show that the Willow Creek Hills have impeded some, but not all ruptures. We speculate that rupture termination or penetration is controlled by the history of LRFZ moment release, displacement, and rupture direction. Our results inform the interpretation of paleoseismic data from near zones of normal-fault structural complexity and demonstrate that these zones may modulate rather than impede rupture displacement.

Idaho↗

USGS2019

Surface rupture and distributed deformation revealed by optical satellite imagery: The intraplate 2016 Mw 6.0 Petermann Ranges earthquake, Australia

High-resolution optical satellite imagery is used to quantify vertical surface deformation associated with the intraplate 20 May 2016 Mw 6.0 Petermann Ranges earthquake, Northern Territory, Australia. The 21 ╓ 1 km long NW-trending rupture resulted from reverse motion on a northeast-dipping fault. Vertical surface offsets of up to 0.7 ╓ 0.1 m distributed across a 0.5-to-1 km wide deformation zone are measured using the Iterative Closest Point (ICP) algorithm to compare pre- and post-earthquake digital elevation models (DEMs) derived from Worldview imagery. The results are validated by comparison with field-based observations and interferometric synthetic aperture radar (InSAR). The pattern of surface uplift is consistent with distributed shear above the propagating tip of a reverse fault, leading to both an emergent fault and folding proximal to the rupture. This study demonstrates the potential for quantifying modest (<1 m) vertical deformation on a reverse fault using optical satellite imagery.

Geophysical Research Letters↗

USGS2019

Vertical coseismic offsets from differential high-resolution stereogrammetric DSMs: The 2013 Baluchistan, Pakistan earthquake

The recent proliferation of high-resolution (< 3-m spatial resolution) digital topography datasets opens a spectrum of geodetic applications in differential topography, including the quantification of coseismic vertical displacement fields. Most investigations of coseismic vertical displacements to date rely, in part, on pre- or post-event lidar surveys that are intractable or non-existent in many locales. Stereogrammetric digital surface models (DSMs) derived from high-resolution satellite optical imagery provide a new avenue for the retrieval of spatially-dense vertical coseismic displacements on a global scale. In this study, we generated 2-m resolution pre- and post-seismic DSMs from satellite optical imagery spanning the 2013 Mw7.7 Baluchistan strike-slip earthquake that occurred on the Hoshab fault in southern Pakistan. We applied the Iterative Closest Point algorithm to the DSMs to quantify the coseismic vertical displacement field at a spatial resolution of 10-30 m and to generate 3D coseismic strain tensors. We found that across-fault vertical offsets alternated between uplift and subsidence and varied between ~1-3 m in a non-systematic manner along the Hoshab fault. We show that the pre-existing topography and near-fault geomorphology are variably consistent and inconsistent with the displacement kinematics of the 2013 earthquake, and we argue that these relationships highlight varied slip sense history along the Hoshab fault. Notably, topography along the southern extents of the Hoshab fault requires different surface displacement kinematics than occurred in the 2013 earthquake, suggesting that the Hoshab fault accommodates varying senses of slip (bimodal slip) through time.

Baluchistan Province↗

USGS2019

Four major Holocene earthquakes on the Reelfoot fault recorded by sackungen in the New Madrid seismic zone, USA

Three sequences of well-documented, major ~M7+ earthquakes (1811-1812 CE, ~1450 CE, and ~900 CE) in the New Madrid seismic zone, USA, contribute significantly to seismic hazard in the region. However, it is unknown whether this <550 yr recurrence interval has been constant throughout the Holocene given limited geomorphic evidence of prior earthquakes. We extend the record of paleoearthquakes along the Reelfoot fault via investigation of ridge-top gravitational failure features, interpreted as sackungen. The sackungen occur in bluffs along the eastern margin of the Mississippi River floodplain and are concentrated near (<15 km) the southwest-dipping Reelfoot reverse fault. A paleoseismic trench excavated across sackungen at the Paw Paw site exposed four packages of colluvial sediment that postdate 30-11 ka Peoria loess. We interpret the colluvial packages to have been deposited following episodic failure of the sackungen as a result of strong ground motions from the following sequence of earthquakes: event 4, 1640 ± 1730 BCE; event 3, 340 ± 670 CE; event 2, 1430 ± 380 CE; and event 1, 1810 ± 50 CE (2-sigma). Event timing corresponds to previously documented earthquakes and represents the longest archive of paleoearthquakes on the Reelfoot fault. If the trenched sackungen record all major Reelfoot fault earthquakes, our observations in combination with prior investigations indicate a period of quiescence from at least 11 – 4.7 ka, followed by four major seismic events culminating in the 1811-1812 CE sequence. This clustered earthquake recurrence helps place bounds on seismic-hazard and geodynamic models in the New Madrid seismic zone.

Arkansas, Missouri↗

USGS2019

A constant slip rate for the western Qilian Shan frontal thrust during the last 200 ka consistent with GPS-derived and geological shortening rates

Active thrust faulting at the front of the Qilian Shan accommodates the northeastward growth of the Tibetan Plateau, however, the lifespan of individual faults and their slip history on different timescales remain largely unknown. Here, we show that the main range-bounding thrust fault of the western Qilian Shan has accrued tectonic slip at an almost constant rate during the last ∼200 ka, and possibly since fault initiation in the mid-Miocene. Our finding is based on 10Be ages from a flight of five deformed fluvial terraces along the Hongshuiba river, which constrain the vertical slip rate of the Qilian Shan frontal thrust to be m/ka during the last 200 ka. With a fault dip of ° constrained by seismic reflection data, we obtain a horizontal shortening rate of m/ka. This value is consistent with both the short-term shortening rate of mm/a derived from GPS data and the long-term shortening rate of m/ka, which is based on a balanced geological cross-section. The latter provides a total shortening estimate of km since the thrust fault initiated Ma ago. The agreement between the shortening rates on the range of timescales between 100 and 107 years suggests that the western Qilian Shan frontal thrust has slipped at a steady rate since its initiation and implies that this fault is the main structure responsible for the growth of the western Qilian Shan.

Tibetan Plateau, Qilian Shan↗

USGS2019

Preliminary report on engineering and geological effects of the July 2019 Ridgecrest earthquake sequence

The Ridgecrest Earthquake sequence included a foreshock event on July 4 2019 (M6.4) and a M7.1 mainshock event on July 5 2019. These events occurred in the Eastern California Shear Zone, near Indian Wells Valley, south of China Lake and west of Searles Valley. GEER has partnered with several organizations to collect perishable data and document the important impacts of these events, including the US Geological Survey, the California Geological Survey, the US Navy, the Southern California Earthquake Center, and local utilities. Critical geotechnical features of this event are extensive left-lateral (M6.4 event) and right-lateral (M7.1 event) surface ruptures over fault segments of variable complexity and width as well as across extensional and compressive step-over zones. We also document lifeline performance at fault crossings (gas, water, electrical), mainshock slip and afterslip, liquefaction and lateral spreading features, and liquefaction effects on structures. These effects are documented using field (ground) mapping and aerial imagery that will support subsequent development of high-resolution digital elevation models. Over 750 ground motions were recorded from the foreshock and mainshock alone, with many additional aftershock records. The data demonstrate significant impacts of site response and rupture directivity on ground motion attributes.

Report↗

USGS2018

Coseismic sackungen in the New Madrid seismic zone, USA

High‐resolution lidar reveals newly recognized evidence of strong shaking in the New Madrid seismic zone in the central United States. We mapped concentrations of sackungen (ridgetop spreading features) on bluffs along the eastern Mississippi River valley in northwestern Tennessee that likely form or are reactivated during large earthquakes. These sackungen are concentrated on the hanging wall of the Reelfoot reverse fault and show a preferential orientation indicating ground failure normal to fault strike. These observations suggest that the sackungen record one or more earthquakes on the southern Reelfoot fault since the deposition of the ~30‐ to 11‐ka Peoria Loess and potentially constrain the minimum intensity of near‐fault ground motion. This study demonstrates that sackungen can be used to infer fault source and mechanism and, in combination with field‐based techniques, improve paleoseismic records and seismic hazard models.

New Madrid Seismic Zone↗

USGS2018

Paleoseismic results from the Alpine site, Wasatch fault zone: Timing and displacement data for six holocene earthquakes at the Salt Lake City–Provo segment boundary

To improve the characterization of Holocene earthquakes on the Wasatch fault zone (WFZ), we conducted light detection and ranging (lidar)‐based neotectonic mapping and excavated a paleoseismic trench across an 8‐m‐high fault scarp near Alpine, Utah, located < 1    km "> < 1 km <1 km south of the boundary between the Salt Lake City and Provo segments (SLCS and PS). We document evidence for six paleoearthquakes (AL6–AL1) from scarp‐derived colluvial wedges and crosscutting relations. A ground‐penetrating radar survey across the scarp resolved fault‐zone width, but not paleoearthquake stratigraphy. Bayesian (OxCal) modeling of 13 radiocarbon and 13 optically stimulated luminescence ages indicates that six earthquakes occurred ∼ 6.2 – 0.4    ka "> ∼ 6.2 – 0.4 ka ∼6.2–0.4 ka . Interseismic recurrence ranges from 0.2 to 1.8 ky (mean 1.2 ky). We estimate 6.5 ± 0.7    m "> 6.5 ± 0.7 m 6.5±0.7 m of cumulative vertical tectonic displacement across the ≥ 14 ‐ m ‐ wide "> ≥ 14 ‐ m ‐ wide ≥14‐m‐wide fault zone used near‐field observations of scarp‐derived colluvial‐wedge thicknesses, antithetic faulting, and graben backtilting. This is similar to our independent estimate of 6.5 ± 0.5    m "> 6.5 ± 0.5 m 6.5±0.5 m using far‐field observations of the offset ground surface and correlation of alluvial‐fan stratigraphy across the WFZ. These results suggest that colluvial‐wedge thickness at the Alpine site approximates one‐half the original west‐facing fault scarp height. Per‐event vertical displacements range from 0.8 to 1.2 m (mean 1.1 m), which we use to estimate surface rupture lengths that may exceed 50 km from earthquakes as large as moment magnitude ( M w "> M w Mw ) ∼ 7.0 "> ∼ 7.0 ∼7.0 . The late Holocene average vertical slip rate is 0.9    mm / yr "> 0.9 mm / yr 0.9 mm/yr ( 0.7 – 1.2    mm / yr "> 0.7 – 1.2 mm / yr 0.7–1.2 mm/yr range). Earthquake frequency has increased in the past ∼ 1    ky "> ∼ 1 ky ∼1 ky , whereas displacement per event has been similar for the past ∼ 6    ky "> ∼ 6 ky ∼6 ky , suggesting that strain accumulation is not the sole factor that controls the frequency and size of earthquakes at the SLCS–PS segment boundary. These findings can be used for a more nuanced characterization of earthquakes at the SLCS–PS boundary and improve earthquake hazard assessments along the Wasatch Front.

Utah↗

USGS2018

Modeling the Holocene slip history of the Wasatch fault (Utah): Coseismic and postseismic Coulomb stress changes and implications for paleoseismicity and seismic hazard

The Wasatch fault zone defines the eastern boundary of the actively extending Basin and Range Province (Utah, western United States) and poses a significant seismic hazard to the metropolitan areas along the Wasatch Range. A wealth of paleoseismological data documents ∼24 surface-rupturing Mw ≥ 7 earthquakes along the Wasatch fault during the past 6400 yr. Here, we simulated the Holocene earthquake sequence on the Wasatch, Oquirrh−Great Salt Lake, and West Valley faults using three-dimensional finite-element forward modeling with the goal to calculate coseismic and postseismic Coulomb stress changes and to evaluate the slip and magnitude of hypothetical present-day and future earthquakes. Our results show that a good fit between modeled and observed paleoevents and time-integrated slip rates can be achieved within the uncertainties of the paleoseismological record and model parameters like the fault geometry. The Coulomb stress change analysis for selected paleoearthquakes showed that maximum positive stress changes are induced on faults located along strike of the source fault, while faults parallel to the source fault are generally located in stress shadow zones. Postseismic viscoelastic relaxation considerably modifies the coseismic stress changes; the resulting transient stress changes are recognizable for more than 100 yr after an earthquake. The modeled present-day state of Coulomb stress changes shows that the Brigham City, Salt Lake City, and Provo segments of the Wasatch fault are prone to failure in a Mw ≥ 6.8 earthquake. Our study shows that simulation of an entire earthquake sequence based on a paleoseismological record is feasible and facilitates identification of possible gaps and inconsistencies in the paleoseismological record. Therefore, forward modeling of earthquake sequences may ultimately contribute to improved seismic hazard estimates.

Utah↗

USGS2018

Combining conflicting Bayesian models to develop paleoseismic records—An example from the Wasatch Fault Zone, Utah

Bayesian statistical analyses of paleoseismic data result in the probabilistic determination of earthquake times using geochronological data evaluated in the context of a stratigraphic model. However, a fundamental problem in paleoseismology is how to use the Bayesian approach to model sparse and/or conflicting geochronological datasets, such as those derived from sites exhibiting episodic sedimentary and pedogenic processes in moderate‐ to high‐energy environments (e.g., a normal‐faulted alluvial fan). Using paleoseismic data for the Corner Canyon site on the Salt Lake City segment of the Wasatch fault zone (Utah), we develop an approach by which multiple Bayesian models are combined to generate an earthquake history at a site. This approach accommodates mutually exclusive interpretations of the geochronological data and thereby limits the influence of sparse data, stratigraphically inconsistent ages, or a single, subjective model interpretation. For the Corner Canyon site, we integrate four OxCal Bayesian models to generate a chronology of six events between ∼ 4.8 "> ∼ 4.8 and ∼ 0.5    ka "> ∼ 0.5 ka ⁠ . Late Holocene (post‐5 ka) mean recurrence and vertical slip‐rate estimates are ∼ 0.9    ky "> ∼ 0.9 ky (0.7–1.0 ky; 95% confidence) and 1.1    mm / yr "> 1.1 mm / yr ( ⁠ 0.8 – 1.7    mm / yr "> 0.8 – 1.7 mm / yr range), respectively. Although our method increases the uncertainty in the timing of individual earthquakes, it more objectively accounts for potential geochronological errors and different interpretations of stratigraphic age control. By relaxing the need to select a single age model, our approach yields more accurate earthquake‐timing results that will better facilitate evaluations of along‐fault event correlation and earthquake rupture length.

Utah↗

USGS2018

An updated stress map of the continental U.S. reveals heterogeneous intraplate stress

Knowledge of the state of stress in the Earth’s crust is key to understanding the forces and processes responsible for earthquakes. Historically, low rates of natural seismicity in the central and eastern United States have complicated efforts to understand intraplate stress, but recent improvements in seismic networks and the spread of human-induced seismicity have greatly improved data coverage. Here we compile a nationwide stress map based on formal inversions of focal mechanisms that challenges the idea that deformation in continental interiors is driven primarily by broad, uniform stress fields derived from distant plate boundaries. Despite plate-boundary compression, extension dominates roughly half of the continent, and second-order forces related to lithospheric structure appear to control extension directions. We also show that the states of stress in several active eastern United States seismic zones differ significantly from surrounding areas and that these anomalies cannot be explained by transient processes, suggesting that earthquakes are focused by persistent, locally-derived sources of stress. Such spatially variable intraplate stress appears to justify the current, spatially variable estimates of seismic hazard. Future work to quantify sources of stress, stressing-rate magnitudes, and their relationship with strain and earthquake rates could allow prospective mapping of intraplate hazard.

Continental United States↗

USGS2018

Integrate urban‐scale seismic hazard analyses with the U.S. National Seismic Hazard Model

For more than 20 yrs, damage patterns and instrumental recordings have highlighted the influence of the local 3D geologic structure on earthquake ground motions (e.g., M "> M M 6.7 Northridge, California, Gao et al. , 1996 ; M "> M M 6.9 Kobe, Japan, Kawase, 1996 ; M "> M M 6.8 Nisqually, Washington, Frankel, Carver, and Williams, 2002 ). Although this and other local‐scale features are critical to improving seismic hazard forecasts, historically they have not been explicitly incorporated into the U.S. National Seismic Hazard Model (NSHM, national model and maps), primarily because the necessary basin maps and methodologies were not available at the national scale. Instead, the U.S. Geological Survey (USGS), its partners, and external groups developed urban seismic hazard maps (urban models and maps) that consider detailed site effects in local areas (e.g., Wong et al. , 2002 ; Cramer et al. , 2006 ; Frankel et al. , 2007 ; Graves et al. , 2011 ). The disconnect between the urban and national hazard models, however, means that the national models, which underlie U.S. building codes and other applications, do not make use of all of the scientific results informing earthquake ground‐shaking hazards. We recommend that future U.S. national seismic hazard assessment make use of all available regional information, including that in urban models. In this column, we describe the roles of and differences between the urban and national models, and discuss the obstacles to and benefits of integrating the urban models with the national model. Future progress on this issue will require further coordination and implementation efforts between the USGS and external partners.

Seismological Research Letters↗

USGS2017

Apparent late Quaternary fault slip rate increase in the southwestern Lower Rhine Graben, central Europe

In regions of low strain, long earthquake recurrence intervals (10 4 –10 6 yrs) and erosive processes limit preservation of Quaternary markers suitable for distinguishing whether faults slip at uniform or secularly varying rates. The Lower Rhine graben in the border region of Germany, The Netherlands, and Belgium provides a unique opportunity to explore Quaternary slip‐rate variations in a region of low strain using the basal (2.29±0.29 Ma) and surface (700±80 ka) contacts of the regionally extensive main terrace (“Hauptterrasse”), deposited by the Rhine and Maas Rivers. These surfaces are vertically offset 3–140 m and 0–68 m, respectively, across individual fault strands within a distributed network of northwest‐trending, slow‐slipping (<0.1 mm/yr) normal faults. In this investigation, we construct Quaternary slip histories for the southern Lower Rhine graben faults using new main terrace surface vertical offset measurements made from light detection and ranging (lidar)‐derived bare‐earth digital terrain models, which we synthesize with existing constraints on the offset basal contact of this fluvial deposit ( n =91 collocated sites with displacement constraints). We find that >80% of the sites record an apparent increase in slip rate for the more recent interval from 700 ka to present, which corresponds to a period of increased uplift of the nearby Rhenish Massif and regional volcanism. However, the apparent increase in slip rate could result, in part, from erosion of the footwall surface below the main terrace, leading to an apparent displacement that is smaller than the total vertical offset since the start of the Quaternary. Prior work focused on characterization of these faults as seismic sources in the Lower Rhine graben has preferentially relied on the average fault‐slip rate constrained using the base of the main terrace. We suggest that average fault‐slip rates calculated using the ∼700 ka main terrace surface are subjected to fewer uncertainties and sample a time interval that is more relevant for seismic‐hazard analysis.

Lower Rhine graben↗

USGS2017

Tectonic summaries of magnitude 7 and greater earthquakes from 2000 to 2015

This paper describes the tectonic summaries for all magnitude 7 and larger earthquakes in the period 2000–2015, as produced by the U.S. Geological Survey National Earthquake Information Center during their routine response operations to global earthquakes. The goal of such summaries is to provide important event-specific information to the public rapidly and concisely, such that recent earthquakes can be understood within a global and regional seismotectonic framework. We compile these summaries here to provide a long-term archive for this information, and so that the variability in tectonic setting and earthquake history from region to region, and sometimes within a given region, can be more clearly understood.

Open-File Report↗

USGS2017

Pulsed strain release on the Altyn Tagh fault, northwest China

Earthquake recurrence models assume that major surface-rupturing earthquakes are followed by periods of reduced rupture probability as stress rebuilds. Although purely periodic, time- or slip-predictable rupture models are known to be oversimplifications, a paucity of long records of fault slip clouds understanding of fault behavior and earthquake recurrence over multiple ruptures. Here, we report a 16 kyr history of fault slip—including a pulse of accelerated slip from 6.4 to 6.0 ka—determined using a Monte Carlo analysis of well-dated offset landforms along the central Altyn Tagh strike-slip fault (ATF) in northwest China. This pulse punctuates a median rate of 8.1 +1.2 / −0.9 mm/a and likely resulted from either a flurry of temporally clustered ∼Mw 7.5 ground-rupturing earthquakes or a single large >Mw 8.2 earthquake. The clustered earthquake scenario implies rapid re-rupture of a fault reach >195 km long and indicates decoupled rates of elastic strain energy accumulation versus dissipation, conceptualized as a crustal stress battery. If the pulse reflects a single event, slip-magnitude scaling implies that it ruptured much of the ATF with slip similar to, or exceeding, the largest documented historical ruptures. Both scenarios indicate fault rupture behavior that deviates from classic time- or slip-predictable models.

Altyn Tagh fault↗

USGS2017

Middle and Late Pleistocene glaciations in the southwestern Pamir and their effects on topography

Glacial chronologies provide insight into the evolution of paleo-landscapes, paleoclimate, topography, and the erosion processes that shape mountain ranges. In the Pamir of Central Asia, glacial morphologies and deposits indicate extensive past glaciations, whose timing and extent remain poorly constrained. Geomorphic data and 15 new 10 Be exposure ages from moraine boulders and roches moutonnées in the southwestern Pamir document multiple Pleistocene glacial stages. The oldest exposure ages, , underestimate the age of the earliest preserved glacial advance and imply that the modern relief of the southwestern Pamir (peaks at ∼5000–6000 m a.s.l.; valleys at ∼2000–3000 m a.s.l.) already existed in the late Middle Pleistocene. Younger exposure ages (∼40–80 ka, ∼30 ka) complement the existing Central Asian glacial chronology and reflect successively less extensive Late Pleistocene glaciations. The topography of the Pamir and the glacial chronologies suggest that, in the Middle Pleistocene, an ice cap or ice field occupied the eastern Pamir high-altitude plateau, whereas westward flowing valley glaciers incised the southwestern Pamir. Since the Late Pleistocene deglaciation, the rivers of the southwestern Pamir adjusted to the glacially shaped landscape. Localized rapid fluvial incision and drainage network reorganization reflect the transient nature of the deglaciated landscape.

Earth and Planetary Science Letters↗

USGS2017

Refining fault slip rates using multiple displaced terrace risers-An example from the Honey Lake fault, NE California, USA

Faulted terrace risers are semi-planar features commonly used to constrain Quaternary slip rates along strike-slip faults . These landforms are difficult to date directly and therefore their ages are commonly bracketed by age estimates of the adjacent upper and lower terrace surfaces. However, substantial differences in the ages of the upper and lower terrace surfaces (a factor of 2.4 difference observed globally) produce large uncertainties in the slip-rate estimate. In this investigation, we explore how the full range of displacements and bounding ages from multiple faulted terrace risers can be combined to yield a more accurate fault slip rate. We use 0.25-m cell size digital terrain models derived from airborne lidar data to analyze three sites where terrace risers are offset right-laterally by the Honey Lake fault in NE California, USA. We use ages for locally extensive subhorizontal surfaces to bracket the time of riser formation: an upper surface is the bed of abandoned Lake Lahontan having an age of 15.8 ± 0.6 ka and a lower surface is a fluvial terrace abandoned at 4.7 ± 0.1 ka. We estimate lateral offsets of the risers ranging between 6.6 and 28.3 m (median values), a greater than fourfold difference in values. The amount of offset corresponds to the riser's position relative to modern stream meanders: the smallest offset is in a meander cutbank position, whereas the larger offsets are in straight channel or meander point-bar positions. Taken in isolation, the individual terrace-riser offsets yield slip rates ranging from 0.3 to 7.1 mm/a. However, when the offset values are collectively assessed in a probabilistic framework, we find that a uniform (linear) slip rate of 1.6 mm/a (1.4–1.9 mm/a at 95% confidence) can satisfy the data, within their respective uncertainties. This investigation demonstrates that integrating observations of multiple offset elements (crest, midpoint, and base) from numerous faulted and dated terrace risers at closely spaced sites can refine slip-rate estimates on strike-slip faults.

California↗