Simulating the Effect of Climate Extremes on Groundwater Flow Through a Lakebed

Groundwater exchanges with lakes resulting from cyclical wet and dry climate extremes maintain lake levels in the environment in ways that are not well understood, in part because they remain difficult to simulate. To better understand the atypical groundwater interactions with lakes caused by climatic extremes, an original conceptual approach is introduced using MODFLOW‐2005 and a kinematic‐wave approximation to variably saturated flow that allows lake size and position in the basin to change while accurately representing the daily lake volume and three‐dimensional variably saturated groundwater flow responses in the basin. Daily groundwater interactions are simulated for a calibrated lake basin in Florida over a decade that included historic wet and dry departures from the average rainfall. The divergent climate extremes subjected nearly 70% of the maximum lakebed area and 75% of the maximum shoreline perimeter to both groundwater inflow and lake leakage. About half of the lakebed area subject to flow reversals also went dry. A flow‐through pattern present for 73% of the decade caused net leakage from the lake 80% of the time. Runoff from the saturated lake margin offset the groundwater deficit only about half of that time. A centripetal flow pattern present for 6% of the decade was important for maintaining the lake stage and generated 30% of all net groundwater inflow. Pumping effects superimposed on dry climate extremes induced the least frequent but most cautionary flow pattern with leakage from over 90% of the actual lakebed area.     

A strategy for modeling ground water rebound in abandoned deep mine systems

Discharges of polluted water from abandoned mines are a major cause of degradation of water resources worldwide. Pollution arises after abandoned workings flood up to surface level, by the process termed ground water rebound. As flow in large, open mine voids is often turbulent, standard techniques for modeling ground water flow (which assume laminar flow) are inappropriate for predicting ground water rebound. More physically realistic models are therefore desirable, yet these are often expensive to apply to all but the smallest of systems. An overall strategy for ground water rebound modeling is proposed, with models of decreasing complexity applied as the temporal and spatial scales of the systems under analysis increase. For relatively modest systems (area < 200 km2), a physically based modeling approach has been developed, in which 3-D pipe networks (representing major mine roadways, etc.) are routed through a variably saturated, 3-D porous medium (representing the country rock). For systems extending more than 100 to 3000 km2, a semidistributed model (GRAM) has been developed, which conceptualizes extensively interconnected volumes of workings as ponds, which are connected to other ponds only at discrete overflow points, such as major inter-mine roadways, through which flow can be efficiently modeled using the Prandtl-Nikuradse pipe-flow formulation. At the very largest scales, simple water-balance calculations are probably as useful as any other approach, and a variety of proprietary codes may be used for the purpose.     

A multiscale investigation of ground water flow at Clear Lake, Iowa

Ground water flow was investigated at Clear Lake, a 1468-ha glacial lake in north-central Iowa, as part of a comprehensive water quality study. A multiscale approach, consisting of seepage meters (and a potentiomanometer), Darcy's law, and an analytic element (AE) model, was used to estimate ground water inflow to and outflow from the lake. Estimates from the three methods disagreed. Seepage meters recorded a median-specific discharge of 0.25 mum/s, which produced a lake inflow rate between 90,750 and 138,200 m3/d, but no detectable outflow. A wave-induced Bernoulli effect probably compromised both inflow and outflow measurements. Darcy's law was applied to 11 zones around the lake, producing inflow and outflow values of 10,500 and 5000 m3/d, respectively. The AE model, GFLOW, coupled with the parameter estimation model, UCODE, simulated ground water flow in a 700-km2 region using 31 hydraulic head and base flow measurements as calibration targets. The model produced ground water inflow and outflow rates of 14,300 and 9200 m3/d, respectively. Although not a substitute for field data, the model's ability to simulate ground water flow to the lake and the region, estimate uncertainty for model parameters, and calculate a lake stage and associated lake water balance makes it a powerful tool for water quality management and an attractive alternative to the traditional methods of ground water/lake investigation.     

Sensitivity and moment analyses of head in variably saturated regimes

A numerical approach for approximating statistical moments of hydraulic heads of variably saturated flows in multi-dimensional porous media is developed. The approximation relies on a first-order Taylor series expansion of a finite element flow model and an adjoint state numerical method for variably saturated flows to evaluate sensitivities. This approach can be employed to analyze uncertainties associated with predictions of head of steady-state or transient flows in variably saturated porous media, with any type of boundary and initial conditions. Limitations of stochastic analytical methods such as spectral/perturbation approaches and the time-consuming Monte Carlo simulation technique are thus alleviated. An example is given to demonstrate the utility of the approach and to investigate the temporal evolution of head variances in a variably saturated flow regime. Results show that the fluctuation of the water table can have significant impacts on the propagation of the head variance.     

VS2DRTI: Simulating Heat and Reactive Solute Transport in Variably Saturated Porous Media

Variably saturated groundwater flow, heat transport, and solute transport are important processes in environmental phenomena, such as the natural evolution of water chemistry of aquifers and streams, the storage of radioactive waste in a geologic repository, the contamination of water resources from acid‐rock drainage, and the geologic sequestration of carbon dioxide. Up to now, our ability to simulate these processes simultaneously with fully coupled reactive transport models has been limited to complex and often difficult‐to‐use models. To address the need for a simple and easy‐to‐use model, the VS2DRTI software package has been developed for simulating water flow, heat transport, and reactive solute transport through variably saturated porous media. The underlying numerical model, VS2DRT, was created by coupling the flow and transport capabilities of the VS2DT and VS2DH models with the equilibrium and kinetic reaction capabilities of PhreeqcRM. Flow capabilities include two‐dimensional, constant‐density, variably saturated flow; transport capabilities include both heat and multicomponent solute transport; and the reaction capabilities are a complete implementation of geochemical reactions of PHREEQC. The graphical user interface includes a preprocessor for building simulations and a postprocessor for visual display of simulation results. To demonstrate the simulation of multiple processes, the model is applied to a hypothetical example of injection of heated waste water to an aquifer with temperature‐dependent cation exchange. VS2DRTI is freely available public domain software.     

Laboratory and numerical investigation of flow and transport near a seepage-face boundary

Laboratory and numerical modeling investigations were completed to study the unconfined ground water flow and transport processes near a seepage-face boundary. The laboratory observations were made in a radial sand tank and included measurements of the height of the seepage face, flow velocity near the seepage face, travel time distribution of multiple tracer slugs, and streamlines. All the observations were reliably reproduced with a three-dimensional, axi-symmetric, variably saturated ground water flow model. Physical data presented in this work demonstrate and quantify the importance of three-dimensional transport patterns within a seepage-face zone. The results imply that vertically averaged flow models that employ Dupuit approximations might introduce error in the analysis of localized solute transport near a seepage-face boundary. The experimental dataset reported in this work will also be of interest for those who are attempting to validate a numerical algorithm for solving ground water and contaminant discharge patterns near a surface-water boundary.     

From rainfall to spring discharge: Coupling conduit flow,subsurface matrix flow and surface flow in karst systems using a discrete–continuum model

Physics-based distributed models for simulating flow in karst systems are generally based on the discrete–continuum approach in which the flow in the three-dimensional fractured limestone matrix continuum is coupled with the flow in discrete one-dimensional conduits. In this study we present a newly designed discrete–continuum model for simulating flow in karst systems. We use a flexible spatial discretization such that complicated conduit networks can be incorporated. Turbulent conduit flow and turbulent surface flow are described by the diffusion wave equation whereas laminar variably saturated flow in the matrix is described by the Richards equation. Transients between free-surface and pressurized conduit flow are handled by changing the capacity term of the conduit flow equation. This new approach has the advantage that the transients in mixed conduit flow regimes can be handled without the Preissmann slot approach. Conduit–matrix coupling is based on the Peaceman’s well-index such that simulated exchange fluxes across the conduit–matrix interface are less sensitive to the spatial discretization. Coupling with the surface flow domain is based on numerical techniques commonly used in surface–subsurface models and storm water drainage models. Robust algorithms are used to simulate the non-linear flow processes in a coupled fashion. The model is verified and illustrated with simulation examples.     

An assembly model for simulation of large-scale ground water flow and transport

When managing large-scale ground water contamination problems, it is often necessary to model flow and transport using finely discretized domains--for instance (1) to simulate flow and transport near a contamination source area or in the area where a remediation technology is being implemented; (2) to account for small-scale heterogeneities; (3) to represent ground water-surface water interactions; or (4) some combination of these scenarios. A model with a large domain and fine-grid resolution will need extensive computing resources. In this work, a domain decomposition-based assembly model implemented in a parallel computing environment is developed, which will allow efficient simulation of large-scale ground water flow and transport problems using domain-wide grid refinement. The method employs common ground water flow (MODFLOW) and transport (RT3D) simulators, enabling the solution of almost all commonly encountered ground water flow and transport problems. The basic approach partitions a large model domain into any number of subdomains. Parallel processors are used to solve the model equations within each subdomain. Schwarz iteration is applied to match the flow solution at the subdomain boundaries. For the transport model, an extended numerical array is implemented to permit the exchange of dispersive and advective flux information across subdomain boundaries. The model is verified using a conventional single-domain model. Model simulations demonstrate that the proposed model operated in a parallel computing environment can result in considerable savings in computer run times (between 50% and 80%) compared with conventional modeling approaches and may be used to simulate grid discretizations that were formerly intractable.     

Comparison of local grid refinement methods for MODFLOW

Many ground water modeling efforts use a finite-difference method to solve the ground water flow equation, and many of these models require a relatively fine-grid discretization to accurately represent the selected process in limited areas of interest. Use of a fine grid over the entire domain can be computationally prohibitive; using a variably spaced grid can lead to cells with a large aspect ratio and refinement in areas where detail is not needed. One solution is to use local-grid refinement (LGR) whereby the grid is only refined in the area of interest. This work reviews some LGR methods and identifies advantages and drawbacks in test cases using MODFLOW-2000. The first test case is two dimensional and heterogeneous; the second is three dimensional and includes interaction with a meandering river. Results include simulations using a uniform fine grid, a variably spaced grid, a traditional method of LGR without feedback, and a new shared node method with feedback. Discrepancies from the solution obtained with the uniform fine grid are investigated. For the models tested, the traditional one-way coupled approaches produced discrepancies in head up to 6.8% and discrepancies in cell-to-cell fluxes up to 7.1%, while the new method has head and cell-to-cell flux discrepancies of 0.089% and 0.14%, respectively. Additional results highlight the accuracy, flexibility, and CPU time trade-off of these methods and demonstrate how the new method can be successfully implemented to model surface water-ground water interactions.     

Modelling lake stage and water balance of Lake Tana,Ethiopia

The level of Lake Tana, Ethiopia, fluctuates annually and seasonally following the patterns of changes in precipitation. In this study, a mass balance approach is used to estimate the hydrological balance of the lake. Water influx from four major rivers, subsurface inflow from the floodplains, precipitation, outflow from the lake constituting river discharge and evapotranspiration from the lake are analysed on monthly and annual bases. Spatial interpolation of precipitation using rain gauge data was conducted using kriging. Outflow from the lake was identified as the evaporation from the lake's surface as well as discharge at the outlet where the Blue Nile commences. Groundwater inflow is estimated using MODular three‐dimensional finite‐difference ground‐water FLOW model software that showed an aligned flow pattern to the river channels. The groundwater outflow is considered negligible based on the secondary sources that confirmed the absence of lake water geochemical mixing outside of the basin. Evaporation is estimated using Penman's, Meyer's and Thornwaite's methods to compare the mass balance and energy balance approaches. Meteorological data, satellite images and temperature perturbation simulations from Global Historical Climate Network of National Oceanographic and Atmospheric Administration are employed for estimation of evaporation input parameters. The difference of the inflow and outflow was taken as storage in depth and compared with the measured water level fluctuations. The study has shown that the monthly and annually calculated lake level replicates the observed values with root mean square error value of 0·17 and 0·15 m, respectively. Copyright © 2009 John Wiley & Sons, Ltd.     

Hydrological processes of a closed catchment‐lake system in a semi‐arid Mediterranean environment

Data collected in 4 years of field observations were used in conjunction with continuous simulation models to study, at the small‐basin scale, the water balance of a closed catchment‐lake system in a semi‐arid Mediterranean environment. The open water evaporation was computed with the Penman equation, using the data set collected in the middle of the lake. The surface runoff was partly measured at the main tributary and partly simulated using a distributed, catchment, hydrological model, calibrated with the observed discharge. The simplified structure of the developed modelling mainly concerns soil moisture dynamics and bedrock hydraulics, whereas the flow components are physically based. The calibration produced high efficiency coefficients and showed that surface runoff is greatly affected by soil water percolation into fractured bedrock. The bedrock reduces the storm‐flow peaks and the interflow and has important multi‐year effects on the annual runoff coefficients. The net subsurface outflow from the lake was calculated as the residual of the lake water balance. It was almost constant in the dry seasons and increased in the wet seasons, because of the moistening of the unsaturated soil. During the years of observation, rainfall 30% higher than average caused abundant runoff and a continuous rise in the lake water levels. The analysis allows to predict that, in years with lower than the average rainfall, runoff will be drastically reduced and will not be able to compensate for negative balance between precipitation and lake evaporation. Such highly unsteady situations, with great fluctuations in lake levels, are typical of closed catchment‐lake systems in the semi‐arid Mediterranean environment. Copyright © 2012 John Wiley & Sons, Ltd.     

Visualization and modeling of the colonization dynamics of a bioluminescent bacterium in variably saturated,translucent quartz sand

An experimental and numerical investigation was conducted to study the colonization dynamics of a bioluminescent bacterium, Pseudomonas fluorescens HK44, during growth in a porous medium under steady, variably saturated flow conditions. Experiments were conducted in a thin-slab light transmission chamber filled with uniform, translucent quartz sand. Steady, variably saturated flow conditions were established using drip emitters mounted on the top of the chamber, with glucose applied through a central dripper located directly above an inoculated region of the chamber. Periodic pulses of salicylate and a dye tracer were applied to induce bioluminescence of the bacterium to monitor colony expansion and to track changes in the hydraulic and transport properties of the sand. Changes in the apparent water saturation of the sand were quantified by monitoring light transmission through the chamber with a CCD camera. The colonized region expanded laterally by about 15 cm, and upward against the flow by 7–8 cm during the 6-day experiment while apparent saturations in the colonized region decreased by 7–9% and the capillary fringe dropped by ∼5 cm. The observed data were reproduced approximately using a numerical model that accounted for the processes of water flow, solute and bacterial transport, cell growth and accumulation, glucose and oxygen consumption, and gas diffusion and exchange. The results of this study illustrate some of the complexities associated with coupled flow, reactive transport, and biological processes in variably saturated porous media, such as localized desaturation, capillary fringe lowering effects, and upstream movement of bacterial colonization, that may not readily observable using other experimental techniques.     

Modeling ground water flow in alluvial mountainous catchments on a watershed scale

In large mountainous catchments, shallow unconfined alluvial aquifers play an important role in conveying subsurface runoff to the foreland. Their relatively small extent poses a serious problem for ground water flow models on the river basin scale. River basin scale models describing the entire water cycle are necessary in integrated water resources management and to study the impact of global climate change on ground water resources. Integrated regional-scale models must use a coarse, fixed discretization to keep computational demands low and to facilitate model coupling. This can lead to discrepancies between model discretization and the geometrical properties of natural systems. Here, an approach to overcome this discrepancy is discussed using the example of the German-Austrian Upper Danube catchment, where a coarse ground water flow model was developed using MODFLOW. The method developed uses a modified concept from a hydrological catchment drainage analysis in order to adapt the aquifer geometry such that it respects the numerical requirements of the chosen discretization, that is, the width and the thickness of cells as well as gradients and connectivity of the catchment. In order to show the efficiency of the developed method, it was tested and compared to a finely discretized ground water model of the Ammer subcatchment. The results of the analysis prove the applicability of the new approach and contribute to the idea of using physically based ground water models in large catchments.     

Understanding the water balance paradox in the Athabasca River Basin,Canada

This study demonstrates the importance of the including and appropriately parameterizing peatlands and forestlands for basin‐scale integrated surface–subsurface models in the northern boreal forest, with particular emphasis on the Athabasca River Basin (ARB). With a long‐term water balance approach to the ARB, we investigate reasons why downstream mean annual stream flow rates are consistently higher than upstream, despite the subhumid water deficit conditions in the downstream regimes. A high‐resolution 3D variably saturated subsurface and surface water flow and evapotranspiration model of the ARB is constructed based on the bedrock and surficial geology and the spatial distribution of peatlands and their corresponding eco‐regions. Historical climate data were used to drive the model for calibration against 40‐year long‐term average surface flow and groundwater observations during the historic instrumental period. The simulation results demonstrate that at the basin‐scale, peatlands and forestlands can have a strong influence on the surface–subsurface hydrologic systems. In particular, peatlands in the midstream and downstream regimes of the ARB increase the water availability to the surface–subsurface water systems by reducing water loss through evapotranspiration. Based on the comparison of forestland evapotranspiration between observation and simulation, the overall spatial average evapotranspiration in downstream forestlands is larger than that in peatlands and thus the water contribution to the stream flow in downstream areas is relatively minor. Therefore, appropriate representation of peatlands and forestlands within the basin‐scale hydrologic model is critical to reproduce the water balance of the ARB.     

The mechanism of Lake Kinneret salinization as a linear reservoir

The salinity of Lake Kinneret, Israel, is significantly higher than the salinity of the water from surface streams that flow to the lake. The relatively high salinity is a result of the activity of saline springs located at the bottom of the lake.The purpose of this work is to establish a general model for the salinization mechanism of Lake Kinneret. The model is based on the main components of the annual water and solute balance. Changes in time of the solute mass of the lake were described as a differential equation of a linear reservoir on an annual time scale. The model assumes that under any long-term operation policy of the lake, the components of the annual solute and water balance stay nearly constant in time.The model was tested for both steady-state conditions, and during changes in time, against measured lake salinity over the years 1968-2000. It was found that the major changes of lake salinity throughout the years were described well, despite the variety of rainfall amounts. Predictions of the expected lake salinity changes were proposed for the cases of controlled increase or decrease of saline springs discharge to the lake; for the changes of water quantity allowed to flow into or pumped out of the lake; and for various initial salinities. Predictions agree well with previous predictions made by statistical models.     

Modeling Variably Saturated Multispecies Reactive Groundwater Solute Transport with MODFLOW‐UZF and RT3D

A numerical model was developed that is capable of simulating multispecies reactive solute transport in variably saturated porous media. This model consists of a modified version of the reactive transport model RT3D (Reactive Transport in 3 Dimensions) that is linked to the Unsaturated‐Zone Flow (UZF1) package and MODFLOW. Referred to as UZF‐RT3D, the model is tested against published analytical benchmarks as well as other published contaminant transport models, including HYDRUS‐1D, VS2DT, and SUTRA, and the coupled flow and transport modeling system of CATHY and TRAN3D. Comparisons in one‐dimensional, two‐dimensional, and three‐dimensional variably saturated systems are explored. While several test cases are included to verify the correct implementation of variably saturated transport in UZF‐RT3D, other cases are included to demonstrate the usefulness of the code in terms of model run‐time and handling the reaction kinetics of multiple interacting species in variably saturated subsurface systems. As UZF1 relies on a kinematic‐wave approximation for unsaturated flow that neglects the diffusive terms in Richards equation, UZF‐RT3D can be used for large‐scale aquifer systems for which the UZF1 formulation is reasonable, that is, capillary‐pressure gradients can be neglected and soil parameters can be treated as homogeneous. Decreased model run‐time and the ability to include site‐specific chemical species and chemical reactions make UZF‐RT3D an attractive model for efficient simulation of multispecies reactive transport in variably saturated large‐scale subsurface systems.     

A two‐dimensional model of hydraulic performance of stormwater infiltration systems

Stormwater infiltration systems are a popular method for urban stormwater control. They are often designed using an assumption of one‐dimensional saturated outflow, although this is not very accurate for many typical designs where two‐dimensional (2D) flows into unsaturated soils occur. Available 2D variably saturated flow models are not commonly used for design because of their complexity and difficulties with the required boundary conditions. A purpose‐built stormwater infiltration system model was thus developed for the simulation of 2D flow from a porous storage. The model combines a soil moisture–based model for unsaturated soils with a ponded storage model and uses a wetting front‐tracking approach for saturated flows. The model represents the main physical processes while minimizing input data requirements. The model was calibrated and validated using data from laboratory 2D stormwater infiltration trench experiments. Calibrations were undertaken using five different combinations of calibration data to examine calibration data requirements. It was found that storage water levels could be satisfactorily predicted using parameters calibrated with either data from laboratory soils tests or observed water level data, whereas the prediction of soil moistures was improved through the addition of observed soil moisture data to the calibration data set. Copyright © 2012 John Wiley & Sons, Ltd.     

Predicting Soil-Water and Soil-Air Transport Properties and Their Effects on Soil-Vapor Extraction Efficiency

Accurate prediction of water and air Iran sport parameters in variably saturated soil is necessary for modeling of soil-vapor extraction (SVE) at soil sites contaminated with volatile organic chemicals (VOCs). An expression for predicting saturated water permeability (k_l,s) in undisturbed soils from the soil total porosity and the field capacity soil-water content was developed by fitting a tortuous-tube fluid flow model to measured water permeability and gas diffusivity data. The new k_l,s expression gave accurate predictions when tested against independent k_l,s data. The k_l,s expression was implemented in the Campbell relative water permeability model to yield a predictive model for water permeability in variably saturated, undisturbed soil. The water permeability model, together with recently developed predictive equations for gas permeability and gas diffusivity, was used in a two-dimensional numerical SVE model that also included non-equilibrium mass transfer of VOC from a separate phase (nonaqueous phase liquid [NAPL]) to the air phase. SVE: calculations showed that gas permeability is likely the most important factor controlling VOC migration and vapor extraction efficiency. Water permeability and gas diffusivity effects became significant at water contents near and above field capacity. The NAPL-air mass transfer coefficient also had large impacts on simulated vapor extraction efficiency. The calculations suggest that realistic SVE models need to include predictive expressions for both conveciive, diffusive. and phase-partitioning processes in natural, undisturbed soils.     

Simulation of field‐scale water flow and bromide transport in a cracked clay soil

Single domain models may seriously underestimate leaching of nutrients and pesticides to groundwater in clay soils with shrinkage cracks. Various two‐domain models have been developed, either empirical or physically based, which take into account the effects of cracks on water flow and solute transport. This paper presents a model concept that uses the clay shrinkage characteristics to derive crack volume and crack depth under transient field conditions. The concept has been developed to simulate field average behaviour of a field with cracks, rather than flow and transport at a small plot. Water flow and solute transport are described with basic physics, which allow process and scenario analysis. The model concept is part of the more general agrohydrological model SWAP, and is applied to a field experiment on a cracked clay soil, at which water flow and bromide transport were measured during 572 days. A single domain model was not able to mimic the field‐average water flow and solute transport. Incorporation of the crack concept considerably improved the simulation of water content and bromide leaching to the groundwater. Still deviations existed between the measured and simulated bromide concentration profiles. The model did not reproduce the observed bromide retardation in the top layer and the high bromide dispersion resulting from water infiltration at various soil depths. A sensitivity analysis showed that the amounts of bromide leached were especially sensitive to the saturated hydraulic conductivity of the top layer, the solute transfer from the soil matrix to crack water flow and the mean residence time of rapid drainage. The shrinkage characteristic and the soil hydraulic properties of the clay matrix showed a low sensitivity. Copyright © 2000 John Wiley & Sons, Ltd.     

Application of a Discrete-Continuum Model to Karst Aquifers in North China

A generalized discrete-continuum model is developed to simulate ground water flow in the karst aquifers of North China. The model is a hybrid numerical flow model, which takes into account both quick conduit flow and diffusive fissure flow. The conduit flow is represented by a discrete network model, and the fissure flow is modeled by a continuum approach. The developed model strongly emphasizes the function of the conduits in the flow fields. They control the general drainage pattern, as demonstrated in the simulation of a complex karst aquifer in North China. The model reproduces reasonably well the flow field in response to an unanticipated discharge of ground water from the karst aquifer into an underground mine based on the aquifer parameters that are manually calibrated from a multiple-well pumping test. Sensitivity of the model to the aquifer parameters was evaluated in the context of the case study.