Bootstrap Calibration and Uncertainty Estimation of Downhole NMR Hydraulic Conductivity Estimates in an Unconsolidated Aquifer

Characterization of hydraulic conductivity (K) in aquifers is critical for evaluation, management, and remediation of groundwater resources. While estimates of K have been traditionally obtained using hydraulic tests over discrete intervals in wells, geophysical measurements are emerging as an alternative way to estimate this parameter. Nuclear magnetic resonance (NMR) logging, a technology once largely applied to characterization of deep consolidated rock petroleum reservoirs, is beginning to see use in near‐surface unconsolidated aquifers. Using a well‐known rock physics relationship—the Schlumberger Doll Research (SDR) equation—K and porosity can be estimated from NMR water content and relaxation time. Calibration of SDR parameters is necessary for this transformation because NMR relaxation properties are, in part, a function of magnetic mineralization and pore space geometry, which are locally variable quantities. Here, we present a statistically based method for calibrating SDR parameters that establishes a range for the estimated parameters and simultaneously estimates the uncertainty of the resulting K values. We used co‐located logging NMR and direct K measurements in an unconsolidated fluvial aquifer in Lawrence, Kansas, USA to demonstrate that K can be estimated using logging NMR to a similar level of uncertainty as with traditional direct hydraulic measurements in unconsolidated sediments under field conditions. Results of this study provide a benchmark for future calibrations of NMR to obtain K in unconsolidated sediments and suggest a method for evaluating uncertainty in both K and SDR parameter values.     

NMR Logging to Estimate Hydraulic Conductivity in Unconsolidated Aquifers

Nuclear magnetic resonance (NMR) logging provides a new means of estimating the hydraulic conductivity (K) of unconsolidated aquifers. The estimation of K from the measured NMR parameters can be performed using the Schlumberger‐Doll Research (SDR) equation, which is based on the Kozeny–Carman equation and initially developed for obtaining permeability from NMR logging in petroleum reservoirs. The SDR equation includes empirically determined constants. Decades of research for petroleum applications have resulted in standard values for these constants that can provide accurate estimates of permeability in consolidated formations. The question we asked: Can standard values for the constants be defined for hydrogeologic applications that would yield accurate estimates of K in unconsolidated aquifers? Working at 10 locations at three field sites in Kansas and Washington, USA, we acquired NMR and K data using direct‐push methods over a 10‐ to 20‐m depth interval in the shallow subsurface. Analysis of pairs of NMR and K data revealed that we could dramatically improve K estimates by replacing the standard petroleum constants with new constants, optimal for estimating K in the unconsolidated materials at the field sites. Most significant was the finding that there was little change in the SDR constants between sites. This suggests that we can define a new set of constants that can be used to obtain high resolution, cost‐effective estimates of K from NMR logging in unconsolidated aquifers. This significant result has the potential to change dramatically the approach to determining K for hydrogeologic applications.     

Hydraulic Conductivity Calibration of Logging NMR in a Granite Aquifer,Laramie Range,Wyoming

In granite aquifers, fractures can provide both storage volume and conduits for groundwater. Characterization of fracture hydraulic conductivity (K) in such aquifers is important for predicting flow rate and calibrating models. Nuclear magnetic resonance (NMR) well logging is a method to quickly obtain near-borehole hydraulic conductivity (i.e., K_NMR) at high-vertical resolution. On the other hand, FLUTe flexible liner technology can produce a K profile at comparable resolution but requires a fluid driving force between borehole and formation. For three boreholes completed in a fractured granite, we jointly interpreted logging NMR data and FLUTe K estimates to calibrate an empirical equation for translating borehole NMR data to K estimates. For over 90% of the depth intervals investigated from these boreholes, the estimated K_NMR are within one order of magnitude of K_FLUTe. The empirical parameters obtained from calibrating the NMR data suggest that “intermediate diffusion” and/or “slow diffusion” during the NMR relaxation time may occur in the flowing fractures when hydraulic aperture are sufficiently large. For each borehole, “intermediate diffusion” dominates the relaxation time, therefore assuming “fast diffusion” in the interpretation of NMR data from fractured rock may lead to inaccurate K_NMR estimates. We also compare calibrations using inexpensive slug tests that suggest reliable K_NMR estimates for fractured rock may be achieved using limited calibration against borehole hydraulic measurements.     

Hydraulic Conductivity from Nuclear Magnetic Resonance Logs in Sediments with Elevated Magnetic Susceptibilities

This study examined the application of slim-hole nuclear magnetic resonance (NMR) tools to estimate hydraulic conductivity (K_NMR) in an unconsolidated aquifer that contains a range of grain sizes (silt to gravel) and high and variable magnetic susceptibilities (MS) (10⁻⁴ to 10⁻² SI). A K calibration dataset was acquired at 1-m intervals in three fully screened wells, and compared to K_NMR estimates using the Schlumberger-Doll research (SDR) equation with published empirical constants developed from previous studies in unconsolidated sediments. While K_NMR using published constants was within an order of magnitude of K, the agreement, overprediction, or underprediction of K_NMR varied with the MS distribution in each well. An examination of the effects of MS on NMR data and site-specific empirical constants indicated that the exponent on T_2ML (n-value in the SDR equation, representing the diffusion regime) was found to have the greatest influence on K_NMR estimation accuracy, while NMR porosity did not improve the prediction of K. K_NMR was further improved by integrating an MS log into the NMR analyses. A first approach detrended T_2ML for the effects of MS prior to calculating K_NMR, and a second approach introduced an MS term into the SDR equation. Both were found to produce similar refinements of K_NMR in intervals of elevated MS. This study found that low frequency NMR logging with short echo times shows promise for sites with moderate to elevated MS levels, and recommends a workflow that examines parameter relationships and integrates MS logs into the estimation of K_NMR.     

Effects of variations of river stage and hydraulic conductivity on temporal scaling of groundwater levels: numerical simulations

It was found in previous studies that groundwater levels may fluctuate as a temporal fractal. In this study numerical simulations of groundwater level fluctuations in an unconfined aquifer near a river were conducted to investigate the effects of aquifer heterogeneity and river stage variations on the fractal behavior of the water levels, h(t). Groundwater recharge was taken to be a white-noise process. The aquifer heterogeneity was simulated with a second-order stationary field of hydraulic conductivity (K) with an exponential variogram model. The results showed that groundwater levels fluctuate as a temporal fractal in both homogeneous and heterogeneous aquifers as long as K is less than 10 m/d. Most aquifers may indeed act as a fractal filter which takes a random non-fractal recharge inputs and produces a fractal responses of groundwater level fluctuations. A crossover in temporal scaling of h(t) may appear in more permeable aquifers. Fluctuations of the groundwater level in a homogeneous aquifer are dominated by the recharge process when the river stage is constant or by the river stage variations when the river stage varies in highly permeable aquifers. Heterogeneity plays an important role in the temporal scaling of h(t) in more permeable aquifers: the stronger the heterogeneity, the stronger the temporal scaling of h(t).     

Numerical Simulations of Seasonally Oscillated Groundwater Dynamics in Coastal Confined Aquifers

Studies investigating the effects of inland recharge on coastal groundwater dynamics were carried out typically in unconfined aquifers, with few in confined aquifers. This study focused on the groundwater dynamics in confined aquifers with seasonally sinusoidally fluctuated inland groundwater head and constant sea level by numerical simulations. It is known that the mixing zone (MZ) of saltwater wedge in response to the seasonal oscillations of inland groundwater head swings around the steady-state MZ. However, our simulation results indicate that even the most landward freshwater-saltwater interface over a year is seaward from the steady-state location when the hydraulic conductivity K is ≤10⁻⁴ m/s under certain boundary conditions with given parameter values. That is, seasonal oscillations of inland groundwater head may reduce seawater intrusion in confined coastal aquifers when K ≤ 10⁻⁴ m/s. Sensitivity analysis indicates that for aquifers of K ≤ 10⁻⁴ m/s, the larger the inland head fluctuation amplitude is, the less the seawater intrudes. This is probably due to the reason that the seawater intrusion time decreases with the increase of fluctuation amplitude when K ≤ 10⁻⁴ m/s. Numerical simulations demonstrate that seasonal inland groundwater head oscillations promote the annual averaged recirculated seawater discharge across the seaward boundary.     

Assimilation of historical head data to estimate spatial distributions of stream bed and aquifer hydraulic conductivity fields

Management of water resources in alluvial aquifers relies mainly on understanding interactions between hydraulically connected streams and aquifers. Numerical models that simulate this interaction often are used as decision support tools for water resource management. However, the accuracy of numerical predictions relies heavily on unknown system parameters (e.g., streambed conductivity and aquifer hydraulic conductivity), which are spatially heterogeneous and difficult to measure directly. This paper employs an ensemble smoother to invert groundwater level measurements to jointly estimate spatially varying streambed and alluvial aquifer hydraulic conductivity along a 35.6‐km segment of the South Platte River in Northeastern Colorado. The accuracy of the inversion procedure is evaluated using a synthetic experiment and historical groundwater level measurements, with the latter constituting the novelty of this study in the inversion and validation of high‐resolution fields of streambed and aquifer conductivities. Results show that the estimated streambed conductivity field and aquifer conductivity field produce an acceptable agreement between observed and simulated groundwater levels and stream flow rates. The estimated parameter fields are also used to simulate the spatially varying flow exchange between the alluvial aquifer and the stream, which exhibits high spatial variability along the river reach with a maximum average monthly aquifer gain of about 2.3 m³/day and a maximum average monthly aquifer loss of 2.8 m³/day, per unit area of streambed (m²). These results demonstrate that data assimilation inversion provides a reliable and computationally affordable tool to estimate the spatial variability of streambed and aquifer conductivities at high resolution in real‐world systems.     

Sulfur Hexafluoride and Potassium Bromide as Groundwater Tracers for Managed Aquifer Recharge

Sulfur hexafluoride (SF₆) is an established tracer for use in managed aquifer recharge projects. SF₆ exsolves from groundwater when it encounters trapped air according to Henry's law. This results in its retardation relative to groundwater flow, which can help determine porous media saturation and flow dynamics. SF₆ and the conservative, nonpartitioning tracer, bromide (Br⁻ added as KBr), were introduced to recharge water infiltrated into stacked glacial aquifers in Thurston County, Washington, providing the opportunity to observe SF₆ partitioning. Br⁻, which is assumed to travel at the same velocity as the groundwater, precedes SF₆ at most monitoring wells (MWs). Average groundwater velocity in the unconfined aquifer in the study area ranges from 3.9 to 40 m/d, except in the southwestern corner where it is slower. SF₆ in the shallow aquifer exhibits an average retardation factor of 2.5 ± 3.8, suggesting an air-to-water ratio on the order of 10⁻³ to 10⁻² in the pore space. Notable differences in tracer arrival times at adjacent wells indicate very heterogeneous conductivity. One MW exhibits double peaks in concentrations of both tracers with different degrees of retardation for the first and second peaks. This suggests multiple flowpaths to the well with variable saturation. The confining layer between the upper two aquifers appears to allow intermittent connection between aquifers but serves as an aquitard in most areas. This study demonstrates the utility of SF₆ partitioning for evaluating hydrologic conditions at prospective recharge sites.     

Evaluation of Hydraulic Properties of Aquitards Using Earthquake‐Triggered Groundwater Variation

The hydraulic properties of aquitards are not easily obtained because monitoring wells are usually installed in aquifers for groundwater resources management. Earthquake‐induced crust stress (strain) triggers groundwater level variations over a short period of time in a large area. These groundwater anomalies can be used to investigate aquifer systems. This study uses a poroelastic model to fit the postseismic variations of groundwater level triggered by the Chi‐Chi earthquake to evaluate the hydraulic properties of aquitards in the Jhoushuei River alluvial fan (JRAF), Taiwan. Six of the adopted eight wells with depths of 70 to 130 m showed good agreement with the recovery theory. The mean hydraulic conductivities (K) of the aquifers for the eight wells are 1.62 × 10^?4 to 9.06 × 10^?4 m/s, and the thicknesses are 18.8 to 46.1 m. The thicknesses of the aquitards are 11.3 to 42.0 m. Under the isotropic assumption for K, the estimated values of K for the aquitards are 3.0 × 10^?8 to 2.1 × 10^?6 m/s, corresponding to a silty medium. The results match the values obtained for the geological material of the drilling core and those reported in previous studies. The estimated values were combined with those given in previous studies to determine the distribution of K in the first two aquitards in the JRAF. The distribution patterns of the aquitards reflect the sedimentary environments and fit the geological material. The proposed technique can be used to evaluate the K value of aquitards using inverse methods. The inversion results can be used in hydrogeological analyses, contaminant modeling, and subsidence evaluation.     

Transient Recharge Estimability Through Field‐Scale Groundwater Model Calibration

The estimation of recharge through groundwater model calibration is hampered by the nonuniqueness of recharge and aquifer parameter values. It has been shown recently that the estimability of spatially distributed recharge through calibration of steady‐state models for practical situations (i.e., real‐world, field‐scale aquifer settings) is limited by the need for excessive amounts of hydraulic‐parameter and groundwater‐level data. However, the extent to which temporal recharge variability can be informed through transient model calibration, which involves larger water‐level datasets, but requires the additional consideration of storage parameters, is presently unknown for practical situations. In this study, time‐varying recharge estimates, inferred through calibration of a field‐scale highly parameterized groundwater model, are systematically investigated subject to changes in (1) the degree to which hydraulic parameters including hydraulic conductivity (K) and specific yield (S_y) are constrained, (2) the number of water‐level calibration targets, and (3) the temporal resolution (up to monthly time steps) at which recharge is estimated. The analysis involves the use of a synthetic reality (a reference model) based on a groundwater model of Uley South Basin, South Australia. Identifiability statistics are used to evaluate the ability of recharge and hydraulic parameters to be estimated uniquely. Results show that reasonable estimates of monthly recharge (<30% recharge root‐mean‐squared error) require a considerable amount of transient water‐level data, and that the spatial distribution of K is known. Joint estimation of recharge, S_y and K, however, precludes reasonable inference of recharge and hydraulic parameter values. We conclude that the estimation of temporal recharge variability through calibration may be impractical for real‐world settings.     

Groundwater mixing between different aquifer types in a complex structural setting discerned by elemental and stable isotope geochemistry

Karst aquifers are well known for their intricate stratigraphy and geologic structures, which make groundwater characterization challenging because flowpaths and recharge sources are complex and difficult to evaluate. Geochemical data, collected from ten closely spaced production wells constructed in two karst aquifers (Bangor Limestone (Mb) and Tuscumbia Limestone/Fort Payne Chert (Mftp)) in Trussville, north‐central Alabama, illustrate two distinctive groundwater end‐members: (1) higher major ion, dissolved inorganic carbon, conductivity, alkalinity concentrations, heavier δ¹³C ratios (max: −10.2 ± 0.2‰ Vienna Pee Dee Belemnite (PDB)) and lower residence times (mean: 19.5 ± 2 years, n = 2) of groundwater in the Mb aquifer and (2) lower constituent concentrations, lighter δ¹³C ratios (min: −13.4 ± 0.2‰ PDB) and longer residence times of groundwater (mean: 23.6 ± 2 years, n = 4) in the Mftp aquifer. Summer and fall data and the binary mixing model show aquifer inter‐flow mixing along solution fractures and confirms the distinctive groundwater geochemistry of the two aquifers. Lowering of static water levels over the summer (drawdown from 2 to 5.2 m) leads to more reducing groundwater conditions (lower Eh values) and slightly enriched δ¹⁸O and δD ratios during the fall [δ¹⁸O: −4.8 ± 0.1 to −5.4 ± 0.1‰ Vienna Standard Mean Oceanic Water (VSMOW), n = 9; δD: −25.4 ± 1 to −27.4 ± 1‰ VSMOW, n = 9] when compared with summer season samples (δ¹⁸O: −5.1 ± 0.1 to −5.7 ± 0.1‰ VSMOW, n = 11; δD: −25.0 ± 1 to −30.6 ± 1‰ VSMOW, n = 11). GIS analyses confirm the localized origin of recharge to the investigated aquifers. The combination of GIS, field parameters and geochemistry analyses can be successfully used to identify recharge sources, evaluate groundwater flow and transport pathways and to improve understanding of how groundwater withdrawals impact the sustainability and susceptibility to contamination of karst aquifers. Copyright © 2015 John Wiley & Sons, Ltd.     

Estimation of the percentage of annual groundwater recharge with bomb tritium using a cumulative mass balance method

The bomb tritium (³H) distribution patterns in the aquifer beneath an abandoned landfill at the Canadian Forces Base (CFB) Borden, Ontario, and in a sandy aquifer at Whiteshell Nuclear Research Establishments (WNRE) Pinawa, Manitoba, all in Canada, were delineated in great detail. A sampling and monitoring network of multilevel samplers and bundle piezometers were used. The directions of groundwater flow were established, and the boundary between the tritiated and non-tritiated zones of the two aquifers were closely demarcated. Using a cumulative mass balance method, the³H input mass into the aquifers was compared with the³H mass in groundwater storage to estimate the percentages of annual groundwater recharge from 1953 to 1978. Two recharge calculations for theeffective recharge zone and thetotal recharge area of the aquifers as established from the flownet analysis, and the distributions of dissolved geochemical constitutents showed that theeffective recharge zone calculations gave higher values of 30.6 cm/yr for CFB Borden and 20.1 cm/yr for WNRE while thetotal recharge areas gave lower values of 19.1 and 10.1 cm/yr for the Borden and WNRE aquifers respectively. The two recharge values provide possible minimum and maximum recharge estimates for the two study areas.     

Homotopy Perturbation Method‐Based Analytical Solution for Tide‐Induced Groundwater Fluctuations

The groundwater variations in unconfined aquifers are governed by the nonlinear Boussinesq's equation. Analytical solution for groundwater fluctuations in coastal aquifers under tidal forcing can be solved using perturbation methods. However, the perturbation parameters should be properly selected and predefined for traditional perturbation methods. In this study, a new dimensional, higher‐order analytical solution for groundwater fluctuations is proposed by using the homotopy perturbation method with a virtual perturbation parameter. Parameter‐expansion method is used to remove the secular terms generated during the solution process. The solution does not require any predefined perturbation parameter and valid for higher values of amplitude parameter A/D, where A is the amplitude of the tide and D is the aquifer thickness.     

Relationships between microbial communities and groundwater chemistry in two pristine confined groundwater aquifers in central China

This study explores linkages between the microbial composition and hydrochemical variables of pristine groundwater to identify active redox conditions and processes. Two confined aquifers underlying the city of Qianjiang in the Jianghan Plain in China were selected for this study, having different recharge sources and strong hydrochemical gradients. Typical methods for establishing redox processes according to threshold concentration criteria for geochemical parameters suggest iron or sulphate reduction processes. High‐throughput 16S rRNA sequencing was used to obtain diversity and taxonomic information on microbial communities. Instead of revealing iron‐ and sulphate‐reducing bacteria, salt‐ and alkali‐tolerant bacteria, such as the phylum Firmicutes and the class Gammaproteobacteria, and in particular, the family Bacillaceae, were dominant in the downstream groundwater of the first aquifer that had high ion concentrations caused by the dissolution of calcite and dolomite; meanwhile, the heterotrophic microaerophilic families Comamonadaceae and Rhodocyclaceae prevailed in the upstream groundwater of the first aquifer. Sulphate‐reducing bacteria were extremely abundant in the upstream groundwater of the second aquifer, as the SO₄^2? concentration was especially high. Methanogens and methanotrophs were predominant in the downstream groundwater of the second aquifer even though the concentration of SO₄^2? was much higher than 0.5 mg L^?1. The microbial communities, together with the geochemical parameters, indicated that the upstream region of the first aquifer was suboxic, that Fe(III) and Mn(IV) reductions were not the main redox processes in the downstream groundwater of the first aquifer with high Fe and Mn concentrations, and that the redox processes in the upstream and downstream regions of the second confined aquifer were SO₄^2? reduction and methanogenesis, respectively. This study expands understanding of the linkages between microbial communities and hydrogeochemistry in pristine groundwaters and provides more evidence for identifying active redox conditions and processes.     

Assessing hydrogeologic controls on dynamic groundwater storage using long‐term instrumental records of water table levels

This study analyzes a long‐term regional compilation of water table response to climate variability based on 124 long‐term groundwater wells distributed across New England, USA, screened in a variety of geologic materials. The New England region of the USA is located in a humid‐temperature climate underlain by low‐storage‐fractured metamorphic and crystalline bedrock dissected by north–south trending valleys filled with glacial and post‐glacial valley fill sediments. Uplands are covered by thin glacial till that comprises more than 60% of the total area. Annual and multi‐annual responses of the water table to climate variability are assessed to understand how local hydraulic properties and hydrogeologic setting (located in recharge/discharge region) of the aquifer influence the hydrologic sensitivity of the aquifer system to climate variability. This study documents that upland aquifer systems dominated by thin deposits of surface till comprise ~70% of the active and dynamic storage of the region. Total aquifer storage changes of +5 to ?7 km³ occur over the region during the study interval. The storage response is dominated by thin and low permeability surficial till aquifer that fills and drains on a multi‐annual basis and serves as the main mechanism to deliver water to valley fill aquifers and underlying bedrock aquifers. Whereas the till aquifer system is traditionally neglected as an important storage reservoir, this study highlights the importance of a process‐based understanding of how different landscape hydrogeologic units contribute to the overall hydrologic response of a region.     

Cost-efficient management of coastal aquifers via recharge with treated wastewater and desalination of brackish groundwater: application to the Akrotiri basin and aquifer,Cyprus

Abstract

We investigate the general methodology for an intensive development of coastal aquifers, described in a companion paper, through its application to the management of the Akrotiri aquifer, Cyprus. The Zakaki area of that aquifer, adjacent to Lemessos City, is managed such that it permits a fixed annual agricultural water demand to be met, as well as and a fraction of the water demand of Lemessos, which varies according to available surface water. Effluents of the Lemessos wastewater treatment plant are injected into the aquifer to counteract the seawater intrusion resulting from the increased pumping. The locations of pumping and injection wells are optimized based on least-cost, subject to meeting the demand. This strategy controls sea intrusion so effectively that desalting of only small volumes of slightly brackish groundwater is required over short times, while ～2.3 m³ of groundwater is produced for each 1 m³ of injected treated wastewater. The cost over the 20-year period 2000–2020 of operation is ～40 M€ and the unit production cost of potable water is under 0.2 €/m³. The comparison between the deterministic and stochastic analyses of the groundwater dynamics indicates the former as conservative, i.e. yielding higher groundwater salinity at the well. The Akrotiri case study shows that the proposed aquifer management scheme yields solutions that are preferable to the widely promoted seawater desalination, also considering the revenues from using the treated wastewater for irrigation.

Citation Koussis, A. D., Georgopoulou, E., Kotronarou, A., Mazi, K., Restrepo, P., Destouni, G., Prieto, C., Rodriguez, J. J., Rodriguez-Mirasol, J., Cordero, T., Ioannou, C., Georgiou, A., Schwartz, J. & Zacharias, I. (2010) Cost-efficient management of coastal aquifers via recharge with treated wastewater and desalination of brackish groundwater: application to the Akrotiri basin and aquifer, Cyprus. Hydrol. Sci. J. 55(7), 1234–1245.     

Isotopic tracing for conceptual models of groundwater hydrodynamics in multilayer aquifer systems of central and southern Tunisia

Abstract

New and previously published data sets including stable and radiogenic isotope measurements (¹⁸O, ²H, ³H, ¹³C and ¹⁴C) were used to investigate, conceptualize and compare groundwater hydrodynamics within three major multilayer aquifer systems located in central and southern Tunisia. It has been demonstrated that the investigated aquifer systems contain modern and palaeoclimatic waters. Modern groundwaters, which refer to contemporaneous and post-nuclear recharge waters, are characterized by enriched stable isotope contents, high carbon-14 activities and high to moderate tritium concentration. While palaeoclimatic groundwaters, which refer to Late Pleistocene and Early Holocene recharge waters, are distinguished by their depleted stable isotope contents, low carbon-14 activities and insignificant tritium concentrations. Established conceptual models have elucidated the groundwater hydrodynamics within the studied aquifer systems. They show that groundwater mixing occurs between end-members from the shallow and deep aquifers that migrate by downward and upward leakage towards the intermediate aquifer.

Editor D. Koutsoyiannis; Associate editor S. Faye

Citation Dassi, L. and Tarki, M., 2014. Isotopic tracing for conceptual models of groundwater hydrodynamics in multilayer aquifer systems of central and southern Tunisia. Hydrological Sciences Journal, 59 (6), 1240–1258. http://dx.doi.org/10.1080/02626667.2014.892206     

Groundwater potential estimation using vertical electrical sounding measurements in the semi-arid Khanasser Valley region,Syria

Abstract

Traditional and polar vertical electrical sounding (VES) techniques were used for computing and characterizing the transmissivity of the Quaternary and Palaeogene aquifers in the Khanasser Valley, northern Syria. The VES technique with its different applied approaches is inexpensive, easy and gives rapid results with high spatial resolution. The Dar-Zarrouk parameters were determined to estimate the groundwater potential of the aquifers. The results show the mean transmissivity of the Quaternary aquifer to be 49 m²/d (minimum: 10 m²/d; maximum: 131 m²/d; standard deviation: 32 m²/d), and that of the Palaeogene aquifer is 0.94 m²/d (minimum: 0.31 m²/d; maximum: 1.54 m²/d; standard deviation: 0.33 m²/d).

Editor D. Koutsoyiannis

Citation Asfahani, J., 2013. Groundwater potential estimation using vertical electrical sounding measurements in the semi-arid Khanasser Valley region, Syria. Hydrological Sciences Journal, 58 (2), 468--482.     

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.     

Palaeosol Control of Arsenic Pollution: The Bengal Basin in West Bengal,India

Groundwater in the Bengal Basin is badly polluted by arsenic (As) which adversely affects human health. To provide low‐As groundwater for As mitigation, it was sought across 235 km² of central West Bengal, in the western part of the basin. By drilling 76 boreholes and chemical analysis of 535 water wells, groundwater with <10 µg/L As in shallow aquifers was found under one‐third of a study area. The groundwater is in late Pleistocene palaeo‐interfluvial aquifers of weathered brown sand that are capped by a palaeosol of red clay. The aquifers form two N‐S trending lineaments that are bounded on the east by an As‐polluted deep palaeo‐channel aquifer and separated by a shallower palaeo‐channel aquifer. The depth to the top of the palaeo‐interfluvial aquifers is mostly between 35 and 38 m below ground level (mbgl). The palaeo‐interfluvial aquifers are overlain by shallow palaeo‐channel aquifers of gray sand in which groundwater is usually As‐polluted. The palaeosol now protects the palaeo‐interfluvial aquifers from downward migration of As‐polluted groundwater in overlying shallow palaeo‐channel aquifers. The depth to the palaeo‐interfluvial aquifers of 35 to 38 mbgl makes the cost of their exploitation affordable to most of the rural poor of West Bengal, who can install a well cheaply to depths up to 60 mbgl. The protection against pollution afforded by the palaeosol means that the palaeo‐interfluvial aquifers will provide a long‐term source of low‐As groundwater to mitigate As pollution of groundwater in the shallower, heavily used, palaeo‐channel aquifers. This option for mitigation is cheap to employ and instantly available.