Simulations of water movement and solute transport through different soil texture configurations under negative‐pressure irrigation

This study examined the effects of different soil texture configurations on water movement and solute transport to provide a reliable scientific basis for the application of negative‐pressure irrigation (NPI) technology. HYDRUS‐2D was used to analyse water movement and solute transport under NPI. The main results are as follows: (a) HYDRUS‐2D can be used to simulate water movement and solute transport under NPI, as there was good agreement between the simulated and measured values for water contents, NaCl concentrations, cumulative water infiltration, and wetting distances in the horizontal and vertical directions; the Nash–Sutcliffe efficiency coefficients were in the range of 0.94–0.97. (b) Layered soils have obvious effects on water movement under NPI. With the emitter position in the loam layer, when a coarse texture of loamy sand was present below the loam layer (namely, L‐LS), irrigation water accumulated in the topsoil, and this led to an increase in evaporation compared with the homogeneous loam profile. However, fine texture silty loam or silty clay loam layers beneath the loam layer (namely, L‐SiL or L‐SiCL, respectively) was more conducive to water infiltration into the lower layer, and this increased the amount of water infiltration and simultaneously reduced the surface evaporation effectively. (c) Layered soils have obvious effects on solute transport under NPI, and salt accumulation will readily occur in the clay‐rich soil layer at the interface. The maximum soil salt accumulation of L‐LS occurred above the soil interface between the two soil layers with a value of 21.80 g/kg; however, for L‐SiCL and L‐SiL, the maximum salt accumulation occurred below the soil interface between the two soil layers, with values of 23.80 g/kg and 20.08 g/kg, respectively. (d) Interlayered soils showed remarkable changes in the water infiltration characteristics and salt‐leaching intensities under NPI, and the properties for the soil profile with a silty loam interlayer were better than those for the soil profile with a silty clay loam interlayer. The soil profile with a loamy sand interlayer had the lowest amount of water infiltration, which resulted in reductions of the salt‐leaching intensities. Thus, NPI is clearly not suitable for loamy sand soil. Overall, the results demonstrated that soil texture configurations affected water movement and solute transport under NPI. Therefore, careful consideration should be given to the use of NPI to achieve target soil water and solution conditions and reduce water loss.     

Agricultural management impacts on groundwater: simulations of existing and alternative management options in Peninsular India

Understanding the principal causes and possible solutions for groundwater depletion in India is important for its water security, especially as it relates to agriculture. A study was conducted in an agricultural watershed in Andhra Pradesh, India to assess the impacts on groundwater of current and alternative agricultural management. Hydrological simulations were used as follows: (1) to evaluate the recharge benefits of water‐harvesting tillage through a modified Soil and Water Assessment Tool (SWAT) model and (2) to predict the groundwater response to changing extent and irrigation management of rice growing areas. The Green–Ampt infiltration routine was modified in SWAT was modified to represent water‐harvesting tillage using maximum depression storage parameter. Water‐harvesting tillage in rainfed croplands was shown to increase basin‐scale groundwater recharge by 3% and decrease run‐off by 43% compared with existing conventional tillage. The groundwater balance (recharge minus irrigation withdrawals), negative 11 mm/year under existing management changed to positive (18–45 mm/year) when rice growing areas or irrigation depths were reduced. Groundwater balance was sensitive to changes in rice cropland management, meaning even small changes in rice cropland management had large impacts on groundwater availability. The modified SWAT was capable of representing tillage management of varying maximum depression storage, and tillage for water‐harvesting was shown to be a potentially important strategy for producers to enhance infiltration and groundwater recharge, especially in semi‐arid regions where rainfall may be becoming increasingly variable. This enhanced SWAT could be used to evaluate the landscape‐scale impacts of alternative tillage management in other regions that are working to develop strategies for reducing groundwater depletion. Copyright © 2013 John Wiley & Sons, Ltd.     

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.     

Impact of water transfer on interaction between surface water and groundwater in the lowland area of North China Plain

The contradiction between the freshwater shortage and the large demand of freshwater by irrigation was the key point in cultivated lowland area of North China Plain. Water transfer project brings fresh water from water resource‐rich area to water shortage area, which can in turn change the hydrological cycle in this region. Major ions and stable isotopes were used to study the temporal variations of interaction between surface water and groundwater in a hydrological year after a water transfer event in November 2014. Irrigation canal received transferred Yellow River, with 2.9% loss by evaporation during water transfer process. The effect of transferred water on shallow groundwater decreased with increasing distance from the irrigation canal. Pit pond without water transfer receives groundwater discharge. During dry season after water transfer event, shallow groundwater near the irrigation canal was recharged by lateral seepage and deep percolation of irrigation, whereas shallow groundwater far from irrigation canal was recharged by deep percolation of deep groundwater irrigation. Canal water lost by evaporation was 2.7–17.4%. Influence of water transfer gradually disappeared until March as the water usage of agricultural irrigation increased. In the dry season, groundwater discharged to irrigation canal and pond; 2.2–31.6% canal water and 11.3–20.0% pond water were lost by evaporation. In the rainy season (June to September), surface water was fed mainly by precipitation and surface run‐off, whereas groundwater was recharged by infiltration of precipitation. The two‐end member mix model showed that the mixing ratio of precipitation in pond and irrigation canal were 73–83.4% (except one pond with 28.1%) and 77.3–99.9%, respectively. Transferred water and precipitation were the important recharge sources for shallow groundwater, which decreased groundwater salinity in cultivated lowland area of North China Plain. With the temporary and spatial limitation of water transfer effects, increased water transfer amounts and frequency may be an effective way of mitigating regional water shortage. In addition, reducing the evaporation of surface water is also an important way to increase the utilization of transfer water.     

A soil‐water‐balance approach to quantify groundwater recharge from irrigated cropland in the North China Plain

Rapidly depleting unconfined aquifers are the primary source of water for irrigation on the North China Plain. Yet, despite its critical importance, groundwater recharge to the Plain remains an enigma. We introduce a one‐dimensional soil‐water‐balance model to estimate precipitation‐ and irrigation‐generated areal recharge from commonly available crop and soil characteristics and climate data. To limit input data needs and to simplify calculations, the model assumes that water flows vertically downward under a unit gradient; infiltration and evapotranspiration are separate, sequential processes; evapotranspiration is allocated to evaporation and transpiration as a function of leaf‐area index and is limited by soil‐moisture content; and evaporation and transpiration are distributed through the soil profile as exponential functions of soil and root depth, respectively. For calibration, model‐calculated water contents of 11 soil‐depth intervals from 0 to 200 cm were compared with measured water contents of loam soil at four sites in Luancheng County, Hebei Province, over 3 years (1998–2001). Each 50‐m² site was identically cropped with winter wheat and summer maize, but received a different irrigation treatment. Average root mean‐squared error between measured and model‐calculated water content of the top 180 cm was 4·2 cm, or 9·3% of average total water content. In addition, model‐calculated evapotranspiration compared well with that measured by a large‐scale lysimeter. To test the model, 12 additional sites were simulated successfully. Model results demonstrate that drainage from the soil profile is not a constant fraction of precipitation and irrigation inputs, but rather the fraction increases as the inputs increase. Because this drainage recharges the underlying aquifer, improving irrigation efficiency by reducing seepage will not reverse water‐table declines. Copyright © 2003 John Wiley & Sons, Ltd.     

Soil macroporosity and infiltration characteristics of a forest podzol

Spatial distribution of soil macroporosity was determined for a forest podzol from tension infiltrometer measurements at the soil surface. Surface‐derived macroporosity values were compared with point infiltration characteristics obtained from soil water content and soil water chemistry measurements during an experimental irrigation, and with parameters of a kinematic wave model applied to soil water content data. Macroporosity estimated by the tension infiltrometer ranged from 0·00087 to 0·0219% of soil volume, and infiltration at these two sites was dominated by propagation of a well‐defined wetting front through the soil profile and bypass flow via soil macropores, respectively. Infiltration at sites with intermediate macroporosities reflected a combination of these two processes, although results were inconclusive at one site owing to lateral flow at the base of the soil profile. There was no agreement between macroporosities estimated by the tension infiltrometer and the kinematic wave model. The maximum soil conductance parameter within the profile at a site, however, was related directly to the surface‐derived macroporosity. The partial agreement between surface‐derived macroporosity estimates and point infiltration characteristics shown here supports the use of tension infiltrometry as a rapid, non‐destructive method of assessing spatial variations in the relative contribution of macropore flow to the infiltration process. Copyright © 2000 John Wiley & Sons, Ltd.     

Genesis of arsenic in groundwater of North Bengal Plain using PCA: A case study of English Bazar Block,Malda District,West Bengal,India

The study area is located on the western part of the alluvium‐filled gap between the Rajmahal hills on the west and the Garo hills on the east. Groundwater occurs under unconfined condition in a thick zone of saturation within the Quaternary alluvial sediments. Three hydrochemical facies with distinct characteristics have been identified which are dominated in general by alkaline earths and weak acids. The major‐ion chemistry of the area is controlled by weathering of silicate minerals, rainfall recharge, ion‐exchange processes and anthropogenic activities such as irrigation return flow and the application of inorganic fertilizers and pesticides. A stoichiometric approach suggests that mineral dissolution and anthropogenic activities contribute 79% and 21% of the total cations dissolved in groundwater. Principal component analysis (PCA) of 42 groundwater samples using 13 chemical parameters indicates that the combined processes of recharge of groundwater from rainfall, sediment water interaction, groundwater flow, infiltration of irrigation return water (which is arsenic rich due to the use of arsenic‐bearing pesticides, wood preservatives, etc. and the pumping of arsenic‐rich groundwater for agriculture purpose), oxidation of natural or anthropogenic organic matter and the reductive dissolution of ferric iron and manganese oxides play a key role in the evolution of groundwater in the study area. Factor 2 scores, associated with the infiltration of irrigation return water and spatial distribution of arsenic concentration reveal that the groundwater of the municipal area will not be affected by arsenic in the future in spite of heavy groundwater abstraction. Another PCA with geologic, geomorphic, anthropogenic, geochemical and landuse factors indicates that arsenic concentration in groundwater increases with increasing area of mango orchards, sand lithofacies and nitrate and decreases with increasing distance of paleochannel from the monitored well and depth of bore wells. High loading on nitrate may be attributed to the use of fertilizer, pesticides, etc. in mango orchards and agricultural land. High loadings on log pCO₂, mango orchards (with negative sign) and phosphate (with positive sign) indicate that mango orchards provide the organic waste material which is decomposed to form organic carbon. The organic carbon undergoes oxidative carbon degeneration by different oxidants and increases the concentration of CO₂ in the aquifer. The reducing condition thus developed in the aquifer helps to dissolve the arsenic adsorbed on iron hydroxide or oxy‐hydroxide coated margins of sand, iron rich heavy mineral grain margins, clay minerals and Fe–Mn concretions present in the aquifer matrix. Copyright © 2007 John Wiley & Sons, Ltd.     

Fast modelling of marine surface multiples

In the development and testing of water‐surface multiple‐removal algorithms, it is valuable to have accurate synthetic seismograms which exhibit multiples, for which the multiple‐free solution is known. A method is presented for constructing 2D and 3D solutions of the acoustic wave equation in water, by combining the solution from a primary source with other scaled solutions of secondary sources, which simulate diffractors. The computation involves function evaluation rather than numerical solution of differential equations and is consequently accurate and comparatively fast. The analytic formulae on which the method is based give insights into methods for multiple removal. Generalized reflection coefficients, defined on a horizontal plane above the diffractors, are derived and used to construct the integral equations which are the basis for many multiple‐removal schemes.     

Importance of river water recharge to the San Joaquin Valley groundwater system

Groundwater is not a sustainable resource, unless abstraction is balanced by recharge. Identifying the sources of recharge in a groundwater basin is critical for sustainable groundwater management. We studied the importance of river water recharge to groundwater in the south‐eastern San Joaquin Valley (24,000 km², population 4 million). We combined dissolved noble gas concentrations, stable isotopes, tritium, and carbon‐14 analyses to analyse the sources, mechanisms, and timescales of groundwater recharge. Area‐representative groundwater sampling and numerical model input data enabled a stable isotope mass balance and quantitative estimates of river and local recharge. River recharge, identified by a lighter stable isotope signature, represents 47 ± 4% of modern groundwater in the San Joaquin Valley (recharged after 1950) but only 26 ± 4% of premodern groundwater (recharged before 1950). This implies that the importance of river water recharge in the San Joaquin valley has nearly doubled and is likely the result of a 40% increase in total recharge, caused by river water irrigation return flows and increased stream depletion and river recharge due to groundwater pumping. Compared with the large and long‐duration capacity for water storage in the subsurface, storage of water in rivers is limited in time and volume, as evidenced by cold river recharge temperatures resulting from fast infiltration and recharge. Groundwater banking of seasonal surface water flows and expansion of managed aquifer recharge practices therefore appear to be a natural and promising method for increasing the resilience of the San Joaquin Valley water supply system.     

Hydrogeological and hydrogeochemical control of groundwater salinity in an arid inland basin: Dunhuang Basin,northwestern China

High groundwater salinity has become a major concern in the arid alluvial plain of the Dunhuang Basin in northwestern China because it poses a significant challenge to water resource management. Isotopic and geochemical analyses were conducted on 55 water samples from springs, boreholes and surface water to identify potential sources of groundwater salinity and analyse the processes that control increasing salinity. The total dissolved solid (TDS) content in the groundwater ranged from 400 to 41 000 mg/l, and high TDS values were commonly associated with shallow water tables and flow‐through and discharge zones in unconfined aquifers. Various groundwater contributions from rainwater, agricultural irrigation, river water infiltration and lateral inflows from mountains were identified by major ions and δD and δ¹⁸O. In general, HCO₃^? and SO₄^2? were the dominant anions in groundwater with a salinity of <2500 mg/l, whereas Cl^? and SO₄^2? were the dominant anions in groundwater with a salinity of >2500 mg/l. The major ion concentrations indicated that mineral weathering, including carbonate and evaporite dissolution, primarily affected groundwater salinity in recharge areas. Evapotranspiration controlled the major ion concentration evolution and salinity distribution in the unconfined groundwaters in the flow‐through and discharge areas, although it had a limited effect on groundwater in the recharge areas and confined aquifers. Agricultural irrigation increased the water table and enhanced evapotranspiration in the oasis areas of the basin. TDS and Cl became more concentrated, but H and O isotopes were not enriched in the irrigation district, indicating that transpiration dominated the increasing salinity. For other places in the basin, as indicated by TDS, Cl, δD and δ¹⁸O characteristics, evaporation, transpiration and water–rock interactions dominated at different hydrogeological zones, depending on the plant coverage and hydrogeological conditions. Groundwater ages of ³H, and δD and δ¹⁸O compositions and distributions suggest that most of the groundwaters in Dunhuang Basin have a paleometeoric origin and experienced a long residence time. These results can contribute to groundwater management and future water allocation programmes in the Dunhuang Basin. Copyright © 2015 John Wiley & Sons, Ltd.     

The hydrogeomorphic influences on alluvial gully erosion along the Mitchell River fluvial megafan

Hydrogeomorphic processes influencing alluvial gully erosion were evaluated at multiple spatial and temporal scales across the Mitchell River fluvial megafan in tropical Queensland, Australia. Longitudinal changes in floodplain inundation were quantified using river gauge data, local stage recorders and HEC‐RAS modelling based on LiDAR topographic data. Intra‐ and interannual gully scarp retreat rates were measured using daily time‐lapse photographs and annual GPS surveys. Erosion was analysed in response to different water sources and associated erosion processes across the floodplain perirheic zone, including direct rainfall, infiltration‐excess runoff, soil‐water seepage, river backwater and overbank flood inundation. The frequency of river flood inundation of alluvial gullies changed longitudinally according to river incision and confinement. Near the top of the megafan, flood water was contained within the macrochannel up to the 100‐year recurrence interval, but river backwater still partially inundated adjacent gullies eroding into Pleistocene alluvium. In downstream Holocene floodplains, inundation of alluvial gullies occurred beyond the 2‐ to 5‐year recurrence interval and contributed significantly to total annual erosion. However, most gully scarp retreat at all sites was driven by direct rainfall and infiltration‐excess runoff, with the 24‐h rainfall total being the most predictive variable. The remaining variability can be explained by seasonal vegetative conditions, complex cycles of soil wetting and drying, tension crack development, near‐surface pore‐water pressure, soil block undermining from spalling and overland flow, and soil property heterogeneity. Implications for grazing management impacts on soil surface and perennial grass conditions include effects on direct rainfall erosion, water infiltration, runoff volume, water concentration along tracks, and the resistance of highly dispersible soils to gully initiation or propagation under intense tropical rainfall. Copyright © 2012 John Wiley & Sons, Ltd.     

A pattern-based approach for multiple removal applied to a 3D Gulf of Mexico data set

Surface‐related multiples are attenuated for one sail line and one streamer of a 3D data set (courtesy of Compagnie Générale de Géophysique). The survey was carried out in the Gulf of Mexico in the Green Canyon area where salt intrusions close to the water‐bottom are present. Because of the complexity of the subsurface, a wavefield method incorporating the full 3D volume of the data for multiple removal is necessary. This method comprises modelling of the multiples, where the data are used as a prediction operator, and a subtraction step, where the model of the multiples is adaptively removed from the data with matching filters. The accuracy of the multiple model depends on the source/receiver coverage at the surface. When this coverage is not dense enough, the multiple model contains errors that make successful subtraction more difficult. In these circumstances, one can either (1) improve the modelling step by interpolating the missing traces, (2) improve the subtraction step by designing methods that are less sensitive to modelling errors, or (3) both. For this data set, the second option is investigated by predicting the multiples in a 2D sense (as opposed to 3D) and performing the subtraction with a pattern‐based approach. Because some traces and shots are missing for the 2D prediction, the data are interpolated in the in‐line direction using a hyperbolic Radon transform with and without sparseness constraints. The interpolation with a sparseness constraint yields the best multiple model. For the subtraction, the pattern‐based technique is compared with a more standard, adaptive‐subtraction scheme. The pattern‐based approach is based on the estimation of 3D prediction‐error filters for the primaries and the multiples, followed by a least‐squares estimation of the primaries. Both methods are compared before and after prestack depth migration. These results suggest that, when the multiple model is not accurate, the pattern‐based method is more effective than adaptive subtraction at removing surface‐related multiples while preserving the primaries.     

How soil texture,channel shape and cross-sectional area affect moisture dynamics and water loss in irrigation channels

To enhance the utilization efficiency of farmland irrigation water and reduce the leakage of water conveyance channels, the leakage process of channels was simulated dynamically. The simulated results were compared with data measured in laboratory experiments, and the performance of the model was evaluated. The results indicated that the simulated values of the model were consistent with the observation values, and the R² values varied between 0.91 and 0.99. In addition, based on the laboratory experiments, a water supply system (Mariotte bottles) and soil box were built using plexiglass. Three influencing factors, namely, the channel form, soil texture and channel cross-sectional area, were varied to observe and calculate the resulting cumulative infiltration amount, infiltration rate and wetting front migration distance. HYDRUS-3D software was used to solve the three-dimensional soil water movement equation under different initial conditions. The results demonstrated that the U-shaped channel was more effective than the trapezoidal channel in increasing the utilization efficiency of the water resources. A U-shaped channel with a small channel cross-sectional area should be adopted and the soil particle size should be prioritized in the construction of water conveyance channels for farmlands. The simulation results were in agreement with the observed results, which indicates that HYDRUS-3D is a reliable tool that can accurately simulate the soil moisture movement in water conveyance channels. The research results can provide a reference for the design and operation of farmland irrigation systems.     

Portable irrigation system for studying hillslope and wetland runoff generation processes

The advantage of an irrigation system is that experiments of varying precipitation duration and intensity can be performed in a controlled situation. An effective portable irrigation system is described for use in experimental plot hillslope and wetland runoff studies. The system consists of four parts:(1) a water pump; (2) a tracer reservoir; (3) a chemical feed pump; and (4) distribution hoses. Relatively uniform water application is achieved by a series of manual valves controlling water flow from the main carrier hose to the distribution hoses. The irrigation system materials are inexpensive and installation and operation costs are minor. The system was used to study runoff generation from a small‐saturated area in a spring‐fed swamp for a range of precipitation intensities and durations. The irrigation system applied a maximum intensity of 14·0mm h^?1, for a maximum duration of 180 min, to a 190m² area. This range of application incorporated all storms up to a one in three year event. The variance of the tracer load was almost three times greater for natural (60%) than with the irrigation system (22%). The irrigation system reduced the uncertainty in the pre‐event water fraction (using a two‐component hydrograph separation) from 3·4% in natural events to only 0·7%. The irrigation system design, operation, calibration and cost are presented. Copyright © 2001 John Wiley & Sons, Ltd.     

Two‐dimensional physically based finite element runoff model for small agricultural watersheds: I. Model development

Soil erosion by water is the root cause of ecological degradation in the Shiwalik foothills of Northern India. Simulation of runoff and its component processes is a pre‐requisite to develop the management strategies to tackle the problem, successfully. A two‐dimensional physically based distributed numerical model, ROMO2D has been developed to simulate runoff from small agricultural watersheds on an event basis. The model employs the 2‐D Richards equation with sink term to simulate infiltration and soil moisture dynamics in the vadoze zone under variable rainfall conditions, and 2‐D Saint‐Venant equations under the kinematic wave approximation along with Manning's equation as the stage‐discharge equation for runoff routing. The various flow‐governing equations have been solved numerically by employing a Galerkin finite element method for spatial discretization using quadrilateral elements and finite difference techniques for temporal solutions. The ROMO2D computer program has been developed as a class‐based program, coded in C + + in such a way that with minor modifications, the model can be used to simulate runoff on a continuous basis. The model writes output for a runoff hydrograph of each storm. Model development is described in this paper and the results of model testing and field application are to be presented in a subsequent paper. Copyright © 2008 John Wiley & Sons, Ltd.     

A statistically based runoff‐yield model coupling infiltration excess and saturation excess mechanisms

A statistically based runoff‐yield model is proposed in this paper. The model considers spatial heterogeneities of rainfall, soil infiltration capacity and soil water storage capacity that are main factors controlling runoff‐yield process. It assumes that the spatial variation of rainfall intensity at each time step can be characterized by a probability density function, which is estimated by matching the hyetograph through goodness‐of‐fit measure, whereas the spatial heterogeneities of soil infiltration capacity and soil water storage capacity are described by parabola‐type functions. Surface runoff is calculated according to infiltration excess mechanism; the statistical distribution of surface runoff rate can be deduced with the joint distribution of rainfall intensity and soil infiltration rate, thus obtaining a quasi‐analytical solution for surface runoff. Based on saturation excess mechanism, the groundwater flow (flows below the ground are collectively referred to as groundwater flow) is calculated by infiltration and the probability distribution of soil water storage capacity. Consequently, the total runoff is composed of infiltration excess and saturation excess runoff components. As an example, this model is applied to flood event simulation in Dongwan catchment, a semi‐humid region and a tributary of Yellow River in China. It indicates that the proposed runoff‐yield model could achieve acceptable accuracy. Copyright © 2011 John Wiley & Sons, Ltd.     

Identification of groundwater recharge sources and processes in a heterogeneous alluvial aquifer: results from multi‐level monitoring of hydrochemistry and environmental isotopes in a riverside agricultural area in Korea

We evaluated sources and pathways of groundwater recharge for a heterogeneous alluvial aquifer beneath an agricultural field, based on multi‐level monitoring of hydrochemistry and environmental isotopes of a riverside groundwater system at Buyeo, Korea. Two distinct groundwater zones were identified with depth: (1) a shallow oxic groundwater zone, characterized by elevated concentrations of NO₃^? and (2) a deeper (>10–14 m from the ground surface) sub‐oxic groundwater zone with high concentrations of dissolved Fe, silica, and HCO₃^?, but little nitrate. The change of redox zones occurred at a depth where the aquifer sediments change from an upper sandy stratum to a silty stratum with mud caps. The δ¹⁸O and δ²H values of groundwater were also different between the two zones. Hydrochemical and δ¹⁸O? δ²H data of oxic groundwater are similar to those of soil water. This illustrates that recharge of oxic groundwater mainly occurs through direct infiltration of rain and irrigation water in the sandy soil area where vegetable cropping with abundant fertilizer use is predominant. Oxic groundwater is therefore severely contaminated by agrochemical pollutants such as nitrate. In contrast, deeper sub‐oxic groundwater contains only small amounts of dissolved oxygen (DO) and NO₃^?. The ³H contents and elevated silica concentrations in sub‐oxic groundwater indicate a somewhat longer mean residence time of groundwater within this part of the aquifer. Sub‐oxic groundwater was also characterized by higher δ¹⁸O and δ²H values and lower d‐excess values, indicating significant evaporation during recharge. We suggest that recharge of sub‐oxic groundwater occurs in the areas of paddy rice fields where standing irrigation and rain water are affected by strong evaporation, and that reducing conditions develop during subsequent sub‐surface infiltration. This study illustrates the existence of two groundwater bodies with different recharge processes within an alluvial aquifer. Copyright © 2009 John Wiley & Sons, Ltd.     

Analysing tension infiltrometer data from sloped surface using two‐dimensional approximation

Tension infiltrometers allow water to infiltrate into the soils at various specified pressure heads. The resulting infiltration rates can then be analysed for soil hydraulic properties by either analytical or inverse numerical methods. Tension infiltrometers however are primarily designed to be deployed on horizontal land surfaces, and their applications have been studied widely using an inverse numerical tool HYDRUS‐2D. However, natural landscapes are often nonhorizontal, and infiltration through tension infiltrometers on sloped surfaces is no longer an axisymmetrical two‐dimensional (2D) process but a fully three‐dimensional (3D) one. In addition, minimal research has examined the effect of simplifying the 3D problem to a 2D one on the hydraulic conductivity estimated using tension infiltrometer data from different land slopes of various soil types. Therefore, in this study, tension infiltrometer data on different slopes have been obtained from a catchment located at National University of Singapore. In addition, tension infiltrometer data of six soil types on different slopes and with different initial water content were simulated using HYDRUS‐3D. Combining field measurements, forward and inverse modelling, the influence of applying a 2D approximation on hydraulic property estimations using tension infiltrometer data was examined. The results show that the estimations for soils with high infiltration rates are more sensitive to application of the 2D approximation. The maximum allowable slopes for employing 2D approximation on clayey and sandy soils are 25° and less than 3°, respectively. Furthermore, the maximum allowable slope decreases with decreasing initial water content. Copyright © 2012 John Wiley & Sons, Ltd.     

Numerical modelling of stream–aquifer interaction: Quantifying the impact of transient streambed permeability and aquifer heterogeneity

Stream–aquifer interaction plays a vital role in the water cycle, and a proper study of this interaction is needed for understanding groundwater recharge, contaminants migration, and for managing surface water and groundwater resources. A model‐based investigation of a field experiment in a riparian zone of the Schwarzbach river, a tributary of the Rhine River in Germany, was conducted to understand stream–aquifer interaction under alternative gaining and losing streamflow conditions. An equivalent streambed permeability, estimated by inverting aquifer responses to flood waves, shows that streambed permeability increased during infiltration of stream water to aquifer and decreased during exfiltration. Aquifer permeability realizations generated by multiple‐point geostatistics exhibit a high degree of heterogeneity and anisotropy. A coupled surface water groundwater flow model was developed incorporating the time‐varying streambed permeability and heterogeneous aquifer permeability realizations. The model was able to reproduce varying pressure heads at two observation wells near the stream over a period of 55 days. A Monte Carlo analysis was also carried out to simulate groundwater flow, its age distribution, and the release of a hypothetical wastewater plume into the aquifer from the stream. Results of this uncertainty analysis suggest (a) stream–aquifer exchange flux during the infiltration periods was constrained by aquifer permeability; (b) during exfiltration, this flux was constrained by the reduced streambed permeability; (c) the effect of temporally variable streambed permeability and aquifer heterogeneity were found important to improve the accurate capture of the uncertainty; and (d) probabilistic infiltration paths in the aquifer reveal that such pathways and the associated prediction of the extent of the contaminant plume are highly dependent on aquifer heterogeneity.     

Drip irrigation enhances shallow groundwater contribution to crop water consumption in an arid area

Shallow groundwater plays a key role in agro‐hydrological processes of arid areas. Groundwater often supplies a necessary part of the water requirement of crops and surrounding native vegetation, such as groundwater‐dependent ecosystems. However, the impact of water‐saving irrigation on cropland water balance, such as the contribution of shallow groundwater to field evapotranspiration, requires further investigation. Increased understanding of quantitative evaluation of field‐scale water productivity under different irrigation methods aids policy and decision‐making. In this study, high‐resolution water table depth and soil water content in field maize were monitored under conditions of flood irrigation (FI) and drip irrigation (DI), respectively. Groundwater evapotranspiration (ET_g) was estimated by the combination of the water table fluctuation method and an empirical groundwater–soil–atmosphere continuum model. The results indicate that daily ET_g at different growth stages varies under the two irrigation methods. Between two consecutive irrigation events of the FI site, daily ET_g rate increases from zero to greater than that of the DI site. Maize under DI steadily consumes more groundwater than FI, accounting for 16.4% and 14.5% of ET_a, respectively. Overall, FI recharges groundwater, whereas DI extracts water from shallow groundwater. The yield under DI increases compared with that under FI, with less ET_a (526 mm) compared with FI (578 mm), and irrigation water productivity improves from 3.51 kg m^?3 (FI) to 4.58 kg m^?3 (DI) through reducing deep drainage and soil evaporation by DI. These results highlight the critical role of irrigation method and groundwater on crop water consumption and productivity. This study provides important information to aid the development of agricultural irrigation schemes in arid areas with shallow groundwater.