Groundwater recharge sources in semiarid irrigated mountain fronts

High-elevation mountains often constitute for basins important groundwater recharge sources through mountain-front recharge processes. These processes include streamflow losses and subsurface inflow from the mountain block. However, another key recharge process is from irrigation practices, where mountain streamflow is distributed across the irrigated piedmont. In this study, coupled groundwater fluctuation measurements and environmental tracers (¹⁸O, ²H, and major ions) were used to identify and compare the natural mountain-front recharge to the anthropogenically induced irrigation recharge. Within the High Atlas mountain front of the Ourika Basin, Central Morocco, the groundwater fluctuation mapping from the dry to wet season showed that recharge beneath the irrigated area was higher than the recharge along the streambed. Irrigation practices in the region divert more than 65% of the stream water, thereby reducing the potential for in-stream groundwater recharge. In addition, the irrigation areas close to the mountain front had greater water table increases (up to 3.5 m) compared with the downstream irrigation areas (<1 m increase). Upstream crops have priority to irrigation with stream water over downstream areas. The latter are only irrigated via stream water during large flood events and are otherwise supplemented by groundwater resources. These changes in water resources used for irrigation practices between upstream and downstream areas are reflected in the spatiotemporal evolution of the stable isotopes of groundwater. In the upstream irrigation area, the groundwater stable isotope values (δ¹⁸O: −8.4‰ to −7.4‰) reflect recharge by the diverted stream water. In the downstream irrigation area, the groundwater isotope values are lower (δ¹⁸O: −8.1‰ to −8.4‰) due to recharge via the flood water. In the nonirrigation area, the groundwater has the highest stable isotope values (δ¹⁸O: −6.8‰ to −4.8‰). This might be due to recharge via subsurface inflow from the mountain block to the mountain front and/or recharge via local low altitude rainfall. These findings highlight that irrigation practices can result in the dominant mountain-front recharge process for groundwater.     

Influence of a coarse interlayer on seawater intrusion and contaminant migration in coastal aquifers

Vertical 2D slice laboratory experiments were carried out in homogenous and layered sand tanks to elucidate the effects of a highly permeable (coarse‐grained sand) interlayer on seawater intrusion and transport of contaminants to a coastal sea. Tidal fluctuations produced oscillations in the seawater–freshwater transition zone, fluctuations of the contaminant infiltration rate and a zigzag contaminant plume outline. The seawater wedge became discontinuous at the (vertical) edges of the interlayer because of increased lateral movement of the seawater–freshwater interface within the interlayer. The contaminant plume formed a tail within the interlayer depending on the tidal stage, and similar to the wedge, its movement was accentuated. A simple analytical model that neglected vertical flow reliably predicted steady‐state seawater intrusion into the coastal aquifer. Numerical modeling was used to gain insight into the groundwater hydrodynamics and contaminant migration. The numerical results confirmed the experimental findings, i.e. that a highly permeable interlayer can provide a rapid transit path for contaminants to reach the seaward boundary and that the interlayer amplifies the effects of tidal fluctuations, resulting in wider transition zones for the seawater wedge and contaminant plume. Numerical simulations further showed that, with increasing interlayer hydraulic conductivity, the maximum seawater intrusion distance inside the interlayer increases approximately linearly. For the fixed‐head contaminant injection condition used, the model showed that contaminant infiltration increases approximately logarithmically with increasing interlayer hydraulic conductivity (other factors held fixed). Copyright © 2013 John Wiley & Sons, Ltd.     

A novel approach for estimating karst groundwater recharge in mountainous regions and its application in Switzerland

A pragmatic and simple approach for estimating the groundwater recharge of karst aquifers in mountainous regions by extrapolation of the hydrological regimes of gauged and well‐documented systems is presented. Specific discharge rates are derived using annual precipitation and spring measurements by taking into account catchment size and elevation, which are assumed to be the dominant factors. Reference sites with high data reliability are used for calibration and regional extrapolation. This is performed with normalized values employing spatial precipitation deviations and correlation with the elevation of the catchment areas. A tiered step procedure provides minimum and maximum normalized gradients for the relationship between recharge quantity and elevation for karst regions. The normalized recharge can therefore be obtained and extrapolated for any location using the spatial precipitation variability to provide an estimate of annual groundwater recharge. The approach was applied to Switzerland (approximately 7500 km² of karst terrain situated between 200 and over 4000 m a.s.l.) using annual precipitation data from meteorological stations for the years 2000 to 2011. Results show that the average recharge rates of different Swiss karst domains range from 20 to 46 L/km²s, which corresponds to an infiltration ratio between 0.6 and 0.9 of total precipitation. Despite uncertainties inherent in the approach, these results provide a benchmark for renewable karst groundwater resources in Switzerland of about 8.4 km³/year. The approach can be applied to any other mountainous karst region, that is, where a clear relationship between elevation, precipitation and recharge can be assumed. Copyright © 2015 John Wiley & Sons, Ltd.     

Complex hydraulic interactions between compartmentalized carbonate aquifers and heterogeneous siliciclastic successions: a case study in southern Italy

The sedimentary sequences containing lithologic units with low permeability represent hydrogeologic systems, which, as of now, have been little studied despite their diffusion worldwide. A hydrogeologic study, aimed to assess the main factors controlling the groundwater flow dynamics in such systems and their hydraulic interactions with nearby carbonate aquifers, has been carried out in Longano (Isernia, Southern Italy). The analysis of the hydraulic heads, combined with the regimes of the springs and the electric conductivity of the groundwater, mainly reflect vertical and lateral heterogeneities of the media in terms of hydraulic properties. In particular, the flow system is controlled by lateral heterogeneities, which characterize a surficial horizon made up of clayey colluviums and talus deposits, separated from the deeper saturated, fissured bedrock. One‐to‐ten relationships in hydraulic heads, monitored in piezometers crossing the fissured media, further uphold the crucial role played by the lateral contrasts of permeability in controlling the flow dynamics. On the whole, significant interactions with the nearby carbonate hydrostructure take place. Nevertheless, the heterogeneities of the siliciclastic succession and surficial horizon, coupled with the compartmentalization of the carbonate system, lead to a complex hydrogeological scenario. In a wider perspective, this study gives information of utmost importance in order to improve the implementation of mathematical models and configuration of tapping works within these heterogeneous and complex settings. Copyright © 2015 John Wiley & Sons, Ltd.     

Groundwater mixing between different aquifer types in a complex structural setting discerned by elemental and stable isotope geochemistry

Karst aquifers are well known for their intricate stratigraphy and geologic structures, which make groundwater characterization challenging because flowpaths and recharge sources are complex and difficult to evaluate. Geochemical data, collected from ten closely spaced production wells constructed in two karst aquifers (Bangor Limestone (Mb) and Tuscumbia Limestone/Fort Payne Chert (Mftp)) in Trussville, north‐central Alabama, illustrate two distinctive groundwater end‐members: (1) higher major ion, dissolved inorganic carbon, conductivity, alkalinity concentrations, heavier δ¹³C ratios (max: −10.2 ± 0.2‰ Vienna Pee Dee Belemnite (PDB)) and lower residence times (mean: 19.5 ± 2 years, n = 2) of groundwater in the Mb aquifer and (2) lower constituent concentrations, lighter δ¹³C ratios (min: −13.4 ± 0.2‰ PDB) and longer residence times of groundwater (mean: 23.6 ± 2 years, n = 4) in the Mftp aquifer. Summer and fall data and the binary mixing model show aquifer inter‐flow mixing along solution fractures and confirms the distinctive groundwater geochemistry of the two aquifers. Lowering of static water levels over the summer (drawdown from 2 to 5.2 m) leads to more reducing groundwater conditions (lower Eh values) and slightly enriched δ¹⁸O and δD ratios during the fall [δ¹⁸O: −4.8 ± 0.1 to −5.4 ± 0.1‰ Vienna Standard Mean Oceanic Water (VSMOW), n = 9; δD: −25.4 ± 1 to −27.4 ± 1‰ VSMOW, n = 9] when compared with summer season samples (δ¹⁸O: −5.1 ± 0.1 to −5.7 ± 0.1‰ VSMOW, n = 11; δD: −25.0 ± 1 to −30.6 ± 1‰ VSMOW, n = 11). GIS analyses confirm the localized origin of recharge to the investigated aquifers. The combination of GIS, field parameters and geochemistry analyses can be successfully used to identify recharge sources, evaluate groundwater flow and transport pathways and to improve understanding of how groundwater withdrawals impact the sustainability and susceptibility to contamination of karst aquifers. Copyright © 2015 John Wiley & Sons, Ltd.     

Site characterization model using least‐square support vector machine and relevance vector machine based on corrected SPT data (Nc)

Statistical learning algorithms provide a viable framework for geotechnical engineering modeling. This paper describes two statistical learning algorithms applied for site characterization modeling based on standard penetration test (SPT) data. More than 2700 field SPT values (N) have been collected from 766 boreholes spread over an area of 220 sqkm area in Bangalore. To get N corrected value (N_c), N values have been corrected (N_c) for different parameters such as overburden stress, size of borehole, type of sampler, length of connecting rod, etc. In three‐dimensional site characterization model, the function N_c=N_c (X, Y, Z), where X, Y and Z are the coordinates of a point corresponding to N_c value, is to be approximated in which N_c value at any half‐space point in Bangalore can be determined. The first algorithm uses least‐square support vector machine (LSSVM), which is related to a ridge regression type of support vector machine. The second algorithm uses relevance vector machine (RVM), which combines the strengths of kernel‐based methods and Bayesian theory to establish the relationships between a set of input vectors and a desired output. The paper also presents the comparative study between the developed LSSVM and RVM model for site characterization. Copyright © 2009 John Wiley & Sons, Ltd.     

一种采用基础矩阵描述的核线影像生成方法

设计了一种采用基础矩阵的立体像对核线影像制作方法,由基础矩阵来描述影像间的相对方位关系,避免了相对定向过程中复杂的迭代过程。通过稀疏匹配获取高精度的同名点,由基础矩阵确定每组同名点对应的同名核线,根据核线必相交于核点的原理,采用最小二乘法确定左右影像的核点坐标;然后根据核点坐标完成影像间核线的快速映射,沿核线方向采用双线性内插法进行核线重采样,完成核线影像制作。采用多组影像进行核线影像制作实验,对生成的核线影像进行同名特征点提取与匹配,结果表明采用该方法制作的核线影像上下视差较小,验证了方法的有效性。     

PEIV模型参数估计新算法

PEIV(Partial Errors-In-Variables)模型是EIV模型的扩展,它能解决系数矩阵含有非随机元素或存在结构特性的问题。针对常规PEIV模型算法的复杂性,提出了一种PEIV模型参数估计的新算法。该算法将系数矩阵含误差的元素看成是一类观测值,与平差模型原观测值构成两类观测值,将PEIV平差模型表示为类似于传统的最小二乘间接平差模型,再通过非线性最小二乘平差理论,推导出了算法的迭代公式和精度评定公式。算法迭代格式与间接平差类似,通过算例验证了算法的可行性和正确性。     

求解GM(1,1)模型新总体最小二乘算法

在GM(1,1)模型中系数矩阵和观测向量都是由原始序列组成的。系数矩阵中同样是有误差的,与观测向量中的误差一样,亦来源于原始序列,即它们误差同源。不同位置的相同元素应该有相同的改正数,采用传统总体最小二乘求解则不能达到此目的。针对这一缺陷,推导了一种新的总体最小二乘算法;并且通过算例验证了新方法的可行性和有效性。