Use of daily precipitation uncertainties in streamflow simulation and forecast

Among other sources of uncertainties in hydrologic modeling, input uncertainty due to a sparse station network was tested. The authors tested impact of uncertainty in daily precipitation on streamflow forecasts. In order to test the impact, a distributed hydrologic model (PRMS, Precipitation Runoff Modeling System) was used in two hydrologically different basins (Animas basin at Durango, Colorado and Alapaha basin at Statenville, Georgia) to generate ensemble streamflows. The uncertainty in model inputs was characterized using ensembles of daily precipitation, which were designed to preserve spatial and temporal correlations in the precipitation observations. Generated ensemble flows in the two test basins clearly showed fundamental differences in the impact of input uncertainty. The flow ensemble showed wider range in Alapaha basin than the Animas basin. The wider range of streamflow ensembles in Alapaha basin was caused by both greater spatial variance in precipitation and shorter time lags between rainfall and runoff in this rainfall dominated basin. This ensemble streamflow generation framework was also applied to demonstrate example forecasts that could improve traditional ESP (Ensemble Streamflow Prediction) method.     

Improving interpolation of daily precipitation for hydrologic modelling: spatial patterns of preferred interpolators

Detailed hydrologic models require high‐resolution spatial and temporal data. This study aims at improving the spatial interpolation of daily precipitation for hydrologic models. Different parameterizations of (1) inverse distance weighted (IDW) interpolation and (2) A local weighted regression (LWR) method in which elevation is the explanatory variable and distance, elevation difference and aspect difference are weighting factors, were tested at a hilly setting in the eastern Mediterranean, using 16 years of daily data. The preferred IDW interpolation was better than the preferred LWR scheme in 27 out of 31 validation gauges (VGs) according to a criteria aimed at minimizing the absolute bias and the mean absolute error (MAE) of estimations. The choice of the IDW exponent was found to be more important than the choice of whether or not to use elevation as explanatory data in most cases. The rank of preferred interpolators in a specific VG was found to be a stable local characteristic if a sufficient number of rainy days are averaged. A spatial pattern of the preferred IDW exponents was revealed. Large exponents (3) were more effective closer to the coast line whereas small exponents (1) were more effective closer to the mountain crest. This spatial variability is consistent with previous studies that showed smaller correlation distances of daily precipitation closer to the Mediterranean coast than at the hills, attributed mainly to relatively warm sea‐surface temperature resulting in more cellular convection coastward. These results suggest that spatially variable, physically based parameterization of the distance weighting function can improve the spatial interpolation of daily precipitation. Copyright © 2009 John Wiley & Sons, Ltd.     

Estimation and spatial interpolation of rainfall intensity distribution from the effective rate of precipitation

Great emphasis is being placed on the use of rainfall intensity data at short time intervals to accurately model the dynamics of modern cropping systems, runoff, erosion and pollutant transport. However, rainfall data are often readily available at more aggregated level of time scale and measurements of rainfall intensity at higher resolution are available only at limited stations. A distribution approach is a good compromise between fine-scale (e.g. sub-daily) models and coarse-scale (e.g. daily) rainfall data, because the use of rainfall intensity distribution could substantially improve hydrological models. In the distribution approach, the cumulative distribution function of rainfall intensity is employed to represent the effect of the within-day temporal variability of rainfall and a disaggregation model (i.e. a model disaggregates time series into sets of higher solution) is used to estimate distribution parameters from the daily average effective precipitation. Scaling problems in hydrologic applications often occur at both space and time dimensions and temporal scaling effects on hydrologic responses may exhibit great spatial variability. Transferring disaggregation model parameter values from one station to an arbitrary position is prone to error, thus a satisfactory alternative is to employ spatial interpolation between stations. This study investigates the spatial interpolation of the probability-based disaggregation model. Rainfall intensity observations are represented as a two-parameter lognormal distribution and methods are developed to estimate distribution parameters from either high-resolution rainfall data or coarse-scale precipitation information such as effective intensity rates. Model parameters are spatially interpolated by kriging to obtain the rainfall intensity distribution when only daily totals are available. The method was applied to 56 pluviometer stations in Western Australia. Two goodness-of-fit statistics were used to evaluate the skill—daily and quantile coefficient of efficiency between simulations and observations. Simulations based on cross-validation show that kriging performed better than other two spatial interpolation approaches (B-splines and thin-plate splines).     

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.     

Addressing uncertainty in reflectivity‐rainfall relations in mountain watersheds during summer convection

The use of precipitation estimates from weather radar reflectivity has become widespread in hydrologic predictions. However, uncertainty remains in the use of the nonlinear reflectivity–rainfall (Z‐R) relation, in particular for mountainous regions where ground validation stations are often lacking, land surface data sets are inaccurate and the spatial variability in many features is high. In this study, we assess the propagation of rainfall errors introduced by different Z‐R relations on distributed hydrologic model performance for four mountain basins in the Colorado Front Range. To do so, we compare spatially integrated and distributed rainfall and runoff metrics at seasonal and event time scales during the warm season when convective storms dominate. Results reveal that the basin simulations are quite sensitive to the uncertainties introduced by the Z‐R relation in terms of streamflow, runoff mechanisms and the water balance components. The propagation of rainfall errors into basin responses follows power law relationships that link streamflow uncertainty to the precipitation errors and streamflow magnitude. Overall, different Z‐R relations preserve the spatial distribution of rainfall relative to a reference case, but not the precipitation magnitude, thus leading to large changes in streamflow amounts and runoff spatial patterns at seasonal and event scales. Furthermore, streamflow errors from the Z‐R relation follow a typical pattern that varies with catchment scale where higher uncertainties exist for intermediate‐sized basins. The relatively high error values introduced by two operational Z‐R relations (WSR‐57 and NEXRAD) in terms of the streamflow response indicate that site‐specific Z‐R relations are desirable in the complex terrain region, particularly in light of other uncertainties in the modelling process, such as model parameter values and initial conditions. Copyright © 2012 John Wiley & Sons, Ltd.     

Evaluation of a water equivalent of snow cover map in a small catchment area using a geostatistical approach

We describe an objective method for evaluating the spatial distribution of water equivalents of the snow cover within a small catchment. Regression analysis is used to quantify the relationship between elevation, presence or absence of forest, and potential direct solar radiation as independent variables and water equivalent as measured at a number of sites. First, this regression relationship is used to interpolate water equivalent data all over the basin area. Then we interpolate the residuals of the regression using a geostatistical approach. Superimposing the results obtained by interpolating the regression relationship and the interpolated residuals eventually yields the water equivalent distribution over the test area. The advantages of the interpolation method used lie in the optimal (effective, unbiased) estimation of the interpolated values as well as in the possibility to quantify the associated estimation variances.     

Regional regression models for hydro‐climate change impact assessment

Hydro‐climatic impacts in water resources systems are typically assessed by forcing a hydrologic model with outputs from general circulation models (GCMs) or regional climate models. The challenges of this approach include maintaining a consistent energy budget between climate and hydrologic models and also properly calibrating and verifying the hydrologic models. Subjective choices of loss, flow routing, snowmelt and evapotranspiration computation methods also increase watershed modelling uncertainty and thus complicate impact assessment. An alternative approach, particularly appealing for ungauged basins or locations where record lengths are short, is to predict selected streamflow quantiles directly from meteorological variable output from climate models using regional regression models that also include physical basin characteristics. In this study, regional regression models are developed for the western Great Lakes states using ordinary least squares and weighted least squares techniques applied to selected Great Lakes watersheds. Model inputs include readily available downscaled GCM outputs from the Coupled Model Intercomparison Project Phase 3. The model results provide insights to potential model weaknesses, including comparatively low runoff predictions from continuous simulation models that estimate potential evapotranspiration using temperature proxy information and comparatively high runoff projections from regression models that do not include temperature as an explanatory variable. Copyright © 2014 John Wiley & Sons, Ltd.     

Long‐term trend of precipitation and runoff in Korean river basins

The spatial and temporal variations of precipitation and runoff for 139 basins in South Korea were investigated for 34 years (1968–2001). The Precipitation‐Runoff Modelling System (PRMS) was selected for the assessment of basin hydrologic response to varying climates and physiology. A non‐parametric Mann–Kendall's test and regression analysis are used to detect trends in annual, seasonal, and monthly precipitation and runoff, while Moran's I is adapted to determine the degree of spatial dependence in runoff trend among the basins. The results indicated that the long‐term trends in annual precipitation and runoff were increased in northern regions and decreased in south‐western regions of the study area during the study period. The non‐parametric Mann–Kendall test showed that spring streamflow was decreasing, while summer streamflow was increasing. April precipitation decreased between 15% and 74% for basins located in south‐western part of the Korean peninsula. June precipitation increased between 18% and 180% for the majority of the basins. Trends in seasonal and monthly streamflow show similar patterns compared to trends in precipitation. Decreases in spring runoff are associated with decreases in spring precipitation which, accompanied by rising temperatures, are responsible for reducing soil moisture. The regional patterns of precipitation and runoff changes show a strong to moderate positive spatial autocorrelation, suggesting that there is a high potential for severe spring drought and summer flooding in some parts of Korea if these trends continue in the future. Copyright © 2008 John Wiley & Sons, Ltd.     

Spatial and temporal variability of temperature,precipitation, and streamflow in upper Sang‐kan basin,China

Temporal streamflow variability in an inland hydrologic station and temporal trends and frequency changes at three weather stations in a semiarid river basin located in Loess Plateau, China, were detected by using linear regression, Mann–Kendall analysis, and wavelet transform methods. Double cumulative curve and ordered clustering were used to identify the hydrological periods of upper Sang‐kan (USK) basin between 1957 and 2012. The results indicate that (1) precipitation in the USK basin over the study period did not show any trend, while the temperature showed a significant increase; (2) streamflow flowing out of the USK basin indicated a significant decrease; (3) two distinct hydrological periods – the ‘natural period’ from 1957 to 1984 and the ‘human impact period’ from 1985 to 2012 – were present; and (4) the contributions of climate change and human activities to reduce the streamflow were 36.9% and 63.1% respectively. The results indicate that human activities may be contributing to a decrease in streamflow in the USK basin. Copyright © 2016 John Wiley & Sons, Ltd.     

Estimation of regional meteorological and hydrological drought characteristics: a case study for Denmark

Information on regional drought characteristics provides critical information for adequate water resource management. This study introduces a method to calculate the probability of a specific area to be affected by a drought of a given severity and demonstrates its potential for calculating both meteorological and hydrological drought characteristics. The method is demonstrated using Denmark as a case study. The calculation procedure was applied to monthly precipitation and streamflow series separately, which were linearly transformed by the Empirical Orthogonal Functions (EOF) method. Denmark was divided into 260 grid-cells of 14×17 km, and the monthly mean and the EOF-weight coefficients were interpolated by kriging. The frequency distributions of the first two (streamflow) or three (precipitation) amplitude functions were then derived. By performing Monte Carlo simulations, amplitude functions corresponding to 1000 years of data were generated. Based on these simulated functions as well as interpolated mean and weight coefficients, long time series of precipitation and streamflow were simulated for each grid-cell. The probability distribution functions of the area covered by a drought and the drought deficit volumes were then derived and combined to produce drought severity-area-frequency curves. These curves allowed an estimation of the probability of an area of a certain extent to have a drought of a given severity, and thereby return periods could be assigned to historical drought events. A comparison of drought characteristics showed that streamflow droughts are less homogeneous over the region, less frequent and last for longer time periods than precipitation droughts.     

Spatial variability and interpolation of daily precipitation amount

Daily precipitation amounts show spatial variation over sub-continential regions. Point measurements, represntative for regions of land, have to be interpolated towards unobserved locations. In this study four days in 1984 were selected to investigate the spatial variability of daily precipitation amount in north-western Europe in relation to the meteorological conditions. Data were interpolated using kriging. Crossvalidation was used to compare interpolated values with measured values. Large differences in the spatial structure of daily precipitation amount are observed as a result of different meteorological conditions. Stratification of the study area into a coast, a mountain and an interior stratum proved to be successful, reducing the Mean Squared Error of Prediction with up to 55%.This article was inadvertently printed in SHH 6(3) 1992 without figures and figure legends. The article is being reprinted in this issue in complete form. The editor apologizes for this error in publication.     

人工神经网络模型预测气候变化对博斯腾湖流域径流影响

温室气体排放量增加造成气候变化,对全球资源环境产生重要影响．本文利用人工神经网络模型建立月降水、气温与径流关系,利用开都河流域降水、气温、径流资料对模型进行训练和验证,通过试算法确定网络模型结构,气温升高和降水量增加对径流影响的敏感程度分析表明,气温升高和降水增加对该区域径流影响较大,且气温升高的影响更为显著,径流增加主要集中在夏季,根据区域气候模型（RCMs）推算的CO2加倍情况下西北地区气候的可能变化,预测位于博斯腾湖流域的开都河大山口站年径流量增加38．6％,其中夏季增加71．8％,冬季增加11．4％。     

A synthetic model for daily streamflow

This paper presents a model for synthesising daily average streamflow data that is suitable for most rivers in Great Britain. The method is based on a linear interpolation of the logorithms of 5-day average flows. The 5-day average flows are produced using N.T. Kottegoda's statistical model (Thesis, Univ. of Birmingham, 1970). The 5-day model preserves the long-term statistical characteristics of the daily data, while the short-term characteristics such as hydrograph shape are imposed by the interpolation method.

A stochastic error term is superimposed on the interpolated daily flows. This term represents the non-deterministic component of the daily time series. The analysis of the observed error terms represents an important part of this paper.

The riverflow in the Severn at Bewdley is used to demonstrate both the analysis of actual data and the generation of synthetic data. The technique is then applied to data from two other rivers with widely differing characteristics to demonstrate the range of the method.     

The Use of Radar and Gauge Measurements to Estimate Areal Precipitation for Several Czech River Basins

Two methods estimating areal precipitation for selected river basins in the Czech Republic are compared. The methods use radar precipitation (the radar-derived precipitation estimate based on column maximum reflectivity) and data from 81 on-line rain gauges routinely provided by the Czech Hydrometeorological Institute. Data from a dense network of climatological rain gauges (the average inter-station distance is approximately 8 km), the measurements of which are not available in real time, are utilized for the verification. The mean areal precipitation, which is used as the ground truth, is obtained by the weighted interpolation of the dense rain gauge network. The accuracy of the methods is evaluated by the root-mean-square-error.The first, pixel-related method merges radar precipitation with rain gauge data to obtain adjusted pixel values. The adjusting procedure combines radar and gauge values in one variable that is interpolated into all radar pixels. The adjusted pixel precipitation is calculated from radar precipitation and from the value of the combined variable. The areal estimates are determined by adding the corresponding pixel values. The second method applies a linear regression model to describe the relationship between the areal precipitation (dependent variable) and its estimates, which are determined from (i) non-adjusted radar precipitation and (ii) on-line rain gauge measurements interpolated into pixels. Classical linear regression, ridge regression and robust regression models are tested.Both the methods decrease the average areal error in comparison with the reference method, which uses the on-line rain gauge data only. The decrease is about 10% and 15% for the pixel-related and regression methods, respectively. When the estimates of the pixel-related method are included as predictors into the regression method then the improvement of accuracy is almost 25%.     

Combining regression and spatial proximity for catchment model regionalization: a comparative study

Abstract

Spatial error regression is employed to regionalize the parameters of a rainfall–runoff model. The approach combines regression on physiographic watershed characteristics with a spatial proximity technique that describes the spatial dependence of model parameters. The methodology is tested for the monthly abcd model at a network of gauges in southeast United States and compared against simpler regression and spatial proximity approaches. Unlike other comparative regionalization studies that only evaluate the skill of regionalized streamflow predictions in ungauged catchments, this study also examines the fit between regionalized parameters and their optimal (i.e. calibrated) values. Interestingly, the spatial error model produces parameter estimates that better resemble the optimal parameters than either of the simpler methods, but the spatial proximity method still yields better hydrologic simulations. The analysis suggests that the superior streamflow predictions of spatial proximity result from its ability to better preserve correlations between compensatory hydrological parameters.

Editor D. Koutsoyiannis; Associate editor Y. Gyasi-Agyei     

Accuracy assessment of optimum interpolation in multichannel precipitation analysis

The accuracy of an optimum interpolation technique in filling missing values in multichannel (or multisite) hydrologic series containing time-coincident data gaps is examined. The applied methodology is based on the maximum entropy method (MEM) of spectral estimation or multivariate autoregressive modeling and heavily depends upon the properties of multichannel prediction error filter (PEF). Six precipitation time series spatially located within a hydrologic basin are used and time-coincident artificial gaps are created in all six series. The performance of the technique is assessed by comparing the filled-in series to the observed and by employing spectral analysis. The results reveal the usefulness of the method in multichannel hydrologic analysis.     

Accuracy assessment of optimum interpolation in multichannel precipitation analysis

The accuracy of an optimum interpolation technique in filling missing values in multichannel (or multisite) hydrologic series containing time-coincident data gaps is examined. The applied methodology is based on the maximum entropy method (MEM) of spectral estimation or multivariate autoregressive modeling and heavily depends upon the properties of multichannel prediction error filter (PEF). Six precipitation time series spatially located within a hydrologic basin are used and time-coincident artificial gaps are created in all six series. The performance of the technique is assessed by comparing the filled-in series to the observed and by employing spectral analysis. The results reveal the usefulness of the method in multichannel hydrologic analysis.     

Estimation of rainfall from infrared‐microwave satellite data for basin‐scale hydrologic modelling

The infrared‐microwave rainfall algorithm (IMRA) was developed for retrieving spatial rainfall from infrared (IR) brightness temperatures (TBs) of satellite sensors to provide supplementary information to the rainfall field, and to decrease the traditional dependency on limited rain gauge data that are point measurements. In IMRA, a SLOPE technique (ST) was developed for discriminating rain/no‐rain pixels through IR image cloud‐top temperature gradient, and 243K as the IR threshold temperature for minimum detectable rainfall rate. IMRA also allows for the adjustment of rainfall derived from IR‐TB using microwave (MW) TBs. In this study, IMRA rainfall estimates were assessed on hourly and daily basis for different spatial scales (4, 12, 20, and 100 km) using NCEP stage IV gauge‐adjusted radar rainfall data, and daily rain gauge data. IMRA was assessed in terms of the accuracy of the rainfall estimates and the basin streamflow simulated by the hydrologic model, Sacramento soil moisture accounting (SAC‐SMA), driven by the rainfall data. The results show that the ST option of IMRA gave accurate satellite rainfall estimates for both light and heavy rainfall systems while the Hessian technique only gave accurate estimates for the convective systems. At daily time step, there was no improvement in IR‐satellite rainfall estimates adjusted with MW TBs. The basin‐scale streamflow simulated by SAC‐SMA driven by satellite rainfall data was marginally better than when SAC‐SMA was driven by rain gauge data, and was similar to the case using radar data, reflecting the potential applications of satellite rainfall in basin‐scale hydrologic modelling. Copyright © 2010 John Wiley & Sons, Ltd.     

Mediating stream baseflow response to climate change: The role of basin storage

Inter‐basin differences in streamflow response to changes in regional hydroclimatology may reflect variations in storage characteristics that control the retention and release of water inputs. These aspects of storage could mediate a basin's sensitivity to climate change. The hypothesis that temporal trends in stream baseflow exhibit a more muted reaction to changes in precipitation and evapotranspiration for basins with greater storage was tested on the Oak Ridges Moraine (ORM) in Southern Ontario, Canada. Long‐term (>25 years) baseflow trends for 16 basins were compared to corresponding trends in precipitation amount and type and in potential evapotranspiration as well as shorter trends in groundwater levels for monitoring wells on the ORM. Inter‐basin differences in storage properties were characterized using physiographic, hydrogeologic, land use/land cover, and streamflow metrics. The latter included the slope of the basin's flow duration curve and basin dynamic storage. Most basins showed temporal increases in baseflow, consistent with limited evidence of increases and decreases in regional precipitation and snowfall: precipitation ratio, respectively, and recent increases in groundwater recharge along the crest of the ORM. Baseflow trend magnitude was uncorrelated to basin physiographic, hydrogeologic, land use/land cover, or flow duration curve characteristics. However, it was positively related to a basin's dynamic storage, particularly for basins with limited coverage of open water and wetlands. The dynamic storage approach assumes that a basin behaves as a first‐order dynamical system, and extensive open water and wetland areas in a basin may invalidate this assumption. Previous work suggested that smaller dynamic storage was linked to greater damping of temporal variations in water inputs and reduced interannual variability in streamflow regime. Storage and release of water inputs to a basin may assist in mediating baseflow response to temporal changes in regional hydroclimatology and may partly account for inter‐basin differences in that response. Such storage characteristics should be considered when forecasting the impacts of climate change on regional streamflow.     

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.