Geology Reports⌕ Search

FIND YOUR NEXT DISCOVERY

Results for “Hydrological Processes”

Original records, connected by a shared subject.

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

392 records · Page 2Linked to original sources
USGS2026

A targeted approach for mapping groundwater discharge to surface water and fish thermal refuge in four Lake Ontario tributaries

The duration, magnitude, and frequency of heatwaves are predicted to increase in the coming decades, a combination that can reduce the survival of many fish species. Across the world, there is broad interest in identifying thermal refuge for heat-intolerant fish species and exploring opportunities to enhance or protect these areas. Because deeper groundwater maintains a relatively constant temperature, groundwater-influenced areas along streams can provide cool-water refuge for fish during periods of extreme heat. A targeted approach was developed for identifying existing cold-water zones and areas of substantial groundwater discharge in four high priority Lake Ontario tributaries. Our approach included: (1) predicting where groundwater discharge is most likely with a simple geospatial model and (2) using model predictions to select field sites for intensive high-resolution study, including ground-based mapping of groundwater features (springs, seeps, tributaries) as well as drone-based optical and thermal infrared surveys. Results from field sites were used to both verify model performance and map different types and aerial extents of thermal anomalies. Geospatial modelling successfully predicted regions of widespread groundwater upwelling, later verified and mapped by field and drone surveys. Comparison of model and field survey results further highlighted specific geospatial layers, such as soil/bedrock types and topographic wetness index, as being particularly useful for predicting groundwater influence on streams in the study area. In addition, a comparison of geospatial model results with a model of fish abundances along the studied streams showed significant positive correlations for many heat-intolerant fish species over a wide geographic area. The approach developed in this study can be applied to other watersheds to highlight areas of probable groundwater discharge and could be used by fishery and water resource managers to support cold-water fish habitat management decision-making and resource conservation.

New York↗
USGS2025

ARCHI: A new R package for automated imputation of regionally correlated hydrologic records

Missing data in hydrological records can limit resource assessment, process understanding, and predictive modeling. Here, we present ARCHI (Automated Regional Correlation Analysis for Hydrologic Record Imputation), a new, open-source software package in R designed to aggregate, impute, cluster, and visualize regionally correlated hydrologic records. ARCHI imputes missing data in “target” records by linear regression using more complete “reference” records as predictors. Automated imputation is implemented using a novel, iterative algorithm that allows each site to be considered a target or reference for regression, growing the pool of complete references with each imputed record until viable gap-filling ceases. Users can limit artifacts from spurious correlations by specifying model-acceptance criteria and applying geospatial, correlation, and group-based filters to control reference selection. ARCHI provides additional functions for visualizing results, clustering records with similar correlation structures, evaluating holdout data, and interactive parameterization with an accessible and intuitive graphical user interface (GUI). This methods brief provides an overview of the ARCHI package, modeling guidelines, and benchmarking on two regional groundwater-level datasets from the Central Valley, CA and Long Island, NY. We evaluate ARCHI alongside widely used multivariate imputation software to highlight and contextualize its computational efficiency, imputation accuracy, and model transparency when applied to large, groundwater-level datasets.

California, New York↗
USGS2025

Hydrologic response of groundwater and streamflow to natural and anthropogenic drivers of change in headwaters of the upper Colorado River basin during recent wet (1982–1999) and drought (2000–2022) conditions

Study region: Headwaters of the upper Colorado River basin (UCOL), USA Study focus: Surface-water and groundwater numerical models incorporating water-use information were used to investigate changes in climate, water use, and simulated hydrologic responses of snow processes, evapotranspiration, groundwater, and streamflow during recent wet (1982–1999) and drought (2000–2022) periods in the headwater subregions of the upper Colorado River basin. New hydrologic insights for the region: Decreases in average streamflow between wet and drought periods ranged from 20 % in the Colorado River headwaters subregion to 23 % in the Gunnison River headwaters subregion. Like streamflow, average surface runoff was statistically less during the drought than the wet period, with decreases from 24–31 % in the headwaters. On a volume basis, runoff decreases were greater than streamflow decreases in both the Colorado River and Gunnison River headwaters. Although the amount of water-year groundwater discharge to streams remained nearly the same between the wet and drought periods, groundwater as a percentage of streamflow increased between the wet and drought periods, highlighting the importance of groundwater in sustaining streamflow during drought conditions. Multiple linear regression analyses revealed that snowmelt-only models were better than the best precipitation and temperature models at explaining streamflow variability from all headwater subregions for both the wet and drought periods.

Arizona, Colorado, New Mexico, Utah, Wyoming↗
USGS2000

Biodegradation of chlorinated ethenes at a karst site in middle Tennessee

This report presents results of field and laboratory investigations examining the biodegradation of chlorinated ethenes in a karst aquifer contaminated with trichloroethylene (TCE). The study site, located in Middle Tennessee, was selected because of the presence of TCE degradation byproducts in the karst aquifer and available site hydrologic and chlorinated-ethene information. Additional chemical, biological, and hydrologic data were gathered to evaluate whether the occurrence of TCE degradation byproducts in the karst aquifer was the result of biodegradation within the aquifer or simply transport into the aquifer. Geochemical analysis established that sulfate-reducing conditions, essential for reductive dechlorination of chlorinated solvents, existed in parts of the contaminated karst aquifer. Other areas of the aquifer fluctuated between anaerobic and aerobic conditions and contained compounds associated with cometabolism, such as ethane, methane, ammonia, and dissolved oxygen. A large, diverse bacteria population inhabits the contaminated aquifer. Bacteria known to biodegrade TCE and other chlorinated solvents, such as sulfate-reducers, methanotrophs, and ammonia-oxidizers, were identified from karst-aquifer water using the RNA-hybridization technique. Results from microcosms using raw karst-aquifer water found that aerobic cometabolism and anaerobic reductive-dechlorination degradation processes were possible when appropriate conditions were established in the microcosms. These chemical and biological results provide circumstantial evidence that several biodegradation processes are active in the aquifer. Additional site hydrologic information was developed to determine if appropriate conditions persist long enough in the karst aquifer for these biodegradation processes to be significant. Continuous monitoring devices placed in four wells during the spring of 1998 indicated that pH, specific conductance, dissolved oxygen, and oxidation-reduction potentials changed very little in areas isolated from active ground-water flow paths. These stable areas in the karst aquifer had geochemical conditions and bacteria conducive to reductive dechlorination of chlorinated ethenes. Other areas of the karst aquifer were associated with active ground-water flow paths and fluctuated between anaerobic and aerobic conditions in response to rain events. Associated with this dynamic environment were bacteria and geochemical conditions conducive to cometabolism. In summary, multiple lines of evidence developed from chemical, biological, and hydrologic data demonstrate that a variety of biodegradation processes are active in this karst aquifer.

Tennessee↗
USGS2025

A new groundwater energy transport model for the MODFLOW hydrologic simulator

Heat transport in the subsurface is an important aspect of research related to the effects of a warming climate on ecological services (i.e., cold-water refugia); the development of geothermal resources for energy banking schemes (i.e., aquifer thermal energy storage [ATES]); and the effects of temperature on other aspects of groundwater quality, such as nutrient cycling. Historically, simulation of heat transport using the MODFLOW groundwater simulator and related codes was performed by scaling the input parameters of a solute-transport model to emulate heat transport. However, that approach required additional pre- and post-processing of input and output and could not account for the variation in effective thermal storage and transport properties during transient, unsaturated flow, for example. True heat-transport capabilities in the context of MODFLOW were first introduced in a variant called USG-Transport. More recently, a new groundwater energy-transport (GWE) model type has been added to MODFLOW 6, the core version of the MODFLOW hydrologic simulator. GWE supports the simulation of heat transport on structured or unstructured grids as well as within and between features of advanced packages that represent streams, lakes, multi-aquifer wells, and the unsaturated zone. GWE is integrated within MODFLOW 6 and is accessible through the FloPy Python package and the MODFLOW 6 application programming interface (API). An example simulation demonstrates conduction between grid cells through both the water and the solid aquifer material, including thermal bleeding from saturated overburden cells into a groundwater flow field.

Groundwater↗
USGS2025

Near-surface geophysics: Environmental applications

The field of geophysics encompasses a broad and diverse compilation of methodologies that employs principles of physics to characterize properties of earth materials within the subsurface. While geophysical methods have a long history in resource exploration and studies of Earth’s interior, the subdiscipline of “near-surface geophysics” has evolved in recent decades for examination of the shallow, near-surface environment for a range of purposes ranging from archaeological or forensic investigations to assessment of geologic, hydrologic, biologic, and geochemical properties and processes. “Environmental geophysics” are near-surface geophysical studies and methods that focus on understanding natural systems (e.g., watershed hydrology, groundwater–surface water connections, biophysical processes) as well as research pertaining to anthropogenic impacts and land management, (e.g., contamination and remediation, saltwater intrusion, agricultural practices). This field can be further subdivided into subdisciplines focused on specific topics and applications, such as water resources and hydrology (hydrogeophysics) or biologic and microbial processes (biogeophysics). Studies in environmental geophysics span a range of scales, from pore-scale laboratory tests to watershed-scale or regional field experiments. Methods vary by the nature of physics employed, the specific measurement acquired, and how that data is ultimately processed and analyzed to produce interpretable results. There exists further diversity in the acquisition logistics, geometry, and timing of data collection. Geophysical data can be collected in boreholes (one-dimensional, 1-D, vertical profiles), along survey lines (two-dimensional, 2-D, cross-sections), or in dense sensor arrays or gridded profiles (three-dimensional, 3-D, models). Regarding the temporal aspect, studies can conduct one-time geophysical surveys to obtain detailed imaging of subsurface structure or use timelapse and continuous monitoring to investigate variations in subsurface properties over time. The cumulation of all possible permutations of these factors (method, acquisition geometry, survey design, and target application) results in an immense diversity among environmental geophysical studies. Nevertheless, this field remains unified in the pursuit of understanding natural and human-impacted near-surface environments through geophysical investigations. Here we highlight some key references within environmental geophysics. Resources on geophysical theory, acquisition logistics, processing and inversion workflows, and example case studies are categorized into the most common geophysical classes within Geophysical Methods. Lastly, example references for the dominant types of applications in environmental geophysical studies are catalogued in Environmental Applications.

Book chapter↗
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↗
USGS2025

Forecasting water levels using the ConvLSTM algorithm in the Everglades, USA

Forecasting water levels in complex ecosystems like wetlands can support effective water resource management, ecological conservation, and understanding surface and groundwater hydrology. Predictive models can be used to simulate the complex interactions among natural processes, hydrometeorological factors, and human activities. The Greater Everglades in the USA is a well-known example of an ecosystem where complexity has motivated adoption of machine learning algorithms in water level prediction studies. This paper aims to contribute to extending existing machine learning algorithms by integrating spatiotemporal data with deep-learning algorithms in the forecasting process. In this study, a deep-learning model is developed to predict water levels on a regional scale, covering a large area of approximately 9,138 square kilometers in the Everglades ecosystem. This model has the architecture of Convolutional Long Short-Term Memory which can deal with spatiotemporal data by capturing both spatial and temporal dependencies in the training data. The forecasting capabilities of this model (referred to as the global model) are assessed by comparing the global model to two Artificial Neural Networks developed at two different gaging stations, referred to here as local models. One local model is developed at a gaging station directly influenced by nearby water control structures, whereas the other is developed at a gaging station located farther away from these structures. By leveraging data from the Everglades Depth Estimation Network spanning from January 2002 to May 2023, the global and local models were trained to forecast water levels with a two-day lead time. Our findings suggest that both the global and local models perform with approximately the same level of accuracy, with Mean Absolute Relative Error values ranging from 0.38% to 1.4% at the selected stations. The developed global model has demonstrated strong potential as a standalone forecasting tool for the entire study area in the Everglades and could eliminate the need for developing multiple local models. This finding also highlights how machine learning can capture complex spatial and temporal relationships to generate accurate water level predictions on a regional scale.

Florida↗
USGS2025

Wetlands, groundwater and seasonality influence the spatial distribution of stream chemistry in a low-relief catchment

Evaluating stream water chemistry patterns provides insight into catchment ecosystem and hydrologic processes. Spatially distributed patterns and controls of stream solutes are well-established for high-relief catchments where solute flow paths align with surface topography. However, the controls on solute patterns are poorly constrained for low-relief catchments where hydrogeologic heterogeneities and river corridor features, like wetlands, may influence water and solute transport. Here, we provide a data set of solute patterns from 58 synoptic surveys across 28 sites and over 32 months in a low-relief wetland-rich catchment to determine the major surface and subsurface controls along with wetland influence across the catchment. In this low-relief catchment, the expected wetland storage, processing, and transport of solutes is only apparent in solute patterns of the smallest subcatchments. Meanwhile, downstream seasonal and wetland influence on observed chemistry can be masked by large groundwater contributions to the main stream channel. These findings highlight the importance of incorporating variable groundwater contributions into catchment-scale studies for low-relief catchments, and that understanding the overall influence of wetlands on stream chemistry requires sampling across various spatial and temporal scales. Therefore, in low-relief wetland-rich catchments, given the mosaic of above and below ground controls on stream solutes, modeling efforts may need to include both surface and subsurface hydrological data and processes.

Michigan↗
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↗
USGS2025

A review of abrupt permafrost thaw: Definitions, usage, and a proposed conceptual framework

Purpose of Review We review how ‘abrupt thaw’ has been used in published studies, compare these definitions to abrupt processes in other Earth science disciplines, and provide a definitive framework for how abrupt thaw should be used in the context of permafrost science. Recent Findings We address several aspects of permafrost systems necessary for abrupt thaw to occur and propose a framework for classifying permafrost processes as abrupt thaw in the future. Based on a literature review and our collective expertise, we propose that abrupt thaw refers to thaw processes that lead to a substantial persistent environmental change within a few decades. Abrupt thaw typically occurs in ice-rich permafrost but may be initiated in ice-poor permafrost by external factors such as hydrologic change (i.e., increased streamflow, soil moisture fluctuations, altered groundwater recharge) or wildfire. Summary Permafrost thaw alters greenhouse gas emissions, soil and vegetation properties, and hydrologic flow, threatening infrastructure and the cultures and livelihoods of northern communities. The term ‘abrupt thaw’ has emerged in scientific discourse over the past two decades to differentiate processes that rapidly impact large depths of permafrost, such as thermokarst, from more gradual, top-down thaw processes that impact centimeters of near-surface permafrost over years to decades. However, there has been no formal definition for abrupt thaw and its use in the scientific literature has varied considerably. Our standardized definition of abrupt thaw offers a path forward to better understand drivers and patterns of abrupt thaw and its consequences for global greenhouse gas budgets, impacts to infrastructure and land-use, and Arctic policy- and decision-making.

Current Climate Change Reports↗
USGS2025

A review of post-wildfire adaptations of surface-water-quality models: Synthesis, gaps, and opportunities

As wildfires increasingly affect water-supply watersheds, the demand for models to predict water-quality responses is increasing. This work reviews and synthesizes existing post-wildfire applications of water-quality models in the context of geographic and ecohydrological distribution, hydrologic and water-quality response process representation, model parameterization, model and input data scales, model calibration data availability, as well as calibration and performance evaluation approaches. Emphasis is placed on models that simulate water-quality output, rather than sediment and erosional response as the primary focus. Here, identified gaps and opportunities to advance the post-wildfire application of water-quality models include: 1. applying models in under-represented geographic and ecohydrologic regions, 2. simulating multiple streamflow generation mechanisms, including groundwater, with an emphasis on shifting dominant flow pathways as the landscape recovers following wildfire, 3. adding studies that include the simulation of metals, 4. incorporating more biogeochemical and in-stream processes to model applications, 5. applying finer spatial and temporal resolution of precipitation data input as well as finer spatial resolution hydrologic response units, 6. implementing fully distributed grid or element models or finer resolution response units to capture burn severity heterogeneity, 7. collecting enhanced water-quality data for model calibration and validation, 8. conducting model-intercomparison studies, and 9. developing model parameter value guidance in post-wildfire applications. These identified gaps and opportunities may assist users in deciding on key processes and approaches to consider in modeling post-wildfire water-quality conditions.

Science of the Total Environment↗
USGS2025

Linking permafrost to the abundance, biomass, and energy density of fish in Arctic headwater streams

Permafrost thaw alters groundwater flow, river hydrology, stream-catchment interactions, and the availability of carbon and nutrients in headwater streams. The impact of permafrost on watershed hydrology and biogeochemistry of headwater streams has been demonstrated, but there is little understanding of how permafrost influences fish in these ecosystems. We examined relations among permafrost characteristics, the resulting changes in water temperature, stream hydrology (e.g., discharge flashiness), and macroinvertebrates, with the abundance, biomass, and energy density of juvenile Dolly Varden ( Salvelinus malma ) and Arctic Grayling ( Thymallus arcticus ) across 10 headwater streams in northwestern Alaska. Macroinvertebrate density was driven by concentrations of dissolved carbon and nutrients supporting stream food webs. Dolly Varden abundance was primarily related to water temperature with fewer fish in warmer streams, whereas Dolly Varden energy density decreased with the flashiness of the headwater streams. Dolly Varden biomass was related to both temperature and bottom-up food web effects. The energy density of Arctic Grayling decreased with warmer temperatures and discharge flashiness. These relations demonstrate the importance of terrestrial–aquatic connections in permafrost landscapes and indicate the complexity of landscape effects on fish. Because permafrost thaw is one of the most impactful changes occurring as the Arctic warms, an improved understanding of how stream temperature, hydrology, and bottom-up food web processes influence fish populations can aid forecasting of future conditions across the Arctic.

Alaska↗
USGS2025

Per- and polyfluoroalkyl substances (PFAS) mass flux and mass balance at an aqueous film-forming foam release site in semiarid eastern New Mexico, USA

Passive flux meters (PFMs) directly measure groundwater chemistry mass flux and Darcy flux, providing insight into contaminant source-zone architecture and transport properties. This study uses PFMs to characterize PFAS flux in groundwater at a semiarid site with a thick (greater than 90-m) unsaturated zone where groundwater has been contaminated with per - and polyfluoroalkyl substances (PFAS) related to the use of aqueous film-forming foam (AFFF) for fire training and fire suppression. PFAS mass discharge (PFAS mass flux integrated over a control plane) in groundwater downgradient from several PFAS release areas is calculated using PFM results. In groundwater downgradient from fire-training areas, total PFAS mass discharge (summed across 14 compounds) was estimated to be between 6.0 and 31 g per day in 2020 and between 5.9 and 23 g per day in 2021. Site-specific documentation, generic information on AFFF properties, and literature values of PFAS concentration in AFFF are used to estimate site-specific PFAS-application rates at fire-training areas. These PFAS-application rates are compared to groundwater PFAS-discharge rates. Results suggest that transformation processes (exact pathways unknown) have led to increased discharge of measured PFAS in groundwater relative to initial AFFF formulations. The mass balance approach has broad applicability as a high-level approach that can provide insight into PFAS transport at AFFF sites.

New Mexico↗
USGS2025

Multi-lithofacies alluvial characterization via airborne electromagnetic-borehole fusion using ordinary interval kriging and geologic constraints

Airborne electromagnetic (AEM) survey provides extensive spatial coverage and detailed resolution at the near surface and can be used to develop hydrogeological models. However, utilization of AEM data is not straightforward because AEM resistivity is an indirect measurement for inferring sediment types. This study develops an ordinary interval kriging (OIK) algorithm and a resistivity-to-multi-lithofacies (R2ML) data fusion workflow for multi-lithofacies alluvial characterization. OIK utilizes irregular interval data to construct three-dimensional (3D) resistivity fields from one-dimensional inverted AEM resistivity models. The R2ML workflow maps the resistivity field generated from OIK into a multi-facies lithological model, incorporating geologic constraints derived from well logs and geological observations. The numerical and real-world cases demonstrate that OIK is computationally efficient, accounts for 3D anisotropy, and minimizes the smoothing effect, thereby preserving resistivity contrasts and reducing interpolation uncertainty. The methodology is applied to lithologic characterization of the Mississippi River Valley alluvial aquifer (MRVA) in the Shellmound area, Mississippi, U.S. A frequency-domain AEM survey was conducted to support groundwater studies for the managed aquifer recharge (MAR) to the MRVA. The resulting lithological model, including four types of lithofacies—clay, very fine sand, fine-medium sands, and graveliferous sands, illustrates the geomorphological processes of the MRVA and implies potential MAR. The alignment between the lithological model and existing geological and hydrogeological investigations demonstrates that OIK and R2ML workflow effectively capture the subsurface architecture of the MRVA. The methods have broad applicability for characterizing alluvial aquifers through AEM-borehole data fusion, supporting sustainable groundwater management.

Mississippi↗
USGS1998

Characterization of hydrogeologic units using matrix properties, Yucca Mountain, Nevada

Determination of the suitability of Yucca Mountain, in southern Nevada, as a geologic repository for high-level radioactive waste requires the use of numerical flow and transport models. Input for these models includes parameters that describe hydrologic properties and the initial and boundary conditions for all rock materials within the unsaturated zone, as well as some of the upper rocks in the saturated zone. There are 30 hydrogeologic units in the unsaturated zone, and each unit is defined by limited ranges where a discrete volume of rock contains similar hydrogeologic properties. These hydrogeologic units can be easily located in space by using three-dimensional lithostratigraphic models based on relation- ships of the properties with the lithostratigraphy. Physical properties of bulk density, porosity, and particle density; flow properties of saturated hydraulic conductivity and moisture-retention characteristics; and the state variables (variables describing the current state of field conditions) of saturation and water potential were determined for each unit. Units were defined using (1) a data base developed from 4,892 rock samples collected from the coring of 23 shallow and 8 deep boreholes, (2) described lithostratigraphic boundaries and corresponding relations to porosity, (3) recognition of transition zones with pronounced changes in properties over short vertical distances, (4) characterization of the influence of mineral alteration on hydrologic properties such as permeability and moisture-retention characteristics, and (5) a statistical analysis to evaluate where boundaries should be adjusted to minimize the variance within layers. This study describes the correlation of hydrologic properties to porosity, a property that is well related to the lithostratigraphy and depositional and cooling history of the volcanic deposits and can, therefore, be modeled to be distributed laterally. Parameters of the hydrogeologic units developed in this study and the relation of flow properties to porosity that are described can be used to produce detailed and accurate representations of the core-scale hydrologic processes ongoing at Yucca Mountain.

Nevada↗
USGS2026

Deep critical zone controls on shallow landslides

The deep critical zone (CZ) has long been recognized for its importance in influencing shallow landslides but was not considered feasible to include in slope stability models at the watershed scale. In this study, we demonstrate that simple approximations of the CZ in a fully coupled hydrologic and soil slope stability model can effectively capture the location, timing, and likely size of shallow landslides. To achieve this, we use coupled, process-based models that incorporate the effects of 1) deep CZ structures, 2) three-dimensional transient hydrology, and 3) multidimensional slope stability, calibrated with data from an intensively monitored field site. Our results show that the hydrologically active deep CZ guides groundwater flow, influencing where it drains from or exfiltrates to the soil mantle, producing distinct patterns of soil saturation and seepage forces at the soil-bedrock boundary. Deep conductive weathered critical zone drains the soil mantle, reducing the likelihood of destabilizing pore pressures, while the downslope thinning of the CZ forces groundwater to the surface. This creates localized instability and a tendency for similar-sized landslides across the landscape. In contrast, the absence of conductive weathered bedrock results in more widespread destabilizing pore pressures, leading to larger landslides and the likelihood of landslides earlier in a storm than in landscapes underlain by a deep CZ. Our findings suggest that first-order variations of deep CZ can provide physical explanations for variations observed in the susceptibility, magnitude, and timing of shallow landslides, and that CZ structure may be inferred from patterns and timing of landsliding.

Proceedings of the National Academy of Sciences↗
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↗