Origin of sulphate in the unsaturated zone and groundwater of a loess aquifer

The use of the sulphate mass balance (SMB) between precipitation and soil water as a supplementary method to estimate the diffuse recharge rate assumes that the sulphate in soil water originated entirely from atmospheric deposition; however, the origin of sulphate in soil and groundwater is often unclear, especially in loess aquifers. This study analysed the sulphur (δ³⁴S-SO₄) and oxygen (δ¹⁸O-SO₄) isotopes of sulphate in precipitation, water-extractable soil water, and shallow groundwater samples and used these data along with hydrochemical data to determine the sources of sulphate in the thick unsaturated zone and groundwater of a loess aquifer. The results suggest that sulphate in groundwater mainly originated from old precipitation. When precipitation percolates through the unsaturated zone to recharge groundwater, sulphates were rarely dissolved due to the formation of CaCO₃ film on the surface of sulphate minerals. The water-extractable sulphate in the deep unsaturated zone (>10 m) was mainly derived from the dissolution of evaporite minerals and there was no oxidation of sulphide minerals during the extraction of soil water by elutriating soil samples with deionized water. The water-extractable concentration of SO₄ was not representative of the actual SO₄ concentration in mobile soil water. Therefore, the recharge rate cannot be estimated by the SMB method using the water-extractable concentration of SO₄ in the loess areas. This study is important for identifying sulphate sources and clarifying the proper method for estimating the recharge rate in loess aquifers.     

Clogging of saturated porous media by silt‐sized suspended solids under varying physical conditions during managed aquifer recharge

Managed aquifer recharge is an effective method for utilizing excess flood flows, but clogging of porous media is a limiting factor in the implementation of this water storage technique. In recent years, much research on the physical clogging of porous media during artificial recharge has been conducted. However, the understanding of clogging due to silt‐sized suspended solids (SS) is still inadequate, especially under varying physical conditions. Here, we subjected sand columns to controlled rates of flow and SS suspensions to investigate the influence of media size, SS size, SS concentration, and flow velocity on the clogging of porous media by silt‐sized SS. The results show that the diameter ratio of SS particles to sand grains is the dominant factor influencing the position of physical clogging. As pore velocity increased, the mobility of silt‐sized SS was enhanced and retention in the porous media decreased noticeably. The spatial retention profiles in the porous media were found to vary greatly at different flow velocities. The SS concentration of the infiltrating suspension also dramatically influenced the mobility and deposition of silt‐sized SS particles, such that high concentrations accelerated the clogging process. As the different physical factors changed, the breakthrough curves and retention profiles of silt‐sized SS particles changed obviously and the mechanisms of retention differed. On the whole, clogging position is mainly determined by particle size ratio, but clogging rate is dominated by a variety of factors including particle size ratio, SS concentration, and flow velocity.     

Quantifying terrain controls on runoff retention and routing in the Northern Prairies

The role of hummocky terrain in governing runoff routing and focussing groundwater recharge in the Northern Prairies of North America is widely recognised. However, most hydrological studies in the region have not effectively utilised information on the surficial geology and associated landforms in large-scale hydrological characterization. The present study uses an automated digital elevation model (DEM) analysis of a 6500-km² area in the Northern Prairies to quantify hydrologically relevant terrain parameters for the common types of terrains in the prairies with different surficial deposits widespread in the prairies, namely, moraines and glaciolacustrine deposits. Runoff retention (and storage) capacity within depressions varies greatly between different surficial deposits and is comparable in magnitude with a typical amount of seasonal snowmelt runoff generation. The terrain constraint on potential runoff retention varies from a few millimetres in areas classified as moraine to tens of millimetres in areas classified as stagnant ice moraine deposits. Fluted moraine and glaciolacustrine deposits have intermediate storage capacity values. The study also identified the probability density function describing a number of immediate upstream neighbours for each depression in a fill-and-spill network. A relationship between depression parameters and surficial deposits, as well as identified depression network structure, allows parametrisation of hydrologic models outside of the high-resolution DEM coverage, which can still account for terrain variation in the Prairies.     

A comparison of methods to estimate groundwater recharge from bare soil based on data observed by a large-scale lysimeter

Effectively estimating groundwater recharge is critical to manage water resources, especially in arid and semi-arid regions as impacted by intensive human activities and climate changes. Rare insights have been gained into groundwater recharge since direct observation is hard to carry out. Although several methods are currently available to estimate groundwater recharge, the estimated results may cover noticeable bias. The behaviours of different methods based on different conceptual frameworks and exhibiting different levels of complexity should be examined to estimate actual groundwater recharge. This study aims to assess the performance of four common methods to estimate groundwater recharge. For this end, large-scale lysimeters equipped with soil water content sensors and water table sensors were set up at a research site established in Guanzhong Basin of China. The data achieved by 1-year observation were employed to compare four estimation methods. As revealed from the results, the following findings are drawn. (a) Groundwater level fluctuation (GLF) method is simple, whereas its accuracy is determined by specific yield, and adopting a water balance method to estimate specific yield can considerably enhance the accuracy of GLF. (b) The calibrated numerical model can obtain the optimal result compared with the other methods, whereas long-term observation data are required for parameter calibration. (c) In the water balance method, the maximum entropy production (MEP) model and a practical method (estimating evaporation between two rainfall events) were used to calculate evaporation. As indicated by the results, water balance method combined with MEP is capable of obtaining more reliable results of groundwater recharge compared with the practical method. (d) With an analytical model based on linearized Richards' equation, accurate results can be achieved. What is more, the analytical model only needs the measurement of soil moisture near the surface. The limitation of this method is that it is difficult to determine the maximal water flux. The mentioned findings are of critical implications to the management and sustainable development of groundwater.     

Estimating groundwater recharge for a freshwater lens in an arid region: Formative and stability assessment

The formation of subsurface freshwater lenses on top of brackish groundwater is a fascinating hydrologic phenomenon that creates groundwater supplies of great potential value in arid regions. Information on the recharge quantity and mechanism of these lenses is both scarce and uncertain. This study examines the formation and macroscale stability of the Rawdatain freshwater lens in Kuwait, for which significant pre-development data are available. The Rawdatain is a large (150 million m³) subsurface freshwater lens overlying brackish groundwater compared to the other freshwater lenses in the Arabian Peninsula. In this study, a three-dimensional (3D) density-dependent groundwater flow model is tested against the following data targets to estimate long-term diffuse and focused groundwater recharge: (i) groundwater head, (ii) total dissolved solids (TDS) groundwater concentration, (iii) volume and vertical thickness of stored groundwater of three different water quality TDS ranges (0–700, 700–1000 and 1000–2000 mg/L) and (iv) geometrical shape features of the lens along cross-sections. To better represent the spatial variation in TDS, six different recharge zones were assigned to allocate diffuse and focused recharge conditions. Twelve recharge rate scenarios, encompassing a wide range of feasible long-term average annual recharge values (200,000–5,000,000 m³/year), were tested against the multiple targets and compared with the groundwater age of the Rawdatain lens. Based on comparison with data targets, the long-term average annual recharge is estimated to be 500,000 m³/year. Scenarios of reduced recharge, which may occur due to changes in land-use or climate, demonstrate the extremely slow response of the lens, which is in agreement with the slow development and formation of the lens (>2,000 years). Within a 100-year time frame, a 50% reduction in annual recharge reduces the lens volumes by 21, 17 and 9% for the three water quality categories, respectively. This study demonstrates the stability of freshwater lenses in arid regions and also provides methodology for similar focused rainfall recharge freshwater lenses.     

Using chlorofluorocarbons (CFCs) and tritium (3H) to estimate groundwater age and flow velocity in Hohhot Basin,China

The concentrations of chlorofluorocarbons (CFC‐11, CFC‐12 and CFC‐113) and tritium (³H) content in groundwater were used to date groundwater age, delineate groundwater flow systems and estimate flow velocity in the Hohhot basin. The estimated young groundwater age is fallen in the bracket of 21 ~ 50 a and indicates the presence of two different age profiles and flow systems in the shallow groundwater system. Older age waters occur under the topographically low areas, where the aquifer is double‐layer aquifer system consisting of shallow unconfined‐semi‐confined aquifer and deep confined aquifer. This reflects long flow paths associated with regional flow. Groundwater (range from 21 to 34 years) in the north piedmont and east hilly areas, where the aquifer is a single‐layer aquifer consisting of alluvial fans, are typically younger than those in the low areas. The combination of CFCs dating with hydrogeological information indicates that both local and regional flow systems are present at the basin. The regional groundwater flow mainly flows from the north and east to the southwest, the local groundwater flow system occurs nearby the Hohhot city. The mean regional groundwater flow velocity of the shallow groundwater is estimated about 0.73 km/a. These findings can aid in refining hydrogeological conceptual model of the study area. Copyright © 2012 John Wiley & Sons, Ltd.     

Stable isotopes of water show deep seasonal recharge in northern bogs and fens

Ground water recharge is assumed to occur primarily at raised bog crests in northern peatlands, which are globally significant terrestrial carbon reservoirs. We synoptically surveyed vertical profiles of peat pore water δ¹⁸O and δ²H from a range of bog and fen landforms across the Glacial Lake Agassiz Peatlands, northern Minnesota. Contrary to our expectations, we find that local‐scale recharge penetrates to not only the basal peat at topographically high bog crests but also transitional Sphagnum lawns and low‐lying fen water tracks. Surface landscape characteristics appear to control the isotopic composition of the deeper pore waters (depths ≥0.5 m), which are partitioned into discrete ranges of δ¹⁸O on the basis of landform type (mean ± standard deviation for bog crests = ?11.9 ± 0.4‰, lawns = ?10.6 ± 0.1‰, fen water tracks = ?8.8 ± 1.0‰). Fen water tracks have a shallow free‐water surface that is seasonally enriched by isotope fractionating evaporation, fingerprinting recharge to underlying pore waters at depths ≥3 m. Isotope mass balance calculations indicate on average 12% of the waters we sampled from the basal peat of the fen water tracks was lost to surface evaporation, which occurred prior to advection and dispersion into the underlying formation. These new data provide direct support for the hypothesis that methane production in deeper peat strata is fuelled by the downward transport of labile carbon substrates from the surface of northern peat basins. Copyright © 2013 John Wiley & Sons, Ltd.     

Near‐surface hydraulic conductivity of northern hemisphere glaciers

The hydrology of near‐surface glacier ice remains a neglected aspect of glacier hydrology despite its role in modulating meltwater delivery to downstream environments. To elucidate the hydrological characteristics of this near‐surface glacial weathering crust, we describe the design and operation of a capacitance‐based piezometer that enables rapid, economical deployment across multiple sites and provides an accurate, high‐resolution record of near‐surface water‐level fluctuations. Piezometers were employed at 10 northern hemisphere glaciers, and through the application of standard bail–recharge techniques, we derive hydraulic conductivity (K) values from 0.003 to 3.519 m day^?1, with a mean of 0.185 ± 0.019 m day^?1. These results are comparable to those obtained in other discrete studies of glacier near‐surface ice, and for firn, and indicate that the weathering crust represents a hydrologically inefficient aquifer. Hydraulic conductivity correlated positively with water table height but negatively with altitude and cumulative short‐wave radiation since the last synoptic period of either negative air temperatures or turbulent energy flux dominance. The large range of K observed suggests complex interactions between meteorological influences and differences arising from variability in ice structure and crystallography. Our data demonstrate a greater complexity of near‐surface ice hydrology than hitherto appreciated and support the notion that the weathering crust can regulate the supraglacial discharge response to melt production. The conductivities reported here, coupled with typical supraglacial channel spacing, suggest that meltwater can be retained within the weathering crust for at least several days. Not only does this have implications for the accuracy of predictive meltwater run‐off models, but we also argue for biogeochemical processes and transfers that are strongly conditioned by water residence time and the efficacy of the cascade of sediments, impurities, microbes, and nutrients to downstream ecosystems. Because continued atmospheric warming will incur rising snowline elevations and glacier thinning, the supraglacial hydrological system may assume greater importance in many mountainous regions, and consequently, detailing weathering crust hydraulics represents a research priority because the flow path it represents remains poorly constrained.     

黄土高原丘陵区小流域地下水补给特征

地下水补给量反映了含水层的可更新能力,是地下水资源管理与合理开发利用的关键参数之一。为定量研究黄土高原丘陵沟壑区小流域地下水的补给特征,基于岔巴沟流域上游、中游、下游3个水文站1959-1969年降水与日径流观测资料,通过退水分析法估算了流域地下水补给量,并分析了与降水量和基流量的关系及其在年内的补给变化过程。结果表明：岔巴沟流域年平均基流量为13.09 mm·a^-1,更新时间为124 d,补给量为11.46 mm·a^-1,降水入渗补给率为0.025,基流补给率为0.89。从上游到下游地下水补给量与入渗补给率逐渐增大,且上游与下游之间有显著差异（P<0.05）;基流补给率逐渐减小,各集水区之间差异均显著。地下水补给量与降水量呈线性正相关（R²>0.40）,在下游集水区内随降水量变化的增幅较大。基流量与降水量也呈正相关关系（R²>0.77）,干流基流80%以上源于降水补给转化。以5月份为节点可将地下水补给过程分为"一次补给"和"二次补给"2个主要阶段,其分别占全年总补给量约30%和70%,并且"二次补给"是造成岔巴沟流域不同集水区地下水补给量差异的主要阶段,并且为无资料地区小流域地下水资源的评价提供借鉴。     

Modelling of water table variation in response to time-varying recharge from multiple basins using the linearised Boussinesq equation

A mathematical model is presented to describe the variations of the water table in an unconfined aquifer due to time-varying recharge applied from four rectangular basins. The model is developed by solving the linearised Boussinesq equation using the extended finite Fourier cosine transform. The time-varying recharge rate is approximated by a number of piecewise linear elements of different lengths and slopes depending on the nature of the variation in recharge rate. Application of this model for the prediction of water table fluctuations and in the sensitivity analysis of various controlling parameters on the aquifer response is demonstrated in an example.