The influence of water retention curve hysteresis on the stability of unsaturated soil slopes

Hysteresis is a common feature exhibited in hydraulic properties of an unsaturated soil. The movement of wetting front and the hysteresis effect are important factors which impact the shear strength of the unsaturated soil and the mechanics of shallow landslides. These failures are mainly triggered by the deepening of the wetting front accompanied by a decrease in matric suction induced by infiltration. This research establishes a method for determining a stability analysis of unsaturated infinite soil slopes, integrating the influence of infiltration and the water retention curve hysteresis. Furthermore, the present stability analysis method including the infiltration model and the advanced Mohr–Coulomb failure criterion calculates the variations of the safety factor (FS) in accordance with different slope angle, depth and hydrological processes. The experimentally measured data on the effect of hysteresis are also carried out for comparison. Numerical analyses, employing both wetting and drying hydraulic behaviour of unsaturated soil, are performed to study the difference in soil‐water content as observed in the experiments. The simulating approximations also fully responded to the experimental data of sand box. The results suggest that the hysteresis behaviour affect the distribution of soil‐water content within the slope indeed. The hysteresis made the FS values a remarkable recovery during the period of non‐rainfall in a rainfall event. The appropriate hydraulic properties of soil (i.e. wetting or drying) should be used in accordance with the processes that unsaturated soil actually experience. This method will enable us to acquire more accurate matric suction head and the unsaturated soil‐shear strength as it changes with the hysteretic flow, in order to calculate into the stability analysis of shallow landslides. An advanced understanding of the process mechanism afforded by this method is critical to realizing a reliable and appropriate design for slope stabilization. It also offers some immediate reference information to the disaster reduction department of the government. Copyright © 2011 John Wiley & Sons, Ltd.     

Influence of heterogeneity on unsaturated hydraulic properties (2) – percentage and shape of heterogeneity

In subsurface porous media, the soil water retention curve (WRC) and unsaturated hydraulic conductivity curve (UHC) are two important soil hydraulic property curves. Spatial heterogeneity is ubiquitous in nature, which may significantly affect soil hydraulic property curves. The main theme of this paper is to investigate how spatial heterogeneities, including their arrangements and amounts in soil flumes, affect soil hydraulic property curves. This paper uses a two‐dimensional variably saturated flow and solute transport finite element model to simulate variations of pressure and moisture content in soil flumes under a constant head boundary condition. To investigate the behavior of soil hydraulic property curves owing to variations of heterogeneities and their arrangements as well, cases with different proportions of heterogeneities are carried out. A quantitative evaluation of parameter variations in the van Genuchten model (VG model) resulting from heterogeneity is presented. Results show that the soil hydraulic properties are strongly affected by variations of heterogeneities and their arrangements. If the pressure head remains at a specific value, the soil moisture increases when heterogeneities increase in the soil flumes. On the other hand, the unsaturated hydraulic conductivity decreases when heterogeneities increase in the soil flumes under a constant pressure head. Moreover, results reveal that parameters estimated from both WRC and UHC also are affected by shapes of heterogeneity; this indicates that the parameters obtained from the WRC are not suitable for predicting the UHC of different shapes in heterogeneous media. Copyright © 2011 John Wiley & Sons, Ltd.     

Inverse upscaling of hydraulic parameters during constant flux infiltration using borehole radar

Modeling unsaturated flow in porous media requires constitutive relations that describe the soil water retention and soil hydraulic conductivity as a function of either potential or water content. Often, the hydraulic parameters that describe these relations are directly measured on small soil cores, and many cores are needed to upscale to the entire heterogeneous flow field. An alternative to the forward upscaling method using small samples are inverse upscaling methods that incorporate soft data from geophysical measurements observed directly on the larger flow field. In this paper, we demonstrate that the hydraulic parameters can be obtained from cross borehole ground penetrating radar by measuring the first arrival travel time of electromagnetic waves (represented by raypaths) from stationary antennae during a constant flux infiltration experiment. The formulation and coupling of the hydrological and geophysical models rely on a constant velocity wetting front that causes critical refraction at the edge of the front as it passes by the antennae. During this critical refraction period, the slope of the first arrival data can be used to calculate (1) the wetting velocity and (2) the hydraulic conductivity of the wet (or saturated) soil. If the soil is undersaturated during infiltration, then an estimate of the saturated water content is needed before calculating the saturated hydraulic conductivity. The hydraulic conductivity value is then used in a nonlinear global optimization scheme to estimate the remaining two parameters of a Broadbridge and White soil.     

Influence of heterogeneity on unsaturated hydraulic properties: (1) local heterogeneity and scale effect

Spatial heterogeneity is ubiquitous in nature, which may significantly affect the soil hydraulic property curves. The models of a closed‐form functional relationship of soil hydraulic property curves (e.g. VG model or exponential model) are valid at point or local scale based on a point‐scale hydrological process, but how do scale effects of heterogeneity have an influence on the parameters of these models when the models are used in a larger scale process? This paper uses a two‐dimensional variably saturated flow and solute transport finite element model (VSAFT2) to simulate variations of pressure and moisture content in the soil flume under a constant head boundary condition. By changing different numerical simulation block sizes, a quantitative evaluation of parameter variations in the VG model, resulting from the scale effects, is presented. Results show that the parameters of soil hydraulic properties are independent of scale in homogeneous media. Parameters of α and n in homogeneous media, which are estimated by using the unsaturated hydraulic conductivity curve (UHC) or the soil water retention curve (WRC), are identical. Variations of local heterogeneities strongly affect the soil hydraulic properties, and the scale affects the results of the parameter estimations when numerical experiments are conducted. Furthermore, the discrepancy of each curve becomes considerable when moisture content becomes closer to a dry situation. Parameters estimated by UHC are totally different from the ones estimated by WRC. Copyright © 2012 John Wiley & Sons, Ltd.     

Hydraulic Properties of Stony Vadose Zonesa

Hydraulic properties of saturated and unsaturated stony soils were studied on a 3.35 m long column, 1.24 m in diameter, filled with alternating sand and boulder layers. The boulders averaged 6.2 × 15 × 20 cm in size and were laid down on their flat side. Tensiometers and a neutron probe access tube were placed in the column for measuring pressure heads and water contents, respectively. Saturated conditions were obtained by ponding the column. The resulting hydraulic conductivity K was 5.1 m/day. This value could also be calculated from the measured K for the sand alone on separate samples, using a simple equation that takes into account the void ratio of the sand alone and that of the boulder-sand mixture. Unsaturated K was determined by applying water at less than ponded infiltration rates. Resulting relations between the unsaturated K and water content or negative pressure head could also be estimated from the relation between unsaturated K and pressure head for the sand alone and the calculated saturated K of the boulder-sand mixture. The method of Millington and Quirk for calculating the relation between unsaturated K and water content also gave reasonable results. The dispersivity of the boulder-sand column was 18 times that of the sand alone. Pore velocity was accurately estimated as the Darcy velocity divided by the volumetric water content. Hydraulic properties of stony vadose zones are difficult to determine. This work shows that they can be estimated from K relations measured in the laboratory on samples of the soil between the rocks. Knowledge of hydraulic properties of vadose zones is important in predicting movement of water and pollutants to the underlying ground water.     

A physical–chemical model for the static water retention characteristic of unsaturated porous media

Water retention characteristic or water retention curve (WRC) is an important constitutive feature of porous media, and also meanwhile is an indispensable requirement in hydraulic transport modelling. Previous experiments have indicated that the specific surface area of porous media has effects on the WRC. It has also been observed that a linear relationship generally exists between the air–water interface area and the water saturation within unsaturated porous media. However it seems that no detailed study on their internal linkage with the WRC has been reported yet. This paper, at first gives a review of the development of WRC modelling, then it tries to explain the water retention characteristic according to the physical and chemical behaviours of the phases involved in unsaturated porous media. Using the traditional capillary theory, the volume averaging theorem and the advances in physical chemistry of interfacial surfaces, this paper then derives out a formula which represents the water retention characteristic of porous media. This formula demonstrates the internal linkage of the WRC to the specific surface area of porosities. It also shows agreements with the experimental observations mathematically. Based on this formula, a fitting model is proposed for the static WRC of porous media. Finally, this model is tested to fit the WRC measurements of a wide range of porous materials. Comparison with other main models is presented.     

Impact of land use on the hydraulic properties of the topsoil in a small French catchment

The hydraulic properties of the topsoil control the partition of rainfall into infiltration and runoff at the soil surface. They must be characterized for distributed hydrological modelling. This study presents the results of a field campaign documenting topsoil hydraulic properties in a small French suburban catchment (7 km²) located near Lyon, France. Two types of infiltration tests were performed: single ring infiltration tests under positive head and tension‐disk infiltration using a mini‐disk. Both categories were processed using the BEST—Beerkan Estimation of Soil Transfer parameters—method to derive parameters describing the retention and hydraulic conductivity curves. Dry bulk density and particle size data were also sampled. Almost all the topsoils were found to belong to the sandy loam soil class. No significant differences in hydraulic properties were found in terms of pedologic units, but the results showed a high impact of land use on these properties. The lowest dry bulk density values were obtained in forested soils with the highest organic matter content. Permanent pasture soils showed intermediate values, whereas the highest values were encountered in cultivated lands. For saturated hydraulic conductivity, the highest values were found in broad‐leaved forests and small woods. The complementary use of tension‐disk and positive head infiltration tests highlighted a sharp increase of hydraulic conductivity between near saturation and saturated conditions, attributed to macroporosity effect. The ratio of median saturated hydraulic conductivity to median hydraulic conductivity at a pressure of − 20 mm of water was about 50. The study suggests that soil texture, such as used in most pedo‐transfer functions, might not be sufficient to properly map the variability of soil hydraulic properties. Land use information should be considered in the parameterizations of topsoil within hydrological models to better represent in situ conditions, as illustrated in the paper. Copyright © 2010 John Wiley & Sons, Ltd.     

Upward infiltration–evaporation method to estimate soil hydraulic properties

Determination of saturated hydraulic conductivity, K_s, and the van Genuchten water retention curve θ(h) parameters is crucial in evaluating unsaturated soil water flow. The aim of this work is to present a method to estimate K_s, α and n from numerical analysis of an upward infiltration process at saturation (Cap0), with (Cap0 + h) and without (Cap0) an overpressure step (h) at the end of the wetting phase, followed by an evaporation process (Evap). The HYDRUS model as well as a brute-force search method were used for theoretical loam soil parameter estimation. The uniqueness and the accuracy of solutions from the response surfaces, K_s–n, α–n and K_s–α, were evaluated for different scenarios. Numerical experiments showed that only the Cap0 + Evap and Cap0 + h + Evap scenarios were univocally able to estimate the hydraulic properties. The method gave reliable results in sand, loam and clay-loam soils.     

Using a pore‐scale model to quantify the effect of particle re‐arrangement on pore structure and hydraulic properties

A pore‐scale model based on measured particle size distributions has been used to quantify the changes in pore space geometry of packed soil columns resulting from a dilution in electrolyte concentration from 500 to 1 mmol l^?1 NaCl during leaching. This was applied to examine the effects of particle release and re‐deposition on pore structure and hydraulic properties. Two different soils, an agricultural soil and a mining residue, were investigated with respect to the change in hydraulic properties. The mining residue was much more affected by this process with the water saturated hydraulic conductivity decreasing to 0·4% of the initial value and the air‐entry value changing from 20 to 50 cm. For agricultural soil, there was little detectable shift in the water retention curve but the saturated hydraulic conductivity decreased to 8·5% of the initial value. This was attributed to localized pore clogging (similar to a surface seal) affecting hydraulic conductivity, but not the microscopically measured pore‐size distribution or water retention. We modelled the soil structure at the pore scale to explain the different responses of the two soils to the experimental conditions. The size of the pores was determined as a function of deposited clay particles. The modal pore size of the agricultural soil as indicated by the constant water retention curve was 45 µm and was not affected by the leaching process. In the case of the mining residue, the mode changed from 75 to 45 µm. This reduction of pore size corresponds to an increase of capillary forces that is related to the measured shift of the water retention curve. Copyright © 2007 John Wiley & Sons, Ltd.     

In situ Mixing of Organic Matter Decreases Hydraulic Conductivity of Denitrification Walls in Sand Aquifers

In a previous study, a denitrification wall was constructed in a sand aquifer using sawdust as the carbon substrate. Ground water bypassed around this sawdust wall due to reduced hydraulic conductivity. We investigated potential reasons for this by testing two new walls and conducting laboratory studies. The first wall was constructed by mixing aquifer material in situ without substrate addition to investigate the effects of the construction technique (mixed wall). A second, biochip wall, was constructed using coarse wood chips to determine the effect of size of the particles in the amendment on hydraulic conductivity. The aquifer hydraulic conductivity was 35.4 m/d, while in the mixed wall it was 2.8 m/d and in the biochip wall 3.4 m/d. This indicated that the mixing of the aquifer sands below the water table allowed the particles to re-sort themselves into a matrix with a significantly lower hydraulic conductivity than the process that originally formed the aquifer. The addition of a coarser substrate in the biochip wall significantly increased total porosity and decreased bulk density, but hydraulic conductivity remained low compared to the aquifer. Laboratory cores of aquifer sand mixed under dry and wet conditions mimicked the reduction in hydraulic conductivity observed in the field within the mixed wall. The addition of sawdust to the laboratory cores resulted in a significantly higher hydraulic conductivity when mixed dry compared to cores mixed wet. This reduction in the hydraulic conductivity of the sand/sawdust cores mixed under saturated conditions repeated what occurred in the field in the original sawdust wall. This indicated that laboratory investigations can be a useful tool to highlight potential reductions in field hydraulic conductivities that may occur when differing materials are mixed under field conditions.     

Water flow path interactions with soil hydraulic properties in till soil at Gårdsjön, Sweden

In 1989, in a hydrological research programme within a deacidification project in the Gårdsjön area in southwest Sweden, flow paths and residence times of soil water and groundwater in microcatchments were examined to support the interpretation of the hydrochemical changes. Saturated hydraulic conductivity and soil water retention were analysed on more than 100 cylinder samples. The catchments have shallow sandy-silty till soil with a mean depth in the main catchment of 43 cm. Porosity of the mineral soil in the main catchment was high and ranged from 38 to 85%. The samples from the B-horizon had generally higher porosity. Porosity and the content of organic matter were correlated. The soil water retention was relatively high at all tensions, likely owing to the high content of organic matter. Dissolved organic substances were most probably transported from the shallow soil on the steep sides of the catchment down to the valley where it precipitated. The high porosities could be a consequence of long-term weathering, provided that the organic substances present have increased the leaching of the weathering products. Measured values of saturated hydraulic conductivity were close to log-normally distributed with a mean for all samples of 3 × 10⁻⁵ m s⁻¹. There was a significant increase in conductivity toward the ground surface with the mean conductivity of the samples in the uppermost 10 cm of the mineral soil of 4 × 10⁻⁵ m s⁻¹, which was about 13 times higher than the conductivity of 3 × 10⁻⁶ m s⁻¹ at 1 m depth. From the relationship between runoff at the catchment outlet and groundwater levels, the conductivity was estimated to be 15–200 times higher in the upper soil layer than in the deeper ones. In one profile, 44–64% of the yearly lateral flow was estimated to occur above 30 cm depth. The conductivity was correlated with the content of drainable water, which indicated the importance of the largest pores for the saturated hydraulic conductivity.     

Test of rock physics models for prediction of seismic velocities in shallow unconsolidated sands: a well log data case‡

This paper tests the ability of various rock physics models to predict seismic velocities in shallow unconsolidated sands by comparing the estimates to P and S sonic logs collected in a shallow sand layer and ultrasonic laboratory data of an unconsolidated sand sample. The model fits are also evaluated with respect to the conventional model for unconsolidated sand. Our main approach is to use Hertz‐Mindlin and Walton contact theories, assuming different weight fractions of smooth and rough contact behaviours, to predict the elastic properties of the high porosity point. Using either the Hertz‐Mindlin or Walton theories with rough contact behaviour to define the high porosity endpoint gives an over‐prediction of the velocities. The P‐velocity is overpredicted by a factor of ～1.5 and the S‐velocity by a factor of ～1.8 for highly porous gas‐sand. The degree of misprediction decreases with increasing water saturation and porosity.Using the Hertz‐Mindlin theory with smooth contact behaviour or weighted Walton models gives a better fit to the data, although the data are best described using the Walton smooth model. To predict the properties at the lower porosities, the choice of bounding model attached to the Walton Smooth model controls the degree of fit to the data, where the Reuss bound best captures the porosity variations of dry and wet sands in this case since they are caused by depositional differences. The empirical models based on lab experiments on unconsolidated sand also fit the velocity data measured by sonic logs in situ, which gives improved confidence in using lab‐derived results.     

Estimation of soil hydraulic properties and their uncertainty through the Beerkan infiltration experiment

The Beerkan method based on in situ single‐ring water infiltration experiments along with the relevant specific Beerkan estimation of soil transfer parameters (BEST) algorithm is attractive for simple soil hydraulic characterization. However, the BEST algorithm may lead to erroneous or null values for the saturated hydraulic conductivity and sorptivity especially when there are only few infiltration data points under the transient flow state, either for sandy soil or soils in wet conditions. This study developed an alternative algorithm for analysis of the Beerkan infiltration experiment referred to as BEST‐generalized likelihood uncertainty estimation (GLUE). The proposed method estimates the scale parameters of van Genuchten water retention and Brooks–Corey hydraulic conductivity functions through the GLUE methodology. The GLUE method is a Bayesian Monte Carlo parameter estimation technique that makes use of a likelihood function to measure the goodness‐of‐fit between modelled and observed data. The results showed that using a combination of three different likelihood measurements based on observed transient flow, steady‐state flow and experimental steady‐state infiltration rate made the BEST‐GLUE procedure capable of performing an efficient inverse analysis of Beerkan infiltration experiments. Therefore, it is more applicable for a wider range of soils with contrasting texture, structure, and initial and saturated water content. Copyright © 2015 John Wiley & Sons, Ltd.     

Influence of antecedent soil moisture on hydraulic conductivity in a series of texture‐contrast soils

Antecedent soil moisture significantly influenced the hydraulic conductivity of the A1, A2e and B21 horizons in a series of strong texture‐contrast soils. Tension infiltration at six supply potentials demonstrated that in the A1 horizon, hydraulic conductivity was significantly lower in the ‘wet’ treatment than in the ‘dry’ treatment. However in the A2e horizon, micropore and mesopore hydraulic conductivity was lower in the ‘dry’ treatment than the ‘wet’ treatment, which was attributed to the precipitation of soluble amorphous silica. In the B21 horizon, desiccation of vertic clays resulted in the formation of shrinkage cracks which significantly increased near‐saturated hydraulic conductivity and prevented the development of subsurface lateral flow in the ‘dry’ treatment. In the ‘wet’ treatment, the difference between the hydraulic conductivity of the A1 and B21 horizons was reduced; however, lateral flow still occurred in the A1 horizon due to difficulty displacing existing soil water further down the soil profile. Results demonstrate the need to account for temporal variation in soil porosity and hydraulic conductivity in soil‐water model conceptualisation and parameterisation. Copyright © 2011 John Wiley & Sons, Ltd.     

Pedotransfer functions estimating soil hydraulic properties using different soil parameters

Estimates of soil hydraulic properties using pedotransfer functions (PTF) are useful in many studies such as hydrochemical modelling and soil mapping. The objective of this study was to calibrate and test parametric PTFs that predict soil water retention and unsaturated hydraulic conductivity parameters. The PTFs are based on neural networks and the Bootstrap method using different sets of predictors and predict the van Genuchten/Mualem parameters. A Danish soil data set (152 horizons) dominated by sandy and sandy loamy soils was used in the development of PTFs to predict the Mualem hydraulic conductivity parameters. A larger data set (1618 horizons) with a broader textural range was used in the development of PTFs to predict the van Genuchten parameters. The PTFs using either three or seven textural classes combined with soil organic mater and bulk density gave the most reliable predictions of the hydraulic properties of the studied soils. We found that introducing measured water content as a predictor generally gave lower errors for water retention predictions and higher errors for conductivity predictions. The best of the developed PTFs for predicting hydraulic conductivity was tested against PTFs from the literature using a subdata set of the data used in the calibration. The test showed that the developed PTFs gave better predictions (lower errors) than the PTFs from the literature. This is not surprising since the developed PTFs are based mainly on hydraulic conductivity data near saturation and sandier soils than the PTFs from the literature. Copyright © 2007 John Wiley & Sons, Ltd.     

Imbibition of saline solutions into dry and prewetted porous media

Infiltration of saline solutions and pure water into homogeneous packs of prewetted and air-dry silica sands was investigated using a light transmission system. Four sand grades and five solutions were considered. Narrow fingers with a sharp, almost saturated, wetting front were observed in the air-dry sands. The water content left behind the fingertip of saline solutions was higher than for pure water, resulting in a greater lateral expansion of the saline fingers over time. The rate of lateral expansion scaled with the square root of time, likely due to classic liquid sorption with the possible addition of water vapor diffusion. At early time, the salty fingers moved faster, but were ultimately overtaken by the pure water fingers. In prewetted sand, the wetting fronts were diffuse and never exceeded 26% saturation, less than third that seen in the initially air-dry media. The plumes in the prewetted sand were also much wider and their shape varied. In the prewetted sand the elevated surface tension of the saline solutions was the major cause for the observed differences in finger width and velocity, yet appeared to be insignificant in air-dry sand. Here, in addition to the density effect, absorption of the saline solution to the silica sand influenced the depth of wetting, finger velocity, and subsequent lateral expansion.     

Evolution of Heterogeneous Mixed Siliciclastic/Carbonate Aquifers Containing Metastable Sediments

Mixed carbonate and siliciclastic marine sediments commonly become freshwater aquifers in eastern coastal regions of the United States and many other global locations. As these deposits age, the carbonate fraction of the sediment is commonly removed by dissolution and the aquifer can become a solely siliciclastic system or contain zones or beds of pure quartz sand. During aquifer evolution, the sediment grain size characteristics, hydraulic conductivity, and porosity change. An investigation of these changes using mixed carbonate/siliciclastic sediment samples collected from a modern barrier island beach in southern Florida showed that the average mean grain diameter decreased with removal of the carbonate fraction, but the average hydraulic conductivity and porosity increased slightly, but not to statistical significance. This counterintuitive result occurs because of the change in the pore types from a combined shelter and intergranular pore system producing a dual porosity system in the mixed sediments to a single intergranular pore system in the siliciclastic sediment fraction. Within the mixed carbonate/siliciclastic sediment, in the pure carbonate fraction, large shell fractions form grain‐supported large pores, which become filled with sand‐sized quartz as the shell fragments decrease in size or as the sediment becomes compacted. The hydraulic conductivity increases because the shell fragments that were oriented perpendicular to flow caused an increase in the length of the flow path, or a larger scale tortuosity, compared with the flow through pure quartz sand.     

Geochemical and Physical Controls on Vadose Zone Hydrology of Holocene Carbonate Sands,Grand Bahama Island

This paper explores the relationship between vadose zone hydrology and geochemical changes in mixed mineralogy carbonate sands from a Bahamian coastal dune of Holocene age. Cores were taken from two sites: at site A, a shallow humic Entisol is developed beneath open scrub vegetation, while at site B a deeper, more organic-rich Inceptisol has formed beneath a mature hardwood coppice. X-ray diffraction analysis reveals significant contrasts in mineralogy both within and between the two sites, with partial stabilization of high-Mg calcite and aragonite, to low-Mg calcite. Stabilization is greater at site B, and is accompanied by a significant increase in total porosity. Diagenetic changes in pore-size distribution have implications for residence times of percolating water, as determined using measurements of moisture retention characteristics using pressure plate apparatus, and hydrological models of unsaturated zone moisture flux. The diagenetically more mature sands from site B have a 50–100 per cent higher moisture retention, although unsaturated hydraulic conductivity is also higher, particularly at greater suctions. The increase in water retention is likely to enhance further rates of mineral-controlled reactions, while development of an organic-rich soil also enhances the geochemical drive for dissolution. Carbonate diagenesis thus appears to be strongly linked to vadose zone hydrology, and the interactions identified here have important consequences for the nature and long-term rates of mineral stabilization. © 1997 by John Wiley & Sons, Ltd.