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.     

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.     

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.     

Subsurface flow measurements using passive flux meters in variably-saturated cold-regions landscapes

To date, passive flux meters have predominantly been applied in temperate environments for tracking the movement of contaminants in groundwater. This study applies these instruments to reduce uncertainty in (typically instantaneous) flux measurements made in a low-gradient, wetland dominated, discontinuous permafrost environment. This method supports improved estimation of unsaturated and over-winter subsurface flows which are very difficult to quantify using hydraulic gradient-based approaches. Improved subsurface flow estimates can play a key role in understanding the water budget of this landscape.     

PS测井在机场隧道地基勘查中的应用

PS测井技术是地震勘探方法之一,也是一种简便、快速、准确的原位测试技术。该文介绍了单孔法PS测井的原理及其在铁路工程勘察中的应用,主要包括根据等效剪切波,进行建筑抗震场地类别划分,铁路工程抗震场地类别划分。根据岩土动力学参数,达到评价岩体质量和划分围岩类别目的,利用剪切波速法估算岩土的承载力基本值。通过青岛机场隧道场地的工程实例说明了PS测井技术在岩土工程勘察设计中的应用情况及应用效果。     

Effects of the geometry of two‐dimensional fractures on their hydraulic aperture and on the validity of the local cubic law

Flow through rough fractures is investigated numerically in order to assess the validity of the local cubic law for different fracture geometries. Two‐dimensional channels with sinusoidal walls having different geometrical properties defined by the aperture, the amplitude, and the wavelength of the walls' corrugations, the corrugations asymmetry, and the phase shift between the two walls are considered to represent different fracture geometries. First, it is analytically shown that the hydraulic aperture clearly deviates from the mean aperture when the walls' roughness, the phase shift, and/or the asymmetry between the fracture walls are relatively high. The continuity and the Navier–Stokes equations are then solved by means of the finite element method and the numerical solutions compared to the theoretical predictions of the local cubic law. Reynolds numbers ranging from 0.066 to 66.66 are investigated so as to focus more particularly on the effect of flow inertial effects on the validity of the local cubic law. For low Reynolds number, typically less than 15, the local cubic law properly describes the fracture flow, especially when the fracture walls have small corrugation amplitudes. For Reynolds numbers higher than 15, the local cubic law is valid under the conditions that the fracture presents a low aspect ratio, small corrugation amplitudes, and a moderate phase lag between its walls.     

Revisiting the existence of an effective stress for wet granular soils with micromechanics

A possible effective stress variable for wet granular materials is numerically investigated based on an adapted discrete element method (DEM) model for an ideal three‐phase system. The DEM simulations consider granular materials made of nearly monodisperse spherical particles, in the pendular regime with the pore fluid mixture consisting of distinct water menisci bridging particle pairs. The contact force‐related stress contribution to the total stresses is isolated and tested as the effective stress candidate for dense or loose systems. It is first recalled that this contact stress tensor is indeed an adequate effective stress that describes stress limit states of wet samples with the same Mohr‐Coulomb criterion associated with their dry counterparts. As for constitutive relationships, it is demonstrated that the contact stress tensor used in conjunction with dry constitutive relations does describe the strains of wet samples during an initial strain regime but not beyond. Outside this so‐called quasi‐static strain regime, whose extent is much greater for dense than loose materials, dramatic changes in the contact network prevent macroscale contact stress‐strain relationships to apply in the same manner to dry and unsaturated conditions. The presented numerical results also reveal unexpected constitutive bifurcations for the loose material, related to stick‐slip macrobehavior.     

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.     

全球气候变化下南海诸岛保护优先区识别分析

全球变化下,珊瑚礁保护区是保护生物多样性、增强珊瑚礁对气候变暖抵抗力的有效方式,而维持珊瑚礁弹性是其核心内容。针对珊瑚礁最具有威胁性的热压力因子,基于南海1982—2009年卫星观测海表面温度（SST）数据和CMIP5加拿大地球系统模式CanESM2模型预估的2006—2100年南海SST数据构建热压力强度模型,从维持珊瑚礁弹性的角度识别IPCC RCP 4.5和RCP 8.5情景下南海诸岛保护优先区。结果表明：RCP 4.5和RCP 8.5情景下13%左右的南海诸岛珊瑚礁识别为保护优先区。根据热压力强度与珊瑚抵抗力及避难所关系,西沙群岛七连屿和晋卿岛近年观测与未来预估的热压力强度均比较低,在保障其服务功能的基础上建议实施完全保护;东沙群岛东沙环礁和中沙环礁排洪滩近年观测急性热压力强度较高但未来预估热压力强度较低,建议实施50%禁止利用保护;中沙群岛黄岩岛近年观测和未来预估的急性热压力强度均比较低,建议实施50%多用途保护。南沙群岛有14%左右的珊瑚礁识别为保护优先区,根据其热压力强度可实施30%~100%禁止利用保护或30%~50%多用途保护。RCP 4.5和RCP 8.5情景下的南海诸岛保护优先区及保护对策,可为维持珊瑚礁生态弹性及应对全球气候变化提供重要的参考价值。     

Parallelization of the flow-path network model using a particle-set strategy

High-performance simulation of flow dynamics remains a major challenge in the use of physical-based, fully distributed hydrologic models. Parallel computing has been widely used to overcome efficiency limitation by partitioning a basin into sub-basins and executing calculations among multiple processors. However, existing partition-based parallelization strategies are still hampered by the dependency between inter-connected sub-basins. This study proposed a particle-set strategy to parallelize the flow-path network (FPN) model for achieving higher performance in the simulation of flow dynamics. The FPN model replaced the hydrological calculations on sub-basins with the movements of water packages along the upstream and downstream flow paths. Unlike previous partition-based task decomposition approaches, the proposed particle-set strategy decomposes the computational workload by randomly allocating runoff particles to concurrent computing processors. Simulation experiments of the flow routing process were undertaken to validate the developed particle-set FPN model. The outcomes of hourly outlet discharges were compared with field gauged records, and up to 128 computing processors were tested to explore its speedup capability in parallel computing. The experimental results showed that the proposed framework can achieve similar prediction accuracy and parallel efficiency to that of the Triangulated Irregular Network (TIN)-based Real-Time Integrated Basin Simulator (tRIBS).