Analytical investigation of the exact groundwater divide between rivers beyond the Dupuit–Forchheimer approximation

The groundwater divide is a key feature of river basins and significantly influenced by subsurface hydrological processes. For an unconfined aquifer between two parallel rivers or ditches, it has long been defined as the top of the water table based on the Dupuit–Forchheimer approximation. However, the exact groundwater divide is subject to the interface between two local flow systems transporting groundwater to rivers from the infiltration recharge. This study contributes a new analytical model for two-dimensional groundwater flow between rivers of different water levels. The flownet is delineated in the model to identify groundwater flow systems and the exact groundwater divide. Formulas with two dimensionless parameters are derived to determine the distributed hydraulic head, the top of the water table and the groundwater divide. The locations of the groundwater divide and the top of the water table are not the same. The distance between them in horizontal can reach up to 8.9% of the distance between rivers. Numerical verifications indicate that simplifications in the analytical model do not significantly cause misestimates in the location of the groundwater divide. In contrast, the Dupuit–Forchheimer approximation yields an incorrect water table shape. The new analytical model is applied to investigate groundwater divides in the Loess Plateau, China, with a Monte Carlo simulation process taking into account the uncertainties in the parameters.     

A biodiversity evaluation framework for restoration of aquatic macrophyte communities in shallow lakes driven by hydrological process management

Macrophyte community diversity and composition respond to ecosystem conservation and local environmental factors. In this study, we developed a multidimensional diversity framework for macrophyte communities, including the taxonomic and functional alpha and beta diversity. We used the framework to explore the relationships among water level regimes and these diversity parameters in a case study of China's Baiyangdian Lake. Analysis of indicators of hydrologic alteration divided the water level from 1959 to 2019 into four regimes (dry, <6.42 m; low, 6.42–7.23 m; medium, 7.23–8.19 m; high, >8.19 m). Alpha and beta diversity were significantly higher in the medium regime than in the low and high regimes. Redundancy analysis indicated that the maximum water depth significantly affected taxonomic alpha diversity, and total nitrogen (TN) and chemical oxygen demand (COD) concentration significantly affected functional alpha diversity, respectively. Mantel tests showed that TN, Secchi depth (SD), and water depth in the high water level regime significantly increased the total beta diversity and turnover components. TN was the main factor that increased total taxonomic beta diversity. Water level regime mainly influenced interspecific relationships by changing the TN and COD concentration. The water level should be maintained between the medium and high water level regimes to promote restoration of the macrophyte community and improve ecosystem stability. The biodiversity evaluation framework would provide a deeper insight into the hydrological process management for restoration of aquatic macrophyte communities in shallow lakes.     

Rainfall intensity affects runoff responses in a semi-arid catchment

The source and hydrochemical makeup of a stream reflects the connectivity between rainfall, groundwater, the stream, and is reflected to water quantity and quality of the catchment. However, in a semi-arid, thick, loess covered catchment, temporal variation of stream source and event associated behaviours are lesser known. Thus, the isotopic and chemical hydrographs in a widely distributed, deep loess, semi-arid catchment of the northern Chinese Loess Plateau were characterized to determine the source and hydrochemical behaviours of the stream during intra-rainfall events. Rainfall and streamflow were sampled during six hydrologic events coupled with measurements of stream baseflow and groundwater. The deuterium isotope (²H), major ions (Cl⁻, SO₄²⁻, NO₃⁻, Ca²⁺, K⁺, Mg²⁺, and Na⁺) were evaluated in water samples obtained during rainfall events. Temporal variation of ²H and Cl⁻ measured in the groundwater and stream baseflow prior to rainfall was similar; however, the isotope compositions of the streamflow fluctuated significantly and responded quickly to rainfall events, likely due to an infiltration excess, overland dominated surface runoff during torrential rainfall events. Time source separation using ²H demonstrated greater than 72% on average, the stream composition was event water during torrential rainfall events, with the proportion increasing with rainfall intensity. Solutes concentrations in the stream had loglinear relationships with stream discharge, with an outling anomaly with an example of an intra-rainfall event on Oct. 24, 2015. Stream Cl⁻ behaved nonconservative during rainfall events, temporal variation of Cl⁻ indicated a flush and washout at the onset of small rainfall events, a dilution but still high concentration pattern in high discharge and old water dominated in regression flow period. This study indicates rainfall intensity affects runoff responses in a semi-arid catchment, and the stored water in the thick, loess covered areas was less connected with stream runoff. Solute transport may threaten water quality in the area, requiring further analysis of the performance of the eco-restoration project.     

Modelling the effect of changing precipitation inputs on deep soil water utilization

Forests in the Southeastern United States are predicted to experience future changes in seasonal patterns of precipitation inputs as well as more variable precipitation events. These climate change‐induced alterations could increase drought and lower soil water availability. Drought could alter rooting patterns and increase the importance of deep roots that access subsurface water resources. To address plant response to drought in both deep rooting and soil water utilization as well as soil drainage, we utilize a throughfall reduction experiment in a loblolly pine plantation of the Southeastern United States to calibrate and validate a hydrological model. The model was accurately calibrated against field measured soil moisture data under ambient rainfall and validated using 30% throughfall reduction data. Using this model, we then tested these scenarios: (a) evenly reduced precipitation; (b) less precipitation in summer, more in winter; (c) same total amount of precipitation with less frequent but heavier storms; and (d) shallower rooting depth under the above 3 scenarios. When less precipitation was received, drainage decreased proportionally much faster than evapotranspiration implying plants will acquire water first to the detriment of drainage. When precipitation was reduced by more than 30%, plants relied on stored soil water to satisfy evapotranspiration suggesting 30% may be a threshold that if sustained over the long term would deplete plant available soil water. Under the third scenario, evapotranspiration and drainage decreased, whereas surface run‐off increased. Changes in root biomass measured before and 4 years after the throughfall reduction experiment were not detected among treatments. Model simulations, however, indicated gains in evapotranspiration with deeper roots under evenly reduced precipitation and seasonal precipitation redistribution scenarios but not when precipitation frequency was adjusted. Deep soil and deep rooting can provide an important buffer capacity when precipitation alone cannot satisfy the evapotranspirational demand of forests. How this buffering capacity will persist in the face of changing precipitation inputs, however, will depend less on seasonal redistribution than on the magnitude of reductions and changes in rainfall frequency.     

An Early Cretaceous slab window beneath central Tibet,SW China: Evidence from OIB‐like alkaline gabbros in the Duolong area

The geodynamic mechanism of the late Early Cretaceous magmatic flare‐up in the collisional zone between the Lhasa and Qiangtang terranes in Tibet is controversial because of a scarcity of robust evidence. To address this problem, we report geochronological, geochemical and Hf isotopic data for the newly discovered Gufeng gabbros from the Duolong Cu–Au mineral district of the western Bangong–Nujiang Suture Zone (BNSZ). The gabbro samples, dated at 126.3 ± 1.8 Ma, show geochemical similarities to typical ocean island basalt (OIB) and have positive εHf(t) values of +3.3 to +6.9. The gabbros were generated by decompression melting of deep upwelling asthenosphere. This event is best explained by slab break‐off and the resultant development of a slab window beneath central Tibet.     

Variations in glacier volume and snow cover and their impact on lake storage in the Paiku Co Basin,in the Central Himalayas

Glaciers and snow cover are important constituents of the surface of the Tibetan Plateau. The responses of these phenomena to global environmental changes are sensitive, rapid and intensive due to the high altitudes and arid cold climate of the Tibetan Plateau. Based on multisource remote sensing data, including Landsat images, MOD10A2 snow product, ICESat, Cryosat-2 altimetry data and long-term ground climate observations, we analysed the dynamic changes of glaciers, snow melting and lake in the Paiku Co basin using extraction methods for glaciers and lake, the degree-day model and the ice and lake volume method. The interaction among the climate, ice-snow and the hydrological elements in Paiku Co is revealed. From 2000 to 2018, the basin tended to be drier, and rainfall decreased at a rate of −3.07 mm/a. The seasonal temperature difference in the basin increased, the maximum temperature increased at a rate of 0.02°C/a and the minimum temperature decreased at a rate of −0.06°C/a, which accelerated the melting from glaciers and snow at rates of 0.55 × 10⁷ m³/a and 0.29 × 10⁷ m³/a, respectively. The rate of contribution to the lake from rainfall, snow and glacier melted water was 55.6, 27.7 and 16.7%, respectively. In the past 18 years, the warmer and drier climate has caused the lake to shrink. The water level of the lake continued to decline at a rate of −0.02 m/a, and the lake water volume decreased by 4.85 × 10⁸ m³ at a rate of −0.27 × 10⁸ m³/a from 2000 to 2018. This evaluation is important for understanding how the snow and ice melting in the central Himalayas affect the regional water cycle.     

A novel coupled FDM-DEM modelling method for flexible membrane boundary in laboratory tests

Difficulties are involved in discrete element method (DEM) modelling of the flexible boundary, that is, the membranes covering the soil sample, which can be commonly found in contemporary laboratory soil tests. In this paper, a novel method is proposed wherein the finite difference method (FDM) and DEM are coupled to simulate the rubber membrane and soil body, respectively. Numerical plane strain and triaxial tests, served by the flexible membrane, are implemented and analysed later. The effect of the membrane modulus on the measurement accuracy is considered, with analytical formulae derived to judge the significance of this effect. Based on an analysis of stress-strain responses and the grain rotation field, the mechanical performances produced by the flexible and rigid lateral boundaries are compared for the plane strain test. The results show that (1) the effect of the membrane on the test result becomes more significant at larger strain level because the membrane applies additional lateral confining pressure to the soil body; (2) the tested models reproduce typical stress and volumetric paths for specimens with shear bands; (3) for the plane strain test, the rigid lateral boundary derives a much higher peak strength and larger bulk dilatation, but a similar residual strength, compared with the flexible boundary. The latter produces a more uniform (or ‘diffuse') rotation field and more mobilised local kinematics than does the former. All simulations show that the proposed FDM-DEM coupling method is able to simulate laboratory tests with a flexible boundary membrane.     

Directional interpolation infinite element for dynamic problems in saturated porous media

In this study, a directional interpolation infinite element suited to a saturated porous medium is presented to account for dynamic problems with semi-infi     

Three-dimensional finite element simulation of fracture propagation in rock specimens with pre-existing fissure(s) under compression and their strength analysis

A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials.     

Hydrological management strategies for the control of algal blooms in regulated lowland rivers

Reservoirs of lowland floodplain rivers with eutrophic backgrounds cause variations in the hydrological and hydraulic conditions of estuaries and low-dam reservoir areas, which can promote planktonic algae to proliferate and algal bloom outbreaks. Understanding the ecological effects of variations in hydrological and hydraulic processes in lowland rivers is important for algal bloom control. In this study, the middle and lower reaches of the Han River, China, a typical regulated lowland river with a eutrophic background, are selected. Based on the effect of hydrological and hydraulic variability on algal blooms, a hydrological management strategy for river algal bloom control is proposed. The results showed that (a) differences in river morphology and background nutrient levels cause significant differences in the critical threshold flow velocities for algal bloom outbreaks between natural river and low-dam reservoir sections; there is no uniform threshold flow velocity for algal bloom control. (b) There are significant differences in the river hydrological/hydraulic conditions between years with and without algal blooms. The average river flow, water level and velocity in years with algal blooms are significantly lower than those in years without algal blooms. (c) For different river sections where algal blooms occur and to meet the threshold flow velocities, the joint operation of cascade reservoirs and diversion projects is an effective method to prevent and control algal blooms in regulated lowland rivers. This study is expected to deepen our understanding of the ecological significance of special hydrological processes and guide algal bloom management in regulated lowland rivers.