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.     

Beaver dam analogues drive heterogeneous groundwater–surface water interactions

Beaver dam analogues (BDAs) are a cost-effective stream restoration approach that leverages the recognized environmental benefits of natural beaver dams on channel stability and local hydrology. Although natural beaver dams are known to exert considerable influence on the hydrologic conditions of a stream system by mediating geomorphic processes, nutrient cycling, and groundwater–surface water interactions, the impacts of beaver-derived restoration methods on groundwater–surface water exchange are poorly characterized. To address this deficit, we monitored hyporheic exchange fluxes and streambed porewater biogeochemistry across a sequence of BDAs installed along a central Wyoming stream during the summer of 2019. Streambed fluxes were quantified by heat tracing methods and vertical hydraulic gradients. Biogeochemical activity was evaluated using major ion porewater chemistry and principal component analysis. Vertical fluxes of approximately 1.0 m/day were observed around the BDAs, as was the development of spatially heterogeneous zones of nitrate production, groundwater upwelling, and anaerobic reduction. Strong contrasts in hyporheic zone processes were observed across BDAs of differing sizes. This suggests that structures may function with size-dependent behaviour, only altering groundwater–surface water interactions after a threshold hydraulic step height is exceeded. Patterns of hyporheic exchange and biogeochemical cycling around the studied BDAs resemble those around natural beaver dams, suggesting that BDAs may provide comparable benefits to channel complexity and near-stream function over a 1-year period.     

张掖盆地地下水硝酸盐污染与人体健康风险评价

地下水是张掖盆地的重要水资源，其硝酸盐污染尚未得到足够重视。对张掖盆地2004、2015年地下水硝酸盐浓度进行了系统分析，并采用美国环境保护署（USEPA）推荐的健康风险评价模型评估了地下水硝酸盐的健康风险。结果表明：自2004年以来张掖盆地地下水硝酸盐污染日趋严重。2015年硝酸盐浓度最高已达到283.32 mg·L^-1，17.61%的采样点硝酸盐氮浓度超过GB5749-2006《生活饮用水卫生标准》中饮用地下水限量值（20 mg·L^-1）。研究区人群经皮肤接触途径摄入硝酸盐的健康风险在可接受水平，而饮水摄入硝酸盐的健康风险较高，总风险中饮水途径引起健康风险的贡献率占99.40%，远大于皮肤接触途径。儿童经饮水摄入和皮肤接触两种途径的健康风险均显著高于成人，分别为成人的1.544倍和1.039倍。32.39%的采样点地下水硝酸盐对儿童的健康风险超出了可接受水平，14.79%的采样点地下水硝酸盐对成人的健康风险不可接受。甘州区城区、临泽县北部边缘及高台县城区周围硝酸盐浓度最高，这些区域内所有人群都面临硝酸盐引发的高健康风险，其余区域硝酸盐引发的健康风险相对较低。     

Time–space fractional governing equations of transient groundwater flow in confined aquifers: Numerical investigation

In order to model non‐Fickian transport behaviour in groundwater aquifers, various forms of the time–space fractional advection–dispersion equation have been developed and used by several researchers in the last decade. The solute transport in groundwater aquifers in fractional time–space takes place by means of an underlying groundwater flow field. However, the governing equations for such groundwater flow in fractional time–space are yet to be developed in a comprehensive framework. In this study, a finite difference numerical scheme based on Caputo fractional derivative is proposed to investigate the properties of a newly developed time–space fractional governing equations of transient groundwater flow in confined aquifers in terms of the time–space fractional mass conservation equation and the time–space fractional water flux equation. Here, we apply these time–space fractional governing equations numerically to transient groundwater flow in a confined aquifer for different boundary conditions to explore their behaviour in modelling groundwater flow in fractional time–space. The numerical results demonstrate that the proposed time–space fractional governing equation for groundwater flow in confined aquifers may provide a new perspective on modelling groundwater flow and on interpreting the dynamics of groundwater level fluctuations. Additionally, the numerical results may imply that the newly derived fractional groundwater governing equation may help explain the observed heavy‐tailed solute transport behaviour in groundwater flow by incorporating nonlocal or long‐range dependence of the underlying groundwater flow field.     

Evaluation of the temporal variations of groundwater storage and its interactions with climatic variables using GRACE data and hydrological models: A study from Turkey

Monitoring of the fluctuations of groundwater storage is particularly important in arid and semi-arid regions where water scarcity brings about various challenges. Remote sensing data and techniques play a preponderant role in developing solutions to environmental problems. The launch of Gravity Recovery and Climate Experiment (GRACE) satellites has eased the remote monitoring and evaluation of groundwater resources with an unprecedented precision over large scales. Within the scope of the current study, the latest release (RL06) of GRACE mass concentrations (Mascons) from Jet Propulsion Laboratory (JPL) dataset as well as Global Land Data Assimilation System (GLDAS) models of Noah and Catchment Land Surface Model (CLSM) were used to provide Groundwater Storage Anomalies (GWSA) over Turkey. The temporal interactions of the estimated GWSA with the climatic variables of precipitation and temperature (derived from the reanalysis datasets of CHELSA [Climatologies at High resolution for the Earth's Land Surface Areas] and FLDAS [the Famine Early Warning Systems Network Land Data Assimilation System], respectively) were investigated statistically. The results suggest that there is a descending trend (from 2003 to 2016) for Terrestrial Water Storage Anomalies (TWSA) and GWSA over Turkey with a total loss of 11 and 6 cm of water, respectively. The statistical analysis results also indicate that the monthly variations of GWSA over Turkey are highly correlated with precipitation and temperature at 2-month lag. The analysis of the climatology (long-term) values of monthly GWSA, precipitation and temperature also revealed high agreement between the variables.     

Properties of geogrid-reinforced marine slope due to the groundwater level changes

Evaluation of slope stability, especially in the absence of a proper bed such as marine soils, is one of the most important issues in geotechnical engineering. Using geogrid layers to enhance the strength and stability of embankments is regarded as a commendable stabilization method. On the other hand, groundwater level erratically fluctuates in coastal areas. Therefore, the aim of this research is to study the effects of groundwater level changes on stability of a geogrid-reinforced slope on loose marine soils in Qeshm Island, Iran. At first, geotechnical properties of the site were obtained by comprehensive series of geotechnical laboratory and in situ tests. Then, by simultaneous changes of groundwater level and several parameters such as embankment slope, loading, geogrid length, geogrid number, and tensile strength of geogrid, different characteristics such as embankment safety factor (SF), vertical and horizontal displacements at embankment top and embankment base were studied. It was observed that groundwater level had significant effects on behavior of the embankment. For most of the observations, by decreasing the groundwater level, the displacements decreased and consequently safety factor increased. Increasing the length, number, and tensile strength of geogrid led to the reduction of displacements and an increase in the safety factor.     

Zircon U–Pb age, trace element and Hf isotope composition of Kongling terrane in the Yangtze Craton: refining the timing of Palaeoproterozoic high-grade metamorphism

U–Pb age, trace element and Hf isotope compositions of zircon were analysed for a metasedimentary rock and two amphibolites from the Kongling terrane in the northern part of the Yangtze Craton. The zircon shows distinct morphological and chemical characteristics. Most zircon in an amphibolite shows oscillatory zoning, high Th/U and ¹⁷⁶Lu/¹⁷⁷Hf ratios, high formation temperature, high trace element contents, clear negative Eu anomaly, as well as HREE-enriched patterns, suggesting that it is igneous. The zircon yields a weighted mean ²⁰⁷Pb/²⁰⁶Pb age of 2857 ± 8 Ma, representing the age of the magmatic protolith. The zircon in the other two samples is metamorphic. It has low Th/U ratios, low trace element concentrations, variable HREE contents (33.8 ≥ Lu_N≥2213; 14.7 ≤ Lu_N/Sm_N ≤ 354) and ¹⁷⁶Lu/¹⁷⁷Hf ratios (0.000030–0.001168). The data indicate that the zircon formed in the presence of garnet and under upper amphibolite facies conditions. The metamorphic zircon yields a weighted mean ²⁰⁷Pb/²⁰⁶Pb age of 2010 ± 13 Ma. These results combined with previously obtained Palaeoproterozoic metamorphic ages suggest a c. 2.0 Ga Palaeoproterozoic collisional event in the Yangtze Craton, which may result from the assembly of the supercontinent Columbia. The zircon in two samples yields weighted mean two-stage Hf model ( T _DM2) ages of 3217 ± 110 and 2943 ± 50 Ma, respectively, indicating that their protoliths were mainly derived from Archean crust.     

平凉隐伏岩溶水环境同位素研究

地下水资源可持续利用的一个急待解决的重要问题,是对地下水补给和更新能力的评价.利用环境同位素技术研究地下水的补给和可更新性是当前较为新颖的方法之一.在西北干旱、半干旱的隐伏岩溶地区,地下水埋藏条件复杂,常规的地质勘探方法所能提供的水文地质信息有限,环境同位素方法在研究地下水的补给及可更新能力方面发挥了优势,可对传统方法进行补充和验证.其结果表明,研究区隐伏岩溶水形成较早,且有大量现代水的混入,平均混入量为54%.说明区内隐伏岩溶水的补给和更新能力较好.环境同位素分析结果还显示,大岔河隐伏岩溶水为一相对独立、半开放的水文地质单元,其补给来源部分为流域内大气降水、地表水的补给,部分为东南部三道沟岩溶地下水的补给;根据环境同位素EPM模型计算,地下水的滞留时间为36 a.地下水储存量为1.314×108 m3; 储水系数为7.29×10-3.这一结果与传统勘探方法的计算结果基本吻合,说明环境同位素方法的实用性.     

西藏结则茶卡湖岸沉积物中铀系年龄及意义

通过对藏北高原西北部结则茶卡湖泊及其沿岸地质地貌调查,发现其为一个富含硼、锂、钾、锶的封闭型盐湖,沿岸海拔4 850 m拔湖325 m有一条明显的高位湖岸线,该湖岸线到湖面之间有六级湖积阶地发育,六级以上阶地保存零星。沿湖岸不同高度上的湖积物U系年龄分别为(14.2±1.2)ka BP(T2)、(38.0±3.5)ka BP(T4)和(41.6±3.2)ka BP(T5)。湖面下降的幅度是藏北高原迄今所知最大。根据湖面平均下降速度推算高位湖岸线和高位湖积层的形成年龄在(120~90)ka BP,与东部的纳木错和西部的甜水海基本一致,说明藏北高原在晚更新世初期有一个明显的泛湖期。大约100 ka BP结则茶卡湖面开始下降,晚更新世以来湖泊演化是在封闭体系干旱环境下进行的,盐湖形成于14 ka BP左右,藏北高原在晚更新世以来气候变化可能为自西向东逐渐变为干寒。     

地下水资源开采潜力评价的综合因素浅析

浙江省现有地下水开采潜力的分析与评价方法存在着不足,本文结合乐清市地下水资源开发利用规划实例,浅析浙江省地下水资源开采潜力分析与评价的综合影响因素.