place

United STates

United STates: explore 1092 source-linked works published from 1918 to 2026, with original documents and citations.

This collection is a preview while coverage and quality are evaluated.

Coverage and selection

Includes records with this source-supplied label or an explicit phrase match in their metadata. Matches indicate a mention, not proof that a paper uses a method or tests a material. Source versions are consolidated by DOI.

Sources: usgs. Collection updated 2026-09-14. Counts describe this index, not the complete source archives.

USGS2026

Insights into Mountain Pass carbonatite formation from in-situ sulfur isotopes and geochemistry of sulfate and sulfide minerals

The Mountain Pass carbonatite stock hosts a world-class rare earth element deposit and may be classified as a carbonate-sulfate igneous rock, as it contains on average > 50 volume percent carbonate minerals and 20 to 30 volume percent sulfate minerals. The sulfates range in composition from barite to celestine and locally occur with sparse sulfide minerals. We investigate the origin of sulfur enrichment and the occurrence of sulfur-bearing minerals in the Mountain Pass carbonatite with in-situ sulfur isotope and mineral chemistry. Barite cores with δ 34 S of 1 to 3‰ do not coexist with sulfides, whereas celestine rims with δ 34 S of > 3‰ are associated with sulfides with δ 34 S < -10‰. We propose a model in which sulfur-bearing sediments were subducted during episodes of plate convergence in the Mojave Province that preceded Mountain Pass magmatism. Metasomatism of the overlying mantle by melts derived from the subducted sediments generated an unusually carbon- and sulfur-rich source to yield carbonatite magmas. Sulfur from primary carbonatite magmas and ~ 1 to 7% sulfur from subducted sediment melts yielded a slightly enriched δ 34 S composition (relative to depleted mantle δ 34 S of -1‰) for early crystallizing barite. Celestine rims on magmatic barite cores formed at low, hydrothermal temperatures (< 350 °C) based on S isotope thermometry for equilibrium celestine-galena and celestine-pyrite pairs. The sparse sulfides in the carbonatite stock are not in equilibrium with the primary barite cores and therefore do not permit S isotope thermometry estimates of magmatic temperatures. The S/Se ratios of sulfide minerals (> 3,400) typically exceed primitive mantle values (S/Se of 3,340), also consistent with their derivation from hydrothermal fluids. Trace occurrences of sulfide and sulfate minerals in alkaline silicate stocks related to the carbonatite stock have similar δ 34 S compositions and yield similarly low formation temperatures, suggesting regionally extensive and chemically similar sulfur-bearing hydrothermal fluids that imparted lithologically diverse rocks with a consistent sulfur isotope fingerprint.

California, Nevada↗

USGS2026

Deaf, deafblind, and hard of hearing university student experiences with earthquake early warning in the United States: Evaluating language planning and technology access

The growing literature on deaf and hard of hearing (DHH+) populations and disasters demonstrates that emergency communication (including alerts) is not reaching global DHH + individuals with dangerous impacts for morbidity and mortality. This is the first research study in the U.S. to qualitatively explore the experiences of DHH + persons with earthquake early warning (EEW) through group-based dialogue sessions. The study investigates eight DHH + university students'past earthquake experiences, access to EEW alerts, and perceptions of ShakeAlert Ⓡ , an EEW system for detecting earthquakes and alerting residents of California, Oregon, and Washington. Findings highlight key gaps in disaster alert usability within four thematic areas: lack of messaging in participants' language(s), unclear alert messaging, deficient message delivery mechanisms for deafblind persons, and insufficient access to earthquake information and training that leads to dependence on informal information networks. Weaknesses identified in these four themes reduce DHH + trust in EEW systems and compromise the capacity of alert recipients to take swift protective action or to mentally prepare before shaking starts. The study also underscores structural factors such as insufficient linguistic representation in disaster language planning and technology design, which ignores the linguistic and sensory access needs of DHH + individuals. Building on disaster language planning frameworks, we recommend involving DHH + populations to co-develop EEW alerts. By centering DHH + perspectives, this research contributes to ongoing efforts to ensure that EEW systems reach everyone.

California↗

USGS2026

Constraining source and path effects of large magnitude earthquakes using ground motion simulations

The purpose of this study is to use ground‐motion simulations to investigate ways in which source and path effects for large‐magnitude earthquakes can be represented in nonergodic ground‐motion models (GMMs). To achieve this, we designed a ground‐motion study in the San Francisco Bay Area that includes earthquakes with a broad range of magnitudes distributed uniformly on a fault plane, and sites covering a large range of rupture distances and azimuths. After running a large suite of kinematic simulations (magnitude 4–7), we then develop a nonergodic GMM with the simulated data. We find that trends in the within‐site residuals are affected significantly by the earthquake radiation pattern, rupture directivity, and slip patterns. Next, we modify an existing rupture directivity model to fit and remove the observed radiation pattern and rupture directivity from the residuals. We also minimize the contributions of slip patterns by averaging the within‐site residuals among multiple source realizations. Finally, after removing the source effects from the within‐site residuals, we compare the path effects computed with different magnitude groups using two approaches. The first approach only considers the small events that have the same shortest path to a site as the large events, whereas the second approach considers all small events on the fault plane. The results indicate that it is difficult to satisfactorily approximate the path effects of large events with those of small events using either approach, at least in the case of simulations.

California↗

USGS2026

Bottom-up characterization of geologic methane emissions in the San Juan Basin in the southwestern USA

Methane is a potent greenhouse gas that plays an important role in atmospheric chemistry and global warming. The current global methane budget has large uncertainties, and a better understanding of the budget would help to guide strategies for reducing anthropogenic emissions to fight climate change. Natural geologic methane emissions are a particularly poorly constrained source, with top-down estimates from 14 C in ice cores suggesting much lower geologic emissions than bottom-up scaling of direct flux measurements. Our study aims to contribute to resolving this discrepancy through improved bottom-up characterization of geologic methane seepage in the San Juan Basin in southwestern Colorado and northwestern New Mexico, USA. We performed 983 new flux chamber measurements in this basin during summer 2022 and winter 2023 field campaigns. Our results, in combination with prior measurements, suggest that natural seepage in the San Juan Basin only occurs on or near the Fruitland coal outcrop. Specifically, our new measurements confirm previous measurements of seepage along the northwestern exposure of the Fruitland outcrop in Colorado (a known hydrodynamic overpressure region) and for the first time, identified seepage locations along the southernmost Fruitland outcrop exposure in New Mexico, in association with a coal cleat and a fault. Overall, seepage along the Fruitland coal outcrop is heterogeneously distributed, with both positive and negative (interpreted as microbial soil sink) methane fluxes. Features that are hypothesized to be predictive of seepage (e.g., faults) were not associated with positive methane fluxes in areas outside of the Fruitland outcrop. Our best estimate for total geologic methane seepage in the San Juan Basin from spatial interpolation and statistical upscaling is approximately 0.14 Tg CH 4 /yr, with a range from 0.029 to 0.48 Tg CH 4 /yr. This best-estimate value is lower than a previous bottom-up estimate from a gridded seepage inventory, but higher than a previous top-down estimate.

Colorado, New Mexico↗

USGS2026

Multi-objective optimization of a hydro-economic model in an over-allocated agricultural basin

Groundwater depletion for agricultural irrigation poses significant environmental and economic challenges. This study introduces a proof-of-concept that combines hydro-economic modeling, scenario-based modeling, and multi-objective optimization to manage pumping curtailment in an over-allocated basin in the western United States. Three optimization scenarios were evaluated, each offering different degrees of management flexibility. Results reveal that scenarios with finer spatial resolution achieved greater environmental benefits per unit profit loss. Additionally, strategies allowing fractional reductions in curtailed wells–rather than complete shutdowns based on water rights seniority–substantially improved efficiency, highlighting the value of increased decision-making flexibility. Although scenario testing can aid stakeholder engagement and strategy exploration, multi-objective optimization provides a systematic framework to quantify tradeoffs between competing objectives. This combined approach demonstrates promise for building consensus and supporting the design of sustainable water management strategies that balance agricultural livelihoods with ecosystem preservation.

Oregon↗

USGS2026

New constraints on northeast Seattle basin structure from converted seismic waves

The Seattle basin is a deep sedimentary basin in the Seattle–Bellevue, Washington metropolitan area within the Puget Lowland of Washington State. We determine the structure of a portion of the basin and the underlying basement using analysis of P waves converted from direct S incident from below. A deep local crustal event beneath Monroe, about 35 km northeast of Seattle, was recorded by a 100‐station nodal array deployed in 2019. The event produced a variety of coherent seismic phases, including converted waves from the sediment—basement boundary, internal structure within the basin, and additional crustal discontinuities. Using observed Sp converted waves, we apply an adjoint‐based full waveform inversion (FWI) method to determine the amplitude and extent of seismic discontinuities at depth. We find the strongest source of converted waves for this event lies ∼6 to 7 km depth below northern Lake Washington, interpreted to be the local depth to basement rock. The newly imaged shallow basement structure may be part of a deformation zone associated with the Siletzia eastern boundary. Our results highlight the utility of converted seismic waves recorded by a dense array, combined with an FWI method, to illuminate crustal structure.

Washington↗

USGS2026

A three-dimensional geologic framework model of the northern Great Plains region of Montana, North Dakota, South Dakota, and Wyoming, USA

This report presents a new three-dimensional geologic framework model (GFM) of the northern Great Plains region, encompassing parts of Montana, North Dakota, South Dakota, and Wyoming. The model provides a regionally consistent, geographic information system (GIS)-ready representation of Phanerozoic sedimentary strata, major fault systems, and Precambrian basement geometry across two sedimentary basins and adjacent uplifts. More than 300,000 geologic and geophysical data inputs were synthesized to model 41 stratigraphic horizons and 47 faults, yielding an internally coherent, sealed-volume interpretation of the subsurface. The modeling workflow developed for this study demonstrates an efficient and scalable approach for constructing basin-to regional-scale GFMs in geologically complex and data-variable settings. Although model fidelity varies with data density and quality, the resulting geometry is broadly consistent with 1:500,000-scale geologic mapping and highlights areas where additional geologic study is most needed. The three-dimensional GFM provides a foundational framework to support groundwater, energy, and mineral resource assessments, and offers a transferable methodology for potential future U.S. Geological Survey efforts to build large-area subsurface models in underexplored regions of the United States.

Montana, North Dakota, South Dakota, Wyoming↗

USGS2026

The U.S. Geological Survey 2025 Puerto Rico and U.S. Virgin Islands time-independent earthquake rupture forecast

We present the 2025 U.S. Geological Survey Puerto Rico and U.S. Virgin Islands (PRVI) time‐independent earthquake rupture forecast (ERF), developed for the 2025 update to the National Seismic Hazard Model (NSHM) for PRVI. The updated ERF improves upon a prior model from 2003, including an expanded fault inventory with slip‐rate estimates, updated seismicity catalogs, and refined subduction zone geometries and deformation models. It applies the fault‐system inversion methodology to solve for rates of ruptures on modeled faults, adapted from the 2023 NSHM (NSHM23) for the western United States, including the first application of the inversion to model rates on a U.S. subduction interface. Off‐fault and intraslab seismicity are constrained by observed seismicity and use updated methods developed for NSHM23. Uncertainties in model components are substantial, and the ERF represents epistemic uncertainties through a comprehensive logic tree consisting of 1.7 billion logic‐tree branches combined across all sources.

Puerto Rico, U.S. Virgin Islands↗

USGS2026

Stress states on the eve of past earthquakes inform earthquake rupture through fault complexity along the San Andreas and San Jacinto faults

Estimating the evolving state of stress along active fault systems can provide insight into the conditions that generated past ground‐rupturing earthquakes and influenced their ability to propagate through areas of geometric complexity, such as fault branches and stepovers. We use quasi‐static forward numerical models that incorporate the 3D complex configuration of active faults in southern California to estimate shear tractions on the geometrically complex southern San Andreas and San Jacinto faults from 1000 to 1900 C.E. These tractions include interseismic accumulation of traction due to tectonic loading, viscoelastic relaxation of shear stress within the upper crust between earthquakes, and effects of other earthquakes on the fault network. We simulate ground‐rupturing earthquakes based on the along‐strike earthquake extents modeled by Scharer and Yule (2020) , assuming that stress drop is complete in each earthquake. We use Monte Carlo simulations to estimate uncertainty in evolving shear tractions due to uncertainties in earthquake timing and in upper‐crustal viscosity. Pre‐earthquake shear tractions typically do not exceed ∼2 MPa. Although ruptures with length <200 km have pre‐earthquake shear tractions that range from near zero to ∼1.75 MPa, these tractions are not less than ∼0.4 MPa for earthquakes with rupture length >200 km. Earthquakes with long (>200 km) ruptures occur only in the single‐stranded part of the system, whereas those with short (<125 km) rupture length and high pre‐earthquake shear traction occur near fault stepovers and branches. This suggests that high accumulated shear traction encourages longer rupture propagation, but may not be sufficient to overcome geometric complexities. This modeling approach informs our understanding of rupture propagation and provides estimates of fault shear tractions that are unavailable from direct measurements.

Callifornia↗

USGS2026

Multi-site evaluation of a postfire debris-flow runout forecast method

Postfire debris flows pose a hazard to human life, property, and infrastructure when they travel from steep source areas to urbanized alluvial fans or other developed areas. Existing methods for rapid (<1 week) postfire debris-flow hazard assessment document the increase in the likelihood and size of debris flows as the magnitude of high-intensity rain necessary to initiate debris flows increases but do not indicate the extent of downstream debris-flow runout. Although many models for the simulation of debris-flow motion are available, there is no established approach for using these models to delineate locations susceptible to postfire debris-flow runout that (a) is feasible to use at the spatial scale of an entire fire; (b) is appropriate for runout onto unconfined areas; (c) reproduces observed relations between runout and rainfall intensity; and (d) characterizes inherent uncertainty in runout, even without spatiotemporally variable rainfall. We propose and evaluate a method for generating postfire debris-flow runout hazard maps that has all the above qualities. Selection of case studies prioritized events triggered by a range of rainfall intensities, locations within and outside of southern California, and observed runout onto unconfined topography. Qualitative and quantitative assessment of performance for four events indicate that simulation results broadly match observations albeit with some discrepancies at a scale larger than structure or land parcel level (approximately 20-m by 20-m). The method may be used to identify potentially hazardous areas immediately following a fire and to provide approximate runout forecasts when a storm is imminent.

Arizona, California↗

USGS2026

Potential impacts of groundwater pumping on stream temperature are greatest in streams with substantial cold groundwater inflows

Groundwater pumping-induced reductions in streamflow (known as ‘streamflow depletion’) have been documented worldwide, but potential impacts of streamflow depletion on stream temperature are not well understood. Here, we use two types of models to identify potential impacts of pumping on stream temperature across the conterminous United States (CONUS) to determine which aspects of a stream's annual thermograph (thermal signatures) can be used to monitor and manage streamflow depletion impacts on stream temperature. We used long-term streamflow and stream temperature data from 30 streamgages across CONUS and surrogate models of streamflow depletion to analyse potential stream temperature impacts at each site. We compared two different stream temperature modelling approaches: (i) a process-based energy balance model and (ii) statistical regression models based on air temperature and stream discharge. We calculated a suite of thermal signatures under depleted and non-depleted conditions for each stream and found that maximum annual 7-day temperature and annual temperature range are potentially the most sensitive to streamflow depletion, with potential changes of at least 2°C at > 70% of the sites when using the process-based model. We also found that the regression-based models predicted much less sensitivity of stream temperature to streamflow depletion than the process-based model. This work provides an initial evaluation and sensitivity analysis of the potential impacts of streamflow depletion on stream temperature. We demonstrate that stream temperature may be most sensitive to pumping in streams with a high proportion of flow sourced from relatively cold groundwater inputs, and that regression-based stream temperature models may underpredict stream temperature changes caused by streamflow depletion.

conterminous United States↗

USGS2026

Orientation dependence of probabilistic seismic hazard estimates from CyberShake physics-based simulations

Earthquake ground‐motion intensities, such as pseudospectral accelerations (SAs), can vary significantly with horizontal orientation. However, conducting probabilistic seismic hazard analysis (PSHA) for each horizontal orientation is challenging because current ground‐motion models used in PSHA consider only a single horizontal intensity value, usually the median across all orientations, known as RotD50. To address this limitation, we employ physics‐based simulations for PSHA, which contain full waveforms from which ground‐motion intensities can be computed for all horizontal orientations to study directional seismic hazard. We apply our approach to the latest CyberShake study of the Greater Los Angeles metropolitan area, developed by the Statewide California Earthquake Center, finding that seismic hazard at a 2475‐yr return period, a common value used for earthquake‐resistant design, varies significantly with horizontal orientation. For instance, for SAs at 3 s, the maximum seismic hazard across all horizontal orientations is, on average, 15% higher than the median RotD50 hazard, with these differences becoming more pronounced at longer periods. These observed variations can generally be attributed to physical mechanisms that polarize seismic waves, such as the radiation pattern of the earthquake source and the influence of the subsurface structure. These results may have important implications for earthquake engineering applications, particularly for long‐period structures in areas with substantial horizontal variations in seismic hazard.

California↗

USGS2026

Evidence for organic and igneous carbon sources in the Mesoproterozoic Adirondack Highlands graphite district, New York

This study examined the origin of graphite mineralization in the Mesoproterozoic Adirondack Highlands. Flake graphite was mined from the historic mining district around Ticonderoga, New York, USA, from the 1850s to the 1920s. Most of the production and reserves from the district are disseminated graphite from feldspathic quartz schists (e.g., Dixon-American Graphite Co. mine, Faxon property, Flake Graphite Company). Protolith sediments of the schists were deposited in the Trans-Adirondack backarc basin at ca. 1.25 Ga and were metamorphosed to high grade at 1.18–1.15 Ga during the accretionary Shawinigan orogeny and at 1.09–1.02 Ga during the Ottawan phase of the Grenvillian orogeny. Major elements and carbon isotopes of graphite from schists (δ 13 C = −28.1‰ to −5.4‰ relative to Vienna Peedee belemnite) are consistent with derivation from organic carbon in protolith sandstones that devolatilized during metamorphism. Some small mines and prospects in the district mined graphitic calcite marbles and clinopyroxene ± scapolite skarns associated with a pegmatite suite dated by U-Pb in zircon to 1.04–1.01 Ga (e.g., Lead Hill and Crown Point deposits), and rare graphite veins are observed at several deposits. Hydrothermal flake graphite from pegmatites, skarns, and veins are distinct from graphite in schist-hosted deposits and have δ 13 C values that cluster around −7‰ to −5‰, consistent with deposition from igneous CO 2 . Covariation of carbon and oxygen stable isotopes and U-Pb zircon geochronology are not compatible with precipitation from fluids derived from country rock metapelites or marbles. Some graphite districts elsewhere show compelling evidence for remobilization and redeposition of sedimentary-derived carbon as hydrothermal graphite, sometimes as thick veins producing world-class, high-grade deposits. In the Adirondacks, there is evidence for two generations of graphite mineralization with distinct carbon sources: syngenetic graphite with an organic origin in Shawinigan quartz schists and epigenetic graphite with an igneous origin in Ottawan skarns and pegmatites.

New York↗

USGS2026

Small-volume tephra deposits of the May 1924 explosions from Halemaʻumaʻu, Kīlauea volcano, and their origin

More than 50 explosive eruptions occurred from Halemaʻumaʻu at Kīlauea volcano over 17 days from May 11 to 27, 1924. Ballistics weighing as much as 14,000 kg were ejected and most landed within 2 km of the vent. Fine ash made up a major component of the tephra and was dispersed tens of kilometers downwind. Draining of the Halemaʻumaʻu lava lake occurred in late February 1924, with the crater floor eventually subsiding by a further ∼70 m (to ∼180 m below the crater rim) by the time the first explosions took place during the night of May 10–11. The largest explosions occurred on May 17–18 and smaller explosions continued until May 27, at which point Halemaʻumaʻu had more than doubled in width and depth. The explosions generated plumes reaching up to ∼10 km high with ballistics ejected up to 2 km from the crater. Almost 100 years later, we investigate and characterize the preserved tephra deposits within ∼3 km of the 1924 crater rim. Grain size and shape analyses were performed on 202 samples collected from 34 tephra profiles using dynamic image analysis, with a subset of layers from nine tephra profiles used for componentry (200 grains per layer in the 0.5–1 mm size fraction). Additionally, we characterize the average diameters (using the five largest clasts) at 216 locations and measure the average diameters of 2291 ballistics (largest per ∼100 m 2 area). Physical descriptions from fieldwork and grain size distributions were used to subdivide the tephra layers into five lithofacies: coarse homogeneous, fine homogenous, red ash, accretionary lapilli-bearing, and finely laminated. Grain size versus shape data show a range of values that demonstrate most grains are dense, smooth, and equant, in alignment with lithic clasts dominating the tephra componentry. The fine grained and accretionary lapilli-bearing nature of some of these lithofacies confirms that water influenced the style of the explosions. However, we also note juvenile clasts within many of the tephra layers, indicating that many of the layers were formed during phreatomagmatic explosions (sensu stricto), despite the eruptive mechanism being dominantly phreatic. Juvenile clasts are more abundant higher in the tephra profiles, suggesting that juvenile magma was more involved later in the explosive sequence. Thermal and hydrologic modeling indicate that groundwater inflow into a short-lived, small-diameter volcanic conduit (10-m to 120-m-diameter used for modeling) during the 78–85 days preceding the first explosion provides a physically plausible mechanism for this eruptive sequence.

Hawaii↗

USGS2026

Groundwater drought in the United States: Spatial and temporal variability

Many communities and ecosystems in the United States that are dependent on groundwater are potentially adversely affected by groundwater drought. We computed yearly groundwater-drought metrics and mean groundwater levels at well locations across the conterminous United States (CONUS), using data from wells and remotely sensed and modeled Gravity Recovery and Climate Experiment Drought Monitor Data Assimilation (GRACE-DADM). We also modeled the probability of low or high human impact at each well location. The spatial distribution of groundwater-drought duration and severity from 2001 to 2020 for 1,510 wells shows longer maximum duration and higher maximum severity events in drier regions like the Southwest than in wetter regions like the Northeast. Based on 613 wells in CONUS from 1981 to 2020, there are many significant decreases in drought duration and severity in the Northeast and many significant increases in annual-mean groundwater levels. In contrast, there are many significant increases in drought metrics and decreases in mean water levels in parts of the Southeast. There are major differences in trends from 2001 to 2020 between well-based and GRACE-DADM-based groundwater metrics in some CONUS regions and a very low correlation between trends at individual locations across CONUS. A potential reason for this disparity is the low GRACE-DADM resolution (∼12 km) and the potential for a large amount of groundwater variation at the local scale. Also, GRACE-DADM represents shallow, unconfined aquifers which may not match the screened interval of the monitoring wells we evaluated. Large spatial gaps in long-term, high frequency, and quality-assured groundwater-well monitoring data present a challenge for understanding groundwater-drought variability across CONUS. Remote sensing tools such as GRACE can help but cannot fully replace well monitoring, as highlighted by our study results. Substantially more long-term monitoring wells would more accurately represent groundwater-drought trends and spatial variability across CONUS, particularly in western regions.

conterminous United States↗

USGS2026

Seasonal and hydrologic variation influences habitat and functional structure of stream fish assemblages

Introduction: Hydrologic variability is a key driver of ecological structure in lotic systems, shaping habitat conditions, taxonomic diversity, and the functional traits that mediate species’ persistence and performance (e.g., reproductive success). While many studies examine taxonomic responses to variation in flows, few evaluate how spatiotemporal hydrologic variation influences the functional organization within stream fish communities. Methods: We quantified seasonal habitat structure and functional trait diversity of fish assemblages across six Ozark Plateau headwater streams representing two contrasting flow regimes: Groundwater Flashy and Runoff/Intermittent Flashy. Fish and habitat data were collected seasonally during a dry year (2002) and a wet year (2003). Functional space was constructed using PCoA of morphological, ecological, and life-history traits, and functional diversity was measured using community weighted means (CWMs), functional richness (FRic), functional evenness (FEve), and functional divergence (FDiv). Results: We found that habitat structure differed strongly by flow regime and season, with Runoff/Intermittent streams exhibiting pronounced reductions in depth, area, and velocity, while groundwater streams remained structurally stable. Functional identity of assemblages was similar across flow regimes, dominated by benthic, hydrodynamic taxa with opportunistic and periodic life-history strategies. However, functional structure differed significantly: FEve and FDiv were consistently lower in Runoff/Intermittent Flashy streams in both years, indicating assemblage dominance of species with similar trait combinations and reduced trait partitioning under variable flow. FRic and taxonomic richness remained stable across seasons and flow regimes, suggesting high functional redundancy despite species turnover. Discussion: Together, results show that flow regime mediates both habitat structural stability and functional organization. As climatic warming and extreme drought increase hydrologic instability in headwaters, functional trait approaches provide a sensitive tool for detecting losses of functional roles that may not be evident by using taxonomic metrics alone.

Arkansas, Kansas, Missouri, Oklahoma↗

USGS2026

Transit time modeling framework for predicting freshwater salinization in urban catchments

The salinity of inland freshwaters is rising globally, particularly in urban watersheds where winter road deicers are widely applied. Attributing stream salinity dynamics to specific sources and transport pathways remains challenging due to episodic salt inputs, engineered drainage, and strong coupling between hydrology and subsurface storage. We present a modeling framework that couples climate-driven deicer build-up and wash-off with transient transit time distribution theory to simulate salt transport through drainage, interflow, and groundwater pathways. Applied to an urban watershed in Northern Virginia (USA), the model reproduces ten years of high-frequency stream salinity measurements across daily-to-decadal timescales. The calibrated model implies an average deicer application of 206 tonnes Cl −1 yr −1 , or roughly one 20 kg bag of rock salt person −1 yr −1 when normalized by the 20,000 people living in the watershed. In winter months, higher infiltration routes a large fraction of snowmelt and deicers into shallow subsurface pathways, enhancing vadose-zone and interflow contributions to stream salinity. Limited subsurface storage capacity and seasonal hydrologic turnover flush excess chloride from the vadose zone and groundwater during subsequent summer storms. By linking climate-driven deicer inputs, hydrologic connectivity, and stream water age, the framework provides a transferable basis for diagnosing and managing freshwater salinization in urban watersheds.

Virginia↗

USGS2026

Water use permits as of July 2024 and reported water use near the North Unit of Theodore Roosevelt National Park, North Dakota, 1980–2023

Starting in the early 2000s, increasing oil and gas development in western North Dakota created a need for additional water resources from surface-water and groundwater sources near the North Unit of Theodore Roosevelt National Park. To summarize the use of water in that area, the U.S. Geological Survey, in cooperation with the National Park Service, developed a map of surface-water and groundwater resources, aquifers, and water-use diversions, and plotted water-use trends from 1980 to 2023. Reported water used from permits in the map area has more than doubled since 2020, increasing from about 750 acre-feet in 2020 to about 2,300 acre-feet in 2022 and 2,000 acre-feet in 2023. Surface water provided the primary source of reported water used for the study period with an average of about 410 acre-feet per year from 1980 through 2017 and about 1,330 acre-feet per year from 2018 through 2023. After 2011, groundwater sourced from the Little Missouri River, Tobacco Garden Creek, Fox Hills, Fort Union, and Dakota aquifers became a larger portion of total annual reported water use from permits in the map area. From 1980 through 2015, water use for irrigation averaged 86 percent of the total annual reported surface-water and groundwater use in the map area. Starting in 2011, however, industrial uses became a proportionally larger total use of water, and in 2015, became the highest reported volume of water use in the map area. From 2011 to 2023, industrial use designated for water depots increased from 50 acre-feet to about 1,370 acre-feet, accounting for about 70 percent of total reported water use in the map area in 2023.

North Dakota↗

Compare source metadata on this page

Work	Published	Source identifier	Source
Insights into Mountain Pass carbonatite formation from in-situ sulfur isotopes and geochemistry of sulfate and sulfide minerals	2026-03-21	70274287	usgs
Deaf, deafblind, and hard of hearing university student experiences with earthquake early warning in the United States: Evaluating language planning and technology access	2026-03-20	70276811	usgs
Constraining source and path effects of large magnitude earthquakes using ground motion simulations	2026-03-19	70274581	usgs
Bottom-up characterization of geologic methane emissions in the San Juan Basin in the southwestern USA	2026-03-19	70274659	usgs
Multi-objective optimization of a hydro-economic model in an over-allocated agricultural basin	2026-03-19	70276384	usgs
New constraints on northeast Seattle basin structure from converted seismic waves	2026-03-18	70276913	usgs
A three-dimensional geologic framework model of the northern Great Plains region of Montana, North Dakota, South Dakota, and Wyoming, USA	2026-03-18	sir20265127	usgs
The U.S. Geological Survey 2025 Puerto Rico and U.S. Virgin Islands time-independent earthquake rupture forecast	2026-03-16	70274265	usgs
Stress states on the eve of past earthquakes inform earthquake rupture through fault complexity along the San Andreas and San Jacinto faults	2026-03-16	70275742	usgs
Multi-site evaluation of a postfire debris-flow runout forecast method	2026-03-12	70274552	usgs
Potential impacts of groundwater pumping on stream temperature are greatest in streams with substantial cold groundwater inflows	2026-03-12	70275104	usgs
Orientation dependence of probabilistic seismic hazard estimates from CyberShake physics-based simulations	2026-03-12	70276903	usgs
Evidence for organic and igneous carbon sources in the Mesoproterozoic Adirondack Highlands graphite district, New York	2026-03-12	70280076	usgs
Small-volume tephra deposits of the May 1924 explosions from Halemaʻumaʻu, Kīlauea volcano, and their origin	2026-03-11	70274211	usgs
Groundwater drought in the United States: Spatial and temporal variability	2026-03-11	70274216	usgs
Seasonal and hydrologic variation influences habitat and functional structure of stream fish assemblages	2026-03-09	70274636	usgs
Transit time modeling framework for predicting freshwater salinization in urban catchments	2026-03-09	70279328	usgs
Water use permits as of July 2024 and reported water use near the North Unit of Theodore Roosevelt National Park, North Dakota, 1980–2023	2026-03-09	sim3545	usgs

These are bibliographic comparisons, not experimental rankings. Follow the original document for methods and conditions.