USGS Science⌕ Search

SEARCH · USGS Science

Results for “Science of the Total Environment”

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

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

1,687 records · Page 78Linked to original sources
USGS2025

Niobium-rich minerals from the Sheep Creek carbonatite-related deposits, Montana, USA

The Sheep Creek area in southwest Montana contains a number of small, tabular, carbonatite-related deposits that are rich in Nb- and REE-minerals. Two of the “vein-dikes” were mined at a small scale for Nb in the late 1950s. The most abundant Nb-bearing minerals are columbite-(Fe), nioboaeschynite-(Ce), baotite, and two varieties of Ti-Nb-Fe-oxide. An early oxide with Ti/Nb (atomic) near unity later broke down to a mixture of columbite + TiO2 with distinctly lower Nb content. Zircon intergrown with the Nb-REE minerals was dated by U-Pb methods to 116±2 Ma. This age is in agreement with previous dating (Th-Pb on monazite, Re-Os on molybdenite), and indicates that the Sheep Creek deposits are early Cretaceous in age.

Montana↗
USGS2024

Subducting plate structure and megathrust morphology from deep seismic imaging linked to earthquake rupture segmentation at Cascadia

The origin of rupture segmentation along subduction zone megathrusts and linkages to the structural evolution of the subduction zone are poorly understood. Here, regional-scale seismic imaging of the Cascadia margin is used to characterize the megathrust spanning ~900 km from Vancouver Island to the California border, across the seismogenic zone to a few tens of kilometers from the coast. Discrete domains in lower plate geometry and sediment underthrusting are identified, not evident in prior regional plate models, which align with changes in lithology and structure of the upper plate and interpreted paleo-rupture patches. Strike-slip faults in the lower plate associated with oblique subduction mark boundaries between regions of distinct lower plate geometry. Their formation may be linked to changes in upper plate structure across long-lived upper plate faults. The Juan de Fuca plate is fragmenting within the seismogenic zone at Cascadia as the young plate bends beneath the heterogeneous upper plate resulting in structural domains that coincide with paleo-rupture segmentation.

Oregon, Washington↗
USGS2026

Rapid seismic and infrasound assessment of large landslides: A case study from Denali National Park and Preserve (Alaska)

Large, rapid landslides are a global hazard that can occur in remote, mountainous areas. Eyewitness reports of landslides and satellite imagery can often be limited or delayed, particularly during inclement weather. However, landslide-generated seismic and infrasound (low-frequency atmospheric sound) waves can be remotely detected in near real-time. This information can significantly expedite characterization and possible landslide response activities. Here, we highlight these capabilities using a > 4 million m³ ice–rock avalanche in Denali National Park and Preserve (Alaska). This event was detected via a landslide-specific seismic location and volume estimation algorithm deployed in Alaska, and — notably — by standard earthquake monitoring systems. Following rapid detection of this event, we combined its seismic and infrasound dataset with optical, synthetic aperture radar, and oblique aerial imagery, multitemporal digital elevation models, and a numerical flow model to reconstruct its failure timeline and dynamics. We apply array processing to infrasound signals traveling > 250 km and find that two precursory events occurred minutes prior to the main failure. We use long-period seismic signals to infer the force exerted by the landslide on the Earth and constrain the rheological parameters of our numerical flow simulation with this result and deposit morphology. The main failure produced a steeply-dipping impulsive initial downward force and reached speeds exceeding 60 m/s. This impulsive force generated relatively strong seismic body waves, which contributed to the earthquake system detection. This large, remote Alaska landslide underscores the key value of seismic and infrasound analysis for rapid landslide assessment and motivates efforts to further operationalize these approaches.

Alaska↗
USGS2024

Twentieth century extreme precipitation detected in a high-resolution, coastal lake-sediment record from California

California faces increasing economic and societal risks from extreme precipitation and flooding associated with atmospheric rivers (ARs) under projected twenty-first century climate warming. Lake sediments can retain signals of past extreme precipitation events, allowing reconstructions beyond the period of instrumental records. Here, we calibrate AR-related extreme precipitation from the last century to proxy data from lake sediments collected in the latitudinal zone of the highest frequency landfall for modern ARs in California. Excursions in erosional proxy data (Ti/Al) are positively and significantly correlated (r median = 0.45, p median = 0.04) with modern records of integrated vapor transport (IVT, kg m −1 s −1 ), a key metric of AR intensity, using correlations that incorporate age-model uncertainty. Despite the land-use change near the study site, the data suggest intense and long-lasting AR storms are identifiable in this sedimentary record. These results allow conservative inferences concerning past extreme hydrology at this site.

California↗
USGS2024

Advancing sustainable groundwater management with a hydro-economic system model: Investigations in the Harney Basin, Oregon

Groundwater resources frequently trend toward unsustainable levels because, absent effective institutions, individual water users generally act independently without considering the impacts on other users. Hydro-economic models (HEMs) of human-natural systems can play a positive role toward successful groundwater management by yielding valuable knowledge and insight. The current study explores how an HEM that captures essential physical and economic characteristics of a system can shed light on the system's processes and dynamics to benefit stakeholders, managers, and also researchers. These propositions are illustrated using the Harney Basin, Oregon, which has seen large groundwater declines in the past 20 years. The HEM shows that: (a) although current groundwater pumping rates will gradually raise costs and reduce well yields, irrigators gain the highest aggregate economic return by continuing current pumping; (b) lowland areas of the basin are hydrologically connected, which limits the efficacy of remedies focused on regulations only in some portions of the basin; (c) community expectations regarding the efficacy of several proposed solutions are overly optimistic; and (d) the study's scenarios identify interventions that would stabilize the groundwater system and prevent additional adverse impacts on residential and livestock wells and groundwater-dependent ecosystems. These interventions would require limiting groundwater pumping by nearly half and reducing annual profits by $7.5–$9.0M. The HEM also demonstrated its value to researchers: its insights shifted attention toward questions about Oregon's existing groundwater institutions and their inability to adaptively manage the transition from abundant groundwater to scarce groundwater in a timely manner.

Oregon↗
USGS2024

The 3D Elevation Program—Supporting Ohio's economy

Introduction High-quality elevation data are proving to be a resource of great economic value in dealing with many important issues in Ohio. Current and accurate high-resolution elevation data support flood risk management, water quantity and quality assessment, precision farming, conservation planning, impervious-surface modeling, forest and other natural resources management, abandoned mine and geologic hazard assessment, karst mapping, and siting of wellhead pads for horizontal drilling. These data also support coastal zone management, traffic safety and preliminary engineering site-selection studies for transportation infrastructure, solar potential and other renewable energy planning, aviation safety, and identification of features of interest or concern such as archaeological sites and orphan oil and gas wells. Critical applications that meet the State’s management needs depend on light detection and ranging (lidar) data that provide a highly detailed three-dimensional (3D) model of the Earth’s surface and aboveground features.

Ohio↗
USGS2025

U.S. Geological Survey National Groundwater Climate Response Network

What is the U.S. Geological Survey National Groundwater Climate Response Network? The U.S. Geological Survey’s (USGS) National Water Monitoring Network is a network of networks that includes the Climate Response Network (CRN). The CRN is a network of wells selected to monitor the effects of climate variability, including droughts, on groundwater systems within the United States, Puerto Rico, and the U.S. Virgin Islands. The primary purpose of the CRN is to monitor these effects on groundwater levels in unconfined or near-surface confined aquifers that are minimally affected by pumping or other human-induced stresses.

Fact Sheet↗
USGS2026

Cajon Pass and the southern San Andreas Fault System: Earthquake cycle stress accumulation and present-day loading

With over a century since the last major rupture affecting the wider Los Angeles region, tectonic stress has steadily built along the southern San Andreas and San Jacinto fault systems, raising concerns of an imminent large earthquake. Cajon Pass, located at the junction of these faults, represents a critical site for potential through-going ruptures in Southern California. We constructed new 4D earthquake cycle simulations using a 1000-year paleoseismic rupture history of the San Andreas Fault System (SAFS) to assess spatial and temporal variations in stress. A semi-analytic Fourier transform model was used to compute stress from 3D dislocations in an elastic plate overlying a Maxwell viscoelastic half-space, assuming a complete coseismic reset of resolved shear stress on ruptured elements. Results show highest stress accumulation north of Cajon Pass (∼1.8 MPa/100 years) due to greater slip rates, and lower rates south of Cajon Pass (∼1.0–1.5 MPa/100 years). By 2025, Coulomb stress is estimated at 2.8 MPa on the Mojave South (MOS) segment, 1.8 MPa on the North San Bernardino (NSB1) segment and 3.6 MPa on the San Jacinto Bernardino (SJB) segment. Segments accumulate stresses with characteristic ranges of pre-event stress interpreted as failure thresholds: 1.2–2.7 MPa for MOS, 0.4–1.6 MPa for NSB1, and 1.2–2.9 MPa for SJB. When the stress disparity between segments SJB and MOS narrows, the faults appear to rupture jointly, suggesting that stress levels may control how Cajon Pass acts as an earthquake gate. These results may inform seismic hazard assessments by linking stress evolution to fault interactions.

California↗
USGS2025

Re-Os geochronology of molybdenite by LA-ICP-MS/MS

Re-Os geochronology is a powerful tool to directly date sulphide mineralisation using phases such as molybdenite, pyrite and chalcopyrite. We present a snapshot of a method in development utilizing laser ablation inductively coupled-plasma ‘triple quad’ mass spectrometry (LA-ICP-MS/MS) for Re-Os geochronology of molybdenite. While traditional LA-(MC)-ICP-MS measurements are limited by the >90% isobaric interference correction of parent radioactive 187 Re on radiogenic daughter 187 Os, the MS/MS technology minimizes this correction through careful instrument tuning and use of a (CH 4 + H 2 + He) reaction gas mixture to mass shift 187 Os away from 187 Re. We demonstrate the potential of the MS/MS technique by presenting robust Re-Os dates for three different molybdenite samples ranging in age from ~28 to 1520 Ma.

Conference Paper↗
USGS2025

Ubiquity and causes of soil water preferential flow across 17 ecoregions

Preferential flow (PF) in soil causes the rapid transport of water, nutrients, and contaminants into the subsurface, influencing groundwater recharge and streamflow. Data scarcity has hindered the quantification of PF occurrence and the identification of its drivers across diverse ecoregions. We address this gap by analyzing high-frequency, multi-depth soil moisture data across 17 ecoregions in the USA, using ∼1,500 sensors at 40 sites. We discovered that PF is widespread, with sites experiencing PF in up to 60% of rainfall events ≥2 mm. Multiple approaches consistently show that PF is more likely to occur with increased peak rainfall intensity, finer textured material, low soil moisture variability, humid climate, and higher net primary productivity. This suggests that PF patterns could shift with projected climate changes, increasing uncertainty in predictions of groundwater recharge, water quality, and streamflow generation.

conterminous United States and Puerto Rico↗
USGS1976

Heat-flow data from southeastern Oregon

With the exception of values from two holes drilled within 2 km of Mickey Hot Springs, 17 new heat-flow values in southeastern Oregon are within or somewhat below the range one would normally expect in non-anomalous parts of the North American Cordillera. This is not surprising for a region in which most igneous rocks on the surface are 5 m.y. old or more. There is a suggestion of a thermal anomaly associated with the very young (late Pleistocene or Holocene) Diamond Craters lava field, and the thermal regime on both sides of Steens Mountain seems to be controlled, to some degree, by lateral and vertical movement of water.

Oregon↗
USGS2025

Controls on water quality below a reclaimed surface coal mine, southeastern Montana

Coal mining and reclamation can have a profound influence on hydrogeologic systems, with clear consequences for groundwater quality, yet their long-term influence on downgradient water quality over time following reclamation is less well documented. Geochemical trends were evaluated in water quality downgradient of a fully reclaimed landscape at the former Big Sky Mine in the Rosebud Creek watershed (southeastern Montana, USA), over a 3-year period (2020–2022), including bond release in 2022. Within 6 km downgradient from the reclaimed area, sulfate concentrations decreased from approximately 3500 to 1800 mg l −1 within the Miller Coulee alluvial aquifer. Major ions, δ 34 S SO4 values, and residence time tracers suggest that the observed decreases in sulfate concentration result from a combination of dilution by mixed-age inflows and incomplete transit of the high salinity plume from the mine boundary. Both bedrock and alluvial aquifers of the Rosebud Creek corridor contained contributions of millennia-old regional groundwater, which may serve to mitigate mine-derived high salinity waters. Rosebud Creek, which traverses the outflow zone of Miller Coulee in the study area, exhibited high sulfate concentrations during low flows and consistent downgradient increases in sulfate concentration. The possibility of plume dynamics in Miller Coulee suggests that the greatest water quality impacts may not yet have reached Rosebud Creek.

Montana↗
USGS2024

Changes in soil erosion caused by wildfire: A conceptual biogeographic model

Soil erosion rates after wildfire are strongly controlled by intrinsic properties such as topography, weather, climate, soil, and vegetation. These landscape and hydroclimatic properties are important in determining post-fire erosion rates; however, their influence on post-fire erosion and their interaction with the intensity of a wildfire remains uncertain. A key limitation in resolving this uncertainty is the lack of conceptual models and frameworks for organising data related to the geomorphic sensitivity of landscapes to wildfire. Our aim is to develop a framework for consolidating understanding of post-fire erosion in the context of hydroclimatic conditions which contribute to system states, for example soil and vegetation properties, and wildfire regime. The framework is developed around a simple conceptual model where the change in erosion due to wildfire is a product of change in runoff generation and sediment supply, which is strongly related to landscape net primary productivity (NPP). We hypothesised that geomorphic sensitivity to wildfire should vary as a unimodal humped relationship across a gradient of NPP, peaking at an intermediate level. To develop this framework and to test the hypothesis, we first review intrinsic soil and vegetation properties related to the supply and transport of sediment from burned and unburned hillslopes. Net primary productivity is systematically related to these intrinsic properties because it integrates many processes involved in soil and vegetation development. Empirical data indicate a trend in the change in surface runoff generation with NPP after wildfire, peaking at an NPP of approximately 15 Mg C ha −1 y −1 . A simple model of fuel availability and soil heating are correlated with a similar “humped” trend in sediment supply. These results are consistent with our conceptual model, which indicates that sediment supply and runoff contribute towards a distinct peak in wildfire effects on erosion at an intermediate level of NPP. We propose that landscapes of intermediate NPP typically have the highest quantity of fuel available to burn, which cause large changes to the soil surface properties. Landscapes at intermediate NPP also tend to produce intrinsic soil and vegetation properties that promote erosion after wildfire. The interplay between these short and long-term landscape characteristics is strongest at intermediate levels of NPP. Our proposed biogeographic model of geomorphic sensitivity to wildfire was supported by erosion data from burned hillslope and zero-order catchments studies from a range fire-prone landscapes in Australia and North America. Our proposed conceptual model will help identify areas most vulnerable to post-fire erosion changes.

Geomorphology↗
USGS2025

Carbonatite-hosted residual REE deposits

Rare earth elements (REEs) occur in magmatic rocks but are especially enriched in carbonatite and alkaline silicates. If these rocks are chemically weathered, then the REEs may become further enriched within the regolith developed from these rocks. Primary magmatic REE minerals, as well as the various carbonate minerals and apatite, provide the REEs which, under pervasive chemical weathering, are incorporated within low-temperature REE minerals forming within the regolith. Many of these minerals, as well as their textures, are characteristic of this mode of formation. Lateritic conditions of weathering are instrumental in producing a thick, weathered, or regolith, profile, and the roles of sulfide oxidation, fluctuating groundwater tables, and downward mass wasting due to carbonate dissolution are identified as the most important controls on REE enrichment in the regolith.

Book chapter↗
USGS2025

Physicochemical properties and bioreactivity of sub-10 µm geogenic particles: Comparison of volcanic ash and desert dust

Exposure to ambient particulate matter (PM) with an aerodynamic diameter of <10 μm (PM 10 ) is a well-established health hazard. There is increasing evidence that geogenic (Earth-derived) particles can induce adverse biological effects upon inhalation, though there is high variability in particle bioreactivity that is associated with particle source and physicochemical properties. In this study, we investigated physicochemical properties and biological reactivity of volcanic ash from the April 2021 eruption of La Soufrière volcano, St. Vincent, and two desert dust samples: a standardized test dust from Arizona and an aeolian Gobi Desert dust sampled in China. We determined particle size, morphology, mineralogy, surface texture and chemistry in sub-10 μm material to investigate associations between particle physicochemical properties and observed bioreactivity. We assessed cellular responses (cytotoxic and pro-inflammatory effects) to acute particle exposures (24 hr) in monocultures at the air-liquid interface using two types of cells of the human airways: BEAS-2B bronchial epithelial cells and A549 alveolar type II epithelial cells. In acellular assays, we also assessed particle oxidative potential and the presence of microorganisms. The results showed that volcanic ash and desert dust exhibit intrinsically different particle morphology, surface textures and chemistry, and variable mineralogical content. We found that Gobi Desert dust is more bioreactive than freshly erupted volcanic ash and Arizona test dust, which is possibly linked to the presence of microorganisms (bacteria) and/or nanoscale elongated silicate minerals (potentially clay such as illite or vermiculite) on particle surfaces.

GeoHealth↗
USGS2024

Comparing subduction ground-motion models to observations for Cascadia

We evaluate Cascadia subduction ground-motion models (GMMs), considered for the 2023 US National Seismic Hazard Model (NSHM) update, by comparing observations to model predictions. The observations comprise regional recordings from intraslab earthquakes, including contributions from 2021 and 2022 events in southern Cascadia and global records from interface earthquakes. Since the 2018 NSHM update, new GMMs for Cascadia have been published by the Next Generation Attenuation (NGA)-Subduction Project that require independent evaluation. In the regional intraslab comparisons, we highlight a characteristic frequency dependence for Cascadia data, with short periods having lower ground motions and longer periods being comparable to other subduction zones. We evaluate differences in northern and southern Cascadia and find that the NGA-Subduction GMMs developed using southern Cascadia data perform better in this region than the model that did not consider these data. We compare ground-motion variability in Cascadia with the NGA-Subduction model predictions and find differences at short periods ( T = 0.1 s) due to the use of global versus regional data in the development of these models. Moreover, the within-event component of aleatory variability from the GMMs overpredicts the standard deviation of Cascadia recordings at very short periods ( T < 0.05 s). Using global interface earthquakes as a proxy to evaluate the Cascadia GMMs, we find long-period overprediction from a simulation-based GMM and some of the empirical GMMs. When comparing recent observations, we find a similar misfit to GMMs and the 2010 and 2022 Ferndale earthquakes. Finally, we observe different basin amplification factors arising in different subsets of the data, which indicate that differences in basin factors between empirical GMMs could arise from the data selection choices by the developers. As part of evaluating the regional basin terms, we apply basin amplification factors from the magnitude 9 Cascadia earthquake simulations to the empirical GMMs for interface earthquakes. The comparisons presented in this study indicate that the NGA-Subduction GMMs for Cascadia perform well relative to observations and older subduction GMMs.

British Columbia, California, Oregon, Washington↗
USGS2026

Only very strong shaking can break a tree

A handful of credible accounts, both recent and historical, describe trees being broken during earthquakes. Apart from landslides and tsunamis, such accounts are uncommon, suggesting that trees break only when ground acceleration is extreme or when trees are especially vulnerable to breakage. Here, I consider the question: what ground acceleration is required to break a living tree? I present simplified first‐principles calculations, considering the properties of common trees, focusing on tall, thin trees that can be approximated as uniform vertical cantilevers. The results suggest that, possibly excepting tall trees with heavy canopies, standing trees can only be broken by shaking approaching or exceeding 1 g , as anticipated by the rarity of documented accounts of snapped trees. I briefly consider several well‐documented instances of tree damage during strong earthquakes, drawing inferences about local ground motions.

The Seismic Record↗
USGS2026

Modeling legacy nitrogen transport under instantaneous, steady-state, and transient groundwater flow conditions

In hydrologic settings where groundwater discharge contributes substantially to surface waters, legacy nitrogen in groundwater can confound surface water nitrogen loads estimated exclusively from current terrestrial sources. Additionally, legacy nitrogen in groundwater can contribute to lagged responses to nitrogen management efforts. Some methods of estimating groundwater contributions to surface water nitrogen loads account for legacy nitrogen, while others do not. The resulting differences are rarely quantified. We used a numerical modeling framework to compare three methods of estimating time-varying annual groundwater nitrogen loads to surface water receptors on eastern Long Island, New York. The instantaneous load method used steady-state contributing areas and includes no temporal groundwater lag. The second method used numerical simulations of nitrogen loads under steady-state flow, which captures groundwater transport lags but omits the annual variability in transient hydrologic stresses. The third method numerically simulated both transient groundwater flow and nitrogen transport to explicitly capture the effects of legacy nitrogen in groundwater. Depending on antecedent nitrogen and hydrologic conditions, historical nitrogen loads estimated from the numerical simulations were sometimes similar (<10% difference) and other times substantially different (±100%) from the instantaneous load estimates. Additionally, simulated future surface water nitrogen loads responded asymptotically over several decades following reductions in terrestrial nitrogen sources, further highlighting the effect of groundwater transport lag times. The comparison of the three methods, quantification of historical interannual variability, and prediction of lagged responses to nitrogen source reductions provide important context for decision makers using estimated groundwater nitrogen loads to help evaluate nitrogen management efficacy.

New York↗