Geostatistical interpolation of space–time rainfall on Tamshui River basin,Taiwan

Taiwan suffers from heavy storm rainfall during the typhoon season. This usually causes large river runoff, overland flow, erosion, landslides, debris flows, loss of power, etc. In order to evaluate storm impacts on the downstream basin, a real‐time hydrological modelling is used to estimate potential hazard areas. This can be used as a decision‐support system for the Emergency Response Center, National Fire Agency Ministry, to make ‘real‐time’ responses and minimize possible damage to human life and property. This study used 34 observed events from 14 telemetered rain‐gauges in the Tamshui River basin, Taiwan, to study the spatial–temporal characteristics of typhoon rainfall. In the study, regionalized theory and cross‐semi‐variograms were used to identify the spatial‐temporal structure of typhoon rainfall. The power form and parameters of the cross‐semi‐variogram were derived through analysis of the observed data. In the end, cross‐validation was used to evaluate the performance of the interpolated rainfall on the river basin. The results show the derived rainfall interpolator represents the observed events well, which indicates the rainfall interpolator can be used as a spatial‐temporal rainfall input for real‐time hydrological modelling. Copyright © 2007 John Wiley & Sons, Ltd.     

Continuous daily hydrograph simulation using duration curves of a precipitation index

The paper describes a parsimonious approach for generating continuous daily stream‐flow time‐series from observed daily rainfall data in a catchment. The key characteristic in the method is a duration curve. It is used to convert the daily rainfall information from source rain gauges into a continuous daily hydrograph at the destination river site. For each source rain gauge a time‐series of rainfall related ‘current precipitation index’ is generated and its duration curve is established. The current precipitation index reflects the current catchment wetness and is defined as a continuous function of precipitation, which accumulates on rainy days and exponentially decays during the periods of no rainfall. The process of rainfall‐to‐runoff conversion is based on the assumption that daily current precipitation index values at rainfall site(s) in a catchment and the destination site's daily flows correspond to similar probabilities on their respective duration curves. The method is tested in several small catchments in South Africa. The method is designed primarily for application at ungauged sites in data‐poor regions where the use of more complex and information consuming techniques of data generation may not be justified. Copyright © 2000 John Wiley & Sons, Ltd.     

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.     

Cokriging for enhanced spatial interpolation of rainfall in two Australian catchments

Rainfall data in continuous space provide an essential input for most hydrological and water resources planning studies. Spatial distribution of rainfall is usually estimated using ground‐based point rainfall data from sparsely positioned rain‐gauge stations in a rain‐gauge network. Kriging has become a widely used interpolation method to estimate the spatial distribution of climate variables including rainfall. The objective of this study is to evaluate three geostatistical (ordinary kriging [OK], ordinary cokriging [OCK], kriging with an external drift [KED]), and two deterministic (inverse distance weighting, radial basis function) interpolation methods for enhanced spatial interpolation of monthly rainfall in the Middle Yarra River catchment and the Ovens River catchment in Victoria, Australia. Historical rainfall records from existing rain‐gauge stations of the catchments during 1980–2012 period are used for the analysis. A digital elevation model of each catchment is used as the supplementary information in addition to rainfall for the OCK and kriging with an external drift methods. The prediction performance of the adopted interpolation methods is assessed through cross‐validation. Results indicate that the geostatistical methods outperform the deterministic methods for spatial interpolation of rainfall. Results also indicate that among the geostatistical methods, the OCK method is found to be the best interpolator for estimating spatial rainfall distribution in both the catchments with the lowest prediction error between the observed and estimated monthly rainfall. Thus, this study demonstrates that the use of elevation as an auxiliary variable in addition to rainfall data in the geostatistical framework can significantly enhance the estimation of rainfall over a catchment.     

Calibration of an oscillating nozzle‐type rainfall simulator

Nozzle‐type rainfall simulators are commonly used in hydrologic and soil erosion research. Simulated rainfall intensity, originating from the nozzle, increases as the distance between the point of measurement and the source is decreased. Hence, rainfall measured using rain gauges would systematically overestimate the rainfall received at the ground level. A simple model was developed to adjust rainfall measured anywhere under the simulator to plot‐wide average rainfall at the ground level. Nozzle height, plot width, gauge diameter and height, and gauge location are required to compute this adjustment factor. Results from 15 runs at different rain intensities and durations, and with different rain gauge layouts, showed that a simple average of measured rain would overestimate the plot‐wide rain by about 20 per cent. Using the adjustment factor to convert measured rainfall for individual gauges before averaging improved the estimate of plot‐wide rainfall considerably. For the 15 runs considered, overall discrepancy between actual and measured rain is reduced to less than 1 per cent with a standard error of 0·97 mm. This model can be easily tested in the ?eld by comparing rainfall depths of different sized gauges. With the adjustment factor they should all give very similar values. Copyright © 2003 John Wiley & Sons, Ltd.     

Occurrence,shape and movement of daily rainfall in the Huai Basin,China

The upper reaches of the Huai River in Central China are located in the East Asian monsoon region. Strong seasonality, as well as large interannual variability of rainfall, causes floods and an uneven supply of water. In order to conserve the water and mitigate the floods, dams and flood protection structures are constructed. Their design requires information about the rainfall. Daily observations from 1957 to 1986 from 78 rain gauges were used to study shape, orientation, movement and geographical and seasonal occurrence of storms in the 79 000 km² study area. The rainfall characteristics were described using graphical plots, cross‐ and autocorrelation. Storms larger than 50 mm/day were found to occur from February to November, whereas storms exceeding 350 mm/day were confined to the main rainfall season from late June to mid‐August. The southern part of the study area experienced a break in the rainfall season in late July, corresponding to the seasonal northward shift of the rain belt. A weekly periodicity of 7–8 days for rainfall was found during June–July, but not during August–September. During the whole period June–September, the spatial pattern of daily rainfall revealed an elongated shape, more pronounced during June–July than August–September. The rainfall area was orientated approximately from WSW to ENE during the whole period, and showed an anticlockwise rotation of about 16° per day during June–July. The cross‐correlation analysis revealed that the rainfall area moved about 100 km/day eastward. These results and an investigation of meteorological maps indicate that the spatial correlation pattern of daily rainfall is produced by cold fronts on the Mei‐Yu front. Suggestions are made as to how to use the results for the construction of design rainfalls in the study area. Copyright © 1999 John Wiley & Sons, Ltd.     

Seasonal snowpack dynamics and runoff in a cool temperate forest: lysimeter experiment in Niigata,Japan

Seasonal snowpack dynamics are described through field measurements under contrasting canopy conditions for a mountainous catchment in the Japan Sea region. Microclimatic data, snow accumulation, albedo and lysimeter runoff are given through the complete winter season 2002–03 in (1) a mature cedar stand, (2) a larch stand, and (3) a regenerating cedar stand or opening. The accumulation and melt of seasonal snowpack strongly influences streamflow runoff during December to May, including winter baseflow, mid‐winter melt, rain on snow, and diurnal peaks driven by radiation melt in spring. Lysimeter runoff at all sites is characterized by constant ground melt of 0·8–1·0 mm day⁻¹. Rapid response to mid‐winter melt or rainfall shows that the snowpack remains in a ripe or near‐ripe condition throughout the snow‐cover season. Hourly and daily lysimeter discharge was greatest during rain on snow (e.g. 7 mm h⁻¹ and 53 mm day⁻¹ on 17 December) with the majority of runoff due to rainfall passing through the snowpack as opposed to snowmelt. For both rain‐on‐snow and radiation melt events lysimeter discharge was generally greatest at the open site, although there were exceptions such as during interception melt events. During radiation melt instantaneous discharge was up to 4·0 times greater in the opening compared with the mature cedar, and 48 h discharge was up to 2·5 times greater. Perhaps characteristic of maritime climates, forest interception melt is shown to be important in addition to sublimation in reducing snow accumulation beneath dense canopies. While sublimation represents a loss from the catchment water balance, interception melt percolates through the snowpack and contributes to soil moisture during the winter season. Strong differences in microclimate and snowpack albedo persisted between cedar, larch and open sites, and it is suggested further work is needed to account for this in hydrological simulation models. Copyright © 2005 John Wiley & Sons, Ltd.     

A multiplier-based method of generating stochastic areal rainfall from point rainfalls

Catchment modelling for water resources assessment is still mainly based on rain gauge measurements as these are more easily available and cover longer periods than radar and satellite-based measurements. Rain gauges however measure the rain falling on an extremely small proportion of the catchment and the areal rainfall obtained from these point measurements are consequently substantially uncertain. These uncertainties in areal rainfall estimation are generally ignored and the need to assess their impact on catchment modelling and water resources assessment is therefore imperative. A method that stochastically generates daily areal rainfall from point rainfall using multiplicative perturbations as a means of dealing with these uncertainties is developed and tested on the Berg catchment in the Western Cape of South Africa. The differences in areal rainfall obtained by alternately omitting some of the rain gauges are used to obtain a population of plausible multiplicative perturbations. Upper bounds on the applicable perturbations are set to prevent the generation of unrealistically large rainfall and to obtain unbiased stochastic rainfall. The perturbations within the set bounds are then fitted into probability density functions to stochastically generate the perturbations to impose on areal rainfall. By using 100 randomly-initialized calibrations of the AWBM catchment model and Sequent Peak Analysis, the effects of incorporating areal rainfall uncertainties on storage-yield-reliability analysis are assessed. Incorporating rainfall uncertainty is found to reduce the required storage by up to 20%. Rainfall uncertainty also increases flow-duration variability considerably and reduces the median flow-duration values by an average of about 20%.     

Which catchment characteristics control the temporal dependence structure of daily river flows?

Hydrological classification systems seek to provide information about the dominant processes in the catchment to enable information to be transferred between catchments. Currently, there is no widely agreed‐upon system for classifying river catchments. This paper develops a novel approach to classifying catchments based on the temporal dependence structure of daily mean river flow time series, applied to 116 near‐natural ‘benchmark’ catchments in the UK. The classification system is validated using 49 independent catchments. Temporal dependence in river flow data is driven by the flow pathways, connectivity and storage within the catchment and can thus be used to assess the influence catchment characteristics have on moderating the precipitation‐to‐flow relationship. Semi‐variograms were computed for the 116 benchmark catchments to provide a robust and efficient way of characterising temporal dependence. Cluster analysis was performed on the semi‐variograms, resulting in four distinct clusters. The influence of a wide range of catchment characteristics on the semi‐variogram shape was investigated, including: elevation, land cover, physiographic characteristics, soil type and geology. Geology, depth to gleyed layer in soils, slope of the catchment and the percentage of arable land were significantly different between the clusters. These characteristics drive the temporal dependence structure by influencing the rate at which water moves through the catchment and/or the storage in the catchment. Quadratic discriminant analysis was used to show that a model with five catchment characteristics is able to predict the temporal dependence structure for un‐gauged catchments. This method could form the basis for future regionalisation strategies, as a way of transferring information on the precipitation‐to‐flow relationship between gauged and un‐gauged catchments. © 2014 The Authors. Hydrological Processes by published by John Wiley & Sons, Ltd.     

Root‐zone moisture replenishment in a native vegetated catchment under Mediterranean climate

The root‐zone moisture replenishment mechanisms are key unknowns required to understand soil hydrological processes and water sources used by plants. Temporal patterns of root‐zone moisture replenishment reflect wetting events that contribute to plant growth and survival and to catchment water yield. In this study, stable oxygen and hydrogen isotopes of twigs and throughfall were continuously monitored to characterize the seasonal variations of the root‐zone moisture replenishment in a native vegetated catchment under Mediterranean climate in South Australia. The two studied hillslopes (the north‐facing slope [NFS] and the south‐facing slope [SFS]) had different environmental conditions with opposite aspects. The twig and throughfall samples were collected every ~20 days over 1 year on both hillslopes. The root‐zone moisture replenishment, defined as percentage of newly replenished root‐zone moisture as a complement to antecedent moisture for plant use, calculated by an isotope balance model, was about zero (±25% for the NFS and ± 15% for the SFS) at the end of the wet season (October), increased to almost 100% (±26% for the NFS and ± 29% for the SFS) after the dry season (April and May), then decreased close to zero (±24% for the NFS and ± 28% for the SFS) in the middle of the following wet season (August). This seasonal pattern of root‐zone moisture replenishment suggests that the very first rainfall events of the wet season were significant for soil moisture replenishment and supported the plants over wet and subsequent dry seasons, and that NFS completed replenishment over a longer time than SFS in the wet season and depleted the root zone moisture quicker in the dry season. The stable oxygen isotope composition of the intraevent samples and twigs further confirms that rain water in the late wet season contributed little to root‐zone moisture. This study highlights the significant role of the very first rain events in the early wet season for ecosystem and provides insights to understanding ecohydrological separation, catchment water yield, and vegetation response to climate changes.     

Rain‐gauge network evaluation and augmentation using geostatistics

Rain‐gauge networks are often used to provide estimates of area average rainfall or point rainfalls at ungauged locations. The level of accuracy a network can achieve depends on the total number and locations of gauges in the network. A geostatistical approach for evaluation and augmentation of an existing rain‐gauge network is proposed in this study. Through variogram analysis, hourly rainfalls are shown to have higher spatial variability than annual rainfalls, with hourly Mei‐Yu rainfalls having the highest spatial variability. A criterion using ordinary kriging variance is proposed to assess the accuracy of rainfall estimation using the acceptance probability defined as the probability that estimation error falls within a desired range. Based on the criterion, the percentage of the total area with acceptable accuracy A_p under certain network configuration can be calculated. A sequential algorithm is also proposed to prioritize rain‐gauges of the existing network, identify the base network, and relocate non‐base gauges. Percentage of the total area with acceptable accuracy is mostly contributed by the base network. In contrast, non‐base gauges provide little contribution to A_p and are subject to removal or relocation. Using a case study in northern Taiwan, the proposed approach demonstrates that the identified base network which comprises of approximately two‐thirds of the total rain‐gauges can achieve almost the same level of performance (expressed in terms of percentage of the total area with acceptable accuracy) as the complete network for hourly Mei‐Yu rainfall estimation. The percentage of area with acceptable accuracy can be raised from 56% to 88% using an augmented network. A threshold value for the percentage of area with acceptable accuracy is also recommended to help determine the number of non‐base gauges which need to be relocated. Copyright © 2007 John Wiley & Sons, Ltd.     

A model to simulate recharge processes of karst massifs

The recharge processes have been evaluated for two karst massifs of southern Italy, the Mt Terminio and Mt Cervialto, characterized by wide endorheic areas. The annual mean recharge has been estimated by Geographic Information System (GIS) tools, from regression of annual mean values of different ground‐elevated rain gauges and thermometers. The recharge has been distinguished for endorheic areas and the other areas of spring catchment, and the ratio between the output spring and input rainfall has been also estimated (recharge coefficient). The annual recharge has been used to calibrate a daily scale model, which allows to estimate the amount of effective rainfall, which is retained as soil moisture; the amount reaching the water table (recharge s.s.); and the amount of rainfall, which develops the runoff and leaves the catchment. All these amounts vary through the hydrological year, in function of soil moisture deficit and daily rainfall intensity. The model allows estimating the recharge conditions through the hydrological year, and it is a useful tool for water management. Copyright © 2014 John Wiley & Sons, Ltd.     

POTENTIAL OF EXTENDING THE RAINFALL INTENSITY–DURATION–FREQUENCY RELATIONSHIP TO NON-RECORDING RAIN GAUGES

Maximum rainfall intensity–duration–frequency (IDF) curves are commonly applied to determine the design rainfall in water resource projects. Normally, the IDF relationship is derived from recording rain gauges. As the network of non-recording rain gauges (daily rainfall) in Taiwan has a higher density than recording rain gauges, attempts were made in this study to extend the IDF relationship to non-recording rain gauges. Eighteen recording rain gauges and 99 non-recording rain gauges over the Chi-Nan area in Southern Taiwan provide the data sets. The regional IDF formulae were generated for ungauged areas to estimate rainfall intensity for various return periods and rainfall durations larger than or equal to one hour. For rainfall durations less than one hour, a set of adjustment formulae were applied to modify the regional IDF formulae. The method proposed in this study had reasonable application to non-recording rain gauges, which was concluded from the verification of four additional recording rain gauges. © 1997 John Wiley & Sons, Ltd.     

Recent climatic trends and linkages to river discharge in Central Vietnam

In this study, change in rainfall, temperature and river discharge are analysed over the last three decades in Central Vietnam. Trends and rainfall indices are evaluated using non‐parametric tests at different temporal levels. To overcome the sparse locally available network, the high resolution APHRODITE gridded dataset is used in addition to the existing rain gauges. Finally, existing linkages between discharge changes and trends in rainfall and temperature are explored. Results are indicative of an intensification of rainfall (+15%/decade), with more extreme and longer events. A significant increase in winter rainfall and a decrease in consecutive dry days provides strong evidence for a lengthening wet season in Central Vietnam. In addition, trends based on APHRODITE suggest a strong orographic signal in winter and annual trends. These results underline the local variability in the impacts of climatic change at the global scale. Consequently, it is important that change detection investigations are conducted at the local scale. A very weak signal is detected in the trend of minimum temperature (+0.2°C/decade). River discharge trends show an increase in mean discharge (31 to 35%/decade) over the last decades. Between 54 and 74% of this increase is explained by the increase in precipitation. The maximum discharge also responds significantly to precipitation changes leading to a lengthened wet season and an increase in extreme rainfall events. Such trends can be linked with a likely increase in floods in Central Vietnam, which is important for future adaptation planning and management and flood preparedness in the region. Copyright © 2012 John Wiley & Sons, Ltd.     

Correlation functions of rainfall field and their application in network design in the tropics

The study examines the correlation function of tropical monsoon rainfall on monthly, seasonal and annual time scales and obtains the relationship between this function and the distance. The area selected for study is Vidarbha with a fairly dense network of rain gauges. Vidarbha is a meteorological sub-division of the state of Maharashtra in India. Utilizing the relationship between the correlation function of the rainfall field and the distance, the errors of optimum interpolation of rainfall at a point have been computed by applying the method of optimum interpolation byGandin (1970). Relationships between the errors of interpolation and distance have been evaluated and from this the maximum spacing allowed between rain gauges for a specified tolerable error in interpolation has been estimated for each of the periods.     

Flood simulation using the gauge‐adjusted radar rainfall and physics‐based distributed hydrologic model

This paper reports the results of an investigation into flood simulation by areal rainfall estimated from the combination of gauged and radar rainfalls and a rainfall–runoff model on the Anseong‐cheon basin in the southern part of Korea. The spatial and temporal characteristics and behaviour of rainfall are analysed using various approaches combining radar and rain gauges: (1) using kriging of the rain gauge alone; (2) using radar data alone; (3) using mean field bias (MFB) of both radar and rain gauges; and (4) using conditional merging technique (CM) of both radar and rain gauges. To evaluate these methods, statistics and hyetograph for rain gauges and radar rainfalls were compared using hourly radar rainfall data from the Imjin‐river, Gangwha, rainfall radar site, Korea. Then, in order to evaluate the performance of flood estimates using different rainfall estimation methods, rainfall–runoff simulation was conducted using the physics‐based distributed hydrologic model, Vflo^?. The flood runoff hydrograph was used to compare the calculated hydrographs with the observed one. Results show that the rainfall field estimated by CM methods improved flood estimates, because it optimally combines rainfall fields representing actual spatial and temporal characteristics of rainfall. Copyright © 2008 John Wiley & Sons, Ltd.     

Optimal design of rain gauge network in the Middle Yarra River catchment,Australia

Rainfall data are a fundamental input for effective planning, designing and operating of water resources projects. A well‐designed rain gauge network is capable of providing accurate estimates of necessary areal average and/or point rainfall estimates at any desired ungauged location in a catchment. Increasing network density with additional rain gauge stations has been the main underlying criterion in the past to reduce error and uncertainty in rainfall estimates. However, installing and operation of additional stations in a network involves large cost and manpower. Hence, the objective of this study is to design an optimal rain gauge network in the Middle Yarra River catchment in Victoria, Australia. The optimal positioning of additional stations as well as optimally relocating of existing redundant stations using the kriging‐based geostatistical approach was undertaken in this study. Reduction of kriging error was considered as an indicator for optimal spatial positioning of the stations. Daily rainfall records of 1997 (an El Niño year) and 2010 (a La Niña year) were used for the analysis. Ordinary kriging was applied for rainfall data interpolation to estimate the kriging error for the network. The results indicate that significant reduction in the kriging error can be achieved by the optimal spatial positioning of the additional as well as redundant stations. Thus, the obtained optimal rain gauge network is expected to be appropriate for providing high quality rainfall estimates over the catchment. The concept proposed in this study for optimal rain gauge network design through combined use of additional and redundant stations together is equally applicable to any other catchment. © 2014 The Authors. Hydrological Processes published by John Wiley & Sons Ltd.     

Statistical corrections of spatially interpolated missing precipitation data estimates

Spatial interpolation methods used for estimation of missing precipitation data generally under and overestimate the high and low extremes, respectively. This is a major limitation that plagues all spatial interpolation methods as observations from different sites are used in local or global variants of these methods for estimation of missing data. This study proposes bias‐correction methods similar to those used in climate change studies for correcting missing precipitation estimates provided by an optimal spatial interpolation method. The methods are applied to post‐interpolation estimates using quantile mapping, a variant of equi‐distant quantile matching and a new optimal single best estimator (SBE) scheme. The SBE is developed using a mixed‐integer nonlinear programming formulation. K‐fold cross validation of estimation and correction methods is carried out using 15 rain gauges in a temperate climatic region of the U.S. Exhaustive evaluation of bias‐corrected estimates is carried out using several statistical, error, performance and skill score measures. The differences among the bias‐correction methods, the effectiveness of the methods and their limitations are examined. The bias‐correction method based on a variant of equi‐distant quantile matching is recommended. Post‐interpolation bias corrections have preserved the site‐specific summary statistics with minor changes in the magnitudes of error and performance measures. The changes were found to be statistically insignificant based on parametric and nonparametric hypothesis tests. The correction methods provided improved skill scores with minimal changes in magnitudes of several extreme precipitation indices. The bias corrections of estimated data also brought site‐specific serial autocorrelations at different lags and transition states (dry‐to‐dry, dry‐to‐wet, wet‐to‐wet and wet‐to‐dry) close to those from the observed series. Bias corrections of missing data estimates provide better serially complete precipitation time series useful for climate change and variability studies in comparison to uncorrected filled data series. Copyright © 2013 John Wiley & Sons, Ltd.     

Characteristics of rain events at an inland locality in northeastern Borneo,Malaysia

Understanding the intensity and duration of tropical rain events is critical to modelling the rate and timing of wet‐canopy evaporation, the suppression of transpiration, the generation of infiltration‐excess overland flow and hence to erosion, and to river responsiveness. Despite this central role, few studies have addressed the characteristics of equatorial rainstorms. This study analyses rainfall data for a 5 km² region largely comprising of the 4 km² Sapat Kalisun Experimental Catchment in the interior of northeastern Borneo at sampling frequencies from 1 min^?1 to 1 day^?1. The work clearly shows that most rainfall within this inland, forested area is received during regular short‐duration events (<15 min) that have a relatively low intensity (i.e. less than two 0·2 mm rain‐gauge tips in almost all 5 min periods). The rainfall appears localized, with significant losses in intergauge correlations being observable in minutes in the case of the typical mid‐afternoon, convective events. This suggests that a dense rain‐gauge network, sampled at a high temporal frequency, is required for accurate distributed rainfall‐runoff modelling of such small catchments. Observed rain‐event intensity is much less than the measured infiltration capacities, and thus supports the tenet of the dominance of quick subsurface responses in controlling river behaviour in this small equatorial catchment. Copyright © 2006 John Wiley & Sons, Ltd.