USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research: Earth Surface”

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

Status and understanding of groundwater quality in the San Joaquin Valley Kern County subbasin domestic-supply aquifer study unit, 2022—California GAMA Priority Basin Project

The quality of water accessed by domestic wells (here referred to as domestic groundwater resources) in the San Joaquin Valley Kern County subbasin (basin number 5-022.14) was assessed as part of the California Groundwater Ambient Monitoring and Assessment (GAMA) Program Priority Basin Project (GAMA-PBP), in cooperation with the California State Water Resources Control Board. Kern County is at the southern end of the San Joaquin Valley in California, and about 30,000 residents are estimated to use privately owned domestic wells for drinking water. Domestic wells typically draw from shallower parts of the aquifer system than public-supply wells and can be more vulnerable to effects from surface activities. Kern County is host to a highly productive agricultural industry, with Bakersfield as the main urban center. The Kern River runs through Bakersfield from the southern Sierra Nevada and intersects the Kern Water Bank, one of the largest groundwater banking operations in California, at the Kern River Intertie. The section of the Kern River running through the Kern Water Bank is dry most years. Kern County also encompasses some of the most productive oil and gas basins in California, with extensive underground and surface disposal of oil-field wastewater. This study was based on data collected from 33 sites sampled by the U.S. Geological Survey for the GAMA-PBP in 2022. To provide context for the water quality assessment, measured concentrations were compared to regulatory and non-regulatory health-based and aesthetic benchmarks. A grid-based method was used to estimate the proportions of the groundwater resources used for domestic-supply wells that have water-quality constituents below (low relative concentration), approaching (moderate relative concentration), or above (high relative concentration) benchmark concentrations. At least one measured constituent with a regulatory benchmark was categorized as having a high relative concentration in 72 percent of the aquifer area used for domestic groundwater resources. Inorganic constituents were detected at high concentrations in 45 percent of the domestic groundwater resources, and the constituents detected above regulatory benchmarks were arsenic, nitrate, and uranium. At least one organic constituent was detected at high concentrations in 41 percent of the domestic groundwater resources, and the constituents exceeding regulatory benchmarks were the fumigants 1,2,3-trichloropropane (1,2,3-TCP), 1,2-dibromo-3-chloropropane (dibromochloropropane [DBCP]), 1,2-dibromoethane (EDB), and the per-and polyfluoroalkyl substance (PFAS) perfluorooctanesulfonate. The disinfection by-product chloroform, the fumigant 1,2-dichloropropane, the herbicides atrazine and hexazinone, and the herbicide degradates 2-chloro-6-ethylamino-4-amino-s-triazine, 2-chloro-4,6-diamino-s-triazine, 4-hydroxychlorothalonil, and metolachlor sulfonic acid were detected in more than 10 percent of domestic groundwater resources, but concentrations did not exceed regulatory benchmarks. Land use, groundwater age (fraction of modern water and mean age), and geochemical environment (oxic or anoxic conditions, pH, alkalinity) were associated with the distribution of high relative concentrations of inorganic and organic constituents. Young, oxygenated water is recharged along the Kern River and adjacent recharge ponds, or as irrigation water in the agricultural areas. High concentrations of nitrate and volatile organic compounds occurred in the oxic water in urban and agricultural areas. The fumigants 1,2,3-TCP, DBCP, and EDB were reported throughout the agricultural areas, whereas chloroform, tetrachloroethene, and PFAS were associated with urban land use. High uranium concentrations were associated with young, modern groundwater in agricultural areas with low pH and high bicarbonate. Total dissolved solids increased with distance from the Kern River, as the contributions of fresh, oxic water decreased. High concentrations of arsenic were present in older anoxic or alkaline groundwater away from areas of recharge. Overall, groundwater age, redox conditions, and the source of recharge as a result of different land uses contribute to large aquifer-scale portions of domestic groundwater resources that exceed health-based benchmarks for nitrate, uranium, and fumigant concentrations.

California↗
USGS2024

Numerical model of the groundwater-flow system near the southeastern part of Puget Sound, Washington

Groundwater flow in the active model area (AMA) was simulated using a groundwater-flow model. A steady-state model version of the model simulates equilibrium conditions, and a transient model version simulates monthly variability. The model corresponds to the physical and temporal dimensions of the conceptual model and groundwater budget. The steady-state model version represents average conditions for an 11-year period (January 1, 2005–December 31, 2015), and the transient model represents monthly hydrologic variability within that period. The 13-layer model was constructed using MODFLOW-NWT with a uniformly spaced grid consisting of 416 rows, 433 columns, and cells with a horizontal dimension of 500 feet (ft) on a side. The model was calibrated to measured values of water levels in wells and lakes and estimated base flow for selected streamflow measurement stations, commonly referred to as streamgages. Model calibration was accomplished using a combination of manual and automatic methods, including the Model-Independent Parameter Estimation (PEST) program that adjusted model input parameters with the aim of minimizing the difference between estimated and model-simulated values of hydraulic head and base flow. Model boundary conditions consist of all simulated groundwater inflow to and outflow from the AMA. For example, a stream reach that simulates a gain from or loss to groundwater is a boundary condition that allows water to exit or enter, respectively, the groundwater system. Other boundary conditions include springs, seeps, precipitation recharge, groundwater exchange with lakes and Puget Sound, and groundwater pumping. A comparison of the estimated groundwater budget to that simulated by the steady-state model version indicates that the relative percentages of total inflow or total outflow for six major categories of boundary conditions are similar for the two budgets. The model was used to simulate three suites of scenarios of potential drought and water-use changes. Scenario 1 suite consisted of the steady-state model version that was run with 0, 15, 20, and 25 percent reduction of precipitation recharge to assess the corresponding reductions in base flow with decreasing recharge. The last simulation for the scenario 1 suite consisted of the transient model version simulating 3 years of consecutive seasonal drought, defined by the months of May through September, to assess the corresponding base-flow reductions. Scenario 2 suite consisted of the steady-state model version with all simulated groundwater use removed, compared with a simulation that includes current groundwater use to evaluate changes to potentiometric surfaces and base flows. Scenario 3 suite consisted of a transient model version of the model that simulated pumping increases for four different categories of water-supply wells (compared to no pumping increases) to evaluate resulting reductions in base flow. Although, these scenarios provide examples of model applications and useful insights, many other scenarios could be simulated. A description of how to download the model is described in the body of this report. Uncertainty is associated with most model inputs. Groundwater levels, lake levels, and land-surface altitudes are relatively certain; other model inputs are far less certain, including precipitation recharge, base flow, hydraulic properties, water use, and the three-dimensional structure of subsurface hydrogeologic units. Models are useful not because of high levels of accuracy of all model inputs, but because they combine the best information and estimates available, thereby providing the best predictions available related to physical processes. The model described in this report simulates groundwater flow on a regional scale, which has inherent limitations for simulating hydrologic scenarios at local scales. Model structures and inputs were generalized to be consistent with this regional scale. For example, the actual groundwater system has much greater heterogeneity of hydraulic conductivity than is possible within the model’s degrees of freedom. Variations in hydraulic gradients over distances less than 500 ft cannot be simulated. The distances between model features, such as a pumping well and a stream, must be placed at 500-ft intervals and are co-located if both features are within the same model cell.

Washington↗
USGS2025

A 10-year continuous daily simulation of chloride flux from a suburban watershed in Fairfax County, Virginia, USA

Increasing levels of chloride in surface water are associated with detrimental effects on water quality, aquatic ecosystems, infrastructure, and human health. Numerous mass-balance studies have inferred watershed transport processes by interpreting chloride inputs and outputs, but few represent internal dynamics explicitly. We constructed a coupled water/chloride mass balance model to gain insights into storage, residence time, and transport processes in a 10-km 2 urban watershed. The model, which operates over a 10-year period at a daily time scale, represents storage in a dynamic soil-moisture reservoir, quick-flow runoff from storm events, and slow-flow runoff that sustains streamflow in dry weather. The calibrated model accurately represented (a)the observed transition from a streamflow enrichment regime in cold months to a dilution regime in warmer months, (b) the observed tendency for late-summer concentrations to be higher after winters with heavy snowfall, and (c) a period-of-record downward trend in chloride concentration likely associated with a downward trend in annual snowfall. Estimated chloride inputs averaged 195 metric tons per year, while the average output was 270 metric tons per year. In contrast, estimated storage was only 107 metric tons. The estimated mean residence time in groundwater was 1.27 years. This short residence time indicates that efforts to reduce inputs will manifest as decreased concentrations in streamflow on a management-relevant time scale of several years. The coupled mass balance model yielded insights into internal watershed dynamics that would not be possible from simple input/output analysis; such models can be useful tools for gaining insight into small watershed hydrology and pollutant transport.

Virginia↗
USGS1992

Public water-supply systems and water use in Tennessee, 1988

This report summarizes the results of a study conducted by the U.S. Geological Survey, in cooperation with the Tennessee Department of Environment and Conservation (TDEC), Division of Water Supply in 1988. Data gathered during an inventory by the TDEC were collated to determine water use, supply sources, population served, and design and storage capacities of the systems. The inventory was limited to systems that were active on June 30, 1988. Results of a survey of the systems conducted by the Tennessee Department of Health and Environment during 1988 were a primary source of data for this report. Data from computer and manual files maintained by the Tennessee Department of Health and Environment and the U.S. Geological Survey also were used. The Division of Water Supply, TDEC, surveyed 541 public water-supply systems. These systems served 81 percent of the population of the State, or 3.95 million people. The gross per capita use statewide for public-supplied water was 179 gallons per day. Total water withdrawals for public supply increased about 39 percent from 510 million gallons per day (Mgal/d) in 1980, to 708 Mgalld in 1988. During the same period, the population increased about 7 percent. Surface-water withdrawals accounted for 63 percent (446 Mgal/d) of the total water withdrawn in the State. All of these withdrawals occurred in the Tennessee (56 percent or 249 Mgal/d) and the Ohio (44 percent or 197 Mgalld) hydrologic regions. Ground water supplied 262 Mgal/d or 37 percent of the total water withdrawn by public-supply systems statewide. Of that amount, 79 percent, or 208 Mgalld, was used in western Tennessee.

Tennessee↗
USGS2025

“Leaky weirs” capture alluvial deposition and enhance seasonal mountain-front recharge in dryland streams

“Leaky weirs” are rock structures installed in dryland streams, which are anchored into exposed bedrock, loosely cemented, and designed to allow water to slowly pass through. They are being tested at a ranch in southeastern Arizona, USA, to restore and conserve the historic range and desert wetlands. Data are collected to assess how leaky weirs impact surface water, subsurface water, and groundwater recharge—including stream discharge, timing, and depth of infiltration, and groundwater elevations. Three adjacent watersheds, two with outlets just below leaky weirs and one with leaky weirs farther upstream, were instrumented with water-level loggers, wildlife cameras, and crest stage instruments with temperature sensors in the soil. As most groundwater recharge is assumed to be focused along the mountain fronts in this region, mountain-block recharge is also evaluated to differentiate between the two using isotope analyses. Finally, a single, late-season flood event is scrutinized to consider the leaky weir effect on all monitored components in the water budget. Results indicated groundwater flow is primarily from the mountains to the east via older, regional mountain-block recharge. However, the development of shallow alluvial aquifers is supported by the leaky weirs, that slow flows, capture permeable sediments, and allow infiltration, thus enhancing mountain-front recharge. In turn, these new pockets of water help support the restoration of historic wetlands. Sediment accumulates where leaky weirs are installed, reducing flashy peak flows, and resulting in a series of infiltration ponds along the channel that support vegetation during growing seasons and recharge the shallow aquifer during non-growing seasons.

Arizona↗
USGS2024

Volcanic gas monitoring

Introduction As magma rises through the crust, decreasing pressure conditions allow volatiles to exsolve from the magma. These volatiles then migrate upward through the crust, where they can be stored at shallower levels or escape to the atmosphere. Rising magma also heats rock masses beneath volcanic centers, causing water in shallow aquifers and hydrothermal systems to boil and release additional gases and steam (see chapter F , this volume; Ingebritsen and Hurwitz, 2024). The chemistry and quantity of gases that reach the surface during periods of quiescence or volcanic unrest can reveal that gas-rich magma is ascending, crystallizing, or alternatively stalling, with important implications for volcanic hazard (for example, Sutton and others, 1992; Aiuppa and others, 2007, 2021; Werner and others, 2009, 2011, 2012; Moretti and others, 2013; de Moor and others, 2016; Lewicki and others, 2019; Edmonds and others, 2022; Kern and others, 2022; Kunrat and others, 2022). Most volcanoes in Alaska and the western United States are characterized by weak degassing, with one or more low-temperature fumaroles (typically near the local boiling temperature of water) and connect to a deeper and sometimes extensive hydrothermal system (for example, McGee and others, 2001; Symonds and others, 2003a, b). Hydrothermal systems will affect the chemistry of rising gases exsolved from deeper magma (Symonds and others, 2001), including sulfur dioxide (SO 2 ), hydrogen chloride (HCl), and water vapor (for example, Doukas and Gerlach, 1995; Gerlach and others, 1998, 2008; Symonds and others, 2001; Werner and others, 2013). As an example, depending on factors such as temperature, pressure, and oxidation state, rising SO 2 will react with groundwater to form hydrogen sulfide (H 2 S) gas, dissolved sulfate (SO 4 2− ), or elemental sulfur (Christenson, 2000; Symonds and others, 2001; Werner and others, 2008). The reaction and dissolution of SO 2 into shallow groundwater is commonly referred to as scrubbing, and can reduce the likelihood that ascending, degassing magma can be detected. Carbon dioxide, however, in addition to exsolving from magma early in the ascent process, is not easily removed by hydrothermal fluids (Lowenstern, 2001). As scrubbing and other processes take place, the SO 2 /H 2 S, CO 2 /SO 2 , and CO 2 /H 2 S ratios may change. High rates of SO 2 emission indicate that magma has moved to relatively shallow levels in the volcano and that the system has heated up enough to establish dry pathways from depth to the surface. Monitoring multiple gas species and the total output of those species is thereby useful for volcano monitoring during both periods of quiescence, to establish background degassing conditions, and during unrest, when gas geochemistry and emission rates can provide information on changing conditions, such as magma ascent. To provide context for multidisciplinary volcano forecasts, we focus on the following two key required capabilities: (1) characterizing baseline geochemistry and gas discharge from volcanoes and volcanic regions and (2) monitoring changes in gas geochemistry and discharge to inform forecasts of volcanic eruptions and their effects. Sufficient baseline data must be collected to identify and interpret anomalous degassing associated with volcanic unrest (for example, Sorey and others, 1998; Rouwet and others, 2014). Differences in volcano type, baseline degassing rates, local hydrology, and geography (for example, high versus low latitude) will result in a different baseline for each volcano. Volcanoes of any threat level that exhibit one or more degassing phenomena would ideally be monitored by techniques needed to establish baseline degassing data, with the sampling frequency of baseline data dictated by the threat level (table E1). Additional monitoring techniques become necessary during periods of unrest. In general, three of the most important techniques for gas monitoring are (1) direct sampling of fumarole, spring, and soil gases for laboratory geochemical measurements, (2) measurements of the chemical composition of the volcanic plume and emission rates of major gas species (for example, H 2 O, CO 2 , SO 2 , and H 2 S) by satellite, airborne, or ground-based techniques, and (3) measurements of diffuse emissions of CO 2 and other gases through soils. Various methods and instruments may be useful both for baseline studies and during unrest.

Scientific Investigations Report↗
USGS2025

Structural controls on splay fault rupture dynamics during Cascadia megathrust earthquakes

Great subduction earthquakes ( M w ≥ 8.0) can generate devastating tsunamis by rapidly displacing the seafloor and overlying water column. These potentially tsunamigenic seafloor offsets result from coseismic fault slip and deformation beneath or within the accretionary wedge. The mechanics of these shallow rupture phenomena and their dependence on subduction zone properties remain unresolved, partly due to the sparsity of offshore observations of shallow megathrust earthquake deformation. Here, we analyze how offshore structure influences shallow rupture mechanics and slip partitioning using 3D dynamic earthquake simulations of the Cascadia subduction zone (CSZ) megathrust with and without variably dipping seaward- or landward-vergent splay faults in the wedge that sole into the megathrust. Resulting tradeoffs between splay and megathrust slip reveal structural controls on rupture partitioning, with greater splay slip leading to less shallow megathrust slip updip. Gently dipping and seaward-vergent splays host more slip than those with steeper, landward-vergent splays. To isolate the underlying mechanisms, we compare models with Andersonian and plunging principal stresses. Results suggest distinct static and dynamic processes control the dip- and vergence-dependence of splay rupture: static (mis)alignment relative to far-field tectonic loading favors slip on more optimally oriented, shallowly dipping splay faults. In contrast, dynamic stress interactions of an updip-propagating megathrust rupture front with the free surface and potential branch faults favor forward branching onto seaward-vergent splays and inhibit backward branching onto landward-vergent splays. Resulting seafloor displacements suggest splay fault structure may influence coseismic tsunami source processes, highlighting the importance of dynamically viable rupture scenarios in subduction hazard assessments.

Cascadia subduction zone↗
USGS2026

Integrating detrital magnetite geochemistry and (U-Th)/He chronometry as a sediment provenance tool in geologic and metallogenic terranes

Magnetite is ubiquitous in porphyry Cu systems and in sediment sourcing both barren and mineralized regions, with potential as an indicator mineral in concealed and coarsely-mapped terranes. We develop and test a workflow for integrated geochemistry and (U–Th)/He (He) dating for inferring detrital magnetite (DMt) provenance in these settings. The ca. 70 Ma Taurus porphyry Cu–Mo(–Au) district in eastern interior Alaska serves as a test case. DMt from streams draining porphyry-related mineralization was characterized by geochemistry and mineral inclusion and microstructure observations, complemented by similar data for potential porphyry and host rock sources. Principal component analysis and clustering of DMt geochemical data resolve multiple populations in our samples geochemically and texturally compatible with derivation from metamorphic, porphyry-related hydrothermal, and igneous sources. Hydrothermal magnetite comprises ∼16–50% of DMt nearest porphyry mineralization but diminishes to ∼4% ∼15 km downstream. Subsampled grains within populations yield ∼160–110 Ma, ∼70 Ma, ∼55 Ma, and ∼20 Ma magnetite He date modes. Combined with provenance, He dates capture Early Cretaceous regional exhumation of metamorphic host rock and Late Cretaceous porphyry Cu mineralization. DMt grains showing partial hematite replacement yield ca. 55–20 Ma dates regardless of source, overlapping regional warm/wet climatic intervals. We interpret Cenozoic dates to reflect exhumation to near-surface oxidizing conditions and(or) supergene weathering. Magnetite is thus a promising target phase for (1) tracking the spatiotemporal distribution of porphyry systems, and (2) linking the formation and exhumation of these systems to a regional geologic history, both in Alaska and globally.

Alaska, Yukon↗
USGS2026

Quantifying methane emissions from a rich fen with uncrewed aircraft systems in boreal Alaska

Thawing of permafrost in northern latitudes is accelerating, potentially releasing substantial amounts of methane (CH 4 ) as forested permafrost plateaus transition into wetlands. This ecosystem shift alters the carbon exchange between the soil and atmosphere, influencing the permafrost-carbon feedback. Monitoring these changes may require measurement platforms operating across varied spatial and temporal scales. Recent advancements in small uncrewed aircraft systems (sUAS) enable high resolution CH 4 flux quantification in remote, complex terrains; however, comparisons with established methods such as eddy covariance flux towers remain limited. We used a hexacopter sUAS to quantify CH 4 emissions from the Alaska Peatland Experiment, a wetland within the Bonanza Creek Experimental Forest. Using an ensemble of methods to define the background CH 4 concentration, along with near surface emissions from soil chambers, helped constrain our flux estimates. The sUAS method yielded an average flux of 0.0077 ± 0.0019 mol s −1 CH 4 , within a factor of two concurrent tower-derived total source flux estimates (0.0036 ± 0.00042 mol s −1 CH 4 ). To assess spatial drivers of observed fluxes, we conducted a 2D footprint analysis and overlaid the results with high-resolution hyperspectral land cover classification, quantifying vegetative contributions within each footprint. This revealed higher fen representation in sUAS measurements (73.8%) than in tower footprints (58.8%), and lower tussock meadow representation (15.6% and 30.3%, respectively). These differences were consistent with known variation in vegetation-specific CH 4 emissions. Our results highlight that combining footprint modeling with land cover characterization can enhance interpretations of CH 4 fluxes and guide cross-platform comparisons.

JGR Atmospheres↗
USGS2025

A framework for understanding the effects of subsurface agricultural drainage on downstream flows

Understanding controls on streamflow volume and magnitude is important to water resource management applications, such as critical water and transportation structure design and floodplain mapping. Changes in land use and agricultural practices, such as subsurface agricultural drainage, may be contributing to changes in streamflow characteristics. Subsurface agricultural drainage, also known as tile drainage, is the practice of installing drains in the subsurface of agricultural fields to improve productivity. Because of the complex interactions between subsurface drainage systems, precipitation, local soil conditions, and land management practices, it is difficult to determine how subsurface agricultural drainage affects downstream flow. Previously developed subsurface agricultural drainage conceptual models under dry, saturated, and winter conditions are summarized, and current literature on the effects of subsurface agricultural drainage on downstream flows, focusing on peak flow, non-event flow, and total flow to develop frameworks for discussing these systems is compiled. The effects that subsurface drainage has on hydrologic systems are expected to vary by site and are seasonally based on system design, soil type, moisture conditions, precipitation characteristics, and land conditions. Subsurface drainage can affect the magnitude of peak flow by converting surface runoff from a storm event to subsurface runoff. By increasing hydrologic connectivity of a catchment, subsurface drainage can increase non-event flow or the flow between two storm events, typically dependent on lateral flow through the subsurface and groundwater. Theoretically, by diverting water from groundwater recharge or by reducing water available for evapotranspiration, subsurface drainage may increase the total volume of flow. Precipitation changes may increase infiltration, excess overland flow, and flood risk regardless of the presence or absence of subsurface drainage.

Illinois, Iowa, Michigan, Minnesota, Missouri, Mon↗
USGS2024

Using a time-of-travel sampling approach to quantify per- and polyfluoroalkyl substances (PFAS) stream loading and source inputs in a mixed-source, urban catchment

Understanding per- and polyfluoroalkyl substances (PFAS) mass distribution in surface and groundwater systems can support source prioritization, load reduction, and water management. Thirteen sites within an urban catchment were sampled utilizing a time-of-travel sampling approach to minimize the influence of subdaily fluctuations in mass from PFAS point sources and to quantify PFAS and ancillary chemical loads from various PFAS sources. A larger increase in perfluoroalkyl sulfonate (PFSA) loads (8 to 11 μg/s, up to 618%) than in perfluoroalkyl carboxylate (PFCA) loads (no change to 3.4 μg/s, up to 122%) was observed at sites below tributaries influenced by military bases with known groundwater discharge. Point discharges from two sewage treatment plants (STPs) resulted in increases in PFCA and PFSA loads that were similar (6 and 10 μg/s respectively) below the first STP and greater for PFCA compared to PFSA loads (23 and 13 μg/s respectively) below the second STP. Overall, percent increases in total PFAS load ranged from 20 to 277% for military base inputs and 44 to 77% for STP inputs. A focus catchment that represents only 14% (76.9 km 2 ) of the drainage area at the most downstream site (544 km 2 ) accounted for about 70% of PFSA and 40% of PFCA loads observed at the most downstream site. Results show that by using a time-of-travel sampling approach in mixed, urban settings with several PFAS sources, it is possible to quantify stream loads from individual PFAS sources, thereby improving source attribution and providing actionable data for water-resource managers.

ACS ES&T Water↗
USGS2025

High frequency and region-scale simulations of large (Mw7+) earthquakes on the southern Whidbey Island fault, Washington, USA

We simulate ground shaking in western Washington State from hypothetical M w 7.0–7.5 earthquakes on the southern Whidbey Island fault (SWIF). Ground motions are modeled considering kinematic source distributions on a complex fault plane, a 3D seismic velocity model, and region‐specific soil velocity models. We run simulations with varying model resolutions, including regional‐scale simulations with a maximum‐modeled frequency of ∼1 Hz and local‐scale simulations with a maximum‐modeled frequency of ∼2.5 Hz. Additional local‐scale simulations are run considering high‐resolution surface topography. We explore how source parameters (i.e., magnitude, hypocenter location, and dip direction) and 3D velocity structure impact peak shaking intensity and its variability. In particular, we find that earthquakes on the SWIF would likely produce strong shaking throughout the populated Puget Lowland, including in the cities of Everett, Seattle, Bellevue, and Tacoma, Washington. Simulated short‐period ( T ≤ 2 s) spectral accelerations are strong throughout the Puget Lowland, and long‐period shaking ( T ≥ 5 s) is strong in the deep regional sedimentary basins, especially the Everett and Seattle basins. Source parameters strongly influence intra‐ and interevent variability in response, primarily through changes in source and site geometry, as well as rupture directivity. We also note a potential coupling between rupture directivity and basin effects, wherein directivity pulses are seemingly guided through the region’s deep, interconnected sedimentary basins. Overall, this work highlights the impacts of 3D source, path, and site effects on seismic hazard in the U.S. Pacific Northwest and substantially expands the catalog of simulated ground motions for Puget Sound area crustal faults.

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

Land-based nutrient flux to a fringing reef: Insights from Ofu Island, American Samoa

Submarine groundwater discharge (SGD) is a critical driver of nutrient transport in coral reef ecosystems, shaping water quality, primary productivity, and overall reef health. This study quantifies SGD fluxes and associated nutrient dynamics in two reef flat pools within the Ofu Unit of the National Park of American Samoa: Papaloloa and Fatuana. A multi-method approach integrating unoccupied aerial system-based thermal infrared (UAS-TIR) surveys, radon-based SGD measurements, multichannel electrical resistivity tomography (ERT), and discrete water sampling was used to assess SGD rates and nutrient contributions. UAS-TIR imagery revealed cooler sea surface temperatures in both pools, indicative of SGD, with the higher fluxes observed in Papaloloa. Radon measurements revealed a strong inverse correlation between SGD rates and tidal stage, with a more immediate SGD response at Papaloloa due to its highly permeable calcareous sand and gravel substrate. In contrast, a 2–3-hour lag in SGD response at Fatuana suggests discharge from a more inland aquifer that has lower diffusivity. Nutrient concentrations correlated with temperature and salinity, confirming SGD as the dominant nutrient transport mechanism, whereas isotopic analyses indicated inputs from both groundwater and potential anthropogenic sources. Despite lower SGD flux at Fatuana, higher algal cover suggests additional factors influencing algal proliferation, including substrate availability and hydrodynamic conditions. Excess nutrient inputs from SGD may contribute to algal overgrowth, which threatens Ofu’s thermally tolerant corals by increasing competition for space and light. These findings underscore the complexity of SGD-mediated nutrient dynamics in reef environments and emphasize the need for integrated hydrological and ecological assessments to support effective reef conservation and management strategies.

American Samoa, Ofu Island, Olosega Island↗
USGS2026

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

Computation of detailed groundwater flow budgets for subdivisions of the Virginia Coastal Plain aquifer system has enabled quantification and more thorough understanding of groundwater flow within this important water resource. A zone budget analysis based on previously published groundwater models of the Virginia Coastal Plain and Virginia Eastern Shore indicates that groundwater conditions vary substantially throughout the Coastal Plain aquifer system because of 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 altered groundwater flow from predevelopment conditions. Rates of sustainable withdrawal are limited because the downward groundwater flow rate into confined aquifers is a relatively small part of the total groundwater budget for the aquifer system compared to the rate of recharge at the land surface. Analyses of groundwater budgets from the Virginia Coastal Plain model indicate that groundwater flow is generally outward from the surficial aquifer to rivers and coastal waterbodies 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 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 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 2023, rates of withdrawal substantially lower than the historical maximum 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 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. For the most recent simulated conditions (2023) downward groundwater flow continues, but because vertical flow rates are a function of the difference between water pressure in the upper surficial systems and lower confined units, rates of downward flow 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. Groundwater inflow from coastal regions could be contributing to saltwater intrusion, even though that was not measured in this study. Analyses of groundwater budgets from the Virginia Eastern Shore peninsula, a geographic region of the Virginia Coastal Plain, indicate 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 more than 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. For the most recent simulated period (2023), estimated storage loss accounts for slightly under 7 percent of withdrawals from the confined aquifer system. The reported withdrawal rate for this period from the confined Yorktown-Eastover system is near the highest reported rate for the Virginia 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 low rates of inflow from coastal water bodies underscore ongoing concerns about up-coning and lateral intrusion of salty groundwater.

Virginia↗
USGS2026

Satellite time series analysis to quantify changing climax ciénegas using a state and transition model approach

Ciénegas are rare wetlands in arid landscapes of the North American Southwest, historically providing critical ecological and hydrological functions but increasingly threatened by changing climate and land use pressures. This study quantifies changes in ciénega condition and floodplain dynamics using a state-and-transition model (STM) informed by expert knowledge and remote sensing. Key factors include woody plant encroachment, water availability, and soil aggradation. We mapped 31 ciénegas with high-resolution imagery and analyzed Landsat data (1985–2023) to assess vegetation health and moisture using the Normalized Difference Vegetation Index (NDVI) and Normalized Difference Infrared Index (NDII). Results show substantial interannual variability in phenology, water stress, and soil moisture, with regional drying and elevation strongly influencing ciénega resilience. We classified ciénegas into three functional states—healthy, desiccated, and dormant—and mapped their 2023 condition. Trend analyses indicate most ciénegas exhibit greening despite drought, though localized variability underscores the need for site-specific management. None are in a stable climax (reference) state; rather, they transition among states in response to external drivers. Increasing woody plant cover and surface drying, likely linked to declining regional water tables, favor deep-rooted species over wetland grasses—a pattern mirrored in adjacent control plots. Spatially explicit analysis revealed intra-ciénega variability often masked by aggregated data, highlighting the importance of high-resolution monitoring. Seasonal and long-term trends provide context for understanding ciénega dynamics, including degradation and restoration pathways. This study emphasizes the importance of groundwater conservation and demonstrates how remote sensing supports long-term monitoring. The STM framework offers a practical tool for adaptive management to sustain freshwater resources in arid environments.

Arizona, New Mexico↗
USGS2024

Bioconcentration of per- and polyfluoroalkyl substances and precursors in fathead minnow tissues environmentally exposed to aqueous film-forming foam-contaminated waters

Exposure to per- and polyfluoroalkyl substances (PFAS) has been associated with toxicity in wildlife and negative health effects in humans. Decades of fire training activity at Joint Base Cape Cod (MA, USA) incorporated the use of aqueous film-forming foam (AFFF), which resulted in long-term PFAS contamination of sediments, groundwater, and hydrologically connected surface waters. To explore the bioconcentration potential of PFAS in complex environmental mixtures, a mobile laboratory was established to evaluate the bioconcentration of PFAS from AFFF-impacted groundwater by flow-through design. Fathead minnows ( n = 24) were exposed to PFAS in groundwater over a 21-day period and tissue-specific PFAS burdens in liver, kidney, and gonad were derived at three different time points. The ∑PFAS concentrations in groundwater increased from approximately 10,000 ng/L at day 1 to 36,000 ng/L at day 21. The relative abundance of PFAS in liver, kidney, and gonad shifted temporally from majority perfluoroalkyl sulfonamides (FASAs) to perfluoroalkyl sulfonates (PFSAs). By day 21, mean ∑PFAS concentrations in tissues displayed a predominance in the order of liver > kidney > gonad. Generally, bioconcentration factors (BCFs) for FASAs, perfluoroalkyl carboxylates (PFCAs), and fluorotelomer sulfonates (FTS) increased with degree of fluorinated carbon chain length, but this was not evident for PFSAs. Perfluorooctane sulfonamide (FOSA) displayed the highest mean BCF (8700 L/kg) in day 21 kidney. Suspect screening results revealed the presence of several perfluoroalkyl sulfinate and FASA compounds present in groundwater and in liver for which pseudo-bioconcentration factors are also reported. The bioconcentration observed for precursor compounds and PFSA derivatives detected suggests alternative pathways for terminal PFAS exposure in aquatic wildlife and humans.

Environmental Toxicology and Chemistry↗
USGS2025

Wet meadow regeneration through restoration of biophysical feedbacks

Wet meadows are globally significant ecosystems that provide critical hydrological, ecological, and biogeochemical functions, yet their extent has declined dramatically due to land use changes and hydrologic alteration. These sedge-dominated wetlands exist at the drier end of the wetland gradient, maintained by shallow groundwater and periodic inundation. This paper is a global synthesis of the ecological, geomorphic, and hydrological dynamics of wet meadows, with an emphasis on alluvial systems, to inform effective restoration strategies. We compare wet meadows to other wetlands, classify them into palustrine, lacustrine, and alluvial types, then focus on alluvial wet meadows and discuss how their formation and persistence depend on ground and surface water interactions, sediment deposition and flow obstructions, all mediated by biological processes. In particular, we highlight the role of hydric graminoids in resisting erosion and maintaining soil cohesion, how beaver promote meadow persistence, and the significance of wet meadows as carbon sinks. We also present stratigraphic evidence demonstrating that incision, often triggered by anthropogenic activity or changing climate, is the primary mechanism of alluvial wet meadow degradation, resulting in water table decline and shifts in vegetation composition. Restoration requires reversing these incisional processes through techniques that elevate water tables, disperse flow and retain sediment—methods traditionally associated with either soil conservation or stream restoration. These include nature-based solutions that create obstructions such as beaver dams and their analogues, rock and wood-based obstructions and incision trench or gully filling and grading. Given their multifunctional value—including but not limited to flood attenuation, biodiversity support, and carbon sequestration—wet meadows warrant a focused restoration framework. This review advocates for a valley-floor scale restoration paradigm that integrates hydrological reconnection, sediment retention, and biological reinforcement to ensure long-term resilience of these systems in the face of changing climate and land use pressures.

Frontiers in Environmental Science↗