Improving global rainfall forecasting with a weather type approach in Japan

An automated version of the weather type classification scheme was performed over Japan to characterize daily circulation conditions. A daily gridded field of mean sea-level pressure (MSLP) from the European Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis dataset (ERA-interim) and the THORPEX Interactive Grand Global Ensemble (TIGGE) daily forecast dataset were used. The weather type is advantageous as it provides an opportunity to improve global rainfall prediction by refining statistical bias correction. We distinguished 11 weather types: anticyclone, cyclone, hybrid and eight purely wind directions. The results indicate that the main weather types contributing to the total volume of rainfall are cyclone, hybrid, purely westerly and northwest winds. A gamma-based bias correction decreases the global rainfall forecast root mean square by 10%, while specific weather type gamma bias correction accounts for 5–10% root mean square error reduction, with a total decrease of errors up to a maximum of 20%. Both global and weather type bias corrections improve the extreme dependency scores (EDS), but for different extreme rainfall thresholds. The study advocates the use of weather type bias-correction methods for extreme event rainfall intensity corrections higher than 100 mm/d.

EDITOR

A. Castellarin

ASSOCIATE EDITOR

A. Jain     

Integration of gauge and radar rainfall to enable best simulation of hydrological parameters

This study is about use of spatially distributed rain in physically based hydrological models. In recent years, spatially distributed radar rainfall data have become available. The distributed radar rain is used to precisely model hydrologic processes and it is more realistic than the past practice of distribution methods like Thiessen polygons. Radar provides a highly accurate spatial distribution of rainfall and greatly improves the basin average rainfall estimates. However, quantification of the exact amount of rainfall from radar observation is relatively difficult. Thus, the fundamental idea of this study is to apply hourly gauge and radar rainfall data in a distributed hydrological model to simulate hydrological parameters. Hence the comparison is made between the outcomes of the WetSpa model from radar rainfall distribution and gauge rainfall distributed by the Thiessen polygon technique. The comparative plots of the hydrograph and the results of hydrological components such as evapotranspiration, surface runoff, soil moisture, recharge and interflow, reflect the spatially distributed radar input performing well for model outflow simulation.

EDITOR D. Koutsoyiannis

ASSOCIATE EDITOR F. Pappenberger     

Assessment of radar-based locally varying anisotropy on daily rainfall interpolation

ABSTRACT

Spatial variability of rainfall has been recognised as an important factor controlling the hydrological response of catchments. However, gauged daily rainfall data are often available at scattered locations over the catchments. This paper looks into how to capitalise on the spatial structure of radar rainfall data for improving kriging interpolation of limited gauge data over catchments at the 1-km² grid scale, using for the case study 117 gauged stations within the 128 km × 128 km region of the Mt Stapylton weather radar field (near Brisbane, Australia). Correlograms were developed using a Fast Fourier Transform method on the Gaussianised radar and gauged data. It is observed that the correlograms vary from day to day and display significant anisotropy. For the radar data, locally varying anisotropy (LVA) was examined by developing the correlogram centred on each pixel and for different radial distances. Cross-validation was carried out using the empirical correlogram tables, as well as different fitting strategies of a two-dimensional exponential distribution for both the gauged and the radar data. The results indicate that the correlograms based on the radar data outperform the gauged ones as judged by statistical measures including root mean square error, mean bias, mean absolute bias, mean standard deviation and mean inter-quartile range. While the radar data display significant LVA, it was observed that LVA did not significantly improve the estimates compared with the global anisotropy. This was also confirmed by conditional simulation of 120 rainfields using different options of correlogram development.

EDITOR M.C. Acreman; ASSOCIATE EDITOR Q. Zhang     

Modelling radar-rainfall estimation uncertainties using elliptical and Archimedean copulas with different marginal distributions

Abstract

Given that radar-based rainfall has been broadly applied in hydrological studies, quantitative modelling of its uncertainty is critically important, as the error of input rainfall is the main source of error in hydrological modelling. Using an ensemble of rainfall estimates is an elegant solution to characterize the uncertainty of radar-based rainfall and its spatial and temporal variability. This paper has fully formulated an ensemble generator for radar precipitation estimation based on the copula method. Each ensemble member is a probable realization that represents the unknown true rainfall field based on the distribution of radar rainfall (RR) error and its spatial error structure. An uncertainty model consisting of a deterministic component and a random error factor is presented based on the distribution of gauge rainfall conditioned on the radar rainfall (GR|RR). Two kinds of copulas (elliptical and Archimedean copulas) are introduced to generate random errors, which are imposed by the deterministic component. The elliptical copulas (e.g. Gaussian and t-copula) generate the random errors based on the multivariate distribution, typically of decomposition of the error correlation matrix using the LU decomposition algorithm. The Archimedean copulas (e.g. Clayton and Gumbel) utilize the conditional dependence between different radar pixels to obtain random errors. Based on those, a case application is carried out in the Brue catchment located in southwest England. The results show that the simulated uncertainty bands of rainfall encompass most of the reference raingauge measurements with good agreement between the simulated and observed spatial dependences. This indicates that the proposed scheme is a statistically reliable method in ensemble radar rainfall generation and is a useful tool for describing radar rainfall uncertainty.

Editor D. Koutsoyiannis; Associate editor S. Grimaldi     

Transfer function noise modelling of groundwater level fluctuation using threshold rainfall-based binary-weighted parameter estimation approach

Considerable uncertainty occurs in the parameter estimates of traditional rainfall–water level transfer function noise (TFN) models, especially with the models built using monthly time step datasets. This is due to the equal weights assigned for rainfall occurring during both water level rise and water level drop events while estimating the TFN model parameters using the least square technique. As an alternative to this approach, a threshold rainfall-based binary-weighted least square method was adopted to estimate the TFN model parameters. The efficacy of this binary-weighted approach in estimating the TFN model parameters was tested on 26 observation wells distributed across the Adyar River basin in Southern India. Model performance indices such as mean absolute error and coefficient of determination values showed that the proposed binary-weighted approach of fitting independent threshold-based TFN models for water level rise and water level drop scenarios considerably improves the model accuracy over other traditional TFN models.

EDITOR D. Koutsoyiannis

ASSOCIATE EDITOR A. Fiori     

Spatial disaggregation and intensity correction of TRMM-based rainfall time series for hydrological applications in dryland catchments

Abstract

A novel approach is presented for combining spatial and temporal detail from newly available TRMM-based data sets to derive hourly rainfall intensities at 1-km spatial resolution for hydrological modelling applications. Time series of rainfall intensities derived from 3-hourly 0.25° TRMM 3B42 data are merged with a 1-km gridded rainfall climatology based on TRMM 2B31 data to account for the sub-grid spatial distribution of rainfall intensities within coarse-scale 0.25° grid cells. The method is implemented for two dryland catchments in Tunisia and Senegal, and validated against gauge data. The outcomes of the validation show that the spatially disaggregated and intensity corrected TRMM time series more closely approximate ground-based measurements than non-corrected data. The method introduced here enables the generation of rainfall intensity time series with realistic temporal and spatial detail for dynamic modelling of runoff and infiltration processes that are especially important to water resource management in arid regions.

Editor D. Koutsoyiannis

Citation Tarnavsky, E., Mulligan, M. and Husak, G., 2012. Spatial disaggregation and intensity correction of TRMM-based rainfall time series for hydrological applications in dryland catchments. Hydrological Sciences Journal, 57 (2), 248–264.     

Simultaneous hydrological prediction at multiple gauging stations using the NARX network for Kemaman catchment,Terengganu, Malaysia

ABSTRACT

This paper presents a neural network model capable of catchment-wide simultaneous prediction of river stages at multiple gauging stations. Thirteen meteorological parameters are considered in the input, which includes rainfall, temperature, mean relative humidity and evaporation. The NARX model is trained with a representative set of hourly data, with optimal time delay for both the input and output. The network trained using 120-day data is able to produce simulations that are in excellent agreement with field observations. We show that for application with one-step-ahead predictions, the loss in network performance is marginal. Inclusion of additional tidal observations does not improve predictions, suggesting that the river stage stations under consideration are not sensitive to tidal backwater effects despite the claim commonly made.

EDITOR D. Koutsoyiannis ASSOCIATE EDITOR F. Pappenberger     

Hydrological system analysis and modelling of the Kara River basin (West Africa) using a lumped metric conceptual model

This paper discusses the analysis and modelling of the hydrological system of the basin of the Kara River, a transboundary river in Togo and Benin, as a necessary step towards sustainable water resources management. The methodological approach integrates the use of discharge parameters, flow duration curves and the lumped conceptual model IHACRES. A Sobol sensitivity analysis is performed and the model is calibrated by applying the shuffled complex evolution algorithm. Results show that discharge generation in three nested catchments of the basin is affected by landscape physical characteristics. The IHACRES model adequately simulates the rainfall–runoff dynamics in the basin with a mean modified Nash-Sutcliffe efficiency measure of 0.6. Modelling results indicate that parameters controlling rainfall transformation to effective rainfall are more sensitive than those routing the streamflow. This study provides insights into understanding the catchment’s hydrological system. Nevertheless, further investigations are required to better understand detailed runoff generation processes.

EDITOR M.C. Acreman; ASSOCIATE EDITOR N Verhoest     

Evaluation of the climate prediction center (CPC) morphing technique (CMORPH) rainfall product on hourly time scales over the source of the Blue Nile River

Limited availability of surface‐based rainfall observations constrains the evaluation of satellite rainfall products over many regions. Observations are also often not available at time scales to allow evaluation of satellite products at their finest resolutions. In the present study, we utilized a 3‐month rainfall data set from an experimental network of eight automatic gauges in Gilgel Abbay watershed in Ethiopia to evaluate the 1‐hourly, 8 × 8‐km Climate Prediction Center morphing technique (CMORPH) rainfall product. The watershed is situated in the Lake Tana basin which is the source of the Blue Nile River. We applied a suite of statistical metrics that included mean difference, bias, standard deviation of differences and measures of association. Our results indicate that the accuracy of the CMORPH product shows a significant variation across the basin area. Its estimates are mostly within ±10 mm h^?1 of the gauge rainfall observations; however, the product does not satisfactorily capture the rainfall temporal variability and is poorly correlated (<0.27) to gauge observations. Its poor rain detection capability led to significant underestimation of the seasonal rainfall depth (total bias reaches up to ?52%) with large amounts of hit rain bias as well as missed rain and false rain biases. In the future refinement of CMORPH algorithm, more attention should be given to reducing missed rain bias over the mountains of Gilgel Abbay, whereas equal attention should be given to hit, missed rain and false rain biases over other parts of the watershed. Copyright © 2012 John Wiley & Sons, Ltd.     

On the potential of Bartlett Lewis rectangular pulse models for simulating rainfall in Peninsular Malaysia

Abstract

The applicability of two versions of the Bartlett Lewis rectangular pulse model, the original and the modified model, is discussed for describing the temporal and spatial variation of rainfall patterns observed at 15 raingauge stations in Peninsular Malaysia over the period 1971–2008; 17 different sets of moment combinations are fitted to these models based on the generalized method of moments approach. The common statistics included in all sets are the mean, variance, lag-1 autocorrelation and the probability of dry based on the hourly rainfall data. The analysis was carried out on hourly rainfall data from all 15 stations for all months of the year. Two stations, Petaling Jaya and Kemaman, located on the west and east coasts of the Peninsula, respectively, are considered for illustration of the results, taking the months of July and November, which correspond to the driest and wettest months, corresponding to the southwest monsoon (May–August) and northeast monsoon (November–February), respectively. The best moment combination found for the illustrative results is based on the common statistics, as well as the mean and variance based on 24-h aggregated rainfall data, the inclusion of which successfully improved the model performance; the errors were significantly reduced. It was also found that the performance of the fitted models based on the mean absolute deviate error varies according to the type of Bartlett Lewis model applied: errors are much smaller for the fitted model based on the modified model as compared to the original model. In addition, the fitted statistics: mean, lag-1 autocorrelation and probability of dry are quite well fitted for several aggregated time scales; however, the variances are underestimated in both models for all aggregated time scales, particularly in the case of the original model. The results of extreme value analysis indicate that the modified model failed to reproduce the annual hourly and daily rainfall extremes satisfactorily.

Editor D. Koutsoyiannis; Associate editor C. Onof

Citation Hanaish, I.S., Ibrahim, K., and Jemain, A.A., 2013. On the potential of Bartlett Lewis rectangular pulse models for simulating rainfall in Peninsular Malaysia. Hydrological Sciences Journal, 58 (8), 1690–1703.     

Correcting the TRMM rainfall product for hydrological modelling in sparsely-gauged mountainous basins

This study developed a correction approach to improve the rainfall field estimation using the TRMM rainfall product in a sparsely-gauged mountainous basin. First, Thiessen polygons were generated to define the measurement domain of each raingauge. Second, the rainfall of TRMM pixels in each Thiessen polygon was corrected using a benchmark method based on the difference between the monthly rainfall estimated by a raingauge and the TRMM pixel that possessed a gauge station (referred to as a gauged pixel). Third, the rainfall values in the gauged pixels were adjusted to the weighted average value of the gauge rainfall and corrected pixel rainfall. Finally, the rainfall in the non-gauged TRMM pixels was corrected as the sum of two terms. The first term is the adjusted rainfall in the corresponding gauged pixel in the same Thiessen polygon, and the second term is the rainfall (after benchmark correction) difference between the current pixel and the gauged pixel. Our results indicate that the corrected rainfall data outperforms the original TRMM product in the simulations of moderate and low flows and outperforms the sparse raingauges in the simulations of both peak and low flows.

EDITOR A. Castellarin; ASSOCIATE EDITOR S. Huang     

A new magnitude category disaggregation approach for temporal high‐resolution rainfall intensities

Simulation of quick runoff components such as surface runoff and associated soil erosion requires temporal high‐resolution rainfall intensities. However, these data are often not available because such measurements are costly and time consuming. Current rainfall disaggregation methods have shortcomings, especially in generating the distribution of storm events. The objectives of this study were to improve point rainfall disaggregation using a new magnitude category rainfall disaggregation approach. The procedure is introduced using a coupled disaggregation approach (Hyetos and cascade) for multisite rainfall disaggregation. The new procedure was tested with ten long‐term precipitation data sets of central Germany using summer and winter precipitation to determine seasonal variability. Results showed that dividing the rainfall amount into four daily rainfall magnitude categories (1–10, 11–25, 26–50, >50 mm) improves the simulation of high rainfall intensity (convective rainfall). The Hyetos model category approach (Hyetos_Cat) with seasonal variation performs representative to observed hourly rainfall compared with without categories on each month. The mean absolute percentage accuracy of standard deviation for hourly rainfall is 89.7% in winter and 95.6% in summer. The proposed magnitude category method applied with the coupled Hyetos_Cat–cascade approach reproduces successfully the statistical behaviour of local 10‐min rainfall intensities in terms of intermittency as well as variability. The root mean square error performance statistics for disaggregated 10‐min rainfall depth ranges from 0.20 to 2.38 mm for summer and from 0.12 to 2.82 mm for the winter season in all categories. The coupled stochastic approach preserves the statistical self‐similarity and intermittency at each magnitude category with a relatively low computational burden. Copyright © 2014 John Wiley & Sons, Ltd.     

Sensitivity of hydrological models to uncertainty in rainfall input

Abstract

Different approaches used in hydrological modelling are compared in terms of the way each one takes the rainfall data into account. We examine the errors associated with accounting for rainfall variability, whether in hydrological modelling (distributed vs lumped models) or in computing catchment rainfall, as well as the impact of each approach on the representativeness of the parameters it uses. The database consists of 1859 rainfall events, distributed on 500 basins, located in the southeast of France with areas ranging from 6.2 to 2851 km². The study uses as reference the hydrographs computed by a distributed hydrological model from radar rainfall. This allows us to compare and to test the effects of various simplifications to the process when taking rainfall information (complete rain field vs sampled rainfall) and rainfall–runoff modelling (lumped vs distributed) into account. The results appear to show that, in general, the sampling effect can lead to errors in discharge at the outlet that are as great as, or even greater than, those one would get with a fully lumped approach. We found that small catchments are more sensitive to the uncertainties in catchment rainfall input generated by sampling rainfall data as seen through a raingauge network. Conversely, the larger catchments are more sensitive to uncertainties generated when the spatial variability of rainfall events is not taken into account. These uncertainties can be compensated for relatively easily by recalibrating the parameters of the hydrological model, although such recalibrations cause the parameter in question to completely lose physical meaning.

Citation Arnaud, P., Lavabre, J., Fouchier, C., Diss, S. & Javelle, P. (2011) Sensitivity of hydrological models to uncertainty of rainfall input. Hydrol. Sci. J. 56(3), 397–410.     

Regionalisation of hydrological responses under land-use change and variable data quality

ABSTRACT

Estimating river flows at ungauged sites is generally recognised as an important area of research. In countries or regions with rapid land development and sparse hydrological gauging networks, three particular challenges may arise—data scarcity, data quality, and hydrological non-stationarity. Using data from 44 gauged sub-catchments of the upper Ping catchment in northern Thailand from the period 1995–2006, three relevant flow response indices (runoff coefficient, base flow index and seasonal elasticity of flow) were regionalised by regression against available catchment properties. The runoff coefficient was the most successfully regionalised, followed by base flow index and lastly the seasonal elasticity. The non-stationarity (represented by the differences between two 6-year sub-periods) was significant both in the flow response indices and in land use indices; however relationships between the two sets of indices were weak. The regression equations derived from regionalisation were not helpful in predicting the non-stationarity in the flow indices except somewhat for the runoff coefficient. A partly subjective data quality scoring system was devised, and showed the clear influence of rainfall and flow data quality on regionalisation uncertainty. Recommendations towards improving data support for hydrological regionalisation in Thailand include more relevant soils databases, improved records of abstractions and investment in the gauge network. Widening of the regionalisation beyond the upper Ping and renewed efforts at using remotely sensed rainfall data are other possible ways forward.

EDITOR Z.W. Kundzewicz ASSOCIATE EDITOR T. Wagener     

Real-time correction of water stage forecast during rainstorm events using combination of forecast errors

This study proposes a real-time error correction method for the forecasted water stage using a combination of forecast errors estimated by the time series models, AR(1), AR(2), MA(1) and MA(2), and the average deviation model to update the water stage forecast during rainstorm events. During flood forecasting and warning operations, the proposed real-time error correction method takes advantage of being individually and continuously implemented and the results not being updated to the hydrological model and hydraulic routings so as to save computational time by recalibrating the parameters of the proposed methods with real-time observation. For model validation, the current study adopts the observed and forecasted data on a severe typhoon, Morakot, collected at eight water level gauges in Southern Taiwan and provided by the flood forecast system FEWS_Taiwan, which is linked with the reliable quantitative precipitation forecast (QPF) at 3 h of lead time provided by the Center Weather Bureau in Taiwan, as the model validation. The results of numerical experiments indicate that the proposed real-time error correction method can effectively reduce the errors of forecasted water stages at the 1-, 2-, and 3-h lead time and so enhance the reliability of forecast information issued by the FEWS_Taiwan. By means of real-time estimating potential forecast error, the uncertainties in hydrology, modules as well as associated parameters, and physiographical features of the river can be reduced.     

Investigation of the groundwater modelling component of the Integrated Water Flow Model (IWFM)

ABSTRACT

The Integrated Water Flow Model (IWFM), developed by the California Department of Water Resources, is an integrated hydrological model that simulates key flow processes including groundwater flows, streamflow, stream–aquifer interactions, rainfall–runoff and infiltration. It also simulates the agricultural water demand as a function of soil, crop and climatic characteristics, as well as irrigation practices, and allows the user to meet these demands through pumping and stream diversions. This study investigates the modelling performance of the groundwater module of IWFM using several hypothetical test problems that cover a wide range of settings and boundary conditions, by comparing the simulation results with analytical solutions, field and laboratory observations, or with results from MODFLOW outputs. The comparisons demonstrate that IWFM is capable of simulating various hydrological processes reliably.

EDITOR M.C. Acreman; ASSOCIATE EDITOR A. Efstratiadis     

Prediction of climate change effects on the runoff regime of a forested catchment in northern Iran

ABSTRACT

The southern coast of the Caspian Sea in northern Iran is bordered by a mountain range with forested catchments which are susceptible to droughts and floods. This paper examines possible changes to runoff patterns from one of these catchments in response to climate change scenarios. The HEC-HMS rainfall–runoff model was used with downscaled future rainfall and temperature data from 13 global circulation models, and meteorological and hydrometrical data from the Casilian (or “Kassilian”) Catchment. Annual and seasonal predictions of runoff change for three future emissions scenarios were obtained, which suggest significantly higher spring rainfall with increased risk of flooding and significantly lower summer rainfall leading to a higher probability of drought. Flash floods arising from extreme rainfall may become more frequent, occurring at any time of year. These findings indicate a need for strategic planning of water resource management and mitigation measures for increasing flood hazards.

EDITOR M.C. Acreman ASSOCIATE EDITOR not assigned     

Assessment of hydrological extremes in the Kamo River Basin,Japan

A suite of extreme indices derived from daily precipitation and streamflow was analysed to assess changes in the hydrological extremes from 1951 to 2012 in the Kamo River Basin. The evaluated indices included annual maximum 1-day and 5-day precipitation (RX1day, RX5day), consecutive dry days (CDD), annual maximum 1-day and 5-day streamflow (SX1day, SX5day), and consecutive low-flow days (CDS). Sen’s slope estimator and two versions of the Mann-Kendall test were used to detect trends in the indices. Also, frequency distributions of the indices were analysed separately for two periods: 1951–1981 and 1982–2012. The results indicate that quantiles of the rainfall indices corresponding to the 100-year return period have decreased in recent years, and the streamflow indices had similar patterns. Although consecutive no rainfall days represented by 100-year CDD decreased, continuous low-flow days represented by 100-year CDS increased. This pattern change is likely associated with the increase in temperature during this period.

EDITOR D. Koutsoyiannis

ASSOCIATE EDITOR E. Gargouri     

基于面平均雨量误差修正的实时洪水预报修正方法

空间集总式水文模型的洪水预报精度会受到面平均雨量估计误差的严重影响.点雨量监测值的误差类型、误差大小以及流域的雨量站点密度和站点的空间分布都会影响到面平均雨量的计算.为提高实时洪水预报精度,本文提出了一种基于降雨系统响应曲线洪水预报误差修正方法.通过此方法估计降雨输入项的误差,从而提高洪水预报精度.此方法将水文模型做为输入和输出之间的响应系统,用实测流量和计算流量之间的差值做为信息,通过降雨系统响应曲线,使用最小二乘估计原理,对面平均雨量进行修正,再用修正后的面平均雨量重新计算出流过程.将此修正方法结合新安江模型使用理想案例进行检验,并应用于王家坝流域的16场历史洪水以及此流域不同雨量站密度的情况下,结果证明均有明显修正效果,且在雨量站密度较低时修正效果更加明显.该方法是一种结构简单且不增加模型参数和复杂度的实时洪水修正的新方法.