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.     

Water uptake of riparian plants in the lower Lhasa River Basin,South Tibetan Plateau using stable water isotopes

Riparian plants can adapt their water uptake strategies based on climatic and hydrological conditions within a river basin. The response of cold-alpine riparian trees to changes in water availability is poorly understood. The Lhasa River is a representative cold-alpine river in South Tibet and an under-studied environment. Therefore, a 96 km section of the lower Lhasa River was selected for a study on the water-use patterns of riparian plants. Plant water, soil water, groundwater and river water were measured at three sites for δ¹⁸O and δ²H values during the warm-wet and cold-dry periods in 2018. Soil profiles differed in isotope values between seasons and with the distance along the river. During the cold-dry period, the upper parts of the soil profiles were significantly affected by evaporation. During the warm-wet period, the soil profile was influenced by precipitation infiltration in the upper reaches of the study area and by various water sources in the lower reaches. Calculations using the IsoSource model indicated that the mature salix and birch trees (Salix cheilophila Schneid. and Betula platyphylla Suk.) accessed water from multiple sources during the cold-dry period, whereas they sourced more than 70% of their requirement from the upper 60–80 cm of the soil profile during the warm-wet period. The model indicated that the immature rose willow tree (Tamarix ramosissima Ledeb) accessed 66% of its water from the surface soil during the cold-dry period, but used the deeper layers during the warm-wet period. The plant type was not the dominant factor driving water uptake patterns in mature plants. Our findings can contribute to strategies for the sustainable development of cold-alpine riparian ecosystems. It is recommended that reducing plantation density and collocating plants with different rooting depths would be conducive to optimal plant growth in this environment.     

Discharge‐sediment processes of the Zhadang glacier on the Tibetan Plateau measured with a high frequency data acquisition system

In high elevation cold regions of the Tibetan Plateau, suspended sediment transfer from glacier meltwater erosion is one of the important hydrological components. The Zhadang glacier is a typical valley‐type glacier in the Nyainqentanglha Mountains on the Tibetan Plateau. To make frequent and long period records of meltwater runoff and sediment processes in the very high elevation and isolated regions, an automatic system was installed near the glacier snout (5400 m a.s.l) in August 2013, to measure the transient discharge and sediment processes at 5‐min interval, which is shorter than the time span for the water flow to traverse the catchment from the farthest end to the watershed outlet. Diurnal variations of discharge, and suspended sediment concentration (SSC) were recorded at high frequency for the Zhadang glacier, before suspended sediment load (SSL) was computed. Hourly SSC varied from the range of 0.2 kg/m³ to 0.5 kg/m³ (at 8:00–9:00) to the range of 2.0 kg/m³ to 4.0 kg/m³ (at 17:00–18:00). The daily SSL was 32.24 t during the intense ablation period. Hourly SSC was linearly correlated with discharge (r = 0.885**, n = 18, p < 0.01). A digit‐eight hysteresis loop was observed for the sediment transport in the glacier area. Air temperature fluctuations influence discharge, and then result in the sediment variations. The results of this study provide insight into the responses of suspended sediment delivery processes with a high frequency data in the high elevation cold regions. Copyright © 2016 John Wiley & Sons, Ltd.     

Towards more systematic perceptual model development: a case study using 3 Luxembourgish catchments

The synthesis of experimental understanding of catchment behaviour and its translation into qualitative perceptual models is an important objective of hydrological sciences. We explore this challenge by examining the cumulative understanding of the hydrology of three experimental catchments and how it evolves through the application of different investigation techniques. The case study considers the Huewelerbach, Weierbach and Wollefsbach headwater catchments of the Attert basin in Luxembourg. Subsurface investigations including bore holes and pits, analysis of soil samples and Electrical Resistivity Tomography measurements are presented and discussed. Streamflow and tracer data are used to gain further insights into the streamflow dynamics of the catchments, using end‐member mixing analysis and hydrograph separation based on dissolved silica and electrical conductivity. We show that the streamflow generating processes in all three catchments are controlled primarily by the subsolum and underlying bedrock. In the Huewelerbach, the permeable sandstone formation supports a stable groundwater component with little seasonality, which reaches the stream through a series of sources at the contact zone with the impermeable marls formation. In the Weierbach, the schist formation is relatively impermeable and supports a ‘fill and spill’‐type of flow mechanism; during wet conditions, it produces a delayed response dominated by pre‐event water. In the Wollefsbach, the impermeable marls formation is responsible for a saturation‐excess runoff generating process, producing a fast and highly seasonal response dominated by event water. The distinct streamflow generating processes of the three catchments are represented qualitatively using perceptual models. The perceptual models are in turn translated into quantitative conceptual models, which simulate the hydrological processes using networks of connected reservoirs and transfer functions. More generally, the paper illustrates the evolution of perceptual models based on experimental fieldwork data, the translation of perceptual models into conceptual models and the value of different types of data for processes understanding and model representation. Copyright © 2014 John Wiley & Sons, Ltd.     

Tensionless–frictionless interaction of flexible annular foundation with a transversely isotropic multi‐layered half‐space

A transversely isotropic multi‐layered half‐space, with axis of material symmetry perpendicular to the free surface, supports a flexible either annular or solid circle foundation. The contact area of the foundation and the half‐space is considered to be both frictionless and tensionless. The foundation is assumed to be affected by a vertical static axisymmetric load. Detailed analysis of the interaction of these two systems with different thickness of layers is the target of this paper. With the use of ring load Green's functions for both the foundation and the continuum half‐space, an integral equation accompanied with some inequalities is introduced to model the complex BVP. With the incorporation of ring‐shape FEM, we are capable of capturing both regular and singular solution smoothly. The validity of the combination of the analytical and numerical method is proved with comparing the results of this paper with a number of benchmark cases of both linear and nonlinear interaction of circular and annular foundation with half‐space. Some new illustrations are presented to portray the aspect of the anisotropy and layering of the half‐space. Copyright © 2014 John Wiley & Sons, Ltd.     

Stable staurolite–cordierite assemblages in K‐poor metapelitic schists in Aston and Hospitalet gneiss domes of the central Pyrenees (France,Andorra)

Staurolite–cordierite assemblages are common in mica schists of the Aston and Hospitalet gneiss domes of the central Axial Zone, Pyrenees (France, Andorra). Within a 200 m wide zone, staurolite, cordierite and andalusite porphyroblasts contain inclusion trails that preserve the same stage of development of a crenulation cleavage, strongly suggesting that all three phases are contemporaneous. Their syntectonic growth occurred during a short period at the beginning of the formation of the dominant schistosity (S₂) of the domes. Staurolite and cordierite touching each other further indicates an equilibrium relationship. Whole‐rock analyses show that some staurolite–cordierite schists are depleted in K₂O compared to post‐Archean shales (PAAS) and amphibolite facies pelites. Analysis of the st‐crd paragenesis in K‐poor schists without muscovite using KFMASH and MnNCKFMASH petrogentic grids, pseudosections and AFM compatibility diagrams predicts stable conditions at pressures of ~3.5 kbar at 575 °C. For metapelites with intermediate X_Mg values (0.7 >  X_Mg >0.48) a ‘muscovite‐out window’ exists from 550–650 °C at 3.5 kbar in the KFMASH system. Conventional thermobarometry (GB‐GASP, AvT‐AvP) and petrogenetic grids show an isobaric P–T path to peak temperatures of ~650 °C, supported by the presence of sillimanite‐K‐feldspar gneiss and migmatites. LP‐HT metamorphism in the Aston dome is related to early Carboniferous (c. 339 Ma) granitic intrusions into the dome core. As metamorphism is directly linked with the formation of the main S₂ schistosity, the temporal relations demonstrated in this study conflict with previous studies which constrained LP‐HT metamorphism and the development of flat‐lying schistosity to the late Carboniferous (315–305 Ma) – at least in the eastern Axial Zone.     

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 nonstationary extreme value distribution for analysing the cessation of karst spring discharge

The effects of climate change and population growth in recent decades are leading us to consider their combined and potentially extreme consequences, particularly regarding hydrological processes, which can be modeled using a generalized extreme value (GEV) distribution. Most of the GEV models were based on a stationary assumption for hydrological processes, in contrast to the nonstationary reality due to climate change and human activities. In this paper, we present the nonstationary generalized extreme value (NSGEV) distribution and use it to investigate the risk of Niangziguan Springs discharge decreasing to zero. Rather than assuming the location, scale, and shape parameters to be constant as one might do for a stationary GEV distribution analysis, the NSGEV approach can reflect the dynamic processes by defining the GEV parameters as functions of time. Because most of the GEV model is designed to evaluate maxima (e.g. flooding, represented by positive numbers), and spring discharge cessation is a ?minima’, we deduced an NSGEV model for minima by applying opposite numbers, i.e. negative instead of positive numbers. The results of the model application to Niangziguan Springs showed that the probability of zero discharge at Niangziguan Springs will be 1/80 in 2025, and 1/10 in 2030. After 2025, the rate of decrease in spring discharge will accelerate, and the probability that Niangziguan Springs will cease flowing will dramatically increase. The NSGEV model is a robust method for analysing karst spring discharge. Copyright © 2013 John Wiley & Sons, Ltd.     

Piecewise prediction model for watershed‐scale erosion and sediment yield of individual rainfall events on the Loess Plateau,China

Establishing a universal watershed‐scale erosion and sediment yield prediction model represents a frontier field in erosion and soil/water conservation. The research presented here was conducted on the Chabagou watershed, which is located in the first sub‐region of the hill‐gully area of the Loess Plateau, China. A back‐propagation artificial neural model for watershed‐scale erosion and sediment yield was established, with the accuracy of the model, then compared with that of multiple linear regression. The sensitivity degree of various factors to erosion and sediment yield was quantitatively analysed using the default factor test. On the basis of the sensitive factors and the fractal information dimension, the piecewise prediction model for erosion and sediment yield of individual rainfall events was established and further verified. The results revealed the back‐propagation artificial neural network model to perform better than the multiple linear regression model in terms of predicting the erosion modulus, with the former able to effectively characterize dynamic changes in sediment yield under comprehensive factor conditions. The sensitivity of runoff erosion power and runoff depth to the erosion and sediment yield associated with individual rainfall events was found to be related to the complexity of surface topography. The characteristics of such a hydrological response are thus closely related to topography. When the fractal information dimension is greater than the topographic threshold, the accuracy of prediction using runoff erosion power is higher than that of using runoff depth. In contrast, when the fractal information dimension is smaller than the topographic threshold, the accuracy of prediction using runoff depth is higher than that of using runoff erosion power. The developed piecewise prediction model for watershed‐scale erosion and sediment yield of individual rainfall events, which introduces runoff erosion power and runoff depth using the fractal information dimension as a boundary, can be considered feasible and reliable and has a high prediction accuracy. Copyright © 2013 John Wiley & Sons, Ltd.     

Water level observations from unmanned aerial vehicles for improving estimates of surface water–groundwater interaction

Integrated hydrological models are usually calibrated against observations of river discharge and piezometric head in groundwater aquifers. Calibration of such models against spatially distributed observations of river water level can potentially improve their reliability and predictive skill. However, traditional river gauging stations are normally spaced too far apart to capture spatial patterns in the water surface, whereas spaceborne observations have limited spatial and temporal resolution. Unmanned aerial vehicles can retrieve river water level measurements, providing (a) high spatial resolution; (b) spatially continuous profiles along or across the water body, and (c) flexible timing of sampling. A semisynthetic study was conducted to analyse the value of the new unmanned aerial vehicle‐borne datatype for improving hydrological models, in particular estimates of groundwater–surface water (GW–SW) interaction. Mølleåen River (Denmark) and its catchment were simulated using an integrated hydrological model (MIKE 11–MIKE SHE). Calibration against distributed surface water levels using the Differential Evolution Adaptive Metropolis algorithm demonstrated a significant improvement in estimating spatial patterns and time series of GW–SW interaction. After water level calibration, the sharpness of the estimates of GW–SW time series improves by ~50% and root mean square error decreases by ~75% compared with those of a model calibrated against discharge only.