SEARCH · USGS Science

Results for “Environmental Modelling & Software”

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

SourceResearch connectionFrom yearThrough year

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

633 records · Page 6Linked to original sources

USGS2025

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

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

Arkansas, Missouri, Oklahoma↗

USGS2024

SUTRA— A code for simulation of saturated-unsaturated, variable-density groundwater flow with solute or energy transport—Documentation of the version 4.0 enhancements—Freeze-thaw capability, saturation and relative-permeability relations, spatially varying properties, and enhanced budget and velocity outputs

Version 4.0 of the Saturated-Unsaturated Transport (SUTRA) software code provides the capability to simulate the freezing and thawing of groundwater during energy transport simulations under saturated and unsaturated conditions. In addition to the types of hydrogeologic processes that SUTRA has been able to simulate in the past, this version can be used to study the effects of the freeze-thaw process on the flow and energy dynamics of hydrogeologic systems. The freeze-thaw simulation capability accounts for the latent heat of fusion and allows thermal property values to vary with changing total-water saturation, liquid-water saturation, and ice saturation. It allows the effective permeability of the porous medium to change as a result of freezing and thawing. This version also provides several user-selectable relations for the dependence of total-water saturation on fluid pressure, the dependence of liquid-water saturation on temperature during freezing and thawing, and the dependence of relative permeability on liquid saturation, as well as three user-selectable formulae for defining the bulk thermal conductivity of a mixture of solid grains, liquid water, ice, and air. For unsaturated simulations without freezing, the selectable total-water saturation relations eliminate the need for the user to program these and their associated relative-permeability functions, as had been required in previous SUTRA versions. Optional nonlinear dependence of fluid density on temperature, which covers the range from supercooled (about −50 degrees Celsius) to superheated (about 400 degrees Celsius), is also provided. Additionally, this version makes it possible to spatially vary parameters that, in previous versions of SUTRA, were required to be spatially uniform: solid-matrix properties, adsorption parameters, and parameters for production of solute mass or energy. Spatial variation is also allowed for the newly included freeze-thaw process parameters. Additional enhancements provide (1) output of water-mass and energy budgets that include values of all component terms in the governing balance equations, and (2) output of Darcy velocities (fluid fluxes), in addition to the velocity output provided by previous SUTRA versions. These enhanced outputs allow fuller interpretation of simulation results, especially for freeze-thaw phenomena. The set of processes simulated by this version of SUTRA are useful for studying a wide range of hydrogeologic system types, conditions, and questions. For cryohydrogeologic simulations, however, this version of the code is limited in that (1) it does not simulate thermomechanical effects of freeze-thaw, (2) pressure changes due to water density change during freezing are neglected, (3) ice saturation cannot exceed the initial porosity of the simulated medium, and (4) cryosuction, the migration of liquid water toward freezing fronts, is neglected. Furthermore, this version does not account for air flow or for water vaporization and sublimation under unsaturated conditions.

Techniques and Methods↗

USGS2026

Assessing future hydrologic extremes using an integrated hydrology and river operations model in the Russian River watershed

Study region The Russian River watershed, situated in coastal, northern California, experiences hydrologic extremes, including periodic droughts and flooding. Water managers are working to maintain sustainable water supplies and environmental flows, while mitigating flood risks. Study focus This paper introduces an integrated hydrology and river operations model for the Russian River watershed. This model is distinct from models in previous studies because it represents surface-groundwater interactions and uses climate forcings to estimate dynamic water use demands that are superimposed onto both reservoir operations and water supply constraints. The model was used to examine three historical (1990–2015) and eight future (2016–2099) water use and climate change scenarios. New hydrological insights for the region The direct connection between streams and aquifers facilitated both annual aquifer replenishment by high winter streamflows and streamflow depletion by groundwater wells (19 % of pumped groundwater in alluvial aquifers from stream leakage) during critical low flow periods. Simulated streamflow changes included 59 % longer and 54 % more severe streamflow droughts, 26 % lower seasonal low streamflows, and up to 125 % higher peak streamflows, averaged over future climate and water use scenarios, suggesting increased future flood and water availability risks. Results showed the importance of reservoir operations for mitigating the impacts of increased hydroclimatic volatility, despite a decrease in reservoir reliability at Lake Mendocino, suggesting that reservoir management may be used to decrease future risks.

California↗

USGS2024

Parsimonious high-resolution landslide susceptibility modeling at continental scales

Landslide susceptibility maps are fundamental tools for risk reduction, but the coarse resolution of current continental-scale models is insufficient for local application. Complex relations between topographic and environmental attributes characterizing landslide susceptibility at local scales are not transferrable across areas without landslide data. Existing maps with multiple susceptibility classifications under-represent landslide potential in moderate and gently sloping terrain. We leverage an extensive landslide database ( N = 613,724), a high-resolution digital elevation model (10-m), and high-performance computing resources, to develop a new nationwide susceptibility map for the contiguous United States, Hawaii, Alaska, and Puerto Rico. We calculate four alternative linear and nonlinear thresholds of topographic slope and relief using an objective split-sample calibration. We down-sample our results to a 90-m grid to account for uncertainty in the digital elevation model and landslide position, and evaluate these thresholds' ability to differentiate areas of greater susceptibility. The less conservative nonlinear model optimally balances our priorities of capturing observed landslides (99%) while minimizing area covered by susceptible terrain (43%). Independent evaluation with four statewide landslide inventories ( N = 172,367) reinforces our model selection but highlights spatially variable performance. Therefore, we propose a novel approach to susceptibility classification using the concentration of landslide-prone terrain within each down-sampled grid. While landslides are possible within any cells containing susceptible terrain, those with the highest concentration capture the majority of observed landslides. Our new map characterizes landside susceptibility more consistently than prior models; our transparent classification approach also provides flexibility for accommodating different tolerances in risk reduction measures.

AGU Advances↗

USGS2025

Component identification of solid biomass fuels using reflected light microscopy: Interlaboratory study 2

As nations transition toward sustainable energy systems , biomass has become a vital component of global energy portfolios. Derived from organic materials such as wood, agricultural residues, forestry byproducts, and organic waste, biomass is a renewable energy source with significant environmental and economic benefits. Responsible biomass energy production can improve waste management, reduce emissions of greenhouse gases, and mitigate environmental pollution. However, as the diversity of biomass-derived fuels increases, robust quality assessment methods are essential to ensure their efficiency, safety, and minimal environmental impact. Reflected light microscopy (RLM) is one such technique with the potential to complement conventional physico-chemical analyses by enabling a rapid identification of material constituents and impurities. To refine this methodology and evaluate the reproducibility of solid biomass component identification using RLM, an interlaboratory study (ILS) was conducted. The study involved the recognition of 58 components across 45 photomicrographs, with the participation of 65 scientists and students from 25 countries. The participants faced high difficulty identifying some of the marked components, and as a result, the percentage of correct answers ranged from 19.0 % to 98.3 %, with an average correct identification rate of 62.7 %. The most challenging aspects of the identification process included distinguishing between woody and non-woody (agro) biomass, accurately identifying petroleum-derived materials, and differentiating agro biomass from inorganic matter. The results suggest that while RLM is an important tool for characterizing solid biomass, further development of methodology guidelines and training are necessary to enhance its effectiveness. Future research should prioritize preparing detailed, image-rich, microscopic morphological descriptions of biomass fuel components, which could improve the accuracy and reliability of using RLM in biomass fuel characterization.

International Journal of Coal Geology↗

USGS2026

Decoding the oxidative digestion mechanism for polystyrene nanoplastic detection in the Great Lakes using a customizable Raman spectral processing algorithm

Despite the concerns surging around nanoplastics (NPs) regarding their prevalence and bioavailability in freshwater systems, robust detection of NPs in complex environmental matrices is hindered by the lack of standardized sample pretreatment and a mechanistic understanding of oxidative digestion. Here, we systematically investigate the interaction between hydrogen peroxide (H 2 O 2 ) and polystyrene (PS) NPs during digestion in deionized (DI) water and four environmental matrices from in and around the Great Lakes Basin. To facilitate high-throughput analysis, we develop Pre_peak, a customizable Raman spectral processing algorithm that achieves >99% accuracy for both NP identification and interference rejection, allowing reliable NP quantification via pixel counting and systematic decoding of the oxidative digestion mechanisms. In DI water, varying H 2 O 2 doses from 0 to 30% has negligible effects on the recovery and Raman signal intensity of PS NPs over 24 hours of digestion. However, morphological changes and aggregation of PS NPs are observed when the H 2 O 2 dose exceeds 20%. Prolonged digestion further leads to progressive NP loss. In natural waters, the optimal dosage and digestion duration depend on matrix characteristics, including dissolved organic matter (DOM) and ion composition. This study provides mechanistic insights into NP–oxidant interactions and underscores the need for matrix-tailored digestion protocols to advance standardized NP detection in freshwater environments.

Great Lakes↗

USGS1998

Design, revisions, and considerations for continued use of a ground-water-flow model of the coastal plain aquifer system in Virginia

A digital numerical model of the ground-water-flow system in the Virginia Coastal Plain was constructed as part of the Regional Aquifer System Analysis program of the U.S. Geological Survey and subsequently revised based on reinterpretations of the hydrogeologic framework. The revised model has been incorporated by the Virginia Department of Environmental Quality as a means to evaluate the effects of existing and proposed ground-water withdrawals. The revised model, which is documented here, includes many design features of the original model. Principal modifications are (1) reassignment of some model layers to represent different aquifers, (2) replacement of applied recharge and the uppermost surface-water layer by constant-head cells to represent the unconfined aquifer, and (3) revision of the lateral extents and transmissivity and vertical leakance values assigned to model layers. Ground-water levels simulated by the revised model under steady-state conditions, using withdrawal rates representing the period 1978-80, generally are similar to the water levels simulated by the original model under transient conditions using the same withdrawal rates. Differences probably result largely from changes in transmissivity and vertical leakance values. T he revised model enables the continued evaluation of development of the large and complex aquifer system. Improved information on parts of some aquifers has been incorporated. Only large, regional-scale trends in ground-water water levels, flow directions, and flow rates are simulated, however, and not local trends and (or) short-term changes in water level and flow. Changing hydraulic stresses, knowledge of hydrogeologic conditions, and resource-management needs require periodic revision of the model to maintain its usefulness for the management of ground-water resources in the Virginia Coastal Plain.

Virginia↗

USGS2025

Shallow geologic framework of the Mississippi Sound and the potential for sediment resources

The Mississippi Sound, an estuarine environment located between the mainland and barrier islands bordering the northern Gulf of America (formerly the Gulf of Mexico), serves as a vital ecosystem for the States of Mississippi and Alabama. Spanning approximately 100 kilometers from east to west and covering 1,400 square kilometers, the sound is home to marine industry and ports, and its shallow and brackish waters sustain a diverse array of marine life. Barrier islands along the southern edge of the sound separate the microtidal estuary from the Gulf of America. This protection from gulf wave action mediates current flow within the sound, resulting in predominantly fine-grained sediment deposition along the seafloor. This study, conducted by the U.S. Geological Survey in cooperation with the U.S. Army Corps of Engineers, provides insight on fluvial and tidal processes spanning the past 5,000 years. The report synthesizes existing research to provide a comprehensive overview of the sound geology, from Pleistocene origins to present-day morphology, and utilizes high-resolution single channel seismic profiles and sediment data to identify and map sedimentary deposits and morphologic features at and below the seafloor. Despite its ecological significance, the Mississippi Sound faces environmental challenges, including water-quality issues, habitat degradation, storm-induced erosion, and the ongoing threats of sea-level rise and environmental changes. This study uses the present-day understanding of the sound's geology to inform coastal management decisions, hazard assessment, and potential mineral resources.

Louisiana, Mississippi↗

USGS2024

RegionGrow3D: A deterministic analysis for characterizing discrete three-dimensional landslide source areas on a regional scale

Regional-scale characterization of shallow landslide hazards is important for reducing their destructive impact on society. These hazards are commonly characterized by (a) their location and likelihood using susceptibility maps, (b) landslide size and frequency using geomorphic scaling laws, and (c) the magnitude of disturbance required to cause landslides using initiation thresholds. Typically, this is accomplished through the use of inventories documenting the locations and triggering conditions of previous landslides. In the absence of comprehensive landslide inventories, physics-based slope stability models can be used to estimate landslide initiation potential and provide plausible distributions of landslide characteristics for a range of environmental and forcing conditions. However, these models are sometimes limited in their ability to capture key mechanisms tied to discrete three-dimensional (3D) landslide mechanics while possessing the computational efficiency required for broad-scale application. In this study, the RegionGrow3D (RG3D) model is developed to broadly simulate the area, volume, and location of landslides on a regional scale (≥1,000 km 2 ) using 3D, limit-equilibrium (LE)-based slope stability modeling. Furthermore, RG3D is incorporated into a susceptibility framework that quantifies landsliding uncertainty using a distribution of soil shear strengths and their associated probabilities, back-calculated from inventoried landslides using 3D LE-based landslide forensics. This framework is used to evaluate the influence of uncertainty tied to shear strength, rainfall scenarios, and antecedent soil moisture on potential landsliding and rainfall thresholds over a large region of the Oregon Coast Range, USA.

Journal of Geophysical Research: Earth Surface↗

USGS2026

Regional models for postfire debris-flow likelihood and rainfall thresholds across the western United States

The U.S. Geological Survey (USGS) uses an empirical model developed with logistic regression (the ‘M1’ model) to rapidly assess debris-flow likelihood and to identify quantitative rainfall thresholds for debris flows after wildfire in the western United States. The M1 model was calibrated to a debris-flow inventory from southern California (United States) and has been applied throughout the western United States. Limited spatial coverage in the calibration dataset has motivated evaluation of M1 model accuracy outside the calibration region (e.g., the Sierra Nevada or the eastern Cascade Range, United States). Previous test cases showed that M1 overpredicts debris-flow likelihood and underpredicts rainfall thresholds for some locations (e.g., Arizona, northern California, Colorado, New Mexico, United States). We sought to improve the regional applicability of a debris-flow likelihood model by expanding the debris-flow inventory used for calibration, testing multiple potential models and generating an updated model framework. The updated inventory includes 3788 observations from 67 burned areas paired with short duration rainfall ratios. The updated model framework consists of a modified model structure and sets of coefficients calibrated separately to the entire updated inventory and to subsets of the inventory that intersect three Environmental Protection Agency (EPA) Level 2 ecoregions (Mediterranean California, Upper Gila Mountains and Western Cordillera). Comparisons of predictions from the updated models with observed rainfall and debris-flow activity show that the updated models outperform the M1 model by ~15%–60% and improve the uniformity of predictive performance across the western United States. The updated models also reduce false positive rates relative to M1 and generate rainfall thresholds that are better aligned with relative differences in regional climatology and debris-flow activity.

Arizona, California, Colorado, Idaho, Montana, Nev↗

USGS1999

The mineral economy of Brazil--Economia mineral do Brasil

This study depicts the Brazilian government structure, mineral legislation and investment policy, taxation, foreign investment policies, environmental laws and regulations, and conditions in which the mineral industry operates. The report underlines Brazil's large and diversified mineral endowment. A total of 37 mineral commodities, or groups of closely related commodities, is discussed. An overview of the geologic setting of the major mineral deposits is presented. This report is presented in English and Portuguese in pdf format.

Data Series↗

USGS2025

Hydrogeologic characterization of the Cahuilla Valley and Terwilliger Valley Groundwater Basins, Riverside County, California

The relation between the groundwater and the amount of natural recharge to the Cahuilla Valley and Terwilliger Valley groundwater basins is not well understood. During the 20th century, the reliance on groundwater near Anza, California, used for agricultural, domestic, and municipal reasons has increased, and there is the potential for changes in groundwater availability related to climate change. Several types of existing data were evaluated, and new data were collected for this study, with the goal of characterizing the region’s hydrogeology. The study’s scope included constructing a geologic framework model to show where the groundwater-bearing units are present and their relation to each other, estimating the major components of the groundwater budget, and understanding local short-term and regional long-term groundwater flow and how that has changed since the early 1900s. Two electrical resistivity tomography surveys were done in the Durasno Valley about 2,150 feet apart to identify the thickness of the alluvium, its horizontal extent, and the depth-to-basement along two profiles perpendicular to Cahuilla Creek. The subsurface sediments were mostly horizontally layered and the transitional boundary between the alluvium and basement was thinner and shallower along the upgradient profile where the depth-to-basement was about 70 feet below land surface; the depth-to-basement at the downgradient profile was more than about 140 feet below land surface. The results from the surveys were used to place four monitoring wells at two sites along the survey profiles. Artesian flow from the deepest well at the downgradient site indicated that the decomposed and competent basement likely contributed some groundwater to the overlying alluvium, laterally, from below, or both. A digital three-dimensional geologic framework model was constructed using EarthVision software to represent the subsurface geometry of the alluvium, decomposed basement, and competent basement. Maps and cross sections of the modeled thicknesses of the alluvium and decomposed basement, and the modeled elevation of the top of the competent basement, were made to show the subsurface geometry of vertical faults, selected wells, and the groundwater-bearing units. Because natural recharge is related to the variable cycles of precipitation, estimates are difficult to quantify. Recharge and runoff have extreme interannual variability in the study area; recharge and runoff can be sporadic, and a substantive amount may not occur in some years. Estimates of recharge from a previous study and the regional-scale Basin Characterization Model for California for four different periods ranged from 3,800 acre-feet/year for 1897–1947 to 5,900 acre-feet/year for 1971–2000. Potential recharge from the disposal of domestic septic systems may have been as much as 500 acre-feet in 2020. It was estimated that between about 400 and 2,400 acre-feet/year of groundwater is lost through evapotranspiration by vegetation and evaporation from open water bodies, but the main source of discharge is through pumpage, mainly used for agriculture from the alluvium in the Cahuilla Valley and Terwilliger Valley groundwater basins. The estimated total pumpage for 1991–2021 ranged from about 1,140 acre-feet in 2019 to about 3,450 acre-feet in 1994. When summed, the cumulative amount of estimated pumpage between 1991 and 2021 was about 81,400 acre-feet. The general direction of groundwater flow is from the northeast along the San Jacinto fault zone at the headwaters of Cahuilla and Hamilton Creeks, to the surface-water outlets at the west and southeast parts of the study area. Groundwater-level data from the 1950s and earlier indicate that there was a natural groundwater divide between the Cahuilla Valley and Terwilliger Valley groundwater basins, but the changing magnitude and extent of the groundwater depressions caused by pumping since about 1950 indicate that the location of the natural groundwater boundary between the Cahuilla Valley and Terwilliger Valley groundwater basins has migrated over time. Flow from the upper to the lower parts of the Cahuilla Valley groundwater basin roughly follows the course of Cahuilla Creek through the narrow Durasno Valley where an estimated volume of flow in April 2019 was about 10–150 acre-feet/year. Short-term trends in groundwater levels, particularly in wells where groundwater is shallow and in the basement unit, show how some areas respond quickly to recharge and discharge. Wells located further to the east within the Cahuilla Valley groundwater basin in the alluvium show much less of a response to recharge events; areas of sustained pumpage from the alluvium, primarily for agriculture, show long-term declines in groundwater levels and generally do not show the effects of storm events or recent runoff. Groundwater levels in wells that are farthest from where most of the recharge occurs and where pumping has been the greatest, had some of the largest long-term groundwater-level declines at a rate of about 0.8 foot/year between 1971 and 2021.

California↗

USGS2025

Scoping decision-maker needs and science availability to support regional natural capital accounting in the U.S. Colorado River Basin

Natural capital accounting has the potential to yield important policy insights at multiple scales, but there remains a disconnect between regional-scale natural capital accounts and their use for informing policy. In this paper, we propose a roadmap that could lead to the creation of policy-relevant regional accounts, with steps split across an initial scoping phase and a subsequent development phase. We demonstrate the scoping steps in action with an application to the Colorado River Basin (“Basin”), a large watershed in the southwestern United States (U.S.) that has faced aridification and substantial high-profile tradeoffs around the use of its water and other natural resources. Drawing on prior U.S. Geological Survey science co-production efforts, we conducted a series of eight discussion sessions with 41 scientists and science representatives whose work is relevant to Basin water, riparian and riverine ecosystems, upland ecosystems and energy and minerals. We summarise participants' thoughts on key topics and economic linkages, their insights and questions of interest and their recommendations on existing scientific data sources and gaps. We evaluate the suitability of the available data for construction of System of Environmental-Economic Accounting (SEEA) Central Framework and SEEA Ecosystem Accounting accounts, including those for land, water, forests, energy and minerals and ecosystems (covering extent, condition and ecosystem services). We present a series of lessons learned during the scoping phase, as well as lessons that could be relevant for future practitioners engaging in the development phase. The information can help guide the development of timely and relevant regional-scale environmental-economic accounts in the U.S. and beyond.

Arizona, California, Colorado, Nevada, New Mexico,↗

USGS2025

Natural source zone depletion of crude oil in the subsurface: Processes controlling mass losses of individual compounds

At many petroleum hydrocarbon spill sites, residual spilled product forms a long-term source of groundwater contamination. The phrase source zone natural depletion is used to refer to the mass loss rates. Overall mass lost under environmental conditions was analyzed using conservative biomarker concentrations for a 1979 oil spill in northern Minnesota, USA. After 40–41 years, an average of 50% of the mass was lost with values ranging from 22% to 57% depending on location. It is also important to understand the composition changes in the source. To understand controls on the losses of individual compounds, concentrations of volatile hydrocarbons in oil samples were compared with aqueous solubilities, and pore-space oil saturations. The results of the comparison show that losses of the oil compounds were controlled by pore-space oil saturations, solubility, and susceptibility to degradation under methanogenic conditions. Compounds that degrade under methanogenic conditions, including toluene, o -xylene, and n -alkanes are more depleted compared to benzene, ethylbenzene, and m - and p -xylene for which losses are dominated by dissolution. These rates and compound-specific behaviors form a foundation for improved modeling approaches and risk analyses.

Minnesota↗

USGS2026

Riverscape heterogeneity shapes population diversity for a migratory fish

Habitat patch dynamics can scale up to influence population demography and diversity with implications for resilience to environmental stochasticity. But how the spatial arrangement and size of habitat patches interact with other components of habitat heterogeneity to shape population diversity at larger spatial scales is not well understood. For riverine fishes, there is increasing evidence that tributary streams provide critical demographic support to main stem rivers. However, the extent to which main stem rivers rely on demographic contributions from tributaries, and the factors underlying this dependence, have not been assessed. Here, we used genetic stock identification to evaluate the effect of tributaries on population diversity of Yellowstone cutthroat trout ( Oncorhynchus virginalis bouvieri ) occupying the main stem Snake River, Wyoming, USA. We found that the main stem relied almost entirely on tributaries for demographic support, but main stem composition varied spatially among river sections. Distance between habitat patches, catchment area, and groundwater availability acted in concert to determine the contribution of specific tributaries to the main stem, but contributions were ultimately modulated by habitat connectivity. We also found evidence for multi-scale spatial structure in tributary contributions, providing insight into untested drivers of main stem river population diversity. Our results demonstrate how spatially discrete and distributed riverscape attributes influence population diversity at broader spatial scales, illustrating how ecosystem resilience emerges from the dynamic, two-way exchange of individuals and energy across habitat networks. Management plans for large rivers that address the ecological contributions of tributaries may be needed to achieve optimal outcomes. Similarly, conservation strategies that exclusively focus on headwater streams may fail to capture the broader habitat requirements necessary to maintain robust cold-water fish populations and associated recreational fisheries, particularly under global environmental change.

Wyoming↗

USGS2026

Simulation of groundwater flow to evaluate hydrogeologic controls on a PFAS plume, Coakley Landfill Superfund Site, Rockingham County, New Hampshire

Per- and polyfluoroalkyl substances (PFAS), including perfluorooctanoic acid (PFOA) and perfluorooctanesulfonic acid (PFOS), have been detected at combined concentrations above 2,000 nanograms per liter (ng/L) at groundwater seep locations near the Coakley Landfill Superfund site, in North Hampton, New Hampshire. The landfill was active from 1972 to 1985. An impermeable cap was placed on the landfill in 1998. The adjacent area to the Coakley Landfill has many water supply wells, and transport of PFAS compounds to the wells is a concern. Fracture anisotropy in the underlying bedrock aquifer complicates the understanding of PFAS transport because groundwater preferentially travels along fractures that may not align with the prevailing groundwater flow direction. In 2018, the U.S Environmental Protection Agency and the U.S. Geological Survey began an investigation of the groundwater flow from the Coakley Landfill site. This report describes the modification of a numerical groundwater-flow model for the local area around the Coakley Landfill and summarizes findings of the investigation. In addition, this report includes a brief description of PFOA and PFOS occurrence, a discussion of model construction, evaluation of model performance through calibration, and discussion of simulation results for two periods (before and after capping). Limitations are also discussed. Results show that simulated groundwater flow moves from the Coakley Landfill to the west and north. Advective transport modeling using particle tracking shows that groundwater from the landfill discharges primarily to streams to the west and north, and a small amount is transported to distal wells. Dilution of contaminants through advection and dispersion likely plays a role in whether PFAS compounds from the landfill will be detected above laboratory reporting levels at distal wells.

EarthArXiv↗

USGS2025

Speleothem evidence for Late Miocene extreme Arctic amplification – An analogue for near-future anthropogenic climate change?

The Miocene provides an excellent climatic analogue for near-future runaway anthropogenic warming, with atmospheric CO 2 concentrations and global average temperatures similar to those projected for the coming century under extreme-emissions scenarios. However, the magnitude of Miocene Arctic warming remains unclear due to the scarcity of reliable proxy data. Here we use stable oxygen isotope and trace element analyses, alongside clumped isotope and fluid inclusion palaeothermometry of speleothems to reconstruct palaeo-environmental conditions near the Siberian Arctic coast during the Tortonian (8.68 ± 0.09 Ma ). Stable oxygen isotope records suggest warmer-than-present temperatures. This is supported by temperature estimates based on clumped isotopes and fluid inclusions giving mean annual air temperatures between + 6.6 and + 11.1 °C , compared with − 12.3 °C today. Trace elements records reveal a highly seasonal hydrological environment. Our estimate of > 18 °C of Arctic warming supports the wider consensus of a warmer-than-present Miocene and provides a rare palaeo-analogue for future Arctic amplification under high-emissions scenarios. The reconstructed increase in mean surface temperature far exceeds temperatures projected in fully coupled global climate models, even under extreme-emissions scenarios. Given that climate models have consistently underestimated the extent of recent Arctic amplification, our proxy data suggest Arctic warming may exceed current projections.

Lena River delta↗

USGS2026

Shallow hydrogeologic framework of the Tully Valley mudboil area, Onondaga County, New York

Mudboils have been documented in the Tully Valley in southern Onondaga County, New York, since the late 1890s. Sediment-laden water from the mudboils flows into Onondaga Creek, which empties into Onondaga Lake at Syracuse 15 miles to the north. Turbidity from the mudboils has degraded the water quality of Onondaga Creek despite a series of mitigation efforts that began in the early 1990s. Turbidity mitigation actions presently (2025) being considered include creek relocation and offline sediment settling. In support of these proposed actions during 2021–23, the U.S. Geological Survey, in cooperation with the New York State Department of Environmental Conservation, U.S. Environmental Protection Agency, Onondaga Nation, Onondaga Environmental Institute, and Central New York Regional Planning and Development Board, collected and analyzed geologic, hydrologic, geophysical, and geotechnical data to characterize the shallow hydrogeology along four proposed creek-relocation paths and in the proposed offline settling basin area. The investigation indicated that the four proposed creek-relocation paths, two east of Onondaga Creek and two west of Onondaga Creek, are underlain by sediments including muck, alluvium, mudboil deposits, alluvial-fan sand and gravel, and lacustrine fines. The proposed excavations would penetrate partially to fully saturated conditions: generally, the water table is shallow near the creek and deep on the alluvial fans. The shallowest excavation, about 5 feet below land surface, would be near the creek and primarily in alluvium, and the deepest excavation, as much as 30 feet below land surface, would be in the alluvial-fan deposits. Brackish waters would be penetrated by proposed channel excavations on the eastern side of Onondaga Creek in an area downgradient from a potentially leaking historical salt-exploration borehole and near the main mudboil area. Excavation in these areas likely would provide a continuous source of brackish groundwater to the relocated creek. Proposed channel excavations of muck, soft to very soft lacustrine fines, and mudboil-type sediments in mudboil and suspected mudboil areas would pose an excavation and slope stability challenge and would have the greatest potential to create new mudboils. Proposed channel excavations below the water table on the Rattlesnake Gulf and Rainbow Creek alluvial fans would intercept groundwater and make the constructed streambank susceptible to seepage-induced slope instability. The substantial water-level fluctuation in the sediments of both alluvial fans would aggravate the stability condition. In addition, excavation on the Rattlesnake Gulf alluvial fan would have the potential to affect water-supply springs at the toe of the fan. The proposed offline settling basin area is in the northern part of the Rattlesnake Gulf alluvial fan. Natural and man-made diversions of Rattlesnake Gulf have resulted in saturated conditions in the general area of the proposed basin. The proposed offline settling basin would be excavated in, and berms would be constructed on, alluvial-fan deposits and lacustrine fines. In the proposed basin area, the alluvial deposits overlying the lacustrine fines are less than 10 feet thick. Excavation, berm construction, and loading of the saturated, soft to very soft lacustrine fines may be problematic and require soil strengthening.

New York↗