USGS Science⌕ Search

SEARCH · USGS Science

Results for “Hydrological Processes”

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.

At least 1,243 records · Page 69Linked to original sources
USGS2018

Hydrologic conditions and simulation of groundwater and surface water in the Great Dismal Swamp of Virginia and North Carolina

The U.S. Geological Survey (USGS), in cooperation with the U.S Fish and Wildlife Service, has investigated the hydrology of the Great Dismal Swamp (Swamp) National Wildlife Refuge (Refuge) in Virginia and North Carolina and developed a three-dimensional numerical model to simulate groundwater and surface-water hydrology. The model was developed with MODFLOW-NWT, a USGS numerical groundwater flow modeling program, in combination with the Surface-Water Routing Process, a software package that simulates dynamic surface-water flows, water control structure management, and groundwater/surface-water interactions. The steady-state model was calibrated to average spring conditions by using automated parameter estimation software (PEST) to reduce simulation errors and assess model parameter sensitivity. The model was then used to simulate wet and dry climatic conditions and a variety of hypothetical scenarios in which water levels in the Swamp were raised and lowered by simulated management of water control structures. Results of the model simulations indicate that, under average spring conditions, precipitation is the primary water input (92%); surface-water (5%) and groundwater (3%) inflows make up the remainder. The primary outflow (or loss) is evapotranspiration (55%), with surface outflows (about 41%) and groundwater outflow (about 4%) making up the remainder. Simulated adjustment of water control structure weir levels demonstrates that groundwater levels are affected by water levels in adjacent ditches and that surface-water and groundwater levels can be controlled through management of water control structures, allowing the Refuge to better manage fire risks and preserve forested-wetland ecosystems in the Refuge. The 13 water control structures proposed in the simulated scenario representing possible future conditions effectively raised simulated water levels in the northeastern corner of the study area, a goal of the Refuge management. Results of this study demonstrate use of MODFLOW with the Surface-Water Routing Process for simulating water management options in peat wetlands and will help Refuge managers to better understand existing hydrologic conditions, assess the hydrologic effects of planned changes to water control structures, and apply the new simulation tool to guide water management on the Refuge.

North Carolina↗
USGS2011

Spatial interpolation schemes of daily precipitation for hydrologic modeling

Distributed hydrologic models typically require spatial estimates of precipitation interpolated from sparsely located observational points to the specific grid points. We compare and contrast the performance of regression-based statistical methods for the spatial estimation of precipitation in two hydrologically different basins and confirmed that widely used regression-based estimation schemes fail to describe the realistic spatial variability of daily precipitation field. The methods assessed are: (1) inverse distance weighted average; (2) multiple linear regression (MLR); (3) climatological MLR; and (4) locally weighted polynomial regression (LWP). In order to improve the performance of the interpolations, the authors propose a two-step regression technique for effective daily precipitation estimation. In this simple two-step estimation process, precipitation occurrence is first generated via a logistic regression model before estimate the amount of precipitation separately on wet days. This process generated the precipitation occurrence, amount, and spatial correlation effectively. A distributed hydrologic model (PRMS) was used for the impact analysis in daily time step simulation. Multiple simulations suggested noticeable differences between the input alternatives generated by three different interpolation schemes. Differences are shown in overall simulation error against the observations, degree of explained variability, and seasonal volumes. Simulated streamflows also showed different characteristics in mean, maximum, minimum, and peak flows. Given the same parameter optimization technique, LWP input showed least streamflow error in Alapaha basin and CMLR input showed least error (still very close to LWP) in Animas basin. All of the two-step interpolation inputs resulted in lower streamflow error compared to the directly interpolated inputs.

Colorado, Georgia↗
USGS2022

Small anthropogenic landforms from past charcoal production control moisture dynamics and chemistry in northcentral Appalachian soils

Throughout the northeastern United States (U.S.) and Europe, relict charcoal hearths (RCHs) are regularly being discovered in proximity to furnaces once used for the extraction of metal from ore or quick-lime production; charcoal produced in hearths was used as a furnace fuel. Given previous research has shown that topographic and subsurface disturbance can be great when a hearth is constructed, we hypothesize that hearth construction alters surface hydrology and soil chemistry in environments in and near hearths. We used a landscape classification process to identify 6,758 hearths near furnaces at Greenwood and Pine Grove Furnace State Park, central and southcentral Pennsylvania, U.S. Two types of digital elevation model wetness indexes were used to quantify surface hydrology effects in and around hearths. Modeled wetness conditions were compared to field soil volumetric water content in RCHs near Greenwood Furnace State Park. Modeled wetness indexes indicate that RCH interiors are significantly wetter than RCH rim areas; RCHs are acting as a landscape moisture sink. Results also indicate that RCHs on slopes result in downslope drier conditions below RCHs. Field measured volumetric water content indicates that as distance from the center of the hearth increases, soil moisture significantly decreases. Geomorphic position was found to not be related to RCH wetness. Soil from RCHs, compared to nearby native soils, has significantly higher total C, a lower Mehlich 3 extractable acidity, higher Ca and P. No trend was evident with RCH soil chemistry and geomorphic position. The high frequency of RCH occurrence, in proximity to the furnace’s RCHs supported, suggests that RCHs today could locally be an important niche for understory flora and fauna. Further research could explore how RCHs might be affecting surrounding plant populations and how within RCH patterns, especially on hillslopes, might represent a distinctly different scale of physical and chemical variability.

Pennsylvania↗
USGS1991

Geochemistry of ground water in alluvial basins of Arizona and adjacent parts of Nevada, New Mexico, and California

Chemical and isotope analyses of ground water from 28 basins in the Basin and Range physiographic province of Arizona and parts of adjacent States were used to evaluate ground-water quality, determine processes that control ground-water chemistry, provide independent insight into the hydrologic flow system, and develop information transfer. The area is characterized by north- to northwest-trending mountains separated by alluvial basins that form a regional topography of alternating mountains and valleys. On the basis of ground-water divides or zones of minimal basin interconnection, the area was divided into 72 basins, each representing an individual aquifer system. These systems are joined in a dendritic pattern and collectively constitute the major water resource in the region. Geochemical models were developed to identify reactions and mass transfer responsible for the chemical evolution of the ground water. On the basis of mineralogy and chemistry of the two major rock associations of the area, a felsic model and a mafic model were developed to illustrate geologic, climatic, and physiographic effects on ground-water chemistry. Two distinct hydrochemical processes were identified: (1) reactions of meteoric water with minerals and gases in recharge areas and (2) reactions of ground water as it moves down the hydraulic gradient. Reactions occurring in recharge and downgradient areas can be described by a 13-component system. Major reactions are the dissolution and precipitation of calcite and dolomite, the weathering of feldspars and ferromagnesian minerals, the formation of montmorillonite, iron oxyhydroxides, and probably silica, and, in some basins, ion exchange. The geochemical modeling demonstrated that relatively few phases are required to derive the ground-water chemistry; 14 phases-12 mineral and 2 gas-consistently account for the chemical evolution in each basin. The final phases were selected through analysis of X-ray diffraction and fluorescence data, aqueous speciation and saturation data, and mass-balance and isotopic constraints and through chemical models developed from mineral combinations among the 27 phases that were considered realistic in these geologically and mineralogically complex basins. X-ray diffraction of basin-fill sediments confirm the presence of the postulated minerals and their weathering sequences. High partial pressures of soil CO2 and large concentrations of dissolved CO2 in recharge areas, and the rapid depletion of CO2 downgradient, accompanied by high weathering rates of the silicates which also decrease downgradient, indicate that carbonic acid is the impetus in the weathering process. Reactions in the soil zone and the unsaturated zone are influential and, in some instances, are as important as the mineralogy of the source rock in determining ground-water compositions. The basins can be divided geochemically into two general categories-closed systems, which evolve under closed hydrologic conditions, and open systems, which are open to CO2 and other constituents along the flow path. The ground-water chemistry of the unconfined aquifers in the eastern part of the study area and of the aquifers underlying the flood plain along the Colorado River generally evolves under open conditions. The ground-water chemistry of most basins in the central and western parts and of the confined aquifers in the eastern part evolves under closed conditions. The factors that determine whether a basin is an open or closed system are the amount of and the spatial and seasonal distribution of annual precipitation and the presence or absence of fine-grained confining units. The basins along the Colorado River are unique among basins in the region. Virtually all ground water underlying the flood plain originated as seepage or overbank flow from the Colorado River. Initial deuterium content of about -120 per mil is indicative of precipitation from the central part of Colorado. Using chemical m

Professional Paper↗
USGS2016

Pleistocene Lake Bonneville as an analog for extraterrestrial lakes and oceans

Geomorphic confirmation for a putative ancient Mars ocean relies on analog comparisons of coastal-like features such as shoreline feature attributes and temporal scales of process formation. Pleistocene Lake Bonneville is one of the few large, geologically young, terrestrial lake systems that exemplify well-preserved shoreline characteristics that formed quickly, on the order of a thousand years or less. Studies of Lake Bonneville provide two essential analog considerations for interpreting shorelines on Mars: (1) morphological variations in expression depend on constructional vs erosional processes, and (2) shorelines are not always correlative at an equipotential elevation across a basin due to isostasy, heat flow, wave setup, fetch, and other factors. Although other large terrestrial lake systems display supporting evidence for geomorphic comparisons, Lake Bonneville encompasses the most integrated examples of preserved coastal features related to basin history, sediment supply, climate, and fetch, all within the context of a detailed hydrograph. These collective terrestrial lessons provide a framework to evaluate possible boundary conditions for ancient Mars hydrology and large water body environmental feedbacks. This knowledge of shoreline characteristics, processes, and environments can support explorations of habitable environments and guide future mission explorations.

Book chapter↗
USGS1998

Hydrology of the Bonneville Salt Flats, northwestern Utah, and simulation of ground-water flow and solute transport in the shallow-brine aquifer

This report describes the hydrologic system of the Bonneville Salt Flats with emphasis on the mechanisms of solute transport. Variable-density, three-dimensional computer simulations of the near-surface part of the ground-water system were done to quantify both the transport of salt dissolved in subsurface brine that leaves the salt-crust area and the salt dissolved and precipitated on the land surface. The study was designed to define the hydrology of the brine ground-water system and the natural and anthropogenic processes causing salt loss, and where feasible, to quantify these processes. Specific areas of study include the transport of salt in solution by ground-water flow and the transport of salt in solution by wind-driven ponds and the subsequent salt precipitation on the surface of the playa upon evaporation or seepage into the subsurface. In addition, hydraulic and chemical changes in the hydrologic system since previous studies were documented.

Utah↗
USGS2006

Solute transport processes

Soils and aquifers are major compartments of the subsurface environment, which together control the terrestrial hydrological cycle. This subsurface is important for water resources and also as repository for municipal, industrial, and government waste. Aquifers are typically recharged by natural rainfall entering the soil profile and leaching into deeper soil layers. Due to intensive agricultural or industrial activities the leachate leaving the soil profile and entering the aquifer may contain concentrations of toxic substances such as agrochemicals, heavy metals, and polycyclic aromatic hydrocarbons. At contaminated industrial sites light and dense non-aqueous phase liquids (LNAPLs, DNAPLs) may be transported as a separate phase to the underlying aquifer systems. Once any of these chemicals have entered the aquifer they can be transported over large horizontal distances thus contaminating large parts of the aquifer and threateningwater supplywells. Remediation of highly contaminated aquifer systems is commonly a long-term and expensive proposition. As safe and effective use of the subsurface environment is a major challenge facing our society, there is a great need to improve our understanding of the shallow subsurface and the groundwater systems. This particularly includes the understanding of transport processes, which are responsible for the fate of contaminants.

Book chapter↗
USGS2013

Remote sensing of soil moisture using airborne hyperspectral data

Landscape assessment of soil moisture is critical to understanding the hydrological cycle at the regional scale and in broad-scale studies of biophysical processes affected by global climate changes in temperature and precipitation. Traditional efforts to measure soil moisture have been principally restricted to in situ measurements, so remote sensing techniques are often employed. Hyperspectral sensors with finer spatial resolution and narrow band widths may offer an alternative to traditional multispectral analysis of soil moisture, particularly in landscapes with high spatial heterogeneity. This preliminary research evaluates the ability of remotely sensed hyperspectral data to quantify soil moisture for the Little River Experimental Watershed (LREW), Georgia. An airborne hyperspectral instrument with a short-wavelength infrared (SWIR) sensor was flown in 2005 and 2007 and the results were correlated to in situ soil moisture values. A significant statistical correlation ( R 2 value above 0.7 for both sampling dates) for the hyperspectral instrument data and the soil moisture probe data at 5.08 cm (2 inches) was determined. While models for the 20.32 cm (8 inches) and 30.48 cm (12 inches) depths were tested, they were not able to estimate soil moisture to the same degree.

GIScience and Remote Sensing↗
USGS2014

Characterizing the primary material sources and dominant erosional processes for post-fire debris-flow initiation in a headwater basin using multi-temporal terrestrial laser scanning data

Wildfire dramatically alters the hydrologic response of a watershed such that even modest rainstorms can produce hazardous debris flows. Relative to shallow landslides, the primary sources of material and dominant erosional processes that contribute to post-fire debris-flow initiation are poorly constrained. Improving our understanding of how and where material is eroded from a watershed during a post-fire debris-flow requires (1) precise measurements of topographic change to calculate volumetric measurements of erosion and deposition, and (2) the identification of relevant morphometrically defined process domains to spatially constrain these measurements of erosion and deposition. In this study, we combine the morphometric analysis of a steep, small (0.01 km 2 ) headwater drainage basin with measurements of topographic change using high-resolution (2.5 cm) multi-temporal terrestrial laser scanning data made before and after a post-fire debris flow. The results of the morphometric analysis are used to define four process domains: hillslope-divergent, hillslope-convergent, transitional, and channelized incision. We determine that hillslope-divergent and hillslope-convergent process domains represent the primary sources of material over the period of analysis in the study basin. From these results we conclude that raindrop-impact induced erosion, ravel, surface wash, and rilling are the primary erosional processes contributing to post-fire debris-flow initiation in the small, steep headwater basin. Further work is needed to determine (1) how these results vary with increasing drainage basin size, (2) how these data might scale upward for use with coarser resolution measurements of topography, and (3) how these results change with evolving sediment supply conditions and vegetation recovery.

California↗
USGS2014

Collaborative modelling and integrated decision support system analysis of a developed terminal lake basin

A terminal lake basin in west-central Nevada, Walker Lake, has undergone drastic change over the past 90 yrs due to upstream water use for agriculture. Decreased inflows to the lake have resulted in 100 km 2 decrease in lake surface area and a total loss of fisheries due to salinization. The ecologic health of Walker Lake is of great concern as the lake is a stopover point on the Pacific route for migratory birds from within and outside the United States. Stakeholders, water institutions, and scientists have engaged in collaborative modeling and the development of a decision support system that is being used to develop and analyze management change options to restore the lake. Here we use an integrated management and hydrologic model that relies on state-of-the-art simulation capabilities to evaluate the benefits of using integrated hydrologic models as components of a decision support system. Nonlinear feedbacks among climate, surface-water and groundwater exchanges, and water use present challenges for simulating realistic outcomes associated with management change. Integrated management and hydrologic modeling provides a means of simulating benefits associated with management change in the Walker River basin where drastic changes in the hydrologic landscape have taken place over the last century. Through the collaborative modeling process, stakeholder support is increasing and possibly leading to management change options that result in reductions in Walker Lake salt concentrations, as simulated by the decision support system.

Nevada↗
USGS2025

Investigating hurricane-induced salt variation across the Land-Estuary-Ocean continuum using a dynamically coupled hydrological-ocean model

Salinity variations across the Land‐Estuary‐Ocean (LEO) continuum are critical for coastal ecosystems and impact the socioeconomic benefits to local communities. However, evaluating these variations is challenging due to the complex interactions of terrestrial and oceanic processes, including river discharge, winds, tides, sea level rise, and storms. This study incorporates a salinity module into a newly developed dynamically coupled hydrological‐ocean model on the COAWST platform, which is, to the best of our knowledge, the first process‐based modeling that achieves seamless two‐way land‐ocean water and salt exchanges across the LEO continuum, facilitating simulations of landward salt input and transport. Using Hurricane Florence (2018) as a case study, we examined how various physical processes influence salinity dynamics and freshwater plume development in the Cape Fear River Estuary, North Carolina. The model simulated saltwater intrusion into freshwater wetlands upstream of the estuary and demonstrated that salinity in the estuary was initially regulated by wind‐driven water level gradients, followed by a dominant influence from significant river runoff. In the coastal ocean, runoff created a large freshwater plume that moved westward, driven by the interplay between runoff, winds, and the estuary's geomorphology. The excellent performance of the coupling approach for salinity simulation underscores the importance of the seamless transport of water and salt at the land‐ocean interface. This study demonstrates that the coupled model is a valuable tool for representing saltwater intrusion, tracking coastal pollutants, and understanding water and material exchange across the LEO continuum.

North Carolina, South Carolina↗
USGS2021

Unsaturated flow processes and the onset of seasonal deformation in slow-moving landslides

Predicting rainfall-induced landslide motion is challenging because shallow groundwater flow is extremely sensitive to the preexisting moisture content in the ground. Here, we use groundwater hydrology theory and numerical modeling combined with five years of field monitoring to illustrate how unsaturated groundwater flow processes modulate the seasonal pore water pressure rise and therefore the onset of motion for slow-moving landslides. The onset of landslide motion at Oak Ridge earthflow in California’s Diablo Range occurs after an abrupt water table rise to near the landslide surface 52–129 days after seasonal rainfall commences. Model results and theory suggest that this abrupt rise occurs from the advection of a nearly saturated wetting front, which marks the leading edge of the integrated downward flux of seasonal rainfall, to the water table. Prior to this abrupt rise, we observe little measured pore water pressure response within the landslide due to rainfall. However, once the wetting front reaches the water table, we observe nearly instantaneous pore water pressure transmission within the landslide body that is accompanied by landslide acceleration. We cast the timescale to reach a critical pore water pressure threshold using a simple mass balance model that considers variable moisture storage with depth and explains the onset of seasonal landslide motion with a rainfall intensity-duration threshold. Our model shows that the seasonal response time of slow-moving landslides is controlled by the dry season vadose zone depth rather than the total landslide thickness.

California↗
USGS2018

Rio Grande transboundary integrated hydrologic model and water-availability analysis, New Mexico and Texas, United States, and Northern Chihuahua, Mexico

Errata **September 28, 2018: The purpose of a USGS Open-file report (OFR) is dissemination of information that must be released immediately to fill a public need or information that is not sufficiently refined to warrant publication in one of the other USGS series. As part of that refinement process, an error was discovered in one of the input data sets of the Rio Grande Transboundary Integrated Hydrologic Model (RGTIHM) that this OFR was based upon. The error involved the assignment of storage properties to “phantom cells.” Phantom cells are required for most variants of MODFLOW that use a structured finite-difference grid when individual stratigraphic layers are represented as separate layers. Using phantom cells is a common practice that allows separate model layers to be maintained without having to combine stratigraphic layers into equivalent model layers or to use an unstructured grid. Typically, phantom cell horizontal hydraulic conductivities and storage properties are set to a small number and vertical hydraulic conductivities are set to a number large enough to allow vertical flow between the vertically adjacent layers. In the RGTIHM, the specific storage properties of the phantom cells for the upper (RGTIHM layers 3 and 4), middle (RGTIHM layers 5 and 6), and lower (RGTIHM layers 7 and 8) members of the Santa Fe Group were inadvertently assigned a value of 1 feet -1 . The revision of these specific storage values to a small number (1.0 x 10 -09 feet -1 ) required additional trial-and-error model calibration and a new sensitivity analysis. After calibration, the overall model fit remained similar to the fit described in the OFR, but the fit for many individual features such as project water available for diversions at the American Canal and Acequia Madre improved due to the reduction in flow coming from lower layers. Overall, there is still an average net depletion of groundwater flow, and the conclusions of the report are not changed. The revised average annual groundwater flow depletion simulated for the period 1953-2014 is -1,480 acre-feet/year for the entire model region, and -3,660 acre-feet/year for the portion of the model in the United States. The final version of the model will be the basis of the USGS Scientific Investigations Report that will supersede this OFR. An updated Model Archive of RGTIHM is available upon request to the USGS California Water Science Center. The corrected version of the model WAS the basis for the USGS Scientific Investigations Report that SUPERSEDED this Open-File Report.** Abstract Changes in population, agricultural development and practices (including shifts to more water-intensive crops), and climate variability are increasing demands on available water resources, particularly groundwater, in one of the most productive agricultural regions in the Southwest—the Rincon and Mesilla Valley parts of Rio Grande Valley, Doña Ana and Sierra Counties, New Mexico, and El Paso County, Texas. The goal of this study was to produce an integrated hydrological simulation model to help evaluate water-management strategies, including conjunctive use of surface water and groundwater for historical conditions, and to support long-term planning for the Rio Grande Project. This report describes model construction and applications by the U.S. Geological Survey, working in cooperation and collaboration with the Bureau of Reclamation. This model, the Rio Grande Transboundary Integrated Hydrologic Model, simulates the most important natural and human components of the hydrologic system, including selected components related to variations in climate, thereby providing a reliable assessment of surface-water and groundwater conditions and processes that can inform water users and help improve planning for future conditions and sustained operations of the Rio Grande Project (RGP) by the Bureau of Reclamation. Model development included a revision of the conceptual model of the flow system, construction of a Transboundary Rio Grande Watershed Model (TRGWM) water-balance model using the Basin Characterization Model (BCM), and construction of an integrated hydrologic flow model with MODFLOW-One-Water Hydrologic Flow Model (referred to as One Water). The hydrologic models were developed for and calibrated to historical conditions of water and land use, and parameters were adjusted so that simulated values closely matched available measurements (calibration). The calibrated model was then used to assess the use and movement of water in the Rincon Valley, Mesilla Basin, and northern part of the Conejos-Médanos Basin, with the entire region referred to as the “Transboundary Rio Grande” or TRG. These tools provide a means to understand hydrologic system response to the evolution of water use in the region, its availability, and potential operational constraints of the RGP. The conceptual model identified surface-water and groundwater inflows and outflows that included the movement and use of water both in natural and in anthropogenic systems. The groundwater-flow system is characterized by a layered geologic sedimentary sequence combined with the effects of groundwater pumping, operation of the RGP, natural runoff and recharge, and the application of irrigation water at the land surface that is captured and reused in an extensive network of canals and drains as part of the conjunctive use of water in the region. Historical groundwater-level fluctuations followed a cyclic pattern that were aligned with climate cycles, which collectively resulted in alternating periods of wet or dry years. Periods of drought that persisted for one or more years are associated with low surface-water availability that resulted in higher rates of groundwater-level decline. Rates of groundwater-level decline also increased during periods of agricultural intensification, which necessitated increasing use of groundwater as a source of irrigation water. Agriculture in the area was initially dominated by alfalfa and cotton, but since 1970 more water-intensive pecan orchards and vegetable production have become more common. Groundwater levels substantially declined in subregions where drier climate combined with increased demand, resulting in periods of reduced streamflows. Most of the groundwater was recharged in the Rio Grande Valley floor, and most of the pumpage and aquifer storage depletion was in Mesilla Basin agricultural subregions. A cyclic imbalance between inflows and outflows resulted in the modeled cyclic depletion (groundwater withdrawals in excess of natural recharge) of the groundwater basin during the 75-year simulation period of 1940–2014. Changes in groundwater storage can vary considerably from year to year, depending on land use, pumpage, and climate conditions. Climatic drivers of wet and dry years can greatly affect all inflows, outflows, and water use. Although streamflow and, to a minor extent, precipitation during inter-decadal wet-year periods replenished the groundwater historically, contemporary water use and storage depletion could have reduced the effects of these major recharge events. The average net groundwater flow-rate deficit for 1953–2014 was estimated to be about 8,990 acre-feet per year.

New Mexico, Northern Chihuahua, Texas↗
USGS2020

DTSGUI: A python program to process and visualize fiber‐optic distributed temperature sensing data

Fiber‐optic distributed temperature sensing (FO‐DTS) has proven to be a transformative technology for the hydrologic sciences, with application to diverse problems including hyporheic exchange, groundwater/surface‐water interaction, fractured‐rock characterization, and cold regions hydrology. FO‐DTS produces large, complex, and information‐rich datasets. Despite the potential of FO‐DTS, adoption of the technology has been impeded by lack of tools for data processing, analysis, and visualization. New tools are needed to efficiently and fully capitalize on the information content of FO‐DTS datasets. To this end, we present DTSGUI, a public‐domain Python‐based software package for editing, parsing, processing, statistical analysis, georeferencing, and visualization of FO‐DTS data.

Groundwater↗
USGS2024

Biogeochemical and community ecology responses to the wetting of non-perennial streams

Transitions between dry and wet hydrologic states are the defining characteristic of non-perennial rivers and streams, which constitute the majority of the global river network. Although past work has focused on stream drying characteristics, there has been less focus on how hydrology, ecology and biogeochemistry respond and interact during stream wetting. Wetting mechanisms are highly variable and can range from dramatic floods and debris flows to gradual saturation by upwelling groundwater. This variation in wetting affects ecological and biogeochemical functions, including nutrient processing, sediment transport and the assembly of biotic communities. Here we synthesize evidence describing the hydrological mechanisms underpinning different types of wetting regimes, the associated biogeochemical and organismal responses, and the potential scientific and management implications for downstream ecosystems. This combined multidisciplinary understanding of wetting dynamics in non-perennial streams will be key to predicting and managing for the effects of climate change on non-perennial ecosystems.

Nature Water↗
USGS2006

User's guide to the Variably Saturated Flow (VSF) process to MODFLOW

A new process for simulating three-dimensional (3-D) variably saturated flow (VSF) using Richards' equation has been added to the 3-D modular finite-difference ground-water model MODFLOW. Five new packages are presented here as part of the VSF Process--the Richards' Equation Flow (REF1) Package, the Seepage Face (SPF1) Package, the Surface Ponding (PND1) Package, the Surface Evaporation (SEV1) Package, and the Root Zone Evapotranspiration (RZE1) Package. Additionally, a new Adaptive Time-Stepping (ATS1) Package is presented for use by both the Ground-Water Flow (GWF) Process and VSF. The VSF Process allows simulation of flow in unsaturated media above the ground-water zone and facilitates modeling of ground-water/surface-water interactions. Model performance is evaluated by comparison to an analytical solution for one-dimensional (1-D) constant-head infiltration (Dirichlet boundary condition), field experimental data for a 1-D constant-head infiltration, laboratory experimental data for two-dimensional (2-D) constant-flux infiltration (Neumann boundary condition), laboratory experimental data for 2-D transient drainage through a seepage face, and numerical model results (VS2DT) of a 2-D flow-path simulation using realistic surface boundary conditions. A hypothetical 3-D example case also is presented to demonstrate the new capability using periodic boundary conditions (for example, daily precipitation) and varied surface topography over a larger spatial scale (0.133 square kilometer). The new model capabilities retain the modular structure of the MODFLOW code and preserve MODFLOW's existing capabilities as well as compatibility with commercial pre-/post-processors. The overall success of the VSF Process in simulating mixed boundary conditions and variable soil types demonstrates its utility for future hydrologic investigations. This report presents a new flow package implementing the governing equations for variably saturated ground-water flow, four new boundary condition packages unique to unsaturated flow, the Adaptive Time-Stepping Package for use with both the GWF Process and the new VSF Process, detailed descriptions of the input and output files for each package, and six simulation examples verifying model performance.

Techniques and Methods↗
USGS2018

A physically based daily simulation of the glacier-dominated hydrology of the Copper River basin, Alaska

The large, highly glacierized Copper River basin is an important water resource for the south-central region of Alaska. Thus, information is needed on the reaction of its hydrologic timing and streamflow volumes to historical changes in climate, in order to assess the possible impact of future changes. However, the basin is remote, and therefore, it has proved difficult to collect field data in a frequent temporal and spatial manner. An extension of the distributed-parameter, physical-process code Precipitation Runoff Modeling System, PRMSglacier, has been specifically developed to simulate daily hydrology without requiring extensive input data. In this study, PRMSglacier was used to characterize the hydrology of the Copper River basin from 1959 to 2015. The basin was split into subbasins for specific regional climatic calibrations and finer resolution characterization. The model was calibrated and performed well against data of glacier mass balance, glacier area change, snow cover, gaged streamflow, evapotranspiration, and solar radiation. Ice melt contributed 26% of the total basin streamflow, with differences temporally from climate oscillations. Furthermore, differences were seen geographically in subbasins depending on the state of the glaciers in each subbasin. Decreasing trends in ice melt volume were mostly seen on smaller steeper glaciers responding to a critical level of glacier recession, while increasing trends in ice melt volume were mostly seen on larger valley glaciers responding to increasing temperature. The areas with substantially decreasing ice melt had decreasing streamflow, possibly indicating health concerns for the ecosystems therein.

Alaska↗
USGS2025

Disentangling geomorphic equifinality in sediment and hydrologic connectivity through the analyses of landscape drivers of hysteresis

Sources, transport mechanisms and pathways of fine sediment in river systems are dependent on a multitude of climatic, geomorphic and anthropogenic factors, resulting in geomorphic equifinality, in which it is difficult to parse how different landscape processes affect sediment transport across different spatiotemporal scales. The objectives of this study are to 1) provide a conceptual model to consider how differing spatial distributions and hydrologic timing of sediment sources, both upland and in-channel, can result in equifinal sediment transport outcomes, and 2) utilize analytical methods with widely available environmental datasets to infer sediment processes from stream gaging data. Hysteretic patterns of observed storm events were classified based on their direction and timing of peak sediment concentration, relative to streamflow, using records from 35 U.S. Geological Survey stream gages in the period between 2007 and 2023 within two different physiographic regions: the Mid-Atlantic Delaware River Basin (DRB) and the Midwestern Illinois River Basin (IRB). The DRB contains mixed forest, urban, suburban and agricultural watersheds over diverse topography, and the IRB is primarily an intensively managed agricultural watershed on flat terrain. We use principal component analysis and linear discriminant analysis to infer regional hydrologic relations with turbidity dynamics, and to identify the primary hydrologic and land surface characteristics most effective at distinguishing between hysteretic classes in each region. These analyses reveal underlying regional relations in storm event hydrodynamics and landscape characteristics that contribute to varying patterns in sediment dynamics. Incorporating these sediment dynamic relations with spatial distributions and hydrologic timing of sediment sources could help to improve process understanding and predictive capability of fine sediment transport in watersheds.

Delaware, Illinois, Indiana, New Jersey, New York,↗