Abstract. Landscape evolution models (LEM) allow studying the earth surface response to a changing climatic and tectonic forcing. While much effort has been devoted to the development of LEMs that simulate a wide range of processes, the numerical accuracy of these models has received much less attention. Most LEMs use first order accurate numerical methods that suffer from substantial numerical diffusion. Numerical diffusion particularly affects the solution of the advection equation and thus the simulation of retreating landforms such as cliffs and river knickpoints with potential unquantified consequences for the integrated response of the simulated landscape. Here we present TTLEM, a spatially explicit, raster based LEM for the study of fluvially eroding landscapes in TopoToolbox 2. TTLEM prevents numerical diffusion by implementing a higher order flux limiting total volume method that is total variation diminishing (TVD-TVM) and solves the partial differential equations of river incision and tectonic displacement. We show that the choice of the TVD-TVM to simulate river incision significantly influences the evolution of simulated landscapes and the spatial and temporal variability of catchment wide erosion rates. Furthermore, a 2D TVD-TVM accurately simulates the evolution of landscapes affected by lateral tectonic displacement, a process whose simulation is hitherto largely limited to LEMs with flexible spatial discretization. By providing accurate numerical schemes on rectangular grids, TTLEM is a widely accessible LEM that is compatible with GIS analysis functions from the TopoToolbox interface.