USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Water Process Engineering”

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,278 records · Page 8Linked to original sources
USGS2025

Multi-scale predictors of Northern Long-eared Bat (Myotis septentrionalis) occupancy in the United States

Historically, Myotis septentrionalis (Northern Long eared Bat) was among the most common forest-interior species in North America. Largely due to high mortality from white-nose syndrome, this species has experienced severe population declines across its range. To create an updated species distribution map representing summer occupancy probabilities from 2017 to 2022, we integrated stationary acoustic data with live-capture data from the database of the North American Bat Monitoring Program into a multi-scale, multi-method occupancy modeling framework. Our results provide data-driven predictions with quantified uncertainty for summer occupancy probabilities for Northern Long-eared Bats at 2 spatial scales across the range of the species, while also accounting for inherent observation biases (e.g., imperfect detection).

Journal of North American Bat Research↗
USGS2026

Evolution of arsenic speciation during thermal treatment simulating wildfire heating in arsenic-rich sediments

Understanding arsenic transformations during wildfires can help better constrain the environmental impacts of increased wildfire intensity, frequency, and burned area. Previous studies have monitored the evolution of metal(loid) speciation at specific time/temperature endpoints, hampering the comprehension of the dynamic transformations of metal(loid)s during wildfires. Here, a novel approach, namely in situ time-resolved X-ray absorption near edge structure (TR-XANES) spectroscopy, was used for the first time to monitor the heat-induced redox transformations of arsenic (As) in two As-rich sediments with different characteristics. No arsenic transformations were observed in either sediment at room temperature or when the sediments were heated to a maximum temperature of 400 ºC. When heated to a maximum temperature of 600 or 700 ºC, As underwent complex, dynamic, and partially reversible redox transformations, the extent of which varied with sediment properties, initial As speciation, and heating temperature and duration. In the case of sediments initially containing As(V) and a low sulfide/sulfate ratio, the As(V) was reduced to As(III) over a short period of time, followed by immediate reoxidation of As(III) to As(V). The extent of the transient As reduction increased with increasing maximum temperature. The final As speciation at the end of the heating experiment was not substantially different from that at the beginning of the experiment. In contrast, in sediments that were best described by a more complex combination of reference compounds (As(V) coprecipitated with or adsorbed on iron oxide, arseniosiderite, arsenopyrite) with a high sulfide/sulfate ratio, As(V) phases were reduced to As(III) and As(III) was further reduced to As(-I) (at 700 ºC, only), followed by delayed reoxidation of As(-I) to As(III) (at 700 ºC, only) and reoxidation of As(III) to As(V). Approximately 9% and 26% of all As remained as As(III) ( i.e., was not reoxidized to As(V)) at the end of the experiment when heated to a maximum temperature of 600 and 700 ºC, respectively, which is higher than the relative abundance of As(III) in the native sediment. These fire-induced As transformations can potentially enhance As mobilization to surface water and groundwater, posing risks to environmental and human health.

California↗
USGS2026

Statewide surficial geologic map of Nebraska underscores Quaternary landscape evolution from the High Plains to the Central Lowland

Surficial geologic mapping in Nebraska has been conducted primarily at the 1:24,000 scale since the mid-1990s, although there have also been limited efforts to map generalized Quaternary and surficial geology within the state and the region. We compiled and evaluated disparate maps (1:24,000–1:1,000,000) and datasets—including geologic, soil and soil parent material, and geomorphic maps as well as LiDAR derivatives—to produce a single 1:500,000 scale surficial geologic map of Nebraska that is the first of its kind. This new map provides a coarse-scale surficial geologic map that will be incorporated into a nationwide U.S. Geological Survey Quaternary geologic map. It also reflects the variation and uniqueness of physical landscapes in the state, where the Great Plains and Central Lowland physiographic provinces meet, further developing a richer interdisciplinary understanding of regional geomorphology in the heart of North America.

Nebraska↗
USGS2026

Synthesizing a twelve-year sediment trap time series of planktic foraminiferal flux in the Gulf of America (Mexico)

Sediment trap time series provide powerful frameworks for testing hypotheses about planktic foraminiferal assemblage composition, seasonality, and geochemical responses to environmental variability, all of which are central to improving paleoceanographic reconstructions. We present high resolution foraminiferal assemblage data from a long-running sediment trap (2008–2020) in the northern Gulf of America (Mexico). This study summarizes the species composition, seasonality, and size distribution of the fifteen most abundant species of planktic foraminifera, which account for 98% of total flux. Foraminiferal flux peaks in winter and reaches a minimum in summer, following the seasonal pattern of primary production in the northern Gulf. Winter assemblages are dominated by non-spinose taxa, whereas spinose taxa prevail during summer. Across nearly all species, average monthly test size covaries with temperature, independent of seasonal flux trends. Notably, Trilobatus sacculifer and Neogloboquadrina dutertrei show a significant decline in relative abundance and flux during 2017–2020, nearly disappearing from the assemblage.

Journal of Foraminiferal Research↗
USGS2025

National population exposure and evacuation potential in the United States to earthquake-generated tsunami threats

Previous efforts to characterize tsunami threats to people have focused primarily on individual scenarios in specific areas but have not recognized multiple scenarios across an entire country. This study addresses this gap by quantifying population exposure and evacuation potential in the United States to 102 earthquake-related, tsunami-hazard zones, including 92 local scenarios, 8 distant scenarios, and 2 probabilistic products. Geospatial path-distance modeling quantified evacuation potential and the influence of departure delays. We focused on residents to support other national, multi-hazard risk analyses. Millions of residents are in distant-tsunami zones, and hundreds of thousands of residents are in local-tsunami zones. In 41 scenarios, there is at least one resident that may have insufficient time to evacuate before wave arrival. Tens of thousands of residents may have insufficient time to evacuate from local tsunamis that impact the U.S. Pacific Northwest or Puerto Rican coastlines. The largest improvements in evacuation potential may come from reducing departure delays in some areas but may involve vertical-evacuation structures or changing land use in other areas.

International Journal of Disaster Risk Reduction↗
USGS2024

Self-guided decision support groundwater modelling with Python

The GMDSI tutorial notebooks repository provides learners with a comprehensive set of tutorials for self-guided training on decision-support groundwater modelling using Python-based tools. Although targeted at groundwater modelling, they are based around model-agnostic tools and readily transferable to other environmental modelling workflows. The tutorials are divided into three parts. The first covers fundamental theoretical concepts. These are intended as background reading for reference on an as-needed basis. Tutorials in the second part introduce learners to some of the core concepts parameter estimation in a groundwater modelling context, as well as providing a gentle introduction to the PEST, PEST++ and pyEMU software. Lastly, the third part demonstrates how to implement highly-parameterized applied decision-support modelling workflows. The tutorials aim to provide examples of both “how to use” the software as well as “how to think” about using the software. A key advantage to using notebooks in this context is that the workflows described run the same code as practitioners would run on a large-scale real- world application. Using a small synthetic model facilitates rapid progression through the workflow.

Journal of Open Source Education↗
USGS2025

Disturbance is the primary determinant of food chain length when the top predator is constant

Food chain length (FCL) is a primary determinant of food web structure and is hypothesized to be influenced by habitat size, productivity, and disturbance. Understanding the environmental characteristics that determine food chain length can assist in understanding how food webs may be impacted due to changes in habitats and environmental characteristics. This study examines the impact of hydrologic disturbance on stream food webs when the top predator is constant. We analyzed FCL in less disturbed groundwater flashy streams and more disturbed runoff flashy streams using stable isotopes. Despite no difference in species richness or fish density, food chains in more disturbed streams had a lower FCL compared to food chains in more stable streams. Assemblage analysis showed that flow regime and drainage area significantly impacted individual species abundances. The more disturbed runoff flashy streams had higher proportions of primary consumer fish, such as the algivorous Campostoma sp. (Stonerollers), which likely drives the reduced FCL. Drainage area and land cover had non-significant relationships with FCL. Shifting community structure due to hydrologic variability likely leads to differences in diet of Micropterus dolomieu (Smallmouth Bass), and thus a difference in FCL.

Arkansas, Missouri, Oklahoma↗
USGS2025

Total uncertainty quantification in inverse solutions with deep learning surrogate models

We propose an approximate Bayesian method for quantifying the total uncertainty in inverse partial differential equation (PDE) solutions obtained with machine learning surrogate models, including operator learning models. The proposed method accounts for uncertainty in the observations, PDE, and surrogate models. First, we use the surrogate model to formulate a minimization problem in the reduced space for the maximum a posteriori (MAP) inverse solution. Then, we randomize the MAP objective function and obtain samples of the posterior distribution by minimizing different realizations of the objective function. We test the proposed framework by comparing it with the iterative ensemble smoother and deep ensembling methods for a nonlinear diffusion equation with an unknown space-dependent diffusion coefficient. Among other applications, this equation describes the flow of groundwater in an unconfined aquifer. Depending on the training dataset and ensemble sizes, the proposed method provides similar or more descriptive posteriors of the parameters and states than the iterative ensemble smoother method. Deep ensembling underestimates uncertainty and provides less-informative posteriors than the other two methods. Our results show that, despite inherent uncertainty, surrogate models can be used for parameter and state estimation as an alternative to the inverse methods relying on (more accurate) numerical PDE solvers.

Journal of Computational Physics↗
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

Detecting snow avalanche activity using infrasound: Hooker Valley, New Zealand

Snow avalanches pose considerable hazards to people and infrastructure in alpine environments. Traditional avalanche monitoring relies on meteorological data and visual observations, which can be limited in scope and timeliness. Infrasound offers a promising complementary monitoring tool by detecting the low-frequency sound waves generated by avalanches. Here, we present infrasound and camera observations during a 50-day field campaign in the Hooker Valley of Aoraki/Mount Cook National Park, New Zealand. Our study detected seven avalanches with the cameras, whereas the infrasound system identified only one of these events, which was the largest and occurred under conditions that likely favoured infrasound propagation. The infrasound system recorded numerous other events not captured by the cameras, indicating the benefit of further investigation to determine their sources. These findings highlight the potential of infrasound technology for detecting avalanches and providing broad spatial coverage, capturing events in areas not monitored by cameras, while also showcasing limitations in infrasound capabilities. The limited detection of smaller avalanches underscores the opportunity for further research to enhance detection capabilities and understand environmental influences such as snow cover and wind noise. Overall, this study emphasises the utility of multidisciplinary monitoring techniques to improve avalanche detection in alpine environments.

Hooker Valley↗
USGS2026

Identifying overlap between native fish movements and a possible seasonal Grass Carp deterrent in the Sandusky River, Ohio, USA

Riverine fishes are vulnerable to habitat fragmentation caused by interrupted access to vital habitats. Fragmentation using artificial structures can be used to limit the spread of invasive fishes by inhibiting movement between critical habitats. Thus, tradeoffs exist between maintaining connectivity for native species and restricting movement of invasive fishes. A nonphysical deterrent has been proposed to prevent invasive Grass Carp ( Ctenopharyngodon idella ) from reaching spawning habitats in the Sandusky River, a tributary to Lake Erie, during late spring–summer. However, this timing overlaps with spawning periods of a suite of native fishes. We used acoustic telemetry and ichthyoplankton surveys to examine movement and spawning activity of multiple native species and to assess potential impacts of deterrent operation on reproduction. All species used the proposed deterrent location during the Grass Carp spawning period, and their movements corresponded with larval fish collections, indicating upstream spawning. These findings suggest the deterrent could unintentionally restrict native species’ access to essential habitats. Balancing habitat connectivity and invasive species control will likely depend on species-specific responses and careful timing of deterrent operation.

Ohio↗
USGS2026

Opportunities for the U.S. Geological Survey’s National Seismic Hazard Model to improve seismic risk assessment of critical infrastructure.

As fragility and risk modeling techniques and computational capabilities evolve, complemented by moving toward more routine and systematic seismic risk assessment of all buildings and critical infrastructure, the authors pose a few critical questions to investigate how the U.S. Geological Survey (USGS) National Seismic Hazard Models (NSHMs) can be used and enhanced further to serve such issues. In this paper, we use three examples from multiple sectors to (1) identify the role of USGS NSHMs in evaluating seismic risks to critical infrastructure, (2) quantify potential impacts from NSHM enhancements (i.e., [i] hazard curves for the vertical component of ground motion, [ii] stochastic event sets, and [iii] maps of probabilistic ground failure hazards), and (3) clarify the feasibility of relevant NSHM improvements. We illuminate that NSHMs are commonly used in location-specific performance assessments, whereas earthquake effects on critical infrastructure can be widespread across large geospatial regions. Further, we found that without the NSHM extensions considered here, risk can be severely underestimated, e.g., neglecting ground failure hazards can underestimate regional loss by a factor of two or more. Although many challenges remain, we developed example prototypes to clarify the feasibility of the NSHM extensions, which can facilitate improved management of risks to critical infrastructure.

Earthquake Spectra Journal↗
USGS2026

Variability and consistency in wildfire susceptibility: Insights from a national compilation

Background Wildfire risk in the United States is rising and remains a land management priority. The quantitative wildfire risk assessment (QWRA) framework integrates fuels, topography, weather and values at risk to estimate the potential change in value from wildfire. Within this, response functions (RFs) represent how values respond to fire intensity. These are often based on expert judgment, but variation across assessments is unclear. Aims This study uses data from the US Geological Survey (USGS) Wildfire Hazard and Risk Assessment Clearinghouse to characterize consistency and variation across categories and contexts. Methods We applied descriptive statistics to summarize RFs, using tables, box-and-whisker plots and heat maps stratified by highly valued resource or asset (HVRA) category and spatial scale. Key results RFs and value definitions vary, especially for ecosystem-related resources. Some functions, such as for buildings in the wildland–urban interface (WUI), translate well across contexts, while others require more input. Conclusions Some functions are broadly transferable, while others need customization. This analysis provides references and starting points for improvement to RFs in QWRAs. Interpretations Expanding the clearinghouse and dataset and building more transparency in expert elicitation can build trust among communities, agencies and end-users, and can support efficient use of limited resources to mitigate wildfire risk.

International Journal of Wildland Fire↗
USGS2026

Post-fire soil hydrologic response and recovery in northern California (USA)

Background Wildfires abruptly change landscapes by altering soil properties and vegetation cover. These changes are thought to reduce soil infiltration capacity, making landscapes susceptible to runoff and erosion. However, post-fire soil response is complex and likely varies across locations and time. Aims Here, we aim to understand regional post-fire soil response and recovery by tracking changes across different northern California (USA) lithology and vegetation types. Methods We conducted repeat in situ soil infiltration tests for 3 years post-fire at 31 burned and 10 unburned sites spanning the 2021 Dixie, 2020 LNU Lightning Complex, 2020 Walbridge and 2020 Glass fires. Key results Our two main findings are: (1) burned chaparral soils have increased hydraulic conductivity compared with unburned sites, and (2) infiltration rates return to pre-fire conditions within 3 years across most lithologies and vegetations. Conclusions Recovery might be generalizable by vegetation and lithology but differ regionally, making it important to identify meaningful hydrologic response units (HRUs). Multi-year studies with paired burned and unburned measurements can constrain the recovery timeline and provide information missed by observations solely of burned soils. Implications Understanding where, and for how long, soil remains susceptible to runoff and erosion can help prioritize areas and time periods most in need of mitigation.

California↗
USGS2025

3D viscoelastic models of slip-deficit rate along the Cascadia subduction zone

Interseismic deformation in the Pacific Northwest is constrained by the horizontal crustal velocity field derived from the Global Positioning System (GPS) in addition to vertical rates derived from GPS, leveling, and tide gauge measurements. Such measurements were folded in to deformation models of fault slip rates as part of the 2023 National Seismic Hazard Model (NSHM) update. Here I build upon one of the contributing models, the viscoelastic earthquake-cycle model of Pollitz [2022]. This model permits inclusion of effects of time-dependent viscoelastic relaxation within earthquake cycles (i.e., ‘ghost transients’) and laterally variable elastic and/or ductile material properties. I lever-age these capabilities to incorporate the Cascadia megathrust into Western U.S.-wide deformation models in which crustal fault slip rates are estimated simultaneously with slip deficit rates along the interplate boundary between the descending Juan de Fuca plate and North American plate. This effort includes construction of a margin-wide model of viscoelastic structure founded on the Slab 2.0 model and probes different models of the ductile properties of the surrounding oceanic asthenosphere, continental lower crust, and mantle asthenosphere. This results in new estimates of the distribution of slip deficit rate along the ∼ 1000 km long margin, highlights the importance of correcting for glacial-isostatic adjustment effects, and permits assessment of sensitivity of results to assumed ductile properties.

California, Oregon, Washington↗
USGS2025

Integrating marine historical ecology into management of Alaska’s Pacific cod fishery for climate readiness

The Pacific cod ( Gadus macrocephalus ) fishery was closed in 2020 after a rapid decline in biomass caused by the marine heat waves of 2014–2019. Pacific cod are exceptionally thermally sensitive and management of this fishery is now challenged by increasingly unpredictable climate conditions. Fisheries monitoring is critical for climate readiness, but short-term monitoring data may be inadequate for recognizing and anticipating change under rapid climate changes. We propose an interdisciplinary, marine historical ecology framework that looks to long-term records (local and traditional knowledge, history, archaeology, and paleoclimatology) to capture a long range of ecological variability and provide historical context for management. In order to connect to contemporary fisheries management, this framework must be built on a common vocabulary and an understanding of the key metrics used in fisheries stock assessments. Here, we propose metrics derived from Pacific cod stock assessment and synthesize information relevant to understanding the effects of past warming periods on cod populations across the Gulf of Alaska and Bering Sea. This case study provides a framework for thinking about how to use these historical records in the context of fisheries management under rapidly changing climate conditions.

Alaska↗
USGS2025

Borehole geophysical time-series logging to monitor passive ISCO treatment of residual chlorinated-ethenes in a confining bed, NAS Pensacola, Florida

In-situ chemical oxidation (ISCO) is a common method to remediate chlorinated ethene contaminants in groundwater. Monitoring the effectiveness of ISCO can be hindered because of insufficient observations to assess oxidant delivery. Advantageously, potassium permanganate, one type of oxidant, provides the opportunity to use its strong electrical signal as a surrogate to track oxidant delivery using time-series borehole geophysical methods, like electromagnetic (EM) induction logging. Here we report a passive ISCO (P-ISCO) experiment, using potassium permanganate cylinders emplaced in boreholes, at a chlorinated ethene contamination site, Naval Air Station Pensacola, Florida. The contaminants are found primarily at the base of a shallow sandy aquifer in contact with an underlying silty-clay confining bed. We used results of the time-series borehole logging collected between 2017 and 2022 in 4 monitoring wells to track oxidant delivery. The EM-induction logs from the monitoring wells showed an increase in EM response primarily along the contact, likely from pooling of the oxidant, during P-ISCO treatment in 2021. Interestingly, concurrent natural gamma-ray (NGR) logging showed a decrease in NGR response at 3 of the 4 wells possibly from the formation of manganese precipitates coating sediments. The coupling of time-series logging and well-chemistry data allowed for an improved assessment of passive ISCO treatment effectiveness.

Florida↗
USGS2025

Preface to the focus section on earthquake-induced ground deformation

Ground failure due to earthquake-induced liquefaction, surface fault rupture, and landsliding has caused substantial damage in the past and continues to pose risks to critical infrastructure in the future. The consideration of ground failure is central to engineering analysis and design, ranging from site-specific engineering assessments of ground failure potential to broader regional-scale ground failure assessments of near-real-time earthquake impacts or longer-term portfolio risk. Uncertainties in all aspects of the problem complicate the prediction of earthquake-induced ground failure, including the identification of hazard-susceptible geomaterials, characterizing their spatial extent and physical properties, spatio-temporal variability in groundwater conditions, characterization of earthquake ground motions, assessing ground failure severity, and linking ground failure severity to infrastructure damage and loss.

Seismological Research Letters↗