New graphical methods for estimating aquifer hydraulic parameters using pumping tests with exponentially decreasing rates

Actual pumping tests may involve continuously decreasing rates over a certain period of time, and the hydraulic conductivity (K) and specific storage (S_s) of the tested confined aquifer cannot be interpreted from the classical constant‐rate test model. In this study, we revisit the aquifer drawdown characteristics of a pumping test with an exponentially decreasing rate using the dimensionless analytical solution for such a variable‐rate model. The drawdown may decrease with time for a short period of time at intermediate pumping times for such pumping tests. A larger ratio of initial to final pumping rate and a smaller radial distance of the observation well will enhance the decreasing feature. A larger decay constant results in an earlier decrease, but it weakens the extent of such a decrease. Based on the proposed dimensionless transformation, we have proposed two graphical methods for estimating K and S_s of the tested aquifer. The first is a new type curve method that does not employ the well function as commonly done in standard type curve analysis. Another is a new analytic method that takes advantage of the decreasing features of aquifer drawdown during the intermediate pumping stage. We have demonstrated the applicability and robustness of the two new graphical methods for aquifer characterization through a synthetic pumping test.     

Size and spatial distribution of loess slides on the Chinese Loess Plateau

ABSTRACT

The size and spatial distribution of loess slides are important for estimating the yield of eroded materials and determining the landslide risk. While previous studies have investigated landslide size distributions, the spatial distribution pattern of landslides at different spatial scales is poorly understood. The results indicate that the loess slide distribution exhibits a power-law scaling across a range of the size distribution. The mean landslide size and size distribution in the different geomorphic types are different. The double Pareto and inverse gamma functions can coincide well with the empirical probability distribution of the loess slide areas and can quantitatively reveal the rollover location, maximum probability, and scaling exponents. The frequency of loess slides increases with mean monthly precipitation. Moreover, point distance analysis showed that > 80% of landslides are located < 3 km from other loess slides. We found that the loess slides at the two study sites (Zhidan and Luochuan County) in northern Shaanxi Province, China show a significant clustered distribution. Furthermore, analysis results of the correlated fractal dimension show that the landslides exhibit a dispersed distribution at smaller spatial scales and a clustered distribution at larger spatial scales.     

Parallelized combined finite-discrete element (FDEM) procedure using multi-GPU with CUDA

This paper focuses on the efficiency of finite discrete element method (FDEM) algorithmic procedures in massive computers and analyzes the time-consuming part of contact detection and interaction computations in the numerical solution. A detailed operable GPU parallel procedure was designed for the element node force calculation, contact detection, and contact interaction with thread allocation and data access based on the CUDA computing. The emphasis is on the parallel optimization of time-consuming contact detection based on load balance and GPU architecture. A CUDA FDEM parallel program was developed with the overall speedup ratio over 53 times after the fracture from the efficiency and fidelity performance test of models of in situ stress, UCS, and BD simulations in Intel i7-7700K CPU and the NVIDIA TITAN Z GPU. The CUDA FDEM parallel computing improves the computational efficiency significantly compared with the CPU-based ones with the same reliability, providing conditions for achieving larger-scale simulations of fracture.     

Temporal variation of soil erosion resistance on sloping farmland during the growth stages of maize (Zea mays L.)

Maize growth has great effects on soil properties and thus likely induces the changes in soil erosion resistance on sloping farmland. However, temporal variation of soil erosion resistance during the growth stages of maize is still unclear in the mountainous yellow soil area where maize is the dominant crop. In this study, four maize plots (MP) and four bare land plots (CK) were conducted to investigate soil erosion resistance, and multiple indicators of soil erosion resistance were measured including the total soil anti-scourability (TAS), mean weight diameter (MWD), soil erodibility K factor and soil shear strength (SH). A comprehensive soil erosion resistance index (CSERI) was employed to quantify the temporal variation of soil erosion resistance during the growth stages of maize (seedling stage, SS; jointing stage, JS; tasselling stage, TS; maturing stage, MS). The results showed that TAS, MWD, SH increased significantly with maize growth and SH decreased when at MS. But K factor decreased significantly over time. CSERI increased significantly during the growth stages of maize and the CSERI of JS, TS, MS increased on average by 74.72, 180.68 and 234.57% than that of SS. Compared to CK, CSERI of MP increased by 49.90, 66.82, 55.60 and 38.61% during the growth stages of maize. The temporal variation of soil erosion resistance was closely related to the changes in maize cover, maize roots and soil organic carbon. The findings demonstrated that it is necessary to consider the temporal variation of soil erosion resistance in the mountainous yellow soil area.     

Constitutive modelling of fine-grained gassy soil: A composite approach

Fine-grained marine sediments containing large undissolved gas bubbles are widely distributed around the world. Presence of the bubbles could degrade the undrained shear strength (s_u ) of the soil, when the gas pressure u_g is relatively high as compared with the effective stress in the saturated soil matrix. Meanwhile, the addition of bubbles may also increase s_u when the difference between u_g and pore water pressure u_w becomes smaller than the water entry value, causing partial water drainage from the saturated matrix into the bubbles (bubble flooding) during globally undrained shearing. A new constitutive model for describing the two competing effects on the stress-strain relationship of fine-grained gassy soil is proposed within the framework of critical state soil mechanics. The gassy soil is considered as a three-phase composite material with compressible cavities, which allows water entry from the saturated matrix. Bubble flooding is modelled by introducing an additional positive volumetric strain increment of the saturated clay matrix, which is dependent on the difference between pore gas and pore water pressure based on experimental observations. A modified hardening law based on that of the modified Cam clay model is employed, which in conjunction with the expression for bubble flooding, can describe both the detrimental and beneficial effects of gas bubbles on soil strength and plastic hardening in shear. Only two extra parameters in addition to those in the modified Cam clay model are used. It is shown that the key features of the stress-strain relationship of three fine-grained gassy soils can be reproduced satisfactorily.     

Evaluating suspended sediment and turbidity reduction projects in a glacially conditioned catchment through dynamic regression and fluvial process-based modelling

Elevated turbidity (T_n) and suspended sediment concentrations (SSC) during and following flood events can degrade water supply quality and aquatic ecosystem integrity. Streams draining glacially conditioned mountainous terrain, such as those in the Catskill Mountains of New York State, are particularly susceptible to high levels of T_n and SSC sourced from erosional contact with glacial-related sediment. This study forwards a novel approach to evaluate the effectiveness of stream restoration best management practices (BMPs) meant to reduce stream T_n and SSC, and demonstrates the approach within the Stony Clove sub-basin of the Catskills, a water supply source for New York City. The proposed approach is designed to isolate BMP effects from natural trends in T_n and SSC caused by trends in discharge and shifts in average T_n or SSC per unit discharge (Q) following large flood events. We develop Dynamic Linear Models (DLMs) to quantify how T_n-Q and SSC-Q relationships change over time at monitoring stations upstream and downstream of BMPs within the Stony Clove and in three other sub-basins without BMPs, providing observational evidence of BMP effectiveness. A process-based model, the River Erosion Model, is then developed to simulate natural, hydrology-driven SSC-Q dynamics in the Stony Clove sub-basin (absent of BMP effects). We use DLMs to compare the modelled and observed SSC-Q dynamics and isolate the influence of the BMPs. Results suggest that observed reductions in SSC and T_n in the Stony Clove sub-basin have been driven by a combination of declining streamflow and the installed BMPs, confirming the utility of the BMPs for the monitored hydrologic conditions.     

Variability in the vertical hyporheic water exchange affected by hydraulic conductivity and river morphology at a natural confluent meander bend

River confluences and their associated tributaries are key morphodynamic nodes that play important roles in controlling hydraulic geometry and hyporheic water exchange in fluvial networks. However, the existing knowledge regarding hyporheic water exchange associated with river confluence morphology is relatively scarce. On January 14 and 15, 2016, the general hydraulic and morphological characteristics of the confluent meander bend (CMB) between the Juehe River and the Haohe River in the southern region of Xi'an City, Shaanxi Province, China, were investigated. The patterns and magnitudes of vertical hyporheic water exchange (VHWE) were estimated based on a one‐dimensional heat steady‐state model, whereas the sediment vertical hydraulic conductivity (K_v) was calculated via in situ permeameter tests. The results demonstrated that 6 hydrodynamic zones and their extensions were observed at the CMB during the test period. These zones were likely controlled by the obtuse junction angle and low momentum flux ratio, influencing the sediment grain size distribution of the CMB. The VHWE patterns at the test site during the test period mostly showed upwelling flow dominated by regional groundwater discharging into the river. The occurrence of longitudinal downwelling and upwelling patterns along the meander bend at the CMB was likely subjected to the comprehensive influences of the local sinuosity of the meander bend and regional groundwater discharge and finally formed regional and local flow paths. Additionally, in dominated upwelling areas, the change in VHWE magnitudes was nearly consistent with that in K_v values, and higher values of both variables generally occurred in erosional zones near the thalweg paths of the CMB, which were mostly made up of sand and gravel. This was potentially caused by the erosional and depositional processes subjected to confluence morphology. Furthermore, lower K_v values observed in downwelling areas at the CMB were attributed to sediment clogging caused by local downwelling flow. The confluence morphology and sediment K_v are thus likely the driving factors that cause local variations in the VHWE of fluvial systems.