Extraction of information content from stochastic disaggregation and bias corrected downscaled precipitation variables for crop simulation

We applied a simple statistical downscaling procedure for transforming daily global climate model (GCM) rainfall to the scale of an agricultural experimental station in Katumani, Kenya. The transformation made was two-fold. First, we corrected the rainfall frequency bias of the climate model by truncating its daily rainfall cumulative distribution into the station’s distribution based on a prescribed observed wet-day threshold. Then, we corrected the climate model rainfall intensity bias by mapping its truncated rainfall distribution into the station’s truncated distribution. Further improvements were made to the bias corrected GCM rainfall by linking it with a stochastic disaggregation scheme to correct the time structure problem inherent with daily GCM rainfall. Results of the simple and hybridized GCM downscaled precipitation variables (total, probability of occurrence, intensity and dry spell length) were linked with a crop model for a more objective evaluation of their performance using a non-linear measure based on mutual information based on entropy. This study is useful for the identification of both suitable downscaling technique as well as the effective precipitation variables for forecasting crop yields using GCM’s outputs which can be useful for addressing food security problems beforehand in critical basins around the world.     

Impact of the spatial variability of daily precipitation on hydrological projections: A comparison of GCM‐ and RCM‐driven cases in the Han River basin,Korea

In this study, we investigate the impact of the spatial variability of daily precipitation on hydrological projections based on a comparative assessment of streamflow simulations driven by a global climate model (GCM) and two regional climate models (RCMs). A total of 12 different climate input datasets, that is, the raw and bias‐corrected GCM and raw and bias‐corrected two RCMs for the reference and future periods, are fed to a semidistributed hydrological model to assess whether the bias correction using quantile mapping and dynamical downscaling using RCMs can improve streamflow simulation in the Han River basin, Korea. A statistical analysis of the daily precipitation demonstrates that the precipitation simulated by the GCM fails to capture the large variability of the observed daily precipitation, in which the spatial autocorrelation decreases sharply within a relatively short distance. However, the spatial variability of precipitation simulated by the two RCMs shows better agreement with the observations. After applying bias correction to the raw GCM and raw RCMs outputs, only a slight change is observed in the spatial variability, whereas an improvement is observed in the precipitation intensity. Intensified precipitation but with the same spatial variability of the raw output from the bias‐corrected GCM does not improve the heterogeneous runoff distributions, which in turn regulate unrealistically high peak downstream streamflow. GCM‐simulated precipitation with a large bias correction that is necessary to compensate for the poor performance in present climate simulation appears to distort streamflow patterns in the future projection, which leads to misleading projections of climate change impacts on hydrological extremes.     

Uncertainty analysis of statistical downscaling methods using Canadian Global Climate Model predictors

Three downscaling models, namely the Statistical Down‐Scaling Model (SDSM), the Long Ashton Research Station Weather Generator (LARS‐WG) model and an artificial neural network (ANN) model, have been compared in terms of various uncertainty attributes exhibited in their downscaled results of daily precipitation, daily maximum and minimum temperature. The uncertainty attributes are described by the model errors and the 95% confidence intervals in the estimates of means and variances of downscaled data. The significance of those errors has been examined by suitable statistical tests at the 95% confidence level. The 95% confidence intervals in the estimates of means and variances of downscaled data have been estimated using the bootstrapping method and compared with the observed data. The study has been carried out using 40 years of observed and downscaled daily precipitation data and daily maximum and minimum temperature data, starting from 1961 to 2000. In all the downscaling experiments, the simulated predictors of the Canadian Global Climate Model (CGCM1) have been used. The uncertainty assessment results indicate that, in daily precipitation downscaling, the LARS‐WG model errors are significant at the 95% confidence level only in a very few months, the SDSM errors are significant in some months, and the ANN model errors are significant in almost all months of the year. In downscaling daily maximum and minimum temperature, the performance of all three models is similar in terms of model errors evaluation at the 95% confidence level. But, according to the evaluation of variability and uncertainty in the estimates of means and variances of downscaled precipitation and temperature, the performances of the LARS‐WG model and the SDSM are almost similar, whereas the ANN model performance is found to be poor in that consideration. Further assessment of those models, in terms of skewness and average dry‐spell length comparison between observed and downscaled daily precipitation, indicates that the downscaled daily precipitation skewness and average dry‐spell lengths of the LARS‐WG model and the SDSM are closer to the observed data, whereas the ANN model downscaled precipitation underestimated those statistics in all months. Copyright © 2006 John Wiley & Sons, Ltd.     

Climate change impact assessment of a river basin using CMIP5 climate models and the VIC hydrological model

Climate change has significant impacts on water availability in larger river basins. The present study evaluates the possible impacts of projected future daily rainfall (2011–2099) on the hydrology of a major river basin in peninsular India, the Godavari River Basin, (GRB), under RCP4.5 and RCP8.5 scenarios. The study highlights a criteria-based approach for selecting the CMIP5 GCMs, based on their fidelity in simulating the Indian Summer Monsoon rainfall. The nonparametric kernel regression based statistical downscaling model is employed to project future daily rainfall and the variable infiltration capacity (VIC) macroscale hydrological model is used for hydrological simulations. The results indicate an increase in future rainfall without significant change in the spatial pattern of hydrological variables in the GRB. The climate-change-induced projected hydrological changes provide a crucial input to define water resource policies in the GRB. This methodology can be adopted for the climate change impacts assessment of larger river basins worldwide.     

Downscaling of precipitation for climate change scenarios: A support vector machine approach

The Climate impact studies in hydrology often rely on climate change information at fine spatial resolution. However, general circulation models (GCMs), which are among the most advanced tools for estimating future climate change scenarios, operate on a coarse scale. Therefore the output from a GCM has to be downscaled to obtain the information relevant to hydrologic studies. In this paper, a support vector machine (SVM) approach is proposed for statistical downscaling of precipitation at monthly time scale. The effectiveness of this approach is illustrated through its application to meteorological sub-divisions (MSDs) in India. First, climate variables affecting spatio-temporal variation of precipitation at each MSD in India are identified. Following this, the data pertaining to the identified climate variables (predictors) at each MSD are classified using cluster analysis to form two groups, representing wet and dry seasons. For each MSD, SVM- based downscaling model (DM) is developed for season(s) with significant rainfall using principal components extracted from the predictors as input and the contemporaneous precipitation observed at the MSD as an output. The proposed DM is shown to be superior to conventional downscaling using multi-layer back-propagation artificial neural networks. Subsequently, the SVM-based DM is applied to future climate predictions from the second generation Coupled Global Climate Model (CGCM2) to obtain future projections of precipitation for the MSDs. The results are then analyzed to assess the impact of climate change on precipitation over India. It is shown that SVMs provide a promising alternative to conventional artificial neural networks for statistical downscaling, and are suitable for conducting climate impact studies.     

A hybrid stochastic-weather-generation method for temporal disaggregation of precipitation with consideration of seasonality and within-month variations

Linking atmospheric and hydrological models is challenging because of a mismatch of spatial and temporal resolutions in which the models operate: dynamic hydrological models need input at relatively fine temporal (daily) scale, but the outputs from general circulation models are usually not realistic at the same scale, even though fine scale outputs are available. Temporal dimension downscaling methods called disaggregation are designed to produce finer temporal-scale data from reliable larger temporal-scale data. Here, we investigate a hybrid stochastic weather-generation method to simulate a high-frequency (daily) precipitation sequence based on lower frequency (monthly) amounts. To deal with many small precipitation amounts and capture large amounts, we divide the precipitation amounts on rainy days (with non-zero precipitation amounts) into two states (named moist and wet states, respectively) by a pre-defined threshold and propose a multi-state Markov chain model for the occurrences of different states (also including non-rain days called dry state). The truncated Gamma and censored extended Burr XII distributions are then employed to model the precipitation amounts in the moist and wet states, respectively. This approach avoids the need to deal with discontinuity in the distribution, and ensures that the states (dry, moist and wet) and corresponding amounts in rainy days are well matched. The method also considers seasonality by constructing individual models for different months, and monthly variation by incorporating the low-frequency amounts as a model predictor. The proposed method is compared with existing models using typical catchment data in Australia with different climate conditions (non-seasonal rainfall, summer rainfall and winter rainfall patterns) and demonstrates better performances under several evaluation criteria which are important in hydrological studies.     

Assessment of hydrologic impacts of climate change in Tunga–Bhadra river basin,India with HEC‐HMS and SDSM

Climate change would significantly affect many hydrologic systems, which in turn would affect the water availability, runoff, and the flow in rivers. This study evaluates the impacts of possible future climate change scenarios on the hydrology of the catchment area of the Tunga–Bhadra River, upstream of the Tungabhadra dam. The Hydrologic Engineering Center's Hydrologic Modeling System version 3.4 (HEC‐HMS 3.4) is used for the hydrological modelling of the study area. Linear‐regression‐based Statistical DownScaling Model version 4.2 (SDSM 4.2) is used to downscale the daily maximum and minimum temperature, and daily precipitation in the four sub‐basins of the study area. The large‐scale climate variables for the A2 and B2 scenarios obtained from the Hadley Centre Coupled Model version 3 are used. After model calibration and testing of the downscaling procedure, the hydrological model is run for the three future periods: 2011–2040, 2041–2070, and 2071–2099. The impacts of climate change on the basin hydrology are assessed by comparing the present and future streamflow and the evapotranspiration estimates. Results of the water balance study suggest increasing precipitation and runoff and decreasing actual evapotranspiration losses over the sub‐basins in the study area. Copyright © 2012 John Wiley & Sons, Ltd.     

Validation of hydrological models for climate scenario simulation: the case of Saguenay watershed in Quebec

This paper presents the results of an investigation into the problems associated with using downscaled meteorological data for hydrological simulations of climate scenarios. The influence of both the hydrological models and the meteorological inputs driving these models on climate scenario simulation studies are investigated. A regression‐based statistical tool (SDSM) is used to downscale the daily precipitation and temperature data based on climate predictors derived from the Canadian global climate model (CGCM1), and two types of hydrological model, namely the physically based watershed model WatFlood and the lumped‐conceptual modelling system HBV‐96, are used to simulate the flow regimes in the major rivers of the Saguenay watershed in Quebec. The models are validated with meteorological inputs from both the historical records and the statistically downscaled outputs. Although the two hydrological models demonstrated satisfactory performances in simulating stream flows in most of the rivers when provided with historic precipitation and temperature records, both performed less well and responded differently when provided with downscaled precipitation and temperature data. By demonstrating the problems in accurately simulating river flows based on downscaled data for the current climate, we discuss the difficulties associated with downscaling and hydrological models used in estimating the possible hydrological impact of climate change scenarios. Copyright © 2007 John Wiley & Sons, Ltd.     

Projection of climate change for daily precipitation: a case study in Shih‐Men reservoir catchment in Taiwan

A statistical downscaling model, based on the outputs of general circulation models (GCMs) as predictors, was proposed to simulate the daily precipitations in the Shih‐Men reservoir catchment in Taiwan. The structure of the proposed downscaling model is composed of two parts: classification and regression. Predictors of classification and regression models were selected from the large‐scale weather variables in the National Centers for Environmental Prediction (NECP) reanalysis data based on statistical tests. Discriminant analysis and support vector regression (SVR) were applied to build the classification and regression models. The outputs of five atmosphere‐ocean GCMs, which are reported to have properly considered tropical cyclone information and East Asian Monsoon modelling, were used for projecting future precipitations. Data from four grids covering Taiwan were used for developing the downscaling model. The potential of the downscaling models in simulating local precipitations was evaluated, and downscaling results reveal that the proposed downscaling model can reproduce local daily precipitations from large‐scale weather variables. Projected local precipitations under two emission scenarios show that the precipitations in the study area tend to decrease. Copyright © 2010 John Wiley & Sons, Ltd.     

Large scale climate and rainfall seasonality in a Mediterranean Area: Insights from a non‐homogeneous Markov model applied to the Agro‐Pontino plain

In the context of climate change and variability, there is considerable interest in how large scale climate indicators influence regional precipitation occurrence and its seasonality. Seasonal and longer climate projections from coupled ocean–atmosphere models need to be downscaled to regional levels for hydrologic applications, and the identification of appropriate state variables from such models that can best inform this process is also of direct interest. Here, a Non‐Homogeneous Hidden Markov Model (NHMM) for downscaling daily rainfall is developed for the Agro‐Pontino Plain, a coastal reclamation region very vulnerable to changes of hydrological cycle. The NHMM, through a set of atmospheric predictors, provides the link between large scale meteorological features and local rainfall patterns. Atmospheric data from the NCEP/NCAR archive and 56‐years record (1951–2004) of daily rainfall measurements from 7 stations in Agro‐Pontino Plain are analyzed. A number of validation tests are carried out, in order to: 1) identify the best set of atmospheric predictors to model local rainfall; 2) evaluate the model performance to capture realistically relevant rainfall attributes as the inter‐annual and seasonal variability, as well as average and extreme rainfall patterns. Validation tests show that the best set of atmospheric predictors are the following: mean sea level pressure, temperature at 1000 hPa, meridional and zonal wind at 850 hPa and precipitable water, from 20°N to 80°N of latitude and from 80°W to 60°E of longitude. Furthermore, the validation tests show that the rainfall attributes are simulated realistically and accurately. The capability of the NHMM to be used as a forecasting tool to quantify changes of rainfall patterns forced by alteration of atmospheric circulation under climate change and variability scenarios is discussed.     

Comparison of daily and sub-daily SWAT models for daily streamflow simulation in the Upper Huai River Basin of China

Despite the significant role of precipitation in the hydrological cycle, few studies have been conducted to evaluate the impacts of the temporal resolution of rainfall inputs on the performance of SWAT (soil and water assessment tool) models in large-sized river basins. In this study, both daily and hourly rainfall observations at 28 rainfall stations were used as inputs to SWAT for daily streamflow simulation in the Upper Huai River Basin. Study results have demonstrated that the SWAT model with hourly rainfall inputs performed better than the model with daily rainfall inputs in daily streamflow simulation, primarily due to its better capability of simulating peak flows during the flood season. The sub-daily SWAT model estimated that 58 % of streamflow was contributed by baseflow compared to 34 % estimated by the daily model. Using the future daily and 3-h precipitation projections under the RCP (Representative Concentration Pathways) 4.5 scenario as inputs, the sub-daily SWAT model predicted a larger amount of monthly maximum daily flow during the wet years than the daily model. The differences between the daily and sub-daily SWAT model simulation results indicated that temporal rainfall resolution could have much impact on the simulation of hydrological process, streamflow, and consequently pollutant transport by SWAT models. There is an imperative need for more studies to examine the effects of temporal rainfall resolution on the simulation of hydrological and water pollutant transport processes by SWAT in river basins of different environmental conditions.     

Evaluating spatio-temporal representations in daily rainfall sequences from three stochastic multi-site weather generation approaches

Many hydrological and agricultural studies require simulations of weather variables reflecting observed spatial and temporal dependence at multiple point locations. This paper assesses three multi-site daily rainfall generators for their ability to model different spatio-temporal rainfall attributes over the study area. The approaches considered consist of a multi-site modified Markov model (MMM), a reordering method for reconstructing space–time variability, and a nonparametric k-nearest neighbour (KNN) model. Our results indicate that all the approaches reproduce adequately the observed spatio-temporal pattern of the multi-site daily rainfall. However, different techniques used to signify longer time scale observed temporal and spatial dependences in the simulated sequences, reproduce these characteristics with varying successes. While each approach comes with its own advantages and disadvantages, the MMM has an overall advantage in offering a mechanism for modelling varying orders of serial dependence at each point location, while still maintaining the observed spatial dependence with sufficient accuracy. The reordering method is simple and intuitive and produces good results. However, it is primarily driven by the reshuffling of the simulated values across realisations and therefore may not be suited in applications where data length is limited or in situations where the simulation process is governed by exogenous conditioning variables. For example, in downscaling studies where KNN and MMM can be used with confidence.     

Effects of temporal scales and space mismatches on the TRMM 3B42 v7 precipitation product in a remote mountainous area

The accurate measurement of precipitation is essential to understanding regional hydrological processes and hydrological cycling. Quantification of precipitation over remote regions such as the Tibetan Plateau is highly unreliable because of the scarcity of rain gauges. The objective of this study is to evaluate the performance of the satellite precipitation product of tropical rainfall measuring mission (TRMM) 3B42 v7 at daily, weekly, monthly, and seasonal scales. Comparison between TRMM grid precipitation and point‐based rain gauge precipitation was conducted using nearest neighbour and bilinear weighted interpolation methods. The results showed that the TRMM product could not capture daily precipitation well due to some rainfall events being missed at short time scales but provided reasonably good precipitation data at weekly, monthly, and seasonal scales. TRMM tended to underestimate the precipitation of small rainfall events (less than 1 mm/day), while it overestimated the precipitation of large rainfall events (greater than 20 mm/day). Consequently, TRMM showed better performance in the summer monsoon season than in the winter season. Through comparison, it was also found that the bilinear weighted interpolation method performs better than the nearest neighbour method in TRMM precipitation extraction.     

Impacts of climate change on the Qingjiang Watershed’s runoff change trend in China

Qingjiang River, the second largest tributary of the Yangtze River in Hubei Province, has taken on the important tasks for power generation and flood control in Hubei Province. The Qingjiang River watershed has a subtropical monsoon climate and, as a result, has dramatic diversity in its water resources. Recently, global warming and climate change have seriously affected the Qingjiang watershed’s integrated water resources management. In this article, general circulation model (GCM) and watershed hydrological models were applied to analyze the impacts of climate change on future runoff of Qingjiang Watershed. To couple the scale difference between GCM and watershed hydrological models, a statistical downscaling method based on the smooth support vector machine was used to downscale the GCM’s large-scale output. With the downscaled precipitation and evaporation, the Xin-anjiang hydrological model and HBV model were applied to predict the future runoff of Qingjiang Watershed under A2 and B2 scenarios. The preformance of the one-way coupling approach in simulating the hydrological impact of climate change in the Qingjiang watershed is evaluated, and the change trend of the future runoff of Qingjiang Watershed under the impacts of climate change is presented and discussed.     

Assessing regression‐based statistical approaches for downscaling precipitation over North America

This paper assesses linear regression‐based methods in downscaling daily precipitation from the general circulation model (GCM) scale to a regional climate model (RCM) scale (45‐ and 15‐km grids) and down to a station scale across North America. Traditional downscaling experiments (linking reanalysis/dynamical model predictors to station precipitation) as well as nontraditional experiments such as predicting dynamic model precipitation from larger‐scale dynamic model predictors or downscaling dynamic model precipitation from predictors at the same scale are conducted. The latter experiments were performed to address predictability limit and scale issues. The results showed that the downscaling of daily precipitation occurrence was rarely successful at all scales, although results did constantly improve with the increased resolution of climate models. The explained variances for downscaled precipitation amounts at the station scales were low, and they became progressively better when using predictors from a higher‐resolution climate model, thus showing a clear advantage in using predictors from RCMs driven by reanalysis at its boundaries, instead of directly using reanalysis data. The low percentage of explained variances resulted in considerable underestimation of daily precipitation mean and standard deviation. Although downscaling GCM precipitation from GCM predictors (or RCM precipitation from RCM predictors) cannot really be considered downscaling, as there is no change in scale, the exercise yields interesting information as to the limit in predictive ability at the station scale. This was especially clear at the GCM scale, where the inability of downscaling GCM precipitation from GCM predictors demonstrates that GCM precipitation‐generating processes are largely at the subgrid scale (especially so for convective events), thus indicating that downscaling precipitation at the station scale from GCM scale is unlikely to be successful. Although results became better at the RCM scale, the results indicate that, overall, regression‐based approaches did not perform well in downscaling precipitation over North America. Copyright © 2013 John Wiley & Sons, Ltd.     

Statistical downscaling of GCM simulations to streamflow using relevance vector machine

General circulation models (GCMs), the climate models often used in assessing the impact of climate change, operate on a coarse scale and thus the simulation results obtained from GCMs are not particularly useful in a comparatively smaller river basin scale hydrology. The article presents a methodology of statistical downscaling based on sparse Bayesian learning and Relevance Vector Machine (RVM) to model streamflow at river basin scale for monsoon period (June, July, August, September) using GCM simulated climatic variables. NCEP/NCAR reanalysis data have been used for training the model to establish a statistical relationship between streamflow and climatic variables. The relationship thus obtained is used to project the future streamflow from GCM simulations. The statistical methodology involves principal component analysis, fuzzy clustering and RVM. Different kernel functions are used for comparison purpose. The model is applied to Mahanadi river basin in India. The results obtained using RVM are compared with those of state-of-the-art Support Vector Machine (SVM) to present the advantages of RVMs over SVMs. A decreasing trend is observed for monsoon streamflow of Mahanadi due to high surface warming in future, with the CCSR/NIES GCM and B2 scenario.     

A spatial temporal downscaling approach to development of IDF relations for Perth airport region in the context of climate change

ABSTRACT

Downscaling of climate projections is the most adapted method to assess the impacts of climate change at regional and local scales. This study utilized both spatial and temporal downscaling approaches to develop intensity–duration–frequency (IDF) relations for sub-daily rainfall extremes in the Perth airport area. A multiple regression-based statistical downscaling model tool was used for spatial downscaling of daily rainfall using general circulation models (GCMs) (Hadley Centre’s GCM and Canadian Global Climate Model) climate variables. A simple scaling regime was identified for 30 minutes to 24 hours duration of observed annual maximum (AM) rainfall. Then, statistical properties of sub-daily AM rainfall were estimated by scaling an invariant model based on the generalized extreme value distribution. RMSE, Nash-Sutcliffe efficiency coefficient and percentage bias values were estimated to check the accuracy of downscaled sub-daily rainfall. This proved the capability of the proposed approach in developing a linkage between large-scale GCM daily variables and extreme sub-daily rainfall events at a given location. Finally IDF curves were developed for future periods, which show similar extreme rainfall decreasing trends for the 2020s, 2050s and 2080s for both GCMs.

Editor M.C. Acreman; Associate editor S. Kanae     

Stormwater quantity and quality response to climate change using artificial neural networks

This research investigates the potential impacts of climate change on stormwater quantity and quality generated by urban residential areas on an event basis in the rainy season. An urban residential stormwater drainage area in southeast Calgary, Alberta, Canada is the focus of future climate projections from general circulation models (GCMs). A regression‐based statistical downscaling tool was employed to conduct spatial downscaling of daily precipitation and daily mean temperature using projection outputs from the coupled GCM. Projected changes in precipitation and temperature were applied to current climate scenarios to generate future climate scenarios. Artificial neural networks (ANNs) developed for modelling stormwater runoff quantity and quality used projected climate scenarios as network inputs. The hydrological response to climate change was investigated through stormwater runoff volume and peak flow, while the water quality responses were investigated through the event mean value (EMV) of five parameters: turbidity, conductivity, water temperature, dissolved oxygen (DO) and pH. First flush (FF) effects were also noted. Under future climate scenarios, the EMVs of turbidity increased in all storms except for three events of short duration. The EMVs of conductivity were found to decline in small and frequent storms (return period < 5 years); but conductivity EMVs were observed to increase in intensive events (return period ≥ 5 years). In general, an increasing EMV was observed for water temperature, whereas a decreasing trend was found for DO EMV. No clear trend was found in the EMV of pH. In addition, projected future climate scenarios do not produce a stronger FF effect on dissolved solids and suspended solids compared to that produced by the current climate scenario. Copyright © 2010 John Wiley & Sons, Ltd.     

Statistically downscaled probabilistic multi‐model ensemble projections of precipitation change in a watershed

This paper presents the development of a probabilistic multi‐model ensemble of statistically downscaled future projections of precipitation of a watershed in New Zealand. Climate change research based on the point estimates of a single model is considered less reliable for decision making, and multiple realizations of a single model or outputs from multiple models are often preferred for such purposes. Similarly, a probabilistic approach is preferable over deterministic point estimates. In the area of statistical downscaling, no single technique is considered a universal solution. This is due to the fact that each of these techniques has some weaknesses, owing to its basic working principles. Moreover, watershed scale precipitation downscaling is quite challenging and is more prone to uncertainty issues than downscaling of other climatological variables. So, multi‐model statistical downscaling studies based on a probabilistic approach are required. In the current paper, results from the three well‐reputed statistical downscaling methods are used to develop a Bayesian weighted multi‐model ensemble. The three members of the downscaling ensemble of this study belong to the following three broad categories of statistical downscaling methods: (1) multiple linear regression, (2) multiple non‐linear regression, and (3) stochastic weather generator. The results obtained in this study show that the new strategy adopted here is promising because of many advantages it offers, e.g. it combines the outputs of multiple statistical downscaling methods, provides probabilistic downscaled climate change projections and enables the quantification of uncertainty in these projections. This will encourage any future attempts for combining the results of multiple statistical downscaling methods. Copyright © 2011 John Wiley & Sons, Ltd.     

Observational and predictive uncertainties for multiple variables in a spatially distributed hydrological model

In this study, uncertainty in model input data (precipitation) and parameters is propagated through a physically based, spatially distributed hydrological model based on the MIKE SHE code. Precipitation uncertainty is accounted for using an ensemble of daily rainfall fields that incorporate four different sources of uncertainty, whereas parameter uncertainty is considered using Latin hypercube sampling. Model predictive uncertainty is assessed for multiple simulated hydrological variables (discharge, groundwater head, evapotranspiration, and soil moisture). Utilizing an extensive set of observational data, effective observational uncertainties for each hydrological variable are assessed. Considering not only model predictive uncertainty but also effective observational uncertainty leads to a notable increase in the number of instances, for which model simulation and observations are in good agreement (e.g., 47% vs. 91% for discharge and 0% vs. 98% for soil moisture). Effective observational uncertainty is in several cases larger than model predictive uncertainty. We conclude that the use of precipitation uncertainty with a realistic spatio‐temporal correlation structure, analyses of multiple variables with different spatial support, and the consideration of observational uncertainty are crucial for adequately evaluating the performance of physically based, spatially distributed hydrological models.