Integrating observations and models to determine the effect of seasonally frozen ground on hydrologic partitioning in alpine hillslopes in the Colorado Rocky Mountains,USA

This study integrated spatially distributed field observations and soil thermal models to constrain the impact of frozen ground on snowmelt partitioning and streamflow generation in an alpine catchment within the Niwot Ridge Long-Term Ecological Research site, Colorado, USA. The study area was comprised of two contrasting hillslopes with notable differences in topography, snow depth and plant community composition. Time-lapse electrical resistivity surveys and soil thermal models enabled extension of discrete soil moisture and temperature measurements to incorporate landscape variability at scales and depths not possible with point measurements alone. Specifically, heterogenous snowpack thickness (~0–4 m) and soil volumetric water content between hillslopes (~0.1–0.45) strongly influenced the depths of seasonal frost, and the antecedent soil moisture available to form pore ice prior to freezing. Variable frost depths and antecedent soil moisture conditions were expected to create a patchwork of differing snowmelt infiltration rates and flowpaths. However, spikes in soil temperature and volumetric water content, as well as decreases in subsurface electrical resistivity revealed snowmelt infiltration across both hillslopes that coincided with initial decreases in snow water equivalent and early increases in streamflow. Soil temperature, soil moisture and electrical resistivity data from both wet and dry hillslopes showed that initial increases in streamflow occurred prior to deep soil water flux. Temporal lags between snowmelt infiltration and deeper percolation suggested that the lateral movement of water through the unsaturated zone was an important driver of early streamflow generation. These findings provide the type of process-based information needed to bridge gaps in scale and populate physically based cryohydrologic models to investigate subsurface hydrology and biogeochemical transport in soils that freeze seasonally.     

Geophysical experiments to image the shallow internal structure and the moisture distribution of a mine waste rock pile

Several field surveys of a waste rock pile were carried out during the summers of 2002 and 2003 using ground-penetrating radar, electromagnetic conductivity and DC resistivity imaging. The waste rock deposit is prone to generate acid mine drainage (AMD) due to the oxidation of sulphidic minerals. One of the most critical factors that lead to the production of AMD is unsaturated water flow and the ensuing moisture distribution in the waste rock. This geophysical characterization study, performed over a 30 m × 30 m test zone, was designed to image the internal structure controlling the water flux at shallow depth. The subsurface was found to consist of three zones for the first 6 m of the pile, mainly based on electrical resistivities: a thin superficial conductive material, an intermediate 2 to 3 m thick highly resistive zone, and a lower, more conductive medium. With the help of hydrogeological tests, chemical analyses and two 2.5 m-deep trenches, it is shown that the two conductive zones are correlated with fine-grained waste rock and the resistive zone correlates with a coarser material. In the two deeper zones, the contact between the two types of waste rock is typically highlighted by a sharp resistive/conductive boundary. An increase of conductance in the relatively thin upper layer towards the edge of the pile appears to be caused by an increase in thickness of the fine-grained material. Additional geophysical surveys carried out on a profile along the flank of the upper bench of the pile show that the main features of the internal structure are sub-parallel to the slope, at least for the first 3 m in depth. The data also show an increase in resistivity from the top to bottom of the slope, in accordance with expected particle segregation, from fine-grained material at the top to coarser material at the bottom. Wide-angle reflection GPR monitoring during large scale infiltration tests seems to indicate preferential flow paths towards the direction of coarser, more pervious material (which also appears to be less oxidized). Water preferentially flows through the coarse-grained material, but it is stored by capillary forces in the fine-grained material. Apart from the deposition methods, the results strongly suggest that factors such as machinery-induced mechanical alteration, construction history of the pile, and increased oxidization near the edges could explain the resistivity model. The model interpreted from geophysical imaging agrees well with the conceptual model of the rock pile. The resistivity and GPR methods appear to be efficient geophysical methods to characterize the internal structure and preferential flow patterns within unsaturated waste rock piles.     

Seasonal transition of hydrological processes in a slow‐moving landslide in a snowy region

A comprehensive understanding of seasonal hydrological dynamics is required to describe the influence of pore‐water pressure on the stability of landslides in snowy regions. This study reports on the results of continuous meteorological and hydrological observations over 2 years on a landslide body comprising Neogene sedimentary rocks in northern Japan, where a thick (3–5 m) seasonal snowpack covers the land surface. Monitoring of the volumetric water content in shallow unsaturated zones (<0.8 m depth) and pore‐water pressure in saturated bedrock at depths of 2.0 and 5.2 m revealed clear seasonality in hydrological responses to rainfall and meltwater supply. During snow‐free periods, both the shallow soil moisture and deep pore‐water pressure responded rapidly to intense rainwater infiltration. In contrast, during snowmelt, the deep pore pressure fluctuated in accordance with the daily cycle of meltwater input, without notable changes in shallow moisture conditions. During occasional foehn events that cause intense snow melting in midwinter, meltwater flows preferentially through the layered snowpack, converging to produce a localized water supply at the ground surface. This episodically triggers a significant rise in pore‐water pressure. The seasonal differences in hydrological responses were characterized by a set of newly proposed indices for the magnitude and quickness of increases in the pressure head near the sliding surface. Under snow‐covered conditions, the magnitude of the pressure increase tends to be suppressed, probably owing to a reduction in infiltration caused by a seasonal decrease in the permeability of surface soils, and effective pore‐water drainage through the highly conductive colluvial layer. Deep groundwater flow within bedrock remained in a steady upwelling state, enhanced by increasing moisture in shallow soils under snow cover, reflecting the convergence of subsurface water from surrounding hillslopes.     

Influence of soil moisture hysteresis on the functioning of capillary barriers

Previous experimental studies of capillary barriers have identified highly hysteretic soil moisture retention characteristics in the materials used. In this study, numerical modelling is used to analyse the role of soil moisture hysteresis in capillary barrier functioning. Comparisons between published experimental results and model simulations indicate that soil moisture hysteresis was a necessary inclusion in the modelling approach to adequately reproduce pore water pressure distributions and the timing of breakthrough occurrences. Under hypothetical intermittent infiltration and evaporation conditions, the predicted volumetric water content in the moisture retention layer was significantly different for hysteretic and non‐hysteretic models. The hysteresis effect was found to be dependent on the nature of infiltration–evaporation cycling, although the predicted volume of flow through the hysteretic barrier was lower than that of the non‐hysteretic case, regardless of the nature of the cyclic upper boundary conditions. For practical engineering designs, where the water leakage through the barrier is the primary concern, the inclusion of soil moisture hysteresis in numerical modelling is needed. Copyright © 2009 John Wiley & Sons, Ltd.     

Soil water content in southern England derived from a cosmic‐ray soil moisture observing system – COSMOS‐UK

Cosmic‐ray soil moisture sensors have the advantage of a large measurement footprint (approximately 700 m in diameter) and are able to operate continuously to provide area‐averaged near‐surface (top 10–20 cm) volumetric soil moisture content at the field scale. This paper presents the application of this technique at four sites in southern England over almost 3 years. Results show the soil moisture response to contrasting climatic conditions during 2011–2014 and are the first such field‐scale measurements made in the UK. These four sites are prototype stations for a UK COsmic‐ray Soil Moisture Observing System, and particular consideration is given to sensor operating conditions in the UK. Comparison of these soil water content observations with the Joint UK Land Environment Simulator 10‐cm soil moisture layer shows that these data can be used to test and diagnose model performance and indicate the potential for assimilation of these data into hydro‐meteorological models. The application of these large‐area soil water content measurements to evaluate remotely sensed soil moisture products is also demonstrated. Numerous applications and the future development of a national COsmic‐ray Soil Moisture Observing System network are discussed. Copyright © 2016 John Wiley & Sons, Ltd.     

An analysis of infiltration with moisture content distribution in a two‐dimensional discretized water content domain

On the basis of unsaturated Darcy's law, the Talbot–Ogden method provides a fast unconditional mass conservative algorithm to simulate groundwater infiltration in various unsaturated soil textures. Unlike advanced reservoir modelling methods that compute unsaturated flow in space, it only discretizes the moisture content domain into a suitable number of bins so that the vertical water movement is estimated piecewise in each bin. The dimensionality of the moisture content domain is extended from one dimensional to two dimensional in this study, which allows us to distinguish pore shapes within the same moisture content range. The vertical movement of water in the extended model imitates the infiltration phase in the Talbot–Ogden method. However, the difference in this extension is the directional redistribution, which represents the horizontal inter‐bin flow and causes the water content distribution to have an effect on infiltration. Using this extension, we mathematically analyse the general relationship between infiltration and the moisture content distribution associated with wetting front depths in different bins. We show that a more negatively skewed moisture content distribution can produce a longer ponding time, whereas a higher overall flux cannot be guaranteed in this situation. It is proven on the basis of the water content probability distribution independent of soil textures. To illustrate this analysis, we also present numerical examples for both fine and coarse soil textures. Copyright © 2014 John Wiley & Sons, Ltd.     

降雨条件下边坡体积含水率变化规律研究

持续降雨是边坡发生失稳破坏的主要诱因之一,基于饱和—非饱和渗流理论,对梅州市大埔县某边坡的渗流场进行模拟,研究在不同降雨工况下该边坡土体体积含水率的时空变化规律。研究结果表明:相同条件下,降雨强度越大(降雨历时越长),边坡表层土体体积含水率变化越大;降雨强度60 mm/d历时1 d的暴雨对边坡表层土体体积含水率的增幅作用存在着一定的滞后性,其余工况未表现出滞后现象;降雨强度为120mm/d和300 mm/d的两种工况各研究点任意时段体积含水率较为接近;当降雨强度达到60 mm/d以上时,边坡内部体积含水率空间变化主要受降雨历时影响,降雨历时越长,降雨入渗深度和体积含水率变化越大。     

The influence of vegetation cover on evapotranspiration atop waste rock piles,Elk Valley,British Columbia

Surface mining in the Elk Valley, British Columbia, involves removing vegetation, soil, and rock to access underlying metallurgical coal. Subsequent waste rock is placed into adjacent valleys, frequently burying headwater streams. Due to their coarse texture, waste rock piles increase infiltration and percolation, increasing solute transport and concentration of geochemicals in downstream surface waters. Previous research suggests that weathering solutes are transport limited, and it is hypothesized that revegetation will enhance evapotranspiration (ET) and reduce percolation through the waste rock, potentially reducing loading. This study examined the surface‐atmosphere water and energy exchanges using the eddy covariance technique for three waste rock surfaces with different levels of reclamation: (a) an ~25‐year‐old mixed coniferous forest, (b) a grass site, and (c) bare waste rock. Measurements were taken from May to October in 2013 and 2014. Soil moisture and matric suction were measured to 1‐m depth. Sap flow at the forested site was measured to partition transpiration from total ET. In all years, ET rates were greatest at the forested site, followed by the grass cover and lowest at the bare waste rock site. Growing season ET rates at the forest were 56% higher than grass in 2013 and 35% higher in 2014. At the vegetated sites, climate was the main driver of ET, with high radiation, and warm and dry conditions enhancing fluxes. Maximum ET at these sites corresponded with peak growing season, with vegetation increasing both transpiration and rainfall interception. At the bare rock site, ET was weakly related to atmospheric conditions, and ET rates briefly increased during periods following rainfall when near‐surface soil moisture was enhanced. Transpiration comprised 29% of overall ET at the forest site from late July to early October. Results suggest that vegetation establishment can be incorporated into mine reclamation plans to enhance ET rates and limit percolation, potentially reducing downstream geochemical loads.     

Infiltration into frozen soil: From core‐scale dynamics to hillslope‐scale connectivity

Infiltration into frozen soil is a key hydrological process in cold regions. Although the mechanisms behind point‐scale infiltration into frozen soil are relatively well understood, questions remain about upscaling point‐scale results to estimate hillslope‐scale run‐off generation. Here, we tackle this question by combining laboratory, field, and modelling experiments. Six large (0.30‐m diameter by 0.35‐m deep) soil cores were extracted from an experimental hillslope on the Canadian Prairies. In the laboratory, we measured run‐off and infiltration rates of the cores for two antecedent moisture conditions under snowmelt rates and diurnal freeze–thaw conditions observed on the same hillslope. We combined the infiltration data with spatially variable data from the hillslope, to parameterise a surface run‐off redistribution model. We used the model to determine how spatial patterns of soil water content, snowpack water equivalent (SWE), and snowmelt rates affect the spatial variability of infiltration and hydrological connectivity over frozen soil. Our experiments showed that antecedent moisture conditions of the frozen soil affected infiltration rates by limiting the initial soil storage capacity and infiltration front penetration depth. However, shallow depths of infiltration and refreezing created saturated conditions at the surface for dry and wet antecedent conditions, resulting in similar final infiltration rates (0.3 mm hr^?1). On the hillslope‐scale, the spatial variability of snowmelt rates controlled the development of hydrological connectivity during the 2014 spring melt, whereas SWE and antecedent soil moisture were unimportant. Geostatistical analysis showed that this was because SWE variability and antecedent moisture variability occurred at distances shorter than that of topographic variability, whereas melt variability occurred at distances longer than that of topographic variability. The importance of spatial controls will shift for differing locations and winter conditions. Overall, our results suggest that run‐off connectivity is determined by (a) a pre‐fill phase, during which a thin surface soil layer wets up, refreezes, and saturates, before infiltration excess run‐off is generated and (b) a subsequent fill‐and‐spill phase on the surface that drives hillslope‐scale run‐off.     

A land surface soil moisture data assimilation framework in consideration of the model subgrid-scale heterogeneity and soil water thawing and freezing

The Ensemble Kalman Filter (EnKF) is well known and widely used in land data assimilation for its high precision and simple operation. The land surface models used as the forecast operator in a land data assimilation system are usually designed to consider the model subgrid-heterogeneity and soil water thawing and freezing. To neglect their effects could lead to some errors in soil moisture assimilation. The dual EnKF method is employed in soil moisture data assimilation to build a soil moisture data as- similation framework based on the NCAR Community Land Model version 2.0 (CLM 2.0) in considera- tion of the effects of the model subgrid-heterogeneity and soil water thawing and freezing: Liquid volumetric soil moisture content in a given fraction is assimilated through the state filter process, while solid volumetric soil moisture content in the same fraction and solid/liquid volumetric soil moisture in the other fractions are optimized by the parameter filter. Preliminary experiments show that this dual EnKF-based assimilation framework can assimilate soil moisture more effectively and precisely than the usual EnKF-based assimilation framework without considering the model subgrid-scale heteroge- neity and soil water thawing and freezing. With the improvement of soil moisture simulation, the soil temperature-simulated precision can be also improved to some extent.     

The effect of biological soil crusts on soil moisture dynamics under different rainfall conditions in the Tengger Desert,China

Rainfall is considered as the dominant water replenishment in desert ecosystems, and the conversion of rainfall into soil water availability plays a central role in sustaining the ecosystem function. In this study, the role of biological soil crusts (BSCs), typically formed in the revegetated desert ecosystem in the Tengger Desert of China, in converting rainfall into soil water, especially for the underlying soil moisture dynamics, was clarified by taking into account the synthetic effects of BSCs, rainfall characteristics, and antecedent soil water content on natural rainfall conditions at point scale. Our results showed that BSCs retard the infiltration process due to its higher water holding capacity during the initial stage of infiltration, such negative effect could be offset by the initial wet condition of BSCs. The influence of BSCs on infiltration amount was dependent on rainfall regime and soil depth. BSCs promoted a higher infiltration through the way of prolonged water containing duration in the ground surface and exhibited a lower infiltration at deep soil layer, which were much more obvious under small and medium rainfall events for the BSCs area compared with the sand area. Generally, the higher infiltration at top soil layer only increased soil moisture at 0.03 m depth; in consequence, there was no water recharge for the deep soil, and thus, BSCs had a negative effect on soil water effectiveness, which may be a potential challenge for the sustainability of the local deep‐rooted vegetation under the site specific rainfall conditions in northwestern China.     

Effects of antecedent moisture and macroporosity on infiltration and water flow in frozen soil

Infiltration into frozen soil plays an important role in soil freeze–thaw and snowmelt-driven hydrological processes. To better understand the complex thermal energy and water transport mechanisms involved, the influence of antecedent moisture content and macroporosity on infiltration into frozen soil was investigated. Ponded infiltration experiments on frozen macroporous and non-macroporous soil columns revealed that dry macroporous soil produced infiltration rates reaching 10³ to 10⁴ mm day⁻¹, two to three orders of magnitude larger than dry non-macroporous soil. Results suggest that rapid infiltration and drainage were a result of preferential flow through initially air-filled macropores. Using recorded flow rates and measured macropore characteristics, calculations indicated that a combination of both saturated flow and unsaturated film flow likely occurred within macropores. Under wet conditions, regardless of the presence of macropores, infiltration was restricted by the slow thawing rate of pore ice, producing infiltration rates of 2.8 to 5.0 mm day⁻¹. Reduced preferential flow under wet conditions was attributed to a combination of soil swelling, due to smectite-rich clay (that reduced macropore volume), and pore ice blockage within macropores. In comparison, dry soil column experiments demonstrated that macropores provided conduits for water and thermal energy to bypass the frozen matrix during infiltration, reducing thaw rates compared with non-macroporous soils. Overall, results showed the dominant control of antecedent moisture content on the initiation, timing, and magnitude of infiltration and flow in frozen macroporous soils, as well as the important role of macropore connectivity. The study provides an important data set that can aid the development of hydrological models that consider the interacting effects of soil freeze–thaw and preferential flow on snowmelt partitioning in cold regions.     

Characterising deep vadose zone water movement and solute transport under typical irrigated cropland in the North China Plain

Understanding the dynamics and mechanisms of soil water movement and solute transport is essential for accurately estimating recharge rates and evaluating the impacts of agricultural activities on groundwater resources. In a thick vadose zone (0–15 m) under irrigated cropland in the piedmont region of the North China Plain, soil water content, matric potential, and solute concentrations were measured. Based on these data, the dynamics of soil water and solutes were analysed to investigate the mechanisms of soil water and solute transport. The study showed that the 0–15‐m vadose zone can be divided into three layers: an infiltration and evaporation layer (0–2 m), an unsteady infiltration layer (2–6 m), and a quasi‐steady infiltration layer (6–15 m). The chloride, nitrate, and sulphate concentrations all showed greater variations in the upper soil layer (0–1 m) compared to values in the deep vadose zone (below 2 m). The average concentrations of these three anions in the deep vadose zone varied insignificantly with depth and approached values of 125, 242, and 116 mg/L. The accumulated chloride, sulphate, and nitrate were 2,179 ± 113, 1,760 ± 383, and 4,074 ± 421 kg/ha, respectively. The soil water potential and solute concentrations indicated that uniform flow and preferential flow both occurred in the deep vadose zone, and uniform flow was the dominant mechanism of soil water movement in this study. The piston‐like flow velocity of solute transport was 1.14 m per year, and the average value of calculated leached nitrate nitrogen was 107 kg/ha?year below the root zone. The results can be used to better understand recharge processes and improve groundwater resources management.     

Soil moisture variability of root zone profiles within SMEX02 remote sensing footprints

Remote sensing of soil moisture effectively provides soil moisture at a large scale, but does not explain highly heterogeneous soil moisture characteristics within remote sensing footprints. In this study, field scale spatio-temporal variability of root zone soil moisture was analyzed. During the Soil Moisture Experiment 2002 (SMEX02), daily soil moisture profiles (i.e., 0–6, 5–11, 15–21, and 25–31 cm) were measured in two fields in Walnut Creek watershed, Ames, Iowa, USA. Theta probe measurements of the volumetric soil moisture profile data were used to analyze statistical moments and time stability and to validate soil moisture predicted by a simple physical model simulation. For all depths, the coefficient of variation of soil moisture is well explained by the mean soil moisture using an exponential relationship. The simple model simulated very similar variability patterns as those observed.As soil depth increases, soil moisture distributions shift from skewed to normal patterns. At the surface depth, the soil moisture during dry down is log-normally distributed, while the soil moisture is normally distributed after rainfall. At all depths below the surface, the normal distribution captures the soil moisture variability for all conditions. Time stability analyses show that spatial patterns of sampling points are preserved for all depths and that time stability of surface measurements is a good indicator of subsurface time stability. The most time stable sampling sites estimate the field average root zone soil moisture value within ±2.1% volumetric soil moisture.     

STUDY OF INFILTRATION INTO A HETEROGENEOUS SOIL USING MAGNETIC RESONANCE IMAGING

In this paper, the feasibility of using magnetic resonance imaging (MRI) to study water infiltration into a heterogeneous soil is examined, together with its difficulties and limitations. MRI studies of ponded water infiltration into an undisturbed soil core show that the combination of one- and two-dimensional imaging techniques provides a visual and non-destructive means of monitoring the temporal changes of soil water content and the moisture profile, and the movement of the wetting front. Two-dimensional images show air entrapment in repetitive ponded infiltration experiments. During the early stages of infiltration, one-dimensional images of soil moisture profiles clearly indicate preferential flow phenomena. The observed advance of wetting fronts can be described by a linear relationship between the square root of infiltration time (√t) and the distance of the wetting front from the soil surface. Similarly, the cumulative infiltration is also directly proportional to √t. Furthermore, from the MRI infiltration moisture profiles, it is possible to estimate the parameters that feature in infiltration equations. © 1997 by John Wiley & Sons, Ltd.     

On the estimation of antecedent wetness conditions in rainfall–runoff modelling

The analysis of the physical processes involved in a conceptual model of soil water content balance is addressed with the objective of its application as a component of rainfall–runoff modelling. The model uses routinely measured meteorological variables (rainfall and air temperature) and incorporates a limited number of significant parameters. Its performance in estimating the soil moisture temporal pattern was tested through local measurements of volumetric water content carried out continuously on an experimental plot located in central Italy. The analysis was carried out for different periods in order to test both the representation of infiltration at the short time‐scale and drainage and evapotranspiration processes at the long time‐scale. A robust conceptual model was identified that incorporated the Green–Ampt approach for infiltration and a gravity‐driven approximation for drainage. A sensitivity analysis was performed for the selected model to assess the model robustness and to identify the more significant parameters involved in the principal processes that control the soil moisture temporal pattern. The usefulness of the selected model was tested for the estimation of the initial wetness conditions for rainfall–runoff modelling at the catchment scale. Specifically, the runoff characteristics (runoff depth and peak discharge) were found to be dependent on the pre‐event surface soil moisture. Both observed values and those estimated by the model gave good results. On the contrary, with the antecedent wetness conditions furnished by two versions of the antecedent precipitation index (API), large errors were obtained. Copyright © 2007 John Wiley & Sons, Ltd.     

Spatial variation of effective porosity and its implications for discharge in an upland headwater catchment in Scotland

Spatial and temporal measurements of shallow sub-surface soil physical properties were made within a 1 km² upland catchment. The surface soil layer of the catchment was organic rich (>70% organic matter) with a corresponding total porosity of 81%. Monthly point observations of volumetric water content (θ) were combined with point estimates of total porosity () and the porosity <50 μm (_residual), to define the ratio of water filled pore volume:pore volume in pores <50 μm (=θ/_residual). Values of θ/_residual were compared with discharge to test whether mass flow occurred when θ/_residual>1. A correlation between water content and discharge was found, with discharge increasing rapidly when θ/_residual approached unity. Similar relationships between water content and catchment discharge were identified for soil units adjacent to the stream when θ/_residual approached unity. These data suggest that soil pores >50 μm are of crucial importance in determining catchment discharge. Spatial and temporal variations in soil properties related to moisture content of the soil were also observed. Under dry conditions, a clear division based on aspect was noted, the west-facing side of the catchment being wettest. In wetter months, total porosity and soil water content were significantly affected by soil type and the spatial pattern of soil water content was more variable than in the dryer months. The physical quantification of soil properties in the shallow sub-surface layer proved important in explaining different initial changes in discharge from the catchment in response to a rainfall event.     

Effects of differing coverage of moss‐dominated soil crusts on hydrological processes and implications for disturbance in the Mu Us Sandland,China

To study the effects of biological soil crusts (BSCs) on hydrological processes and their implications for disturbance in the Mu Us Sandland, the water infiltration, evaporation and soil moisture of high coverage (100% BSCs), middle coverage (40% BSCs) and low coverage (0% BSCs, bare sand) of moss‐dominated crusts were conducted in this study, respectively. The conclusions are as follows: (1) the main effects of moss‐dominated crusts in the Mu Us Sandland on the infiltration of rainwater were to reduce the infiltration depths and to retain the limited rainwater in shallow soil; (2) moss‐dominated crusts have no significant effects on daily evaporation when the volumetric water content at 4 cm depth in 100% BSCs (VWC4) was over 24.7%, on enhanced daily evaporation when the VWC4 ranged from 6.5% to 24.7% and on reduced daily evaporation when the VWC4 was less than 6.5%; and (3) decreasing the coverage of moss‐dominated crusts (from 100% to 40%) did not significantly change its effects on infiltration, evaporation and soil moisture. Our results demonstrated that for the growth and regeneration of shrubs, which were dominated by Artemisia ordosica in the Mu Us Sandland, high coverage of moss‐dominated crusts has negative effects on hydrological processes, and these negative effects could not be significantly reduced by decreasing the coverage of moss‐dominated crusts from 100% to 40%. Therefore, for the sustained and healthy development of shrub communities in the Mu Us Sandland, it is necessary to take appropriate measures for the well‐developed BSCs in the sites with high vegetation coverage in the rainy season. Copyright © 2015 John Wiley & Sons, Ltd.     

Impact of subsurface rock fragments on runoff and interrill soil loss from cultivated soils

Field and laboratory studies have indicated that rock fragments in the topsoil may have a large impact on soil properties, soil quality, hydraulic, hydrological and erosion processes. In most studies, the rock fragments investigated still remain visible at the soil surface and only properties of these visible rock fragments are used for predicting runoff and soil loss. However, there are indications that rock fragments completely incorporated in the topsoil could also significantly influence the percolation and water distribution in stony soils and therefore, also infiltration, runoff and soil loss rates. Therefore, in this study interrill laboratory experiments with simulated rainfall for 60 min were conducted to assess the influence of subsurface rock fragments incorporated in a disturbed silt loam soil at different depths below the soil surface (i.e. 0.001, 0.01, 0.05 and 0.10 m), on infiltration, surface runoff and interrill erosion processes for small and large rock fragment sizes (i.e. mean diameter 0.04 and 0.20 m, respectively). Although only small differences in infiltration rate and runoff volume are observed between the soil without rock fragments (control) and the one with subsurface rock fragments, considerable differences in total interrill soil loss are observed between the control treatment and both contrasting rock fragments sizes. This is explained by a rapid increase in soil moisture in the areas above the rock fragments and therefore a decrease in topsoil cohesion compared with the control soil profile. The observed differences in runoff volume and interrill soil loss between the control plots and those with subsurface rock fragments is largest after a cumulative rainfall (P_cum) of 11 mm and progressively decreases with increasing P_cum. The results highlight the impacts and complexity of subsurface rock fragments on the production of runoff volume and soil loss and requires their inclusion in process‐based runoff and erosion models. Copyright © 2011 John Wiley & Sons, Ltd.