Controls on groundwater flow in the Bengal Basin of India and Bangladesh: regional modeling analysis

Groundwater for domestic and irrigation purposes is produced primarily from shallow parts of the Bengal Basin aquifer system (India and Bangladesh), which contains high concentrations of dissolved arsenic (exceeding worldwide drinking water standards), though deeper groundwater is generally low in arsenic. An essential first step for determining sustainable management of the deep groundwater resource is identification of hydrogeologic controls on flow and quantification of basin-scale groundwater flow patterns. Results from groundwater modeling, in which the Bengal Basin aquifer system is represented as a single aquifer with higher horizontal than vertical hydraulic conductivity, indicate that this anisotropy is the primary hydrogeologic control on the natural flowpath lengths. Despite extremely low hydraulic gradients due to minimal topographic relief, anisotropy implies large-scale (tens to hundreds of kilometers) flow at depth. Other hydrogeologic factors, including lateral and vertical changes in hydraulic conductivity, have minor effects on overall flow patterns. However, because natural hydraulic gradients are low, the impact of pumping on groundwater flow is overwhelming; modeling indicates that pumping has substantially changed the shallow groundwater budget and flowpaths from predevelopment conditions.     

含水层层状非均质对地下水流系统的影响

区域尺度上含水层非均质具有复杂的结构性和随机性,难以准确刻画,造成非均质对区域地下水流系统的影响机制研究不够深入。本文以鄂尔多斯盆地白垩系地下水流系统为研究实例,选择典型剖面,采用剖面二维随机数值模拟方法,通过对比不同非均质刻画方法下地下水流场的变化,探讨含水层层状非均质对地下水流系统的影响机制。结果显示,均质条件下模型各向异性（含水层水平和垂向渗透系数比值Kh/Kv）取值为1000时,地下水流场与实际条件较为接近;非均质条件下,渗透系数方差取值0.91,水平相关长度取值5000 m,Kh/Kv取值150时,接近实际条件。研究表明,在大尺度地下水流模拟研究中,采用水平相关长度、渗透系数方差和各向异性值三个变量生成的随机场能很好地刻画含水层的层状非均质特征及其对水流系统的影响控制作用。由于含水层不同尺度层状非均质的叠加效应,采用均质各向异性介质等效概化含水层层状非均质性会造成等效各向异性值偏大失真的效应。     

Ambient vertical flow in long-screen wells: a case study in the Fontainebleau Sands Aquifer (France)

A tritium (³H) profile was constructed in a long-screened well (LSW) of the Fontainebleau Sands Aquifer (France), and the data were combined with temperature logs to gain insight into the potential effects of the ambient vertical flow (AVF) of water through the well on the natural aquifer stratification. AVF is commonly taken into account in wells located in fracture aquifers or intercepting two different aquifers with distinct hydraulic heads. However, due to the vertical hydraulic gradient of the flow lines intercepted by wells, AVF of groundwater is a common process within any type of aquifer. The detection of ³H in the deeper parts of the studied well (approximate depth 50 m), where ³H-free groundwater is expected, indicates that shallow young water is being transported downwards through the well itself. The temperature logs show a nearly zero gradient with depth, far below the mean geothermal gradient in sedimentary basins. The results show that the age distribution of groundwater samples might be biased in relation to the age distribution in the surroundings of the well. The use of environmental tracers to investigate aquifer properties, particularly in LSWs, is then limited by the effects of the AVF of water that naturally occurs through the well.     

Stochastic groundwater modeling of a sedimentary aquifer: evaluation of the impacts of abstraction scenarios under conditions of reduced recharge

A transient finite difference groundwater flow model has been calibrated for the Nasia sub-catchment of the White Volta Basin. This model has been validated through a stochastic parameter randomization process and used to evaluate the impacts of groundwater abstraction scenarios on resource sustainability in the basin. A total of 1500 equally likely model realizations of the same terrain based on 1500 equally likely combinations of the data of the key aquifer input parameters were calibrated and used for the scenario analysis. This was done to evaluate model non-uniqueness arising from uncertainties in the key aquifer parameters especially hydraulic conductivity and recharge by comparing the realizations and statistically determining the degree to which they differ from each other. Parameter standard deviations, computed from the calibrated data of the key parameters of hydraulic conductivity and recharge, were used as a yardstick for evaluating model non-uniqueness. All model realizations suggest horizontal hydraulic conductivity estimates in the range of 0.03–78.4 m/day, although over 70 % of the area has values in the range of 0.03–14 m/day. Low standard deviations of the horizontal hydraulic conductivity estimates from the 1500 solutions suggest that this range adequately reflects the properties of the material in the terrain. Lateral groundwater inflows and outflows appear to constitute significant components of the groundwater budgets in the terrain, although estimated direct vertical recharge from precipitation amounts to about 7 % of annual precipitation. High potential for groundwater development has been suggested in the simulations, corroborating earlier estimates of groundwater recharge. Simulation of groundwater abstraction scenarios suggests that the domain can sustain abstraction rates of up to 200 % of the current estimated abstraction rates of 12,960 m³/day under the current recharge rates. Decreasing groundwater recharge by 10 % over a 20-year period will not significantly alter the results of this abstraction scenario. However, increasing abstraction rates by 300 % over the period with decreasing recharge by 10 % will lead to drastic drawdowns in the hydraulic head over the entire terrain by up to 6 m and could cause reversals of flow in most parts of the terrain.     

Constraining flow paths of saline groundwater at basin margins using hydrochemistry and environmental isotopes: Lake Cooper,Murray Basin,Australia

Solutes in saline groundwater (total dissolved solids up to 37 000 mg/L) in the Lake Cooper region in the southern margin of the Riverine Province of the Murray Basin are derived by evapotranspiration of rainfall with minor silicate, carbonate and halite dissolution. The distribution of hydraulic heads, salinity, percentage modern carbon (pmc) contents, and Cl/Br ratios imply that the groundwater system is complex with vertical flow superimposed on lateral flow away from the basin margins. Similarities in major ion composition, stable (O, H, and C) isotope, and ⁸⁷Sr/⁸⁶Sr ratios between groundwater from the shallower Shepparton Formation and the deeper Calivil – Renmark aquifer also imply that these aquifers are hydraulically interconnected. Groundwater in the deeper Calivil – Renmark aquifer in the Lake Cooper region has residence times of up to 25 000 years, implying that pre-land-clearing recharge rates were <1 mm/y. As in other regions of the Murray Basin, the low recharge rates account for the occurrence of high-salinity groundwater. Shallow (<20 m) groundwater yields exclusively modern ¹⁴C ages and shows a greater influence of evaporation over transpiration. Both these observations reflect the rise of the regional water-table following land clearing over the last 200 years and a subsequent increase in recharge to 10 – 20 mm/y. The rise of the regional water-table also has increased vertical and horizontal hydraulic gradients that may ultimately lead to the export of salt from the Lake Cooper embayment into the adjacent fresher groundwater resources.     

Modeling of groundwater flow for Mujib aquifer, Jordan

Jordan is an arid country with very limited water resources. Groundwater is the main source for its water supply. Mujib aquifer is located in the central part of Jordan and is a major source of drinking water for Amman, Madaba and Karak cities. High abstraction rates from Mujib aquifer during the previous years lead to a major decline in water levels and deterioration in groundwater quality. Therefore, proper groundwater management of Mujib aquifer is necessary; and groundwater flow modeling is essential for proper management. For this purpose, Modflow was used to build a groundwater flow model to simulate the behavior of the flow system under different stresses. The model was calibrated for steady state condition by matching observed and simulated initial head counter lines. Drawdown data for the period 1985–1995 were used to calibrate the transient model by matching simulated drawdown with the observed one. Then, the transient model was validated by using drawdown data for the period 1996–2002. The results of the calibrated model showed that the horizontal hydraulic conductivity of the B2/A7 aquifer ranges between 0.001 and 40m/d. Calibrated specific yield ranges from 0.0001 to 0.15. The water balance for the steady state condition of Mujib aquifer indicated that the total annual direct recharge is 20.4 × 10⁶m³, the total annual inflow is 13.0 × 10⁶ m³, springs discharge is 15.3 × 10⁶ m³, and total annual outflow is 18.7 × 10⁶ m³. Different scenarios were considered to predict aquifer system response under different conditions. The results of the sensitivity analysis show that the model is highly sensitive to horizontal hydraulic conductivity and anisotropy and with lower level to the recharge rates. Also the model is sensitive to specific yield     

Field characterization and data integration to define the hydraulic heterogeneity of a shallow granular aquifer at a sub-watershed scale

Providing a sound basis for aquifer management or remediation requires that hydrogeological investigations carried out to understand groundwater flow and contaminant transport be based on representative data that capture the heterogeneous spatial distribution of aquifer hydraulic properties. This paper describes a general workflow allowing the characterization of the heterogeneity of the hydraulic properties of granular aquifers at an intermediate scale of a few km². The workflow involves characterization and data integration steps that were applied on a 12-km² study area encompassing a decommissioned landfill emitting a leachate plume and its main surface water receptors. The sediments composing the aquifer were deposited in a littoral–sublittoral environment and show evidence of small-scale transitional heterogeneities. Cone penetrometer tests (CPT) combined with soil moisture and electrical resistivity (SMR) measurements were thus used to identify and characterize spatial heterogeneities in hydraulic properties over the study area. Site-specific statistical relationships were needed to infer hydrofacies units and to estimate hydraulic properties from high-resolution CPT/SMR soundings distributed all over the study area. A learning machine approach was used due to the complex statistical relationships between colocated hydraulic and CPT/SMR data covering the full range of aquifer materials. Application of this workflow allowed the identification of hydrofacies units and the estimation of horizontal hydraulic conductivity, vertical hydraulic conductivity and porosity over the study area. The paper describes and discusses data acquisition and integration methodologies that can be adapted to different field situations, while making the aquifer characterization process more time-efficient and less labor-intensive.     

Analysis of Steady-State Flow To Multiscreened Wells Under Natural Vertical Flow Conditions

Water levels measured at multiscreened wells in unconfined aquifers may not coincide, in general, with the elevation of the water table. The presence of vertical gradients (as often is the case in recharge areas) or the existence of confining layers may cause the water levels to differ from local hydraulic heads in the aquifer. In these cases, a misinterpretation of water levels may lead to the erroneous conclusion that observed drawdowns are provoked by overpumping. In this paper, we analyze the effect that a natural vertical gradient has on water levels in wells screened over their entire saturated thickness. As one would expect, it is observed that, even without pumping, the water level in the wells lies below the water table. Type curves relating the steady-state drawdown to the vertical gradient and to the hydraulic conductivity anisotropy are presented. These curves were obtained using a groundwater flow numerical model (FREESURF: Neuman and Witherspoon, 1970). The theoretical results are checked with field data from deep wells in the detrital Madrid aquifer. In this particular aquifer, it is observed that the effect of vertical gradients is important both in terms of drawdowns and flow rates.

---

     

The hydraulic conductivity field and groundwater flow in the unconfined aquifer system of the Keta Strip,Ghana

This study investigates the hydraulic conductivity field and the groundwater flow pattern as predicted by a calibrated steady state groundwater flow model for the Keta Strip, southeastern Ghana. The hydraulic conductivity field is an important parameter in evaluating aquifer properties in space, and in general basin-wide groundwater resources evaluation and management. This study finds that the general hydraulic conductivity of the unconsolidated unconfined aquifer system of the Keta Strip ranges between 2 m/d and 20 m/d, with an average of 15 m/d. The spatial variation in horizontal hydraulic conductivity appears to take the trend in the variations in the nature of the material in space. Calibrated groundwater recharge suggests that 6.9–34% of annual precipitation recharges the shallow aquifer system. This amount of recharge is significant and suggests high fortunes in terms of groundwater resources development for agriculture and industrial activities in the area. A spatial distribution of groundwater recharge from precipitation is presented in this study. The spatial pattern appears to take the form of the distribution in horizontal hydraulic conductivity, and suggests that the vertical hydraulic conductivity takes the same pattern of spatial variation as the horizontal hydraulic conductivity. This is consistent with observations in other areas. The resulting groundwater flow is dominated by local flow systems as the unconfined system is quite shallow. A general northeast – southwest flow pattern has been observed in the study area.     

Evaluation of heterogeneity and field-scale groundwater flow regime in a leaky till aquitard

Recent work in southern Ontario, Canada, demonstrates anomalously high vertical groundwater flow velocities (>1 m/year) through a thick (as much as 60 m), sandy silt till aquitard (Northern till), previously assumed to be of very low permeability (hydraulic conductivity <10^–10 m/s). Rapid recharge is attributed to the presence of fractures and sedimentary heterogeneities within the till, but the field-scale flow regime is poorly understood. This study identifies the nature of physical groundwater pathways through the till and provides estimates of the associated groundwater fluxes. The aquitard groundwater flow system is characterized by integrating details of the outcrop and subsurface sedimentary characteristics of the till with field-based hydrogeologic investigation and numerical modeling. Outcrop and subsurface data identify a composite internal aquitard stratigraphy consisting of tabular till beds (till elements) separated by laterally continuous sheet-like sands and gravels (interbeds) and boulder pavements. Individual till elements contain sedimentary heterogeneities, including discontinuous sand and gravel lenses, vertical sand dikes, and zones of horizontal and vertical fractures. Hydrogeologic field investigations indicate a three-layer aquitard flow system, consisting of upper and lower zones of more hydraulically active and heterogeneous till separated by a middle unit of relatively lower hydraulic conductivity. Groundwater pathways and fluxes in the till were evaluated using a two-dimensional aquitard/aquifer flow model which indicates a step-wise flow mechanism whereby groundwater moves alternately downward along vertical pathways (fractures, sedimentary dikes) and laterally along horizontal sand interbeds within the till. This model is consistent with observed hydraulic-head and isotope profiles, and the presence of tritiated pore waters at various depths throughout the till. Simulations suggest that a bulk aquitard vertical hydraulic conductivity on the order of 1×10^–9 m/s is required to reproduce observed hydraulic-head and tritium profiles. Electronic Publication     

Re-evaluating the US Geological Survey’s pumping tests (1967) in the Punjab region of Pakistan for use in groundwater studies

In 1967, the US Geological Survey (USGS) published the results of 141 pumping tests carried out throughout the Pakistani Punjab to establish representative hydraulic parameters of its large aquifer. Many authors have since concluded that the USGS had over-estimated the horizontal hydraulic conductivity (k _r) by 25–100 %, leaving vertical anisotropy and aquifer depth unresolved. No test wells have ever been drilled below 450 m to reach the base of the aquifer, although petroleum explorations mention depths between 1,500 and 4,500 m. After comparison and re-evaluation of all related papers, this study concludes that the USGS interpretation was correct, that its hydraulic values still stand without change, and that the USGS’s applied distance drawdown interpretation is valid to prevent influence of partial penetration on the results. This study also uniquely resolved vertical anisotropy and aquifer thickness by using early- and late-time drawdowns separately and proper scaling of the coordinates, which has often been omitted. With appropriate scaling, all interpretations match the data. The representative hydraulic aquifer values are: k _r?=?65 m/d, vertical anisotropy k _r/k _z?=?25 and aquifer depth 500–1,500 m. The conclusion is that these values can be used, at least as first estimates, for groundwater studies in the Pakistani Punjab.     

Correlation of hydraulic conductivity with stratigraphy in a fractured-dolomite aquifer, northeastern Wisconsin, USA

This study integrated surface and subsurface stratigraphic data with geophysical logs and hydrogeologic data in order to characterize the hydraulic properties of the Silurian dolomite in northeastern Wisconsin. Silurian stratigraphy consists of predictable alternations of characteristic facies associations. A vertical profile of hydraulic conductivity, obtained from short-interval packer tests in a core hole that penetrates a majority of the Silurian section, indicates that hydraulic conductivity ranges over five orders of magnitude (10^–1 to 10^–6 cm/s). Matrix conductivity is generally low and varies with texture; the finer-grained restricted-marine and transitional facies being less permeable than the coarser-grained open-marine facies. High-conductivity values are generally associated with bedding-plane fractures, and fracture frequency is greater in the restricted-marine facies. High-flow features in 16 wells were identified using fluid-temperature/resistivity and heat-pulse flowmeter logs. Natural-gamma logs were used to identify the stratigraphic position of flow features in each well and to correlate flow features to specific stratigraphic horizons. By combining stratigraphic, geophysical, and hydrogeologic data, 14 high-permeability zones within the Silurian aquifer have been identified and correlated in wells more than 16 km apart. These zones parallel bedding, appear most pronounced at contacts between contrasting lithologies, and are most abundant in restricted-marine lithologies. Electronic Publication     

Structural controls on groundwater flow in a fractured bedrock aquifer underlying an agricultural region of northwestern New Brunswick,Canada

A hydrogeological study was conducted in northwestern New Brunswick, Canada, to improve the predictability of fracture-dominated groundwater flow within folded bedrock composed of fine-grained turbidites. Borehole televiewer logging and outcrop mapping, integrated with hydraulic packer tests revealed enhanced hydraulic conductivity associated with northeasterly striking bedding-plane fractures formed during folding and flexural slip. These fractures impart azimuthal anisotropy to the aquifer because of moderately dipping fold limbs. High-angle fractures form a well-developed non-stratabound network, comprising two open fracture sets striking NNE parallel to the current direction of principal stress, and WNW parallel to the direction of principal stress that dominated during the Acadian orogeny. The subset of fractures showing significant oxidation, deemed most important to the groundwater flow system, is dominated by bedding-plane and high-angle fractures striking near-parallel to the maximum principal stress direction, resulting in extensional opening and enhanced hydraulic conductivities. An equivalent porous media model, incorporating anisotropy and varying hydraulic conductivity with depth, indicates that horizontal flow dominates the aquifer with relatively minor exchange between different model layers. These findings have implications for understanding flow directions in the Black Brook Watershed and elsewhere in the Matapédia Basin where fractures formed under similar stress conditions.     

Migration of arsenic in multi-aquifer system of southern Bengal Basin: analysis via numerical modeling

A heterogeneous anisotropic steady-state groundwater flow model for the multi-aquifer system of a part of southern Bengal Basin shows that human intervention has changed the natural groundwater flow system. At present, the shallow groundwater flow is restricted within the aquifer, with very short travel time of tens of years and vertical path length. The deep aquifer is fed by surface water or rainwater from distant locations with travel time of thousands of years and has no hydraulic connection with the arsenic-rich shallow aquifer. Numerical simulations indicate that the future pumping of deep groundwater is not likely to drive in arsenic from the shallow aquifer. Therefore, new wells may be installed in the deep aquifer. High pumping of shallow unpolluted aquifer consisting of brown sand will drive in groundwater containing organic matter from the post-Last Glacial Maximum aquifer-aquitard system. The organic matter drives reduction of manganese oxides at strip interfaces between palaeo-channel and palaeo-interfluve. After the completion of manganese reduction, FeOOH reduction may take place in the marginal palaeo-interfluvial aquifer and release sorbed arsenic. Arsenic then moves into the interior of palaeo-interfluvial aquifer polluting its fresh groundwater. Arsenic migration rates ranges between 0.21 and 6.3 and 1.3 × 10^?2 and 0.4 m/year in horizontal and vertical directions, respectively. Therefore, palaeo-interfluvial aquifer will remain arsenic-free for hundreds to thousands of years to supply safe drinking water.     

Aquifer hydraulic conductivity estimation from surface geoelectrical measurements for Krauthausen test site, Germany

The physical explanation is presented for two types of theoretical equations which encompass a direct and an inverse relationship between the hydraulic conductivity and electrical resistivity of an aquifer. Using the theory of direct current (DC) resistivity prospecting, it is shown that a direct relation exists if the bulk of the current flow is in the vertical direction, whereas the inverse relation exists if the current is in the horizontal direction. Consequently, the relation in question is direct in the case of a relatively conducting basement and inverse in the case of a comparatively resistive basement. Then, the aquifer hydraulic conductivity for the hydrogeologically well-studied area of Krauthausen in Germany is estimated by two different techniques using surface geoelectrical data. The estimated values are compared with those conventionally determined by a pump test experiment and it can be concluded that the estimation is better whenever geoelectrical data are used.     

Role of high-elevation groundwater flows in the hydrogeology of the Cimino volcano (central Italy) and possibilities to capture drinking water in a geogenically contaminated environment

Origin, yield and quality of the groundwater flows at high elevation in the Cimino volcano (central Italy) were examined. In this area, groundwater is geogenically contaminated by arsenic and fluoride, yet supplies drinking water for approximately 170,000 inhabitants. The origin of the high-elevation groundwater flows is strictly related to vertical and horizontal variability of the rock types (lava flows, lava domes and ignimbrite) in an area of limited size. In some cases, groundwater circuits are related to perched aquifers above noncontinuous aquitards; in other cases, they are due to flows in the highly fractured dome carapace, limited at the bottom by a low-permeability dome core. The high-elevation groundwater outflow represents about 30% of the total recharge of Cimino’s hydrogeological system, which has been estimated at 9.8 L/s/km². Bicarbonate alkaline-earth, cold, neutral waters with low salinity, and notably with low arsenic and fluoride content, distinguish the high-elevation groundwaters from those of the basal aquifer. Given the quantity and quality of these resources, approaches in the capture and management of groundwater in this hydrogeological environment should be reconsidered. Appropriate tapping methods such as horizontal drains, could more efficiently capture the high-elevation groundwater resources, as opposed to the waters currently pumped from the basal aquifer which often require dearsenification treatments.

     

Evaluation of methods for determination of hydraulic properties in an aquifer–aquitard system hydrologically connected to a river

Field and laboratory methods have been used to determine the hydraulic properties in a multiple-layer aquifer–aquitard system that is hydrologically connected to a river. First, hypothetical pumping tests in aquifer–aquitard systems were performed to evaluate the feasibility of MODFLOW-PEST in determining these parameters. Sensitivity analyses showed that: the horizontal hydraulic conductivity in the aquifer has the highest composite sensitivity; the vertical hydraulic conductivity has higher composite sensitivity than the horizontal hydraulic conductivity in the aquitard; and a partial penetration pumping well in an aquifer layer can improve the quality of the estimated parameters. This inverse approach was then used to analyze a pumping-recovery test conducted near the Platte River in southeastern Nebraska, USA. The hydraulic conductivities and specific yield were calculated for the aquitard and aquifer. The direct-push technique was used to generate sediment columns; permeameter tests on these columns produced the vertical hydraulic conductivities that are compatible with those obtained from the pumping-recovery test. Thus, the combination of the direct-push technique with permeameter tests provides a new method for estimation of vertical hydraulic conductivity. The hydraulic conductivity, determined from grain-size analysis, is smaller than the horizontal one but larger than the vertical one determined by the pumping-recovery test.     

应用沉降和水位数据计算上海地区弱透水层的参数

弱透水层参数的确定是准三维或三维地下水流数学模型应用时的难点,而应力-应变图解法是一种估计弱透水层参数的有效方法。在介绍该法的基本原理之后,结合上海沉降和水位观测资料,求出了上海地区弱透水层的弹性贮水率、非弹性贮水率和垂向渗透系数,为上海地面沉降的数值模拟提供了较好的初始参数值。对于弱透水层,抽水试验的方法是不可取的,如果又没有条件取土样作室内试验时,若有沉降资料和相应的水位资料,则应力-应变图解法是确定近似参数值的一种比较好的方法。     

Large-scale hydraulic conductivity distribution in an unconfined carbonate aquifer using the ocean tidal propagation

The hydraulic conductivity of an unconfined carbonate aquifer at the uplifted atoll of Minami-Daito, Japan, was evaluated by a combination of cross-spectral analysis, analytical solution, and density-dependent groundwater modeling based on observed groundwater levels in 15 wells and at sea level. The island area was divided into 10 subregions based on island morphology and on inland propagation of ocean tides. The hydraulic conductivity was obtained for each subregion using analytical solutions based on phase lags of M₂ constituents of ocean tides at each well by assuming two aquifer thicknesses (300 and 1,800 m) and two effective porosities (0.1 and 0.3). The density-dependent groundwater model evaluated the hydraulic conductivity of the subregions by reproducing observed groundwater levels. The hydraulic conductivity in the subregions was estimated as 3.46?×?10^?3 to 6.35?×?10^?2 m/s for aquifer thickness of 300 m and effective porosity of 0.1, and as 1.73?×?10^?3 to 3.17?×?10^?2 m/s for aquifer thickness of 1,800 m and the effective porosity of 0.3. It was higher in southern and northern areas, and higher in interior lowland than in the western and eastern areas. Fissures and dolomite distributions on the island control differences of the omnidirectional ocean tidal propagation and cause these differences in hydraulic conductivity. The method used for this study may also be applicable to other small islands that have few or no data for hydraulic conductivity.