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

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

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.     

Effects of rock fragments incorporated in the soil matrix on concentrated flow hydraulics and erosion

Rock fragments can act as a controlling factor for erosional rates and patterns in the landscape. Thus, the objective of this study is to better understand the role that rock fragments incorporated into the soil matrix play in concentrated flow hydraulics and erosion. Laboratory flume experiments were conducted with soil material that was mixed with rock fragments. Rock fragment content ranged from 0 to 40 per cent by volume. Other treatments were slope (7 and 14%) and flow discharge (5·7 and 11·4 l min^?1). An increase in rock fragment content resulted in lower sediment yield, and broader width of flow. Rock fragment cover at the soil surface, i.e. surface armour, increased with time in experiments with rock fragments. Flow energy was largely dissipated by rock fragment cover. For more turbulent flow conditions, when roughness elements were submerged in the flow, hydraulic roughness was similar for different rock fragment contents. In experiments with few or no rock fragments a narrow rill incised. Flow energy was dissipated by headcuts. Total sediment yield was much larger than for experiments with rock fragments in the soil. Adding just a small number of rock fragments in the soil matrix resulted in a significant reduction of sediment yield. Copyright © 2007 John Wiley & Sons, Ltd.     

Runoff generation and routing on artificial slopes in a Mediterranean–continental environment: the Teruel coalfield,Spain

The aim of this study was to identify the mechanisms of runoff generation and routing and their controlling factors at the hillslope scale, on artificial slopes derived from surface coal mining reclamation in a Mediterranean–continental area. Rainfall and runoff at interrill and microcatchment scales were recorded for a year on two slopes with different substrata: topsoil cover and overburden cover. Runoff coefficient and runoff routing from interrill areas to microcatchment outlets were higher in the overburden substratum than in topsoil, and greater in the most developed rill network. Rainfall volume is the major parameter responsible for runoff response on overburden, suggesting that this substratum is very impermeable—at least during the main rainfall periods of the year (late spring and autumn) when the soil surface is sealed. In such conditions, most rainfall input is converted into runoff, regardless of its intensity. Results from artificial rainfall experiments, conducted 3 and 7 years after seeding, confirm the low infiltration capacity of overburden when sealed. The hydrological response shows great seasonal variability on the overburden slope in accordance with soil surface changes over the year. Rainfall volume and intensities (I₃₀, I₆₀) explain runoff at the interrill scale on the topsoil slope, where rainfall experiments demonstrated a typical Hortonian infiltration curve. However, no correlation was found at the microcatchment level, probably because of the loss of functionality of the only rill as ecological succession proceeded. The runoff generation mechanism on the topsoil slope is more homogeneous throughout the year. Runoff connectivity, defined as the ratio between runoff rates recorded at the rill network scale and those recorded at the interrill area scale in every rainfall event, was also greater on the rilled overburden slope, and in the most developed rill network. The dense rill networks of the overburden slope guarantee very effective runoff drainage, regardless of rainfall magnitude. Rills drain overland flow from interrill‐sealed areas, reducing the opportunity of reinfiltration in areas not affected by siltation. Runoff generation and routing on topsoil slopes are controlled by grass cover and soil moisture content, whereas on overburden slopes rill network density and soil moisture content are the main controlling factors. Copyright © 2002 John Wiley & Sons, Ltd.     

Sediment concentration in interrill flow: interactions between soil surface conditions,vegetation and rainfall

A database composed of 673 natural rainfall events with sediment concentration measurements at the field or plot scale was analysed. Measurements were conducted on similar soil type (loess soils prone to sealing phenomenon) to apprehend the variability and complexity involved in interrill erosion processes attributable to soil surface conditions. The effects of the dominant controlling factors are not described by means of equations; rather, we established a classification of potential sediment concentration domain according to combination of the dominant parameters. Thereby, significant differences and evolution trends of mean sediment concentration between the different parameter categories are identified. Further, when parameter influences interact, it allows us to discern the relative effects of factors according to their respective degree of expression. It was shown that crop cover had a major influence on mean sediment concentration, particularly when soil surface roughness is low and when maximum 6‐min intensity of rainfall events exceeds 10 mm h^?1: mean sediment concentration decreases from 8·93 g l^?1 for 0–20 per cent of coverage to 0·97 g l^?1 for 21–60 per cent of coverage. The established classification also indicates that the increase of the maximum 6‐min intensity of the rainfall factor leads to a linear increase of mean sediment concentration for crop cover over 21 per cent (e.g. from 2·96 g l^?1 to 14·44 g l^?1 for the 1–5 cm roughness class) and to an exponential increase for low crop cover (e.g. from 3·92 g l^?1 to 58·76 g l^?1 for the 1–5 cm roughness class). The implication of this work may bring perspective for erosion prediction modelling and give references for the development of interrill erosion equation. Copyright © 2002 John Wiley & Sons, Ltd.     

Effects of rock fragment size and cover on overland flow hydraulics,local turbulence and sediment yield on an erodible soil surface

The interactions between overland flow hydraulics and sediment yield were studied in flume experiments on erodible soil surfaces covered by rock fragments. The high erodibility of a non-cohesive fine sediment (D₅₀ + 0·09mm) permitted the effects of local turbulence and scour on sediment yield to be examined. Overland flow hydraulics and sediment yield were compared for experiments with pebble (D₅₀ + 1·5cm) and cobble (D₅₀ + 8·6cm) rock fragment covers. Cover percentages range from 0 to 99 per cent. Rock fragment size strongly affects the relations between flow hydraulics and rock fragment cover. For pebbles spatially-averaged hydraulic parameters (flow velocity, flow depth, effective flow width, unit discharge, total shear stress, Darcy-Weisbach friction factor, percentage grain friction and grain shear stress) vary most rapidly within cover percentages at low covers (power functions). In contrast, for cobbles these parameters vary most rapidly within cover percentages at high covers (exponential functions). As the type of the function that describes the relation between flow hydraulics and cover percentage can be deduced from the ratio of rock fragment height to flow depth, the continuity equation can be employed to determine the actual coefficients of the functions, provided the regression of one hydraulic parameter (e.g. flow velocity) with cover percentage is known and a good estimate exists for two values of another hydraulic variable for a low and a high cover percentage. The variation of sediment yield with cover percentage is also strongly dependent on rock fragment size, but neither the convex-upward relation for pebbles, nor the positive relation for cobbles can be solely attributed to the spatially averaged hydraulics of sheet-flow. Rock fragments induce local turbulence that leads to scour hole development on the stoss side of the rock fragments while deposition commonly occurs in the wake. This local scour and deposition substantially affects sediment yield. However, scour dimensions cannot be predicted by spatially averaged flow hydraulics. An adjustment of existing scour formulas that predict scour around bridge piers is suggested. Sediment yield from non-cohesive soils might then be estimated by a combination of sediment transport and scour formulas.     

Fire severity and surface rock fragments cause patchy distribution of soil water repellency and infiltration rates after burning

Although fire‐induced soil water repellency (SWR) and its effects on soil hydrology and geomorphology have been studied in detail, very few studies have considered the effect of rock fragments resting on the soil surface or partly embedded in soil. In this research, we have studied the effect of rock fragments on the strength and spatial distribution of fire‐induced SWR at different fire severities. A fire‐affected area was selected for this experiment and classified into different zones according to fire severity (unburned, low, moderate and high) and rock fragment cover (low, <20% and high, >60%). During 7 days after fire, SWR and infiltration rates were assessed in the soil surface covered by individual rock fragments and in the midpoint between two adjacent rock fragments (with maximum spacing of 20 cm). SWR increased with fire severity. Rock fragments resting on the soil surface increased the heterogeneity of the spatial distribution of fire‐induced SWR. SWR increased significantly with rock fragment cover in bare areas under moderate and high fire severity, but quantitatively important changes were only observed under high fire severity. In areas with a low rock fragment cover, water repellency from soil surfaces covered by rock fragments increased relative to bare soil surfaces, with increasing SWR. In areas with a high rock fragment cover, SWR increased significantly from non‐covered to covered soil surfaces only after low‐severity burning. Rock fragment cover did not affect infiltration rates, although it decreased significantly in soil surfaces after high‐severity burning in areas under low and high rock fragment cover. Copyright © 2013 John Wiley & Sons, Ltd.     

The hydrological response of soil surfaces to rainfall as affected by cover and position of rock fragments in the top layer

Rainfall experiments have been conducted in the laboratory in order to assess the hydrological response of top soils very susceptible to surface sealing and containing rock fragments in different positions with respect to the soil surface. For a given cover level, rock fragment position in the top soil has an ambivalent effect on water intake and runoff generation. Compared to a bare soil surface rock fragments increase water intake rates as well as time of runoff concentration and decrease runoff volume if they rest on the soil surface. For the same cover level, rock fragments reduce infiltration rate and enhance runoff generation if they are well embedded in the top layer. The effects of rock fragment position on infiltration rate and runoff generation are proportional to cover percentage. Micromorphological analysis and measurements of the saturated hydraulic conductivity of bare top soils and of the top layer underneath rock fragments resting on the soil surface reveal significant differences supporting the mechanism proposed by Poesen (1986): i.e. runoff generated as rock flow or as Horton overland flow can (partly) infiltrate into the unsealed soil surface under the rock fragments, provided that they are not completely embedded in the top layer. Hence, rock fragment position, beside other rock fragment properties, should be taken into account when assessing the hydrological response of soils susceptible to surface sealing and containing rock fragments in their surface layers. A simple model, based on the proportions of bare soil surface, soil surface occupied by embedded rock fragments, and soil surface covered with rock fragments resting on the soil surface, describes the runoff coefficient data relatively well.     

Relevant effects of vegetal cover and litter on the soil hydrological response of two contrasting Mediterranean hillslopes at the end of the dry season (south of Spain)

In Mediterranean regions, hillslopes are generally considered to be a mosaic of sink and source areas that control runoff generation and water erosion processes. These hillslopes used to be characterized by a complex hydrological and erosive response combining Hortonian and saturation excess overland flows. The hydrological response of soils is highly dependent on the soil surface components (e.g. vegetation patches, bare soil, rock fragment cover, crusts), which each one of them is dominated by a certain hydrological process. One of these soil surface components, not widely considered in studies of soil hydrology under Mediterranean conditions, is the accumulation of litter beneath shrubs enhancing water repellency in soils. This study investigates the influence of soil surface components, especially the litter accumulated beneath Cistus spp., in the hydrological and erosive responses of soils on two Mediterranean hillslopes having different exposures. The study was performed by means of rainfall simulation experiments and the Water Drop Penetration Time for measuring water repellency of soils, both techniques being carried out at the end of summer (September 2010) with very dry soils. The results indicate that (i) soil surface components from the north facing hillslope are characterized by a more uniform hydrological and erosive response than those from the south‐facing ones; (ii) the water repellency is more influential on the hydrological response of the north‐facing hillslope due to a greater accumulation of organic rest on the soils as the vegetation cover is also higher; (iii) the south‐facing hillslope seemed to follow the fertility island theory with very degraded bare soil areas, which are the most generated areas of runoff and mobilized sediments; (iv) the experimental area can be considered as a threshold area between the semiarid and subhumid Mediterranean environments, with the south‐facing hillslope being comparable with the former and the north facing one with the latter. Copyright © 2011 John Wiley & Sons, Ltd.     

Sensitivity analysis of EUROSEM using Monte Carlo simulation II: the effect of rills and rock fragments

A sensitivity analysis of the surface and catchment characteristics in the European soil erosion model (EUROSEM) was carried out with special emphasis on rills and rock fragment cover. The analysis focused on the use of Monte Carlo simulation but was supplemented by a simple sensitivity analysis where input variables were increased and decreased by 10%. The study showed that rock fragments have a significant effect upon the static output parameters of total runoff, peak flow rate, total soil loss and peak sediment discharge, but with a high coefficient of variation. The same applied to the average hydrographs and sedigraphs although the peak of the graphs was associated with a low coefficient of variation. On average, however, the model was able to simulate the effect of rock fragment cover quite well. The sensitivity analysis through the Monte Carlo simulation showed that the model is particularly sensitive to changes in parameters describing rills and the length of the plane when no rock fragments are simulated but that the model also is sensitive to changes in the fraction of non‐erodible material and interrill slope when rock fragments were embedded in the topsoil. For rock fragments resting on the surface, changes in parameter values did not affect model output significantly. The simple sensitivity analysis supported the findings from the Monte Carlo simulation and illustrates the importance when choosing input parameters to describe both rills and rock fragment cover when modelling with EUROSEM. Copyright © 2000 John Wiley & Sons, Ltd.     

Spatial distribution of surface rock fragment on hill-slopes in a small catchment in wind-water erosion crisscross region of the Loess Plateau

The cover and size distributions of surface rock fragment in hillslopes were investigated by using digital photographing and treating technique in a small catchment in wind-water erosion crisscross region of the Loess Plateau. The results indicated that the maximal cover of rock fragment was pre-sented at mid-position in steep hillslope. Rock fragment presented a general decreasing-trend along the hillslope in gentle hillslope. Rock fragment cover was positively related to gradient, rock fragment size decreased generally along the hillslope, and the size reduced with the gradient. The mean size of rock fragment was at a range of 6―20 mm in the steep hillslope, rock fragment size > 50 mm was rarely presented. The covers of rock fragment at different positions were markedly related to the quantities of rock fragment < 40 mm. The area of rock fragment of 2―50 mm accounted for 60% or more of the total area, dominating the distribution of rock fragment in the hillslopes.     

Fire effects on rangeland hydrology and erosion in a steep sagebrush‐dominated landscape

Post‐fire runoff and erosion from wildlands has been well researched, but few studies have researched the degree of control exerted by fire on rangeland hydrology and erosion processes. Furthermore, the spatial continuity and temporal persistence of wildfire impacts on rangeland hydrology and erosion are not well understood. Small‐plot rainfall and concentrated flow simulations were applied to unburned and severely burned hillslopes to determine the spatial continuity and persistence of fire‐induced impacts on runoff and erosion by interrill and rill processes on steep sagebrush‐dominated sites. Runoff and erosion were measured immediately following and each of 3 years post‐wildfire. Spatial and temporal variability in post‐fire hydrologic and erosional responses were compared with runoff and erosion measured under unburned conditions. Results from interrill simulations indicate fire‐induced impacts were predominantly on coppice microsites and that fire influenced interrill sediment yield more than runoff. Interrill runoff was nearly unchanged by burning, but 3‐year cumulative interrill sediment yield on burned hillslopes (50 g m^?2) was twice that of unburned hillslopes (25 g m^?2). The greatest impact of fire was on the dynamics of runoff once overland flow began. Reduced ground cover on burned hillslopes allowed overland flow to concentrate into rills. The 3‐year cumulative runoff from concentrated flow simulations on burned hillslopes (298 l) was nearly 20 times that measured on unburned hillslopes (16 l). The 3‐year cumulative sediment yield from concentrated flow on burned and unburned hillslopes was 20 400 g m^?2 and 6 g m^?2 respectively. Fire effects on runoff generation and sediment were greatly reduced, but remained, 3 years post‐fire. The results indicate that the impacts of fire on runoff and erosion from severely burned steep sagebrush landscapes vary significantly by microsite and process, exhibiting seasonal fluctuation in degree, and that fire‐induced increases in runoff and erosion may require more than 3 years to return to background levels. Published in 2008 by John Wiley & Sons, Ltd.     

Mining soil databases for landscape‐scale patterns in the abundance and size distribution of hillslope rock fragments

Landscape‐scale variation in rock fragments on soil‐mantled hillslopes is poorly understood, despite the potential importance of rock fragments in soil weathering and coarse sediment supply to river networks. We explored the utility of soil survey databases for data mining, with the goals of identifying landscape‐scale patterns in the abundance and size distribution of rock fragments (diameter D > 2 mm) and potential controls on grain size production. We focus on data from three regions: the Hawaiian Islands, and the Sierra Nevada and Cascade Mountains, where elevation transects span a range of environmental conditions. We selected pedons from pits dug on hillslopes with active soil production and transport. For the 27 pedons selected, we constructed depth‐averaged grain size distributions and calculated the mass fraction of rock fragments (F_RF) and the median rock fragment grain size (D_50RF). We also categorized as bimodal, size distributions with a clear ‘breakpoint’ between fine and coarse modes. Several strong patterns emerge from the data. We find rock fragments in 85% of the pedons, primarily in distinct coarse modes within bimodal size distributions. Values of F_RF and D_50RF are strongly correlated, although the best‐fit power law scaling between F_RF and D_50RF differs between the warmer Hawaiian, and colder Sierra Nevada and Cascade Mountain sites. We also find a regional contrast in the variation in F_RF with elevation; F_RF declines with elevation in Hawaii, but increases in the mainland sites. Although this contrast could be an artifact of variable lithology, precipitation may influence many patterns in the data. Lower mean‐annual precipitation correlates with higher F_RF, dominantly bimodal distributions and surface enrichment in the vertical distribution of rock fragments. These observations may be useful in refining models of coarse sediment supply to rivers, and suggest opportunities for future work to test mechanistic hypotheses for rock fragment production on soil‐mantled hillslopes. Copyright © 2011 John Wiley & Sons, Ltd.     

Impact of rainfall pattern on interrill erosion process

The impact of rainfall pattern on the interrill erosion process is not fully understood despite its importance. Systematic rainfall simulation experiments involving various rainfall intensities, stages, intensity sequences, and surface cover conditions were conducted in this study to investigate their effects on the interrill erosion process. Five rainfall patterns designed with the same total kinetic energy/precipitation (increasing, decreasing, rising–falling, falling–rising and constant patterns) were randomly delivered to a pre‐wet clay loam soil surface at a 10° slope gradient. Significant differences in soil losses were observed among the different rainfall patterns and stages, but there was no obvious difference in runoff. Kinetic energy flux (KE_r) was a governing factor for interrill erosion, and constant rainfall pattern (CST) produced nine times greater soil loss than runs with no KE_r. Varied‐intensity patterns had a profound effect on raindrop‐induced sediment transport processes; path analysis results indicated that said effect was complex, interactive and intensity‐dependent. Low hydraulic parameter thresholds further indicated that KE_r was the dominant factor in detaching soil particles, while overland flow mainly contributed to transporting the pre‐detached particles. This study not only sheds light on the mechanism of interrill sediment transport capacity and detachability, but also may provide a useful database for developing event‐based interrill erosion prediction models. Copyright © 2017 John Wiley & Sons, Ltd.     

Changes of soil surface roughness under water erosion process

The objective of this study was to determine the changing characteristics of soil surface roughness under different rainfall intensities and examine the interaction between soil surface roughness and different water erosion processes. Four artificial management practices (raking cropland, artificial hoeing, artificial digging, and contour tillage) were used according to the local agriculture customs of the Loess Plateau of China to simulate different types of soil surface roughness, using an additional smooth slope for comparison purposes. A total of 20 rainfall simulation experiments were conducted in five 1 m by 2 m boxes under two rainfall intensities (0.68 and 1.50 mm min^?1) on a 15° slope. During splash erosion, soil surface roughness decreased in all treatments except raking cropland and smooth baseline under rainfall intensity of 0.68 mm min^?1, while increasing for all treatments except smooth baseline under rainfall intensity of 1.50 mm min^?1. During sheet erosion, soil surface roughness decreased for all treatments except hoeing cropland under rainfall intensity of 0.68 mm min^?1. However, soil surface roughness increased for the artificial hoeing and raking cropland under rainfall intensity of 1.50 mm min^?1. Soil surface roughness has a control effect on sheet erosion for different treatments under two rainfall intensities. For rill erosion, soil surface roughness increased for raking cropland and artificial hoeing treatments, and soil surface roughness decreased for artificial digging and the contour tillage treatments under two rainfall intensities. Under rainfall intensity of 0.68 mm min^?1, the critical soil surface roughness was 0.706 cm for the resistance control of runoff and sediment yield. Under rainfall intensity of 1.50 mm min^?1, the critical soil surface roughness was 1.633 cm for the resistance control of runoff, while the critical soil surface roughness was 0.706 cm for the resistance control of sediment yield. These findings have important implications for clarifying the erosive nature of soil surface roughness and harnessing sloped farmland. Copyright © 2013 John Wiley & Sons, Ltd.     

Evaluation of an interrill soil erosion model using laboratory catchment data

Physically based soil erosion simulation models require input parameters of soil detachment and sediment transport owing to the action and interactions of both raindrops and overland flow. A simple interrill soil water transport model is applied to a laboratory catchment to investigate the application of raindrop detachment and transport in interrill areas explicitly. A controlled laboratory rainfall simulation study with slope length simulation by flow addition was used to assess the raindrop detachment and transport of detached soil by overland flow in interrill areas. Artificial rainfall of moderate to high intensity was used to simulate intense rain storms. However, experiments were restricted to conditions where rilling and channelling did not occur and where overland flow covered most of the surface. A simple equation with a rainfall intensity term for raindrop detachment, and a simple sediment transport equation with unit discharge and a slope term were found to be applicable to the situation where clear water is added at the upper end of a small plot to simulate increased slope length. The proposed generic relationships can be used to predict raindrop detachment and the sediment transport capacity of interrill flow and can therefore contribute to the development of physically‐based erosion models. Copyright © 1999 John Wiley & Sons, Ltd.     

Spatial patterns of soil water balance on intensively cultivated hillslopes in a semi‐arid environment: the impact of rock fragments and soil thickness

During past decades, a diverse system of subsistence agriculture in south‐east Spain (annual rainfall of less than 300 mm) has been overturned in favour of large‐scale plantations of almond trees without consideration for topography and related spatial patterns in soil hydrological properties. The objective of this paper is to investigate the spatial pattern in soil physical properties induced by this cultivation system, and to highlight its impact on the water balance. Soil properties were recorded along hillslopes with shallow soils developed on slates and greywackes in the upper part of the Guadelentin drainage basin (Murcia region). Frequent tillage of these almond plantations covering entire hillslopes has resulted in denudation by tillage erosion on the topographic convexities, as well as transport of rock fragments and fine earth along the slopes. These processes have created a systematic spatial pattern of soil thickness and rock fragment content: shallow and stony soils on the topographic convexities and deep soils with a rock fragment mulch in the concavities at the foot of the slopes. At the same time, a negative relationship between rock fragment content and fine earth bulk density was observed. The impact of this spatial pattern in soil properties on the water balance was evaluated using the PATTERN one‐dimensional hydrological and plant growth model. The model simulates the water balance of soil profiles covering the observed variation in soil thickness, stoniness and bulk density. The model results indicate that the highest rates of infiltration, evaporation and drainage, as well as the lowest rates of overland flow are restricted to shallow soils on the hilltops. In contrast, the deeper soils in the valley bottoms produce a more stable moisture regime than shallower soils, which tend to saturate and dry out quickly. These model results are in agreement with the spatial patterns of almond productivity: an asymptotic increase with soil thickness. Copyright © 2000 John Wiley & Sons, Ltd.     

Thresholds for gully initiation and sedimentation in Mediterranean Europe

In Mediterranean areas the dynamics of gully development act as an important indicator of desertification. However, little is known about the influence of climate and land‐use changes, and almost no field data exist to assess the sensitivity of a landscape to gully erosion. Two important components of gully erosion studies are the prediction of where gullies begin and where they end. To address some of these issues, topographical thresholds for gully initiation and sedimentation in six different Mediterranean study areas were established. Field measurements of local soil surface slope (S) and drainage‐basin area (A) at the point of initiation of ephemeral gullies in intensively cultivated fields (five datasets) and permanent gullies in rangelands (three datasets) were carried out. A negative power relationship of the form S = aA^−b was fitted through all datasets, and defined as the mean topographical threshold for gullying in the respective area. Topographically controlled slopes of sedimentation at the gully bottom were also measured. Compared to theoretical relationships for channel initiation by overland flow, relatively low values for b are obtained, suggesting a dominance of overland flow and an influence of subsurface flow. The influence of landsliding at steeper slopes appeared from the flattening of the overall negative trend in the higher slope range (S > 0·30) of the integrated dataset. Comparing the threshold lines of our datasets to the average trend lines through data found in literature revealed that vegetation type and cover could better explain differences in topographical thresholds level than climatic conditions. In cultivated fields, soil structure and moisture conditions, as determined by the rainfall distribution, are critical factors influencing topographical thresholds rather than daily rainfall amounts of the gully‐initiating events. In rangelands, vegetation cover at the time of incision appears to be the most important factor differentiating between topographical thresholds, overruling the effect of average annual rainfall amounts. Soil texture and rock fragment cover contributed little to the explanation of the relative threshold levels. Differences in regression slopes (b) between the S–A relationships found in this study have been attributed to the soil characteristics in the different study areas, determining the relative importance of subsurface flow and Hortonian overland flow. Sedimentation slopes where both ephemeral and permanent gullies end were generally high because of the high rock fragment content of the transported sediment. A positive relationship was found between the rock fragment content at the apex of the sedimentation fan and the slope of the soil surface at this location. Copyright © 2000 John Wiley & Sons, Ltd.     

Runoff production processes in small alpine catchments within the unconsolidated Pleistocene sediments of the Lainbach area (Upper Bavaria)

In the Lainbach catchment, unconsolidated Pleistocene moraine sediments are widely distributed. Because of the great natural risk of floods, together with extreme loads of sediments, investigations of runoff production processes have been conducted in this area. At hillslope scale three test sites with different states of soil development and vegetation cover were instrumented with V‐shaped weirs, precipitation gauges and measurement devices for electrical conductivity (EC) of discharge water. The EC has been used as a geochemical tracer for hydrograph separation, since the statistical relationship between content of dissolved Ca²⁺, Mg²⁺ cations and EC is highly significant for different stages of runoff. This method allows hydrograph separation at high temporal resolution for both the rising and falling limb of the hydrograph. The following results of the investigations can be resumed. If relief conditions are similar, the effectiveness of runoff production decreases with an increasing density of vegetation cover. The runoff delivery ratio decreases as well as the peaks of runoff. In contrast, concentration times of hillslope catchments are equal, even if vegetation cover is of great density and soils are well developed. As a reason for the short reaction times, different runoff production processes have been detected. On bare ground, infiltration excess overland flow intensified by surface sealing processes is the main source for quick runoff. On hillslopes well covered by vegetation, translatory flow processes indicated by soil water with high solute contents force a rapid runoff reaction only a few minutes after rainfall has begun. It is to be assumed that translatory flow is a runoff production process typical for hillslopes covered by vegetation in a steep alpine relief. By means of the areal distribution of the topographic index, concentration of runoff production on a small part of the catchment has been demonstrated for hillslopes densely covered by vegetation. The investigations have shown that there is a lack of studies on runoff production processes in steep alpine relief, as well as a deficit of methods to quantify hydraulic properties of coarse‐grained soils with a wide grain size distribution. Copyright © 2003 John Wiley & Sons, Ltd.     

A simulation model for unified interrill erosion and rill erosion on hillslopes

A mathematical model was developed for simulating runoff generation and soil erosion on hillslopes. The model is comprised of three modules: one for overland flow, one for soil infiltration, and one for soil erosion including rill erosion and interrill erosion. Rainfall and slope characteristics affecting soil erosion on hillslopes were analysed. The model results show that the slope length and gradient, time distribution rainfall, and distribution of rills have varying influence on soil erosion. Erosion rate increases nonlinearly with increase in the slope length; a long slope length leads to more serious erosion. The effect of the slope gradient on soil erosion can be both positive and negative. Thus, there exists a critical slope gradient for soil erosion, which is about 45° for the rate of erosion at the end of the slope and about 25° for the accumulated erosion. Copyright © 2005 John Wiley & Sons, Ltd.