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.     

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.     

Variable power‐law scaling of hillslope Hortonian rainfall–runoff processes

Hydrological studies focused on Hortonian rainfall–run‐off scaling have found that the run‐off depth generally declines with the plot length in power‐law scaling. Both the power‐law proportional coefficient and the scaling exponent show great variability for specific conditions, but why and how they vary remain unclear. In the present study, the scaling of hillslope Hortonian rainfall–run‐off processes is investigated for different rainfall, soil infiltration, and hillslope surface characteristics using the physically based cell‐based rainfall‐infiltration‐run‐off model. The results show that both temporally intermittent and steady rainfalls can result in prominent power‐law scaling at the initial stage of run‐off generation. Then, the magnitude of the power‐law scaling decreases gradually due to the decreasing run‐on effect. The power‐law scaling is most sensitive to the rainfall and soil infiltration parameters. When the ratio of rainfall to infiltration exceeds a critical value, the magnitude of the power‐law scaling tends to decrease notably. For different intermittent rainfall patterns, the power‐law exponent varies in the range of ?1.0 to ?0.113, which shows an approximately logarithmic increasing trend for the proportional coefficient as a function of the run‐off coefficient. The scaling is also sensitive to the surface roughness, soil sealing, slope angle, and hillslope geometry because these factors control the run‐off routing and run‐on infiltration processes. These results provide insights into the variable scaling of the Hortonian rainfall–run‐off process, which are expected to benefit modelling of large‐scale hydrological and ecological processes.     

Assessment of transmission loss in a Mediterranean karstic watershed (Wadi Natuf,West Bank)

Run‐off transmission loss into karstified consolidated aquifer bedrock below ephemeral streams (wadis) has rarely been described nor quantified. This study presents unique data of long‐term high‐resolution field measurements and field observations in a semiarid to subhumid Mediterranean carbonatic mountainshed. The catchment with a 103 km² surface area is subdivided into 5 subcatchments. Coupled run‐off measurements were made in the different stream sections (reaches), and transmission loss calculated from differences in discharge. Rainfall and run‐off observations from 9 automated precipitation gauging stations and 5 pressure transducers for automatic water level recording are complemented by manual measurements during 34 run‐off events covering a total measurement period of 8 consecutive years. Run‐off generation is strongly event based depending on rainfall intensities and depths. Both, run‐off generation and transmission losses are related to spatial patterns of bedrock lithologies (and hydrostratigraphy). Transmission losses range between 62% and 80% of generated run‐off, with most of the smaller events showing 100% transmission loss. Therefore, although event run‐off coefficients in the mountains can reach up to 22%, only 0.11% of total annual precipitation leaves the catchment as run‐off. Most run‐off infiltrates directly into the regional karst aquifers (Upper Cretaceous carbonate series), with transmission loss intensities of up to 40 mm/h below the stream channels. The factors determining run‐off—such as geology, pedology, vegetation cover and land use, relief and morphology, the semiarid to subhumid Mediterranean climate with a strong elevation gradient, and the patchiness of individual storm events distributed over the winter seasons—as well as the lithology and epikarst features of the bedrock are all characteristic for larger areas in the Mediterranean region. Therefore, we expect that our findings can be generalized to a large extent.     

Spatially distributed long‐term hydrologic simulation using a continuous SCS CN method‐based hybrid hydrologic model

A continuous Soil Conservation Service (SCS) curve number (CN) method that considers time‐varied SCS CN values was developed based on the original SCS CN method with a revised soil moisture accounting approach to estimate run‐off depth for long‐term discontinuous storm events. The method was applied to spatially distributed long‐term hydrologic simulation of rainfall‐run‐off flow with an underlying assumption for its spatial variability using a geographic information systems‐based spatially distributed Clark's unit hydrograph method (Distributed‐Clark; hybrid hydrologic model), which is a simple few parameter run‐off routing method for input of spatiotemporally varied run‐off depth, incorporating conditional unit hydrograph adoption for different run‐off precipitation depth‐based direct run‐off flow convolution. Case studies of spatially distributed long‐term (total of 6 years) hydrologic simulation for four river basins using daily NEXRAD quantitative precipitation estimations demonstrate overall performances of Nash–Sutcliffe efficiency (E_NS) 0.62, coefficient of determination (R²) 0.64, and percent bias 0.33% in direct run‐off and E_NS 0.71, R² 0.72, and percent bias 0.15% in total streamflow for model result comparison against observed streamflow. These results show better fit (improvement in E_NS of 42.0% and R² of 33.3% for total streamflow) than the same model using spatially averaged gauged rainfall. Incorporation of logic for conditional initial abstraction in a continuous SCS CN method, which can accommodate initial run‐off loss amounts based on previous rainfall, slightly enhances model simulation performance; both E_NS and R² increased by 1.4% for total streamflow in a 4‐year calibration period. A continuous SCS CN method‐based hybrid hydrologic model presented in this study is, therefore, potentially significant to improved implementation of long‐term hydrologic applications for spatially distributed rainfall‐run‐off generation and routing, as a relatively simple hydrologic modelling approach for the use of more reliable gridded types of quantitative precipitation estimations.     

Effective run‐off flow length over biological soil crusts on silty loam soils in drylands

Surface hydrological behaviour is important in drylands because it affects the distribution of soil moisture and vegetation and the hydrological functioning of slopes and catchments. Microplot scale run‐off can be relatively easily measured, i.e. by rainfall simulations. However, slope or catchment run‐off cannot be deduced from microplots, requiring long‐time monitoring, because run‐off coefficients decrease with increasing drainage area. Therefore, to determine the slope length covered by run‐off (run‐off length) is crucial to connect scales. Biological soil crusts (BSCs) are good model systems, and their hydrology at slope scale is insufficiently known. This study provides run‐off lengths from BSCs, by field factorial experiments using rainfall simulation, including two BSC types, three rain types, three antecedent soil moistures and four plot lengths. Data were analysed by generalized linear modelling, including vascular plant cover as covariates. Results were the following: (i) the real contributing area is almost always much smaller than the topographical contributing area; (ii) the BSC type is key to controlling run‐off; run‐off length reached 3 m on cyanobacterial crust, but hardly over 1 m on lichen crust; this pattern remained through rain type or soil moisture; (iii) run‐off decreased with BSC development because soil sealing disappears; porosity, biomass and roughness increase and some changes occur in the uppermost soil layer; and (iv) run‐off flow increased with both rain type and soil moisture but run‐off coefficient only with soil moisture (as larger rains increased both run‐off and infiltration); vascular plant cover had a slight effect on run‐off because it was low and random. Copyright © 2014 John Wiley & Sons, Ltd.     

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.     

Estimating daily inflows of large lakes using a water‐balance‐based runoff coefficient scaling approach

Daily river inflow time series are highly valuable for water resources and water environment management of large lakes. However, the availability of continuous inflow data for large lakes is still relatively limited, especially for large lakes situated within humid plain regions with tens or even hundreds of tributaries. In this study, we choose the fifth largest freshwater Lake Chaohu in China as our study area to introduce a new approach to reconstruct historical daily inflows at ungauged subcatchments of large lakes. This approach makes use of water level, lake surface rainfall, evaporation from the lake, and catchment rainfall observations. Rainfall–runoff relationship at a reference catchment was analysed to select rainfall input and estimate run‐off coefficient firstly, and the run‐off coefficient was then transferred to ungauged subcatchments to initially estimate daily inflows. Run‐off coefficient was scaled to adjust daily inflows at ungauged subcatchments according to water balance of the lake. This approach was evaluated using sparsely measured inflows at eight subcatchments of Lake Chaohu and compared with the commonly used drainage area ratio method. Results suggest that the inflow time series reconstructed from this approach consistent well to corresponding observations, with mean R² and Nash–Sutcliffe efficiency values of 0.69 and 0.6, respectively. This approach outperforms drainage area ratio method in terms of mean R² and Nash–Sutcliffe efficiency values. Accuracy of this approach holds well when the number of water‐level station being used decreased from four to one.     

Quantifying the role of the snowpack in generating water available for run‐off during rain‐on‐snow events from snow pillow records

Rain‐on‐snow events have generated major floods around the world, particularly in coastal, mountainous regions. Most previous studies focused on a limited number of major rain‐on‐snow events or were based primarily on model results, largely due to a lack of long‐term records from lysimeters or other instrumentation for quantifying event water balances. In this analysis, we used records from five automated snow pillow sites in south coastal British Columbia, Canada, to reconstruct event water balances for 286 rain‐on‐snow events over a 10‐year period. For large rain‐on‐snow events (event rainfall >40 mm), snowmelt enhanced the production of water available for run‐off (WAR) by approximately 25% over rainfall alone. For smaller events, a range of antecedent and meteorological factors influenced WAR generation, particularly the antecedent liquid water content of the snowpack. Most large events were associated with atmospheric rivers. Rainfall dominated WAR generation during autumn and winter events, whereas snowmelt dominated during spring and summer events. In the majority of events, the sensible heat of rain contributed less than 10% of the total energy consumed by snowmelt. This analysis illustrated the importance of understanding the amount of rainfall occurring at high elevations during rain‐on‐snow events in mountainous regions.     

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.     

Long‐term and seasonal hydrologic performance of an extensive green roof

Sustainable strategies such as green roofs have been implemented as stormwater management tools to mitigate disturbance of the hydrologic cycle resulting from urbanization. Green roofs, also referred to as vegetated roofs, can improve the urban landscape by reducing heat island effects, providing ecosystem services, and facilitating the retention and treatment of stormwater. Green roofs have received particular attention because they do not require acquisition and development of land and represent an application of biomimicry in design and construction. In this paper, we evaluate the effects of precipitation, evapotranspiration (ET), antecedent dry period (ADP), and seasonal variation on the run‐off quantity and distribution of an extensive, sedum covered, green roof on a commercial building in Syracuse, NY, USA. The green roof greatly facilitated retention of precipitation events without significant changes over the 4‐year study. The green roof retained on average 95.9 ± 3.6% (6.5 ± 5.6 mm) per rainfall event, with a range from 75% to 99.6% (33.2 to 3.3 mm). However, as precipitation quantity increased, the retention of water decreased. This high water retention capacity was the result of the combined effects of ET, stormwater storage (plants, growth media, and stormwater retention layer), and limited surface run‐off from the roof deck due to variation in the sloping of the green roof and the tapered insulation to the deck drains. The water retention capacity of the green roof did not change significantly between growing and nongrowing seasons. Slightly greater precipitation during the growing season coincided with increased ET. Average potential ET during the growing season was approximately 3 times greater than during the nongrowing season. The hydrologic performance of the green roof was not significantly impacted by an ADP greater than 2 days.     

An analytical approach to ascertain saturation‐excess versus infiltration‐excess overland flow in urban and reference landscapes

Uncontrolled overland flow drives flooding, erosion, and contaminant transport, with the severity of these outcomes often amplified in urban areas. In pervious media such as urban soils, overland flow is initiated via either infiltration‐excess (where precipitation rate exceeds infiltration capacity) or saturation‐excess (when precipitation volume exceeds soil profile storage) mechanisms. These processes call for different management strategies, making it important for municipalities to discern between them. In this study, we derived a generalized one‐dimensional model that distinguishes between infiltration‐excess overland flow (IEOF) and saturation‐excess overland flow (SEOF) using Green–Ampt infiltration concepts. Next, we applied this model to estimate overland flow generation from pervious areas in 11 U.S. cities. We used rainfall forcing that represented low‐ and high‐intensity events and compared responses among measured urban versus predevelopment reference soil hydraulic properties. The derivation showed that the propensity for IEOF versus SEOF is related to the equivalence between two nondimensional ratios: (a) precipitation rate to depth‐weighted hydraulic conductivity and (b) depth of soil profile restrictive layer to soil capillary potential. Across all cities, reference soil profiles were associated with greater IEOF for the high‐intensity set of storms, and urbanized soil profiles tended towards production of SEOF during the lower intensity set of storms. Urban soils produced more cumulative overland flow as a fraction of cumulative precipitation than did reference soils, particularly under conditions associated with SEOF. These results will assist cities in identifying the type and extent of interventions needed to manage storm water produced from pervious areas.     

Application of the WEPP model to determine sources of run‐off and sediment in a forested watershed

This study investigates critical run‐off and sediment production sources in a forested Kasilian watershed located in northern Iran. The Water Erosion Prediction Project (WEPP) watershed model was set up to simulate the run‐off and sediment yields. WEPP was calibrated and validated against measured rainfall–run‐off–sediment data. Results showed that simulated run‐off and sediment yields of the watershed were in agreement with the measured data for the calibration and validation periods. While low and medium values of run‐off and sediment yields were adequately simulated by the WEPP model, high run‐off and sediment yield values were underestimated. Performance of the model was evaluated as very good and satisfactory during the calibration and validation stages, respectively. Total soil erosion and sediment load of the study watershed during the study period were determined to be 10 108 t yr^?1 and 8735 t yr^?1, respectively. The northern areas of the watershed with dry farming were identified as the critical erosion prone zones. To prioritize the subwatersheds based on their contribution to the run‐off and sediment production at the watershed's main outlet, unit response approach (URA) was applied. In this regard, subwatersheds close to the main outlet were found to have the highest contribution to sediment yield of the whole watershed. Results indicated that depending on the objective of land and water conservation practices, particularly, for controlling sediment yield at the main outlet, critical areas for implementing the best management practices may be identified through conjunctive application of WEPP and URA. 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.     

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.     

Run‐off formation in a humid,temperate headwater catchment using a combined hydrological,hydrochemical and isotopic approach (Jizera Mountains,Czech Republic)

The formation of baseflow and stormflow was examined in the 1.18 km² part of the headwater catchment Uhlí?ská, Jizera Mountains, Czech Republic, over the period 2007–2011, by means of run‐off data and environmental tracers ¹⁸O and SiO₂. The baseflow, computed using the digital filter approach BFLOW, contributes 67% to total streamflow and has a mean residence time of 12.3 months. It is formed by groundwater discharge from the valley deluviofluvial granitic sediments, in combination with soil water in weathered layers on hillslopes during rainfall and snowmelt periods. The prevailing source of the groundwater is the infiltration of snowmelt water. Analysis of 20 run‐off events and their hysteretic patterns demonstrated that the stormflow water has a residence time of about 4 months and is generated by preferential flow on hillslopes combined by soil matrix drainage. Because of slower flow in the soil matrix, the enrichment of pore water in SiO₂ is more pronounced. The stormflow and snowmelt water flowing via preferential pathways of upslope minerals soils pushes the pre‐event groundwater through the pathways in wetlands to the stream, and the wetland can be therefore considered as groundwater supplied. This mechanism has been found to be typical for the groundwater‐supplied headwater catchments of the Jizera Mountains and can be also assumed in other mountainous headwaters of the granitic massif in Central Europe. The main methodological contribution of this study are the residence time calculations stratified by baseflow and event flow, identifying run‐off components of different travel times to streams and linking them with geochemical run‐off sources. This achievement was possible because of a comprehensive dataset on hydrology, stable isotopes and silica hydrochemistry in all relevant run‐off generation components. This concept indicates that a possible long‐term change in snowmelt may affect the run‐off regime of headwater catchments to climate or land‐use changes. Copyright © 2013 John Wiley & Sons, Ltd.     

Sources of uncertainty in the NRCS CN model: Recognition and solutions

The Natural Resources Conservation Service Curve Number model is one of the most recognizable procedures in the field of rainfall–run‐off estimation. It has been widely applied for different purposes in hydrological models. In spite of its widespread use, some uncertainties have not even clarified and must be examined for its proper application. Particularly, choosing the most representative rainfall–run‐off events, and the coefficient λ that relates the parameters of the model (curve number CN and initial abstraction I_a). In this research, an advanced analysis is developed to evaluate the influence of λ for a set of representative watersheds of the Agricultural Research Service of the United Stated Department of Agriculture. They are characterized by different soil properties, land uses, and climatic conditions. Finally, 2 novel methodologies for the selection of the most representative rainfall–run‐off events and for the adaptation of coefficient λ are included, based on the pattern of rainfall distribution.     

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.     

Effects of rock fragment cover on hydrological processes under rainfall simulation in a semi‐arid region of China

Rock fragment cover has long been an important agricultural crop production technique on the Loess Plateau, China. Although this approach plays an important role in controlling hydrological processes and preventing soil erosion, inconsistent results have been recovered in this field. In this study, we investigated the effects of rock fragment cover on infiltration, run‐off, soil erosion, and hydraulic parameters using rainfall simulation in the field in a semi‐arid region of China. Two field plots encompassing 6 rock fragment coverages (0%, 10%, 20%, 25%, 30%, and 40%), as well as 2 rock fragment positions and sizes were exposed to rainfall at a particular intensity (60 mm h^?1). The results of this study showed that increasing the rock fragment coverage with rock fragments resting on the soil surface increased infiltration but decreased run‐off generation and sediment yield. A contrasting result was found, however, when rock fragments were partially embedded into the soil surface; in this case, a positive relationship between rock fragment coverage and run‐off rate as well as a nonmonotonic relationship with respect to soil loss rate was recovered. The size of rock fragments also exerted a positive effect on run‐off generation and sediment yield but had a negative effect on infiltration. At the same time, both mean flow velocity and Froude number decreased with increasing rock fragment coverage regardless of rock fragment position and size, whereas both Manning roughness and Darcy–Weisbach friction factor were positively correlated. Results show that stream power is the most sensitive hydraulic parameter affecting soil loss. Combined with variance analysis, we concluded that the order of significance of rock fragment cover variables was position followed by coverage and then size. We also quantitatively incorporated the effects of rock fragment cover on soil loss via the C and K factors in the Revised Universal Soil Loss Equation. Overall, this study will enable the development of more accurate modelling approaches and lead to a better understanding of hydrological processes under rock fragment cover conditions.     

Downslope soil detachment–transport on steep slopes via rain splash

This study developed a one‐dimensional model of downslope rain splash transport based on field experiments and previous studies. The developed model considers soil detachment processes, ground cover, probability densities, and the effect of overland run‐off in preventing detachment. Field monitoring was conducted to observe precipitation run‐off, ground cover, and sediment production on steep hillslopes. Field‐observed data were used to develop the splash detachment rate equation, probability densities for splash transport, and the maximum splash transport distance. Observed and estimated splash transport showed overall agreement, with some differences for small storm events or events with relatively low intensity, probably caused by variation of overland run‐off depth and connectivity as well as differences in soil surface cohesion at various degrees of wetness. Our model can provide insights on the interactions among rainfall intensity, soil surface condition, soil wetness, and splash transport on forested hillslopes. Copyright © 2011 John Wiley & Sons, Ltd.