Effects of sediment load on hydraulics of overland flow on steep slopes

Eroded sediment may have significant effects on the hydraulics of overland flow, but few studies have been performed to quantify these effects on steep slopes. This study investigated the potential effects of sediment load on Reynolds number, Froude number, flow depth, mean velocity, Darcy–Weisbach friction coefficient, shear stress, stream power, and unit stream power of overland flow in a sand‐glued hydraulic flume under a wide range of hydraulic conditions and sediment loads. Slope gradients were varied from 8·7 to 34·2%, unit flow rates from 0·66 to 5·26×10^?3 m² s^?1, and sediment loads from 0 to 6·95 kg m^?1 s^?1. Both Reynolds number (Re) and Froude number (Fr) decreased as sediment load increased, implying a decrease in flow turbulence. This inverse relationship should be considered in modeling soil erosion processes. Flow depth increased as sediment load increased with a mean value of 1·227 mm, caused by an increase in volume of sediment‐laden flow (contribution 62·4%) and a decrease in mean flow velocity (contribution 37·6%). The mean flow velocity decreased by up to 0·071 m s^?1 as sediment load increased. The Darcy–Weisbach friction coefficient (f) increased with sediment load, showing that the total energy consumption increased with sediment load. The effects of sediment load on f depended on flow discharge: as flow discharge increased, the influence of sediment load on f decreased due to increased flow depth and reduced relative roughness. Flow shear stress and stream power increased with sediment load, on average, by 80·5% and 60·2%, respectively; however, unit stream power decreased by an average of 11·1% as sediment load increased. Further studies are needed to extend and apply the insights obtained under these controlled conditions to real‐world overland flow conditions. Copyright © 2010 John Wiley & Sons, Ltd.     

Global heat flow: A new look

A global heat flow map has been derived from existing observations supplemented in areas without data by an empirical predictor based on tectonic setting and age. In continental areas the predictor is based on the observed correlation of heat flow with age of last tectono-thermal event, and in oceanic regions on the observed relation of heat flow to age of ocean floor. The predictor was used to assign mean heat flow values to 5° × 5° grid areas on the globe, weighted according to the relative area of tectonic provinces represented. A spherical harmonic analysis to degree 12 of the heat flow field yields a mean value of 59 mW m^?2, a rms residual of 13 mW m^?2, and an amplitude spectrum which decreases gradually and almost monotonically fromn = 1. The spherical harmonic representation of the heat flow field is free of the unreal distortions which have characterized earlier analyses based on a geographically sparse data set. Areas with residuals greater than 15 mW m^?2 comprise less than 19% of the area of the globe, thus indicating that most heat flow provinces have characteristic dimensions adequately represented in a 12-degree analysis.     

An experimental study on semi-brittle and plastic rheology of Panzhihua gabbro

Rheological properties of the crust and upper mantle are essential data that are needed in modelling the mechanical behaviour of the shallow part of the earth. The importance of such data has been seen in discussions about the strength profile of continen…     

Thick lava flows of Karisimbi Volcano, Rwanda: insights from SIR-C interferometric topography

 We use a digital elevation model (DEM) derived from interferometrically processed SIR-C radar data to estimate the thickness of massive trachyte lava flows on the east flank of Karisimbi Volcano, Rwanda. The flows are as long as 12 km and average 40–60 m (up to >140 m) in thickness. By calculating and subtracting a reference surface from the DEM, we derived a map of flow thickness, which we used to calculate the volume (up to 1 km³ for an individual flow, and 1.8 km³ for all the identified flows) and yield strength of several flows (23–124 kPa). Using the DEM we estimated apparent viscosity based on the spacing of large folds (1.2×10¹² to 5.5×10¹² Pa s for surface viscosity, and 7.5×10¹⁰ to 5.2×10¹¹ Pa s for interior viscosity, for a strain interval of 24 h). We use shaded-relief images of the DEM to map basic flow structures such as channels, shear zones, and surface folds, as well as flow boundaries. The flow thickness map also proves invaluable in mapping flows where flow boundaries are indistinct and poorly expressed in the radar backscatter and shaded-relief images. Received: 6 September 1997 / Accepted: 15 May 1998     

The significance of channel recharge rates for estimating debris‐flow magnitude and frequency

The quantification of debris‐flow hazard requires estimates of debris‐flow frequency and magnitude. Several methods have been proposed to determine the probable volume of future debris flows from a given basin, but most have neglected to account for debris recharge rates over time, which may lead to underestimation of debris‐flow volumes in basins with rare debris flows. This paper deals with the determination of debris recharge rates in debris‐flow channels based on knowledge of debris storage and the elapsed time since the last debris flow. Data are obtained from coastal British Columbia and a relation is obtained across a sample of basins with similar terrain and climatic conditions. For Rennell Sound on the west coast of the Queen Charlotte Islands, the power‐law relation for area‐normalized recharge rate, R_t, versus elapsed time, t_e was R_t = 0·23t_e^?0·58 with an explained variance of 75 per cent. A difference in recharge rates may exist between creeks in logged and unlogged forested terrain. The power function for undisturbed terrain was R_t = 0·20t_e^?0·49, while the function for logged areas was R_t = 0·30t_e^?0·77. This result suggests that for the same elapsed time since the last debris flow, clearcut gullies tend to recharge at a slower rate than creeks in old growth forest. This finding requires verification, particularly for longer elapsed times since debris flow, but would have important implications for forest resource management in steep coastal terrain. This study demonstrates that commonly used encounter probability equations are inappropriate for recharge‐limited debris flow channels. Copyright © 2005 John Wiley & Sons, Ltd.     

The linkage between velocity patterns and sediment entrainment in a forced‐pool and riffle unit

A field‐based project was initiated in order to characterize velocities and sediment entrainment in a forced‐pool and riffle sequence. Three‐dimensional velocities and turbulence intensities were measured with an acoustic Doppler velocimeter at 222 different points at three similar flows that averaged approximately 4·35 m³ s⁻¹ within a large pool–riffle unit on North Saint Vrain Creek, Colorado. Sediment‐sorting patterns were observed with the introduction of 500 tracer particles painted according to initial seeding location. Tracer particles moved sporadically during a 113 day period in response to the annual snowmelt peak flow, which reached a maximum level of 14·8 m³ s⁻¹. Velocity data indicate high instantaneous velocities and turbulence levels in the centre of pools. Patterns of sediment deposition support the notion that stream competence is higher in the pool than the downstream riffle. Flow convergence around a large channel constriction appears to play a major role in multiple processes that include helical flow development and sediment routing, and backwater development with low velocities and turbulence levels above the constriction that may locally limit sediment supply. Jet flow, flow separation, vortex scour and turbulence generation enhance scour in the centre of pools. Ultimately, multiple processes appear to play some role in maintenance of this forced pool and the associated riffle. Copyright © 2008 John Wiley & Sons, Ltd.     

A spreadsheet-based technique (Lotus 1-2-3) for separating tropical forest storm hydrographs using Hewlett and Hibbert's slope

The inclined line separation technique of Hewlett and Hibbert has been widely adopted to separate delayed flow from the total stream storm runoff. Presented here is the application of the technique to highly responsive storm hydrographs using a personal computer method based on a Lotus 1-2-3 spreadsheet. Using discharge measurements (in m³ s⁻¹), catchment area (in km²) and time (in Julian days), the separation slope is adjusted on the monitor screen until the precise time at which the end of quickflow as storm runoff gives way to delayed flow may be established. The application of the inclined line method is compared with other separation techniques applied to the same dataset. The annual stream quickflow runoff for the study catchment was calculated by the four different separating lines — (i) best-fit curve, (ii) N-day after peak, (iii) inclined line and (iv) horizontal line — was 250, 312, 368, and 588 mm, amounting to 33, 31, 51 and 78 per cent respectively of the annual total stream runoff. Separation of flow by computer spreadsheet methods may be consistently applied throughout a dataset and therefore have a comparative advantage over more arbitrary techniques. © 1997 John Wiley & Sons, Ltd.     

The dynamics of a channel-fed lava flow on Pico Partido volcano, Lanzarote

A short length of channel on Pico Partido volcano, Lanzarote, provides us the opportunity to examine the dynamics of lava flowing in a channel that extends over a sudden break in slope. The 1–2-m-wide, 0.5–2-m-deep channel was built during the 1730–1736 eruptions on Lanzarote and exhibits a sinuous, well-formed channel over a steep (11° slope) 100-m-long proximal section. Over-flow units comprising smooth pahoehoe sheet flow, as well as evidence on the inner channel walls for multiple (at least 11) flow levels, attest to unsteady flow in the channel. In addition, superelevation is apparent at each of the six bends along the proximal channel section. Superelevation results from banking of the lava as it moves around the bend thus causing preferential construction of the outer bank. As a result, the channel profile at each bend is asymmetric with an outer bank that is higher than the inner bank. Analysis of superelevation indicates flow velocities of ~8 m s^–1. Our analysis of the superelevation features is based on an inertia-gravity balance, which we show is appropriate, even though the down-channel flow is in laminar flow. We use a viscosity-gravity balance model, together with the velocities calculated from superelevation, to obtain viscosities in the range 25–60 Pa s (assuming that the lava behaved as a Newtonian liquid). Estimated volume fluxes are in the range 7–12 m³ s^–1. An apparent down-flow increase in derived volume flux may have resulted from variable supply or bulking up of the flow due to vesiculation. Where the channel moves over a sharp break in slope and onto slopes of ~6°, the channel becomes less well defined and widens considerably. At the break of slope, an elongate ridge extends across the channel. We speculate that this ridge was formed as a result of a reduction in velocity immediately below the break of slope to allow deposition of entrained material or accretion of lava to the channel bed as a result of a change in flow regime or depth.     

Terrestrial heat flow in Japan

Terrestrial heat flow, Q=K×ΔT/ΔZ cal/cm² sec has been determined at 51 localities (39 on land and 12 in the sea) in and around the Japanese Islands. The average values of observed heat flow in land and sea are 1.53µ cal/cm²sec and 1.48µcal/cm²sec respectively. These value do not differ greatly from the world’s averages. The outstanding features of the heat flow distribution are as follows:a) High heat flow region (Q>2.0µcal/cm²sec) exists in the Inner Zone of the Honshu Arc. This region of high heat flow is more distinct in the northeastern Japan than in the southwestern Japan.b) The High heat flow region seems to extend, through the Fossa Magna area, down to the Izu-Mariana Arc.c) It is also probable that a similar high heat flow zone exists in the inner side of the Kurile Arc.d) These zones of high heat flow precisely coincide with the zones of the Cenozoic orogeny in the area concerned.e) Far off the coast of the northeastern Japan, the area at about 150° E may be a high heat flow region.f) Low heat flow (Q<1.0µcal/cm²sec) prevails in the Pacific coast side of the northeastern Japan and in the oceanic area directly east of it, including the area of the Japan Trench.g) The region bounded by the above mentioned high and low heat flow regions has heat flow which is more or less normal. Based on these measurements, a « steady state ” temperature distribution in the crust has been calculated for each of the above regions of high, low and intermediate heat flow, and it was found that there is a large temperature differences between the bottom of the crust of the high and low heat flow regions: the temperature at the Moho boundary in the high heat flow regions should be as high as some 800～1000°C (d=27 km), whereas that under the low heat flow region should be only about 200°C (d=23 km). The high general temperature at the Moho under the high heat flow region seems to favor a production of magma in the upper mantle. Calculated Moho temperatures disfavor the hypothesis that the Moho boundary is due to phase transition.     

GIS‐assisted modelling for debris flow hazard assessment based on the events of May 1998 in the area of Sarno,Southern Italy: Part I. Maximum run‐out

Based on the debris flow events that occurred in May 1998 in the area of Sarno, Southern Italy, this paper presents an approach to simulate debris flow maximum run‐out. On the basis of the flow source areas and an average thickness of 1·2 m of the scarps, we estimated debris flow volumes of the order of 10⁴ and 10⁵ m³. Flow mobility ratios (ΔH/L) derived from the x, y, z coordinates of the lower‐most limit of the source areas (i.e. apex of the alluvial fan) and the distal limit of the flows ranged between 0·27 and 0·09. We performed regression analyses that showed a good correlation between the estimated flow volumes and mobility ratios. This paper presents a methodology for predicting maximum run‐out of future debris flow events, based on the developed empirical relationship. We implemented the equation that resulted from the calibration as a set of GIS macros written in Visual Basic for Applications (VBA) and running within ArcGIS. We carried out sensitivity analyses and observed that hazard mapping with this methodology should attempt to delineate hazard zones with a minimum horizontal resolution of 0·4 km. The developed procedure enables the rapid delineation of debris flow maximum extent within reasonable levels of uncertainty, it incorporates sensitivities and it facilitates hazard assessments via graphic user interfaces and with modest computing resources. Copyright © 2007 John Wiley & Sons, Ltd.     

Effect of bed topography on soil aggregates transport by rill flow

After its formation, a rill may remain in the field for months, often receiving lower flow rates than the formative discharge. The objective of this work was to evaluate the rill flow transport capacity of soil aggregates at discharges unable to erode the rill, and to analyse the influence of the rill macro‐roughness on this transport process. A non‐erodible rill was built in which roughness was reproduced in detail. In order to assess only the rill macro‐roughness, a flat channel with a similar micro‐roughness to that in the rill replica was built. Rill and channel experiments were carried out at a slope of 8 and at six discharges (8·3 × 10^?5 to 5·2 × 10^?4 m³ s^?1) in the rill, and eight discharges (1·6 × 10^?5 to 5·2 × 10^?4 m³ s^?1) in the channel. Non‐erodible aggregates of three sizes (1–2, 3–5 and 5–10 mm) were released at the inlet of the rill/channel. The number of aggregates received at the outlet was registered. The number and position of the remaining aggregates along the rill/channel were also determined. The rill flow was a major sediment transport mechanism only during the formation of the rill, as during that period the power of the flow was great enough to overcome the influence of the macro‐roughness of the rill bed. At lower discharges the transport capacity in the previously formed rill was significantly less than that in the flat channel under similar slope and discharge. This was determined to be due to local slowing of flow velocities at the exit of rill pools. Copyright © 2006 John Wiley & Sons, Ltd.     

Assessment of transmission losses and groundwater recharge from runoff events in a wadi under shortage of data on lateral inflow,Negev, Israel

A hydrological–lithostratigraphical model was developed for assessment of transmission losses and groundwater recharge from runoff events in arid water courses where hydrological and meteorological records are incomplete. Water balance equations were established for reaches between hydrometric stations. Because rainfall and tributary flow data are scarce, lateral inflow, which is an essential component of the water balance equation, could not be estimated directly. The solution was obtained by developing a method which includes a hydrological–lithostratigraphical analogy. This is based on the following assumptions: (a) runoff resulting from a given rainfall event is related to the watershed surface lithology; (b) for a given event, the spatial distribution of runoff reflects the distribution of rainfall: and (c) transmission losses are uniquely related to the total inflow to the reach. The latter relationship, called the loss function, and the water balance equation comprise a model which simultaneously assesses lateral inflow and transmission losses for runoff events recorded at the terminal stations. The model was applied to three reaches of the arid Nahal Tsin in Israel. In this case study, the transmission losses were of the same order of magnitude as the flow at the major hydrometric stations. The losses were subdivided into channel moistening, which subsequently evaporates, and deep percolation, which recharges groundwater. For large runoff events, evaporation was substantially smaller than the losses. The mean annual recharge of groundwater from runoff events in the Tsin watershed was 4·1×10⁶ m³, while the mean annual flow volume at the major stations ranged from 0·6 to 1·5×10⁶ m³. Once in 100 years, the annual recharge may be seven times higher than the mean annual value, but the recharge during most years is very small. Copyright © 1999 John Wiley & Sons, Ltd.     

Heat flow from eastern Panama and northwestern Colombia

Heat flows were determined at 12 sites in four distinct areas between longitude 77° and 80°W in eastern Panama and northwestern Colombia. Evidently, most of the region is underlain by mafic oceanic crust so that the crustal radiogenic component of heat flow is very small (～ 0.1 μcal cm^?2 sec^?1). Low heat-flow values (～ 0.7 μcal cm^?2 sec^?1) in northwestern Colombia may reflect thermal transients associated with shallow subduction. The normal values (～ 1) at about 78°W are consistent with the mean heat flow from the western Caribbean and the Gulf of Mexico. At 80°W, a fairly high value of 1.8 may define the easterly limit of thermal transients due to Cenozoic volcanic activity in Central America.     

One‐ and two‐dimensional modelling of overland flow in semiarid shrubland,Jornada basin,New Mexico

Two distributed parameter models, a one‐dimensional (1D) model and a two‐dimensional (2D) model, are developed to simulate overland flow in two small semiarid shrubland watersheds in the Jornada basin, southern New Mexico. The models are event‐based and represent each watershed by an array of 1‐m² cells, in which the cell size is approximately equal to the average area of the shrubs. Each model uses only six parameters, for which values are obtained from field surveys and rainfall simulation experiments. In the 1D model, flow volumes through a fixed network are computed by a simple finite‐difference solution to the 1D kinematic wave equation. In the 2D model, flow directions and volumes are computed by a second‐order predictor–corrector finite‐difference solution to the 2D kinematic wave equation, in which flow routing is implicit and may vary in response to flow conditions. The models are compared in terms of the runoff hydrograph and the spatial distribution of runoff. The simulation results suggest that both the 1D and the 2D models have much to offer as tools for the large‐scale study of overland flow. Because it is based on a fixed flow network, the 1D model is better suited to the study of runoff due to individual rainfall events, whereas the 2D model may, with further development, be used to study both runoff and erosion during multiple rainfall events in which the dynamic nature of the terrain becomes an important consideration. Copyright © 2006 John Wiley & Sons, Ltd.     

Characteristics of concentrated flow hydraulics for rangeland ecosystems: implications for hydrologic modeling

Concentrated flow is often the dominant source of water erosion following disturbance on rangelands. Because of the lack of studies that explain the hydraulics of concentrated flow on rangelands, cropland‐based equations have typically been used for rangeland hydrology and erosion modeling, leading to less accurate predictions due to different soil and vegetation cover characteristics. This study investigates the hydraulics of concentrated flow using unconfined field experimental data over diverse rangeland landscapes within the Great Basin Region, United States. The results imply that the overall hydraulics of concentrated flow on rangelands differ significantly from those of cropland rills. Concentrated flow hydraulics on rangelands are largely controlled by the amount of cover or bare soil and hillslope angle. New predictive equations for concentrated flow velocity (R² = 0·47), hydraulic friction (R² = 0·52), and width (R² = 0·4) representing a diverse set of rangeland environments were developed. The resulting equations are applicable across a wide span of ecological sites, soils, slopes, and vegetation and ground cover conditions and can be used by physically‐based rangeland hydrology and erosion models to estimate rangeland concentrated flow hydraulic parameters. Published in 2011. This article is a US Government work and is in the public domain in the USA.     

Estimation of water velocities in boreholes using a new passive flowmeter

Abstract

Preferential flow pathways in a fractured aquifer may yield abrupt reductions of the water velocity in a well. We propose a new device for measuring low (5–13 cm d^-1) velocities in wells originating from fractures at different depths. The presented flowmeter has been applied in a well in the Bari (southern Italy) fractured aquifer. In the same well, the horizontal flowmeter velocity (9.6 cm d^-1) at 0.5 m depth was compared with velocity (8 cm d^-1) derived from a field tracer test, providing a value 16.5% higher. Moreover, the flowmeter measurements at 1.5 m depth gave a horizontal velocity of 7.2 cm d^-1, which is 11% less than water flow velocity estimated from the field test. The new flowmeter implements the tracer point-dilution method in a plastic (PVC) pipe by causing the water flow to pass through an artificial filter. Laboratory calibration tests have confirmed the good performance of the proposed flowmeter technique, even for water flow up to 300 cm d^-1. The flowmeter was sensitive to 0.1 cm d^-1, with a detection limit of 1.5 cm d^-1, i.e. half the measurable flow velocity of existing flowmeters in wells.

Editor D. Koutsoyiannis; Associate editor S. Grimaldi     

Landslide and tsunami hazard at Yate volcano,Chile as an example of edifice destruction on strike-slip fault zones

The edifice of Yate volcano, a dissected stratocone in the Andean Southern Volcanic Zone, has experienced multiple summit collapses throughout postglacial time restricted to sectors NE and SW of the summit. The largest such historic event occurred on 19th February 1965 when ～6.1–10?×?10⁶ m³ of rock and ice detached from 2,000-m elevation to the SW of the summit and transformed into a debris flow. In the upper part of the flow path, velocities are estimated to have reached 40 m s^?1. After travelling 7,500 m and descending 1,490 m, the flow entered an intermontane lake, Lago Cabrera. A wavemaker of estimated volume 9?±?3?×?10⁶ m³ generated a tsunami with an estimated amplitude of 25 m and a run-up of ～60 m at the west end of the lake where a settlement disappeared with the loss of 27 lives. The landslide followed 15 days of unusually heavy summer rain, which may have caused failure by increasing pore water pressure in rock mechanically weathered through glacial action. The preferential collapse directions at Yate result from the volcano’s construction on the dextral strike-slip Liquiñe-Ofqui fault zone. Movement on the fault during the lifetime of the volcano is thought to have generated internal instabilities in the observed failure orientations, at ～10° to the fault zone in the Riedel shear direction. This mechanically weakened rock may have led to preferentially orientated glacial valleys, generating a feedback mechanism with collapse followed by rapid glacial erosion, accelerating the rate of incision into the edifice through repeated landslides. Debris flows with magnitudes similar to the 1965 event are likely to recur at Yate, with repeat times of the order of 10² years. With a warming climate, increased glacial meltwater due to snowline retreat and increasing rain, at the expense of snow, may accelerate rates of edifice collapse, with implications for landslide hazard and risk at glaciated volcanoes, in particular those in strike-slip tectonic settings where orientated structural instabilities may exist.     

Flow separation on Zongo Glacier,Cordillera Real,Bolivia

Meltwaters collected from the proglacial stream escaping from Zongo Glacier (2·1 km²), Bolivia (16°S), have been monitored in order to analyse the internal drainage system of an Andean glacier. Electrical conductivity has been measured sporadically between February 1995 and March 1996, during 16 one-day field surveys, under various meteorological conditions in summer and winter. The mixing-model technique based on the electrical conductivity is used for a quantitative separation of discharge which is derived from continuous water level registration. Tracer experiments (mainly uranine dye and NaCl salt) have been carried out from March to June 1997 to obtain information about the internal drainage system. In the tropical Andes, accumulation only occurs in austral summer, whereas ablation occurs throughout the year and is higher during the accumulation season, between November and March. The assumptions involved in the use of mixing models for analysis of glacial drainage structure are applicable for tropical glaciers because glacial conduits do not suffer complete closure, and are permanently supplied by meltwaters, even in wintertime. Two components of discharge are separated: an englacial flow originating from surface meltwater which is routed without chemical enrichment, and offering low electrical conductivity; and a subglacial one routed in contact with bedrock or sediments showing high ionic concentrations. Electrical conductivity of meltwater varies diurnally, inversely to discharge fluctuations. According to this behaviour, total discharge is mainly formed by the englacial component. The drainage structures for englacial and subglacial flow have to be widely interconnected, as indicated by diurnal variations of the subglacial discharge. Comparison of hydrograph separation based on conductivity and on ¹⁸O isotope confirms that the subglacial flow is influenced by surface melting. A hydrograph separation of the subglacial flow is proposed, between a diurnal variable component, composed of water coming from the englacial network, and a base flow, which may vary seasonally. The dye tracing experiments confirm the drainage complexity of Zongo Glacier and demonstrate the interest of identifying three main drainage components. © 1998 John Wiley & Sons, Ltd.     

Nondimensional sediment transport capacity of sand soils and its response to parameter in the Loess Plateau of China

Soil erosion is a major contributor to land degradation in the Loess Plateau in China. To clarify the sediment transport capacity of overland flow influenced by hydraulic parameters, such as shear stress, sand shear stress (hydraulic gradient partition method and hydraulic radius partition method), mean flow velocity, Froude number, stream power, and unit stream power, indoor experiments with eight-unit-width flow discharges from 0.0667 × 10⁻³ to 0.3333 × 10⁻³ m²·s⁻¹, six slope gradients from 3.49 to 20.79%, and two kinds of sand soils (d₅₀ = 0.17 and 0.53 mm) were systematically investigated. A nondimensional method was adopted in data processing. Results showed that there was a partition phenomenon of relation curves because of the different median grain diameters. The correlation between the nondimensional stream power and nondimensional sediment transport capacity was the highest, followed by the correlation between the nondimensional unit stream power and nondimensional sediment transport capacity. However, there was a poor correlation between the flow intensity indices of velocity category and nondimensional sediment transport capacity. Nondimensional stream power, nondimensional unit stream power, and nondimensional shear stress could predict sediment transport capacity well. Ignoring the partition phenomenon of the relation curves, stream power could be used to predict sediment transport capacity, with a coefficient of determination of .85. Furthermore, a general flow intensity index was obtained to predict sediment transport capacity of overland flow. Finally, an empirical formula for predicting sediment transport capacity with a coefficient of determination of .90 was established by multiple regression analyses based on the general flow intensity index. During the analysis between measured sediment transport capacities in present study and predicted values based on Zhang model, Mahmoodabadi model, and Wu model, it was found that these three models could not accurately predict sediment transport capacities of this study because different models are estimated on the basis of different experimental conditions.     

Evolution of overland flow connectivity in bare agricultural plots

Soil surface roughness not only delays overland flow generation but also strongly affects the spatial distribution and concentration of overland flow. Previous studies generally aimed at predicting the delay in overland flow generation by means of a single parameter characterizing soil roughness. However, little work has been done to find a link between soil roughness and overland flow dynamics. This is made difficult because soil roughness and hence overland flow characteristics evolve differently depending on whether diffuse or concentrated erosion dominates. The present study examined whether the concept of connectivity can be used to link roughness characteristics to overland flow dynamics. For this purpose, soil roughness of three 30‐m² tilled plots exposed to natural rainfall was monitored for two years. Soil micro‐topography was characterized by means of photogrammetry on a monthly basis. Soil roughness was characterized by the variogram, the surface stream network was characterized by network‐based indices and overland flow connectivity was characterized by Relative Surface Connection function (RSCf) functional connectivity indicator. Overland flow hydrographs were generated by means of a physically‐based overland flow model based on 1‐cm resolution digital elevation models. The development of eroded flow paths at the soil surface not only reduced the delay in overland flow generation but also resulted in a higher continuity of high flow velocity paths, an increase in erosive energy and a higher rate of increase of the overland flow hydrograph. Overland flow dynamics were found to be highly correlated to the RSCf characteristic points. By providing information regarding overland flow dynamics, the RSCf may thus serve as a quantitative link between soil roughness and overland flow generation in order to improve the overland flow hydrograph prediction. Copyright © 2016 John Wiley & Sons, Ltd.