Numerical Modeling of Nitrate Removal in Anoxic Groundwater during River Flooding of Riparian Zones

A numerical study demonstrates the effects of flooding on subsurface hydrological flowpaths and nitrate removal in anoxic groundwater in riparian zones with a top peat layer. A series of two-dimensional numerical simulations with changing conditions for flow (steady state or transient with flooding), hydrogeology, denitrification, and duration of flooding demonstrate how flowpaths, residence times, and nitrate removal are affected. In periods with no flooding groundwater flows horizontally and discharges to the river through the riverbed. During periods with flooding, shallow groundwater is forced upwards as discharge through peat layers that often have more optimal conditions for denitrification caused by the presence of highly reactive organic matter. The contrast in hydraulic conductivity between the sand aquifer and the overlying peat layer, as well as the flooding duration, have a significant role in determining the degree of nitrate removal.     

An Application of HFE-D for Evaluating Sea Water Intrusion in Coastal Aquifers of Southern Vietnam

The coastal aquifers and inland waters of the Long Xuyen Quadrangle and Ca Mau Peninsula of southern Vietnam have been significantly impacted by sea water intrusion (SI) as a result of recent anthropogenic activities. This study identified the evolution and spatial distribution of hydrochemical conditions in coastal aquifers at this region using Hydrochemical Facies Evolution Diagram (HFE-D) and Geographical Information System mapping. Hydraulic heads and water chemistry were measured at 31 observation wells in four layered aquifers during dry and rainy seasons in early (2005), and more recent (2016), stages of agricultural development. Hydrochemical facies associated with intrusion or freshening stages were mapped in each aquifer after assigning mixing index values to each facies. The position of groundwater freshening and SI phases differed in Holocene, Upper Pleistocene, Middle Pleistocene, and Lower Pleistocene aquifers. The geographic position of freshening and intrusion fronts differ in dry and rainy seasons, and shifted after 11 years of groundwater abstraction in all four aquifers. The spatial and temporal differences in hydrochemical facies distributions according to HFE-D reflect the relative impact of SI in the four aquifers. The study results provide a better understanding of the evolution of groundwater quality associated with SI in a peninsular coastal aquifer system, and highlight the need for improving groundwater quality and management in similar coastal regions.     

Geomorphologic control on pollutant retardation at the groundwater–surface water interface

The results of research on the pollutant retardation potential of permeable riverbed sediments in catchments with significant groundwater–surface water (GW-SW) interaction are presented. The fraction of organic carbon and cation exchange capacity of fluvial sediments in various geomorphologic environments have been quantified. Sediments in selected reaches of the rivers Tern and Leith (UK), from the underlying Permian sandstone aquifers, and from along the length of the rivers Severn and Eden into which the Tern and Leith discharge have been investigated. Statistical analyses show significant variation in the geochemistry and pollutant retardation potential of sediments from different geomorphologic features, and between upland and lowland rivers. The sorption potential of fine-grained sediments deposited in pools was greater than sand in runs and coarser deposits in riffles. Similarly, sediments in lowland rivers were found to have a greater retardation potential than those in upland rivers. There was generally greater retardation potential in fluvial sediments of all types than in the underlying aquifers, and in lowland rivers the fluvial sediment retardation potential greatly dominated that of the aquifer. The findings demonstrate the potential for pollutant retardation processes in riverbed sediments of sandstone catchments, and suggest that consideration of retardation processes at the groundwater–surface water interface should be included into environmental risk-assessment studies, in order to better assess and manage the effects of contaminated groundwater discharges to rivers, particularly in lowland catchments. Copyright © 2008 John Wiley & Sons, Ltd.     

River-aquifer interactions, geologic heterogeneity, and low-flow management

Low river flows are commonly controlled by river-aquifer exchange, the magnitude of which is governed by hydraulic properties of both aquifer and aquitard materials beneath the river. Low flows are often important ecologically. Numerical simulations were used to assess how textural heterogeneity of an alluvial system influences river seepage and low flows. The Cosumnes River in California was used as a test case. Declining fall flows in the Cosumnes River have threatened Chinook salmon runs. A ground water-surface water model for the lower river basin was developed, which incorporates detailed geostatistical simulations of aquifer heterogeneity. Six different realizations of heterogeneity and a homogenous model were run for a 3-year period. Net annual seepage from the river was found to be similar among the models. However, spatial distribution of seepage along the channel, water table configuration and the level of local connection, and disconnection between the river and aquifer showed strong variations among the different heterogeneous models. Most importantly, the heterogeneous models suggest that river seepage losses can be reduced by local reconnections, even when the regional water table remains well below the riverbed. The percentage of river channel responsible for 50% of total river seepage ranged from 10% to 26% in the heterogeneous models as opposed to 23% in the homogeneous model. Differences in seepage between the models resulted in up to 13 d difference in the number of days the river was open for salmon migration during the critical fall months in one given year.     

Assimilation of historical head data to estimate spatial distributions of stream bed and aquifer hydraulic conductivity fields

Management of water resources in alluvial aquifers relies mainly on understanding interactions between hydraulically connected streams and aquifers. Numerical models that simulate this interaction often are used as decision support tools for water resource management. However, the accuracy of numerical predictions relies heavily on unknown system parameters (e.g., streambed conductivity and aquifer hydraulic conductivity), which are spatially heterogeneous and difficult to measure directly. This paper employs an ensemble smoother to invert groundwater level measurements to jointly estimate spatially varying streambed and alluvial aquifer hydraulic conductivity along a 35.6‐km segment of the South Platte River in Northeastern Colorado. The accuracy of the inversion procedure is evaluated using a synthetic experiment and historical groundwater level measurements, with the latter constituting the novelty of this study in the inversion and validation of high‐resolution fields of streambed and aquifer conductivities. Results show that the estimated streambed conductivity field and aquifer conductivity field produce an acceptable agreement between observed and simulated groundwater levels and stream flow rates. The estimated parameter fields are also used to simulate the spatially varying flow exchange between the alluvial aquifer and the stream, which exhibits high spatial variability along the river reach with a maximum average monthly aquifer gain of about 2.3 m³/day and a maximum average monthly aquifer loss of 2.8 m³/day, per unit area of streambed (m²). These results demonstrate that data assimilation inversion provides a reliable and computationally affordable tool to estimate the spatial variability of streambed and aquifer conductivities at high resolution in real‐world systems.     

Surface water-groundwater exchange dynamics in buried-valley aquifer systems

Groundwater is a primary source of drinking water worldwide, but excess nutrients and emerging contaminants could compromise groundwater quality and limit its usage as a drinking water source. As such contaminants become increasingly prevalent in the biosphere, a fundamental understanding of their fate and transport in groundwater systems is necessary to implement successful remediation strategies. The dynamics of surface water-groundwater (hyporheic) exchange within a glacial, buried-valley aquifer system are examined in the context of their implications for the transport of nutrients and contaminants in riparian sediments. High conductivity facies act as preferential flow pathways which enhance nutrient and contaminant delivery, especially during storm events, but transport throughout the aquifer also depends on subsurface sedimentary architecture (e.g. interbedded high and low conductivity facies). Temperature and specific conductance measurements indicate extensive hyporheic mixing close to the river channel, but surface water influence was also observed far from the stream-aquifer interface. Measurements of river stage and hydraulic head indicate that significant flows during storms (i.e., hot moments) alter groundwater flow patterns, even between consecutive storm events, as riverbed conductivity and, more importantly, the hydraulic connectivity between the river and aquifer change. Given the similar mass transport characteristics among buried-valley aquifers, these findings are likely representative of glacial aquifer systems worldwide. Our results suggest that water resources management decisions based on average (base) flow conditions may inaccurately represent the system being evaluated, and could reduce the effectiveness of remediation strategies for nutrients and emerging contaminants.     

Depth‐dependent hydraulic conductivity distribution patterns of a streambed

Characterization of streambed hydraulic conductivity from the channel surface to a great depth below the channel surface can provide needed information for the determination of stream‐aquifer hydrologic connectedness, and it is also important to river restoration. However, knowledge on the streambed hydraulic conductivity for sediments 1 m below the channel surface is scarce. This study describes a method that was used to determine the distribution patterns of streambed hydraulic conductivity for sediments from channel surface to a depth of 15 m below. The method includes Geoprobe's direct‐push techniques and Permeameter tests. Direct‐push techniques were used to generate the electrical conductivity (EC) logs and to collect sequences of continuous sediment cores from river channels, as well as from the alluvial aquifer connected to the river. Permeameter tests on these sediment cores give the profiles of vertical hydraulic conductivity (K_v) of the channel sediments and the aquifer materials. This method was applied to produce K_v profiles for a streambed and an alluvial aquifer in the Platte River Valley of Nebraska, USA. Comparison and statistical analysis of the K_v profiles from the river channel and from the proximate alluvial aquifer indicates a special pattern of K_v in the channel sediments. This depth‐dependent pattern of K_v distribution for the channel sediments is considered to be produced by hyporheic processes. This K_v‐distribution pattern implied that the effect of hyporheic processes on streambed hydraulic conductivity can reach the sediments about 9 m below the channel surface. Copyright © 2010 John Wiley & Sons, Ltd.     

Environmental reconstructions of a late devensian terrace sequence. Some preliminary findings

Facies analysis of Severn Main terrace sediments at Eardington, Shropshire is instructive in determining aspects of the sedimentary environment at a ‘proximal’ location in the terrace. Evaluation of the structural and directional properties of the gravels indicates that the terrace is a composite feature comprising two units. The lower unit is notable for the widespread occurrence of large sandy facies, a near absence of massive gravel facies, and the presence of large, locally derived, lithoclasts in excess of 1 m in diameter. Facies associations observed indicate the existence of in-channel bars and large channels with bedforms at the dune-plane bed transition. The directional properties of the unit indicate flows parallel to the main valley axis and although the unit is unlike previously reported examples it is interpreted as being the product of a proximal low sinuosity environment. The upper unit is composed mainly of multistorey, coarse, massive gravel units with limited interbedded sand facies. Fabric analysis indicates progradation of the unit from the west. The facies sequence is similar to those encompassed by facies models proposed for proximal alluvial fan deposits and therefore, on structural and directional grounds, the upper unit is interpreted as an alluvial fan deposit which prograded into the main Severn valley from the adjacent Mor Brook tributary. The implications of the composite nature of the terrace at this point are considered both in terms of terrace correlation and palaeohyrological estimation.     

Heterogeneity influences on stream water–groundwater interactions in a gravel-dominated floodplain

ABSTRACT

Floodplains are composed of complex depositional patterns of ancient and recent stream sediments, and research is needed to address the manner in which coarse floodplain materials affect stream–groundwater exchange patterns. Efforts to understand the heterogeneity of aquifers have utilized numerous techniques typically focused on point-scale measurements; however, in highly heterogeneous settings, the ability to model heterogeneity is dependent on the data density and spatial distribution. The objective of this research was to investigate the correlation between broad-scale methodologies for detecting heterogeneity and the observed spatial variability in stream/groundwater interactions of gravel-dominated alluvial floodplains. More specifically, this study examined the correlation between electrical resistivity (ER) and alluvial groundwater patterns during a flood event at a site on Barren Fork Creek, in the Ozark ecoregion of Oklahoma, USA, where chert gravels were common both as streambed and as floodplain material. Water table elevations from groundwater monitoring wells for a flood event on 1–5 May 2009 were compared to ER maps at various elevations. Areas with high ER matched areas with lower water table slope at the same elevation. This research demonstrated that ER approaches were capable of indicating heterogeneity in surface water–groundwater interactions, and that these heterogeneities were present even in an aquifer matrix characterized as highly conductive. Portions of gravel-dominated floodplain vadose zones characterized by high hydraulic conductivity features can result in heterogeneous flow patterns when the vadose zone of alluvial floodplains activates during storm events.

EDITOR D. Koutsoyiannis; ASSOCIATE EDITOR X. Chen     

Simulation of less‐mobile porosity dynamics in contrasting sediment water interface porous media

Considering heterogeneity in porous media pore size and connectivity is essential to predicting reactive solute transport across interfaces. However, exchange with less‐mobile porosity is rarely considered in surface water/groundwater recharge studies. Previous research indicates that a combination of pore‐fluid sampling and geoelectrical measurements can be used to quantify less‐mobile porosity exchange dynamics using the time‐varying relation between fluid and bulk electrical conductivity. For this study, we use macro‐scale (10 s of cm) advection–dispersion solute transport models linked with electrical conduction in COMSOL Multiphysics to explore less‐mobile porosity dynamics in two different types of observed sediment water interface porous media. Modeled sediment textures contrast from strongly layered streambed deposits to poorly sorted lakebed sands and cobbles. During simulated ionic tracer perturbations, a lag between fluid and bulk electrical conductivity, and the resultant hysteresis, is observed for all simulations indicating differential loading of pore spaces with tracer. Less‐mobile exchange parameters are determined graphically from these tracer time series data without the need for inverse numerical model simulation. In both sediment types, effective less‐mobile porosity exchange parameters are variable in response to changes in flow direction and fluid flux. These observed flow‐dependent effects directly impact local less‐mobile residence times and associated contact time for biogeochemical reaction. The simulations indicate that for the sediment textures explored here, less‐mobile porosity exchange is dominated by variable rates of advection through the domain, rather than diffusion of solute, for typical low‐to‐moderate rate (approximately 3–40 cm/day) hyporheic fluid fluxes. Overall, our model‐based results show that less‐mobile porosity may be expected in a range of natural hyporheic sediments and that changes in flowpath orientation and magnitude will impact less‐mobile exchange parameters. These temporal dynamics can be assessed with the geoelectrical experimental tracer method applied at laboratory and field scales.     

Numerical modelling of riverbed grain size stratigraphic evolution

For several decades,quantification of riverbed grain size stratigraphic evolution has been based upon the active layer formulation（ALF）,which unfortunately involves considerable uncertainty.While it is the sediment exchange across the bed surface that directly affects the riverbed stratigraphy,it has been assumed in the ALF that the sediment fraction at the lower interface of the active layer is a linear function of the sediment fraction in the flow.Here it is proposed that the sediment fraction of the sediment exchange flux is used directly in estimating the sediment fraction at the lower surface of the active layer.Together with the size-specific mass conservation for riverbed sediment,the modified approach is referred to as the surface-based formulation（SBF）.When incorporated into a coupled non-capacity modelling framework for fluvial processes,the SBF leads to results that agree as well or better than those using ALF with laboratory and field observations.This is illustrated for typical cases featuring bed aggradation and degradation due to graded bed-load sediment transport.Systematic experiments on graded sediment transport by unsteady flows are warranted for further testing the modified formulation.     

A Search for Freshwater in the Saline Aquifer of Coastal Bangladesh

In the polder region of coastal Bangladesh, shallow groundwater is primarily brackish with unpredictable occurrence of freshwater pockets. Delta building processes, including the codeposition of fresh-to-saline porewater and sediments, have formed the shallow aquifer. Impermeable clay facies and the lack of a topographical gradient limit the flow of groundwater and its mixing with surface water so controls on spatial variability of salinity are not obvious. By characterizing groundwater-surface water (GW-SW) interactions, this study attempted to identify areas of potable groundwater for the polder communities. We used transects of piezometers, cores, electromagnetic induction, and water chemistry surveys to explore two sources of potential fresh groundwater: (1) tidal channel-aquifer exchange and (2) meteoric recharge. Fresh groundwater proved difficult to find due to heterogeneous subsurface lithology, asymmetrical tidal dynamics, extreme seasonal fluctuations in rainfall, and limited field data. Geophysical observations suggest substantial lateral variability in shallow subsurface conductivity profiles. Piezometers show varying degrees of tidal pressure attenuation away from the channels. Nevertheless, the active exchange of freshwater appears to be limited due to low permeability of banks and surface sediments. Results indicate that pockets of fresh groundwater cannot be identified using readily available hydrogeological methods, so alternative drinking water sources should be pursued. By better understanding the hydrogeology of the system, however, communities will be better equipped to redirect water management resources to more feasible and sustainable drinking water options.     

Groundwater flow modelling of Kwa Ibo River watershed,southeastern Nigeria

Groundwater flow modelling of the Kwa Ibo River watershed in Abia State of Nigeria is presented in this paper with the aim of assessing the degree of interaction between the Kwa Ibo River and the groundwater regime of the thick sandy aquifer. The local geology of the area comprises the Quaternary to recent Benin Formation. Potential aquifer zones that were delineated earlier using geoelectrical resistivity soundings and borehole data for the area formed the basis for groundwater flow modelling. The watershed has been modelled with a grid of 65 rows by 43 columns and with two layers. Lateral inflow from the north has been simulated with constant heads at the Government College, Umuahia, and outflow at Usaka Elegu in the south. The Kwa Ibo River traverses the middle of the watershed from north to south. The river‐stage data at Umudike, Amawom, Ntalakwu and Usaka Elegu have been used for assigning surface water levels and riverbed elevations in the model. Permeability distribution was found to vary from 3 to 14·5 m day^?1. Natural recharge due to rainfall formed the main input to the aquifer system, and abstraction from wells was the main output. A steady‐state groundwater flow simulation was carried out and calibrated against the May 1980 water levels using 26 observation wells. The model computations have converged after 123 iterations. Under the transient‐state calibration, the highest rainfall (and hence groundwater recharge) over the 10‐year study period was recorded in 1996, whereas the lowest was recorded in 1991. The computed groundwater balance of 55 274 m³ day^?1 was comparable to that estimated from field investigations. Results from the modelling show that abstraction is much less than groundwater recharge. Hence there is the possibility for additional groundwater exploitation in the watershed through drilling of boreholes. Copyright © 2008 John Wiley & Sons, Ltd.     

Groundwater flow path determination during riverbank filtration affected by groundwater exploitation: a case study of Liao River,Northeast China

The groundwater flow path plays an important role in maintaining hydrological and ecological quality and security, which are important in the comprehensive management and use of both groundwater and surface water. In this study, an integrated multi-tracer-constrained framework was used to determine the groundwater flow path. The results show that there are shallow and deep flow paths in riverbank filtration, controlled by the different permeabilities of riverbed sediments and aquifers at different depths. The contribution of river water to shallow groundwater is less than that to deep groundwater because of the low permeability of the riverbed sediment in the dense muddy layer in the shallow slope of the river valley. This contribution decreases with increasing distance from the Liao River. The shallow groundwater quality is better than the deep groundwater quality because of its longer residence time.     

Patterns and dynamics of river–aquifer exchange with variably-saturated flow using a fully-coupled model

The parallel physically-based surface–subsurface model PARFLOW was used to investigate the spatial patterns and temporal dynamics of river–aquifer exchange in a heterogeneous alluvial river–aquifer system with deep water table. Aquifer heterogeneity at two scales was incorporated into the model. The architecture of the alluvial hydrofacies was represented based on conditioned geostatistical indicator simulations. Subscale variability of hydraulic conductivities (K) within hydrofacies bodies was created with a parallel Gaussian simulation. The effects of subscale heterogeneity were investigated in a Monte Carlo framework. Dynamics and patterns of river–aquifer exchange were simulated for a 30-day flow event. Simulation results show the rapid formation of saturated connections between the river channel and the deep water table at preferential flow zones that are characterized by high conductivity hydrofacies. Where the river intersects low conductivity hydrofacies shallow perched saturated zones immediately below the river form, but seepage to the deep water table remains unsaturated and seepage rates are low. Preferential flow zones, although only taking up around 50% of the river channel, account for more than 98% of total seepage. Groundwater recharge is most efficiently realized through these zones. Subscale variability of K_sat slightly increased seepage volumes, but did not change the general seepage patterns (preferential flow zones versus perched zones). Overall it is concluded that typical alluvial heterogeneity (hydrofacies architecture) is an important control of river–aquifer exchange in rivers overlying deep water tables. Simulated patterns and dynamics are in line with field observations and results from previous modeling studies using simpler models. Alluvial heterogeneity results in distinct patterns and dynamics of river–aquifer exchange with implications for groundwater recharge and the management of riparian zones (e.g. river channel-floodplain connectivity via saturated zones).     

Coupled groundwater flow and heat transport simulation for estimating transient aquifer–stream exchange at the lowland River Spree (Germany)

Subsurface flow and heat transport near Freienbrink, NE Germany, was simulated in order to study groundwater–surface water exchange between a floodplains aquifer and a section of the lowland River Spree and an adjacent oxbow. Groundwater exfiltration was the dominant process, and only fast surface water level rises resulted in temporary infiltration into the aquifer. The main groundwater flow paths are identified based on a 3D groundwater flow model. To estimate mass fluxes across the aquifer–surface water interfaces, a 2D flow and heat transport modelling approach along a transect of 12 piezometers was performed. Results of steady‐state and transient water level simulations show an overall high accuracy with a Spearman coefficient ρ = 0.9996 and root mean square error (RMSE) = 0.008 m. Based on small groundwater flow velocities of about 10^?7 to 10^?6 ms^?1, mean groundwater exfiltration rates of 233 l m^?2 d^?1 are calculated. Short periods of surface water infiltration into the aquifer do not exceed 10 days, and the infiltration rates are in the same range. The heat transport was modelled with slightly less accuracy (ρ = 0.8359 and RMSE = 0.34 °C). In contrast to the predominant groundwater exfiltration, surface water temperatures determine the calculated temperatures in the upper aquifer below both surface water bodies down to 10 m during the whole simulation period. These findings emphasize prevailing of heat conduction over advection in the upper aquifer zones, which seems to be typical for lowland streams with sandy aquifer materials and low hydraulic gradients. Moreover, this study shows the potential of coupled numerical flow and heat transport modelling to understand groundwater–surface water exchange processes in detail. Copyright © 2013 John Wiley & Sons, Ltd.     

Scaling procedures for heterogeneous unconfined aquifers

Complex flows in heterogeneous confined and unconfined aquifers is a phenomenon that continues to present difficulties in flow mapping and modelling in the field, laboratory, and through numerical simulations. It is often the case with complicated phenomena that transformative scaling and reduction of the problem through symmetry is of great efficacy in the formation of predictive models in both the laboratory and computational settings. A detailed a study of the application of a broad class of Lie scaling transformations on a set of equations representing the groundwater flows in heterogeneous confined and unconfined aquifers has produced a set of scaling relationships between the spatial variables, hydrologic variables, and parameters. The set of scaling transformations preserve the structure of the equations in the sense that the scaling transformations leave the initial‐boundary value system representing the invariant groundwater flows. This theoretical approach elucidates not only the scaling relationships but also the properties that hydrologic variables and parameters must satisfy in order for calling to be possible. Validation of the theory developed is carried out through a series of four numerical simulations using the USGS modflow ‐2005 software package. The results of these experiments demonstrate that the derived scaling transformations can effectively form predictive models of large‐scale phenomena at small scales with negligible error in many cases. Comments on the limitations of the approach and directions for future research are made in the closing sections. Copyright © 2016 John Wiley & Sons, Ltd.     

Investigating Groundwater‐Lake Interactions by Hydraulic Heads and a Water Balance

Discharge of groundwater into lakes (lacustrine groundwater discharge, LGD) can play a major role in water balances of lakes. Unfortunately, studies often neglect this input path because of methodological difficulties in its determination. Direct measurements of LGD are labor‐consuming and prone to error. The present study uses both spatially variable hydraulic‐head data and meteorological data to estimate groundwater input by LGD and lake water output through infiltration. The study sites are two shallow, groundwater‐fed lakes without any surface inflows or outflows. Horizontally interpolated groundwater heads were combined with lake water levels to obtain vertical hydraulic gradients between the aquifer and the lake, which are separated by a thick layer of lake bed sediment which has a much lower hydraulic conductivity than the underlying aquifer. By fitting the hydraulic gradient to the results of a simple mass balance and considering the process of clogging, we were able to estimate the hydraulic conductivity of the lake bed sediments. We calculated groundwater inputs by LGD and lake water outputs by infiltration on an annual basis. Although our method requires several assumptions, the results are reasonable and provide useful information about the exchange between the aquifer and the lake, which can, for example, be used for the calculation of nutrient mass balances.