Saturated and unsaturated River Nile/groundwater aquifer interaction systems in the Nile Valley, Egypt

Surface water bodies interact with underlying aquifer systems, creating a complex flow system and flow paths. In general, a surface water body may be classified as gaining, losing, or flow through on the basis of its interaction with the surrounding aquifer. In the Nile Valley, the quaternary aquifer system is in a direct hydraulic interaction with the River Nile, canals, and drains. In this study, a regional numerical model was developed and used to evaluate the interaction between surface water bodies and the quaternary aquifer system in the Nile Valley. The solution is considered for a quasi three-dimensional, steady-state groundwater flow. The model used simulates the interaction between surface water bodies and groundwater for saturated and unsaturated flow conditions. In addition, a hydrodynamic model was used to simulate different extreme (high and low) scenarios for Nile surface water levels along the distance between Old Aswan Dam and Delta Barrages. Model calibration shows close results, and the model was used to simulate surface water levels. Results indicate that the Nile River acts as a drain for the quaternary aquifer (gaining water from the aquifer), although in the reaches upstream of the main barrages, the Nile loses the water, recharging the aquifer. All other main canals are recharging the aquifer system. The seepage rate depends mainly on the difference in piezometric head between the aquifer system and surface water bodies, as well as the hydraulic conductance of the base layer sediments of the surface water body. The model was used to evaluate the regional water balance for the Nile Valley and to estimate the surface water bodies' gains and losses.     

Elevated arsenic in deeper groundwater of the western Bengal basin, India: Extent and controls from regional to local scale

The deeper groundwater (depending on definition) of the Bengal basin (Ganges-Brahmaputra delta) has long been considered as an alternate, safe drinking-water source in areas with As-enrichment in near-surface groundwater. The present study provides the first collective discussion on extent and controls of elevated As in deeper groundwater of a regional study area in the western part of the Bengal basin. Deeper groundwater is defined here as non-brackish, potable (Cl⁻ ? 250 mg/L) groundwater available at the maximum accessed depth (∼80-300 m). The extent of elevated As in deeper groundwater in the study area seems to be largely controlled by the aquifer-aquitard framework. Arsenic-enriched deeper groundwater is mostly encountered north of 22.75°N latitude, where an unconfined to semi-confined aquifer consisting of Holocene- to early Neogene-age gray sand dominates the hydrostratigraphy to 300 m depth below land surface. Aquifer sediments are not abnormally enriched in As at any depth, but sediment and water chemistry are conducive to As mobilization in both shallow and deeper parts of the aquifer(s). The biogeochemical triggers are influenced by complex redox disequilibria. Results of numerical modeling and profiles of environmental tracers at a local-scale study site suggest that deeper groundwater abstraction can draw As-enriched water to 150 m depth within a few decades, synchronous with the advent of wide-scale irrigational pumping in West Bengal (India).     

Geologic controls on porosity development in the Maynardville Limestone, Oak Ridge, Tennessee

 Understanding the geologic controls of porosity development and their relationship to the karst aquifer system in the Cambrian Maynardville Limestone is important in determining possible contaminant transport pathways and provides essential data for hydrologic models within the Oak Ridge Reservation of east Tennessee. In the Maynardville Limestone, several important factors control porosity development: (1) lithologic controls on secondary microporosity and mesoporosity are related to dissolution of evaporite minerals and dedolomitization in supratidal facies; (2) depth below the ground surface controls the formation of karst features because the most active portion of the groundwater system is at shallow depths, and karst features are rare below ≈35 m; and (3) structural controls are related to solution enlargement of fractures and faults. Received: 21 May 1996 · Accepted: 30 August 1996     

Hydrogeochemical stratification of the unconfined Samobor aquifer (Zagreb,Croatia)

In the western part of the city of Zagreb, Croatia, hydrogeochemical and isotopic investigations of the Samobor aquifer were carried out with the aim of determining the differences in hydrogeochemical characteristics at increasing aquifer depths. The aquifer comprises 40-m thick gravelly–sandy deposits, with lenses and interlayers of silt and clay. The analyses have proven that with increasing aquifer depth, there are decreases in groundwater temperature and the values of electrical conductivity and increases in the sodium, iron and manganese concentrations. The δ¹³C distribution shows an evident increase in biogenic carbon concentrations with increasing aquifer depth. The measured specific ¹⁴C activities showed that the deeper part of the aquifer is characterized by slow water exchange, while the shallower part is influenced by current recharge, although the pumping wells located on the well-field downstream penetrate the aquifer fully. A direct exchange of water from the Sava River and groundwater occurs in the near vicinity of the river. This exchange weakens further away, while the difference in hydrogeochemical characteristics between the Sava River water and groundwater increases.     

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.     

Depth of the base of the Jackson aquifer, based on geophysical exploration, southern Jackson Hole, Wyoming, USA

 A geophysical survey was conducted to determine the depth of the base of the water-table aquifer in the southern part of Jackson Hole, Wyoming, USA. Audio-magnetotellurics (AMT) measurements at 77 sites in the study area yielded electrical-resistivity logs of the subsurface, and these were used to infer lithologic changes with depth. A 100–600 ohm-m geoelectric layer, designated the Jackson aquifer, was used to represent surficial saturated, unconsolidated deposits of Quaternary age. The median depth of the base of the Jackson aquifer is estimated to be 200 ft (61 m), based on 62 sites that had sufficient resistivity data. AMT-measured values were kriged to predict the depth to the base of the aquifer throughout the southern part of Jackson Hole. Contour maps of the kriging predictions indicate that the depth of the base of the Jackson aquifer is shallow in the central part of the study area near the East and West Gros Ventre Buttes, deeper in the west near the Teton fault system, and shallow at the southern edge of Jackson Hole. Predicted, contoured depths range from 100 ft (30 m) in the south, near the confluences of Spring Creek and Flat Creek with the Snake River, to 700 ft (210 m) in the west, near the town of Wilson, Wyoming. Received, May 1997 · Revised, February 1998 · Accepted, April 1998     

Alluvial aquifer recharge enhanced by a natural dam: feasibility assessment based on multidisciplinary characterization (Khuzestan,Southwest Iran)

In the Kushkak Valley (Khuzestan, Southwest Iran) an anticlinal structure has partially impounded an ephemeral stream. This natural impounded area has been chosen for an artificial recharge site due to its current rate of recharge, capability to store water and favorable situation for the construction of man-made barriers to stream flow. The aquifer to be recharged is the Kushkak unconfined aquifer which consists of medium to coarse-grained alluvial deposits that overlie consolidated conglomerate rock. In this semi-arid area with infrequent relatively heavy falls of rain, alluvial aquifer recharge can be an important process that sustains shallow, over-exploited groundwater bodies. In this investigation a multidisciplinary approach including: hydrometerological studies, and a detailed hydrogeochemical survey, have been carried out. Other essential prerequisite parameters for the scheme were also taken into account to determine the suitability of this location for groundwater artificial recharge. The assessment has brought out that (1) the proposed reservoir will conserve a major part of the water being lost, (2) annual runoff of about 0.27 MCM can be injected into the aquifer through recharge from impounded water, (3) hydrochemical data from surface water and from the Kushkak aquifer water demonstrates that dilution and change in compositional trend in the groundwater proximal to the impounded alluvial bed areas would be expected based on the infiltration capacity of this site, and (4) cost–benefit ratio of the project is 1:2 and it is assumed to recover the investment within six years.     

The relationship between the water table and the surface flow of a losing stream, lower Medano Creek, Great Sand Dunes National Monument, Colorado

 Proposed groundwater withdrawals in the San Luis Valley of Colorado may lower the water table in Great Sand Dunes National Monument. In response, the National Park Service initiated a study that has produced a generalized conceptual model of the hydrologic system in order to assess whether a lowering of the water table might decrease the surface flow of lower Medano Creek. Based upon information obtained during the drilling of several boreholes, there appear to be five important hydrostratigraphic units underlying lower Medano Creek within the upper 30 m of the ground surface: 1. a perched aquifer overlying an aquitard located between about 5 and 6 m below the ground surface; 2. the aquitard itself; 3. an unconfined aquifer located between the upper and lower aquitards; 4. an aquitard located between about 27 and 29 m below the ground surface; and 5. a confined underlying the lower aquitard. Because the areal extent of the aquitards cannot be determined from the borehole data, a detailed conceptual model of the hydrogeologic system underlying lower Medano Creek cannot be developed. However, a generalized conceptual model can be envisioned that consists of a complex system of interlayered aquifers and leaky aquitards, with each aquifer having a unique hydraulic head. Water levels in the perched aquifer rise rapidly to their annual maximum levels in response to the arrival of the flow terminus of Medano Creek during the spring runoff event, and the location of the flow terminus is directly dependent upon the discharge of the creek. Water levels in the deeper, non-perched aquifers do not appear to fluctuate significantly in response to the arrival of the flow terminus, demonstrating that it is unlikely that the proposed groundwater withdrawals will decrease the surface flow of lower Medano Creek. Received: 27 December 1995 · Accepted: 20 February 1996     

Hydrogeology of the Ordos Basin,China

The Ordos Basin is located in the east of NW China that is composed of different aquifer systems. Karst groundwater is stored in the Carmbrian-Ordovician carbonates along the margins of the basin. Fissured-pore water is present in the Cretaceous strata in the central-western basin and pore water is stored in the overlying Quaternary deposits discontinuously. The main origin of groundwater in the basin is direct or indirect infiltration of precipitation. Groundwater flows from recharge areas to adjacent local discharge areas. Besides evaporation and abstraction, groundwater feeds springs and rivers, such as the Yellow River and its tributaries. According to the karst aquifer lithologic structure, the features of karst development and circulation, the karst aquifer is divided into three structural and circulation patterns. Based on the control of Cretaceous sedimentary environment, lithologic structure, lithofacies, and palaeogeographic characteristics, the Cretaceous system is divided into the northern desert simple plateau aquifer system and the southern loess plateau aquifer system. PACKER was used to obtain temperature, hydrogeochemical and isotope data at specific depths. Groundwater circulation is studied using hydrodynamic fields, temperature fields, isotopes, hydrogeochemical data and numerical simulations. According to the result, it is divided into local, intermediate and regional systems.     

Apportioning deformation among depth intervals in an aquifer system using InSAR and head data

Land surface subsidence due to excessive groundwater pumping is an increasing concern in California, USA. Interferometric Synthetic Aperture Radar (InSAR) is a remote sensing technique for measuring centimeter-to-millimeter surface deformation at 10–100 m spatial resolution. Here, a data-driven approach that attributes deformation to individual depth intervals within an aquifer system by integrating head data acquired from each of three screened intervals in a monitoring well with InSAR surface deformation measurements was developed. The study area was the Colusa Basin in northern Central Valley. To reconstruct the surface deformation history over the study area, 13 ALOS-PALSAR scenes acquired between 2006 and 2010 were processed. Up to ~3-cm year^?1 long-term subsidence and up to ~6 cm seasonal subsidence were observed using the InSAR technique. The technique developed in this paper integrates the InSAR-observed seasonal deformation rate and the co-located head measurements in multiple depth intervals to estimate the elastic skeletal storage coefficient, the time delay between the head change and the observed deformation, and subsequently the deformation of each depth interval. This technique can be implemented when hydraulic head measurements within each depth interval are not correlated with each other. Using this approach, the depth interval that contributed the most to the total subsidence, as well as storage parameters for all intervals, are estimated. The technique can be used for identification of the depth interval within the aquifer system responsible for deformation.

     

Disposition of aquifer system, geo-electric characteristics and gamma-log anomaly in the Kuttanad alluvium of Kerala

Kuttanad near the western coast of Kerala is a low lying area below mean sea level. The hydrogeological scenario in this wet land dominated tract, known as kole lands, was investigated by using the lithological and geophysical data generated under the R&D exploration programme of Central Ground Water Board (CGWB) during the year 2000, as also other related data generated till the year 2009. The hydrogeology of the tract is complex, and the presence of carbonaceous and organic matter in the soil and sub-soil influences the geochemical and geo-electric characteristics of the sub-surface geological formations as reflected by SP and resistivity logs. The lithologs and geophysical logs of bore wells were used to define the aquifer geometry. While the lateral and vertical variations of the aquifer systems identified were high, the bottom aquifer was relatively brackish over a large area. The gamma log behavior in these granular zones is anomalous and misguiding in nature. The radioactivities recorded against clay horizons were in the range of 130 to 200cps and in the granular zones 150 to 350 cps. The presence of monazite bearing sands at certain depths was responsible for anomalously high gamma counts against granular zones. The water quality in the granular zones was inferred from the electrical logs. The depth of occurrence of high radiation sand beds was also demarcated.     

Water availability and land subsidence in the Central Valley,California, USA

The Central Valley in California (USA) covers about 52,000 km² and is one of the most productive agricultural regions in the world. This agriculture relies heavily on surface-water diversions and groundwater pumpage to meet irrigation water demand. Because the valley is semi-arid and surface-water availability varies substantially, agriculture relies heavily on local groundwater. In the southern two thirds of the valley, the San Joaquin Valley, historic and recent groundwater pumpage has caused significant and extensive drawdowns, aquifer-system compaction and subsidence. During recent drought periods (2007–2009 and 2012-present), groundwater pumping has increased owing to a combination of decreased surface-water availability and land-use changes. Declining groundwater levels, approaching or surpassing historical low levels, have caused accelerated and renewed compaction and subsidence that likely is mostly permanent. The subsidence has caused operational, maintenance, and construction-design problems for water-delivery and flood-control canals in the San Joaquin Valley. Planning for the effects of continued subsidence in the area is important for water agencies. As land use, managed aquifer recharge, and surface-water availability continue to vary, long-term groundwater-level and subsidence monitoring and modelling are critical to understanding the dynamics of historical and continued groundwater use resulting in additional water-level and groundwater storage declines, and associated subsidence. Modeling tools such as the Central Valley Hydrologic Model, can be used in the evaluation of management strategies to mitigate adverse impacts due to subsidence while also optimizing water availability. This knowledge will be critical for successful implementation of recent legislation aimed toward sustainable groundwater use.     

Spatio-temporal variation of groundwater recharge in response to variability in precipitation, land use and soil in Yanqing Basin, Beijing, China

It is important to understand how groundwater recharge responds to precipitation variability in space and time, especially in those areas such as Yanqing Basin (China), where groundwater represents the sole water resource. A simple soil-water balance method is applied for spatio-temporal simulation of groundwater recharge in Yanqing Basin from 1981 to 2000. It was implemented on a monthly time step considering the effects of land use and soil texture. The area-average recharge associated with various land uses and soil textures was then compared with zonal analysis using a geographic information system (GIS). The main findings include: (1) the variation of groundwater recharge follows precipitation changes, either at yearly or seasonal intervals, (2) land use plays a more influential role in groundwater recharge than soil texture in this area, and (3) the water table quickly rises in response to recharge in the shallow parts of the aquifer, while there is a delay of 0.5–1.0?years where the groundwater level is at depth 4–10?m. The application demonstrates how spatio-temporal analysis can be utilized for groundwater-recharge estimation through distributed modeling and GIS.     

Identification of cross-valley faults in the Maynardville Limestone,Oak Ridge Reservation,Tennessee, using seismic refraction tomography

First arrival times from P-wave refraction and reflection seismic surveys along Bear Creek Valley on the Oak Ridge Reservation, Tennessee, were inverted to produce refraction tomographic velocity images showing seismic velocity variations within thinly mantled karstic bedrock to a depth of approximately 20 m. Inverted velocities are consistent with two distinct bedrock groups: the Nolichucky Shale (2,730–5,150 m/s) and Maynardville Limestone (3,940–7,575 m/s). Low-velocity zones (2,700–4,000 m/s) in the tomographic images correspond to previously inferred cross-valley strike-slip faults; in places, these faults create permeability barriers that offset or block groundwater flowing along Bear Creek Valley. These faults may also force groundwater contaminants, such as dense non-aqueous phase liquids, to migrate laterally or downward, spreading contamination throughout the groundwater system. Other, previously unmapped cross-valley faults may also be visible in the tomographic images. Borehole logs suggest the low-velocity values are caused by low rigidity fractured and vuggy rock, water zones, cavities and collapse features. Surface streams, including Bear Creek, tend to lie directly above these low-velocity zones, suggesting fault and fracture control of surface drainage, in addition to the subsurface flow system. In some cases, fault zones are also associated with bedrock depressions and thicker accumulations of unconsolidated sediment.     

Hydrogeological implications of paleo-fluvial architecture for the Paskapoo Formation, SW Alberta, Canada: a stochastic analysis

Fluvial systems tend to deposit sediment in well-defined relational geometries and in vertically and laterally repeating patterns. These sedimentary deposits are preserved to varying degrees depending on how much the fluvial system reworks the deposits. The Paskapoo bedrock aquifer system in southern Alberta, Canada, was deposited in a foreland depositional basin during uplift of the Rocky Mountains, and both the geomorphic model and field evidence indicate that the upper 100 m of the local aquifer system contains well-preserved, highly connected paleo-channels and associated overbank deposits. In order to evaluate the value of different types of data, a simplified stochastic-numerical groundwater flow model was developed to examine the sensitivity of results to model parameters. Parameters examined include: fraction of the formation made up of channel sands; meander and sinuosity factors; width-to-depth ratios of preserved channels; and crevasse splay conductivity. In all cases examined, the system exhibited anisotropic behavior with the along-channel flow direction being the most permeable and the vertical direction being least permeable. In general, the strongest control on the resulting effective anisotropic hydraulic conductivities was channel fraction, but geometric factors that control between-channel connectivity (e.g., channel sinuosity) had an appreciable effect on the across-channel flow direction effective permeability.     

A quantitative, three-dimensional depositional model of gravelly braided rivers

A quantitative, three‐dimensional depositional model of gravelly, braided rivers has been developed based largely on the deposits of the Sagavanirktok River in northern Alaska. These deposits were described using cores, wireline logs, trenches and ground‐penetrating radar profiles. The origin of the deposits was inferred from observations of: (1) channel and bar formation and migration and channel filling, interpreted from aerial photographs; (2) water flow during floods; and (3) the topography and texture of the river bed at low‐flow stage. This depositional model quantitatively represents the geometry of the different scales of strataset, the spatial relationships among them and their sediment texture distribution. Porosity and permeability in the model are related to sediment texture. The geometry of a particular type and scale of strataset is related to the geometry and migration of the bedform type (e.g. ripples, dunes, bedload sheets, bars) associated with deposition of the strataset. In particular, the length‐to‐thickness ratio of stratasets is similar to the wavelength‐to‐height ratio of associated bedforms. Furthermore, the wavelength and height of bedforms such as dunes and bars are related to channel depth and width. Therefore, the thickness of a particular scale of strataset (i.e. medium‐scale cross‐sets and large‐scale sets of inclined strata) will vary with river dimensions. These relationships between the dimensions of stratasets, bedforms and channels mean that this depositional model can be applied to other gravelly fluvial deposits. The depositional model can be used to interpret the origin of ancient gravelly fluvial deposits and to aid in the characterization of gravelly fluvial aquifers and hydrocarbon reservoirs.     

Logging in deep water wells in central Jutland, Denmark

During a search for unpolluted resources, new wells have been drilled to a deep confined aquifer in central Jutland. Since little is known about the extent and vulnerability of the deep reservoir, geophysical logging was used to supplement hydraulic well-testing. Gamma-ray, SP and resistivity logs were recorded in the boreholes. Gamma-logs were also recorded in two cased wells. Analysis of the logs, combined with information from samples, shows an unconfined aquifer of quartz and mica sand to a depth of about 55 m at Grindsted. Below this level there is a 40 m thick aquitard with an apparent southerly dip. The aquitard is composed of interbedded sand, silt and clay. The deep confined aquifer extends from 95 to 110 m below the surface at Grindsted. Both the aquitard and the deep aquifer are provisionally dated as Early/Middle Miocene. The electrical logs discriminate well between sand, silt and clay beds. If a suitable mud is used in the well it is possible to calculate the formation water resistivities in good agreement with water samples obtained after completion. From the resistivity logs formation factors were calculated in the range 3–5 using no-invasion departure curves. The high content of heavy mineralsin the Tertiary deposits makes it impossible to obtain precise lithological information from the gamma-ray log. On the other hand this log provides clear markers for identifying the hydrogeological units in the well-bores.     

Case history of one of the few successful Superfund remediation sites: A site at Salinas, California, USA

In 1985, a former tire manufacturing plant surrounded by agricultural fields in the Salinas Valley was designated a Superfund Site by the US Environmental Protection Agency. The plant had been operating for seventeen years, from 1963 until 1980. When dismantling of the plant was started, it was determined that toxic hydrocarbon solvents and oils from the plant had contaminated soil and groundwater in alluvial deposits alongside the plant. It was determined later that the groundwater contamination also lay beneath the agricultural fields in a narrow groundwater plume that extends about 4.3 km downgradient from the plant. Because of the complex architecture of the aquifer system, the gradient, and extensive pumping of agricultural wells, the contaminants migrated northwestward and downward to deeper levels away from the plant.

The agricultural fields are underlain by an unconfined shallow aquifer and by a system of confined aquifers that extend to more than 180 m below surface. Aquifers are discontinuous beds of channel sand and gravel; confining beds are overbank clay and silt, and estuarine clay.

Geophysical data, logs of existing agricultural and other wells, and careful consideration of the stratigraphic architecture of the depositional environment provided the basis for a conceptual hydrogeologic model and for locating characterization wells for detailed visual and geophysical logging and hydrologic testing. Successive refinements of the characterization by sequential installation of wells indicated optimal locations for installation of extraction and monitoring wells. Validity of the concept of the hydrogeologic regime was verified by close match of predictions made by modeling with the later results of pumping from the extraction wells in a pump-and-treat system.

Successful remediation was accomplished by analyzing data from 110 agricultural wells, the few domestic water wells, nearly 200 sequentially installed stratigraphic-characterization and monitoring wells, 25 extraction wells, and by close cooperation among federal, state and local agencies, and the ranchers and growers.

Total contaminants recovered from activated-charcoal strippers of the treatment system totalled < 230 kg. Large quantities were harmlessly volatilized and dispersed into the atmosphere by air strippers and by agricultural sprinkling systems spraying water onto the fields. Crop testing showed no contamination of food crops. The activity has taken seven years and has cost more than US$22 million.     

Mid-Holocene Hydrologic Model of the Shingobee Watershed, Minnesota

A hydrologic model of the Shingobee Watershed in north-central Minnesota was developed to reconstruct mid-Holocene paleo-lake levels for Williams Lake, a surface-water body located in the southern portion of the watershed. Hydrologic parameters for the model were first estimated in a calibration exercise using a 9-yr historical record (1990–1998) of climatic and hydrologic stresses. The model reproduced observed temporal and spatial trends in surface/groundwater levels across the watershed. Mid-Holocene aquifer and lake levels were then reconstructed using two paleoclimatic data sets: CCM1 atmospheric general circulation model output and pollen-transfer functions using sediment core data from Williams Lake.Calculated paleo-lake levels based on pollen-derived paleoclimatic reconstructions indicated a 3.5-m drop in simulated lake levels and were in good agreement with the position of mid-Holocene beach sands observed in a Williams Lake sediment core transect. However, calculated paleolake levels based on CCM1 climate forcing produced only a 0.05-m drop in lake levels. We found that decreases in winter precipitation rather than temperature increases had the largest effect on simulated mid-Holocene lake levels. The study illustrates how watershed models can be used to critically evaluate paleoclimatic reconstructions by integrating geologic, climatic, limnologic, and hydrogeologic data sets.