Geology topics

Daniel F. Stockli

Publications and source records attributed to Daniel F. Stockli.

8 recordsLinked to original sources

USGS2026

Depositional interplay between the Ancestral Rocky Mountains and Ouachita–Marathon–Sonora orogenies: Insights from provenance records in the late Paleozoic Marfa Basin, West Texas, U.S.A.

The Marfa Basin in West Texas is a late Palaeozoic synorogenic depocenter associated with regional deformation linked to the Ancestral Rocky Mountains (ARM) and Ouachita–Marathon–Sonora (OMS) orogenies in southwestern Laurentia. Basin strata range in age from Middle Pennsylvanian to the middle Permian and include the Cieneguita, Alta, Pinto Canyon, Rose Mine and Mina Grande Formations. Sandstone petrography and detrital zircon (DZ) U–Pb and (U–Th)/He double dating data from these strata reveal three tectonically driven sedimentation stages: syntectonic ARM deposition, progressive OMS foredeep deposition and an orogenic transition. The Cieneguita and lower part of the Alta Formations exhibit a Mesoproterozoic DZ age signature (~1318 and ~1076 Ma age peaks) and quartzo-feldspathic sandstone compositions sourced from the adjacent ARM-related Diablo Platform basement uplift in the Middle Pennsylvanian to earliest Permian. In contrast, the upper part of the Alta Formation, as well as the Pinto Canyon and Rose Mine Formations, have peri-Gondwanan DZ age signatures, with Mesoproterozoic (~1069–1036 Ma age peaks), Neoproterozoic–Cambrian (~700–490 Ma) and Palaeozoic (~490–300 Ma) age modes and litho-quartzose sandstone compositions derived from the OMS fold-and-thrust belt and orogenic hinterland during the early to middle Permian. The lower to middle parts of the Alta Formation have alternating DZ age signatures and sandstone compositions from both ARM and OMS sources, revealing that the transition in the sediment supply occurred during the middle Wolfcampian. This transition was not characterised by source mixing, but rather by sediment interfingering alternately sourced from the Diablo Platform uplift and the advancing OMS belt. These observations are confirmed by the DZ He ages, which reveal distinct cooling histories for both source terranes. These results document a switch from ARM- to OMS-related syntectonic deposition in southwestern Laurentia during the early Permian, demonstrating that ARM-driven deformation largely preceded the continental collision along the Marathon segment of the OMS orogen.

Texas↗

USGS2024

Tracking cycles of Phanerozoic opening and closing of ocean basins using detrital rutile and zircon geochronology and geochemistry

Sedimentary basins provide a deep time archive of tectonic and Earth-surface processes that can be leveraged by detrital mineral U-Pb dating and geochemistry to track paleogeography, magmatism, and crustal evolution. Zircon preserves the long-term (billions of years) record of supercontinent cycles; however, it is biased toward preserving felsic crustal records. Detrital rutile complements the detrital zircon record by providing constraints on the time and temperature of rifting and mafic magmatism, metamorphism, exhumation of the middle and lower crust, subduction, and amagmatic orogenesis. We use detrital zircon U-Pb and detrital rutile U-Pb geochronology and trace element analysis of Permian to Eocene siliciclastic rocks in the southern Pyrenees to capture supercontinent cycles of ocean basins opening and closing. Detrital rutile age spectra show peaks at ca. 100 Ma associated with rifting and hyperextension in the Pyrenean realm, 200 Ma associated with the Central Atlantic Magmatic Province, and 330 Ma, 375 Ma, and 400 Ma associated with subduction and Rheic Ocean crust formation. Zr-in-rutile thermometry and rutile Cr-Nb systematics provide further insight into metamorphic facies (peak metamorphic temperatures) and source rock lithology (mafic versus felsic affinity). Detrital zircon age spectra have peaks at ca. 300 Ma, 450 Ma, and 600 Ma associated with major orogenic events and felsic magmatism, and Th/U ratios provide information on relative zircon formation temperatures. Comparison of these independent records shows that detrital rutile reflects rifting, magma-poor orogenesis, and oceanic lithospheric processes, while detrital zircon detects continental lithospheric processes. Integrated detrital zircon and rutile data sets archive past geological events across multiple Wilson cycles.

Geology↗

USGS2023

Evolution of Miocene normal and dextral faulting in the lower Colorado River region near Blythe, California, USA

The evolution of strain in nascent continental plate boundaries commonly involves distributed deformation and transitions between different styles of deformation as the plate boundary matures. Distributed NW-striking faults, many with km-scale right-lateral separation, are prevalent near Blythe, California, and have been variably interpreted to have accommodated either Middle Miocene NE-SW extension as normal faults or Late Miocene to Pliocene dextral shear as strike-slip faults. However, with poor timing and kinematic constraints, it is unclear how these faults relate to known domains of Neogene deformation and the evolution of the Pacific–North America plate boundary. We present kinematic data (n = 642 fault planes, n = 512 slickenlines) that demonstrate that these faults dominantly dip steeply northeast; ~96% of measured faults record normal, dextral, or oblique dextral-normal kinematics that likely reflect a gradational transition between normal and dextral oblique kinematic regimes. We constrain fault timing with 11.7 Ma and 7.0 Ma 40 Ar/ 39 Ar dates of rocks cut by faults, and laser ablation–inductively coupled plasma–mass spectrometry U-Pb dating of calcite mineralized during oblique dextral faulting that demonstrates fault slip at ca. 10–7 Ma and perhaps as late as ca. 4 Ma. This Late Miocene dextral oblique faulting is best compatible with a documented regional transition from Early to Middle Miocene NE-directed extension during detachment fault slip to subsequent NW-directed dextral shear. We estimate 11–38 km of cumulative dextral slip occurred across a 50-km-wide zone from the Palen to Riverside mountains, including up to 20 km of newly documented dextral shear that may partly alleviate the regional discrepancy of cumulative dextral shear along this part of the Late Miocene Pacific–North America plate boundary.

California↗

USGS2022

Tectonic influence on axial-transverse sediment routing in the Denver Basin

Detrital zircon U-Pb and (U-Th)/He ages from latest Cretaceous–Eocene strata of the Denver Basin provide novel insights into evolving sediment sourcing, recycling, and dispersal patterns during deposition in an intracontinental foreland basin. In total, 2464 U-Pb and 78 (U-Th)/He analyses of detrital zircons from 21 sandstone samples are presented from outcrop and drill core in the proximal and distal portions of the Denver Basin. Upper Cretaceous samples that predate uplift of the southern Front Range during the Laramide orogeny (Pierre Shale, Fox Hills Sandstone, and Laramie Formation) contain prominent Late Cretaceous (84–77 Ma), Jurassic (169–163 Ma), and Proterozoic (1.69–1.68 Ga) U-Pb ages, along with less abundant Paleozoic through Archean zircon grain ages. These grain ages are consistent with sources in the western U.S. Cordillera, including the Mesozoic Cordilleran magmatic arc and Yavapai-Mazatzal basement, with lesser contributions of Grenville and Appalachian zircon recycled from older sedimentary sequences. Mesozoic zircon (U-Th)/He ages confirm Cordilleran sources and/or recycling from the Sevier orogenic hinterland. Five of the 11 samples from syn-Laramide basin fill (latest Cretaceous–Paleocene D1 Sequence) and all five samples from the overlying Eocene D2 Sequence are dominated by 1.1–1.05 Ga zircon ages that are interpreted to reflect local derivation from the ca. 1.1 Ga Pikes Peak batholith. Corresponding late Mesoproterozoic to early Neoproterozoic zircon (U-Th)/He ages are consistent with local sourcing from the southern Front Range that underwent limited Mesozoic–Cenozoic unroofing. The other six samples from the D1 Sequence yielded detrital zircon U-Pb ages similar to pre-Laramide units, with major U-Pb age peaks at ca. 1.7 and 1.4 Ga but lacking the 1.1 Ga age peak found in the other syn-Laramide samples. One of these samples yielded abundant Mesozoic and Paleozoic (U-Th)/He ages, including prominent Early and Late Cretaceous peaks. We propose that fill of the Denver Basin represents the interplay between locally derived sediment delivered by transverse drainages that emanated from the southern Front Range and a previously unrecognized, possibly extraregional, axial-fluvial system. Transverse alluvial-fluvial fans, preserved in proximal basin fill, record progressive unroofing of southern Front Range basement during D1 and D2 Sequence deposition. Deposits of the upper and lower D1 Sequence across the basin were derived from these fans that emanated from the southern Front Range. However, the finer-grained, middle portion of the D1 Sequence that spans the Cretaceous-Paleogene boundary was deposited by both transverse (proximal basin fill) and axial (distal basin fill) fluvial systems that exhibit contrasting provenance signatures. Although both tectonic and climatic controls likely influenced the stratigraphic development of the Denver Basin, the migration of locally derived fans toward and then away from the thrust front suggests that uplift of the southern Front Range may have peaked at approximately the Cretaceous-Paleogene boundary.

Colorado↗

USGS2021

Detrital signals of coastal erosion and fluvial sediment supply during glacio-eustatic sea-level rise, Southern California, USA

Coastal erosion, including sea-cliff retreat, represents both an important component of some sediment budgets and a significant threat to coastal communities in the face of rising sea level. Despite the importance of predicting future rates of coastal erosion, few prehistoric constraints exist on the relative importance of sediment supplied by coastal erosion versus rivers with respect to past sea-level change. We used detrital zircon U-Pb geochronology as a provenance tracer of river and deep-sea fan deposits from the Southern California Borderland (United States) to estimate relative sediment contributions from rivers and coastal erosion from late Pleistocene to present. Mixture modeling of submarine canyon and fan samples indicates that detrital zircon was dominantly (55%–86%) supplied from coastal erosion during latest Pleistocene (ca. 13 ka) sea-level rise, with lesser contributions from rivers, on the basis of unique U-Pb age modes relative to local Peninsular Ranges bedrock sources. However, sediment that was deposited when sea level was stable at its highest and lowest points since the Last Glacial Maximum was dominantly supplied by rivers, suggesting decreased coastal erosion during periods of sea-level stability. We find that relative sediment supply from coastal erosion is strongly dependent on climate state, corroborating predictions of enhanced coastal erosion during future sea-level rise.

California↗

USGS2021

Extreme Quaternary plate boundary exhumation and strike slip localized along the southern Fairweather fault, Alaska, USA

The Fairweather fault (southeastern Alaska, USA) is Earth’s fastest-slipping intracontinental strike-slip fault, but its long-term role in localizing Yakutat–(Pacific–)North America plate motion is poorly constrained. This plate boundary fault transitions northward from pure strike slip to transpression where it comes onshore and undergoes a <25°, 30-km-long restraining double bend. To the east, apatite (U-Th)/He (AHe) ages indicate that North America exhumation rates increase stepwise from ~0.7 to 1.7 km/m.y. across the bend. In contrast, to the west, AHe age-depth data indicate that extremely rapid 5–10 km/m.y. Yakutat exhumation rates are localized within the bend. Further northwest, Yakutat AHe and zircon (U-Th)/He (ZHe) ages gradually increase from 0.3 to 2.6 Ma over 150 km and depict an interval of extremely rapid >6–8 km/m.y. exhumation rates that increases in age away from the bend. We interpret this migration of rapid, transient exhumation to reflect prolonged advection of the Cenozoic–Cretaceous sedimentary cover of the eastern Yakutat microplate through a stationary restraining bend along the edge of the North America plate. Yakutat cooling ages imply a long-term strike-slip rate (54 ± 6 km/m.y.) that mimics the millennial (53 ± 5 m/k.y.) and decadal (46 mm/yr) rates. Fairweather fault slip can account for all Pacific–North America relative plate motion throughout Quaternary time and indicates stability of highly localized plate boundary strike slip on a single fault where extreme rock uplift rates are persistently localized within a restraining bend.

Alaska↗

USGS2013

New thermochronometric constraints on the Tertiary landscape evolution of the central and eastern Grand Canyon, Arizona

Thermal histories are modeled from new apatite (U-Th)/He and apatite fission-track data in order to quantitatively constrain the landscape evolution of the Grand Canyon region. Fifty new samples and their associated thermochronometric ages are presented here. Samples span from Lee’s Ferry in the east to Quartermaster Canyon in the west and include four age-elevation transects into Grand Canyon and borehole samples from the Coconino Plateau. Twenty-seven samples are inversely modeled to provide continuous thermal histories. This represents the most extensive and complete dataset on patterns of long-term exhumation in the Grand Canyon region, and it enables us to constrain the timing and magnitude of erosion and also discriminate between canyon incision and broader planation. The new data suggest that the early Cenozoic landscape in eastern Grand Canyon was low in relief and does not indicate the presence of an early Cenozoic precursor to the modern Grand Canyon. However, there is evidence for the incision of a smaller-scale canyon across the Kaibab Uplift at 28–20 Ma. This middle-Cenozoic denudation event was accompanied by the removal of a majority of remaining Mesozoic strata west of the Kaibab Uplift. In contrast, just upstream in the area of Lee’s Ferry, ∼2 km of Mesozoic strata remained over the middle Cenozoic and were removed after 10 Ma.

Arizona↗

USGS2011

Rapid cooling rates at an active mid-ocean ridge from zircon thermochronology

Oceanic spreading ridges are Earth's most productive crust generating environment, but mechanisms and rates of crustal accretion and heat loss are debated. Existing observations on cooling rates are ambiguous regarding the prevalence of conductive vs. convective cooling of lower oceanic crust. Here, we report the discovery and dating of zircon in mid-ocean ridge dacite lavas that constrain magmatic differentiation and cooling rates at an active spreading center. Dacitic lavas erupted on the southern Cleft segment of the Juan de Fuca ridge, an intermediate-rate spreading center, near the intersection with the Blanco transform fault. Their U–Th zircon crystallization ages (29.3 − 4.6 + 4.8 ka; 1σ standard error s.e.) overlap with the (U–Th)/He zircon eruption age (32.7 ± 1.6 ka) within uncertainty. Based on similar 238 U− 230 Th disequilibria between southern Cleft dacite glass separates and young mid-ocean ridge basalt (MORB) erupted nearby, differentiation must have occurred rapidly, within ~ 10–20 ka at most. Ti-in-zircon thermometry indicates crystallization at 850–900 °C and pressures > 70–150 MPa are calculated from H 2 O solubility models. These time-temperature constraints translate into a magma cooling rate of ~ 2 × 10 − 2 °C/a. This rate is at least one order-of-magnitude faster than those calculated for zircon-bearing plutonic rocks from slow spreading ridges. Such short intervals for differentiation and cooling can only be resolved through uranium-series ( 238 U– 230 Th) decay in young lavas, and are best explained by dissipating heat convectively at high crustal permeability.

Juan de Fuca ridge and the Blanco transform fault ↗