Modelling and experimental study of coupled exchange of p,p′‐DDE and colloids in a flume simulated stream–streambed system

Stream–subsurface exchange plays a significant role in the fate and transport of contaminants in streams. It has been modelled explicitly by considering fundamental processes such as hydraulic exchange, colloid filtration, and contaminant interactions with streambed sediments and colloids. The models have been successfully applied to simulate the transport of inorganic metals and nutrients. In this study, laboratory experiments were conducted in a recirculating flume to investigate the exchange of a hydrophobic organic contaminant, p,p′‐dichloro‐diphenyl‐dichloroethane (DDE), between a stream and a quartz sand bed. A previously developed process‐based multiphase exchange model was modified by accounting for the p,p′‐DDE kinetic adsorption to and desorption from the bed sediments/colloids and was applied to interpret the experimental results. Model input parameters were obtained by conducting independent small‐scale batch experiments. Results indicate that the immobilization of p,p′‐DDE in the quartz sand bed can occur under representative natural stream conditions. The observed p,p′‐DDE exchange was successfully simulated by the process‐based model. The model sensitivity analysis results show that the exchange of p,p′‐DDE can be sensitive to either the sediment sorption/desorption parameters or colloidal parameters depending on the experimental conditions tested. For the experimental conditions employed here, the effect of colloids on contaminant transport is expected to be minimal, and the stream–subsurface exchange of p,p′‐DDE is dominated by the interaction of p,p′‐DDE with bed sediment. The work presented here contributes to a better mechanistic understanding of the complex transport process that hydrophobic organic contaminants undergo in natural streams and to the development of reliable, predictive models for the assessment of impacted streams. Copyright © 2015 John Wiley & Sons, Ltd.     

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.     

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.     

Numerical and theoretical research on stress distribution in the loosening zone of the trapdoor problem

The traditional theory of soil arching effect was developed on the assumption that stress distribution in the loosening zone is uniform. However, because of the deflection of principal stress' direction, the stress distribution in the loosening zone is actually ununiform. For the evaluation of principal stress axis deflection and stress redistribution, a discrete element method numerical model of trapdoor problem is established for the simulation of soil arching effect. Based on the numerical results, an arc shape of major principal stress trajectory and uniform horizontal stress distribution at the same depth of the loosening zone are adopted. An analytical model is raised to estimate the average loosening earth pressure acting on the trapdoor and stress distribution in the loosening zone at a limit state. In addition, comparison studies are carried out between the predictions of the proposed solutions and discrete element method numerical results as well as available model test results, thereby validating the accuracy of the proposed theoretical model. Both numerical and theoretical results indicate that the vertical stress distribution in the loosening zone is obviously ununiform. The load acting in the middle of loosening zone is transferred toward two sides so that the vertical stress distribution in loosening zone is concave.     

Water mining from the deep critical zone by apple trees growing on loess

There have been significant recent advances in understanding the ecohydrology of deep soil. However, the links between root development and water usage in the deep critical zone remains poorly understood. To clarify the interaction between water use and root development in deep soil, we investigated soil water and root profiles beyond maximum rooting depth in five apple orchards planted on farmland with stand ages of 8, 11, 15, 18, and 22 years in a subhumid region on the Chinese Loess Plateau. Apple trees rooted progressively deeper for water with increasing stand age and reached 23.2 ± 0.8 m for the 22‐year‐old trees. Soil water deficit in deep soil increased with tree age and was 1,530 ± 43 mm for a stand age of 22 years. Measured root deepening rate was far great than the reported pore water velocity, which demonstrated that trees are mining resident old water. The deficits are not replenished during the life‐span of the orchard, showing a one‐way mining of the critical zone water. The one‐way root water mining may have changed the fine root profile from an exponential pattern in the 8‐year‐old orchard to a relative uniform distribution in older orchards. Our findings enhance our understanding of water‐root interaction in deep soil and reveal the unintended consequences of critical zone dewatering during the lifespan of apple trees.     

Root‐zone moisture replenishment in a native vegetated catchment under Mediterranean climate

The root‐zone moisture replenishment mechanisms are key unknowns required to understand soil hydrological processes and water sources used by plants. Temporal patterns of root‐zone moisture replenishment reflect wetting events that contribute to plant growth and survival and to catchment water yield. In this study, stable oxygen and hydrogen isotopes of twigs and throughfall were continuously monitored to characterize the seasonal variations of the root‐zone moisture replenishment in a native vegetated catchment under Mediterranean climate in South Australia. The two studied hillslopes (the north‐facing slope [NFS] and the south‐facing slope [SFS]) had different environmental conditions with opposite aspects. The twig and throughfall samples were collected every ~20 days over 1 year on both hillslopes. The root‐zone moisture replenishment, defined as percentage of newly replenished root‐zone moisture as a complement to antecedent moisture for plant use, calculated by an isotope balance model, was about zero (±25% for the NFS and ± 15% for the SFS) at the end of the wet season (October), increased to almost 100% (±26% for the NFS and ± 29% for the SFS) after the dry season (April and May), then decreased close to zero (±24% for the NFS and ± 28% for the SFS) in the middle of the following wet season (August). This seasonal pattern of root‐zone moisture replenishment suggests that the very first rainfall events of the wet season were significant for soil moisture replenishment and supported the plants over wet and subsequent dry seasons, and that NFS completed replenishment over a longer time than SFS in the wet season and depleted the root zone moisture quicker in the dry season. The stable oxygen isotope composition of the intraevent samples and twigs further confirms that rain water in the late wet season contributed little to root‐zone moisture. This study highlights the significant role of the very first rain events in the early wet season for ecosystem and provides insights to understanding ecohydrological separation, catchment water yield, and vegetation response to climate changes.     

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.     

New P–T constraints on the Tamayen glaucophane‐bearing rocks,eastern Taiwan: Perple_X modelling results and geodynamic implications

New pseudosection modelling was applied to better constrain the P–T conditions and evolution of glaucophane‐bearing rocks in the Tamayen block of the Yuli belt, recognized as the world's youngest known blueschist complex. Based on the predominant clinoamphibole, textural relationships, and mineral compositions, these glaucophane‐bearing high‐P rocks can be divided into four types. We focused on the three containing garnet. The chief phase assemblages are (in decreasing mode): amphibole + quartz + epidote + garnet + chlorite + rutile/titanite (Type‐I), phengite + amphibole + quartz + garnet + chlorite + epidote + titanite + biotite + magnetite (Type‐II), and amphibole + quartz + albite + epidote + garnet + rutile + hematite + titanite (Type‐III). Amphibole exhibits compositional zoning from core to rim as follows: glaucophane → pargasitic amphibole → actinolite (Type‐I), barroisite → Mg‐katophorite/taramite → Fe‐glaucophane (Type‐II), glaucophane → winchite (Type‐III). Using petrographic data, mineral compositions and Perple_X modelling (pseudosections and superimposed isopleths), peak P–T conditions were determined as 13 ± 1 kbar and 550 ± 40 °C for Type‐I, 10.5 ± 0.5 kbar and 560 ± 30 °C for Type‐II (thermal peak) and 11 ± 1 kbar and 530 ± 30 °C for Type‐III. The calculations yield higher pressures and temperatures than previously thought; the difference is ~1–6 kbar and 50–200 °C. The three rock types record similar P–T retrograde paths with clockwise trajectories; all rocks followed trajectories with substantial pressure decrease under near‐isothermal conditions (Type‐I and Type‐III), with the probable exception of Type‐II where decompression followed colder geotherms. The P–T paths suggest a tectonic environment in which the rocks were exhumed from maximum depths of ~45 km within a subduction channel along a relative cold geothermal gradient of ~11–14 °C km^?1.     

Monitoring over‐winter soil water dynamics in a freezing and snow‐covered environment using a thermally insulated tensiometer

Water potential below a frozen soil layer was continuously monitored over an entire winter period (using thermally insulated tensiometers sheltered in a heated chamber) along with other soil, snow and atmospheric variables. In early winter, the freezing front advanced under a thin snow cover, inducing upward soil water flow in the underlying unfrozen soil. The freezing front started to retreat when the snow cover became thick enough to insulate the soil, resulting in the reversal of the flow direction in the unfrozen zone. These data provide a clear illustration of soil water dynamics, which have rarely been monitored with a tensiometer. Copyright © 2005 John Wiley & Sons, Ltd.