Multi-phase water-rhyolite interaction and ore fluid generation at Aljustrel, Portugal

Detailed study of the Aljustrel volcanic rocks, through petrography, mineral chemistry, whole-rock geochemistry, and previously published stable isotope data, show that the rocks experienced extreme metasomatism. The metasomatic phenomena presented are explained accurately by a sea-water convective model. The evolution of the sea-water-derived fluid as it progresses downwards through the rocks is evident, from cold, oxidising, Mg-bearing, to hot, reducing, Mg-free, and enriched in Si, Fe, Cu, Zn. Several lines of evidence (feldspar mineralogy, paragenetic relations between regional and ore-zone alteration, and stable isotope data) show that the temperature of the hydrothermal circulation increased through time, and that the water-to-rock ratio evolved from very high to low values. According to our data, about 30 ppm Cu and some 100 ppm Zn were removed from the Aljustrel Volcanics. This implies, for both metals, the same amount of leached rock to produce the 2 Mt Cu and 6 Mt Zn contained in the massive sulphide deposits of the Aljustrel area. With regard to Fe, about 1% was on average removed from the rocks, producing enough metal to account for all the Fe in the MS deposits and overlying chemical sediments. This study documents semi-quantitatively the provenance of the main metals in major massive sulphide deposits via sea water leaching of the immediate host rocks.

---

Resumen (translated by E. Pascual) Un estudio detallado de las rocas volcánicas de Aljustrel, mediante petrografia, química mineral, geoquímica de roca total y datos de isótopos estables previamente publicados muestra que estas rocas han sufrido metasomatismo en grado extremo. Los fenómenos metasomáticos presentados se explican con precisión mediante un modelo convectivo con agua de mar. Es evidente la evolución del fluido derivado del agua de mar conforme progresa en profundidad a través de las rocas, desde un fluido frío, oxidante y con magnesio a otro reductor, sin magnesio y enriquecido en Si, Fe, Cu y Zn. Varias líneas de evidencia (mineralogía de los feldespatos, relaciones paragenéticas entre la alteración regional y de las zonas mineralizadas, isótopos estables) muestran que la temperatura y la circulación hidrotermal se incrementaron con el tiempo y que la relación agua-roca evolucionó de valores muy altos a bajos. De acuerdo con nuestros datos, de las rocas volcánicas de Aljustrel se han extraído unas 30 ppm de Cu y unas 100 ppm de Zn. Esto implica, para ambos metales, que la misma roca lixiviada ha producido los 2 Mt de Cu y los 6 Mt de Zn contenidos en los depósitos de sulfuros masivos del área de Aljustrel. En cuanto al Fe, una media de alrededor de 1% fue extraída de las rocas, produciendo metal suficiente para explicar todo el Fe de los depósitos de sulfuros masivos y de los sedimentos químicos suprayacentes. Este estudio documenta de forma semicuantitativa la procedencia de los principales metales en depósitos mayores de sulfuros masivos, mediante lixiviado por agua de mar de sus rocas huéspedes inmediatas.

---

---

Received: 1 February 1997 / Accepted: 30 July 1997     

Magmatism in the Iberian Pyrite Belt: petrological constraints on a metallogenic model

Igneous formations associated with massive sulphide deposits in the Iberian Pyrite Belt (IPB) are essentially composed of basic lavas and dolerites, and dacitic to rhyolitic volcanites; intermediate lavas are subordinate. The basic rocks show variable geochemical characteristics: lavas and dolerites comparable to recent within-plate alkaline basalts seem restricted to the western and southern parts of the IPB, whereas basic rocks comparable to continental tholeiites or arc-related basalts occur across the whole belt. The felsic rocks are classified as calc-alkaline and belong to the “low-Al₂O₃ and high-Yb type”. At given SiO₂, Al₂O₃ and TiO₂ contents, they show variable Zr, Nb, and HREE contents. Heavy-rare-earth element fractionation decreases from the dacites to the rhyolites ([Gd/Yb]_N ∼ 1), whereas the negative Eu-anomaly becomes more pronounced. The characteristics of the rhyolites are typical of sulphide-fertile volcanic packages. Trace-element modelling suggests that the felsic rocks evolved from a dacitic parent magma through fractional crystallization of hornblende and plagioclase. Partial melting of an amphibolite protolith, which appears as the most probable model for the origin of this dacitic magma, requires a high T/P gradient in the crust. The occurrence of alkaline basalts and continental tholeiites is consistent with formation of the IPB in a tensional tectonic setting. However, the associated island-arc tholeiites suggest a location in a domain of plate convergence. Emplacement in a fore-arc basin over a recently accreted crustal segment is envisaged as a possible hypothesis to account for the geological and petrological constraints. A high geothermal gradient and eruption in a submarine tensional basin could have been two key ingredients for the development of massive sulphide deposits within the IPB.

---

Resumen (translated by E. Pascual) Las formaciones ígneas asociadas con los depósitos de sulfuros masivos de la Faja Pirítica Ibérica (IPB) se componen esencialmente de doleritas y lavas básicas y de rocas volcánicas dacíticas a riolíticas; las lavas intermedias son poco abundantes. Las rocas básicas muestran caracteres geoquímicos variables: lavas y doleritas comparables a basaltos alcalinos intraplaca recientes parecen hasta ahora restringidos a las partes W y S de la IPB, mientras que rocas básicas de caracteres comparables a las de toleítas continentales o basaltos relacionados con arcos aparecen en toda la zona. Las rocas ácidas se clasifican como calcoalcalinas del tipo “low-Al₂O₃, high-Yb”. Para un contenido dado en SiO₂, Al₂O₃ y TiO₂, muestran contenidos variables en Zr, Nb y REE. El contenido en tierras raras pesadas decrece de dacitas a riolitas ([Gd/Yb]_N ∼ 1), al tiempo que la anomalía de Eu se hace más pronunciada. La modelización de elementos trazas sugiere que las rocas evolucionaron a partir de un magma parental dacítico mediante cristalización fraccionada de hornblenda y plagioclasa. La fusión parcial de un protolito anfibolítico, que parece el modelo más plausible para el origen del magma dacítico, requiere un elevado gradiente T/P en la corteza. La existencia de basaltos alcalinos y de toleítas continentales es congruente con la formación de la IPB en un entorno tectónico distensivo. El emplazamiento en una cuenca “fore-arc”, en un segmento cortical de reciente acreción, se contempla como una hipótesis posible para explicar los caracteres geológicos y petrológicos. Un alto gradiente térmico, junto con la erupción en una cuenca submarina extensional, pueden haber sido los dos ingredientes clave en el desarrollo de los depósitos de sulfuros masivos en la IPB.

---

---

Received: 3 March 1996 / Accepted: 7 April 1997     

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.

---

     

Geomechanical model for fracture deformation under hydraulic, mechanical and thermal loads

Hydraulic flow and transport (heat and solute) within crystalline rocks is dominated by the fracture systems found within them. In situ stress conditions have a significant impact on the hydraulic, mechanical and thermal coupled processes, and quantification of these processes provides a key to understanding the often transient time-dependent behaviour of crystalline rocks. In this paper, a geomechanical model is presented which describes fracture closure as a function of effective stress and the changes in parameters such as storage, permeability, porosity and aperture. Allowing the fracture closure to be defined by the change in normal effective stress provides a link to the numerical consideration of parametrical changes due to rock stress alterations caused for example by changes in fracture fluid pressure, stress release, tectonic stress, thermal stress, orientation of the natural fracture in the pervasive stress system and local changes in a rock mass due to stress alteration. The model uses geometrical considerations based on a fractal distribution of apertures on the fracture surface, and applies well-established analytical elastic deformation solutions to calculate the deformation response to changes in effective stress. Analysis of the fractal generation method allows a standard normal distribution of fracture apertures to be predicted for all common fractal dimensions relating to a 2D surface. Changes in the fracture aperture are related to hydraulic functions such as permeability, storage and porosity of the fracture. The geomechanical model is experimentally validated against laboratory scale experimental data gained from the closure of a fractured sample recovered at a depth of 3,800 m from the KTB pilot borehole. Parameters for matching the experimental data were established externally, the only fitting parameters applied were the minimum and maximum contact area between the surfaces and the number of allowable contacts. The model provides an insight into the key processes determining the closure of a fracture, and can act as a material input function for numerical models linking the effects of changes in the stress field, hydraulic or thermal conditions, to the flow and transport parameters of a fractured system.

---

Résumé L’écoulement et le transport (chaleur et soluté) dans les roches cristallines sont dominés par les systèmes de fracture. Les conditions de stress in-situ ont un impact significatif sur l’hydraulique, les processus couplés de mécanique et thermique et la quantification de ces processus apportent une clé pour comprendre le comportement transitoire des roches cristallines. Dans cet article un modèle géomécanique est présenté, modèle qui décrit la fermeture des fractures comme une fonction de la contrainte effective et des changements de paramètres tels le coefficient d’emmagasinement, la perméabilité, la porosité et l’ouverture. En s’accordant que la fermeture des fractures est définit par les changements de la contrainte effective normale, on apporte le lien avec la considération numérique des changements paramétriques dus aux altérations de la contrainte des roches, causés par exemple par des variations de la pression des fluides dans les fractures, du dégagement de la contrainte, des contraintes tectoniques et thermiques, des orientations des fractures naturelles dans le système de contraintes pénétrantes, et des changements locaux dans un massif de roches dus à l’altération des contraintes. Le modèle utilise des considérations géométriques basées sur une distribution fractale des ouvertures à la surface des fractures, et permet d’établir des solutions analytiques de la déformation élastique pour calculer la réponse de la déformation à la contrainte effective. L’analyse de la méthode par génération fractale permet de prédire une distribution normale standard de l’ouverture des fractures, pour toutes les dimensions fractales en relation avec les surfaces 2D. Les changements dans l’ouverture des fractures sont mis en relation avec les fonctions hydrauliques tels la perméabilité, l’emmagasinement et la porosité de la fracture. Le modèle géoméchanique est expérimentalement validé à l’échelle du laboratoire sur un échantillon fracturé récupéré à une profondeur de 3,800 mètres sur le puits du site pilote KTB. Les paramètres du calibrage des données expérimentales ont été établies extérieurement, les seuls paramètres utilisés étant les surfaces de contact minimum et maximum, et le nombre de contacts permis. Le modèle apporte une connaissance perspicace sur le processus clé déterminant la fermeture des fractures, et peut servir de fonction input dans les modèles numériques reliant les effets des variations de la contrainte du terrain, les conditions hydrauliques ou thermales, les paramètres de l’écoulement et du transport et les systèmes de fracture.

---

Resumen El flujo hidráulico y transporte (de calor y solutos) dentro de rocas cristalinas está dominado por los sistemas de fracturas que se encuentran en ellas. Las condiciones de esfuerzos in-situ tienen un impacto significativo en los procesos aparejados termales, mecánicos e hidráulicos y la cuantificación de estos procesos aporta una clave para entender el frecuente comportamiento transitorio dependiente de las rocas cristalinas. En este artículo se presenta un modelo geomecánico que describe el cierre de fracturas en función del esfuerzo efectivo y los cambios en parámetros tal como almacenamiento, permeabilidad, porosidad y apertura. El definir el cierre de fractura mediante el cambio en esfuerzo normal efectivo aporta un vínculo con la consideración numérica de cambios paramétricos ocasionados por alteraciones de esfuerzos en la roca causadas, por ejemplo, por cambios en presión de fluidos en fractura, liberación de esfuerzo, esfuerzo tectónico, esfuerzo termal, orientación de fracturas naturales en el sistema de esfuerzos penetrante, y cambios locales en una masa rocosa ocasionados por alteración de esfuerzos. El modelo utiliza consideraciones geométricas basadas en la distribución fractal de aperturas en la superficie de fractura y aplica soluciones analíticas bien establecidas de deformación elástica para calcular la respuesta de deformación a cambios en el esfuerzo efectivo. Los análisis del método de generación fractal permiten predecir una distribución normal standard para la distribución de aperturas de fracturas para todas las dimensiones fractales comunes que se relacionan con una superficie 2D. Los cambios en la apertura de fractura se relacionan con funciones hidráulicas tal como permeabilidad, almacenamiento y porosidad de la fractura. El modelo geomecánico se ha validado experimentalmente en contra de datos experimentales a escala de laboratorio obtenidos a partir del cierre de una muestra fracturada recuperada a una profundidad de 3,800 m en el pozo piloto KTB. Se establecieron externamente parámetros que se ajustan a los datos experimentales, con los parámetros de ajuste aplicados que fueron el área máxima y mínima de contacto entre las superficies y el número de contactos permisibles. El modelo arroja luz sobre los procesos clave que determinan el cierre de una fractura y puede actuar como un material de función de entrada para modelos numéricos que vinculan los efectos de cambios en el campo de esfuerzos, condiciones termales o hidráulicas, con los parámetros de flujo y transporte de un sistema fracturado.

     

From well-test data to input to stochastic continuum models: effect of the variable support scale of the hydraulic data

The question of how well the true underlying hydraulic conductivity statistics of heterogeneous media are captured by well tests is addressed. The hydraulic conductivity value and the corresponding support volume associated with a theoretical well are correlated, causing a bias in the statistics derived from well-test analyses. Statistics derived from numerically simulated well tests are compared with the known underlying conductivity statistics and the results indicate an under-prediction by simulations at higher hydraulic conductivities. The deviation starts at about mean conductivity and can be as large as an order of magnitude, with the conductivity in the vicinity of the well defining the upper boundary. In other words, the conductivity value interpreted from the well test cannot be larger than the value that the well test first encounters. Consequently, for data in this simulation exercise, the standard deviation, if only determined for the upper range of the conductivity values, would be underestimated by a factor of 1.6–2. While this specific range is likely to depend on the scale and degree of the underlying heterogeneity as well as the duration of the test, the results should be indicative of a more general behaviour and are likely to occur in other heterogeneous data as well.

---

Resumen Se plantea la pregunta de qué tan bien son representadas en las pruebas de pozo, las estadísticas reales de conductividad hidráulica subyacente de medios heterogéneos. Son correlacionados el valor de conductividad hidráulica y el volumen de apoyo correspondiente asociado con un pozo teórico, causando una distorsión en las estadísticas derivadas del análisis de la prueba de pozo. Las estadísticas derivadas de las pruebas de pozo simuladas numéricamente son comparadas con las estadísticas de conductividad subyacente conocidas, y los resultados indican una sub-predicción por las simulaciones hechas con conductividades hidráulicas más altas. La desviación empieza casi con la conductividad media y puede ser tan grande como un orden de magnitud, con la conductividad en la vecindad del pozo definiendo el límite superior. En otras palabras, el valor de conductividad interpretado a partir de la prueba del pozo no puede ser más grande que el valor que la prueba de pozo encuentre primero. Por consiguiente, para los datos en este ejercicio de simulación, la desviación estándar, si solamente fue determinada para el rango superior de los valores de conductividad, se subestimaría en un factor de 1.6–2. Mientras es probable que este rango específico dependa de la escala y del grado de la heterogeneidad subyacente, así como de la duración de la prueba, los resultados deben ser indicativos de un comportamiento más general y son probables también de ocurrir en otros datos heterogéneos.

---

Résumé Nous posons ici la question de savoir dans quelle mesure les statistiques de la conductivité hydraulique des milieux hétérogènes pourrait être révélée par des essais de puits. La valeur de la conductivité hydraulique et le volume capté correspondant sont corrélés, créant un biais dans l’analyse des statistiques dérivées des essais de puits. Les statistiques en provenance de simulations numériques d’essais de puits sont comparées avec les statistiques de conductivités connues et les résultats indiquent une sous-évaluation par les simulations, pour les conductivités hydrauliques les plus élevées: la déviation commence à partir de la valeur moyenne de la conductivité et peut atteindre la magnitude d’un ordre de grandeur en considérant la conductivité mesurée au voisinage du puits. Autrement dit, la valeur de la conductivité interprétée via l’essais de pompage ne peut être plus importante que les premières valeurs rencontrées. Par conséquence, pour les données de cet exercice de simulation, la déviation standard sera sous-estimée d’un facteur compris entre 1.6–2 pour les valeurs les plus élevées. Tandis que l’échelle spécifique de valeurs est dépendante de l’échelle et du degré de l’hétérogénéité souterraine, de même que de la durée du test, les résultats pourraient être indicatifs d’un comportement plus général et seraient sans doute observables dans d’autres cas de données hétérogènes.

     

Lake bank filtration at Nainital,India: water-quality evaluation

There are different water-supply schemes in Uttarakhand, India to tap the water from streams, rivers and lakes. At Nainital, seven tube-wells (depths 22.6–36.7 m), located at a distance of <100 m from the lake, are being used to abstract (1) lake water after passage through the soil and (2) subsurface water/groundwater flowing towards the lake. Water samples from the lake and tube-wells were analyzed in monsoon and non-monsoon periods from 1997 to 2006. Total dissolved solids, EC, alkalinity and hardness were found to be marginally greater in tube-well waters. The difference in hydrochemistry of tube-well water was mainly due to variation in flow regimes during monsoon and non-monsoon periods. Results clearly indicate that lake water as such is not potable as it contains unacceptable levels of organic matter in terms of COD (~44 mg/L), coliforms (~15.6 × 10⁴ MPN/100 mL) and nutrients. Coliform bacteria and COD have not been detected in any of the tube-well water samples over the years. Lake water, treated by sand filters did not conform to drinking water standards. These investigations have led to the closure of the treatment facility and installation of two tube-wells in addition to the existing five tube-wells.

---

Résumé  Il existe divers projets d’alimentation en eau dans l’état d’Uttarakhand, Inde, afin de capter l’eau de ruisseaux, de rivières et de lacs. A Nainital, sept puits tubés (profondeur de 22.6–36.7 m), situés à une distance < à 100 m du lac, sont utilisés pour prélever (1) de l’eau du lac après transit à travers le sol et (2) de l’eau de sub-surface/eau souterraine s’écoulant vers le lac. Des échantillons d’eau du lac et des puits tubés ont été analysés en périodes de mousson et de celles sans mousson 1997 à 2006. Résidu sec, C.E., alcalinité et dureté ont été trouvés marginalement supérieurs dans l’eau des puits tubés. La différence d’hydrochimie de l’eau des puits tubés était surtout due à la variation des régimes d’écoulement pendant les périodes de mousson et de celles sans mousson. Les résultats indiquent clairement que l’eau du lac en tant que telle n’est pas potable car elle contient des teneurs inacceptables de matière organique en termes de COD (~44 mg/L), de coliformes (~15.6 × 10⁴ MPN/100 mL) et d’éléments nutritifs. Des bactéries coliformes et du COD n’ont été détectés dans aucun des échantillons d’eau de puits tubés au fil des années. L’eau du lac traitée par des filtres à sable ne se conformait pas aux normes de l’eau potable. Ces recherches ont conduit à la fermeture de l’installation de traitement et à l’implantation de deux puits tubés en plus des cinq puits existants.

---

Resumen  En el estado de Uttarakhand, India, existen diferentes esquemas de abastecimiento de agua que explotan agua de arroyos, ríos y lagos. En Nainital, siete pozos (profundidades entre 22.6–36.7 m), ubicados a una distancia de <100 m del lago, se usan para extraer (1) agua del lago luego de su pasaje a través del suelo y (2) agua superficial y subterránea que fluye hacia el lago. En períodos de monzón y de no monzón de 1997 a 2006 se han analizado muestras de agua del lago y de las captaciones. Se halló que el agua de los pozos es ligeramente mayor en términos del total de sólidos disueltos, la conductividad eléctrica, la alcalinidad y la dureza. La diferencia en la hidroquímica del agua de las perforaciones se debe principalmente a la variación de los regímenes de flujo durante los períodos de monzón y de no monzón. Los resultados claramente indican que el agua del lago no es potable por su contenido inaceptable de materia orgánica medida como demanda de carbono orgánico/oxígeno –DCO- (~44 mg/L), coliformes (~15.6 × 10⁴ NMP/100 mL) y nutrientes. En el período, no se han detectado bacterias coliformes ni DCO en las muestras de agua de las captaciones. El agua del lago, tratada con filtros de arena, no conformó los estándares de agua para bebida. Estas investigaciones han demostrado la necesidad de clausurar las instalaciones de tratamiento y la adición de dos captaciones a las cinco ya existentes.

     

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.

     

Mineralogy of suspended sediment in three karst springs

Springs in karstic carbonate rocks frequently carry a sediment load as well as a dissolved load. Analysis of morphology and mineralogy of suspended sediment from three contrasting karst springs reveals a suite of clastic particles that reflect both source areas and processes that take place within the aquifer. Nolte Spring in Lancaster County, Pennsylvania, USA, discharges sediment of apparently precipitated calcite, indicating that at some point in the aquifer or vadose zone, water exceeds saturation with respect to calcite. Sediment morphologies and chemical conditions in the aquifer point to two different scenarios for this precipitation. The other two springs, Arch Spring in Blair County, Pennsylvania and Bushkill Spring in Northampton County, Pennsylvania, show no evidence of calcite precipitation. Arch Spring discharges mainly layer silicates, while Bushkill Spring discharges mainly silica.

---

Resumen Los manantiales en rocas carbonatadas karstficadas tienen frecuentemente una carga de sedimentos y de sólidos disueltos. Los análisis de la morfología y mineralogía de los sedimentos suspendidos de tres manantiales kársticos contrastados revela un conjunto de partículas clásticas que reflejan ambas áreas fuente y los procesos que tienen lugar dentro del acuífero. Nolte Spring en el Condado de Lancaster, Pennsylvania, USA, descarga sedimento de calcita aparentemente precipitada indicando que en cierto punto del acuífero o en la zona vadosa, el agua supera la saturación con respecto a la calcita. La morfología de los sedimentos y las condiciones químicas en el acuífero apuntan a dos escenarios diferentes para esta precipitación. Los otras dos manantiales, Arch Spring en el Condado de Blair, Pennsylvania y Bushkill Spring en el Condado de Norhtampton, Pennsylvania, no muestran evidencia de precitación de calcita. Arch Spring descarga principalmente filosilicatos mientras que Bushkill Spring descarga fundamentalmente Sílice.

---

Résumé Les sources situées en contexte carbonaté karstique charrient fréquemment, outre les matières dissoutes, une charge solide. La morphologie et la minéralogie des sédiments en suspension issus de trois sources karstiques de caractéristiques distinctes ont été analysées : une série de particules clastiques, reflétant à la fois les secteurs des sources et les processus internes de l’aquifère, en ressort. Les sédiments issus de Nolte Spring (Comté de Lancaster, Pennsylvanie, Etats-Unis) sont constitués de calcite apparemment précipitée: ceci indique que l’eau atteint la sursaturation vis-à-vis de la calcite à un certain point de l’aquifère ou de la zone non-saturée. La morphologie des sédiments et les conditions chimiques au sein de l’aquifère mènent à deux scénarii différents pour expliquer cette précipitation. Les deux autres sources, Arch Spring (Comté de Blair, Pennsylvanie) et Bushkill Spring (Comté de Northampton, Pennsylvanie), ne présentent aucun indice de précipitation de calcite. Arch Spring exporte essentiellement des phyllosilicates, alors que la silice est majoritaire à Bushkill Spring.

     

Hydrogeological environment and groundwater occurrences of the alluvial aquifers in parts of the Central Ganga Plain, Uttar Pradesh, India

A detailed hydrogeological investigation was carried out in parts of the Central Ganga Plain, India, with the objective of assessing the aquifer framework and its resource potential. The area was studied because of its dual hydrogeological situation, that is water logging and soil salinization in the canal command areas and depletion of aquifers in the western part of the basin. A comprehensive investigation of the aquifer system between the Ganga River and Kali River revealed its lateral and vertical dimensions and hydrogeological characteristics. Moreover, study of the groundwater occurrences, movements and behaviour, in terms of water level fluctuation with time and space, confirms the coexistence of over exploitation as well as water logging in the area.

---

Electronic Supplementary Material Supplementary material is available for this article at

---

Resumen Una investigación hidrogeológica detallada se llevó a cabo en partes de la Llanura de Ganga Central, India, con el objetivo de evaluar la estructura del acuífero y su potencial del recurso. El área fue estudiada por su doble situación hidrogeológica, es decir la saturación con agua y salinización de suelos en las áreas dominadas por el canal y vaciamiento de acuíferos en el la parte occidental de la cuenca. Una investigación completa del sistema acuífero entre el Río Ganga y el Río Kali, reveló sus dimensiones verticales y laterales y las características hidrogeológicas. Además, los estudios sobre la ocurrencia del agua subterránea, sus movimientos y comportamiento, en lo que se refiere al nivel de agua, y a su fluctuación en el tiempo y el espacio, confirma la co-existencia en el área de sobre- explotación así como de saturación con agua.

---

Résumé L'objectif de cette étude hydrogéologique détaillée de portions de la Plaine Centrale du Gange est de déreminer la structure aquifère et la ressource potentielle. L'intérêt de la zone repose sur sa dualité du point de vue hydrogéologique, les zones influencées par le canal présentent une remontée de la nappe avec une salinisation des sols, la portion Ouest du bassin présente une baisse du niveau des aquifères. Par cette étude, le système aquifère compris entre la rivière du Gange et la rivière Kali révèle ses dimensions latérales et verticales ainsi que ses caractéristiques hydrogéologiques. De plus, l'étude des événements, des mouvements et du régime hydrogéologique affectant le niveau phréatique confirme la co-existence de surexploitation et de saturation des sols dans la région.

     

On the solutions of modeling a slug test performed in a two-zone confined aquifer

The slug test is a popular field method for estimating aquifer parameters. As a finite-thickness skin is present at a wellbore, the result of test-data analysis may represent average hydraulic properties of the skin and formation zones. Moench and Hsieh (Int. Assoc. of Hydrogeology, The IAH 17th International Congress on Hydrology of Rocks of Low Permeability, Int. Assoc. of Hydrogeology. Tucson, Arizona, 7–12 January) provided the Laplace-domain solution for a confined aquifer with a finite-thickness skin that satisfies the governing equations, boundary conditions, and continuity requirements at the interface of the skin and formation zones. Two new formulae presented here relate the modified zero-order Bessel functions of the first and second kinds to their first and second derivatives. Their solution is evaluated using the modified Crump approach for the Laplace inversion. The numerical values obtained by the modified Crump approach are compared with others based on finite-element modeling. The results indicate that the values based on finite-element modeling exhibit poor accuracy, attributed to errors in the mass balance equation and interface condition in their mathematical model.

---

Resumen Los ensayos slug o tipo pulso son métodos populares de campo para estimar parámetros de acuíferos. Debido a que en los sondeos existe una membrana de espesor finito, el resultado del análisis de los datos del ensayo pueden representar las propiedades hidráulicas promedio de las zonas de la formación y dicha membrana. Moench AF y Hsieh PA (Asociación Internacional de Hidrogeólogos, 7–12 Enero 1985) Análisis de datos de prueba slug en un pozo con membrana de espesor finito. 17 Congreso Internacional de la Asociación Internacional de Hidrogeólogos sobre Hidrología de Rocas de Baja Permeabilidad. Asociación Internacional de Hidrogeólogos. Tucson, Arizona. 7–12 Enero. han aportado la solución en el ámbito de Laplace para un acuífero confinado con una membrana de espesor finito que satisface los requerimientos de las ecuaciones dominantes, las condiciones restrictivas, y de continuidad en la interfase de las zonas de formación y membrana. Las dos nuevas fórmulas que se presentan aquí relacionan las funciones modificadas Bessel de orden cero de la primera y segunda clase con su primera y segunda derivadas. Se evalúa su solución usando el enfoque modificado de Crump para la inversión Laplace. Los valores numéricos obtenidos mediante el enfoque modificado de Crump se comparan con otros valores numéricos basados en modelos de elementos finitos. Los resultados indican que los valores basados en el modelo de elementos finitos presentan poca precisión, que se atribuye a errores en la ecuación de balance de masa y a las condiciones de interfase en su modelo matemático.

---

Résumé Le slug-test représente une méthode populaire de terrain pour estimer les paramètres d’un aquifère. Du fait de la présence d’une couche à l’épaisseur limitée sur un site de forages, le résultat de l’analyse des données du test pourrait représenté le terrain mais aussi cette couche. Moench AF et Hsieh PA (Association Internationale d’Hydrogéologie, 7–12 Jan 1985) Analyses de données de slug test dans un puits présentant une couche d’épaisseur finie. 17ème congrès international de l’AIH sur l’Hydrologie des roches de faible perméabilité, Association Internationale d’Hydrogéologie. Tucson, Arizona, 7–12 Jan. apportent une solution du domaine de Laplace pour un aquifère confiné avec une épaisseur limitée qui satisfait l’équation en question, les conditions aux limites, et les besoins de continuité à l’interface de la couche et de l’encaissant. Deux nouvelles formules présentées dans cet article met en relation les fonctions de Bessel modifiées d’ordre zéro du premier et second type, avec leur dérivées première et deuxième. Leurs solutions sont évaluées en utilisant une approche de Crump modifiée pour l’inversion de Laplace. Les valeurs numériques obtenues avec l’approche de Crump sont comparées avec d’autres approches utilisant les modèles aux éléments finis. Les résultats indiquent que les valeurs basées sur les modèles aux éléments finis montrent une faible précision, du fait d’erreurs dans le bilan hydrologique et les conditions aux limites dans le modèle.

     

Multiple-method estimation of recharge rates at diverse locations in the North Carolina Coastal Plain, USA

Recharge rates determined at diverse study sites in a shallow, unconfined aquifer differed from one another depending on the analytical method used and on each method’s applicability and limitations. Total recharge was quantified with saturated-zone methods using water-table fluctuations at seven sites in North Carolina, USA and using groundwater-age dating at three of the seven sites; at two of the sites, potential recharge was quantified with an unsaturated-zone method using Darcy’s law; and at five of the sites, net recharge was quantified with a stream hydrograph-separation method using streamflow-recession curves. Historical mean net recharge was 25 to 69% of the historical total recharge rates. The large disparity is attributed to groundwater losses between recharge and discharge areas, primarily by evapotranspiration and seepage to underlying aquifers. The spatial distribution of historical mean annual total recharge did not vary between landscape units, as suggested in a previous study. Similarly, total recharge did not correlate significantly with mean annual rainfall, mean annual water table depth, or the surficial soil properties of percent clay and bulk density. Total recharge did correlate significantly with the surficial soil properties of percent sand and percent silt.

---

Résumé Les taux de recharge déterminés sur divers sites d’études dans un aquifère phréatique libre, varient les uns par rapport aux autres selon la méthode analytique utilisée et les limites et l’applicabilité de chaque méthode. La recharge totale a été quantifiée avec des méthodes en zone saturée utilisant les fluctuations des nappes sur sept sites de la Caroline du Nord, USA, et en utilisant les datations des eaux souterraines sur trois des sept sites ; sur deux des sites, la recharge potentielle a été quantifiée avec une méthode en zone non-saturée utilisant la loi de Darcy, et sur deux sites, la recharge nette a été quantifiée suivant une méthode de séparation des hydrographes de cours d’eau utilisant les courbes de récession du débit. La moyenne historique de la recharge nette est comprise entre 25 et 69% des taux historiques de recharge totale. L’importante disparité est attribuée aux pertes en eaux souterraines entre les zones de recharge et les zones d’exutoire, d’abord par évapotranspiration et ensuite par infiltration vers les aquifères sous-jacents. La distribution spatiale de la recharge historique moyenne annuelle et totale ne varie pas selon les unités paysagères, comme cela était suggéré dans une étude antérieure. De même, la recharge totale n’est pas corrélée significativement avec la moyenne annuelle des précipitations, la moyenne annuelle de la profondeur de la nappe, ou les propriétés de surface des sols que sont le pourcentage d’argile ou la densité apparente. La recharge totale est corrélée significativement avec les propriétés de surface du sol et les teneurs en sable et en silt.

---

Resumen Las tasas de recarga determinadas en diversos sitios de estudio en un acuífero somero no confinado difieren entre sí dependiendo del método analítico usado y de las limitaciones y aplicabilidad de cada método. Se cuantificó la recarga total con métodos de zona saturada usando fluctuaciones del nivel freático en siete sitios en Carolina del Norte, Estados Unidos, y usando datación de edades de agua subterránea en tres de los siete sitios; en dos de los sitios se cuantificó la recarga potencial con un método de zona no saturada usando la Ley de Darcy y en cinco de los sitios se cuantificó la recarga neta con el método de separación hidrográfica usando curvas de recesión de flujo en arroyos. La recarga neta media histórica varió entre 25 y 69% de las tasas de recarga total histórica. La diferencia tan grande se atribuye a pérdidas de agua subterránea entre áreas de descarga y recarga, principalmente por evapotranspiración y escurrimiento hacia acuíferos subyacentes. La distribución espacial de la recarga total anual media histórica no varió entre las unidades de paisaje, como ya se había sugerido en un estudio previo. De modo similar, la recarga total no guarda correlación significativa con la lluvia media anual, la profundidad media anual del nivel freático, o las propiedades superficiales del suelo, particularmente porcentaje de arcilla y densidad volumétrica. La recarga total si tuvo una correlación significativa con las propiedades superficiales del suelo tal como porcentaje de arena y porcentaje de limo.

     

Estimating groundwater recharge through tills: a sensitivity analysis of soil moisture budgets and till properties in Ireland

The study highlights the dangers of limiting recharge calculations in humid climates to meteorological and soil considerations. It highlights the importance of developing a conceptual understanding of the influence of geology on recharge mechanisms and recharge rates. Uncertainties in recharge estimates through tills are examined by a study of the available literature, combined with sensitivity analyses of soil moisture budget parameters and of hypothetical scenarios of till properties and hydraulic gradients that are realistic for Irish conditions. The sensitivity analyses demonstrate that till properties have a greater influence over recharge than the soil moisture budgeting parameters. Results are considered in terms of recharge coefficients, representing the proportion of effective precipitation that becomes actual recharge to the aquifer. The literature review identifies a range in recharge coefficients from 4% to 90% in selected field and catchment scale studies from bedrock aquifers in Ireland and the U.K. The sensitivity analyses are able to simulate a similar range in recharge coefficient of 2% to 80% by varying only till permeability, thickness and vertical hydraulic gradients. In the scenarios examined, the greatest sensitivity to vertical till permeability lies between 0.001 m/day and 0.01 m/day. The sensitivity to soil moisture budgeting parameters is much less apparent. It is concluded that this is due to the dominance of grassland in Ireland and the relatively wet summers, resulting in the limited development of soil moisture deficits by comparison with some parts of the U.K.

---

Résumé Cette étude met en lumière les dangers résultants de la pratique consultant à limiter les calculs de recharge aux données météorologiques ainsi qu’aux considérations liées au sol. Elle met en lumière l’importance de développer des modèles conceptuels comprenant l’influence de la géologie sur les mécanismes de la recharge et les taux de recharge. Les incertitudes sur les estimations de la recharge à travers les argiles sont examinées à travers les données disponibles dans la littérature, combinées avec une analyse de sensibilité des paramètres du bilan hydrique et des scénarios hypothétiques reprenant les propriétés des argiles et les gradients hydrauliques qui sont réalistes dans le contexte irlandais de l’étude. Les analyses de sensibilité démontrent que les propriétés des argiles ont une plus grande influence sur la recharge que les paramètres du bilan hydrique des sols. Les résultats sont considérés en terme de coefficients de recharge, représentant la proportion de précipitation effective qui devient la recharge réelle de l’aquifère. La littérature identifie des coefficients de recharge allant de 4 à 90% sur des terrains sélectionnés à l’échelle de bassins-versants, sur les aquifères de socle en Irlande et au Royaume-Uni. Les analyses de sensibilité sont capables de simuler une variation aussi large de coefficients de recharge, en ne changeant rien que la perméabilité des argiles, leur épaisseur et les gradients hydrauliques verticaux. Dans les scénarios étudiés, la meilleure sensibilité à la perméabilité verticale se situe entre 0.001 et 0.01 m/jour. La sensibilité au bilan hydrique est beaucoup moins apparente. On en conclut que ceci est du à la dominance des prairies en Ireland et aux étés relativement humides, résultant des déficits en eau des sols relativement modérés en comparaison avec certaines zones du Royaume-uni.

---

Resumen El estudio resalta los peligros de limitar los cálculos de recarga en climas húmedos a las consideraciones meteorológicas y del suelo. Se destaca la importancia de desarrollar un entendimiento conceptual de la influencia de la geología en los mecanismos de recarga y los ritmos de recarga. Se examinan las incertidumbres en los estimados de recarga a través de tills mediante un estudio de la literatura disponible combinado con análisis de sensibilidad de parámetros de balance de humedad del suelo y de escenarios hipotéticos de las propiedades del till y gradientes hidráulicos que son realistas para condiciones típicas de Irlanda. Los análisis de sensibilidad demuestran que las propiedades del till tienen una influencia más grande en la recarga que los parámetros del balance de humedad del suelo. Se consideran los resultados en términos de coeficientes de recarga, representando la proporción de precipitación efectiva que se convierte en recarga real al acuífero. La revisión de literatura identifica un rango de coeficientes de recarga que varía de 4% a 90% en estudios seleccionados de cuencas y escala de campo realizados en acuíferos de macizo rocoso en Irlanda y el Reino Unido. Los análisis de sensibilidad pueden simular un rango similar en los coeficientes de recarga de 2% a 80% al variar únicamente la permeabilidad del till, el espesor y los gradientes verticales hidráulicos. En los escenarios examinados, la relación de sensibilidad más grande con la permeabilidad vertical del till se encuentra entre 0.001 m/día y 0.01 m/día. La sensibilidad en relación a los parámetros del balance de humedad del suelo es mucho menos aparente. Se concluye que esto se debe a la predominancia de tierras con hierba en Irlanda y los veranos relativamente húmedos, lo que resulta en el desarrollo limitado de déficit de humedad de suelo en comparación con algunas partes del Reino Unido.

     

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.

     

Roughness control on hydraulic conductivity in fractured rocks

The influence of joint roughness on the typologies of fluid flow inside fractures is well known and, thanks to experiences in the field of hydraulics, it has been studied from both a physical and mathematical point of view. Nevertheless, the formulations adopted by traditional hydraulic models are hardly applicable in the geological field, because of the difficulty encountered in the roughness parameter estimation. Normally this parameter can be estimated using the joint roughness coefficient (JRC), which considers both the asperity height and its regularity and directional trend. The main advantage in using the JRC arises from the fact that it can easily be obtained from geological-technical surveys and from comparison with the standard Barton profiles. Some relationships have been built up that allow for the estimation of the hydraulic conductivity tensor (an essential parameter for understanding water flow in fractured rock masses), not only as a function of traditional parameters like aperture, spacing, dip and dip direction, etc., but also of joint roughness, precisely expressed in terms of the roughness coefficient. These relationships have been studied initially from a theoretical point of view and then practically, through laboratory investigations.

---

Resumen  Se conoce muy bien la influencia de la rugosidad de las grietas en las tipologías del flujo de fluidos a lo interior de las fracturas y gracias a las experiencias en el campo de hidráulica ha sido posible estudiarla desde puntos de vista matemáticos y físicos. Sin embargo, las formulaciones adoptadas por los modelos hidráulicos tradicionales tienen poca aplicabilidad en el campo geológico debido a la dificultad relacionada con la estimación del parámetro de rugosidad. Normalmente este parámetro puede estimarse usando el coeficiente de rugosidad de grieta (JRC) el cual considera tanto la altura de la aspereza como su regularidad y tendencia direccional. La principal ventaja de utilizar el JRC se deriva del hecho que puede obtenerse fácilmente de levantamientos técnico-geológicos y de la comparación con los perfiles Standard Barton. Se han construido algunas relaciones que permiten la estimación del tensor de conductividad hidráulica (un parámetro esencial para el entendimiento del flujo de agua en masas de roca fracturadas), no solo en función de parámetros tradicionales como apertura, espaciado, buzamiento y dirección de buzamiento, etc., sino también en función de la rugosidad de la grieta estimada con precisión en términos del coeficiente de rugosidad. Estas relaciones se han estudiado inicialmente desde un punto de vista teórico y luego de modo práctico a través de investigaciones de laboratorio.

---

Résumé  L’influence de la rugosité des joints sur les types d’écoulement de fluide dans les fractures est bien connue et a été étudiée aussi bien du point de vue physique que mathématique grace à des expériences menées dans le domaine de l’hydraulique. Cependant les formulations adoptées dans les modèles hydrauliques traditionnels sont difficilement applicables dans le domaine de la géologie à cause de la difficulté rencontrée pour estimer la rugosité. Ce paramètre peut normalement être apprécié grace au coefficient de rugosité du joint (JRC), lequel prend en compte à la fois la hauteur de l’aspérité ainsi que sa régularité et sa direction. Le principal avantage dans l’utilisation du JRC réside dans le fait qu’il peut facilement être obtenu à partir d’études techniques-géologiques et par comparaison avec la classification de Barton. Des relations qui permettent une estimation du tenseur de conductivité hydraulique (un paramètre essentiel pour comprendre l’écoulement de l’eau dans les masses rocheuses fracturées) ont été élaborées, pas seulement en fonction de paramètres traditionnels tels que l’ouverture, l’espacement, l’inclinaison et la direction d’inclinaison, etc , mais aussi en prenant en compte la rugosité des joints à travers le coefficient de rugosité. Ces relations ont initialement été étudiées d’un point de vue théorique puis expérimentalement à travers des recherches en laboratoire.

     

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.

     

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.

     

Dryland salinity processes within the discharge zone of a local groundwater system,southeastern Australia

Detailed study of a localised saline discharge zone in southeastern Australia shows that the salinisation is mostly due to the shallow water table (<1–2 m from the surface). Direct evaporation, particularly in summer, leads to extremely high soil–water salinities at the surface, even though the underlying groundwater is moderately fresh. Groundwater discharge is localised at a break of slope, where the water table intersects the surface, and where the transition from permeable sands to clay-rich sediments inhibits lateral groundwater flow. Higher salt concentrations build up in the clays because of the long residence times during which soil-waters are exposed to evapotranspiration and the reduced potential for salts to be flushed from the sediments. As a result the saline discharge area does not correspond to the part of the site with the largest salt store. Results of the study demonstrate that for dryland salinisation to occur, the groundwater beneath the discharge zone need not be saline, and the presence of a large salt store does not necessarily lead to problems of dryland salinisation if, as in the clay-rich sediments at the site, the salt lies below the pasture root zone. Furthermore, mobilisation of salt stores within low permeability sediments by rising groundwater may be minor.

---

Résumé L’étude détaillée d’une zone localisée d’émergence salée au sud-est de l’Australie montre que la salinisation est surtout due à la nappe phréatique (<1–2 m sous la surface). L’évaporation directe, particulièrement durant l’été, conduit à des salinités de l’eau du sol extrêmement élevées à la surface, même si l’eau souterraine sous-jacente est modérément douce. L’émergence de l’eau souterraine est localisée à la rupture de pente, là où la nappe phréatique rencontre la surface du sol et où la transition entre sables perméables et sédiments riches en argiles inhibe les écoulements d’eau souterraine latéraux. Les plus fortes concentrations en sel s’accumulent dans les argiles du fait de temps de résidence élevés, durant lesquels les eaux du sol sont exposées à l’évapotranspiration et à un lessivage réduit des sédiments. Il en résulte que l’aire d’émergence des eaux salées ne correspond pas à la partie du site rencontrant la réserve de sel la plus importante. Les résultats de l’étude démontrent que pour que la salinisation d’une zone aride devienne effective, l’eau souterraine sous la zone ne doit pas être forcément salée, et la présence d’une zone étendue de réserve de sel ne conduit pas forcément à des problèmes de salinisation de zones arides si, comme dans les sédiments argileux du site, le sel ne repose pas sous la zone de pature. De plus, la mobilisation des réserves de sel dans les sédiments peu perméables par la montée du niveau de l’eau souterraine devrait être mineure.

---

Resumen El estudio detallado de una zona local de descarga salina en el sudeste de Australia, muestra que la salinización es principalmente debido a un nivel freático poco profundo (<1–2 m de la superficie). La evaporación Directa, particularmente en verano, conlleva a salinidades de suelo-agua sumamente altas en la superficie, aunque el agua subterránea subyacente es moderadamente dulce. La descarga de Agua subterránea se localiza en una interrupción de la ladera, dónde el nivel freático intercepta la superficie, y donde la transición de las arenas permeables a los sedimentos ricos en arcilla inhibe el flujo lateral del agua subterránea. Las concentraciones de sal más altas se forman en las arcillas debido a los tiempos de residencia largos durante los cuales se exponen el conjunto suelo-agua a la evapotranspiración y también por el potencial reducido para las sales de ser expulsadas de los sedimentos. Como resultado el área de la descarga salina no corresponde a la parte del sitio con el contenido de sal más grande. Los resultados del estudio demuestran que para que ocurra la salinización en terrenos secos, el agua subterránea bajo la de zona de descarga no necesita ser salina, y que la presencia de un almacenamiento de sal grande, no necesariamente lleva a los problemas de salinización en terrenos secos si, como en los sedimentos ricos en arcilla del sitio, la sal yace debajo de la zona de raíz de la pastura. Además, la movilidad de depósitos de sal dentro de los sedimentos de permeabilidad baja pueden ser menores, por causa del agua subterránea ascendente.

     

Conceptual hydrodynamic model of the Pamukkale hydrothermal field, southwestern Turkey, based on hydrochemical and isotopic data

The results of study on the hydrochemical and isotope characteristics of shallow and deep waters at Pamukkale hydrothermal field Turkey are described in order to obtain a better understanding of the hydrological circulation. The field can be grouped into two groundwater sub-systems; cold water springs of Ca–HCO₃ type (10–12 °C), and CO₂-rich thermal waters of Ca–HCO₃–SO₄ type (25–58 °C). The occurrence of these water types is closely related to the morphology of the region, where intense tectonism formed horst and graben structures. Hence, two hydrogeological systems were defined: a deep geothermal system which is related to extensive and deep circulation of meteoric water in the regional flow system, and a shallow system which is related to local groundwater flow through sedimentary strata. The meteoric water falling at higher elevations percolates to the local groundwater system at a shallow level and flows to the deep geothermal system. During a deep convection cycle from a recharge to discharge area, the cold water attains heat from the asthenospheric intrusions, causing it to ascend. Variations of chemical and isotopic composition of thermal waters result from their mixing with cool groundwater in a shallow aquifer during their ascent to the surface.

---

Résumé Les résultats d'une étude portant sur les caractéristiques hydrochimiques et isotopiques de puits phréatiques et profonds situés dans le champ hydrothermal de Pamukkale, sont décrits de telle manière à éclairer le fonctionnement des circulations hydrologiques. Le champ peut étre divisé en deux sous-systèmes d'eaux souterraines, l'un avec des eaux de sources froides (10–12 °C) de type Ca–HCO3, et les eaux thermales (25–58 °C) riches en CO₂ et de type Ca–HCO₃–SO4. L'occurrence de ces types d'eaux est fermement liée à la morphologie de la région, oùne tectonique intense a engendré des structures en horsts et en grabens. Dés lors deux systèmes hydrogéologiques ont été définis : un système profond, qui est lié à la circulation extensive et profonde des eaux météoritiques dans le système régional d'écoulement, et un système phréatique lié aux écoulements locaux des eaux souterraines à travers les strates sédimentaires. Les eaux météoritiques aux altitudes élevées, percolent jusqu'aux systèmes locaux phréatiques, puis coulent jusqu'aux systèmes géothermaux plus profonds. Durant le cycle de convection profond des zones de recharge jusqu'aux zones de décharge, l'eau froide atteint les zones chaudes liées aux intrusions athenosphériques, provoquant la remontée. Les variations de la composition chimique et isotopique des eaux thermales résultent dans leurs mélanges avec des eaux souterraines froides dans les aquifères phréatiques durant leur remontée jusqu'à la surface.

---

Resumen Se describen los resultados del estudio de las características isotópicas e hidroquímicas de las aguas someras y profundas para obtener un mejor entendimiento de la circulación hidrológica del campo hidrotermal Pamukkale. El campo puede agruparse en dos sub-sistemas de agua subterránea: manantiales de agua fría del tipo Ca–HCO₃ (10–12°C) y aguas termales ricas en CO₂ del tipo Ca–HCO₃–SO₄(25–58°C). El ambiente de estos tipos de aguas se relaciona estrechamente con la morfología de la región donde el tectonismo intenso ha formado estructuras extensionales tipo graben y horst. Se definieron dos sistemas hidrogeológicos: un sistema geotermal profundo que se relaciona con la circulación profunda y extensa de agua meteórica en el sistema regional de flujo y un sistema somero el cual se relaciona con flujo local de agua subterránea a través de estratos sedimentarios. El agua meteórica que cae en altas elevaciones percola al sistema local de agua subterránea en un nivel somero y fluye hacia el sistema geotermal profundo. Durante un ciclo de convección del área de recarga hacia la zona de descarga, el agua fría se calienta a partir de los intrusivos astenosféricos lo que ocasiona que asciendan. Como resultado de la mezcla de las aguas recalentadas, con agua subterránea fría en un acuífero somero durante el ascenso hacia la superficie, se derivan variaciones en la composición química e isotópica de las aguas termales.

     

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.

     

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.