An effective depression filling algorithm for DEM-based 2-dimensional surface flow modelling
Abstract. The surface runoff process in fluvial/pluvial flood modelling is often simulated employing a two-dimensional (2-D) diffusive wave approximation to described by grid based digital elevation models (DEMs). However, a serious problem of this approach may arise when using a 2-D surface flow model which exchanges flows through adjacent cells, or conventional rink removal algorithms which also allow flow to be exchanged along diagonal directions, due to the existence of artificial depression in DEMs. This study firstly analyses the two types of depressions in DEMs and reviews the current depression filling algorithms with a medium sized basin in South-East England, the Upper Medway Catchment (220 km2) used to demonstrate the depression issue in 2-D surface runoff simulation by MIKE SHE with different DEM resolutions (50 m, 100 m and 200 m). An alternative depression-filling algorithm for 2-D overland flow modelling is developed and evaluated by comparing the simulated flows at the outlet of the catchment. This result suggests that the depression estimates at different grid resolution of DEM highly influences overland flow estimation and the new depression filling algorithm is shown to be effective in tackling this issue when comparing simulations in sink-dominated and sink-free digital elevation models, especially for depressions in relatively flat areas on digital land surface models.