Experimental and analytical studies of solute transport during run‐off over vegetated surfaces

Solute transport in overland flow is considered as one of the main contributors to water pollution. Although many models of pollutant transport mechanism from soil to run‐off water have been proposed, the characteristics of solute transport accompanying the water run‐off over vegetated surface have not been well studied. In this study, a series of laboratory experiments were conducted to study the solute transport over vegetated surfaces. Based on the experimental results, an idea of the “stationary water layer” in run‐off was proposed. Applying the complete mixing theory in the stationary water layer, an analytical solute transport model was developed with the assumption that the upper run‐off completely mixes with the underlying water in the stationary water layer for each site. The results show that the predictions made by the present model are in good agreement with the measured experimental data. For the vegetated surfaces, the depth of stationary water layer is related to the rainfall intensity, bed slope, and vegetation density. The analytical solution shows that the maximum solute transport occurs at the time of concentration. This study advances our understanding of the mechanisms of solute transport over vegetated areas.     

Hillslope run‐off thresholds with shrink–swell clay soils

Irrigation experiments on 12 instrumented field plots were used to assess the impact of dynamic soil crack networks on infiltration and run‐off. During applications of intensity similar to a heavy rainstorm, water was seen being preferentially delivered within the soil profile. However, run‐off was not observed until soil water content of the profile reached field capacity, and the apertures of surface‐connected cracks had closed >60%. Electrical resistivity measurements suggested that subsurface cracks persisted and enhanced lateral transport, even in wet conditions. Likewise, single‐ring infiltration measurements taken before and after irrigation indicated that infiltration remained an important component of the water budget at high soil water content values, despite apparent surface sealing. Overall, although the wetting and sealing of the soil profile showed considerable complexity, an emergent property at the hillslope scale was observed: all of the plots demonstrated a strikingly similar threshold run‐off response to the cumulative precipitation amount. Copyright © 2014 John Wiley & Sons, Ltd.     

A finite element model for simulating surface run‐off and unsaturated seepage flow in the shallow subsurface

This work introduces water–air two‐phase flow into integrated surface–subsurface flow by simulating rainfall infiltration and run‐off production on a soil slope with the finite element method. The numerical model is formulated by partial differential equations for hydrostatic shallow flow and water–air two‐phase flow in the shallow subsurface. Finite element computing formats and solution strategies are presented to obtain a numerical solution for the coupled model. An unsaturated seepage flow process is first simulated by water–air two‐phase flow under the atmospheric pressure boundary condition to obtain the rainfall infiltration rate. Then, the rainfall infiltration rate is used as an input parameter to solve the surface run‐off equations and determine the value of the surface run‐off depth. In the next iteration, the pressure boundary condition of unsaturated seepage flow is adjusted by the surface run‐off depth. The coupling process is achieved by updating the rainfall infiltration rate and surface run‐off depth sequentially until the convergence criteria are reached in a time step. A well‐conducted surface run‐off experiment and traditional surface–subsurface model are used to validate the new model. Comparisons with the traditional surface–subsurface model show that the initiation time of surface run‐off calculated by the proposed model is earlier and that the water depth is larger, thus providing values that are closer to the experimental results.     

Simulating saturation‐excess surface run‐off in a semi‐distributed hydrological model

The level of complexity, and the number of parameters, to include in a hydrological model is a relatively contentious issue in hydrological modelling. However, it can be argued that explicitly representing important run‐off generation processes can improve the practical value of a model's outputs. This paper explores the benefits of including a new function into an existing semi‐distributed hydrological model (the Pitman model) that is widely used in the sub‐Saharan Africa region. The new function was designed to represent saturation‐excess surface run‐off processes at subcatchment scales and was motivated by the evidence of dambo (low topography riparian areas) type features in many sub‐Saharan river basins. The results for uncertainty versions of the model, with and without the new function, were compared for 25 catchments, which were divided up into those where evidence of dambos exists and those where there is no such evidence. The results suggest that the new function certainly improves the model results for the catchments where dambos exist, but not in situations where saturation‐excess surface run‐off is not expected to occur. The overall conclusion is therefore that the addition of the new function is justified.     

Identifying run‐off contributions during melt‐induced run‐off events in a glacierized alpine catchment

We analysed contributions to run‐off using hourly stream water samples from seven individual melt‐induced run‐off events (plus one rainfall event) during 2011, 2012 and 2013 in two nested glacierized catchments in the Eastern Italian Alps. Electrical conductivity and stable isotopes of water were used for mixing analysis and two‐component and three‐component hydrograph separation. High‐elevation snowmelt, glacier melt and autumn groundwater were identified as major end‐members. Discharge and tracers in the stream followed the diurnal variations of air temperature but markedly reacted to rainfall inputs. Hysteresis patterns between discharge and electrical conductivity during the melt‐induced run‐off events revealed contrasting loop directions at the two monitored stream sections. Snowmelt contribution to run‐off was highest in June and July (up to 33%), whereas the maximum contribution of glacier melt was reached in August (up to 65%). The maximum groundwater and rainfall contributions were 62% and 11%, respectively. Run‐off events were generally characterized by decreasing snowmelt and increasing glacier melt fractions from the beginning to the end of the summer 2012, while run‐off events in 2013 showed less variable snowmelt and lower glacier melt contributions than in 2012. The results provided essential insights into the complex dynamics of melt‐induced run‐off events and may be of further use in the context of water resource management in alpine catchments. Copyright © 2015 John Wiley & Sons, Ltd.     

Modelling soil solute release and transport in run‐off on a loessial slope with and without surface stones

Stone covers on loessial slopes can increase the time of infiltration by slowing the velocity of the overland flow, which reduces the transport of solutes, but few mechanistic models have been tested under water‐scouring conditions. We carried out field experiments to test a previously proposed, physically based model of water and solute transport. The area of soil infiltration was calculated from the uncovered surface area, and Richards' equation and the kinematic wave equation were used to describe water infiltration and flow along slopes with stone covers. The transport of chemicals into the run‐off from the surface soil, presumably by diffusion, and their movement in the soil profile could be described by the convection–diffusion equations of the model. The simulated and measured data correlated well. The stones on the soil surface reduced the area available for infiltration but increased the Manning coefficient, eventually leading to increased water infiltration and decreased solute loss with run‐off. Our results indicated that the traditional model of water movement and solute migration could be used to simulate water transport and solute migration for stone‐covered soil on loessial slopes.     

Run‐off processes from mountains to foothills: The role of soil stratigraphy and structure in influencing run‐off characteristics across high to low relief landscapes

The critical zone features that control run‐off generation, specifically at the regional watershed scale, are not well understood. Here, we addressed this knowledge gap by quantitatively and conceptually linking regional watershed‐scale run‐off regimes with critical zone structure and climate gradients across two physiographic provinces in the Southeastern United States. We characterized long‐term (~20 years) discharge and precipitation regimes for 73 watersheds with United States Geological Survey in‐stream gaging stations across the Appalachian Mountain and Piedmont physiographic provinces of North Carolina. Watersheds included in this analysis had <10% developed land and ranged in size from 14.1–4,390 km². Thirty‐four watersheds were located in the Piedmont physiographic province, which is typically classified as a low relief landscape with deep, highly weathered soils and regolith. Thirty‐nine watersheds were located in the Appalachian Mountain physiographic province, which is typically classified as a steeper landscape with highly weathered, but shallower soils and regolith. From the United States Geological Survey daily mean run‐off time series, we calculated annual and seasonal baseflow indices (BFI), minimum, mean, and maximum daily run‐off, and Pearson's correlation coefficients between precipitation and baseflow. Our results showed that Appalachian Mountain watersheds systematically had higher minimum daily flows and BFI values. Piedmont watersheds displayed much larger deviations from mean annual BFI in response to year‐to‐year variability in precipitation. A series of linear regression models between 21 landscape metrics and annual BFIs showed non‐linear and complex terrestrial–hydrological relationships across the two provinces. From these results, we discuss how distinct features of critical zone architecture, with specific focus on soil depth and stratigraphy, may be dominating the regulation of hydrological processes and run‐off regimes across these provinces.     

Thresholds for run‐off generation in a drained closed depression

Hydrological threshold behaviour has been observed across hillslopes and catchments with varying characteristics. Few studies, however, have evaluated rainfall–run‐off response in areas dominated by agricultural land use and artificial subsurface drainage. Hydrograph analysis was used to identify distinct hydrological events over a 9‐year period and examine rainfall characteristics, dynamic water storage, and surface and subsurface run‐off generation in a drained and farmed closed depression in north‐eastern Indiana, USA. Results showed that both surface flow and subsurface tile flow displayed a threshold relationship with the sum of rainfall amount and soil moisture deficit (SMD). Neither surface flow nor subsurface tile flow was observed unless rainfall amount exceeded the SMD. Timing of subsurface tile flow relative to soil moisture response on the shoulder slope of the depression indicated that the formation and drainage of perched water tables on depression hillslopes were likely the main mechanism that produced subsurface connectivity. Surface flow generation was delayed compared with subsurface tile flow during rainfall events due to differences in soil water storage along depression hillslopes and run‐off generation mechanisms. These findings highlight the substantial impact of subsurface tile drainage on the hydrology of closed depressions; the bottom of the depression, the wettest area prior to drainage installation, becomes the driest part of the depression after installation of subsurface drainage. Rapid connectivity of localized subsurface saturation zones during rainfall events is also greatly enhanced because of subsurface drainage. Thus, less fill is required to generate substantial spill. Understanding hydrologic processes in drained and farmed closed depressions is a critical first step in developing improved water and nutrient management strategies in this landscape.     

Surface moisture area during rainfall–run‐off events to understand the hydrological dynamics of a basin in a plain region

The understanding of the hydrology of plain basins may be improved by the combined analysis of rainfall–run‐off records and remote sensed surface moisture data. Our work evaluates the surface moisture area (SMA) produced during rainfall–run‐off events in a plain watershed of the Argentine Pampas Region, and studies which hydrological variables are related to the generated SMA. The study area is located in the upper and middle basins of the Del Azul stream, characterized by the presence of small gently hilly areas surrounded by flat landscapes. Data from 9 rainfall–run‐off events were analysed. MODIS surface reflectance data were processed to calculate SMA subsequent to the peak discharge (post‐SMA), and previous to the rainfall events (prev‐SMA), to consider the antecedent wetness. Rainfall–run‐off data included total precipitation depth (P), maximum intensity of rainfall over 6 hr (I6max), surface run‐off registered between the beginning of the event and the day previous to the analysed MODIS scene (R), peak flow (Qp), and flood intensity (IF). In contrast with other works, post‐SMA showed a negative relationship with the R. Three groups of cases were identified: (a) Events of low I6max, high prev‐SMA, and low R were associated with slow and weakly channelized flow over plain areas, leading to saturated overland flow (SOF), with large SMA; (b) events of high I6max, low prev‐SMA, and medium to high R were rapidly transported along the gentle slopes of the basin, related to Hortonian overland flow (HOF) and low post‐SMA; and (c) events of medium to high I6max and prev‐SMA with medium R were related to heterogeneous input‐antecedent‐run‐off conditions combined: Local spatial conditions may have produced HOF or SOF, leading to an averaged response with medium SMA. The interactions between the geomorphology of the basin, the characteristics of the events, and the antecedent conditions may explain the obtained results. This analysis is relevant for the general knowledge of the hydrology of large plains, whose functioning studies are still in their early stages.     

Tracing variability of run‐off generation in mountainous permafrost of semi‐arid north‐eastern Mongolia

The headwaters of mountainous, discontinuous permafrost regions in north‐eastern Mongolia are important water resources for the semi‐arid country, but little is known about hydrological processes there. Run‐off generation on south‐facing slopes, which are devoid of permafrost, has so far been neglected and is totally unknown for areas that have been affected by recent forest fires. To fill this knowledge gap, the present study applied artificial tracers on a steppe‐vegetated south‐facing and on two north‐facing slopes, burned and unburned. Combined sprinkling and dye tracer experiments were used to visualize processes of infiltration and water fluxes in the unsaturated zone. On the unburned north‐facing slope, rapid and widespread infiltration through a wet organic layer was observed down to the permafrost. On the burned profile, rapid infiltration occurred through a combusted organic and underlying mineral layer. Stained water seeped out at the bottom of both profiles suggesting a general tendency to subsurface stormflow (SSF). Ongoing SSF could directly be studied 24 h after a high‐intensity rainfall event on a 55‐m hillslope section in the burned forest. Measurements of water temperature proved the role of the permafrost layer as a base horizon for SSF. Repeated tracer injections allowed direct insights into SSF dynamics: A first injection suggested rather slow dispersive subsurface flow paths; whereas 18 h later, a second injection traced a more preferential flow system with 20 times quicker flow velocities. We speculate that these pronounced SSF dynamics are limited to burned slopes where a thermally insulating organic layer is absent. On three south‐facing soil profiles, the applied tracer remained in the uppermost 5 cm of a silt‐rich mineral soil horizon. No signs of preferential infiltration could be found, which suggested reduced biological activity under a harsh, dry and cold climate. Instead, direct observations, distributed tracers and charcoal samples provided evidence for the occurrence of overland flow. Copyright © 2014 John Wiley & Sons, Ltd.     

Flow threads in surface run‐off: implications for the assessment of flow properties and friction coefficients in soil erosion and hydraulics investigations

Soil surface microtopography produces non‐uniform surface run‐off, in which narrow threads of relatively deep and fast ?ow move within broader, shallower, slower‐moving regions. This kind of ?ow is probably widespread, given that microtopography is itself common. Methods used to record the properties of surface run‐off include grid‐ or transect‐based depth observations, with a single mean ?ow speed derived by calculation from V = Q/WD, and the use of dye timing to estimate velocity, with an effective mean depth calculated from D = Q/WV. Because these methods allow only single, ?ow‐?eld mean values to be derived for V or D, neither is well suited to non‐uniform ?ows. The use of depth data to derive a ?ow‐?eld mean V furthermore implicitly applies area weighting to the depth data; likewise, the use of dye speeds for V inherently overestimates mean V because dye dominantly follows the faster ?ow threads. The associated errors in derived parameters such as friction coef?cients are not readily quanti?ed and appear not to have been addressed previously. New ?eld experiments made on untilled soil surfaces in arid western NSW, Australia, explore these circumstances and the implications for deriving meaningful measures of ?ow properties, including friction coef?cients. On surfaces deliberately chosen for their very subtle microtopography, average thread velocities are shown to be commonly 2·5 times greater than the ?ow‐?eld mean, and locally 6–7 times greater. On the other hand, non‐thread ?ow speeds lie below the ?ow‐?eld mean, on average reaching only 84 per cent of this value, and often considerably less. Flow‐?eld means conceal the existence of regions of the ?ow ?eld whose properties are statistically distinct. Results con?rm that a reliance on ?ow‐?eld average depths yields estimates of friction coef?cients that are biased toward the shallower, high‐roughness parts of the ?ow, while if dye speeds are relied upon the results are biased toward the deeper, smoother threads of ?ow. A new approach to the evaluation of friction coef?cients in non‐uniform ?ows is advanced, involving the determination of separate coef?cients for threads and non‐thread zones of the ?ow ?eld. In contrast, ?ow‐?eld friction coef?cients as they are customarily derived in run‐off plot experiments subsume these distinct coef?cients in proportions that are generally unknown. The value of such coef?cients is therefore questionable. Copyright © 2004 John Wiley & Sons, Ltd.     

Water repellency as run‐off and soil detachment controlling factor in a dry‐Mediterranean hillslope (South of Spain)

This study examines the effect of water repellency on controlling temporal variability of runoff generation mechanisms and soil detachment on metamorphic derived soils under dry‐Mediterranean climate. The research is carried out in an unburnt Mediterranean hillslope in souther Spain characterized by a patchy vegetation pattern and shallow soils. The Water Drop Penetration Time test (WDPT) is applied to measure water repellency at the end of summer (Sep‐2008), mid autumn (Nov‐2008) and mid winter (Feb‐2009). Rainfall simulations were used to obtain runoff generation and soil detachment in the same periods of time. The main shrub specie is Cistus monspeliensis which leaves a load of litter during the summer due to the lack of water. This great amount of organic material is accumulated under the shrubs triggering an extreme water repellence (WDPT > 6,000 s) that limits infiltration processes. This process is enforced due to the low soil water content at the end of dry season. Certain water repellency (WDPT > 1,500 s) is also observed on bare soil as consequence of their sandier texture and the accumulation of annual plants which die at the end of the wet season. Soil moisture increases during the autumn and water repellency disappears in both shrub and bare soil at the middle of the wet season (WDPT < 5 s). The main consequence is that the temporal trend of water repellency controls the mechanism and frequency of runoff generation and, hence, soil detachment. At the end of the summer, Hortonian mechanisms predominates when water repellency is extreme, even in soils under Cistus monspeliensis where runoff generation can reach higher peaks of overland flow and sediment concentration. Conversely, only the saturation of soil could generate runoff during the wet season being this quite less frequent in bare soil and absent in shrub. Copyright © 2010 John Wiley & Sons, Ltd.     

Long‐term hydrological simulation based on the Soil Conservation Service curve number

Presenting a critical review of daily flow simulation models based on the Soil Conservation Service curve number (SCS‐CN), this paper introduces a more versatile model based on the modified SCS‐CN method, which specializes into seven cases. The proposed model was applied to the Hemavati watershed (area = 600 km²) in India and was found to yield satisfactory results in both calibration and validation. The model conserved monthly and annual runoff volumes satisfactorily. A sensitivity analysis of the model parameters was performed, including the effect of variation in storm duration. Finally, to investigate the model components, all seven variants of the modified version were tested for their suitability. Copyright © 2004 John Wiley & Sons, Ltd.     

Micro‐geomorphology determines community structure of biological soil crusts at small scales

Biological soil crusts (BSCs) are ubiquitous communities of diminutive organisms such as cyanobacteria, green algae, lichens, mosses and others associated closely with particles of surface soil, forming a cohesive thin horizontal layer. The ecological roles of BSCs affecting soil nutrient cycling, stability and hydrological processes, influencing the germination and establishment of vascular plants, and serving as habitats for numerous arthropods and microorganism have been well documented. We tested the hypothesis that micro‐geomorphological features determine the spatial distribution of BSCs by reallocating related abiotic resources at small‐ and medium‐scales in the Tengger Desert. Our results showed that higher soil pH and higher total potassium content in topsoil positively correlated with the colonization of cyanobacteria and algae in the earliest successional stages of BSCs, while increasing dust deposition onto the topsoil enhanced the development of lichen and mosses in the later stages of BSCs. Increasing soil moisture raised the proportion of mosses and lichen in BSCs, this will possibly change the ecological functions of BSCs, such as nitrogen‐fixation by cyanobacteria, due to the conversion from a complex to relative simple type of BSC. Micro‐geomorphology has created various habitats at a small‐scale affecting colonization and development of cryptogams. This paper considers the contribution of micro‐geomorphology to biodiversity in the extreme arid desert systems. Copyright © 2010 John Wiley & Sons, Ltd.     

Analytical and experimental study on dissolved pollutant wash‐off over impervious surfaces

Non‐point source pollution in the impervious surface of city, which including dissolved and particulate pollutants, is a significant source of water pollution. Simple first‐order decay models can generally simulate the cumulative wash‐off process of the particulate pollutants. There is inadequate knowledge as to whether or not they are suitable for dissolved pollutants. This study presents a mathematical wash‐off model for dissolved pollutants, which combines analytical equations for overland flows and the exponential equation for the pollutant wash‐off. A series of laboratory experiments have been conducted to verify this wash‐off model. It shows that the pollutant concentration and pollutant transport rate can be predicted well by the newly developed equations. It is found that the pollutant concentration monotonically decreases to zero as the accumulated pollutants are washed off, whereas the pollutant transport rate first increases to the maximum value and then decreases to zero. The maximum pollutant transport rate is found to increase with the decrease of the arrival time of the maximum value. The difference between the simplified exponential model and the amended wash‐off equation depends on the initial residual percentage (P_c), but the present equation generally provides a more accurate representation of the wash‐off process of dissolved pollutants.     

Development of a soil moisture‐based distributed hydrologic model for determining hydrologically based critical source areas

A simple grid cell‐based distributed hydrologic model was developed to provide spatial information on hydrologic components for determining hydrologically based critical source areas. The model represents the critical process (soil moisture variation) to run‐off generation accounting for both local and global water balance. In this way, it simulates both infiltration excess run‐off and saturation excess run‐off. The model was tested by multisite and multivariable evaluation on the 50‐km² Little River Experimental Watershed I in Georgia, U.S. and 2 smaller nested subwatersheds. Water balance, hydrograph, and soil moisture were simulated and compared to observed data. For streamflow calibration, the daily Nash‐Sutcliffe coefficient was 0.78 at the watershed outlet and 0.56 and 0.75 at the 2 nested subwatersheds. For the validation period, the Nash‐Sutcliffe coefficients were 0.79 at the watershed outlet and 0.85 and 0.83 at the 2 subwatersheds. The per cent bias was less than 15% for all sites. For soil moisture, the model also predicted the rising and declining trends at 4 of the 5 measurement sites. The spatial distribution of surface run‐off simulated by the model was mainly controlled by local characteristics (precipitation, soil properties, and land cover) on dry days and by global watershed characteristics (relative position within the watershed and hydrologic connectivity) on wet days when saturation excess run‐off was simulated. The spatial details of run‐off generation and travel time along flow paths provided by the model are helpful for watershed managers to further identify critical source areas of non‐point source pollution and develop best management practices.     

Understanding land use and cover change impacts on run‐off and sediment load at flood events on the Loess Plateau,China

The Loess Plateau has been experiencing large‐scale land use and cover changes (LUCCs) over the past 50 years. It is well known about the significant decreasing trend of annual streamflow and sediment load in the catchments in this area. However, how surface run‐off and sediment load behaved in response to LUCC at flood events remained a research question. We investigated 371 flood events from 1963 to 2011 in a typical medium‐sized catchment within the Plateau in order to understand how LUCC affected the surface run‐off generation and sediment load and their behaviours based on the analysis of return periods. The results showed that the mean annual surface run‐off and sediment load from flood events accounted for 49.6% and 91.8% of their mean annual totals. The reduction of surface run‐off and associated sediment yield in floods explained about 85.0% and 89.2% of declines in the total annual streamflow and sediment load, respectively. The occurrences of flood events and peak sediment concentrations greater than 500 kg/m³ showed a significantly downward trend, yet the counterclockwise loop events still dominated the flood event processes in the catchment. The results suggest that LUCC over the past 50 years resulted in significant changes in the water balance components and associated soil erosion and sediment transportation in the catchment. This was achieved mainly by reducing surface run‐off and sediment yield during floods with return period of less than 5 years. Run‐off–sediment load behaviour during the extreme events with greater than 10‐year return periods has not changed. Outcomes from this study are useful in understanding the eco‐hydrological processes and assisting the sustainable catchment management and land use planning on the Loess Plateau, and the methodologies are general and applicable to similar areas worldwide.     

Rainfall,run‐off,and sediment transport dynamics in a humid mountain badland area: Long‐term results from a small catchment

In this study, we investigated rainfall, run‐off, and sediment transport dynamics (414 run‐off events and 231 events with sediment information) of a humid mountain badland area—the Araguás catchment (Central Pyrenees, Spain)—from October 2005 to September 2016. Use of this long‐term database allows characterization of the hydrological response, which consist of low‐magnitude/high‐frequency events and high‐magnitude/low‐frequency events, and identification of seasonal dynamics and rainfall‐run‐off thresholds. Our results indicate that the Araguás catchment, similarly to other humid badlands, had high hydrological responsiveness (mean annual run‐off coefficient: 0.52), a non‐linear relationship of rainfall with run‐off (common in Mediterranean environments), and seasonal hydrological and sedimentological dynamics. We created and validated a multivariate regression model to characterize the hydrological variables (stormflow and peak discharge) and sedimentological variables (mean and maximum suspended sediment concentrations and total suspended sediment load). In summer and at the beginning of autumn, the response was mainly related to rainfall intensity, suggesting a predomination of Hortonian flows. In contrast, in spring and winter, the responses were mainly related to the antecedent conditions (previous rainfall and baseflow), suggesting the occurrence of saturated excess flow processes, and the contribution of neighbouring vegetated areas. The multivariate analysis also showed that total sediment load is better predicted by a multivariate regression model that integrates pre‐event, rainfall, and run‐off variables. In general, our models provided more accurate predictions of small‐magnitude/high‐frequency events than high‐magnitude/low‐frequency events. This study highlights the high inter‐ and intra‐annual variability response in humid badland areas and that long‐term records are needed to reduce the uncertainty of hydrological and sedimentological responses in Mediterranean badland areas.     

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.