Inverse groundwater modelling in the Willunga Basin,South Australia

A new inverse technique for modelling groundwater flow, based on a functional minimization technique, has been used to calibrate a groundwater flow model of a subregion of the Port Willunga aquifer within the Willunga Basin in South Australia. The Willunga Basin is the location of extensive viticulture, irrigated primarily by groundwater, the levels and quality of which have declined significantly over the last 40 years. The new method is able to generate estimates of transmissivity, storativity and groundwater recharge over the whole subregion as a time-varying continuous surface; previous methods estimate local discrete parameter values at specific times. The new method has also been shown to produce accurate head values for the subregion and very good estimates of groundwater recharge. Its ultimate goal will be to provide a new and invaluable tool for significantly improved groundwater resource management. Supported in part by US National Science Foundation grants, DMS-0107492 and DMS-0079478.     

Aquifer development planning to supply a seaside resort: a case study in Goa,India

Using the hydrogeological and socio-economic data derived from a European Commission research project on the measurement, monitoring and sustainability of the coastal environment, two optimization models have been applied to satisfy the future water resources needs of the coastal zone of Bardez in Goa, India. The number of tourists visiting Goa since the 1970s has risen considerably, and roughly a third of them go to Bardez taluka, prompting growth in the tourist-related infrastructure in the region. The optimization models are non-linear mixed integer models that have been solved using GAMS/DICOPT++ commercial software. Optimization models were used, firstly, to indicate the most suitable zones for building seaside resorts and wells to supply the tourist industry with an adequate amount of water, and secondly, to indicate the best location for wells to adequately supply pre-existing hotels. The models presented will help to define the optimal locations for the wells and the hydraulic infrastructures needed to satisfy demand at minimum cost, taking into account environmental constraints such as the risk of saline intrusion.

---

Résumé A l’aide de données hydrogéologiques et socio-économiques, dérivées d’un projet de recherche de la Commission Européenne sur l’étude, la surveillance et la durabilité de l’environnement c?tier, deux modèles d’optimisation ont été appliqués en vue de satisfaire les besoins futurs en ressources en eau de la zone c?tière de Bardez au Goa en Inde. Le nombre de touristes qui visitent le Goa depuis les années 70 a considérablement augmenté et environ un tiers d’entre eux se rend à Bardez taluka, ce qui encourage la croissance des infrastructures liées au tourisme dans la région. Les modèles d’optimisation sont des modèles non linéaires mixtes en nombres entiers qui ont été résolus en utilisant le logiciel commercial GAMS/DICOPT++. Les modèles d’optimisation ont été utilisés pour indiquer premièrement les zones les plus appropriées pour construire des stations balnéaires et des puits pour approvisionner en quantité d’eau suffisante l’industrie touristique, et deuxièmement pour indiquer la meilleure localisation de puits pour l’alimentation des h?tels préexistants. Les modèles présentés aideront à définir les localisations optimales des puits et des infrastructures hydrauliques nécessaires à la satisfaction de la demande à un co?t minimum, tout en prenant en compte les contraintes environnementales, tel que le risque d’intrusion saline.

---

Resumen Se han aplicado dos modelos de optimización en base a datos socio-económicos e hidrogeológicos derivados de un proyecto de investigación de la Comisión Europea sobre medición, monitoreo y sostenibilidad del ambiente costero para satisfacer las necesidades futuras de recursos hídricos de la zona costera de Bardez en Goa, India. El número de turistas que visita Goa desde la década de 1970′s ha subido considerablemente y cerca de un tercio de ellos van a Bardez taluka impulsando el crecimiento en la infraestructura turística de la región. Los modelos de optimización son modelos de números enteros mixtos no-lineales que se han resuelto usando el programa comercial GAMS/DICOPT++. Los modelos de optimización se usaron, primero para indicar las zonas más adecuadas para la construcción de centros de diversión en las márgenes del océano y pozos para abastecer la industria turística con una adecuada cantidad de agua, y segundo, para indicar la mejor localización de pozos para el abastecimiento de hoteles pre-existentes. Los modelos que se presentan ayudarán a definir las localizaciones óptimas para los pozos y las infraestructuras hidráulicas necesarias para satisfacer la demanda al mínimo costo, tomando en consideración restricciones ambientales tal como el riesgo de intrusión salina.

     

Effects of climate change on the groundwater system in the Grote-Nete catchment,Belgium

The effects of climate change on the groundwater systems in the Grote-Nete catchment, Belgium, covering an area of 525 km², is modeled using wet (greenhouse), cold or NATCC (North Atlantic Thermohaline Circulation Change) and dry climate scenarios. Low, central and high estimates of temperature changes are adopted for wet scenarios. Seasonal and annual water balance components including groundwater recharge are simulated using the WetSpass model, while mean annual groundwater elevations and discharge are simulated with a steady-state MODFLOW groundwater model. WetSpass results for the wet scenarios show that wet winters and drier summers are expected relative to the present situation. MODFLOW results for wet high scenario show groundwater levels increase by as much as 79 cm, which could affect the distribution and species richness of meadows. Results obtained for cold scenarios depict drier winters and wetter summers relative to the present. The dry scenarios predict dry conditions for the whole year. There is no recharge during the summer, which is mainly attributed to high evapotranspiration rates by forests and low precipitation. Average annual groundwater levels drop by 0.5 m, with maximum of 3.1 m on the eastern part of the Campine Plateau. This could endanger aquatic ecosystem, shrubs, and crop production.     

The geochemical study of fluvio-lacustrine aquifers in the Kathmandu Basin (Nepal) and the implications for the mobilization of arsenic

The geochemical analyses of fluvial-lacustrine aquifer sediments of the Kathmandu Valley have been made as a step in assessing the environment for the mobilization of arsenic in groundwater. Elements measured by X-ray fluorescence (XRF) include 4 major oxides (Fe₂O₃, TiO₂, CaO, P₂O₅) and 14 trace elements (As, Pb, Zn, Cu, Ni, Cr, V, Sr, Y, Nb, Zr, Th, Sc and TS). Elution tests of 15 selected core samples were also carried out to determine the potential leaching of arsenic from the sediments. The XRF results show that average bulk concentrations of the major oxides and trace elements are similar to modern unconsolidated sediments and average upper continental crust. However, the abundance of elements varies with grain size, with higher concentrations in finer-size grades. Variations in elemental abundances within the basin are strongly controlled by sediment facies. The elution tests show that greater amounts of arsenic are generally eluted from the fine sediments, although the rates are variable. The results overall suggest that As concentrations in the bulk sediments are not a controlling factor for elevated As in the Kathmandu Basin groundwater, and the roles of other factors such as redox conditions and organic matter contents are likely to be more significant.     

Effects of intraborehole flow on groundwater age distribution

Environmental tracers are used to estimate groundwater ages and travel times, but the strongly heterogeneous nature of many subsurface environments can cause mixing between waters of highly disparate ages, adding additional complexity to the age-estimation process. Mixing may be exacerbated by the presence of wells because long open intervals or long screens with openings at multiple depths can transport water and solutes rapidly over a large vertical distance. The effect of intraborehole flow on groundwater age was examined numerically using direct age transport simulation coupled with the Multi-Node Well Package of MODFLOW. Ages in a homogeneous, anisotropic aquifer reached a predevelopment steady state possessing strong depth dependence. A nonpumping multi-node well was then introduced in one of three locations within the system. In all three cases, vertical transport along the well resulted in substantial changes in age distributions within the system. After a pumping well was added near the nonpumping multi-node well, ages were further perturbed by a flow reversal in the nonpumping multi-node well. Results indicated that intraborehole flow can substantially alter groundwater ages, but the effects are highly dependent on local or regional flow conditions and may change with time.     

Spatial assessment of groundwater use potential for irrigation in Teesta Barrage Project in Bangladesh

In Bangladesh, development of the groundwater resource for irrigation is a vital component of the government’s agricultural strategy to attain food self-sufficiency. Amidst reports of falling groundwater levels in many parts of Bangladesh, the potential of groundwater use has been investigated in the Teesta Barrage Project (TBP) in which large-scale groundwater development for dry-season irrigation has taken place in the recent past. Several techniques and tools have been applied such as the combination of analysis of groundwater hydrographs and mathematical modelling to derive key hydrogeological variables, calculation of net irrigation requirement and the use of geographical information systems. The results show that the economically attractive high-yielding variety (HYV) Boro (dry season) rice cultivation during the groundwater irrigation season may not be sustained in large parts of the project area if the current trends in abstraction are continued. However, due to spatial variation in abstraction, nine thanas (sub-districts)—out of a total of 21 in the project area—may still be able to expand groundwater-irrigated cropland and a groundwater-use potential of 40 mm/year may be created if deep-set shallow tubewells are used by the farmers to abstract groundwater. A structured approach, based on zoning of potential areas, is recommended for groundwater development and use.

---

Resumen En Bangladesh, el desarrollo de las aguas subterráneas para riego es un componente vital del la estrategia agrícola del Gobierno para obtener la auto-suficiencia alimentaria. Un informe intermedio evidencia un descenso de los niveles de aguas subterráneas en muchas partes de Bangladesh, el uso potencial de aguas subterránea ha sido investigado en el Proyecto Tessta Barrage (TBP) en el cual ha tenido lugar el desarrollo de las aguas subterráneas a gran escala para el riego en la estación seca en un pasado reciente. Se ha aplicado algunas técnicas y herramientas, como la combinación del análisis de datos gráficos de aguas subterráneas y modelización matemática para obtener variables hidrogeológicas clave, el cálculo de la necesidad neta de riego y el uso de Sistemas de Información Geográfica. El resultado muestra que el cultivo, económicamente atractivo, de arroz HYV Boro (estación seca) durante la estación de riego con aguas subterráneas puede no ser sostenible en gran parte del área del proyecto si continúan las tendencias actuales en las extracciones. Sin embargo, debido a la variación espacial en la extracción, nueve thanas (subdistritos)—de un total de 21 en el área del Proyecto—pueden todavía sufrir un aumento de la tierra cultivable regada con aguas subterráneas y que se puede asumir un uso potencial de agua subterránea de 40 mm/estación si los granjeros utilizan sondeos someros instalados en fosas para la extracción. Se recomienda una aproximación estructurada para el desarrollo y uso de las aguas subterráneas, basada en la zonificación de áreas potenciales.

---

Résumé Au Bangladesh, l’exploitation des ressources en eau souterraines pour l’irrigation est une composante vitale de la stratégie agricole gouvernementale pour atteindre l’auto-suffisance alimentaire. Parmi les baisses de niveaux piézométriques rapportées sur de nombreuses régions du Bangladesh, le potentiel d’utilisation des eaux souterraines a été étudié sur le Projet de Barrage de Teesta (TBP), où, dans un passé récent, les eaux souterraines ont été massivement exploitées pour l’irrigation en saison sèche. De nombreux outils et techniques ont été utilisés, comme l’analyse combinée des chroniques piézométriques et des modèles mathématiques pour dégager les variables hydrogéologiques clés, le calcul des besoins nets pour l’irrigation, et l’utilisation des systèmes d’information géographique. Les résultats démontrent que la culture en saison d’irrigation de la variété de riz HYV Boro (saison sèche), économiquement attrayante, n’est potentiellement pas viable sur de nombreux secteurs de la zone étudiée si les prélèvements perdurent selon la tendance actuelle. Cependant, du fait de la variation spatiale des prélèvements, neuf thanas (sous-districts)—sur un total de 21 sur le secteur d’étude—pourraient encore augmenter les surfaces irriguées, et une disponibilité en eau de 40 mm/saison pourrait être générée si les agriculteurs utilisaient des puits peu profonds pour exploiter les ressources souterraines. Une approche structurée, basée sur le zonage des secteurs potentiels, est recommandée pour l’exploitation et l’utilisation des eaux souterraines.

     

2D magnetic resonance tomography applied to karstic conduit imaging

Karstic conduits play a crucial role for water supply in many parts of the world. However, the imaging of such targets is generally a difficult task for most geophysical methods. Magnetic Resonance Sounding (MRS) is a geophysical method designed for imaging of water bearing structures. Initially, MRS was developed for characterizing horizontally stratified aquifers. However, when applying a 1D MRS measuring setup to the imaging of 2D–3D targets, the size of which may be much smaller than the loop, the accuracy and the lateral resolution may not be sufficient. We have studied the possibility of simultaneously processing several MRS aligned along a profile to perform a Magnetic Resonance Tomography (MRT). This work emphasizes the gain of resolution for 2D–3D imagery of MRT versus the interpolation of 1D inversion results of MRS along the same profile. Numerical modelling results show that the MRT response is sensitive to the size and location of the 2D target in the subsurface. Sensitivity studies reveal that by using the coincident transmitting/receiving (TX/RX) setup and shifting the loop around the anomaly area, the depth, section and position of a single karstic conduit with a size smaller than the MRS loop size can be resolved. The accuracy of the results depends on the noise level and signal level, the latter parameter being linked to the depth and volume of the karstic conduit and the water content in the limestone matrix. It was shown that when applying MRT to the localization of 2D anomalies such as karstic conduits, the inclination of the geomagnetic field, the orientation of the MRT profile and the angle of crossover of the conduit by the MRT profile must be taken into account. Otherwise additional errors in interpretation should be expected. A 2D inversion scheme was developed and tested. Both numerical and experimental results confirm the efficiency of the developed approach.     

Comments on “Steady- and transient-state inversion in hydrogeology by successive flux estimation” by P. Pasquier and D. Marcotte

Pasquier and Marcotte [Pasquier P, Marcotte D. Steady- and transient-state inversion in hydrogeology by successive flux estimation. Adv Wat Res 2006;29:1934–52] propose some modifications to the Comparison Model Method (CMM), in order to apply it to transient 3D ground water flow data for conductivity identification. We present some remarks on that paper to improve the comprehension of the basic features of the CMM and of the real value of the novelties introduced by Pasquier and Marcotte.