On the dependence of the saturated hydraulic conductivity upon the effective porosity through a power law model at different scales

We study the scale dependence of the saturated hydraulic conductivity K_s through the effective porosity n_e by means of a newly developed power‐law model (PLM) which allows to use simultaneously measurements at different scales. The model is expressed as product between a single PLM (capturing the impact of the dominating scale) and a characteristic function κ^? accounting for the correction because of the other scale(s). The simple (closed form) expression of the κ^?‐function enables one to easily identify the scales which are relevant for K_s. The proposed model is then applied to a set of real data taken at the experimental site of Montalto Uffugo (Italy), and we show that in this case two (i.e. laboratory and field) scales appear to be the main ones. The implications toward an important application (solute transport) in Hydrology are finally discussed. Copyright © 2016 John Wiley & Sons, Ltd.     

Investigation of the elevation of saltwater wedge due to subsurface dams

Subsurface dams are rather effective and used for the prevention of saltwater intrusion in coastal regions around the world. We carried out the laboratory experiments to investigate the elevation of saltwater wedge after the construction of subsurface dams. The elevation of saltwater wedge refers to the upward movement of the downstream saltwater wedge because the subsurface dams obstruct the regional groundwater flow and reduce the freshwater discharge. Consequently, the saltwater wedge cannot further extend in the longitudinal direction but rises in the vertical profile resulting in significant downstream aquifer salinization. In order to quantitatively address this issue, field-scale numerical simulations were conducted to explore the influence of various dam heights, distances, and hydraulic gradients on the elevation of saltwater wedge. Our investigation shows that the upward movement of the saltwater wedge and its areal extension in the vertical domain of the downstream aquifer become more severe with a higher dam and performed a great dependence on the freshwater discharge. Furthermore, the increase of the hydraulic gradient and the dam distance from the sea boundary leads to a more pronounced wedge elevation. This phenomenon comes from the variation of the freshwater discharge due to the modification of dam height, location, and hydraulic gradient. Large freshwater discharge can generate greater repulsive force to restrain the elevation of saltwater wedge. These conclusions provide theoretical references for the behaviour of the freshwater–seawater interface after the construction of subsurface dams and help optimize the design strategy to better utilize the coastal groundwater resources.     

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.     

Lateral extrusion along the Altyn Tagh Fault,Qilian Shan (NE Tibet): insight from a 3D crustal budget

The lithospheric strike‐slip Altyn Tagh Fault has accommodated hundreds of kilometres of displacement between the Qaidam and Tarim blocks since its Eocene reactivation. However, the way the deformation is accommodated in the Qilian Shan and further east remains uncertain. Based on 360 km of north‐eastward migration of the relatively rigid Qaidam block along the Altyn Tagh Fault and 3D isovolumetric balancing of the crustal deformation within the Altyn Tagh Fault–Qilian Shan system, we demonstrate that 250 ± 28 km (43.8–49.4%) of N20E directed crustal shortening and an additional ~250–370 km of eastward motion of the Qilian Shan crust must be accounted for by strike‐slip faulting in the Qilian Shan and crustal thickening in the Qinling area, as well as by extension in the adjoining North China block graben systems.     

A discontinuous finite element suspended sediment transport model for water quality assessments in river networks

Suspended sediment plays an important role in the distribution and transport of many pollutants (such as radionuclides) in rivers. Pollutants may adsorb on fine suspended particles (e.g. clay) and spread according to the suspended sediment movement. Hence, the simulation of the suspended sediment mechanism is indispensable for realistic transport modelling. This paper presents and tests a simple mathematical model for predicting the suspended sediment transport in river networks. The model is based on the van Rijn suspended load formula and the advection–diffusion equation with a source or sink term that represents the erosion or deposition fluxes. The transport equation is solved numerically with the discontinuous finite element method. The model evaluation was performed in two steps, first by comparing model simulations with the measured suspended sediment concentrations in the Grote Nete–Molse Nete River in Belgium, and second by a model intercomparison with the sediment transport model NST MIKE 11. The simulations reflect the measurements with a Nash‐Sutcliffe model efficiency of 0.6, while the efficiency between the proposed model and the NST MIKE 11 simulations is 0.96. Both evaluations indicate that the proposed sediment transport model, that is sufficiently simple to be practical, is providing realistic results.     

Simulation of site-scale water fluxes in desert and natural oasis ecosystems of the arid region in Northwest China

The study of water fluxes is important to better understand hydrological cycles in arid regions. Data-driven machine learning models have been recently applied to water flux simulation. Previous studies have built site-scale simulation models of water fluxes for individual sites separately, requiring a large amount of data from each site and significant computation time. For arid areas, there is no consensus as to the optimal model and variable selection method to simulate water fluxes. Using data from seven flux observation sites in the arid region of Northwest China, this study compared the performance of random forest (RF), support vector machine (SVM), back propagation neural network (BPNN), and multiple linear regression (MLR) models in simulating water fluxes. Additionally, the study investigated inter-annual and seasonal variation in water fluxes and the dominant drivers of this variation at different sites. A universal simulation model for water flux was built using the RF approach and key variables as determined by MLR, incorporating data from all sites. Model performance of the SVM algorithm (R² = 0.25–0.90) was slightly worse than that of the RF algorithm (R² = 0.41–0.91); the BPNN algorithm performed poorly in most cases (R² = 0.15–0.88). Similarly, the MLR results were limited and unreliable (R² = 0.00–0.66). Using the universal RF model, annual water fluxes were found to be much higher than the precipitation received at each site, and natural oases showed higher fluxes than desert ecosystems. Water fluxes were highest during the growing season (May–September) and lowest during the non-growing season (October–April). Furthermore, the dominant drivers of water flux variation were various among different sites, but the normalized difference vegetation index (NDVI), soil moisture and soil temperature were important at most sites. This study provides useful insights for simulating water fluxes in desert and oasis ecosystems, understanding patterns of variation and the underlying mechanisms. Besides, these results can make a contribution as the decision-making basis to the water management in desert and oasis ecosystems.     

Lake sedimentary processes and vegetation changes over the last 45k cal a bp in the uplands of south-eastern Amazonia

This study addresses paleoclimate influences in a southern Amazonia ecotone based on multiproxy records from lakes of the Carajás region during the last 45k cal a bp. Wet and cool environmental conditions marked the initial deposition in shallow depressions with detrital sediments and high weathering rates until 40k cal a bp. Concomitantly, forest and C3 canga plants, along with cool-adapted taxa, developed; however, short drier episodes enabled expansion of C4 plants and diagenetic formation of siderite. A massive event of siderite formation occurred approximately 30k cal a bp under strong drier conditions. Afterwards, wet and cool environmental conditions returned and persisted until the Last Glacial Maximum (LGM). The LGM was marked by lake-level lowstands and subaerial exposure. The transition from the LGM to the Holocene is marked by the onset of oscillations in temperature and humidity, with an expansion of forest and canga plants. Cool taxa were present for the last time in the Carajás region ~ 9.5–9k cal a bp. After 10k cal a bp , shallow lakes became upland swamps due to natural infilling processes, but the current vegetation types and structures of the plateaus were acquired only after 3k cal a bp under wetter climatic conditions.