Prediction of groundwater‐induced flooding in a chalk aquifer for future climate change scenarios

Current climate change models for the southeast UK predict changing rainfall patterns, with increased incidence of extreme events. The chalk aquifer in the UK and northern France is susceptible to groundwater‐induced flooding under such conditions. In this methodological study we apply a frequency domain analysis approach to the chalk aquifer to derive a transfer function between effective rainfall and groundwater level from 7 years of monitoring data from the North Heath Barn site, near Brighton. The derived transfer function was calibrated and validated against monitoring data and then used to predict groundwater level for rainfall models for high, medium and low emission scenarios from the UKCP09 database. The derived transfer function is most closely comparable to the linear aquifer model, despite evidence for both matrix and fracture or karst water flow in the chalk, with transmissivity and unconfined storativity at the catchment scale of 1548 m² day^?1 and 1.6 × 10^?2. The application of the transfer function to UKCP09 rainfall data suggests that groundwater‐induced flooding may be about four times more frequent by 2040–2069 compared with 1961–1990 and seven times more frequent by 2070–2099. The model data also suggest an increase in the duration of groundwater minima relative to the reference period. Compared to deterministic modelling which requires detailed knowledge of aquifer heterogeneity and processes, the transfer function approach, although with limitations, is simpler, incorporating these factors into the analysis through frequency and phase coefficients, and thus may have the potential for groundwater risk assessment in other areas. Copyright © 2015 John Wiley & Sons, Ltd.     

Simulation of the spatio‐temporal extent of groundwater flooding using statistical methods of hydrograph classification and lumped parameter models

This article presents the development of a relatively low cost and rapidly applicable methodology to simulate the spatio‐temporal occurrence of groundwater flooding in chalk catchments. In winter 2000/2001 extreme rainfall resulted in anomalously high groundwater levels and groundwater flooding in many chalk catchments of northern Europe and the southern United Kingdom. Groundwater flooding was extensive and prolonged, occurring in areas where it had not been recently observed and, in places, lasting for 6 months. In many of these catchments, the prediction of groundwater flooding is hindered by the lack of an appropriate tool, such as a distributed groundwater model, or the inability of models to simulate extremes adequately. A set of groundwater hydrographs is simulated using a simple lumped parameter groundwater model. The number of models required is minimized through the classification and grouping of groundwater level time‐series using principal component analysis and cluster analysis. One representative hydrograph is modelled then transposed to other observed hydrographs in the same group by the process of quantile mapping. Time‐variant groundwater level surfaces, generated using the discrete set of modelled hydrographs and river elevation data, are overlain on a digital terrain model to predict the spatial extent of groundwater flooding. The methodology is applied to the Pang and Lambourn catchments in southern England for which monthly groundwater level time‐series exist for 52 observation boreholes covering the period 1975–2004. The results are validated against observed groundwater flood extent data obtained from aerial surveys and field mapping. The method is shown to simulate the spatial and temporal occurrence of flooding during the 2000/2001 flood event accurately. British Geological Survey © NERC 2011. Hydrological Processes © 2011 John Wiley & Sons, Ltd.     

Modelling of the groundwater flow and of tracer movement in the porous and fissured media: Chalk Aquifer (Northern part of Paris Basin,France)

A groundwater flow model has been developed in order to study the chalk aquifer of Paris Basin, based on most of the geological and hydrological available data. The numerical processes are intended to modelling the groundwater flow in the Senonian (Late Cretaceous) formations and to visualize the tracer movement in groundwater resources in the experimental site of LaSalle Beauvais (northern part Paris Basin). Both objectives were achieved as follows: (i) the comprehension of the spatial distribution of the hydraulic conductivity in the chalk aquifer taking into account the characteristics of the hydrogeological system and (ii) the use of the analytical solution for describing one‐dimensional to two‐dimensional solute transport in a unidirectional steady‐state flow tracer with scale‐dependent dispersion. Advection and diffusion mechanisms are taken into account. Comparison between the breakthrough curves of the analytical and the numerical solutions provided an excellent agreement for various ranges of scale‐related transport parameters of interest. The developed power series solution facilitates fast prediction of the breakthrough curves at each observation point. Thus, the derived new solutions are widely applicable and are very useful for the validation of numerical transport. The numerical approach is carried out by MT3DMS, a Modular 3‐D Multi‐Species Transport Model for Simulation of Advection, Dispersion, and Chemical Reactions of Contaminants in Groundwater Systems, and based on total variation‐diminishing method using the ULTIMATE algorithm. The estimation of the infected surface could constitute an approach in water management and allows to prevent the risks of pollution and to manage the groundwater resource from a durable development perspective. Copyright © 2016 John Wiley & Sons, Ltd.     

A proposed structure–groundwater model: application to chalky media of the Paris Basin

The groundwater flow in a fissured chalky environment at the northern border of the Paris Basin depends on several geological and hydrogeological parameters. Although the studied sector of the basin presents a homogeneous rock type, it is affected by a fracture network. In this type of environment, in which the permeability is low, the groundwater flow displays significant disruption, which is localized in the Fruges region (northern France). The interconnection of the discontinuities (network of fault and/or joints) is reliant on the structural control of groundwater flow through increases in the hydraulic connection between the unsaturated and the saturated zone. The methodology developed herein makes use of microstructural and regional analysis of the fracture patterns, and allowed consideration of the piezometric variations of the chalk aquifer during periods of low and high groundwater levels (April and October 2001) and a diagraphic representation of the estimated physical parameters (electrical resistivity). This enabled us to construct a ‘flow structure’ conceptual model in which we identify two types of faults: tight walls and flow paths that control the piezometric heads and the flow rate. Model validation was carried out on a similar sector. Copyright © 2007 John Wiley & Sons, Ltd.     

Modelling groundwater/surface water interaction in a managed riparian chalk valley wetland

Understanding hydrological processes in wetlands may be complicated by management practices and complex groundwater/surface water interactions. This is especially true for wetlands underlain by permeable geology, such as chalk. In this study, the physically based, distributed model MIKE SHE is used to simulate hydrological processes at the Centre for Ecology and Hydrology River Lambourn Observatory, Boxford, Berkshire, UK. This comprises a 10‐ha lowland, chalk valley bottom, riparian wetland designated for its conservation value and scientific interest. Channel management and a compound geology exert important, but to date not completely understood, influences upon hydrological conditions. Model calibration and validation were based upon comparisons of observed and simulated groundwater heads and channel stages over an equally split 20‐month period. Model results are generally consistent with field observations and include short‐term responses to events as well as longer‐term seasonal trends. An intrinsic difficulty in representing compressible, anisotropic soils limited otherwise excellent performance in some areas. Hydrological processes in the wetland are dominated by the interaction between groundwater and surface water. Channel stage provides head boundaries for broad water levels across the wetland, whilst areas of groundwater upwelling control discrete head elevations. A relic surface drainage network confines flooding extents and routes seepage to the main channels. In‐channel macrophyte growth and its management have an acute effect on water levels and the proportional contribution of groundwater and surface water. The implications of model results for management of conservation species and their associated habitats are discussed. Copyright © 2015 John Wiley & Sons, Ltd.     

Dynamics of river–aquifer interactions along a chalk stream: the River Lambourn,UK

This paper describes the spatial and temporal pattern of groundwater flow accretion to the River Lambourn, a 234 km² chalk catchment of the West Berkshire Downs, UK, which has been largely unaffected by groundwater abstraction. Variations in the discharge measured at four fixed gauges in the catchment, coupled with information on the length of flowing channel over the period 1983–2001, are used to describe regional patterns in flow accretion. Mean catchment accretion generally exceeds 0·15 m³ s⁻¹ km⁻¹, but there are significant differences between perennial reaches indicating how the combination of local structural controls and seasonal changes in the drainage net affect flow accretion. Data from current meter surveys were used to determine the spatial variability in flow accretion: 505 paired observations along 12 reaches between 1 and 2·95 km in length indicated a consistent spatial trend in accretion. Accretion was high in upstream and downstream channel reaches, and in middle reaches where dry valleys intersected the main valley. A flow accretion index was developed to describe the relationship of flow accretion in each of the 12 study reaches to catchment discharge. The relationship varied from a strong positive correlation with catchment discharge (two reaches), a weak positive correlation (three reaches), a strong negative correlation (two reaches), to no relationship to catchment discharge (four reaches). The results highlight the need to reconsider the usual assumption of uniform, or uniformly increasing, flow accretion in chalk catchments. Moreover, they emphasize the importance of catchment topography, and illustrate how flow accretion in individual reaches may vary between high and low groundwater levels. Copyright © 2005 John Wiley & Sons, Ltd.     

Numerical Modeling of Nitrate Removal in Anoxic Groundwater during River Flooding of Riparian Zones

A numerical study demonstrates the effects of flooding on subsurface hydrological flowpaths and nitrate removal in anoxic groundwater in riparian zones with a top peat layer. A series of two-dimensional numerical simulations with changing conditions for flow (steady state or transient with flooding), hydrogeology, denitrification, and duration of flooding demonstrate how flowpaths, residence times, and nitrate removal are affected. In periods with no flooding groundwater flows horizontally and discharges to the river through the riverbed. During periods with flooding, shallow groundwater is forced upwards as discharge through peat layers that often have more optimal conditions for denitrification caused by the presence of highly reactive organic matter. The contrast in hydraulic conductivity between the sand aquifer and the overlying peat layer, as well as the flooding duration, have a significant role in determining the degree of nitrate removal.     

Fault zone hydrogeology in arid environments: The origin of cold springs in the Wadi Araba Basin,Egypt

The role of faults in controlling groundwater flow in the Sahara and most of the hyper-arid deserts is poorly understood due to scarcity of hydrological data. The Wadi Araba Basin (WAB), in the Eastern Sahara, is highly affected by folds and faults associated with Senonian tectonics and Paleogene rifting. Using the WAB as a test site, satellite imagery, aeromagnetic maps, field observations, isotopic and geochemical data were examined to unravel the structural control on groundwater flow dynamics in the Sahara. Analysis of satellite imagery indicated that springs occur along structurally controlled scarps. Isotopic data suggested that cold springs in the WAB showed a striking similarity with the Sinai Nubian aquifer system (NAS) water and the thermal springs along the Gulf of Suez (e.g., δ¹⁸O = −8.01‰ to −5.24‰ and δD = −53.09‰ to −31.12‰) demonstrating similar recharge sources. The findings advocated that cold springs in the WAB represent a natural discharge from a previously undefined aquifer in the Eastern Desert of Egypt rather than infiltrated precipitation over the plateaus surrounding the WAB or through hydrologic windows from deep crystalline basement flow. A complex role of the geological structures was inferred including: (1) channelling of the groundwater flow along low-angle faults, (2) compartmentalization of the groundwater flow upslope from high-angle faults, and (3) reduction of the depth to the main aquifer in a breached anticline setting, which resulted in cold spring discharge temperatures (13–22°C). Our findings emphasize on the complex role of faults and folds in controlling groundwater flow, which should be taken into consideration in future examination of aquifer response to climate variability in the Sahara and similar deserts worldwide.     

Modelling of nitrate leaching from arable land into unsaturated soil and chalk 1. Development of a management model for applications of sewage sludge and fertilizer

A mathematical model has been produced to examine the impact of sewage sludge and fertilizer application to arable land and the effect of different crop regimes on the amount of nitrate leached to chalk groundwater. Previous work on nitrate leaching has concentrated on either a soil science or a hydrogeological approach with little overlap between the two. This study considered both fields to obtain an overall picture of the nitrate leaching process. IMPACT is a layered deterministic N-leaching model which predicts the nitrogen loads entering groundwater daily from arable land, and can be used as a management tool in development of sludge application and agricultural policy. The model relates nitrogen species movement resulting from the application of sewage sludge and fertilizer to differing vegetation-soil-hydrogeological conditions. Field data collected at three sites on the unconfined chalk aquifer of East Anglia, England over a two and a half year period was used to produce an initial conceptual model and to constrain the mathematical model during development. IMPACT simulates nitrogen and transport processes in the soil and unsaturated zone of the chalk. The nitrogen processes include: mineralisation of soil organic-N and sewage sludge organic-N, nitrification; crop uptake; volatilization; denitrification; and N inputs from fertilizers and precipitation. A mixing cell method is used to model solute transport in both the soil and chalk. Matrix flow and combined fissure-matrix flow are considered for the chalk. The model enables examination of the relationship between the arable/hydrogeological systems and the environmental implications of sludge application and of different arable regimes. Results are of use in developing strategies for arable farming and sludge application in areas sensitive to nitrate leaching. This Part 1 paper describes the model development approach. Results of associated modelling scenarios are presented separately in the associated Part 2 paper.     

Hydrodynamics of floodplain wetlands in a chalk catchment: The River Lambourn, UK

Variations in floodplain channel water levels and valley floor groundwater levels (measured in piezometers and boreholes) are examined at selected points along the course of the River Lambourn, a chalk river in southern England. A local alluvial gravel aquifer in the valley bottom is associated with numerous small wetlands that extend over much of the river's perennial profile. Variations in hydraulic gradient between local borehole levels and/or floodplain channel water levels are described for three sites in the seasonal section of the channel at Bockhampton, East Garston and West Shefford. The results indicate that observed groundwater levels are closely associated with flows from discrete springs at the margins of the channel and floodplain. However, as the floodplain widens and the alluvial gravel aquifer increases in size, the gravel aquifer accounts for a substantial down-valley component of groundwater flow with a diffuse vertical water flux. In the lower catchment, the exchange of flows between the gravel aquifer and the river enables some attenuation of floodplain water-table variability, providing a stable hydrological regime for valley-bottom wetlands. Catchment controls upon the local, valley-bottom, wetland regime are demonstrated with the application of a simple groundwater model developed using MODFLOW. The model is used to simulate groundwater discharge to the river in the upper and lower catchment, in addition to the water level regime at selected points in the valley bottom in the lower catchment. The results demonstrate the importance of taking catchment-scale water flow into account when managing isolated wetlands in a permeable catchment.     

Tracing groundwater flow using chemical data

Abstract

In the southern part of the Great Hungarian Plain there are two different groundwater flow systems located in the study area. Between these two systems there is a narrow groundwater divide where groundwater flow cannot be detected. In the western part, groundwater of greater hardness moves from the northwest, west and southwest. In the eastern part, groundwater flow is from southeast to northwest. The directions of groundwater flow have been established on the basis of dissolved mono- and di-valent cation concentrations. The major direction of groundwater flow detected by a statistical evaluation of water chemical data agrees with previous geological investigations.     

Methods of the Theory of Water Flow in Porous Media in Problems Arising in Bog Watering

The analysis of hydrogeological problems arising in bog watering by methods enabling the depth to groundwater table to be controlled was used to develop a method for estimating the efficiency of such control when data for directly solving groundwater flow problems are lacking. The effect of nonmonotone dependence of the rate of groundwater level rise on the hydraulic conductivity of peat underlain by high-permeability sand rocks is described. This effect is shown to be significant in assessing the potentialities of the operational control of groundwater level in bog massifs with this type of geological structure.     

Groundwater dynamics in a till hillslope: flow directions,gradients and delay

Knowledge of groundwater dynamics is important for the understanding of hydrological controls on chemical processes along the water flow pathways. To increase our knowledge of groundwater dynamics in areas with shallow groundwater, the groundwater dynamics along a hillslope were studied in a boreal catchment in Southern Sweden. The forested hillslope had a 1‐ to 2‐m deep layer of sandy till above bedrock. The groundwater flow direction and slope were calculated under the assumption that the flow followed the slope of the groundwater table, which was computed for different triangles, each defined by three groundwater wells. The flow direction showed considerable variations over time, with a maximum variation of 75°. During periods of high groundwater levels the flow was almost perpendicular to the stream, but as the groundwater level fell, the flow direction became gradually more parallel to the stream, directed in the downstream direction. These findings are of importance for the interpretation of results from hillslope transects, where the flow direction usually is assumed to be invariable and always in the direction of the hillslope. The variations in the groundwater flow direction may also cause an apparent dispersion for groundwater‐based transport. In contrast to findings in several other studies, the groundwater level was most responsive to rainfall and snowmelt in the upper part of the hillslope, while the lower parts of the slope reached their highest groundwater level up to 40 h after the upper parts. This can be explained by the topography with a wetter hollow area in the upper part. Copyright © 2011 John Wiley & Sons, Ltd.     

Spatial delineation of riparian groundwater within alluvium deposit of mountainous region using Laplace equation

The interaction between surface and groundwater plays a key role in a riparian ecosystem while the size of riparian groundwater has not been typically incorporated into hydrological modelling systems. An extensive geophysical survey composed of 25 individual DC electrical resistivity profiles was conducted at the Blair–Wallis site in Wyoming. The observed resistivity images show a near‐surface aquifer interpreted as the saturated alluvium deposit along the channel, rather than the geological bedrock. Based on the electrical resistivity images, it can be inferred that only the near‐surface portion of the groundwater actively interacts with the stream flow in the mountainous and hilly watershed. This study attempted the spatial extrapolation of the measured riparian aquifer depths by means of fitting functions based on the surface topography. The analysis indicated that the boundary of the riparian aquifer well corresponds to the topographical inflexion point of the hill slope profile. It was also demonstrated that the extent of alluvium deposit, where the area of riparian aquifer is indicated, can be delineated using the slope and curvature maps in the geographic information system. Then, the parabolic and biharmonic functions were tested for the groundwater depth estimation using the developed alluvium deposit map. The proposed methodology was effective if geological diffusion processes by wind and water dominated the topography. The spatial map of the active aquifer will be useful in hydrological drought analysis because it is considered to be a main source of baseflow during dry seasons.     

Groundwater occurrence and flow pattern in the Enugu coal-mine area,Anambra State,Nigeria

Abstract

The distribution of major geological units, static water level data, water chemistry data, and observations of surface features influenced by groundwater seepage were used to ascertain the nature of groundwater occurrence and flow pattern in the Enugu coal field, Nigeria. Considerations of the geological units, the static groundwater levels and groundwater seepages in the mines indicate that the coal sequence is a multiaquifer system in which sandstone and coal aquifers alternate with shale aquitards. Based on the hydraulic head data, the groundwater flow is predominantly downwards. Groundwater velocity calculation across the multiaquifer system using the Darcy equation gave a flow velocity of about 1 m day^?1. For groundwater systems, such a calculated velocity is considered high. The high velocity is most probably due to the high fracture porosity as well as the presence of other stratigraphic and structural features such as alluvial fills that provide high hydraulic conductivity pathways across the aquifer system. The pattern of groundwater inflow into the mines is also influenced by these stratigraphie and structural features.     

Assessing the Transport of Pharmaceutical Compounds in a Layered Aquifer Discharging to a Stream

A groundwater plume containing high concentrations of pharmaceutical compounds, mainly sulfonamides, barbiturates, and ethyl urethane, in addition to chlorinated ethenes and benzene was investigated. The contamination originating from a former pharmaceutical industry discharges into a multilayered aquifer system and a downgradient stream. In this study, geological and hydrogeological data were integrated into a numerical flow model to examine identified trends using statistical approaches, including principal component analysis and hierarchal cluster analysis. A joint interpretation of the groundwater flow paths and contaminant concentrations in the different compartments (i.e., groundwater and hyporheic zone) provided insight on the transport processes of the different contaminant plumes to the stream. The analysis of historical groundwater concentrations of pharmaceutical compounds at the site suggested these compounds are slowly degrading. The pharmaceutical compounds migrate in both a deep semiconfined aquifer, as well as in the shallow unconfined aquifer, and enter the stream along a 2-km stretch. This contrasted with the chlorinated ethenes, which mainly discharge to the stream as a focused plume from the unconfined aquifer. The integrated approach developed here, combining groundwater flow modeling and statistical analyses of the contaminant concentration data collected in groundwater and the hyporheic zone, lead to an improved understanding of the observed distribution of contaminants in the unconfined and semiconfined aquifers, and thus to their discharge to the stream. This approach is particularly relevant for large and long-lasting contaminant sources and plumes, such as abandoned landfills and industrial production sites, where field investigations may be very expensive.     

Combination of CFCs and stable isotopes to characterize the mechanism of groundwater–surface water interactions in a headwater basin of the North China Plain

Mountainous areas are characterized by steep slopes and rocky landforms, with hydrological conditions varying rapidly from upstream to downstream, creating variable interactions between groundwater and surface water. In this study, mechanisms of groundwater–surface water interactions within a headwater catchment of the North China Plain were assessed along the stream length and during different seasons, using hydrochemical and stable isotope data, and groundwater residence times estimated using chlorofluorocarbons. These tracers indicate that the river is gaining, due to groundwater discharge in the headwater catchment both in the dry and rainy seasons. Residence time estimation of groundwater using chlorofluorocarbons data reveals that groundwater flow in the shallow sedimentary aquifer is dominated by the binary mixing of water approximating a piston flow model along 2 flow paths: old water, carried by a regional flow system along the direction of river flow, along with young water, which enters the river through local flow systems from hilly areas adjacent to the river valley (particularly during the rainy season). The larger mixing ratio of young water from lateral groundwater recharge and return flow of irrigation during the rainy season result in higher ion concentrations in groundwater than in the dry season. The binary mixing model showed that the ratio of young water versus total groundwater ranged from 0.88 to 0.22 and 1.0 to 0.74 in the upper and lower reaches, respectively. In the middle reach, meandering stream morphology allows some loss of river water back into the aquifer, leading to increasing estimates of the ratio of young water (from 0.22 to 1). This is also explained by declining groundwater levels near the river, due to groundwater extraction for agricultural irrigation. The switch from a greater predominance of regional flow in the dry season, to more localized groundwater flow paths in the wet season is an important groundwater–surface water interactions mechanism, with important catchment management implications.     

Slopes of an airborne electromagnetic resistivity model interpolated jointly with borehole data for 3D geological modelling

We investigate a novel way to introduce resistivity models deriving from airborne electromagnetic surveys into regional geological modelling. Standard geometrical geological modelling can be strengthened using geophysical data. Here, we propose to extract information contained in a resistivity model in the form of local slopes that constrain the modelling of geological interfaces. The proposed method is illustrated on an airborne electromagnetic survey conducted in the region of Courtenay in France. First, a resistivity contrast corresponding to the clay/chalk interface was interpreted confronting the electromagnetic soundings to boreholes. Slopes were then sampled on this geophysical model and jointly interpolated with the clay/chalk interface documented in boreholes using an implicit 3D potential‐field method. In order to evaluate this new joint geophysical–geological model, its accuracy was compared with that of both pure geological and pure geophysical models for various borehole configurations. The proposed joint modelling yields the most accurate clay/chalk interface whatever the number and location of boreholes taken into account for modelling and validation. Compared with standard geological modelling, the approach introduces in between boreholes geometrical information derived from geophysical results. Compared with conventional resistivity interpretation of the geophysical model, it reduces drift effects and honours the boreholes. The method therefore improves what is commonly obtained with geological or geophysical data separately, making it very attractive for robust 3D geological modelling of the subsurface.     

Groundwater recharge estimation in arid hardrock‐alluvium aquifers using combined water‐table fluctuation and groundwater balance approaches

This paper proposes an approach to estimate groundwater recharge using an optimization‐based water‐table fluctuation method combined with a groundwater balance model in an arid hardrock‐alluvium region, located at the Oman–United Arab Emirates border. We introduce an “effective hardrock thickness” term to identify the percentage of the considered hardrock thickness in which effective groundwater flow takes place. The proposed method is based upon a Thiessen polygon zoning approach. The method includes subpolygons to represent specific geologic units and to enhance the confidence of the estimated groundwater recharge. Two linear and 1 nonlinear submodels were developed to evaluate the model components for the calibration (October 1996 to September 2008) and validation (October 2008 to September 2013) periods. Long‐term annual groundwater recharge from rainfall and return flow over the model domain are estimated as 24.62 and 5.71 Mm³, respectively, while the effective groundwater flow circulation is found to occur in the upper 7% of the known hardrock thickness (42 m), confirming conclusions of previous field studies. Considering a total difference in groundwater levels between eastern and western points of the study area of the order of 220 m and a 12‐year monthly calibration period, a weighted root mean squared error in predicted groundwater elevation of 2.75 m is considered quite reasonable for the study area characterized by remarkable geological and hydrogeological diversity. The proposed approach provides an efficient and robust method to estimate groundwater recharge in regions with a complex geological setting in which interaction between fractured and porous media cannot be easily assessed.     

Groundwater recharge and discharge in a hyperarid alluvial plain (Akesu,Taklimakan Desert,China)

Groundwater recharge and discharge in the Akesu alluvial plain were estimated using a water balance method. The Akesu alluvial plain (4842 km²) is an oasis located in the hyperarid Tarim River basin of central Asia. The land along the Akesu River has a long history of agricultural development and the irrigation area is highly dependent on water withdrawals from the river. We present a water balance methodology to describe (a) surface water and groundwater interaction and (b) groundwater interaction between irrigated and non‐irrigated areas. Groundwater is recharged from the irrigation system and discharged in the non‐irrigated area. Uncultivated vegetation and wetlands are supplied from groundwater in the hyperarid environment. Results show that about 90% of groundwater recharge came from canal loss and field infiltration. The groundwater flow from irrigated to non‐irrigated areas was about 70% of non‐irrigated area recharge and acted as subsurface drainage for the irrigation area. This desalinated the irrigation area and supplied water to the non‐irrigated area. Salt moved to the non‐irrigation area following subsurface drainage. We conclude that the flooding of the Akesu River is a supplemental groundwater replenishment mechanism: the river desalinates the alluvial plain by recharging fresh water in summer and draining saline regeneration water in winter. Copyright © 2007 John Wiley & Sons, Ltd.