USGS Science⌕ Search

SEARCH · USGS Science

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 559 records · Page 31Linked to original sources
USGS2012

Exploring changes in the spatial distribution of stream baseflow generation during a seasonal recession

Relating watershed structure to streamflow generation is a primary focus of hydrology. However, comparisons of longitudinal variability in stream discharge with adjacent valley structure have been rare, resulting in poor understanding of the distribution of the hydrologic mechanisms that cause variability in streamflow generation along valleys. This study explores detailed surveys of stream base flow across a gauged, 23 km 2 mountain watershed. Research objectives were (1) to relate spatial variability in base flow to fundamental elements of watershed structure, primarily topographic contributing area, and (2) to assess temporal changes in the spatial patterns of those relationships during a seasonal base flow recession. We analyzed spatiotemporal variability in base flow using (1) summer hydrographs at the study watershed outlet and 5 subwatershed outlets and (2) longitudinal series of discharge measurements every ~100 m along the streams of the 3 largest subwatersheds (1200 to 2600 m in valley length), repeated 2 to 3 times during base flow recession. Reaches within valley segments of 300 to 1200 m in length tended to demonstrate similar streamflow generation characteristics. Locations of transitions between these segments were consistent throughout the recession, and tended to be collocated with abrupt longitudinal transitions in valley slope or hillslope-riparian characteristics. Both within and among subwatersheds, correlation between the spatial distributions of streamflow and topographic contributing area decreased during the recession, suggesting a general decrease in the influence of topography on stream base flow contributions. As topographic controls on base flow evidently decreased, multiple aspects of subsurface structure were likely to have gained influence.

Water Resources Research↗
USGS2022

Integrated modeling of dynamic marsh feedbacks and evolution under sea-level rise in a mesotidal estuary (Plum Island, MA, USA)

Around the world, wetland vulnerability to sea-level rise (SLR) depends on different factors including tidal regimes, topography, creeks and estuary geometry, sediment availability, vegetation type, etc. The Plum Island estuary (PIE) is a mesotidal wetland system on the east coast of the United States. This research applied a newly updated Hydro-MEM (integrated hydrodynamic-marsh) model to assess the impacts of intermediate-low (50 cm), intermediate (1 m), and intermediate-high (1.5 m) SLR on marsh evolution by the year 2100. Model advancements include capturing vegetation change, inorganic and below and aboveground organic matter portion of marsh platform accretion, and mudflat creation. Although the results indicate a low vulnerability marsh at the PIE, the vegetation changes from high to low marsh under all SLR scenarios (2%–22%), with the higher bounds belonging to higher rise scenarios. Lower SLR produces more productive marsh (13% gain in high productivity regions), whereas the highest SLR scenario causes increased tidal inundation, which leads to loss in productivity (12% change from high to low productivity regions), generation of mudflats (17% of the domain land), and marsh migration to higher lands. Sensitive nonlinear tidal flow changes, which may be increased or decreased with SLR as a result of mudflat creation, marsh migration, and bottom friction change, emphasize the importance of integrated modeling approaches that include dynamic marsh feedbacks in hydrodynamic modeling and varying hydrodynamic effects on the marsh system.

Massachusetts↗
USGS1998

Interactions between pool geometry and hydraulics

An experimental and computational research approach was used to determine interactions between pool geometry and hydraulics. A 20-m-long, 1.8-m-wide flume was used to investigate the effect of four different geometric aspects of pool shape on flow velocity. Plywood sections were used to systematically alter constriction width, pool depth, pool length, and pool exit-slope gradient, each at two separate levels. Using the resulting 16 unique geometries with measured pool velocities in four-way factorial analyses produced an empirical assessment of the role of the four geometric aspects on the pool flow patterns and hence the stability of the pool. To complement the conclusions of these analyses, a two-dimensional computational flow model was used to investigate the relationships between pool geometry and flow patterns over a wider range of conditions. Both experimental and computational results show that constriction and depth effects dominate in the jet section of the pool and that pool length exhibits an increasing effect within the recirculating-eddy system. The pool exit slope appears to force flow reattachment. Pool length controls recirculating-eddy length and vena contracta strength. In turn, the vena contracta and recirculating eddy control velocities throughout the pool.

Water Resources Research↗
USGS2010

Modeling biotic uptake by periphyton and transient hyporrheic storage of nitrate in a natural stream

To a convection-dispersion hydrologic transport model we coupled a transient storage submodel (Bencala, 1984) and a biotic uptake submodel based on Michaelis-Menten kinetics (Kim et al., 1990). Our purpose was threefold: (1) to simulate nitrate retention in response to change in load in a third-order stream, (2) to differentiate biotic versus hydrologie factors in nitrate retention, and (3) to produce a research tool whose properties are consistent with laboratory and field observations. Hydrodynamic parameters were fitted from chloride concentration during a 20-day chloride-nitrate coinjection (Bencala, 1984), and biotic uptake kinetics were based on flume studies by Kim et al. (1990) and Triska et al. (1983). Nitrate concentration from the 20-day coinjection experiment served as a base for model validation. The complete transport retention model reasonably predicted the observed nitrate concentration. However, simulations which lacked either the transient storage submodel or the biotic uptake submodel poorly predicted the observed nitrate concentration. Model simulations indicated that transient storage in channel and hyporrheic interstices dominated nitrate retention within the first 24 hours, whereas biotic uptake dominated thereafter. A sawtooth function for V max ranging from 0.10 to 0.17 μg NO 3 -N s −1 gAFDM −1 (grams ash free dry mass) slightly underpredicted nitrate retention in simulations of 2–7 days. This result was reasonable since uptake by other nitrate-demanding processes were not included. The model demonstrated how ecosystem retention is an interaction between physical and biotic processes and supports the validity of coupling separate hydrodynamic and reactive submodels to established solute transport models in biological studies of fluvial ecosystems.

California↗
USGS2008

Analytical and numerical analyses of an unconfined aquifer test considering unsaturated zone characteristics

A 7-d, constant rate aquifer test conducted by University of Waterloo researchers at Canadian Forces Base Borden in Ontario, Canada, is useful for advancing understanding of fluid flow processes in response to pumping from an unconfined aquifer. Measured data include not only drawdown in the saturated zone but also volumetric soil moisture measured at various times and distances from the pumped well. Analytical analyses were conducted with the model published in 2001 by Moench and colleagues, which allows for gradual drainage but does not include unsaturated zone characteristics, and the model published in 2006 by Mathias and Butler, which assumes that moisture retention and relative hydraulic conductivity (RHC) in the unsaturated zone are exponential functions of pressure head. Parameters estimated with either model yield good matches between measured and simulated drawdowns in piezometers. Numerical analyses were conducted with two versions of VS2DT: one that uses traditional Brooks and Corey functional relations and one that uses a RHC function introduced in 2001 by Assouline that includes an additional parameter that accounts for soil structure and texture. The analytical model of Mathias and Butler and numerical model of VS2DT with the Assouline model both show that the RHC function must contain a fitting parameter that is different from that used in the moisture retention function. Results show the influence of field-scale heterogeneity and suggest that the RHC at the Borden site declines more rapidly with elevation above the top of the capillary fringe than would be expected if the parameters were to reflect local- or core-scale soil structure and texture.

Water Resources Research↗
USGS2004

Estimating accumulation rates and physical properties of sediment behind a dam: Englebright Lake, Yuba River, northern California

Studies of reservoir sedimentation are vital to understanding scientific and management issues related to watershed sediment budgets, depositional processes, reservoir operations, and dam decommissioning. Here we quantify the mass, organic content, and grain-size distribution of a reservoir deposit in northern California by two methods of extrapolating measurements of sediment physical properties from cores to the entire volume of impounded material. Englebright Dam, completed in 1940, is located on the Yuba River in the Sierra Nevada foothills. A research program is underway to assess the feasibility of introducing wild anadromous fish species to the river upstream of the dam. Possible management scenarios include removing or lowering the dam, which could cause downstream transport of stored sediment. In 2001 the volume of sediments deposited behind Englebright Dam occupied 25.5% of the original reservoir capacity. The physical properties of this deposit were calculated using data from a coring campaign that sampled the entire reservoir sediment thickness (6–32 m) at six locations in the downstream ∼3/4 of the reservoir. As a result, the sediment in the downstream part of the reservoir is well characterized, but in the coarse, upstream part of the reservoir, only surficial sediments were sampled, so calculations there are more uncertain. Extrapolation from one-dimensional vertical sections of sediment sampled in cores to entire three-dimensional volumes of the reservoir deposit is accomplished via two methods, using assumptions of variable and constant layer thickness. Overall, the two extrapolation methods yield nearly identical estimates of the mass of the reservoir deposit of ∼26 × 10 6 metric tons (t) of material, of which 64.7–68.5% is sand and gravel. Over the 61 year reservoir history this corresponds to a maximum basin-wide sediment yield of ∼340 t/km 2 /yr, assuming no contribution from upstream parts of the watershed impounded by other dams. The uncertainties and limitations of the estimates of overall sediment quantities are discussed. Implications for watershed management and future reservoir sedimentation studies are also presented.

California↗
USGS2000

Regional interdisciplinary paleoflood approach to assess extreme flood potential

In the past decade, there has been a growing interest of dam safety officials to incorporate a risk‐based analysis for design‐flood hydrology. Extreme or rare floods, with probabilities in the range of about 10 −3 to 10 −7 chance of occurrence per year, are of continuing interest to the hydrologic and engineering communities for purposes of planning and design of structures such as dams [ National Research Council , 1988]. The National Research Council stresses that as much information as possible about floods needs to be used for evaluation of the risk and consequences of any decision. A regional interdisciplinary paleoflood approach was developed to assist dam safety officials and floodplain managers in their assessments of the risk of large floods. The interdisciplinary components included documenting maximum paleofloods and a regional analyses of contemporary extreme rainfall and flood data to complement a site‐specific probable maximum precipitation study [ Tomlinson and Solak , 1997]. The cost‐effective approach, which can be used in many other hydrometeorologic settings, was applied to Elkhead Reservoir in Elkhead Creek (531 km 2 ) in northwestern Colorado; the regional study area was 10,900 km 2 . Paleoflood data using bouldery flood deposits and noninundation surfaces for 88 streams were used to document maximum flood discharges that have occurred during the Holocene. Several relative dating methods were used to determine the age of paleoflood deposits and noninundation surfaces. No evidence of substantial flooding was found in the study area. The maximum paleoflood of 135 m 3 s −1 for Elkhead Creek is about 13% of the site‐specific probable maximum flood of 1020 m 3 s −1 . Flood‐frequency relations using the expected moments algorithm, which better incorporates paleoflood data, were developed to assess the risk of extreme floods. Envelope curves encompassing maximum rainfall (181 sites) and floods (218 sites) were developed for northwestern Colorado to help define maximum contemporary and Holocene flooding in Elkhead Creek and in a regional frequency context. Study results for Elkhead Reservoir were accepted by the Colorado State Engineer for dam safety certification.

Water Resources Research↗
USGS2009

Recent experimental data may point to a greater role for osmotic pressures in the subsurface

Uncertainty about the origin of anomalous fluid pressures in certain geologic settings has caused researchers to take a second look at osmosis, or flow driven by chemical potential differences, as a pressure‐generating process in the subsurface. Interest in geological osmosis has also increased because of an in situ experiment by Neuzil (2000) suggesting that Pierre Shale could generate large osmotic pressures when highly compacted. In the last few years, additional laboratory and in situ experiments have greatly increased the number of data on osmotic properties of argillaceous formations, but they have not been systematically examined. In this paper we compile these data and explore their implications for osmotic pressure generation in subsurface systems. Rather than base our analysis on osmotic efficiencies, which depend strongly on concentration, we calculated values of a quantity we term osmotic specific surface area ( A so ) that, in principle, is a property of the porous medium only. The A so values are consistent with a surprisingly broad spectrum of osmotic behavior in argillaceous formations, and all the formations tested exhibited at least a modest ability to generate osmotic pressure. It appears possible that under appropriate conditions some formations can be highly effective osmotic membranes able to generate osmotic pressures exceeding 30 MPa (3 km of head) at porosities as high as ∼0.1 and pressures exceeding 10 MPa at porosities as high as ∼0.2. These findings are difficult to reconcile with the lack of compelling field evidence for osmotic pressures, and we propose three explanations for the disparity: (1) Our analysis is flawed and argillaceous formations are less effective osmotic membranes than it suggests; (2) the necessary subsurface conditions, significant salinity differences within intact argillaceous formations, are rare; or (3) osmotic pressures are unlikely to be detected and are not recognized when encountered. The last possibility, that osmotic pressures routinely escape detection or are attributed to other mechanisms, has important implications for understanding subsurface flow regimes.

Water Resources Research↗
USGS2010

Comparison of episodic acidification of Mid-Atlantic Upland and Coastal Plain streams

Episodic acidification was examined in five mid-Atlantic watersheds representing three physiographic provinces: Coastal Plain, Valley and Ridge, and Blue Ridge. Each of the watersheds receives a similar loading of atmospheric pollutants (SO 4 2− and NO 3 − ) and is underlain by different bedrock type. The purpose of this research was to quantify and compare the episodic variability in storm flow chemistry in Reedy Creek, Virginia (Coastal Plain), Mill Run and Shelter Run, Virginia (Valley and Ridge), and Fishing Creek Tributary and Hunting Creek, Maryland (Blue Ridge). Because episodic responses were similar from storm to storm in each of the watersheds, a representative storm from each watershed was discussed. Acidification, defined as the loss of acid-neutralizing capacity (ANC), was observed in all streams except Mill Run. Mill Run chemistry showed little episodic variability. During storms in the other streams, p H decreased while SO 4 2− , NO 3 − , and K + concentrations increased. Concentrations of Mg 2+ and Ca 2+ increased in Reedy Creek and Fishing Creek Tributary, but decreased in Shelter Run and Hunting Creek. Therefore the net effect of episodic changes on the acid-base status differed among the streams. In general, greater losses of ANC were observed during storms at Shelter Run and Hunting Creek, watersheds underlain by reactive bedrock (carbonate, metabasalt); comparatively smaller losses in ANC were observed at Reedy Creek and Fishing Creek Tributary, watersheds underlain by quartzites and unconsolidated quartz sands and cobbles. Increased SO 4 2− concentrations were most important during storms at Reedy Creek and Fishing Creek Tributary, but organic anions (inferred by anion deficit) were also a factor in causing the loss of ANC. Dilution of base cations was the most important factor in the loss of ANC at Shelter Run. Both increased sulfate and dilution of base flow were important in causing the episodic acidification at Hunting Creek. The role of SO 4 2− in contributing to episodic acidification in these watersheds is similar to that documented in studies conducted in other regions of the United States, Scandinavia, Canada, and the United Kingdom. The importance of SO 4 2− in mid-Atlantic United States streams contrasts with northeastern United States streams, in which increased NO 3 - derived from snowpack is more important in causing episodic acidification. Results support the importance of shallow subsurface processes in determining storm flow chemistry, regional climatic characteristics in determining the different sources of acidity during episodes, and the importance of bedrock geology in determining the amount of ANC loss.

Maryland, Virginia, West Virginia↗
USGS2011

Nutrient loadings to streams of the continental United States from municipal and industrial effluent?

Data from the United States Environmental Protection Agency Permit Compliance System national database were used to calculate annual total nitrogen (TN) and total phosphorus (TP) loads to surface waters from municipal and industrial facilities in six major regions of the United States for 1992, 1997, and 2002. Concentration and effluent flow data were examined for approximately 118,250 facilities in 45 states and the District of Columbia. Inconsistent and incomplete discharge locations, effluent flows, and effluent nutrient concentrations limited the use of these data for calculating nutrient loads. More concentrations were reported for major facilities, those discharging more than 1 million gallons per day, than for minor facilities, and more concentrations were reported for TP than for TN. Analytical methods to check and improve the quality of the Permit Compliance System data were used. Annual loads were calculated using "typical pollutant concentrations" to supplement missing concentrations based on the type and size of facilities. Annual nutrient loads for over 26,600 facilities were calculated for at least one of the three years. Sewage systems represented 74% of all TN loads and 58% of all TP loads. This work represents an initial set of data to develop a comprehensive and consistent national database of point-source nutrient loads. These loads can be used to inform a wide range of water-quality management, watershed modeling, and research efforts at multiple scales.

Journal of the American Water Resources Associatio↗
USGS2014

Coupled hydrological and biogeochemical processes controlling variability of nitrogen species in streamflow during autumn in an upland forest

Autumn is a season of dynamic change in forest streams of the northeastern United States due to effects of leaf fall on both hydrology and biogeochemistry. Few studies have explored how interactions of biogeochemical transformations, various nitrogen sources, and catchment flow paths affect stream nitrogen variation during autumn. To provide more information on this critical period, we studied (1) the timing, duration, and magnitude of changes to stream nitrate, dissolved organic nitrogen (DON), and ammonium concentrations; (2) changes in nitrate sources and cycling; and (3) source areas of the landscape that most influence stream nitrogen. We collected samples at higher temporal resolution for a longer duration than typical studies of stream nitrogen during autumn. This sampling scheme encompassed the patterns and extremes that occurred during base flow and stormflow events of autumn. Base flow nitrate concentrations decreased by an order of magnitude from 5.4 to 0.7 µmol L −1 during the week when most leaves fell from deciduous trees. Changes to rates of biogeochemical transformations during autumn base flow explained the low nitrate concentrations; in-stream transformations retained up to 72% of the nitrate that entered a stream reach. A decrease of in-stream nitrification coupled with heterotrophic nitrate cycling were primary factors in the seasonal nitrate decline. The period of low nitrate concentrations ended with a storm event in which stream nitrate concentrations increased by 25-fold. In the ensuing weeks, peak stormflow nitrate concentrations progressively decreased over closely spaced, yet similarly sized events. Most stormflow nitrate originated from nitrification in near-stream areas with occasional, large inputs of unprocessed atmospheric nitrate, which has rarely been reported for nonsnowmelt events. A maximum input of 33% unprocessed atmospheric nitrate to the stream occurred during one event. Large inputs of unprocessed atmospheric nitrate show direct and rapid effects on forest streams that may be widespread, although undocumented, throughout nitrogen-polluted temperate forests. In contrast to a week-long nitrate decline during peak autumn litterfall, base flow DON concentrations increased after leaf fall and remained high for 2 months. Dissolved organic nitrogen was hydrologically flushed to the stream from riparian soils during stormflow. In contrast to distinct seasonal changes in base flow nitrate and DON concentrations, ammonium concentrations were typically at or below the detection limit, similar to the rest of the year. Our findings reveal couplings among catchment flow paths, nutrient sources, and transformations that control seasonal extremes of stream nitrogen in forested landscapes.

Vermont↗
USGS2014

Assessment of parametric uncertainty for groundwater reactive transport modeling,

The validity of using Gaussian assumptions for model residuals in uncertainty quantification of a groundwater reactive transport model was evaluated in this study. Least squares regression methods explicitly assume Gaussian residuals, and the assumption leads to Gaussian likelihood functions, model parameters, and model predictions. While the Bayesian methods do not explicitly require the Gaussian assumption, Gaussian residuals are widely used. This paper shows that the residuals of the reactive transport model are non-Gaussian, heteroscedastic, and correlated in time; characterizing them requires using a generalized likelihood function such as the formal generalized likelihood function developed by Schoups and Vrugt (2010). For the surface complexation model considered in this study for simulating uranium reactive transport in groundwater, parametric uncertainty is quantified using the least squares regression methods and Bayesian methods with both Gaussian and formal generalized likelihood functions. While the least squares methods and Bayesian methods with Gaussian likelihood function produce similar Gaussian parameter distributions, the parameter distributions of Bayesian uncertainty quantification using the formal generalized likelihood function are non-Gaussian. In addition, predictive performance of formal generalized likelihood function is superior to that of least squares regression and Bayesian methods with Gaussian likelihood function. The Bayesian uncertainty quantification is conducted using the differential evolution adaptive metropolis (DREAM (zs) ) algorithm; as a Markov chain Monte Carlo (MCMC) method, it is a robust tool for quantifying uncertainty in groundwater reactive transport models. For the surface complexation model, the regression-based local sensitivity analysis and Morris- and DREAM (ZS) -based global sensitivity analysis yield almost identical ranking of parameter importance. The uncertainty analysis may help select appropriate likelihood functions, improve model calibration, and reduce predictive uncertainty in other groundwater reactive transport and environmental modeling.

Water Resources Research↗
USGS2013

Effect of land cover and use on dry season river runoff, runoff efficiency, and peak storm runoff in the seasonal tropics of Central Panama

A paired catchment methodology was used with more than 3 years of data to test whether forests increase base flow in the dry season, despite reduced annual runoff caused by evapotranspiration (the “sponge-effect hypothesis”), and whether forests reduce maximum runoff rates and totals during storms. The three study catchments were: a 142.3 ha old secondary forest, a 175.6 ha mosaic of mixed age forest, pasture, and subsistence agriculture, and a 35.9 ha actively grazed pasture subcatchment of the mosaic catchment. The two larger catchments are adjacent, with similar morphology, soils, underlying geology, and rainfall. Annual water balances, peak runoff rates, runoff efficiencies, and dry season recessions show significant differences. Dry season runoff from the forested catchment receded more slowly than from the mosaic and pasture catchments. The runoff rate from the forest catchment was 1–50% greater than that from the similarly sized mosaic catchment at the end of the dry season. This observation supports the sponge-effect hypothesis. The pasture and mosaic catchment median runoff efficiencies were 2.7 and 1.8 times that of the forest catchment, respectively, and increased with total storm rainfall. Peak runoff rates from the pasture and mosaic catchments were 1.7 and 1.4 times those of the forest catchment, respectively. The forest catchment produced 35% less total runoff and smaller peak runoff rates during the flood of record in the Panama Canal Watershed. Flood peak reduction and increased streamflows through dry periods are important benefits relevant to watershed management, payment for ecosystem services, water-quality management, reservoir sedimentation, and fresh water security in the Panama Canal watershed and similar tropical landscapes.

Water Resources Research↗
USGS2014

Distributed Evaluation of Local Sensitivity Analysis (DELSA), with application to hydrologic models

This paper presents a hybrid local-global sensitivity analysis method termed the Distributed Evaluation of Local Sensitivity Analysis (DELSA), which is used here to identify important and unimportant parameters and evaluate how model parameter importance changes as parameter values change. DELSA uses derivative-based “local” methods to obtain the distribution of parameter sensitivity across the parameter space, which promotes consideration of sensitivity analysis results in the context of simulated dynamics. This work presents DELSA, discusses how it relates to existing methods, and uses two hydrologic test cases to compare its performance with the popular global, variance-based Sobol' method. The first test case is a simple nonlinear reservoir model with two parameters. The second test case involves five alternative “bucket-style” hydrologic models with up to 14 parameters applied to a medium-sized catchment (200 km 2 ) in the Belgian Ardennes. Results show that in both examples, Sobol' and DELSA identify similar important and unimportant parameters, with DELSA enabling more detailed insight at much lower computational cost. For example, in the real-world problem the time delay in runoff is the most important parameter in all models, but DELSA shows that for about 20% of parameter sets it is not important at all and alternative mechanisms and parameters dominate. Moreover, the time delay was identified as important in regions producing poor model fits, whereas other parameters were identified as more important in regions of the parameter space producing better model fits. The ability to understand how parameter importance varies through parameter space is critical to inform decisions about, for example, additional data collection and model development. The ability to perform such analyses with modest computational requirements provides exciting opportunities to evaluate complicated models as well as many alternative models.

Water Resources Research↗
USGS2014

Lithologic influences on groundwater recharge through incised glacial till from profile to regional scales: Evidence from glaciated Eastern Nebraska

Variability in sediment hydraulic properties associated with landscape depositional and erosional features can influence groundwater recharge processes by affecting soil-water storage and transmission. This study considers recharge to aquifers underlying river-incised glaciated terrain where the distribution of clay-rich till is largely intact in upland locations but has been removed by alluvial erosion in stream valleys. In a stream-dissected glacial region in eastern Nebraska (Great Plains region of the United States), recharge estimates were developed for nested profile, aquifer, and regional scales using unsaturated zone profile measurements (matric potentials, Cl- and 3H), groundwater tracers (CFC-12 and SF6), and a remote sensing-assisted water balance model. Results show a consistent influence of till lithology on recharge rates across nested spatial scales despite substantial uncertainty in all recharge estimation methods, suggesting that minimal diffuse recharge occurs through upland glacial till lithology whereas diffuse recharge occurs in river valleys where till is locally absent. Diffuse recharge is estimated to account for a maximum of 61% of total recharge based on comparison of diffuse recharge estimated from the unsaturated zone (0-43 mm yr-1) and total recharge estimated from groundwater tracers (median 58 mm yr-1) and water balance modeling (median 56 mm yr-1). The results underscore the importance of lithologic controls on the distributions of both recharge rates and mechanisms.

Nebraska↗
USGS2014

Dynamic hyporheic exchange at intermediate timescales: testing the relative importance of evapotranspiration and flood pulses

Hyporheic fluxes influence ecological processes across a continuum of timescales. However, few studies have been able to characterize hyporheic fluxes and residence time distributions (RTDs) over timescales of days to years, during which evapotranspiration (ET) and seasonal flood pulses create unsteady forcing. Here we present a data-driven, particle-tracking piston model that characterizes hyporheic fluxes and RTDs based on measured vertical head differences. We used the model to test the relative influence of ET and seasonal flood pulses in the Everglades (FL, USA), in a manner applicable to other low-energy floodplains or broad, shallow streams. We found that over the multiyear timescale, flood pulses that drive relatively deep (∼1 m) flow paths had the dominant influence on hyporheic fluxes and residence times but that ET effects were discernible at shorter timescales (weeks to months) as a break in RTDs. Cumulative RTDs on either side of the break were generally well represented by lognormal functions, except for when ET was strong and none of the standard distributions applied to the shorter timescale. At the monthly timescale, ET increased hyporheic fluxes by 1–2 orders of magnitude; it also decreased 6 year mean residence times by 53–87%. Long, slow flow paths driven by flood pulses increased 6 year hyporheic fluxes by another 1–2 orders of magnitude, to a level comparable to that induced over the short term by shear flow in streams. Results suggest that models of intermediate-timescale processes should include at least two-storage zones with different RTDs, and that supporting field data collection occur over 3–4 years.

Florida↗
USGS2014

Precipitation isotopes link regional climate patterns to water supply in a tropical mountain forest, eastern Puerto Rico

Like many mountainous areas in the tropics, watersheds in the Luquillo Mountains of eastern Puerto Rico have abundant rainfall and stream discharge and provide much of the water supply for the densely populated metropolitan areas nearby. Projected changes in regional temperature and atmospheric dynamics as a result of global warming suggest that water availability will be affected by changes in rainfall patterns. It is essential to understand the relative importance of different weather systems to water supply to determine how changes in rainfall patterns, interacting with geology and vegetation, will affect the water balance. To help determine the links between climate and water availability, stable isotope signatures of precipitation from different weather systems were established to identify those that are most important in maintaining streamflow and groundwater recharge. Precipitation stable isotope values in the Luquillo Mountains had a large range, from fog/cloud water with δ 2 H, δ 18 O values as high as +12 ‰, −0.73 ‰ to tropical storm rain with values as low as −127 ‰, −16.8 ‰. Temporal isotope values exhibit a reverse seasonality from those observed in higher latitude continental watersheds, with higher isotopic values in the winter and lower values in the summer. Despite the higher volume of convective and low-pressure system rainfall, stable isotope analyses indicated that under the current rainfall regime, frequent trade -wind orographic showers contribute much of the groundwater recharge and stream base flow. Analysis of rain events using 20 years of 15 -minute resolution data at a mountain station (643 m) showed an increasing trend in rainfall amount, in agreement with increased precipitable water in the atmosphere, but differing from climate model projections of drying in the region. The mean intensity of rain events also showed an increasing trend. The determination of recharge sources from stable isotope tracers indicates that water supply will be affected if regional atmospheric dynamics change trade- wind orographic rainfall patterns in the Caribbean.

Luquillo Mountains↗
USGS2014

Watershed-scale modeling of streamflow change in incised montane meadows

Land use practices have caused stream channel incision and water table decline in many montane meadows of the Western United States. Incision changes the magnitude and timing of streamflow in water supply source watersheds, a concern to resource managers and downstream water users. The hydrology of montane meadows under natural and incised conditions was investigated using watershed simulation for a range of hydrologic conditions. The results illustrate the interdependence between: watershed and meadow hydrology; bedrock and meadow aquifers; and surface and groundwater flow through the meadow for the modeled scenarios. During the wet season, stream incision resulted in less overland flow and interflow and more meadow recharge causing a net decrease in streamflow and increase in groundwater storage relative to natural meadow conditions. During the dry season, incision resulted in less meadow evapotranspiration and more groundwater discharge to the stream causing a net increase in streamflow and a decrease in groundwater storage relative to natural meadow conditions. In general, for a given meadow setting, the magnitude of change in summer streamflow and long-term change in watershed groundwater storage due to incision will depend on the combined effect of: reduced evapotranspiration in the eroded meadow; induced groundwater recharge; replenishment of dry season groundwater storage depletion in meadow and bedrock aquifers by precipitation during wet years; and groundwater storage depletion that is not replenished by precipitation during wet years.

Water Resources Research↗