USGS Science⌕ Search

SEARCH · USGS Science

Results for “Canadian Journal of Fisheries and Aquatic Sciences”

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,688 records · Page 19Linked to original sources
USGS2026

Evaluating groundwater quality influences from oil field operations and other anthropogenic activities in an urban setting, Santa Fe Springs, California

Groundwater quality is often affected by anthropogenic activities in urban settings. This study examines groundwater quality in and around the Santa Fe Springs Oil Field in Los Angeles County, California, where oil and gas production commonly intersects with high density industrial, commercial and residential land uses. Utilizing a combination of new and historical data, we evaluated potential pathways that would allow for oil field formation fluids to migrate into groundwater and whether mixing may have occurred based on the distribution of groundwater and oil field formation fluid tracers in samples. Samples were analyzed for a wide array of constituents including volatile organic compounds, light hydrocarbons, major ions, and various isotopic compositions. Despite evidence of oil field infrastructure providing potential pathways of migration via uncemented annular spaces, casing breaches and historical disposal of oil field formation water in surface ponds, the distribution and occurrence of stable isotopes of water, chloride, boron, and total dissolved solids do not indicate mixing of oil field formation water and groundwater. However, methane isotopic signatures and the presence of heavier alkanes suggest gas from oil-bearing formations have migrated from depth via oil field well infrastructure. Volatile organic compound detections were mainly from manufactured compounds unrelated to oil and gas production, with a relatively limited number of petroleum hydrocarbons also detected. Volatile organic compounds were generally found in wells tapping shallow, modern aged groundwater, indicating anthropogenic activities occurring at or near land surface as the source. Study results suggest that while oil field infrastructure provides migration pathways for oil field formation fluids to be introduced into groundwater, urban land uses not related to oil and gas production are the primary drivers of groundwater quality degradation.

California↗
USGS2025

Assessing nonpoint-source uranium pollution in an irrigated stream-aquifer system

Uranium (U) in rocks and soils of arid and semi-arid environments can be mobilized by irrigation and fertilization, posing environmental and health risks. Elevated U, along with selenium (Se) and nitrate (NO 3 ) co-constituents, necessitates careful monitoring and management. We developed a distributed-parameter numerical model to assess U pollution in an irrigated stream-aquifer system, applying it to a 552 km 2 region in Colorado's Lower Arkansas River Valley (LARV) over 14 years. A MODFLOW model, describing groundwater and stream flow, was coupled with an RT3D-OTIS model to portray reactive U transport. Calibration using the PESTPP-iES iterative ensemble smoother (iES) software indicated good agreement with observed U concentrations. The model revealed substantial and variable U levels across the LARV, highlighting potential hotspots and possible contributing factors, such as geological composition of the bedrock and near-surface shale and aquifer sediments derived from them, irrigation practices, and riparian landscape. U levels exceed the chronic standard (85th percentile = 30 μg/L, set by the US Environmental Protection Agency), which is the permissible regulatory threshold, in groundwater across 44 % of the region and along the river by an average factor of 2.9. Simulated average U concentrations in the non-riparian aquifer and river are 124 μg/L and 60 μg/L, respectively, compared with 112 μg/L and 62 μg/L for measured values. The average 85th percentile U concentration is 222 μg/L in the aquifer and 82 μg/L in the river. Average simulated U mass loading to the river is 0.17 kg/day per km, compared to an estimated 0.23 kg/day per km. Findings provide a baseline for comparing future simulated outcomes of alternative best management practices (BMPs) for U pollution mitigation and offer a methodology applicable to other irrigated regions.

Colorado↗
USGS2024

Produced water geochemistry from hydraulically stimulated Niobrara Formation petroleum wells: Origin of salinity and temporal perspectives on treatment and reuse

Produced water (i.e., a mixture of returned injection fluids and geologic formation brines) represents the largest volumetric waste stream associated with petroleum production in the United States. As such, produced water has been the focus of intense study with emphasis on understanding the geologic origin of the fluids, environmental impacts of unintended or intentional release, disposal concerns, and their commodity (e.g., lithium) potential. However, produced water geochemistry from many active petroleum plays remain poorly constrained leading to knowledge gaps associated with the origin of brine salinity and parameters (e.g., radium levels) that can impact treatment, disposal, and possible reuse. Here we evaluate the major ion geochemistry, radium concentrations, and stable water isotope composition of ~120 produced water samples collected from 17 producing unconventional petroleum wells in Weld County, Colorado from the Late Cretaceous Niobrara Formation. This sample set encompasses eight produced water time series from four new wells across production days 0 to ~365 and from four established wells across production days ~1000 to ~1700. Additionally, produced water from nine other established Niobrara Formation wells were sampled at discrete time points ranging from day 458 to day 2256, as well as hydraulic fracturing input fluids. These results expand the available Niobrara Formation produced water geochemical data, previously limited to few wells sampled within the first year of production, allowing for the heterogeneity of major ions and radium to be evaluated. Furthermore, we explore the geochemical relationships between major ion ratios and stable water isotope composition to understand the origin of salinity in Niobrara Formation brines from the Denver-Julesburg Basin. These findings are discussed with perspective toward potential treatment and reuse of Niobrara produced water prior to disposal.

Colorado, Nebraska, Wyoming↗
USGS2024

Global variability of the composition and temperature at the 410-km discontinuity from receiver function analysis of dense arrays

Seismic boundaries caused by phase transitions between olivine polymorphs in Earth's mantle provide thermal and compositional markers that inform mantle dynamics. Seismic studies of the mantle transition zone often use either global averaging with sparse arrays or regional sampling from a single dense array. The intermediate approach of this study utilizes many densely spaced seismic arrays distributed around the globe. We systematically compute teleseismic P-to-S receiver functions for each seismic array and invert for the 1-D seismic velocity structure of the mantle transition zone beneath each array to facilitate a comparison between densely sampled regions. We stack 3,600 receiver functions on average at 67 arrays in total. The stack is used in a probabilistic inversion to estimate the mantle transition zone interface depths and velocities beneath each array. We focus on the 410-km discontinuity (410) because it is a prominent seismic interface that is clearly linked to a single mineral phase transition between olivine and wadsleyite. The depths and velocity contrasts of the 410 are mapped to temperatures and compositions using mineral physics constraints. The depth of the 410 ranges from ∼405–440 km, which is consistent with a ∼360 K temperature range in a dry mantle and a ∼260 K temperature range in a wet mantle (2 wt. % water). The Vs contrast across the 410 ranges from ∼2.5–8 %, which is consistent with ∼20–70 vol. % olivine composition in a dry mantle and ∼25–80 vol. % in a wet mantle. The bulk composition of the upper mantle near the 410-km discontinuity is typically considered to be well-mixed because there is no thermodynamic impediment to convection at the olivine to wadsleyite phase transition. However, the wide range of inferred olivine content from our study suggests that there are large lateral variations in the bulk composition of the upper mantle near the 410-km discontinuity.

Earth and Planetary Science Letters↗
USGS2025

System characterization report on the Environmental Mapping and Analysis Program (EnMAP)

This report addresses system characterization of the Environmental Mapping and Analysis Program hyperspectral sensor by the DLR (German Aerospace Center, ground segment project management), GFZ (Deutsches Geoforschungszentrum, science lead) and is part of a series of system characterization reports produced and delivered by the U.S. Geological Survey Earth Resources Observation and Science Cal/Val Center of Excellence. These reports present and detail the methodology and procedures for characterization; present technical and operational information about the EnMAP hyperspectral sensor; and provide a summary of test measurements, data retention practices, data analysis results, and conclusions. The Earth Resources Observation and Science Cal/Val Center of Excellence system characterization team completed data analyses to characterize the geometric (interior and exterior), and radiometric performances of the EnMAP hyperspectral sensor. Results of these analyses indicate that the Environmental Mapping and Analysis Program has a band-to-band geometric performance in the range of −0.135 to 0.15 pixel, geometric performance relative to the Operational Land Imager in the range of −27.716 meters (−0.92 pixel) to 32.892 meters (1.09 pixels) offset in comparison to Landsat 8 Operational Land Imager, offset of a radiometric comparison in the range of −0.012 to 0.020, slope of a radiometric comparison in the range of 0.947 to 1.031.

Open-File Report↗
USGS2024

Potentially toxic elements in wild Agassiz’s desert tortoises: Tissue concentrations and association with disease

Background: Desert tortoise ( Gopherus agassizii ) populations have continued to decline due to infectious and other diseases, predation, and habitat alteration. The potential contribution of minerals and heavy metals to tortoise health and susceptibility to disease remains uncertain. Objective: The objective of this study was to evaluate the results of elemental analysis of trace minerals and macrominerals in scute keratin, kidney, and liver tissue from ill and dying desert tortoises salvaged for necropsy between 1993 and 2000. Methods: Salvaged tortoises were categorized by age (adult, juvenile), geographic location, and primary disease based on necropsy findings. A subset of tortoises that were injured or killed by vehicular trauma or predation but with no notable pathologic abnormalities were used as controls. A panel of 21 trace minerals and 6 macrominerals was analyzed in scute keratin, kidney, and liver tissue samples by inductively-coupled plasma spectrometry and atomic absorption spectrophotometry. Results: Necropsies were done on 46 tortoises, including 9 juveniles salvaged from 5 regions in the Colorado and Mojave deserts of California. Primary diseases were cutaneous dyskeratosis (n=9), infection/ inflammation (n=8), malnutrition (n=7), mycoplasmosis (n=5), and urolithiasis (n=3); 14 tortoises were classified as controls. Concentrations of elements differed significantly by tissue, age, desert region, and disease (P < 0.05). Tortoises with cutaneous dyskeratosis had significantly higher Se concentrations, primarily in keratin and liver, than tortoises with other diseases (P < 0.001). Juveniles were more likely than adults to have high Pb, Sn, and Zn levels (P < 0.05). All tortoises had detectable levels of more than one potentially toxic heavy metal, including As, Cd, Cr, Hg, Ni, Pb, Sn and V. Conclusions: Potentially toxic elements are frequently found in tissue from tortoises in desert regions of California, with significantly higher concentrations in diseased tortoises. Metal exposure from soils, mining, historic and ongoing military activities, and other human activities could increase susceptibility to disease in desert tortoises.

California↗
USGS2026

Importance of fish in the diet of Adélie penguins across multiple life stages

Over recent decades, the Adélie penguin ( Pygoscelis adeliae ) population has grown across most of its range, the exception being the northern coast of the western Antarctic Peninsula (WAP). In the Ross Sea, the very large Cape Crozier colony grew to reach ~9% of the global population. Demographic factors driving the Cape Crozier trend have yet to be identified. However, low chick fledging mass - a negative influence in the northern WAP where the diet is mostly less-energy-dense krill - is not showing an effect in the Cape Crozier colony. Previously, we hypothesized that Cape Crozier fledglings, and post-breeding adults, must be finding sufficient higher-quality prey once free of the prey-depleted colony foraging area, with subsequent higher survival. Here we explore penguin diet within and outside the colony foraging area using stable isotope analysis (SIA) of feathers: those grown by post-breeding adults outside the foraging area, prior to moult, compared to those of near-fledged chicks raised on meals obtained within the foraging area. Second, we explore whether a major difference occurs between breeding and post-breeding diet of adults from a small, nearby colony (Cape Royds) that exploits a smaller, little-depleted preyscape. Finally, we compare the results to SIA analyses published by others to explore whether the pre-moult diet of Ross Island adults, mostly foraging within the north-eastern Ross Sea, compares with that of adults nesting in the southern WAP (Bellingshausen Sea coast). At the latter, populations are not decreasing, with diet containing an appreciable contribution of fish. The results confirmed that fish contributed greatly to chick diets at both Ross Island colonies (δ 15 N 10.6–11.4). Once adults were released from central-place foraging, SIA levels (δ 15 N 9.7–12.9) were similar to those of the southern WAP. Foraging on energy-dense fish in coastal Antarctic waters could be one reason Adélie penguins are among the most abundant and therefore most ecologically successful penguin species.

Antarctic Science↗
USGS2024

The Europa Imaging System (EIS) investigation

The Europa Imaging System (EIS) consists of a Narrow-Angle Camera (NAC) and a Wide-Angle Camera (WAC) that are designed to work together to address high-priority science objectives regarding Europa’s geology, composition, and the nature of its ice shell. EIS accommodates variable geometry and illumination during rapid, low-altitude flybys with both framing and pushbroom imaging capability using rapid-readout, 8-megapixel (4k × 2k) detectors. Color observations are acquired using pushbroom imaging with up to six broadband filters. The data processing units (DPUs) perform digital time delay integration (TDI) to enhance signal-to-noise ratios and use readout strategies to measure and correct spacecraft jitter. The NAC has a 2.3° × 1.2° field of view (FOV) with a 10-μrad instantaneous FOV (IFOV), thus achieving 0.5-m pixel scale over a swath that is 2 km wide and several km long from a range of 50 km. The NAC is mounted on a 2-axis gimbal, ±30° cross- and along-track, that enables independent targeting and near-global (≥90%) mapping of Europa at ≤100-m pixel scale (to date, only ∼15% of Europa has been imaged at ≤900 m/pixel), as well as stereo imaging from as close as 50-km altitude to generate digital terrain models (DTMs) with ≤4-m ground sample distance (GSD) and ≤0.5-m vertical precision. The NAC will also perform observations at long range to search for potential erupting plumes, achieving 10-km pixel scale at a distance of one million kilometers. The WAC has a 48° × 24° FOV with a 218-μrad IFOV, achieving 11-m pixel scale at the center of a 44-km-wide swath from a range of 50 km, and generating DTMs with 32-m GSD and ≤4-m vertical precision. The WAC is designed to acquire three-line pushbroom stereo and color swaths along flyby ground-tracks.

Space Science Reviews↗
USGS2025

Bathymetric contour maps, surface area and capacity tables, and bathymetric change maps for selected water-supply lakes in Missouri, 2022–23

Bathymetric data were collected at 13 water-supply lakes around the periphery of Missouri by the U.S. Geological Survey in cooperation with the Missouri Department of Natural Resources and various local agencies, as part of a multiyear effort to establish or update the surface area and capacity tables for the surveyed lakes. Surveys were carried out during the months of April and May in 2022 and 2023. All but two of the lakes had been surveyed previously by the U.S. Geological Survey, and the recent surveys were compared to the earlier surveys to document the changes in the bathymetric surface and capacity of the lake. Bathymetric data were collected using a high-resolution multibeam mapping system mounted on a boat. Supplemental depth data at three of the lakes were collected in shallow areas with an acoustic Doppler current profiler on a remote-controlled boat. Data points from the various sources were exported at a gridded data resolution appropriate to each lake, either 0.82 foot, 1.64 feet, or 3.28 feet. Data outside the multibeam survey extent and greater than the surveyed water-surface elevation were obtained from data collected using aerial light detection and ranging (lidar) point cloud data. A linear enforcement technique was used to add points to the dataset in areas of sparse data (the upper ends of coves where the water was shallow or aquatic vegetation precluded data acquisition) based on surrounding multibeam and upland data values. The various point datasets were used to produce a three-dimensional triangulated irregular network surface of the lake-bottom elevations for each lake. A surface area and capacity table was produced from the three-dimensional surface for each lake showing surface area and capacity at specified lake water-surface elevations. Various quality-assurance tests were conducted to ensure quality data were collected with the multibeam, including beam angle checks and patch tests. Additional quality-assurance tests were conducted on the gridded bathymetric data from the survey, the bathymetric surface created from the gridded data, and the contours created from the bathymetric survey.

Missouri↗
USGS2025

Groundwater flowpath characteristics drive variability in per- and polyfluoroalkyl substances (PFAS) loading across a stream-wetland system

Groundwater-dependent ecosystems in areas with industrial land use are at risk of exposure to a PFAS chemicals. We investigated one such system with several known PFAS source areas, where high and low permeability sediments (glacial) coupled with groundwater-lake and groundwater/surface-water interactions created complex ‘source to seep’ dynamics. Using heat-tracing and chemical methods, numerous preferential groundwater discharge zones were identified and sampled across the upper Quashnet River stream-wetland system in Mashpee, MA, USA, downgradient of Joint Base Cape Cod (JBCC). Surface-water and groundwater samples were analyzed for 40 PFAS compounds between March and October 2022. Samples were collected from groundwater seeps identified as preferential discharge points (PDPs), wells upgradient of the stream-wetland system, contributing flow-through kettle lakes, and along Quashnet River surface-waters. PFAS from sampled waters contained perfluorinated carboxylates (PFPea, PFHxA, PFNA), perfluorinated sulfonates (PFBS, PFPeS, PFHxS, PFOS), fluorotelomer sulfonates (6:2, 8:2 FtS), and perfluoroalkyl sulfonamides (PFOSA). Samples from PDPs and wells had measured PFAS concentrations ranging from non-detect to 4677 ng/L ng/L (mean = 418 ng/L, std. = 709 ng/L), and a range of deuterium excess values (3.2 to 15.9 per mil) indicative of varying degrees of groundwater-lake interaction prior to groundwater flowpath emergence at PDPs. Correlations ( p < 0.01) between deuterium excess, %PFAS precursors, and terminal PFAS compounds highlighted potential precursor transformations associated with lake-groundwater exchange along flowpaths sourcing PDPs. However, some seepages had higher total PFAS concentrations (>1000 ng/L) than upgradient kettle lakes despite showing lake (evaporative) isotopic signatures, indicating the potential for groundwater flowpath convergence at wetland discharge zones and the influence of lakebed PFAS precursor reactions. Results from these synoptic surveys address gaps in the existing PFAS literature by demonstrating the importance of subsurface fate and transport on PFAS compound concentrations and mass loading in preferential groundwater discharge zones.

Massachusetts↗
USGS2026

Repeated coral bleaching events have eroded reef accretion potential in the Florida Keys

Coral bleaching is now the primary driver of coral-reef degradation globally, threatening the persistence of coral-reef functions and the invaluable ecosystem services they support. The Florida Keys have been impacted by four severe, regional-scale coral-bleaching events, which occurred in 1997–1998, 2005, 2014–2015, and 2023. Our aim in this study was to quantify the effects of those thermal-stress events on the persistence of coral-reef structures. We leveraged long-term coral-reef monitoring data to quantify the impacts of coral bleaching on coral cover, carbonate production, and reef-accretion potential throughout the Florida Keys from 1996–2024. Our study indicates that coral bleaching drove significant declines of reef-building corals, which successively diminished regional-scale reef-accretion potential (from 1.65 ± 0.33 mm y -1 in 1996) and culminated with a regional-scale transition to net reef erosion (-0.12 ± 0.08 mm y -1 ) following the 2023 event. Coral bleaching was the primary driver of declining reef-building capacity, with the cumulative impacts of the four coral-bleaching events accounting for >75% of total declines since 1996. The 1997–1998 event alone drove ~50% of overall declines, largely due to losses of the shallow-water ecosystem engineer, Acropora palmata . More moderate impacts of thermal stress on regional reef-accretion potential in 2005 (6%) and 2014–2015 (9%) could lend some support to the hypothesis that the remaining corals were more thermally tolerant; however, the unprecedented thermal stress event in 2023 caused widespread coral-bleaching and bleaching-related mortality, and caused an additional 12% of the total decline in reef-accretion potential. Stony-coral tissue loss disease accounted for an additional 6% of the decline, with the remaining 19% unexplained. Coral bleaching has already dramatically reshaped the structure and function of coral reefs in the Florida Keys over the last three decades. With the high likelihood that thermal stress will accelerate moving forward, the persistence of Florida’s essential coral-reef functions and services is uncertain.

Florida↗
USGS2025

Submarine groundwater discharge creates cold‐water refugia that can mitigate exposure of heat stress in nearshore corals

Coral reef mortality around the world is accelerating due to human activities and rising sea temperatures that cause bleaching, which is expected to become more frequent. Our ability to predict which corals will be most resilient, however, remains limited due to insufficient information characterizing nearshore temperature and habitat conditions. In this study, we examine how submarine groundwater discharge (SGD) reduces nearshore water temperatures and exposure of corals to heat stress, complementing the understanding that SGD can adversely affect coral when it contains elevated nutrient concentrations. Data from fixed nearshore sensors and vertical depth profiles along ~100 km of the western shoreline of the Island of Hawai’i from 2003 to 2014 demonstrate that submarine groundwater discharge (SGD) can reduce nearshore water temperatures by 1 °C–5°C and create estuarine-like conditions with salinities as low as 20 PSU, where the prevalent coral species, Pocillopora meandrina, Porites lobata , and Montipora capitata , thrive. Time-series temperature records reveal that exposure to high ambient ocean temperatures, which are known to initiate bleaching events, are reduced up to 5%–46% of the time. Coral health surveys indicated coral bleaching in response to moderately high annual temperatures in 2010 and 2011, with more colonies affected farther from cold, SGD-fed waters. Synthesis of these results, along with coral response data following the more extreme marine heat wave of 2014–2015, demonstrates lower coral loss and greater coral recovery near groundwater seeps, particularly those with higher flux and influence on reducing nearshore water temperatures. Our results demonstrate that SGD may therefore provide a beneficial ecosystem service and enhance coral reef resilience, particularly where human-related nutrient additions to groundwater can be mitigated. The implications of our findings are relevant across tropical coasts where groundwater inputs can be substantial, such as the Caribbean and Indo-Pacific, and contribute to improving our understanding of coral sensitivity to gradients in temperature and nutrient stress. Improved management of groundwater resources could thus be vital to local–regional strategies for mitigating future heat stress.

Hawaii↗
USGS2026

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

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

Colorado, New Mexico↗
USGS2024

A benchmark dataset and workflow for landslide susceptibility zonation

Landslide susceptibility shows the spatial likelihood of landslide occurrence in a specific geographical area and is a relevant tool for mitigating the impact of landslides worldwide. As such, it is the subject of countless scientific studies. Many methods exist for generating a susceptibility map, mostly falling under the definition of statistical or machine learning. These models try to solve a classification problem: given a collection of spatial variables, and their combination associated with landslide presence or absence, a model should be trained, tested to reproduce the target outcome, and eventually applied to unseen data. Contrary to many fields of science that use machine learning for specific tasks, no reference data exist to assess the performance of a given method for landslide susceptibility. Here, we propose a benchmark dataset consisting of 7360 slope units encompassing an area of about 4,100 km 2 "> 4,100 km 2 in Central Italy. Using the dataset, we tried to answer two open questions in landslide research: (1) what effect does the human variability have in creating susceptibility models; (2) how can we develop a reproducible workflow for allowing meaningful model comparisons within the landslide susceptibility research community. With these questions in mind, we released a preliminary version of the dataset, along with a “call for collaboration,” aimed at collecting different calculations using the proposed data, and leaving the freedom of implementation to the respondents. Contributions were different in many respects, including classification methods, use of predictors, implementation of training/validation, and performance assessment. That feedback suggested refining the initial dataset, and constraining the implementation workflow. This resulted in a final benchmark dataset and landslide susceptibility maps obtained with many classification methods. Values of area under the receiver operating characteristic curve obtained with the final benchmark dataset were rather similar, as an effect of constraints on training, cross–validation, and use of data. Brier score results show larger variability, instead, ascribed to different model predictive abilities. Correlation plots show similarities between results of different methods applied by the same group, ascribed to a residual implementation dependence. We stress that the experiment did not intend to select the “best” method but only to establish a first benchmark dataset and workflow, that may be useful as a standard reference for calculations by other scholars. The experiment, to our knowledge, is the first of its kind for landslide susceptibility modeling. The data and workflow presented here comparatively assess the performance of independent methods for landslide susceptibility and we suggest the benchmark approach as a best practice for quantitative research in geosciences.

Earth-Science Reviews↗
USGS2025

Groundwater structures fish growth and production across a riverscape

Landscapes are composed of habitat patches and conditions that vary across space and time. While habitat variability and complexity can support important ecological processes and ecosystem services, the dynamic nature of habitats can also constrain organismal growth and production as optimal conditions are fleeting. In riverine ecosystems, groundwater discharge to streams stabilises water temperature and flow regimes, thus mediating how habitat complexity is expressed. Yet, how stable habitats structure growth and production within the broader landscape matrix is not well understood. In this study, we explored the effects of groundwater on spatiotemporal variation in growth and production for juvenile Yellowstone cutthroat trout ( Oncorhynchus virginalis bouvieri ) across the upper Snake River catchment, Wyoming, USA. We combined machine learning techniques and remotely sensed landscape data to estimate groundwater availability across the river network, which we linked to stream temperature regimes and conspecific density. We then used Bayesian hierarchical models to quantify the effects of temperature, density and groundwater on spatiotemporal variation in fish growth and production in 52 focal reaches. Finally, we predicted body size trajectories and trends in total production continuously over both space and time to understand the effect of groundwater at the riverscape scale. Groundwater discharged to streams where topography changes abruptly in valley-bottom areas underlain by coarse glacial deposits. Groundwater stabilised temperature regimes and was associated with high trout densities. Temperature and density, in turn, interacted to influence growth rates: growth increased strongly with temperature, but this effect was reduced when density was high. Accordingly, variation in groundwater availability among stream reaches diversified growth and production regimes. In reaches with low groundwater availability, growth and production declined over time from summer maxima. In contrast, in reaches with high groundwater availability, temporal trends in growth and production were hump-shaped—peaking in autumn—and mean production was greater. At the riverscape scale, temporal asynchrony in growth rates generated convergent spatial variation in growth capacity, but—when combined with density—led to the formation of distinct hotspots of production. Our results demonstrate how groundwater, an important driver of aquatic ecosystem heterogeneity, structures trout growth and production across space and time. Importantly, rare, but stable habitats may disproportionately affect ecological processes and serve as key sources of population diversity at larger spatial scales.

Wyoming↗
USGS2025

Tapwater exposures, residential risk, and mitigation in a PFAS-impacted-groundwater community

Tapwater (TW) safety and sustainability are priorities in the United States. Per/polyfluoroalkyl substance(s) (PFAS) contamination is a growing public-health concern due to prolific use, widespread TW exposures, and mounting human-health concerns. Historically-rural, actively-urbanizing communities that rely on surficial-aquifer private wells incur elevated risks of unrecognized TW exposures, including PFAS, due to limited private-well monitoring and contaminant-source proliferation in urbanizing landscapes. Here, a broad-analytical-scope TW-assessment was conducted in a hydrologically-vulnerable, Mississippi River alluvial-island community, where PFAS contamination of the shallow-alluvial drinking-water aquifer has been documented, but more comprehensive contaminant characterization to inform decision-making is currently lacking. In 2021, we analyzed 510 organics, 34 inorganics, and 3 microbial groups in 11 residential and community locations to assess (1) TW risks beyond recognized PFAS issues, (2) day-to-day and year-to-year risk variability, and (3) suitability of the underlying sandstone aquifer as an alternative source to mitigate TW-PFAS exposures. Seventy-six organics and 25 inorganics were detected. Potential human-health risks of detected TW exposures were explored based on cumulative benchmark-based toxicity quotients ( ∑ TQ ). Elevated risks ( ∑ TQ ≥ 1) from organic and inorganic contaminants were observed in all alluvial-aquifer-sourced synoptic samples but not in sandstone-aquifer-sourced samples. Repeated sampling at 3 sites over 52–55 h indicated limited variability in risk over the short-term. Comparable PFAS-specific ∑ TQ for spatial-synoptic, short-term (3 days) temporal, and long-term (3 years quarterly) temporal samples indicated that synoptic results provided useful insight into the risks of TW-PFAS exposures at French Island over the long-term. No PFAS detections in sandstone-aquifer-sourced samples over a 3 year period indicated no PFAS-associated risk and supported the sandstone aquifer as an alternative drinking-water source to mitigate community TW-PFAS exposures. This study illustrated the importance of expanded contaminant monitoring of private-well TW, beyond known concerns (in this case, PFAS), to reduce the risks of a range of unrecognized contaminant exposures.

Wisconsin↗
USGS2026

Inference of pattern-based geological CO2 sequestration and oil recovery potential in a commingled main pay and residual oil zone CO2-EOR flood

Several detailed studies have shown that residual oil zones (ROZs) can present significant resources for additional hydrocarbon recovery as well as subsurface carbon dioxide (CO 2 ) sequestration via enhanced oil recovery by injecting CO 2 (CO 2 -EOR). Field development strategies included new wells drilled dedicated to main pay zones (MPZ) and ROZs, or existing wells in MPZs deepened to ROZs for commingled injection-production using different well patterns. The latter presented a challenge when discerning the injection and production from each of the zones, and for subsequent quantification of CO 2 sequestration and EOR potential from different patterns and from the field. In this paper, an innovative method for analyzing commingled injections and productions from MPZs and ROZs, with application to pattern-based data from four staggered line drive patterns in Wasson Field's Denver Unit, Texas, USA, was developed. Decline curve and ratio-trend methods were used as means of history-matching and forecasting. Cumulative production-time and cumulative production-rate data for oil, gas, and water, as well as water-oil ratio (WOR) and gas-oil ratio (GOR), were analyzed along with injection data for time intervals covering major injection events in MPZ, or MPZ and ROZ combined. A combined analysis enabled inference of allocation of fluids into different zones during WAG (water alternating gas) injection and thereby estimation of CO 2 storage, utilization, and retention in different zones as a function of total injection. Results show that ROZs generally present higher CO 2 sequestration potential compared to MPZs, and a comparable incremental oil recovery factor of ∼20%, on average. Results based on ratio analysis further show that while the WOR trend of the pattern production is mostly dominated and controlled by ROZ, GOR is controlled by both intervals. Although the method relying on decline curves and the approach used in zonal fluid allocations are subject to their limitations, this study presents a practical and innovative well-pattern-based method to infer and forecast CO 2 sequestration and oil recovery quantities and fluid ratios from MPZs and ROZs in commingled operations and highlight the added potential offered by ROZs.

Texas↗
USGS2024

Influence of inherited structure on flexural extension in foreland basin systems: Evidence from the northern Arkoma basin and southern Ozark dome, USA

Extensional faults are key components of foreland basin systems. They form within the upper crust in response to flexure of the lithosphere and accommodate subsidence within the foredeep and forebulge depozones. Such faults are excellent proxies for orogenic system evolution and control the distribution of natural resources and hazards. However, the spatiotemporal evolution of flexural extension has not been documented previously at a regional scale, thereby limiting our understanding of underlying geodynamic controls. Here, we resolve late Paleozoic flexural extension in the northern Arkoma basin and southern Ozark dome, USA. We synthesize a large database of previous mapping, existing research, subsurface data, and geophysical data into 3D geologic and 2D kinematic models. Mesh surfaces representing several key horizons from the Carboniferous Period (ca. 335-306 Ma) were constructed. These surfaces were built from oil and gas well tops (n = ∼10,000) and surface geologic map contacts using an advanced kriging method. The mesh surfaces are offset by a complex 3D fault network, allowing detailed analysis of along-strike and down-dip variations in fault displacement. Analysis of the 3D model reveals a regular and repeated fault segmentation pattern wherein E -W striking, left- and foreland-stepping en échelon normal faults are segmented by inherited NE striking basement faults. Maximum vertical separation along the E -W normal faults is generally focused between the inherited NE-trending faults. This suggests that the inherited basement faults delocalized extensional strain during late Paleozoic normal faulting. Maximum vertical separation and fault localization may correlate to areas with high-amplitude positive magnetic anomalies interpreted as Mesoproterozoic granitic rocks. Speculative covariance of magnetic anomalies and fault displacements implies that the relatively strong basement granite concentrated stress, leading to localized faulting within the relatively thin sedimentary cover. Lastly, we show that flexural extension migrated southeast to northwest from the Chesterian-Morrowan (ca. 335-319 Ma) to the Desmoinesian (ca. 306 Ma). The migratory flexural extension may be explained by diachronous loading during Pangean assembly, or by synchronous loading but variable load compensation due to inherent factors.

Arkansas, Oklahoma↗