Search USGSSearch

SEARCH · Search USGS

Results for “Water Research”

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 883 records · Page 49Linked to original sources

Integrated hydrology and operations modeling to evaluate climate change impacts in an agricultural valley irrigated with snowmelt runoff

Applying models to developed agricultural regions remains a difficult problem because there are no existing modeling codes that represent both the complex physics of the hydrology and anthropogenic manipulations to water distribution and consumption. We apply an integrated groundwater – surface water and hydrologic river operations model to an irrigated river valley in northwestern Nevada/northern California, United States to evaluate the impacts of climate change on snow-fed agricultural systems that use surface water and groundwater conjunctively. We explicitly represent individual surface water rights within the hydrologic model and allow the integrated code to change river diversions in response to earlier snowmelt runoff and water availability. Historically under-used supplemental groundwater rights are dynamically activated within the model to offset diminished surface water deliveries. The model accounts for feedbacks between the natural hydrology and anthropogenic stresses, which is a first-of-its-kind assessment of the impacts of climate change on individual water rights, and more broadly on river basin operations. Earlier snowmelt decreases annual surface water deliveries to all water rights, not just the junior water rights, owing to a lack of surface water storage in the upper river basin capable of capturing earlier runoff. Conversely, downstream irrigators with access to reservoir storage benefit from earlier runoff flowing past upstream points of diversion prior to the start of the irrigation season. Despite regional shifts toward greater reliance on groundwater for irrigation, crop consumption (a common surrogate for crop yield) decreases due to spatiotemporal changes in water supply that preferentially impact a subset of growers in the region.

California

Unifying advective and diffusive descriptions of bedform pumping in the benthic biolayer of streams

Many water quality and ecosystem functions performed by streams occur in the benthic biolayer, the biologically active upper (~5 cm) layer of the streambed. Solute transport through the benthic biolayer is facilitated by bedform pumping, a physical process in which dynamic and static pressure variations over the surface of stationary bedforms (e.g., ripples and dunes) drive flow across the sediment-water interface. In this paper we derive two predictive modeling frameworks, one advective and the other diffusive, for solute transport through the benthic biolayer by bedform pumping. Both frameworks closely reproduce patterns and rates of bedform pumping previously measured in the laboratory, provided that the diffusion model's dispersion coefficient declines exponentially with depth. They are also functionally equivalent, such that parameter sets inferred from the 2D advective model can be applied to the 1D diffusive model, and vice versa. The functional equivalence and complementary strengths of these two models expand the range of questions that can be answered, for example, by adopting the 2D advective model to study the effects of geomorphic processes (such as bedform adjustments to land use change) on flow-dependent processes and the 1D diffusive model to study problems where multiple transport mechanisms combine (such as bedform pumping and turbulent diffusion). By unifying 2D advective and 1D diffusive descriptions of bedform pumping, our analytical results provide a straightforward and computationally efficient approach for predicting, and better understanding, solute transport in the benthic biolayer of streams and coastal sediments.

Water Resources Research

Characterizing physical properties of streambed interface sediments using in situ complex electrical conductivity measurements

Streambed sediment physical properties such as surface area, are difficult to quantify in situ but exert a high‐level control on a wide range of biogeochemical processes and sorption of contaminants. We introduce the use of complex electrical conductivity (CC) methods (also known as spectral induced polarization (SIP)) that measure both real and imaginary conductivity to non‐invasively and efficiently characterize shallow streambed sediments. We explore the method through synthetic modeling, laboratory, and field measurements to demonstrate the sensitivity of imaginary conductivity to sediment surface area, controlled in part by fine‐grained iron oxides produced by anoxic groundwater discharge. Laboratory measurements verify expected relationships between CC parameters and sediment properties. Synthetic modeling using a 1D analytical model illustrates the influence of water layer depth and conductivity on the field CC measurements made at the streambed‐stream water interface. Specifically, the inverted sediment imaginary conductivity is less impacted by uncertainty in the water layer depth and conductivity relative to the real conductivity and phase shift. Field CC measurements along a landfill‐impacted river reveal discrete streambed zones with enhanced bulk surface area generally corresponding to anoxic groundwater discharges zones with high concentrations of fine‐grained iron oxide precipitates.

Water Resources Research

Improved prediction of management-relevant groundwater discharge characteristics throughout river networks

Groundwater discharge zones connect aquifers to surface water, generating baseflow and serving as ecosystem control points across aquatic ecosystems. The influence of groundwater discharge on surface flow connectivity, fate and transport of contaminants and nutrients, and thermal habitat depends strongly on hydrologic characteristics such as the spatial distribution, age, and depth of source groundwater flow paths. Groundwater models have the potential to predict spatial discharge characteristics within river networks, but models are often not evaluated against these critical characteristics and model equifinality with respect to discharge processes is a known challenge. We quantify discharge characteristics across a suite of groundwater models with commonly used frameworks and calibration data. We developed a base model (MODFLOW‐NWT) for a 1,570‐km 2 watershed in the northeastern United States and varied the calibration data, control of river‐aquifer exchange directionality, and resolution. Most models ( n = 11 of 12) fit similarly to calibration metrics, but patterns in discharge location, flow path depth, and subsurface travel time varied substantially. We found (1) a 15% difference in the percent of discharge going to first‐order streams, (2) threefold variations in flow path depth, and (3) sevenfold variations in the subsurface travel times among the models. We recalibrated three models using a synthetic discharge location data set. Calibration with discharge location data reduced differences in simulated discharge characteristics, suggesting an approach to improved equifinality based on widespread field‐based mapping of discharge zones. Our work quantifying variation across common modeling approaches is an important step toward characterizing and improving predictions of groundwater discharge characteristics.

Connecticut, Massachusetts

Precipitation characteristics and land cover control wet season runoff source and rainfall partitioning in three humid tropical catchments in central Panama

Mechanisms of runoff generation in the humid tropics are poorly understood, particularly in the context of land-use/land cover change. This study analyzed the results of 124 storm hydrographs from three humid tropical catchments of markedly different vegetation cover and land-use history in central Panama during the 2017 wet season: actively grazed pasture, young secondary succession, and near-mature forest. We used electrical conductivity to separate baseflow (old water) from storm-event water (new-water). In all three land covers, new-water dominated storm runoff generation in 44% of the sampled storm events, indicating the dominance of fast shallow flow paths in the landscape. Activation of these flow paths was found to depend on a combination of maximum rainfall intensity and total storm rainfall, which, in turn, relates to markedly contrasting hydrograph separation results among land covers. Relationships between these rainfall characteristics and storm runoff generation were nonlinear, producing a threshold response with the exceedance of specific rainfall volumes and/or intensities. The pastoral catchment delivered order of magnitude more new-water during storm events than the two forested catchments. Although new-water contributed minimally (<10%) to total wet season runoff in the forested catchments, 43% of runoff generation in the pasture came from five large rainfall events where a threshold response produced substantial increases in total runoff and new-runoff efficiency. Based on our results, we propose a conceptual model of hydrologic flow paths in humid tropical systems that can explain previously observed disparities in seasonal storage and runoff with respect to land use/land cover.

Agua Salud Watershed

Assessment of peak flow scaling and Its effect on flood quantile estimation in the United Kingdom

Regional flood frequency analysis (RFFA) methods are essential tools to assess flood hazard and plan interventions for its mitigation. They are used to estimate flood quantiles when the at‐site record of streamflow data is not available or limited. One commonly used RFFA method is the index flood method (IFM), which assumes that peak floods satisfy the simple scaling hypothesis. In this work we present an integrated approach to assess the spatial scaling behavior of floods in the United Kingdom (UK) for 540 catchments, where the IFM is currently used operationally. This assessment employs product moments, probability weighted moments, and quantile analysis, and is applied to two different types of “hydrologically homogeneous” UK regions: geographical regions as defined in the Flood Studies Report (NERC, 1975) and pooling‐groups as defined in the updated Flood Estimation Handbook (FEH; Institute of Hydrology, 1999). To understand which variables play a significant role in the flood‐peak generating mechanism, the assessment approach considers scaling not only of drainage area alone but also of other hydro‐geomorphological variables. Results provided by the different methodologies consistently showed that only part (ranging from 30% to 70%) of the peak flow variability is explained by drainage area alone; this fraction increases (up to 80%–95%) when multiple regression is used. Supported by the peak flow spatial scaling assessment, we compared the proposed approach for peak flow quantile estimation with the current FEH method in ungauged catchments. The quantile regression method based on the pooling‐group outperforms the current FEH‐ungauged method, providing a 14% relative improvement in root mean square error over the entire country.

Water Resources Research

Using turbulence to identify preferential areas for grass carp (Ctenopharyngodon idella) larvae in streams: A laboratory study

In this experimental series, we studied the swimming capabilities and response of grass carp ( Ctenopharyngodon idella ) larvae to flow turbulence in a laboratory flume. We compared three different experimental configurations, representing in‐stream obstructions commonly found in natural streams (e.g., a gravel bump, a single vertical cylinder, and patches of submerged rigid vegetation). Grass carp larvae (postgas bladder emergence) were introduced to each experimental configuration and subjected to a variety of hydrodynamic forces of different magnitudes and scales. We varied the flow velocities and water depths and found ranges of turbulent kinetic energy and Reynolds stresses that triggered a response in larval trajectories, identified by measured horizontal and vertical swimming speeds for each flow condition. Larvae apparently actively avoided areas with increased levels of turbulence by swimming away, moving faster in short bursts, and expending more energy. In addition to the magnitude of turbulent kinetic energy, the length scale and time scale of turbulent eddies also influenced the larvae response. These findings support the development of new strategies for controlling the spread of grass carp larvae in rivers, as well as the development of numerical tools incorporating active swimming capabilities to predict larval transport in streams.

Water Resources Research

Baseflow age distributions and depth of active groundwater flow in a snow‐dominated mountain headwater basin

Deeper flows through bedrock in mountain watersheds could be important, but lack of data to characterize bedrock properties limits understanding. To address data scarcity, we combine a previously published integrated hydrologic model of a snow‐dominated, headwater basin of the Colorado River with a new method for dating baseflow age using dissolved gas tracers SF 6 , CFC‐113, N 2 , and Ar. The original flow model predicts the majority of groundwater flow through shallow alluvium (<8 m) sitting on top of less permeable bedrock. The water moves too quickly and is unable to reproduce observed SF 6 concentrations. To match gas data, bedrock permeability is increased to allow a larger fraction of deeper and older groundwater flow (median 112 m). The updated hydrologic model indicates interannual variability in baseflow age (3–12 years) is controlled by the volume of seasonal interflow and tightly coupled to snow accumulation and monsoon rain. Deeper groundwater flow remains stable (11.7 ± 0.7 years) as a function mean historical recharge to bedrock hydraulic conductivity (R/K). A sensitivity analysis suggests that increasing bedrock K effectively moves this alpine basin away from its original conceptualization of hyperlocalized groundwater flow (high R/K) with groundwater age insensitive to changes in water inputs. Instead, this basin is situated close to the precipitation threshold defining recharge controlled groundwater flow conditions (low R/K) in which groundwater age increases with small reductions in precipitation. Work stresses the need to explore alternative methods characterizing bedrock properties in mountain basins to better quantify deeper groundwater flow and predict their hydrologic response to change.

Colorado

Measuring and interpreting multilayer aquifer-system compactions for a sustainable groundwater-system development

Ever decreasing water resources and climate change have driven the increasing use of groundwater causing land subsidence in many countries. Geodetic sensors such as InSAR, GPS and leveling can detect surface deformation but cannot measure subsurface deformation. A single‐well, single‐depth extensometer can be used to measure subsurface deformation, but it cannot delineate the depths of major compaction and provide insight about the deformation mechanism throughout a complex aquifer system, unless man extensometers at different depths are used. We present a multilayer compaction well (MLCW), installed in a borehole, that uses magnetic rings to detect stratum compaction at 25 depths as deep as 300 m below land surface. Our laboratory and field assessments indicate 1 mm precision and accuracy for one single‐depth magnetic reading. We tested the performance of MLCW by measuring aquifer‐system compaction over the proximal, middle, and distal fans of the Choushui River Alluvial Fan (CRAF) that has long experienced severe land subsidence. The MLCW measurements were used to create time‐depth diagrams of compaction, showing different compaction rates at different layers of aquifers and aquitards to identify the depths of major compactions. The elastic (reversible) and inelastic (irreversible) compactions from MLCW were used in stress‐strain analyses to estimate skeletal specific storages and the safe groundwater levels, below which groundwater extractions have caused irreversible compactions. The hydrogeological parameters derived from MLCW measurements can help governmental agencies to determine effective land‐use and water‐use policies, and ascertain the best strategy for utilizing artificial recharge to prevent land subsidence and achieve sustainable groundwater management.

Water Resources Research

Machine learning predicted redox conditions in the glacial aquifer system, northern continental United States

Groundwater supplies 50% of drinking water worldwide and 30% in the United States. Geogenic and anthropogenic contaminants can, however, compromise water quality, thus limiting groundwater availability. Reduction/oxidation (redox) processes and redox conditions affect groundwater quality by influencing the mobility and transport of common geogenic and anthropogenic contaminants. In the glacial aquifer system, northern United States (GLAC, 1.87 million km 2 ), groundwater with high arsenic or manganese concentration is associated with reducing conditions and high nitrate with oxidizing conditions. This study uses machine learning to identify the relative influence of drivers of redox conditions (e.g., residence time vs. reactivity) across the glacial landscape. We developed three‐dimensional boosted regression tree models to predict redox conditions using the likelihood of low dissolved oxygen or high iron as indicators of anoxic conditions. Results indicate that variation in redox condition is controlled primarily by residence time (e.g., groundwater age and relative depth) and to a lesser extent by geochemical reactivity (e.g., subsurface contact time, soil carbon). Older water and deeper wells, along with more water storage or slower water movement was associated with higher probability of anoxic conditions. Mapped model results illustrate regions where anoxic redox conditions may mobilize geogenic contaminants or oxic conditions may limit denitrification potential. Results may also provide simplified redox input for process or predictive models of, for example, arsenic, manganese, or nitrate. Machine learning modeling methods can lead to improved understanding of contaminant occurrence and what drives redox conditions, and the methods may be transferable to other settings.

Northern Continental United States

Identification of seasonal streamflow regimes and streamflow drivers for daily and peak flows in Alaska

Alaska is among northern high‐latitude regions where accelerated climate change is expected to impact streamflow properties, including seasonality and primary flow drivers. Evaluating changes to streamflow, including flood characteristics, across this large and diverse environment can be improved by identifying the distribution and influence of flow drivers. Using metrics of mean monthly streamflow data from 253 streamgages, seasonal flow regimes were clustered to guide identification of seasonal‐flow drivers and form hydrologic groups for identification of peak‐flow populations. Nine seasonally distinct subclasses described variability within three classes dominated by (mostly fall) rainfall, (spring) snowmelt, and (summer) high‐elevation melt. The most glacierized basins exclusively grouped into high‐elevation melt subclasses, and less glacierized basins sometimes exhibited seasonal patterns aligned with rainfall‐ and snowmelt‐dominated regimes. Peak‐flow populations varied by subclass from dominant rainfall or dominant snowmelt to mixed rainfall‐snowmelt or mixed rainfall, snowmelt, and high‐elevation melt. Within subclasses, rainfall generated higher mean peak flows (relative to mean annual flow) than snowmelt or high‐elevation melt. Seasonal flow regimes showed clear but complex associations with basin characteristics, primarily elevation and winter temperature, and with geographic location. These dependencies provided elevation‐based analogies for changes associated with warming and insights for seasonal flow regime prediction and hydrologic region delineation. These results provide a spatially comprehensive perspective on seasonal streamflow drivers across Alaska from historical data and serve as an important historical basis for analysis.

Alaska

Characterization of groundwater recharge and flow in California's San Joaquin Valley from InSAR-observed surface deformation

Surface deformation in California's Central Valley (CV) has long been linked to changes in groundwater storage. Recent advances in remote sensing have enabled the mapping of CV deformation and associated changes in groundwater resources at increasingly higher spatiotemporal resolution. Here, we use interferometric synthetic aperture radar (InSAR) from the Sentinel‐1 missions, augmented by continuous Global Positioning System (cGPS) positioning, to characterize the surface deformation of the San Joaquin Valley (SJV, southern two‐thirds of the CV) for consecutive dry (2016) and wet (2017) water years. We separate trends and seasonal oscillations in deformation time series and interpret them in the context of surface and groundwater hydrology. We find that subsidence rates in 2016 (mean −42.0 mm/yr; peak −345 mm/yr) are twice that in 2017 (mean −20.4 mm/yr; peak −177 mm/yr), consistent with increased groundwater pumping in 2016 to offset the loss of surface‐water deliveries. Locations of greatest subsidence migrated outwards from the valley axis in the wetter 2017 water year, possibly reflecting a surplus of surface‐water supplies in the lowest portions of the SJV. Patterns in the amplitude of seasonal deformation and the timing of peak seasonal uplift reveal entry points and potential pathways for groundwater recharge into the SJV and subsequent groundwater flow within the aquifer. This study provides novel insight into the SJV aquifer system that can be used to constrain groundwater flow and subsidence models, which has relevance to groundwater management in the context of California's 2014 Sustainable Groundwater Management Act (SGMA).

California

Exploring the factors controlling the error characteristics of the Surface Water and Ocean Topography mission discharge estimates

The Surface Water and Ocean Topography (SWOT) satellite mission will measure river width, water surface elevation, and slope for rivers wider than 50-100 m. SWOT observations will enable estimation of river discharge by using simple flow laws such as the Manning-Strickler equation, complementing in-situ streamgages. Several discharge inversion algorithms designed to compute unobserved flow law parameters ( e.g. friction coefficient, bathymetry) have been proposed, but to date, a systematic assessment of factors controlling algorithm performance has not been conducted. Here, we assess the performance of the five algorithms that are expected to be used in the construction of the SWOT product. To perform this assessment, we used synthetic SWOT observations created with hydraulic model output corrupted with SWOT-like error. Prior information provided to the algorithms was purposefully limited to an estimate of mean annual flow (MAF), designed to produce a “worst case” benchmark. Prior MAF error was an important control on algorithm performance, but discharge estimates produced by the algorithms are less biased than the MAF; thus, the discharge algorithms improve on the prior. We show for the first time that accuracy and frequency of remote sensing observations are less important than prior bias, hydraulic variability among reaches, and flow law accuracy in governing discharge algorithm performance. The discharge errors and error sensitivities reported herein are a bounding benchmark, representing worst possible expected errors and error sensitivities. This study lays the groundwork to develop predictive power of algorithm performance, and thus map the global distribution of worst-case SWOT discharge accuracy.

Water Resources Research

Direct observation of the depth of active groundwater circulation in an alpine watershed

The depth of active groundwater circulation is a fundamental control on stream flows and chemistry in mountain watersheds, yet it remains challenging to characterize and is rarely well constrained. We collected hydraulic conductivity, hydraulic head, temperature, chemical, noble gas, and 3 H/ 3 He groundwater age data from discrete levels in two boreholes 46 and 81 m deep in an alpine watershed, in combination with chemical and age data from shallow groundwater discharge, to discern groundwater flow rates at different depths and directly observe active and inactive groundwater. Vertical head gradients are steep (average of 0.4) and thermal profiles are consistent with typical linear conductive continental geotherms. Groundwater deeper than ∼20 m is distinct from shallow groundwater and creek water in its chemistry, noble gas signature, and age (dominantly >65 years compared to <9 years). Together these results suggest low vertical groundwater flow velocities and a relatively shallow active circulation depth of ∼20 m. This hypothesis is tested with a simple 2‐D numerical fluid flow and heat transport model representing a hillslope transect through the two boreholes. The modeling indicates that the subhorizontally bedded sedimentary rocks underlying the basin are highly anisotropic with low vertical hydraulic conductivity, and at most ∼10% of bedrock recharge (equivalent to <2% of stream baseflow) flows below a depth of 20 m. The study demonstrates the considerable value of discrete‐depth hydrogeologic, chemical, and age data for determining active circulation depth, and illustrates an approach for maximizing the utility of individual boreholes drilled for mountain bedrock aquifer characterization.

Colorado

Substantial declines in salinity observed across the Upper Colorado River Basin during the 20th century, 1929 to 2019

Salinity in the Colorado River Basin causes an estimated $300 to $400 million per year in economic damages in the U.S. To inform and improve salinity‐control efforts, this study quantifies long‐term trends in salinity (dissolved solids) across the Upper Colorado River Basin (UCRB), including time periods prior to the construction of large dams and preceding the implementation of salinity‐control projects. Weighted Regressions on Time, Discharge, and Season was used with datasets of dissolved‐solids and specific‐conductance measurements, collected as early as 1929, to evaluate long‐term trends in dissolved‐solids loads and concentrations in streams from 1929 to 2019 ( n =14). Results indicate that large, widespread, and sustained downward trends in dissolved‐solids concentrations and loads occurred over the last 50 to 90 years. For 12 of the 14 stream sites with significant downward change, median declines of ‐38% (range of ‐14 to ‐57%) and ‐40% (range of ‐9 to ‐65%) were observed for flow‐normalized concentration and load, respectively. Steepest rates of decline occurred from 1980 to 2000, coincident with the initiation of salinity‐control efforts in the 1980s. However, there was a consistent slowing or reversing of downward trends after 2000 even though salinity‐control efforts continued. Significant decreases in salinity occurred as early as the 1940s at some streams, indicating that, in addition to salinity‐control projects, changes in land cover, land use, and/or climate substantially affect salinity transport in the UCRB. Observed dissolved‐solids trends are likely the result of changes to watershed‐related processes, not due to changes in the streamflow regime.

Arizona, Colorado, New Mexico, Utah, Wyoming

Influence of invasive submerged aquatic vegetation (E. densa) on currents and sediment transport in a freshwater tidal system

We present a field study combining measurements of vegetation density, vegetative drag, and reduction of suspended-sediment concentration (SSC) within patches of the invasive submerged aquatic plant Egeria densa . Our study was motivated by concern that sediment trapping by E. densa , which has proliferated in the Sacramento–San Joaquin Delta, is impacting marsh accretion and reducing turbidity. In the freshwater tidal Delta, E. densa occupies shallow regions, frequently along channel margins. We investigated two sites: Lindsey Slough, a muddy low-energy backwater, and the lower Mokelumne River, with stronger currents and sandy bed sediments. At the two sites biomass density, frontal area, and areal density of the submerged aquatic vegetation (SAV) were similar. Current attenuation within E. densa exceeded 90% and the vegetative drag coefficient followed , where is stem Reynolds number. The SAV reduced SSC by an average of 18% in Lindsey Slough. At Mokelumne River the reduction ranged 0–40%, with greatest trapping when discharge and SSC were elevated. This depletion of SSC decreases the transport of sediment to marshes by the same percentage, as the rising tide must pass through fringing SAV before reaching marshes. Sediment trapping in E. densa in the Delta is limited by low flux through the canopy and low settling velocity of suspended sediment (mostly flocculated mud). Sediment trapping by SAV has the potential to reduce channel SSC, but the magnitude and sign of the effect can vary with local factors including vegetative coverage and the depositional or erosional nature of the setting.

California

Estimating and forecasting time-varying groundwater recharge in fractured rock: A state-space formulation with preferential and diffuse flow to the water table

Rapid infiltration following precipitation may result in groundwater contamination from surface contaminants or pathogens. In fractured rock, contaminants can migrate rapidly to points of groundwater withdrawals. In contrast to the temporal availability of groundwater quality chemical indicators, meteorological and groundwater level observations are available in real-time to estimate time-varying recharge, which can act as a surrogate to identify periods of rapid infiltration that may indicate contamination susceptibility. Estimating recharge using methods, such as base-flow recession, unsaturated infiltration models, or Water-Table Fluctuations (WTF), cannot capitalize on currently available technologies and telecommunication infrastructure to conduct real-time recharge estimation at scales relevant to characterizing rapid infiltration. We present a linear, physics-based State-Space (SS) model of one-dimensional infiltration to estimate recharge, which includes preferential and diffuse-flow to the water table. The model can take advantage of real-time data for water-table altitude, precipitation, and evapotranspiration. Model parameters are calibrated over an observation period, and the Kalman Filter (KF) is subsequently applied to continuously update the observed (water-table altitude) and unobserved (groundwater recharge) system states and predict future states as new data become available. The SS/KF algorithm is demonstrated at the Masser Groundwater Recharge Site in Pennsylvania, USA and comparisons are made with recharge estimates from WTF methods. Model results indicate that the frequency of observations (daily versus sub-daily) dictates the allocation between preferential and diffuse flow. Additionally, because infiltration processes encompass many nonlinearities, model parameters estimated from observation periods need to be updated at least seasonally to account for changing recharge conditions.

Water Resources Research

Climate impacts on source contributions and evaporation to flow in the Snake River Basin using surface water isoscapes (δ2H and δ18O)

Rising global temperatures are expected to decrease the precipitation amount that falls as snow, causing greater risk of water scarcity, groundwater overdraft, and fire in areas that rely on mountain snowpack for their water supply. Streamflow in large river basins varies with the amount, timing, and type of precipitation, evapotranspiration, and drainage properties of watersheds; however, these controls vary in time and space making it difficult to identify the areas contributing most to flow and when. In this study, we separate the evaporative influences from source values of water isotopes from the Snake River Basin in the western United States (US) to relate source area to flow dynamics. We developed isoscapes (δ 2 H and δ 18 O) for the basin and found that isotopic composition of surface water in small watersheds is primarily controlled by longitude, latitude, and elevation. To examine temporal variability in source contributions to flow, we present a six-year record of Snake River water isotopes from King Hill, Idaho after removing evaporative influences. During periods of low flow, source water values were isotopically lighter indicating a larger contribution to flow from surface waters in the highest elevation, eastern portion of the basin. River evaporation increases were evident during summer likely reflecting climate, changing water availability, and management strategies within the basin. Our findings present a potential tool for identifying critical portions of basins contributing to river flow as climate fluctuations alter flow dynamics. This tool can be applied in other continental-interior basins where evaporation may obscure source water isotopic signatures.

Idaho, Montana, Oregon, Wyoming