The use of a landscape simulator in the validation of the SIBERIA landscape evolution model: transient landforms

SIBERIA is a physically based model for the geomorphic evolution of landforms. It is essential that the SIBERIA model be tested or validated against controlled landform development. Previous studies have demonstrated that SIBERIA is able to simulate declining equilibrium landforms and in this paper we examine SIBERIA's ability to simulate landforms as they evolve to their declining equilibrium form. These landscapes are termed transient landforms. Landscapes generated by SIBERIA were compared to those produced by a physical model (experimental model landforms) at stages of evolution. Comparison of the experimental landscapes with the simulated landscapes using total mass, hypsometric curve, width function, cumulative area distribution and area–slope demonstrate that SIBERIA can simulate the experimental model landscape during development (i.e. transient landscapes). Copyright © 2002 John Wiley & Sons, Ltd.     

Smoothing of digital elevation models and the alteration of overland flow path length distributions

DEM smoothing is a common pre-processing technique used to remove undesirable roughness from a DEM. However, it is hypothesized that smoothing straightens and reduces the length of overland flow paths, which is an important factor controlling modelled time-to-peak flow. Currently, there is a lack of research comparing how different smoothing techniques alter the distribution of overland flow path length. Four low-pass filtering techniques were applied to three fine-resolution LiDAR DEMs of varying relief: the mean filter, the median filter, the Gaussian filter, and the feature-preserving DEM smoothing (FPDEMS) filter, each with different degrees of smoothing. Downslope-distance-to-stream distributions were then derived using D8 and D∞ flow directions and statistically compared to distributions derived from the unsmoothed DEM for each study site. The results indicate that the alteration of flow path length distributions as a result of smoothing is complex. Mean flow path lengths may decrease or increase in response to smoothing, depending on landscape relief and the derivation of flow directions, and generalized flow paths may become longer. The largest increase in mean flow path lengths was 19.2 m using the 21 × 21 median filter and D8 flow directions in the high-relief study site, relative to an unsmoothed mean length of 138.6 m in this site. The largest decrease in mean flow path length was 48.9 m using the 21 × 21 mean filter and D∞ flow directions in the low-relief study site, relative to an unsmoothed mean length of 290.9 m in this site. Furthermore, minimal flow path length alterations were achieved with the Gaussian filter when gentle smoothing is required, and with the FPDEMS filter when moderate to aggressive smoothing is required. These results suggest that an appropriate smoothing method should be chosen based on the relief of the landscape and the degree of smoothing required.     

The production of digital elevation models for experimental model landscapes

Automated digital photogrammetry was used to produce digital elevation models of experimental model landscapes under controlled laboratory conditions as part of a series of rainfall erosion experiments looking at the evolution of landforms in response to erosion. The method allowed the elevations of the experimental landscapes to be studied in great detail on a regular grid digital terrain map with relatively very little effort. Digital photogrammetry produced elevation data at a resolution of 6 mm with a standard deviation of 2·0 mm over an experimental catchment relief of approximately 200 mm; this resolution is considerably better than that achievable by conventional manual photogrammetry. The density of grid points was sufficiently high that small‐scale details such as knickpoints developing in channels were represented. The method can facilitate the study of both experimental and natural landscapes in great detail. Copyright © 2001 John Wiley & Sons, Ltd.     

Evaluating preprocessing methods of digital elevation models for hydrological modelling

With the introduction of high‐resolution digital elevation models, it is possible to use digital terrain analysis to extract small streams. In order to map streams correctly, it is necessary to remove errors and artificial sinks in the digital elevation models. This step is known as preprocessing and will allow water to move across a digital landscape. However, new challenges are introduced with increasing resolution because the effect of anthropogenic artefacts such as road embankments and bridges increases with increased resolution. These are problematic during the preprocessing step because they are elevated above the surrounding landscape and act as artificial dams. The aims of this study were to evaluate the effect of different preprocessing methods such as breaching and filling on digital elevation models with different resolutions (2, 4, 8, and 16 m) and to evaluate which preprocessing methods most accurately route water across road impoundments at actual culvert locations. A unique dataset with over 30,000 field‐mapped road culverts was used to assess the accuracy of stream networks derived from digital elevation models using different preprocessing methods. Our results showed that the accuracy of stream networks increases with increasing resolution. Breaching created the most accurate stream networks on all resolutions, whereas filling was the least accurate. Burning streams from the topographic map across roads from the topographic map increased the accuracy for all methods and resolutions. In addition, the impact in terms of change in area and absolute volume between original and preprocessed digital elevation models was smaller for breaching than for filling. With the appropriate methods, it is possible to extract accurate stream networks from high‐resolution digital elevation models with extensive road networks, thus providing forest managers with stream networks that can be used when planning operations in wet areas or areas near streams to prevent rutting, sediment transport, and mercury export.     

A combined algorithm for automated drainage network extraction from digital elevation models

Drainage networks are the basis for segmentation of watersheds, an essential component in hydrological modelling, biogeochemical applications, and resource management plans. With the rapidly increasing availability of topographic information as digital elevation models (DEMs), there have been many studies on DEM‐based drainage network extraction algorithms. Most of traditional drainage network extraction methods require preprocessing of the DEM in order to remove “spurious” sink, which can cause unrealistic results due to removal of real sinks as well. The least cost path (LCP) algorithm can deal with flow routing over sinks without altering data. However, the existing LCP implementations can only simulate either single flow direction or multiple flow direction over terrain surfaces. Nevertheless, terrain surfaces in the real world are usually very complicated including both convergent and divergent flow patterns. The triangular form‐based multiple flow (TFM) algorithm, one of the traditional drainage network extraction methods, can estimate both single flow and multiple flow patterns. Thus, in this paper, it is proposed to combine the advantages of the LCP algorithm and the TFM algorithm in order to improve the accuracy of drainage network extraction from the DEM. The proposed algorithm is evaluated by implementing a data‐independent assessment method based on four mathematical surfaces and validated against “true” stream networks from aerial photograph, respectively. The results show that when compared with other commonly used algorithms, the new algorithm provides better flow estimation and is able to estimate both convergent and divergent flow patterns well regarding the mathematical surfaces and the real‐world DEM.     

Hydrological response characteristics of Mediterranean catchments at different time scales: a meta-analysis

ABSTRACT

This work examines 140 hydrological studies conducted in the Mediterranean region. It identifies key characteristics of the hydrological responses of Mediterranean catchments at various time scales and compares different methods and modelling approaches used for individual-catchment studies. The study area is divided into the northwestern (NWM), eastern (EM) and southern (SM) Mediterranean. The analysis indicates regional discrepancies in which the NWM shows the most extreme rainfall regime. A tendency for reduced water resources driven by both anthropogenic and climatic pressures and a more extreme rainfall regime are also noticeable. Catchments show very heterogeneous responses over time and space, resulting in limitations in hydrological modelling and large uncertainties in predictions. However, few models have been developed to address these issues. Additional studies are necessary to improve the knowledge of Mediterranean hydrological features and to account for regional specificities.

Editor D. Koutsoyiannis Associate editor A. Efstratiadis     

Testing of the SIBERIA landscape evolution model using the Tin Camp Creek,Northern Territory,Australia, field catchment

The SIBERIA landscape evolution model was used to simulate the geomorphic development of the Tin Camp Creek natural catchment over geological time. Measured hydrology, erosion and geomorphic data were used to calibrate the SIBERIA model, which was then used to make independent predictions of the landform geomorphology of the study site. The catchment, located in the Northern Territory, Australia is relatively untouched by Europeans so the hydrological and erosion processes that shaped the area can be assumed to be the same today as they have been in the past, subject to the caveats regarding long‐term climate fluctuation. A qualitative, or visual comparison between the natural and simulated catchments indicates that SIBERIA can match hillslope length and hillslope profile of the natural catchments. A comparison of geomorphic and hydrological statistics such as the hypsometric curve, width function, cumulative area distribution and area–slope relationship indicates that SIBERIA can model the geomorphology of the selected Tin Camp Creek catchments. Copyright 2002 © Environmental Research Institute of the Supervising Scientist, Commonwealth of Australia.     

The concept of effective length in hillslopes: assessing the influence of climate and topography on the contributing areas of catchments

This paper introduces a concept of ‘effective length’ in hillslopes to define the effective area influencing the runoff‐producing saturated zones of a hillslope or catchment. This effective area of a catchment usually is less than that given by its physical boundaries, particularly in regions where the total potential evaporation exceeds total rainfall on an annual basis. In this paper, expressions for effective lengths in hillslopes with different scale, shape and soil properties are derived for given climatic conditions. The influence of these variables on effective length is investigated. It is shown that, for a given rainfall frequency and soil parameters, the effective length changes with the planform geometry and profile shape of a hillslope; it is also a function of the ratio of available travel opportunity time to the hillslope's scale response time. The application of the concept to three natural catchments, subdivided into a number of simple hillslopes, is described. It is shown that, for these three test catchments and over 24 years of record, rarely would the entire catchment areas contribute to flow at the respective outlets. The implications of the concept of effective length for several land‐use practices, such as clearing for forest for greater water yield, and planting trees for salinity control, are discussed. Copyright © 2002 John Wiley & Sons, Ltd.     

A comparision of interpolation methods for producing digital elevation models at the field scale

Digital elevation models have been used in many applications since they came into use in the late 1950s. It is an essential tool for applications that are concerned with the Earth's surface such as hydrology, geology, cartography, geomorphology, engineering applications, landscape architecture and so on. However, there are some differences in assessing the accuracy of digital elevation models for specific applications. Different applications require different levels of accuracy from digital elevation models. In this study, the magnitudes and spatial patterning of elevation errors were therefore examined, using different interpolation methods. Measurements were performed with theodolite and levelling. Previous research has demonstrated the effects of interpolation methods and the nature of errors in digital elevation models obtained with indirect survey methods for small‐scale areas. The purpose of this study was therefore to investigate the size and spatial patterning of errors in digital elevation models obtained with direct survey methods for large‐scale areas, comparing Inverse Distance Weighting, Radial Basis Functions and Kriging interpolation methods to generate digital elevation models. The study is important because it shows how the accuracy of the digital elevation model is related to data density and the interpolation algorithm used. Cross validation, split‐sample and jack‐knifing validation methods were used to evaluate the errors. Global and local spatial auto‐correlation indices were then used to examine the error clustering. Finally, slope and curvature parameters of the area were modelled depending on the error residuals using ordinary least regression analyses. In this case, the best results were obtained using the thin plate spline algorithm. Copyright © 2009 John Wiley & Sons, Ltd.     

A tracer‐based assessment of hydrological pathways at different spatial scales in a mesoscale Scottish catchment

Geochemically based hydrograph separation techniques were used in a preliminary assessment to infer how runoff processes change with landscape characteristics and spatial scale (1–233 km²) within a mesoscale catchment in upland Scotland. A two‐component end‐member mixing analysis (EMMA) used Gran alkalinity as an assumed conservative tracer. Analysis indicated that, at all scales investigated, acidic overland flow and shallow subsurface storm flows from the peaty soils covering the catchment headwaters dominated storm runoff generation. The estimated groundwater contribution to annual runoff varied from 30% in the smallest (ca 1 km²) peat‐dominated headwater catchment with limited groundwater storage, to >60% in larger catchments (>30 km²) with greater coverage of more freely draining soils and more extensive aquifers in alluvium and other drift. This simple approach offers a useful, integrated conceptualization of the hydrological functioning in a mesoscale catchment, which can be tested and further refined by focused modelling and process‐based research. However, even as it stands, the simple conceptualization of system behaviour will have significant utility as a tool for communicating hydrological issues in a range of planning and management decisions. Copyright © 2002 John Wiley & Sons, Ltd.     

A comparison of drainage networks derived from digital elevation models at two scales

Automated generation of drainage networks has become increasingly popular with powerful analytical functions in geographic information systems (GIS) and with the increased availability of digital elevation models (DEMs). This paper compares drainage networks derived from DEMs at two scales, 1:250 000 (250K) and 1:24 000 (24K), using various drainage parameters common in hydrology and geomorphology. The comparison of parameters derived from the 250K DEMs with those from the 24K DEMs in 20 basins ranging from 150 to 1000 km² in West Virginia shows that the goodness-of-fit between parameter estimates based on the DEMs varies. Results clearly show that superior estimations are produced from the 24K DEMs. Better estimates can be obtained from the 250K DEMs for stream length and frequency parameters than for gradient parameters. However, the estimation of the mean gradient parameters based on the 250K DEMs seems to improve with increasing terrain complexity. Finally, basin size does not strongly affect the accuracy of parameter estimates based on the 250K DEMs.     

Advancing understanding of runoff and sediment transfers in agricultural catchments through simultaneous observations across scales

Our understanding of the effect of scale on runoff and sediment transfers within catchments is currently limited by a lack of available data. A multi‐scale dataset of 17 rainfall events collected simultaneously at four spatial scales within a small agricultural catchment in 2005–2006 is presented. Analysis using exploratory techniques and a two‐step, zero‐inflated lognormal mixed‐effects regression model, has demonstrated that event responses, and event response characteristics representing runoff and sediment peaks and area‐normalized yields, are scale dependent, and hence cannot be transferred directly between scales. Runoff and sediment yields increase as scale increases, and it is proposed that this effect, which differs from that observed in the few other studies of scale effects undertaken, is due to increasing connectivity within the catchment, and the dominance of preferential flow pathways including through macropores and field drains. The processes contributing to scale dependence in the data, and the possibility that certain processes dominate at particular scales, are discussed. The data presented here help to improve our spatial understanding of runoff and sediment transport in small agricultural catchments, and provide examples of the type of spatial dataset and the type of analysis that are essential if we are to develop models which are able to predict runoff and soil erosion accurately, and allow us to manage runoff and sediment transport effectively across scales. Copyright © 2011 John Wiley & Sons, Ltd.     

Gully erosion in Mangatu Forest,New Zealand,estimated from digital elevation models

The methodology and errors involved in determining the amount of sediment produced during two (19·5 and 33·2 year) periods by 11 (c. 0·01 − >0·20 km²) gullies within a 4 km² area in the headwaters of the Waipaoa River basin, New Zealand, using sequential digital elevation models are described. Sediment production from all gullies within the study area was 0·99 ± 0·03 × 10⁶ t a⁻¹ (2480 ± 80 t ha⁻¹ a⁻¹) during the period from 1939 to 1958. It declined to 0·62 ± 0·02 × 10⁶ t a⁻¹ (1550 ± 50 t ha⁻¹ a⁻¹) during the period from 1958 to 1992, when many of the smaller gullies were stabilized by a programme of afforestation, which commenced in 1960. Both figures are very high by global standards. The two largest (the Tarndale and Mangatu) gully complexes together generated 73 and 95 per cent of the sediment in the specified time periods, but the latter amount is equivalent to only c. 5 per cent of the total annual sediment load of the Waipaoa River. © 1998 John Wiley & Sons, Ltd.     

A qualitative and quantitative evaluation of experimental model catchment evolution

Due to the geological time scales required for observation of catchment evolution, surrogates or analogues of field data are necessary to understand long‐term processes. To investigate long‐term catchment behaviour, two experimental model catchments that developed without rigid boundaries under controlled conditions are examined and a qualitative and quantitative analysis of their evolution is presented. Qualitatively, the experimental catchments have the visual appearance of field scale data. Observation demonstrates that changes in catchment shape and network form are conservative. Quantitative analysis suggests that the catchments reach an equilibrium form while a reduction in the channel network occurs. While the catchments are laboratory scale models, the results provide insights into field scale behaviour. Copyright © 2003 John Wiley & Sons, Ltd.