Adam S. Ward

A field comparison of multiple techniques to quantify groundwater-surface-water interactions

Groundwater–surface-water (GW-SW) interactions in streams are difficult to quantify because of heterogeneity in hydraulic and reactive processes across a range of spatial and temporal scales. The challenge of quantifying these interactions has led to the development of several techniques, from centimeter-scale probes to whole-system tracers, including chemical, thermal, and electrical methods. We co-applied conservative and smart reactive solute-tracer tests, measurement of hydraulic heads, distributed temperature sensing, vertical profiles of solute tracer and temperature in the stream bed, and electrical resistivity imaging in a 450-m reach of a 3rd-order stream. GW-SW interactions were not spatially expansive, but were high in flux through a shallow hyporheic zone surrounding the reach. NaCl and resazurin tracers suggested different surface–subsurface exchange patterns in the upper ⅔ and lower ⅓ of the reach. Subsurface sampling of tracers and vertical thermal profiles quantified relatively high fluxes through a 10- to 20-cm deep hyporheic zone with chemical reactivity of the resazurin tracer indicated at 3-, 6-, and 9-cm sampling depths. Monitoring of hydraulic gradients along transects with MINIPOINT streambed samplers starting ∼40 m from the stream indicated that groundwater discharge prevented development of a larger hyporheic zone, which progressively decreased from the stream thalweg toward the banks. Distributed temperature sensing did not detect extensive inflow of ground water to the stream, and electrical resistivity imaging showed limited large-scale hyporheic exchange. We recommend choosing technique(s) based on: 1) clear definition of the questions to be addressed (physical, biological, or chemical processes), 2) explicit identification of the spatial and temporal scales to be covered and those required to provide an appropriate context for interpretation, and 3) maximizing generation of mechanistic understanding and reducing costs of implementing multiple techniques through collaborative research.

How does subsurface characterization affect simulations of hyporheic exchange

We investigated the role of increasingly well-constrained geologic structures in the subsurface (i.e., subsurface architecture) in predicting streambed flux and hyporheic residence time distribution (RTD) for a headwater stream. Five subsurface realizations with increasingly resolved lithological boundaries were simulated in which model geometries were based on increasing information about flow and transport using soil and geologic maps, surface observations, probing to depth to refusal, seismic refraction, electrical resistivity (ER) imaging of subsurface architecture, and time-lapse ER imaging during a solute tracer study. Particle tracking was used to generate RTDs for each model run. We demonstrate how improved characterization of complex lithological boundaries and calibration of porosity and hydraulic conductivity affect model prediction of hyporheic flow and transport. Models using hydraulic conductivity calibrated using transient ER data yield estimates of streambed flux that are three orders of magnitude larger than uncalibrated models using estimated values for hydraulic conductivity based on values published for nearby hillslopes (10(-4) vs. 10(-7) m(2)/s, respectively). Median residence times for uncalibrated and calibrated models are 10(3) and 10(0) h, respectively. Increasingly well-resolved subsurface architectures yield wider hyporheic RTDs, indicative of more complex hyporheic flowpath networks and potentially important to biogeochemical cycling. The use of ER imaging to monitor solute tracers informs subsurface structure not apparent from other techniques, and helps to define transport properties of the subsurface (i.e., hydraulic conductivity). Results of this study demonstrate the value of geophysical measurements to more realistically simulate flow and transport along hyporheic flowpaths.

Hydrogeomorphic controls on hyporheic and riparian transport in two headwater mountain streams during base flow recession

Solute transport along riparian and hyporheic flow paths is broadly expected to respond to dynamic hydrologic forcing by streams, aquifers, and hillslopes. However, direct observation of these dynamic responses is lacking, as is the relative control of geologic setting as a control on responses to dynamic hydrologic forcing. We conducted a series of four stream solute tracer injections through base flow recession in each of two watersheds with contrasting valley morphology in the H.J. Andrews Experimental Forest, monitoring tracer concentrations in the stream and in a network of shallow riparian wells in each watershed. We found hyporheic mean arrival time, temporal variance, and fraction of stream water in the bedrock-constrained valley bottom and near large roughness elements in the wider valley bottom were not variable with discharge, suggesting minimal control by hydrologic forcing. Conversely, we observed increases in mean arrival time and temporal variance and decreasing fraction stream water with decreasing discharge near the hillslopes in the wider valley bottom. This may indicate changes in stream discharge and valley bottom hydrology control transport in less constrained locations. We detail five hydrogeomorphic responses to base flow recession to explain observed spatial and temporal patterns in the interactions between streams and their valley bottoms. Models able to account for the transition from geologically dominated processes in the near-stream subsurface to hydrologically dominated processes near the hillslope will be required to predict solute transport and fate in valley bottoms of headwater mountain streams.

Antecedent Moisture Controls on Stream Nitrate Flux in an Agricultural Watershed

Evaluating nitrate-N fluxes from agricultural landscapes is inherently complex due to the wide range of intrinsic and dynamic controlling variables. In this study, we investigate the influence of contrasting antecedent moisture conditions on nitrate-N flux magnitude and dynamics in a single agricultural watershed on intra-annual and rainfall-event temporal scales. High temporal resolution discharge and nitrate concentration data were collected to evaluate nitrate-N flux magnitude associated with wet (2009) and dry (2012) conditions. Analysis of individual rainfall events revealed a marked and consistent difference in nitrate-N flux response attributed to wet/dry cycles. Large-magnitude dilutions (up to 10 mg N L) persisted during the wet antecedent conditions (2009), consistent with a dominant baseflow contribution and excess groundwater release in relation to precipitation volume (discharge > > precipitation). Smaller-magnitude concentrations (<7 mg N L) were observed during the drought conditions of 2012, consistent with a quickflow-dominated response to rain events and infiltration/storage of precipitation resulting in discharge < precipitation. Nitrate-N loads and yields from the watershed were much higher (up to an order of magnitude) in the wet year vs. the dry year. Our results suggest that the response of nitrate-N loading to rain events is highly dependent on intra-annual antecedent moisture conditions and subsurface hydrologic connectivity, which together dictate the dominant hydrologic pathways for stream recharge. Additionally, the results of our study indicate that continued pronounced wet/dry cycles may become more dominant as the short-term driver of future nitrate-N exports.

Response of the hyporheic zone to transient groundwater fluctuations on the annual and storm event time scales

The volume of the water stored in and exchanged with the hyporheic zone is an important factor in stream metabolism and biogeochemical cycling. Previous studies have identified groundwater direction and magnitude as one key control on the volume of the hyporheic zone, suggesting that fluctuation in the riparian water table could induce large changes under certain seasonal conditions. In this study, we analyze the transient drivers that control the volume of the hyporheic zone by coupling the Brinkman-Darcy equation to the Navier-Stokes equations to simulate annual and storm induced groundwater fluctuations. The expansion and contraction of the hyporheic zone was quantified based on temporally dynamic scenarios simulating annual groundwater fluctuations in a humid temperate climate. The amplitude of the groundwater signal was varied between scenarios to represent a range of annual hydrologic forcing. Storm scenarios were then superimposed on the annual scenario to simulate the response to short-term storm signals. Simulations used two different groundwater storm responses; one in-phase with the surface water response and one 14 h out-of-phase with the surface water response to represent our observed site conditions. Results show that annual groundwater fluctuation is a dominant control on the volume of the hyporheic zone, where increasing groundwater fluctuation increases the amount of annual variation. Storm responses depended on the antecedent conditions determined by annual scenarios, where the time of year dictated the duration and magnitude of the storm induced response of the hyporheic zone.

Stream solute tracer timescales changing with discharge and reach length confound process interpretation

Improved understanding of stream solute transport requires meaningful comparison of processes across a wide range of discharge conditions and spatial scales. At reach scales where solute tracer tests are commonly used to assess transport behavior, such comparison is still confounded due to the challenge of separating dispersive and transient storage processes from the influence of the advective timescale that varies with discharge and reach length. To better resolve interpretation of these processes from field-based tracer observations, we conducted recurrent conservative solute tracer tests along a 1 km study reach during a storm discharge period and further discretized the study reach into six segments of similar length but different channel morphologies. The resulting suite of data, spanning an order of magnitude in advective timescales, enabled us to (1) characterize relationships between tracer response and discharge in individual segments and (2) determine how combining the segments into longer reaches influences interpretation of dispersion and transient storage from tracer tests. We found that the advective timescale was the primary control on the shape of the observed tracer response. Most segments responded similarly to discharge, implying that the influence of morphologic heterogeneity was muted relative to advection. Comparison of tracer data across combined segments demonstrated that increased advective timescales could be misinterpreted as a change in dispersion or transient storage. Taken together, our results stress the importance of characterizing the influence of changing advective timescales on solute tracer responses before such reach-scale observations can be used to infer solute transport at larger network scales.

A Comparison of Hyporheic Transport at a Cross‐Vane Structure and Natural Riffle

While restoring hyporheic flowpaths has been cited as a benefit to stream restoration structures, little documentation exists confirming that constructed restoration structures induce comparable hyporheic exchange to natural stream features. This study compares a stream restoration structure (cross-vane) to a natural feature (riffle) concurrently in the same stream reach using time-lapsed electrical resistivity (ER) tomography. Using this hydrogeophysical approach, we were able to quantify hyporheic extent and transport beneath the cross-vane structure and the riffle. We interpret from the geophysical data that the cross-vane and the natural riffle induced spatially and temporally unique hyporheic extent and transport, and the cross-vane created both spatially larger and temporally longer hyporheic flowpaths than the natural riffle. Tracer from the 4.67-h injection was detected along flowpaths for 4.6 h at the cross-vane and 4.2 h at the riffle. The spatial extent of the hyporheic zone at the cross-vane was 12% larger than that at the riffle. We compare ER results of this study to vertical fluxes calculated from temperature profiles and conclude significant differences in the interpretation of hyporheic transport from these different field techniques. Results of this study demonstrate a high degree of heterogeneity in transport metrics at both the cross-vane and the riffle and differences between the hyporheic flowpath networks at the two different features. Our results suggest that restoration structures may be capable of creating sufficient exchange flux and timescales of transport to achieve the same ecological functions as natural features, but engineering of the physical and biogeochemical environment may be necessary to realize these benefits.

Dynamic hyporheic and riparian flow path geometry through base flow recession in two headwater mountain stream corridors

The hydrologic connectivity between streams and their valley bottoms (stream corridor) is a critical determinant of their ecological function. Ecological functions are known to be spatially and temporally variable, but spatial dimensions of the problem are not easily quantified and thus they are usually overlooked. To estimate the spatial patterns of connectivity, and how connectivity varies with changes in discharge, we developed the hyporheic potential model. We used the model to interpret a series of solute tracer injections in two headwater mountain streams with contrasting valley bottom morphologies to estimate connectivity in the stream corridor. The distributions of flow path origination locations and the lengths of hyporheic flow paths appear to vary with base flow recession, even in cases where transport timescales are apparently unchanged. The modeled distribution of origination locations further allowed us to define a spatial analog to the temporal window of detection associated with solute tracer studies, and enables assessment of connectivity dynamics between streams and their corridors. Altogether, the reduced complexity hyporheic potential model provides an easy way to anticipate the spatial distribution and origination locations of hyporheic flow paths from a basic understanding of the valley bottom characteristics and solute transport timescales. Plain Language Summary The manuscript details a simple method to assess the spatial connectivity of streams and their riparian zones. While the timescales of exchange in the river corridor have been broadly studied, the complimentary spatial dimension (i.e., the geometry of exchange flowpaths) remains largely unknown. The major challenge in assessing the spatial dimensions of exchange is the limited information available in the subsurface. Here, we develop a reduced complexity model of valley bottom transport to overcome these information limitations. With this model, relatively simple field site characterization and solute tracer data are combined to assess the spatial distribution of downwelling along a headwater mountain stream. We validate the model with a numerical experiment, and demonstrate its application in two watersheds of contrasting geology, repeated through baseflow recession.

How does reach‐scale stream‐hyporheic transport vary with discharge? Insights from rSAS analysis of sequential tracer injections in a headwater mountain stream

The models of stream reach hyporheic exchange that are typically used to interpret tracer data assume steady-flow conditions and impose further assumptions about transport processes on the interpretation of the data. Here we show how rank Storage Selection (rSAS) functions can be used to extract ‘‘process-agnostic’’ information from tracer breakthrough curves about the time-varying turnover of reach storage. A sequence of seven slug injections was introduced to a small stream at base flow over the course of a diel fluctuation in stream discharge, providing breakthrough curves at discharges ranging from 0.7 to 1.2 L/s. Shifted gamma distributions, each with three parameters varying stepwise in time, were used to model the rSAS function and calibrated to reproduce each breakthrough curve with Nash-Sutcliffe efficiencies in excess of 0.99. Variations in the fitted parameters over time suggested that storage within the reach does not uniformly increase its turnover rate when discharge increases. Rather, changes in transit time are driven by both changes in the average rate of turnover (external variability) and changes in the relative rate that younger and older water contribute to discharge (internal variability). Specifically, at higher discharge, the turnover rate increased for the youngest part of the storage (corresponding to approximately 5 times the volume of the channel), while discharge from the older part of the storage remained steady, or declined slightly. The method is shown to be extensible as a new approach to modeling reach-scale solute transport that accounts for the time-varying, discharge-dependent turnover of reach storage.

Optimizing Sampling Strategies for Riverine Nitrate Using High-Frequency Data in Agricultural Watersheds

Understanding linked hydrologic and biogeochemical processes such as nitrate loading to agricultural streams requires that the sampling bias and precision of monitoring strategies be known. An existing spatially distributed, high-frequency nitrate monitoring network covering ∼40% of Iowa provided direct observations of in situ nitrate concentrations at a temporal resolution of 15 min. Systematic subsampling of nitrate records allowed for quantification of uncertainties (bias and precision) associated with estimates of various nitrate parameters, including: mean nitrate concentration, proportion of samples exceeding the nitrate drinking water standard (DWS), peak (>90th quantile) nitrate concentration, and nitrate flux. We subsampled continuous records for 47 site-year combinations mimicking common, but labor-intensive, water-sampling regimes (e.g., time-interval, stage-triggered, and dynamic-discharge storm sampling). Our results suggest that time-interval sampling most efficiently characterized all nitrate parameters, except at coarse frequencies for nitrate flux. Stage-triggered storm sampling most precisely captured nitrate flux when less than 0.19% of possible 15 min observations for a site-year were used. The time-interval strategy had the greatest return on sampling investment by most precisely and accurately quantifying nitrate parameters per sampling effort. These uncertainty estimates can aid in designing sampling strategies focused on nitrate monitoring in the tile-drained Midwest or similar agricultural regions.