Simulations of freshwater lens recharge and salt/freshwater interfaces using the HYDRUS and SWI2 packages for MODFLOW

Abstract The paper presents an evaluation of the combined use of the HYDRUS and SWI2 packages for MODFLOW as a potential tool for modeling recharge in coastal aquifers subject to saltwater intrusion. The HYDRUS package for MODFLOW solves numerically the one-dimensional form of the Richards equation describing water flow in variablysaturated media. The code computes groundwater recharge to or capillary rise from the groundwater table while considering weather, vegetation, and soil hydraulic property data. The SWI2 package represents in a simplified way variable-density flow associated with saltwater intrusion in coastal aquifers. Combining these two packages within the MODFLOW framework provides a more accurate description of vadose zone processes in subsurface systems with shallow aquifers, which strongly depend upon infiltration. The two packages were applied to a two-dimensional problem of recharge of a freshwater lens in a sandy peninsula, which is a typical geomorphologic form along the Baltic and the North Sea coasts, among other places. Results highlighted the sensitivity of calculated recharge rates to the temporal resolution of weather data. Using daily values of precipitation and potential evapotranspiration produced average recharge rates more than 20% larger than those obtained with weekly or monthly averaged weather data, leading to different trends in the evolution of freshwater-saltwater interfaces. Root water uptake significantly influenced both the recharge rate and the position of the freshwater-saltwater interface. The results were less sensitive to changes in soil hydraulic parameters, which in our study were found to affect average yearly recharge rates by up to 13%.

[1]  Jeffrey G. Arnold,et al.  Soil and Water Assessment Tool Theoretical Documentation Version 2009 , 2011 .

[2]  R. Winston.  ModelMuse - A Graphical User Interface for MODFLOW-2005 and PHAST , 2009 .

[3]  Ling Li,et al.  Capillary effect on flow in the drainage layer of highway pavement , 2012 .

[4]  G. O. Essink,et al.  Effects of climate change on coastal groundwater systems: A modeling study in the Netherlands , 2010 .

[5]  S. P. Neuman,et al.  Finite element simulation of flow in saturated - unsaturated soils considering water uptake by plants , 1974 .

[6]  Miroslav Šejna,et al.  Recent Developments and Applications of the HYDRUS Computer Software Packages , 2016 .

[7]  G. O. Essink,et al.  Shallow rainwater lenses in deltaic areas with saline seepage , 2011 .

[8]  Hamed Ketabchi,et al.  Conceptualization of a fresh groundwater lens influenced by climate change: A modeling study of an arid-region island in the Persian Gulf, Iran , 2014 .

[9]  Philip D. Meyer,et al.  Uncertainty analyses of infiltration and subsurface flow and transport for SDMP sites , 1997 .

[10]  T. Prabhakar Clement,et al.  Impacts of the 2004 tsunami on groundwater resources in Sri Lanka , 2006 .

[11]  R. Healy.  Simulating Water, Solute, and Heat Transport in the Subsurface with the VS2DI Software Package , 2008 .

[12]  C. Voss,et al.  SUTRA, a model for saturated-unsaturated, variable-density groundwater flow with solute or energy transport—Documentation of generalized boundary conditions, a modified implementation of specified pressures and concentrations or temperatures, and the lake capability , 2019, Techniques and Methods.

[13]  M. V. Genuchten,et al.  Review and comparison of models for describing non-equilibrium and preferential flow and transport in the vadose zone , 2003 .

[14]  D. A. Barry,et al.  Seawater intrusion processes, investigation and management: Recent advances and future challenges , 2013 .

[15]  Bernward Hölting,et al.  Hydrogeologie : Einführung in die Allgemeine und Angewandte Hydrogeologie , 2013 .

[16]  G. Destouni,et al.  The influence of temporal hydrological randomness on seawater intrusion in coastal aquifers , 2006 .

[17]  G.,et al.  Consequences of varied soil hydraulic and meteorological complexity on unsaturated zone time lag estimates , 2017 .

[18]  Alyssa M. Dausman,et al.  Documentation of the Seawater Intrusion (SWI2) Package for MODFLOW , 2014 .

[19]  J. Join,et al.  Modelling the response of fresh groundwater to climate and vegetation changes in coral islands , 2014, Hydrogeology Journal.

[20]  Christian D. Langevin,et al.  A comparison between SWI and SEAWAT - the importance of dispersion, inversion and vertical anisotropy , 2010 .

[21]  R. L. Johnson,et al.  User's guide to the Variably Saturated Flow (VSF) process to MODFLOW , 2006 .

[22]  Dinshaw N. Contractor,et al.  Recharge and aquifer response: Northern Guam Lens Aquifer, Guam, Mariana Islands , 2002 .

[23]  K. Sudheer,et al.  Implementation of Solute Transport in the Vadose Zone into the “HYDRUS Package for MODFLOW” , 2017, Ground water.

[24]  S. van der Zee,et al.  Cation exchange in a temporally fluctuating thin freshwater lens on top of saline groundwater , 2017, Hydrogeology Journal.

[25]  I. Neuweiler,et al.  Experimental and Numerical Analysis of Air Trapping in a Porous Medium with Coarse Textured Inclusions , 2016, Acta Geophysica.

[26]  A. W. Harbaugh.  MODFLOW-2005 : the U.S. Geological Survey modular ground-water model--the ground-water flow process , 2005 .

[27]  B. Faybishenko,et al.  Effect of temporal averaging of meteorological data on predictions of groundwater recharge , 2018 .

[28]  T. Clement,et al.  Impacts of Climate Change and Urbanization on Groundwater Resources in a Barrier Island , 2016 .

[29]  Warren J. Busscher,et al.  Simulation of Field Water Use and Crop Yield , 1980 .

[30]  D. E. Prudic,et al.  Documentation of the Unsaturated-Zone Flow (UZF1) Package for modeling Unsaturated Flow Between the Land Surface and the Water Table with MODFLOW-2005 , 2006 .

[31]  P. Hsieh,et al.  VS2DI-A graphical software package for simulating fluid flow and solute or energy transport in variably saturated porous media , 2000 .

[32]  M. Persson,et al.  Effects of hysteresis, rainfall dynamics, and temporal resolution of rainfall input data in solute transport modelling in uncropped soil , 2016 .

[33]  B. Scanlon,et al.  Estimating Groundwater Recharge , 2010 .

[34]  Jirka Šimůnek,et al.  Evaluating Interactions between Groundwater and Vadose Zone Using the HYDRUS‐Based Flow Package for MODFLOW , 2008 .

[35]  T. Gates,et al.  Modeling Variably Saturated Multispecies Reactive Groundwater Solute Transport with MODFLOW‐UZF and RT3D , 2013, Ground water.

[36]  G. Destouni,et al.  Solute transport through a heterogeneous coupled vadose‐saturated zone system with temporally random rainfall , 2001 .

[37]  Gerald N. Flerchinger,et al.  Intercode comparisons for simulating water balance of surficial sediments in semiarid regions , 2002 .

[38]  T. Günther,et al.  Numerical modelling of climate change impacts on freshwater lenses on the North Sea Island of Borkum using hydrological and geophysical methods , 2012 .

[39]  A. Verruijt.  A NOTE ON THE GHYBEN-HERZBERG FORMULA , 1968 .

[40]  N. Greggio,et al.  Seasonal dynamic of a shallow freshwater lens due to irrigation in the coastal plain of Ravenna, Italy , 2014, Hydrogeology Journal.

[41]  S. Luoma,et al.  Impacts of Future Climate Change and Baltic Sea Level Rise on Groundwater Recharge, Groundwater Levels, and Surface Leakage in the Hanko Aquifer in Southern Finland , 2014 .

[42]  Groundwater Recharge Modeling of the Nete Catchment ( Belgium ) Using the HYDRUS-1 D – MODFLOW Package , 2013 .

[43]  Dinshaw N. Contractor,et al.  Simulated effect of vadose infiltration on water levels in the Northern Guam Lens Aquifer. , 2000 .

[44]  R. Carsel,et al.  Developing joint probability distributions of soil water retention characteristics , 1988 .

[45]  D. Allen,et al.  From days to decades: numerical modelling of freshwater lens response to climate change stressors on small low-lying islands , 2014 .

[46]  M. Antonellini,et al.  Seasonal variation in natural recharge of coastal aquifers , 2013, Hydrogeology Journal.

[47]  Miroslav Šejna,et al.  Development and Applications of the HYDRUS and STANMOD Software Packages and Related Codes , 2008 .

[48]  B. Hasholt.  Bernward Hölting: HYDROGEOLOGIE: Einführung in die Allgemeine und Angewandte Hydrogeologie. Stuttgart 1980, X, 340 s., ill., 20 cm, DM26,80. Ferd. Enke Verlag. , 1981 .

[49]  Mary P Anderson,et al.  Importance of Unsaturated Zone Flow for Simulating Recharge in a Humid Climate , 2008, Ground water.

[50]  A. Elshorbagy,et al.  Water availability and forest growth in coarse-textured soils , 2011, Canadian Journal of Soil Science.

[51]  P. D. Louw,et al.  Response to recharge variation of thin rainwater lenses and their mixing zone with underlying saline groundwater , 2012 .

[52]  Brian D. Smerdon,et al.  Influence of subhumid climate and water table depth on groundwater recharge in shallow outwash aquifers , 2008 .