Soil Hydraulic Properties Research Articles

New Method for Modelling Seasonal Variation in Resistance and Performance of Earthing Systems

The current techniques utilized for estimating seasonal fluctuations in earthing system resistance, including artificial neural networks (ANNs) and correlation/correction factors, rely on resistance records, soil resistivity measurements, and meteorological data collected across broad areas. However, they frequently fail to consider the impact of soil conditions and properties at the actual earthing location. As a solution, this research introduces a new method that models atmospheric conditions as soil suction and incorporates hydraulic soil properties (soil water retention characteristics and hydraulic conductivity) to estimate the seasonal changes in earthing resistance and performance. To illustrate this approach, this study constructs geometric models of vertical earthing rods for three homogeneous soil textures (clayey, silty, and sandy) utilizing COMSOL Multiphysics software. By coupling the differential equations governing electric current and water flow using Archie’s formula and solving numerically with the finite element method (FEM) for various soil suctions, this research reveals that soil water retention and resistivity variations are notably influenced by soil texture. Sandy soil displays higher variability, silt soil demonstrates moderate changes, while clayey soil exhibits lower fluctuation. By linking soil resistivity changes to soil suction and hydraulic properties, this innovative method predicts seasonal trends in earthing resistance and performance.

Soil water retention and hydraulic conductivity measured in a wide saturation range

Abstract. Soil hydraulic properties (SHPs), particularly soil water retention capacity and hydraulic conductivity of unsaturated soils, are among the key properties that determine the hydrological functioning of terrestrial systems. Some large collections of SHPs, such as the UNSODA and HYPRES databases, have already existed for more than 2 decades. They have provided an essential basis for many studies related to the critical zone. Today, sample-based SHPs can be determined in a wider saturation range and with higher resolution by combining some recently developed laboratory methods. We provide 572 high-quality SHP data sets from undisturbed, mostly central European samples covering a wide range of soil texture, bulk density and organic carbon content. A consistent and rigorous quality filtering ensures that only trustworthy data sets are included. The data collection contains (i) SHP data, which consist of soil water retention and hydraulic conductivity data, determined by the evaporation method and supplemented by retention data obtained by the dewpoint method and saturated conductivity measurements; (ii) basic soil data, which consist of particle size distribution determined by sedimentation analysis and wet sieving, bulk density and organic carbon content; and (iii) metadata, which include the coordinates of the sampling locations. In addition, for each data set, we provide soil hydraulic parameters for the widely used van Genuchten–Mualem model and for the more advanced Peters–Durner–Iden model. The data were originally collected to develop and test SHP models and associated pedotransfer functions. However, we expect that they will be very valuable for various other purposes such as simulation studies or correlation analyses of different soil properties to study their causal relationships. The data are available at https://doi.org/10.5880/fidgeo.2023.012 (Hohenbrink et al., 2023).

Study on soil hydraulic properties of slope farmlands with different degrees of erosion degradation in a typical black soil region.

In order to explore the impact of soil erosion degradation on soil hydraulic properties of slope farmland in a typical black soil region, typical black soils with three degrees of erosion degradation (light, moderate and heavy) were selected as the research objects. The saturated hydraulic conductivity, water holding capacity and water supply capacity of the soils were analyzed, as well as their correlations with soil physicochemical properties. The results showed that the saturated hydraulic conductivity of black soils in slope farmlands decreased with erosion degradation degree, which was higher in 0-10 cm soil layer than in 10-20 cm soil layer. The water holding capacity and water supplying capacity of typical black soils also decreased with the increase of erosion degradation degree, and both of them were stronger in the upper soil than in the lower soil. With the aggravation of erosion degradation of black soils, soil organic matter content decreased while soil bulk density increased, leading to the decline of soil hydraulic conductivity. The increase of soil bulk density and the decrease of contents of organic matter and >0.25 mm water stable aggregates were the main factors leading to the decrease of soil water holding capacity. These findings provide scientific basis and basic data for rational utilization of soil water, improvement of land productivity and prevention of soil erosion.

Changes in physical and hydraulic properties in sandy soils of the Pampa Biome under different uses

ABSTRACT The naturally fragile sandy soils of the Pampa Biome (PB) may be degraded with the introduction of poorly managed agricultural crops. Anthropic use can markedly decrease vegetation cover on sandy soils, leaving them more exposed to erosive agents. Decreases in organic matter content, biodiversity, and nutrient availability, increased soil compaction, and decreased water availability are also some of the impacts caused on PB soils by implementing poorly managed agricultural crops. In Rio Grande do Sul, in areas with sandier soils, the intense replacement of PB with commercial crops that has occurred in recent years (2000-2020) may be starting a disastrous cycle of degradation of these soils. However, it is not yet known how much these soils are degraded by recent anthropic use. There are no local scientific publications dedicated to the diagnosis of the degradation of these soils by anthropic use. Therefore, the need for this study was based on the objective of evaluating the physical and hydraulic properties of sandy PB soils with the insertion of agricultural crops. The study was conducted on three different soils, where soil samples were collected under three systems of use (extensive cattle raising on native grassland, eucalyptus afforestation, and soybean crop). Our results show that it can take more than nine years for soybean and ten years for eucalyptus land-use change to indicate some level of degradation in soil physical and hydraulic properties after replacing PB with these cultivated crops.

Development and functional evaluation of pedotransfer functions for estimating soil hydraulic properties in Portuguese soils: Implications for soil water dynamics

Accurate knowledge of soil hydraulic properties is essential for effectively simulating soil water dynamics in field plots. However, determining these properties directly can be challenging, especially when a large amount of data is required. Therefore, the primary objective of this study was to develop four hierarchical pedotransfer functions (PTFs) to estimate the soil hydraulic parameters of the Mualem-van Genuchten model. The PTFs that were developed ranged from utilizing a limited set of variables to those comprising a more comprehensive set of variables. These included particle size distribution, soil bulk density, soil depth, and soil water contents at matric potentials of −33 and − 1580 kPa. The functional performance of these PTFs was also assessed in simulations of the soil water balance in ten agricultural fields over two growing seasons using the HYDRUS-1D model. The simpler PTFs exhibited moderate accuracy in estimating soil water contents at different matric potentials, while the more complex PTFs with a greater number of variables yielded increasing accuracy. However, estimation errors were more noticeable in the drier range of the soil water retention curve. The inclusion of bulk density as a predictor only affected the wetter range of the same curve. Estimation of log Ks improved from the simpler to the more complex PTFs, but the statistical indicators were consistently modest. The functional performance of the PTFs showed increased errors in simulating soil water contents in the ten study sites compared to reference simulations. The more complex PTFs with the greatest number of variables provided the best predictions in most cases. The PTFs primarily affected flow conditions in the soil profile, resulting in increased soil water storage and impacting predictions of deep percolation the most. Water availability in the soil profile then influenced actual transpiration, alleviating stress in some cases but causing waterlogging in one situation. The study emphasizes the cautious use of the PTFs and the need for additional information on clay-rich soils and hydraulic conductivity data to enhance their performance.

A Novel Analytical Solution for Ponded Infiltration With Consideration of a Developing Saturated Zone

AbstractPonding at the soil surface exerts profound impacts on infiltration. However, the effects of ponding depth on infiltration, especially the development of a saturated zone below the soil surface, have yet to be considered in present infiltration models. A new general Green‐Ampt model solution (GAMS) was derived for a one‐dimensional vertical infiltration problem under a uniform initial moisture distribution with ponding on its surface. An expression was included in the new solution for simulating the saturated layer developed below the soil surface as long as the pressure head at the surface is sufficiently high to saturate the soil. The GAMS simulates the infiltration processes closer to the numerical solution by HYDRUS‐1D than the traditional and the recently improved Green‐Ampt model. Moreover, an inversion method to improve the estimates of soil hydraulic parameters from one‐dimensional vertical infiltration experiments that is based on the GAMS was suggested. The effect of ponding depth (hp), initial soil moisture content, soil texture, and hydraulic soil properties (saturated hydraulic conductivity Ks, water‐entry suction hd and shape coefficient n of soil water retention curve) in the saturated zone was also evaluated. The results indicate that the saturated zone length increased at a comparable rate with the unsaturated wetted zone length during infiltration. Generally, a larger saturated zone was found for soils with higher initial soil moisture contents, coarser texture, higher Ks values, greater n, and lower −hd. Our findings reveal that including the saturated zone in the infiltration model yields a better estimate of the soil hydraulic parameters. The proposed GAMS model can improve irrigation design and rainfall‐runoff simulations.

Divergent effects of food waste derived hydrochar on hydraulic properties and infiltration in a sandy soil

Hydrothermal carbonization of food waste (FW) and application of the solid product (hydrochar, HC) in agriculture, is expected to be an effective way to promote the circulation in food-energy-water nexus. However, little is known on the effects of applying FW-derived HC as a soil amendment on soil properties and functions. The current study investigated the effects of FW-derived HC on soil hydraulic properties and infiltration process. Two types of FW-derived HC were mixed into the top 0–20 cm soil layer at a rate of 1 % and 2 % (w/w) in a soil column experiment with grass growth. In-situ infiltration experiment was conducted within the soil column using a disc infiltrometer at different negative pressure (−20 cm, −6 cm, and −2 cm water head). Unsaturated hydraulic conductivity (K0) and sorptivity (S0) at each negative pressure head was estimated by fitting experimental cumulative infiltration data to a 4-terms infiltration model. Saturated hydraulic conductivity (Ks) and soil water retention curve (SWRC) were measured using undisturbed soil cores taken from the soil column. Soil pore size distribution and soil water constants were estimated from the SWRCs. The results indicate that application of FW-derived HC significantly (p < 0.05) affected soil structure by increasing meso (0.03–0.0002 mm) and macro (> 0.03 mm) pores, which can be attributed to the promoted soil aggregation. Saturated hydraulic conductivity was significantly increased and water retention characteristics were changed with significantly (p < 0.05) increased soil water holding capacity and plant available water content. Besides, FW-derived HCs were found hydrophobic and can render wettable sandy soil water repellent. The K0 and S0 were significantly (p < 0.05) decreased and cumulative infiltration was reduced when soil was wetted from an initial dry condition. Negative effect was found on grass growth with significantly (p < 0.05) decreased aboveground biomass after HCs application. Our study indicates that agriculture application of FW-derived HC as a soil amendment could effectively improve soil fertility by recycling nutrient elements and improving soil structure. However, attention should be paid to the hydrophobicity and phototoxicity of FW-derived HC to avoid the related detrimental influences on plant growth and environment.

Semi-analytical solution of pore-water pressure in unsaturated ground and infinite slope considering highly nonlinear soil hydraulic properties

A multi-exponential function (ME) is proposed to depict the highly nonlinear soil water retention curve and hydraulic conductivity function. Then, semi-analytical solution of pore-water pressure in an unsaturated infinite slope at steady state is derived using ME to depict soil hydraulic properties. The solution is verified with experimental and numerical results. It is found that ME can better describe the highly nonlinear hydraulic properties of soil than the single exponential function (SE), which is commonly-used previously. At steady state, the adoption of SE results in overestimating negative pore-water pressure (PWP; i.e., matric suction) by about 25 kPa and slope’s factor of safety (FOS) by 100% under rainfall, but underestimating negative PWP by about 50 kPa under evaporation. The difference between PWP and FOS calculated by using SE and ME becomes more significant for soil with bimodal hydraulic conductivity function than that with unimodal one. This difference reduces at a larger slope angle and deeper depth. Moreover, this difference generally increases at a lower SME/SSE, where SME and SSE represent the suction at which the unsaturated hydraulic conductivity expressed by ME and SE equals the applied rainfall intensity or evaporation rate, respectively. When SME/SSE is greater than 0.1, the difference basically diminishes.

The appraisal of pedotransfer functions with legacy data; an example from southern Africa

Predictions of soil hydraulic properties by pedotransfer functions (PTFs) must be treated with caution when they are used in an application domain which differs from the domain of their original development and calibration. However, in some settings, scientists may have little alternative but to use PTFs calibrated elsewhere. In this paper we consider how legacy data can be used to evaluate PTFs in new regions, paying particular attention to the challenges that arise when, as is often the case, the legacy data are not obtained by independent random sampling, and may be clustered at multiple scales. We undertook this work in southern Africa (Zimbabwe, Zambia and Malawi) where PTFs have been little-used, despite the scarcity of direct measurements of the soil properties of interest. We evaluated the extent to which existing PTFs provide a useful tool for the prediction of soil moisture content at field-capacity (−33 kPa) and permanent wilting-point (−1500 kPa) at different spatial scales. Soil legacy data for Zambia, Zimbabwe and Malawi were collated from various sources and PTFs from temperate and tropical domains were evaluated. We examined error variance components of predictions at within-profile, within-site and between-site scales; and estimated their mean errors. In general the better-performing PTFs (with respect to bias and the size of the error variance components) were ones calibrated with data from a tropical domain. This was most apparent at −1500 kPa. However, not all PTFs calibrated with data on tropical soils performed well, and predictions from some PTFs calibrated over a temperate domain were better at −33 kPa. The observations were spatially clustered, with data from different depth intervals in the same profile, from profiles in the same experimental site or farm, and from clusters across the region. This enabled us to show, with an appropriate mixed model analysis, that PTFs which effectively capture regional-scale variation may be less useful for predicting variation within a profile. We propose that such studies, based on legacy data, and with a suitable linear mixed model, should be used to screen PTFs of any provenance before their wider application.

CLSoilMaps: A national soil gridded database of physical and hydraulic soil properties for Chile

Spatially explicit soil information is crucial for comprehending and managing many of Earth´s processes related to carbon, water, and other biogeochemical cycles. We introduced a gridded database of soil physical properties and hydraulic parameters at 100 meters spatial resolution. It covers the continental area of Chile and binational basins shared with Argentina for six standardized depths following the specifications of the GlobalSoilMap project. We generated soil maps based on digital soil mapping techniques based on more than 4000 observations, including unpublished data from remote areas. These maps were used as input for the pedotransfer function Rosetta V3 to obtain predictions of soil hydraulic properties, such as field capacity, permanent wilting point, total available water capacity, and other parameters of the water retention curve. The trained models outperformed several other DSM studies applied at the national and regional scale for soil physical properties (nRMSE ranging from 6.93% to 15.7%) and delivered acceptable predictions (nRMSE ranging from 10.4% to 15.6%) for soil hydraulic properties, making them suitable for countless environmental studies.

Saturated Hydraulic Conductivity of a Sandy Loam under No-Till and Intensive Tillage in a Corn–Soybean Rotation

Tillage management practices have a dynamic impact on soil hydraulic properties and processes. There is a need for information about the effect of tillage practices on soil hydraulic properties for crops growing under sprinkler irrigation in the northern Great Plains. A long-term study was conducted from 2014 to 2018 to examine the effect of no tillage (NT) and conventional tillage (CT) on the saturated hydraulic conductivity (Ks) of a sandy loam soil in a two-year corn (Zea mays L.)–soybean (Glycine max L.) rotation. In situ Ks measurements were taken in the center of crop rows within NT and CT plots using a pressure ring infiltrometer at the soil surface (0–15 cm) and a constant head well permeameter at the subsurface (15–30 cm). Results indicated that Ks values were well described by a log-normal distribution at both depths. Results showed that logarithmic Ks (log Ks) was not significantly impacted by tillage. Averaged over the five-year study, the log-transformed Ks of 100 measurements was not significantly affected by tillage in the surface layer under either corn or soybean nor in the subsurface layer under soybean. However, the mean soil log Ks in CT plots (1.784 mm h−1) was significantly greater than that in NT plots (1.186 mm h−1) in the 15–30 cm layer under corn, while Ks was nearly 50% greater in CT than in NT. Large values for the coefficient of variation (CV%) of Ks measurements exhibited significant spatial variations of Ks among plots within each tillage treatment at both the soil surface and subsurface layers under corn and soybean. Thus, more studies under different soils and cropping systems with a larger sample size per treatment are needed to lower spatial variability within treatments and validate the effect of tillage on soil hydraulic properties.

Integrative effect of activated biochar to reduce water stress impact and enhance antioxidant capacity in crops

Biochar is a soil amendment that can change soil's physical and hydraulic properties. However, biochar application is far from being a ‘one size fits-all’ approach. The impact of the management practices is dependent on biochar type (feedstock and production conditions), application depth and method, climate and site characteristics. Hence, this study aims to enrich the available inconclusive information on how biochar could affect clay loamy soil and to assess the potential impact of the induced change on water stress mitigation of rain-fed durum wheat under the specific condition of the semi-arid environment of North West of Tunisia. A field experiment was investigated in which three biochar rates 0 (B0), 10 (equivalent to 0.5% of weight) (B1) and 20 t/ha (equivalent to 1% of weight), (B2), were tested. Other laboratory analysis allowed the evaluation of soil water retention curve (SWRC), saturated hydraulic conductivity (Ks), dry density (ρb) and biostress biomarkers such as glutathione-S-transferase (GST), catalase activities (CAT) and malondialdehyde content (MDA) as well as yield attributes. Results showed that treatment B2 significantly decreased ρb and Ks with relative change values of about −3.1% and −19%. Consequently, SWRC showed a better water retention capacity, mostly from saturation to matric potential value (h) of 33 kPa. Total (TAWC), plant (PAWC) and readily (RAWC) available water contents, significantly increased under B2 with relative changes of +6%, +44% and +44% respectively. Moreover, GST and CAT were also boosted under B2. Consequently, biological and grain yields as well as grain water use efficiency (GWUE) significantly increased. GWUE increased from 0.81 ± 0.04 in B0 to 1.09 ± 0.01 kg/m3 in B2. The correlation analysis showed a significant and positive correlation, between GWUE and soil water parameters (θs, θfc and θmre) suggesting the indirect effect of biochar on water-use efficiency for grain yield of wheat. Therefore, among the tested rates 20 t/ha could be suggested to improve plant soil water availability and reduce the harmful impact of drought stress on rain-fed durum wheat.

Impacts of Land Use Changes on Soil Functions and Water Security: Insights from a Three-Year-Long Study in the Cantareira System, Southeast of Brazil

Maintaining soil functions is crucial for human well-being, but there is a lack of integration between soil, water security, ecosystem services, and climate change. To bridge this knowledge gap and address erosion-induced soil and water losses and considering intrinsic impacts of soil structure, a three-year-long study was conducted focused on three dominant soil types (Typic Hapludult, Typic Dystrudept, and Typic Usthortent) combined with different land uses (native forest, eucalyptus plantation, rotational grazing, and extensive grazing) in a critical water supply region for the São Paulo metropolitan area in Southeastern Brazil. Surface runoff, evaluated for erosion resistance, was measured using the Cornell infiltrometer, and soil electrical resistivity tomography estimated soil water content to a depth of 1.5 m for groundwater recharge analysis. Soil hydraulic properties were also measured. The results revealed that native forest soils had higher hydraulic conductivity, particularly in the surface layer, compared to eucalyptus and pastures. Native forests in Typic Hapludult showed a higher runoff rate (200 to 250 mm h−1) due to a naturally dense subsoil layer that negatively impacted water infiltration and recharge with a high erosion potential, therefore reducing the amount of water stored. Typic Usthortent maintained a higher soil water content in pastures than in other land uses and also showed a low rate of water infiltration, resulting in perched water in the surface layer. In Typic Dystrudept, the native forest presented higher hydraulic conductivity (0–5 cm: 115.9 cm h−1) than eucalyptus (0–5 cm: 36.4 cm h−1), rotational grazing (0–5 cm: 19.4 cm h−1), and extensive grazing (0–5 cm: 2.6 cm h−1), but there were no significant differences in soil water content among land uses. This work illustrates the crucial role of native forests in affecting deep water recharge, reducing the soil surface erosion, mainly in soils without naturally subsoil layer, maintaining recharge potential. For Ultisols, pastures preserved soil structure and are therefore less impacted by soil management. With these results, a contribution is made to soil and water conservation, providing support for sustainable management practices in erosion-prone areas.

Analyzing the role of soil and vegetation spatial variability in modelling hydrological processes for irrigation optimization at large scale

The main purpose of this paper was to study the effect of spatial variability of soil hydraulic properties (HP) and vegetation parameters (VP) (e.g., leaf-area index, LAI, and crop coefficient, Kc) on modelling agro-hydrological processes and optimising irrigation volumes at large scale. Based on this analysis, the effect of partly overlooking the spatial variability of soil HP and/or VP inputs was verified on a 140 ha irrigation sector in “Sinistra Ofanto” irrigation system in Apulia Region, Southern Italy. Five soil profiles were excavated and the HP were measured in all the soil horizons. Additionally, measurements of soil HP were taken in the surface soil layer in ninety sites distributed over the whole irrigation sector. All the HP measurements were carried out using tension infiltrometer. Remote sensing applications were used to obtain LAI and Kc using European Space Agency’s (ESA) Sentinel-2 images with 10 m resolutions. First, distributed (on ninety polygons with an average area of about 1.5 ha) optimal irrigation volumes and related deep percolation volumes at a depth of 80 cm, were computed using an agro-hydrological model and accounting for the actual observed variability of soil HP and VP inputs. The sector scale irrigation and deep percolation volumes were obtained by aggregating the distributed irrigation volumes. This was considered as the reference scenario (hereafter DVS—Detailed Variability Scenario). Then, reduced variability scenarios (hereafter RVS—Reduced Variability Scenario) were considered, where the information on the actual spatial variability of the soil HP and VP was gradually overlooked to find the minimum data set needed to still have sector scale irrigation volumes and related deep percolation volumes comparable to those obtained under the DVS. Results showed that overlooking VP (RVS-VP) variability did not significantly change the optimal irrigation volumes and the deep percolation fluxes. By contrast, neglecting the HP variability (RVS-HP) showed significant effects on both the irrigation and percolation volumes compared to the DVS. The main practical finding was that, at least for the area investigated in this study, hydraulic characterization of one soil profile in an area of approximately 30 ha provides sector scale irrigation volumes and percolation fluxes comparable to those obtained under the DVS, thus by accounting for all the observed local variability.

Predicting soil hydraulic properties for binary mixtures – concept and application for constructed Technosols

Abstract. Constructed Technosols are an important means of substituting natural soil material, such as peat and geogenic material, for use in urban green infrastructure. One characteristic of Technosols important to their role in urban green infrastructure, specifically with respect to urban water management, is their soil hydraulic properties (SHPs). The SHPs depend on the composition of the constructed Technosols (e.g. their components and their mixing ratio). The diversity of possible components and the infinite number of mixing ratios practically prohibit the experimental identification of the composition needed to achieve suitable soil hydrological functions. In this study, we propose a compositional model for predicting the water retention curves (WRCs) of any binary mixture based on the measured WRCs of its two pure components only (basic scheme) or with one additional mixture (extended scheme). The unsaturated hydraulic conductivity curves (HCCs) are predicted based on the modelled WRCs. The compositional model is developed from existing methods for estimating the porosity of binary mixtures. The model was tested on four data sets of measured WRCs of different binary mixtures. The distribution of water and air in 50 cm high soil columns filled with these mixtures was predicted under hydrostatic conditions in order to assess their suitability for typical urban applications. The difference between the maxima of the pore size distributions ΔPSDmax (m) of the components indicates the applicability of the compositional approach. For binary mixtures with small ΔPSDmax, the water content deviations between the predicted and the measured WRCs range from 0.004 to 0.039 cm3 cm−3. For mixtures with a large ΔPSDmax, the compositional model is not applicable. The prediction of the soil hydraulic properties of any mixing ratio facilitates the simulation of flow and transport processes in constructed Technosols before they are produced (e.g. for specific urban water management purposes).

A modified approach to model the hydrological effects of living moss and fallen leaves on carbonate-derived laterite in karst forests

Due to recent vegetation restoration projects in karst mountainous areas of China, the living moss and litter layer play important roles in controlling soil water dynamics in karst forests. However, the underlying mechanisms of the effects of moss and the litter layer on soil water dynamics are still controversial. To explore the underlying mechanisms, we employ the CoreKarHydro model, capturing soil water content dynamics driven by evapotranspiration and influenced by living moss and fallen leaves. This model integrates surface changes' impact on evapotranspiration, categorizing the domain into: no cover (S), moss (SM), pine needles (SP), and both (SMP). Model validation employs the MOCOM-UA method, using diverse metrics (RMSE, R2, ACC, ED, KGE and SA) against lysimeter data with 4 moss biomass levels and 3 pine needle biomass levels. The results indicated: 1) The CoreKarHydro model exhibited reliable performance in simulating the soil water content curves, with an RMSE ranging from 1.49% to 4.21% and an R2 value ranging from 0.60 to 0.77. 2) Among the factors that influenced evaporation processes, soil hydraulic properties emerged as the most sensitive parameters, whereas the properties of the moss and pine needles played a relatively lesser role. 3) The impact of moss cover on evaporation processes varied depending on the soil texture and the stage of evaporation. Moss cover reduced evaporation loss in 8 out of 12 soil textures (loamy sand, loam, silt, silt loam, sandy clay loam, silty clay loam, sandy clay, and silty clay). On the other hand, moss cover increased evaporation loss in 4 out of 12 soil textures (sand, sandy loam, clay loam, and clay). This study contributes to the understanding of the mechanisms of coupled soil-moss(litter)-hydrological processes in karst areas.

Microplastics affect C, N, and P cycling in natural environments: Highlighting the driver of soil hydraulic properties

As microplastics (MPs) are organic polymers with a carbon-based framework, they may affect nutrient cycling. Information regarding how MPs influence N, P, and C cycling and the underlying driving force remains lacking. N, P, and C cycling induced by soil hydraulic properties under MPs exposure (including polyethylene (PE), polyvinyl chloride (PVC), polystyrene (PS), polypropylene (PP)) in the natural environment were investigated in this study. MPs exposure increased the soil water content (11.2–84.5%) and reduced bulk density (11.4–42.8%); soil saturated hydraulic conductivity increased by 7.3–69.4% under PP and PE exposure. MPs exposure led to increases in available phosphorus, NO3--N, NH4+-N, and soil organic matter; the bacterial communities related to N and C cycling were significantly changed. Expression levels of soil N and C cycling-related genes were enhanced under low concentrations (0.5% and 2%) of MPs, except PVC; consequently, soil nitrogen storage and organic carbon storage increased by 12–75% and 6.7–93%, respectively. Correlation analyses among soil hydraulic properties, bacterial communities, and functional genes related to nutrient cycling revealed that soil hydraulic properties (including soil water content, saturated water capacity, and soil saturated hydraulic conductivity) were the dominant factors affecting soil N and C storage under MPs exposure.

Effect of fines loss on the microstructure and shear behaviors of gap-graded soils: A multiscale perspective

Natural sand is typically composed of both coarse and fine particles and are susceptible to suffusion erosion under seepage effect. Fines loss caused by suffusion is a critical factor affecting the mechanical and hydraulic properties of soils. This paper investigates the effect of fines loss on the microstructure and shear behaviors of gap-graded soils using a three-dimensional finite element and discrete element coupling method. A modified particle removal scheme based on the internal moment tensor is first presented, which is utilized to generate several specimens with varying degrees of fines loss. The macroscopic shear responses of both eroded and non-eroded specimens are benchmarked using the responses of representative volume elements (RVEs) in pure discrete element method (DEM). The results reveal that the eroded specimen exhibits a reduced peak stress ratio yet almost unchanged critical state values compared to the non-eroded specimen. Fines loss restrains the formation of shear band and results in more uniform stress and strain distribution inside the specimens. In addition, fine particles in gap-graded soils only sustain a small proportion of external loads, particularly for the specimens with more fines loss. The RVEs at different locations of the specimens are found to exhibit distinct packing assemblies and force networks at the end of shearing.

Gas permeability and emission in unsaturated vegetated landfill cover with biochar addition

Plant–biochar interaction has been recognized to affect the hydraulic properties of landfill cover soils, while its influence on landfill gas emission is rarely studied. This study investigated the coupled effects of biochar and vegetation on gas permeability and emission in unsaturated landfill cover through an integrated theoretical modelling and laboratory investigation. First, a gas permeability model was developed for vegetated coarse-grained soils with biochar addition. Then, a well-instrumented laboratory column test and two tests from the literature, considering bare, grass, biochar and grass + biochar conditions, were used for model validation. Finally, a numerical parametric study was conducted to investigate the influence of root growth and drought conditions on the gas emission rate. Results showed that the developed model can satisfactorily capture the gas permeability of unsaturated soils at various degrees of saturation. The lowest water retention capacity, the highest gas permeability and gas emission rate after 24 months of growth were observed in the grassed column. However, adding biochar in vegetated soils can maximize the water retention capacity and decrease the gas permeability, resulting in the lowest gas emission rate. The measured gas emission rates for the four cases meet the recommended value by the design guideline. The parametric study showed that the increased root depth from 0.2 m to 0.4 m improved the gas emission rate by 170% in the grass case but decreased by 97% in the grass + biochar case. Under the severe drought condition with soil suction around 500 kPa, the gas emission rate in the grassed case exceeded the design value by 18%, while those in the biochar cases were far below the allowable value. Therefore, peanut shell biochar should be considered to amend the grassed landfill cover using coarse-grained soils as it can significantly improve engineering performance in reducing gas emissions under extreme drought conditions.Graphical abstract

Reducing uncertainties in hydromechanical modeling with a recently developed Rosetta 3 podeotransfer function

Stability analysis of unsaturated landslide deposits requires reliable estimates of soil moisture and pore water pressure. However, modeled soil moisture and pore water pressure contain substantial uncertainties due to imperfect information on soil hydraulic properties. Due to the relatively high dimensionality, commonly used parameter optimization strategies can be significantly affected by equifinality problems. This study investigates the effectiveness of reducing parameter estimation dimensionality using soil pedo-transfer functions. Specifically, we first estimated soil hydraulic parameters using the traditional Generalized Likelihood Uncertainty Estimation (GLUE) method, with parameters randomly drawn from the entire space (refer to as GLUE-random). In a second strategy, we use the Rosetta 3 pedotransfer function to constrain soil hydraulic parameters (refer to as GLUE-Rosetta). The two methods were tested in a typical landslide deposit with in-situ measured soil moisture dynamics for inverse modeling. The GLUE-random estimated soil hydraulic parameters contained substantial uncertainties –resulting in poorly constrained soil water retention curves (SWCC) and hydraulic conductivity functions (HCF). As a result, the uncertainty bands of pore water pressure and slope stability can cross values with several orders of magnitudes. In contrast, GLUE-Rosetta provided well-constrained SWCC and HCF, which significantly reduce the uncertainties in pore water pressure and slope stability estimates. These results suggest that the Rosetta 3 pedotransfer function can significantly improve the reliability of soil hydraulic parameters by reducing the dimensionality of the optimization problem and high-quality prior information of soil hydraulic properties. In conclusion, Rosetta 3 can enhance the reliability of soil parameters estimates and the reliability of subsurface hydrology, which may benefit the development of landslide early-warning systems.