USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Waterway, Port, Coastal, and Ocean 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,722 records · Page 90Linked to original sources
USGS2025

PFAS sampling activities in the U.S. Geological Survey national networks

Per- and polyfluoroalkyl substances (PFAS), frequently called “forever chemicals,” are used for a wide variety of industrial purposes and are often found in common household and industrial items such as firefighting foams, non-stick cookware, and water-resistant materials. The contamination of water, air, and soil by PFAS is a national and global issue due to their widespread occurrence in multiple applications and resistance to biodegradation and other traditional treatment processes. Research indicates that many PFAS can be emitted to the atmosphere and transported and deposited long distances from the source. The U.S. Geological Survey (USGS) Water Resources Mission Area received funding to implement a national-scale sampling effort to assess PFAS occurrence. To follow agency directives, the National Water Quality Network (NWQN) added PFAS sample monitoring for both surface water and groundwater, and also added PFAS monitoring to selected sites in the National Atmospheric Deposition Program (NADP).

General Information Product↗
USGS2025

Observing northern high-latitude river systems to understand changes in a warming Arctic

Purpose of Review Streams and rivers are undergoing rapid change as the Arctic warms and thaws. We review recent observations in Arctic stream systems to identify ubiquitous changes and the most useful tools for observing change and exploring the underlying processes. Recent Findings Recent literature indicates increasingly significant trends in river hydrology and chemistry due to persistent warming in the Arctic and longer observational records for analysis. However, regional differences in the magnitude and direction of these trends persist. We also observe thresholds in ground thaw and surface–groundwater interactions that can impact river hydrology and chemistry. Summary Warming and thaw are occurring rapidly at high latitudes, resulting in increasing, yet variable responses in stream systems across regions and scales. These differences highlight the need for long-term records and an interdisciplinary approach to explain trends and predict future states. Stream systems respond to multiple landscape changes related to hydrology (changing precipitation and subsurface flow), geology (ground thaw dynamics), and ecology (vegetation change).

Current Climate Change Reports↗
USGS2025

Evidence for nonlocal sediment transport on hillslopes from fault scarp morphology

Hillslope sediment transport processes such as bioturbation, rainsplash, and granular mechanics occur across the entire planet. Yet, it remains uncertain how these small-scale processes act together to shape landscapes. Longstanding hillslope diffusion theory posits that hillslope processes are spatially limited, whereas new concepts of nonlocal sediment transport argue otherwise. However, each theory produces subtly different, but distinct, predictions for the evolution of fault scarps. We use the topographic change of fault scarps to demonstrate that hillslope processes produce nonlocal sediment transport. Analysis of a global compilation of 340 dated single-earthquake scarp profiles reveals a statistically significant ( p < 0.05 ) relationship between scarp age and scarp asymmetry, here defined as the ratio of imaginary to real components of the Fourier transform of absolute slope. Numerical simulations show that nonlocal models predict this relationship, whereas hillslope diffusion models do not. To further investigate this result, we examined the depositional geometry of a well-exposed colluvial wedge along the Wasatch fault in central Utah, United States. Our quantitative comparison between the exposure and numerical simulations reveals better agreement with the nonlocal model. Nonlocal sediment transport theory appears to best capture the physics of how hillslope processes shape fault scarps, yet hillslope diffusion provides a useful approximation in many cases. As the processes that act on fault scarps are nearly identical to those acting on hillslopes, our results provide evidence supporting nonlocality as a generalized model of hillslope sediment transport.

Utah↗
USGS2026

Is satellite-derived bathymetry vertical accuracy dependent on satellite mission and processing method?

This research focusses on three satellite-derived bathymetry methods and optical satellite instruments: (1) a stereo photogrammetry bathymetry module (SaTSeaD) developed for the NASA Ames stereo pipeline open-source software (version 3.6.0) using stereo WorldView data; (2) physics-based radiative transfer equations (PBSDB) using Landsat data; and (3) a modified composite band-ratio method for Sentinel-2 (SatBathy) with an initial simplified calibration, followed by a more rigorous linear regression against in situ bathymetry data. All methods were tested in three different areas with different geological and environmental conditions, Cabo Rojo, Puerto Rico; Key West, Florida; and Cocos Lagoon and Achang Flat Reef Preserve, Guam. It is demonstrated that all satellite derived bathymetry (SDB) methods have increased accuracy when the results are aligned with higher-accuracy ICESat-2 ATL24 track bathymetry data using the iterative closest point (ICP). SDB vertical accuracy depends more on location characteristics than the method or optical satellite instrument used. All error metrics considered (mean absolute error, median absolute deviation, and root mean square error) can be less than 5% of the maximum bathymetry depth penetration for at least one method, although not necessarily for the same method for all sites. The SDB error distribution tends to be bimodal irrespective of method, satellite instrument, alignment, site, or maximum bathymetry depth, leading to the potential ineffectiveness of traditional error metrics, such as the root mean square error. However, our analysis demonstrates that performing detrending where possible can achieve an error distribution as close to normality as possible for which error metrics are more diagnostic.

Florida↗
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↗
USGS2026

Roadway runoff induced acute mortality in juvenile coho salmon during spring storm events

Extensive mortalities of adult coho salmon ( Oncorhynchus kisutch ), often called “Urban Runoff Mortality Syndrome” (URMS), have been documented during the fall in creeks where water quality has been degraded by roadway runoff. The primary cause of mortality is 6PPD-quinone (6PPDQ; N -(1,3-dimethylbutyl)- N ′-phenyl- p -phenylenediamine quinone)–an ozone transformation product that forms on all vehicle tires. Laboratory studies have shown that juvenile coho salmon are highly sensitive to 6PPDQ exposure. Unlike adults, juveniles reside in impacted watersheds year-round, including during the spring when 6PPDQ concentrations can frequently exceed lethal thresholds during storms. To assess the potential incidence of URMS in springtime rearing habitats for juvenile coho salmon, we conducted a paired water quality and toxicology study at Miller Creek, a runoff-impacted watershed in Normandy Park, WA, USA. Using a small field facility, three naïve groups of juvenile coho salmon ( N = 720) were exposed to either creek water or groundwater ( N = 120 per treatment per storm), across three spring storms while comparing water quality and mortality end points. In creek water during exposures, peak 6PPDQ concentrations reached 73–110 ng/L, exceeding reported median lethal concentrations (LC 50 ) for coho salmon. Over each 24–73 h storm exposure period, ∼80% of Miller Creek-exposed juvenile salmon died. No mortality occurred among control fish exposed to groundwater. These results indicate previously unidentified mortality risks for juvenile life stages of coho salmon during spring storms, suggesting substantial and year-round water quality impediments to coho salmon health and recovery across roadway runoff-impacted spawning, rearing, and migratory habitats.

Washington↗
USGS2025

Formation of the Mount Weld rare earth deposit, Western Australia: Geochronology constraints

Constraining the age of protracted chemical weathering in stable cratonic areas that may form thick regoliths and the potential enrichment of various elements is challenging. Economic deposits of aluminium, iron, copper, nickel, cobalt, niobium, and rare earth elements (REEs) form in this manner. Determining the age of formation can provide key information for exploration of similar deposits as well as to better constrain paleoclimatic conditions. This study describes our effort to constrain the age of formation of the Mount Weld deposit, a world-class carbonatite-derived REE laterite deposit. We utilize multiple geochronological techniques on different minerals. The oldest dates, ranging from ca. 100 to 50 Ma, were from laser ablation ICPMS, Lu-Hf dating of churchite, a heavy REE phosphate mineral formed by mineral saturation in groundwater. Growth bands on individual minerals show a younging outwards. 40 Ar/ 39 Ar geochronology of cryptomelane, a potassium-bearing manganese oxide mineral, yielded ages ranging from ca. 40 to 26 Ma. Similarly, (U–Th)/He geochronology of goethite yielded ages from ca. 45 to 19 Ma.

Mount Weld deposit↗
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↗
USGS2026

Estimating aftershock risk for entry into earthquake-damaged buildings

We present a simple method to estimate the risk of experiencing strong shaking from aftershocks during entry into earthquake-damaged buildings. We compute wait times until the probability of strong ground shaking from aftershocks reaches a predefined risk threshold; for example, a 0.4 percent probability of experiencing Modified Mercalli Intensity 7 or greater shaking during the planned building entry. We also develop a relation between aftershock probability and the probability of strong shaking, so that users can reference the U.S. Geological Survey aftershock forecast during an ongoing aftershock sequence to determine if the risk threshold has been met. We apply our method to active continental regions (for example, the Western United States), stable continental regions (for example, the Central and Eastern United States), and subduction zones (for example, Cascadia or Alaska).

Open-File Report↗
USGS2025

Aftershocks in stress shadows are inconsistent with modeled static Coulomb stress changes

Aftershock triggering is commonly attributed to increases in static Coulomb stress. In some areas, termed "stress shadows", a decrease in Coulomb stress is predicted to suppress earthquake occurrence. However, aftershocks are often observed in the modeled stress shadows. We examine several hypotheses that attempt to reconcile these shadow aftershocks with the static Coulomb stress change model: (1) they appear to be in shadows because of inaccuracy in the stress change calculations, (2) they occur on faults of unusual orientation which actually experienced increased Coulomb stress, (3) they occur on faults with different frictional properties, not modeled well by Coulomb stress, and (4) they are secondary aftershocks triggered by prior aftershocks or afterslip. When tested on the 2016 Mw7.0 Kumamoto, Japan, and 2019 Mw7.1 Ridgecrest, California, aftershock sequences, none of these hypotheses can explain the majority of the shadow aftershocks, and taken together these hypotheses can explain only about half of these aftershocks. This implies that Coulomb stress modeling that lacks small-scale fault zone heterogeneity might be inadequate to fully capture the true static stress changes and/or that other physical triggering models are needed, for example transient processes such as delayed triggering by dynamic stress changes from the passing seismic waves.

California↗
USGS2024

Distinguishing natural sources from anthropogenic events in seismic data

As seismic data are increasingly used to investigate a diverse range of subsurface phenomena beyond regular fast-rupturing earthquakes (Peng and Gomberg, 2010; Beroza and Ide, 2011), it is important to acknowledge that human-generated ground vibrations may be mistaken for naturally generated subsurface processes (Larose et al., 2015; Li et al., 2018). Correct discrimination of natural processes from anthropogenic noise is especially pressing given the trend in seismic detection research toward automated algorithms and machine learning methods (Yoon et al., 2015; Kong et al., 2019;Mousavi and Beroza, 2022) and the growth in seismic data collection in new environments such as urban and industry settings (e.g., Díaz et al.,2017).

Seismological Research Letters↗
USGS2026

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

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

Seismological Research Letters↗
USGS1976

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↗
USGS2024

Uncertainty in ground-motion-to-intensity conversions significantly affects earthquake early warning alert regions

We examine how the choice of ground‐motion‐to‐intensity conversion equations (GMICEs) in earthquake early warning (EEW) systems affects resulting alert regions. We find that existing GMICEs can underestimate observed shaking at short rupture distances or overestimate the extent of low‐intensity shaking. Updated GMICEs that remove these biases would improve the accuracy of alert regions for the ShakeAlert EEW system for the West Coast of the United States. ShakeAlert uses ground‐motion prediction equations (GMPEs), which calculate spatial distributions of peak ground acceleration (PGA) and peak ground velocity (PGV) from earthquake source estimates, combined with GMICEs to translate GMPE output into modified Mercalli intensity (MMI). We find significant epistemic uncertainty in alert distances; near‐source MMI estimates from different GMICEs can differ by over 1 MMI unit, and MMI extents used for public EEW alerts can differ by hundreds of kilometers for larger magnitude earthquakes ( M ∼6.5+). We use a catalog of “Did You Feel It?” shaking reports to evaluate how well GMICEs predict observed shaking. Our preferred GMICE is the one that computes MMI using PGV for high intensities and transitions to using PGA for nondamaging intensities. These results motivate updating GMICE relationships more generally, including in ShakeMap applications.

The Seismic Record↗
USGS2026

Magnitude conversion relations create substantial differences in seismic hazard models

Earthquake catalogs are essential data inputs for seismic hazard modeling. Because earthquake magnitudes are reported in a variety of types (e.g., local magnitudes and moment magnitudes), magnitude conversion relationships must be used to convert the different magnitude types present in a catalog to a uniform magnitude type to avoid biases in the hazard computation. However, these conversion relationships are often uncertain and have been shown to sometimes perform poorly. Here, we investigate the sensitivity of the gridded seismicity component of the National Seismic Hazard Model (NSHM) to the catalog conversion equations in the Eastern United States. In the 2023 NSHM, magnitudes of various types were converted to moment magnitudes using equations developed by the Central and Eastern United States Seismic Source Characterization for Nuclear Facilities (CEUS‐SSCn), based on least‐squares (LS) regressions made using data from a catalog containing events up through 2008. We recompute these equations using events in the Advanced National Seismic System Comprehensive Earthquake Catalog with multiple magnitudes from 2000 to 2023. Although we prefer the use of orthogonal regressions for our datasets, LS regressions produce broadly similar results, with both approaches exhibiting large deviations from the CEUS‐SSCn conversions, especially at smaller magnitudes. We compare the spatial distribution of annual rates using three different models: (1) the 2023 NSHM conversions, (2) our updated conversions, and (3) no conversions. We find that the choice of conversions leads to substantial differences in the rate forecasts, which can greatly impact the seismic hazard model, particularly in regions with low‐seismicity rates such as the Eastern United States, where the hazard is dominated by gridded seismicity rather than a fault model.

Seismological Research Letters↗
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↗
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↗
USGS2026

Surface variable‐based machine learning for scalable arsenic prediction in undersampled areas

In the United States, private wells are not federally regulated, and many households do not test for Arsenic (As). Chronic exposure is linked with multiple health outcomes, and risk can change sharply over short distances and with well depth. Coarse maps or sparse sampling often miss exceedances. Most existing models operate at ∼1 km resolution and use groundwater chemistry or detailed geologic logs, which limits their use in undersampled areas where improved guidance is most needed. We overcome these limitations by developing a machine learning model for Minnesota, USA, that predicts As exposure risk using only surficial variables from remote sensing and global data sets. Variables related to surface water hydrology and geomorphology are selected based on mechanistic links that control redox conditions and As mobilization. Local training was essential, and surficial geology variables that are more sensitive to local conditions were needed to maximize model accuracy. The resulting complete model was sufficiently sensitive to generate accurate and detailed risk maps and depth profiles of As concentrations above the 10 μg/L maximum contaminant level. Accuracy depended on local training data density. We identified a training data density of 0.07 wells/km 2 as a practical target for stable county-level performance. Maps of exceedance probabilities highlight priority areas for testing that are particularly important in rural communities that have received less sampling. These results support public health action by guiding where to install wells and where to test them, how much new sampling is needed, and where treatment outreach is most urgent.

Minnesota↗