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.     

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.     

Comparison of nutrient loss pathways: Run‐off and seepage flow in Vertisols

Recent studies have highlighted the agronomic and environmental importance of phosphorus (P) movement through the soil profile. Thus, faced with challenges such as high‐profile cases of P enrichment of surface water, better understanding of nutrient movement through soil is needed to better manage agricultural fertilizers and manures and their contribution to water quality degradation. In particular, field‐scale research is especially needed in soils with preferential flow transport pathways. Thus, we collected nitrogen (N) and P transport data in run‐off and seepage (lateral subsurface return flow) from 13 field‐ and farm‐scale watersheds on Vertisols in Central Texas for a 14‐year period. For 2004–2017, seepage accounted for ~20% of the total surface flow, and nutrient concentrations were generally similar in run‐off and seepage. As surface run‐off contributed ~80% of the flow, it follows that median annual N and P loads in run‐off were significantly greater than in seepage for every watershed. N loads in both run‐off and seepage flow from cultivated land were an order of magnitude greater than in native prairie and improved pasture, and the highest run‐off and seepage P loads both occurred on cultivated land with organic fertilizer sources. Increasing watershed scale (size) did not to produce consistent patterns in N or P loss in run‐off or seepage. Land use and watershed scale produced significant differences in seepage volume but did not affect run‐off volumes or total surface flow/rainfall. Although less significant in terms of total offsite flux, nutrient movement in vadose zones has important agronomic and environmental implications as considerable N and P are transported through and within the root zone and eventually offsite. And in terms of P, this contradicts the traditionally held scientific viewpoint that P movement through the vadose zone is unimportant agronomically and environmentally.     

Effects of Coal Gangue Content on Water Movement and Solute Transport in a China Loess Plateau Soil

The mining industry has grown strongly in China in recent decades, resulting in large amounts of coal gangues, which cause water and soil pollution, soil erosion, and various other environmental problems. They are often used in reclamation projects in attempts to restore land damaged by mining, hence they are frequently present (in widely varying proportions) in the topsoil in areas around mines. Their presence can strongly affect key soil variables, including its bulk density, structure, water retention, water movement, and solute transport rates. In the study presented here, the effects of gangue contents on infiltration, saturated hydraulic conductivity, and solute transport parameters of a Chinese Loess plateau soil were examined. The results show that infiltration rates and saturated hydraulic conductivity decreased with increasing gangue content. The Peck–Watson equation modeled these relationships well, but Bouwer–Rice equations provided poorer matches with the acquired data. Cumulative infiltration over time was described well by both the Philip equation and Kostiakov equation. Both the simplified convection–dispersion equation and a two‐region model described the solute transport processes well. In addition, the dispersion increased, while both the Peclet number and mobile water fraction decreased, with increases in gangue contents.     

Water and solute movement in a coarse-textured water-repellent field soil

Unstable water flow in water-repellent unsaturated soils can significantly affect the processes of infiltration and soil water redistribution. A field experiment was carried out to study the effect of water-repellency on water and bromide movement in a coarse-textured soil in the southwestern part of The Netherlands. The field data were analyzed using a relatively simple numerical model based on the standard Richards' equation for unsaturated water flow and the Fickian-based convection-dispersion equation for solute transport. Water-repellency was accounted for by multiplying the water content and the unsaturated hydraulic conductivity of the soil with F, a factor equal to the volumetric fraction of soil occupied by preferential flow paths resulting from the unstable flow process. The good comparison of simulated and measured bromide concentrations suggests that the model provides a viable method for simulating unstable water flow in water-repellent soils.     

Simulations of water movement and solute transport through different soil texture configurations under negative‐pressure irrigation

This study examined the effects of different soil texture configurations on water movement and solute transport to provide a reliable scientific basis for the application of negative‐pressure irrigation (NPI) technology. HYDRUS‐2D was used to analyse water movement and solute transport under NPI. The main results are as follows: (a) HYDRUS‐2D can be used to simulate water movement and solute transport under NPI, as there was good agreement between the simulated and measured values for water contents, NaCl concentrations, cumulative water infiltration, and wetting distances in the horizontal and vertical directions; the Nash–Sutcliffe efficiency coefficients were in the range of 0.94–0.97. (b) Layered soils have obvious effects on water movement under NPI. With the emitter position in the loam layer, when a coarse texture of loamy sand was present below the loam layer (namely, L‐LS), irrigation water accumulated in the topsoil, and this led to an increase in evaporation compared with the homogeneous loam profile. However, fine texture silty loam or silty clay loam layers beneath the loam layer (namely, L‐SiL or L‐SiCL, respectively) was more conducive to water infiltration into the lower layer, and this increased the amount of water infiltration and simultaneously reduced the surface evaporation effectively. (c) Layered soils have obvious effects on solute transport under NPI, and salt accumulation will readily occur in the clay‐rich soil layer at the interface. The maximum soil salt accumulation of L‐LS occurred above the soil interface between the two soil layers with a value of 21.80 g/kg; however, for L‐SiCL and L‐SiL, the maximum salt accumulation occurred below the soil interface between the two soil layers, with values of 23.80 g/kg and 20.08 g/kg, respectively. (d) Interlayered soils showed remarkable changes in the water infiltration characteristics and salt‐leaching intensities under NPI, and the properties for the soil profile with a silty loam interlayer were better than those for the soil profile with a silty clay loam interlayer. The soil profile with a loamy sand interlayer had the lowest amount of water infiltration, which resulted in reductions of the salt‐leaching intensities. Thus, NPI is clearly not suitable for loamy sand soil. Overall, the results demonstrated that soil texture configurations affected water movement and solute transport under NPI. Therefore, careful consideration should be given to the use of NPI to achieve target soil water and solution conditions and reduce water loss.     

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.     

Catchment‐scale conceptual modelling of water and solute transport in the dual flow system of the karst critical zone

Hydrological and hydrochemical processes in the critical zone of karst environments are controlled by the fracture‐conduit network. Modelling hydrological and hydrochemical dynamics in such heterogeneous hydrogeological settings remains a research challenge. In this study, water and solute transport in the dual flow system of the karst critical zone were investigated in a 73.5‐km² catchment in southwest China. We developed a dual reservoir conceptual run‐off model combined with an autoregressive and moving average model with algorithms to assess dissolution rates in the “fast flow” and “slow flow” systems. This model was applied to 3 catchments with typical karst critical zone architectures, to show how flow exchange between fracture and conduit networks changes in relation to catchment storage dynamics. The flux of bidirectional water and solute exchange between the fissure and conduit system increases from the headwaters to the outfall due to the large area of the developed conduits and low hydraulic gradient in the lower catchment. Rainfall amounts have a significant influence on partitioning the relative proportions of flow and solutes derived from different sources reaching the underground outlet. The effect of rainfall on catchment function is modulated by the structure of the karst critical zone (e.g., epikarst and sinkholes). Thin epikarst and well‐developed sinkholes in the headwaters divert more surface water (younger water) into the underground channel network, leading to a higher fraction of rainfall recharge into the fast flow system and total outflow. Also, the contribution of carbonate weathering to mass export is also higher in the headwaters due to the infiltration of younger water with low solute concentrations through sinkholes.     

Experimental study and mathematical modelling of soluble chemical transfer from unsaturated/saturated soil to surface runoff

A one‐dimensional, two‐layer solute transport model is developed to simulate chemical transport process in an initially unsaturated soil with ponding water on the soil surface before surface runoff starts. The developed mathematical model is tested against a laboratory experiment. The infiltration and diffusion processes are mathematically lumped together and described by incomplete mixing parameters. Based on mass conservation and water balance equations, the model is developed to describe solute transport in a two‐zone layer, a ponding runoff zone and a soil mixing zone. The two‐zone layer is treated as one system to avoid describing the complicated chemical transport processes near the soil surface in the mixing zone. The proposed model was analytically solved, and the solutions agreed well with the experimental data. The developed experimental method and mathematical model were used to study the effect of the soil initial moisture saturation on chemical concentration in surface runoff. The study results indicated that, when the soil was initially saturated, chemical concentration in surface runoff was significantly (two orders of magnitude) higher than that with initially unsaturated soil, while the initial chemical concentrations at the two cases were of the same magnitude. The soil mixing depth for the initially unsaturated soil was much larger than that for the initially saturated soil, and the incomplete runoff mixing parameter was larger for the initially unsaturated soil. The higher the infiltration rate of the soil, the greater the infiltration‐related incomplete mixing parameter. According to the quantitative analysis, the soil mixing depth was found to be sensitive for both initially unsaturated and saturated soils, and the incomplete runoff mixing parameter was sensitive for initially saturated soil but not for the initially unsaturated soil; the incomplete infiltration mixing parameter behaved just the opposite. Some suggestions are made for reducing chemical loss from runoff. Copyright © 2010 John Wiley & Sons, Ltd.     

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.     

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.     

Solute transport by surface runoff from low‐angle slopes: theory and application

The removal of chemicals in solution by overland flow from agricultural land has the potential to be a significant source of chemical loss where chemicals are applied to the soil surface, as in zero tillage and surface‐mulched farming systems. Currently, we lack detailed understanding of the transfer mechanism between the soil solution and overland flow, particularly under field conditions. A model of solute transfer from soil solution to overland flow was developed. The model is based on the hypothesis that a solute is initially distributed uniformly throughout the soil pore space in a thin layer at the soil surface. A fundamental assumption of the model is that at the time runoff commences, any solute at the soil surface that could be transported into the soil with the infiltrating water will already have been convected away from the area of potential exchange. Solute remaining at the soil surface is therefore not subject to further infiltration and may be approximated as a layer of tracer on a plane impermeable surface. The model fitted experimental data very well in all but one trial. The model in its present form focuses on the exchange of solute between the soil solution and surface water after the commencement of runoff. Future model development requires the relationship between the mass transfer parameters of the model and the time to runoff to be defined. This would enable the model to be used for extrapolation beyond the specific experimental results of this study. The close agreement between experimental results and model simulations shows that the simple transfer equation proposed in this study has promise for estimating solute loss to surface runoff. Copyright © 2000 John Wiley & Sons, Ltd.     

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.     

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.     

Hydrologic flowpaths during snowmelt in forested headwater catchments under differing winter climatic and soil frost regimes

Changes in hydrologic flowpaths have important impacts on the timing, magnitude and hydrochemistry of run‐off during snowmelt in forested catchments, but how flowpaths are affected by variation in winter climate and the irregular presence of soil frost remains poorly understood. The depth and extent of soil frost may be expected to increase as snowpack decreases or develops later because of climate change. In this study, we used end‐member mixing analysis to determine daily contributions of snow, forest floor soil water and groundwater to stream run‐off during snowmelt under different soil frost regimes resulting from interannual and elevational variation at the Hubbard Brook Experimental Forest in New Hampshire, USA. We observed greater routing of run‐off through forest floor flowpaths during early snowmelt in 2011, when the snowpack was deep and soil frost was minimal, compared with the early snowmelt in 2012 under conditions of deep and extensive soil frost. The results indicate that widespread soil frost that penetrated the depth of the forest floor decreased the flow signal through the shallowest subsurface flowpaths, but did not reduce overall infiltration of melt waters, as the contribution from the snow‐precipitation end‐member was similar under both conditions. These results are consistent with development of granular soil frost which permits vertical infiltration of melt waters, but either reduces lateral flow in the forest floor or prevents the solute exchange that would produce the typical chemical signature of shallow subsurface flowpaths in streamwater. Copyright © 2016 John Wiley & Sons, Ltd.     

Simulation of field‐scale water flow and bromide transport in a cracked clay soil

Single domain models may seriously underestimate leaching of nutrients and pesticides to groundwater in clay soils with shrinkage cracks. Various two‐domain models have been developed, either empirical or physically based, which take into account the effects of cracks on water flow and solute transport. This paper presents a model concept that uses the clay shrinkage characteristics to derive crack volume and crack depth under transient field conditions. The concept has been developed to simulate field average behaviour of a field with cracks, rather than flow and transport at a small plot. Water flow and solute transport are described with basic physics, which allow process and scenario analysis. The model concept is part of the more general agrohydrological model SWAP, and is applied to a field experiment on a cracked clay soil, at which water flow and bromide transport were measured during 572 days. A single domain model was not able to mimic the field‐average water flow and solute transport. Incorporation of the crack concept considerably improved the simulation of water content and bromide leaching to the groundwater. Still deviations existed between the measured and simulated bromide concentration profiles. The model did not reproduce the observed bromide retardation in the top layer and the high bromide dispersion resulting from water infiltration at various soil depths. A sensitivity analysis showed that the amounts of bromide leached were especially sensitive to the saturated hydraulic conductivity of the top layer, the solute transfer from the soil matrix to crack water flow and the mean residence time of rapid drainage. The shrinkage characteristic and the soil hydraulic properties of the clay matrix showed a low sensitivity. Copyright © 2000 John Wiley & Sons, Ltd.     

LISEM: A SINGLE-EVENT PHYSICALLY BASED HYDROLOGICAL AND SOIL EROSION MODEL FOR DRAINAGE BASINS. I: THEORY,INPUT AND OUTPUT

A new physically based hydrological and soil erosion model has been developed, which can be used for planning and conservation purposes: the LImburg Soil Erosion Model (LISEM). The LISEM model is one of the first examples of a physically based model that is completely incorporated in a raster Geographical Information System. This incorporation facilitates easy application in larger catchments, improves the user friendliness by avoiding conversion routines and allows remotely sensed data to be used. Processes incorporated in the model are rainfall, interception, surface storage in micro-depressions, infiltration and vertical movement of water in the soil, overland flow, channel flow, detachment by rainfall and throughfall, detachment by overland flow and transport capacity of the flow. Special attention has been given to the influence of tractor wheelings, small roads and surface sealing. Vertical movement of water in the soil is simulated using the Richard's equation. Optionally, the user can choose the Holtan or the Green–Ampt infiltration model. For the distribution flow routing, a four-point finite-difference solution of the kinematic wave is used together with Manning's equation.     

A quasi three-dimensional model for flow and transport in unsaturated and saturated zones: 1. Implementation of the quasi two-dimensional case

A quasi three-dimensional (QUASI 3-D) model is presented for simulating the subsurface water flow and solute transport in the unsaturated and in the saturated zones of soil. The model is based on the assumptions of vertical flow in the unsaturated zone and essentially horizontal groundwater flow. The 1-D Richards equation for the unsaturated zone is coupled at the phreatic surface with the 2-D flow equation for the saturated zone. The latter was obtained by averaging 3-D flow equation in the saturated zone over the aquifer thickness. Unlike the Boussinesq equation for a leaky-phreatic aquifer, the developed model does not contain a storage term with specific yield and a source term for natural replenishment. Instead it includes a water flux term at the phreatic surface through which the Richards equation is linked with the groundwater flow equation. The vertical water flux in the saturated zone is evaluated on the basis of the fluid mass balance equation while the horizontal fluxes, in that equation, are prescribed by Darcy law. A 3-D transport equation is used to simulate the solute migration. A numerical algorithm to solve the problem for the general quasi 3-D case was developed. The developed methodology was exemplified for the quasi 2-D cross-sectional case (QUASI2D). Simulations for three synthetic problems demonstrate good agreement between the results obtained by QUASI2D and two fully 2-D flow and transport codes (SUTRA and 2DSOIL). Yet, simulations with the QUASI2D code were several times faster than those by the SUTRA and the 2DSOIL codes.     

Temporal variability of solute transport under vadose zone conditions

The heterogeneity of the solute flux field in the horizontal plane at the field scale has been documented in several field studies. On the other hand, little information is available on the persistence of certain solute transport scenarios over consecutive infiltration cycles. This study was initiated to analyse the recurrence of solute leaching behaviour as estimated in two soil column tests emphasizing the preferential flow phenomenon. Twenty-four small-sized soil samples were subjected to two consecutive unsaturated steady-state flow leaching experiments with bromide as tracer. Observed breakthrough curves (BTCs) were analysed by the method of moments and by the advection–dispersion equation (ADE) to classify solute behaviour. Frequency distributions of the parameters indicating the solute velocity were heavily skewed or bimodal, reflecting the broad variability of the leaching scenarios, including some with pronounced preferential solute breakthrough. Exclusion of the preferential flow columns from our calculations revealed an average amount of 37% of immobile water. The large-scale BTCs derived from assembling the individual concentration courses of each run showed similar features, such as an early bromide breakthrough. However, two distinct apices, viz. one preferential and one matrix, were observed only in the first run, whereas the concentration decrease between the peaks was missing from the second run. A change in soil structure with continuous leaching was presumed to modify the interplay of the various flow domains, thereby altering the spreading of the BTCs. Correlation analysis between parameters of both tests suggests that preferential transport conditions are likely to occur at the same locations in the field over several infiltration cycles, whereas the ‘classical’ or expected matrix flow is time variant and therefore seems to be hardly predictable. © 1998 John Wiley & Sons, Ltd.     

Relation of sediment transport capacity to stone cover and size in rain‐impacted interrill flow

This study examines how the sediment transport capacity of interrill overland flow varies with stone cover and stone size at two flow intensities. Six series of flume experiments were conducted on two slopes (2° and 10°) with stones of three sizes (28·0, 45·5 and 91·3 mm) serving as roughness elements. Bed sediment size, water discharge and simulated rainfall intensity were the same in all experiments. It was found (1) that transport capacity is positively related to stone size, with the relation becoming stronger as stone cover increases and flow intensity decreases; and (2) that transport capacity is negatively related to stone cover at the high flow intensity and curvilinearly related to stone cover at the low flow intensity. The curvilinear relations are concave‐upward with the lowest transport capacities occurring at stone covers between 0·40 and 0·60. The highest transport capacities are found at stone covers of 0 and 1, with the transport capacity being greater at the former stone cover than at the latter. Copyright © 2000 John Wiley & Sons, Ltd.