Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Water Resources 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 145 records · Page 8Linked to original sources
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↗
USGS2003

An analytical formulation of two‐dimensional groundwater dispersion induced by surficial recharge variability

A predominant cause of dispersion in groundwater is advective mixing due to variability in seepage rates. Hydraulic conductivity variations have been extensively researched as a cause of this seepage variability. In this paper the effect of variations in surface recharge to a shallow surficial aquifer is investigated as an important additional effect. An analytical formulation has been developed that relates aquifer parameters and the statistics of recharge variability to increases in the dispersivity. This is accomplished by solving Fourier transforms of the small perturbation forms of the groundwater flow equations. Two field studies are presented in this paper to determine the statistics of recharge variability for input to the analytical formulation. A time series of water levels at a continuous groundwater recorder is used to investigate the temporal statistics of hydraulic head caused by recharge, and a series of infiltrometer measurements are used to define the spatial variability in the recharge parameters. With these field statistics representing head fluctuations due to recharge, the analytical formulation can be used to compute the dispersivity without an explicit representation of the recharge boundary. Results from a series of numerical experiments are used to define the limits of this analytical formulation and to provide some comparison. A sophisticated model has been developed using a particle‐tracking algorithm (modified to account for temporal variations) to estimate groundwater dispersion. Dispersivity increases of 9 percent are indicated by the analytical formulation for the aquifer at the field site. A comparison with numerical model results indicates that the analytical results are reasonable for shallow surficial aquifers in which two‐dimensional flow can be assumed.

Water Resources Research↗
USGS2004

Scale effects on headwater catchment runoff timing, flow sources, and groundwater‐streamflow relations

The effects of catchment size and landscape organization on runoff generation are poorly understood. Little research has integrated hillslope and riparian runoff investigation across catchments of different sizes to decipher first‐order controls on runoff generation. We investigated the role of catchment sizes on riparian and hillslope dynamics based on hydrometric and tracer data observed at five scales ranging from trenched hillslope sections (55–285 m 2 ) to a 280‐ha catchment at Maimai on the west coast of the South Island, New Zealand. The highly organized landscape is comprised of similar headwater catchments, regular geology, steep highly dissected topography, relatively consistent soil depths, and topographically controlled shallow through flow. We found a strong correlation between riparian zone groundwater levels and runoff for the headwaters, whereas the water tables in the valley bottom of the larger catchments were uncorrelated to runoff for 14 months of record. While there was no clear relationship between catchment size and new water contribution to runoff in the two storms analyzed in detail, lag times of tracer responses increased systematically with catchment size. The combination of hydrometric and tracer data allowed assessment of the runoff contributions from different parts of the landscape. Runoff was generated consistently in headwater riparian zones. This agreed also with the observed variations of tracer ( 18 O and silica) responses for the different catchments. During wetter antecedent conditions or during larger events (>30 mm under dry antecedent conditions) hillslope and valley bottom floodplains did contribute to event runoff directly. We propose that analysis of landscape‐scale organization and the distribution of dominant landscape features provide a structure for investigation of runoff production and solute transport, especially as catchment‐scale increases from headwaters to the mesoscale.

Water Resources Research↗
USGS2005

CO2 dynamics in the Amargosa Desert: Fluxes and isotopic speciation in a deep unsaturated zone

Natural unsaturated-zone gas profiles at the U.S. Geological Survey's Amargosa Desert Research Site, near Beatty, Nevada, reveal the presence of two physically and isotopically distinct CO 2 sources, one shallow and one deep. The shallow source derives from seasonally variable autotrophic and heterotrophic respiration in the root zone. Scanning electron micrograph results indicate that at least part of the deep CO 2 source is associated with calcite precipitation at the 110-m-deep water table. We use a geochemical gas-diffusion model to explore processes of CO 2 production and behavior in the unsaturated zone. The individual isotopic species 12 CO 2 , 13 CO 2 , and 14 CO 2 are treated as separate chemical components that diffuse and react independently. Steady state model solutions, constrained by the measured δ 13 C (in CO 2 ), and δ 14 C (in CO 2 ) profiles, indicate that the shallow CO 2 source from root and microbial respiration composes ∼97% of the annual average total CO 2 production at this arid site. Despite the small contribution from deep CO 2 production amounting to ∼0.1 mol m −2 yr −1 , upward diffusion from depth strongly influences the distribution of CO 2 and carbon isotopes in the deep unsaturated zone. In addition to diffusion from deep CO 2 production, 14 C exchange with a sorbed CO 2 phase is indicated by the modeled δ 14 C profiles, confirming previous work. The new model of carbon-isotopic profiles provides a quantitative approach for evaluating fluxes of carbon under natural conditions in deep unsaturated zones.

Water Resources Research↗
USGS2019

Hillslope hydrology in global change research and earth system modeling

Earth System Models (ESMs) are essential tools for understanding and predicting global change, but they cannot explicitly resolve hillslope‐scale terrain structures that fundamentally organize water, energy, and biogeochemical stores and fluxes at subgrid scales. Here we bring together hydrologists, Critical Zone scientists, and ESM developers, to explore how hillslope structures may modulate ESM grid‐level water, energy, and biogeochemical fluxes. In contrast to the one‐dimensional (1‐D), 2‐ to 3‐m deep, and free‐draining soil hydrology in most ESM land models, we hypothesize that 3‐D, lateral ridge‐to‐valley flow through shallow and deep paths and insolation contrasts between sunny and shady slopes are the top two globally quantifiable organizers of water and energy (and vegetation) within an ESM grid cell. We hypothesize that these two processes are likely to impact ESM predictions where (and when) water and/or energy are limiting. We further hypothesize that, if implemented in ESM land models, these processes will increase simulated continental water storage and residence time, buffering terrestrial ecosystems against seasonal and interannual droughts. We explore efficient ways to capture these mechanisms in ESMs and identify critical knowledge gaps preventing us from scaling up hillslope to global processes. One such gap is our extremely limited knowledge of the subsurface, where water is stored (supporting vegetation) and released to stream baseflow (supporting aquatic ecosystems). We conclude with a set of organizing hypotheses and a call for global syntheses activities and model experiments to assess the impact of hillslope hydrology on global change predictions.

Water Resources Research↗
USGS2011

Effects of uncertain topographic input data on two-dimensional flow modeling in a gravel-bed river

Many applications in river research and management rely upon two-dimensional (2D) numerical models to characterize flow fields, assess habitat conditions, and evaluate channel stability. Predictions from such models are potentially highly uncertain due to the uncertainty associated with the topographic data provided as input. This study used a spatial stochastic simulation strategy to examine the effects of topographic uncertainty on flow modeling. Many, equally likely bed elevation realizations for a simple meander bend were generated and propagated through a typical 2D model to produce distributions of water-surface elevation, depth, velocity, and boundary shear stress at each node of the model's computational grid. Ensemble summary statistics were used to characterize the uncertainty associated with these predictions and to examine the spatial structure of this uncertainty in relation to channel morphology. Simulations conditioned to different data configurations indicated that model predictions became increasingly uncertain as the spacing between surveyed cross sections increased. Model sensitivity to topographic uncertainty was greater for base flow conditions than for a higher, subbankfull flow (75% of bankfull discharge). The degree of sensitivity also varied spatially throughout the bend, with the greatest uncertainty occurring over the point bar where the flow field was influenced by topographic steering effects. Uncertain topography can therefore introduce significant uncertainty to analyses of habitat suitability and bed mobility based on flow model output. In the presence of such uncertainty, the results of these studies are most appropriately represented in probabilistic terms using distributions of model predictions derived from a series of topographic realizations. Copyright 2011 by the American Geophysical Union.

Water Resources Research↗
USGS1998

Base cation concentrations in subsurface flow from a forested hillslope: The role of flushing frequency

A 20-m-wide trench was excavated to bedrock on a hillslope at the Panola Mountain Research Watershed in the Piedmont region of Georgia to determine the effect of upslope drainage area from the soil and bedrock surfaces on the geochemical evolution of base cation concentrations in subsurface flow. Samples were collected from ten 2-m sections and five natural soil pipes during three winter rainstorms in 1996. Base cation concentrations in hillslope subsurface flow were generally highest early and late in the storm response when flow rates were low, but during peak flow, concentrations varied little. Base cation concentrations in matrix flow from the 10 trench sections were unrelated to the soil surface drainage area and weakly inversely related to the bedrock surface drainage area. Base cation concentrations in pipe flow were lower than those in matrix flow and were also consistent with the inverse relation to bedrock surface drainage area found in matrix flow. The left side of the trench, which has the highest bedrock surface drainage area, had consistently lower mean base cation concentrations than the right side of the trench, which has the lowest bedrock surface drainage area. During moderate size rain events of about 20–40 mm, subsurface flow occurred only on the left side of the trench. The greater volume of water that has flowed through the left side of the trench appears to have resulted in greater leaching of base cations from soils and therefore lower base cation concentrations in subsurface flow than in flow from the right side of the trench. Alternatively, a greater proportion of flow that bypasses the soil matrix may have occurred through the hillslope on the left side of the trench than on the right side. Flushing frequency links spatial hillslope water flux with the evolution of groundwater and soil chemistry.

Georgia↗
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↗
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↗