Seasonal variations in infiltrability of moss‐dominated biocrusts on aeolian sand and loess soil in the Chinese Loess Plateau

Biocrust effects on soil infiltration have attracted increasing attention in dryland ecosystems, but their seasonal variations in infiltrability have not yet been well understood. On the Chinese Loess Plateau, soil infiltrability indicated by saturated hydraulic conductivity (K_s) of biocrusts and bare soil, both on aeolian sand and loess soil, was determined by disc infiltrometer in late spring (SPR), midsummer (SUM), and early fall (FAL). Then their correlations with soil biological and physiochemical properties and water repellency index (RI) were analysed. The results showed that the biocrusts significantly decreased K_s both on sand during SPR, SUM, and FAL (by 43%, 66%, and 35%, respectively; P < .05) and on loess (by 42%, 92%, and 10%, respectively; P <.05). As compared with the bare soil, the decreased K_s in the biocrusted surfaces was mostly attributed to the microorganism biomass and also to the increasing content of fine particles and organic matter. Most importantly, both the biocrusts and bare soil exhibited significant (F ≥ 11.89, P ≤ .003) seasonal variations in K_s, but their patterns were quite different. Specifically, the K_s of bare soil gradually decreased from SPR to SUM (32% and 42% for sand and loess, respectively) and FAL (29% and 39%); the K_s of biocrusts also decreased from SPR to SUM (59% and 92%) but then increased in FAL (36% and 588%). Whereas the seasonal variations in K_s of the biocrusts were closely correlated with the seasonal variations in RI, the RI values were not high enough to point at hydrophobicity. Instead of that, the seasonal variations of K_s were principally explained by the changes in the crust biomass and possibly by the microbial exopolysaccharides. We conclude that the biocrusts significantly decreased soil infiltrability and exhibited a different seasonal variation pattern, which should be carefully considered in future analyses of hydropedological processes.     

Bound and mobile soil water isotope ratios are affected by soil texture and mineralogy,whereas extraction method influences their measurement

Questions persist about interpreting isotope ratios of bound and mobile soil water pools, particularly relative to clay content and extraction conditions. Interactions between pools and resulting extracted water isotope composition are presumably related to soil texture, yet few studies have manipulated the bound pool to understand its influence on soil water processes. Using a series of drying and spiking experiments, we effectively labelled bound and mobile water pools in soils with varying clay content. Soils were first vacuum dried to remove residual water, which was then replaced with heavy isotope-enriched water prior to oven drying and spiking with heavy isotope-depleted water. Water was extracted via centrifugation or cryogenic vacuum distillation (at four temperatures) and analysed for oxygen and hydrogen isotope ratios via isotope ratio mass spectrometry. Water from centrifuged samples fell along a mixing line between the two added waters but was more enriched in heavy isotopes than the depleted label, demonstrating that despite oven drying, a residual pool remains and mixes with the mobile water. Soils with higher clay + silt content appeared to have a larger bound pool. Water from vacuum distillation samples have a significant temperature effect, with high temperature extractions yielding progressively more heavy isotope-enriched values, suggesting that Rayleigh fractionation occurred at low temperatures in the vacuum line. By distinctly labelling bound and mobile soil water pools, we detected interactions between the two that were dependent on soil texture. Although neither extraction method appeared to completely extract the combined bound and mobile (total water) pool, centrifugation and high temperature cryogenic vacuum distillations were comparable for both δ²H and δ¹⁸O of soil water isotope ratios.     

The exponential decline in saturated hydraulic conductivity with depth: a novel method for exploring its effect on water flow paths and transit time distribution

The strong vertical gradient in soil and subsoil saturated hydraulic conductivity is characteristic feature of the hydrology of catchments. Despite the potential importance of these strong gradients, they have proven difficult to model using robust physically based schemes. This has hampered the testing of hypotheses about the implications of such vertical gradients for subsurface flow paths, residence times and transit time distribution. Here we present a general semi‐analytical solution for the simulation of 2D steady‐state saturated‐unsaturated flow in hillslopes with saturated hydraulic conductivity that declines exponentially with depth. The grid‐free solution satisfies mass balance exactly over the entire saturated and unsaturated zones. The new method provides continuous solutions for head, flow and velocity in both saturated and unsaturated zones without any interpolation process as is common in discrete numerical schemes. This solution efficiently generates flow pathlines and transit time distributions in hillslopes with the assumption of depth‐varying saturated hydraulic conductivity. The model outputs reveal the pronounced effect that changing the strength of the exponential decline in saturated hydraulic conductivity has on the flow pathlines, residence time and transit time distribution. This new steady‐state model may be useful to others for posing hypotheses about how different depth functions for hydraulic conductivity influence catchment hydrological response. Copyright © 2016 John Wiley & Sons, Ltd.     

Micromechanical approach to swelling behavior of capillary-porous media with coupled physics

A detailed multiscale analysis is presented of the swelling phenomenon in unsaturated clay-rich materials in the linear regime through homogenization. Herein, the structural complexity of the material is formulated as a three-scale, triple porosity medium within which microstructural information is transmitted across the various scales, leading ultimately to an enriched stress-deformation relation at the macroscopic scale. As a side note, such derived relationship leads to a tensorial stress partitioning that is reminiscent of a Terzaghi-like effective stress measure. Otherwise, a major result that stands out from previous works is the explicit expression of swelling stress and capillary stress in terms of micromechanical interactions at the very fine scale down to the clay platelet level, along with capillary stress emerging due to interactions between fluid phases at the different scales, including surface tension, pore size, and morphology. More importantly, the swelling stress is correlated with the disjoining forces due to electrochemical effects of charged ions on clay minerals and van der Waals forces at the nanoscale. The resulting analytical expressions also elucidate the role of the various physics in the deformational behavior of clayey material. Finally, the capability of the proposed formulation in capturing salient behaviors of unsaturated expansive clays is illustrated through some numerical examples.     

ISA model: A constitutive model for soils with yield surface in the intergranular strain space

In this article, a new constitutive model for soils is proposed. It is formulated by means of plasticity, but in contrast to the precedent works, it presents a yield function describing a surface within the intergranular strain space. This latter is a state variable providing information of the recent strain history. An expression for the plastic strain rate has been proposed to guarantee the stress rate continuity. Under the application of medium or large strain amplitudes, the constitutive equation becomes independent of the intergranular strain and delivers a mathematical structure similar to some Karlsruhe hypoplastic models. Some simulations of monotonic and cyclic triaxial test are provided to evaluate and analyze the model performance. Copyright © 2015 John Wiley & Sons, Ltd.     

Diffusive mass exchange of non‐reactive substances in dual‐porosity porous systems – column experiments under saturated conditions

Diffusive mass exchange into immobile water regions within heterogeneous porous aquifers influences the fate of solutes. The percentage of immobile water is often unidentified in natural aquifers though. Hence, the mathematical prediction of solute transport in such heterogeneous aquifers remains challenging. The objective of this study was to find a simple analytical model approach that allows quantifying properties of mobile and immobile water regions and the portion of immobile water in a porous system. Therefore, the Single Fissure Dispersion Model (SFDM), which takes into account diffusive mass exchange between mobile and immobile water zones, was applied to model transport in well‐defined saturated dual‐porosity column experiments. Direct and indirect model validation was performed by running experiments at different flow velocities and using conservative tracer with different molecular diffusion coefficients. In another column setup, immobile water regions were randomly distributed to test the model applicability and to determine the portion of immobile water. In all setups, the tracer concentration curves showed differences in normalized maximum peak concentration, tailing and mass recovery according to their diffusion coefficients. These findings were more pronounced at lower flow rates (larger flow times) indicating the dependency of diffusive mass exchange into immobile water regions on tracers' molecular diffusion coefficients. The SFDM simulated all data with high model efficiency. Successful model validation supported the physical meaning of fitted model parameters. This study showed that the SFDM, developed for fissured aquifers, is applicable in porous media and can be used to determine porosity and volume of regions with immobile water. Copyright © 2015 John Wiley & Sons, Ltd.     

北部湾沉积物粘土矿物分布特征及其环境意义

本文以南海北部湾SO-31沉积柱为研究对象,研究了14C年代学和粘土矿物学特征,并对部分全球气候事件进行了对比,为古环境、古气候的恢复提供基础资料,也为全球重大气候事件在该区域的响应提供信息。结果显示全新世以来地层沉积正常,平均沉积速率为0.57mm/a。粘土成分主要由蒙脱石、伊利石、高岭石和绿泥石组成,组合类型为蒙脱石-伊利石-高岭石-绿泥石型。全新世以来环境气候演变可划分为五个阶段:低温期阶段、干湿交替的寒冷气候阶段、逐渐升温阶段、干旱温暖气候阶段、湿热阶段。气候在每个阶段背景下还存在一些次级波动,总体趋势为干湿交替,温度逐渐上升。由于海域环境及矿物指标的影响,北部湾SO-31沉积柱粘土矿物记录的降温事件时间比其他指标记录的新仙女木降温事件发生时间滞后500~800a。     

贵州威宁麻窝山岩溶盆地沉积物中粘土矿物特征及其古气候指示意义

利用X射线衍射等分析方法,对贵州威宁麻窝山岩溶盆地沉积物中粘土矿物进行研究,获得了粘土矿物的组成、相对含量、结晶度指数、化学指数以及相对含量比值等方面的信息,探讨该盆地的沉积环境演化特征。研究表明,粘土矿物以伊/蒙混层矿物(30%~75%)和绿泥石(10%~45%)为主,高岭石(5%~20%)和伊利石(10%)次之,基本不含蒙脱石。伊/蒙混层矿物和绿泥石含量呈负相关性。高岭石含量在剖面底部较稳定,在484cm上部出现较大的波动。伊利石结晶度和化学指数与伊/蒙混层矿物含量具有较好的正相关性,且在第Ⅰ阶段相对于其它指标较敏感,出现波动。通过综合对比分析,将研究区剖面划分为4个阶段,其是3个冷暖干湿的大循环和若干个小循环气候环境的物质记录。     

Axial friction response of full-scale pipes in soft clays

The axial friction response of subsea pipelines in soft clays is a very important aspect for designers of subsea pipelines but the response is not well understood so far. There is a pressing need for the comprehension of the response. In this paper, model tests are performed using full-scale pipes coated with polyethylene (PE) to study the effects of the set-up period, the pipe diameter, the buried depth of the pipe, the shear strength of soft clays and the loading rate on the axial friction response of pipelines in soft clays. The variations of the axial friction coefficient are analyzed using the effective stress method based on model test results. The results show that the axial friction resistance increases with the increasing pipe diameter but the effect of the pipe diameter on the axial friction coefficient can be neglected. The ultimate axial resistance also increases with the increase of the buried depth of pipelines, the undrained shear strength of soft clays and the loading rate. The axial friction coefficient increases with the increasing loading rate. However, the axial friction coefficient decreases with the increasing buried depth. The method to determine the axial friction coefficient is developed by analyzing model test results, which considers the effects of the diameter, the buried depth, the undrained shear strength of soft clays and the loading rate. The study results not only extend the industry data base but also supply a basis to determine the axial friction coefficient of PE-coated pipes in soft clays for ocean engineering geological investigations.     

Incremental elastoplastic FEM for simulating the deformation process of suction caissons subjected to cyclic loads in soft clays

This paper presents an incremental elastoplastic finite element method (FEM) to simulate the undrained deformation process of suction caisson foundations subjected to cyclic loads in soft clays. The method is developed by encoding the total-stress-based bounding surface model proposed by the authors in the ABAQUS software package. According to the model characteristics, elastoplastic stress states associated with the incremental strains of each iteration are determined using the sub-incremental explicit Euler algorithm, and the state parameters describing the cyclic accumulative rates of strains are updated by setting state variables during the calculations. The radial fallback method is also proposed to modify the stress states outside the bounding surface to the surface during determination of the elastoplastic stress states. The stress reversals of soil elements are judged by the angle between the incremental deviatoric stress and the exterior normal vector at the image stress point on the bounding surface to update the mapping centre and state variables during cyclic loading. To assess the general validity of the method, the reduced scale model tests and centrifuge tests of suction caissons subjected to cyclic loads are simulated using the method. Predictions are in relative good agreement with test results. Compared with the limit equilibrium and quasi-static methods, the method can not only determine the cyclic bearing capacity, but can also analyse the deformation process and the failure mechanisms of suction caisson under cyclic loads in soft clays.