Geology Reports⌕ Search

SEARCH · Geology Reports

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

Climate-driven sulfate export in alpine watersheds may stimulate methylmercury production

Climate change is increasing sulfate export and changing wetland extent in mountain regions. These changes may increase microbially mediated production of the neurotoxic substance methylmercury due to enhanced sulfate metabolism in mountain environments. Here, we assess methylmercury concentrations and formation rates across high-elevation wetlands in the Colorado Rocky Mountains. We also investigate sulfate controls on methylmercury production within subalpine peatlands by amending soils with sulfate to mimic increased stream export of sulfate from the alpine zone and measuring methylmercury formation rates for different sulfate treatments. We found that subalpine peatlands have statistically significant higher methylmercury concentrations and formation rates compared to alpine, mineral-soil wetlands. Methylmercury production in subalpine peatlands also increased significantly ( p < 0.05) following sulfate additions; the highest rates occurred in sediments with intermediate extractable sulfate concentrations (∼0.60–1.4 mg sulfate g −1 dry soil). Our study is the first to identify soil sulfate-related thresholds for methylmercury production and sulfate-limitation of methylmercury production in subalpine peatlands. These findings highlight important linkages between climate-driven mineral weathering and mercury cycling in mountain regions globally.

Colorado↗
USGS2025

A soil velocity model for improved ground motion simulations in the U. S. Pacific Northwest

Near-surface seismic velocity structure may significantly impact the intensity, duration, and frequency content of ground shaking during an earthquake. In this study, we compile 649 shear wave velocity (Vs) profiles throughout the U.S. Pacific Northwest and southern British Columbia (PNW) and use these measured profiles to develop a representative soil velocity model for four major Holocene soil provinces: Puget Lowlands, Willamette Valley, fill and alluvium, and `other' soils. The resulting soil velocity model shows good agreement to measured data for a wide range of site conditions, with variability between different geologic domains reflecting fundamental differences in depositional environments. We then show that using this regional soil velocity model in simulations of the 2001 M6.8 Nisqually, Washington earthquake improves the fit to observed high-frequency (≥ 0.5 Hz) ground motions in the Puget Sound region compared to simulations that do not incorporate shallow (≤ 200 m) seismic velocity structure. Overall, this work shows that incorporating localized soil velocity profiles into seismic velocity models is important for accurately estimating high-frequency ground motion and regional seismic hazard in earthquake simulations. Future earthquake simulations and hazard studies in the PNW could incorporate these soil velocity profiles to capture the region's distinct site response characteristics.

Washington↗
USGS2026

The 2025 Puerto Rico and Virgin Islands U.S. National Seismic Hazard Model Update: Ground motion model selection and comparison

We evaluate, select, and describe the ground-motion models (GMMs) used in the 2025 update of the U.S. National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI). We identify the most appropriate models that align with GMM selection criteria for use in the PRVI region to improve the accuracy of seismic hazard assessments. The update incorporates globally applicable GMMs suited for the active crustal and subduction earthquakes in the region. We include region-specific adjustments to these GMMs derived from local site response analyses derived from ground motion records. The unadjusted and regionally-corrected GMMs are combined to create a robust model for predicting median ground motion. The model integrates epistemic uncertainty through a median ground motion logic tree that accounts for variations in magnitude and distance. This study compares the GMMs selected for the 2025 PRVI NSHM, including both as-provided and regionally adjusted NGA-West2 and NGA-Subduction models, with those used in the 2003 PRVI NSHM. We evaluate how changes in model selection, weighting, aleatory variability, and epistemic uncertainty influence seismic hazard estimates. Trends with distance, magnitude, and spectral period are analyzed to evaluate how the scaling behavior of the newer GMMs differs from that of earlier models. Relative to the GMMs used in the 2003 NSHM for this region, the 2025 models generally predict lower ground motions. Comparisons with additional GMMs indicate that the adjustments applied for PRVI are consistent with regional-specific modifications developed elsewhere globally. The increase in aleatory variability and epistemic uncertainty in the 2025 update results in a notable increase in hazard levels from these wider uncertainty bounds. These changes can result in as much as a 10%–20% variation in probabilistic ground motion at the 2% in 50 years exceedance level for hazard maps computed across the region for representative site classes and periods.

Puerto Rico, Virgin Islands↗
USGS2026

The damability function: A probabilistic approach to regional landslide dam susceptibility analysis applied to the Oregon Coast Range, USA

Landslides can dam rivers and require rapid response to mitigate catastrophic outburst floods. Here, we present a workflow to map landslide dam formation susceptibility at a regional scale. We define a probabilistic function that combines river valley width and landslide volume to efficiently determine the likelihood of a landslide dam or “damability”. We combine damability values with landslide susceptibility to estimate landslide dam susceptibility. The valley width measurements are automated using a new elevation threshold-based algorithm. Landslide volume is represented as a statistical distribution from mapped landslides. We validate and apply our approach to the Oregon Coast Range, USA and find that 36 % of river stretches exceed a dam potential threshold; these are in river headwaters and steeper terrain, which in this case correlate with more resistant lithologies. We also estimate volumes of the potential dammed lakes and find that most rivers with high dam susceptibility are less likely to impound large lakes because they have low drainage areas. However, widespread susceptibility, and the potential impacts from exceptionally large landslides, suggest that this hazard should be considered in the Pacific Northwest. The damability function workflow can ingest new data and be applied more broadly to assess future landslide dam hazards.

Oregon↗
USGS2026

Characterizing changes in postfire debris-flow hazard as burned areas recover

Emergency assessments of postfire debris-flow hazards that are performed by the U.S. Geological Survey (USGS) provide estimates of debris-flow likelihood and rainfall triggering conditions that are used for evaluating and managing runoff-generated debris-flow hazards in recently burned areas throughout the western United States. Although the immediate postfire period, within roughly one year after fire, is typically the most susceptible to runoff-generated debris flows, the hazard evolves in time and space as the burned area recovers. The recovery trajectory a given burned area will take depends on local climate and weather and can be difficult to predict. Some burned areas recover quickly, whereas others experience debris flows for multiple years after fire. As a result, extending our ability to update debris-flow likelihood estimates and rainfall thresholds based on observed recovery of the burned area would be beneficial. We present a method for multi-year runoff-generated debris-flow hazard assessment that leverages the USGS “M1” debris-flow likelihood model and integrates updated, satellite-derived, normalized burn ratio data to estimate vegetation recovery. We predict recovery-aware rainfall thresholds and validate them against a multi-year debris-flow hazard prediction and could be adapted for use with other debris-flow models that incorporate burn severity data.

Arizona, California, Colorado, New Mexico, Washing↗
USGS2025

No evidence for an active margin-spanning megasplay fault at the Cascadia Subduction Zone

It has been previously proposed that a megasplay fault within the Cascadia accretionary wedge, spanning from offshore Vancouver Island to Oregon, has the potential to slip during a future Cascadia subduction zone earthquake. This hypothetical fault has major implications for tsunami size and arrival times and is included in disaster-planning scenarios currently in use in the region. This hypothesis is evaluated in this study using CASIE21 deep-penetrating and U.S. Geological Survey high-resolution seismic reflection profiles. We map changes in wedge structural style and seismic character to identify the inner-outer wedge transition zone where a megasplay fault has been previously hypothesized to exist and evaluate evidence for active faulting within this zone. Our results indicate that there is not an active, through-going megasplay fault in Cascadia, but instead, the structure and activity of faulting at the inner-outer wedge transition zone is highly variable and segmented along strike, consistent with the segmentation of other physical and mechanical properties in Cascadia. Wedge sedimentation, plate dip, and subducting topography are proposed to play a major role in controlling megasplay fault development and evolution. Incorporating updated megasplay fault location, geometry, and activity into modeling of Cascadia earthquakes and tsunamis could help better constrain associated hazards.

British Columbia, Oregon, Washington↗
USGS2026

Variable pH habitats could help prepare stone crabs for coastal acidification

Coastal acidification is being exacerbated by terrestrial organic inputs, especially after high precipitation events. Florida’s rainy season coincides with stone crab ( Menippe mercenaria) reproduction, and pH extremes could limit future harvests by reducing reproductive output. Populations that experience pH variability can serve as “natural laboratories” for estimating the tolerance of a species to coastal acidification. Here, we conducted a series of experiments to determine if ovigerous stone crabs conditioned in more variable pH habitats (seagrass) would result in faster embryonic development, greater hatching success, and higher larval survival relative to crabs conditioned in habitats with lower pH variability (sandy habitats). After field conditioning, crabs were transported to the laboratory and randomly acclimated to either a control pH (pH = ~ 7.90) or a reduced pH condition (pH = ~ 7.60) until larval release. The rate of embryo development was slower in the laboratory reduced pH treatment, however, there were no observable field condition effects on embryo development rate. Crabs conditioned in the more pH variable seagrass habitat did have greater hatching success and higher larval survival than crabs in less pH variable sandy habitats; however, larval survival was low across all treatments. These results suggest that the pH variability experienced in seagrass habitats during brooding may serve as a mechanism for stone crabs to acclimatize to extremes in seawater pH.

Marine Biology↗
USGS2025

2023 Earthquake Ground-Motion Workshop for the Central and Eastern United States, with a focus on the Gulf and Atlantic Coastal Plains—Agenda and abstracts

The U.S. Geological Survey held a virtual workshop December 7–8, 2023, to share research and ideas about earthquake ground motions in the Central and Eastern United States, with a focus on the Atlantic and Gulf Coastal Plains. The workshop was organized to learn about potential regionalization of ground-motion characteristics (source, path, and site), consider new explanatory variables for site response, and hear and discuss updates on ground-motion research on the Atlantic and Gulf Coastal Plains. The workshop was organized into a series of contributed presentations and three panel discussions held during 2 days. This report documents the agenda, contributed abstracts, and panel summaries.

Scientific Investigations Report↗
USGS2026

Age and colony variation in Adélie penguin metapopulation vital rates: Insights from a 25-year mark–recapture study

Understanding how vital rates vary with age, life-history stage, and among populations is fundamental for predicting the demographic consequences of environmental change, especially in longer-lived species with complex life histories. These species often exhibit delayed maturity and iteroparity, making nuanced demographic insights critical for assessing their long-term viability. This study investigated age- and colony-related variation in Adélie penguin vital rates including survival, recruitment, and breeding propensity. We used mark–recapture data collected over 25 years (1996–2020) from three Adélie penguin breeding colonies that differed in population sizes and trends but comprised a metapopulation located at capes Royds, Bird, and Crozier on Ross Island, Antarctica. We used multi-state models to estimate survival and detection rates relative to reproductive state and breeding colony and estimated transition probabilities reflecting movements between reproductive states and colonies. Apparent survival varied by reproductive state, colony, and age and averaged 0.80 (SD = 0.02) at Bird, and 0.72 (SD = 0.03) and 0.73 (SD = 0.03) at Crozier and Royds, respectively, for pre-breeders age 2–7 years with strong declines in pre-breeder survival after age 8. We observed less age-related variation in survival of breeders and non-breeders, but we observed differences between colonies with lower survival for breeders (0.72 to 0.80) compared to non-breeders (0.75 to 0.82). The average probability of surviving the first 2 years after fledging ranged from 0.43 (SD = 0.14) at Royds and Crozier (0.43, SD = 0.10) to 0.55 (SD = 0.15) at Bird. Movement between colonies was highest for pre-breeders (0.00%–12.00% depending on age and colony) and lowest for breeders (<0.20%). We observed the lowest age-related recruitment rates at Royds, with recruitment at Crozier almost twice as high, and intermediate at Bird. Breeding propensity was highest at Crozier and lowest at Bird. Colony-specific variation in vital rates likely contributed differently to population trajectories, suggesting that care must be taken to extrapolate vital rate estimates across colonies even within a metapopulation. These findings also highlight the importance of considering age, life-history stage, and geographic variation when assessing population vital rates in species with complex life histories.

Frontiers in Ecology and Environment↗
USGS2024

The 2023 U.S. 50-state National Seismic Hazard Model: Changes in 2023 compared to 2018 ground motions

We present the 2023 U.S. National Seismic Hazard Model (NSHM) for all 50 states that applies new smoothed seismicity, fault rupture, and ground motion models. New data and methods are introduced in the 2023 earthquake rupture forecast that include: new earthquake catalogs - excluding induced earthquakes, alternative declustering methods, spatially smoothed seismicity distributions, full-catalog scaled rates to account for aftershocks, updated CEUS-WUS attenuation boundary, new magnitude-scaling equations, new geodetic and geologic deformation models, and alternative fault system solutions accounting for a more complete representation of epistemic uncertainty potential for earthquake generation in Alaska, Hawaii, and the conterminous U.S. Improved ground motion models consider new Next Generation Attenuation NGASubduction, modified NGA-East, and adjustments to account for regional biases in ground shaking observations. Semi-empirical and 3D simulations of ground motion are applied to account for shaking at 21 oscillator periods, 2 peak motions, and 8 site conditions. Site effects models are constructed for western U.S. basins (Seattle, Portland/Tualatin, San Francisco, Central Valley of California, Los Angeles, and Salt Lake City) and for sites with deep sedimentary wedges found across the central and eastern U.S. Gulf Coast and Atlantic coastal plain regions. These models result in substantial changes compared to the older NSHMs and are differentiated for the earthquake rupture forecast and ground motion model changes to display sensitivities and impacts.

Conference Paper↗
USGS2026

Cascadia Subduction Zone science: Call for the next generation community seismic velocity model

The Cascadia subduction zone (CSZ) hosts major seismic and tsunami hazards, yet key questions persist about the relationship between margin structure, fluid distribution, episodic tremor and slip, shallow megathrust behavior, shaking and tsunamigenesis, and the resulting hazard estimates. Addressing these problems requires an empirically grounded, three‐dimensional seismic velocity model to illuminate subsurface structure and properties and to provide a basis for geophysical studies such as earthquake simulations and ground‐motion estimation. In May 2024, the National Science Foundation‐funded Cascadia Region Earthquake Science Center (CRESCENT) community velocity model (CVM) working group, with U.S. Geological Survey and regional partners, convened a workshop to identify priorities for such a model. Participants emphasized the features necessary for addressing key science questions, including implementing findability, accessibility, interoperability, and reusability (FAIR) access, capturing along‐strike and along‐dip structural heterogeneity, resolving shallow offshore–onshore structure, constraining elastic properties and quantifying their uncertainties for numerical wave propagation simulations, their validation benchmarks, and supporting associated accurate earthquake ground‐motion simulations and hazard assessments. This article describes the priorities defined in the workshop, and a description of how, guided by these needs, CRESCENT plans to develop multiple generations of a CVM to advance CSZ science and improve seismic and tsunami hazard modeling across the Pacific Northwest. The CVM will span the CSZ from the surface to ∼100 km depth, offshore and east of the Cascades into Idaho (∼132°–110° W) and the southern and northern tectonic regime transitions (∼36°–52° N) to capture the entire tectonic system as well as its surroundings.

Cascadia Subduction Zone↗
USGS2026

Asynchronous landslide seasonality across the United States

Mid-range landslide outlooks can facilitate weather-related landslide preparedness and disaster response planning, but seasonal landslide activity remains poorly quantified at continental scales. Leveraging >55,000 reported landslides from across the United States (U.S.), we used circular statistics to quantify landslide seasonality in 67 National Weather Service County Warning Areas (CWAs). We found regional differences in landslide season timing and duration, with transitions between domains variably corresponding to climate class or river basin. We assessed differences in seasonality by movement type for slides, flows, and falls, detecting apparent, but uncertain, differences between slide and fall seasonalities in 27 of 35 (77%) of CWAs with both types reported. In the Pacific Northwest, where long records exist, we found a credible shift toward a later mean landslide season in western Washington from 1990 to 2020, but no trend in western Oregon. Our results can provide emergency planners a resource to assess seasonal landslide probability nationwide.

Geophysical Research Letters↗
USGS2026

Critical minerals memory match game

An educational information packet about the 2025 List of Critical Minerals, which includes a memory match game about select critical minerals and how they are used.

General Information Product↗
USGS2026

Temperate wetlands lose climate-cooling capacity under warming

Temperate inland wetlands have long been viewed as climate regulators, with sustained CO₂ sequestration offsetting continuous CH₄ emissions over centennial–millennial timescales. However, the net climatic effect of these opposing fluxes remains uncertain under anthropogenic warming as they differ in both magnitude and atmospheric lifetime. We combine 19,000 CH₄ chamber measurements from 202 temperate inland wetlands with carbon burial rates from 16 dated sediment cores in North America’s Prairie Pothole Region to estimate greenhouse gas dynamics over two centuries. Here we show that CH₄ emissions rise by ~1.39 kg CH₄ ha⁻¹ yr⁻¹ over the historical period and are projected to approach 950 kg CH₄ ha⁻¹ yr⁻¹ under high-forcing scenarios by 2100, while CO₂ sequestration is expected to decline by ~0.3 kg CO₂ ha⁻¹ yr⁻¹, driving a shift toward persistent net CO₂-equivalent source behaviour and a weakening of the long-term cooling effect that wetlands have provided under natural conditions.

Prairie Pothole Region↗
USGS2026

Bayesian belief network model to predict human-wildlife conflict in protected areas

Human-wildlife conflict (HWC) poses a pervasive global challenge, affecting livelihoods and threatening biodiversity. To better anticipate and mitigate HWC risk, we developed a large-scale predictive model using a Bayesian Belief Network (BBN). We surveyed 1,011 park rangers across 135 terrestrial protected areas in three Andean countries, documenting recent HWC incidents involving wildlife persecution or killing, livestock depredation, crop damage, or threats to human safety and property. We identified key drivers of HWC risk, including governance, wildlife acceptance, participation, and habitat quality. A sensitivity analysis revealed that enhancing governance and improving wildlife acceptance could reduce HWC risk by > 85%. The BBN model demonstrated scalability, effectively identifying strategies to reduce HWC risk at multiple scales, from individual protected areas to national networks. Our findings highlight the importance of strengthening governance, increasing wildlife acceptance, and enhancing community participation in conservation efforts. BBNs provide a flexible, cost-effective, and data-driven tool to guide protected areas and wildlife managers in monitoring, anticipating, and making informed decisions to mitigate conflict and promote coexistence.

Scientific Reports↗
USGS2026

A 481 m-high landslide-tsunami in a cruise ship-frequented Alaska fjord

Early in the morning of 10 August 2025, a >64 × 10 6 –cubic meter landslide struck Tracy Arm fjord in Alaska. The landslide was preconditioned by glacial retreat caused by climate change. The resulting 481-meter runup megatsunami followed an initial 100-meter-high breaking wave traveling at >70 meters per second. The landslide was preceded by several days of microseismicity, which increased in rate and magnitude until ~1 hour before failure. The landslide produced globally observed long-period seismic waves equivalent in size to a moment magnitude 5.4 earthquake. A long-period (~66 second) global seismic signal, produced by a landslide-induced seiche trapped within the fjord, persisted for up to 36 hours, the second time a days-long seiche had thus been observed. With fjord regions increasingly visited by cruise ships, and climate change making similar events more likely, this unanticipated, near-miss event highlights the growing risk from landslides and tsunamis in coastal environments.

Alaska↗
USGS2026

Unprecedented burning in tropical peatlands during the 20th century compared to the previous two millennia

Tropical peatland wildfire incidence has risen in recent decades, driven by drainage for land use and intensified by severe droughts with global climate change. These disturbances have altered vegetation structure, disrupted ecosystem functioning, and increased carbon emissions, particularly in Southeast Asia. However, the long-term history and characteristics of wildfires in tropical peatlands remain largely unknown. Here, we compiled fifty-eight macro-charcoal records from peatlands across the tropics, ranging from lowland forested to montane peatlands, to assess millennia-scale changes and controlling factors of tropical peatland burning. We divided the datasets into four main sub-regions: Neotropical, Afrotropical, Indomalayan and Australasian ecoregions to explore regional variability. Tropical peatlands had high burning levels between 0 and 850 ce , followed by a relatively low and stable period until a marked increase during the 20th century. The general trend in tropical peatland burning follows changes in global temperature, and climate variables that control the length and severity of drought events have a notable influence on peat burning before 1900 ce . During the 20th century, regional differences were observed, with declining fire trends in the Neotropical and Afrotropical regions and increasing fire trends in the Indomalayan and Australasian regions. This difference is likely attributable to human activities, and such intervention is also evident in palm swamps and hardwood swamps under similar wet, weakly seasonal climates. With the increase in anthropogenic pressures on peatlands and greater climate variability, future wildfires in peatlands are likely to become more frequent and widespread across all tropical ecoregions. Conservation and sustainable land-use practices could be used to mitigate and control peatland burning and protect these carbon-rich sinks.

Global Change Biology↗
USGS2026

Rapid earthquake magnitude classification via P-wave strains from borehole strainmeters and Distributed Acoustic Sensing

Distributed Acoustic Sensing (DAS) offers a promising approach for earthquake early warning (EEW) in settings where seismic networks are costly to maintain. By repurposing fiber-optic cables as dense strainmeter arrays, DAS enables real-time earthquake detection wherever those fibers are accessible. However, poor azimuthal coverage and challenges in estimating magnitude from strain measurements remain key hurdles in applying for earthquake monitoring. Here, we develop a machine learning method to distinguish large (M≥5.4) earthquakes from smaller ones within the first 4 seconds of a strain waveform after a P-wave arrival without determining location. Using ensemble decision tree models trained on borehole strainmeter data (3.5≤M≤7.1) and tested on onshore DAS waveforms (including the 2024 M7 Offshore Cape Mendocino earthquake), we find that low-frequency (0.2–0.5 Hz) continuous wavelet transform coefficients are the strongest predictors of magnitude, in addition to strain amplitude. Both DAS and borehole strainmeters effectively capture long-period strain signals, making these findings valuable for EEW systems. Our method shows high precision compared to the real-time EEW system, ShakeAlert®, supporting the position that DAS is a viable technology for earthquake monitoring and magnitude classification.

California↗