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.     

P‐wave transmission across fractures with nonlinear deformational behaviour

Stress wave attenuation across fractured rock masses is a great concern of underground structure safety. When the wave amplitude is large, fractures experience nonlinear deformation during the wave propagation. This paper presents a study on normal transmission of P‐wave across parallel fractures with nonlinear deformational behaviour (static Barton–Bandis model). The results show that the magnitude of transmission coefficient is a function of incident wave amplitude, nondimensional fracture spacing and number of fractures. Two important indices of nondimensional fracture spacing are identified, and they divide the area of nondimensional fracture spacing into three parts (individual fracture area, transition area and small spacing area). In the different areas, the magnitude of transmission coefficient has different trends with nondimensional fracture spacing and number of fractures. In addition, the study reveals that under some circumstances, the magnitude of transmission coefficient increases with increasing number of fractures, and is larger than 1. Copyright © 2006 John Wiley & Sons, Ltd.     

Seasonal variations in prey selection by estuarine black-headed gulls (Larus ridibundus)

The number of black-headed gulls (Larus ridibundus) in the Clyde Estuary is large. In summer the average density has reached 1350 gulls km^?2 and in winter 180 gulls km^?2. This paper compares prey selection and feeding efficiency in gulls during summer and winter on tidal flats, and considers how seasonal differences may be adaptations to cope with seasonal changes in prey availability.Gross and net rates of energy intake were highest in summer because gulls captured more of the polychaete N. diversicolor than the amphipod C. volutator. In winter, gulls selected for C. volutator and therefore an energetically less profitable diet. Throughout the year gulls selected more C. volutator relative to N. diversicolor than expected on energetic grounds and so apparently did not maximize potential net rate of energy intake.Gulls used three techniques to capture prey and made most intensive use of the ‘crouch’ technique. Crouching gulls attained a much higher net rate of energy intake than ‘upright’ or ‘paddling’ gulls.A log-linear model showed that (a) season, water depth and gull density determined feeding technique and (b) feeding technique and season independently determined foraging success and prey selection. Thus gull density and water depth acted on prey selection through imposed variations in feeding technique.Reasons for gulls selecting energetically unprofitable C. volutator and for the use of several distinct feeding techniques are discussed.     

Generalized strain probing of constitutive models

Advanced material constitutive models are used to describe complex soil behaviour. These models are often used in the solution of boundary value problems under general loading conditions. Users and developers of constitutive models need to methodically investigate the represented soil response under a wide range of loading conditions. This paper presents a systematic procedure for probing constitutive models. A general incremental strain probe, 6D hyperspherical strain probe (HSP), is introduced to examine rate‐independent model response under all possible strain loading conditions. Two special cases of HSP, the true triaxial strain probe (TTSP) and the plane‐strain strain probe (PSSP), are used to generate 3‐D objects that represent model stress response to probing. The TTSP, PSSP and general HSP procedures are demonstrated using elasto‐plastic models. The objects resulting from the probing procedure readily highlight important model characteristics including anisotropy, yielding, hardening, softening and failure. The PSSP procedure is applied to a Neural Network (NN) based constitutive model. It shows that this probing is especially useful in understanding NN constitutive models, which do not contain explicit functions for yield surface, hardening, or anisotropy. Copyright © 2004 John Wiley & Sons, Ltd.     

A micromechanical finite element model for linear and damage‐coupled viscoelastic behaviour of asphalt mixture

This study presents a finite element (FE) micromechanical modelling approach for the simulation of linear and damage‐coupled viscoelastic behaviour of asphalt mixture. Asphalt mixture is a composite material of graded aggregates bound with mastic (asphalt and fine aggregates). The microstructural model of asphalt mixture incorporates an equivalent lattice network structure whereby intergranular load transfer is simulated through an effective asphalt mastic zone. The finite element model integrates the ABAQUS user material subroutine with continuum elements for the effective asphalt mastic and rigid body elements for each aggregate. A unified approach is proposed using Schapery non‐linear viscoelastic model for the rate‐independent and rate‐dependent damage behaviour. A finite element incremental algorithm with a recursive relationship for three‐dimensional (3D) linear and damage‐coupled viscoelastic behaviour is developed. This algorithm is used in a 3D user‐defined material model for the asphalt mastic to predict global linear and damage‐coupled viscoelastic behaviour of asphalt mixture. For linear viscoelastic study, the creep stiffnesses of mastic and asphalt mixture at different temperatures are measured in laboratory. A regression‐fitting method is employed to calibrate generalized Maxwell models with Prony series and generate master stiffness curves for mastic and asphalt mixture. A computational model is developed with image analysis of sectioned surface of a test specimen. The viscoelastic prediction of mixture creep stiffness with the calibrated mastic material parameters is compared with mixture master stiffness curve over a reduced time period. In regard to damage‐coupled viscoelastic behaviour, cyclic loading responses of linear and rate‐independent damage‐coupled viscoelastic materials are compared. Effects of particular microstructure parameters on the rate‐independent damage‐coupled viscoelastic behaviour are also investigated with finite element simulations of asphalt numerical samples. Further study describes loading rate effects on the asphalt viscoelastic properties and rate‐dependent damage behaviour. Copyright © 2006 John Wiley & Sons, Ltd.     

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.     

Hydrochemical characteristics and the effects of irrigation on groundwater quality in Harran Plain,GAP Project,Turkey

Improper design, faulty planning, mismanagement and incorrect operation of irrigation schemes are the principle reasons for the deterioration of groundwater quality in a large number of countries, in particular in semi-arid and arid regions. The aim of this study is to determine the dimensions of groundwater quality after surface irrigation was begun in the semi-arid Harran Plain. Physical and chemical parameters of the groundwater including pH, temperature, electrical conductivity (EC), sodium, potassium, calcium, magnesium, chloride, bicarbonate, sulphate, nitrate, nitrite, ammonium, total phosphorus, total organic carbon and turbidity were determined monthly during the 2006 water year. The quality of the groundwater in the study area was assessed hydrochemically in order to determine its suitability for human consumption and agricultural purposes. In the general plain, the EC values measured were considerably above the guide level of 650 μS/cm, while nitrate in particular was found in almost all groundwater samples to be significantly above the maximum admissible concentration of 50 mg/l for the quality of water intended for human consumption as per the international and national standards. Total hardness reveals that a majority of the groundwater samples fall in the very hard water category. Interpretation of analytical data shows that Ca–HCO₃ and Ca–SO₄ are the dominant hydrochemical facies in the study area.