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 km2) 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 km2 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.

[1]  Jeremy B. Jones,et al.  Surface-subsurface interactions in stream ecosystems. , 1996, Trends in ecology & evolution.

[2]  P. Vervier,et al.  Hydromorphological control of phosphorus in a large free‐flowing gravel bed river: the Garonne River (France) , 2001 .

[3]  F. Triska,et al.  RETENTION AND TRANSPORT OF NUTRIENTS IN A THIRD-ORDER STREAM IN NORTHWESTERN CALIFORNIA: HYPORHEIC PROCESSES' , 1989 .

[4]  M. Trimmer,et al.  Nutrient dynamics in relation to surface-subsurface hydrological exchange in a groundwater fed chalk stream , 2006 .

[5]  S. Findlay.  Importance of surface‐subsurface exchange in stream ecosystems: The hyporheic zone , 1995 .

[6]  E. Swain,et al.  Simulation of Integrated Surface-Water/Ground-Water Flow and Salinity for a Coastal Wetland and Adjacent Estuary , 2005 .

[7]  Luc Lambs,et al.  Interactions between groundwater and surface water at river banks and the confluence of rivers , 2004 .

[8]  A. Boudou,et al.  Genotoxic and stress inductive potential of cadmium in Xenopus laevis larvae. , 2006, Aquatic toxicology.

[9]  Angela M. Gurnell,et al.  Linking hydrology and ecology , 2000 .

[10]  Roger Alexander Falconer,et al.  Coupling surface and subsurface flows in a depth averaged flood wave model , 2007 .

[11]  R. Lancaster.  Fluvial evolution of the Garonne River, France: integrating field data with numerical simulations , 2005 .

[12]  A. Bronstert,et al.  The impact of groundwater–surface water interactions on the water balance of a mesoscale lowland river catchment in northeastern Germany , 2007 .

[13]  Klement Tockner,et al.  A landscape perspective of surface-subsurface hydrological exchanges in river corridors , 2002 .

[14]  J. Stanford,et al.  An Ecosystem Perspective of Alluvial Rivers: Connectivity and the Hyporheic Corridor , 1993, Journal of the North American Benthological Society.

[15]  Emily H. Stanley,et al.  THE FUNCTIONAL SIGNIFICANCE OF THE HYPORHEIC ZONE IN STREAMS AND RIVERS , 1998 .

[16]  L. Ridolfi,et al.  Bedform-induced hyporheic exchange with unsteady flows , 2007 .

[17]  J. Sánchez-Pérez,et al.  Differentiated free-living and sediment-attached bacterial community structure inside and outside denitrification hotspots in the river–groundwater interface , 2008, Hydrobiologia.

[18]  Doerthe Tetzlaff,et al.  Connectivity between landscapes and riverscapes—a unifying theme in integrating hydrology and ecology in catchment science? , 2007 .

[19]  V. Zlotnik,et al.  Impact of heterogeneity, bed forms, and stream curvature on subchannel hyporheic exchange , 2004 .

[20]  J. Nash,et al.  River flow forecasting through conceptual models part I — A discussion of principles☆ , 1970 .

[21]  M. Dobson,et al.  Role of interaction zones between surface and ground waters in DOC transport and processing: considerations for river restoration , 1993 .

[22]  S. Recous,et al.  STICS : a generic model for the simulation of crops and their water and nitrogen balances. I. Theory, and parameterization applied to wheat and corn , 1998 .

[23]  D. M. Ely,et al.  A method for evaluating the importance of system state observations to model predictions, with application to the Death Valley regional groundwater flow system , 2004 .

[24]  S. Sauvage,et al.  Assessment of the quantitative and qualitative buffer function of an alluvial wetland: hydrological modelling of a large floodplain (Garonne River, France) , 2003 .

[25]  M. Thoms.  Floodplain–river ecosystems: lateral connections and the implications of human interference , 2003 .

[26]  Sabine Sauvage,et al.  A numerical tool to integrate biophysical diversity of a large regulated river: hydrobiogeochemical bases. The case of the Garonne River (France) , 2003 .

[27]  M. Dole‐Olivier.  Surface water-groundwater exchanges in three dimensions on a backwater of the Rhône river , 1998 .

[28]  J. M. Sánchez-Péreza,et al.  Change in groundwater chemistry as a consequence of suppression of floods : the case of the Rhine floodplain , 2002 .

[29]  G. Pinay,et al.  Change in Groundwater Nitrate Concentration in a Large River Floodplain: Denitrification, Uptake, or Mixing? , 1998, Journal of the North American Benthological Society.

[30]  M. Sophocleous.  Interactions between groundwater and surface water: the state of the science , 2002 .

[31]  Mustafa M. Aral,et al.  Simultaneous Solution Of Coupled Surface Water/groundwater Flow Systems , 2003 .

[32]  David S. White,et al.  Perspectives on Defining and Delineating Hyporheic Zones , 1993, Journal of the North American Benthological Society.

[33]  J. Ward,et al.  The Four-Dimensional Nature of Lotic Ecosystems , 1989, Journal of the North American Benthological Society.

[34]  Donald I. Siegel,et al.  Modeling surface and ground water mixing in the hyporheic zone using MODFLOW and MT3D , 2006 .

[35]  Keith Loague,et al.  Further testing of the Integrated Hydrology Model (InHM): event‐based simulations for a small rangeland catchment located near Chickasha, Oklahoma , 2005 .

[36]  Mustafa M. Aral,et al.  River networks and groundwater flow: a simultaneous solution of a coupled system , 2005 .

[37]  P. E. O'connell,et al.  River flow forecasting through conceptual models part III - The Ray catchment at Grendon Underwood , 1970 .

[38]  K. Bencala,et al.  Hyporheic zone hydrological processes , 2000 .

[39]  J. Bear.  Dynamics of Fluids in Porous Media , 1975 .

[40]  A. Howard.  Simulation modeling and statistical classification of escarpment planforms , 1995 .

[41]  K. Bencala,et al.  Modeling surface-subsurface hydrological interactions , 2000 .

[42]  R. Maxwell,et al.  Integrated surface-groundwater flow modeling: A free-surface overland flow boundary condition in a parallel groundwater flow model , 2006 .

[43]  M. Baker,et al.  Hydrological variability, organic matter supply and denitrification in the Garonne River ecosystem , 2004 .

[44]  Michael E. Campana,et al.  Seasonal variation in surface‐subsurface water exchange and lateral hyporheic area of two stream‐aquifer systems , 1998 .

[45]  A. I. Delis,et al.  Numerical solution of the two-dimensional shallow water equations by the application of relaxation methods , 2005 .

[46]  K. Loague,et al.  Hydrologic‐Response simulations for the R‐5 catchment with a comprehensive physics‐based model , 2001 .

[47]  M. Brunke,et al.  The ecological significance of exchange processes between rivers and groundwater , 1997 .

[48]  D. Corenblit,et al.  Les ajustements morphologiques contemporains du lit mineur de la Garonne, France et leurs effets sur l'hydrosystème fluvial , 2000 .

[49]  Harihar Rajaram,et al.  Modeling hyporheic zone processes , 2003 .

[50]  Ben Chie Yen,et al.  Modeling of conjunctive two-dimensional surface-three-dimensional subsurface flows , 2002 .

[51]  S. Sauvage,et al.  Nitrogen dynamics in the shallow groundwater of a riparian wetland zone of the Garonne, SW France: nitrate inputs, bacterial densities, organic matter supply and denitrification measurements , 2003 .

[52]  William E. Dietrich,et al.  Modeling fluvial erosion on regional to continental scales , 1994 .

[53]  Timothy A. Davis,et al.  A column pre-ordering strategy for the unsymmetric-pattern multifrontal method , 2004, TOMS.

[54]  P. Huyakorn,et al.  A fully coupled physically-based spatially-distributed model for evaluating surface/subsurface flow , 2004 .

[55]  L. Lautz,et al.  Impact of debris dams on hyporheic interaction along a semi‐arid stream , 2006 .