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.     

Implementation of a coupled plastic damage distinct lattice spring model for dynamic crack propagation in geomaterials

This paper studies dynamic crack propagation by employing the distinct lattice spring model (DLSM) and 3‐dimensional (3D) printing technique. A damage‐plasticity model was developed and implemented in a 2D DLSM. Applicability of the damage‐plasticity DLSM was verified against analytical elastic solutions and experimental results for crack propagation. As a physical analogy, dynamic fracturing tests were conducted on 3D printed specimens using the split Hopkinson pressure bar. The dynamic stress intensity factors were recorded, and crack paths were captured by a high‐speed camera. A parametric study was conducted to find the influences of the parameters on cracking behaviors, including initial and peak fracture toughness, crack speed, and crack patterns. Finally, selection of parameters for the damage‐plasticity model was determined through the comparison of numerical predictions and the experimentally observed cracking features.     

Assessment of alternative adsorption models and global sensitivity analysis to characterize hexavalent chromium loss from soil to surface runoff

We investigate our ability to assess transfer of hexavalent chromium, Cr(VI), from the soil to surface runoff by considering the effect of coupling diverse adsorption models with a two‐layer solute transfer model. Our analyses are grounded on a set of two experiments associated with soils characterized by diverse particle size distributions. Our study is motivated by the observation that Cr(VI) is receiving much attention for the assessment of environmental risks due to its high solubility, mobility, and toxicological significance. Adsorption of Cr(VI) is considered to be at equilibrium in the mixing layer under our experimental conditions. Four adsorption models, that is, the Langmuir, Freundlich, Temkin, and linear models, constitute our set of alternative (competing) mathematical formulations. Experimental results reveal that the soil samples characterized by the finest grain sizes are associated with the highest release of Cr(VI) to runoff. We compare the relative abilities of the four models to interpret experimental results through maximum likelihood model calibration and four model identification criteria (i.e., the Akaike information criteria [AIC and AIC_C] and the Bayesian and Kashyap information criteria). Our study results enable us to rank the tested models on the basis of a set of posterior weights assigned to each of them. A classical variance‐based global sensitivity analysis is then performed to assess the relative importance of the uncertain parameters associated with each of the models considered, within subregions of the parameter space. In this context, the modelling strategy resulting from coupling the Langmuir isotherm with a two‐layer solute transfer model is then evaluated as the most skilful for the overall interpretation of both sets of experiments. Our results document that (a) the depth of the mixing layer is the most influential factor for all models tested, with the exception of the Freundlich isotherm, and (b) the total sensitivity of the adsorption parameters varies in time, with a trend to increase as time progresses for all of the models. These results suggest that adsorption has a significant effect on the uncertainty associated with the release of Cr(VI) from the soil to the surface runoff component.     

Predicting quantiles of water quality from catchment characteristics

Water quality is often highly variable both in space and time, which poses challenges for modelling the more extreme concentrations. This study developed an alternative approach to predicting water quality quantiles at individual locations. We focused on river water quality data that were collected over 25 years, at 102 catchments across the State of Victoria, Australia. We analysed and modelled spatial patterns of the 10th, 25th, 50th, 75th and 90th percentiles of the concentrations of sediments, nutrients and salt, with six common constituents: total suspended solids (TSS), total phosphorus (TP), filterable reactive phosphorus (FRP), total Kjeldahl nitrogen (TKN), nitrate-nitrite (NO_x), and electrical conductivity (EC). To predict the spatial variation of each quantile for each constituent, we developed statistical regression models and exhaustively searched through 50 catchment characteristics to identify the best set of predictors for that quantile. The models predict the spatial variation in individual quantiles of TSS, TKN and EC well (66%–96% spatial variation explained), while those for TP, FRP and NO_x have lower performance (37%–73% spatial variation explained). The most common factors that influence the spatial variations of the different constituents and quantiles are: annual temperature, percentage of cropping land area in catchment and channel slope. The statistical models developed can be used to predict how low- and high-concentration quantiles change with landscape characteristics, and thus provide a useful tool for catchment managers to inform planning and policy making with changing climate and land use conditions.     

Solution for a plane strain rough‐walled hydraulic fracture driven by turbulent fluid through impermeable rock

The impact of turbulent flow on plane strain fluid‐driven crack propagation is an important but still poorly understood consideration in hydraulic fracture modeling. The changes that hydraulic fracturing has experienced over the past decade, especially in the area of fracturing fluids, have played a major role in the transition of the typical fluid regime from laminar to turbulent flow. Motivated by the increasing preponderance of high‐rate, water‐driven hydraulic fractures with high Reynolds number, we present a semianalytical solution for the propagation of a plane strain hydraulic fracture driven by a turbulent fluid in an impermeable formation. The formulation uses a power law relationship between the Darcy‐Weisbach friction factor and the scale of the fracture roughness, where one specific manifestation of this generalized friction factor is the classical Gauckler‐Manning‐Strickler approximation for turbulent flow in a rough‐walled channel. Conservation of mass, elasticity, and crack propagation are also solved simultaneously. We obtain a semianalytical solution using an orthogonal polynomial series. An approximate closed‐form solution is enabled by a choice of orthogonal polynomials embedding the near‐tip asymptotic behavior and thus giving very rapid convergence; a precise solution is obtained with 2 terms of the series. By comparison with numerical simulations, we show that the transition region between the laminar and turbulent regimes can be relatively small so that full solutions can often be well approximated by either a fully laminar or fully turbulent solution.     

Modelling the poroelastoplastic behaviour of soils subjected to internal erosion by suffusion

Granular soils subjected to seepage flow may suffer suffusion, ie, a selective internal erosion. Extending the classical approach of poromechanics, we deduce a new form of the Clausius-Duhem inequality accounting for dissipation due to suffusion, and we deduce restrictions on the constitutive laws of the soil. We suggest (a) a possible coupling between the seepage forces and the suffusion kinetics and (b) an extension of an existing elastoplastic model for the skeleton mechanical behaviour. Numerical integrations of the elastoplastic model are carried out under drained axisymmetric triaxial and oedometric conditions. As a result, we prove that the extended model is able to qualitatively reproduce the suffusion induced strains and the strength reduction experimentally observed. Predictions on the oedometric behaviour of suffusive soils are also provided.     

黑龙江哈尔滨白渔泡国家湿地公园沼泽、林地和农田土壤物理、化学和生物性质差异

为了揭示黑龙江哈尔滨白渔泡国家湿地公园沼泽、林地和农田土壤物理、化学和生物性质的差异,于2018年7月25日~8月2日,在湿地公园内,在天然芦苇(Phragmites australis)沼泽、林地、旱田和水田中设置采样地,采集不同深度(0~10 cm、10~20 cm和20~30 cm)的土壤样品,测定土壤样品的物理、化学和生物指标。研究结果表明,白渔泡国家湿地公园不同采样地土壤指标存在差异;与天然芦苇沼泽土壤相比,其它采样地土壤的含水量明显偏低,土壤全氮、全磷、碱解氮和有机质含量都明显偏小,水田土壤速效磷含量偏大;天然芦苇沼泽土壤脲酶、硝酸还原酶、纤维素酶、蛋白酶和β-葡萄糖苷酶活性都高于林地和农田土壤,水田0~10 cm和10~20 cm深度土壤的硝酸还原酶活性显著高于旱田和林地;与天然芦苇沼泽土壤相比,旱田土壤小于0.25 mm的小团聚体含量偏大,而其它采样地土壤的各粒级团聚体的比例变化较小,水田土壤团聚体平均重量直径比天然芦苇沼泽和旱田土壤低。     

A modified numerical model for moisture-salt transport in unsaturated sandy soil under evaporation

Soil salinization, caused by salt migration and accumulation underneath the soil surface, will corrode structures. To analyze the moisture-salt migration and salt precipitation in soil under evaporation conditions, a mathematical model consisting of a series of theoretical equations is briefly presented. The filling effect of precipitated salts on tortuosity factor and evaporation rate are taken into account in relevant equations. Besides, a transition equation to link the solute transport equation before and after salt precipitation is proposed. Meanwhile, a new relative humidity equation deduced from Pitzer ions model is used to modify the vapor transport flux equation. The results show that the calculated values are in good agreement with the published experimental data, especially for the simulation of volume water content and evaporation rate of Toyoura sand, which confirm the reliability and applicability of the proposed model.     

自然和自然贡献情景模型研究综述

本文回顾了自然和自然贡献情景模型发展的背景、历史和内容,概括总结了自然和自然贡献情景模型的发展进程以及联合国生物多样性与生态系统服务政府间科学—政策平台(IPBES)情景模型的概念框架,讨论了自然和自然贡献情景模型存在的问题和发展方向。为了在全球层面解决现有综合集成模型存在的问题,根据地球表层建模基本定理和生态环境曲面建模基本定理,提出了具有中国原创特点的自然与自然贡献情景模型概念框架。     

Overview of In - Situ Biodegradation and Enhancement

Microbial degradation technologies have been developed to restore ground water quality in aquifers polluted by organic contaminants effectively in recent years. However, in course of the degradation, the formation of biofilms in ground water remediation technology can be detrimental to the effectiveness of a ground water remediation project. Several alternatives are available to a remedial design engineer, such as Permeable Reactive Barriers (PRBs) and in -situ bioremediation, Hydrogen Releasing Compounds (HRCs) barrier, Oxygen Releasing Compounds (ORCs) barrier etc. which are efficient and cost- effective technologies. Excessive biomass formation renders a barrier ineffective in degrading the contaminants, Efforts are made to develop kinetics models which accurately determine bio - fouling and bio - filn formation and to control excessive biomass formation.