Stochastic modelling of groundwater flow and solute transport in aquifers

This paper presents an introductory overview of recently developed stochastic theories for tackling spatial variability problems in predicting groundwater flow and solute transport. Advantages and limitations of the theories are discussed. Lastly, strategies based on the stochastic approaches to predict solute transport in aquifers are recommended.     

WEKU – a GIS-Supported stochastic model of groundwater residence times in upper aquifers for the supraregional groundwater management

 The supraregional GIS-supported stochastical model, WEKU, for the determination of groundwater residence times in the upper aquifers of large groundwater provinces is presented. Using a two-dimensional analytical model of groundwater flow, groundwater residence times are determined within two extreme cases. In the first case, maximal groundwater residence times are calculated, representing the part of groundwater, that is drained by the main surface water of a groundwater catchment area. In the second case, minimal groundwater residence times for drainage into the nearest surface water are determined. Using explicit distribution functions of the input parameters, mean values as well as potential ranges of variations of the groundwater residence times are derived. The WEKU model has been used for the determination of groundwater residence times throughout Germany. The model results – mean values and deviations of the groundwater velocity and the maximal and minimal groundwater residence times in the upper aquifers – are presented by general maps and discussed in detail. It is shown that the groundwater residence times in the upper aquifer vary regionally, differentiated between less than 1 year and more than 2000 years. Using this information, the time scales can be specified, until measures to remediate polluted groundwater resources may lead to a substantial groundwater quality improvement in the different groundwater provinces of Germany. With respect to its supraregional scale of application, the WEKU model may serve as a useful tool for the supraregional groundwater management on a state, federal or international level. Received: 15 August 1995 · Accepted: 15 October 1995     

Integrated approach for identification of potential groundwater zones in Seethanagaram Mandal of Vizianagaram District,Andhra Pradesh,India

Identifying a good site for groundwater exploration in hard rock terrain is a challenging task. In hard rocks, groundwater occurs in secondary porosity developed due to weathering, fracturing, faulting, etc., which is highly variable within short distance and contributing to near-surface inhomogeneity. In such situations topographic, hydrogeological and geomorphological features provide useful clues for the selection of suitable sites. Initially, based on satellite imagery, topographical, geomorphological and hydrogeological features, an area of about 149 km² was demarcated as a promising zone for groundwater exploration in the hard rock tract of Seethanagaram Mandal, Vizianagaram District, Andhra Pradesh, India. A total of 50 Vertical Electrical Soundings (VES) were carried out using Wenner electrode configuration. An interactive interpretation of the VES data sharpened the information inferred from geomorphological and hydrogeological reconnaissance. Ten sites were recommended for drilling. Drilling with Down-The-Hole Hammer (DTH) was carried out at the recommended sites down to 50 to 70 m depths. The interpreted VES results matched well with the drilled bore well lithologs. The yields of bore wells vary from 900 to 9000 liters per hour (lph).     

Implications of Kali-Hindon inter-stream aquifer water balance for groundwater management in western Uttar Pradesh

The Kali-Hindon inter-stream region extends over an area of 395 km² within the Ganga-Yamuna interfluve. It is a fertile tract for sugarcane cultivation. Groundwater is a primary resource for irrigation and industrial purposes. In recent years, over-exploitation has resulted in an adverse impact on the groundwater regime. In this study, an attempt has been made to calculate a water balance for the Kali-Hindon inter-stream region. Various inflows and outflows to and from the aquifer have been calculated. The recharge due to rainfall and other recharge parameters such as horizontal inflow, irrigation return flow and canal seepage were also evaluated. Groundwater withdrawals, evaporation from the water table, discharge from the aquifer to rivers and horizontal subsurface outflows were also estimated. The results show that total recharge into the system is 148.72 million cubic metres (Mcum), whereas the total discharge is 161.06 Mcum, leaving a deficit balance of −12.34 Mcum. Similarly, the groundwater balance was evaluated for the successive four years. The result shows that the groundwater balance is highly sensitive to variation in rainfall followed by draft through pumpage. The depths to water level are shallow in the canal-irrigated northern part of the basin and deeper in the southern part. The pre-monsoon and post-monsoon water levels range from 4.6 to 17.7 m below ground level (bgl) and from 3.5 to 16.5 m bgl respectively. It is concluded that the groundwater may be pumped in the canal-irrigated northern part, while withdrawals may be restricted to the southern portion of the basin, where intense abstraction has led to rapidly falling water table levels.     

Travel time based approach for the assessment of vulnerability of karst groundwater: the Transit Time Method

Karst aquifers represent important water resources in many parts of the world. Unfortunately, karst aquifers are characterised by high contamination risks. This paper presents a travel time based method for the estimation of karst groundwater vulnerability. It considers (1) physics-based lateral flow within the uppermost weathered zone (epikarst) in a limestone-dominated region and (2) high velocities of vertical infiltration at discrete infiltration points (e.g. sinkholes) or lines (e.g. dry valleys, faults). Consequently, the Transit Time Method honours the actual flow path within the unsaturated zone of a karst aquifer system. A test site in Northern Jordan was chosen for the demonstration of the assessment technique, i.e. the catchment area of the Qunayyah Spring north of the capital Amman. The results demonstrate that zones of highest vulnerability lie within valleys and nearby main fault zones. It also reveals that regions, categorised as protected areas by other methods due to thick unsaturated zones, contribute to a major degree to the total risk.     

Hydrogeological framework and water balance studies in parts of Krishni–Yamuna interstream area, Western Uttar Pradesh, India

The Krishni–Yamuna interstream area is a micro-watershed in the Central Ganga Plain and a highly fertile track of Western Uttar Pradesh. The Sugarcane and wheat are the major crops of the area. Aquifers of Quaternary age form the major source of Irrigation and municipal water supplies. A detailed hydrogeological investigation was carried out in the study area with an objective to assess aquifer framework, groundwater quality and its resource potential. The hydrogeological cross section reveals occurrence of alternate layers of clay and sand. Aquifer broadly behaves as a single bodied aquifer down to the depth of 100 m bgl (metre below ground level) as the clay layers laterally pinch out. The depth to water in the area varies between 5 and 16.5 m bgl. The general groundwater flow direction is from NE to SW with few local variations. An attempt has been made to evaluate groundwater resources of the area. The water budget method focuses on the various components contributing to groundwater flow and groundwater storage changes. Changes in ground water storage can be attributed to rainfall recharge, irrigation return flow and ground water inflow to the basin minus baseflow (ground water discharge to streams or springs), evapotranspiration from ground water, pumping and ground water outflow from the basin. The recharge is obtained in the study area using Water table fluctuation and Tritium methods. The results of water balance study show that the total recharge in to the interstream region is of the order of 185.25 million m³ and discharge from the study area is of the order of 203.24 million m³, leaving a deficit balance of −17.99 million m³. Therefore, the present status of groundwater development in the present study area has acquired the declining trend. Thus, the hydrogeological analysis and water balance studies shows that the groundwater development has attained a critical state in the region.     

Application of MODFLOW and geographic information system to groundwater flow simulation in North China Plain, China

MODFLOW is a groundwater modeling program. It can be compiled and remedied according to the practical applications. Because of its structure and fixed data format, MODFLOW can be integrated with Geographic Information Systems (GIS) technology for water resource management. The North China Plain (NCP), which is the politic, economic and cultural center of China, is facing with water resources shortage and water pollution. Groundwater is the main water resource for industrial, agricultural and domestic usage. It is necessary to evaluate the groundwater resources of the NCP as an entire aquifer system. With the development of computer and internet information technology it is also necessary to integrate the groundwater model with the GIS technology. Because the geological and hydrogeological data in the NCP was mainly in MAPGIS format, the powerful function of GIS of disposing of and analyzing spatial data and computer languages such as Visual C and Visual Basic were used to define the relationship between the original data and model data. After analyzing the geological and hydrogeological conditions of the NCP, the groundwater flow numerical simulation modeling was constructed with MODFLOW. On the basis of GIS, a dynamic evaluation system for groundwater resources under the internet circumstance was completed. During the process of constructing the groundwater model, a water budget was analyzed, which showed a negative budget in the NCP. The simulation period was from 1 January 2002 to 31 December 2003. During this period, the total recharge of the groundwater system was 49,374 × 10⁶ m³ and the total discharge was 56,530 × 10⁶ m³ the budget deficit was −7,156 × 10⁶ m³. In this integrated system, the original data including graphs and attribution data could be stored in the database. When the process of evaluating and predicting groundwater flow was started, these data were transformed into files that the core program of MODFLOW could read. The calculated water level and drawdown could be displayed and reviewed online.     

Problems of soil and groundwater pollution in the disposal of “marble” slurries in NW Sicily

This work deals with disposal of slurries generated during the cutting and polishing processes of slabs of decorative sedimentary carbonate rocks in the north western Sicily. At present, they are used as fillers of dismantled quarries near the sawmills and, as a final step of reclamation, are covered with earth layers. In spite of such inexpensive solution, there is lack of knowledge about the composition of the waste. In order to assess if there is any threat for the environment and to suggest indications for alternative solutions, such as recycling or inactivation processes, the slurries were analysed by XR diffraction, simultaneous thermal analysis, ICP/MS, ionic chromatography, FTIR, UV-Vis, COD and TOC measurements, grain size analysis. Results indicate that the slurries can threaten the groundwater, because of the high chemical oxygen demand; furthermore they can modify the mechanism of groundwater recharge, because of their grain size distribution. Some laboratory tests show that, even in very aggressive conditions, the solid pollutants persist in the waste and slowly release into water the products of their degradation. The slurry therefore should be subjected to inactivation treatment before disposal or, alternatively, recycled as secondary raw material for a suitable process.     

Impacts of at-site wastewater disposal systems on the groundwater aquifer in arid regions: case of Sfax City, Southern Tunisia

Groundwater in Sfax City (Tunisia) has been known since the beginning of the century for its deterioration in quality, as a result of wastewater recharge into the aquifer. An average value of 12 × 10⁶ m³ of untreated wastewater reaches the groundwater aquifer each year. This would result not only in a chemical and biological contamination of the groundwater, but also in an increase of the aquifer piezometric level. Quantitative impacts were evaluated by examining the groundwater piezometric level at 57 surface wells and piezometers. The survey showed that, during the last two decades, the groundwater level was ever increasing in the urban area with values reaching 7 m in part; and decreasing in Sidi Abid (agricultural area) with values exceeding −3 m. Groundwater samples for chemical and microbial analysis were collected from 41 wells spread throughout the study area. Results showed significantly elevated levels of sodium, chlorides, nitrates and coliform bacteria all over the urban area. High levels (NO₃: 56–254 mg/l; Na >1,500 mg/l; Coliforms >30/100 ml) can be related to more densely populated areas with a higher density of pit latrine and recharge wells. Alternatively results showed a very variable chemical composition of groundwater, e.g. electrical conductivity ranges from 4,040 to19,620 μs/cm and the dry residual varies between 1.4 and 14 g/l with concentrations increasing downstream. Furthermore a softening of groundwater in Set Ezzit (highly populated sector) was observed.