skip to main content


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
NSF-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. Abstract The Wrangell Arc in Alaska (USA) and adjacent volcanic fields in the Yukon provide a long-term record of interrelations between flat-slab subduction of the Yakutat microplate, strike-slip translation along the Denali–Totschunda–Duke River fault system, and magmatism focused within and proximal to a Cretaceous suture zone. Detrital zircon (DZ) U-Pb (n = 2640) and volcanic lithic (DARL) 40Ar/39Ar dates (n = 2771) from 30 modern river sediment samples document the spatial-temporal evolution of Wrangell Arc magmatism, which includes construction of some of the largest Quaternary volcanoes on Earth. Mismatches in DZ and DARL date distributions highlight the impact of variables such as mineral fertility and downstream mixing/dilution on resulting provenance signatures. Geochronologic data document the initiation of Wrangell Arc magmatism at ca. 30–17 Ma along both sides of the Totschunda fault on the north flank of the Wrangell–St. Elias Mountains in Alaska, followed by southeastward progression of magmatism at ca. 17–10 Ma along the Duke River fault in the Yukon. This spatial-temporal evolution is attributable to dextral translation along intra-arc, strike-slip faults and a change in the geometry of the subducting slab (slab curling/steepening). Magmatism then progressed generally westward outboard of the Totschunda and Duke River faults at ca. 13–6 Ma along the southern flank of the Wrangell–St. Elias Mountains in Alaska and then northwestward from ca. 6 Ma to present in the western Wrangell Mountains. The 13 Ma to present spatial-temporal evolution is consistent with dextral translation along intra-arc, strike-slip faults and previously documented changes in plate boundary conditions, which include an increase in plate convergence rate and angle at ca. 6 Ma. Voluminous magmatism is attributed to shallow subduction-related flux melting and slab edge melting that is driven by asthenospheric upwelling along the lateral edge of the Yakutat flat slab. Magmatism was persistently focused within or adjacent to a remnant suture zone, which indicates that upper plate crustal heterogeneities influenced arc magmatism. Rivers sampled also yield subordinate Paleozoic–Mesozoic DZ and DARL age populations that reflect earlier episodes of magmatism within underlying accreted terranes and match magmatic flare-ups documented along the Cordilleran margin. 
    more » « less
  2. The Alaska Range suture zone exposes Cretaceous to Quaternary marine and nonmarine sedimentary and volcanic rocks sandwiched between oceanic rocks of the accreted Wrangellia composite terrane to the south and older continental terranes to the north. New U-Pb zircon ages, 40Ar/39Ar, ZHe, and AFT cooling ages, geochemical compositions, and geological field observations from these rocks provide improved constraints on the timing of Cretaceous to Miocene magmatism, sedimentation, and deformation within the collisional suture zone. Our results bear on the unclear displacement history of the seismically active Denali fault, which bisects the suture zone. Newly identified tuffs north of the Denali fault in sedimentary strata of the Cantwell Formation yield ca. 72 to ca. 68 Ma U-Pb zircon ages. Lavas sampled south of the Denali fault yield ca. 69 Ma 40Ar/39Ar ages and geochemical compositions typical of arc assemblages, ranging from basalt-andesite-trachyte, relatively high-K, and high concentrations of incompatible elements attributed to slab contribution (e.g., high Cs, Ba, and Th). The Late Cretaceous lavas and bentonites, together with regionally extensive coeval calc-alkaline plutons, record arc magmatism during contractional deformation and metamorphism within the suture zone. Latest Cretaceous volcanic and sedimentary strata are locally overlain by Eocene Teklanika Formation volcanic rocks with geochemical compositions transitional between arc and intraplate affinity. New detrital-zircon data from the modern Teklanika River indicate peak Teklanika volcanism at ca. 57 Ma, which is also reflected in zircon Pb loss in Cantwell Formation bentonites. Teklanika Formation volcanism may reflect hypothesized slab break-off and a Paleocene–Eocene period of a transform margin configuration. Mafic dike swarms were emplaced along the Denali fault from ca. 38 to ca. 25 Ma based on new 40Ar/39Ar ages. Diking along the Denali fault may have been localized by strike-slip extension following a change in direction of the subducting oceanic plate beneath southern Alaska from N-NE to NW at ca. 46–40 Ma. Diking represents the last recorded episode of significant magmatism in the central and eastern Alaska Range, including along the Denali fault. Two tectonic models may explain emplacement of more primitive and less extensive Eocene–Oligocene magmas: delamination of the Late Cretaceous–Paleocene arc root and/or thickened suture zone lithosphere, or a slab window created during possible Paleocene slab break-off. Fluvial strata exposed just south of the Denali fault in the central Alaska Range record synorogenic sedimentation coeval with diking and inferred strike-slip displacement. Deposition occurred ca. 29 Ma based on palynomorphs and the youngest detrital zircons. U-Pb detrital-zircon geochronology and clast compositional data indicate the fluvial strata were derived from sedimentary and igneous bedrock presently exposed within the Alaska Range, including Cretaceous sources presently exposed on the opposite (north) side of the fault. The provenance data may indicate ~150 km or more of dextral offset of the ca. 29 Ma strata from inferred sediment sources, but different amounts of slip are feasible. Together, the dike swarms and fluvial strata are interpreted to record Oligocene strike-slip movement along the Denali fault system, coeval with strike-slip basin development along other segments of the fault. Diking and sedimentation occurred just prior to the onset of rapid and persistent exhumation ca. 25 Ma across the Alaska Range. This phase of reactivation of the suture zone is interpreted to reflect the translation along and convergence of southern Alaska across the Denali fault driven by highly coupled flat-slab subduction of the Yakutat microplate, which continues to accrete to the southern margin of Alaska. Furthermore, a change in Pacific plate direction and velocity at ca. 25 Ma created a more convergent regime along the apex of the Denali fault curve, likely contributing to the shutting off of near-fault extension- facilitated arc magmatism along this section of the fault system and increased exhumation rates. 
    more » « less
  3. The uplift history of the Sierra Nevada, California, is a topic of long-standing disagreement with much of it centered on the timing and nature of slip along the range-bounding normal fault along the east flank of the southern Sierra Nevada. The history of normal fault slip is important for characterizing the uplift history of the Sierra Nevada, as well as for characterizing the geologic and geodynamic factors that drove, and continue to drive, normal faulting. To address these issues, we completed new structural studies and extensive apatite (U-Th)/He (AHe) thermochronometry on samples collected from three vertical transects in the footwall to the east-dipping southern Sierra Nevada normal fault (SNNF). Our structural studies on bedrock fault planes show that the SNNF is a steeply (~70°) east-dipping normal fault. The new AHe data reveal two elevation-invariant AHe age arrays, indicative of two distinct periods of cooling and exhumation, which we interpret as initiation of normal faulting along the SNNF at ca. 28–27 Ma with a second phase of normal faulting at ca. 17–13 Ma. We argue that beginning in the late Oligocene, the SNNF marked the now long-standing stable western limit, or break-away zone, of the Basin and Range. Slip along SNNF, and the associated unloading of the footwall, likely resulted in two periods of uplift of Sierra Nevada during the late Cenozoic. Trench retreat, driven by westward motion of the North American plate, along the Farallon–North American subduction zone boundary, as well as the gravitationally unstable northern and southern Basin and Range pushing on the cold Sierra Nevada, likely drove the late Oligocene- aged normal slip along the SNNF and the similar-aged but generally local and minor extension within the Basin and Range. We posit that the thick proto–Basin and Range lithosphere was primed for late Oligocene extension by replacement of the steepening Farallon slab with hot and buoyant asthenosphere. While steepening of the Farallon slab had not yet reached the southern Sierra Nevada by late Oligocene time, we speculate that late Oligocene slip along the SNNF reactivated a late Cretaceous dextral shear zone as the Sierra Nevada block was pulled and pushed westward in response to trench retreat and gravitational potential energy. The dominant middle Miocene normal fault-slip history along the SNNF is contemporaneous with high-magnitude slip recorded along range-bounding normal faults across the Basin and Range, including the east-adjacent Inyo and White mountains, indicating that this period of extension was a major regional tectonic event. We infer that a combination of slab-driven trench retreat along the Juan de Fuca–North America subduction zone boundary and clockwise rotation of the southern ancestral Cascade Range superimposed on continental lithosphere pre-conditioned for extension drove this episode of middle Miocene normal slip along the SNNF and extension to the east across the Basin and Range. Transtensional plate motion along the Pacific–North America plate boundary, and likely a growing slab window, continued to drive extension along the SNNF and the western Basin and Range, but not until ca. 11 Ma when the Mendocino triple junction reached the latitude of our northernmost (U-Th)/He transect. 
    more » « less
  4. Wöfler, Andreas (Ed.)
    Abstract Classically held mechanisms for removing mountain topography (e.g., erosion and gravitational collapse) require 10-100 Myr or more to completely remove tectonically generated relief. Here, we propose that mountain ranges can be completely and rapidly (<2 Myr) removed by a migrating hotspot. In western North America, multiple mountain ranges, including the Teton Range, terminate at the boundary with the relatively low relief track of the Yellowstone hotspot. This abrupt transition leads to a previously untested hypothesis that preexisting mountainous topography along the track has been erased. We integrate thermochronologic data collected from the footwall of the Teton fault with flexural-kinematic modeling and length-displacement scaling to show that the paleo-Teton fault and associated Teton Range was much longer (min. original length 190-210 km) than the present topographic expression of the range front (~65 km) and extended across the modern-day Yellowstone hotspot track. These analyses also indicate that the majority of fault displacement (min. 11.4-12.6 km) and the associated footwall mountain range growth had accumulated prior to Yellowstone encroachment at ~2 Ma, leading us to interpret that eastward migration of the Yellowstone hotspot relative to stable North America led to removal of the paleo-Teton mountain topography via posteruptive collapse of the range following multiple supercaldera (VEI 8) eruptions from 2.0 Ma to 600 ka and/or an isostatic collapse response, similar to ranges north of the Snake River plain. While this extremely rapid removal of mountain ranges and adjoining basins is probably relatively infrequent in the geologic record, it has important implications for continental physiography and topography over very short time spans. 
    more » « less
  5. We interpret the kinematics of the Tangra Yumco (TYC) rift by evaluating spatiotemporal trends in fault displacement, extension onset, and exhumation rates. We present new geologic mapping, U-Pb geochronology, zircon (U-Th)/He (ZHe) thermochronology, and HeFTy thermal modeling results that are critical to testing dynamic models of extension in Tibet. The TYC rift is bounded by two NNE striking (~N10°E-N35°E) high angle (~45-70°) active normal faults that alternate dominance along strike. Footwall granodiorites show foliation, slip lineation, and fault plane striation measurements indicative of northeast directed oblique sinistral-normal slip. In North and South TYC, hanging wall deposits are cut by a series of active high-angle normal faults which likely sole into a master fault at depth, while in central TYC, hanging wall deposits display synthetic graben structures potentially indicative of low-angle faulting. Analysis of ~50 samples collected across key structural relationships in and around TYC yield 14 mean U-Pb dates between ~59-49 Ma and ~190 single-grain ZHe dates between ~60-4 Ma with spatial trends in ZHe data correlating strongly with latitude. Samples from Gangdese latitudes show a concentration of ~28-15 Ma ages, while those north of ~29.8° latitude yield both younger (~9-4 Ma) and older (~59-45 Ma) ages. We interpret (1) Gangdese Range samples reflect exhumation during contraction and uplift along the GCT peaking at ~21-20 Ma, (2) ~9-4 Ma ages reveal extension timing along fault segments experiencing significant rift-related exhumation, and (3) ~59-45 Ma ages represent un-reset or partially-reset samples from fault segments that have experienced lesser magnitudes of rift exhumation. HeFTy thermal models indicate a two-stage cooling history with initial slow cooling followed by accelerated cooling rates in Late Miocene-Pliocene time (~13-4 Ma) consistent with prior results from TYC and other Tibetan rifts. Our data are consistent with a segment linkage fault evolution model for the TYC rift, with underthrusting of Indian lithosphere likely related to the northward acceleration of rifting. Future work will utilize advanced HeFTy modeling including U-Pb and apatite fission track data to further constrain the exhumation history of TYC and test dynamic models of extension for southern Tibet. 
    more » « less