Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Hydrology and Earth Systems Science”

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

Recent scientific contributions by the U.S. Geological Survey in the San Francisco Bay and Sacramento–San Joaquin Delta Estuary

Introduction The San Francisco Bay and Sacramento–San Joaquin Delta Estuary (hereafter, Bay-Delta) is the largest estuary on the West Coast of the United States. The Bay-Delta covers more than 1,600 square miles and drains a watershed of more than 75,000 square miles, which is greater than 40 percent of California. The region surrounding the Bay-Delta is home to about 10 million people, and its habitats (fig. 1) support more than 800 plant and animal species. The waterways of the Bay-Delta are the central hub of California’s extensive freshwater delivery system, supplying water to more than 27 million Californians and 4 million acres of farmland in the Central Valley. The U.S. Geological Survey (USGS) is the primary science agency of the U.S. Department of the Interior, providing a broad range of Earth, water, biological, and mapping data and expertise to inform natural resource management across the country. This fact sheet focuses on research conducted by the USGS in the Bay-Delta region, mostly within the past 5 years. The fact sheet is organized across five major themes: water flow and water quality, fish and wildlife, wetland restoration, invasive species, and hazards. In each of these areas, the USGS works closely with Federal, State, and local agencies; academia; and non-governmental organizations to provide objective, evidence-based science. The data and knowledge gained from USGS research and monitoring are publicly available and are routinely used by agencies, including the U.S. Bureau of Reclamation, U.S. Army Corps of Engineers, U.S. Fish and Wildlife Service, California Department of Water Resources, California Department of Fish and Wildlife, California State Water Resources Control Board, Delta Stewardship Council, and many additional organizations and stakeholders to inform management of the Bay-Delta region and its resources.

California↗
USGS2025

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

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

Current Climate Change Reports↗
USGS1998

Continuous hydrologic simulation of runoff for the Middle Fork and South Fork of the Beargrass Creek basin in Jefferson County, Kentucky

The Hydrological Simulation Pro-gram-FORTRAN (HSPF) was applied to an urban drainage basin in Jefferson County, Ky to integrate the large amounts of information being collected on water quantity and quality into an analytical framework that could be used as a management and planning tool. Hydrologic response units were developed using geographic data and a K-means analysis to characterize important hydrologic and physical factors in the basin. The Hydrological Simulation Program FORTRAN Expert System (HSPEXP) was used to calibrate the model parameters for the Middle Fork Beargrass Creek Basin for 3 years (June 1, 1991, to May 31, 1994) of 5-minute streamflow and precipitation time series, and 3 years of hourly pan-evaporation time series. The calibrated model parameters were applied to the South Fork Beargrass Creek Basin for confirmation. The model confirmation results indicated that the model simulated the system within acceptable tolerances. The coefficient of determination and coefficient of model-fit efficiency between simulated and observed daily flows were 0.91 and 0.82, respectively, for model calibration and 0.88 and 0.77, respectively, for model confirmation. The model is most sensitive to estimates of the area of effective impervious land in the basin; the spatial distribution of rain-fall; and the lower-zone evapotranspiration, lower-zone nominal storage, and infiltration-capacity parameters during recession and low-flow periods. The error contribution from these sources varies with season and antecedent conditions.

Kentucky↗
USGS2026

Birds, breakpoints, and baselines: How citizen science data can reveal ecological boundaries in Kenya’s Upper Tana watershed

Tropical watersheds are increasingly threatened by climate change, land-use conversion, and resource extraction, yet conventional biodiversity monitoring in these systems is often spatially and temporally limited. Citizen science offers a complementary approach, enabling biodiversity data collection over large areas that can supplement professional scientific surveys. We analyzed 10 years (2012–2022) of Upper Tana Watershed bird data from the Kenya Bird Map project, covering 114 pentads (9 × 9 km) within a 17,000 km 2 watershed to assess patterns of bird community composition and distribution across this watershed and to also evaluate the effects of environmental variables, seasonality, and sampling effort to help inform improvements in future citizen science projects. Citizen (or community) scientists recorded 575 species (>50% of Kenya’s total avifauna) in 74 families. Asymptotic species accumulation indicates that most probable species present in the watershed were detected. Threshold indicator taxa analysis revealed distinct ecological boundaries along elevation (∼1,500 m), precipitation (∼1,100 mm), and mean temperature (∼19°C) gradients, corresponding to a turnover from xeric savanna to mesic montane forest assemblages. Notably, bird communities showed little seasonal differentiation between wet and dry periods, consistent with dominance by resident year-round species. Data limitations including uneven survey distribution and frequency, absence of abundance metrics, and coarse representation of local environmental conditions that likely reduced our ability to detect fine-scale species–habitat relationships. Addressing these gaps through spatially balanced sampling at greater resolution, greater survey frequency in underrepresented areas, and improved capture of habitat metrics could strengthen the use of citizen-science bird data for watershed bioassessment. Our findings demonstrate that structured citizen-science initiatives can identify ecological boundaries and inform adaptive management of tropical socio-ecological systems under rapid environmental change.

Upper Tana watershed↗
USGS2024

Aboveground carbon stocks across a hydrological gradient: Ghost forests to non-tidal freshwater forested wetlands

Upper estuarine forested wetlands (UEFWs) play an important role in the sequestration of atmospheric carbon (C), which is facilitated by their position at the boundary of terrestrial and maritime environments but threatened by sea level rise. This study assessed the change in aboveground C stocks along the estuarine–riverine hydrogeomorphic gradient spanning salt-impacted freshwater tidal forested wetlands to freshwater forested wetlands in seasonally tidal and nontidal landscape positions. Standing stocks of C in forested wetlands were measured along two major coastal river systems, the Winyah Bay in South Carolina and the Savannah River in Georgia (USA), replicating and expanding a previous study to allow the assessment of change over time. Aboveground C stocks on these systems averaged 172.9 Mg C ha −1 , comparable to those found in UEFWs across the globe and distinct from the terrestrial forested ecosystems they are often considered to be a part of during large-scale C inventory efforts. Groundwater salinity conditions as low as 1.3 ppt were observed in conjunction with losses of aboveground C. When viewed in context alongside expected sea level rise and corresponding saltwater intrusion estimates, these data suggest a marked decrease in aboveground C stocks in forested wetlands situated in and around tidal estuaries.

Georgia, South Carolina↗
USGS2024

Repeated coseismic uplift of coastal lagoons above the Patton Bay Splay Fault System, Montague Island, Alaska, USA

Coseismic slip on the Patton Bay splay fault system during the 1964 M w 9.2 Great Alaska Earthquake contributed to local tsunami generation and vertically uplifted shorelines as much as 11 m on Montague Island in Prince William Sound (PWS). Sudden uplift of 3.7–4.3 m caused coastal lagoons along the island's northwestern coast to gradually drain. The resulting change in depositional environment from marine lagoon to freshwater muskeg created a sharp, laterally continuous stratigraphic contact between silt and overlying peat. Here, we characterize the geomorphology, sedimentology, and diatom ecology across the 1964 earthquake contact and three similar prehistoric contacts within the stratigraphy of the Hidden Lagoons locality. We find that the contacts signal instances of abrupt coastal uplift that, within error, overlap the timing of independently constrained megathrust earthquakes in PWS—1964 Common Era, 760–870 yr BP, 2500–2700 yr BP, and 4120–4500 yr BP. Changes in fossil diatom assemblages across the inferred prehistoric earthquake contacts reflect ecological shifts consistent with repeated draining of a lagoon system caused by >3 m of coseismic uplift. Our observations provide evidence for four instances of combined megathrust-splay fault ruptures that have occurred in the past ∼4,200 years in PWS. The possibility that 1964-style combined megathrust-splay fault ruptures may have repeated in the past warrants their consideration in future seismic and tsunami hazards assessments.

Alaska↗
USGS2024

Correlation analysis of groundwater and hydrologic data, Kaloko-Honokōhau National Historical Park, Hawai‘i

Designated in 1978, Kaloko-Honokōhau National Historical Park is located on the west coast of the Island of Hawaiʻi. The Kaloko-Honokōhau National Historical Park encompasses about 1,200 acres of coastal land and nearshore ecosystems, which include wetlands, anchialine pools (landlocked bodies of brackish water with hydrologic connections to the ocean), fishponds, a fishtrap, and coral reefs. These nearshore ecosystems are dependent on groundwater discharge with a freshwater component and provide habitat for threatened and endangered, endemic species, such as the orangeblack Hawaiian damselfly ( Megalagrion xanthomelas ) and the Hawaiian coot (ʻAlae keʻokeʻo, Fulica alai ). The populations of these native species, however, are threatened because of habitat loss related to urban development and environmental changes. Kaloko-Honokōhau National Historical Park is within the Keauhou aquifer system and the North Kona District, which experienced a 52 percent resident-population increase between 2000 and 2020 and a 41 percent visitor increase between 2008 and 2019. To support the current water demand associated with this growing population, groundwater is the primary source of freshwater used in the North Kona District, with about 15 million gallons of groundwater withdrawn from the Keauhou aquifer system per day since 2009. With anticipated development, future (2015–35) groundwater withdrawal from the Keauhou aquifer system is projected to be about 55 percent greater than recent (2012–14) withdrawal. Because Kaloko-Honokōhau National Historical Park is located within a coastal aquifer, natural and human-induced changes can affect the quality and quantity of groundwater, which can threaten groundwater-dependent ecosystems. To improve understanding of recent groundwater conditions, the U.S. Geological Survey, in cooperation with the National Park Service, undertook this study to document correlations between hydrologic time-series datasets from sites in and near Kaloko-Honokōhau National Historical Park using the nonparametric (distribution-free) Kendall’s tau statistical test. For the statistical analyses, dependent variables representing the groundwater system include groundwater level, the groundwater-level difference between pairs of sites, and specific conductance, and independent variables include datasets of sea level, rainfall, and groundwater withdrawal. About 34 percent of the 140 non-time-lagged Kendall’s tau statistical tests evaluated in this report are statistically significant ( p -value ≤ 0.050) with generally weak (0.1 ≤ tau ≤ 0.2) to moderate (0.2 ≤ tau ≤ 0.3) correlations. Groundwater levels measured at monitoring sites have the strongest correlation with the multivariate El Niño–Southern Oscillation index and withdrawal from production wells at the nearby Kohanaiki Private Club Community. Specific conductance is not consistently and significantly correlated with the independent hydrologic variables investigated in this report. Because the relations between hydrologic variables are commonly not instantaneous, a second set of correlations was evaluated after applying a range of time lags to the independent variable datasets. Relative to the non-time-lagged case (the set of correlations that did not use time-lagged independent variables), some of the time-lagged independent variables improved correlations with some of the dependent variables. For a particular independent variable, similar time lags were expected between the independent variable and dependent variable at all four monitoring sites. However, different time lags among the four sites sometimes produced the strongest correlations. This study identified several correlations that are statistically significant and hydrologically plausible, but the correlations could indicate that multiple concurrent factors are controlling the observed groundwater-system response, which might be better addressed using multivariate analyses. This study only investigates bivariate correlations, which may not explain all the variance in the data. The correlations analyzed in this report are limited by the quantity of available hydrologic data in the area near Kaloko-Honokōhau National Historical Park and are based on 14 years of time-series data, which were aggregated to a relatively coarse monthly temporal resolution that represents the minimum resolution common to all datasets.

Hawaii↗
USGS2026

Cajon Pass and the southern San Andreas Fault System: Earthquake cycle stress accumulation and present-day loading

With over a century since the last major rupture affecting the wider Los Angeles region, tectonic stress has steadily built along the southern San Andreas and San Jacinto fault systems, raising concerns of an imminent large earthquake. Cajon Pass, located at the junction of these faults, represents a critical site for potential through-going ruptures in Southern California. We constructed new 4D earthquake cycle simulations using a 1000-year paleoseismic rupture history of the San Andreas Fault System (SAFS) to assess spatial and temporal variations in stress. A semi-analytic Fourier transform model was used to compute stress from 3D dislocations in an elastic plate overlying a Maxwell viscoelastic half-space, assuming a complete coseismic reset of resolved shear stress on ruptured elements. Results show highest stress accumulation north of Cajon Pass (∼1.8 MPa/100 years) due to greater slip rates, and lower rates south of Cajon Pass (∼1.0–1.5 MPa/100 years). By 2025, Coulomb stress is estimated at 2.8 MPa on the Mojave South (MOS) segment, 1.8 MPa on the North San Bernardino (NSB1) segment and 3.6 MPa on the San Jacinto Bernardino (SJB) segment. Segments accumulate stresses with characteristic ranges of pre-event stress interpreted as failure thresholds: 1.2–2.7 MPa for MOS, 0.4–1.6 MPa for NSB1, and 1.2–2.9 MPa for SJB. When the stress disparity between segments SJB and MOS narrows, the faults appear to rupture jointly, suggesting that stress levels may control how Cajon Pass acts as an earthquake gate. These results may inform seismic hazard assessments by linking stress evolution to fault interactions.

California↗
USGS2026

Development of a hydrologic flow model of the Zuñi Salt Lake and surrounding areas, west-central New Mexico

The terminal Zuñi Salt Lake is in a maar in west-​central New Mexico and contains hypersaline water that has long been used by Native Americans for religious purposes and the collection of salt. The U.S. Geological Survey (USGS), in cooperation with the Zuni Tribe of the Zuni Reservation, New Mexico, and Bureau of Reclamation, developed the Zuñi Salt Lake hydrologic flow model to simulate the steady-state conditions of the lake and surrounding groundwater-flow system. The model is a tool that can be used to analyze the potential hydrologic responses to resource development in the region under different water-use scenarios, which can support negotiations that could help protect future use of the lake. The USGS modular hydrologic model, MODFLOW 6, was used to simulate the hydrologic system of the lake and surrounding areas. To define the layering of the model, a three-dimensional hydrogeologic framework was constructed, which included seven informal hydrogeologic units and the locations of faults that truncate and offset the units. Calibration of the model was performed using the iterative ensemble smoother included in a parameter estimation software suite (PEST++). The iterative ensemble smoother approach resulted in a posterior parameter ensemble of model realizations with parameter and simulated values that show similar success in matching observations. Groundwater flux to the Zuñi Salt Lake simulated by using the posterior parameter distribution ranged from 26 to 730 acre-feet per year with a mean, median, and standard deviation of 530, 570, and 170 acre-feet per year, respectively. The relative contribution of groundwater discharge to the lake was simulated using a particle-tracking model for MODFLOW (MODPATH) to be from, in order of decreasing contribution, the Mesaverde, Cenozoic, Mancos, and Dakota hydrogeologic units. Simulations of future conditions using the ensemble of models can provide an empirical estimate of forecast uncertainty without substantial additional analysis.

New Mexico↗
USGS2025

Multi-scale geophysical imaging of a hydrothermal system in Yellowstone National Park, USA

Little is known about the local plumbing systems that fuel Yellowstone’s famous hot springs, geysers and mud pots. A multi-method, multi-scale geophysical investigation was carried out in the Obsidian Pool Thermal Area (OPTA) to: (i) delineate the lateral extent of the hydrothermal area and associated surface features; (ii) estimate the dimensions of the upflow zone and identify its main controlling structures; (iii) assess fluids circulation pathways from depth to surface. Ground and airborne geophysical data were acquired to connect local and regional scales, from shallow to large depths. Maps of surface electrical resistivity show a strong correlation with hydrothermal features. At in-termediate depths, electrical resistivity permits delineating the upper limit of the upflow zone, while Poisson’s ratio highlights differences in subsurface fluid content. Combining these results with surface observations and topographic information, we speculate that differential mixing of hydrothermal and fresh water could explain the wide diversity of features observed at OPTA. Low electrical resistivity observed at large depths also suggest that a vast upflow zone, controlled by rhyolite flows and conjugate faults, underlies the OPTA. We speculate that hydrothermal fluids rise along fractures and reach the surface in topographic lows to form hydrothermal features. Our results show that synoptic, multi-scale geophysical measurements provide a roadmap for understanding where and how geologic heterogeneity, topography, fluid-gas separation, and the mixing of thermal and meteoric waters conspire to produce the wide variety of Yellowstone’s renowned hydrothermal features.

Wyoming↗
USGS2026

The use of fluorite geochemistry and machine learning to identify critical mineral systems

Fluorite (CaF 2 ) is a potential pathfinder to critical mineral and rare earth element (REE) deposits but its application has been limited to a narrow range of mineralization types. I show that fluorite is a robust recorder of mineralization fertility by applying statistical and machine-learning methods to a new global fluorite geochemical database. Distinct median rare earth and trace element patterns are observed among deposit types and genetic environments. Fluorite associated with carbonatites and REE deposits are relatively enriched in Sr and have minimal Eu anomalies. These characteristics define new bivariate discrimination diagrams that correctly identify 78% of carbonatite-related fluorite and 88% of fluorite from REE deposits. Random forest classifiers were developed for a wide range of mineralization types and genetic settings. Trained solely on rare earth element patterns, these models achieve accuracies of 77–79%. Higher classification accuracies (up to 88–96%) are obtained when including elements such as Sr, highlighting the significance of trace elements for optimal fluorite classification. The recognition of diagnostic fluorite compositional fingerprints, particularly in REE-fertile systems, underscores its potential as a pathfinder and indicator for critical mineral exploration in F-bearing environments.

Mineralium Deposita↗
USGS2024

Factors contributing to pesticide contamination in riverine systems: The role of wastewater and landscape sources

Wastewater treatment plant (WWTP) discharges can be a source of organic contaminants, including pesticides, to rivers. An integrated model was developed for the Potomac River watershed (PRW) to determine the amount of accumulated wastewater percentage of streamflow (ACCWW) and calculate predicted environmental concentrations (PECs) for 14 pesticides in non-tidal National Hydrography Dataset Plus Version 2.1 stream segments. Predicted environmental concentrations were compared to measured environmental concentrations (MECs) from 32 stream sites that represented a range of ACCWW and land use to evaluate model performance and to assess possible non-WWTP loading sources. Statistical agreement between PECs and MECs was strongest for insecticides, followed by fungicides and herbicides. Principal component analysis utilizing optical fluorescence and ancillary water quality data identified wastewater and urban runoff sources. Pesticides that indicated relatively larger sources from WWTPs included dinotefuran, fipronil, carbendazim, thiabendazole, and prometon whereas imidacloprid, azoxystrobin, propiconazole, tebuconazole, and diuron were more related to urban runoff. In addition, PECs generally comprised a low proportion of MECs, which indicates possible dominant loading sources beyond WWTP discharges. Cumulative potential toxicity was higher for sites with greater ACCWW and/or located in developed areas. Imidacloprid, fipronil, and carbendazim accounted for the largest portion of predicted potential toxicity across sites. The chronic aquatic life toxicity benchmarks for freshwater invertebrates were exceeded for 82 % of the imidacloprid detections ( n = 28) and 47 % of the fipronil detections ( n = 19). These results highlight the ecological implications of pesticide contamination from WWTP discharges and also the potential legacy effects from accumulated soil and groundwater sources. Pesticide management strategies that mitigate both current and historical impacts may improve the health of aquatic ecosystems.

Potomac River watershed↗
USGS2025

Water supply in the conterminous United States, Alaska, Hawaii, and Puerto Rico, water years 2010–20

We present an assessment of water supply across the conterminous United States (CONUS), Alaska, Hawaii, and Puerto Rico covering water years 2010–20. Our analysis drew on two national hydrologic models, the National Hydrologic Model Precipitation-Runoff Modeling System and the Weather Research and Forecasting model hydrologic modeling system. Both models produced estimates of streamflow, evapotranspiration, soil moisture, snow water equivalent, and other hydrologic states and fluxes. The models were driven by the bias-adjusted 4-kilometer-resolution, long-term regional hydroclimate simulation over the conterminous United States dataset (CONUS404). We assessed spatial and temporal error distributions by comparing monthly simulations at the 12-digit hydrologic unit code and regional scale from both models against external benchmarking datasets. Results showed that average annual rainfall across the CONUS was 857 millimeters per year for the period of analysis, with water year 2012 the driest year (729 millimeters) and water year 2019 the wettest year (995 millimeters). Key interannual variability results included the following: (1) the California–Nevada hydrologic region had the highest variability in precipitation and snow accumulation, and (2) the Texas hydrologic region was among hydrologic regions with the highest variability in precipitation. We related interannual variability in precipitation to storage volumes in soil moisture, snow water equivalent, and lakes and reservoirs to highlight areas with little storage and large year-to-year variability in precipitation. These areas included the Southern High Plains, Central High Plains, Texas, Souris–Red–Rainy, Mississippi Embayment, and Midwest regions. Our analysis of groundwater-level data showed that several of these areas overlap aquifers where groundwater levels were considerably lower than historical averages, including the Colorado Plateaus aquifers, the Rio Grande aquifer system, and the Central and Southern regions of the High Plains aquifer. Many of these lowered groundwater levels are continuations of decades-long declines from overpumping that started well before the assessment period. The resulting water budgets and their analyses provide a high-resolution foundational assessment of the mean state and variability of the terrestrial hydrologic cycle across the CONUS and Alaska, Hawaii, and Puerto Rico to support a wide range of water resource management applications.

Professional Paper↗
USGS2025

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

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

Massachusetts↗
USGS2024

Thermo-hydrologic processes governing supra-permafrost talik dynamics in discontinuous permafrost near Umiujaq (Québec, Canada)

Widespread supra-permafrost talik formation is currently recognized as a critical mechanism that could accelerate permafrost thaw in the Arctic (e.g., Connon et al. 2018; Farquharson et al. 2022). However, the trajectory of permafrost dynamics following talik formation may prove difficult to predict. Physically-based cryohydrogeologic models provide a powerful tool for understanding processes and factors controlling talik dynamics and, ultimately, how permafrost will respond to climate change. Such models are typically used to represent multiple non-linear processes relevant for groundwater systems in cold regions, such as coupled heat and groundwater movement, including freeze-thaw dynamics and the effects on the surface energy balance and the subsurface thermal and hydraulic properties (Lamontagne-Hallé et al. 2020). Though cryohydrogeologic modeling advances have been made in simulating talik dynamics, few applications have been tested against robust long-term hydrometeorological and subsurface observations.

Quebec↗
USGS2025

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

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

Colorado↗
USGS2021

Remote sensing-based actual evapotranspiration assessment in a data-scarce area of Brazil: A case study of the Urucuia Aquifer System

The large groundwater reserves of the Urucuia Aquifer System (UAS) enabled agricultural development and economic growth in the western Bahia State, in northeastern Brazil. Over the last several years, concern has grown around the aquifer’s diminishing water levels, and water balance (WB) studies are in demand. Considering the lack of measured actual evapotranspiration (ET a ), a major component of the water cycle, this work uses the Operational Simplified Surface Energy Balance (SSEBop) model to estimate ET a , and compares it to basin-scale estimates from the Soil Moisture Accounting Procedure (SMAP) monthly model and from an annual WB closure method, based on gridded meteorological data and the Gravity Recovery and Climate Experiment (GRACE) product. Additionally, a comparative assessment of different versions of the SSEBop parameterization was performed. Moderate Resolution Imaging Spectroradiometer (MODIS) imagery was used to implement eight different versions of the SSEBop algorithm over the UAS between 2000 and 2013. SSEBop and SMAP ET a yielded similar seasonal patterns, with correlation coefficient (r) up to 0.65, mean difference (MD) of 0.8 mm/month and mean absolute difference (MAD) of 18.5 mm/month. Comparison of SSEBop annual ET a estimates to annual SMAP and WB closure estimates yielded low MD (12.1 and −7.3 mm/year, respectively) and MAD (82.5 and 82.8 mm/year, respectively), but also low r values (0.00 and 0.37, respectively). The comparison of the different SSEBop versions indicated the need to incorporate a calibration step of the aerodynamic heat resistance (r ah ) parameter. SSEBop results were also used for land cover and drought monitoring. Analysis indicates that agriculture, associated with an increasing trend of atmospheric evaporative demand, is responsible for the decrease in groundwater levels and streamflow in the studied time period.

Urucuia Aquifer System↗
USGS2025

Fiber-optic distributed temperature sensing of hydrologic processes—Diverse deployments and new applications by the U.S. Geological Survey

Fiber-optic distributed temperature sensing instruments harness the temperature-dependent properties of glass to measure temperature continuously along optical fibers by using precise pulses of laser light. In the mid-2000s, this technology was refined for environmental monitoring purposes such as snowpack-air exchange, groundwater/surface-water exchange, and lake-water stratification. Fiber-optic distributed temperature sensing has revealed unprecedented details about preferential flow processes; however, the method is labor intensive and requires specific training, resulting in limited use by the broader water community. With the ongoing national implementation of the U.S. Geological Survey Next Generation Water Observing System, there has been renewed interest in harnessing the unique spatiotemporal monitoring capabilities of fiber-optic distributed temperature sensing. This fact sheet briefly describes this technology, highlights uses by the U.S. Geological Survey, and discusses current applications and future opportunities.

Fact Sheet↗