skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


  • NSF Public Access
  • Search Results
  • Quantifying Normal Fault Evolution from River Profile Analysis in the Northern Basin and Range Province, Southwest Montana, USA
Title: Quantifying Normal Fault Evolution from River Profile Analysis in the Northern Basin and Range Province, Southwest Montana, USA
Abstract Over the past few decades, tectonic geomorphology has been widely implemented to constrain spatial and temporal patterns of fault slip, especially where existing geologic or geodetic data are poor. We apply this practice along the eastern margin of Bull Mountain, Southwest Montana, where 15 transient channels are eroding into the flat, upstream relict landscape in response to an ongoing period of increased base level fall along the Western North Boulder fault. We aim to improve constraints on the spatial and temporal slip rates across the Western North Boulder fault zone by applying channel morphometrics, cosmogenic erosion rates, bedrock characteristics, and calibrated reproductions of the modern river profiles using a 1-dimensional stream power incision model that undergoes a change in the rate of base level fall. We perform over 104 base level fall simulations to explore a wide range of fault slip dynamics and stream power parameters. Our best fit simulations suggest that the Western North Boulder fault started as individual fault segments along the middle to southern regions of Bull Mountain that nucleated around 6.2 to 2.5 Ma, respectively. This was followed by the nucleation of fault segments in the northern region around 1.5 to 0.4 Ma. We recreate the evolution of the Western North Boulder fault to show that through time, these individual segments propagate at the fault tips and link together to span over 40 km, with a maximum slip of 462 m in the central portion of the fault. Fault slip rates range from 0.02 to 0.45 mm/yr along strike and are consistent with estimates for other active faults in the region. We find that the timing of fault initiation coincides well with the migration of the Yellowstone hotspot across the nearby Idaho-Montana border and thus attribute the initiation of extension to the crustal bulge from the migrating hotspot. Overall, we provide the first quantitative constraints on fault initiation and evolution of the Western North Boulder fault, perhaps the farthest north basin in the Northern Basin and Range province that such constraints exist. We show that river profiles are powerful tools for documenting the spatial and temporal patterns of normal fault evolution, especially where other geologic/geodetic methods are limited, proving to be a vital tool for accurate tectonic hazard assessments.  more » « less
Award ID(s):
1727139
PAR ID:
10311498
Author(s) / Creator(s):
; ; ; ;
Editor(s):
Valla, Pierre
Date Published:
Journal Name:
Lithosphere
Volume:
2021
Issue:
1
ISSN:
1941-8264
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. (, Geomorphology)
    Advancements in the quantitative investigation of fluvial topography in tectonically active regions account for the emerging numerical modelling of river profiles and their linear inversions. Applications of a geomorphic approach strive for the reconstruction of long-term tectonic deformation histories by decoding base-level fall signals transiently embedded in the geomorphic record. I present integrated results from river profile inversions and marine terraces analyses, here used to outline the deformation model associated to a debated crustal fault system at the southern termination of the Calabrian Forearc High (Central Mediterranean). The study aims at constraining spatial and temporal variations in geometry, strain partitioning, slip rate, and time-transgressive propagation of the tectonic deformation associated with a fault system. In particular, I systematically analyze river profiles draining the eastern flank of the Peloritani Mts. in northeastern Sicily (southern Italy), a NNE-SSW- trending mountain ridge thought as located at the footwall of an active, ESE-dipping normal fault. I perform the linear inversions of fluvial topography by applying recently available MATLAB scripts, and I carry out the analysis of terraced surfaces by both GIS tools and MATLAB-based software packages. The results I obtain suggest that the eastern flank of the Peloritani Mts. have been deformed according to at least three main stages of uplift accommodated along distinct, ~10–15 km-long en-échelon arranged fault segments. The reconstructed evolution of tectonic deformation unveils time-transgressive, southward propagation since the last ~600 kyr, inset in a general increase, through time, of the regional component of uplift. The results of this study contribute to address the issues of the deformation style and strain partitioning along complex and/or debated fault systems. These results also demonstrate the potential of the geomorphic approach in defining the spatial and temporal tectonic evolutionary model of a region. 
    more » « less
  2. ; (, Earth and Space Science)
    Abstract Geologic and geodetic observations provide constraints on tectonic and earthquake cycle kinematics. Block models offer one approach to integrating the effects of plate rotations, elastic strain accumulation, applied basal displacements, internal block strain, and idealized pressure sources. Here, we describe the construction of block models where spatially variable slip rates are parameterized by distance‐weighted eigenmodes operating over meshes of triangular dislocation elements. This dimensionally reduced model is recast as a quadratic programming problem with upper and lower bounds on both geologic fault slip rates and spatially variable slip deficit rates. We propose iterating over successive quadratic programming estimates with evolving slip rate bounds to find a solution consistent with specified coupling at all points on geometrically complex fault surfaces. 
    more » « less
  3. ; ; ; ; ; ; ; (, Journal of Geophysical Research: Solid Earth)
    Abstract The potential for future earthquakes on faults is often inferred from inversions of geodetically derived surface velocities for locking on faults using kinematic models such as block models. This can be challenging in complex deforming zones with many closely spaced faults or where deformation is not readily described with block motions. Furthermore, surface strain rates are more directly related to coupling on faults than surface velocities. We present a methodology for estimating slip deficit rate directly from strain rate and apply it to New Zealand for the purpose of incorporating geodetic data in the 2022 revision of the New Zealand National Seismic Hazard Model. The strain rate inversions imply slightly higher slip deficit rates than the preferred geologic slip rates on sections of the major strike‐slip systems including the Alpine Fault, the Marlborough Fault System and the northern part of the North Island Fault System. Slip deficit rates are significantly lower than even the lowest geologic estimates on some strike‐slip faults in the southern North Island Fault System near Wellington. Over the entire plate boundary, geodetic slip deficit rates are systematically higher than geologic slip rates for faults slipping less than one mm/yr but lower on average for faults with slip rates between about 5 and 25 mm/yr. We show that 70%–80% of the total strain rate field can be attributed to elastic strain due to fault coupling. The remaining 20%–30% shows systematic spatial patterns of strain rate style that is often consistent with local geologic style of faulting. 
    more » « less
  4. (, Seismological Research Letters)
    Abstract Seismic hazard assessment, such as the U.S. Geological Survey (USGS) National Seismic Hazard Model (NSHM), relies on estimates of fault slip rate based on geology and/or geodetic observations such as the Global Navigation Satellite System (GNSS), including the Global Positioning System. Geodetic fault slip rates may be estimated within a 3D spherical block model, in which the crust is divided into microplates bounded by mapped faults; fault slip rates are determined by the relative rotations of adjacent microplates. Uncertainty in selecting appropriate block-bounding faults and in forming closed microplates has limited the interpretability of block models for seismic hazard modeling. By introducing an automated block closure algorithm and regularizing the resulting densely spaced block model with total variation regularization, I develop the densest and most complete block model of the western continental United States to date. The model includes 853 blocks bounded by 1017 geologically identified fault sections from the USGS NSHM Fault Sections database. Microplate rotations and fault slip rates are constrained by 4979 GNSS velocities and 1243 geologic slip rates. I identify a regularized solution that fits the GNSS velocity field with a root mean square misfit of 1.9 mm/yr and reproduces 57% of geologic slip rates within reported geologic uncertainty and model sensitivity, consistent with other geodetic-based models in this Focus Section. This block model includes slip on faults that are not included in the USGS NSHM Fault sections database (but are required to form closed blocks) for an estimate of “off-fault” deformation of 3.62×1019  N·m/yr, 56% of the total calculated moment accumulation rate in the model. 
    more » « less
  5. (, Journal of Geophysical Research: Solid Earth)
    Abstract I employ an elasticity‐based method to invert a geodetically derived surface velocity field in the western US using for present‐day surface strain rate fields with uncertainties. The method uses distributed body forces in a thin elastic sheet and allows for discontinuities in velocity across creeping faults using the solution for dislocations in a thin elastic plate. I compare the strain rate fields with previously published stress orientations and moment rates from geological slip rate data and previous geodetic studies. Geologic and geodetic moment rates are calculated using slip rate and off‐fault strain rates from the 2023 US National Seismic Hazard Model (NSHM) deformation models. I find that computed total geodetic moment rates are higher than NSHM summed moment rates on faults for all regions of the western US except the highest deforming rate regions including the Western Transverse Ranges and the northern and southern San Andreas Fault (SAF) system in California. Computed geodetic moment rates are comparable to the moment rates derived from the geodetically based NSHM deformation models in all regions. I find systematic differences in orientations of maximum horizontal shortening rate and maximum horizontal compressive stress in the Pacific Northwest region and along much of the SAF system. In the Pacific Northwest, the maximum horizontal stress orientations are rotated counterclockwise 40–90° relative to the maximum horizontal strain rate directions. Along the SAF system, the maximum horizontal stresses are rotated systematically 25–40° clockwise (closer to fault normal) relative to the strain rates. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email