Salinization of groundwater in the Nefzawa oases region,Tunisia: results of a regional-scale hydrogeologic approach

Groundwater pumped from the semi-confined Complexe Terminal (CT) aquifer is an important production factor in irrigated oases agriculture in southern Tunisia. A rise in the groundwater salinity has been observed as a consequence of increasing abstraction from the aquifer during the last few decades. All sources of contamination were investigated using hydrochemical data available from the 1980s. Water samples were taken from drains and observation wells tapping both the CT and the phreatic aquifers and analyzed with regard to chemistry, temperature, isotopes and other environmental tracers. Local salinization mechanisms are suggested, i.e. the upwelling of saline water from the underlying, confined Continental Intercalaire (CI) aquifer, as well as backflow of agricultural drainage water. At this stage, the main salt pan, the Chott el Djerid, is not a contamination source. A finite difference model was also developed to simulate groundwater flow and contaminant transport in the oases. Calibration for the period 1950–2000 was carried out in order to adjust geological and chemical system parameters. The simulation of planned extraction projects predicts a worsening of the present situation. Maintenance of the present abstraction regime will not reduce or stop the salinity increase.

---

Résumé L’eau souterraine pompée dans l’aquifère semi-captif du Complexe Terminal (CT) est un facteur de production important pour l’agriculture des oasis irriguées du sud de la Tunisie. Une augmentation de la salinité de l’eau souterraine a été considérée comme la conséquence de l’augmentation des prélèvements dans l’aquifère au cours des dernières décades. Toutes les sources de contamination ont été étudiées à l’aide de données hydro-chimiques disponibles depuis les années 80. Des échantillons d’eau ont été prélevés dans des drains et des puits d’observation qui captent à la fois le CT et les aquifères phréatiques; la chimie, la température, les isotopes ainsi que d’autres traceurs environnementaux ont été analysés. Des processus locaux de salinisation sont proposés, comme par exemple la remontée d’eau salée en provenance de l’aquifère captif du Continental Intercalaire (CI) sous-jacent, ainsi que le reflux des eaux du drainage agricole. A ce stade, le principal marais salé, le Chott el Djerid, n’est pas une source de contamination. Un modèle aux différences finies a également été élaboré pour simuler les écoulements souterrains et le transport de contaminants dans les oasis. Une calibration pour la période 1950–2000 a été effectuée afin d’ajuster les paramètres des systèmes géologique et chimique. La simulation de projets planifiés d’extraction prédit une aggravation de la situation actuelle. La conservation du régime d’extraction actuel ne réduira ou ne stoppera pas l’augmentation de la salinité.

---

Resumen El agua subterránea que se bombea del acuífero semi-confinado Terminal Complejo (CT) es un importante factor de producción en la agricultura de riego con oasis en el sur de Túnez. Se ha observado un incremento en la salinidad del agua subterránea como consecuencia de la abstracción creciente del acuífero durante las últimas décadas. Se investigaron todas las fuentes de contaminación usando datos hidroquímicos disponibles de la década de 1980s. Las muestras de agua se tomaron de drenajes y pozos de observación alojados tanto en el CT como los acuíferos freáticos y se analizaron en relación con química, temperatura, isótopos y otros trazadores ambientales. Se sugieren mecanismos de salinización local, i.e. el ascenso de agua salada a partir del acuífero confinado subyacente Intercalado Continental (CI), así como también retorno de flujo de agua de drenaje agrícola. En esta etapa, el pan salado principal, el Chott el Djerid, no es una fuente de contaminación. También se desarrolló un modelo de diferencia finita para simular el flujo de agua subterránea y el transporte de contaminantes en los oasis. Se calibró el modelo para el periodo 1950–2000 para de este modo ajustar parámetros del sistema químico y geológico. La simulación de los proyectos de extracción que se han planeado predice un empeoramiento de la situación actual. El mantenimiento del régimen de extracción actual no reducirá o detendrá el incremento de salinidad.

     

Simulating the impacts of groundwater pumping on stream–aquifer dynamics in semiarid northwestern Oklahoma,USA

Visual MODFLOW, a numerical groundwater flow model, was used to evaluate the impacts of groundwater exploitation on streamflow depletion in the Alluvium and Terrace aquifer of the Beaver-North Canadian River (BNCR) in northwestern Oklahoma, USA. Water demand in semi-arid northwestern Oklahoma is projected to increase by 53% during the next five decades, driven primarily by irrigation, public water supply, and agricultural demand. Using MODFLOW’s streamflow routing package, pumping-induced changes in baseflow and stream leakage were analyzed to estimate streamflow depletion in the BNCR system. Simulation results indicate groundwater pumping has reduced baseflow to streams by approximately 29% and has also increased stream leakage into the aquifer by 18% for a net streamflow loss of 47%. The magnitude and intensity of streamflow depletion, however, varies for different stream segments, ranging from 0 to 20,804 m³/d. The method provides a framework for isolating and quantifying impacts of aquifer pumping on stream function in semiarid alluvial environments.

---

Résumé  Visual MODFLOW, un modèle numérique d’écoulement des eaux souterraines, a été utilisé pour évaluer les impacts de l’exploitation de l’eau souterraine sur la réduction du régime de la rivière dans un aquifère d’alluvions et de terrasse de la Rivière canadienne Beaver-North (BNCR) au nord-ouest de l’Oklahoma, aux USA. La demande d’eau dans le nord-ouest semi-aride de l’Oklahoma devrait augmenter de 53% durant les cinq prochaines décennies, essentiellement du fait de l’irrigation, la demande en eau potable et la demande agricole. En utilisant le module numérique traitant les rivières dans MODFLOW, les changements sur le régime de base des cours d’eau induits par le pompage et la recharge par la rivière ont été analysés pour estimer la réduction du régime de la rivière du système BNCR. Les résultats de la simulation indiquent que le pompage de l’eau souterraine réduit le débit de base de la rivière de 18% pour une perte nette du débit de la rivière de 47%. L’importance et l’intensité de la réduction du débit de base, toutefois, varie selon le bief, de 0 à 20,840 m³/d. La méthode apporte un cadre pour isoler et quantifier les impacts sur le pompage de l’aquifère sur le fonctionnement de la rivière dans un environnement alluvial semi-aride.

---

Resumen  El Visual MODFLOW, que es un modelo numérico de flujo de agua subterránea, fue usado para evaluar los impactos de la explotación del agua subterránea, en la disminución de la cantidad de flujo de una corriente superficial, relacionada al acuífero del aluvión y de la terraza del Río canadiense Beaver-North (BNCR) en Oklahoma del noroeste, EUA. Se proyecta que la demanda de agua aumentará en un 53%, en el sector semiárido de Oklahoma del noroeste, durante las próximas cinco décadas, principalmente debido a la irrigación, suministro de agua público, y la demanda agrícola. Mediante el uso del paquete especial para relaciones corriente superficial- acuífero (SFR) de MODFLOW, se analizaron cambios inducidos por el bombeo, tanto en el flujo base, como en el goteo del río, para estimar la disminución del flujo de la corriente en el sistema de BNCR. Los resultados de la simulación indican, que el bombeo de agua subterránea ha reducido el flujo base a los arroyos aproximadamente en un 29% y también ha aumentado el goteo del arroyo hacia el acuífero en 18%, resultando una pérdida neta del flujo del río de un 47%. La magnitud e intensidad de la disminución del flujo del río, sin embargo, varía en sus diferentes segmentos, oscilando de 0 a 20,804 m³/d. El método suministra un marco, para aislar y cuantificar los impactos del bombeo del acuífero, sobre el funcionamiento de una corriente superficial en los ambientes aluviales semiáridos.

     

Analysis of radial movement of an unconfined leaky aquifer due to well pumping and injection

Radial movement of an unconfined leaky aquifer was studied with respect to hydraulic forces that are induced by well recharge and discharge. New analytic solutions in the velocity and displacement fields were found and applied to describe transient movement in an unconfined leaky aquifer. Linear momentum and mass balance of saturated porous sediments, the Darcy-Gersevanov law, and the analytic solution of hydraulic drawdown for unsteady flow within the unconfined leaky aquifer were introduced to find the new solutions. Analytic results indicate that the nonlinear relation between the initial hydraulic head (h₀) and the well function has an insignificant effect on the aquifer transient movement when the drawdown s<0.02h ₀. When the well function is simplified with different assumptions and pumping conditions, the new solutions correspondingly reduce to cases that are similar to the Hantush-Jacob, Muskat, and Theis transient movement of a confined leaky aquifer. It was found that large leakance is important in slowing radial movement and reducing aquifer deformation. Flow velocity in the aquifer is more responsive to leakance than to cumulative displacement within the aquifer. The zones and boundary with tensile stress can be located using the same approach applied to a confined aquifer for risk assessment of earth fissuring.

---

Résumé Le mouvement radial dans un aquifère semi-perméable libre a été étudié en considérant les forces hydrauliques induites par la recharge et le prélèvement au niveau d’un puits. De nouvelles solutions analytiques des champs de vitesse et de déplacement ont été trouvées et appliquées pour décrire le mouvement transitoire dans un aquifère semi-perméable libre. Pour trouver les nouvelles solutions, le moment linéaire et le bilan de masse dans les sédiments poreux saturés, la loi de Darcy-Gersevanov, et la solution analytique du rabattement hydraulique pour un écoulement transitoire dans un aquifère semi-perméable libre sont pris en compte. Les résultats analytiques montrent que la relation non-linéaire entre le potentiel hydraulique initial (h₀) et la fonction de puits a un effet insignifiant sur le mouvement transitoire de l’aquifère quand le rabattements s <0.02h₀. Si la fonction de puits est simplifiée en tenant compte de différentes hypothèses et conditions de pompage, les nouvelles solutions se réduisent en conséquence à des cas similaires au mouvement transitoire dans un aquifère semi-perméable captif selon Hantush-Jacob, Muskat, et Theis. Il s’est avéré qu’une drainance élevée est importante pour le ralentissement du mouvement radial et la réduction de la déformation de l’aquifère. La vitesse d’écoulement dans l’aquifère est plus sensible à la drainance qu’au déplacement cumulé dans l’aquifère. Les zones et la frontière présentant des contraintes de traction peuvent être localisées en utilisant la même approche appliquée à un aquifère captif pour l’évaluation des risques de fissuration du sol.

---

Resumen Se ha estudiado el movimiento radial en un acuífero no confinado con fuga en relación a las fuerzas hidráulicas que son inducidas por la descarga y recarga del pozo. Se encontraron nuevas soluciones analíticas en los campos de desplazamiento y velocidad, las cuales se aplicaron para describir el movimiento transitorio en un acuífero no confinado con fuga. Se han introducido cálculos de balance de masa y momentum lineal de los sedimentos porosos saturados, la Ley Darcy-Gersevanov, y la solución analítica de descenso hidráulico para flujo sin régimen permanente dentro del acuífero no confinado con goteo para encontrar las nuevas soluciones. Los resultados analíticos indican que la relación no lineal entre la presión hidráulica inicial (h₀) y la función del pozo tiene efecto insignificante en el movimiento transitorio del acuífero cuando el descenso s < 0.02h₀. Cuando se simplifica la función del pozo en base a diferentes supuestos y condiciones de bombeo se reducen las nuevas soluciones a casos que son similares a el movimiento transitorio Hantush-Jacob, Muskat, y Theis de un acuífero confinado con goteo. Se encontró que las fugas grandes son importantes en la disminución del movimiento radial y en la reducción de la deformación del acuífero. La velocidad de flujo en el acuífero responde más a las fugas que al desplazamiento acumulativo dentro del acuífero. Pueden localizarse zonas y límites con esfuerzos tensionales usando el mismo enfoque aplicado a acuíferos confinados para evaluación de riesgo de generación de fracturas en el terreno.

---

     

Hydrogeology of the Winnipeg Formation in Manitoba, Canada

The Winnipeg Formation is the basal sedimentary unit throughout much of southern and central Manitoba, Canada, where it forms a regional aquifer over most of its extent. This aquifer is an important source of water in southeastern Manitoba and in Manitoba’s Interlake area, but in most other areas, groundwater within the aquifer is saline. Chemical and isotopic evidence indicate the presence of groundwaters of three different origins: (1) basin brines; (2) modern meteoric recharge; and (3) subglacial recharge that occurred during the late Pleistocene. Hydraulic head and sedimentary facies distributions indicate that the flow system in parts of the area is not in a state of equilibrium and saline waters will encroach on areas currently occupied by freshwater in some areas, while in other areas, freshwater will replace saline water. These features must be considered in groundwater resource management, as groundwater withdrawals will likely hasten these processes.

---

Resumen La Formación Winnipeg es la unidad sedimentaria basal en la mayor parte de Manitoba central, Canadá, donde forma un acuífero regional en la mayor parte de su extensión. Este acuífero es una fuente importante de agua en el Sureste de Manitoba y el área de entrelagos de Manitoba, pero en la mayoría de las otras zonas del acuífero, el agua es salina. Las evidencias químicas e isotópicas indican que existen aguas subterráneas de tres orígenes diferentes: (1) salmueras de cuenca; (2) recarga meteórica actual; y (3) recarga subglacial ocurrida durante el Pleistoceno Superior. Los niveles piezométricos y la distribución de las facies sedimentarias indican que el sistema de flujo no se encuentra en estado de equilibrio en parte del área y las aguas salinas irán invadiendo áreas actualmente ocupadas con aguas dulces, mientras que en otras zonas el agua dulce está reemplazando al agua salina. Estos hechos deben ser considerados en la gestión de las aguas subterráneas como recurso, ya que las extracciones de agua acelerarán probablemente estos procesos.

---

Résumé La Formation de Winnipeg est l’unité sédimentaire de base sur la plus grande partie du Sud et du centre du Manitoba au Canada, où elle forme un aquifère régional sur pratiquement toute son extension. Cet aquifère représente une ressource en eau importante dans le Sud-Est du Manitoba et dans les zones d’entre les lacs, mais salée dans la plus part des autres zones. Les indications isotopiques et chimiques permettent de distinguer trois différentes origines des eaux souterraines: (1) les saumures de bassin; (2) la recharge météoritique moderne; (3) la recharge sub-glaciaire qui est apparue durant le Pléistocène récent. Les charges hydrauliques et la distribution des faciès sédimentaires indiquent que le système d’écoulement dans certaines zones n’est pas dans un état d’équilibre et que les eaux salées empièteront sur des zones d’eau douce, tandis que dans d’autres zones l’eau douce remplacera les eaux salées. Ces aspects doivent être considérés dans la gestion des ressources en eau souterraine, car le prélèvement des eaux souterraines pourrait accentuer ces processus.

     

Evaluation of the ability of an artificial neural network model to simulate the input-output responses of a large karstic aquifer: the La Rochefoucauld aquifer (Charente,France)

The ability of artificial neural networks (ANN) to model the rainfall-discharge relationships of karstic aquifers has been studied in the La Rochefoucauld karst system, south-west France, which supplies water to the city of Angoulême. A neural networks model was developed based on MLP (multi-layer perceptron) networks and the Levenberg-Marquardt optimization algorithm. Raw rainfall data were used without transformation into effective rainfall. This allowed for the elimination of certain non-verifiable simplifying assumptions and their subsequent introduction into the modeling procedure. The raw rainfall and discharge data were divided into three groups for the training, the validation and the prediction test of the ANN model. The training and validation phases led to a very satisfactory calibration of the ANN model. The attempt to predict discharges showed that the ANN model is able to simulate the karstic aquifer discharges. The shape of the simulated hydrographs was found to be similar to that of the actual hydrographs. These encouraging results make it possible to consider interesting and new prospects for the modeling of karstic aquifers, which are highly non-linear systems.

---

Resumen Se ha estudiado la capacidad de las redes artificiales neurales (ANN) para modelizar las relaciones de lluvia-descarga de acuíferos kársticos en el sistema kárstico La Rochefocauld, al suroeste de Francia, el cual abastece de agua a la ciudad de Angoulême. Se desarrolló un modelo de redes neurales en base a redes MLP (Perceptron Multi-Capas) y el algoritmo de optimización Levenberg-Marquardt. Se utilizaron datos de lluvia sin la transformación hacia lluvia efectiva. Esto permitió la eliminación de ciertos supuestos simplificadores no verificables y su subsiguiente introducción en el procedimiento de modelizado. Los datos brutos de descarga y lluvia se dividieron en 3 grupos para la preparación, validación y la prueba de predicción del modelo ANN. Las fases de preparación y validación llevaron a una calibración muy satisfactoria del modelo ANN. El intento por predecir descargas mostró que el modelo ANN es capaz de simular las descargas del acuífero kárstico. Se encontró que la forma de los hidrogramas sintéticos es similar a la de los hidrogramas reales. Estos resultados alentadores hacen posible considerar prospectos nuevos e interesantes para el modelizado de acuíferos kársticos los cuales son sistemas altamente no-lineares.

---

Résumé L’aptitude des réseaux de neurones artificiels (RNA) à modéliser les relations pluie-débit des aquifères karstiques a été évaluée sur le karst de La Rochefoucauld (Sud-Ouest de la France), qui fournit l’eau potable à la capitale régionale Angoulême. Un modèle RNA a été développé à cet effet, basé sur les réseaux PMC (Perceptron Multicouche) et l’algorithme d’optimisation de Levenberg-Marquardt. Les données de pluie utilisées concernent la pluie brute, sans transformation en pluie efficace, ce qui permet de s’affranchir de certaines hypothèses simplificatrices non vérifiables. Les données de pluie brute et de débit ont été réparties en 3 groupes pour l’apprentissage, la validation et le test de prédiction du RNA. Les phases d’apprentissage et de validation ont permis d’aboutir à une calibration très satisfaisante du modèle RNA. La tentative de prédiction a montré que le RNA est apte à simuler les débits de l’aquifère karstique à partir de la pluie brute. La forme des hydrogrammes simulés est semblable à celle des hydrogrammes réels. Les résultats obtenus sont très encourageants et permettent d’envisager des perspectives intéressantes et nouvelles de modélisation des aquifères karstiques, qui sont des systèmes hautement non-linéaires.

     

Remote sensing and GIS for mapping groundwater recharge and discharge areas in salinity prone catchments,southeastern Australia

Identifying groundwater recharge and discharge areas across catchments is critical for implementing effective strategies for salinity mitigation, surface-water and groundwater resource management, and ecosystem protection. In this study, a synergistic approach has been developed, which applies a combination of remote sensing and geographic information system (GIS) techniques to map groundwater recharge and discharge areas. This approach is applied to an unconfined basalt aquifer, in a salinity and drought prone region of southeastern Australia. The basalt aquifer covers ~11,500 km² in an agriculturally intensive region. A review of local hydrogeological processes allowed a series of surface and subsurface indicators of groundwater recharge and discharge areas to be established. Various remote sensing and GIS techniques were then used to map these surface indicators including: terrain analysis, monitoring of vegetation activity, and mapping of infiltration capacity. All regions where groundwater is not discharging to the surface were considered potential recharge areas. This approach, applied systematically across a catchment, provides a framework for mapping recharge and discharge areas. A key component in assigning surface and subsurface indicators is the relevance to the dominant recharge and discharge processes occurring and the use of appropriate remote sensing and GIS techniques with the capacity to identify these processes.

---

Résumé Sur un bassin versant, l’identification des zones d’alimentation et des aires d’émergence des eaux souterraines est capitale pour mettre en œuvre des stratégies de réduction de la salinité, de gestion des eaux souterraines et superficielles, et de protection des écosystèmes. Une approche synergique a été développée, sur la base de recherches existantes pour cartographier de manière systématique les zones de recharge et d'émergence. Elle combine des techniques de télédétection et de Systèmes d’Information Géographique, dans le but de cartographier les zones d’alimentation et d’émergence. Ces techniques ont été appliquées à un aquifère basaltique libre, dans un secteur sujet à la sécheresse et à la salinité du Sud-Est australien. L’aquifère considéré couvre environ 11,500 km² d’une région essentiellement agricole. Un inventaire des processus hydrogéologiques locaux a conduit à l’établissement d’une série d’indicateurs superficiels et subsuperficiels des zones d’alimentation et d’émergence. Plusieurs méthodes de télédétection et de SIG furent alors utilisées pour cartographier ces indicateurs de surface, y compris les analyses de terrain, le suivi de l’activité de la végétation et la cartographie des capacités d’infiltration. Toutes les régions en dehors des zones d'émergence sont considérées commes des zones d'alimentation potentielles. Un des composants clés de l’attribution des indicateurs superficiels et subsuperficiels est la concordance avec les processus majeurs d’alimentation et de prélèvement entrant en jeu, et l’utilisation de méthodes de télédétection et SIG adaptées permettant d’identifier ces processus.

---

Resumen La identificación de áreas de descarga y recarga de agua subterránea es crítica para la implementación de estrategias efectivas relacionadas con mitigación de salinidad, gestión de recursos de agua subterránea y agua superficial, y protección de ecosistemas. En este estudio se ha desarrollado un enfoque sinérgico apoyándose en investigación existente el cual aplica una combinación de técnicas de sistemas de información geográfico (SIG) y sensores remotos para cartografiar áreas de recarga y descarga de agua subterránea. Este enfoque se aplica a un acuífero de basalto no confinado en una región del sureste de Australia propensa a sequía y salinidad. El acuífero de basalto cubre aproximadamente 11,500 km² en una región con agricultura intensiva. Una revisión de los procesos hidrogeológicos locales permitió el establecimiento de una serie de indicadores superficiales y subterráneos de áreas de recarga y descarga de agua subterránea. Luego se utilizaron varias técnicas de SIG y sensores remotos para cartografiar estos indicadores superficiales incluyendo: análisis del terreno, monitoreo de la actividad de la vegetación, y cartografiado de la capacidad de infiltración. Todas las regiones donde el agua subterránea no descarga en la superficie fueron consideradas áreas potenciales de recarga. Este enfoque que se ha aplicado sistemáticamente a través de la cuenca aporta un marco para la cartografía de áreas de recarga y descarga. Un componente clave en la asignación de indicadores superficiales y subterráneos lo constituye la relevancia a los procesos dominantes de recarga y descarga que ocurren y el uso de técnicas apropiadas de SIG y sensores remotos con la capacidad para identificar esos procesos.

     

Aquifer storage capacity and maximum annual yield from long-term aquifer fluxes

Long-term time series data of aquifer recharge, groundwater extraction, and discharge are used to estimate aquifer storage capacity and maximum annual yield. Aquifer storage capacity is defined as the maximum volume of water that can be stored in an aquifer. It is estimated using a transient water-balance approach. The maximum annual yield is defined as the maximum combined groundwater extraction plus discharge that can be sustained in an aquifer judged by the historical record of recharge. It is determined according to a graphical mass-curve method. These two quantities are useful in aquifer characterization and groundwater management, the apportionment of groundwater rights and aquifer storage and recovery operations being two frequent applications. Time series data from the Edwards Aquifer, Texas, USA, illustrate the application of the methods presented.

---

Résumé Des données sur le long terme de la recharge d’aquifère, de l’exploitation de l’eau souterraine et de la vidange sont utilisées pour estimer la capacité d’emmagasinement de l’aquifère et la capacité annuelle maximum. La capacité d’emmagasinement de l’aquifère est définit comme le volume maximum de l’eau qui peut être emmagasinée dans le réservoir. Elle est estimée en utilisant une approche par bilan hydrologique en transitoire. La capacité annuelle maximum est définit comme la combinaison de l’extraction de l’eau souterraine et de la vidange, combinaison qui peut être jugée durable ou non au regard de l’examen des valeurs de recharge. Elle est déterminée au moyen d’une méthode de courbe des valeurs cumulées. Ces deux quantités sont très utiles pour caractériser les aquifères et pour gérer les eaux souterraines, l’imputation des droits aux eaux souterraines, et du fait que les opérations de recharge des aquifères et de récupération sont très fréquentes. Les données temporelles provenant de l’aquifère Edwards, Texas, USA, illustrent l’application de la méthode présentée.

---

Resumen Para estimar la capacidad de almacenamiento de un acuífero y su rendimiento máximo anual, se han utilizado series temporales largas de datos de recarga al acuífero, de extracciones de agua subterránea y de descarga. La capacidad de almacenamiento de un acuífero se define como el volumen máximo de agua que puede ser almacenado en el mismo. Se estima utilizando una aproximación con un balance de agua en régimen transitorio. El rendimiento máximo anual se define como el máximo resultante de la combinación de la extracción de aguas subterráneas y la descarga que puede ser sostenido por un acuífero según los registros históricos de recarga y se define según un método gráfico masa-curva. Estos dos valores son útiles en la caracterización de los acuíferos y en la gestión de las aguas subterráneas, el reparto de los derechos sobre las aguas subterráneas y las operaciones de almacenaje y recuperación de los acuíferos, que son dos aplicaciones frecuentes. Las series de tiempo para el Acuífero Edwards, Texas, USA, ilustra la aplicación de los métodos presentados.

     

Comparative analysis between analytical approximations and numerical solutions describing recession flow in unconfined hillslope aquifers

Recession flow of aquifers from a hillslope can be described by the non-linear Boussinesq equation. Under strong assumptions and for specific conceptual formulations, different authors derived analytical approximations or linearized versions to this partial differential equation. A comparative analysis between some analytical approximations of the Boussinesq equation and the numerical solution of the recession flow of an unconfined homogeneous aquifer (horizontal, inclined and concave aquifer floor) was carried out. The objective was to define the range where the analytical solutions approximate the numerical solution. The latter was considered in this study as the reference method, because it requires fewer assumptions. From the considered analytical approximations, exponential decay relationships were found to be mainly valid for fine domain materials when horizontal, mild slopes (less than 2%) and concave aquifer floors were considered, but failed to reproduce coarse aquifer numerical model outflows, in contrast to the quadratic decay relationship, which better reproduce outflows in such domains. On the basis of the comparative analysis, it has been found that recession flows obtained with the considered analytical approximations yield similar values only for certain ranges of aquifer properties and geometries.

---

Résumé L’écoulement de décrue dans les aquifères de versant peut être décrit par l’équation non-linéaire de Boussinesq. En considérant d’importantes hypothèses et des formulations conceptuelles spécifiques, différents auteurs ont dérivé des approximations analytiques ou des versions linéarisées de cette équation partielle différentielle. Une analyse comparative entre des approximations analytiques de l’équation de Boussinesq et une solution numérique de la décrue d’un aquifère homogène libre (limite inférieure de l’aquifère horizontale, inclinée et concave) a été effectuée. L’objectif était de définir la plage pour laquelle les solutions analytiques approchent la solution numérique. Cette dernière a été considérée dans cette étude comme la méthode de référence car elle demande moins d’hypothèses. Parmi les approximations analytiques considérées, les relations de décroissance exponentielle se sont avérées être principalement valides pour des matériaux fins lorsque des couches aquifères horizontales, concaves et de pentes douces (moins de 2%) étaient considérées; en revanche, celles-ci n’ont pas réussi à reproduire les écoulements sortants d’un modèle numérique pour un aquifère à matériaux grossiers, contrairement à la relation de décroissance quadratique qui reproduit mieux les écoulements sortants dans de tels milieux. Sur la base de l’analyse comparative, il s’est avéré que les écoulements en décrue obtenus avec les approximations analytiques considérées présentent des valeurs similaires seulement pour certains critères géométriques et propriétés de l’aquifère.

---

Resumen Puede describirse el flujo de recesión de acuíferos de una ladera mediante la ecuación no-lineal de Boussinesq. Bajo fuertes supuestos y para formulaciones conceptuales específicas, distintos autores han derivado aproximaciones analíticas o versiones lineales para esta ecuación parcial diferencial. Se realizó un análisis comparativo entre algunas aproximaciones analíticas de la ecuación de Boussinesq y la solución numérica del flujo de recesión de un acuífero homogéneo no confinado (piso de acuífero cóncavo, inclinado, y horizontal). El objetivo consistió en definir el rango donde las soluciones analíticas aproximan la solución numérica. En este estudio se consideró la solución numérica como el método de referencia debido a que requiere menores supuestos. De las aproximaciones analíticas consideradas se encontró que las relaciones de desintegración exponencial eran principalmente válidas para materiales de ámbitos finos donde se consideraron pisos de acuíferos horizontal, pendientes medianas (menores a 2%) y cóncavos pero fallaron en reproducir flujos de salida de modelos numéricos de ámbito grueso, en contraste con la relación de desintegración cuadrática la cual reproduce mejor los flujos de salida en tales ámbitos. En base al análisis comparativo se encontró que los flujos de recesión obtenidos con las aproximaciones analíticas consideradas aportaron valores similares solo para ciertos rangos de propiedades y geometrías de acuíferos.

     

Sources of dissolved solids and water in Wadi Al Bih aquifer,Ras Al Khaimah Emirate,United Arab Emirates

Regional brines that underlie the potable groundwater appear to be responsible for the increase in dissolved solids in the Wadi Al Bih aquifer in the Ras Al Khaimah Emirate, United Arab Emirates. In this karstic carbonate aquifer, groundwater extraction exceeds recharge and the reduced heads can induce transport of underlying brines into the potable water aquifer. Increasing dissolved solids with time threatens the continued use of groundwater for agricultural and domestic uses. The potential of intrusion of seawater, dissolution of minerals, or intrusion of regional brines as a source of these solutes were evaluated based on groundwater samples collected in April and September 1996 from the Wadi Al Bih well field and isotope data from previously collected samples. Hydrogeologic conditions and solute modeling suggest the dominant source of solutes in the Wadi Al Bih aquifer is most likely an underlying regional brine. The pervasive presence of tritium is consistent with recent (since 1960s) recharge in the drainage basin, thus, the solutes and water appear to be from different sources. The chemical and isotopic data are also consistent with enhanced groundwater recharge associated with dams constructed to reduce flooding.

---

Résumé Dans l’aquifère de Wadi Al Bih, dans l’Emirat de Ras Al Khaimah (Emirats Arabes Unis), les eaux sursalées présentes sous les eaux souterraines potables paraissent responsables de l’augmentation en solides dissous. Dans cet aquifère carbonaté karstique, l’exploitation des eaux souterraines est supérieure à la recharge, et la baisse des niveaux piézométriques peut engendrer une remontée des eau sursalées vers la ressource en eau potable. L’augmentation des solides dissous dans le temps met en danger la pérennité de l’utilisation des eaux souterraines pour l’agriculture et les besoins domestiques. Les potentiels de l’intrusion d’eau de mer, de la dissolution des minéraux et de l’intrusion d’eaux sursalées régionales en tant que sources de ces solutés ont été évalués sur la base d’échantillons prélevés en avril 1996 sur le champ captant de Wadi Al Bih, et de données isotopiques acquises antérieurement. Les conditions hydrogéologiques et les modélisations de solutés suggèrent que la source majeure des solutés la plus vraisemblable est une eau sursalée sous-jacente à l’échelle régionale, dans le cas de l’aquifère de Wadi Al Bih. L’omniprésence du tritium concorde avec une réalimentation récente (depuis les années 1960) sur le bassin versant, c’est pourquoi l’eau et les solutés paraissent provenir de sources différentes. Les données chimiques et isotopiques concordent également avec une réalimentation accrue des eaux souterraines, liée à la construction récente de digues pour réduire les inondations.

---

Resumen El incremento de sólidos disueltos en el acuífero Wadi Al Bih, en el Emirato Ras Al Khaimah, Emiratos árabes Unidos, parece deberse a la existencia de salmueras regionales, subyacentes a las aguas subterráneas potables. En este acuífero carbonatado kárstico, la extracción de agua subterránea supera la recarga y el descenso de niveles puede provocar el transporte de las salmueras hacia las aguas potables del acuífero. El incremento temporal en sólidos disueltos amenaza el uso continuo de agua subterránea para usos agrícolas y domésticos. Se ha evaluado el potencial de intrusión marina, disolución de minerales o intrusión de salmueras regionales como una fuente de estos solutos, a partir de muestras de agua subterránea recogidas en abril y septiembre de 1996 en el campo de pozos de Wadi Al Bih y de datos isotópicos de muestreos anteriores. Las condiciones hidrogeológicas y el modelo de solutos sugiere que la fuente dominante de solutos en el acuífero Wadi Al Bih son probablemente las salmueras subyacentes. Los datos químicos e isotópicos son también consistentes con un aumento de la recarga de agua subterránea asociado con la construcción de presas para reducir inundaciones.

     

Use of geochemical modeling to evaluate the hydraulic connection of aquifers: a case study from Chianan Plain, Taiwan

Aquifers are generally composed of highly permeable layers that can conduct a considerable amount of groundwater. Traditionally, aquifer units are correlated through the concept of lithostratigraphy. For low-permeable aquifers, it is difficult to define the spatial distribution of hydrogeological units, and this study attempts to use geochemical modeling to identify the groundwater flow paths in an area of Taiwan. Multiple geochemical analyses, including groundwater chemistry; stable isotopic compositions of hydrogen, oxygen and carbon; and radiocarbon contents were performed. Using these parameters as the constraints of geochemical models, the hydraulic connection was examined between pairs of possibly interlinked wells along four selected cross sections, and the conceptual groundwater model was accordingly established. The resultant model suggests that the hydraulic connection between aquifers should be correlated with the concept of chronological stratigraphy, especially for low-permeable, unconsolidated aquifers. Using Darcy’s law, the hydraulic conductivities of the fine-sand aquifers were estimated to be between 3.14 × 10⁻⁵ and 1.80 × 10⁻⁴ m/s, which are roughly one order of magnitude higher than those derived by in situ pumping tests. The substantial extraction of groundwater over a long period in the studied area could accelerate groundwater flow, leading to an overestimation of the aquifer permeability.

---

Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.

---

Résumé En règle générale, les aquifères sont constitués de couches très perméables capables de conduire des quantités d’eau considérables. Les unités aquifères sont traditionnellement corrélées d’un point de vue lithostratigraphique. Dans le cas des aquifères à faible perméabilité, il s’avère difficile de définir une distribution spatiale des unités hydrogéologiques, et cette étude tente d’utiliser la modélisation géochimique pour identifier les cheminements des eaux souterraines dans un secteur de Taiwan. Plusieurs analyses géochimiques ont été réalisées ; elles incluent la chimie des eaux souterraines, les isotopes stables de l’hydrogène, de l’oxygène et du carbone et l’abondance en carbone 14. En utilisant ces paramètres comme contraintes pour les modèles géochimiques, les liaisons hydrauliques entre les puits potentiellement interconnectés ont été étudiées deux par deux, le long de quatre coupes présélectionnées, et le modèle conceptuel a été établi en conséquence. Ce modèle résultant suggère que les connexions hydrauliques inter-aquifères devraient être corrélées avec les concepts stratigraphiques, surtout pour les aquifères non consolidés à faible perméabilité. Les perméabilités dans les aquifères de sables fins ont été estimées par la loi de Darcy entre 3.14 × 10⁻⁵ et 1.80 × 10⁻⁴ m/s, soit environ un ordre de grandeur au-dessus de celles issues des tests de pompage sur site. L’exploitation substantielle des eaux souterraines sur une longue période dans la zone d’étude peut accélérer les écoulements souterrains, menant à une surestimation de la perméabilité de l’aquifère.

---

Resumen Los acuíferos están generalmente compuestos de capas muy permeables que pueden conducir una cantidad considerable de agua subterránea. Tradicionalmente, se ponen en correlación las unidades acuíferas a través del concepto de litoestratigrafía. Para los acuíferos de baja-permeabilidad, es difícil de definir la distribución espacial de unidades hidrogeológicas y este estudio intenta usar modelamiento geoquímico para identificar las direcciones del flujo de agua subterránea en un área de Taiwán. Se realizó el análisis geoquímico múltiple, incluyendo la química del agua subterránea, la composición de isótopos estables de hidrógeno, oxígeno y carbono, y el contenido de radiocarbono. Usando estos parámetros como limitantes de modelos geoquímicos, la conexión hidráulica se examinó entre pares de pozos posiblemente interrelacionados, a lo largo de cuatro cortes transversales seleccionados y de acuerdo con esto se estableció el modelo conceptual del agua subterránea. El modelo resultante sugiere que la conexión hidráulica entre los acuíferos deba interrelacionarse con el concepto de estratigrafía cronológica, sobre todo para los acuíferos sin consolidar de baja permeabilidad. Usando la ley de Darcy, se estimaron las conductividades hidráulicas de los acuíferos de arena fina obteniendo valores entre 3.14 × 10⁻⁵ y 1.80 × 10⁻⁴ m/s, las cuales son aproximadamente un orden de magnitud mayor, que aquéllos obtenidos in situ por las pruebas de bombeo. La extracción considerable de agua subterránea durante un periodo largo en el área estudiada, podría acelerar el flujo de agua subterránea, llevando a una sobrestimación de la permeabilidad en los acuíferos.

     

A simple daily soil–water balance model for estimating the spatial and temporal distribution of groundwater recharge in temperate humid areas

Quantifying the spatial and temporal distribution of natural groundwater recharge is usually a prerequisite for effective groundwater modeling and management. As flow models become increasingly utilized for management decisions, there is an increased need for simple, practical methods to delineate recharge zones and quantify recharge rates. Existing models for estimating recharge distributions are data intensive, require extensive parameterization, and take a significant investment of time in order to establish. The Wisconsin Geological and Natural History Survey (WGNHS) has developed a simple daily soil–water balance (SWB) model that uses readily available soil, land cover, topographic, and climatic data in conjunction with a geographic information system (GIS) to estimate the temporal and spatial distribution of groundwater recharge at the watershed scale for temperate humid areas. To demonstrate the methodology and the applicability and performance of the model, two case studies are presented: one for the forested Trout Lake watershed of north central Wisconsin, USA and the other for the urban-agricultural Pheasant Branch Creek watershed of south central Wisconsin, USA. Overall, the SWB model performs well and presents modelers and planners with a practical tool for providing recharge estimates for modeling and water resource planning purposes in humid areas.

---

Resumen La cuantificación de la distribución espacial y temporal de la recarga natural de agua subterránea es un requisito previo para una modelación y una gestión efectivas de las aguas subterráneas. Dado que se está incrementando el uso de los modelos de flujo para la toma de decisiones de gestión, existe una necesidad creciente de métodos simples y prácticos para delimitar las zonas de recarga y cuantificar los rangos de la misma. Los modelos existentes para la estimación de la distribución de la recarga requieren datos intensivos, una parametrización extensiva y una inversión de tiempo significativa para ser establecidos. El Servicio Geológico y de Historioa Natural de Wisconsin (WGNHS) ha desarrollado un modelo simple de balance diario de agua en el suelo (SWB) que usa de forma sencilla datos disponibles de suelo, de cobertera, topográficos y climáticos conjuntamente con un Sistema de Información Geográfica (GIS) para estimar la distribución espacial y temporal de la recarga de aguas subterráneas a escala de cuenca para zonas templadas húmedas. Para demostrar la metodología, aplicabilidad y el comportamiento del modelo, se presentan dos casos: uno en la cuenca boscosa de Trout Lake, en la zona Norte-Central de Wisconsin, USA y el otro en la Cuenca urbano-agrícola de Pheasant Branch Creek, Sur-Centro de Winconsin, USA. En conjunto, el modelo SWB funciona bien y presenta a los modeladores y planificadores una herramienta práctica para llevar a cabo una estimación de la recarga para propósitos de modelación y planificación de los recursos de agua en zonas húmedas.

---

Résumé Quantifier les distributions temporelle et spatiale de la réalimentation naturelle des eaux souterraines est en règle générale un préalable à une modélisation et une gestion efficaces des eaux souterraines. Etant donné que les modèles numériques sont utilisés de manière croissante dans les prises de décisions en matière de gestion, il existe un besoin accru pour des méthodes simples et pratiques, afin de délimiter les zones d’alimentation et de quantifier les recharges associées. Les modèles existants, destinés à l’estimation de la répartition des réalimentations, demandent énormément de données, un paramétrage long, et un investissement conséquent en temps de mise en œuvre. Le Wisconsin Geologic and Natural History Survey (WGNHS) a développé un modèle simple basé sur un bilan en eau quotidien dans les sols (SWB); il utilise les données directement disponibles sur les sols, l’occupation des sols, la topographie et le climat, en conjonction avec un Système d’Information Géographique, afin d’estimer les distributions temporelle et spatiale de la réalimentation des eaux souterraines à l’échelle du bassin versant, pour les zones humides tempérées. Afin de démontrer la méthodologie, l’applicabilité et les performances du modèle, deux applications sont présentées: la première sur le bassin versant boisé de Trout Lake au centre-nord du Wisconsin (Etats-Unis), et le second sur le bassin versant agricole et urbanisé de Pheasant Branch Creek au centre-sud du Wisconsin. Le modèle SWB se comporte globalement bien, et offre aux modélisateurs un outil fonctionnel pour estimer les réalimentations, dans le cadre de modélisations et de plans de gestion des ressources en eau souterraine dans les zones humides.

     

Development of a GIS-based model for extrapolating mesoscale groundwater recharge estimates using integrated geospatial data sets

A geographic information systems (GIS)-based model was used to derive spatially explicit estimates of recharge using the elevation-dependent chloride mass-balance (EDCMB) approach in a 14-basin area in southwestern Nevada, USA. For each of the 14 basins, a non-linear regression equation was developed relating chloride enrichment to the elevation of spring watersheds. The ratio of chloride enrichment as a function of elevation was used in a GIS framework to transform continuous precipitation data to recharge. Spatial masks that represented two definitions of the lower limits of recharge—one definition based on alluvium/non-alluvium boundaries, the other based on both alluvium/non-alluvium boundaries and an elevation threshold—were applied to each basin. Resultant recharge estimates were then summed for each basin. When compared to summaries of previous recharge estimates for the study area, the EDCMB approach produced results that were within 14 and 3% of two studies, but were significantly greater (31%) than a third. GIS proved to be a very effective tool for combining large spatial data with widely different resolutions into an integrated data set, and also proved to be an efficient mechanism for implementing robust statistical models to estimate recharge.

---

Résumé Un modèle basé sur les systèmes d’informations géographiques (SIG) a été utilisé pour obtenir des estimations explicites de la recharge dans l’espace, en utilisant l’approche du bilan des chlorures en fonction de l’altitude (EDCMB en anglais) dans une zone de 14 bassins, dans le Sud-Ouest du Nevada aux Etats-Unis. Pour chacun des 14 bassins, une régression non-linéaire a été établie, reliant l’enrichissement des chlorures à l’altitude des bassins d’alimentation des sources. Le rapport d’enrichissement des chlorures en fonction de l’altitude a été utilisé dans une approche SIG pour transformer des données pluviométriques continues en terme de recharge. Des masques spatiaux représentant deux définitions des limites de recharge les plus basses ont été appliquées à chaque bassin ; l’une basée sur les limites alluvions/autres formations, l’autre sur ces mêmes limites ainsi que sur un seuil altitudinal. Les recharges estimées résultantes ont ensuite été additionnées pour chaque bassin. Comparée aux précédents bilans de recharge estimés dans la zone d’étude, l’approche EDCMB conduit à des résultats se situant entre 14 et 3% de ceux de deux autres études et largement supérieurs (31%) à ceux d’une troisième étude. Le SIG s’est révélé être un outil très efficace pour combiner dans une base de données intégrée, des données spatiales étendues avec des résolutions largement différentes, ainsi que pour implémenter des modèles statistiques conséquents permettant d’estimer la recharge.

---

Resumen Se utilizó un modelo basado en un Sistema de Información Geográfico (SIG) para derivar espacialmente estimados explícitos de recarga usando el enfoque de balance de masa de cloruro dependiente de elevación (EDCMB) en un área de 14 cuencas en el suroeste de Nevada, Estados Unidos. Para cada una de las 14 cuencas se desarrolló una ecuación de regresión no linear que relaciona el enriquecimiento de cloruro con la elevación de las cuencas con manantiales. La relación de enriquecimiento de cloruro en función de la elevación se usó en un marco SIG para transformar datos continuos de precipitación a recarga. Se aplicaron máscaras espaciales a cada cuenca las cuales representaron dos definiciones de los límites inferiores de recarga, una definición basada en límites de aluvión/sin aluvión, el otro basado en ambos límites aluvión/sin aluvión y una elevación normal de fondo. Los estimados resultantes de recarga fueron sumados para cada cuenca. Al compararse con estimados de recarga previos para el área de estudio se encontró que el enfoque EDCMB produjo resultados que estaban entre el 14 y 3% de dos estudios pero fueron significativamente mayores (31%) que un tercer estudio. El SIG ha demostrado ser una herramienta muy efectiva para combinar datos espaciales grandes con resoluciones muy distintas en un grupo integrado de datos y también ha mostrado ser un mecanismo eficiente para implementar modelos estadísticos robustos para estimar recarga.

     

Solute and isotope constraint of groundwater recharge simulation in an arid environment,Abu Dhabi Emirate,United Arab Emirates

A simple, well-constrained simulation of solute increases in a downgradient direction was used in the shallow unconfined aquifer of eastern Abu Dhabi Emirate (United Arab Emirates). The simulation indicates that the observed exponential increase in solute concentrations results from a combination of upward transport of solutes from underlying mudstones and evaporites, and groundwater losses by evaporation. Groundwater recharge and discharge flux in unconfined regional aquifers in arid regions commonly are difficult to estimate because there are few constraints on the flux of water lost or gained from the system. Total dissolved solids (TDS) and deuterium isotopes (δ²H) in groundwater are used to constrain estimated fluxes to the shallow aquifer of eastern Abu Dhabi Emirate. Vertical upward transport of solutes from underlying mudstones and evaporites accounts for solute increases along approximately the first 80 km of the simulated flow path, but a combination of upward solute transport and evaporation is necessary to explain observed solute concentrations beyond 80 km. Mobilization and transport of solutes in the unsaturated zone by recharging precipitation is not a significant factor.

---

Résumé Une simulation simple et bien contrainte de l’augmentation des solutés dans la direction de l’écoulement a été réalisée sur l’aquifère phréatique et libre de l’Est de l’émirat d’Abu Dhabi (émirats Arabes Unis). La simulation indique que l’augmentation observée de solutés, de forme exponentielle, est le produit de la combinaison du transport ascendant de solutés provenant des mudstones et des évaporites sous-jacents, et de l’évaporation des eaux souterraines. La recharge des eaux souterraines et les flux de vidange dans les aquifères régionaux libres des régions arides sont habituellement difficiles à estimer du fait du peu de contraintes agissant sur la perte ou le gain d’eau dans le système. La charge dissoute totale (TDS, en anglais) et les isotopes deutérium (δ²H) des eaux souterraines sont utilisés pour contraindre l’estimation des flux en direction de l’aquifère phréatique de l’Est de l’émirat d’Abu Dhabi. Le transport vertical ascendant des solutés à partir des mudstones et des évaporites fait parti des processus qui concentrent les solutés approximativement le long des 80 premiers kilomètres du trajet des écoulements simulés, mais une combinaison du transport ascendant des solutés et de l’évaporation est nécessaire pour expliquer les concentrations observées en soluté au delà des 80 km. La mobilisation et le transport des solutés par les précipitations efficaces dans la zone non saturée ne représentent pas des facteurs significatifs.

---

Resumen Se utilizó una simulación simple y muy restringida, de incrementos de soluto en la dirección del flujo, en el acuífero libre poco profundo del Emirato oriental de Abu Dhabi, (Emiratos árabes Unidos). La simulación indica que el aumento exponencial observado en las concentraciones del soluto, resulta de una combinación de transporte ascendente de solutos desde las evaporitas y lodolitas subyacentes, y por pérdidas del agua subterránea por evaporación. Los flujos de recarga y descarga de agua subterránea, en los acuíferos libres regionales en regiones áridas, normalmente son difíciles estimar, porque hay algunas restricciones en el flujo de agua perdido o ganado por el sistema. Los sólidos disueltos totales (TDS) y los isótopos de deuterio (δ²H) en el agua subterránea, se usan para forzar los flujos estimados hacia el acuífero poco profundo del Emirato oriental de Abu Dhabi. El transporte ascendente vertical de solutos desde las lodolitas y evaporitas subyacentes, involucra los aumentos del soluto a lo largo de aproximadamente los primeros 80 Km. de la dirección de flujo simulada, pero es necesaria una combinación de transporte ascendente del soluto y la evaporación para explicar las concentraciones del soluto observadas más allá de 80 Km. La movilización y transporte de solutos en la zona no saturada, por recarga debida a precipitación, no es un factor significante.

     

Groundwater flow modelling of the regional aquifer of the Pampa del Tamarugal,northern Chile

The Pampa del Tamarugal Aquifer (PTA) is an important source of groundwater in northern Chile. In this study, a groundwater flow model of this aquifer is developed and calibrated for the period 1983–2004. The model reproduces the observed flow-field and the water balance components reasonably well. Five scenarios are defined to evaluate the response to different pumping situations. These scenarios show that groundwater heads will continue to decrease with the present pumping discharge rates. To account for variations in the model results due to uncertainties in average recharge rates, randomly generated recharge realizations with different levels of uncertainty are simulated. Evaporation flow rates and groundwater flowing out of the modelled area seem unaffected by the recharge uncertainty, whereas the storage terms can vary considerably. For the most intensive pumping scenario under the generated random recharge rates, it is unlikely that the cumulative discharged volume from the aquifer, at the end of the simulation period, will be larger than 12% of the estimated groundwater reserve. Fluctuations in simulated groundwater heads due to uncertainties in the average recharge values are more noticeable in certain areas. These fluctuations could explain unusual behaviour in the observed groundwater heads in these areas.

---

Resumen El Acuífero de la Pampa del Tamarugal (PTA) es una fuente importante de agua subterránea en el norte de Chile. En este estudio se desarrolla y calibra un modelo de flujo de agua subterránea para el periodo 1983–2004. El modelo reproduce razonablemente bien el campo de flujo observado y los componentes del balance hídrico. Se definen cinco escenarios para evaluar la respuesta a diferentes situaciones de bombeo. Estos escenarios muestran que con las tasas de descarga de bombeo actuales las presiones de agua subterránea continuarán en descenso. Para explicar las variaciones en los resultados del modelo debidas a incertidumbres en tasas de recarga promedio se han simulado realizaciones de recarga generadas aleatoriamente con diferentes niveles de incertidumbre. Las tasa de flujo por evaporación y el agua subterránea que fluye fuera del área modelizada parecen no haber sido afectadas por la incertidumbre en recarga mientras que los términos de almacenamiento pueden variar considerablemente. Para el escenario de bombeo más intenso, bajo las tasas de recarga generadas aleatoriamente, no es probable que el volumen descargado acumulado del acuífero al final del periodo de simulación sea mayor del 12% de la reserva de agua subterránea estimada. Las fluctuaciones en las presiones simuladas de agua subterránea debido a las incertidumbres en los valores de recarga promedio son más notorios en ciertas áreas. Las fluctuaciones podrían explicar el comportamiento inusual en las presiones de agua subterránea observadas en estas áreas.

---

Résumé L’aquifère de la Pampa del Tamarugal (PTA) représente une importante source d’eau souterraine dans le Nord du Chili. Dans cette étude, un modèle d’écoulement souterrain de cet aquifère est réalisé et calibré pour la période 1983–2004. Celui-ci reproduit le champs d’écoulement observé et les composantes du bilan d’eau raisonnablement bien. Cinq scénarios sont définis pour évaluer la réponse du système à des conditions de pompage différentes. D’après ces scénarios, les niveaux piézomètriques vont continuer à baisser sous l’effet des taux de pompage actuels. Pour tenir compte de la variabilité des résultats du modèle due aux incertitudes sur les taux de recharge moyens, des épisodes de recharge, générés aléatoirement et présentant différents niveaux d’incertitude sont simulés. Les taux d’évaporation et l’eau souterraine s’écoulant hors de la zone modélisée semblent non affectés par l’incertitude sur la recharge, les termes d’emmagasinement en revanche peuvent varier considérablement. Dans le cas du scénario avec le pompage le plus important et des taux de recharge aléatoires, il est peu probable que le volume cumulé pompé hors de l’aquifère à la fin de la simulation atteigne une valeur supérieure à 12% de la réserve d’eau souterraine estimée. Les fluctuations des niveaux piézomètriques simulés liées aux incertitudes sur les recharges moyennes sont plus perceptibles dans certaines zones. Ces fluctuations pourraient expliquer un comportement inhabituel des niveaux piézomètriques observés dans ces zones.

     

Climate-induced hydrological impacts on the groundwater system of the Pingtung Plain,Taiwan

The Pingtung Plain is one of the most important groundwater-resource areas in southwestern Taiwan. The overexploitation of groundwater in the last two decades has led to serious deterioration in the quantity and quality of groundwater resources in this area. Furthermore, the manifestation of climate change tends to induce the instability of surface-water resources and strengthen the importance of the groundwater resources. Southwestern Taiwan in particular shows decreasing tendencies in both the annual amount of precipitation and annual precipitation days. To effectively manage the groundwater resources of the Pingtung Plain, a numerical modeling approach is adopted to investigate the response of the groundwater system to climate variability. A hydrogeological model is constructed based on the information from geology, hydrogeology, and geochemistry. Applying the linear regression model of precipitation to the next two decades, the modeling result shows that the lowering water level in the proximal fan raises an alarm regarding the decrease of available groundwater in the stress of climate change, and the enlargement of the low-groundwater-level area on the coast signals the deterioration of water quantity and quality in the future. Suitable strategies for water-resource management in response to hydrological impacts of future climatic change are imperative.

---

Resumen La llanura Pingtung constituye uno de los recursos de aguas subterráneas más importantes en el Suroeste de Taiwan. La sobreexplotación de las aguas subterráneas en las dos últimas décadas ha dado lugar a un serio deterioro de la calidad y la cantidad de los recursos subterráneos en esta área. En particular el Suroeste de Taiwan muestra una tendencia decreciente en las cantidades de precipitación y en los días anuales de lluvia. Para gestionar de forma efectiva los recursos subterráneos en la Llanura Pingtung, se ha utilizado un modelo numérico aproximado para investigar la respuesta de las aguas subterráneas a la variabilidad climática. Un modelo hidrogeológico se construye a partir de la información geológica, hidrogeológica y geoquímica. Aplicando el modelo de regresión linear de la precipitación para las próximas dos décadas, el modelo resultante muestra que el descenso de los niveles de agua en el abanico proximal es alarmante, observándose el descenso del agua subterránea disponible en la presión del cambio climático, y el crecimiento del área de descensos de niveles de agua subterránea en la costa apunta a un deterioro de la cantidad y calidad del agua subterránea en el futuro. Se imponen pues, estrategias apropiadas para la gestión de los recursos hídricos en respuesta a los impactos hidrológicos del futuro cambio climático.

---

Résumé La plaine de Pingtung est l’une des plus importantes zones de ressource en eau souterraine du Sud-Ouest de Taiwan. La surexploitation de l’eau souterraine durant les deux dernières décennies a conduit à une sérieuse détérioration de la quantité et de la qualité des ressources en eau souterraine dans cette zone. De plus, la manifestation des changements climatiques tend à induire une instabilité des ressources en eau souterraine et renforce l’importance des ressources en eau souterraine. Le Sud-Ouest de Taiwan montre, en particulier, des précipitations annuelles et des nombres annuels de jour de pluie à la baisse. Pour gérer efficacement les ressources en eau souterraine de la Plaine de Pingtung, une approche par modélisation numérique est adoptée pour étudier la réponse du système hydrogéologique aux variabilités climatiques. La construction du modèle hydrogéologique est basée sur les informations géologiques, hydrogéologiques et géochimiques. En appliquant le modèle de régression linéaire aux précipitations pour les deux prochaines décades, le résultat de la modélisation montre que la baisse du niveau d’eau atteint un état alarmant au regard de la décroissance des eaux souterraines disponibles et la contrainte du changement climatique, tandis que l’extension de la zone de niveau bas des eaux souterraines à la cote indique une détérioration de la quantité et de la qualité de l’eau dans le futur. Des stratégies convenables pour la gestion des ressources en eau souterraine en réponse aux impacts hydrologiques de futurs changements climatiques sont impératifs.

     

Determination of specific yield using a water balance approach – case study of Torla Odha watershed in the Deccan Trap province, Maharastra State, India

Specific yield is an essential parameter for any groundwater management plan. Volumetric analysis in the domain of groundwater budgeting for the non-monsoon months has been undertaken for a typical watershed of the Deccan basalt province. The Torla Odha watershed covers an area of over 22 km² on a third-order tributary of the westerly flowing Bhima River. The watershed receives a normal annual rainfall of 643 mm. The entire water demand is supplied by dug wells, which penetrate a shallow aquifer. The specific yield was estimated by comparing the monthly net volume of water removed from the aquifer, with the volume of aquifer de-saturated, based on monthly water level data. The estimated specific yield ranges from 0.0019 in May to 0.0173 in November with an average value of 0.0093. A correlation of the groundwater levels with the detailed geology suggests that the higher specific yield value (0.017) corresponds to dewatering of the weathered zone within the shallow aquifer. The specific yield of the massive basalt immediately below the weathered zone varies from 0.0089 to 0.0103. The underlying vesicular basalt, which is dissected by sheet joints, has a relatively higher specific yield (0.0121). The massive basalt, which forms the base of the shallow aquifer system, has a lower specific yield from 0.0019 to 0.0022.

---

Résumé Le débit spécifique est un paramètre essentiel pour tout plan de gestion des eaux souterraines. Les analyses volumétriques, dans le cadre des bilans hydriques des eaux en dehors des mois de mousson, ont été entreprises pour un bassin-versant typique de la province basaltique du Deccan. Le bassin-versant du Torla Odha couvre une superficie de 22 km², et alimente l’affluent du troisième ordre de la rivière Bhima, qui coule vers l’Ouest. La pluviométrie annuelle atteint 643 mm. Toute la demande en eau es assurée par des puits foncés pénétrant dans l’aquifère phréatique. Le débit spécifique a été estimé en comparant le volume net mensuel d’eau captée dans l’aquifère, avec le volume de l’aquifère dé-saturé, basé sur les données des niveaux piézométriques mensuels. Le débit spécifique estimé s’étend entre 0,0019 en Mai et 0.173 en Novembre; la moyenne se situe à 0,0093. Une corrélation entre les niveaux des eaux souterraines et la géologie, suggère que les débits spécifiques les plus importants (0,017) correspondent aux zones altérées de l’aquifère phréatique. Le débit spécifique du massif basaltique, immédiatement sous la zone altérée, varie entre 0,0089 et 0,0103. Le basalte vésiculaire, situé juste en dessous et traversé par des diaclases parallèles, possède un débit spécifique sensiblement plus élevé (0,0121). Le basalte massif, qui forme la base de l’aquifère phréatique, possède un débit spécifique moins important, compris entre 0,0019 et 0,0022.

---

Resumen El rendimiento específico es un parámetro esencial para cualquier plan de manejo de aguas subterráneas. Se ha llevado a cabo el análisis volumétrico, en el entorno de balance de aguas subterráneas, para los meses sin monzón de una cuenca típica de la provincia de basaltos Deccan. La cuenca Rorla Odha cubre un área de 22 km² en un tributario de tercer orden del Río Bhima que fluye al oeste. La cuenca capta una lluvia anual normal de 643 mm. La totalidad de la demanda de agua es abastecida por pozos manuales que penetran un acuífero somero. Se estimó el rendimiento específico al comparar el volumen neto mensual de agua removido del acuífero con el volumen de agua de-saturado estimado a partir de datos de niveles de agua mensuales. Los valores estimados de rendimiento específico varían de 0.0019 en mayo a 0.0173 en noviembre con un valor promedio de 0.0093. La correlación de niveles de agua subterránea con la geología de detalle sugieren que el valor más alto (0.017) de rendimiento específico corresponden con el desaguado de la zona de intemperismo dentro del acuífero somero. El rendimiento específico del basalto masivo que se encuentra inmediatamente debajo de la zona de intemperismo varía de 0.0089 a 0.0103. El basalto vesicular subyacente, el cual está disectado por fracturas laminares, tiene un rendimiento específico relativamente más alto (0.0121). El basalto masivo, que forma la base del sistema de acuífero somero, tiene un rendimiento específico más bajo el cual varía de 0.0019 a 0.0022. Palabras clave: basalto Deccan. India. Rendimiento específico. Recarga de agua subterránea. Balance hídrico.

     

Flowpath structure in a limestone aquifer: multi-borehole logging investigations at the hydrogeological experimental site of Poitiers, France

Since 2002, a new hydrogeological experimental site (HES) has been developed in Poitiers, France. The overall research objective related to this site is to improve the understanding of flow and solute transport in carbonate aquifers. The benchmarking of various types of numerical models against the HES field data is one of the main research projects supported by three national scientific programs. Within this framework, the purpose of this report is to synthesize existing knowledge about both aquifer geology and flowpath structure, based on core analysis and well logging. The combined use of core-hole data, borehole logs and outcrop data provide valuable information about lithostratigraphy and fracturing. The comparison of flow-meter data with borehole images indicates that flowpaths in the HES aquifer are strongly constrained within (1) subhorizontal karstic structures and (2) subvertical fractures. The presence of karstic levels appears to be conditioned by the stratigraphy and are unevenly developed in the HES. The vertical interconnectivity between the three karstic levels seems to result from fractures occurrence in this limestone formation. Following the observations, data and interpretations, a conceptual model of the Dogger aquifer is proposed.

---

Résumé  La vocation première du Site Expérimental Hydrogéologique (SEH) de Poitiers, France, est de fournir des données pertinentes-y compris des chroniques ou expériences long terme-pour la caractérisation, la quantification et la modélisation des transferts d’eau et de solutés dans les aquifères carbonatés. Depuis 2002, les investigations menées sur le SEH ont permis de recueillir une importante quantité de données concernant un aquifère calcaire captif de 100 m d’épaisseur. Ces données servent actuellement de support à un exercice de modélisation comparée (benchmark) soutenu par trois programmes scientifiques nationaux. Dans le cadre de ce projet, l’objectif ici est de synthétiser les connaissances actuelles relatives à la géologie et la structure des écoulements dans l’aquifère du Dogger. La structure lithostratigraphique et la fracturation de l’aquifère sont analysées à partir des carottes de forages, des affleurements, et des diagraphies de puits. L’analyse croisée des données de débitmétrie verticale et des imageries de parois des forages indique que les écoulements sont principalement associés à des structures karstiques subhorizontales et des fractures subverticales. Les niveaux karstiques sont irrégulièrement répartis sur le SEH mais apparaissent fortement corrélés avec des niveaux stratigraphiques précis. La connectivité hydraulique des différents niveaux karstiques semble liée à la fracturation des calcaires.

---

Resumen  Desde 2002, se desarrolló en Potiers, Francia, un sitio experimental hidrogeológico (SEH). El objetivo general de la investigación relacionada a este sitio es mejorar el entendimiento del flujo y transporte de solutos en acuíferos calcáreos. El punto de referencia de varios tipos de modelos numéricos utilizando los datos del SEH es uno de los principales proyectos de investigación financiados por tres programas científicos nacionales. Dentro de este marco, el propósito de este reporte es sintetizar el conocimiento existente tanto acerca de la geología como de la estructura de los caminos de flujo, basado en el análisis de testigos y perfiles de pozos. El uso combinado de testigos, perfiles de pozos e datos de afloramientos provee información valiosa acerca de la litoestratigrafía y el fracturamiento. La comparación de datos de medidores de flujo con imágenes de pozos indica que los caminos de flujo en el acuífero del SEH están fuertemente condicionados por (1) estructuras cársticas subhorizontales, y (2) fracturas subverticales. La presencia de niveles cársticos pareciera estar condicionada por la estratigrafía y estar desarrollada en forma no igualmente distribuida dentro del SHE. La interconectividad vertical entre tres niveles cársticos sería el resultado de la ocurrencia de fracturas en la formación de las calizas. En base a las observaciones, los datos e interpretaciones, se propone un modelo conceptual del acuífero Dogger.

     

Interactions Between An Alluvial-Aquifer Wellfield And The Scioto River,Ohio, USA

Data from wellfields in the Scioto River valley aquifer near Piketon, Ohio, USA, were used to evaluate groundwater quantity and quality and to formulate wellfield-management concepts for future use of the alluvial aquifer. The percentage of induced river infiltration pumped from wells next to the river was estimated using a numerical groundwater-flow model and chemical mass-balance techniques. These analyses indicate that 50 to 88 percent of the wellfield water is derived from induced infiltration of river water. Model simulations correlate the percentage of well water from induced infiltration to wellfield pumpage, river stage, proximity of wells to the river, riverbed hydraulic conductivity, aquifer hydraulic conductivity, and river-stage level. Induced infiltration of river water has influenced the quality of groundwater, making the wellfield water more similar to the chemical character of the river than to that of the alluvial aquifer.

---

Resumen: Los datos procedentes de un conjunto de pozos en el acuífero de valle del Río Scioto, cerca de Pinkerton, Ohio; se usaron para evaluar la cantidad y calidad de las aguas subterráneas, así como para planificar la gestión de los pozos cercanos al mismo, fue estimado mediante un modelo numérico de flujo y técnicas de balance químico de masa. Los análisis indican que entre el 50 y 88 porciento del agua bombeada procede de infiltración inducida. Simulaciones realizadas con el modelo permitieron correlacionar el porcentaje del agua bombeada que proviene de infiltración inducida, con distintas variables: caudal de bombeo, tramo del río, proximidad de los pozos al río, conductividad hidráulica del lecho del río, conductividad hidráulica del acuífero y nivel del tramo de río. La infiltración inducida ha afectado a la calidad del agua de los pozos, cuyo carácter químico es más parecido al del agua del río que al del acuífero aluvial.

---

Résumé: Les données provenant de champs captants dans l'aquifère de la valée de la Scioto, près de Piketon (Ohio), ont été utilisées pour évaluer la quantité et la qualité de l'eau souterraine et pour définir les conditions de gestion de l'exploitation future de cet aquifère alluvial. La part de la réalimentation par la rivière, sous l'effet de pompages, dans les eaux des puits voisins du cours d'eau, a été estimée au moyen d'un modèle numérique d'écoulement souterrain et par des méthodes chimiques du bilan de matière. Ces analyses indiquent que 50 à 88 pourcent de l'eau pompée provient de la réalimentation par la rivière. Les simultations corrèlent la part pompée provenant de la réalimentation avec les pompages, le niveau de la rivière. La réalimentation provoquée à partir de la rivière a influencé la qualité de l'eau souterraine en conférant à l'eau des puits des caractères plus proches de ceux de la rivière que de ceux de l'aquifère alluvial.

     

Optimization of the DRASTIC method for groundwater vulnerability assessment via the use of simple statistical methods and GIS

The assessment of groundwater vulnerability to pollution has proved to be an effective tool for the delineation of protection zones in areas affected by groundwater contamination due to intensive fertilizer applications. By modifying and optimizing the well known and widely used DRASTIC model it was possible to predict the intrinsic vulnerability to pollution as well as the groundwater pollution risk more accurately. This method incorporated the use of simple statistical and geostatistical techniques for the revision of the factor ratings and weightings of all the DRASTIC parameters under a GIS environment. The criterion for these modifications was the correlation coefficient of each parameter with the nitrates concentration in groundwater. On the basis of their statistical significance, some parameters were subtracted from the DRASTIC equation, while land use was considered as an additional DRASTIC parameter. Following the above-mentioned modifications, the correlation coefficient between groundwater pollution risk and nitrates concentration was considerably improved and rose to 33% higher than the original method. The model was applied to a part of Trifilia province, Greece, which is considered to be a typical Mediterranean region with readily available hydrogeological and hydrochemical data.

---

Resumen La evaluación de vulnerabilidad del agua subterránea a la contaminación ha demostrado ser una herramienta efectiva para la delimitación de zonas de protección en áreas afectadas por contaminación de aguas subterráneas debido a aplicaciones intensivas de fertilizantes. Mediante la modificación y optimización del bien conocido y ampliamente utilizado modelo DRASTIC fue posible predecir la vulnerabilidad intrínseca a la contaminación así como el riesgo a la contaminación del agua subterránea con mayor precisión. Este método incorporó el uso de técnicas estadísticas y geoestadísticas simples para la revisión del pesaje y establecimiento de rangos de factores de todos los parámetros DRASTIC bajo un ambiente SIG. El criterio para estas modificaciones fue el coeficiente de correlación de cada parámetro con las concentraciones de nitraros en agua subterránea. En base al grado significativo estadístico algunos parámetros fueros sustraídos de la ecuación DRASTIC, mientras que se consideró el uso de la tierra como un parámetro adicional de DRASTIC. Siguiendo las modificaciones antes mencionadas se mejoró considerablemente el coeficiente de correlación entre el riesgo a la contaminación del agua subterránea y las concentraciones de nitratos incrementando en 33% su valor en relación al método original. El modelo se aplicó en una parte de la provincia Trifilia, Grecia, la cual se considera ser una región Mediterránea típica con datos hidroquímicos e hidrogeológicos fácilmente disponibles.

---

Résumé L’évaluation de la vulnérabilité des eaux souterraines à la pollution a montré qu’elle est un outil efficace pour délimiter les zones de protection dans les zones affectées par la contamination des eaux souterraines due à l’utilisation intensive de fertilisants. En modifiant et optimisant le modèle DRASTIC, bien connu et souvent utilisé, il a été possible de prédire la vulnérabilité intrinsèque à la pollution, et de définir plus précisément le risque de pollution. Cette méthode incorpore l’utilisation de simples techniques statistiques et géostatistiques, pour la révision des facteurs d’estimation et de pondération de tous les paramètres de DRASTIC sous S.I.G. Le critère de ces modifications était le coefficient de corrélation de chaque paramètre avec la concentration en nitrates dans les eaux souterraines. Sur la base de leur signification statistique, certains paramètres ont été soustraits de l’équation DRASTIC. Suivant les modifications mentionnées ci-dessus, le coefficient de corrélation entre les concentrations en nitrate et le risque de pollution des eaux souterraines a été considérablement amélioré de 33% par rapport à la méthode originale. Le modèle a été appliqué sur une partie de la province de Trifilia en Grèce, qui est considérée comme une région typiquement méditerranéenne avec des données hydrogéologiques et hydrochimiques aisément accessibles.

     

Hydrogeology and hydrochemistry of a shallow alluvial aquifer,western Saudi Arabia

A hydrogeological and hydrochemical study was conducted on a shallow alluvial aquifer, Wadi Wajj, in western Saudi Arabia to assess the influence of protection measures on groundwater quality. The hydrochemistry was assessed up-gradient and down-gradient from potential contamination sources in the main city in dry and wet seasons prior to and after the installation of major drainage and wastewater facilities. Wadi Wajj is an unconfined aquifer where water is stored and transmitted through fractured and weathered bedrock and the overlying alluvial sediments. Natural recharge to the aquifer is about 5% of rainfall-runoff. Hydrochemistry of the aquifer shows temporal and seasonal changes as influenced by protection measures and rainfall runoff. Both groundwater and runoff showed similar chemical signature, which is mostly of chloride-sulfate-bicarbonate and sodium-calcium type. Groundwater downstream of the city, though of poorer quality than upstream, showed significant improvement after the installation of a concrete runoff tunnel and a wastewater treatment plant. Concentrations of many of the groundwater quality indicators (e.g., TDS, coliform bacteria, and nitrate) exceed US Environmental Protection Agency drinking-water standards. Heavy metal content is, however, within allowable limits by local and international standards. The chemical analyses also suggest the strong influence of stream runoff and sewage water on the groundwater quality.

---

Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.

---

Résumé Une étude hydrogéologique et géochimique a été menée sur l’aquifère phréatique alluviale Wadi Wajj dans l’Ouest de l’Arabie Saoudite afin d’évaluer l’influence de mesures de protection sur la qualité de l’eau souterraine. L’hydrogéochimie de l’eau a été étudiée en amont et en aval de sources potentielles de contamination dans la ville principale, pendant les saisons sèches et humides, avant et après l’installation de réseaux majeurs de drainage et d’eaux usées. L’aquifère Wadi Wajj est libre, l’eau est stockée et s’écoule dans les sédiments alluviaux et dans le socle fracturé et altéré sous-jacent. La recharge naturelle de l’aquifère représente 5% des eaux de pluie et de ruissellement. L’hydrogéochimie de l’eau de l’aquifère montre que les changements saisonniers et temporaires sont influencés par les mesures de protection et par le ruissellement des eaux pluviales. L’eau souterraine et l’eau de ruissellement ont présenté le même faciès chimique, de type bicarbonatée-sulfatée-chlorée et calco-sodique. En aval de la ville, l’eau souterraine, bien que de moins bonne qualité qu’en amont, a présenté une nette amélioration après l’installation d’un système de récupération et d’une station de traitement des eaux usées. Les concentrations de plusieurs paramètres indicateurs de la qualité de l’eau (tels que la charge totale dissoute, les coliformes, et les nitrates) dépassent les normes de potabilité de l’eau de consommation de l’agence américaine de la protection de l’environnement. Les teneurs en métaux lourds en revanche n’excèdent pas les normes locales et internationales. Les analyses chimiques indiquent aussi l’influence importante de l’écoulement par ruissellement et des eaux usées sur la qualité de l’eau souterraine.

---

Resumen Un estudio hidrogeológico e hidroquímico fue hecho en un acuífero somero de Wadi Wajj, en Arabia Saudi oeste para evaluar la influecia de medidas de protección en la calidad del agua subterránea. La hidroquímica fue evaluada gradiente-arriba y gradiente-abajo de las fuentes potenciales de contaminación de la ciudad principal durante las estaciones seca y lluviosa, antes y después de la instalación de sistemas principales de drenaje y aguas servidas. Wadi Wajj es un acuífero no-confinado donde el agua es almacenada y transmitida a través de roca fracturada y meteorizada, y los sedimentos aluviales que le sobreyacen. La recarga natural del acuífero es de cerca del 5% de la precipitación-escorrentía. La hidroquímica del acuífero muestra cambios temporales y estacionales influenciados por las medidas de protección y la escorrentía de precipitación. Ambas, agua subterránea y escorrentía mostraron composición química similar, siendo mayoritariamente de tipos cloruro-sulfato-bicarbonato y sodio-calcio. El agua subterránea aguas arriba de la ciudad, aunque de calidad más pobre que aguas abajo, mostró significante mejoría después de la instalación de un tunel de concreto para escorrentía y una planta de tratamiento de aguas servidas. Las concentraciones de muchos de los indicadores de calidad de agua subterránea (e.g., STD, coliformes, y nitrato) exceden los estándares de la Agencia de Protección Ambiental USA para agua potable. El contenido de metales pesados está, sin embargo, dentro de los límites permisibles de los estándares locales e internacionales. Los análisis químicos también sugieren la fuerte influencia de la escorrentía y aguas residuales en la calidad del agua subterránea.