Using high‐resolution isotope data and alternative calibration strategies for a tracer‐aided runoff model in a nested catchment

Testing hydrological models over different spatio‐temporal scales is important for both evaluating diagnostics and aiding process understanding. High‐frequency (6‐hr) stable isotope sampling of rainfall and runoff was undertaken during 3‐week periods in summer and winter within 12 months of daily sampling in a 3.2‐km² catchment in the Scottish Highlands. This was used to calibrate and test a tracer‐aided model to assess the (a) information content of high‐resolution data, (b) effect of different calibration strategies on simulations and inferred processes, and (c) model transferability to <1‐km² subcatchment. The 6‐hourly data were successfully incorporated without loss of model performance, improving the temporal resolution of the modelling, and making it more relevant to the time dynamics of the isotope and hydrometric response. However, this added little new information due to old‐water dominance and riparian mixing in this peatland catchment. Time variant results, from differential split sample testing, highlighted the importance of calibrating to a wide range of hydrological conditions. This also provided insights into the nonstationarity of catchment mixing processes, in relation to storage and water ages, which varied markedly depending on the calibration period. Application to the nested subcatchment produced equivalent parameterization and performance, highlighting similarity in dominant processes. The study highlighted the utility of high‐resolution data in combination with tracer‐aided models, applied at multiple spatial scales, as learning tools to enhance process understanding and evaluation of model behaviour across nonstationary conditions. This helps reveal more fully the catchment response in terms of the different mechanistic controls on both wave celerites and particle velocities.     

On the effects of triangulated terrain resolution on distributed hydrologic model response

Distributed hydrologic models based on triangulated irregular networks (TIN) provide a means for computational efficiency in small to large‐scale watershed modelling through an adaptive, multiple resolution representation of complex basin topography. Despite previous research with TIN‐based hydrology models, the effect of triangulated terrain resolution on basin hydrologic response has received surprisingly little attention. Evaluating the impact of adaptive gridding on hydrologic response is important for determining the level of detail required in a terrain model. In this study, we address the spatial sensitivity of the TIN‐based Real‐time Integrated Basin Simulator (tRIBS) in order to assess the variability in the basin‐averaged and distributed hydrologic response (water balance, runoff mechanisms, surface saturation, groundwater dynamics) with respect to changes in topographic resolution. Prior to hydrologic simulations, we describe the generation of TIN models that effectively capture topographic and hydrographic variability from grid digital elevation models. In addition, we discuss the sampling methods and performance metrics utilized in the spatial aggregation of triangulated terrain models. For a 64 km² catchment in northeastern Oklahoma, we conduct a multiple resolution validation experiment by utilizing the tRIBS model over a wide range of spatial aggregation levels. Hydrologic performance is assessed as a function of the terrain resolution, with the variability in basin response attributed to variations in the coupled surface–subsurface dynamics. In particular, resolving the near‐stream, variable source area is found to be a key determinant of model behaviour as it controls the dynamic saturation pattern and its effect on rainfall partitioning. A relationship between the hydrologic sensitivity to resolution and the spatial aggregation of terrain attributes is presented as an effective means for selecting the model resolution. Finally, the study highlights the important effects of terrain resolution on distributed hydrologic model response and provides insight into the multiple resolution calibration and validation of TIN‐based hydrology models. Copyright © 2005 John Wiley & Sons, Ltd.     

How to improve the representation of hydrological processes in SWAT for a lowland catchment – temporal analysis of parameter sensitivity and model performance

Model diagnostic analyses help to improve the understanding of hydrological processes and their representation in hydrological models. A detailed temporal analysis detects periods of poor model performance and model components with potential for model improvements, which cannot be found by analysing the whole discharge time series. In this study, we aim to improve the understanding of hydrological processes by investigating the temporal dynamics of parameter sensitivity and of model performance for the Soil and Water Assessment Tool model applied to the Treene lowland catchment in Northern Germany. The temporal analysis shows that the parameter sensitivity varies temporally with high sensitivity for three groundwater parameters (groundwater time delay, baseflow recession constant and aquifer fraction coefficient) and one evaporation parameter (soil evaporation compensation factor). Whereas the soil evaporation compensation factor dominates in baseflow and resaturation periods, groundwater time delay, baseflow recession constant and aquifer fraction coefficient are dominant in the peak and recession phases. The temporal analysis of model performance identifies three clusters with different model performances, which can be related to different phases of the hydrograph. The lowest performance, when comparing six performance measures, is detected for the baseflow cluster. A spatially distributed analysis for six hydrological stations within the Treene catchment shows similar results for all stations. The linkage of periods with poor model performance to the dominant model components in these phases and with the related hydrological processes shows that the groundwater module has the highest potential for improvement. This temporal diagnostic analysis enhances the understanding of the Soil and Water Assessment Tool model and of the dominant hydrological processes in the lowland catchment. Copyright © 2013 John Wiley & Sons, Ltd.     

Comparison of stormflow responses of surface‐mined and forested watersheds in the Appalachian Mountains,USA

The results of a hydrological analysis that was conducted as part of a larger, multifaceted, collaborative effort to quantify ecosystem functions in watersheds subjected to land‐use and land‐cover change are presented. The primary goal of the study was to determine whether a small watershed in the Appalachian region (USA) that was recently subjected to surface mining and reclamation practices produces stormflow responses to rain events that are different from those produced by a nearby reference watershed covered by young, second‐growth forest. Water balances indicated that runoff yields did not vary significantly between the two watersheds on an annual basis. Statistically significant differences (p?0·05) in runoff responses were observed on an event basis, however, with the mined/reclaimed watershed producing, on average (a) higher storm runoff coefficients (2·5×), (b) greater total storm runoff (3×), and (c) higher peak hourly runoff rates (2×) when compared with the reference watershed. Results of a unit hydrograph analysis also showed, unexpectedly, that the modelled unit responses of the two watersheds to effective rainfall pulses were similar, despite the noted differences in land cover. Differences in stormflow responses were thus largely explained by dramatic reductions in cumulative rates of rainfall abstraction (measured using infiltrometers) attributable to soil compaction during land reclamation. Additional field hydrological measurements on other mined watersheds will be needed to generalize our results, as well as to understand and predict the cumulative hydrological impacts of widespread surface mining in larger watersheds and river basins. Copyright © 2006 John Wiley & Sons, Ltd.     

Performance of a conceptual and physically based model in simulating the response of a semi‐urbanized watershed in San Antonio,Texas

The need for accurate hydrologic analysis and rainfall–runoff modelling tools has been rapidly increasing because of the growing complexity of operational hydrologic and hydraulic problems associated with population growth, rapid urbanization and expansion of agricultural activities. Given the recent advances in remote sensing of physiographic features and the availability of near real‐time precipitation products, rainfall–runoff models are expected to predict runoff more accurately. In this study, we compare the performance and implementation requirements of two rainfall–runoff models for a semi‐urbanized watershed. One is a semi‐distributed conceptual model, the Hydrologic Engineering Center‐Hydrologic Modelling System (HEC‐HMS). The other is a physically based, distributed‐parameter hydrologic model, the Gridded Surface Subsurface Hydrologic Analysis (GSSHA). Four flood events that took place on the Leon Creek watershed, a sub‐watershed of the San Antonio River basin in Texas, were used in this study. The two models were driven by the Multisensor Precipitation Estimator radar products. One event (in 2007) was used for HEC‐HMS and GSSHA calibrations. Two events (in 2004 and 2007) were used for further calibration of HEC‐HMS. Three events (in 2002, 2004 and 2010) were used for model validation. In general, the physically based, distributed‐parameter model performed better than the conceptual model and required less calibration. The two models were prepared with the same minimum required input data, and the effort required to build the two models did not differ substantially. Copyright © 2012 John Wiley & Sons, Ltd.     

Field investigation and modelling of coupled stream discharge and shallow water‐table dynamics in a small Mediterranean catchment (Sardinia)

In the semi‐arid Mediterranean environment, the rainfall–runoff relationships are complex because of the markedly irregular patterns in rainfall, the seasonal mismatch between evaporation and rainfall, and the spatial heterogeneity in landscape properties. Watersheds often display considerable non‐linear threshold behavior, which still make runoff generation an open research question. Our objectives in this context were: to identify the primary processes of runoff generation in a small natural catchment; to test whether a physically based model, which takes into consideration only the primary processes, is able to predict spatially distributed water‐table and stream discharge dynamics; and to use the hydrological model to increase our understanding of runoff generation mechanisms. The observed seasonal dynamics of soil moisture, water‐table depth, and stream discharge indicated that Hortonian overland‐flow was negligible and the main mechanism of runoff generation was saturated subsurface‐flow. This gives rise to base‐flow, controls the formation of the saturated areas, and contributes to storm‐flow together with saturation overland‐flow. The distributed model, with a 1D scheme for the kinematic surface‐flow, a 2D sub‐horizontal scheme for the saturated subsurface‐flow, and ignoring the unsaturated flow, performed efficiently in years when runoff volume was high and medium, although there was a smoothing effect on the observed water‐table. In dry years, small errors greatly reduced the efficiency of the model. The hydrological model has allowed to relate the runoff generation mechanisms with the land‐use. The forested hillslopes, where the calibrated soil conductivity was high, were never saturated, except at the foot of the slopes, where exfiltration of saturated subsurface‐flow contributed to storm‐flow. Saturation overland‐flow was only found near the streams, except when there were storm‐flow peaks, when it also occurred on hillslopes used for pasture, where soil conductivity was low. The bedrock–soil percolation, simulated by a threshold mechanism, further increased the non‐linearity of the rainfall–runoff processes. Copyright © 2013 John Wiley & Sons, Ltd.     

Effects of area under‐estimations of sloped mountain terrain on simulated hydrological behaviour: a case study using the ACRU model

Sloped areas calculated from a GIS raster file, such as a digital elevation model, are smaller than the true surface area, because they are projected to a planimetric plane. In mountainous regions this sloped area under‐estimation (SAUE) can have significant consequences on hydrological calculations. A sensitivity analysis is conducted, using the ACRU agro‐hydrological modelling system in a small watershed in Glacier National Park, Montana, USA, to investigate the sensitivity of the SAUE on key elements of the hydrological cycle, including precipitation depth, April snow depth, August soil moisture deficit, actual evapotranspiration depth, and runoff depth. The sensitivity analysis is based on 224 unique combinations of slope, soil and land cover types, elevation with associated precipitation depths, and north and south facing radiation regimes. Results revealed an increasing influence of the SAUE on all hydrological processes with increasing slope steepness. Distinct differences and magnitudes between different land cover types, different elevations, and, in particular, different exposition were quantified. Actual evapotranspiration increases with SAUE, while runoff decreases. April soil water is simulated to decrease with an increase in SAUE. Finally, a comparison of a streamflow simulation of a small and steep alpine watershed with and without consideration of the SAUE is carried out. The sloped area of the small watershed is under‐estimated by 20·9%, and the difference in simulated runoff is 12·3%. When the SAUE was not considered, runoff was simulated to be higher, the associated coefficient of determination was slightly lower, and the slope of the regression line was flatter. Copyright © 2010 John Wiley & Sons, Ltd.     

Impact of model simplifications on soil erosion predictions: application of the GLUE methodology to a distributed event‐based model at the hillslope scale

In this paper, we analyse how the performance and calibration of a distributed event‐based soil erosion model at the hillslope scale is affected by different simplifications on the parameterizations used to compute the production of suspended sediment by rainfall and runoff. Six modelling scenarios of different complexity are used to evaluate the temporal variability of the sedimentograph at the outlet of a 60 m long cultivated hillslope. The six scenarios are calibrated within the generalized likelihood uncertainty estimation framework in order to account for parameter uncertainty, and their performance is evaluated against experimental data registered during five storm events. The Nash–Sutcliffe efficiency, percent bias and coverage performance ratios show that the sedimentary response of the hillslope in terms of mass flux of eroded soil can be efficiently captured by a model structure including only two soil erodibility parameters, which control the rainfall and runoff production of suspended sediment. Increasing the number of parameters makes the calibration process more complex without increasing in a noticeable manner the predictive capability of the model. Copyright © 2015 John Wiley & Sons, Ltd.     

Seasonal evolution of meltwater generation,storage and discharge at a non‐temperate glacier in Svalbard

In glacierized catchments, meteorological inputs driving surface melting are translated into runoff outputs mediated by the glacier hydrological system: analysis of the relationship between meteorology and diurnal and seasonal patterns of runoff should reflect the functioning of that system, with the role of meltwater storage likely to be of particular importance. Daily meltwater storage is determined for a glacier at 78 °N in the Svalbard archipelago, by comparing inputs calculated from a surface energy balance model with measured outputs (proglacial discharge). Solar radiation, air temperature, wind speed and proglacial discharge are then analysed by regression and time‐series methods, in order to assess the meteorology–discharge relationship and its variation at diurnal and seasonal time‐scales. The recorded discharge time‐series can be divided into two contrasting intervals: up to early August, proglacial discharge was high and variable, mean hydrographs showed little indication of diurnal cycling, ARIMA models of discharge indicated a non‐seasonal, moving‐average generating process, and there was a net loss of meltwater from storage; from early August, proglacial discharge was low and relatively invariable, but with clearer diurnal cycles, regression models of discharge showed substantially improved correlations with air temperature and solar radiation, ARIMA models indicated a non‐seasonal, autoregressive generating process, and eventually a seasonal component, and there was a net gain in meltwater storage. The transition between the two periods is brief compared with the duration of the melt season. The runoff response to meteorology therefore lacks the strongly progressive element previously identified in mid‐latitude glacierized catchments. In particular, the glacier hydrological system only appears responsive to diurnal forcing following the depletion of the seasonal snowpack meltwater store. Copyright © 2001 John Wiley & Sons, Ltd.     

Use of multi‐platform,multi‐temporal remote‐sensing data for calibration of a distributed hydrological model: an application in the Arno basin,Italy

Images from satellite platforms are a valid aid in order to obtain distributed information about hydrological surface states and parameters needed in calibration and validation of the water balance and flood forecasting. Remotely sensed data are easily available on large areas and with a frequency compatible with land cover changes. In this paper, remotely sensed images from different types of sensor have been utilized as a support to the calibration of the distributed hydrological model MOBIDIC, currently used in the experimental system of flood forecasting of the Arno River Basin Authority. Six radar images from ERS‐2 synthetic aperture radar (SAR) sensors (three for summer 2002 and three for spring–summer 2003) have been utilized and a relationship between soil saturation indexes and backscatter coefficient from SAR images has been investigated. Analysis has been performed only on pixels with meagre or no vegetation cover, in order to legitimize the assumption that water content of the soil is the main variable that influences the backscatter coefficient. Such pixels have been obtained by considering vegetation indexes (NDVI) and land cover maps produced by optical sensors (Landsat‐ETM). In order to calibrate the soil moisture model based on information provided by SAR images, an optimization algorithm has been utilized to minimize the regression error between saturation indexes from model and SAR data and error between measured and modelled discharge flows. Utilizing this procedure, model parameters that rule soil moisture fluxes have been calibrated, obtaining not only a good match with remotely sensed data, but also an enhancement of model performance in flow prediction with respect to a previous calibration with river discharge data only. Copyright © 2006 John Wiley & Sons, Ltd.     

Investigation of the impacts of vegetation distribution and evaporative cooling on synthetic urban daytime climate using a coupled LES—LSM model

Local flow properties and regional weather or climate are strongly affected by land‐atmosphere interactions of momentum and scalars within the daytime convective boundary layer (CBL). In this study, we investigate the impact of green space scale on the daytime atmospheric boundary layer (ABL) over a synthetic urban domain using a recently developed large‐eddy simulation‐land surface model (LES–LSM) framework. With the use of realistic soundings as initial conditions, a series of numerical experiments over synthetic urban surfaces with varied scale of vegetated area is performed. Simulated micrometeorological properties, surface fluxes, basic CBL characteristics, and cloud distribution are analysed. The results show reference‐level air potential temperature and specific humidity as well as surface fluxes over green space are significantly affected by the scale of green space in the urban domain. The surface organization due to vegetated area scale also has impacts on horizontally averaged scalar and momentum profiles; however, the magnitude in this study is smaller than the results of a previous study using a set of offline surface fluxes as the lower boundary condition for LES. In addition, even though this study only performs a daytime diurnal cycle, the impact of green space scale on cloud distribution in simulations is significant. The cases with more organized green space yield lower‐elevated cumulus cloud and larger‐cloud cover fraction, which impacts the energy budget at the top of boundary layer and, in turn, could lead to additional surface cooling with respect to longer‐term weather and climate. Copyright © 2010 John Wiley & Sons, Ltd.     

Effects of hillslope topography on hydrological responses in a weathered granite mountain,Japan: comparison of the runoff response between the valley‐head and the side slope

To evaluate the effects of hillslope topography on storm runoff in a weathered granite mountain, discharge rate, soil pore water pressures, and water chemistry were observed on two types of hillslope: a valley‐head (a concave hillslope) and a side slope (a planar hillslope). Hydrological responses on the valley‐head and side slope reflected their respective topographic characteristics and varied with the rainfall magnitude. During small rainfall events (<35 mm), runoff from the side slope occurred rapidly relative to the valley‐head. The valley‐head showed little response in storm runoff. As rainfall amounts increased (35–60 mm), the valley‐head yielded a higher flow relative to the side slope. For large rainfall events (>60 mm), runoff from both hillslopes increased with rainfall, although that from the valley‐head was larger than that from the side slope. The differences in the runoff responses were caused by differences in the roles of lower‐slope soils and the convergence of the hillslope. During small rainfall events, the side slope could store little water; in contrast, all rainwater could be stored in the soils at the valley‐head hollow. As the amount of rainfall increased, the subsurface saturated area of the valley‐head extended from the bottom to the upper portion of the slope, with the contributions of transient groundwater via lateral preferential flowpaths due to the high concentration of subsurface water. Conversely, saturated subsurface flow did not contribute to runoff responses, and the subsurface saturated area at the side slope did not extend to the upper slope for the same storm size. During large rainfall events, expansion of the subsurface saturated area was observed in both hillslopes. Thus, differences in the concentration of subsurface water, reflecting hillslope topography, may create differences in the extension of the subsurface saturated area, as well as variability in runoff responses. Copyright © 2007 John Wiley & Sons, Ltd.     

Distributed and integrated response of a geographic information system‐based hydrologic model in the eastern Palouse region,Idaho

Verification of distributed hydrologic models is rare owing to the lack of spatially detailed field measurements and a common mismatch between the scale at which soil hydraulic properties are measured and the scale of a single modelling unit. In this study, two of the most commonly calibrated parameters, i.e. soil depth and the vertical distribution of lateral saturated hydraulic conductivity K_s, were eliminated by a spatially detailed soil characterization and results of a hillslope‐scale field experiment. The soil moisture routing (SMR) model, a geographic information system‐based hydrologic model, was modified to represent the dominant hydrologic processes for the Palouse region of northern Idaho. Modifications included K_s as a double exponential function of depth in a single soil layer, a snow accumulation and melt algorithm, and a simple relationship between storage and perched water depth (PWD) using the drainable porosity. The model was applied to a 2 ha catchment without calibration to measured data. Distributed responses were compared with observed PWD over a 3‐year period on a 10 m × 15 m grid. Integrated responses were compared with observed surface runoff at the catchment outlet. The modified SMR model simulated the PWD fluctuations remarkably well, especially considering the shallow soils in this catchment: a 0·20 m error in PWD is equivalent to only a 1·6% error in predicted soil moisture content. Simulations also captured PWD fluctuations during a year with high spatial variability of snow accumulation and snowmelt rates at upslope, mid‐slope, and toe slope positions with errors as low as 0·09 m, 0·12 m, and 0·12 m respectively. Errors in distributed and integrated model simulations were attributed mostly to misrepresentation of rain events and snowmelt timing problems. In one location in the catchment, simulated PWD was consistently greater than observed PWD, indicating a localized recharge zone, which was not identified by the soil morphological survey. Copyright © 2006 John Wiley & Sons, Ltd.     

Daily streamflow modelling and assessment based on the curve‐number technique

A cell‐based long‐term hydrological model (CELTHYM) that can be integrated with a geographical information system (GIS) was developed to predict continuous stream flow from small agricultural watersheds. The CELTHYM uses a cell‐by‐cell soil moisture balance approach. For surface runoff estimation, the curve number technique considering soil moisture on a daily basis was used, and release rate was used to estimate baseflow. Evapotranspiration was computed using the FAO modified Penman equation that considered land‐use‐based crop coefficients, soil moisture and the influence of topography on radiation. A rice paddy field water budget model was also adapted for the specific application of the model to East Asia. Model sensitivity analysis was conducted to obtain operational information about the model calibration parameters. The CELTHYM was calibrated and verified with measured runoff data from the WS#1 and WS#3 watersheds of the Seoul National University, Department of Agricultural Engineering, in Hwaseong County, Kyounggi Province, South Korea. The WS#1 watershed is comprised of about 35·4% rice paddy fields and 42·3% forest, whereas the WS#3 watershed is about 85·0% forest and 11·5% rice paddy fields. The CELTHYM was calibrated for the parameter release rate, K, and soil moisture storage coefficient, STC, and results were compared with the measured runoff data for 1986. The validation results for WS#1 considering all daily stream flow were poor with R², E² and RMSE having values of 0·40, ?6·63 and 9·69 (mm), respectively, but validation results for days without rainfall were statistically significant (R² = 0·66). Results for WS#3 showed good agreement with observed data for all days, and R², E² and RMSE were 0·92, 0·91 and 2·23 (mm), respectively, suggesting potential for CELTHYM application to other watersheds. The direct runoff and water balance components for watershed WS#1 with significant areas of paddy fields did not perform well, suggesting that additional study of these components is needed. Copyright © 2002 John Wiley & Sons, Ltd.     

Influence of soil‐surface types on the overall runoff of the Tabernas badlands (south‐east Spain): field data and model approaches

The Tabernas desert, an extensive badlands area in Almeria province (south‐east Spain), is characterized by a high variability in soil surface cover and soil properties along with important topographical contrasts giving rise to a wide range of hydrological behaviour. A double approach through field monitoring and modelling has been used to ascertain the influence of soil‐surface variability on the overall hydrological response. Small plots were monitored for 3 years to assess runoff from the different surface types. Data provided by the long‐term monitoring of three small catchments formed by different soil surfaces were used to find out the specific contribution of each soil surface to the catchment runoff. A simple spatially distributed model was built to predict runoff generation based on the infiltration rate of each soil‐surface type (defined as terrain units with the same cover, the same soil type and on the same landform). Plot results prove that the soil surface units within the study area behave differently in terms of hydrological response to natural rainfall. These responses are explained by the types of cover, topographical characteristics and soil properties. When runoff events are simple (with one or two runoff peaks), the modelled hydrographs reproduce the hydrographs observed reasonably well, but in complex events (with several runoff peaks) the adjustment is not as good. The model also shows the influence of the spatial distribution of soil surfaces on the overall runoff, aiding exploration of the spatial hydrological relationships among different landscape units. Copyright © 2002 John Wiley & Sons, Ltd.     

Impact of an inter‐basin water transfer and reservoir operation on a karst open streamflow hydrological regime: an example from the Dinaric karst (Croatia)

This paper describes the hydrological changes caused by inter‐basin water transfer and the reservoir development on the hydrological regimes of two rivers. The Sabljaki Reservoir in the Zagorska Mre?nica River and the Bukovik Reservoir in the upper Dobra River began operation in 1959. Both are part of the hydroelectric power plant (HEPP) Gojak, whose installed capacity is 50 m³/s. Their water volumes at the spillway altitudes of 320·10 and 320·15 m a. s. l. are 3·3 × 10⁶ and 0·24 × 10⁶ m³ respectively. Both the Dobra and Mre?nica Rivers are losing, sinking and underground karst rivers. A 9376‐m‐long tunnel provides water from the Sabljaki Reservoir to the HEPP Gojak, which was constructed in the Lower Dobra River. The Sabljaki Reservoir is located in the Pla?ki karst polje, while the Bukovik Reservoir is located in the neighbouring Ogulin karst polje. The consequences of the inter‐basin water transfer are strong and have caused abrupt changes in the hydrological regimes of the downstream sections of both rivers. At the same time, the construction and development of both the reservoirs have also caused hydrological changes to the upstream section of the Upper Dobra River. Copyright © 2010 John Wiley & Sons, Ltd.     

Multi‐scale relief model (MSRM): a new algorithm for the visualization of subtle topographic change of variable size in digital elevation models

Morphological analysis of landforms has traditionally relied on the interpretation of imagery. Although imagery provides a natural view of an area of interest (AOI) images are largely hindered by the environmental conditions at the time of image acquisition, the quality of the image and, mainly, the lack of topographical information, which is an essential factor for a correct understanding of the AOI's geomorphology. More recently digital surface models (DSMs) have been incorporated into the analytical toolbox of geomorphologists. These are usually high‐resolution models derived from digital photogrammetric processes or LiDAR data. However, these are restricted to relatively small areas and are expensive or complex to acquire, which limits widespread implementation. In this paper, we present the multi‐scale relief model (MSRM), which is a new algorithm for the visual interpretation of landforms using DSMs. The significance of this new method lies in its capacity to extract landform morphology from both high‐ and low‐resolution DSMs independently of the shape or scale of the landform under study. This method thus provides important advantages compared to previous approaches as it: (1) allows the use of worldwide medium resolution models, such as SRTM, ASTER GDEM, ALOS, and TanDEM‐X; (2) offers an alternative to traditional photograph interpretation that does not rely on the quality of the imagery employed nor on the environmental conditions and time of its acquisition; and (3) can be easily implemented for large areas using traditional GIS/RS software. The algorithm is tested in the Sutlej‐Yamuna interfluve, which is a very large low‐relief alluvial plain in northwest India where 10 000 km of palaeoriver channels have been mapped using MSRM. The code, written in Google Earth Engine's implementation of JavaScript, is provided as Supporting Information for its use in any other AOI without particular technical knowledge or access to topographical data. © 2017 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.