Stage‐dependent transient storage of phosphorus in alluvial floodplains

Models for contaminant transport in streams commonly idealize transient storage as a well mixed but immobile system. These transient storage models capture rapid (near‐stream) hyporheic storage and transport, but do not account for large‐scale, stage‐dependent interaction with the alluvial aquifer. The objective of this research was to document transient storage of phosphorus (P) in coarse gravel alluvium potentially influenced by large‐scale, stage‐dependent preferential flow pathways (PFPs). Long‐term monitoring was performed at floodplain sites adjacent to the Barren Fork Creek and Honey Creek in northeastern Oklahoma. Based on results from subsurface electrical resistivity mapping which was correlated to hydraulic conductivity data, observation wells were installed both in higher hydraulic conductivity and lower hydraulic conductivity subsoils. Water levels in the wells were monitored over time, and water samples were obtained from the observation wells and the stream to document P concentrations at multiple times during high flow events. Contour plots indicating direction of flow were developed using water table elevation data. Contour plots of total P concentrations showed the alluvial aquifer acting as a transient storage zone, with P‐laden stream water heterogeneously entering the aquifer during the passage of a storm pulse, and subsequently re‐entering the stream during baseflow conditions. Some groundwater in the alluvial floodplains had total P concentrations that mirrored the streams' total P concentrations. A detailed analysis of P forms indicated that particulate P (i.e. P attached to particulates greater than 0·45 µm) was a significant portion of the P transport. This research suggests the need for more controlled studies on stage‐dependent transient storage in alluvial systems. Copyright © 2011 John Wiley & Sons, Ltd.     

Surface water–groundwater exchange dynamics in a tidal freshwater zone

In coastal rivers, tides can propagate for tens to hundreds of kilometres inland beyond the saltwater line. Yet the influence of tides on river–aquifer connectivity and solute transport in tidal freshwater zones (TFZs) is largely unknown. We estimate that along the TFZ of White Clay Creek (Delaware, USA), 11% of river water exchanges through tidal bank storage zones. Additional hyporheic processes such as flow through bedforms likely contribute even more exchange. The turnover length associated with tidal bank storage is 150 km, on the order of turnover lengths for all hyporheic exchange processes in non‐tidal rivers of similar size. Based on measurements at a transect of piezometers located 17 km from the coast, tides exchange 0.36 m³ of water across the banks and 0.86 m³ across the bed per unit river length. Exchange fluxes range from ?1.66 to 2.26 m day^?1 across the bank and ?0.84 to 1.88 m day^?1 across the bed. During rising tide, river water infiltrates into the riparian aquifer, and the downstream transport rate in the channel is low. During falling tide, stored groundwater is released to the river, and the downstream transport rate in the channel increases. Tidal bank storage zones may remove nutrients or other contaminants from river water and attenuate nutrient loads to coasts. Alternating expansion and contraction of aerobic zones in the riparian aquifer likely influence contaminant removal along flow paths. A clear need exists to understand contaminant removal and other ecosystem services in TFZs and adopt best management practices to promote these ecosystem services. Copyright © 2015 John Wiley & Sons, Ltd.     

Temporal Hyporheic Zone Response to Water Table Fluctuations

Expansion and contraction of the hyporheic zone due to temporal hydrologic changes between stream and riparian aquifer influence the biogeochemical cycling capacity of streams. Theoretical studies have quantified the control of groundwater discharge on the depth of the hyporheic zone; however, observations of temporal groundwater controls are limited. In this study, we develop the concept of groundwater‐dominated differential hyporheic zone expansion to explain the temporal control of groundwater discharge on the hyporheic zone in a third‐order stream reach flowing through glacially derived terrain typical of the Great Lakes region. We define groundwater‐dominated differential expansion of the hyporheic zone as: differing rates and magnitudes of hyporheic zone expansion in response to seasonal vs. storm‐related water table fluctuation. Specific conductance and vertical hydraulic gradient measurements were used to map changes in the hyporheic zone during seasonal water table decline and storm events. Planar and riffle beds were monitored in order to distinguish the cause of increasing hyporheic zone depth. Planar bed seasonal expansion of the hyporheic zone was of a greater magnitude and longer in duration (weeks to months) than storm event expansion (hours to days). In contrast, the hyporheic zone beneath the riffle bed exhibited minimal expansion in response to seasonal groundwater decline compared to storm related expansion. Results indicated that fluctuation in the riparian water table controlled seasonal expansion of the hyporheic zone along the planar bed. This groundwater induced hyporheic zone expansion could increase the potential for biogeochemical cycling and natural attenuation.     

Reducing monitoring gaps at the aquifer–river interface by modelling groundwater–surface water exchange flow patterns

This study investigates spatial patterns and temporal dynamics of aquifer–river exchange flow at a reach of the River Leith, UK. Observations of sub‐channel vertical hydraulic gradients at the field site indicate the dominance of groundwater up‐welling into the river and the absence of groundwater recharge from surface water. However, observed hydraulic heads do not provide information on potential surface water infiltration into the top 0–15 cm of the streambed as these depths are not covered by the existing experimental infrastructure. In order to evaluate whether surface water infiltration is likely to occur outside the ‘window of detection’, i.e. the shallow streambed, a numerical groundwater model is used to simulate hydrological exchanges between the aquifer and the river. Transient simulations of the successfully validated model (Nash and Sutcliff efficiency of 0·91) suggest that surface water infiltration is marginal and that the possibility of significant volumes of surface water infiltrating into non‐monitored shallow streambed sediments can be excluded for the simulation period. Furthermore, the simulation results show that with increasing head differences between river and aquifer towards the end of the simulation period, the impact of streambed topography and hydraulic conductivity on spatial patterns of exchange flow rates decreases. A set of peak flow scenarios with altered groundwater‐surface water head gradients is simulated in order to quantify the potential for surface water infiltration during characteristic winter flow conditions following the observation period. The results indicate that, particularly at the beginning of peak flow conditions, head gradients are likely to cause substantial increase in surface water infiltration into the streambed. The study highlights the potential for the improvement of process understanding of hyporheic exchange flow patterns at the stream reach scale by simulating aquifer‐river exchange fluxes with a standard numerical groundwater model and a simple but robust model structure and parameterization. Copyright © 2011 John Wiley & Sons, Ltd.     

From Grain to Floodplain: Evaluating heterogeneity of floodplain hydrostatigraphy using sedimentology,geophysics, and remote sensing

Stratigraphy is a fundamental component of floodplain heterogeneity and hydraulic conductivity and connectivity of alluvial aquifers, which affect hydrologic processes such as groundwater flow and hyporheic exchange. Watershed-scale hydrological models commonly simplify the sedimentology and stratigraphy of floodplains, neglecting natural floodplain heterogeneity and anisotropy. This study, conducted in the upper reach of the East River in the East River Basin, Colorado, USA, combines point-, meander-, and floodplain-scale data to determine key features of alluvial aquifers important for estimating hydrologic processes. We compare stratigraphy of two meanders with disparate geometries to explore floodplain heterogeneity and connectivity controls on flow and transport. Meander shape, orientation, and internal stratigraphy affected residence time estimates of laterally exchanged hyporheic water. Although the two meanders share a sediment source, vegetation, and climate, their divergent river migration histories resulted in contrasting meander hydrofacies. In turn, the extent and orientation of these elements controlled the effective hydraulic conductivity and, ultimately, estimates of groundwater transport and hyporheic residence times. Additionally, the meanders’ orientation relative to the valley gradient impacted the hydraulic gradient across the meanders—a key control of groundwater velocity. Lastly, we combine our field data with remotely sensed data and introduce a potential approach to estimate key hydrostratigraphic packages across floodplains. Prospective applications include contaminant transport studies, hyporheic models, and watershed models. © 2019 John Wiley & Sons, Ltd.     

The effects of floods on the temperature of riparian groundwater

The transition area between rivers and their adjacent riparian aquifers, which may comprise the hyporheic zone, hosts important biochemical reactions, which control water quality. The rates of these reactions and metabolic processes are temperature dependent. Yet the thermal dynamics of riparian aquifers, especially during flooding and dynamic groundwater flow conditions, has seldom been studied. Thus, we investigated heat transport in riparian aquifers during 3 flood events of different magnitudes at 2 sites along the same river. River and riparian aquifer temperature and water‐level data along the Lower Colorado River in Central Texas, USA, were monitored across 2‐dimensional vertical sections perpendicular to the bank. At the downstream site, preflood temperature penetration distance into the bank suggested that advective heat transport from lateral hyporheic exchange of river water into the riparian aquifer was occurring during relatively steady low‐flow river conditions. Although a small (20‐cm stage increase) dam‐controlled flood pulse had no observable influence on groundwater temperature, larger floods (40‐cm and >3‐m stage increases) caused lateral movement of distinct heat plumes away from the river during flood stage, which then retreated back towards the river after flood recession. These plumes result from advective heat transport caused by flood waters being forced into the riparian aquifer. These flood‐induced temperature responses were controlled by the size of the flood, river water temperature during the flood, and local factors at the study sites, such as topography and local ambient water table configuration. For the intermediate and large floods, the thermal disturbance in the riparian aquifer lasted days after flood waters receded. Large floods therefore have impacts on the temperature regime of riparian aquifers lasting long beyond the flood's timescale. These persistent thermal disturbances may have a significant impact on biochemical reaction rates, nutrient cycling, and ecological niches in the river corridor.     

Aerobic respiration in riparian exchange zones of regulated river corridors

River stage fluctuations drive surface water-groundwater exchanges within river corridors. This study evaluates how repeated daily stage fluctuations, representative of hydropeaking conditions, influence aerobic respiration of river-sourced dissolved organic carbon (DOC) in the riparian exchange zone using reactive flow and transport simulations. Over 50 hypothetical scenarios were modelled to evaluate how the duration of the daily flood signal, river DOC concentration, aquifer hydraulic conductivity and ambient groundwater flow condition affect the fate and transport of DOC and DO in the riparian aquifer. Time series subsurface snapshots highlight how the various factors influence the subsurface distribution of DOC and DO. The total mass of DOC respired per meter of river had a wide range depending on the parameters, spanning from 1.4 to 71 g over 24-h, with high hydraulic conductivity and losing ambient groundwater flow conditions favouring the largest amount of DOC respired. The ratio of DOC mass entering the riparian zone with the mass returning to the river showed that as little as 5% to as much as 76% of the DOC that enters the bank during stage fluctuations returns to the river. This return ratio is dependent on river DOC concentration, hydraulic conductivity and ambient groundwater flow. The results illustrate that stage variations due to river regulation can be a significant control on aerobic respiration in riparian exchange zones.     

Divergence and flow direction as indicators of subsurface heterogeneity and stage‐dependent storage in alluvial floodplains

Assuming homogeneity in alluvial aquifers is convenient, but limits our ability to accurately predict stream‐aquifer interactions. Research is needed on (i) identifying the presence of focused, as opposed to diffuse, groundwater discharge/recharge to streams and (ii) the magnitude and role of large‐scale bank and transient storage in alluvial floodplains relative to changes in stream stage. The objective of this research was to document and quantify the effect of stage‐dependent aquifer heterogeneity and bank storage relative to changes in stream stage using groundwater flow divergence and direction. Monitoring was performed in alluvial floodplains adjacent to the Barren Fork Creek and Honey Creek in northeastern Oklahoma. Based on results from subsurface electrical resistivity mapping, observation wells were installed in high and low electrical resistivity subsoils. Water levels in the wells were recorded real time using pressure transducers (August to October 2009). Divergence was used to quantify heterogeneity (i.e. variation in hydraulic conductivity, porosity, and/or aquifer thickness), and flow direction was used to assess the potential for large‐scale (100 m) bank or transient storage. Areas of localized heterogeneity appeared to act as divergence zones allowing stream water to quickly enter the groundwater system, or as flow convergence zones draining a large groundwater area. Maximum divergence or convergence occurred with maximum rates of change in flow rates or stream stage. Flow directions in the groundwater changed considerably between base and high flows, suggesting that the floodplains acted as large‐scale bank storage zones, rapidly storing and releasing water during passage of a storm hydrograph. During storm events at both sites, the average groundwater direction changed by at least 90° from the average groundwater direction during baseflow. Aquifer heterogeneity in floodplains yields hyporheic flows that are more responsive and spatially and temporally complex than would be expected compared to more common assumptions of homogeneity. Copyright © 2012 John Wiley & Sons, Ltd.     

Effect of rapidly changing river stage on uranium flux through the hyporheic zone

Measurement of ground water/surface water interaction within the hyporheic zone is increasingly recognized as an important aspect of subsurface contaminant fate and transport. Understanding the interaction between ground water and surface water is critical in developing a complete conceptual model of contaminant transport through the hyporheic zone. At the Hanford Site near Richland, Washington, ground water contaminated with uranium discharges to the Columbia River through the hyporheic zone. Ground water flux varies according to changes in hydraulic gradient caused by fluctuating river stage, which changes in response to operation of dams on the Columbia River. Piezometers and continuous water quality monitoring probes were installed in the hyporheic zone to provide long-term, high-frequency measurement of hydraulic gradient and estimated uranium concentrations. Subsequently, the flux of water and uranium was calculated for each half-hour time period over a 15-month study period. In addition, measurement of water levels in the near-shore unconfined aquifer enhanced the understanding of the relationship between river stage, aquifer elevation, and uranium flux. Changing river stage resulted in fluctuating hydraulic gradient within the hyporheic zone. Further, influx of river water caused lower uranium concentrations as a result of dilution. The methods employed in this study provide a better understanding of the interaction between surface and ground water in a situation with a dynamically varying vertical hydraulic gradient and illustrate how the combination of relatively standard methods can be used to derive an accurate estimation of water and contaminant flux through the hyporheic zone.     

Analytic solutions to transient groundwater flow under time-dependent sources in a heterogeneous aquifer bounded by fluctuating river stage

Analytical solutions for the water table and lateral discharge in a heterogeneous unconfined aquifer with time-dependent source and fluctuating river stage were derived and compared with those in an equivalent homogeneous aquifer. The heterogeneous aquifer considered consists of a number of sections of different hydraulic conductivity values. The source term and river stage were assumed to be time-dependent but spatially uniform. The solutions derived is useful in studying various groundwater flow problems in a horizontally heterogeneous aquifer since the spatially piecewise-constant hydraulic conductivity and temporally piecewise-constant recharge and lateral discharge can be used to quantify variations in these processes commonly observed in reality. Applying the solutions derived to an aquifer of three sections of different hydraulic conductivity values shown that (1) the aquifer heterogeneity significantly increases the spatial variation of the water table and thus its gradient but it has little effect on lateral discharge in the case of temporally and spatially uniform recharge, (2) the time-dependent but spatially uniform recharge increases the temporal variation of groundwater table over the entire aquifer but its effect on lateral discharge is limited in the zone near the river, and (3) the effect of river stage fluctuation on the water table and lateral discharge is limited in the zone near the river and the effect of the heterogeneity is to increase lateral discharge to or recharge from the river.     

Variation in Saturated Hydraulic Conductivity at the Outcrop Scale,the Whanganui Basin,New Zealand

Groundwater flow and contaminant transport are strongly influenced by hydrogeological spatial variation. Understanding the textural heterogeneity of aquifer and aquitard units is critical for predicting preferential flow pathways, but is often hindered by sparse hydrogeological data, widely spaced data points, and complex stratigraphy. Here, we demonstrate the application of a relatively new air permeameter technology, providing a cost-effective, rapid alternative for characterizing hydrostratigraphic units in the field. The aim of this research is to (1) characterize the variation of saturated hydraulic conductivity across shallow-marine hydrostratigraphic units of the Whanganui Basin, New Zealand, and (2) assess the variation of saturated hydraulic conductivity within individual hydrostratigraphic units and relate these changes to facies and depositional environments. Results suggest heterogeneity within fine-grained aquitard units is controlled by bioturbation, whereby burrowing, ingestion and defecation results in grain size segregation and differential micrite cementation. Coarse-grained heterolithic aquifer facies display sharp changes in permeability across planar to cross-bedded sets, related to current and wave energy fluctuations within shallow-marine depositional settings. Bedding plane orientation creates high permeability zones that promotes down dip subsurface flow. Down dip gradation of coarse-grained nearshore facies into fine-grained shelf facies along the paleo shoreline-shelf transect is suggested to promote lateral and vertical groundwater flow within the basin fill. Air permeameter techniques have potential for application within groundwater basins around the world, providing datasets that facilitate greater understanding of groundwater systems, informing practices and policies for targeted water quality management.     

A geomorphology‐based integrated stream–aquifer interaction model for semi‐gauged catchments

Many of the existing stream–aquifer interaction models available in the literature are very complex with limited applicability in semi‐gauged and ungauged catchments. In this study, to estimate the influent and effluent subsurface water fluxes under limited geo‐hydrometeorological data availability conditions, a simple stream–aquifer interaction model, namely, the variable parameter McCarthy–Muskingum (VPMM) hillslope‐storage Boussinesq (hsB) model, has been developed. This novel model couples the VPMM streamflow transport with the hsB groundwater flow transport modules in online mode. In this integrated model, the surface water–groundwater flux exchange process is modelled by the Darcian approach with the variable hydraulic heads between the river stage and groundwater table accounting for the rainfall forcing. Considering the exchange fluxes in the hyporheic zone and lateral overland flow contribution, this approach is field tested in a typical 48‐km stretch of the Brahmani River in eastern India to simulate the streamflow and its depth with the minimum Nash–Sutcliffe efficiency of 94% and 88%; the maximum root mean square error of 134 m³/s and 0.35 m; and the minimum index of agreement of 98% and 97%, respectively. This modelling approach could be very well utilized in data‐scarce world‐river basins to estimate the stream–aquifer exchange flux due to rainfall forcings.     

Coastal aquifer response to extreme storm events in Emilia‐Romagna,Italy

With global warming and sea level rise, many coastal systems will experience increased levels of inundation and storm flooding, especially along sandy lowland coastal areas, such as the Northern Adriatic coast (Italy). Understanding how extreme events may directly affect groundwater hydrology in shallow unconfined coastal aquifers is important to assess coastal vulnerability and quantify freshwater resources. This study investigates shallow coastal aquifer response to storm events. The transitory and permanent effects of storm waves are evaluated through the real time monitoring of groundwater and soil parameters, in order to characterize both the saturated and unsaturated portions of the coastal aquifer of Ravenna and Ferrara (southern Po Delta, Italy). Results highlight a general increase in hydraulic head and soil moisture, along with a decrease in groundwater salinity and pore water salinity due to rainfall infiltration during the 2 days storm event. The only exceptions are represented by the observation wells in proximity to the coastline (within 100 m), which recorded a temporary increase in soil and water salinity caused by the exceptional high waves, which persist on top of the dune crest during the storm event. This generates a saline plume that infiltrates through the vadose zone down to the saturated portion of the aquifer causing a temporary disappearance of the freshwater lens generally present, although limited in size, below the coastal dunes. Despite the high hydraulic conductivity, the aquifer system does not quickly recover the pre‐storm equilibrium and the storm effects are evident in groundwater and soil parameters after 10 days past the storm overwash recess.     

Geomorphology,hydrology, and aquatic vegetation drive seasonal hyporheic flow patterns across a gravel‐dominated floodplain

Across 1·7 km² of the Umatilla River floodplain (Oregon, USA), we investigated the influences of an ephemeral tributary and perennial ‘spring channel’ (fed only by upwelling groundwater) on hyporheic hydrology. We derived maps of winter and summer water‐table elevations from data collected at 46 monitoring wells and 19 stage gauges and used resulting maps to infer groundwater flow direction. Groundwater flow direction varied seasonally across the floodplain and was influenced by main channel stage, flooding, the tributary creek, and the location and direction of hyporheic exchange in the spring channel. Hyporheic exchange in the spring channel was evaluated with a geochemical mixing model, which confirmed patterns of floodplain groundwater movement inferred from water‐table maps and showed that the spring channel was fed predominantly by hyporheic water from the floodplain aquifer (87% during winter, 80% during summer), with its remaining flow supplied by upslope groundwater from the adjacent catchment aquifer. Summertime growth of aquatic macrophytes in the spring channel also influenced patterns of hyporheic exchange and groundwater flow direction in the alluvial aquifer by increasing flow resistance in the spring channel, locally raising surface water stage and adjacent water‐table elevation, and thereby altering the slope of the water‐table in the hyporheic zone. The Umatilla River floodplain is larger than most sites where hyporheic hydrology has been investigated in detail. Yet, our results corroborate other research that has identified off‐channel geomorphic features as important drivers of hyporheic hydrology, including previously published modeling efforts from a similar river and field observations from smaller streams. Copyright © 2007 John Wiley & Sons, Ltd.     

Impacts of river bed gas on the hydraulic and thermal dynamics of the hyporheic zone

Despite the presence of gas in river beds being a well known phenomenon, its potential feedbacks on the hydraulic and thermal dynamics of the hyporheic zone has not been widely studied. This paper explores hypotheses that the presence of accumulated gas impacts the hydraulic and thermal dynamics of a river bed due to changes in specific storage, hydraulic conductivity, effective porosity, and thermal diffusivity. The hypotheses are tested using data analysis and modelling for a study site on the urban River Tame, Birmingham, UK. Gas, predominantly attributed to microbial denitrification, was observed in the river bed up to around 14% by volume, and to at least 0.8 m depth below river bed. Numerical modelling indicates that, by altering the relative hydraulic conductivity distribution, the gas in the river bed leads to an increase of groundwater discharge from the river banks (relative to river bed) by a factor of approximately 2 during river low flow periods. The increased compressible storage of the gas phase in the river bed leads to an increase in the simulated volume of river water invading the river bed within the centre of the channel during storm events. The exchange volume can be more than 30% greater in comparison to that for water saturated conditions. Furthermore, the presence of gas also reduces the water-filled porosity, and so the possible depth of such invading flows may also increase markedly, by more than a factor of 2 in the observed case. Observed diurnal temperature variations within the gaseous river bed at 0.1 and 0.5 m depth are, respectively, around 1.5 and 6 times larger than those predicted for saturated sediments. Annual temperature fluctuations are seen to be enhanced by around 4 to 20% compared to literature values for saturated sediments. The presence of gas may thus alter the bulk thermal properties to such a degree that the use of heat tracer techniques becomes subject to a much greater degree of uncertainty. Although the likely magnitude of thermal and hydraulic changes due to the presence of gas for this site have been demonstrated, further research is needed into the origins of the gas and its spatial and temporal variability to enable quantification of the significance of these changes for chemical attenuation and hyporheic zone biology.     

Flow‐Through Stream Modeling with MODFLOW and MT3D: Certainties and Limitations

This paper aims to assess MODFLOW and MT3D capabilities for simulating the spread of contaminants from a river exhibiting an unusual relationship with an alluvial aquifer, with the groundwater head higher than the river head on one side and lower on the other (flow‐through stream). A series of simulation tests is conducted using a simple hypothetical model so as to characterize and quantify these limitations. Simulation results show that the expected contaminant spread could be achieved with a specific configuration composed of two sets of parameters: (1) modeled object parameters (hydraulic groundwater gradient, hydraulic conductivity values of aquifer and streambed), and (2) modeling parameters (vertical discretization of aquifer, horizontal refinement of stream modeled with River [RIV] package). The influence of these various parameters on simulation results is investigated, and potential complications and errors are identified. Contaminant spread from stream to aquifer is not always reproduced by MT3D due to the RIV package's inability to simulate lateral exchange fluxes between stream and aquifer. This paper identifies the need for a MODFLOW streamflow package allowing lateral stream‐aquifer interactions and streamflow routine calculations. Such developments could be of particular interest for modeling contaminated flow‐through streams.     

Dynamic river–groundwater exchange in the presence of a saline,semi‐confined aquifer

Understanding groundwater–surface water exchange in river banks is crucial for effective water management and a range of scientific disciplines. While there has been much research on bank storage, many studies assume idealized aquifer systems. This paper presents a field‐based study of the Tambo Catchment (southeast Australia) where the Tambo River interacts with both an unconfined aquifer containing relatively young and fresh groundwater (<500 μS/cm and <100 years old) and a semi‐confined artesian aquifer containing old and saline groundwater (electrical conductivity > 2500 μS/cm and >10 000 years old). Continuous groundwater elevation and electrical conductivity monitoring within the different aquifers and the river suggest that the degree of mixing between the two aquifers and the river varies significantly in response to changing hydrological conditions. Numerical modelling using MODFLOW and the solute transport package MT3DMS indicates that saline water in the river bank moves away from the river during flooding as hydraulic gradients reverse. This water then returns during flood recession as baseflow hydraulic gradients are re‐established. Modelling also indicates that the concentration of a simulated conservative groundwater solute can increase for up to ~34 days at distances of 20 and 40 m from the river in response to flood events approximately 10 m in height. For the same flood event, simulated solute concentrations within 10 m of the river increase for only ~15 days as the infiltrating low‐salinity river water drives groundwater dilution. Average groundwater fluxes to the river stretch estimated using Darcy's law were 7 m³/m/day compared with 26 and 3 m³/m/day for the same periods via mass balance using Radon (²²²Rn) and chloride (Cl), respectively. The study shows that by coupling numerical modelling with continuous groundwater–surface water monitoring, the transient nature of bank storage can be evaluated, leading to a better understanding of the hydrological system and better interpretation of hydrochemical data. Copyright © 2015 John Wiley & Sons, Ltd.     

A coupled vertically integrated model to describe lateral exchanges between surface and subsurface in large alluvial floodplains with a fully penetrating river

This paper presents a vertically averaged model for studying water and solute exchanges between a large river and its adjacent alluvial aquifer. The hydraulic model couples horizontal 2D Saint Venant equations for river flow and a 2D Dupuit equation for aquifer flow. The dynamic coupling between river and aquifer is provided by continuity of fluxes and water level elevation between the two domains. Equations are solved simultaneously by linking the two hydrological system matrices in a single global matrix in order to ensure the continuity conditions between river and aquifer and to accurately model two‐way coupling between these two domains. The model is applied to a large reach (about 36 km²) of the Garonne River (south‐western France) and its floodplain, including an instrumented site in a meander. Simulated hydraulic heads are compared with experimental measurements on the Garonne River and aquifer in the floodplain. Model verification includes comparisons for one point sampling date (27 piezometers, 30 March 2000) and for hydraulic heads variations measured continuously over 5 months (5 piezometers, 1 January to 1 June 2000). The model accurately reproduces the strong hydraulic connections between the Garonne River and its aquifer, which are confirmed by the simultaneous variation of the water level in the river and in piezometers located near the river bank. The simulations also confirmed that the model is able to reproduce groundwater flow dynamics during flood events. Given these results, the hydraulic model was coupled with a solute‐transport component, based on advection‐dispersion equations, to investigate the theoretical dynamics of a conservative tracer over 5 years throughout the 36 km² reach studied. Meanders were shown to favour exchanges between river and aquifer, and although the tracer was diluted in the river, the contamination moved downstream from the injection plots and affected both river banks. Copyright © 2008 John Wiley & Sons, Ltd.     

Dissolution in a variably confined carbonate platform: effects of allogenic runoff,hydraulic damming of groundwater inputs,and surface–groundwater exchange at the basin scale

In variably confined carbonate platforms, impermeable confining units collect rainfall over large areas and deliver runoff to rivers or conduits in unconfined portions of platforms. Runoff can increase river stage or conduit heads in unconfined portions of platforms faster than local infiltration of rainfall can increase groundwater heads, causing hydraulic gradients between rivers, conduits and the aquifer to reverse. Gradient reversals cause flood waters to flow from rivers and conduits into the aquifer where they can dissolve limestone. Previous work on impacts of gradient reversals on dissolution has primarily emphasized individual caves and little research has been conducted at basin scales. To address this gap in knowledge, we used legacy data to assess how a gradient of aquifer confinement across the Suwannee River Basin, north‐central Florida affected locations, magnitudes and processes of dissolution during 2005–2007, a period with extreme ranges of discharge. During intense rain events, runoff from the confining unit increased river stage above groundwater heads in unconfined portions of the platform, hydraulically damming inputs of groundwater along a 200 km reach of river. Hydraulic damming allowed allogenic runoff with SI_CAL < ?4 to fill the entire river channel and flow into the aquifer via reversing springs. Storage of runoff in the aquifer decreased peak river discharges downstream and contributed to dissolution within the aquifer. Temporary storage of allogenic runoff in karst aquifers represents hyporheic exchange at a scale that is larger than found in streams flowing over non‐karst aquifers because conduits in karst aquifers extend the area available for exchange beyond river beds deep into aquifers. Post‐depositional porosity in variably confined carbonate platforms should thus be enhanced along rivers that originate on confining units. This distribution should be considered in models of porosity distribution used to manage water and hydrocarbon resources in carbonate rocks. Copyright © 2013 John Wiley & Sons, Ltd.