Search USGS⌕ Search

SEARCH · Search USGS

Search Search USGS

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.

13 recordsLinked to original sources
USGS2026

Investigation of drivers and sources of potential harmful algal blooms in Lake Anna, Virginia, 2023-24

Algal blooms in freshwater reservoirs are an increasing concern because of the risk to aquatic ecosystem health, recreational use, and water quality. Among the taxa that comprise an algal bloom, cyanobacteria are of particular concern due to the potential to produce toxic compounds which can cause acute and chronic illness in humans if ingested through contaminated water or shellfish. When toxins are above a toxin-specific threshold in the water column, the Commonwealth of Virginia categorizes the bloom as a harmful algal bloom (HAB), distinguishing it from an algal bloom without toxins present. Toxins are of particular concern in bodies of water which serve as major recreational areas, such as Lake Anna in central Virginia. Lake Anna is a 27 km-long and 47 km2 reservoir in central Virginia. Since the Commonwealth of Virginia revised the HAB monitoring framework in 2018, algal bloom advisories for potentially toxigenic blooms have been issued in the western, riverine zone of the reservoir. There were detections of microcystin and anatoxin production but no detections of cylindrospermopsin and saxitoxin production in Lake Anna; however, toxin concentrations were low, at values less than 0.65 µg/L. In 2023, the U.S. Geological Survey, in cooperation with Virginia Department of Environmental Quality, initiated a 19-month intensive and extensive study of western Lake Anna and its tributaries (North Anna River and Pamunkey Creek) to (1) characterize the algal community, biomass, and toxin production; (2) identify potential drivers leading to algal bloom initiation, persistence, and decline; and (3) assess the primary sources of the drivers contributing to algal blooms. Continuous and discrete monitoring of water quality within the lake, watershed inputs, meteorological conditions, and algal communities were monitored. Multiple modeling scenarios identified that the formation and persistence of algal blooms in Lake Anna are driven by interactions among chemical, macro- and micronutrient, and physical factors. Key drivers include total nitrogen, total phosphorus, water temperature, suspended sediment concentration, wind speed, copper, iron, molybdenum, nickel, sodium, and alkalinity. During this study, winter stormflow events delivered the highest loading of macro- and micronutrients from the tributaries to the lake. The lake is a sink for these nutrients in the water column or in the lake-bed sediments. During high-occupancy and recreation periods on the lake in the summer months, boating activity can promote internal loading by the resuspension of the lake-bed sediments, thereby increasing the nutrient availability for algal uptake. Thermal stratification and hypoxic conditions were observed, which can contribute to the passive release of nutrients into the water column, further facilitating algal bloom proliferation. These processes may create a positive feedback loop that supports the formation, persistence, and annual recurrence of algal blooms in the western region of Lake Anna. Continued water-quality monitoring of early warning indicators may improve the detection and forecasting of algal growth in the lake and help managers proactively manage water quality and algal growth.

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

Decadal trends in the quality of groundwater used for public drinking-water supply in California, 2004–2023, California groundwater ambient monitoring and assessment program, priority basin project

This study provides a comprehensive assessment of decadal changes in the quality of groundwater used for public drinking-water supply at 444 monitoring sites across California during 2004–2023. We assessed decadal step trends in groundwater quality for 145 water-quality constituents and geochemical indicators statewide and across geographic and land-use based network groups. We evaluated the statistical significance of directional changes (predominant increase or decrease of constituent concentrations) and the magnitude of those changes across all network groups. Uranium showed the most widespread directional and high-magnitude increases of all constituents with regulatory benchmarks statewide, particularly in the agriculture-dominated Central Valley as well as urban- and desert-dominated regions of Southern California. Fluoride and perchlorate showed the most widespread directional and high-magnitude decreases of all constituents with regulatory benchmarks statewide, which were also most pronounced in Southern California. Although arsenic and nitrate did not often register significant directional changes across network groups, they showed widespread, high-magnitude changes in both directions (increase and decrease) at levels often exceeding 10 percent of respective regulatory benchmarks statewide. Triazine herbicides (atrazine and simazine) and the gasoline oxygenate methyl tert-butyl ether (MTBE) showed significant directional decreases statewide, but not at levels considered to be of high magnitude compared to respective regulatory benchmarks. We observed significant directional and high-magnitude increases of total dissolved solids (TDS) statewide, which were most pronounced in agricultural areas. Analysis of explanatory geochemical indicators indicated that prevalent statewide increases of alkalinity and calcium were the predominant components of the observed statewide increases in TDS by mass. Widespread increases in groundwater alkalinity and calcium across agricultural and urban areas may be related, in part, to warm-season irrigation and other anthropogenic factors that have shifted soil weathering dynamics over the long term. Increasing alkalinity concentrations were related to increasing uranium concentrations, particularly in areas with aquifer materials derived from granitic rocks. Conversely, increasing calcium concentrations were related to decreasing fluoride concentrations, particularly in areas where fluoride occurred naturally at elevated concentrations. Decrease of perchlorate, triazine herbicides, and MTBE are likely related to decreased anthropogenic source inputs over time and natural attenuation in aquifers.

EarthArXiv↗
USGS2026

Effects of groundwater withdrawals for water bottling and municipal use, Wards Brook Valley, Maine and New Hampshire

Hydrologic models for the Wards Brook valley near Fryeburg, Maine were developed for historical (2016 – 2021) and hypothetical future conditions (2046 – 2065 and 2080 – 2099) to understand the effects of groundwater withdrawals for bottled water and municipal use on hydrologic conditions (stream base flows and groundwater levels). Analyses showed that the simulated base flows in Wards Brook were reduced because of pumping for both municipal water supplies and for water bottling, and about half of the total pumping impact on the base flows in Wards Brook was from the bottled water extraction. Simulated flows were greater than the minimum recommended streamflow of 2,180 cubic meters per day (400 gallons per minute) throughout the historical period. Simulated groundwater levels at two of three nearby ponds (Round Pond and Davis Pond) were minimally affected by pumping conditions, and effects were primarily from the municipal well closest to the ponds. Several estimates of future projected recharge were used to understand the potential effects of groundwater withdrawals on hydrologic conditions under multiple hypothetical climate conditions. Annual projected recharge rates in the mid- and late-21st century from two climate scenarios (stabilized greenhouse-gas emissions and high greenhouse-gas emissions) were similar to rates for 2016 – 2021. However, monthly recharge patterns for the future periods shifted toward more recharge in the winter months (December, January, and February) and less recharge in April, May, and October relative to 2016 – 2021. The lowest mean monthly base flows from the future emission scenarios all remain larger than the minimum recommended streamflow and indicate no long-term declines in flow relative to historical conditions. However, simulated base flows during hypothetical 3-year drought scenarios declined below minimum recommended streamflow during the summer months in the stabilized- and high-emission scenarios in the mid-21st century. Although water is generally plentiful in the Wards Brook valley, reduced pumping may be needed to maintain streamflows in Wards Brook under future climate conditions similar to modeled drought scenarios.

EarthArXiv↗
USGS2026

Multiple-well monitoring site adjacent to the Midway- Sunset and Buena Vista Oil Fields, Kern County, California

Groundwater quality in and around oil fields in the Southern San Joaquin Valley is of interest to many California residents that rely heavily on groundwater for domestic, commercial, and agricultural use. To help assess the effects of historical oil-field activities and natural geologic sources on groundwater near the southwest margins of the Kern County Groundwater Subbasin, a multiple-well monitoring site was installed near the administrative boundary between the Midway-Sunset and Buena Vista Oil Fields in Kern County, California. The installation of the Midway-Sunset Buena Vista multiple-well monitoring site (MSBV) supports regional analysis of the relations of oil and gas sources to groundwater quality by providing information about the geology, hydrology, geophysical properties, and water quality of the alluvial and upper Tulare aquifers in areas where groundwater data were limited. Data collected from the site included drill cuttings, whole core samples, sidewall core samples, mud-gas analysis, borehole geophysical logs, depth to water measurements, and water quality samples. Whole cores were scanned using dual energy computed tomography. Subsamples of selected cores were analyzed for density, porosity, specific retention, and bulk minerology. Thin sections of the subsamples were prepared, photographed, and examined. Two samples were analyzed using scanning electron microscope technology to examine the microporosity of diatomite laden sediment. Instrumentation installed in the wells collect hourly depth to water measurements. Analysis of the data show there is 355 feet of alluvium overlying the Tulare Formation at the well site. The contact between the two formations is an aquitard resulting in a perched aquifer in the alluvium and unconfined aquifer in the Tulare Formation. The alluvium is more heterogenous and finer grained than the Tulare Formation resulting in markedly higher porosity in the alluvium compared to the Tulare Formation. Higher specific retention observed in the alluvium is attributed to the finer grained sediment and greater abundance of reworked diatomite (as represented by opal-CT [cristobalite-tridymite]) compared to the Tulare Formation. Total dissolved solids (TDS) approached or exceeded 10,000 milligrams per liter (mg/L) in the alluvium from approximately 176 to 242 feet below land surface and at the top of the Amnicola clay at approximately 670 feet below land surface within the Tulare Formation. Elevated TDS, chloride, and boron concentrations in the alluvium and on top of the Amnicola clay likely reflect groundwater that is mixed with oil-field water. Water chemistry and modern-aged groundwater in the alluvial monitoring well (MSBV #3) are consistent with the oil-field water in the alluvium being derived from documented historical surface disposal of oil-field water upslope (northwest) of the site. Water chemistry and pre-modern groundwater age in the deeper Tulare monitoring well (MSBV #1) on top of the Amnicola clay are consistent with oil-field fluids derived from upslope natural geologic sources or old oil wells that leak in the subsurface. Shallow groundwater in the Tulare (MSBV #2) is not affected by mixing with oil-field sources.

EarthArXiv↗
USGS2026

Computation of regional groundwater budgets for the Virginia Coastal Plain aquifer system

Computation of detailed groundwater flow budgets for subdivisions of Virginia’s Coastal Plain aquifer system has enabled quantification and more thorough understanding of groundwater flow within this important water resource. A zone budget analysis conducted on previously published groundwater models of the Virginia Coastal Plain and Virginia Eastern Shore shows that groundwater conditions vary substantially throughout the Coastal Plain aquifer system due to local variations in hydrogeology and historical and ongoing variations in groundwater use and management. Decades of substantial groundwater withdrawal from the Coastal Plain aquifer system have fundamentally altered groundwater flow from pre-development conditions. Rates of sustainable withdrawal are limited because the downward groundwater flow rate into confined aquifers supplying groundwater is a relatively small portion of the total groundwater water budget for the aquifer system. Analyses of groundwater budgets from the Virginia Coastal Plain model show that groundwater flow is generally outward from the surficial aquifer to rivers and coastal water bodies and downward through a series of underlying aquifers and confining units to the Potomac aquifer, which is the deepest aquifer and the source of most groundwater withdrawals. Downward flow into the Potomac aquifer currently is estimated to be only 7 percent of total net precipitation-derived net recharge at the land surface but makes up about 66 percent of inflow to the aquifer in Virginia, with much of the remaining inflow occurring laterally from areas outside of defined groundwater budget regions in Virginia. For several decades prior to 2010, high rates of withdrawal from the Potomac aquifer resulted in substantial decline in groundwater storage in the aquifer and in most overlying aquifers and confining units. From 2010 to 2025, rates of withdrawal substantially lower than the historical maximum have resulted in small net increases in groundwater storage in the confined aquifer system for most regions of the Virginia Coastal Plain. Nevertheless, for the same period, groundwater storage for the entire model domain continues to incrementally decline, indicating that storage recovery in Virginia is offset by a continued decrease in storage in areas beneath the Chesapeake Bay or in adjacent areas of Maryland and North Carolina. Withdrawals from the Potomac aquifer have induced substantial downward flow which is a large part of groundwater budgets for confined aquifers such as the Potomac. Downward groundwater flow continues under current conditions, but because vertical flow rates are a function of the difference between water pressure in the upper surficial systems and lower confined units, those rates are lower than those in earlier decades as the confined water levels partially recover from larger groundwater withdrawals in the past. Geographically, groundwater flow is generally inward from perimeter regions of the Virginia Coastal Plain toward central regions with the largest withdrawal rates. Estimated groundwater inflow from coastal regions could be contributing to saltwater intrusion, though that was not measured directly in this study. Analyses of groundwater budgets from the Virginia Eastern Shore peninsula, a geographic region of the Virginia Coastal Plain, show that groundwater flow for that isolated aquifer system is generally outward from the surficial aquifer to coastal water bodies and downward into the confined Yorktown-Eastover aquifer system, which is the source of most withdrawals. Downward groundwater flow into the confined Yorktown-Eastover aquifer system is estimated to be less than 2 percent of total recharge and less than 9 percent of net recharge at the water table but makes up over 93 percent of all inflow to the confined aquifer system. Decades of substantial but relatively consistent groundwater withdrawals have induced greater downward flow rates into the confined aquifer system but also have resulted in loss of groundwater from storage. Currently, estimated storage loss accounts for slightly under 7 percent of withdrawals from the confined aquifer system. The current withdrawal rate from the confined Yorktown-Eastover system is near the highest reported rate for the Eastern Shore, which means that the storage depletion is expected to continue, even though groundwater levels appear to be relatively stable. Estimated groundwater flow rates upward from the confining unit underlying the Yorktown-Eastover system and small rates of inflow from coastal water bodies underscore ongoing concerns about up-coning and lateral intrusion of salty groundwater.

EarthArXiv↗
USGS2025

Regional characterization of coal resources in the U.S. Gulf Coast

There is increasing interest in extracting critical minerals (CM), including rare earth elements (REE), from coals in the United States to address the overreliance on imported REE. The U.S. Gulf Coast and the Williston basins are the two major lignite-bearing basins within the country. Recent REE and CM studies of the lignite in these basins have indicated that the coals may be a viable source material for REE and CM extraction. To evaluate in-place coal as a potential source of REE and CM, the coal resources need to be quantified. This study presents the results of a regional analysis of the U.S. Gulf Coast lignite and bituminous coal resources that might be available as potential sources of REE and CM. The resource analysis used kriging methods to develop isopleth maps of cumulative coal thickness throughout the region using data from 31,181 drill holes and other data points. The estimated total coal resource in the Gulf Coast is about 83 billion metric tons in the upper 90 m (~ 300 ft) of the subsurface. Texas accounted for 40 percent (32 billion metric tons) of the total resource, followed by Mississippi (24 %, 20 billion metric tons), Louisiana (14 %, 12 billion metric tons), Tennessee (10 %, 8.5 billion metric tons), and Arkansas (6 %, 5.1 billion metric tons). The remaining states each accounted for less than 5 percent of the total resource. Georgia had the smallest resource estimated at 7 million metric tons. Here we report the first known state-wide lignite resource estimates for Georgia, Kentucky (820 million metric tons), and Missouri (1,800 million metric tons). A comparison of the results of this study with those of previous Gulf Coast and Williston Basin resource studies is difficult because each study used different data sources, assessment methodologies, overburden depths, and qualifying coal thicknesses. Coal-power electric generation has sharply decreased in past decades and mining of these coals for CM and REE could provide additional co-products such as activated carbon and other uses such as fertilizer (soil enhancer).

EarthArXiv↗
USGS2025

A spatial analysis of the groundwater emergence flood hazard in Long Island, New York and near coastal areas surrounding Long Island Sound in New York, Connecticut, and Rhode Island

Long Island, New York and near coastal areas surrounding Long Island Sound are densely populated and, like other coastal areas, are susceptible to flooding from several potential sources, including stormwater from precipitation events, tidal flooding and storm surge, and groundwater inundation or groundwater emergence flooding. The latter refers to the intersection of a rising water table with land surface or critical infrastructure. Many studies of flood drivers either neglect or only briefly discuss how shallow groundwater conditions may contribute to or exacerbate flood conditions. As part of a comprehensive study of compound flood hazards in the near coastal areas surrounding Long Island and Long Island Sound, a spatial analysis was completed, in cooperation with the Environmental Protection Agency’s Long Island Sound Study, using available regional datasets to characterize the potential hazard for groundwater emergence flooding. The approximately 3,100 square mile study area was subdivided into 11,407 900-meter by 900-meter (approximately 3,000-feet by 3,000-feet) grid cells, for the purposes of integrating the spatial datasets to calculate and map the groundwater emergence flood hazard. The depth to the water table, hydrologic soil groups, and National Land Cover Database were harmonized to the common grid. A groundwater emergence flood hazard rank was calculated for each grid cell for current average conditions following a set of rules accounting for the depth to the water table and the percent of area within each cell with slow infiltrating soils. A higher sea level position scenario was also calculated for the Long Island part of the study area. The calculated groundwater emergence flood hazard rank was reviewed in concert with the National Land Cover Data Base to identify developed areas and associated infrastructure that may be at risk to groundwater emergence flooding. Study results indicate that the groundwater emergence flood hazard is highest in coastal areas and near surface water where the water table is close to ground surface. Inland areas away from surface water bodies are not likely to be exposed to groundwater emergence flooding. For Long Island, under a scenario with higher sea level position, a greater groundwater emergence flood hazard is calculated in some locations closer to the coast and where land is submerged. Away from the coast and surface-water drainage, the groundwater emergence flood hazard is similar between the current average sea level condition and a higher sea level position scenario.

EarthArXiv↗
USGS2025

Salinas Valley integrated hydrologic and reservoir operations models, Monterey and San Luis Obispo Counties, California

The area surrounding the Salinas Valley groundwater basin in Monterey and San Luis Obispo Counties of California is a highly productive agricultural area, contributes significantly to the local economy, and provides a substantial portion of vegetables and other agricultural commodities to the Nation. This region of California provides about half of the Nation’s lettuce, celery, broccoli, and spinach each year. Thus, this agricultural area provides significant volumes of agricultural products not just for California but the entire United States. Changes in population and increased agricultural development, which includes a shift toward more water-intensive crops, and climate variability, have put increasing demand on both surface water and groundwater resources in the valley. This has resulted in water management challenges in the Salinas Valley that are predominantly related to distribution of water supply throughout the basin. Where and when the water is present in the surface and subsurface does not coincide with where and when the water is needed. To deal with the distribution issue, historically water has been used conjunctively in the valley. Conjunctive use is a water management strategy that coordinates surface water and groundwater use to maximize water availability. Groundwater is used throughout the Salinas Valley to meet water demands when surface water supplies are insufficient. Availability of surface water is constrained by climate. Precipitation and streamflow vary seasonally and year to year. Although there are two reservoirs in the Salinas Valley to capture and store water during wet periods, the only conveyance of reservoir water to coastal agricultural areas is the Salinas River. Increasing demand on groundwater and surface water resources throughout the Salinas Valley has resulted in undesirable effects of unsustainable water use, such as surface water depletion, groundwater level declines, storage depletion in the principal aquifers, and seawater intrusion. To address these escalating issues, local communities, water management agencies, and groundwater sustainability agencies are evaluating how to sustainably manage both their surface water and groundwater resources. To meet water demands and reduce undesirable effects of unsustainable water use, continued conjunctive management of surface water and groundwater would ideally incorporate strategies to deal with increases in demand and a variable climate. To evaluate the challenging water management issues in the Salinas Valley, the U.S. Geological Survey, Monterey County Water Resource Agency, and the Salinas Valley Basin Groundwater Sustainability Agency developed a comprehensive suite of models that represent the Salinas Valley Hydrogeologic system called the Salinas Valley System Model. The Salinas Valley Geologic Framework was developed to characterize the subsurface using various topographic and geologic data sources, including information on hydrogeologic units, their surfaces and extents, geologic structures, lithology, and elevations from borehole data and cross sections, as well as details on faults and existing models. The Salinas Valley Watershed Model simulates the entire Salinas River watershed. Monthly surface water inflows into the integrated hydrologic model domain were simulated using the Salinas Valley Watershed Model. The historical model uses historical climate data, water and land use data, and reservoir releases to simulate agricultural operations, including landscape water demands, diversions, and reclaimed wastewater. The operational model adds an embedded reservoir operations framework to the simulation of the historical model that allows specified operational rules to simulate reservoir releases and changes in reservoir storage. The operational model assumes current reservoir operations and constant land use, which differs from historical conditions. Thus, the operational model is a hypothetical baseline model that can be used by local water managers to evaluate and quantify potential benefits of water supply projects. Together, the geologic framework, watershed, historical, and operational models form a tool that can be used to simulate irrigated agriculture and associated reservoir operations of the integrated hydrologic system of the Salinas Valley.

EarthArXiv↗
USGS2025

Simulation of the impacts of spring diversions on streamflow in the Strawberry Creek watershed, San Bernardino County, California, using an integrated hydrological model

The Strawberry Creek watershed, situated in the San Bernardino Mountains of southern California, features a group of natural springs known as Arrowhead Springs that have been augmented with diversions in the form of sub-horizontal borings and tunnels. Understanding the impact of these structures on streamflow through groundwater capture is crucial for managing surface-water resources in this watershed. In this study we constructed the Strawberry Creek integrated hydrological model (SCIHM) to increase this understanding. The SCIHM is an integrated surface runoff and groundwater model that uses the coupled groundwater and surface-water flow model (GSFLOW), which is based on the integration of the precipitation-runoff modeling system (PRMS) and the modular groundwater flow model commonly called MODFLOW, version MODFLOW-2005 software to simulate surface runoff and infiltration and groundwater flow. The model has three layers, 263 rows, and 176 columns. The model area includes the Strawberry Creek and four adjacent watersheds. The PRMS model was calibrated using two streamflow gaging stations and the GSFLOW model was calibrated to reported spring diversion discharge and a sparse number of groundwater-level measurements. The SCIHM was run with and without diversions active and simulated streamflow was compared, finding that in the headwaters of Strawberry Creek about 35 percent of the diversion flow was captured from base flow.

EarthArXiv↗
USGS2025

Assessing legacy nitrogen in groundwater using numerical models of the Long Island aquifer system, New York

Nitrogen transported along groundwater flow paths in coastal aquifers can contribute substantially to nitrogen loading into surface water receptors, particularly in hydrologic systems dominated by groundwater discharge. Nitrogen entrained in the aquifer is a function of land use and associated nitrogen sources at the time of groundwater recharge, which may differ considerably from present-day sources. Legacy nitrogen can result in substantial discrepancies between observed present-day nitrogen loading to surface water receptors and loading estimated from present-day sources. Additionally, legacy nitrogen can continue to discharge into surface waters after nitrogen mitigation actions have been undertaken. Here, we use a numerical modeling framework to compare three methods of estimating time-varying historical nitrogen loads to four water bodies (receptors) on eastern Long Island, New York. The methods span a range of data requirements and process complexity, from instantaneous receptor loads calculated from steady-state groundwater contributing areas, to transient loads estimated by explicitly simulating legacy groundwater nitrogen transport over a century with large changes in nitrogen sources and hydrologic conditions. The effects of legacy nitrogen on estimated receptor loads varied temporally and spatially within the study area. Depending on antecedent nitrogen inputs and hydrologic conditions, historical annual nitrogen loads estimated from transient simulations accounting for legacy nitrogen can be quite similar (<10% difference) or substantially different (±100%) from those estimated from simpler instantaneous methods. Continued input of present-day nitrogen sources using methods that account for legacy nitrogen results in asymptotic increases in receptor nitrogen loads over time, indicating that simulated present-day receptor nitrogen loads are not in equilibrium with present-day inputs. For these receptors in disequilibrium, models simulating transient groundwater nitrogen transport could be used to account for legacy nitrogen lag times to help resource managers evaluate the potential effectiveness of proposed nitrogen mitigation actions.

EarthArXiv↗
USGS2025

Wetter winters, drier summers: Quantifying the change in hydrological response around the Puget Sound area using the wflow_sbm hydrological model and CMIP6 projections

Climate change is expected to impact hydrological regimes worldwide, including the Pacific Northwest of the United States. This study investigates how climate change will affect river discharge in the Puget Sound region of the State of Washington, with a focus on King and Pierce Counties. We simulated river discharge under historical and future conditions using the physically based, spatially distributed wflow_sbm hydrological model, which was calibrated and validated against U.S. Geological Survey discharge records. Future forcing was based on an ensemble of six high-resolution CMIP6 climate models, which were bias corrected using empirical quantile mapping. The results indicate a decrease in summer discharges (5–10%) and an increase in winter discharges (5–10%) across the study region. The high discharges (90th percentile) are projected to increase in winter, and the low discharges are projected to decrease in summer, due to shifts in precipitation regimes, snowpack hydrology, and evapotranspiration. However, variability between individual CMIP6 models often exceeds the magnitude of ensemble mean changes, underscoring substantial uncertainty in climate projections and the importance of including multiple climate models in climate change analysis. Furthermore, model consensus increased with elevation, which could be the result of the higher elevation areas being driven by less diverse hydrological processes. These findings highlight potential challenges for regional water management, ecosystem health, and flood risk mitigation in the Puget Sound region under future climate conditions.

EarthArXiv↗
USGS2025

Tracing mercury from land to river: Global sources, retention, and implications for sustainability

Mercury (Hg) pollution in river systems is a global sustainability challenge. Yet the transport, transformation, and retention of Hg within global rivers remain poorly quantified, particularly in regions with sparse observations such as Southeast Asia and Africa, hindering effective pollution mitigation and reinforcing geographic inequities in scientific knowledge and environmental governance. Here, we present the first global, high-resolution simulation of riverine Hg dynamics using a process-based model that traces Hg from land-based sources through river networks to the ocean. Under a realistic scenario, we estimate that ~1,900 megagrams per year (Mg/yr) of Hg enters global rivers, including 1,500 Mg/yr from human-induced sources and 400 Mg/yr from soil erosion. Nearly half of this flux (~1,000 Mg/yr) is retained in reservoirs and dams, which act as major sinks. While such retention limits downstream delivery to the oceans, it also heightens in-reservoir Hg methylation risks. By bridging the gap between Hg releases and observed riverine exports, our framework offers a scalable tool for data-limited regions, promotes data access, and supports global freshwater and pollution-management strategies.

EarthArXiv↗