Numerical modeling groundwater recharge and its implication in water cycles of two interdunal valleys in the Sand Hills of Nebraska

The topography and geomorphology of the sand dunes and interdunal valleys in the Nebraska Sand Hills play important roles in regional water cycle by influencing groundwater recharge and evapotranspiration (ET). In this study, groundwater recharge, associated with precipitation and ET as well as soil hydraulics, and its spatial variations owing to the topography of dunes and valleys are examined. A method is developed to describe the recharge as a function of the storage capacity of dunes of various heights. After the method is tested using observations from a network of wells in the Sand Hills, it is used in the MODFLOW model to simulate and describe recharge effects on groundwater table depth at two different dune-valley sites. Analysis of modeled groundwater budget shows that the groundwater table depth in the interdunal valleys is critically influenced by vertical groundwater flows from surrounding dunes. At the site of higher dunes there are steadier and larger vertical groundwater flows in the dunes from their previous storage of precipitation. These vertical flows change to be horizontal converging groundwater flows and create upwelling in the interdunal valleys, where larger ET loss at the surface further enhances groundwater upwelling. Such interdunal valley is the major concentration area of the surface water and groundwater flow in the Sand Hills. At the site of shallow dunes and a broad interdunal valley the supply of groundwater from the dunes is trivial and inadequate to support upwelling of groundwater in the valley. The groundwater flows downward in the valley, and the valley surface is dry. Weak ET loss at the surface has a smaller effect on the groundwater storage than the precipitation recharge, making such area a source for groundwater.     

Seasonal expansion and contraction of stream networks in shallow groundwater systems

Surface water and groundwater are normally closely connected in areas with shallow aquifer systems. Stream systems can thus be considered as the outcrops of associated groundwater flows in areas with a shallow groundwater table and a previous subsurface. This situation prevails in sandy lowland areas where almost all rainfall percolates into the subsurface so that the surplus over evapotranspiration becomes part of a groundwater drainage system before it reappears at the surface in a stream. The stream network, being the interface with the groundwater system, must have the capacity to release the seasonally dependent precipitation surplus through the continuum of ground and surface waters. A river network therefore consists of a hierarchical system of different order and incision depth, of which the discharge-contributing component contracts and expands with the seasonal fluctuation in recharge and water table depth.

Coupling the mathematical expressions for groundwater drainage and stream flow enables development of a conjunctive model which relates the properties of a seasonally contracting and expanding stream network and related groundwater level fluctuation to the seasonal rainfall character for given geological and geomorphological conditions. This model further allows for assessment of drainage network response to a changing environment.     

Ground water flow in a desert basin: challenges of simulating transport of dissolved chromium

A large chromium plume that evolved from chromium releases in a valley near the Mojave River was studied to understand the processes controlling fate and migration of chromium in ground water and used as a tracer to study the dynamics of a basin and range ground water system. The valley that was studied is naturally arid with high evapotranspiration such that essentially no precipitation infiltrates to the water table. The dominant natural hydrogeologic processes are recharge to the ground water system from the Mojave River during the infrequent episodes when there is flow in the river, and ground water flow toward a playa lake where the ground water evaporates. Agricultural pumping in the valley from the mid-1930s to the 1970s significantly altered ground water flow conditions by decreasing water levels in the valley by more than 20 m. This pumping declined significantly as a result of dewatering of the aquifer, and water levels have since recovered modestly. The ground water system was modeled using MODFLOW, and chromium transport was simulated using MT3D. Several innovative modifications were made to these modeling programs to simulate important processes in this ground water system. Modifications to MODFLOW include developing a new well package that estimates pumping rates from irrigation wells at each time step based on available drawdown. MT3D was modified to account for mass trapped above the water table when the water table declines beneath nonirrigated areas and to redistribute mass to the system when water levels rise.     

RIPGIS-NET: a GIS tool for riparian groundwater evapotranspiration in MODFLOW

RIPGIS-NET, an Environmental System Research Institute (ESRI's) ArcGIS 9.2/9.3 custom application, was developed to derive parameters and visualize results of spatially explicit riparian groundwater evapotranspiration (ETg), evapotranspiration from saturated zone, in groundwater flow models for ecohydrology, riparian ecosystem management, and stream restoration. Specifically RIPGIS-NET works with riparian evapotranspiration (RIP-ET), a modeling package that works with the MODFLOW groundwater flow model. RIP-ET improves ETg simulations by using a set of eco-physiologically based ETg curves for plant functional subgroups (PFSGs), and separates ground evaporation and plant transpiration processes from the water table. The RIPGIS-NET program was developed in Visual Basic 2005, .NET framework 2.0, and runs in ArcMap 9.2 and 9.3 applications. RIPGIS-NET, a pre- and post-processor for RIP-ET, incorporates spatial variability of riparian vegetation and land surface elevation into ETg estimation in MODFLOW groundwater models. RIPGIS-NET derives RIP-ET input parameters including PFSG evapotranspiration curve parameters, fractional coverage areas of each PFSG in a MODFLOW cell, and average surface elevation per riparian vegetation polygon using a digital elevation model. RIPGIS-NET also provides visualization tools for modelers to create head maps, depth to water table (DTWT) maps, and plot DTWT for a PFSG in a polygon in the Geographic Information System based on MODFLOW simulation results.     

Effect of elevation resolution on evapotranspiration simulations using MODFLOW

Surface elevations represented in MODFLOW head-dependent packages are usually derived from digital elevation models (DEMs) that are available at much high resolution. Conventional grid refinement techniques to simulate the model at DEM resolution increases computational time, input file size, and in many cases are not feasible for regional applications. This research aims at utilizing the increasingly available high resolution DEMs for effective simulation of evapotranspiration (ET) in MODFLOW as an alternative to grid refinement techniques. The source code of the evapotranspiration package is modified by considering for a fixed MODFLOW grid resolution and for different DEM resolutions, the effect of variability in elevation data on ET estimates. Piezometric head at each DEM cell location is corrected by considering the gradient along row and column directions. Applicability of the research is tested for the lower Rio Grande (LRG) Basin in southern New Mexico. The DEM at 10 m resolution is aggregated to resampled DEM grid resolutions which are integer multiples of MODFLOW grid resolution. Cumulative outflows and ET rates are compared at different coarse resolution grids. Results of the analysis conclude that variability in depth-to-groundwater within the MODFLOW cell is a major contributing parameter to ET outflows in shallow groundwater regions. DEM aggregation methods for the LRG Basin have resulted in decreased volumetric outflow due to the formation of a smoothing error, which lowered the position of water table to a level below the extinction depth.     

IHMS—Integrated Hydrological Modelling System. Part 1. Hydrological processes and general structure

A newly Integrated Hydrological Modelling System (IHMS) has been developed to study the impact of changes in climate, land use and water management on groundwater and seawater intrusion (SWI) into coastal areas. The system represents the combination of three models, which can, if required, be run separately. It has been designed to assess the combined impact of climate, land use and groundwater abstraction changes on river, drainage and groundwater flows, groundwater levels and, where appropriate, SWI. The approach is interdisciplinary and reflects an integrated water management approach. The system comprises three packages: the Distributed Catchment Scale Model (DiCaSM), MODFLOW (96 and 2000) and SWI models. In addition to estimating all water balance components, DiCaSM, produces the recharge data that are used as input to the groundwater flow model of the US Geological Survey, MODFLOW. The latter subsequently generates the head distribution and groundwater flows that are used as input to the SWI model, SWI. Thus, any changes in land use, rainfall, water management, abstraction, etc. at the surface are first handled by DiCaSM, then by MODFLOW and finally by the SWI. The three models operate at different spatial and temporal scales and a facility (interface utilities between models) to aggregate/disaggregate input/output data to meet a desired spatial and temporal scale was developed allowing smooth and easy communication between the three models. As MODFLOW and SWI are published and in the public domain, this article focuses on DiCaSM, the newly developed unsaturated zone DiCaSM and equally important the interfacing utilities between the three models. DiCaSM simulates a number of hydrological processes: rainfall interception, evapotranspiration, surface runoff, infiltration, soil water movement in the root zone, plant water uptake, crop growth, stream flow and groundwater recharge. Input requirements include distributed data sets of rainfall, land use, soil types and digital terrain; climate data input can be either distributed or non‐distributed. The model produces distributed and time series output of all water balance components including potential evapotranspiration, actual evapotranspiration, rainfall interception, infiltration, plant water uptake, transpiration, soil water content, soil moisture (SM) deficit, groundwater recharge rate, stream flow and surface runoff. This article focuses on details of the hydrological processes and the various equations used in DiCaSM, as well as the nature of the interface to the MODFLOW and SWI models. Furthermore, the results of preliminary tests of DiCaSM are reported; these include tests related to the ability of the model to predict the SM content of surface and subsurface soil layers, as well as groundwater levels. The latter demonstrates how the groundwater recharge calculated from DiCaSM can be used as input into the groundwater model MODFLOW using aggregation and disaggregation algorithms (built into the interface utility). SWI has also been run successfully with hypothetical examples and was able to reproduce the results of some of the original examples of Bakker and Schaars ( 2005 ). In the subsequent articles, the results of applications to different catchments will be reported. Copyright © 2010 John Wiley & Sons, Ltd.     

Recharge estimation for transient ground water modeling

Reliable ground water models require both an accurate physical representation of the system and appropriate boundary conditions. While physical attributes are generally considered static, boundary conditions, such as ground water recharge rates, can be highly variable in both space and time. A practical methodology incorporating the hydrologic model HELP3 in conjunction with a geographic information system was developed to generate a physically based and highly detailed recharge boundary condition for ground water modeling. The approach uses daily precipitation and temperature records in addition to land use/land cover and soils data. The importance of the method in transient ground water modeling is demonstrated by applying it to a MODFLOW modeling study in New Jersey. In addition to improved model calibration, the results from the study clearly indicate the importance of using a physically based and highly detailed recharge boundary condition in ground water quality modeling, where the detailed knowledge of the evolution of the ground water flowpaths is imperative. The simulated water table is within 0.5 m of the observed values using the method, while the water levels can differ by as much as 2 m using uniform recharge conditions. The results also show that the combination of temperature and precipitation plays an important role in the amount and timing of recharge in cooler climates. A sensitivity analysis further reveals that increasing the leaf area index, the evaporative zone depth, or the curve number in the model will result in decreased recharge rates over time, with the curve number having the greatest impact.     

Evaluation of Noble Gas Recharge Temperatures in a Shallow Unconfined Aquifer

Water table temperatures inferred from dissolved noble gas concentrations (noble gas temperatures, NGT) are useful as a quantitative proxy for air temperature change since the last glacial maximum. Despite their importance in paleoclimate research, few studies have investigated the relationship between NGT and actual recharge temperatures in field settings. This study presents dissolved noble gas data from a shallow unconfined aquifer heavily impacted by agriculture. Considering samples unaffected by degassing, NGT calculated from common physically based interpretive gas dissolution models that correct measured noble gas concentrations for "excess air" agreed with measured water table temperatures (WTT). The ability to fit data to multiple interpretive models indicates that model goodness-of-fit does not necessarily mean that the model reflects actual gas dissolution processes. Although NGT are useful in that they reflect WTT, caution is recommended when using these interpretive models. There was no measurable difference in excess air characteristics (amount and degree of fractionation) between two recharge regimes studied (higher flux recharge primarily during spring and summer vs. continuous, low flux recharge). Approximately 20% of samples had dissolved gas concentrations below equilibrium concentration with respect to atmospheric pressure, indicating degassing. Geochemical and dissolved gas data indicate that saturated zone denitrification caused degassing by gas stripping. Modeling indicates that minor degassing (<10% ΔNe) may cause underestimation of ground water recharge temperature by up to 2°C. Such errors are problematic because degassing may not be apparent and degassed samples may be fit by a model with a high degree of certainty.     

Estimating Persistent Mass Flux of Volatile Contaminants from the Vadose Zone to Ground Water

Contaminants may persist for long time periods within low permeability portions of the vadose zone where they cannot be effectively treated and are a potential continuing source of contamination to ground water. Setting appropriate vadose zone remediation goals typically requires evaluating these persistent sources in terms of their impact on meeting ground water remediation goals. Estimating the impact on ground water can be challenging at sites with low aqueous recharge rates where vapor-phase movement is the dominant transport process in the vadose zone. Existing one-dimensional approaches for simulating transport of volatile contaminants in the vadose zone are considered and compared to a new flux-continuity-based assessment of vapor-phase contaminant movement from the vadose zone to the ground water. The flux-continuity-based assessment demonstrates that the ability of the ground water to move contaminant away from the water table controls the vapor-phase mass flux from the vadose zone across the water table. Limitations of these approaches are then discussed with respect to the required assumptions and the need to incorporate three-dimensional processes when evaluating vapor-phase transport from the vadose zone to the ground water. The carbon tetrachloride plume at the U.S. Department of Energy Hanford Site is used as the example site where persistent vadose zone contamination needs to be considered in the context of ground water remediation.     

Importance of unsaturated zone flow for simulating recharge in a humid climate

Transient recharge to the water table is often not well understood or quantified. Two approaches for simulating transient recharge in a ground water flow model were investigated using the Trout Lake watershed in north-central Wisconsin: (1) a traditional approach of adding recharge directly to the water table and (2) routing the same volume of water through an unsaturated zone column to the water table. Areas with thin (less than 1 m) unsaturated zones showed little difference in timing of recharge between the two approaches; when water was routed through the unsaturated zone, however, less recharge was delivered to the water table and more discharge occurred to the surface because recharge direction and magnitude changed when the water table rose to the land surface. Areas with a thick (15 to 26 m) unsaturated zone were characterized by multimonth lags between infiltration and recharge, and, in some cases, wetting fronts from precipitation events during the fall overtook and mixed with infiltration from the previous spring snowmelt. Thus, in thicker unsaturated zones, the volume of water infiltrated was properly simulated using the traditional approach, but the timing was different from simulations that included unsaturated zone flow. Routing of rejected recharge and ground water discharge at land surface to surface water features also provided a better simulation of the observed flow regime in a stream at the basin outlet. These results demonstrate that consideration of flow through the unsaturated zone may be important when simulating transient ground water flow in humid climates with shallow water tables.     

Impact of Recharge Through Residual Oil Upon Sampling of Underlying Ground Water

At an aviation gasoline spill site in Traverse City, Michigan, historical records indicate a positive correlation between significant rainfall events and increased concentrations of slightly soluble organic compounds in the monitoring wells of the site. To investigate the recharge effect on ground water quality due to infiltrating, water percolating past residual oil and into the saturated zone, an in situ infiltration experiment was performed at the site. Sampling cones were set at various depths below a circular test area, 13 feet (4 meters) in diameter. Rainfall was simulated by sprinkling the test area at a rate sufficiently low to prevent runoff. The sampling cones for soil-gas and ground water quality were installed in the unsaturated and saturated zones to observe the effects of the recharge process. At the time of the test, the water table was below the residual oil layer. The responses of the soil-gas and ground water quality were monitored during the recharge and drainage periods, which resulted from the sprinkling.
Infiltrated water was determined to have transported organic constituents of the residual oil, specifically benzene, toluene, ethylbenzene, and ortho-xylene (BTEX), into the ground water beneath the water table, elevating the aqueous concentrations of these constituents in the saturated zone. Soil-gas concentrations of the organic compounds in the unsaturated zone increased with depth and time after the commencement of infiltration. Reaeration of the unconfined aquifer via the infiltrated water was observed. It is concluded that water quality measurements are directly coupled to recharge events for the sandy type of aquifer with an overlying oil phase, which was studied in this work. Ground water sampling strategies and data analysis need to reflect the effect of recharge from precipitation on shallow, unconfined aquifers where an oil phase may be present.     

Dynamics of flood water infiltration and ground water recharge in hyperarid desert

A study on flood water infiltration and ground water recharge of a shallow alluvial aquifer was conducted in the hyperarid section of the Kuiseb River, Namibia. The study site was selected to represent a typical desert ephemeral river. An instrumental setup allowed, for the first time, continuous monitoring of infiltration during a flood event through the channel bed and the entire vadose zone. The monitoring system included flexible time domain reflectometry probes that were designed to measure the temporal variation in vadose zone water content and instruments to concurrently measure the levels of flood and ground water. A sequence of five individual floods was monitored during the rainy season in early summer 2006. These newly generated data served to elucidate the dynamics of flood water infiltration. Each flood initiated an infiltration event which was expressed in wetting of the vadose zone followed by a measurable rise in the water table. The data enabled a direct calculation of the infiltration fluxes by various independent methods. The floods varied in their stages, peaks, and initial water contents. However, all floods produced very similar flux rates, suggesting that the recharge rates are less affected by the flood stages but rather controlled by flow duration and available aquifer storage under it. Large floods flood the stream channel terraces and promote the larger transmission losses. These, however, make only a negligible contribution to the recharge of the ground water. It is the flood duration within the active streambed, which may increase with flood magnitude that is important to the recharge process.     

Comparing MODFLOW simulation options for predicting intra‐meander flux

During the evolution of meander bends, the intra‐meander groundwater head gradients steepen and generate zones of accelerated water and nutrient intra‐meander fluxes important for ecosystem processes. This paper compares and contrasts three MODFLOW groundwater model packages based on their simulation of intra‐meander flux for two stages of meander evolution observed in a sandbox river table and one level of river bed clogging, where the hydraulic conductivity in the river bed is lower than in the adjacent aquifer. These packages are the Time‐Variant Specified Head package [constant head (CHD)], River package (RIV), and Streamflow‐Routing package (SFR2), each controlling the groundwater or river head bounding the intra‐meander region. The RIV and SFR2 packages fix river stage and allow for variation in groundwater head below the river, which is suggested for simulating intra‐meander flux for all sinuosities with and without river bed clogging whenever river bed parameters are available. The CHD package fixes below river groundwater head and fails to simulate intra‐meander head loss and flux in meanders with high sinuosity or river bed clogging. In low sinuosity meanders and in cases without river bed clogging, there were no significant differences between MODFLOW packages for simulating river intra‐meander head loss and flux. This research demonstrates why MODFLOW users need to consider the limitations of each package when simulating intra‐meander flux in reaches with river bed clogging, high sinuosity, or similarly steep hydraulic gradients. Copyright © 2014 John Wiley & Sons, Ltd.     

Salinization of a fresh palaeo-ground water resource by enhanced recharge

Deterioration of fresh ground water resources caused by salinization is a growing issue in many arid and semi-arid parts of the world. We discuss here the incipient salinization of a 10(4) km2 area of fresh ground water (<3,000 mg/L) in the semiarid Murray Basin of Australia caused by widespread changes in land use. Ground water 14C concentrations and unsaturated zone Cl soil water inventories indicate that the low salinity ground water originated mainly from palaeo-recharge during wet climatic periods more than 20,000 years ago. However, much of the soil water in the 20 to 60 m thick unsaturated zone throughout the area is generally saline (>15,000 mg/L) because of relatively high evapotranspiration during the predominantly semiarid climate of the last 20,000 years. Widespread clearing of native vegetation over the last 100 years and replacement with crops and pastures leads to enhancement of recharge rates that progressively displace the saline soil-water from the unsaturated zone into the ground water. To quantify the impact of this new hydrologic regime, a one-dimensional model that simulates projected ground water salinities as a function of depth to ground water, recharge rates, and soil water salt inventory was developed. Results from the model suggest that, in some areas, the ground water salinity within the top 10 m of the water table is likely to increase by a factor of 2 to 6 during the next 100 years. Ground water quality will therefore potentially degrade beyond the point of usefulness well before extraction of the ground water exhausts the resource.     

Extinction depth and evapotranspiration from ground water under selected land covers

In many landscapes, vegetation extracts water from both the unsaturated and the saturated zones. The partitioning of evapotranspiration (ET) into vadose zone evapotranspiration and ground water evapotranspiration (GWET) is complex because it depends on land cover and subsurface characteristics. Traditionally, the GWET fraction is assumed to decay with increasing depth to the water table (DTWT), attaining a value of 0 at what is termed the extinction depth. A simple assumption of linear decay with depth is often used but has never been rigorously examined using unsaturated-saturated flow simulations. Furthermore, it is not well understood how to relate extinction depths to characteristics of land cover and soil texture. In this work, variable saturation flow theory is used to simulate GWET for three land covers and a range of soil properties under drying soil conditions. For a water table within half a meter of the land surface, nearly all ET is extracted from ground water due to the close hydraulic connection between the unsaturated and the saturated zones. For deep-rooted vegetation, the decoupling of ground water and vadose zone was found to begin at water table depths between 30 and 100 cm, depending on the soil texture. The decline of ET with DTWT is better simulated by an exponential decay function than the commonly used linear decay. A comparison with field data is consistent with the findings of this study. Tables are provided to vary the extinction depth for heterogeneous landscapes with different vegetation cover and soil properties.     

Estimating groundwater recharge in lowland watersheds

Little attention has been given to the role of groundwater in the hydrological cycle of lowland watersheds. Our objective in this study was to estimate total recharge to groundwater by analysing water table response to storm events and the rate at which water was transferred into the shallow aquifer. This was conducted at three sites in a rural watershed in the lower Atlantic coastal plain near Charleston, South Carolina, USA. A novel version of the water table fluctuation method was used to estimate total recharge to the shallow aquifer by comparing hourly data of water table position following storm events and measuring water table recession behavior, rather than subjective graphical analysis methods. Also, shallow aquifer recharge rates (vertical fluxes) were estimated using Darcy's Law by comparing static water levels in a water table well and in a shallow piezometer during dry periods. The total annual recharge estimated ranged from 107 ± 39 mm·yr^–1 (5–10% of annual precipitation) at a poorly drained topographic low area to 1140 ± 230 mm·yr^–1 (62–94% of annual precipitation) for a moderately well‐drained upland site. The average aquifer recharge rate was 114 ± 60 mm·yr^–1, which is similar to previous estimations of base flow for the ephemeral third‐order streams in this watershed. The difference in the two methods may have been caused by processes not accounted for in the Darcy flux method, soil moisture deficits, and average evapotranspiration demand, which is about 1100 mm·yr^–1 for this region. Although other factors also can affect partitioning of recharge, an integrated approach to inspecting easily gathered groundwater data can provide information on an often neglected aspect of water budget estimation. We also discuss the effects of land use change on recharge reduction, given a typical development scenario for the region. Copyright © 2011 John Wiley & Sons, Ltd.     

A Simple Method for Simulating Groundwater Interactions with Fens to Forecast Development Effects

Protection of fens–wetlands dependent on groundwater discharge–requires characterization of groundwater sources and stresses. Because instrumentation and numerical modeling of fens is labor intensive, easy-to-apply methods that model fen distribution and their vulnerability to development are desirable. Here we demonstrate that fen areas can be simulated using existing steady-state MODFLOW models when the unsaturated zone flow (UZF) package is included. In cells where the water table is near land surface, the UZF package calculates a head difference and scaled conductance at these “seepage drain” cells to generate average rates of vertical seepage to the land. This formulation, which represents an alternative to blanketing the MODFLOW domain with drains, requires very little input from the user because unsaturated flow-routing is inactive and results are primarily driven by easily obtained topographic information. Like the drain approach, it has the advantage that the distribution of seepage areas is not predetermined by the modeler, but rather emerges from simulated heads. Beyond the drain approach, it takes account of intracell land surface variation to explicitly quantify multiple surficial flows corresponding to infiltration, rejected recharge, recharge and land-surface seepage. Application of the method to a basin in southeastern Wisconsin demonstrates how it can be used as a decision-support tool to first, reproduce fen distribution and, second, forecast drawdown and reduced seepage at fens in response to shallow pumping.     

Simulating Piecewise-Linear Surface Water and Ground Water Interactions with MODFLOW

The standard MODFLOW packages offer limited capabilities to model piecewise-linear boundary conditions to describe ground water–surface water interaction. Specifically, MODFLOW is incapable of representing a Cauchy-type boundary with different resistances for discharge or recharge conditions. Such a more sophisticated Cauchy boundary condition is needed to properly represent surface waters alternatively losing water through the bottom (high resistance) or gaining water mostly near the water surface (low resistance). One solution would be to create a new package for MODFLOW to accomplish this. However, it is also possible to combine multiple instances of standard packages in a single cell to the same effect. In this specific example, the general head boundary package is combined with the drain package to arrive at the desired piecewise-linear behavior. In doing so, the standard USGS MODFLOW version can be used without any modifications at the expense of a minor increase in preprocessing and postprocessing and computational effort. The extra preprocessing for creating the input and extra postprocessing to determine the water balance in terms of the physical entities from the MODFLOW cell fluxes per package can be taken care of by a user interface.     

Evaluation of the soil water balance in an alluvial flood plain with a shallow groundwater table

Abstract

Many of the hydrological and ecological functions of alluvial flood plains within watersheds depend on the water flow exchanges between the vadoze soil zone and the shallow groundwater. The water balance of the soil in the flood plain is investigated, in order to evaluate the main hydrological processes that underlie the temporal dynamics of soil moisture and groundwater levels. The soil moisture and the groundwater level in the flood plain were monitored continuously for a three-year period. These data were integrated with the results derived from applying a physically-based numerical model which simulated the variably-saturated vertical water flow in the soil. The analysis indicated that the simultaneous processes of lateral groundwater flow and the vertical recharge from the unsaturated zone caused the observed water table fluctuations. The importance of these flows in determining the rises in the water table varied, depending on soil moisture and groundwater depth before precipitation. The monitoring period included two hydrological years (September 2009–September 2011). About 13% of the precipitation vertically recharged the groundwater in the first year and about 50% in the second. The difference in the two recharge coefficients was in part due to the lower groundwater levels in the recharge season of the first hydrological year, compared to those observed in the second. In the latter year, the shallow groundwater increased the soil moisture in the unsaturated zone due to capillary rise, and so the mean hydraulic conductivity of the unsaturated soil was high. This moisture state of soil favoured a more efficient conversion of infiltrated precipitation into vertical groundwater recharge. The results show that groundwater dynamics in the flood plain are an important source of temporal variability in soil moisture and vertical recharge processes, and this variability must be properly taken into account when the water balance is investigated in shallow groundwater environments.

Citation Pirastru, M. and Niedda, M., 2013. Evaluation of the soil water balance in an alluvial flood plain with a shallow groundwater table. Hydrological Sciences Journal, 58 (4), 898–911.     

A digital procedure for ground water recharge and discharge pattern recognition and rate estimation

A digital procedure to estimate recharge/discharge rates that requires relatively short preparation time and uses readily available data was applied to a setting in central Wisconsin. The method requires only measurements of the water table, fluxes such as stream baseflows, bottom of the system, and hydraulic conductivity to delineate approximate recharge/discharge zones and to estimate rates. The method uses interpolation of the water table surface, recharge/discharge mapping, pattern recognition, and a parameter estimation model. The surface interpolator used is based on the theory of radial basis functions with thin-plate splines. The recharge/discharge mapping is based on a mass-balance calculation performed using MODFLOW. The results of the recharge/discharge mapping are critically dependent on the accuracy of the water table interpolation and the accuracy and number of water table measurements. The recharge pattern recognition is performed with the help of a graphical user interface (GUI) program based on several algorithms used in image processing. Pattern recognition is needed to identify the recharge/discharge zonations and zone the results of the mapping method. The parameter estimation program UCODE calculates the parameter values that provide a best fit between simulated heads and flows and calibration head-and-flow targets. A model of the Buena Vista Ground Water Basin in the Central Sand Plains of Wisconsin is used to demonstrate the procedure.