Groundwater discharge mechanism in semi‐arid regions and the role of evapotranspiration

The present study examined groundwater recharge/discharge mechanisms in the regional Central Sudan Rift Basins (CSRB). Aquifers in CSRB constitute poorly sorted silisiclastics of sand, clay and gravels deposited in closed hydrologic systems of the Cretaceous–Pleistocene fluviolacustrine environments. CSRB are bounded to the north by the highlands of the Central African Shear Zone (CAZS) that represents the surface and groundwater divides. Sporadic recharge in the peripheries of the basins along the CASZ occurs subsequent to decadal and centennial storm events. Inflow from the Nile into the aquifers represents an additional source of recharge. Thus, groundwater resources cannot be labelled fossil nor can they be readily recharged. Closed hydrologic troughs located adjacent to the influent Nile system mark areas of main groundwater discharge characterized by lower hydraulic heads. This study has examined mechanisms that derive the discharge of the groundwater in these closed basins and concluded that only evapotranspirative discharge can provide a plausible explanation. Groundwater abstraction is mainly through deep‐rooted trees and effective evaporation. The increase of TDS along the flow indicates local recharge at the peripheries of basins and shows the influence of evaporation and rock/water interaction. The decline in groundwater level along a flow path was calculated using Darcy's law to estimate average recharge and evapotranspirative discharge, which are equal under natural equilibrium and make the only fluxes in CSRB. Steady‐state 2D flow modelling has demonstrated that an average recharge of 4–8 mm yr^?1 and evapotranspirative discharge of 1–22 mm yr^?1 will maintain natural equilibrium in CSRB. Sporadic storms provide recharge in the highlands to preserve the current hydraulic gradient and maintain aquifer dynamics. Simulated recharge from the Nile totals about 17·5 mm yr^?1 and is therefore a significant contributor to the water balance. Copyright © 2007 John Wiley & Sons, Ltd.     

Groundwater–surface water interaction in Lake Nasser,Southern Egypt

A cross‐sectional model, based on the two dimensional groundwater flow equation of Edelman, was applied at seven transects distributed over four geological cross sections to estimate groundwater heads and recharge from/or groundwater discharge to Lake Nasser. The lake with a length of 500 km and an average width of 12 km was created over the period 1964–1970, the time for constructing the Aswan High Dam (AHD). The model, constrained by regional‐scale groundwater flow and groundwater head data in the vicinity of the lake, was successfully calibrated to timeseries of piezometeric heads collected at the cross sections in the period 1965–2004. Inverse modeling yielded high values for the horizontal hydraulic conductivity in the range of 6.0 to 31.1 m day^?1 and storage coefficient between 0.01 and 0.40. The results showed the existence of a strong vertical anisotropy of the aquifer. The calibrated horizontal permeability is systematically higher than the vertical permeability (≈1000:1). The calibrated model was used to explore the recharge from/or groundwater discharge to Lake Nasser at the seven transects for a 40‐year period, i.e. from 1965 to 2004. The analysis for the last 20‐year period, 1985–2004, revealed that recharge from Lake Nasser reduced by 37% compared to the estimates for the first 20‐year period, 1965–1984. In the period 1965–2004, seepage of Lake Nasser to the surrounding was estimated at 1.15 × 10⁹ m³ year^?1. This led to a significant rise of the groundwater table. Variance‐based sensitivity and uncertainty analysis on the Edelman results were conducted applying quasi‐Monte Carlo sequences (Latin Hypercube sampling). The maximum standard deviation of the total uncertainty on the groundwater table was 0.88 m at Toshka (west of the lake). The distance from the lake, followed by the storage coefficient and hydraulic conductivity, were identified as the most sensitive parameters. Copyright © 2012 John Wiley & Sons, Ltd.     

Ground water flow analysis of a mid-Atlantic outer coastal plain watershed, Virginia, U.S.A

Models for ground water flow (MODFLOW) and particle tracking (MODPATH) were used to determine ground water flow patterns, principal ground water discharge and recharge zones, and estimates of ground water travel times in an unconfined ground water system of an outer coastal plain watershed on the Delmarva Peninsula, Virginia. By coupling recharge and discharge zones within the watershed, flowpath analysis can provide a method to locate and implement specific management strategies within a watershed to reduce ground water nitrogen loading to surface water. A monitoring well network was installed in Eyreville Creek watershed, a first-order creek, to determine hydraulic conductivities and spatial and temporal variations in hydraulic heads for use in model calibration. Ground water flow patterns indicated the convergence of flow along the four surface water features of the watershed; primary discharge areas were in the nontidal portions of the watershed. Ground water recharge zones corresponded to the surface water features with minimal development of a regional ground water system. Predicted ground water velocities varied between < 0.01 to 0.24 m/day, with elevated values associated with discharge areas and areas of convergence along surface water features. Some ground water residence times exceeded 100 years, although average residence times ranged between 16 and 21 years; approximately 95% of the ground water resource would reflect land use activities within the last 50 years.     

Distinguishing Mountain Front and Mountain Block Recharge in an Intermontane Basin of the Himalayan Region

Intermontane basin aquifers worldwide, particularly in the Himalayan region, are recharged largely by the adjoining mountains. Recharge in these basins can occur either by water infiltrating from streams near mountain fronts (MFs) as mountain front recharge (MFR) or by sub-surface mountain block infiltration as mountain block recharge (MBR). MFR and MBR recharge are challenging to distinguish and are least quantified, considering the lack of extensive understanding of the hydrological processes in the mountains. This study used oxygen and hydrogen isotopes (δ¹⁸O and δ²H), electrical conductivity (EC) data, hydraulic head, and groundwater level data to differentiate MFR and MBR. Groundwater level data provide information about the groundwater-surface water interactions and groundwater flow directions, whereas isotopes and EC data are used to distinguish and quantify different recharge sources. The present methodology is tested in an intermontane basin of the Himalayan region. The results suggest that karst springs (KS) and deep groundwater (DGW) recharge are dominated by snowmelt (47% ± 10% and 46% ± 9%) as MBR from adjacent mountains, insignificantly affected by evaporation. The hydraulic head data and isotopes indicate Quaternary shallow groundwater (SGW) aquifer system recharge as MFR of local meteoric water with significant evaporation. The results indicate several flow paths in the aquifer system, a local flow for KS, intermediate flow for SGW, and regional flow for DGW. The findings will significantly impact water resource management in the area and provide vital baseline knowledge for sustainable groundwater management in other Himalayan intermontane basins.     

Groundwater recharge in an arid grassland as indicated by soil chloride profile and multiple tracers

Previous studies have shown that shallow groundwater in arid regions is often not in equilibrium with near‐surface boundary conditions due to human activities and climate change. This is especially the case where the unsaturated zone is thick and recharge rate is limited. Under this nonequilibrium condition, the unsaturated zone solute profile plays an important role in estimating recent diffuse recharge in arid environments. This paper combines evaluation of the thick unsaturated zone with the saturated zone to investigate the groundwater recharge of a grassland in the arid western Ordos Basin, NW China, using the soil chloride profiles and multiple tracers (²H, ¹⁸O, ¹³C, ¹⁴C, and water chemistry) of groundwater. Whereas conventional water balance and Darcy flux measurements usually involve large errors in recharge estimations for arid areas, chloride mass balance has been widely and generally successfully used. The results show that the present diffuse recharge beneath the grassland is 0.11–0.32 mm/year, based on the chloride mass balance of seven soil profiles. The chloride accumulation age is approximately 2,500 years at a depth of 13 m in the unsaturated zone. The average Cl content in soil moisture in the upper 13 m of the unsaturated zone ranges from 2,842 to 7,856 mg/L, whereas the shallow groundwater Cl content ranges from 95 to 351 mg/L. The corrected ¹⁴C age of shallow groundwater ranges from 4,327 to 29,708 years. Stable isotopes show that the shallow groundwater is unrelated to modern precipitation. The shallow groundwater was recharged during the cold and wet phases of the Late Pleistocene and Holocene humid phase based on palaeoclimate, and consequently, the groundwater resources are nonrenewable. Due to the limited recharge rate and thick unsaturated zone, the present shallow groundwater has not been in hydraulic equilibrium with near‐surface boundary conditions in the past 2,500 years.     

Heuristic numerical and analytical models of the hydrologic controls over vertical solute transport in a domed peat bog,Jura Mountains,Switzerland

We report the results of numerical and analytical simulations to test the hypothesis that downward vertical flow of porewater from the crests of domed alpine and kettle bogs controls vertical porewater distributions of major solutes such as Ca and Mg. The domed Etang de la Gruère bog (EGr), Switzerland, characterized by a vertical downward gradient of 0·04 and stratified layers of peat, is chosen as a field site for the model calibration and evaluation. The middle 4‐m section of the 6·5 m thick bog peat is heavily humified and has a hydraulic conductivity of ～10^?5·6 cm s^?1. Above and below, peat is less humified with a hydraulic conductivity of ～10^?3 cm s^?1. Heuristic finite difference simulations, using Visual MODFLOW, of the bog hydraulics show that the higher conductivity peat at the bog base is critical to create the observed deep, local flow cells that substantively recharge porewater. Model results and Peclet number calculations show that before ～7000 ¹⁴C yr BP diffusion of solutes from underlying mineral soils controlled the vertical distribution of porewater chemistry. From 7000 to ～1250 ¹⁴C BP the porewater chemistry was probably controlled by both upward diffusion and downward advection, and after ～1250 ¹⁴C yr BP porewater chemistry was probably controlled by downward advection. Concentrations of conservative major solutes in the porewaters of alpine, ombrotrophic bogs are the net effect of both downward vertical porewater movement and upward vertical diffusion, the magnitudes of which are delicately poised to the configuration of the bog water table over time and subsurface peat stratigraphy. Copyright © 2002 John Wiley & Sons, Ltd.     

Groundwater–surface water interactions in a large semi‐arid floodplain: implications for salinity management

Flow regulation and water diversion for irrigation have considerably impacted the exchange of surface water between the Murray River and its floodplains. However, the way in which river regulation has impacted groundwater–surface water interactions is not completely understood, especially in regards to the salinization and accompanying vegetation dieback currently occurring in many of the floodplains. Groundwater–surface water interactions were studied over a 2 year period in the riparian area of a large floodplain (Hattah–Kulkyne, Victoria) using a combination of piezometric surface monitoring and environmental tracers (Cl⁻, δ²H, and δ¹⁸O). Despite being located in a local and regional groundwater discharge zone, the Murray River is a losing stream under low flow conditions at Hattah–Kulkyne. The discharge zone for local groundwater, regional groundwater and bank recharge is in the floodplain within ∼1 km of the river and is probably driven by high rates of transpiration by the riparian Eucalyptus camaldulensis woodland. Environmental tracers data suggest that the origin of groundwater is principally bank recharge in the riparian zone and a combination of diffuse rainfall recharge and localized floodwater recharge elsewhere in the floodplain. Although the Murray River was losing under low flows, bank discharge occurred during some flood recession periods. The way in which the water table responded to changes in river level was a function of the type of stream bank present, with point bars providing a better connection to the alluvial aquifer than the more common clay‐lined banks. Understanding the spatial variability in the hydraulic connection with the river channel and in vertical recharge following inundations will be critical to design effective salinity remediation strategies for large semi‐arid floodplains. Copyright © 2005 John Wiley & Sons, Ltd.     

Recharge and mineralization of groundwater of the Upper Cretaceous aquifer in Orontes basin,Syria

Abstract

Chemical and isotopic data of groundwater of the Upper Cretaceous aquifer in the Orontes basin, Syria, have been used to assess the groundwater geochemistry, the origin of groundwater recharge and groundwater residence time. The chemical data indicate that dissolution of evaporite minerals is the main process controlling groundwater mineralization. The composition of stable isotopes δ¹⁸O and δ²H, together with ¹⁴C activity, reflect the existence of three different groups: (a) groundwater in the Coastal Mountains with δ¹⁸O of –6.65‰, quite similar to modern-day precipitation, and high ¹⁴C (>50 pmC); (b) groundwater in the unconfined aquifer of the Hama Uplift, which has δ¹⁸O of –5.52‰ and ¹⁴C near 20 pmC, and is recharged locally; and (c) groundwater from the confined aquifer of the Homs Depression, which is characterized by more depleted δ¹⁸O,, –8.01‰, and low ¹⁴C (<7 pmC), and might be recharged in the northern piedmont of the Anti-Lebanon Mountains. The distinctive isotope signatures of the latter two groups indicate different recharge elevations and palaeoclimatic effects. The low recharge rate of the groundwater in the Hama Uplift aquifer, and the even slower recharge rate in the Homs Depression aquifer, are reflected by groundwater ¹⁴C residence times of 5 and over 22 Ka BP, respectively.

Editor D. Koutsoyiannis

Citation Al-Charideh, A., 2013. Recharge and mineralization of groundwater of the Upper Cretaceous aquifer in Orontes basin (Syria). Hydrological Sciences Journal, 58 (2), 452–467.     

Estimating groundwater recharge from water isotope (δ2H, δ18O) depth profiles in the Densu River basin,Ghana

Abstract

Accurate estimation of groundwater recharge is essential for the proper management of aquifers. A study of water isotope (δ²H, δ¹⁸O) depth profiles was carried out to estimate groundwater recharge in the Densu River basin in Ghana, at three chosen observation sites that differ in their altitude, geology, climate and vegetation. Water isotopes and water contents were analysed with depth to determine water flow in the unsaturated zone. The measured data showed isotope enrichment in the pore water near the soil surface due to evaporation. Seasonal variations in the isotope signal of the pore water were also observed to a depth of 2.75 m. Below that depth, the seasonal variation of the isotope signal was attenuated due to diffusion/dispersion and low water flow velocities. Groundwater recharge rates were determined by numerical modelling of the unsaturated water flow and water isotope transport. Different groundwater recharge rates were computed at the three observation sites and were found to vary between 94 and 182 mm/year (± max. 7%). Further, the approximate peak-shift method was applied to give information about groundwater recharge rates. Although this simple method neglects variations in flow conditions and only considers advective transport, it yielded mean groundwater recharge rates of 110–250 mm/year (± max. 30%), which were in the same order of magnitude as computed numerical modelling values. Integrating these site-specific groundwater recharge rates to the whole catchment indicates that more water is potentially renewed than consumed nowadays. With increases in population and irrigation, more clean water is required, and knowledge about groundwater recharge rates – essential for improving the groundwater management in the Densu River basin – can be easily obtained by measuring water isotope depth profiles and applying a simple peak-shift approach.

Citation Adomako, D., Maloszewski, P., Stumpp, C., Osae, S. & Akiti, T. T. (2010) Estimating groundwater recharge from water isotope (δ²H, δ¹⁸O) depth profiles in the Densu River basin, Ghana. Hydrol. Sci. J. 55(8), 1405–1416.     

Groundwater recharge in arid regions: Review and critique of estimation methods

Arid-site recharge, while generally low, can be highly variable. Recharge under similar climate and soil conditions but with different plant cover and topography can vary from zero to more than the annual precipitation. Simple estimates of recharge based on fixed fractions of annual precipitation are misleading because they do not reflect the plant and soil factors controlling recharge. Detailed water balance models, successful for irrigated agriculture, fail to predict evapotranspiration accurately under conditions where plants suffer seasonal water stress and cover is sparse. Recharge, when estimated as a residual in water balance models, may be in error by as much as an order of magnitude. Similar errors can occur when soil water flow models are used with measured or estimated soil hydraulic conductivities and tension gradients. Lysimetry and tracer tests offer the best hope for evaluating recharge at arid sites, particularly in siting waste disposal facilities, where reliable recharge estimates are needed. Quantification of drainage using lysimetry over several years under a given set of soil, plant, and climate conditions for a specific site can provide a basis for calibrating models for recharge prediction. Tracer tests using such long-lived tracers as ³⁶Cl or perhaps stable isotopes (¹⁸0, deuterium) can provide qualitative estimates of recent recharge at a given site.     

Influence of aquifer properties and the spatial and temporal distribution of recharge and abstraction on sustainable yields in semi‐arid regions

This paper aims to contribute to understanding the importance of four factors on the determination of sustainable yields: (i) aquifer properties; (ii) temporal distribution of recharge; (iii) temporal distribution of groundwater pumping; and (iv) spatial distribution of pumping wells. It is important to comprehend how the present‐day and future vulnerability of groundwater systems to pumping activities depend on these critical factors and what the risks are of considering sustainable yield as a fixed percentage of mean annual recharge (MAR). A numerical model of the Querença–Silves aquifer in Portugal is used to develop hypothetical scenarios with which these factors are studied. Results demonstrate the aquifer properties, particularly the storage coefficient, have an important role in determining the resilience of an aquifer and therefore to which degree it is dependent on the spatial and temporal distribution of abstraction and recharge, as well as the occurrence of extreme events. Sustainable yields are determined for the developed scenarios based on specific criteria rather than a fraction of MAR. Under simplified current recharge and abstraction conditions, the sustainable yield was determined at approximately 73% of MAR or 76 million m³. When considering a concentration of rainfall in time, as predicted by climate scenarios for the region, sustainable yield could drop to ca 70% of MAR. However, a more even distribution of pumping volumes throughout the year could increase this value. The location of the pumping wells is seen to affect the distribution of hydraulic heads in the aquifer, albeit without significant changes in sustainable yield. Copyright © 2011 John Wiley & Sons, Ltd.     

Impacts of human activities on recharge in a multilayered semiarid aquifer (Campo de Cartagena,SE Spain)

The development of intense agriculture in semiarid areas modifies intensity and spatial distribution of groundwater recharge by summing irrigation return flow to limited rainfall infiltration. Environmental tracers provide key information, but their interpretation is complicated by more complex groundwater flow patterns. In multilayered aquifers, the real origin of the groundwater samples is hard to assess because of local mixing processes occurring inside long‐screened boreholes. We use environmental tracers (¹⁴C, ¹³C, ²H, ¹⁸O, ³H) to investigate the long‐term evolution of recharge in the five‐layer Campo de Cartagena aquifer in South‐Eastern Spain, in addition to high‐resolution temperature loggings to identify the depth of origin of groundwater. Despite the complex background, this methodology allowed a reliable interpretation of the geochemistry and provided a better understanding of the groundwater flow patterns. The tritium method did not give good quantitative results because of the high variability of the recharge signal but remained an excellent indicator of recent recharge. Nonetheless, both pre‐anthropization and post‐anthropization recharge regime could be identified and quantified by radiocarbon. Before the development of agriculture, recharge varied from 17 mm.year^‐1 at the mountain ranges to 6 mm.year^‐1 in the plain, whereas the mean annual rainfall is about 300 mm. In response to the increase of agricultural activity, recharge fluxes to the plain were amplified and nowadays reach up to 210 mm.year^‐1 in irrigated areas. These values are strengthened by global water budget and local unsaturated zone studies. Copyright © 2013 John Wiley & Sons, Ltd.     

Channel network extent in the context of historical land use,flow generation processes,and landscape evolution in the North Carolina Piedmont

Intensive agricultural land use in the 18th to early 20th centuries on the southeastern Piedmont resulted in substantial soil erosion and gully development. Today, many historically farmed areas have been abandoned and afforested, and such landscapes are an opportunity to study channel network recovery from disturbance by gullying. Channel initiation mapping, watershed area–slope relationships, and field monitoring of flow generation processes are used to identify channel network extent and place it in hydrologic, historical and landscape evolution context. In six study areas in the North Carolina Piedmont, 100 channel heads were mapped in fully‐forested watersheds, revealing a channel initiation relationship of 380 = AS^1.27, where A is contributing area (m²) and S is local slope (m/m). Flow in these channels is generated by subsurface and overland flow. The measured relative slope exponent is lower than expected based on literature values of ~2 for forested watersheds with subsurface and overland flow, suggesting that the channel network extent may reflect a former hydrological regime. However, geomorphic evidence of recovery in channel heads within fully forested watersheds is greater than those with present day pasture. Present day channel heads lie within hollows or downslope of unchanneled valleys, which may be remnants of historical gullies, and area–slope relationships provide evidence of colluvial aggradation within the valleys. Channel network extent appears to be sensitive to land use change, with recovery beginning within decades of afforestation. Channel initiation mapping and area–slope relationships are shown to be useful tools for interpreting geomorphic effects of land use change. Copyright © 2012 John Wiley & Sons, Ltd.     

Groundwater temperature evolution in the subsurface urban heat island of Cologne,Germany

Long‐term heating of shallow urban aquifers is observed worldwide. Our measurements in the city of Cologne, Germany revealed that the groundwater temperatures found in the city centre are more than 5 K higher than the undisturbed background. To explore the role of groundwater flow for the development of subsurface urban heat islands, a numerical flow and heat transport model is set up, which describes the hydraulic conditions of Cologne and simulates the transient evolution of thermal anomalies in the urban ground. A main focus is on the influence of horizontal groundwater flow, groundwater recharge and trends in local ground warming. To examine heat transport in groundwater, a scenario consisting of a local hot spot with a length of 1 km of long‐term ground heating was set up in the centre of the city. Groundwater temperature‐depth profiles at upstream, central and downstream locations of this hot spot are inspected. The simulation results indicate that the main thermal transport mechanisms are long‐term vertical conductive heat input, horizontal advection and transverse dispersion. Groundwater recharge rates in the city are low (<100 mm a^?1) and thus do not significantly contribute to heat transport into the urban aquifer. With groundwater flow, local vertical temperature profiles become very complex and are hard to interpret, if local flow conditions and heat sources are not thoroughly known. Copyright © 2014 John Wiley & Sons, Ltd.     

A Black Hills-Madison Aquifer Origin for Dakota Aquifer Groundwater in Northeastern Nebraska

Previous studies of the Dakota Aquifer in South Dakota attributed elevated groundwater sulfate concentrations to Madison Aquifer recharge in the Black Hills with subsequent chemical evolution prior to upward migration into the Dakota Aquifer. This study examines the plausibility of a Madison Aquifer origin for groundwater in northeastern Nebraska. Dakota Aquifer water samples were collected for major ion chemistry and isotopic analysis (¹⁸O, ²H, ³H, ¹⁴C, ¹³C, ³⁴S, ¹⁸O-SO₄, ⁸⁷Sr, ³⁷Cl). Results show that groundwater beneath the eastern, unconfined portion of the study area is distinctly different from groundwater sampled beneath the western, confined portion. In the east, groundwater is calcium-bicarbonate type, with δ¹⁸O values (−9.6 to −12.4) similar to local, modern precipitation (−7.4 to −10), and tritium values reflecting modern recharge. In the west, groundwater is calcium-sulfate type, having depleted δ¹⁸O values (−16 to −18) relative to local, modern precipitation, and ¹⁴C ages 32,000 to more than 47,000 years before present. Sulfate, δ¹⁸O, δ²H, δ³⁴S, and δ¹⁸O-SO₄ concentrations are similar to those found in Madison Aquifer groundwater in South Dakota. Thus, it is proposed that Madison Aquifer source water is also present within the Dakota Aquifer beneath northeastern Nebraska. A simple Darcy equation estimate of groundwater velocities and travel times using reported physical parameters from the Madison and Dakota Aquifers suggests such a migration is plausible. However, discrepancies between ¹⁴C and Darcy age estimates indicate that ¹⁴C ages may not accurately reflect aquifer residence time, due to mixtures of varying aged water.     

The Effect of Modeled Recharge Distribution on Simulated Groundwater Availability and Capture

Simulating groundwater flow in basin‐fill aquifers of the semiarid southwestern United States commonly requires decisions about how to distribute aquifer recharge. Precipitation can recharge basin‐fill aquifers by direct infiltration and transport through faults and fractures in the high‐elevation areas, by flowing overland through high‐elevation areas to infiltrate at basin‐fill margins along mountain fronts, by flowing overland to infiltrate along ephemeral channels that often traverse basins in the area, or by some combination of these processes. The importance of accurately simulating recharge distributions is a current topic of discussion among hydrologists and water managers in the region, but no comparative study has been performed to analyze the effects of different recharge distributions on groundwater simulations. This study investigates the importance of the distribution of aquifer recharge in simulating regional groundwater flow in basin‐fill aquifers by calibrating a groundwater‐flow model to four different recharge distributions, all with the same total amount of recharge. Similarities are seen in results from steady‐state models for optimized hydraulic conductivity values, fit of simulated to observed hydraulic heads, and composite scaled sensitivities of conductivity parameter zones. Transient simulations with hypothetical storage properties and pumping rates produce similar capture rates and storage change results, but differences are noted in the rate of drawdown at some well locations owing to the differences in optimized hydraulic conductivity. Depending on whether the purpose of the groundwater model is to simulate changes in groundwater levels or changes in storage and capture, the distribution of aquifer recharge may or may not be of primary importance.     

Radially convergent groundwater flow in sloping terrain

Abstract

The groundwater flow equation governing the elevation (h) of the steady-state phreatic surface in a sloping aquifer fed by constant recharge over a bi-circular sector is rhh′ ? r ² Bh′ + Pr ² ? PR ² = 0, where r is the radial coordinate, P is a constant involving recharge and aquifer properties, and B is the slope of the aquifer—bedrock boundary. The derived flow equation describes radially convergent flow through a sloping aquifer that discharges to a water body of fixed head. One important simplification is that in which the width of the bi-circular sector is constant, and the draining land becomes a rectangular aquifer. The bi-circular sector and rectangular-strip groundwater flow problems are solved in terms of implicit equations. The solutions for the steady-state phreatic surfaces depend on the ratio of recharge to hydraulic conductivity, the slope of the aquifer-bedrock, and the downstream constant-head boundary. Computational examples illustrate the application of the solutions.     

Groundwater evolution and mean water age inferred from hydrochemical and isotopic tracers in a tropical confined aquifer

The coastal confined aquifer in the Gulf of Urabá (Colombia) is an important water source for the banana agro‐industry as well as for urban and rural communities. However, the main processes controlling recharge and mixing in the aquifer are still poorly understood. Hydrochemical analyses and stable isotope monitoring were conducted to (a) determine groundwater recharge origin, mean groundwater age, and the main processes governing groundwater chemistry and the potential mixing of marine water and the influence of diffusive processes from the two surrounding aquitard layers. Hydrochemical data indicate that the main processes affecting the dissolved chemical composition include cation exchange, dissolution of carbonated and CO₂, and silicate weathering. δ¹⁸O and δ²H compositions combined with ¹⁴C data highlight the differences in climatic conditions between the recharge zone and the confined section of the aquifer, which is close to the Atlantic Ocean. Groundwater samples with ¹⁴C ages from recent to 28,300 years BP show a depleted isotopic trend ranging from ?6.43‰ to ?9.14‰ in δ¹⁸O and from ?43.2‰ to ?65.7‰ in δ²H. The most depleted δ¹⁸O and δ²H compositions suggest a cooler recharge climate than the current conditions (corresponding to the last glacial period of the late Pleistocene). Depleted δ¹³C values in the total dissolved inorganic carbon indicate the existence of organic material oxidation processes within the geologic formation. These results can be used or transferred to enhance groundwater modelling efforts in other confined coastal aquifers of South America where scarcity of long‐term monitoring data limits water resources planification under a changing climate.     

Innovative Environmental Tracer Techniques for Evaluating Sources of Spring Discharge from a Carbonate Aquifer Bisected by a River

Littlefield Springs discharge about 1.6 m³/s along a 10‐km reach of the Virgin River in northwestern Arizona. Understanding their source is important for salinity control in the Colorado River Basin. Environmental tracers suggest that Littlefield Springs are a mixture of older groundwater from the regional Great Basin carbonate aquifer and modern (post‐1950s) seepage from the Virgin River. While corrected ¹⁴C apparent ages range from 1 to 9 ka, large amounts of nucleogenic ⁴He and low ³He/⁴He ratios suggest that the carbonate aquifer component is likely even older Pleistocene recharge. Modeled infiltration of precipitation, hydrogeologic cross sections, and hydraulic gradients all indicate recharge to the carbonate aquifer likely occurs in the Clover and Bull Valley Mountains along the northern part of the watershed, rather than in the nearby Virgin Mountains. This high‐altitude recharge is supported by relatively cool noble‐gas recharge temperatures and isotopically depleted δ²H and δ¹⁸O. Excess (crustal) SF₆ and ⁴He precluded dating of the modern component of water from Littlefield Springs using SF₆ and ³H/³He methods. Assuming a lumped‐parameter model with a binary mixture of two piston‐flow components, Cl^?/Br^?, Cl^?/F^?, δ²H, and CFCs indicate the mixture is about 60% Virgin River water and 40% groundwater from the carbonate aquifer, with an approximately 30‐year groundwater travel time for Virgin River seepage to re‐emerge at Littlefield Springs. This suggests that removal of high‐salinity sources upstream of the Virgin River Gorge would reduce the salinity of water discharging from Littlefield Springs into the Virgin River within a few decades.     

Hydraulic conductivity in frozen soils

An apparatus has been developed for investigation of hydraulic conductivity of frozen soils. The test procedure is isothermal and involves the passage of water from one reservoir into the frozen sample and out of the frozen sample into a second reservoir. The water in the reservoirs remains unfrozen because it contains dissolved lactose. The concentration of lactose is such that, initially, the water in the reservoirs is in thermodynamic equilibrium with the water in the soil. On application of pressure to one reservoir a known hydraulic gradient is established and flow takes place. Flow is shown to vary linearly with hydraulic gradient. The hydraulic conductivity coefficient depends on soil type and temperature and is related to the unfrozen water content. At temperatures within a few tenths of 0°C the coefficient apparently ranges from 10^?5 to 10^?9 cm sec^?1, and decreases only slowly below about ?0·5°C. Soils known to be susceptible to frost heave are shown to have significant hydraulic conductivities well below 0°C.