M. Bayani Cardenas

Water table dynamics and groundwater–surface water interaction during filling and draining of a large fluvial island due to dam‐induced river stage fluctuations

[1] Dam‐controlled river stage fluctuations alter groundwater–surface water interaction between persistent bars and islands and the rivers bounding them by rapidly changing hydraulic gradients and expanding hyporheic zones. A 300‐m long and 80‐m wide sand‐ gravel island with established vegetation located on the Colorado River (Austin, Texas, USA) is subjected to >1 m daily river stage variations due to upstream dam operations. Piezometer nests with probes monitored the evolution of the water table and groundwater flow paths through several cycles of dam‐induced stage fluctuations. Results show that hydraulic head and the water table within the island closely track the river stage associated with dam release. Water table mounds and depressions which overlap in time were mapped through the course of one storage‐release cycle over which >4,000 m 3 of water moved in and out of the island. Dam operations have drastically altered groundwater–surface water connectivity between the Colorado River and the fluvial island aquifer by pumping substantial amounts of water in and out of the aquifer during dam release and storage cycles.

Groundwater flow, transport, and residence times through topography-driven basins with exponentially decreasing permeability and porosity.

Received 27 March 2010; revised 29 June 2010; accepted 9 August 2010; published 20 November 2010. [1] In this paper, we investigate the effects of systematic and local heterogeneity on groundwater flow, transport, and residence time distributions (RTDs) of basins where groundwater flow is topography driven. Systematic heterogeneity is represented by an exponentially depth‐decreasing hydraulic conductivity and porosity, and local heterogeneity is represented by the dispersivity. The RTDs for both a simple basin with one flow system and a basin with nested local and regional systems gradually evolve to a power law RTD with more pronounced systematic heterogeneity. Exponential decrease of poromechanical properties enhances shallow circulation and subdues deep and regional flows leading to longer flushing times for the large part of the domain, while the shallower portions flush solutes rapidly. Therefore, deeper basins lead to more persistent and pronounced power law RTDs when the poromechanical properties systematically decrease with depth. Separate contributions to the RTD due to stagnation zones associated with local flow cells and due to deeper immobile zones were identified; each leads to a different tailing behavior. Local heterogeneity slightly enhances the power law RTD by causing the tailing to begin earlier but does not affect the late time portion of the RTD. Systematic depth‐dependent heterogeneity is an important factor controlling the circulation and associated RTDs of subsurface fluids. It contributes significantly to generation of power law RTDs.

Hyporheic zone hydrologic science: A historical account of its emergence and a prospectus

The hyporheic zone, defined by shallow subsurface pathways through river beds and banks beginning and ending at the river, is an integral and unique component of fluvial systems. It hosts myriad hydrologically controlled processes that are potentially coupled in complex ways. Understanding these processes and the connections between them is critical since these processes are not only important locally but integrate to impact increasingly larger scale biogeochemical functioning of the river corridor up to the river network scale. Thus, the hyporheic zone continues to be a growing research focus for many hydrologists for more than half the history of Water Resources Research. This manuscript partly summarizes the historical development of hyporheic zone hydrologic science as gleaned from papers published in Water Resources Research, from the birth of the concept of the hyporheic zone as a hydrologic black box (sometimes referred to as transient storage zone), to its adolescent years of being torn between occasionally competing research perspectives of interrogating the hyporheic zone from a surface or subsurface view, to its mature emergence as an interdisciplinary research field that employs the wide array of state-of-the-art tools available to the modern hydrologist. The field is vibrant and moving in the right direction of addressing critical fundamental and applied questions with no clear end in sight in its growth. There are exciting opportunities for scientists that are able to tightly link the allied fields of geology, geomorphology, hydrology, geochemistry, and ecology to tackle the many open problems in hyporheic zone science.

Geoelectrical imaging of hyporheic exchange and mixing of river water and groundwater in a large regulated river.

Hyporheic mixing and surface water-groundwater interactions are critical processes in aquatic environments. Yet, there is a lack of methods for assessing the spatial extent and distribution of these mixing zones. This study applied time-lapse electrical resistivity (ER) imaging in a 60-m wide and 0.7-m deep alluvial river whose stage periodically varied by 0.7 m due to dam operations to assess dynamic hyporheic mixing and surface water-groundwater interactions. Sixteen channel-spanning repeat ER tomograms (2D sections) over one flood cycle captured the dynamic ER distribution. We mapped a laterally discontinuous hyporheic zone, which had mainly river water circulating through it, several meters into the bed. Underneath the hyporheic zone was a transitional mixing zone intermittently flushed by mixing river water and deep groundwater. Minimally mixed groundwater dominated the deepest areas. ER imaging allows for unraveling hyporheic and deep mixing zone dynamics in large regulated rivers.

Effects of Multiscale Anisotropy on Basin and Hyporheic Groundwater Flow

Various subsurface flow systems exhibit a combination of small-scale to large-scale anisotropy in hydraulic conductivity (K). The large-scale anisotropy results from systematic trends (e.g., exponential decrease or increase) of K with depth. We present a general two-dimensional solution for calculation of topography-driven groundwater flow considering both small- and large-scale anisotropy in K. This solution can be applied to diverse systems with arbitrary head distribution and geometry of the water table boundary, such as basin or hyporheic flow. In a special case, this solution reduces to the well-known Tóth model of uniform isotropic basin. We introduce an integral measure of flushing intensity that quantifies flushing at different depths. Using this solution, we simulate heads and streamlines and provide analyses of flow structure in the flow domain, relevant to basin analyses or hyporheic flow. It is shown that interactions between small-scale anisotropy and large-scale anisotropy strongly control the flow structure. In the classic Tóth flow model, the flushing intensity curves exhibit quasi-exponential decrease with depth. The new measure is capable of capturing subtle changes in the flow structure. Our study shows that both small- and large-scale anisotropy characteristics have substantial effects that need to be integrated into analysis of topography-driven flow.

Ex-Stream: A MATLAB program for calculating fluid flux through sediment-water interfaces based on steady and transient temperature profiles

Temperature is a useful environmental tracer for quantifying movement and exchange of water and heat through and near sediment-water interfaces (SWI). Heat tracing involves analyzing temperature time series or profiles from temperature probes deployed in sediments. Ex-Stream is a MATLAB program that brings together two transient and two steady one-dimensional coupled heat and fluid flux analytical models. The program includes a graphical user interface, a detailed user manual, and postprocessing capabilities that enable users to extract fluid fluxes from time-series temperature observations. Program output is written to comma-separated values files, displayed within the MATLAB command window, and may be optionally plotted. The models that are integrated into Ex-Stream can be run collectively, allowing for direct comparison, or individually.

Transport zonation limits coupled nitrification-denitrification in permeable sediments.

Measurement of biogeochemical processes in permeable sediments (including the hyporheic zone) is difficult because of complex multidimensional advective transport. This is especially the case for nitrogen cycling, which involves several coupled redox-sensitive reactions. To provide detailed insight into the coupling between ammonification, nitrification and denitrification in stationary sand ripples, we combined the diffusion equilibrium thin layer (DET) gel technique with a computational reactive transport biogeochemical model. The former approach provided high-resolution two-dimensional distributions of NO3(-) and (15)N-N2 gas. The measured two-dimensional profiles correlate with computational model simulations, showing a deep pool of N2 gas forming, and being advected to the surface below ripple peaks. Further isotope pairing calculations on these data indicate that coupled nitrification-denitrification is severely limited in permeable sediments because the flow and transport field limits interaction between oxic and anoxic pore water. The approach allowed for new detailed insight into subsurface denitrification zones in complex permeable sediments.

Modification of the Local Cubic Law of fracture flow for weak inertia, tortuosity, and roughness

The classical Local Cubic Law (LCL) generally overestimates flow through real fractures. We thus developed and tested a modified LCL (MLCL) which takes into account local tortuosity and roughness, and works across a low range of local Reynolds Numbers. The MLCL is based on (1) modifying the aperture field by orienting it with the flow direction and (2) correcting for local roughness changes associated with local flow expansion/contraction. In order to test the MLCL, we compared it with direct numerical simulations with the Navier-Stokes equations using real and synthetic three-dimensional rough-walled fractures, previous corrected forms of the LCL, and experimental flow tests. The MLCL performed well and the effective errors (δ) in volumetric flow rate range from −3.4% to 13.4% with an arithmetic mean of |δ| ( ) equal to 3.7%. The MLCL is more accurate than previous modifications of the LCL. We also investigated the error associated with applying the Cubic Law (CL) while utilizing modified aperture field. The δ from the CL ranges from −14.2% to 11.2%, with a slightly higher  = 6.1% than the MLCL. The CL with the modified aperture field considering local tortuosity and roughness may also be sufficient for predicting the hydraulic properties of rough fractures.

Non‐Fickian transport through two‐dimensional rough fractures: Assessment and prediction

Non-Fickian transport ubiquitously occurs across all scales within fractured geological media. Detailed characterization of non-Fickian transport through single fractures is thus critical for predicting the fate of solutes and other fluid-borne entities through fractured media. Our direct numerical simulations of solute transport through two-dimensional rough-walled fractures showed early arrival and heavy tailing in breakthrough curves (BTCs), which are salient characteristics of non-Fickian transport. Analyses for dispersion coefficients (DADE) using the standard advection-dispersion equation (ADE) led to errors which increased linearly with fracture heterogeneity. Estimated Taylor dispersion coefficients deviated from estimated DADE even at higher Peclet numbers. Alternatively, we used continuous time random walk (CTRW) model with truncated power law transition rate probability to characterize the non-Fickian transport. CTRW modeling markedly and consistently improved fits to the BTCs relative to those fitted with ADE solutions. The degree of deviation of transport from Fickian to non-Fickian is captured by the parameter β of the truncated power law. We found that β is proportional to fracture heterogeneity. We also found that the CTRW transport velocity can be predicted based on the flow velocity. Along with the ability to predict β, this is a major step toward prediction of transport through CTRW using measurable physical properties.

Wave‐driven porewater and solute circulation through rippled elastic sediment under highly transient forcing

Waves induce porewater flow and solute transport through permeable marine sediment. However, past studies have ignored high-frequency pressure pulses, under the assumption that the porewater flow field is adequately represented by a time-averaged one or that the saturated sediment is incompressible. We modeled porewater flow and solute transport inside ripples, forced by instantaneous pressure profiles along the sediment-water interface (SWI) with 0.1-s temporal resolution. The transient pressure profiles were taken from a field data–driven large-eddy simulation model of wave-driven oscillatory flow. The simulations suggest that in elastic, permeable, and saturated sediment, a time-averaged representation of the flow field may be inadequate and that this also leads to shortcomings in how transport is modeled. Bursts in fluid flushing occur when high-frequency pressure fluctuations were considered, leading to larger long-term average fluid fluxes compared to a steady flow field driven by a time-averaged pressure profile. The pressure perturbations along the SWI propagate within a few milliseconds to meter depths within the sediment leading to strongly transient porewater velocity fields. This leads to enhanced dispersion of solutes and larger time-averaged solute fluxes. However, enhanced solute flux across the SWI diminished through time with increasing permeability. The high-frequency transient pressures and sediment elastic properties we considered have been largely ignored and unrecognized. Future observational and modeling studies should consider these processes, especially since they mediate timing-sensitive biogeochemical reactions.