Functional behaviour of pedotransfer functions in soil water flow simulation

Soil hydraulic parameter values for large‐scale modelling cannot be obtained by direct methods. Pedotransfer functions (PTFs) that relate soil hydraulic properties (SHPs) to generally available soil texture data may provide an alternative. A considerable number of PTF models has been developed, the application of three recent PTFs is evaluated. As a first step sets of SHPs derived from the PTFs are compared with measured sets of SHPs for three sites. No good agreement was found statistically between measured and PTF results or between PTF results. As a second step and from a practical point of view results for three hydrologically functional variables were compared and evaluated. The three selected functional variables are saturated hydraulic conductivity, k₀, in relation to infiltration excess runoff, available soil water amounts for evapotranspiration and water table depth for a specified upward flux or capillary rise. Derived k₀ distributions from PTFs show substantially less variance than from the measured data at all three sites. This can have a considerable impact on infiltration excess runoff, depending on the actual rainfall regime. Simulated available soil water amounts for evapotranspiration for some combinations of PTFs and sites are close to those obtained for measured SHPs, however, no consistency in results can be detected. Water table depths for specified upward flux densities using PTF derived SHPs are generally deeper than those based on measured SHPs and means a potentially higher water availability. Overall, differences in capillary rise among the selected PTFs and between measured and PTF based results are again inconsistent and show no clear relationship with soil texture. Finally, as a third step, effective SHPs were calculated by using spatially averaged texture as PTF input representing areal average behaviour. For these effective SHPs the calculated effective values for the three selected functional variables appear to be close to the areally averaged values obtained with step 2. The selected functional variables thus appear to depend linearly on the PTFs over the range for which the data are representative. This suggests that for our specific PTFs areal mean or effective values for the three functional variables can be obtained fairly accurately from a single measurement of a bulk collection of soil samples as input. Copyright © 2003 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.     

Hydraulic properties of peat soils along a bulk density gradient—A meta study

Our understanding of hydraulic properties of peat soils is limited compared with that of mineral substrates. In this study, we aimed to deduce possible alterations of hydraulic properties of peat soils following degradation resulting from peat drainage and aeration. A data set of peat hydraulic properties (188 soil water retention curves [SWRCs], 71 unsaturated hydraulic conductivity curves [UHCs], and 256 saturated hydraulic conductivity [K_s] values) was assembled from the literature; the obtained data originated from peat samples with an organic matter (OM) content ranging from 23 to 97 wt% (weight percent; and according variation in bulk density) representing various degrees of peat degradation. The Mualem‐van Genuchten model was employed to describe the SWRCs and UHCs. The results show that the hydraulic parameters of peat soils vary over a wide range confirming the pronounced diversity of peat. Peat decomposition significantly modifies all hydraulic parameters. A bulk density of approximately 0.2 g cm^?3 was identified as a critical threshold point; above and below this value, macroporosity and hydraulic parameters follow different functions with bulk density. Pedotransfer functions based on physical peat properties (e.g., bulk density and soil depth) separately computed for bog and fen peat have significantly lower mean square errors than functions obtained from the complete data set, which indicates that not only the status of peat decomposition but also the peat‐forming plants have a large effect on hydraulic properties. The SWRCs of samples with a bulk density of less than 0.2 g cm^?3 could be grouped into two to five classes for each peat type (botanical composition). The remaining SWRCs originating from samples with a bulk density of larger than 0.2 g cm^?3 could be classified into one group. The Mualem‐van Genuchten parameter values of α can be used to estimate K_s if no K_s data are available. In conclusion, the derived pedotransfer functions provide a solid instrument to derive hydraulic parameter values from easily measurable quantities; however, additional research is required to reduce uncertainty.     

Performance of soil particle-size distribution models for describing deposited soils adjacent to constructed dams in the China Loess Plateau

Soil particle-size distributions (PSD) have been used to estimate soil hydraulic properties. Various parametric PSD models have been proposed to describe the soil PSD from sparse experimental data. It is important to determine which PSD model best represents specific soils. Fourteen PSD models were examined in order to determine the best model for representing the deposited soils adjacent to dams in the China Loess Plateau; these were: Skaggs (S-1, S-2, and S-3), fractal (FR), Jaky (J), Lima and Silva (LS), Morgan (M), Gompertz (G), logarithm (L), exponential (E), log-exponential (LE), Weibull (W), van Genuchten type (VG) as well as Fredlund (F) models. Four-hundred and eighty samples were obtained from soils deposited in the Liudaogou catchment. The coefficient of determination (R ²), the Akaike’s information criterion (AIC), and the modified AIC (mAIC) were used. Based upon R ² and AIC, the three- and four-parameter models were both good at describing the PSDs of deposited soils, and the LE, FR, and E models were the poorest. However, the mAIC in conjunction with R ² and AIC results indicated that the W model was optimum for describing PSD of the deposited soils for emphasizing the effect of parameter number. This analysis was also helpful for finding out which model is the best one. Our results are applicable to the China Loess Plateau.     

Sensitivity of soil parameters in unsaturated zone modelling and the relation between effective,laboratory and in situ estimates

Simulation of soil moisture content requires effective soil hydraulic parameters that are valid at the modelling scale. This study investigates how these parameters can be estimated by inverse modelling using soil moisture measurements at 25 locations at three different depths (at the surface, at 30 and 60 cm depth) on an 80 by 20 m hillslope. The study presents two global sensitivity analyses to investigate the sensitivity in simulated soil moisture content of the different hydraulic parameters used in a one‐dimensional unsaturated zone model based on Richards' equation. For estimation of the effective parameters the shuffled complex evolution algorithm is applied. These estimated parameters are compared to their measured laboratory and in situ equivalents. Soil hydraulic functions were estimated in the laboratory on 100 cm³ undisturbed soil cores collected at 115 locations situated in two horizons in three profile pits along the hillslope. Furthermore, in situ field saturated hydraulic conductivity was estimated at 120 locations using single‐ring pressure infiltrometer measurements. The sensitivity analysis of 13 soil physical parameters (saturated hydraulic conductivity (K_s), saturated moisture content (θ_s), residual moisture content (θ_r), inverse of the air‐entry value (α), van Genuchten shape parameter (n), Averjanov shape parameter (N) for both horizons, and depth (d) from surface to B horizon) in a two‐layer single column model showed that the parameter N is the least sensitive parameter. K_s of both horizons, θ_s of the A horizon and d were found to be the most sensitive parameters. Distributions over all locations of the effective parameters and the distributions of the estimated soil physical parameters from the undisturbed soil samples and the single‐ring pressure infiltrometer estimates were found significantly different at a 5% level for all parameters except for α of the A horizon and K_s and θ_s of the B horizon. Different reasons are discussed to explain these large differences. Copyright © 2004 John Wiley & Sons, Ltd.     

Comparing hydraulic properties of different forest floors

The forest floor plays an important role in runoff rate, soil erosion and soil infiltration capacity by protecting mineral soils from the direct impact of falling raindrops. Forest floor consists of different kinds of litter with different hydraulic properties. In this study, the inverse method was used to estimate the hydraulic properties of three kinds of forest floor (broad‐leaved, needle‐leaved and mixed‐stand) at three replications in a completely random design. Forest floor samples were collected from the Gilan Province, Iran. The samples were piled up to make long columns 40.88 cm high with an inner diameter of 18.1 cm. Artificial rainfall experiments were conducted on top of the columns, and free drainage from the bottom of the columns was measured in the laboratory. Saturated hydraulic conductivity (K_s), saturated water content and water retention curve parameters (van Genuchten equation) were estimated by the inverse method. The results showed that the K_s of needle‐leaved samples differed significantly (p < 0.05) from those of broad‐leaved and mixed‐stand samples, whereas the latter two did not differ in this regard. No significant differences emerged in the water retention function parameters of van Genuchten (θ_r, β and α) in the three forest floor samples. The saturated water content of mixed‐stand samples was significantly different (p < 0.05) from that of broad‐leaved and needle‐leaved treatments with the latter two samples showing no significant difference. The good agreement between simulated and observed free drainage for all forest floor samples in the validation period indicates that the estimated hydraulic properties efficiently characterize the unsaturated water flow in the forest floor. Copyright © 2013 John Wiley & Sons, Ltd.     

Infiltration and runoff generation processes in fire‐affected soils

Post‐wildfire runoff was investigated by combining field measurements and modelling of infiltration into fire‐affected soils to predict time‐to‐start of runoff and peak runoff rate at the plot scale (1 m²). Time series of soil‐water content, rainfall and runoff were measured on a hillslope burned by the 2010 Fourmile Canyon Fire west of Boulder, Colorado during cyclonic and convective rainstorms in the spring and summer of 2011. Some of the field measurements and measured soil physical properties were used to calibrate a one‐dimensional post‐wildfire numerical model, which was then used as a ‘virtual instrument’ to provide estimates of the saturated hydraulic conductivity and high‐resolution (1 mm) estimates of the soil‐water profile and water fluxes within the unsaturated zone. Field and model estimates of the wetting‐front depth indicated that post‐wildfire infiltration was on average confined to shallow depths less than 30 mm. Model estimates of the effective saturated hydraulic conductivity, K_s, near the soil surface ranged from 0.1 to 5.2 mm h^?1. Because of the relatively small values of K_s, the time‐to‐start of runoff (measured from the start of rainfall), t_p, was found to depend only on the initial soil‐water saturation deficit (predicted by the model) and a measured characteristic of the rainfall profile (referred to as the average rainfall acceleration, equal to the initial rate of change in rainfall intensity). An analytical model was developed from the combined results and explained 92–97% of the variance of t_p, and the numerical infiltration model explained 74–91% of the variance of the peak runoff rates. These results are from one burned site, but they strongly suggest that t_p in fire‐affected soils (which often have low values of K_s) is probably controlled more by the storm profile and the initial soil‐water saturation deficit than by soil hydraulic properties. Published 2013. This article is a U.S. Government work and is in the public domain in the USA.     

Effects of biological crust coverage on soil hydraulic properties for the Loess Plateau of China

Biological soil crusts (BSCs) are ubiquitous living covers that have been allowed to grow on abandoned farmlands over the Loess Plateau because the “Grain for Green” project was implemented in 1999 to control serious soil erosion. However, few studies have been conducted to quantify the effects of BSC coverage on soil hydraulic properties. This study was performed to assess the effects of BSC coverage on soil hydraulic properties, which are reflected by the soil sorptivity under an applied pressure of 0 (S ₀ ) and ?3 (S ₃ ) cm, saturated hydraulic conductivity (K _s ), wetting front depth (WFD ), and mean pore radius (λ _m ), for the Loess Plateau of China. Five classes of BSC coverage (i.e., 1–20%, 20–40%, 40–60%, 60–80%, and 80–100%) and a bare control were selected at both cyanobacteria‐ and moss‐covered sites to measure soil hydraulic properties using a disc infiltrometer under 2 consecutive pressure heads of 0 and ?3 cm, allowing the direct calculation of S ₀ , S ₃ , K _s , and λ _m . The WFD was measured onsite using a ruler immediately after the experiments of infiltration. The results indicated that both cyanobacteria and moss crusts were effective in changing the soil properties and impeding soil infiltration. The effects of moss were greater than those of cyanobacteria. Compared to those of the control, the S ₀ , S ₃ , K _s , WFD , and λ _m values of cyanobacteria‐covered soils were reduced by 13.7%, 11.0%, 13.3%, 10.6%, and 12.6% on average, and those of moss‐covered soils were reduced by 27.6%, 22.1%, 29.5%, 22.2%, and 25.9%, respectively. The relative soil sorptivity under pressures of 0 (RS ₀ ) and ?3 (RS ₃ ) cm, the relative saturated hydraulic conductivity (RK _s ), the relative wetting front depth (RWFD ), and the relative mean pore radius (Rλ _m ) decreased exponentially with coverage for both cyanobacteria‐ and moss‐covered soils. The rates of decrease in RS ₀ , RS ₃ , RK _s , RWFD , and Rλ _m of cyanobacteria were significantly slower than those of moss, especially for the coverage of 0–40%, with smaller ranges. The variations of soil hydraulic properties with BSC coverage were closely related to the change in soil clay content driven by the BSC coverage on the Loess Plateau. The results are useful for simulating the hydraulic parameters of BSC‐covered soils in arid and semiarid areas.     

Parameter estimation for multiple post-wildfire hydrologic models

Predictions of post-wildfire flooding and debris flows are needed, typically with short lead times. Measurements of soil-hydraulic properties necessary for model parameterization are, however, seldom available. This study quantified soil-hydraulic properties, soil-water retention, and selected soil physical properties within the perimeter of the 2017 Thomas Fire in California. The Thomas Fire burn scar produced catastrophic debris flows in January 2018, highlighting the need for improved prediction capability. Soil-hydraulic properties were also indirectly estimated using relations tied to soil-water retention. These measurements and estimates are examined in the context of parameterizing post-wildfire hydrologic models. Tension infiltrometer measurements showed significant decreases (p < .05) in field-saturated hydraulic conductivity (K_fs) and sorptivity (S) in burned areas relative to unburned areas. Wildfire effects on soil water-retention were dominated by significant decreases in saturated soil-water content (θ_S). The van Genuchten parameters α, N, and residual water content did not show significant wildfire effects. The impacts of the wildfire on hydraulic and physical soil properties were greatest in the top 1 cm, emphasizing that measurements of post-fire soil properties should focus on the near-surface. Reductions in K_fs, θ_s, and soil-water retention in burned soils were attributed to fire-induced decreases in soil structure evidenced by increases in dry bulk density. Sorptivity reductions in burned soils were attributed to increases in soil-water repellency. Rapid post-fire assessments of flash flood and debris flow hazards using physically-based hydrologic models are facilitated by similarities between K_fs, S, and the Green–Ampt wetting front potential (ψ_f) with measurements at other southern CA burned sites. We suggest that ratios of burned to unburned K_fs (0.37), S (0.36), and ψ_f (0.66) could be used to scale unburned values for model parameterization. Alternatively, typical burned values (K_fs = 20 mm hr⁻¹; S = 6 mm hr^−0.5; ψ_f = 1.6 mm) could be used for model parameterization.     

Effects of biological soil crusts on soil detachment process by overland flow in the Loess Plateau of China

Biological soil crusts (BSCs) cover up to 60 to 70% of the soil surface in grasslands after the ‘Grain for Green’ project was implemented in 1999 to rehabilitate the Loess Plateau. However, few studies exist that quantify the effects of BSCs on the soil detachment process by overland flow in the Loess Plateau. This study investigated the potential effects of BSCs on the soil detachment capacity (D_c), and soil resistance to flowing water erosion reflected by rill erodibility and critical shear stress. Two dominant BSC types that developed in the Loess Plateau (the later successional moss and the early successional cyanobacteria mixed with moss) were tested against natural soil samples collected from two abandoned farmland areas. The samples were subjected to flow scouring under six different shear stresses ranging from 7.15 to 24.08 Pa. The results showed that D_c decreased significantly with crust coverage under both moss and mixed crusts. The mean D_c of bare soil (0.823 kg m^?2 s^?1) was 2.9 to 48.4 times greater than those of moss covered soil (0.017–0.284 kg m^?2 s^?1), while it (3.142 kg m^?2 s^?1) was 4.9 to 149.6 times greater than those of mixed covered soil (0.021–0.641 kg m^?2 s^?1). The relative detachment rate of BSCs compared with bare soils decreased exponentially with increasing BSC coverage for both types of BSCs. The D_c value can be simulated by flow shear stress, cohesion, and BSC coverage using a power function (NSE ≥ 0.59). Rill erodibility also decreased with coverage of both crust types. Rill erodibility of bare soil was 3 to 74 times greater than those of moss covered soil and was 2 to 165 times greater than those of mixed covered soil. Rill erodibility could also be estimated by BSC coverage in the Loess Plateau (NSE ≥ 0.91). The effect of crust coverage on critical shear stress was not significant. Copyright © 2015 John Wiley & Sons, Ltd.     

Hydraulic properties of concentrated flow of a bank gully in the dry‐hot valley region of southwest China

To quantify spatiotemporal variation in hydraulic properties of bank gully concentrated flow, a series of scour experiments were run under water discharge rates ranging from 30 to 120 l min^?1. Concentrated flows were found to be turbulent and supercritical in the upstream catchment area and downstream gully beds. As discharge increased, values of the soil erosion rate, Reynolds number (Re), shear stress, stream power, and flow energy consumption (ΔE) increased while values of the Froude number (Fr) and the Darcy–Weisbach friction factor (resistance f ) did not. With the exception of gully headcut collapse under discharge rates of 60, 90, and 120 l min^?1, a declining power function trend (P < 0.05) in the soil erosion rate developed in the upstream catchment area, headcuts, and downstream gully beds. However, increasing trends were observed in temporal variations of hydraulic properties for downstream gully beds and the upstream catchment area. Despite significant differences in temporal variation between the soil erosion rate and hydraulic property values, relative steady state conditions of the soil erosion rate and ΔE were attained following an initial period of adjustment in the upstream catchment area, headcuts, and downstream gully beds under different discharge rates. A logarithmic growth of flow energy consumption per unit soil loss (ΔE_u) was observed in bank gullies and the upstream catchment area as the experiment progressed, further illustrating the actual reason behind the discrepancy in temporal variation between soil erosion rates and ΔE. Results demonstrate that ΔE can be used to estimate headcut erosion soil loss, but further quantitative studies are required to quantify coupling effects between hydraulic properties and vertical variation in soil mechanical properties on temporal variation for bank gully soil erosion rates. Copyright © 2015 John Wiley & Sons, Ltd.     

Thresholds and relations for soil‐hydraulic and soil‐physical properties as a function of burn severity 4 years after the 2011 Las Conchas Fire,New Mexico,USA

Wildfire effects on soil‐physical and ‐hydraulic properties as a function of burn severity are poorly characterized, especially several years after wildfire. A stratified random sampling approach was used in 2015 to sample seven sites representing a spectrum of remotely sensed burn severity in the area impacted by the 2011 Las Conchas Fire in New Mexico, USA. Replicate samples from each site were analysed in the laboratory. Linear and linear indicator regression were used to assess thresholds in soil‐physical and ‐hydraulic properties and functional relations with remotely sensed burn severity. Significant thresholds were present for initial soil‐water content (θ_i) at 0–6 cm depth between the change in the Normalized Burn Ratio (dNBR) equal to 618–802, for bulk density (ρ_b) at 3–6 cm between dNBR equal to 416–533, for gravel fraction at 0–1 cm between dNBR equal to 416–533, for fines (the silt + clay fraction) at 0–1 cm for dNBR equal to 416–533, and for fines at 3–6 cm for dNBR equal to 293–416. Significant linear relations with dNBR were present between ρ_b at 0–1 cm, loss on ignition (LOI) at 0–1 cm, gravel fraction at 0–1 cm, and the large organic fraction at 1–3 cm. No thresholds or effects on soil‐hydraulic properties of field‐saturated hydraulic conductivity or sorptivity were observed. These results suggest that ρ_b and LOI at 0–1 cm have residual direct impacts from the wildfire heat impulse. The θ_i threshold is most likely from delayed groundcover/vegetation recovery that increases evaporation at the highest burn severity sites. Gravel and silt + clay thresholds at 0–1 cm at the transition to high burn severity suggest surface gravel lag development from hydraulic erosion. Thresholds in ρ_b from 3 to 6 cm and in silt + clay fraction from 3 to 6 cm appear to be the result of soil variability between sites rather than wildfire impacts. This work suggests that gravel‐rich soils may have increased resilience to sustained surface runoff generation and erosion following wildfire, with implications for assessments of postwildfire hydrologic and erosion recovery potential.     

Seepage mechanisms and permeability differences between loess and paleosols in the critical zone of the Loess Plateau

The soil in the Loess Plateau has special permeability characteristics due to the alternating distribution of loess and paleosol layers. Using an analysis of the physical properties, microstructure and thermogravimetric analysis of loess and paleosol, this paper examines strata seepage mechanisms in the loess tableland area and considers the applicability of a hydraulic conductivity empirical formula. The analysis shows that hydraulic conductivity attenuation with depth can be represented by a negative exponential model, while hydraulic conductivity values are not normally distributed. The best-fitting models of hydraulic conductivity in the horizontal (K_H) and vertical (K_V) directions are Gaussian models, and both have strong spatial correlations. This study of the difference in permeability between loess and paleosol found that the depositional environment was the dominant causal factor, making the average hydraulic conductivity of most loess layers greater than that of the underlying paleosol layers. Different microstructures between loess and paleosol also confirmed the microscopic explanation in permeability anisotropy and their permeability difference. Thermogravimetric analysis determined temperature ranges for different types of water lost by heat, and then calculated ratios of bound water mass to liquid limit, with an average of 0.768. A modified formula suitable for loess was obtained by integrating the consistency index method and effective porosity ratio model into the hydraulic conductivity empirical formula. Compared with the results of laboratory tests and uncorrected formulas, the modified formula provides a good estimate of strata hydraulic conductivity. Accurate understandings of seepage mechanisms and permeability differences in the loess area are important, promoting ecological restoration and providing scientific guidance for the sustainable development of the Loess Plateau.     

Vertical variations and transport mechanism of soil moisture in response to vegetation restoration on the Loess Plateau of China

Soil moisture is essential for vegetation restoration in arid and semi-arid regions. Ascertaining the vertical distribution and transportation of soil moisture under different vegetation types has a profound effect on the ecological construction. In this study, the soil moisture at a depth of 500 cm for four typical vegetation types, including Robinia pseudoacacia, Caragana korshinskii, Stipa bungeana, and corn, were investigated and compared in the Zhifanggou watershed of the Loess plateau. Additionally, hydrogen and oxygen stable isotopes were detected to identify the transport mechanism of soil moisture. The results showed vertical distribution and transportation of soil moisture were different under different vegetation types. Depth-averaged soil moisture under S. bungeana and corn generally increased along the profile, while C. korshinskii and R. pseudoacacia showed weakly increasing and relatively stable after an obvious decreasing trend (0–40 cm). The soil moisture under R. pseudoacacia was lower than that under other vegetation types, especially in deep layer. However, the effect of R. pseudoacacia on soil moisture in the topsoil (< 30 cm) could be positive. For R. pseudoacacia (160–500 cm), C. korshinskii (0–500 cm), and S. bungeana (0–100 cm), the soil moisture declined with increased in vegetation age. Planting arbor species such as R. pseudoacacia intensified the decline of soil moisture on the Loess Plateau. The capacity of evaporation fractionation of soil moisture followed the sequence: corn > S. bungeana > R. pseudoacacia > C. korshinskii. The δ¹⁸O values in soil water fluctuated across the profile. The δ¹⁸O values changed sharply in upper layer and generally remained stable in deep layer. However, in middle layer, the vertical distribution characteristics of the δ¹⁸O values were different under different vegetation types. We estimated that piston flow was the main mode of precipitation infiltration, and the occurrence of preferential flow was related to vegetation types. These results were helpful to improve the understanding of the response of deep soil moisture to vegetation restoration and inform practices for sustainable water management.     

Geomorphic controls on the physical and hydrologic properties of soils in a valley floor

The distribution of soil hydraulic and physical properties strongly influences runoff processes in landscapes. Although much work has been done to quantify and predict the properties of hillslope soils, far less is known about the distribution of soil properties in valley floors. A technique that links the estimation and distribution of soil hydraulic properties in valleys, with easily identified geomorphic features, was developed along a 2 km length of a valley at Brooks Creek in New South Wales, Australia. Soil physical and hydraulic property data were collected across a set of floodplain and fan features within the valley and analysed statistically to determine if soil properties varied significantly between geomorphic features and stratigraphic layers. The results show that the depth‐averaged saturated hydraulic conductivity, K_s, of the soil varies significantly with landform: fan units have K_g values that are twice that of floodplains and colluvial toeslope deposits have K_s values four times higher than floodplains. Given the notorious variability of K_s values in space, the strong statistical separation of soil properties by landform, backed up by strong separation of soil particle size by landform, suggests a way forward in understanding the distribution of soil properties in valleys and their influence on catchment hydrology. Copyright © 2000 John Wiley & Sons, Ltd.     

Temporal variation in soil detachment under different land uses in the Loess Plateau of China

Measurements of temporal variations in soil detachability under different land uses are badly needed to develop new algorithms or evaluate the existing ones for temporal adjustment of soil detachability in continuous soil erosion models. Few studies have been conducted in the Loess Plateau to quantify temporal variations in detachment rate of runoff under different land uses. The objectives of this study were to investigate the temporal variations of soil detachment rate under different land uses and to further identify the potential factors causing the change in detachment rate in the Loess Plateau. Undisturbed soil samples were collected in the fields of arable land (millet, soybean, corn, and potato), grassland, shrub land, wasteland, and woodland and tested in a laboratory flume under a constant hydraulic condition. The measurements started in mid‐April and ended in early October, 2006. The results showed that soil detachment rate of each land use fluctuated considerably over time. Distinctive temporal variation in detachment rate was found throughout the summer growing season of measurement in each land use. The maximum detachment rates of different land uses varied from 0·019 to 0·490 kg m^–2 s^–1 and the minimum detachment rates ranged from 0·004 to 0·092 kg m^–2 s^–1. Statistical analysis using a paired‐samples t‐test indicated that variations in soil detachment rate differed significantly at the 0·05 level between land uses in most cases. The major factors responsible for the temporal variation of soil detachment were tillage operations (such as planting, ploughing, weeding, harvesting), soil consolidation, and root growth. The influence of tillage operations on soil detachment depended on the degree of soil disturbance caused by the operations. The consolidation of the topsoil over time after tillage was reflected by increases in soil bulk density and soil cohesion. As soil bulk density and cohesion increased, detachment rate decreased. The impact of root density was inconclusive in this study. Further studies are needed to quantify the effects of root density on temporal variations of soil detachment. This work provides useful information for developing temporal adjustments to soil detachment rate in continuous soil erosion models in the Loess Plateau. Copyright © 2009 John Wiley & Sons, Ltd.     

Carbon isotopic composition of modern soil and paleosol as a response to vegetation change on the Chinese Loess Plateau

The relationship between the carbon isotopic composition of paleosols and pale-ovegetation on the Loess Plateau is still unclear. One of the main reasons is that we are short of knowledge about the characteristics of the carbon isotopic composition of modern soil in this area. A preliminary investigation of the carbon isotopic compositions of the modern soil and the loess/paleosol sequence on the Loess Plateau shows that the carbon isotopic composition of modern soil is consistent with the distribution of modern plants on the Loess Plateau, where the ecosystem is dominated by a mixture of C4 and C3 plants. Comparing the δ13C values of modern soil and loess-paleosol sequences from the Xunyi profile, we conclude that C3 plants dominated the landscape during loess sediment stages, while C4 plants expanded during paleosol stages.     

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.     

Unsaturated hydraulic properties of xerophilous mosses: towards implementation of moss covered soils in hydrological models

Evaporation from mosses and lichens can form a major component of the water balance, especially in ecosystems where mosses and lichens often grow abundantly, such as tundra, deserts and bogs. To facilitate moss representation in hydrological models, we parameterized the unsaturated hydraulic properties of mosses and lichens such that the capillary water flow through moss and lichen material during evaporation could be assessed. We derived the Mualem‐van Genuchten parameters of the drying retention and the hydraulic conductivity functions of four xerophilous moss species and one lichen species. The shape parameters of the retention functions (2.17 < n < 2.35 and 0.08 < α < 0.13 cm^?1) ranged between values that are typical for sandy loam and loamy sand. The shapes of the hydraulic conductivity functions of moss and lichen species diverged from those of mineral soils, because of strong negative pore‐connectivity parameters (?2.840 < l < ?2.175) and low hydraulic conductivities at slightly negative pressure heads (0.016 < K₀ < 0.280 cm/d). These K₀ values are surprisingly low, considering that mosses are very porous. However, during evaporation, large pores and voids were air filled and did not participate in capillary water flow. Small K₀ values cause mosses and lichens to be conservative with water during wet conditions, thus tempering evaporation compared to mineral soils. On the other hand, under dry conditions, mosses and lichens are able to maintain a moisture supply from the soil, leading to a higher evaporation rate than mineral soils. Hence, the modulating effect of mosses on evaporation possibly differs between wet and dry climates. Copyright © 2013 John Wiley & Sons, Ltd.