Electrical-hydraulic conductivity model for a weatheredfractured aquifer system of Olbanita,Lower Baringo Basin,Kenya Rift

Groundwater yield in the Kenya Rift is highly unsustainable owing to geological variability.In this study,field hydraulic characterization was performed by using geoelectric approaches.The relations between electrical-hydraulic(eh)conductivities were modeled hypothetically and calibrated empirically.Correlations were based on the stochastic models and field-scale hydraulic parameters were contingent on pore-level parameters.By considering variation in pore-size distributions over eh conduction interval,the relations were scaled-up for use at aquifer-level.Material-level electrical conductivities were determined by using Vertical Electrical Survey and hydraulic conductivities by analyzing aquifer tests of eight boreholes in the Olbanita aquifer located in Kenya rift.VES datasets were inverted by using the computer code IP2Win.The main result is that ln T=0.537(ln Fa)+3.695;the positive gradient indicating eh conduction through pore-surface networks and a proxy of weathered and clayey materials.An inverse(1/F-K)correlation is observed.Hydraulic parameters determined using such approaches may possibly contribute significantly towards sustainable yield management and planning of groundwater resources.     

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.     

Evaluation of unconfined aquifer parameters from flow to partially penetrating wells in Tailan River basin, China

Effective evaluation, management and abstraction of groundwater resources of any aquifer require accurate and reliable estimates of its hydraulic parameters. This study, therefore, looks at the determination of hydraulic parameters of an unconfined aquifer using both analytical and numerical approaches. A long-duration pumping test data obtained from an unconfined aquifer system within the Tailan River basin in Xinjiang Autonomous Region in the northwest of China is used, in this study, to investigate the best method for estimating the parameters of the aquifer. The pumping test was conducted by pumping from a radial collector well and measuring the response in nine observation wells; all the wells used in the test were partially penetrating. Using two well-known tools, namely AquiferTest and MODFLOW, as an aid for the analytical and numerical approaches, respectively, the parameters of the aquifer were determined and their outputs compared. The estimated horizontal hydraulic conductivity, vertical hydraulic conductivity, and specific yield for the analytical approach are 38.1–50.30 m/day, 3.02–9.05 m/day and 0.204–0.339, respectively, while the corresponding numerical estimates are 20.50–35.24 m/day, 0.10–3.40 m/day, and 0.27–0.31, respectively. Comparing the two, the numerical estimates were found to be more representative of the aquifer in the study area since it simulated the groundwater flow conditions of the pumping test in the aquifer system better than the analytical solution.     

Using lithologic modeling techniques for aquifer characterization and groundwater flow modeling of the Sohag area, Egypt

Three dimensional lithologic modeling techniques have been used for detailed characterization and groundwater flow modeling of the Quaternary aquifer system of the Sohag area, Egypt. Well log data were used for building the lithologic model using RockWorks. A groundwater flow model, facilitated by MODFLOW 2000, was built using results of the lithologic model. The obtained lithologic model honored the well log data and revealed a complex sedimentary system, which is mainly composed of six lithologic categories: clay, clay and sand, fine sand, coarse sand, sand and gravel, and gravel. Inter-fingering and presence of lenses are the main characteristics of the sedimentary basin in the study area. A wide range of hydraulic conductivities is present, which vary spatially and control the groundwater flow. Heterogeneity of the aquifer system is spatially represented where different hydraulic conductivities are found in the different directions. Sandy layers tend to be connected. Hydraulic continuity is represented by inter-fingering and connection of sandy materials within the aquifer system.     

Characterisation of conduction phenomena in soils at the particle-scale: Finite element analyses in conjunction with synthetic 3D imaging

This paper describes an innovative method to characterise conduction parameters in geomaterials at the particle-scale. The technique is exemplified using 3D synthetic grain packing generated with discrete element approaches. This creates a geo-mechanically viable user-defined 3D granular image through which the particle skeleton and the corresponding pore network are constructed. Images are then imported into the finite element analyses to solve the governing equations of hydraulic and thermal conduction. Navier–Stokes equation is uniquely upscaled to Darcy’s law to assess hydraulic conductivity in soils, while a similar approach implements the Fourier equation to evaluate thermal conduction through grain chains and pore network. High performance computing is used to meet demanding numerical calculations of 3D meshed geometries. Packing density (i.e., porosity) and inter-particle contact areas are explored as variables to highlight the effects of pore volume and inter-particle contact condition in hydraulic and thermal conduction. This emerging technique allows not only characterising the macro-scale behaviour of conduction phenomena in soils but also quantifying and visualising the preferential and local conduction behaviour at the particle-scale. Laboratory measurements of hydraulic and thermal conductivities support numerically obtained results and validate the viability of the new methods used herein. This study introduces an alternative way to determine physical parameters of soils using emerging technology of rigorous numerical simulations in conjunction with 3D images, and to enable fundamental observation of particle-scale mechanisms of macro-scale manifestation.     

Estimating hydraulic properties from tidal attenuation in the Northern Guam Lens Aquifer, territory of Guam, USA

Tidal-signal attenuations are analyzed to compute hydraulic diffusivities and estimate regional hydraulic conductivities of the Northern Guam Lens Aquifer, Territory of Guam (Pacific Ocean), USA. The results indicate a significant tidal-damping effect at the coastal boundary. Hydraulic diffusivities computed using a simple analytical solution for well responses to tidal forcings near the periphery of the island are two orders of magnitude lower than for wells in the island’s interior. Based on assigned specific yields of ～0.01–0.4, estimated hydraulic conductivities are ～20–800 m/day for peripheral wells, and ～2,000–90,000 m/day for interior wells. The lower conductivity of the peripheral rocks relative to the interior rocks may best be explained by the effects of karst evolution: (1) dissolutional enhancement of horizontal hydraulic conductivity in the interior; (2) case-hardening and concurrent reduction of local hydraulic conductivity in the cliffs and steeply inclined rocks of the periphery; and (3) the stronger influence of higher-conductivity regional-scale features in the interior relative to the periphery. A simple numerical model calibrated with measured water levels and tidal response estimates values for hydraulic conductivity and storage parameters consistent with the analytical solution. The study demonstrates how simple techniques can be useful for characterizing regional aquifer properties.     

Fracture flow modelling based on satellite images of the Wajid Sandstone, Saudi Arabia

Large-scale geological features have been identified by satellite imagery and global positioning system data in the Wajid Sandstone in Saudi Arabia. The main objective is to evaluate the importance of fractures for the overall flow behaviour in this fractured rock aquifer and to estimate in-situ hydraulic apertures. Data on fractures and lineaments were available for three outcrops. By applying a “cut-out” routine on the fracture endpoint data of these fracture trace windows, three deterministic discrete fracture networks (DFN), with an area of 100 m?×?100 m, could be generated. These were used to simulate the fracture flow and to determine the hydraulic conductivity tensors. Using additional data on hydraulic pumping tests and matrix conductivities, in-situ hydraulic apertures could be determined. Average in-situ hydraulic apertures range from 1,300 to 1,700 µm. Observations from the field support these results. In addition, a hydraulic conductivity ratio between the matrix and fracture system was used to identify the contribution of the DFN to the overall fluid transport. A ratio of 10.4 was determined, which indicates that the effective flow behaviour in the Wajid Sandstone aquifer is not entirely dominated by the fracture system, though evidently strongly controlled by it.     

The typology of Irish hard-rock aquifers based on an integrated hydrogeological and geophysical approach

Groundwater flow in hard-rock aquifers is strongly controlled by the characteristics and distribution of structural heterogeneity. A methodology for catchment-scale characterisation is presented, based on the integration of complementary, multi-scale hydrogeological, geophysical and geological approaches. This was applied to three contrasting catchments underlain by metamorphic rocks in the northern parts of Ireland (Republic of Ireland and Northern Ireland, UK). Cross-validated surface and borehole geophysical investigations confirm the discontinuous overburden, lithological compartmentalisation of the bedrock and important spatial variations of the weathered bedrock profiles at macro-scale. Fracture analysis suggests that the recent (Alpine) tectonic fabric exerts strong control on the internal aquifer structure at meso-scale, which is likely to impact on the anisotropy of aquifer properties. The combination of the interpretation of depth-specific hydraulic-test data with the structural information provided by geophysical tests allows characterisation of the hydrodynamic properties of the identified aquifer units. Regionally, the distribution of hydraulic conductivities can be described by inverse power laws specific to the aquifer litho-type. Observed groundwater flow directions reflect this multi-scale structure. The proposed integrated approach applies widely available investigative tools to identify key dominant structures controlling groundwater flow, characterising the aquifer type for each catchment and resolving the spatial distribution of relevant aquifer units and associated hydrodynamic parameters.     

Estimation of regional-scale groundwater flow properties in the Bengal Basin of India and Bangladesh

Quantitative evaluation of management strategies for long-term supply of safe groundwater for drinking from the Bengal Basin aquifer (India and Bangladesh) requires estimation of the large-scale hydrogeologic properties that control flow. The Basin consists of a stratified, heterogeneous sequence of sediments with aquitards that may separate aquifers locally, but evidence does not support existence of regional confining units. Considered at a large scale, the Basin may be aptly described as a single aquifer with higher horizontal than vertical hydraulic conductivity. Though data are sparse, estimation of regional-scale aquifer properties is possible from three existing data types: hydraulic heads, ¹⁴C concentrations, and driller logs. Estimation is carried out with inverse groundwater modeling using measured heads, by model calibration using estimated water ages based on ¹⁴C, and by statistical analysis of driller logs. Similar estimates of hydraulic conductivities result from all three data types; a resulting typical value of vertical anisotropy (ratio of horizontal to vertical conductivity) is 10⁴. The vertical anisotropy estimate is supported by simulation of flow through geostatistical fields consistent with driller log data. The high estimated value of vertical anisotropy in hydraulic conductivity indicates that even disconnected aquitards, if numerous, can strongly control the equivalent hydraulic parameters of an aquifer system.     

含水层水力参数的尺度效应研究现状

主要论述了国内外目前对弥散率尺度效应以及渗透系数等夺数的空间变异性研究的现状和今后的研究趋势。文中提出了含水层水力参数的尺度效应问题，并指出庄重视地下水驱动力场变化对参数的影响研究。     

Aquifer vulnerability using recharge,depth to groundwater,soil type and slope to classify the vadose zone (Molototsi and Middle Letaba quaternary catchments,Limpopo Province,South Africa)

The aquifer vulnerability of the Molototsi (B81G) and Middle Letaba (B82D) quaternary catchments was assessed to determine the influence of the vadose zone on the groundwater regime. The aquifer vulnerability was assessed by developing a new method, which evaluates the vadose zone as a pathway for pollutants using the following four parameters: recharge, depth to water table, soil type (saturated vertical hydraulic conductivity) and slope (RDSS). Recharge was estimated using the chloride mass balance method, and the depth to the water table was measured in the field using dipmeter. The seepage behaviour (soil type) was determined as hydraulic conductivity from in situ infiltration and percolation testing (SABS 0252-2:1993 and double ring infiltrometer). The slopes were determined with the digital elevation method using ArcGIS software. The four parameters were overlaid using weighted sum, weighted overlay and raster calculator to produce the vulnerability map. Different weightings were attributed in the methods and the best selected. The results obtained indicated high vulnerability on the lower and upper parts of both catchments. The benefits of the method described are (a) the easy quantification of the parameters through fairly simple methods and (b) the exclusion of arbitrary index values.     

Numerical groundwater-flow modeling to evaluate potential effects of pumping and recharge: implications for sustainable groundwater management in the Mahanadi delta region,India

Process-based groundwater models are useful to understand complex aquifer systems and make predictions about their response to hydrological changes. A conceptual model for evaluating responses to environmental changes is presented, considering the hydrogeologic framework, flow processes, aquifer hydraulic properties, boundary conditions, and sources and sinks of the groundwater system. Based on this conceptual model, a quasi-three-dimensional transient groundwater flow model was designed using MODFLOW to simulate the groundwater system of Mahanadi River delta, eastern India. The model was constructed in the context of an upper unconfined aquifer and lower confined aquifer, separated by an aquitard. Hydraulic heads of 13 shallow wells and 11 deep wells were used to calibrate transient groundwater conditions during 1997–2006, followed by validation (2007–2011). The aquifer and aquitard hydraulic properties were obtained by pumping tests and were calibrated along with the rainfall recharge. The statistical and graphical performance indicators suggested a reasonably good simulation of groundwater flow over the study area. Sensitivity analysis revealed that groundwater level is most sensitive to the hydraulic conductivities of both the aquifers, followed by vertical hydraulic conductivity of the confining layer. The calibrated model was then employed to explore groundwater-flow dynamics in response to changes in pumping and recharge conditions. The simulation results indicate that pumping has a substantial effect on the confined aquifer flow regime as compared to the unconfined aquifer. The results and insights from this study have important implications for other regional groundwater modeling studies, especially in multi-layered aquifer systems.     

Surface electrical resistivity tomography and hydrogeological characterization to constrain groundwater flow modeling in an agricultural field site near Ferrara (Italy)

A portion of an unconfined alluvial aquifer located in the Padana Plain (Italy) was characterized following an integrated hydro-geophysical approach. Initially an electrical resistivity tomography (ERT) survey was employed to localize the boundaries of a modest paleo-channel body and to design the installation of a groundwater monitoring network. Multilevel slug-tests were performed to estimate the aquifer’s saturated hydraulic conductivities. Determined permeability values together with electrical resistivity data were correlated. The correlation resulted in a site specific bi-logarithmic linear relationship. Based on this relationship, punctually determined hydraulic conductivities were spatially extended over the studied flow domain. Finally, continuously measured piezometric heads were used to calibrate a 3D flow model. Sensitivity analysis was performed to confirm the reliability of the reconstructed permeability field, as well as, to assess the minimum number of direct measurements needed to safely characterize the selected aquifer portion. The integration of the ERT survey results with the classical hydrogeological tests can be conveniently applied to constrain the permeability field during flow model calibration. Although the applicability of the determined relationship is site specific, the followed procedure is useful especially when there is a need to optimize the available resources and in case of small-scale pilot studies.     

Contaminant Spill in a Sandbox with Non-Gaussian Conductivities: Simultaneous Identification by the Restart Normal-Score Ensemble Kalman Filter

Mathematical Geosciences - The joint identification of the parameters defining a contaminant source and the heterogeneous distribution of the hydraulic conductivities of the aquifer where the...     

Analysis of pumping test data for determining unconfined-aquifer parameters: Composite analysis or not?

Recently, composite analysis (CA), which simultaneously analyzes all drawdown data from multiple observation wells, has been applied to determine the hydraulic parameters of an unconfined aquifer. Moench (1994) claimed that the value of specific yield (S _y) determined from non-composite analysis (nonCA) is sometimes unrealistically low as compared with that obtained by water-balance calculation, and results from CA are better representative of aquifer properties than those from nonCA. To examine the validity of this assertion, the drawdown data from a pumping test conducted at Cape Cod, Massachusetts, USA, were analyzed using both nonCA and CA methods. The results show that the mean estimates of hydraulic conductivity and S _y determined from CA are close to those determined from nonCA. In some cases the analysis based on CA also results in low estimates of S _y as compared with those determined based on nonCA. A hypothetical case study is presented, which examines the effect of measurement errors on the estimated parameters. The results indicate that the CA method also gives poorer estimates of S _y than the nonCA method if the pumping test data contain measurement errors. Moench AF (1994) Specific yield as determined by type-curve analysis of aquifer-test data. Ground Water, 32(6):949–957.     

Pattern Recognition in a Bimodal Aquifer Using the Normal-Score Ensemble Kalman Filter

The ensemble Kalman filter (EnKF) is now widely used in diverse disciplines to estimate model parameters and update model states by integrating observed data. The EnKF is known to perform optimally only for multi-Gaussian distributed states and parameters. A new approach, the normal-score EnKF (NS-EnKF), has been recently proposed to handle complex aquifers with non-Gaussian distributed parameters. In this work, we aim at investigating the capacity of the NS-EnKF to identify patterns in the spatial distribution of the model parameters (hydraulic conductivities) by assimilating dynamic observations in the absence of direct measurements of the parameters themselves. In some situations, hydraulic conductivity measurements (hard data) may not be available, which requires the estimation of conductivities from indirect observations, such as piezometric heads. We show how the NS-EnKF is capable of retrieving the bimodal nature of a synthetic aquifer solely from piezometric head data. By comparison with a more standard implementation of the EnKF, the NS-EnKF gives better results with regard to histogram preservation, uncertainty assessment, and transport predictions.     

Evaluating the salinity distribution of a shallow coastal aquifer by vertical multielectrode profiling (Denmark)

A monitoring system, including five groups of piezometers and five vertical multielectrode profiling probes (VMEP), has been installed in an aquifer beneath a coastal dune in Denmark. In order to assess the salinity distribution within the aquifer, geoelectrical data were gathered in March, June and September 2008, by measuring a dipole-dipole and gradient array using multielectrode profiling. Interpretation of the processed resistivity data was performed by regularized inversion using a one-dimensional, horizontally layered model of formation resistivity. The standard deviation on estimated layer log-resistivity was 0.01–0.03. By estimating two parameters of a power function, observed fluid conductivities derived from samples of porewater were related to corresponding estimated formation resistivities. The conductivity profiles correlate with a winter situation in March with high sea level, active recharge and significant wave activity, causing increased hydraulic heads, a thicker freshwater lens and salt water overlying freshwater close to the sea. In June, the thickness of the freshwater lens is reduced due to less recharge and prevailing offshore winds, imposing density-stable conditions and a sharper transition between fresh and brackish water. During the autumn, aquifer recharge is enhanced and hydraulic heads increase, resulting in a thicker freshwater lens.     

Developing conceptual hydrogeological model for Potsdam sandstones in southwestern Quebec,Canada

A hydrogeological study was conducted in Potsdam sandstones on the international border between Canada (Quebec) and the USA (New York). Two sandstone formations, arkose and conglomerate (base) and well-cemented quartz arenite (upper), underlie the study area and form the major regional aquifer unit. Glacial till, littoral sand and gravel, and marine silt and clay discontinuously overlie the aquifer. In both sandstone formations, sub-horizontal bedding planes are ubiquitous and display significant hydraulic conductivities that are orders of magnitude more permeable than the intact rock matrix. Aquifer tests demonstrate that the two formations have similar bulk hydrologic properties, with average hydraulic conductivities ranging from 2?×?10^?5 to 4?×?10^?5 m/s. However, due to their different lithologic and structural characteristics, these two sandstones impose rather different controls on groundwater flow patterns in the study area. Flow is sustained through two types of fracture networks: sub-horizontal, laterally extensive fractures in the basal sandstone, where hydraulic connectivity is very good horizontally but very poor vertically and each of the water-bearing bedding planes can be considered as a separate planar two-dimensional aquifer unit; and the more fractured and vertically jointed system found in the upper sandstone that promotes a more dispersed, three-dimensional movement of groundwater.     

What can be learned from sequential multi-well pumping tests in fracture-karst media? A case study in Zhangji, China

A fracture-karst aquifer is a karst aquifer with a fractured rock matrix, and its parameters are difficult to determine. Two sequential pumping tests in a fracture-karst aquifer system at the Zhangji well field of China are considered, one carried out before (in 2000) and one after (in 2005) the operation of a pumping station in the well field (2003–2005). The sequential tests serve multiple purposes. First, they provide a cross check of the parameters obtained. Second, they can be used to assess the effect of long-term groundwater exploitation of the aquifer. A three-dimensional finite-element transient flow model has been developed to simulate groundwater flow at the site. Generally good agreement has been found between the simulated and observed hydraulic heads for both tests. The hydraulic parameters obtained from the 2005 test are generally consistent with their counterparts from the 2000 test. However, a small but steady increase of hydraulic conductivities from 2000 to 2005 at the site has been observed. A 10-year prediction of groundwater resources has been made and indicates that the well field can accommodate the proposed 8.0?×?10⁴ m³/day exploitation rate under relative drought conditions without causing a steady decline of groundwater levels.