SEARCH · USGS Science

Results for “Rapid Communications in Mass Spectrometry”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,602 records · Page 82Linked to original sources

USGS2024

3-D mapping of the conterminous U.S. within the USGS National Cooperative Geologic Mapping Program: Progress and future prospects

The U.S. Geological Survey (USGS) National Cooperative Mapping Program (NCGMP) is bringing together subsurface and three-dimensional information at multiple scales for the conterminous United States from data produced throughout the USGS and by federal and state partners. Components of this work include data inventory and catalog development, data integration and database development, and construction of subregional and basin-scale three-dimensional (3-D) geologic models. At small scales, the NCGMP is compiling subsurface data to create a limited number of subsurface layers that map the majority of the lower 48 States. At larger scales, 3-D models of stratigraphic units are intended to be the subsurface analog and extension of the surface geologic maps made within the Program. The Mapping Program is developing methods of storing, visualizing, and distributing 3-D models and subsurface data from multiple sources in easily shareable, queryable, non-proprietary format(s). We report on progress in 3-D mapping of the conterminous U.S. and on future prospects.

conterminous United States↗

USGS2024

Hypothetical CO2 leakage into, and hydrological plume management within, an underground source of drinking water at a proposed CO2 storage facility, Kemper County, Mississippi, USA

A large Geologic Carbon Sequestration (GCS) hub has been proposed in Kemper County, Mississippi. The target injection interval consists of numerous Cretaceous-aged deep saline aquifers overlain by a competent and extensive regional sealing layer. Above the seal, the deepest Underground Source of Drinking Water (USDW) at the site is the Eutaw aquifer of the Eutaw Group and McShan Formation, undifferentiated. To assess potential risks of leakage from the deep sequestration reservoir, a model of a portion of the Cretaceous Eutaw Group was constructed in this study. Simulations tested various permeabilities, hypothetical leakage rates, and plume mitigation strategies utilizing existing wells. Results suggest that, under the influence of regional groundwater flow fields, leaking CO 2 would effectively bypass the existing wells, and to influence this migration would require very large water extraction rates. Therefore, to ensure plume detection, monitoring for leakage at the injection wells themselves is very important.

Mississippi↗

USGS2025

Combining multisite tsunami and deformation modeling to constrain slip distributions for the 1700 C.E. Cascadia earthquake

A major earthquake ruptured the Cascadia subduction zone (CSZ) on 26 January 1700. Key paleoseismic evidence associated with this event include tsunami deposits, stratigraphic evidence of coastal coseismic subsidence, written Japanese records of a tsunami unaccompanied by earthquake shaking, and margin‐wide turbidites found offshore and in lacustrine environments. Despite this wealth of independent clues, important details about this event remain unresolved. Dating uncertainties do not conclusively establish whether the proxies are from one earthquake or a sequence of them, and we have limited knowledge of the likely slip distributions of the event or events. Here, we use a catalog of 37,500 candidate synthetic ruptures between M w 7.8 and 9.2 and simulate their resulting coseismic deformation and tsunami inundation. Each model is then compared against estimated Japan tsunami arrivals, regional coastal subsidence records, and local paleotsunami deposits mapped at six different coastal marshes and one coastal lake along the CSZ. We find that seven full‐margin ruptures with a median magnitude of M w 9.1 satisfy all three constraints. We favor one M w 9.11 model that best matches all site paleoseismic observations and suggests that the Cascadia megathrust slipped up to ∼30 m and must have shallow geodetic coupling. We also find that some sequences composed of three or four ruptures can still satisfy the observations, yet no sequences of two ruptures can. Sequences are differentiated into three groups based on whether they contain a mainshock rupture located in the south (>44° N) or further north. All sequences contain unruptured portions of the megathrust and most contain mainshocks with peak slip above 40 m. The fit of the geologic evidence from sequences is poor in comparison to single‐event models. Therefore, sequences are generally less favored compared to full‐margin events.

British Columbia, California, Oregon, Washington↗

USGS2026

Regional models for postfire debris-flow likelihood and rainfall thresholds across the western United States

The U.S. Geological Survey (USGS) uses an empirical model developed with logistic regression (the ‘M1’ model) to rapidly assess debris-flow likelihood and to identify quantitative rainfall thresholds for debris flows after wildfire in the western United States. The M1 model was calibrated to a debris-flow inventory from southern California (United States) and has been applied throughout the western United States. Limited spatial coverage in the calibration dataset has motivated evaluation of M1 model accuracy outside the calibration region (e.g., the Sierra Nevada or the eastern Cascade Range, United States). Previous test cases showed that M1 overpredicts debris-flow likelihood and underpredicts rainfall thresholds for some locations (e.g., Arizona, northern California, Colorado, New Mexico, United States). We sought to improve the regional applicability of a debris-flow likelihood model by expanding the debris-flow inventory used for calibration, testing multiple potential models and generating an updated model framework. The updated inventory includes 3788 observations from 67 burned areas paired with short duration rainfall ratios. The updated model framework consists of a modified model structure and sets of coefficients calibrated separately to the entire updated inventory and to subsets of the inventory that intersect three Environmental Protection Agency (EPA) Level 2 ecoregions (Mediterranean California, Upper Gila Mountains and Western Cordillera). Comparisons of predictions from the updated models with observed rainfall and debris-flow activity show that the updated models outperform the M1 model by ~15%–60% and improve the uniformity of predictive performance across the western United States. The updated models also reduce false positive rates relative to M1 and generate rainfall thresholds that are better aligned with relative differences in regional climatology and debris-flow activity.

Arizona, California, Colorado, Idaho, Montana, Nev↗

USGS2026

Hyperspectral narrowband imaging spectroscopy: A new paradigm for Earth observation

This editorial introduces the Special Issue entitled “Hyperspectral Narrowband Imaging Spectroscopy: A New Paradigm for Earth Observation” in the August 2026 issue of Photogrammetric Engineering & Remote Sensing (PE&RS), the flagship journal of the American Society for Photogrammetry and Remote Sensing (ASPRS). This volume represents the fourth dedicated hyperspectral special issue published in PE&RS, following earlier contributions by Thenkabail et al. (2025, 2024a, 2024b), and continues ASPRS’s commitment to advancing cutting‑edge imaging spectroscopy research and its applications across Earth system science. Remote sensing is undergoing one of the most profound transformations in its history. The emergence of hyperspectral narrowband (HNB) imaging spectroscopy data, capable of acquiring hundreds of contiguous, narrow spectral bands, has shifted the discipline from observing Earth in a handful of broad spectral windows to capturing continuous spectral signatures of the Planet. This transition marks a decisive break from the multispectral paradigm that has dominated satellite remote sensing for nearly five decades, driven by the advent of new orbital imaging spectrometers such as EnMAP, PRISMA, and NASA’s EMIT, and by the forthcoming Surface Biology and Geology (SBG) mission ((Pires Silva et al., 2026; Bourriz et al., 2025; Thenkabail et al., 2025; Chabrillat et al., 2024; Aneece et al., 2024; Dave et al., 2024; Thenkabail et al., 2024a; Thenkabail et al., 2024b; Thenkabail, 2024a; Thenkabail, 2024b; Thompson et al., 2022; Kokaly et al., 2022; Aneece & Thenkabail, 2022; Cawse Nicholson et al., 2021; Guanter et al., 2021; Vangi et al., 2021; Thenkabail et al., 2021). These missions (e.g., Table 1) deliver unprecedented spectral fidelity, improved signal to noise ratios, and global coverage capabilities, enabling a new era of quantitative, spectroscopy based Earth observation. Where multispectral broadbands (MBBs) provide only a few discrete measurements along the electromagnetic spectrum, HNB systems deliver rich, diagnostic information that enables scientists to characterize Earth’s surface with unprecedented biochemical, biophysical, and structural detail (Figure 1a, 1b). The implications for environmental monitoring, agriculture, water resources, and mineral exploration are profound. Several overarching themes emerge: • Spectral fidelity matters. The ability to preserve subtle absorption features is essential for mineral mapping, vegetation trait retrieval, and biochemical modeling. • AI and deep learning are indispensable. From destriping to classification, modern analytics must be scalable, label‑efficient, and capable of exploiting the full spectral–spatial richness of HNB data. • Physics‑based and data‑driven approaches must converge. Radiative transfer models such as PROSAIL, enhanced with localized soil parameterizations, remain foundational for biophysical retrievals and model‑based inference. • Dimensionality reduction and feature extraction are critical. Techniques such as L1‑ISOMAP demonstrate that intelligent manifold learning can unlock the structure of fused, high‑dimensional datasets. • Next‑generation architectures must be interactive and multimodal. ICTNet exemplifies the future of hyperspectral classification: hybrid, synergistic, and capable of modeling both local textures and global spectral dependencies.

Photogrammetric Engineering and Remote Sensing (PE↗

USGS2025

Reply to, “Comment on ‘The 1886 Charleston, South Carolina, earthquake: Relic railroad offset reveals rupture,’ by Roger Bilham and Susan E. Hough”

We welcome this opportunity to respond to Pratt et al. (2024) (hereinafter P24). Bilham and Hough (2023) proposed a “first-cut” elastic deformation model for the 1886 earthquake, a quantitative source model constrained by identified coseismic constraints. A key observation was the measurement of a lateral offset of a railroad line south of Summerville, leading to a model with predominately dextral slip and minor convergence, from which we concluded that active faulting had raised the Penholoway Marine Terrace >6 m since ∼770 ka. P24 questioned these constraints and proposed an alternative rupture model with predominantly reverse slip. This alternative model is neither consistent with coseismic constraints nor with other geophysical data. In a revised model presented here, we recognize that uplift of the Penholoway Terrace is confined to the eastern edge of the terrace, which we conclude results from active folding and tectonic transpression centered on the dextral fault that offset the railroad in 1886.

South Carolina↗

USGS2025

Ultralong, supershear rupture of the 2025 Mw 7.7 Mandalay earthquake reveals unaccounted risk

The 28 March 2025 moment magnitude ( M w ) 7.7 earthquake in Mandalay, Burma (Myanmar), ruptured 475 kilometers of the Sagaing Fault, which was more than twice the length predicted by magnitude scaling relationships. Kinematic slip models and observation of a Rayleigh Mach wave that passed through parts of Thailand confirmed that rupture occurred at supershear velocities of greater than 5 kilometers per second. The anomalous length exposed a vast population to violent near-fault shaking. The Mandalay earthquake is a modern analog for the M w 7.9 1906 San Francisco earthquake, another atypically long and fast rupture. Probabilistic seismic hazard analyses use scaling relations that do not account for such long ruptures at moderate magnitudes. This limitation, in conjunction with a likely increased population and infrastructure exposure for atypically long ruptures, contributes to a potential mischaracterization of seismic risk.

Mandalay↗

USGS2024

Event classification, seismicity, and eruption forecasting at Great Sitkin Volcano, Alaska: 1999–2023

The frequency content of volcanogenic seismicity is often used to classify events and their spatial and temporal progression is then used to map subsurface volcanic processes. The progression of volcano-seismic events and associated source processes also plays a critical role in eruption forecasting. Here we develop and evaluate a computerized methodology for characterizing volcano-seismic event types using Frequency Index and Average Peak Frequency. We apply and test this technique at Great Sitkin Volcano, Alaska, classifying over 9000 hypocenters between 1999 and 2023. This 24-year time span covers periods of seismic quiescence, earthquake activity on nearby tectonic (bookshelf) faults, precursory unrest from 2016 to 2021, and the explosive onset in May 2021 of the ongoing effusive eruption. We use the spatial and temporal evolution of classified event types to map the active volcanic and tectonic processes, develop a conceptual model of the subsurface magmatic system, and perform a retrospective analysis of eruption forecasts at Great Sitkin Volcano between 2016 and the present. The classification and progression of hypocenters suggests the subsurface Great Sitkin Volcano magmatic system consists of a mid- to lower- crustal source zone between 10 and 40 km depth and an upper crustal magma storage area between −1 and 10 km depth (hypocenter depth is referenced to sea level and negative depths reflect height above sea level). The earliest precursors occurred in July 2016 and consisted of deep long-period and volcano-tectonic earthquakes at mid-crustal depths suggesting the subsequent unrest and eruption were triggered by a deeper intrusion of magma. This mid-crustal seismic activity was immediately followed by the onset upper-crustal long-period events and volcano-tectonic earthquakes VTs suggesting a strong linkage between the shallow and deeper portions of the magmatic system. The upper crustal area was likely capped by the 1974 lava dome until the magmatic explosion on May 26, 2021.

Alaska↗

USGS2025

Zircon as a pathfinder to REE mineralization

Carbonatites and alkaline silicate rocks are major primary sources of the rare earth elements (REE) and other critical metals, such as Nb. Despite the economic significance of these rocks, their formation and the processes of REE enrichment are poorly understood. Here, statistical analysis of a global dataset demonstrates that zircon geochemistry is a powerful recorder of REE metallogenesis and a potential pathfinder for REE deposits. Zircons from REE and Nb fertile intrusions lack Eu anomalies and have elevated Gd/Yb and Th/Yb, indicating they crystallised from magmas that originated from deep, oxidised and enriched mantle sources. Complexes with Nb enrichment have low U/Nb, reflecting an enriched mantle source, whereas high U/Nb in REE-only fertile intrusions suggest a subduction-metasomatised mantle source. Machine learning models demonstrate high accuracy in classifying zircon from barren and fertile deposits. Classification of detrital zircons shows that REE-enriched deposits correlate with supercontinent assembly, whereas Nb fertile complexes are associated with supercontinent breakup. This approach offers a new, mineral to global scale, petrologic and exploration tool that enhances understanding of REE metallogenesis.

Geochemical Perspectives Letters↗

USGS2026

A ground-motion model derived using a generalized mean rupture distance for large slab interface earthquakes

Source–station distance is a central input to ground‐motion models (GMMs) for predicting seismic shaking. GMM development uses distance metrics including the Joyner–Boore distance, which is the shortest distance from an observation point to the surface projection of the earthquake rupture, and R rup the shortest distance to the rupture in three dimensions. Thompson and Baltay (2018) proposed the generalized mean rupture distance R p to address observed near‐fault ground‐motion saturation. R p accounts for the contribution to the shaking of all parts of the rupture and provides a simple method for incorporating spatially variable slip. They used R p to develop a GMM for shallow crustal earthquakes, assuming uniform slip. Here, we investigate the improvement offered by an R p ‐based GMM for large subduction interface earthquakes by recalibrating the path term for a published GMM ( Parker et al. , 2022 ) using R p derived from distributed slip models (DSMs). Inspection of within‐event and total residuals indicates that the recalibrated model fits the data at least as well as an alternative in which the path term was recalibrated using R rup and the same dataset. Incorporating slip information in its entirety or trimming the DSM geometry is preferable to assuming uniform moment release on a prescribed model fault that extends beyond the actual ruptured area. R p tuned to minimize model uncertainty is closer to the maximum distance to rupture for peak ground velocity (PGV) than acceleration, possibly reflecting the contribution to high‐frequency shaking of slip on local asperities. However, when the moment is concentrated far from stations, PGV is fit adequately with R p closer to R rup ⁠ . Our results suggest incorporating slip‐derived moment release through R p could improve GMMs for slab interface events, especially if the implementation of R p ‐based models is refined using a larger dataset of earthquakes with greater geographic diversity to account for regional ground‐motion variations.

Seismological Research Letters↗

USGS1976

Ground magnetic and VLF studies at Midnite uranium mine, Stevens County, Washington

Ground magnetic and E-mode VLF measurements were made at more than 200 stations at the Midnite uranium Mine, Stevens County, Washington. The equipment was found to be light and foolproof, and each set of data was collected in one day, exclusive of surveying and marking the station points. The maps of contoured magnetic and VLF data resembled geologic maps of the area, showing many common features, as well as a few unique features on each map. Each showed some evidence of a small normal fault cutting the region; taken together, the evidence for this fault was judged to be much stronger, so that further structural details were inferred. Because of soil cover and inadequate outcrops, the presence of this fault had hitherto been unsuspected. A target area for further uranium mineralization north of the fault was indicated by this quick, inexpensive geophysical survey. The drilling program which Midnite Mine management had underway at that time had already outlined a mineralized zone in the indicated area.

Washington↗

USGS2024

Field evidence and indicators of rockfall fragmentation and implications for mobility

Rockfall fragmentation can play an important role in hazard studies and the design of protective measures. However, the current lack of modeling tools that incorporate rock fragmentation mechanics is a limitation to enhancing studies and design. This research investigates the fragmentation patterns of rockfalls and analyzes the resulting distribution of fragment sizes within corresponding rockfall deposits. We focus on small rock fragments, which provide insights into the dynamics of the rockfall event and can be used as input for numerical modeling. We analyzed multiple rockfall events from locations worldwide, each exhibiting different degrees of fragmentation. Using image analysis techniques, we mapped all visible blocks, determined their volumes, and measured the distances they travelled from the initial point of impact. A key finding is the identification of three indicators of fragmentation. First, in cases where fragmentation was largely absent, we observed a trend of increasing block size with distance from the impact point or source area, which aligns with previously published findings. However, for energetic rockfall events characterized by intense fragmentation, we observed that small fragments exhibited longer travel distances compared to larger fragments. This distinction allowed us to differentiate blocks primarily resulting from the disaggregation process from those primarily resulting from dynamic fragmentation, with implications for rockfall mobility. Second, although the size distribution of rockfall deposits exhibits a power-law scaling for volumes larger than a minimum size threshold corresponding to a rollover of the distribution, in some case studies a deviation from power-law scaling is observed, indicating a process of larger block comminution due to fragmentation. Third, we found that rockfalls with fragmentation experience reduced mobility, indicated by higher reach angles, and higher lateral dispersion showing a wider distribution of trajectories. We interpret these findings as being directly related to the energy-consuming nature of fragmentation, which prevents farther deposition of fragmented rock blocks.

Albacete province, Lombardy and Aosta Valley, Yose↗

USGS2026

The potential impact of three-dimensional distributed slip models derived from real-time GNSS data on the performance of the ShakeAlert earthquake early warning system for slab interface earthquakes

The ShakeAlert® earthquake early warning (EEW) system is designed to warn users of imminent strong ground motion with sufficient time to take protective actions. ShakeAlert currently uses three algorithms to characterize the earthquake source. One estimates the location and magnitude using the first few seconds of the P wave and, while fast, tends to underestimate magnitude for M w 7.0+ earthquakes. A second estimates the location, orientation, length, and corresponding magnitude of a line source using observed peak ground acceleration and contributes primarily to M w 5.5+ earthquakes. The third infers earthquake magnitude from peak ground displacement measured using Global Navigation Satellite System (GNSS) data and offers nonsaturating magnitudes for M w 7.0+ earthquakes. Other EEW algorithms exist that infer temporally evolving spatially variable slip on a 3D fault surface from real‐time GNSS data, information that might enable more accurate and timely alerts in the event of large‐magnitude subduction interface earthquakes. Here, we evaluate the potential contribution of one such algorithm, BEFORES ( Minson et al. , 2014 ), to improve ShakeAlert performance through a simulated real‐time implementation of Bayesian evidence‐based fault orientation and real‐time earthquake slip (BEFORES) and other ShakeAlert algorithms using data for eight M w 7.6+ earthquakes. The test results demonstrate that BEFORES can produce well‐constrained and accurate magnitude estimates as soon as or sooner than other EEW algorithms, in turn enabling it to increase the amount of warning time users receive in many cases. However, with a modified Mercalli intensity (MMI) threshold of 3.5, which is commonly used for issuing alerts, BEFORES would tend to alert large geographic regions that did not feel strong shaking (MMI 6+). This effect can be mitigated using a higher alert threshold of MMI 4.5 without negative impact on the amount of warning time obtainable with BEFORES.

Bulletin of the Seismological Society of America↗

USGS2026

Groundwater budget for the Mountain Home area, southern Idaho, 2022–23

The U.S. Geological Survey, with funding from the Idaho Department of Water Resources, developed a groundwater budget for the Mountain Home area in southern Idaho for irrigation year 2023 (November 1, 2022–October 31, 2023). This study focused on the water balance across the Cinder Cone Butte Critical Groundwater Area (CGWA), Mountain Home Groundwater Management Area (GWMA), and the rest of the study area (RoSA), compiling data from various sources, including precipitation records, groundwater level measurements, metered groundwater pumpage data, surface water diversions and evapotranspiration (ET) estimates derived from remote sensing satellite imagery, and ground-based reference data. Key inflow components included recharge from applied surface water irrigation (which incorporates incidental recharge from irrigation practices and conveyance losses), estimated tributary streamflow, and estimated mountain block recharge. The key outflow components were groundwater pumpage for irrigation, municipal, industrial, and domestic uses, and ET. Recharge from applied irrigation and mountain block recharge were the largest inflows, and groundwater pumpage for irrigation was the largest outflow. The CGWA had a positive groundwater budget residual of 2,170 acre-feet (acre-ft), which contrasts with observed long-term groundwater level declines and historical trends of storage depletion. This positive residual is likely associated with unquantified outflows, including lateral groundwater flow out of the subregion, or other complexities, such as overestimated tributary contributions relative to the actual recharge for the 2023 water budget. The GWMA exhibited a positive residual of 56,563 acre-ft, primarily owing to recharge from applied surface water irrigation and areal recharge during a wetter-than-average year, which allowed irrigation entities to deliver more water from in-basin and out-of-basin reservoirs. The RoSA showed a large positive residual of 124,933 acre-ft. The interpretation of these positive residuals must account for significant uncertainties, including estimations of areal recharge, tributary streamflow (particularly losses and diversions), ET, the volume of surface water loss to the Snake River, lateral groundwater flows between subregions and across study area boundaries, and the unquantified groundwater discharge to the Snake River. These uncertainties, in combination with the complex hydrogeologic controls on water movement and limitations of remotely sensed data, directly affect the accuracy of water availability assessments. Future data collection efforts would help reduce these uncertainties and support water resource management decisions in the Mountain Home area. Key efforts could include installing additional streamflow gaging stations (particularly to quantify tributary losses and gains and surface water losses to the Snake River), improving groundwater pumpage metering, and validating remotely sensed ET data with ground-based measurements. Furthermore, to better quantify unrepresented or highly uncertain fluxes, focused investigations on groundwater discharge to the Snake River, lateral groundwater flows between subregions and across study area boundaries, and a more robust determination of the actual influence and volume of mountain block recharge would help refine future water availability assessments for the Mountain Home area.

Idaho↗

USGS2025

Mapping a Carrington storm

A map is presented of median 1-min-resolution peak geoelectric-field strength across the United States as would be induced by magnetic storms as intense as the 2 September 1859 Carrington storm. The map is constructed from two data sets: Magnetometer time series from 22 ground-based observatories recording 40 magnetic storms, and surface impedance tensors derived from magnetotelluric measurements acquired at 1616 survey sites across the contiguous United States. Carrington-class storm geoelectric fields are likely to be very strong in the United States East and Midwest; > 5.00 V/km at many places. In Virginia, strengths would likely range from 30.30 V/km, with a 68% confidence interval of [19.44,47.20] V/km, to as low as 0.05 [0.03,0.07] V/km. Comparison of model geopotentials with those measured on 30 long lines, indicates errors of about 18%. A Carrington-class storm would likely induce geoelectric fields with strengths 55% greater than for the 13–14 March 1989 storm.

contiguous United States↗

USGS2022

Waters divided: A history of alluvial fan research and a view of its future

Flows exiting confined valleys tend to deposit sediment in fan-shaped landforms. Where deposition is wholly or largely by the tractive forces of flowing water, these landforms are called alluvial fans. They are the product of the progressive division of water and sediment downfan, from slopes that may exceed 0.10 to distal slopes that may be below 0.01. Channel depths also tend to decline, from values that approach one to several meters at steep fanheads, to a few decimeters at distal fan margins. The result is a radiating, depositional ramp where confined or unconfined flows transport sediment from source basins to bounding streams, subsiding basins, or stable platforms. Where streams or subsiding basins consume the sediment supply from the source basin, the fan may approach a steady form whose extent and distal slope are set by stream location or subsidence rate. Where boundary conditions do not remove sediment, the fan may prograde out to long distances and low slopes (<0.01). Theoretical and experimental work over the past several decades support the notion that alluvial fan long-profiles become steeper as sediment supply increases or transport capacity decreases, and increasingly concave upward as the rate of bed material deposition decreases downfan. Grainsize distributions of alluvial fans seem to span the range observed in alluvial rivers, with no processes that uniquely identify them, apart from the distributary pattern of deposition. Bed sand cover tends to increase downfan in arid-region fans, with an absence of systematic downfan fining of coarser grain sizes. Surficial mapping and geochronology have demonstrated that fan deposition varies greatly through time, arguably from climate variations that alter hillslope sediment supply. The combination of surficial mapping and hydraulic modeling with high-resolution topography can now produce detailed flood susceptibility maps. The effective use of these maps to protect lives and property, however, depends on answering many of the enduring questions about the mechanics of how water and sediment divide down alluvial fans.

Book chapter↗

USGS2026

Archean to Mesozoic–Cenozoic seismic crustal structure: Implications for geological and biological evolution

We use >4500 measurements of crustal structure to investigate the seismic structure of continental crust, Archean to Mesozoic–Cenozoic. The mean crustal thickness of continents, including their margins, is 36.5 km. We find that Archean, Paleoproterozoic, and Mesoproterozoic crust have similar mean crustal thickness (41 km), seismic velocities, Vp/Vs ratio, and density. Crusts of these ages span 3 Ga, from 4.0 to 1.0 Ga, and their similarity of physical properties (thickness, density, and seismic velocities) suggests that the process of crustal formation may have been similar during this time period, which covers 66% of Earth's history. Neoproterozoic crustal properties differ significantly from Mesoproterozoic and older crusts. The mean thickness of Neoproterozoic crust is 32 km, some 9 km thinner than the mean thickness of Archean, Paleoproterozoic, and Mesoproterozoic crust. A lithospheric root with a thickness of 150–200 km underlies Archean, Paleoproterozoic, and Mesoproterozoic crusts, and this root resists lithospheric rifting and crustal extension and thinning. Based on these observations, particularly the lithospheric thickness, we conclude that Archean, Paleoproterozoic, and Mesoproterozoic lithosphere are unique and together form the stable cratonic nuclei, defined as the thick (150–250 km), long-lived cores of continents. Higher mantle temperatures during the Archean, Paleoproterozoic, and Mesoproterozoic may have played a key role in the formation of the thick lithospheric roots. The second key finding is that the similarities of mean seismic properties indicate that the process of crustal formation operating in the Archean continued during the Paleoproterozoic and Mesoproterozoic. The thin (mean value < 135 km) lithospheric root beneath Neoproterozoic and younger crust may be related to the steady decrease in mantle temperature through time. Neoproterozoic and Paleozoic crust have similar physical properties, and these eras are characterized by pronounced biodiversification, including the renowned Garden of Ediacara, the Cambrian Explosion, and the Great Ordovician Biodiversification. Mesozoic–Cenozoic crust is the most diverse and reflects the current tectonic and magmatic processes of crustal formation.

Book chapter↗

USGS2024

Produced water geochemistry from hydraulically stimulated Niobrara Formation petroleum wells: Origin of salinity and temporal perspectives on treatment and reuse

Produced water (i.e., a mixture of returned injection fluids and geologic formation brines) represents the largest volumetric waste stream associated with petroleum production in the United States. As such, produced water has been the focus of intense study with emphasis on understanding the geologic origin of the fluids, environmental impacts of unintended or intentional release, disposal concerns, and their commodity (e.g., lithium) potential. However, produced water geochemistry from many active petroleum plays remain poorly constrained leading to knowledge gaps associated with the origin of brine salinity and parameters (e.g., radium levels) that can impact treatment, disposal, and possible reuse. Here we evaluate the major ion geochemistry, radium concentrations, and stable water isotope composition of ~120 produced water samples collected from 17 producing unconventional petroleum wells in Weld County, Colorado from the Late Cretaceous Niobrara Formation. This sample set encompasses eight produced water time series from four new wells across production days 0 to ~365 and from four established wells across production days ~1000 to ~1700. Additionally, produced water from nine other established Niobrara Formation wells were sampled at discrete time points ranging from day 458 to day 2256, as well as hydraulic fracturing input fluids. These results expand the available Niobrara Formation produced water geochemical data, previously limited to few wells sampled within the first year of production, allowing for the heterogeneity of major ions and radium to be evaluated. Furthermore, we explore the geochemical relationships between major ion ratios and stable water isotope composition to understand the origin of salinity in Niobrara Formation brines from the Denver-Julesburg Basin. These findings are discussed with perspective toward potential treatment and reuse of Niobrara produced water prior to disposal.

Colorado, Nebraska, Wyoming↗