Modelling catchment‐scale shallow landslide occurrence and sediment yield as a function of rainfall return period

A model‐based method is proposed for improving upon existing threshold relationships which define the rainfall conditions for triggering shallow landslides but do not allow the magnitude of landsliding (i.e. the number of landslides) to be determined. The SHETRAN catchment‐scale shallow landslide model is used to quantify the magnitude of landsliding as a function of rainfall return period, for focus sites of 180 and 45 km² in the Italian Southern Alps and the central Spanish Pyrenees. Rainfall events with intensities of different return period are generated for a range of durations (1‐day to 5‐day) and applied to the model to give the number of landslides triggered and the resulting sediment yield for each event. For a given event duration, simulated numbers of landslides become progressively less sensitive to return period as return period increases. Similarly, for an event of given return period, landslide magnitude becomes less sensitive to event duration as duration increases. The temporal distribution of rainfall within an event is shown to have a significant impact on the number of landslides and the timing of their occurrence. The contribution of shallow landsliding to catchment sediment yield is similarly quantified as a function of the rainfall characteristics. Rainfall intensity–duration curves are presented which define different levels of landsliding magnitude and which advance our predictive capability beyond, but are generally consistent with, published threshold curves. The magnitude curves are relevant to the development of guidelines for landslide hazard assessment and forecasting. Copyright © 2011 John Wiley & Sons, Ltd.     

Effects of spatial aggregation of soil spatial information on watershed hydrological modelling

Many researchers have examined the impact of detailed soil spatial information on hydrological modelling due to the fact that such information serves as important input to hydrological modelling, yet is difficult and expensive to obtain. Most research has focused on the effects at single scales; however, the effects in the context of spatial aggregation across different scales are largely missing. This paper examines such effects by comparing the simulated runoffs across scales from watershed models based on two different levels of soil spatial information: the 10‐m‐resolution soil data derived from the Soil‐Land Inference Model (SoLIM) and the 1:24000 scale Soil Survey Geographic (SSURGO) database in the United States. The study was conducted at three different spatial scales: two at different watershed size levels (referred to as full watershed and sub‐basin, respectively) and one at the model minimum simulation unit level. A fully distributed hydrologic model (WetSpa) and a semi‐distributed model (SWAT) were used to assess the effects. The results show that at the minimum simulation unit level the differences in simulated runoff are large, but the differences gradually decrease as the spatial scale of the simulation units increases. For sub‐basins larger than 10 km² in the study area, stream flows simulated by spatially detailed SoLIM soil data do not significantly vary from those by SSURGO. The effects of spatial scale are shown to correlate with aggregation effect of the watershed routing process. The unique findings of this paper provide an important and unified perspective on the different views reported in the literature concerning how spatial detail of soil data affects watershed modelling. Different views result from different scales at which those studies were conducted. In addition, the findings offer a potentially useful basis for selecting details of soil spatial information appropriate for watershed modelling at a given scale. Copyright © 2011 John Wiley & Sons, Ltd.     

Flow network derivation from a high resolution DEM in a low relief,agrarian landscape

Digital flow networks derived from digital elevation models (DEMs) sensitively react to errors due to measurement, data processing and data representation. Since high‐resolution DEMs are increasingly used in geomorphological and hydrological research, automated and semi‐automated procedures to reduce the impact of such errors on flow networks are required. One such technique is stream‐carving, a hydrological conditioning technique to ensure drainage connectivity in DEMs towards the DEM edges. Here we test and modify a state‐of‐the‐art carving algorithm for flow network derivation in a low‐relief, agricultural landscape characterized by a large number of spurious, topographic depressions. Our results show that the investigated algorithm reconstructs a benchmark network insufficiently in terms of carving energy, distance and a topological network measure. The modification to the algorithm that performed best, combines the least‐cost auxiliary topography (LCAT) carving with a constrained breaching algorithm that explicitly takes automatically identified channel locations into account. We applied our methods to a low relief landscape, but the results can be transferred to flow network derivation of DEMs in moderate to mountainous relief in situations where the valley bottom is broad and flat and precise derivations of the flow networks are needed. Copyright © 2013 John Wiley & Sons, Ltd.     

Interactions between subgrid‐scale resolution,feature representation and grid‐scale resolution in flood inundation modelling

Numerical modelling of flood inundation over large and complex floodplains often requires mesh resolutions coarser than the structural features (e.g. buildings) that are known to influence the inundation process. Recent research has shown that this mismatch is not well represented by conventional roughness treatments, but that finer‐scale features can be represented through porosity‐based subgrid‐scale treatments. This paper develops this work by testing the interactions between feature representation, subgrid‐scale resolution and mesh resolution. It uses as the basis for this testing a 2D diffusion‐based flood inundation model which is applied to a 2004 flood event in a topologically complex upland floodplain in northern England. This study formulated simulations with different grid mesh resolution and subgrid mesh ratio. The sensitivity of the model to mesh resolution and roughness specification was investigated. Model validation and verification suggest that the subgrid treatment with higher subgrid mesh ratio can give much improved predictions of flood propagation, in particular, in terms of the predicted water depth. This study also highlighted the limitation of using at‐a‐point in time inundation extent for validation of flood models of this type. Copyright © 2010 John Wiley & Sons, Ltd.     

Geomorphological and historical data in assessing landslide hazard

Traditionally, earth scientists assess landslide occurrence on the basis of geomorphological investigations carried out through aerial photograph interpretation and ?eldwork. Conversely, local administrators primarily evaluate the impact of natural catastrophes, such as landsliding, on the basis of historical records and data. Owing to the substantial difference in the structure and spatial density of these two types of information, it is dif?cult to compare them directly and few investigators have attempted this. We compared landslide information derived from geomorphological mapping and historical data in a pilot area (the Staffora river basin, northern Italy). To do this we generated two multivariate statistical models where the dependent variable was either the mapped landslide deposits (geomorphological model), or the historical sites affected by landslide‐induced damage (historical model). By quantitatively comparing these two model maps, we demonstrate that the geomorphological model performs better in terms of percentage of terrain units correctly predicted as stable or unstable. The historical model underestimates landslide hazard mainly where human structures are lacking. However, it highlights slopes where landslide movements take place with a high frequency at the temporal scale of human life. Hence, the joint use of these two models may facilitate the knowledge of the overall instability conditions of a given region and the identi?cation of the landslides that are most frequently reactivated. Copyright © 2003 John Wiley & Sons, Ltd.     

Hydrological modelling of a small watershed using generated rainfall in the soil and water assessment tool model

An adequately tested soil and water assessment tool (SWAT) model was applied to the runoff and sediment yield of a small agricultural watershed in eastern India using generated rainfall. The capability of the model for generating rainfall was evaluated for a period of 18 years (1981–1998). The watershed and subwatershed boundaries, drainage networks, slope, soil series and texture maps were generated using a geographical information system (GIS). A supervised classification method was used for land‐use/cover classification from satellite imageries. Model simulated monthly rainfall for the period of 18 years was compared with observations. Simulated monthly rainfall, runoff and sediment yield values for the monsoon season of 8 years (1991–1998) were also compared with their observed values. In general monthly average rainfall predicted by the model was in close agreement with the observed monthly average values. Also, simulated monthly average values of surface runoff and sediment yield using generated rainfall compared well with observed values during the monsoon season of the years 1991–1998. Results of this study revealed that the SWAT model can generate monthly average rainfall satisfactorily and thereby can produce monthly average values of surface runoff and sediment yield close to the observed values. Therefore, it can be concluded that the SWAT model could be used for developing a multiple year management plan for the critical erosion prone areas of a small watershed. Copyright © 2004 John Wiley & Sons, Ltd.     

The impact of depression removal on catchment geomorphology,soil erosion and landscape evolution

The current generation of landscape evolution models use a digital elevation model for landscape representation. These programs also contain a hydrological model that defines overland flow with the drainage network routed to an outlet. One of the issues with landscape evolution modelling is the hydrological correctness of the digital elevation model used for the simulations. Despite the wide use and increased quality of digital elevation models, data pits and depressions in the elevation data are a common feature and their removal will remain a necessary step for many data sets. This study examines whether a digital elevation model can be hydrologically correct (i.e. all depressions removed so that all water can run downslope) before use in a landscape evolution model and what effect depression removal has on long‐term geomorphology and hydrology. The impact on sediment transport rates is also examined. The study was conducted using a field catchment and a proposed landform for a post‐mining landscape. The results show that there is little difference in catchment geomorphology and hydrology for the non‐depression removed and depression removed data sets. The non‐depression removed and depression removed digital elevation models were also evaluated as input to a landscape evolution model for a 50 000 year simulation period. The results show that after 1000 years there is little difference between the data sets, although sediment transport rates did vary considerably early on in the simulation. Copyright © 2007 John Wiley & Sons, Ltd.     

Effect of digital elevation model spatial resolution on depression storage

Surface water storage—including wetlands and other small waterbodies—has largely been disregarded in traditional hydrological models. In this paper, the grid resampling method is adopted to study the influence of the digital elevation model (DEM) grid resolution on depression storage (DS) considering different rainfall return periods. It is observed that the DEM grid size highly affects DS, and the higher the grid resolution is, the larger the DS value. However, when the grid resolution reaches a certain value, the maximum DS value decreases. This suggests that a critical grid resolution value exists at which the water storage capacity of depressions is maximized, namely, 20 m in this work (except for the overall area simulation under infiltration). This phenomenon is further verified in two test cases with and without the infiltration process, that is, calculations of the local area and without infiltration area, respectively. This research may facilitate the accurate computation of the DS process, which is greatly affected by the grid resolution, thereby improving the reliability of hydrological models.     

How to evaluate the quality of coarse‐resolution DEM‐derived drainage networks

The quality of digital elevation model (DEM)‐derived river drainage networks (RDNs) is influenced by DEM quality, basin physical characteristics, scale, and algorithms used; these factors should not be neglected. However, few research studies analyse the different evaluation approaches used in the literature with respect to adequacy, meaning of the results, advantages, and limitations. Focusing on coarse‐resolution networks, this paper reviews the use of these techniques and offers new insights on these issues. Additionally, we propose adaptations for selected metrics and discuss distinct interpretations for the evaluation of RDNs derived at different spatial resolutions (1, 5, 10, 20, and 30 km) considering the Uruguay River basin (206,000 km²) as a case study. The results demonstrate that lumped basin/river characteristics and basin delineation analysis should not be used as evaluation criteria for RDN quality; however, some of these metrics offer useful complementary information. Percentage of the DEM‐derived RDN within a uniform buffer placed around a river network considered as reference and mean separation distance between these two networks are more suitable metrics, but the former is insensitive to serious errors. The change in reference from a fine‐scale network to a coarse‐resolution manual tracing network significantly augments the discrepancy of these largest errors when the mean distance metric was applied, and visual comparison analysis is necessary to interpret the results for other metrics. We recommend the use of the mean distance metric in combination with a detailed visual assessment, the importance of which increases as the resolution coarsens. In both cases, the impact of network quality can be further refined by quantifying the basin shape and river length errors.     

Modelling the effects of spatial and temporal resolution of rainfall and basin model on extreme river discharge

Abstract

Important characteristics of an appropriate river basin model, intended to study the effect of climate change on basin response, are the spatial and temporal resolution of the model and the rainfall input. The effects of input and model resolution on extreme discharge of a large river basin are assessed to give some indication on appropriate resolutions. A simple stochastic rainfall model and a river basin model with uniform parameters and multiple rainfall input have been developed and applied to the River Meuse basin in northwestern Europe. The results show that the effect of model resolution on extreme river discharge is much greater than that of input resolution. The highest model resolution seems to be quite accurate in determining extreme discharge. Although the results should be interpreted with caution, they may give some indication of appropriate input and model resolutions for the determination of extreme discharge of a large river basin.     

The influence of elevation error on the morphometrics of channel networks extracted from DEMs and the implications for hydrological modelling

Stream network morphometrics have been used frequently in environmental applications and are embedded in several hydrological models. This is because channel network geometry partly controls the runoff response of a basin. Network indices are often measured from channels that are mapped from digital elevation models (DEMs) using automated procedures. Simulations were used in this paper to study the influence of elevation error on the reliability of estimates of several common morphometrics, including stream order, the bifurcation, length, area and slope ratios, stream magnitude, network diameter, the flood magnitude and timing parameters of the geomorphological instantaneous unit hydrograph (GIUH) and the network width function. DEMs of three UK basins, ranging from high to low relief, were used for the analyses. The findings showed that moderate elevation error (RMSE of 1·8 m) can result in significant uncertainty in DEM‐mapped network morphometrics and that this uncertainty can be expressed in complex ways. For example, estimates of the bifurcation, length and area ratios and the flood magnitude and timing parameters of the GIUH each displayed multimodal frequency distributions, i.e. two or more estimated values were highly likely. Furthermore, these preferential estimates were wide ranging relative to the ranges typically observed for these indices. The wide‐ranging estimates of the two GIUH parameters represented significant uncertainty in the shape of the unit hydrograph. Stream magnitude, network diameter and the network width function were found to be highly sensitive to elevation error because of the difficulty in mapping low‐magnitude links. Uncertainties in the width function were found to increase with distance from outlet, implying that hydrological models that use network width contain greater uncertainty in the shape of the falling limb of the hydrograph. In light of these findings, care should be exercised when interpreting the results of analyses based on DEM‐mapped stream networks. Copyright © 2007 John Wiley & Sons, Ltd.     

Influence of DEM resolution on modelling hydrological connectivity in a complex agricultural catchment with woody crops

Digital elevation model (DEM) resolution influences hydrological simulation. However, its influence when modelling hydrological connectivity (HC) in woody crops remains to be seen. We assessed surface topography, microtopography and HC in an agricultural sub‐catchment (27.4 ha) using six photogrammetry‐derived DEMs with 0.03, 0.05, 0.1, 0.2, 0.5 and 1 m cell sizes. Land uses included vineyards, olive groves, cereal fields, and forestry patches. We ran an updated version of Borselli's HC index (IC) using the D‐Infinity approach. We assessed HC in woody crops at high spatial resolution for the first time. After analysing the differences in the contributing area, the flow width, the soil roughness, the convergence index and the IC (normalised and non‐normalised) at different scales (hillslope, land uses and compartments, ephemeral gullies, depositional areas and the sub‐catchment outlet) and accounting for the field vertical components, we propose an optimum DEM resolution (0.2 m) to improve modelling of structural HC in woody crops. The modelled hydrological features at this threshold resolution matched well with the geomorphic features associated with the short‐ and medium‐term patterns of soil redistribution. Higher DEM resolutions, especially at 0.03 and 0.05 m, introduced bias in the input data and the IC computations. Finally, we observed good agreement between the outputs at the lowest resolution, 1 × 1 m, and the long‐term soil redistribution patterns (functional connectivity). Copyright © 2017 John Wiley & Sons, Ltd.     

Grid‐resolution effects on a model for integrating urban and rural areas

The impact of grid‐cell size on calibrated parameters and on the performance of a variable source area model intended for urbanizing catchments is studied in this research. The model uses TOPMODEL concepts that were modified to consider urban areas in both the topographic index and the mechanism of surface runoff generation. The revised model known as TOPURBAN, was applied to a small catchment of roughly 8 km² in southern Ontario. Ten different grid‐cell sizes ranging from 10 m to 100 m were selected to study scale effects in this catchment with mild to moderate relief, on three separate time periods. The model performed adequately with calibration efficiencies for all three time periods in the range of 0\65 to 0\85. The verification efficiencies were not as high and ranged from 0\4 to 0\6. Larger cell sizes produced higher averages of topographic index, and this resulted in larger calibrated transmissivities. The most important parameter in determining the quantity of urban runoff was slightly affected by grid resolution. During the calibration process, this parameter was also found to interact with important parameters that dealt primarily with rural runoff generation. As cell size increased, contributions from urban areas increased and overland flow contributions from rural areas decreased. Results showed that in this integrated model of urban and rural areas, predicted processes based on calibrated parameters were dependent on grid resolution. Calibration of internal state variables is recommended to draw conclusions on the influences between urban and rural areas on the overall flow. Copyright © 2000 John Wiley & Sons, Ltd.