Construction of design hyetographs for locations without observed data

A design hyetograph which represents the time distribution of design rainfall depth corresponding to a duration and a return period is essential in hydrologic design. However, for locations without observed data (ungauged sites), construction of design hyetographs is a difficult task because of the lack of data. Hence, an approach based on self‐organizing map (SOM) is proposed in this paper to construct design hyetographs at ungauged sites. SOM, which is a special kind of artificial neural networks (ANNs), is a powerful technique for extracting and visualizing salient features of data and for solving classification problems. The proposed approach is composed of three steps: classification, assignment and construction. First, the SOM‐based classification is performed to analyse gauged sites' design hyetographs. Second, based on the concept of indicator kriging, a method is developed to assign an ungauged site of interest to a certain cluster. Third, based on the spatial information, the clustering results, and the design hyetographs of gauged sites, the design hyetograph at the site of interest is constructed using the reciprocal‐distance‐squared method. An application is conducted to assess the advantages of the proposed approach over the conventional approaches. Moreover, cross‐validation tests are applied to evaluate the performance of the accuracy and the robustness of the proposed approach. The results confirm the improvement in performance by using the proposed approach instead of conventional approaches. The proposed approach is useful for constructing design hyetographs at ungauged sites. Copyright © 2009 John Wiley & Sons, Ltd.     

Application of L‐moments and Bayesian inference for low‐flow regionalization in Sefidroud basin,Iran

Reliable estimation of low flows at ungauged catchments is one of the major challenges in water‐resources planning and management. This study aims at providing at‐site and ungauged sites low‐flow frequency analysis using regionalization approach. A two‐stage delineating homogeneous region is proposed in this study. Clustering sites with similar low‐flow L‐moment ratios is initially conducted, and L‐moment‐based discordancy and heterogeneity measures are then used to detect unusual sites. Based on the goodness‐of‐fit test statistic, the best‐fit regional model is identified in each hydrologically homogeneous region. The relationship between mean annual 7‐day minimum flow and hydro‐geomorphic characteristics is also constructed in each homogeneous region associated with the derived regional model for estimating various low‐flow quantiles at ungauged sites. Uncertainty analysis of model parameters and low‐flow estimations is carried out using the Bayesian inference. Applied in Sefidroud basin located in northwestern Iran, two hydrologically homogeneous regions are identified, i.e. the east and west regions. The best‐fit regional model for the east and west regions are generalized logistic and Pearson type III distributions, respectively. The results show that the proposed approach provides reasonably good accuracy for at‐site as well as ungauged‐site frequency analysis. Besides, interval estimations for model parameters and low flows provide uncertainty information, and the results indicate that Bayesian confidence intervals are significantly reduced when comparing with the outcomes of conventional t‐distribution method. Copyright © 2013 John Wiley & Sons, Ltd.     

A regional model of low flow for southern Taiwan

The primary purpose of this study is to develop regional models of the lower part of flow duration curves (LPFDCs) to synthesize low‐flow characteristics at ungauged sites in southern Taiwan. Because of the close relationship between low streamflow regimes and hydrogeological features, the model development first involved delimiting homogeneous hydrogeological regions by using two‐step cluster analysis. Each homogeneous region was then discriminated by an equation developed on the basis of its hydrogeological features, which was then used to determine which of three sets of regional LPFDC models would be appropriate for a particular ungauged site. Each of the three sets of regional LPFDC models were developed using both conventional multivariate statistical regression and fuzzy regression. Thirty‐four stream‐gauged watersheds located in southern Taiwan provide the data set. The study results reveal that the regional LPFDC models developed in this study could be applied reasonably at ungauged sites. Copyright © 2002 John Wiley & Sons, Ltd.     

Using synthetic flow duration curves for rainfall–runoff model calibration at ungauged sites

The primary objective of the study is to propose a strategy for rainfall–runoff model calibration at ungauged sites. This strategy comprises two main components: (1) development of the regional analysis method to synthesize the flow duration curves at ungauged sites; and (2) utilization of the synthetic flow duration curves for model calibration. Since the regional analysis method can synthesize the flow duration curves at ungauged sites, the continuous rainfall–runoff model coupled with a global optimization method were applied in southern Taiwan using the synthetic flow duration curve as an objective for model calibration. The results reveal that the regional flow duration curve and the strategy for model calibration at ungauged sites have good performances in the study area. Copyright © 2000 John Wiley & Sons, Ltd.     

Storm‐event rainfall–runoff modelling approach for ungauged sites in Taiwan

This work develops a top‐down modelling approach for storm‐event rainfall–runoff model calibration at unmeasured sites in Taiwan. Twenty‐six storm events occurring in seven sub‐catchments in the Kao‐Ping River provided the analytical data set. Regional formulas for three important features of a streamflow hydrograph, i.e. time to peak, peak flow, and total runoff volume, were developed via the characteristics of storm event and catchment using multivariate regression analysis. Validation of the regional formulas demonstrates that they reasonably predict the three features of a streamflow hydrograph at ungauged sites. All of the sub‐catchments in the study area were then adopted as ungauged areas, and the three streamflow hydrograph features were calculated by the regional formulas and substituted into the fuzzy multi‐objective function for rainfall–runoff model calibration. Calibration results show that the proposed approach can effectively simulate the streamflow hydrographs at the ungauged sites. The simulated hydrographs more closely resemble observed hydrographs than hydrographs synthesized using the Soil Conservation Service (SCS) dimensionless unit hydrograph method, a conventional method for hydrograph estimation at ungauged sites in Taiwan. Copyright © 2008 John Wiley & Sons, Ltd.     

A feature weighting and selection method for improving the homogeneity of regions in regionalization of watersheds

The method used for feature selection or feature weighting in regionalization of watersheds may affect the results of regionalization methods considerably. It can play a key role in forming hydrologically homogeneous regions for regional flood frequency analysis. In this study, a method based on exploring the nearest and farthest neighbours of data points is introduced for identifying salient features for regionalization of watersheds. The method includes options to relate watershed features to flood data records in order to increase the homogeneity of the regions. The nearest and farthest neighbours are identified based on the criteria such as the mutual information criterion and Spearman's rank correlation coefficient. Then, the watershed features more able to explain the relationships between the nearest and farthest neighbours are identified as salient features to form homogeneous features for regional flood frequency analysis. The results show that the optimum option of the proposed method improves the performances of the hard and fuzzy clustering algorithms in more than half of the cases based on the cluster validity indices. Furthermore, the results reveal that the optimum option can increase the number of the homogeneous regions formed by clustering algorithms to a great extent. By using the optimum option with 5 nearest and 5 farthest neighbours, longitude, drainage area, and run‐off coefficient are identified as the salient features to regionalize Sefidrud basin. The results show that the proposed method can be considered as an efficient method to form homogeneous regions for regional flood frequency analysis.     

Synthetic regional flow duration curve for southern Taiwan

The primary purpose of this study is to develop the regional flow duration curves for southern Taiwan. To define homogeneous regions for developing regional flow duration curves, multivariate statistical analysis (principal component and cluster analysis) was applied to daily flow data from 34 stream-gauged stations in southern Taiwan. Two kinds of clustering variables, the dimensionless flow duration curve and specific flow duration curve, were compared in this study. It was found that three homogeneous regions delineated by specific flow duration curves as clustering variables have more reasonable results. The three homogeneous regions not only have well-defined geographical boundaries, but also correspond to the rainfall and geology characteristics of the regions. It seems that the technique of cluster analysis can reasonably define the homogeneous regions. In each homogeneous region, the synthetic regional flow duration curves were developed by a family of parametric duration curves. This approach has the advantage of being simple and needing only the basin area as an index. The performance of the regional flow duration curve was verified by the comparison of areas under the actual and synthetic flow duration curves; the latter were generated from the regional flow duration curve. Almost all the 34 stream-gauged stations had less than 25% absolute error.     

Design flood estimation in ungauged catchments and statistical characterization using principal components analysis: application of Gradex method in Upper Moulouya

Design flood estimation in ungauged catchments is of great importance in hydrologic practice especially where there is no available data about streamflow. Except the watershed of Anseghmir who is equipped with a gauge station, all the other watersheds are ungauged catchments. The use of frequency analysis of series of rainfall and streamflow is very important for the characterization of the hydrologic resources of the Upper Moulouya. The region has a semiarid climate that requires a good knowledge of the watershed's potential water to assist policy makers in forecasting extreme events, managing water resources and decision making. The frequency analysis was used to determine the design flood of different return periods. The results obtained are used in Gradex method to estimate the hydrologic variables of each subcatchment of the Upper Moulouya. Once the hydrologic study is completed, a principal components analysis was made to highlight the affinities between the different subcatchments and to deduce the hydrologic and hydrographic parameters that better characterize them. Copyright © 2012 John Wiley & Sons, Ltd.     

A canonical correlation approach to the determination of homogeneous regions for regional flood estimation of ungauged basins

Abstract

A canonical correlation method for determining the homogeneous regions used for estimating flood characteristics of ungauged basins is described. The method emphasizes graphical and quantitative analysis of relationships between the basin and flood variables before the data of the gauged basins are used for estimating the flood variables of the ungauged basin. The method can be used for both homogeneous regions, determined a priori by clustering algorithms in the space of the flood-related canonical variables, as well as for “regions of influence” or “neighbourhoods” centred on the point representing the estimated location of the ungauged basin in that space.     

Formation of homogeneous regions for regional frequency analysis of extreme precipitation events in the Czech Republic

Extreme high precipitation amounts are among environmental events with the most disastrous consequences for human society. This paper deals with the identification of ‘homogeneous regions’ according to statistical characteristics of precipitation extremes in the Czech Republic, i.e. the basic and most important step toward the regional frequency analysis. Precipitation totals measured at 78 stations over 1961–2000 are used as an input dataset. Preliminary candidate regions are formed by the cluster analysis of site characteristics, using the average-linkage clustering and Ward’s method. Several statistical tests for regional homogeneity are utilized, based on the 10-yr event and the variation of L-moment statistics. In compliance with results of the tests, the area of the Czech Republic has been divided into four homogeneous regions. The findings are supported by simulation experiments proposed to evaluate stability of the test results. Since the regions formed reflect also climatological differences in precipitation regimes and synoptic patterns causing high precipitation amounts, their future application may not be limited to the frequency analysis of extremes.     

Regional flood frequency analysis using residual kriging in physiographical space

In this article, an approach using residual kriging (RK) in physiographical space is proposed for regional flood frequency analysis. The physiographical space is constructed using physiographical/climatic characteristics of gauging basins by means of canonical correlation analysis (CCA). This approach is a modified version of the original method, based on ordinary kriging (OK). It is intended to handle effectively any possible spatial trends within the hydrological variables over the physiographical space. In this approach, the trend is first quantified and removed from the hydrological variable by a quadratic spatial regression. OK is therefore applied to the regression residual values. The final estimated value of a specific quantile at an ungauged station is the sum of the spatial regression estimate and the kriged residual. To evaluate the performance of the proposed method, a cross‐validation procedure is applied. Results of the proposed method indicate that RK in CCA physiographical space leads to more efficient estimates of regional flood quantiles when compared to the original approach and to a straightforward regression‐based estimator. Copyright © 2010 John Wiley & Sons, Ltd.     

Clustering spatial–temporal precipitation data using wavelet transform and self-organizing map neural network

A data analysis method is proposed to cluster and explore spatio-temporal characteristics of the 22 years of precipitation data (1982–2003) for Taiwan. The wavelet transform self-organizing map (WTSOM) framework combines the wavelet transform (WT) and a self-organizing map (SOM) neural network. WT is used to extract dynamic and multiscale features of the non-stationary precipitation time-series, and SOM is applied to objectively identify spatially homogeneous clusters on the high-dimensional wavelet-transformed feature space. Haar and Morlet wavelets are applied in the data preprocessing stage to preserve the desired characteristics of the precipitation data. A two-level SOM neural network is applied to identify clusters in the wavelet space in the clustering stage. The performance of clustering is evaluated using silhouette coefficients. The results indicate that singularities or sharp transitions are more significant than changes in the periodicity or data structure in the spatial–temporal precipitation data. The WTSOM results show that six clusters are optimal for both Haar and Morlet wavelet functions, but their corresponding geographic locations are different. The geographic locations of clusters based on the Haar wavelet, which captures the occurrence of extreme hydrological events, appear in blocks while those classified by the Morlet wavelet, which indicates periodicity changes and describes fine structures, appear in strips that cross the island of Taiwan. Principal component analysis is applied to the precipitation data of each cluster. The first principal components explain 62–90% of the total variation of data. Characteristics of precipitation data for each cluster are explored using scalogram analysis. The results show that both extreme hydrological events and periodicity changes appear in the spatial and temporal precipitation data but with different characteristics for each cluster. Recognizing homogeneous hydrologic regions and identifying the associated precipitation characteristics improves the efficiency of water resources management in adapting to climate change, preventing the degradation of the water environment, and reducing the impact of climate-induced disasters. Measures for countering the stress of precipitation variation for water resources management are provided.     

Regional frequency analysis of rainfall extremes in Southern Malawi using the index rainfall and L-moments approaches

Rainfall extremes often result in the occurrence of flood events with associated loss of life and infrastructure in Malawi. However, an understanding of the frequency of occurrence of such extreme events either for design or disaster planning purposes is often limited by data availability at the desired temporal and spatial scales. Regionalisation, which involves “trading time for space” by pooling together observations for stations with similar behavior, is an alternative approach for more accurate determination of extreme events even at ungauged areas or sites with short records. In this study, regional frequency analysis of rainfall extremes in Southern Malawi, large parts of which are flood prone, was undertaken. Observed 1-, 3-, 5- and 7-day annual maximum rainfall series for the period 1978–2007 at 23 selected rainfall stations in Southern Malawi were analysed. Cluster analysis using scaled at-site characteristics was used to determine homogeneous rainfall regions. L-moments were applied to derive regional index rainfall quantiles. The procedure also validated the three rainfall regions identified through homogeneity and heterogeneity tests based on Monte Carlo simulations with regional average L-moment ratios fitted to the Kappa distribution. Based on assessments of the accuracy of the derived index rainfall quantiles, it was concluded that the performance of this regional approach was satisfactory when validated for sites not included in the sample data. The study provides an estimate of the regional characteristics of rainfall extremes that can be useful in among others flood mitigation and engineering design.     

Curve‐Number/Green–Ampt mixed procedure for streamflow predictions in ungauged basins: Parameter sensitivity analysis

A procedure combining the Soil Conservation Service‐Curve Number (SCS‐CN) method and the Green–Ampt (GA) infiltration equation was recently developed to overcome some of the drawbacks of the classic SCS‐CN approach when estimating the volume of surface runoff at a sub‐daily time resolution. The rationale of this mixed procedure, named Curve Number for Green–Ampt (CN4GA), is to use the GA infiltration model to distribute the total volume of the net hyetograph (rainfall excess) provided by the SCS‐CN method over time. The initial abstraction and the total volume of rainfall given by the SCS‐CN method are used to identify the ponding time and to quantify the hydraulic conductivity parameter of the GA equation. In this paper, a sensitivity analysis of the mixed CN4GA parameters is presented with the aim to identify conditions where the mixed procedure can be effectively used within the Prediction in Ungauged Basin perspective. The effects exerted by changes in selected input parameters on the outputs are evaluated using rectangular and triangular synthetic hyetographs as well as 100 maximum annual storms selected from synthetic rainfall time series. When applied to extreme precipitation events, which are characterized by predominant peaks of rainfall, the CN4GA appears to be rather insensitive to the input hydraulic parameters of the soil, which is an interesting feature of the CN4GA approach and makes it an ideal candidate for the rainfall excess estimation at sub‐daily temporal resolution at ungauged sites. Copyright © 2013 John Wiley & Sons, Ltd.     

Regionalization of precipitation characteristics in the Canadian Prairie Provinces using large-scale atmospheric covariates and geophysical attributes

Observed data at most stations are often inadequate to obtain reliable estimates of many hydro-meteorological variables that not only define water availability across a region but also the vulnerability of social infrastructure to climatic extremes. To overcome this, data from neighboring sites with similar statistical characteristics are often pooled. The pooling process is based on partitioning of a larger region into smaller sub-regions with homogeneous features of interest. The established approaches rely heavily on statistics computed from observed precipitation data rather than the covariates that play a significant role in modulating the regional and local climate patterns at various temporal and spatial scales. In this study, a new approach for identifying homogeneous regions for regionalization of precipitation characteristics is proposed for the Canadian Prairie Provinces. This approach incorporates information about large-scale atmospheric covariates, teleconnection indices and geographical site attributes that impact spatial patterns of precipitation in order to delineate homogeneous precipitation regions through combined use of multivariate approaches—principal component analysis, canonical correlation analysis and fuzzy C-means clustering. Results of the analyses suggest that the study area can be partitioned into five homogeneous regions. These partitions are validated independently for homogeneity using statistics computed from monthly and seasonal precipitation totals, and seasonal extremes from a network of observation stations. Furthermore, based on the identified regions, precipitation magnitude-frequency relationships of warm and cold season single- and multi-day precipitation extremes, developed through regional frequency analysis, are mapped spatially. Such estimates are important for numerous water resources related activities.     

Regionalisation of the parameters of the log‐Pearson 3 distribution: a case study for New South Wales,Australia

The index flood method is widely used in regional flood frequency analysis (RFFA) but explicitly relies on the identification of ‘acceptable homogeneous regions’. This paper presents an alternative RFFA method, which is particularly useful when ‘acceptably homogeneous regions’ cannot be identified. The new RFFA method is based on the region of influence (ROI) approach where a ‘local region’ can be formed to estimate statistics at the site of interest. The new method is applied here to regionalize the parameters of the log‐Pearson 3 (LP3) flood probability model using Bayesian generalized least squares (GLS) regression. The ROI approach is used to reduce model error arising from the heterogeneity unaccounted for by the predictor variables in the traditional fixed‐region GLS analysis. A case study was undertaken for 55 catchments located in eastern New South Wales, Australia. The selection of predictor variables was guided by minimizing model error. Using an approach similar to stepwise regression, the best model for the LP3 mean was found to use catchment area and 50‐year, 12‐h rainfall intensity as explanatory variables, whereas the models for the LP3 standard deviation and skewness only had a constant term for the derived ROIs. Diagnostics based on leave‐one‐out cross validation show that the regression model assumptions were not inconsistent with the data and, importantly, no genuine outlier sites were identified. Significantly, the ROI GLS approach produced more accurate and consistent results than a fixed‐region GLS model, highlighting the superior ability of the ROI approach to deal with heterogeneity. This method is particularly applicable to regions that show a high degree of regional heterogeneity. Copyright © 2014 John Wiley & Sons, Ltd.     

Estimation of rainfall intensity–duration–frequency curves at ungauged locations using quantile regression methods

Intensity–duration–frequency (IDF) curves of extreme rainfall are used extensively in infrastructure design and water resources management. In this study, a novel regional framework based on quantile regression (QR) is used to estimate rainfall IDF curves at ungauged locations. Unlike standard regional approaches, such as index-storm and at-site ordinary least-squares regression, which are dependent on parametric distributional assumptions, the non-parametric QR approach directly estimates rainfall quantiles as a function of physiographic characteristics. Linear and nonlinear methods are evaluated for both the regional delineation and IDF curve estimation steps. Specifically, delineation by canonical correlation analysis (CCA) and nonlinear CCA (NLCCA) is combined, in turn, with linear QR and nonlinear QR estimation in a regional modelling framework. An exhaustive comparative study is conducted between standard regional methods and the proposed QR framework at sites across Canada. Overall, the fully nonlinear QR framework, which uses NLCCA for delineation and nonlinear QR for estimation of IDF curves at ungauged sites, leads to the best results.     

Estimating minimum streamflow from measurements at ungauged sites in regions with streamflow‐gauging networks

Estimation of low flows in rivers continues to be a vexing problem despite advances in statistical and process‐based hydrological models. We develop a method to estimate minimum streamflow at seasonal to annual timescales from measured streamflow based on regional similarity in the deviations of daily streamflow from minimum streamflow for a period of interest. The method is applied to 1,019 gauged sites in the Western United States for June to December 2015. The gauges were clustered into six regions with distinct timing and magnitude of low flows. A gamma distribution was fit each day to the deviations in specific discharge (daily streamflow divided by drainage area) from minimum specific discharge for gauges in each region. The Kolmogorov–Smirnov test identified days when the gamma distribution was adequate to represent the distribution of deviations in a region. The performance of the gamma distribution was evaluated at gauges by comparing daily estimates of minimum streamflow with estimates from area‐based regression relations for minimum streamflow. Each region had at least 8 days during the period when streamflow measurements would provide better estimates than the regional regression equation, but the number of such days varied by region depending on aridity and homogeneity of streamflow within the region. Synoptic streamflow measurements at ungauged sites have value for estimating minimum streamflow and improving the spatial resolution of hydrological model in regions with streamflow‐gauging networks.     

A Methodology for Clustering Lakes in Alberta on the basis of Water Quality Parameters

In this study, a methodology for clustering 18 lakes in Alberta, Canada using the data of 19 water quality parameters for a period of 11 years (1988–2002) is presented. The methods consist of (i) principal component analysis (PCA) to determine the dominant water quality parameters, (ii) cluster analysis techniques to develop the characteristics of the clusters, and (iii) pattern‐match lakes to determine the appropriate cluster for each of the lakes. The PCA revealed that three principal components (PCs) were able to explain ～88% of the variability and the dominant water quality parameters were total dissolved solids, total phosphorus, and chlorophyll‐a. We obtained five clusters for the period 1994–1997 by using the dominant parameters with water quality deteriorating as the cluster number increased from 1 to 5. Upon matching cluster patterns with the entire dataset, it was observed that some of the lakes belonged to the same cluster all the time (e.g., cluster 1 for lakes Elkwater, Gregg, and Jarvis; cluster 3 for Sturgeon; cluster 4 for Moonshine; and cluster 5 for Saskatoon), while others changed with time. This methodology could be applied in other regions of the world to identify the most suitable source waters and prioritize their management. It could be helpful to analyze the natural controlling processes, pollution types, impact of seasonal changes and overall quality of source waters. This methodology could be used for monitoring water bodies in a cost effective and efficient way by sampling only less number of dominant parameters instead of using a large set of parameters.