Multi-models for SPI drought forecasting in the north of Haihe River Basin,China

Drought is one of the most devastating climate disasters. Hence, drought forecasting plays an important role in mitigating some of the adverse effects of drought. Data-driven models are widely used for drought forecasting such as ARIMA model, artificial neural network (ANN) model, wavelet neural network (WANN) model, support vector regression model, grey model and so on. Three data-driven models (ARIMA model; ANN model; WANN model) are used in this study for drought forecasting based on standard precipitation index of two time scales (SPI; SPI-6 and SPI-12). The optimal data-driven model and time scale of SPI are then selected for effective drought forecasting in the North of Haihe River Basin. The effectiveness of the three data-models is compared by Kolmogorov–Smirnov (K–S) test, Kendall rank correlation, and the correlation coefficients (R²). The forecast results shows that the WANN model is more suitable and effective for forecasting SPI-6 and SPI-12 values in the north of Haihe River Basin.     

<Emphasis Type="Italic">Clostridium difficile</Emphasis> infection incidence prediction in hospitals (CDIIPH): a predictive model based on decision tree and fuzzy techniques

Several risk factors associated with the increased likelihood of healthcare-associated Clostridium difficile infection (CDI) have been identified in the literature. These risk factors are mainly related to age, previous CDI, antimicrobial exposure, and prior hospitalization. No model is available in the published literature that can be used to predict the CDI incidence using healthcare administration data. However, the administrative data can be imprecise and may challenge the building of classical statistical models. Fuzzy set theory can deal with the imprecision inherent in such data. This research aimed to develop a model based on deterministic and fuzzy mathematical techniques for the prediction of hospital-associated CDI by using the explanatory variables controllable by hospitals and health authority administration. Retrospective data on CDI incidence and other administrative data obtained from 22 hospitals within a regional health authority in British Columbia were used to develop a decision tree (deterministic technique based) and a fuzzy synthetic evaluation model (fuzzy technique based). The decision tree model had a higher prediction accuracy than that of the fuzzy based model. However, among the common results predicted by two models, 72 % were correct. Therefore, this relationship was used to combine their results to increase the precision and the strength of evidence of the prediction. These models were further used to develop an Excel-based tool called C. difficile Infection Incidence Prediction in Hospitals (CDIIPH). The tool can be utilized by health authorities and hospitals to predict the magnitude of CDI incidence in the following quarter.     

Fast and direct nonparametric procedures in the L-moment homogeneity test

Regional frequency analysis is an important tool to properly estimate hydrological characteristics at ungauged or partially gauged sites in order to prevent hydrological disasters. The delineation of homogeneous groups of sites is an important first step in order to transfer information and obtain accurate quantile estimates at the target site. The Hosking–Wallis homogeneity test is usually used to test the homogeneity of the selected sites. Despite its usefulness and good power, it presents some drawbacks including the subjective choice of a parametric distribution for the data and a poorly justified rejection threshold. The present paper addresses these drawbacks by integrating nonparametric procedures in the L-moment homogeneity test. To assess the rejection threshold, three resampling methods (permutation, bootstrap and Pólya resampling) are considered. Results indicate that permutation and bootstrap methods perform better than the parametric Hosking–Wallis test in terms of power as well as in time and procedure simplicity. A real-world case study shows that the nonparametric tests agree with the HW test concerning the homogeneity of the volume and the bivariate case while they disagree for the peak case, but that the assumptions of the HW test are not well respected.     

Regularized Bayesian estimation for GEV-B-splines model

Large observed datasets are not stationary and/or depend on covariates, especially, in the case of extreme hydrometeorological variables. This causes the difficulty in estimation, using classical hydrological frequency analysis. A number of non-stationary models have been developed using linear or quadratic polynomial functions or B-splines functions to estimate the relationship between parameters and covariates. In this article, we propose regularised generalized extreme value model with B-splines (GEV-B-splines models) in a Bayesian framework to estimate quantiles. Regularisation is based on penalty and aims to favour parsimonious model especially in the case of large dimension space. Penalties are introduced in a Bayesian framework and the corresponding priors are detailed. Five penalties are considered and the corresponding priors are developed for comparison purpose as: Least absolute shrinkage and selection (Lasso and Ridge) and smoothing clipped absolute deviations (SCAD) methods (SCAD1, SCAD2 and SCAD3). Markov chain Monte Carlo (MCMC) algorithms have been developed for each model to estimate quantiles and their posterior distributions. Those approaches are tested and illustrated using simulated data with different sample sizes. A first simulation was made on polynomial B-splines functions in order to choose the most efficient model in terms of relative mean biais (RMB) and the relative mean-error (RME) criteria. A second simulation was performed with the SCAD1 penalty for sinusoidal dependence to illustrate the flexibility of the proposed approach. Results show clearly that the regularized approaches leads to a significant reduction of the bias and the mean square error, especially for small sample sizes (n < 100). A case study has been considered to model annual peak flows at Fort-Kent catchment with the total annual precipitations as covariates. The conditional quantile curves were given for the regularized and the maximum likelihood methods.     

Seasonal copula models for the analysis of glacier discharge at King George Island,Antarctica

Modelling glacier discharge is an important issue in hydrology and climate research. Glaciers represent a fundamental water resource when melting of ice and snow contributes to runoff. Glaciers are also studied as natural global warming sensors. GLACKMA association has implemented one of their Pilot Experimental Catchment areas at the King George Island in the Antarctica which records values of the liquid discharge from Collins glacier. In this paper, we propose the use of time-varying copula models for analyzing the relationship between air temperature and glacier discharge, which is clearly non constant and non linear through time. A seasonal copula model is defined where both the marginal and copula parameters vary periodically along time following a seasonal dynamic. Full Bayesian inference is performed such that the marginal and copula parameters are estimated in a one single step, in contrast with the usual two-step approach. Bayesian prediction and model selection is also carried out for the proposed model such that Bayesian credible intervals can be obtained for the conditional glacier discharge given a value of the temperature at any given time point. The proposed methodology is illustrated using the GLACKMA real data where there is, in addition, a hydrological year of missing discharge data which were not possible to measure accurately due to problems in the sounding.     

Assessment of the climate change impacts on flood frequency (case study: Bazoft Basin,Iran)

The present study attempts to investigate potential impacts of climate change on floods frequency in Bazoft Basin which is located in central part of Iran. A combination of four general circulation models is used through a weighting approach to assess uncertainty in the climate projections. LARS-WG model is applied to downscale large scale atmospheric data to local stations. The resulting data is in turn used as input of the hydrological model Water and Energy Transfer between Soil, plants and atmosphere (WetSpa) to simulate runoff for present (1971–2000), near future (2020–2049) and far future (2071–2100) conditions. Results demonstrate good performance of both WetSpa and LARS-WG models. In addition in this paper instantaneous peak flow (IPF) is estimated using some empirical equations including Fuller, Sangal and Fill–Steiner methods. Comparison of estimated and observed IPF shows that Fill–Steiner is better than other methods. Then different probability distribution functions are fit to IPF series. Results of flood frequency analysis indicate that Pearson III is the best distribution fitted to IPF data. It is also indicated that flood magnitude will decrease in future for all return periods.     

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.