Mouth of the Mahomet Regional Groundwater Model, Imperial Valley Region of Mason, Tazewell, and Logan Counties, Illinois /

"December 20, 1994."

The cost of uniqueness in groundwater model calibration

Calibration of a groundwater model requires that hydraulic properties be estimated throughout a model domain. This generally constitutes an underdetermined inverse problem, for which a Solution can only be found when some kind of regularization device is included in the inversion process. Inclusion of regularization in the calibration process can be implicit, for example through the use of zones of constant parameter value, or explicit, for example through solution of a constrained minimization problem in which parameters are made to respect preferred values, or preferred relationships, to the degree necessary for a unique solution to be obtained. The cost of uniqueness is this: no matter which regularization methodology is employed, the inevitable consequence of its use is a loss of detail in the calibrated field. This, ill turn, can lead to erroneous predictions made by a model that is ostensibly well calibrated. Information made available as a by-product of the regularized inversion process allows the reasons for this loss of detail to be better understood. In particular, it is easily demonstrated that the estimated value for an hydraulic property at any point within a model domain is, in fact, a weighted average of the true hydraulic property over a much larger area. This averaging process causes loss of resolution in the estimated field. Where hydraulic conductivity is the hydraulic property being estimated, high averaging weights exist in areas that are strategically disposed with respect to measurement wells, while other areas may contribute very little to the estimated hydraulic conductivity at any point within the model domain, this possibly making the detection of hydraulic conductivity anomalies in these latter areas almost impossible. A study of the post-calibration parameter field covariance matrix allows further insights into the loss of system detail incurred through the calibration process to be gained. A comparison of pre- and post-calibration parameter covariance matrices shows that the latter often possess a much smaller spectral bandwidth than the former. It is also demonstrated that, as all inevitable consequence of the fact that a calibrated model cannot replicate every detail of the true system, model-to-measurement residuals can show a high degree of spatial correlation, a fact which must be taken into account when assessing these residuals either qualitatively, or quantitatively in the exploration of model predictive uncertainty. These principles are demonstrated using a synthetic case in which spatial parameter definition is based oil pilot points, and calibration is Implemented using both zones of piecewise constancy and constrained minimization regularization. (C) 2005 Elsevier Ltd. All rights reserved.

Groundwater flow model of the Logan river alluvial aquifer system Josephville, South East Queensland

The study focuses on an alluvial plain situated within a large meander of the Logan River at Josephville near Beaudesert which supports a factory that processes gelatine. The plant draws water from on site bores, as well as the Logan River, for its production processes and produces approximately 1.5 ML per day (Douglas Partners, 2004) of waste water containing high levels of dissolved ions. At present a series of treatment ponds are used to aerate the waste water reducing the level of organic matter; the water is then used to irrigate grazing land around the site. Within the study the hydrogeology is investigated, a conceptual groundwater model is produced and a numerical groundwater flow model is developed from this. On the site are several bores that access groundwater, plus a network of monitoring bores. Assessment of drilling logs shows the area is formed from a mixture of poorly sorted Quaternary alluvial sediments with a laterally continuous aquifer comprised of coarse sands and fine gravels that is in contact with the river. This aquifer occurs at a depth of between 11 and 15 metres and is overlain by a heterogeneous mixture of silts, sands and clays. The study investigates the degree of interaction between the river and the groundwater within the fluvially derived sediments for reasons of both environmental monitoring and sustainability of the potential local groundwater resource. A conceptual hydrogeological model of the site proposes two hydrostratigraphic units, a basal aquifer of coarse-grained materials overlain by a thick semi-confining unit of finer materials. From this, a two-layer groundwater flow model and hydraulic conductivity distribution was developed based on bore monitoring and rainfall data using MODFLOW (McDonald and Harbaugh, 1988) and PEST (Doherty, 2004) based on GMS 6.5 software (EMSI, 2008). A second model was also considered with the alluvium represented as a single hydrogeological unit. Both models were calibrated to steady state conditions and sensitivity analyses of the parameters has demonstrated that both models are very stable for changes in the range of ± 10% for all parameters and still reasonably stable for changes up to ± 20% with RMS errors in the model always less that 10%. The preferred two-layer model was found to give the more realistic representation of the site, where water level variations and the numerical modeling showed that the basal layer of coarse sands and fine gravels is hydraulically connected to the river and the upper layer comprising a poorly sorted mixture of silt-rich clays and sands of very low permeability limits infiltration from the surface to the lower layer. The paucity of historical data has limited the numerical modelling to a steady state one based on groundwater levels during a drought period and forecasts for varying hydrological conditions (e.g. short term as well as prolonged dry and wet conditions) cannot reasonably be made from such a model. If future modelling is to be undertaken it is necessary to establish a regular program of groundwater monitoring and maintain a long term database of water levels to enable a transient model to be developed at a later stage. This will require a valid monitoring network to be designed with additional bores required for adequate coverage of the hydrogeological conditions at the Josephville site. Further investigations would also be enhanced by undertaking pump testing to investigate hydrogeological properties in the aquifer.

Improved modeling of groundwater recharge in agricultural watersheds using a combination of crop model and remote sensing

For improved water management and efficiency of use in agriculture, studies dealing with coupled crop-surface water-groundwater models are needed. Such integrated models of crop and hydrology can provide accurate quantification of spatio-temporal variations of water balance parameters such as soil moisture store, evapotranspiration and recharge in a catchment. Performance of a coupled crop-hydrology model would depend on the availability of a calibrated crop model for various irrigated/rainfed crops and also on an accurate knowledge of soil hydraulic parameters in the catchment at relevant scale. Moreover, such a coupled model should be designed so as to enable the use/assimilation of recent satellite remote sensing products (optical and microwave) in order to model the processes at catchment scales. In this study we present a framework to couple a crop model with a groundwater model for applications to irrigated groundwater agricultural systems. We discuss the calibration of the STICS crop model and present a methodology to estimate the soil hydraulic parameters by inversion of crop model using both ground and satellite based data. Using this methodology we demonstrate the feasibility of estimation of potential recharge due to spatially varying soil/crop matrix.

Groundwater balance research to improve the accuracy of drawdown predictions in the Surat Basin, Queensland

Groundwater modelling studies rely on an accurate determination of inputs and outputs that make up the water balance. Often there is large uncertainty associated with estimates of recharge and unmetered groundwater use. This can translate to equivalent uncertainty in the forecasting of sustainable yields, impacts of extraction, and susceptibility of groundwater dependent ecosystems. In the case of Coal Seam Gas, it is important to characterise the temporal and special distribution of depressurisation in the reservoir and how this may or may not extend to the adjacent aquifers. A regional groundwater flow model has been developed by the Queensland Government to predict drawdown impacts due to Coal Seam Gas activities in the Surat basin. This groundwater model is undergoing continued refinement and there is currently scope to address some of the key areas of uncertainty including better quantification of groundwater recharge and unmetered groundwater extractions. Research is currently underway to improve the accuracy of estimates of both of these components of the groundwater balance in order to reduce uncertainty in predicted groundwater drawdowns due to CSG activities.

Contribution of geophysical surveys to groundwater modelling of a porous aquifer in semiarid Niger: An overview

Subsurface geophysical surveys were carried out using a large range of methods in an unconfined sandstone aquifer in semiarid south-western Niger for improving both the conceptual model of water flow through the unsaturated zone and the parameterization of numerical a groundwater model of the aquifer. Methods included: electromagnetic mapping, electrical resistivity tomography (ERT), resistivity logging, time domain electromagnetic sounding (TDEM), and magnetic resonance sounding (MRS). Analyses of electrical conductivities, complemented by geochemical measurements, allowed us to identify preferential pathways for infiltration and drainage beneath gullies and alluvial fans. The mean water content estimated by MRS (13%) was used for computing the regional groundwater recharge from long-term change in the water table. The ranges in permeability and water content obtained with MRS allowed a reduction of the degree of freedom of aquifer parameters used in groundwater modelling.

Modeling the vulnerability of an urban groundwater system due to the combined impacts of climate change and management scenarios

Climate change impact on a groundwater-dependent small urban town has been investigated in the semiarid hard rock aquifer in southern India. A distributed groundwater model was used to simulate the groundwater levels in the study region for the projected future rainfall (2012-32) obtained from a general circulation model (GCM) to estimate the impacts of climate change and management practices on groundwater system. Management practices were based on the human-induced changes on the urban infrastructure such as reduced recharge from the lakes, reduced recharge from water and wastewater utility due to an operational and functioning underground drainage system, and additional water extracted by the water utility for domestic purposes. An assessment of impacts on the groundwater levels was carried out by calibrating a groundwater model using comprehensive data gathered during the period 2008-11 and then simulating the future groundwater level changes using rainfall from six GCMs Institute of Numerical Mathematics Coupled Model, version 3.0 (INM-CM. 3.0); L'Institut Pierre-Simon Laplace Coupled Model, version 4 (IPSL-CM4); Model for Interdisciplinary Research on Climate, version 3.2 (MIROC3.2); ECHAM and the global Hamburg Ocean Primitive Equation (ECHO-G); Hadley Centre Coupled Model, version 3 (HadCM3); and Hadley Centre Global Environment Model, version 1 (HadGEM1)] that were found to show good correlation to the historical rainfall in the study area. The model results for the present condition indicate that the annual average discharge (sum of pumping and natural groundwater outflow) was marginally or moderately higher at various locations than the recharge and further the recharge is aided from the recharge from the lakes. Model simulations showed that groundwater levels were vulnerable to the GCM rainfall and a scenario of moderate reduction in recharge from lakes. Hence, it is important to sustain the induced recharge from lakes by ensuring that sufficient runoff water flows to these lakes.

Evaluation of effective groundwater recharge of freshwater lens in small islands by the combined modeling of geoelectrical data and water heads

In small islands, a freshwater lens can develop due to the recharge induced by rain. Magnitude and spatial distribution of this recharge control the elevation of freshwater and the depth of its interface with salt water. Therefore, the study of lens morphology gives useful information on both the recharge and water uptake due to evapotranspiration by vegetation. Electrical resistivity tomography was applied on a small coral reef island, giving relevant information on the lens structure. Variable density groundwater flow models were then applied to simulate freshwater behavior. Cross validation of the geoelectrical model and the groundwater model showed that recharge exceeds water uptake in dunes with little vegetation, allowing the lens to develop. Conversely, in the low-lying and densely vegetated sectors, where water uptake exceeds recharge, the lens cannot develop and seawater intrusion occurs. This combined modeling method constitutes an original approach to evaluate effective groundwater recharge in such environments.
[Comte, J.-C., O. Banton, J.-L. Join, and G. Cabioch (2010), Evaluation of effective groundwater recharge of freshwater lens in small islands by the combined modeling of geoelectrical data and water heads, Water Resour. Res., 46, W06601, doi:10.1029/2009WR008058.]

Modelling the response of fresh groundwater to climate and vegetation changes in coral islands

In coral islands, groundwater is a crucial freshwater resource for terrestrial life, including human water supply. Response of the freshwater lens to expected climate changes and subsequent vegetation alterations is quantified for Grande Glorieuse, a low-lying coral island in the Western Indian Ocean. Distributed models of recharge, evapotranspiration and saltwater phytotoxicity are integrated into a variable-density groundwater model to simulate the evolution of groundwater salinity. Model results are assessed against field observations including groundwater and geophysical measurements. Simulations show the major control currently exerted by the vegetation with regards to the lens morphology and the high sensitivity of the lens to climate alterations, impacting both quantity and salinity. Long-term changes in mean sea level and climatic conditions (rainfall and evapotranspiration) are predicted to be responsible for an average increase in salinity approaching 140 % (+8 kg m-3) when combined. In low-lying areas with high vegetation density, these changes top +300 % (+10 kg m-3). However, due to salinity increase and its phytotoxicity, it is shown that a corollary drop in vegetation activity can buffer the alteration of fresh groundwater. This illustrates the importance of accounting for vegetation dynamics to study groundwater in coral islands.

Aluminium smelter wastewater : groundwater fluoride fate

This study assessed the viability of utilising groundwater systems to manage Portland Aluminium's trade waste water stream. A mathematical groundwater model was created to assess laboratory and field results. Fluoride was modelled for spatio-temporal longevity in the system. All models indicate that fluoride will ultimately exit the site to the seaward boundary.

Numerical groundwater flow modeling in the Wakal River basin, India

Increasing dependence on groundwater in the Wakal River basin, India, jeopardizes water supply sustainability. A numerical groundwater model was developed to better understand the aquifer system and to evaluate its potential in terms of quantity and replenishment. Potential artificial recharge areas were delineated using landscape and hydrogeologic parameters, Geographic Information System (GIS), and remote sensing. Groundwater models are powerful tools for recharge estimation when transmissivity is known. Proper recharge must be applied to reproduce field-measured heads. The model showed that groundwater levels could decline significantly if there are two drought years in every four years that result in reduced recharge, and groundwater withdrawal is increased by 15%. The effect of such drought is currently uncertain however, because runoff from the basin is unknown. Remote sensing and GIS revealed areas with slopes less than 5%, forest cover, and Normalized Difference Vegetative Index greater than 0.5 that are suitable recharge sites.

基于盆地地下水模型与山区地表水模型耦合的济源盆地地下水模拟研究

The aim of this study is to increase the precision of groundwater modeling. The way is use the distributed model calculate the mountain basin groundwater lateral discharge and the river runoff. With appropriate technique help, the groundwater model can couple the distributed model results. This paper’s study object is makeing the distributed hydrological model HEC-HMS coupled to the popular groundwater model Visual MODFLOW. The application example is Jiyuan basin which is a typical basin of North China. HEC-HMS can calculate the surface runoff and subsurface runoff at mountain-pass. The subsurface runoff can turn to recharge well straightly. The water level - runoff course and Trial method is used to back analyze the parameters of surface runoff to Visual MODFLOW. So the distributed hydrological model can coupled to the groundwater model. The research proved that base on couple the distributed surface water model the groundwater model’s results are notability improved. The example is Jiyuan basin where use the distributed model coupled to the groundwater model. On the base of the coupled model applied to Jiyuan basin groundwater modeling. The paper estimates the groundwater change in the study area. Then, by use the water resources integrated planning results, the article calculate the basin groundwater can be development and utilization quantity and potential.