Modelling soil erosion with a downscaled landscape evolution model

The measurement and prediction of soil erosion is important for understanding both natural and disturbed landscape systems. In particular numerical models of soil erosion are important tools for managing landscapes as well as understanding how they have evolved over time. Over the last 40 years a variety of methods have been used to determine rates of soil loss from a landscape and these can be loosely categorized into empirical and physically based models. Alternatively, physically based landscape evolution models (LEMs) have been developed that provide information on soil erosion rates at much longer decadal or centennial scales, over large spatial scales and examine how they may respond to environmental and climatic changes. Both soil erosion LEMs are interested in similar outcomes (landscape development and sediment delivery) yet have quite different methodologies and parameterizations. This paper applies a LEM (the CAESAR model) for the first time at time and space scales where soil erosion models have largely been used. It tests the ability of the LEM to predict soil erosion on a 30 m experimental plot on a trial rehabilitated landform in the Northern Territory, Australia. It then continues to discuss the synergies and differences between soil erosion and LEMs. The results demonstrate that once calibrated for the site hydrology, predicted suspended sediment and bedload yields from CAESAR show a close correspondence in both volume and timing of field measured data. The model also predicts, at decadal scales, sediment loads close to that of field measured data. Findings indicate that the small‐scale drainage network that forms within these erosion plots is an important control on the timing and magnitude of sediment delivery. Therefore, it is important to use models that can alter the DEM to reflect changing topography and drainage network as well as having a greater emphasis on channel processes. Copyright © 2012 John Wiley & Sons, Ltd. and Commonwealth of Australia     

A method for estimating spatially variable seepage and hydraulic conductivity in channels with very mild slopes

Infiltration along ephemeral channels plays an important role in groundwater recharge in arid regions. A model is presented for estimating spatial variability of seepage due to streambed heterogeneity along channels based on measurements of streamflow‐front velocities in initially dry channels. The diffusion‐wave approximation to the Saint‐Venant equations, coupled with Philip's equation for infiltration, is connected to the groundwater model MODFLOW and is calibrated by adjusting the saturated hydraulic conductivity of the channel bed. The model is applied to portions of two large water delivery canals, which serve as proxies for natural ephemeral streams. Estimated seepage rates compare well with previously published values. Possible sources of error stem from uncertainty in Manning's roughness coefficients, soil hydraulic properties and channel geometry. Model performance would be most improved through more frequent longitudinal estimates of channel geometry and thalweg elevation, and with measurements of stream stage over time to constrain wave timing and shape. This model is a potentially valuable tool for estimating spatial variability in longitudinal seepage along intermittent and ephemeral channels over a wide range of bed slopes and the influence of seepage rates on groundwater levels. Copyright © 2012 John Wiley & Sons, Ltd.