Assessing the impact of land use change on hydrology by ensemble modelling (LUCHEM) II: Ensemble combinations and predictions

This paper reports on a project to compare predictions from a range of catchment models applied to a mesoscale river basin in central Germany and to assess various ensemble predictions of catchment streamflow. The models encompass a large range in inherent complexity and input requirements. In approximate order of decreasing complexity, they are DHSVM, MIKE-SHE, TOPLATS, WASIM-ETH, SWAT, PRMS, SLURP, HBV, LASCAM and IHACRES. The models are calibrated twice using different sets of input data. The two predictions from each model are then combined by simple averaging to produce a single-model ensemble. The 10 resulting single-model ensembles are combined in various ways to produce multi-model ensemble predictions. Both the single-model ensembles and the multi-model ensembles are shown to give predictions that are generally superior to those of their respective constituent models, both during a 7-year calibration period and a 9-year validation period. This occurs despite a considerable disparity in performance of the individual models. Even the weakest of models is shown to contribute useful information to the ensembles they are part of. The best model combination methods are a trimmed mean (constructed using the central four or six predictions each day) and a weighted mean ensemble (with weights calculated from calibration performance) that places relatively large weights on the better performing models. Conditional ensembles, in which separate model weights are used in different system states (e.g. summer and winter, high and low flows) generally yield little improvement over the weighted mean ensemble. However a conditional ensemble that discriminates between rising and receding flows shows moderate improvement. An analysis of ensemble predictions shows that the best ensembles are not necessarily those containing the best individual models. Conversely, it appears that some models that predict well individually do not necessarily combine well with other models in multi-model ensembles. The reasons behind these observations may relate to the effects of the weighting schemes, non-stationarity of the climate series and possible cross-correlations between models.     

On the importance of soil moisture in calibration of rainfall–runoff models: two case studies

Streamflow modelling results from the GR4H and PDM hydrological models were evaluated in two Australian sub-catchments, using (1) calibration to streamflow and (2) joint-calibration to streamflow and soil moisture. Soil moisture storage in the models was evaluated against soil moisture observations from field measurements. The PDM had the best performance in terms of both streamflow and soil moisture estimations during the calibration period, but was outperformed by GR4H during validation. It was also shown that the soil moisture estimation was improved significantly by joint-calibration for the case where streamflow and soil moisture estimations were poor. In other cases, addition of the soil moisture constraint did not degrade the results. Consequently, it is recommended that GR4H be used, in preference to the PDM, in the foothills of the Murrumbidgee catchment or other Australian catchments with semi-arid to sub-humid climate, and that soil moisture data be used in the calibration process.     

Exploratory analysis of statistical post‐processing methods for hydrological ensemble forecasts

Despite many recent improvements, ensemble forecast systems for streamflow often produce under‐dispersed predictive distributions. This situation is problematic for their operational use in water resources management. Many options exist for post‐processing of raw forecasts. However, most of these have been developed using meteorological variables such as temperature, which displays characteristics very different from streamflow. In addition, streamflow data series are often very short or contain numerous gaps, thus compromising the estimation of post‐processing statistical parameters. For operational use, a post‐processing method has to be effective while remaining as simple as possible. We compared existing post‐processing methods using normally distributed and gamma‐distributed synthetic datasets. To reflect situations encountered with ensemble forecasts of daily streamflow, four normal distribution parameterizations and six gamma distribution parameterizations were used. Three kernel‐based approaches were tested, namely, the ‘best member’ method and two improvements thereof, and one regression‐based approach. Additional tests were performed to assess the ability of post‐processing methods to cope with short calibration series, missing values or small numbers of ensemble members. We thus found that over‐dispersion is best corrected by the regression method, while under‐dispersion is best corrected by kernel‐based methods. This work also shows key limitations associated with short data series, missing values, asymmetry and bias. One of the improved best member methods required longer series for the estimation of post‐processing parameters, but if provided with adequate information, yielded the best improvement of the continuous ranked probability score. These results suggest guidelines for future studies involving real operational datasets. Copyright © 2014 John Wiley & Sons, Ltd.     

Choosing an arbitrary calibration period for hydrologic models: How much does it influence water balance simulations?

The selection of calibration and validation time periods in hydrologic modelling is often done arbitrarily. Nonstationarity can lead to an optimal parameter set for one period which may not accurately simulate another. However, there is still much to be learned about the responses of hydrologic models to nonstationary conditions. We investigated how the selection of calibration and validation periods can influence water balance simulations. We calibrated Soil and Water Assessment Tool hydrologic models with observed streamflow for three United States watersheds (St. Joseph River of Indiana/Michigan, Escambia River of Florida/Alabama, and Cottonwood Creek of California), using time period splits for calibration/validation. We found that the choice of calibration period (with different patterns of observed streamflow, precipitation, and air temperature) influenced the parameter sets, leading to dissimilar simulations of water balance components. In the Cottonwood Creek watershed, simulations of 50-year mean January streamflow varied by 32%, because of lower winter precipitation and air temperature in earlier calibration periods on calibrated parameters, which impaired the ability for models calibrated to earlier periods to simulate later periods. Peaks of actual evapotranspiration for this watershed also shifted from April to May due to different parameter values depending on the calibration period's winter air temperatures. In the St. Joseph and Escambia River watersheds, adjustments of the runoff curve number parameter could vary by 10.7% and 20.8%, respectively, while 50-year mean monthly surface runoff simulations could vary by 23%–37% and 169%–209%, depending on the observed streamflow and precipitation of the chosen calibration period. It is imperative that calibration and validation time periods are chosen selectively instead of arbitrarily, for instance using change point detection methods, and that the calibration periods are appropriate for the goals of the study, considering possible broad effects of nonstationary time series on water balance simulations. It is also crucial that the hydrologic modelling community improves existing calibration and validation practices to better include nonstationary processes.     

Improving streamflow estimation in ungauged basins using a multi-modelling approach

ABSTRACT

In this study, a multi-modelling approach is proposed for improved continuous daily streamflow estimation in ungauged basins using regionalization—the process of transferring hydrological data from gauged to ungauged watersheds. Four regionalization models, two data-driven and two hydrological, were used for continuous daily streamflow estimation. Comparison of the individual models reveals that each of the four models performed well on a limited number of ungauged basins while none of them performed well for the entire 90 selected watersheds. The results obtained from the four models are evaluated and reported in a deterministic way by a model combination approach along with its uncertainty range consisting of 16 ensemble members. It is shown that a combined model of the four individual models performed well on all 90 watersheds and the ensemble range can account for the uncertainty of models. The combined model was more efficient and appeared more robust compared to the individual models. Furthermore, continuous ranked probability scores (CRPS) calculated for the ensemble model outputs indicate better performance compared to individual models and competitive with the combined model.

EDITOR A. Castellarin ASSOCIATE EDITOR G. Di Baldassarre     

Parameter identification and uncertainty analysis for simulating streamflow in a river basin of Eastern India

Uncertainty is inherent in modelling studies. However, the quantification of uncertainties associated with a model is a challenging task, and hence, such studies are somewhat limited. As distributed or semi‐distributed hydrological models are being increasingly used these days to simulate hydrological processes, it is vital that these models should be equipped with robust calibration and uncertainty analysis techniques. The goal of the present study was to calibrate and validate the Soil and Water Assessment Tool (SWAT) model for simulating streamflow in a river basin of Eastern India, and to evaluate the performance of salient optimization techniques in quantifying uncertainties. The SWAT model for the study basin was developed and calibrated using Parameter Solution (ParaSol), Sequential Uncertainty Fitting Algorithm (SUFI‐2) and Generalized Likelihood Uncertainty Estimation (GLUE) optimization techniques. The daily observed streamflow data from 1998 to 2003 were used for model calibration, and those for 2004–2005 were used for model validation. Modelling results indicated that all the three techniques invariably yield better results for the monthly time step than for the daily time step during both calibration and validation. The model performances for the daily streamflow simulation using ParaSol and SUFI‐2 during calibration are reasonably good with a Nash–Sutcliffe efficiency and mean absolute error (MAE) of 0.88 and 9.70 m³/s for ParaSol, and 0.86 and 10.07 m³/s for SUFI‐2, respectively. The simulation results of GLUE revealed that the model simulates daily streamflow during calibration with the highest accuracy in the case of GLUE (R² = 0.88, MAE = 9.56 m³/s and root mean square error = 19.70 m³/s). The results of uncertainty analyses by SUFI‐2 and GLUE were compared in terms of parameter uncertainty. It was found that SUFI‐2 is capable of estimating uncertainties in complex hydrological models like SWAT, but it warrants sound knowledge of the parameters and their effects on the model output. On the other hand, GLUE predicts more reliable uncertainty ranges (R‐factor = 0.52 for daily calibration and 0.48 for validation) compared to SUFI‐2 (R‐factor = 0.59 for daily calibration and 0.55 for validation), though it is computationally demanding. Although both SUFI‐2 and GLUE appear to be promising techniques for the uncertainty analysis of modelling results, more and more studies in this direction are required under varying agro‐climatic conditions for assessing their generic capability. Copyright © 2015 John Wiley & Sons, Ltd.     

Dual state–parameter estimation of hydrological models using ensemble Kalman filter

Hydrologic models are twofold: models for understanding physical processes and models for prediction. This study addresses the latter, which modelers use to predict, for example, streamflow at some future time given knowledge of the current state of the system and model parameters. In this respect, good estimates of the parameters and state variables are needed to enable the model to generate accurate forecasts. In this paper, a dual state–parameter estimation approach is presented based on the Ensemble Kalman Filter (EnKF) for sequential estimation of both parameters and state variables of a hydrologic model. A systematic approach for identification of the perturbation factors used for ensemble generation and for selection of ensemble size is discussed. The dual EnKF methodology introduces a number of novel features: (1) both model states and parameters can be estimated simultaneously; (2) the algorithm is recursive and therefore does not require storage of all past information, as is the case in the batch calibration procedures; and (3) the various sources of uncertainties can be properly addressed, including input, output, and parameter uncertainties. The applicability and usefulness of the dual EnKF approach for ensemble streamflow forecasting is demonstrated using a conceptual rainfall-runoff model.     

Calibration and uncertainty analysis of the SWAT model using Genetic Algorithms and Bayesian Model Averaging

In this paper, the Genetic Algorithms (GA) and Bayesian Model Averaging (BMA) were used to simultaneously conduct calibration and uncertainty analysis for the Soil and Water Assessment Tool (SWAT). In this combined method, several SWAT models with different structures are first selected; next GA is used to calibrate each model using observed streamflow data; finally, BMA is applied to combine the ensemble predictions and provide uncertainty interval estimation. This method was tested in two contrasting basins, the Little River Experimental Basin in Georgia, USA, and the Yellow River Headwater Basin in China. The results obtained in the two case studies show that this combined method can provide deterministic predictions better than or comparable to the best calibrated model using GA. The 66.7% and 90% uncertainty intervals estimated by this method were analyzed. The differences between the percentage of coverage of observations and the corresponding expected coverage percentage are within 10% for both calibration and validation periods in these two test basins. This combined methodology provides a practical and flexible tool to attain reliable deterministic simulation and uncertainty analysis of SWAT.     

Multi‐variable and multi‐site calibration and validation of SWAT in a large mountainous catchment with high spatial variability

Many methods developed for calibration and validation of physically based distributed hydrological models are time consuming and computationally intensive. Only a small set of input parameters can be optimized, and the optimization often results in unrealistic values. In this study we adopted a multi‐variable and multi‐site approach to calibration and validation of the Soil Water Assessment Tool (SWAT) model for the Motueka catchment, making use of extensive field measurements. Not only were a number of hydrological processes (model components) in a catchment evaluated, but also a number of subcatchments were used in the calibration. The internal variables used were PET, annual water yield, daily streamflow, baseflow, and soil moisture. The study was conducted using an 11‐year historical flow record (1990–2000); 1990–94 was used for calibration and 1995–2000 for validation. SWAT generally predicted well the PET, water yield and daily streamflow. The predicted daily streamflow matched the observed values, with a Nash–Sutcliffe coefficient of 0·78 during calibration and 0·72 during validation. However, values for subcatchments ranged from 0·31 to 0·67 during calibration, and 0·36 to 0·52 during validation. The predicted soil moisture remained wet compared with the measurement. About 50% of the extra soil water storage predicted by the model can be ascribed to overprediction of precipitation; the remaining 50% discrepancy was likely to be a result of poor representation of soil properties. Hydrological compensations in the modelling results are derived from water balances in the various pathways and storage (evaporation, streamflow, surface runoff, soil moisture and groundwater) and the contributions to streamflow from different geographic areas (hill slopes, variable source areas, sub‐basins, and subcatchments). The use of an integrated multi‐variable and multi‐site method improved the model calibration and validation and highlighted the areas and hydrological processes requiring greater calibration effort. Copyright © 2005 John Wiley & Sons, Ltd.     

赣江流域多种数据同化方案的径流模拟比较

为了探讨水文模型在不同水文数据同化方案下的径流模拟差异,本文采用集合卡尔曼滤波算法,以遥感蒸散发产品、实测径流为观测数据,构建了基于新安江模型的数据同化框架。基于此框架设计了4种不同同化方案(DA-ET、DAET(K)、DA-ET-Q、DA-ET-Q(K))以及1种对照方案OL,以赣江流域开展实例研究,评估了水文数据同化中遥感蒸散发产品的时间分辨率、模型蒸散发相关参数时变与否以及多源数据同化对径流模拟的影响。结果表明:在DA-ET方案下,同化两种不同时间分辨率的蒸散发产品均能提高模型整体的径流模拟精度,且时间分辨率更高的产品的同化效果更好;在DA-ET方案的基础上,考虑加入实测径流进行同化能够提升模型径流模拟精度,且DA-ET(K)与DA-ET-Q(K)方案所得径流相对误差的减幅均超过了20%,说明在蒸散发同化过程中同时考虑蒸散发参数动态变化的结果更优;相较于OL方案,4种同化方案均能不同程度地提高模型对径流高水部分的模拟能力,但DA-ET-Q(K)方案表现最差,而其余方案差异并不显著。本研究有助于进一步了解不同数据同化方案在径流模拟中的差异,从而为水资源高效利用与科学管理提供科学依据...     

Evaluation of distributed hydrologic impacts of temperature-index and energy-based snow models

Two commonly used strategies in modeling snowmelt are the energy balance and temperature-index methods. Here we evaluate the distributed hydrologic impacts of these two different snowmelt modeling strategies, each in conjunction with a physics-based hydrologic model (PIHM). Results illustrate that both the Isnobal energy-balance and calibrated temperature-index methods adequately reproduce snow depletion at the observation site. However, the models exhibit marked differences in the distribution of snowmelt. When combined with PIHM, both models capture streamflow reasonably during calibration year (WY06), but Isnobal model gives better streamflow results in the validation year (WY07). The uncalibrated temperature-index model predicts streamflow poorly in both years. Differences between distributed snowmelt, as predicted by Isnobal and calibrated temperature-index method, and its consequent effect on predicted hydrologic states suggest the need to carefully calibrate temperature-index models in both time and space. Combined physics-based snow and hydrologic models provide the best accuracy, while a temperature-index model using coefficients from the literature the poorest.     

Leveraging ensemble meteorological forcing data to improve parameter estimation of hydrologic models

As continental to global scale high-resolution meteorological datasets continue to be developed, there are sufficient meteorological datasets available now for modellers to construct a historical forcing ensemble. The forcing ensemble can be a collection of multiple deterministic meteorological datasets or come from an ensemble meteorological dataset. In hydrological model calibration, the forcing ensemble can be used to represent forcing data uncertainty. This study examines the potential of using the forcing ensemble to identify more robust parameters through model calibration. Specifically, we compare an ensemble forcing-based calibration with two deterministic forcing-based calibrations and investigate their flow simulation and parameter estimation properties and the ability to resist poor-quality forcings. The comparison experiment is conducted with a six-parameter hydrological model for 30 synthetic studies and 20 real data studies to provide a better assessment of the average performance of the deterministic and ensemble forcing-based calibrations. Results show that the ensemble forcing-based calibration generates parameter estimates that are less biased and have higher frequency of covering the true parameter values than the deterministic forcing-based calibration does. Using a forcing ensemble in model calibration reduces the risk of inaccurate flow simulation caused by poor-quality meteorological inputs, and improves the reliability and overall simulation skill of ensemble simulation results. The poor-quality meteorological inputs can be effectively filtered out via our ensemble forcing-based calibration methodology and thus discarded in any post-calibration model applications. The proposed ensemble forcing-based calibration method can be considered as a more generalized framework to include parameter and forcing uncertainties in model calibration.     

Improving monthly streamflow prediction in alpine regions: integrating HBV model with Bayesian neural network

Statistical methods have been widely used to build different streamflow prediction models; however, lacking of physical mechanism prevents precise streamflow prediction in alpine regions dominated by rainfall, snow and glacier. To improve precision, a new hybrid model (HBNN) integrating HBV hydrological model, Bayesian neural network (BNN) and uncertainty analysis is proposed. In this approach, the HBV is mainly used to generate initial snow-melt and glacier-melt runoffs that are regarded as new inputs of BNN for precision improvement. To examine model reliability, a hybrid deterministic model called HLSSVM incorporating the HBV model and least-square support vector machine is also developed and compared with HBNN in a typical region, the Yarkant River basin in Central Asia. The findings suggest that the HBNN model is a robust streamflow prediction model for alpine regions and capable of combining strengths of both the BNN statistical model and the HBV hydrological model, providing not only more precise streamflow prediction but also more reasonable uncertainty intervals than competitors particularly at high flows. It can be used in predicting streamflow for similar regions worldwide.     

Post‐processing of ensemble forecasts in low‐flow period

For water supply, navigational, ecological protection or water quality control purposes, there is a great need in knowing the likelihood of the river level falling below a certain threshold. Ensemble streamflow prediction (ESP) based on simulations of deterministic hydrologic models is widely used to assess this likelihood. Raw ESP results can be biased in both the ensemble means and the spreads. In this study, we applied a modified general linear model post‐processor (GLMPP) to correct these biases. The modified GLMPP is built on the basis of regression of simulated and observed streamflow calculated on the basis of canonical events, instead of the daily values as is carried out in the original GLMPP. We conducted the probabilistic analysis of post‐processed ESP results falling below pre‐specified low‐flow levels at seasonal time scale. Raw ESP forecasts from the 1980 to 2006 periods by four different land surface models (LSMs) in eight large river basins in the continental USA are included in the analysis. The four LSMs are Noah, Mosaic, variable infiltration capacity and Sacramento models. The major results from this study are as follows: (1) a modified GLMPP was proposed on the basis of canonical events; (2) post‐processing can improve the accuracy and reduce the uncertainty of hydrologic forecasts; (3) post‐processing can help deal with the effect of human activity; and (4) raw simulation results from different models vary greatly in different basins. However, post‐processing can always remove model biases under different conditions. Copyright © 2014 John Wiley & Sons, Ltd.     

How informative are stream level observations in different geographic regions?

Simple runoff models with a low number of model parameters are generally able to simulate catchment runoff reasonably well, but they rely on model calibration, which makes their use in ungauged basins challenging. In a previous study it has been shown that a limited number of streamflow measurements can be quite informative for constraining runoff models. In practice, however, instead of performing such repeated flow measurements, it might be easier to install a stream level logger. Here, a dataset of 600+ gauged basins in the USA was used to study how well models perform when only stream level data, rather than streamflow data, are available. A runoff model (the HBV model) was calibrated assuming that only stream level observations were available, and the simulations were evaluated on the full observed streamflow record. The results indicate that stream level data alone can already provide surprisingly good model simulation results in humid catchments, whereas in arid catchments some form of quantitative information (e.g. a streamflow observation or a regional average value) is needed to obtain good results. These results are encouraging for hydrological observations in data scarce regions as level observations are much easier to obtain than streamflow measurements. Based on runoff modelling, it might even be possible to derive streamflow time series from the level data obtained from loggers, satellites or community‐based approaches. Copyright © 2016 John Wiley & Sons, Ltd.     

Effect of single and multi-site calibration techniques on hydrological model performance,parameter estimation and predictive uncertainty: a case study in the Logone catchment,Lake Chad basin

Understanding hydrological processes at catchment scale through the use of hydrological model parameters is essential for enhancing water resource management. Given the difficulty of using lump parameters to calibrate distributed catchment hydrological models in spatially heterogeneous catchments, a multiple calibration technique was adopted to enhance model calibration in this study. Different calibration techniques were used to calibrate the Soil and Water Assessment Tool (SWAT) model at different locations along the Logone river channel. These were: single-site calibration (SSC); sequential calibration (SC); and simultaneous multi-site calibration (SMSC). Results indicate that it is possible to reveal differences in hydrological behavior between the upstream and downstream parts of the catchment using different parameter values. Using all calibration techniques, model performance indicators were mostly above the minimum threshold of 0.60 and 0.65 for Nash Sutcliff Efficiency (NSE) and coefficient of determination (R ²) respectively, at both daily and monthly time-steps. Model uncertainty analysis showed that more than 60% of observed streamflow values were bracketed within the 95% prediction uncertainty (95PPU) band after calibration and validation. Furthermore, results indicated that the SC technique out-performed the other two methods (SSC and SMSC). It was also observed that although the SMSC technique uses streamflow data from all gauging stations during calibration and validation, thereby taking into account the catchment spatial variability, the choice of each calibration method will depend on the application and spatial scale of implementation of the modelling results in the catchment.     

A multi-objective calibration approach using in-situ soil moisture data for improved hydrological simulation of the Prairies

ABSTRACT

Traditionally, hydrological models are only calibrated to reproduce streamflow regime without considering other hydrological state variables, such as soil moisture and evapotranspiration. Limited studies have been performed on constraining the model parameters, despite the fact that the presence of a large number of parameters may provide large degree of freedom, resulting in equifinality and poor model performance. In this study, a multi-objective optimization approach is adopted, and both streamflow and soil moisture data are calibrated simultaneously for an experimental study basin in the Saskatchewan Prairies in western Canada. The results of this study show that the multi-objective calibration improves model fidelity compared to the single objective calibration. Moreover, the study demonstrates that single objective calibration performed against only streamflow can fairly mimic the streamflow hydrograph but does not yield realistic estimation of other fluxes such as evapotranspiration and soil moisture (especially in deeper soil layers).     

A hybrid factorial stepwise-cluster analysis method for streamflow simulation – a case study in northwestern China

ABSTRACT

In this study, a hybrid factorial stepwise-cluster analysis (HFSA) method is developed for modelling hydrological processes. The HFSA method employs a cluster tree to represent the complex nonlinear relationship between inputs (predictors) and outputs (predictands) in hydrological processes. A real case of streamflow simulation for the Kaidu River basin is applied to demonstrate the efficiency of the HFSA method. After training a total of 24?108 calibration samples, the cluster tree for daily streamflow is generated based on a stepwise-cluster analysis (SCA) approach and is then used to reproduce the daily streamflows for calibration (1995–2005) and validation (2008–2010) periods. The Nash-Sutcliffe coefficients for calibration and validation are 0.68 and 0.65, respectively, and the deviations of volume are 1.68% and 4.11%, respectively. Results show that: (i) the HFSA method can formulate a SCA-based hydrological modelling system for streamflow simulation with a satisfactory fitting; (ii) the variability and peak value of streamflow in the Kaidu River basin can be effectively captured by the SCA-based hydrological modelling system; (iii) results from 2⁶ factorial experiments indicate that not only are minimum temperature and precipitation key drivers of system performance, but also the interaction between precipitation and minimum temperature significantly impacts on the streamflow. The findings are useful in indicating that the streamflow of the study basin is a mixture of snowmelt and rainfall water.

EDITOR D. Koutsoyiannis; ASSOCIATE EDITOR G. Thirel     

A One-dimensional Ensemble Forecast and Assimilation System for Fog Prediction

A probabilistic fog forecast system was designed based on two high resolution numerical 1-D models called COBEL and PAFOG. The 1-D models are coupled to several 3-D numerical weather prediction models and thus are able to consider the effects of advection. To deal with the large uncertainty inherent to fog forecasts, a whole ensemble of 1-D runs is computed using the two different numerical models and a set of different initial conditions in combination with distinct boundary conditions. Initial conditions are obtained from variational data assimilation, which optimally combines observations with a first guess taken from operational 3-D models. The design of the ensemble scheme computes members that should fairly well represent the uncertainty of the current meteorological regime. Verification for an entire fog season reveals the importance of advection in complex terrain. The skill of 1-D fog forecasts is significantly improved if advection is considered. Thus the probabilistic forecast system has the potential to support the forecaster and therefore to provide more accurate fog forecasts.     

Lateral terrestrial water flow contribution to summer precipitation at continental scale – A comparison between Europe and West Africa with WRF-Hydro-tag ensembles

It is well accepted that summer precipitation can be altered by soil moisture condition. Coupled land surface – atmospheric models have been routinely used to quantify soil moisture – precipitation feedback processes. However, most of the land surface models (LSMs) assume a vertical soil water transport and neglect lateral terrestrial water flow at the surface and in the subsurface, which potentially reduces the realism of the simulated soil moisture – precipitation feedback. In this study, the contribution of lateral terrestrial water flow to summer precipitation is assessed in two different climatic regions, Europe and West Africa, for the period June–September 2008. A version of the coupled atmospheric-hydrological model WRF-Hydro with an option to tag and trace land surface evaporation in the modelled atmosphere, named WRF-Hydro-tag, is employed. An ensemble of 30 simulations with terrestrial routing and 30 simulations without terrestrial routing is generated with random realizations of turbulent energy with the stochastic kinetic energy backscatter scheme, for both Europe and West Africa. The ensemble size allows to extract random noise from continental-scale averaged modelled precipitation. It is found that lateral terrestrial water flow increases the relative contribution of land surface evaporation to precipitation by 3.6% in Europe and 5.6% in West Africa, which enhances a positive soil moisture – precipitation feedback and generates more uncertainty in modelled precipitation, as diagnosed by a slight increase in normalized ensemble spread. This study demonstrates the small but non-negligible contribution of lateral terrestrial water flow to precipitation at continental scale.