Effects of gradient,distance, curvature and aspect on steep burned and unburned hillslope soil erosion and deposition

Wildfire denudes vegetation and impacts chemical and physical soil properties, which can alter hillslope erosion rates. Post‐wildfire erosion can also contribute disproportionately to long‐term erosion rates and landscape evolution. Post‐fire hillslope erosion rates remain difficult to predict and document at the hillslope scale. Here we use ²¹⁰Pb_aex (lead‐210 mineral‐adsorbed excess) inventories to describe net sediment erosion on steep, convex hillslopes in three basins (unburned, moderately and severely burned) in mountainous central Idaho. We analyzed nearly 300 soil samples for ²¹⁰Pb_aex content with alpha spectrometry and related net sediment erosion to burn severity, aspect, gradient, curvature and distance from ridgetop. We also tested our data against models for advective, linear and non‐linear diffusive erosion. Statistically lower net soil losses on north‐ versus south‐facing unburned hillslopes suggest that greater vegetative cover and soil cohesion on north‐facing slopes decrease erosion. On burned hillslopes, erosion differences between aspects were less apparent and net erosion was more variable, indicating that vegetation influences erosion magnitude and fire drives erosion variability. We estimated net soil losses throughout the length of unburned hillslopes, including through a footslope transition to concave form. In contrast, on burned hillslopes, the subtle shift from convex to concave form was associated with deposition of a post‐fire erosion pulse. Such overall patterns of erosion and deposition are consistent with predictions from a non‐linear diffusion equation. This finding also suggests that concave sections of overall convex hillslopes affect post‐disturbance soil erosion and deposition. Despite these patterns, no strong relationships were evident between local net soil losses and gradient, curvature, distance from ridgetop, or erosion predicted with advection or diffusion equations. The observed relationship between gradient and erosion is therefore likely more complex or stochastic than often described theoretically, especially over relatively short timescales (60–100 years). Copyright © 2016 John Wiley & Sons, Ltd.     

Asymmetric hillslope erosion following wildfire in Fourmile Canyon,Colorado

Infrequent, high‐magnitude events cause a disproportionate amount of sediment transport on steep hillslopes, but few quantitative data are available that capture these processes. Here we study the influence of wildfire and hillslope aspect on soil erosion in Fourmile Canyon, Colorado. This region experienced the Fourmile Fire of 2010, strong summer convective storms in 2011 and 2012, and extreme flooding in September 2013. We sampled soils shortly after these events and use fallout radionuclides to trace erosion on polar‐ and equatorial‐facing burned slopes and on a polar‐facing unburned slope. Because these radionuclides are concentrated in the upper decimeter of soil, soil inventories are sensitive to erosion by surface runoff. The polar‐facing burned slope had significantly lower cesium‐137 (¹³⁷Cs) and lead‐210 (²¹⁰Pb) inventories (p < 0.05) than either the polar‐facing unburned slope or equatorial‐facing burned slope. Local slope magnitude does not appear to control the erosional response to wildfire, as relatively gently sloping (~20%) polar‐facing positions were severely eroded in the most intensively burned area. Field evidence and soil profile analyses indicate up to 4 cm of local soil erosion on the polar‐facing burned slope, but radionuclide mass balance indicates that much of this was trapped nearby. Using a ¹³⁷Cs‐based erosion model, we find that the burned polar‐facing slope had a net mean sediment loss of 2 mm (~1 kg m^?2) over a one to three year period, which is one to two orders of magnitude higher than longer‐term erosion rates reported for this region. In this part of the Colorado Front Range, strong hillslope asymmetry controls soil moisture and vegetation; polar‐facing slopes support significantly denser pine and fir stands, which fuels more intense wildfires. We conclude that polar‐facing slopes experience the most severe surface erosion following wildfires in this region, indicating that landscape‐scale aridity can control the geomorphic response of hillslopes to wildfires. Copyright © 2018 John Wiley & Sons, Ltd.     

The sensitivity of hillslope bedrock erosion to precipitation

Decoupling the impacts of climate and tectonics on hillslope erosion rates is a challenging problem. Hillslope erosion rates are well known to respond to changes in hillslope boundary conditions (e.g. channel incision rates) through their dependence on soil thickness, and precipitation is an important control on soil formation. Surprisingly though, compilations of hillslope denudation rates suggest little precipitation sensitivity. To isolate the effects of precipitation and boundary condition, we measured rates of soil production from bedrock and described soils on hillslopes along a semi‐arid to hyperarid precipitation gradient in northern Chile. In each climate zone, hillslopes with contrasting boundary conditions (actively incising channels versus non‐eroding landforms) were studied. Channel incision rates, which ultimately drive hillslope erosion, varied with precipitation rather than tectonic setting throughout the study area. These precipitation‐dependent incision rates are mirrored on the hillslopes, where erosion shifts from relatively fast and biologically‐driven to extremely slow and salt‐driven as precipitation decreases. Contrary to studies in humid regions, bedrock erosion rates increase with precipitation following a power law, from ～1 m Ma^?1 in the hyperarid region to ～40 m Ma^?1 in the semi‐arid region. The effect of boundary condition on soil thickness was observed in all climate zones (thicker soils on hillslopes with stable boundaries compared to hillslopes bounded by active channels), but the difference in bedrock erosion rates between the hillslopes within a climate region (slower erosion rates on hillslopes with stable boundaries) decreased as precipitation decreased. The biotic‐abiotic threshold also marks the precipitation rate below which bedrock erosion rates are no longer a function of soil thickness. Our work shows that hillslope processes become sensitive to precipitation as life disappears and the ability of the landscape to respond to tectonics decreases. Copyright © 2010 John Wiley & Sons, Ltd.     

Quantifying hillslope erosion rates and processes for a coastal California landscape over varying timescales

The Earth's surface erodes by processes that occur over different spatial and temporal scales. Both continuous, low‐magnitude processes as well as infrequent, high‐magnitude events drive erosion of hilly soil‐mantled landscapes. To determine the potential variability of erosion rates we applied three independent, field‐based methods to a well‐studied catchment in the Marin Headlands of northern California. We present short‐term, basin‐wide erosion rates determined by measuring pond sediment volume (40 years) and measured activities of the fallout nuclides ¹³⁷Cs and ²¹⁰Pb (40–50 years) for comparison with long‐term (>10 ka) rates previously determined from in situ‐produced cosmogenic ¹⁰Be and ²⁶Al analyses. In addition to determining basin‐averaged rates, ¹³⁷Cs and ²¹⁰Pb enable us to calculate point‐specific erosion rates and use these rates to infer dominant erosion processes across the landscape. When examined in the context of established geomorphic transport laws, the correlations between point rates of soil loss from ¹³⁷Cs and ²¹⁰Pb inventories and landscape morphometry (i.e. topographic curvature and upslope drainage area) demonstrate that slope‐driven processes dominate on convex areas while overland flow processes dominate in concave hollows and channels. We show a good agreement in erosion rates determined by three independent methods: equivalent denudation rates of 143 ± 41 m Ma^?1 from pond sediment volume, 136 ± 36 m Ma^?1 from the combination of ¹³⁷Cs and ²¹⁰Pb, and 102 ± 25 m Ma^?1 from ¹⁰Be and ²⁶Al. Such agreement suggests that erosion of this landscape is not dominated by extreme events; rather, the rates and processes observed today are indicative of those operating for at least the past 10 000 years. Copyright © 2006 John Wiley & Sons, Ltd.     

210Pb dating of scottish lake sediments,afforestation and accelerated soil erosion

Changes in the rate of soil erosion in lake catchments can be identified from changes in the rate of sediment accumulation in lakes. Here we compare recently afforested sites with non-afforested sites in the Galloway area of Southwest Scotland. We show that lakes with non-afforested catchments have slow, constant sediment accumulation rates, whereas lakes with recently afforested catchments have changes in accumulation that parallel the known history of afforestation. For Loch Grannoch the sediment accumulation rate increases from 0.1 cm yr^?1 to over 2 cm yr^?2 during the disturbance period. Data from L. Skerrow, however, suggest that the rate might decline to predisturbance levels after approximately 10 years as the forest canopy closes and drainage channels stabilize.     

Modelling the effect of soil burn severity on soil erosion at hillslope scale in the first year following wildfire in NW Spain

Fire severity is recognized as a key factor in explaining post‐fire soil erosion. However, the relationship between soil burn severity and soil loss has not been fully established until now. Sediment availability may also affect the extent of post‐fire soil erosion. The objective of this study was to determine whether soil burn severity, estimated by an operational classification system based on visual indicators, can significantly explain soil loss in the first year after wildfire in shrubland and other areas affected by crown fires in northwest (NW) Spain. An additional aim was to establish indicators of sediment availability for use as explanatory variables for post‐fire soil loss. For these purposes, we measured hillslope‐scale sediment production rates and site characteristics during the first year after wildfire in 15 experimental sites using 65 plots. Sediment yields varied from 0.2 Mg ha^?1 to 50.1 Mg ha^?1 and soil burn severity ranged from low (1.8) to very high (4.5) in the study period. A model that included soil burn severity, annual precipitation and a land use factor (as a surrogate for sediment availability) as explanatory variables reasonably explained the erosion losses measured during the first year after fire. Model validation confirmed the usefulness of this empirical model. The proposed empirical model could be used by forest managers to help evaluate erosion risks and to plan post‐fire stabilization activities. Copyright © 2015 John Wiley & Sons, Ltd.     

Post-fire soil erosion mitigation at the scale of swales using forest logging residues at a reduced application rate

Mulching with forest residues has proved to be highly effective in reducing post-fire soil losses at the plot scale. However, its effectiveness has not been quantified at the application rates that are typically used in operational post-fire land management (2–3 Mg ha^-1 using straw), as well as at scales larger than 100 m². The present study compared post-fire erosion rates for six convergent hillslopes or swales of 500 to 800 m², three of which were left untreated while the other three were mulched immediately after the fire with shredded eucalypt bark at a rate of 2.4 Mg ha^-1. Erosion rates were monitored at irregular intervals during the first three post-fire years, whilst ground cover was assessed yearly. Selected topsoil properties (0–2 cm) such as organic matter content and aggregate stability were determined at a single occasion – two years after the wildfire, for three micro-environments separately: bare soil, and under mulch/litter and vegetation. Soil losses on the untreated swales decreased with post-fire year from 2.2 to 0.4 and 0.11 Mg ha^-1 yr^-1 (respectively for the first, second and third post-fire years), while the mulched swales produced 84%, 77% and 38% less soil losses than the untreated swales. Soil losses also depended on slope aspect, with the north-facing swales producing less erosion than the west-facing ones. This could be linked to their significant differences in bare soil, vegetation and stone cover, or a combination thereof. The type of micro-environment also played a significant role in topsoil properties (stone content, bulk density, resistance to penetration/shear stress, porosity and organic matter content). The present results add to the increasing evidence that forest residues should be duly considered for operational post-fire land management. Forest residues were highly effective in reducing erosion from swales at application rates as low as the typical 2 Mg ha^-1 of post-fire straw mulch. © 2019 John Wiley & Sons, Ltd. © 2019 John Wiley & Sons, Ltd.     

Use of fallout tracers 7Be, 210Pb and 137Cs to distinguish the form of sub‐surface soil erosion delivering sediment to rivers in large catchments

Fallout radionuclides (FRNs) ¹³⁷Cs and ²¹⁰Pb are well established as tracers of surface and sub‐surface soil erosion contributing sediment to river systems. However, without additional information, it has not been possible to distinguish sub‐surface soil erosion sources. Here, we use the FRN ⁷Be (half‐life 53 days) in combination with ¹³⁷Cs and excess ²¹⁰Pb to trace the form of erosion contributing sediment in three large river catchments in eastern Australia; the Logan River (area 3700 km²), Bowen River (9400 km²) and Mitchell River (4700 km²). We show that the combination of ¹³⁷Cs, excess ²¹⁰Pb and ⁷Be can discriminate horizontally aligned sub‐surface erosion sources (rilled and scalded hillslopes and the floors of incised drainage lines and gully ‘badland’ areas) from vertical erosion sources (channel banks and gully walls). Specifically, sub‐surface sources of sediment eroded during high rainfall and high river flow events have been distinguished by the ability of rainfall‐derived ⁷Be to label horizontal soil surfaces, but not vertical. Our results indicate that in the two northern catchments, erosion of horizontal sub‐surface soil sources contributed almost as much fine river sediment as vertical channel banks, and several times the contribution of hillslope topsoils. This result improves on source discrimination provided previously and indicates that in some areas erosion of hillslope soils may contribute significantly to sediment yield, but not as topsoil loss. We find that in north‐eastern Australia, scalded areas on hillslopes and incising drainage lines may be sediment sources of comparable importance to vertical channel banks. Previous studies have used the combination of ¹³⁷Cs, excess ²¹⁰Pb and ⁷Be to estimate soils losses at the hillslope scale. Here, we show that with timely and judicious sampling of soil and sediment during and immediately after high flow events ⁷Be measurements can augment fallout ¹³⁷Cs and ²¹⁰Pb to provide important erosion source information over large catchments. Copyright © 2013 John Wiley & Sons, Ltd.     

How the spatial distribution of grass contributes to controlling hillslope erosion

The aim of this study was to investigate how the spatial distribution of grass influenced run-off and erosion from a hillslope with loess and cinnamon soils in the rocky area of Northern China. We set up a trial to test the two soils with different treatments, including bare soil (BS), grass strips on the upper (UGS) and lower (DGS) parts of the slope, grass cover over the entire slope (GS), and a grass carpet on the lower part of the slope (GC), under simulated rainfall conditions. The results showed that the run-off coefficients for the loess and cinnamon soils decreased by between 4% and 20% and by between 2% and 37%, respectively, when covered with grass. Grass spatial distribution had little effect on the run-off, but more effect on erosion than vegetation coverage degree. The most effective location of grass cover for decreasing hillslope erosion was at the foot, and the high efficiency was mainly due to controlling of rill formation and sediment deposition. The soil loss from GS, DGS, and GC on the loess and cinnamon soils was between 77% and 93% less and 55% and 80% less, respectively, compared with the loss from BS. However, the soil characteristics had little effect on soil erosion for well-vegetated slopes. The results highlight the importance of vegetation re-establishment at the foot of hillslope in controlling soil erosion.     

Unsteady soil erosion due to rainfall impact: a model of sediment sorting on the hillslope

A new method is presented for predicting sediment sorting associated with soil erosion by raindrop impact for non-equilibrium conditions. The form of soil erosion considered is that which results from raindrop impact in the presence of shallow overland flow itself where the flow is not capable of eroding sediment. The method specifically considers early time runoff and erosion when sediment leaving an eroding area is generally finer and thus may have a higher potential for transport of sorbed pollutants. The new mechanism described is the formation of a deposited layer on the soil surface, which is shown to lead to sediment sorting during an erosion event. The deposited layer is taken to have two roles in this process: to temporarily store sediment on the surface between successive trajectories, and to shield the underlying soil from erosive stresses. Equations describing the dynamics of the suspended sediment mixture and the deposited layer are developed. By integrating these equations over the length of eroding land element and over the duration of the erosion event, an event-based solution is proposed which predicts total sediment sorting over the event. This solution is shown to be consistent with experimentally observed trends in enrichment of fine sediment. Predictions using this approach are found to only partly explain measured enrichment for sets of experimental data for two quite different soils, but to be in poor agreement for an aridsol of dispersive character. It is concluded that the formation of the deposited layer is a significant mechanism in the enrichment of fine sediment and associated sorbed pollutants, but that processes in the dispersive soil are not as well described by the theory presented.     

Use of magnetic and radiometric measurements to investigate erosion and sedimentation in a British upland catchment

In response to the potential shortcomings of single-technique strategies in the investigation of erosion and sedimentation, a combined magnetic and radiometric (¹³Cs and ²¹⁰Pb) approach has been undertaken in the upland, watershed-lake system of Howden Reservoir, Derbyshire. By combining these techniques, some assessment of sediment sources and the erosion status of the catchment has been achieved. Alone, each approach would have been unable to determine unequivocally sediment provenance. Furthermore, the parallel use of these measurements has highlighted limitations and/or uncertainties in both the magnetic and ¹³⁷Cs techniques. These problems reflect the particular soil characteristics and drainage conditions of this upland catchment. Despite the documented severity of peat erosion in the region, Howden Reservoir has a mixture of sediment sources and a relatively moderate rate of sedimentation. Sediment yields (total 127·7 t km^?2 yr^?1 including organic fraction 31·3 t km^?2 yr^?1) are, however, higher than in other British upland areas.     

Testing the use of an image‐based technique to measure gully erosion at Sparacia experimental area

The first part of this investigation was aimed at testing the use of a three‐dimensional (3D) digital terrain model and a quasi‐tridimensional (2.5D) digital elevation model obtained by a large series of oblique images of eroded channels taken from consumer un‐calibrated and non‐metric cameras. For two closed earth channels having a different sinuosity, the ground measurement of some cross sections by a profilometer (P) was carried out and their real volume was also measured. The comparison among the three methods (3D, 2.5D, and P) pointed out that a limited underestimation of the total volume always occurs and that the 3D method is characterized by the minimum difference between measured and real volume. For this reason, 3D model can be used as benchmark. In the subsequent part of the investigation, the three ground measurement methods were applied for surveying of an ephemeral gully (EG) channel at the Sparacia area. The morphological and hydraulic variable values of the 24 surveyed cross sections determined by both 2.5D model and profilometer were compared. This comparison showed that the estimate error is generally less than ±10%. The EG measurements carried out by the three methods supported the applicability both of the empirical relationship between EG length and its eroded volume and the theoretical dimensionless relationship among the morphological variables describing the channelized erosion process. Finally, it was demonstrated that the effect of the distance interval on the EG volume measurement by 3D and 2.5D models is negligible for the investigated EG.     

Impact of model simplifications on soil erosion predictions: application of the GLUE methodology to a distributed event‐based model at the hillslope scale

In this paper, we analyse how the performance and calibration of a distributed event‐based soil erosion model at the hillslope scale is affected by different simplifications on the parameterizations used to compute the production of suspended sediment by rainfall and runoff. Six modelling scenarios of different complexity are used to evaluate the temporal variability of the sedimentograph at the outlet of a 60 m long cultivated hillslope. The six scenarios are calibrated within the generalized likelihood uncertainty estimation framework in order to account for parameter uncertainty, and their performance is evaluated against experimental data registered during five storm events. The Nash–Sutcliffe efficiency, percent bias and coverage performance ratios show that the sedimentary response of the hillslope in terms of mass flux of eroded soil can be efficiently captured by a model structure including only two soil erodibility parameters, which control the rainfall and runoff production of suspended sediment. Increasing the number of parameters makes the calibration process more complex without increasing in a noticeable manner the predictive capability of the model. Copyright © 2015 John Wiley & Sons, Ltd.     

210Pb和137Cs定年技术在湖泊沉积物中的应用与问题

湖泊沉积记录研究对了解流域侵蚀、揭示区域气候环境变化及人类活动影响、建立外生地球化学循环模式等均具有重要意义.放射性同位素~(210)Pb和~(137)Cs组合定年法是百年尺度湖泊沉积记录研究的首选方法,在过去几十年中得到了广泛应用,大力推动了湖泊沉积记录的高分辨率气候环境变化研究的开展.然而,如何理解和选择适用的计算模式来尽可能地提高定年准确度目前仍存在争议,这在一定程度上限制了该方法的应用.鉴于上述问题,本文基于本实验室近年来从青藏高原获取的若干支代表性湖芯的定年结果,从样品前处理、仪器及测试、数据分析及计算等全过程对~(210)Pb、~(137)Cs定年原理和方法进行了较为系统的阐述;在此基础上,提出了~(210)Pb定年样品封装原则,探讨了高纯锗γ谱仪在放射性活度测量中的影响因素,并建立了准确有效的比活度校正方法;明确指出了恒定初始浓度模式(constant initial concentration,CIC)可用于计算沉积物每层的年龄,而不是部分研究中认为的平均沉积速率,这对于深刻理解该定年模式并提供准确的定年结果具有重要的参考意义.同时,通过对比两种放射性核素的分布特征发现,相对于~(210)Pb而言,~(137)Cs更容易发生垂向迁移进而影响其时标准确性,这需引起研究者的充分重视,因此,当二者定年结果出现明显差异时,建议以~(210)Pb定年结果为准.     

Bank erosion of legacy sediment at the transition from vertical to lateral stream incision

Streambank erosion is a primary source of suspended sediments in many waterways of the US Atlantic Piedmont. This problem is exacerbated where banks are comprised of fine sediment produced by the intensive land use practices of early European settlers. A stream in this region, Richland Creek incises into banks comprised of three stratigraphic layers associated with historic land use: pre‐European settlement, early European agriculture and development, and water‐powered milldam operation. This study aims to identify the bank processes along a reach of Richland Creek that is eroding towards its pre‐disturbance elevation. The volume of material that has eroded along this stream since the milldam breached was calculated by differencing a reconstructed surface of the pond bed and an aerial lidar digital terrain model (DTM). Immediately downstream from the study reach, the channel is floored by bedrock and immediately upstream the rate of channel erosion approximately doubled along the longitudinal profile of Richland Creek, which indicate that the study reach spans the transition from a channel dominated by vertical incision in the upstream direction to horizontal widening in the downstream direction. The combined hydrometeorological conditions and dominant processes causing reach‐scale cut bank erosion were investigated with analyses of stream stage, precipitation, and streambank volumetric and surfaces change that was measured during nine terrestrial lidar surveys in 2010–2012. The spatial variability of erosion during a simulated precipitation event was examined in a field‐based experiment. Erosion was greatest where mill pond sediment columns detached along vertical desiccation and horizontal seepage cracks. This sediment accumulated on the bank toe throughout the study and was a source of readily‐entrained fine sediment contrary to the upper reaches where depositional accommodation space is more limited. Findings suggest that hotspots of sediment excavation progress upstream, indicating that restoration efforts should focus upon stabilizing banks at these locations. Copyright © 2015 John Wiley & Sons, Ltd.     

Numerical modeling of surface runoff and erosion due to moving rainstorms at the drainage basin scale

A physically-based distributed erosion model (MEFIDIS) was applied to evaluate the consequences of storm movement on runoff and erosion from the Alenquer basin in Portugal. Controlled soil flume laboratory experiments were also used to test the model. Nine synthetic circular storms were used, combining three storm diameters (0.5, 1 and 2 times the Alenquer basin’s axial length) with three speeds of storm movement (0.5, 1 and 2 m/s); storm intensities were synthesized in order to maintain a constant rainfall depth of 50 mm. The model was applied to storms moving downstream as well as upstream along the basin’s axis. In all tests, downstream-moving storms caused significantly higher peak runoff (56.5%) and net erosion (9.1%) than did upstream-moving storms. The consequences for peak runoff were amplified as the storm intensity increased. The hydrograph shapes were also different: for downstream-moving storms, runoff started later and the rising limb was steeper, whereas for upstream moving storms, runoff started early and the rising limb was less steep. Both laboratory and model simulations on the Alenquer basin showed that the direction of storm movement, especially in case of extreme rainfall events, significantly affected runoff and soil loss.     

Equifinality and uncertainty in physically based soil erosion models: application of the GLUE methodology to WEPP–the Water Erosion Prediction Project–for sites in the UK and USA

Despite the wealth of soil erosion models available for the prediction of both runoff and soil loss at a variety of scales, little quantification is made of uncertainty and error associated with model output. This in part reflects the need to produce unequivocal or optimal results for the end user, which will often be an unrealistic goal. This paper presents a conceptually simple methodology, Generalized Likelihood Uncertainty Estimation (GLUE), for assessing the degree of uncertainty surrounding output from a physically based soil erosion model, the Water Erosion Prediction Project (WEPP). The ability not only to be explicit about model error but also to evaluate future improvements in parameter estimation, observed data or scientific understanding is demonstrated. This approach is applied to two sets of soil loss/runoff plot replicates, one in the UK and one in the USA. Although it is demonstrated that observations can be largely captured within uncertainty bounds, results indicate that these uncertainty bounds are often wide, reflecting the need to qualify results that derive from ‘optimum’ parameter sets, and to accept the concept of equifinality within soil erosion models. Attention is brought to the problem of under‐prediction of large events/over‐prediction of small events, as an area where model improvements could be made, specifically in the case of relatively dry years. Finally it is proposed that such a technique of model evaluation be employed more widely within the discipline so as to aid the interpretation and understanding of complex model output. Copyright © 2000 John Wiley & Sons, Ltd.     

Application of lattice gas techniques to the study of sediment erosion and transport caused by laminar sheetflow

In order to gain a deeper understanding of the dynamics of erosion and sediment transport on hillslopes, it seems important to clarify the role of some basic mechanisms involved in these processes. While there is evidence that this cannot be done using the theoretical framework of river hydraulics, the use of numerical analysis could be of considerable help. The nature of the problem requires a technique capable of solving Navier–Stokes equations at low Reynolds number, with geometrically complex boundaries and solid particles moving inside the flow field. All these requirements make a novel method, known as lattice gas automaton LGA, a natural candidate for the study of the hydrodynamics of sheetflows. However, due to the recent introduction of this technique, there is a lack of a clear definition of its operational limits. Considering the case of a viscous sheetflow on an erodible rough boundary, we argue that by using LGA the stream Reynolds number can be increased only at the expense of a reduction of the boundary shearing stress. Accordingly, LGA cannot profitably be used to study the beginning of sediment motion and transport. On the other hand, a further evolution of LGA, known as the lattice Boltzmann method, seems highly promising for the numerical study of the erosion processes that eventually lead to drainage network evolution along hillslopes. © 1997 by John Wiley & Sons, Ltd.     

Seepage erosion and its implication to the formation of amphitheatre valley heads: A case study at Obara,Japan

Seepage erosion was investigated in an amphitheatre with a semicircular valley head, steep slopes, and a flat bottom developed in granodiorite hills at Obara, Aichi prefecture, Japan. A high sediment yield occurred where the measuring sites were located at the base of the landslide debris in the base of the convex slopes, whereas sediment outflows were small where the measuring sites were located at the base of the strong convex slopes. This implies that the seepage erosion was an effective agent for removal of debris deposited at the base of the slope. Small landslides can be found at the lower slopes within the area of the observed amphitheatre. The slope stability analysis and subsurface water observation of the lower slope suggest that the small landslides in this amphitheatre are due to over-steepened slopes, and relatively insensitive to subsurface water status. Colluvium in the flat valley bottom thinly covers the bedrock surface. Therefore the topography of the amphitheatre was found to be formed by parallel retreat of slopes by the repetition of basal seepage erosion and subsequent small landslides.