Geology topics

Jens-Erik Lundstern

Publications and source records attributed to Jens-Erik Lundstern.

8 recordsLinked to original sources

USGS2024

Recent advances in characterizing the crustal stress field and future applications of stress data: Perspectives from North America

The stress field controls patterns of crustal deformation, including which faults are likeliest to cause earthquakes or transmit fluids. Since the 1950s, maps of maximum horizontal stress ( S Hmax ) orientations have advanced dramatically, and the style of faulting (relative principal stress magnitudes) has recently been mapped in some regions as well. This perspectives paper summarizes developments in characterizing stress orientations and (relative) magnitudes, including new seismic and borehole methods, as well as progress in identifying the causes of stress variations. Despite these advances, adding far more spatiotemporal detail would allow geoscientists to address many of today's key challenges regarding natural hazards, energy development, and geodynamics. In particular, it is critically important to characterize stress heterogeneity at multiple scales while also recognizing the coherent variability of the stress field. The second part of the paper considers how more detailed stress datasets could prove essential to addressing some of the grand questions in geoscience, including deciphering the poorly understood feedbacks between crustal dynamics and surface processes, improving earthquake and eruption forecasts, and determining the origins and shared properties of plate boundaries.

Geological Society, London, Special Publications↗

USGS2024

Paleogene sedimentary basin development in southern Nevada, USA

The cause of the transition from Mesozoic and early Cenozoic crustal shortening to later extension in the western United States is debated. In many parts of the extant Sevier hinterland, now the Basin and Range Province, the sedimentary sections that provide the most direct record of that transition remain poorly studied and lack meaningful age control. In this paper, we present field characterization supported by U-Pb detrital zircon and 40 Ar/ 39 Ar feldspar ages for ten sections across southern Nevada. We describe a newly identified basin, here named the Fallout Hills basin, which preserves >1.0 km of sedimentary deposits as old as middle Eocene, ca. 48 Ma. Deposition occurred during the 20 m.y. (million years) before the 27.60 ± 0.03 Ma Monotony Tuff blanketed much of south-central Nevada, based on 47.6 Ma and younger detrital zircon maximum depositional ages (MDAs) from near the Pintwater and Spotted Ranges. Elsewhere in southern Nevada, prevolcanic Cenozoic strata commonly form thinner (~100 m), isolated exposures that yield detrital zircon MDAs ≤10 m.y. older than (and in some cases essentially the same age as) the ca. 27–28 Ma ignimbrites that cap the sections. A variable but overall upward-fining facies pattern is observed in both the Fallout Hills basin and the thinner sections. These localized patterns imply topographic changes that are unlikely to reflect plate-scale processes and are not consistent with large-magnitude extension. Instead, variable uplift due to magmatism combined with antecedent topographic relief from thrust faulting and subsequent erosion likely provided accommodation for these deposits.

Nevada↗

USGS2024

The importance of nodal plane orientation diversity for earthquake focal mechanism stress inversions

Inversions of earthquake focal mechanisms are among the most accessible and reliable methods for determining crustal stress. However, the use of this method varies widely, and assumptions that underpin it are often violated, potentially compromising stress estimates. We investigate the consequences of violating the little-studied assumption that the focal mechanisms have diverse orientations. Our approach is to employ data-informed synthetic mechanisms, with nodal plane orientations defined by recent earthquake lineaments in the Midland Basin, western Texas, and rakes consistent with slip in the mapped stress field. Using both the traditional stress inversion method that assumes constant shear stress magnitudes on the causative faults as well as a recently published variable shear stress method, we show that low fault plane diversity can cause maximum horizontal stress ( S Hmax ) orientation and relative principal stress magnitude (faulting regime) estimates to differ markedly from the true values. This problem is compounded for catalogs with even modest amounts of noise (≤15°) or few (e.g., 20) mechanisms. Significantly, traditional approaches for quantifying uncertainty such as the bootstrap can severely underestimate the true uncertainty under these circumstances. To remedy this, we provide simple tools to quantify nodal plane orientation diversity and stress inversion reliability.

Geological Society of London Special Publications↗

USGS2022

Revised age and regional correlations of Cenozoic strata on Bat Mountain, Death Valley region, California, USA, from zircon U-Pb geochronology of sandstones and ash-fall tuffs

Basin analysis and tectonic reconstructions of the Cenozoic history of the Death Valley region, California, USA, are hindered by a lack of volcanic (tuff) age control in many stratigraphic successions exposed in the Grapevine and Funeral Mountains of California, USA. Although maximum depositional ages (MDAs) interpreted from detrital zircon U-Pb data may be a promising alternative to volcanic ages, arguments remain regarding the calculation of MDAs including, but not limited to, the number of “young” grains to consider (i.e., the spectrum of dates used to calculate the MDA); which grains, if any, should be ignored; which approaches yield results that are statistically rigorous; and ultimately, which approaches result in ages that are geologically reasonable. We compare commonly used metrics of detrital zircon MDA for five sandstone samples from the Cenozoic strata exposed on Bat Mountain in the southern Funeral Mountains of California—i.e., the youngest single grain (YSG), the weighted mean of the youngest grain cluster of two or more grains at 1σ uncertainty (YC1σ(2+)) and of three or more grains at 2σ uncertainty (YC2σ(3+)), the youngest graphical peak (YPP), and the maximum likelihood age (MLA). Every sandstone sample yielded abundant Cenozoic zircon U-Pb dates that formed unimodal, near-normal age distributions that were clearly distinguishable from the next-oldest grains in each sample and showed an apparent up-section decrease in peak age. Benchmarked against published K/Ar and 40 Ar/ 39 Ar ages and five new zircon U-Pb ages of ash-fall tuffs, our analysis parallels prior studies and demonstrates that many MDA metrics—YSG, YC1σ(2+), YC2σ(3+), and YPP—drift toward unreasonably young or old values. In contrast, the maximum likelihood estimation approach and the resulting MLA metric consistently produce geologically appropriate estimates of MDA without arbitrary omission of any young (or old) zircon dates. Using the MLAs of sandstones and zircon U-Pb ages of interbedded ash-fall tuffs, we develop a new age model for the Oligocene–Miocene Amargosa Valley Formation (deposited ca. 28.5–18.5 Ma) and the Miocene Bat Mountain Formation (deposited ca. 15.5–13.5 Ma) and revise correlations to Cenozoic strata across the eastern Death Valley region.

California↗

USGS2022

Magmatism, migrating topography, and the transition from Sevier shortening to Basin and Range extension, western United States

The paleogeographic evolution of the western U.S. Great Basin from the Late Cretaceous to the Cenozoic is critical to understanding how the North American Cordillera at this latitude transitioned from Mesozoic shortening to Cenozoic extension. According to a widely applied model, Cenozoic extension was driven by collapse of elevated crust supported by crustal thicknesses that were potentially double the present ~30–35 km. This model is difficult to reconcile with more recent estimates of moderate regional extension (≤50%) and the discovery that most high-angle, Basin and Range faults slipped rapidly ca. 17 Ma, tens of millions of years after crustal thickening occurred. Here, we integrated new and existing geochronology and geologic mapping in the Elko area of northeast Nevada, one of the few places in the Great Basin with substantial exposures of Paleogene strata. We improved the age control for strata that have been targeted for studies of regional paleoelevation and paleoclimate across this critical time span. In addition, a regional compilation of the ages of material within a network of middle Cenozoic paleodrainages that developed across the Great Basin shows that the age of basal paleovalley fill decreases southward roughly synchronous with voluminous ignimbrite flareup volcanism that swept south across the region ca. 45–20 Ma. Integrating these data sets with the regional record of faulting, sedimentation, erosion, and magmatism, we suggest that volcanism was accompanied by an elevation increase that disrupted drainage systems and shifted the continental divide east into central Nevada from its Late Cretaceous location along the Sierra Nevada arc. The north-south Eocene–Oligocene drainage divide defined by mapping of paleovalleys may thus have evolved as a dynamic feature that propagated southward with magmatism. Despite some local faulting, the northern Great Basin became a vast, elevated volcanic tableland that persisted until dissection by Basin and Range faulting that began ca. 21–17 Ma. Based on this more detailed geologic framework, it is unlikely that Basin and Range extension was driven by Cretaceous crustal overthickening; rather, preexisting crustal structure was just one of several factors that that led to Basin and Range faulting after ca. 17 Ma—in addition to thermal weakening of the crust associated with Cenozoic magmatism, thermally supported elevation, and changing boundary conditions. Because these causal factors evolved long after crustal thickening ended, during final removal and fragmentation of the shallowly subducting Farallon slab, they are compatible with normal-thickness (~45–50 km) crust beneath the Great Basin prior to extension and do not require development of a strongly elevated, Altiplano-like region during Mesozoic shortening.

Nevada↗

USGS2022

Downhill from Austin and Ely to Las Vegas: U-Pb detrital zircon suites from the Eocene–Oligocene Titus Canyon Formation and associated strata, Death Valley, California

In a reconnaissance investigation aimed at interrogating the changing topography and paleogeography of the western United States prior to Basin and Range faulting, a preliminary study made use of U-Pb ages of detrital zircon suites from 16 samples from the Eocene–Oligocene Titus Canyon Formation, its overlying units, and correlatives near Death Valley. The Titus Canyon Formation unconformably overlies Neoproterozoic to Devonian strata in the Funeral and Grapevine Mountains of California and Nevada. Samples were collected from (1) the type area in Titus Canyon, (2) the headwaters of Monarch Canyon, and (3) unnamed Cenozoic strata exposed in a klippe of the Boundary Canyon fault in the central Funeral Mountains. Red beds and conglomerates at the base of the Titus Canyon Formation at locations 1 and 2, which contain previously reported 38–37 Ma fossils, yielded mostly Sierran batholith–age detrital zircons (defined by Triassic, Jurassic, and Cretaceous peaks). Overlying channelized fluvial sandstones, conglomerates, and minor lacustrine shale, marl, and limestone record an abrupt change in source region around 38–36 Ma or slightly later, from more local, Sierran arc–derived sediment to extraregional sources to the north. Clasts of red radiolarian-bearing chert, dark radiolarian chert, and quartzite indicate sources in the region of the Golconda and Roberts Mountains allochthons of northern Nevada. Sandstones intercalated with conglomerate contain increasing proportions of Cenozoic zircon sourced from south-migrating, caldera-forming eruptions at the latitude of Austin and Ely in Nevada with maximum depositional ages (MDAs) ranging from 36 to 24 Ma at the top of the Titus Canyon Formation. Carbonate clasts and ash-rich horizons become more prevalent in the overlying conglomeratic Panuga Formation (which contains a previously dated 15.7 Ma ash-flow tuff). The base of the higher, ash-dominated Wahguyhe Formation yielded a MDA of 14.4 Ma. The central Funeral Mountains section exposes a different sequence of units that, based on new data, are correlative to the Titus Canyon, Panuga, and Wahguyhe Formations at locations 1 and 2. An ash-flow tuff above its (unexposed) base provided a MDA of 34 Ma, and the youngest sample yielded a MDA of 12.7 Ma. The striking differences between age-correlative sections, together with map-based evidence for channelization, indicate that the Titus Canyon Formation and overlying units likely represent fluvial channel, floodplain, and lacustrine deposits as sediments mostly bypassed the region, moving south toward the Paleogene shoreline in the Mojave Desert. The profound changes in source regions and sedimentary facies documented in the Titus Canyon Formation took place during ignimbrite flareup magmatism and a proposed eastward shift of the continental divide from the axis of the Cretaceous arc to a new divide in central Nevada in response to thermal uplift and addition of magma to the crust. This uplift initiated south-flowing fluvial systems that supplied sediments to the Titus Canyon Formation and higher units.

Arizona, California, Nevada, Utah↗

USGS2022

State of stress in areas of active unconventional oil and gas development in North America

In this paper, we present comprehensive data on stress orientation and relative magnitude in areas throughout North America where unconventional oil and gas are currently being developed. We find excellent agreement between maximum horizontal principal stress (S Hmax ) orientations over a wide range of depths, using multiple methods. In all basins studied, we observed coherent stress fields that in some cases vary systematically from one part of a basin to another. In the Appalachian Basin in the eastern United States, S Hmax is oriented northeast–southwest to east-northeast–west-southwest and the style of faulting is compressive, transitioning from reverse faulting in eastern Pennsylvania and New York to principally strike-slip faulting in western Pennsylvania, Ohio, and West Virginia. In the midcontinent, central Oklahoma is characterized by an approximately east–west S Hmax direction and strike-slip faulting. The Fort Worth Basin in northeastern Texas is characterized by normal–strike-slip faulting and a north-northeast–south-southwest S Hmax direction. In the Midland subbasin of western Texas, S Hmax is consistently approximately east–west and normal–strike-slip faulting is observed. Farther west, the Delaware subbasin of western Texas and southeastern New Mexico is characterized by normal faulting and S Hmax rotates ∼150° clockwise from north to south. Marked changes in S Hmax direction also occur across the Raton Basin of southern Colorado and northern New Mexico, the Denver-Julesburg Basin in northern Colorado, and the Uinta Basin in northeastern Utah, likely associated with their location near the margins of extensional provinces. The new data sets we present help improve operational efficiency by constraining absolute stress magnitudes and the ideal azimuth to drill horizontal wells (i.e., perpendicular to the local S Hmax orientation) and make it possible to predict which fractures and faults are likely to be activated during hydraulic stimulation.

AAPG Bulletin↗

USGS2021

Pore pressure threshold and fault slip potential for induced earthquakes in the Dallas-Fort Worth area of north central Texas

Earthquakes were induced in the Fort Worth Basin from 2008 through 2020 by increase in pore pressure from injection of oilfield wastewater (SWD). In this region and elsewhere, a missing link in understanding the mechanics of causation has been a lack of comprehensive models of pore pressure evolution (ΔPp) from SWD. We integrate detailed earthquake catalogs, ΔPp, and probabilistic fault slip potential (FSP) and find that faults near large-scale SWD operations became unstable early, when ΔPp reached ∼0.31 MPa and FSP reached 0.24. Faults farther from SWD became unstable later, when FSP reached 0.17 and at much smaller ΔPp. Earthquake sequences reactivated with mean ΔPp of ∼0.05 MPa. The response of faults shows strong variability, with many remaining stable at higher ΔPp and few that became seismogenic at smaller changes. As ΔPp spread regionally, an ever-increasing number of faults were impacted and the most sensitive became unstable.

Texas↗