Modelling long‐term (300 ka) upland catchment response to multiple lava damming events

Landscapes respond in complex ways to external drivers such as base level change due to damming events. In this study, landscape evolution modelling was used to understand and analyse long‐term catchment response to lava damming events. PalaeoDEM reconstruction of a small Turkish catchment (45 km²) that endured multiple lava damming events in the past 300 ka, was used to derive long‐term net erosion rates. These erosion rates were used for parameter calibration and led to a best fit parameter set. This optimal parameter set was used to compare net erosion landscape time series of four scenarios: (i) no uplift and no damming events; (ii) no uplift and three damming events; (iii) uplift and no damming events; and (iv) uplift and three damming events. Spatial evolution of net erosion and sediment storage of scenario (iii) and (iv) were compared. Simulation results demonstrate net erosion differences after 250 000 years between scenarios with and without dams. Initially, trunk gullies show less net erosion in the scenario with damming events compared with the scenario without damming events. This effect of dampened erosion migrates upstream to smaller gullies and local slopes. Finally, an intrinsic incision pulse in the dam scenario results in a higher net erosion of trunk gullies while decoupled local slopes are still responding to the pre‐incision landscape conditions. Sediment storage differences also occur on a 100 ka scale. These differences behaved in a complex manner owing to different timings of the migration of erosion and sediment waves along the gullies for each scenario. Although the specific spatial and temporal sequence of erosion and deposition events is sensitive to local parameters, this model study shows the manner in which past short‐lived events like lava dams have long‐lasting effects on catchment evolution. Copyright © 2014 John Wiley & Sons, Ltd.     

Modelling centennial sediment waves in an eroding landscape – catchment complexity

Sediment flux dynamics in fluvial systems have often been related to changes in external drivers of topography, climate or land cover. It is well known that these dynamics are non‐linear. Recently, model simulations of fluvial activity and landscape evolution have suggested that self‐organization in landscapes can also cause internal complexity in the sedimentary record. In this contribution one particular case of self‐organization is explored in the Sabinal field study area, Spain, where several dynamic zones of sedimentation and incision are observed along the current river bed. Whether these zones can be caused by internal complexity was tested with landscape evolution model (LEM) LAPSUS (Landscape Process Modelling at Multi‐dimensions and Scales). During various 500 year simulations, zones of sedimentation appear to move upstream and downstream in eroding river channels (‘waves’). These waves are visualized and characterized for a range of model settings under constant external forcing, and the self‐organizing process behind their occurrence is analysed. Results indicate that this process is not necessarily related to simplifications in the model and is more generic than the process of bed‐armouring that has recently been recognized as a cause for complexity in LEM simulations. We conclude that autogenic sediment waves are the result of the spatial propagation in time of feedbacks in local transport limited (deposition) and detachment limited (erosion) conditions. Copyright © 2014 John Wiley & Sons, Ltd.     

Modelling the long-term geomorphic response to check dam failures in an alpine channel with CAESAR-Lisflood

Globally, between 1950 and 2011 nearly 80,000 debris flow fatalities occurred in densely populated regions in mountainous terrain. Mitigation of these hazards includes the construction of check dams, which limit coarse sediment transport and in the European Alps number in the 100,000s. Check dam functionality depends on periodic, costly maintenance, but maintenance is not always possible and check dams often fail. As such, there is a need to quantify the long-term (10–100 years) geomorphic response of rivers to check dam failures. Here, for the first time, a landscape evolution model (CAESAR-Lisflood) driven by a weather generator is used to replicate check dam failures due to the lack of maintenance, check dam age, and flood occurrence. The model is applied to the Guerbe River, Switzerland, a pre-Alpine catchment containing 73 check dams that undergo simulated failure. Also presented is a novel method to calibrate CAESAR-Lisflood's hydrological component on this ungauged catchment. Using 100-year scenarios of check dam failure, the model indicates that check dam failures can produce 8 m of channel erosion and a 322% increase in sediment yield. The model suggests that after check dam failure, channel erosion is the remobilization of deposits accumulated behind check dams, and, after a single check dam failure channel equilibrium occurs in five years, but after many check dam failures channel equilibrium may not occur until 15 years. Overall, these findings support the continued maintenance of check dams.     

Fluvial and slope‐wash erosion of soil‐mantled landscapes: detachment‐ or transport‐limited?

Many numerical landform evolution models assume that soil erosion by flowing water is either purely detachment‐limited (i.e. erosion rate is related to the shear stress, power, or velocity of the flow) or purely transport‐limited (i.e. erosion/deposition rate is related to the divergence of shear stress, power, or velocity). This paper reviews available data on the relative importance of detachment‐limited versus transport‐limited erosion by flowing water on soil‐mantled hillslopes and low‐order valleys. Field measurements indicate that fluvial and slope‐wash modification of soil‐mantled landscapes is best represented by a combination of transport‐limited and detachment‐limited conditions with the relative importance of each approximately equal to the ratio of sand and rock fragments to silt and clay in the eroding soil. Available data also indicate that detachment/entrainment thresholds are highly variable in space and time in many landscapes, with local threshold values dependent on vegetation cover, rock‐fragment armoring, surface roughness, soil texture and cohesion. This heterogeneity is significant for determining the form of the fluvial/slope‐wash erosion or transport law because spatial and/or temporal variations in detachment/entrainment thresholds can effectively increase the nonlinearity of the relationship between sediment transport and stream power. Results from landform evolution modeling also suggest that, aside from the presence of distributary channel networks and autogenic cut‐and‐fill cycles in non‐steady‐state transport‐limited landscapes, it is difficult to infer the relative importance of transport‐limited versus detachment‐limited conditions using topography alone. Copyright © 2011 John Wiley & Sons, Ltd.     

The impact of land use change and check‐dams on catchment sediment yield

Extensive land use changes have occurred in many areas of SE Spain as a result of reforestation and the abandonment of agricultural activities. Parallel to this the Spanish Administration spends large funds on hydrological control works to reduce erosion and sediment transport. However, it remains untested how these large land use changes affect the erosion processes at the catchment scale and if the hydrological control works efficiently reduce sediment export. A combination of field work, mapping and modelling was used to test the influence of land use scenarios with and without sediment control structures (check‐dams) on sediment yield at the catchment scale. The study catchment is located in SE Spain and suffered important land use changes, increasing the forest cover 3‐fold and decreasing the agricultural land 2·5‐fold from 1956 to 1997. In addition 58 check‐dams were constructed in the catchment in the 1970s accompanying reforestation works. The erosion model WATEM‐SEDEM was applied using six land use scenarios: land use in 1956, 1981 and 1997, each with and without check‐dams. Calibration of the model provided a model efficiency of 0·84 for absolute sediment yield. Model application showed that in a scenario without check dams, the land use changes between 1956 and 1997 caused a progressive decrease in sediment yield of 54%. In a scenario without land use changes but with check‐dams, about 77% of the sediment yield was retained behind the dams. Check‐dams can be efficient sediment control measures, but with a short‐lived effect. They have important side‐effects, such as inducing channel erosion downstream. While also having side‐effects, land use changes can have important long‐term effects on sediment yield. The application of either land use changes (i.e. reforestation) or check‐dams to control sediment yield depends on the objective of the management and the specific environmental conditions of each area. Copyright © 2008 John Wiley & Sons, Ltd.     

Relations between rainfall–runoff‐induced erosion and aeolian deposition at archaeological sites in a semi‐arid dam‐controlled river corridor

Process dynamics in fluvial‐based dryland environments are highly complex with fluvial, aeolian, and alluvial processes all contributing to landscape change. When anthropogenic activities such as dam‐building affect fluvial processes, the complexity in local response can be further increased by flood‐ and sediment‐limiting flows. Understanding these complexities is key to predicting landscape behavior in drylands and has important scientific and management implications, including for studies related to paleoclimatology, landscape ecology evolution, and archaeological site context and preservation. Here we use multi‐temporal LiDAR surveys, local weather data, and geomorphological observations to identify trends in site change throughout the 446‐km‐long semi‐arid Colorado River corridor in Grand Canyon, Arizona, USA, where archaeological site degradation related to the effects of upstream dam operation is a concern. Using several site case studies, we show the range of landscape responses that might be expected from concomitant occurrence of dam‐controlled fluvial sand bar deposition, aeolian sand transport, and rainfall‐induced erosion. Empirical rainfall‐erosion threshold analyses coupled with a numerical rainfall–runoff–soil erosion model indicate that infiltration‐excess overland flow and gullying govern large‐scale (centimeter‐ to decimeter‐scale) landscape changes, but that aeolian deposition can in some cases mitigate gully erosion. Whereas threshold analyses identify the normalized rainfall intensity (defined as the ratio of rainfall intensity to hydraulic conductivity) as the primary factor governing hydrologic‐driven erosion, assessment of false positives and false negatives in the dataset highlight topographic slope as the next most important parameter governing site response. Analysis of 4+ years of high resolution (four‐minute) weather data and 75+ years of low resolution (daily) climate records indicates that dryland erosion is dependent on short‐term, storm‐driven rainfall intensity rather than cumulative rainfall, and that erosion can occur outside of wet seasons and even wet years. These results can apply to other similar semi‐arid landscapes where process complexity may not be fully understood. Published 2015. This article is a U.S. Government work and is in the public domain in the USA     

Landscapes on the edge: examining the role of climatic interactions in shaping coastal watersheds using a coastal–terrestrial landscape evolution model

Incised coastal gullies (ICGs) are dynamic features found at the terrestrial‐coastal interface. Their geomorphic evolution is driven by the interactions between processes of fluvial knickpoint migration and coastal cliff erosion. Under scenarios of future climate change the frequency and magnitude of the climatological drivers of both terrestrial (fluvial and hillslope) and coastal (cliff erosion) processes are likely to change, with an adjunct impact on these types of coastal features. Here we explore the response of an incised coastal gully to changes in both terrestrial and coastal climate in order to elucidate the key process interactions which drive ICG evolution. We modify an extant landscape evolution model, CHILD, to incorporate processes of soft‐cliff erosion. This modified version, termed the Coastal‐Terrestrial‐CHILD (CT‐CHILD) model, is then employed to explore the interactions between changing terrestrial and coastal driving forces on the future evolution of an ICG found on the south‐west Isle of Wight, UK. It was found that the magnitude and frequency of storm events will play a key role in determining the future trajectory of ICGs, highlighting a need to understand the role of event sequencing in future projections of landscape evolution. Furthermore, synergistic (positive) and antagonistic (negative) interactions were identified between coastal and terrestrial parameters, such as wave height intensity and precipitation duration, which act to modulate the impact of changes in any one parameter. Of note was the role played by wave height intensity in driving coastal erosion, which was found to play a more important role than sea‐level rise in determining rates of coastal erosion. This highlights the need for a greater focus on wave height in studies of soft‐cliff erosion. Copyright © 2014 John Wiley & Sons, Ltd.     

Catchment response to lava damming: integrating field observation,geochronology and landscape evolution modelling

Combining field reconstruction and landscape evolution modelling can be useful to investigate the relative role of different drivers on catchment response. The Geren Catchment (~45 km²) in western Turkey is suitable for such a study, as it has been influenced by uplift, climate change and lava damming. Four Middle Pleistocene lava flows (⁴⁰Ar/³⁹Ar‐ dated from 310 to 175 ka) filled and dammed the Gediz River at the Gediz–Geren confluence, resulting in base‐level fluctuations of the otherwise uplift‐driven incising river. Field reconstruction and luminescence dating suggest fluvial terraces in the Geren Catchment are capped by Middle Pleistocene aggradational fills. This showed that incision of the Geren trunk stream has been delayed until the end of MIS 5. Subsequently, the catchment has responded to base‐level lowering since MIS 4 by 30 m of stepped net incision. Field reconstruction left us with uncertainty on the main drivers of terrace formation. Therefore, we used landscape evolution modelling to investigate catchment response to three scenarios of base‐level change: (i) uplift with climate change (rainfall and vegetation based on arboreal pollen); (ii) uplift, climate change and short‐lived damming events; (iii) uplift, climate and long‐lived damming events. Outputs were evaluated for erosion–aggradation evolution in trunk streams at two different distances from the catchment outlet. Climate influences erosion–aggradation activity in the catchment, although internal feedbacks influence timing and magnitude. Furthermore, lava damming events partly control if and where these climate‐driven aggradations occur. Damming thus leaves a legacy on current landscape evolution. Catchment response to long‐duration damming events corresponds best with field reconstruction and dating. The combination of climate and base level explains a significant part of the landscape evolution history of the Geren Catchment. By combining model results with fieldwork, additional conclusions on landscape evolution could be drawn. Copyright © 2016 John Wiley & Sons, Ltd.     

Drainage basin sensitivity to tectonic and climatic forcing: implications of a stochastic model for the role of entrainment and erosion thresholds

Long‐term average rates of channel erosion and sediment transport depend on the frequency–magnitude characteristics of ?ood ?ows that exceed an erosion threshold. Using a Poisson model for rainfall and runoff, analytical solutions are developed for average rates of stream incision and sediment transport in the presence of such a threshold. Solutions are derived and numerically tested for three erosion/transport formulas: the Howard–Kerby shear‐stress incision model, the Bridge–Dominic sediment transport model, and a generic shear‐stress sediment transport model. Results imply that non‐linearity resulting from threshold effects can have a ?rst‐order impact on topography and patterns of dynamic response to tectonic and climate forcing. This non‐linearity becomes signi?cant when fewer than about half of ?ood events are capable of detaching rock or sediment. Predicted morphology and uplift‐gradient scaling is more closely consistent with observations and laboratory experiments than conventional slope‐linear or shear‐linear erosion laws. These results imply that particle detachment thresholds are not details that can be conveniently ignored in long‐term landscape evolution models. Copyright © 2004 John Wiley & Sons, Ltd.     

A demonstration of the importance of bedload transport for fluvial bedrock erosion and knickpoint propagation

We provide field evidence for the role of bedload in driving fluvial incision and knickpoint propagation. Using aerial photographs, field surveys, and hydrological data, we constrain the incision history of a bedrock gorge 1200 m long and up to 20 m deep cut by Da'an River in western Taiwan. This reach of the river experienced 10 m of uplift during the 1999 Chi‐Chi earthquake. For five years following the earthquake, bedload was prevented from entering the uplift zone, the knickpoint was static and little incision took place. Bedload transport across the uplift zone resumed in 2004, initiating extremely rapid incision, with 620 m of knickpoint propagation and up to 20 m of downcutting by 2008. This change highlights the relative inefficiency of suspended sediment and the dominant role of bedload as a tool for fluvial erosion and knickpoint propagation. Once bedload tools became available, knickpoint propagation was influenced by geological structure, lithology, and drainage organization. In particular, a change in dip of the sandstone beds at the site caused a decrease of knickpoint propagation velocity. Copyright © 2012 John Wiley & Sons, Ltd.     

Raindrop‐impact‐induced erosion processes and prediction: a review

Raindrop‐impact‐induced erosion is initiated when detachment of soil particles from the surface of the soil results from an expenditure of raindrop energy. Once detachment by raindrop impact has taken place, particles are transported away from the site of the impact by one or more of the following transport processes: drop splash, raindrop‐induced flow transport, or transport by flow without stimulation by drop impact. These transport processes exhibit varying efficiencies. Particles that fall back to the surface as a result of gravity produce a layer of pre‐detached particles that provides a degree of protection against the detachment of particles from the underlying soil. This, in turn, influences the erodibility of the eroding surface. Good understanding of rainfall erosion processes is necessary if the results of erosion experiments are to be properly interpreted. Current process‐based erosion prediction models do not deal with the issue of temporal variations in erodibility during a rainfall event or variabilities in erodibility associated with spatial changes in dominance of the transport processes that follow detachment by drop impact. Although more complex erosion models may deal with issues like this, their complexity and high data requirement may make them unsuitable for use as general prediction tools. Copyright © 2005 John Wiley & Sons, Ltd.     

Sediment connectivity: a framework for understanding sediment transfer at multiple scales

A major challenge for geomorphologists is to scale up small‐magnitude processes to produce landscape form, yet existing approaches have been found to be severely limited. New ways to scale erosion and transfer of sediment are thus needed. This paper evaluates the concept of sediment connectivity as a framework for understanding processes involved in sediment transfer across multiple scales. We propose that the concept of sediment connectivity can be used to explain the connected transfer of sediment from a source to a sink in a catchment, and movement of sediment between different zones within a catchment: over hillslopes, between hillslopes and channels, and within channels. Using fluvial systems as an example we explore four scenarios of sediment connectivity which represent end‐members of behaviour from fully linked to fully unlinked hydrological and sediment connectivity. Sediment‐travel distance – when combined with an entrainment parameter reflecting the frequency–magnitude response of the system – maps onto these end‐members, providing a coherent conceptual model for the upscaling of erosion predictions. This conceptual model could be readily expanded to other process domains to provide a more comprehensive underpinning of landscape‐evolution models. Thus, further research on the controls and dynamics of travel distances under different modes of transport is fundamental. Copyright © 2014 John Wiley & Sons, Ltd.     

Modeling the impact of dam removal on channel evolution and sediment delivery in a multiple dam setting

Dam removal can generate geomorphic disturbances, including channel bed and bank erosion and associated abrupt/pulsed release and downstream transfer of reservoir sediment, but the type and rate of geomorphic response often are hard to predict. The situation gets even more complex in systems which have been impacted by multiple dams and a long and complex engineering history. In previous studies one-dimensional (1-D) models were used to predict aspects of post-removal channel change. However, these models do not consider two-dimensional (2-D) effects of dam removal such as bank erosion processes and lateral migration. In the current study the impacts of multiple dams and their removal on channel evolution and sediment delivery were modeled by using a 2-D landscape evolution model (CAESAR-Lisflood) focusing on the following aspects: patterns, rates, and processes of geomorphic change and associated sediment delivery on annual to decadal timescales. The current modeling study revealed that geomorphic response to dam removal (i.e., channel evolution and associated rates of sediment delivery) in multiple dam settings is variable and complex in space and time. Complexity in geomorphic system response is related to differences in dam size, the proximity of upstream dams, related buffering effects and associated rates of upstream sediment supply, and emerging feedback processes as well as to the presence of channel stabilization measures. Modeled types and rates of geomorphic adjustment, using the 2-D landscape evolution model CAESAR-Lisflood, are similar to those reported in previous studies. Moreover, the use of a 2-D method showed some advantages compared to 1-D models, generating spatially varying patterns of erosion and deposition before and after dam removal that provide morphologies that are more readily comparable to field data as well as features like the lateral re-working of past reservoir deposits which further enables the maintenance of sediment delivery downstream.     

Modelling sediment dynamics due to hillslope–river interactions: incorporating fluvial behaviour in landscape evolution model LAPSUS

Landscape evolution models (LEMs) simulate the three‐dimensional development of landscapes over time. Different LEMs have different foci, e.g. erosional behaviour, river dynamics, the fluvial domain, hillslopes or a combination. LEM LAPSUS is a relatively simple cellular model operating on timescales of centuries to millennia and using annual timesteps that has had a hillslope focus. Our objective was to incorporate fluvial behaviour in LAPSUS without changing the existing model equations. The model should be able to reproduce alternating aggradation and incision in the floodplains of catchments, depending on simulated conditions. Testing was done using an artificial digital elevation model (DEM) and a demonstration of the ability for fluvial simulation was performed for a real landscape (Torrealvilla catchment, southeast Spain). Model equations to calculate sediment dynamics and water routing were similar for both hillslope and fluvial conditions, but different parameter values were used for these domains, defined based on annual discharge. Parameters changing between the domains are convergence factor p, which is used in the multiple flow algorithm to route water, and discharge and gradient exponents m and n, used in transport capacity calculations. Erodibility and ‘sedimentability’ factors K and P were changed between cold (little vegetation, high erodibility) and warm conditions (more vegetation, lower erodibility). Results show that the adapted parameters reproduced alternating aggradation – due to divergent flow in the floodplain and sediment supply under cold conditions – and incision due to reduced sediment supply and resulting clean water erosion during simulated warm conditions. The simulated results are due to interactions between hillslopes and floodplains, as the former provide the sediments that are deposited in the latter. Similar behaviour was demonstrated when using the real DEM. Sensitivity and resolution analysis showed that the model is sensitive to changes in m, n and p and that model behaviour is influenced by DEM resolution. Copyright © 2012 John Wiley & Sons, Ltd.     

Hillslope response to knickpoint migration in the Southern Appalachians: implications for the evolution of post‐orogenic landscapes

The southern Appalachians represent a landscape characterized by locally high topographic relief, steep slopes, and frequent mass movement in the absence of significant tectonic forcing for at least the last 200 Ma. The fundamental processes responsible for landscape evolution in a post‐orogenic landscape remain enigmatic. The non‐glaciated Cullasaja River basin of south‐western North Carolina, with uniform lithology, frequent debris flows, and the availability of high‐resolution airborne lidar DEMs, is an ideal natural setting to study landscape evolution in a post‐orogenic landscape through the lens of hillslope–channel coupling. This investigation is limited to channels with upslope contributing areas >2.7 km², a conservative estimate of the transition from fluvial to debris‐flow dominated channel processes. Values of normalized hypsometry, hypsometric integral, and mean slope vs elevation are used for 14 tributary basins and the Cullasaja basin as a whole to characterize landscape evolution following upstream knickpoint migration. Results highlight the existence of a transient spatial relationship between knickpoints present along the fluvial network of the Cullasaja basin and adjacent hillslopes. Metrics of topography (relief, slope gradient) and hillslope activity (landslide frequency) exhibit significant downstream increases below the current position of major knickpoints. The transient effect of knickpoint‐driven channel incision on basin hillslopes is captured by measuring the relief, mean slope steepness, and mass movement frequency of tributary basins and comparing these results with the distance from major knickpoints along the Cullasaja River. A conceptual model of area–elevation and slope distributions is presented that may be representative of post‐orogenic landscape evolution in analogous geologic settings. Importantly, the model explains how knickpoint migration and channel–hillslope coupling is an important factor in tectonically‐inactive (i.e. post‐orogenic) orogens for the maintenance of significant relief, steep slopes, and weathering‐limited hillslopes. Copyright © 2011 John Wiley & Sons, Ltd.     

Predicting reach‐specific properties of fluvial terraces to guide future fieldwork. A case study for the Late Quaternary River Allier (France) with the FLUVER2 model

Numerical models have not yet systematically been used to predict properties of fluvial terrace records in order to guide fieldwork and sampling. This paper explores the potential of the longitudinal profile model FLUVER2 to predict testable field properties of the relatively well‐studied, Late Quaternary Allier system in France. For the Allier terraces an overlapping ¹⁴C and U‐series chronology as well as a record of ¹⁰Be erosion rates exist. The FLUVER2 modelling exercise is focused on the last 50 ka of the upper Allier reach because for this location and period the constraints of the available dating techniques are tightest. A systematic calibration based on terrace occurrence and thicknesses was done using three internal parameters related to (1) the sediment erodibility; (2) the sediment transport distance; and (3) the sediment supply derived from the surrounding landscape. As external model inputs, the best available, reconstructed, tectonic, climatic and base‐level data were used. Calibrated model outputs demonstrate a plausible match with the existing fluvial record. Validation of model output was done by comparing the modelled and measured timing of aggradation and incision phases for the three locations. The modelled range of landscape erosion rates showed a reasonably good match with existing erosion rate estimates derived from ¹⁰Be measurements of fluvial sands. The quasi‐validated model simulation was subsequently used to make new testable predictions about the timing and location of aggradation and erosion phases for three locations along the Allier river. The validated simulations predict that along the Allier, reach‐specific dynamics of incision and aggradation, related to the variations in sediment supply by major tributaries, cause relevant differences in the local fluvial terrace stratigraphy. Copyright © 2016 John Wiley & Sons, Ltd.     

Impact of vegetation on erosion: Insights from the calibration and test of a landscape evolution model in alpine badland catchments

While it is well recognized that vegetation can affect erosion, sediment yield and, over longer timescales, landform evolution, the nature of this interaction and how it should be modeled is not obvious and may depend on the study site. In order to develop quantitative insight into the magnitude and nature of the influence of vegetation on catchment erosion, we build a landscape evolution model to simulate erosion in badlands, then calibrate and evaluate it against sediment yield data for two catchments with contrasting vegetation cover. The model couples hillslope gravitational transport and stream alluvium transport. Results indicate that hillslope transport processes depend strongly on the vegetation cover, whereas stream transport processes do not seem to be affected by the presence of vegetation. The model performance in prediction is found to be higher for the denuded catchment than for the reforested one. Moreover, we find that vegetation acts on erosion mostly by reducing soil erodibility rather than by reducing surface runoff. Finally, the methodology we propose can be a useful tool to evaluate the efficiency of previous revegetation operations and to provide guidance for future restoration work. © 2019 John Wiley & Sons, Ltd.     

Laboratory soil piping and internal erosion experiments: evaluation of a soil piping model for low‐compacted soils

Mechanistic models have been proposed for soil piping and internal erosion on well‐compacted levees and dams, but limited research has evaluated these models in less compacted (more erodible) soils typical of hillslopes and streambanks. This study utilized a soil box (50 cm long, 50 cm wide and 20 cm tall) to conduct constant‐head, soil pipe and internal erosion experiments for two soils (clay loam from Dry Creek and sandy loam from Cow Creek streambanks) packed at uniform bulk densities. Initial gravimetric moisture contents prior to packing were 10, 12 and 14% for Dry Creek soil and 8, 12, and 14% for Cow Creek soil. A 1‐cm diameter rod was placed horizontally along the length of the soil bed during packing and carefully removed after packing to create a continuous soil pipe. A constant head was maintained at the inflow end. Flow rates and sediment concentrations were measured from the pipe outlet. Replicate submerged jet erosion tests (JETs) were conducted to derive erodibility parameters for repacked samples at the same moisture contents. Flow rates from the box experiments were used to calibrate the mechanistic model. The influence of the initial moisture content was apparent, with some pipes (8% moisture content) expanding so fast that limited data was collected. The mechanistic model was able to estimate equivalent flow rates to those observed in the experiments, but had difficulty matching observed sediment concentrations when the pipes rapidly expanded. The JETs predicted similar erodibility coefficients compared to the mechanistic model for the more erodible cases but not for the less erodible cases (14% moisture content). Improved models are needed that better define the changing soil pipe cross‐section during supply‐ and transport‐limited internal erosion, especially for piping through lower compacted (more erodible) soils as opposed to more well‐compacted soils resulting from constructing levees and dams. Copyright © 2013 John Wiley & Sons, Ltd.     

Erosion prediction on unpaved mountain roads in northern Thailand: validation of dynamic erodibility modelling using KINEROS2

The event‐ and physics‐based KINEROS2 runoff/erosion model for predicting overland flow generation and sediment production was applied to unpaved mountain roads. Field rainfall simulations conducted in northern Thailand provided independent data for model calibration and validation. Validation shows that KINEROS2 can be parameterized to simulate total discharge, sediment transport and sediment concentration on small‐scale road plots, for a range of slopes, during simulated rainfall events. The KINEROS2 model, however, did not accurately predict time‐dependent changes in sediment output and concentration. In particular, early flush peaks and the temporal decay in sediment output were not predicted, owing to the inability of KINEROS2 to model removal of a surface sediment layer of finite depth. After 15–20 min, sediment transport declines as the supply of loose superficial material becomes depleted. Modelled erosion response was improved by allowing road erodibility to vary during an event. Changing the model values of erosion detachment parameters in response to changes in surface sediment availability improved model accuracy of predicted sediment transport by 30–40%. A predictive relationship between road erodibility ‘states’ and road surface sediment depth is presented. This relationship allows implementation of the dynamic erodibility (DE) method to events where pre‐storm sediment depth can be estimated (e.g., from traffic usage variables). Copyright © 2001 John Wiley & Sons, Ltd.