Response of soil aggregate disintegration to the different content of organic carbon and its fractions during splash erosion

Aggregate disintegration is a critical process in soil splash erosion. However, the effect of soil organic carbon (SOC) and its fractions on soil aggregates disintegration is still not clear. In this study, five soils with similar clay contents and different contents of SOC have been used. The effects of slaking and mechanical striking on splash erosion were distinguished by using deionized water and 95% ethanol as raindrops. The simulated rainfall experiments were carried out in four heights (0.5, 1.0, 1.5 and 2.0 m). The result indicated that the soil aggregate stability increased with the increases of SOC and light fraction organic carbon (LFOC). The relative slaking and the mechanical striking index increased with the decreases of SOC and LFOC. The reduction of macroaggregates in eroded soil gradually decreased with the increase of SOC and LFOC, especially in alcohol test. The amount of macroaggregates (>0.25 mm) in deionized water tests were significantly less than that in alcohol tests under the same rainfall heights. The contribution of slaking to splash erosion increased with the decrease of heavy fractions organic carbon. The contribution of mechanical striking was dominant when the rainfall kinetic energy increased to a range of threshold between 9 J m⁻² mm⁻¹ and 12 m⁻² mm⁻¹. This study could provide the scientific basis for deeply understanding the mechanism of soil aggregates disintegration and splash erosion.     

A plasticity constitutive model for unsaturated,anisotropic, nonexpansive soils

The paper describes and evaluates an incremental plasticity constitutive model for unsaturated, anisotropic, nonexpansive soils (CMUA). It is based on the modified Cam-Clay (MCC) model for saturated soils and enhances it by introducing anisotropy (via rotation of the MCC yield surface) and an unsaturated compressibility framework describing a double dependence of compressibility on suction and on the degree of saturation of macroporosity. As the anisotropic and unsaturated features can be activated independently, the model is downwards compatible with the MCC model. The CMUA model can simulate effectively: the dependence of compressibility on the level of developed anisotropy, uniqueness of critical state independent of the initial anisotropy, an evolving compressibility during constant suction compression, and a maximum of collapse. The model uses Bishop's average skeleton stress as its first constitutive variable, favouring its numerical implementation in commercial numerical analysis codes (eg, finite element codes) and a unified treatment of saturated and unsaturated material states.     

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.     

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.     

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.     

Enhancement of backtracking search algorithm for identifying soil parameters

In this paper, an enhanced backtracking search algorithm (so-called MBSA-LS) for parameter identification is proposed with two modifications: (a) modifying the mutation of original backtracking search algorithm (BSA) considering the contribution of current best individual for accelerating convergence speed and (b) novelly incorporating an efficient differential evolution (DE) as local search for improving the quality of population. The proposed MBSA-LS is first validated with better performance than the original BSA and some other typical state-of-the-art optimization algorithms on a benchmark of soil parameter identification in terms of effectiveness, efficiency, and robustness. Then, the efficiency of the MBSA-LS is further illustrated by two representative cases: identifying soil parameters from both laboratory tests and field measurements. All comparisons demonstrate that the proposed MBSA-LS algorithm can give accurate results in a short time. Finally, to conveniently solve the problems of parameter identification, a practical tool ErosOpt for parameter identification is developed by integrating the proposed MBSA-LS and some other efficient algorithms for readers to conduct the parameter identification using optimisation algorithms.     

A smoothed cap model for variably saturated soils and its robust numerical implementation

This paper deals with the numerical implementation of a cap model for unsaturated soils. It provides a brief review of existing cap model approaches, based on which an improved model formulated in terms of generalised effective stress and matric suction is derived and described in detail. Although the proposed model is a multisurface plasticity model, it can efficiently be implemented using only single‐surface projections because of the smoothness of the model, which is obtained by construction. Numerical algorithms are provided for these single‐surface stress projections, using a single‐equation approach whenever possible. The robustness of the utilised single‐equation approaches is enhanced by proposing problem‐fitted start‐up procedures based on investigations of the nonlinear projection equations. A comparison of the model response with extensive material test data is used to validate the model and to demonstrate the robust application of the approach to silty sands and low to medium plasticity clays. Copyright © 2015 John Wiley & Sons, Ltd.     

Generation of continuous surface soil moisture dataset using combined optical and thermal infrared images

Surface soil moisture (SSM) is a critical variable for understanding water and energy flux between the atmosphere and the Earth's surface. An easy to apply algorithm for deriving SSM time series that primarily uses temporal parameters derived from simulated and in situ datasets has recently been reported. This algorithm must be assessed for different biophysical and atmospheric conditions by using actual geostationary satellite images. In this study, two currently available coarse‐scale SSM datasets (microwave and reanalysis product) and aggregated in situ SSM measurements were implemented to calibrate the time‐invariable coefficients of the SSM retrieval algorithm for conditions in which conventional observations are rare. These coefficients were subsequently used to obtain SSM time series directly from Meteosat Second Generation (MSG) images over the study area of a well‐organized soil moisture network named REMEDHUS in Spain. The results show a high degree of consistency between the estimated and actual SSM time series values when using the three SSM dataset‐calibrated time‐invariable coefficients to retrieve SSM, with coefficients of determination (R²) varying from 0.304 to 0.534 and root mean square errors ranging from 0.020 m³/m³ to 0.029 m³/m³. Further evaluation with different land use types results in acceptable debiased root mean square errors between 0.021 m³/m³ and 0.048 m³/m³ when comparing the estimated MSG pixel‐scale SSM with in situ measurements. These results indicate that the investigated method is practical for deriving time‐invariable coefficients when using publicly accessed coarse‐scale SSM datasets, which is beneficial for generating continuous SSM dataset at the MSG pixel scale.     

Temporal stability of soil moisture in an oasis of northwestern China

Temporal stability of soil moisture has been widely used in hydrological monitoring since it emerged. However, the spatial analysis of temporal stability at the landscape scale is often limited because of insufficient sampling numbers. This work made an effort to investigate the spatial variations of temporal stability of soil moisture in an oasis landscape. The specific objectives of the study were to explore the spatial patterns of temporal stability and to determine the controlling factors of temporal stability in the desert oasis. A time series of soil moisture measurements were gathered on 23 occasions at 118 locations over 3 years in a rectangular transect of approximately 100 km². The nonparametric Spearman's rank correlation coefficient, standard deviation of relative difference (SDRD), and mean absolute bias error (MABE) were used to quantify the temporal stability of soil moisture. Results showed that the temporal stability of soil moisture was depth dependent and season dependent. The spatial pattern of soil moisture in a deep soil layer and between two same seasons generally had a high temporal stability. SDRD and MABE were spatially autocorrelated and exhibited strong spatial structures in the geographic space. The concept of temporal stability can be extended to describe the time‐stable areas of soil moisture with geostatistics. There were great differences between SDRD and MABE in describing the temporal stability of soil moisture and in identifying the controlling factors of temporal stability. In this case, MABE was a better alternative to estimate the areal mean soil moisture using representative locations than SDRD. Land use type, soil moisture condition, and soil particle composition were the dominant controls of temporal stability in the oasis. These insights could help to better understand the essence of temporal stability of soil moisture in arid regions.     

Groundwater recharge in an arid grassland as indicated by soil chloride profile and multiple tracers

Previous studies have shown that shallow groundwater in arid regions is often not in equilibrium with near‐surface boundary conditions due to human activities and climate change. This is especially the case where the unsaturated zone is thick and recharge rate is limited. Under this nonequilibrium condition, the unsaturated zone solute profile plays an important role in estimating recent diffuse recharge in arid environments. This paper combines evaluation of the thick unsaturated zone with the saturated zone to investigate the groundwater recharge of a grassland in the arid western Ordos Basin, NW China, using the soil chloride profiles and multiple tracers (²H, ¹⁸O, ¹³C, ¹⁴C, and water chemistry) of groundwater. Whereas conventional water balance and Darcy flux measurements usually involve large errors in recharge estimations for arid areas, chloride mass balance has been widely and generally successfully used. The results show that the present diffuse recharge beneath the grassland is 0.11–0.32 mm/year, based on the chloride mass balance of seven soil profiles. The chloride accumulation age is approximately 2,500 years at a depth of 13 m in the unsaturated zone. The average Cl content in soil moisture in the upper 13 m of the unsaturated zone ranges from 2,842 to 7,856 mg/L, whereas the shallow groundwater Cl content ranges from 95 to 351 mg/L. The corrected ¹⁴C age of shallow groundwater ranges from 4,327 to 29,708 years. Stable isotopes show that the shallow groundwater is unrelated to modern precipitation. The shallow groundwater was recharged during the cold and wet phases of the Late Pleistocene and Holocene humid phase based on palaeoclimate, and consequently, the groundwater resources are nonrenewable. Due to the limited recharge rate and thick unsaturated zone, the present shallow groundwater has not been in hydraulic equilibrium with near‐surface boundary conditions in the past 2,500 years.