Regional bivariate modeling of droughts using L-comoments and copulas

The regional bivariate modeling of drought characteristics using the copulas provides valuable information for water resources management and drought risk assessment. The regional frequency analysis (RFA) can specify the similar sites within a region using L-comoments approach. One of the important steps in the RFA is estimating regional parameters of the copula function. In the present study, an optimization-based method along with the adjusted charged system search are introduced and applied to estimate the regional parameters of the copula models. The capability of the proposed methodology is illustrated by copula functions on drought events. Three commonly used copulas containing Clayton, Frank and Gumbel are employed to derive the joint distribution of drought severity and duration. The result of the new method are compared to the method of moments and after applying several goodness-of-fit tests, the results indicate that the new method provides higher accuracy than the classic one. Furthermore, the results of the upper tail dependence coefficient indicate that the Gumbel copula is the best-fitted copula among the other ones for modeling drought characteristics.     

Comparison of parameter uncertainty analysis techniques for a TOPMODEL application

Parameter uncertainty in hydrologic modeling is crucial to the flood simulation and forecasting. The Bayesian approach allows one to estimate parameters according to prior expert knowledge as well as observational data about model parameter values. This study assesses the performance of two popular uncertainty analysis (UA) techniques, i.e., generalized likelihood uncertainty estimation (GLUE) and Bayesian method implemented with the Markov chain Monte Carlo sampling algorithm, in evaluating model parameter uncertainty in flood simulations. These two methods were applied to the semi-distributed Topographic hydrologic model (TOPMODEL) that includes five parameters. A case study was carried out for a small humid catchment in the southeastern China. The performance assessment of the GLUE and Bayesian methods were conducted with advanced tools suited for probabilistic simulations of continuous variables such as streamflow. Graphical tools and scalar metrics were used to test several attributes of the simulation quality of selected flood events: deterministic accuracy and the accuracy of 95 % prediction probability uncertainty band (95PPU). Sensitivity analysis was conducted to identify sensitive parameters that largely affect the model output results. Subsequently, the GLUE and Bayesian methods were used to analyze the uncertainty of sensitive parameters and further to produce their posterior distributions. Based on their posterior parameter samples, TOPMODEL’s simulations and the corresponding UA results were conducted. Results show that the form of exponential decline in conductivity and the overland flow routing velocity were sensitive parameters in TOPMODEL in our case. Small changes in these two parameters would lead to large differences in flood simulation results. Results also suggest that, for both UA techniques, most of streamflow observations were bracketed by 95PPU with the containing ratio value larger than 80 %. In comparison, GLUE gave narrower prediction uncertainty bands than the Bayesian method. It was found that the mode estimates of parameter posterior distributions are suitable to result in better performance of deterministic outputs than the 50 % percentiles for both the GLUE and Bayesian analyses. In addition, the simulation results calibrated with Rosenbrock optimization algorithm show a better agreement with the observations than the UA’s 50 % percentiles but slightly worse than the hydrographs from the mode estimates. The results clearly emphasize the importance of using model uncertainty diagnostic approaches in flood simulations.     

On the criteria of model performance evaluation for real-time flood forecasting

Model performance evaluation for real-time flood forecasting has been conducted using various criteria. Although the coefficient of efficiency (CE) is most widely used, we demonstrate that a model achieving good model efficiency may actually be inferior to the naïve (or persistence) forecasting, if the flow series has a high lag-1 autocorrelation coefficient. We derived sample-dependent and AR model-dependent asymptotic relationships between the coefficient of efficiency and the coefficient of persistence (CP) which form the basis of a proposed CE–CP coupled model performance evaluation criterion. Considering the flow persistence and the model simplicity, the AR(2) model is suggested to be the benchmark model for performance evaluation of real-time flood forecasting models. We emphasize that performance evaluation of flood forecasting models using the proposed CE–CP coupled criterion should be carried out with respect to individual flood events. A single CE or CP value derived from a multi-event artifactual series by no means provides a multi-event overall evaluation and may actually disguise the real capability of the proposed model.     

A big-data spatial,temporal and network analysis of bovine tuberculosis between wildlife (badgers) and cattle

Bovine tuberculosis (TB) poses a serious threat for agricultural industry in several countries, it involves potential interactions between wildlife and cattle and creates societal problems in terms of human-wildlife conflict. This study addresses connectedness network analysis, the spatial, and temporal dynamics of TB between cattle in farms and the European badger (Meles meles) using a large dataset generated by a calibrated agent based model. Results showed that infected network connectedness was lower in badgers than in cattle. The contribution of an infected individual to the mean distance of disease spread over time was considerably lower for badger than cattle; badgers mainly spread the disease locally while cattle infected both locally and across longer distances. The majority of badger-induced infections occurred when individual badgers leave their home sett, and this was positively correlated with badger population growth rates. Point pattern analysis indicated aggregation in the spatial pattern of TB prevalence in badger setts across all scales. The spatial distribution of farms that were not TB free was aggregated at different scales than the spatial distribution of infected badgers and became random at larger scales. The spatial cross correlation between infected badger setts and infected farms revealed that generally infected setts and farms do not coexist except at few scales. Temporal autocorrelation detected a two year infection cycle for badgers, while there was both within the year and longer cycles for infected cattle. Temporal cross correlation indicated that infection cycles in badgers and cattle are negatively correlated. The implications of these results for understanding the dynamics of the disease are discussed.     

Impact of sensor measurement error on sensor positioning in water quality monitoring networks

This paper studies the impact of sensor measurement error on designing a water quality monitoring network for a river system, and shows that robust sensor locations can be obtained when an optimization algorithm is combined with a statistical process control (SPC) method. Specifically, we develop a possible probabilistic model of sensor measurement error and the measurement error model is embedded into a simulation model of a river system. An optimization algorithm is used to find the optimal sensor locations that minimize the expected time until a spill detection in the presence of a constraint on the probability of detecting a spill. The experimental results show that the optimal sensor locations are highly sensitive to the variability of measurement error and false alarm rates are often unacceptably high. An SPC method is useful in finding thresholds that guarantee a false alarm rate no more than a pre-specified target level, and an optimization algorithm combined with the thresholds finds a robust sensor network.     

Simulating and mapping the risk of surge floods in multiple typhoon scenarios: a case study of Yuhuan County,Zhejiang Province,China

We propose a scenario-based method for simulating and mapping the risk of surge floods for use by local authorities concerned with public safety and urban planning in coastal areas. Focusing on the triad of hazard, vulnerability and adaptation capability, we estimate the comprehensive risk and display its spatial distribution using the raster calculation tool in ArcGIS. The detailed methodology is introduced via a case study of Yuhuan, an island county in Zhejiang Province, China, which is frequently affected by typhoon storm surges. First, we designed 24 typhoon scenarios and modeled the flood process in each scenario using the hydrodynamic module of MIKE 21. Second, flood depth and area were used for hazard assessment; an authorized indicator system of land use categories and a survey of emergency shelters were used for vulnerability and adaptation capability assessment, respectively; and a quantified model was used for assessment of the comprehensive risk. Lastly, we used the GIS raster calculation tool for mapping the risk of storm surges in multiple typhoon scenarios. Our principal findings are as follows: (1) Seawalls are more likely to be overtopped or destroyed by more severe storm surges with increasing typhoon intensity. (2) Most of the residential areas with inadequate emergency shelters are highly vulnerable to flood events. (3) As projected in the risk mapping, if an exceptional typhoon with a central pressure of 915 or 925 hPa made a landfall in Yuhuan, a wide range of areas would be flooded and at high risk. (4) Determining optimal strategies based on identification of risk-inducing factors is the most effective way of promoting safe and sustainable development in coastal cities.     

Intercomparison of multiple statistical downscaling methods: multi-criteria model selection for South Korea

A number of statistical downscaling methodologies have been introduced to bridge the gap in scale between outputs of climate models and climate information needed to assess potential impacts at local and regional scales. Four statistical downscaling methods [bias-correction/spatial disaggregation (BCSD), bias-correction/constructed analogue (BCCA), multivariate adaptive constructed analogs (MACA), and bias-correction/climate imprint (BCCI)] are applied to downscale the latest climate forecast system reanalysis (CFSR) data to stations for precipitation, maximum temperature, and minimum temperature over South Korea. All methods are calibrated with observational station data for 19 years from 1973 to 1991 and validated for the more recent 19-year period from 1992 to 2010. We construct a comprehensive suite of performance metrics to inter-compare methods, which is comprised of five criteria related to time-series, distribution, multi-day persistence, extremes, and spatial structure. Based on the performance metrics, we employ technique for order of preference by similarity to ideal solution (TOPSIS) and apply 10,000 different weighting combinations to the criteria of performance metrics to identify a robust statistical downscaling method and important criteria. The results show that MACA and BCSD have comparable skill in the time-series related criterion and BCSD outperforms other methods in distribution and extremes related criteria. In addition, MACA and BCCA, which incorporate spatial patterns, show higher skill in the multi-day persistence criterion for temperature, while BCSD shows the highest skill for precipitation. For the spatial structure related criterion, BCCA and MACA outperformed BCSD and BCCI. From the TOPSIS analysis, we found that MACA is the most robust method for all variables in South Korea, and BCCA and BCSD are the second for temperature and precipitation, respectively. We also found that the contribution of the multi-day persistence and spatial structure related criteria are crucial to ranking the skill of statistical downscaling methods.     

Let-It-Rain: a web application for stochastic point rainfall generation at ungaged basins and its applicability in runoff and flood modeling

We present a web application named Let-It-Rain that is able to generate a 1-h temporal resolution synthetic rainfall time series using the modified Bartlett–Lewis rectangular pulse (MBLRP) model, a type of Poisson stochastic rainfall generator. Let-It-Rain, which can be accessed through the web address http://www.LetItRain.info, adopts a web-based framework combining ArcGIS Server from server side for parameter value dissemination and JavaScript from client side to implement the MBLRP model. This enables any desktop and mobile end users with internet access and web browser to obtain the synthetic rainfall time series at any given location at which the parameter regionalization work has been completed (currently the contiguous United States and Republic of Korea) with only a few mouse clicks. Let-It-Rain shows satisfactory performance in its ability to reproduce observed rainfall mean, variance, auto-correlation, and probability of zero rainfall at hourly through daily accumulation levels. It also shows a reasonably good performance in reproducing watershed runoff depth and peak flow. We expect that Let-It-Rain can stimulate the uncertainty analysis of hydrologic variables across the world.     

Improving at-site flood frequency analysis with additional spatial information: a probabilistic regional envelope curve approach

Extreme flood events have detrimental effects on society, the economy and the environment. Widespread flooding across South East Queensland in 2011 and 2013 resulted in the loss of lives and significant cost to the economy. In this region, flood risk planning and the use of traditional flood frequency analysis (FFA) to estimate both the magnitude and frequency of the 1-in-100 year flood is severely limited by short gauging station records. On average, these records are 42 years in Eastern Australia and many have a poor representation of extreme flood events. The major aim of this study is to test the application of an alternative method to estimate flood frequency in the form of the Probabilistic Regional Envelope Curve (PREC) approach which integrates additional spatial information of extreme flood events. In order to better define and constrain a working definition of an extreme flood, an Australian Envelope Curve is also produced from available gauging station data. Results indicate that the PREC method shows significant changes to the larger recurrence intervals (≥100 years) in gauges with either too few, or too many, extreme flood events. A decision making process is provided to ascertain when this method is preferable for FFA.     

An assessment of using subsampling method in selection of a flood frequency distribution

In flood frequency analysis, a suitable probability distribution function is required in order to establish the flood magnitude-return period relationship. Goodness of fit (GOF) techniques are often employed to select a suitable distribution function in this context. But they have been often criticized for their inability to discriminate between statistical distributions for the same application. This paper investigates the potential utility of subsampling, a resampling technique with the aid of a GOF test to select the best distribution for frequency analysis. The performance of the methodology is assessed by applying the methodology to observed and simulated annual maximum (AM) discharge data series. Several AM series of different record lengths are used as case studies to determine the performance. Numerical analyses are carried out to assess the performance in terms of sample size, subsample size and statistical properties inherent in the AM data series. The proposed methodology is also compared with the standard Anderson–Darling (AD) test. It is found that the methodology is suitable for a longer data series. A better performance is obtained when the subsample size is taken around half of the underlying data sample. The methodology has also outperformed the standard AD test in terms of effectively discriminating between distributions. Overall, all results point that the subsampling technique can be a promising tool in discriminating between distributions.