Spatial measurements of stream baseflow,a relevant method for aquifer characterization and permeability evaluation. Application to a hard‐rock aquifer,the Oman ophiolite

The structure, functioning and hydrodynamic properties of aquifers can be determined from an analysis of the spatial variability of baseflow in the streams with which they are associated. Such analyses are based on simple low‐cost measurements. Through interpreting the hydrological profiles (Q = f(A)) it is possible to locate the aquifer(s) linked to the stream network and to determine the type of interrelated flow, i.e. whether the stream drains or feeds the aquifer. Using an analytical solution developed for situations with a positive linear relationship, i.e. where the baseflow increases linearly with increasing catchment size, it is also possible to estimate the permeability of the aquifer(s) concerned at catchment scale. Applied to the hard‐rock aquifers of the Oman ophiolite, this method shows that the ‘gabbro’ aquifer is more permeable than the ‘peridotite’ aquifer. As a consequence the streams drain the peridotites and ‘leak’ into the gabbro. The hydrological profiles within the peridotite are linear and positive, and indicate homogeneity in the hydrodynamic properties of these formations at the kilometre scale. The permeability of the peridotite is estimated at 5 · 10^?7 to 5 · 10^?8 m/s. Copyright © 2004 John Wiley & Sons, Ltd.     

A methodology for regionalizing 3‐D effective porosity at watershed scale in crystalline aquifers

An innovative approach for regionalizing the 3‐D effective porosity field is presented and applied to two large, overexploited, and deeply weathered crystalline aquifers located in southern India. The method derives from earlier work on regionalizing a 2‐D effective porosity field in that part of an aquifer where the water table fluctuates, which is now extended over the entire aquifer using a 3‐D approach. A method based on geological and geophysical surveys has also been developed for mapping the weathering profile layers (saprolite and fractured layers). The method for regionalizing 3‐D effective porosity combines water table fluctuation and groundwater budget techniques at various cell sizes with the use of satellite‐based data (for groundwater abstraction), the structure of the weathering profile, and geostatistical techniques. The approach is presented in detail for the Kudaliar watershed (983 km²) and tested on the 730 km² Anantapur watershed. At watershed scale, the effective porosity of the aquifer ranges from 0.5% to 2% in Kudaliar and between 0.3% and 1% in Anantapur, which agrees with earlier works. Results show that (a) depending on the geology and on the structure of the weathering profile, the vertical distribution of effective porosity can be very different and that the fractured layers in crystalline aquifers are not necessarily characterized by a rapid decrease in effective porosity and (b) that the lateral variations in effective porosity can be larger than the vertical ones. These variations suggest that within a same weathering profile, the density of open fractures and/or degree of weathering in the fractured zone may significantly vary from a place to another. The proposed method provides information on the spatial distribution of effective porosity that is of prime interest in terms of flux and contaminant transport in crystalline aquifers. Implications for mapping groundwater storage and scarcity are also discussed, which should help in improving groundwater resource management strategies.     

Seismic wave velocity and anisotropy of serpentinized peridotite in the Oman ophiolite

Shallow seismic measurements in harzburgite from the Oman ophiolite performed in a zone where the maximum horizontal anisotropy is expected (vertical foliation and horizontal lineation) point to a dominant dependence of seismic properties on fracturing.

Optical microscopy studies show that microcracks are guided by the serpentine (lizardite) penetrative network oriented subparallel to the harzburgite foliation and subperpendicular to the mineral lineation, and that serpentine (lizardite) vein filling has a maximum concentration of (001) planes parallel to the veins walls. The calculated elastic properties of the oriented alteration veins filled with serpentine in an anisotropic matrix formed by oriented crystals of olivine and orthopyroxene are compared with seismic velocities measured on hand specimens.

Laboratory ultrasonic data indicate that open microcracks are closed at 100 MPa pressure, e.g. (J. Geophys. Res. 65, (1960) 1083) and (Proc. ODP Sci. Results Leg 118, (1990) 227). Above this pressure, laboratory measurements and modeling show that P-compressional and S-shear wave velocities are mainly controlled by the mineral preferred orientation. Veins sealed with serpentine are effective in slightly lowering P and S velocities and increasing anisotropy. The penetrative lizardite network does not affect directly the geometry of seismic anisotropy, but contributes indirectly in the fact that this network controls the microcrack orientations.

Comparison between seismic measurements of peridotite and gabbro in the same conditions suggest that P- and S-waves anisotropies are a possible discriminating factor between the two lithologies in the suboceanic lithosphere.     

A conceptual hydrogeological model of ophiolite hard-rock aquifers in Oman based on a multiscale and a multidisciplinary approach

Ophiolites are found all over the world: from the Alps to the Himalayas, in Cuba, Papua-New Guinea, New Caledonia, Newfoundland, etc. They are composed of hard rocks—basalt, dolerite, gabbro and peridotite, which are formed at the mid-oceanic ridges, with specific ridge-related tectonic fracturing and intense hydrothermal alteration. Their geological and thus their hydrogeological properties differ from those of both granite or classical gabbro and classical basaltic lava. A conceptual hydrogeological model of these hard-rock aquifers was developed based on the convergent results of a multidisciplinary approach at several spatial scales, from rock-sample (centimetre) to catchment (kilometre), on well-preserved ophiolite rocks in Oman. In ophiolite rocks, groundwater circulation takes place mostly in the fissured near-surface horizon (50 m thick), and, to a lesser degree, in the tectonic fractures. Hydrograph analysis (Water Resour Res 34:233–240, 1977), interpretation of numerous pumping tests using both classical Theis and dual porosity models [Water Resour Res 32:2733–2745, 1996; Comput Geosci J (in press)], and mercury porosity and hydraulic conductivity lab-measurements support the aquifer parameter estimates. The hydraulic conductivity K of the fissured horizon is estimated at 10^–5 to 10^–6 m/s for gabbro and dolerite, and 10^–7 m/s for peridotite. The storage coefficient S of the peridotite aquifer is estimated at 10^–3 and appears to be controlled mainly by microcracks (20 to 100 m wide). Tectonic fractures in the ophiolite have similar hydrodynamic properties regardless of lithology (10^–1<T<10^–4 m²/s and 10^–1<S<10^–3) though the probability of obtaining productive wells is two to three times greater in gabbro and dolerite than in peridotite. Some of the tectonic fractures produce small hydrothermal, hyperalkaline springs in the peridotite. The water budget and hydrochemistry of the Oman ophiolite are characterized and support the conceptual hydrogeological model. Despite low annual rainfall, a relatively low hydraulic conductivity and a significant storage coefficient explain why most of the streams in peridotite are perennial.

---

Résumé Les ophiolites se trouvent partout dans le monde : des Alpes à lHymalaya, à Cuba, Papouasie Nouvelle-Guinée, Nouvelle Calédonie, etc. Les ophiolites sont composées de roches dures basaltiques, de dolérites, gabbros et péridotites, qui sont formées aux rides mi-océaniques, avec, en relation avec les rides, une fracturation tectonique et une altération hydrothermale intense. Leurs propriétés géologiques et hydrogéologiques diffèrent de celles des granites ou des gabbros classiques ou des laves basaltiques classiques. Un modèle hydrogéologique conceptuel des aquifères de hard-rock a été développé sur les résultats convergents dune approche multidisciplinaire à différentes échelles spatiales, de léchantillon de roche de quelques centimètres à la taille du bassin versant en kilomètres, des ophiolites bien conservées dOman. Dans les roches ophiolitiques, la circulation des eaux souterraines est surtout localisée dans les fissures dun horizon de surface de plus dune cinquantaine de mètres et à un degré moindre dans les fractures tectoniques. Lanalyse des hydrographes (Water Resour Res 34:233–240, 1977), linterprétation des nombreux essais de pompage utilisant des modèles classiques de Theis et des modèles à double porosité (Water Resour Res 32:2733–2745, 1996; Comput Geosci J (Soumis), la porosité au mercure et la conductivité hydraulique calculée en laboratoire supportent lestimation des paramètres de laquifère. La conductivité hydraulique K de lhorizon fissuré est estimée entre 10^–5 et 10^–6 m/s pour les gabbros et les dolérites, et 10^–7 m/s pour les péridotites. Le coefficient demmagasinement S de laquifère de péridotite est estimé à 10^–3 et semble être contrôlé par les micro-fractures (20 à 100 m de largeur). Les fractures tectoniques dans les ophiolites possèdent des propriétés hydrodynamiques en regard de la lithologie (10-1<T<10-4 m²/s et 10-1/S/10-3) bien que la probabilité dobtenir des puits productifs est trois plus grand dans le grabbro et la dolérite que la péridétite. Quelques unes des fractures tectoniques produisent de petites sources hydrothermales hyperalcalines dans la péridotite. Le bilan hydrologique et hydrochimique de lophiolite dOman caractérise le modèle hydrogéologique conceptuel. Malgré des faibles précipitations annuelles, une relativement basse conductivité hydraulique et un coefficient demmagasinement significatif expliquent pourquoi la plus part des cours deau sur les péridotites ne sont jamais taris.

---

Resumen Las rocas ofiolitas se encuentran en muchas partes del mundo desde los Alpes a los Himilaya, Newfoundland, Nueva Caledonia etc. Son rocas duras compuestas de basalto, dolerita, gabro y peridotita que son formadas por los cordones medio-oceánicos y por fracturación relacionada con los cordones tectónicos específicos y alteración hidrotermal intensa. La geología y los propiedades hidrogeológicas varian entre las de granito y gabro clásico y lava de basalto clásico. Se desarrolló un modelo conceptual de los acuíferos de las rocas duras Basado en los resultados convergentes de un enfoce multidisciplinario a varias escalas a partir de muestras de roca (centímetro) a cuenca kilómetro) enfocado en las rocas ofilíticas bien preservados en Oman. En las rocas ofiolíticas la circulación de agua ocurre principalmente en el horizonte de fisuras cercanas a la superficie (>50 m de espesor) y en un grado menos importante en fracturas tectónicas. Los estimados de los parámetros de estos acuíferos se basan en análisis de hidrógrafos (Water Resour Res 34:233–240, 1977), en la interpretación de pruebas de bombeo usando los modelos clásicos de Theis y de porosidad doble (Hamm y Bidaux 1996; Lods y Gouze, [en prensa]), y porosidad de mecurio y conductividad hidraúlica medida en el laboratorio. La conductividad hidraúlica K del horizonte fisurado se estima en 10E-5 a 10E-6 m-s para gabro y dolerita y 10E-7 para peridotita. El coeficiente de abastecimiento S del acuífero peridotita se estima en 10–3 y parece ser controlado principalmente por microfracturas (20 to 100 m de ancho). Las fracturas tectónicas en la ofiolita tienen propiedades hidrodinámicas similares independientes de la litología (10–1<T<10–4 m²/s and 10–1<S<10–3) pero la probabilidad de obtener pozos productivos es dos o tres veces mayor en gabro que en dolerita o peridotita. Algunas de las fracturas tectónicas producen pequeñas descargas de agua hidrotermal hiperalcalina en la peridotita. El presupuesto e hidroquímica de las ofilitas de Omán se caracteriza y apoya en el modelo conceptual hidrogeológico. La baja lluvia anual, una conductividad hidraúlica relativamente baja y un coeficiente de abastecimiento importante explican porqué la mejoría de los drenajes en peridotita es peremne.

     

Optimized groundwater management applied to a large-yield (10 M m3/year) well field in France

Hydrogeology Journal - The strong impact of population increase and the effects of climate change on drinking-water resources mean that it is essential to optimize the management of groundwater....     

An efficient methodology for estimating irrigation return flow coefficients of irrigated crops at watershed and seasonal scale

Irrigation return flow coefficients, i.e. the ratio between the quantity of water returned from the cultivated area to the groundwater system and the amount of abstraction, vary by more than 50% for rice cultivation using standing water irrigation to 0% in the case of drip irrigation technique. This component of the groundwater budget plays an important role, particularly in intensively irrigated areas. Thus, to avoid any inaccurate aquifer budgeting, modelling and consequently any erroneous watershed management, this component needs to be accurately assessed for a particular time‐step (e.g. weekly, seasonally) onto the studied area. The present paper proposes a cost‐effective and useful methodology for assessing irrigation return flow coefficients (C_f = irrigation return flow/pumping flow) based on (i) basic crops field survey and meteorological data and (ii) the use of a simple hydraulic model that combines both water balance technique and unsaturated/saturated flow theory. An attempt to estimate the uncertainty of irrigation return flow coefficient estimates based on the uncertainty introduced by the pumping and the natural spatial variability of the soil characteristics is also proposed. Results have been compared to real field conditions and allow us to (i) estimate the uncertainty and (ii) validate and demonstrate the robustness of the applied methodology. The proposed methodology allows relatively good estimates of the irrigation return flow coefficients at watershed and seasonal scale. The irrigation return flow coefficients are calculated as: 51 ± 8% in rainy season (Kharif) and 48 ± 4% in summer (Rabi) for rice; 26 ± 11% in rainy season and 24 ± 4% in summer for vegetables; 13 ± 8% in rainy season and 11 ± 3% in summer for flowers. These results were found to be consistent with the existing literature. Copyright © 2007 John Wiley & Sons, Ltd.     

Estimating aquifer thickness using multiple pumping tests

A method to estimate aquifer thickness and hydraulic conductivity has been developed, consisting of multiple pumping tests. The method requires short-duration pumping cycles on an unconfined aquifer with significant seasonal water-table fluctuations. The interpretation of several pumping tests at a site in India under various initial conditions provides information on the change in hydrodynamic parameters in relation to the initial water-table level. The transmissivity linearly decreases compared with the initial water level, suggesting a homogeneous distribution of hydraulic conductivity with depth. The hydraulic conductivity is estimated from the slope of this linear relationship. The extrapolation of the relationship between transmissivity and water level provides an estimate of the aquifer thickness that is in good agreement with geophysical investigations. The hydraulically active part of the aquifer is located in both the shallow weathered and the underlying densely fractured zones of the crystalline basement. However, no significant relationship is found between the aquifer storage coefficient and initial water level. This new method contributes to filling the methodological gap between single pumping tests and hydraulic tomography, in providing information on the variation of the global transmissivity according to depth. It can be applied to any unconfined aquifer experiencing large seasonal water-table fluctuations and short pumping cycles.     

Improving hydrogeological understanding through well-test interpretation by diagnostic plot and modelling: a case study in an alluvial aquifer in France

Hydrogeology Journal - The study of groundwater resources using pumping test data is usually carried out with the Theis solution, which enables the hydraulic parameters of porous aquifers such as...     

High-yielding aquifers in crystalline basement: insights about the role of fault zones,exemplified by Armorican Massif,France

While groundwater constitutes a crucial resource in many crystalline-rock regions worldwide, well-yield conditions are highly variable and barely understood. Nevertheless, it is well known that fault zones may have the capacity to ensure sustainable yield in crystalline media, but there are only a few and disparate examples in the literature that describe high-yield conditions related to fault zones in crystalline rock basements. By investigating structural and hydraulic properties of remarkable yielding sites identified in the Armorican Massif, western France, this study discusses the main factors that may explain such exceptional hydrogeological properties. Twenty-three sites, identified through analysis of databases available for the region, are investigated. Results show that: (1) the highly transmissive fractures are related to fault zones which ensure the main water inflow in the pumped wells; (2) the probability of intersecting such transmissive fault zones does not vary significantly with depth, at least within the range investigated in this study (0–200 m); and (3) high yield is mainly controlled by the structural features of the fault zones, in particular the fault dip and the presence of a connected storage reservoir. Conceptual models that summarize the hydrological properties of high-yield groundwater resources related to fault zones in crystalline basement are shown and discussed.