Multi-criteria evaluation of CMIP5 GCMs for climate change impact analysis

Climate change is expected to have severe impacts on global hydrological cycle along with food-water-energy nexus. Currently, there are many climate models used in predicting important climatic variables. Though there have been advances in the field, there are still many problems to be resolved related to reliability, uncertainty, and computing needs, among many others. In the present work, we have analyzed performance of 20 different global climate models (GCMs) from Climate Model Intercomparison Project Phase 5 (CMIP5) dataset over the Columbia River Basin (CRB) in the Pacific Northwest USA. We demonstrate a statistical multicriteria approach, using univariate and multivariate techniques, for selecting suitable GCMs to be used for climate change impact analysis in the region. Univariate methods includes mean, standard deviation, coefficient of variation, relative change (variability), Mann-Kendall test, and Kolmogorov-Smirnov test (KS-test); whereas multivariate methods used were principal component analysis (PCA), singular value decomposition (SVD), canonical correlation analysis (CCA), and cluster analysis. The analysis is performed on raw GCM data, i.e., before bias correction, for precipitation and temperature climatic variables for all the 20 models to capture the reliability and nature of the particular model at regional scale. The analysis is based on spatially averaged datasets of GCMs and observation for the period of 1970 to 2000. Ranking is provided to each of the GCMs based on the performance evaluated against gridded observational data on various temporal scales (daily, monthly, and seasonal). Results have provided insight into each of the methods and various statistical properties addressed by them employed in ranking GCMs. Further; evaluation was also performed for raw GCM simulations against different sets of gridded observational dataset in the area.     

Revisiting hydrological drought propagation and recovery considering water quantity and quality

Climate extremes, in particular droughts, are significant driving forces towards riverine ecosystem disturbance. Drought impacts on stream ecosystems include losses that can be either direct (e.g., destruction of habitat for aquatic species) or indirect (e.g., deterioration of water quality, soil quality, and increased chance of wildfires). This paper combines hydrologic drought and water quality changes during droughts and represents a multistage framework to detect and characterize hydrological droughts while considering water quality parameters. This method is applied to 52 streamflow stations in the state of California, USA, over the study period of 1950–2010. The framework is assessed and validated based on two drought events declared by the state in 2002 and 2008. Results show that there are two opposite drought propagation patterns in northern and southern California. In general, northern California indicates more frequent droughts with shorter time to recover. Chronology of drought shows that stations located in southern California have not followed a specific pattern but they experienced longer drought episodes with prolonged drought recovery. When considering water quality, results show that droughts either deteriorate or enhance water systems, depending on the parameter of interest. Undesirable changes (e.g., increased temperature and decreased dissolved oxygen) are observed during droughts. In contrast, decreased turbidity is detected in rivers during drought episodes, which is desirable in water systems. Nevertheless, water quality deteriorates during drought recovery, even after drought termination. Depending on climatic and streamflow characteristics of the watersheds, it was found that it would take nearly 2 months on average for water quality to recover after drought termination.     

Radiance data assimilation for operational snow and streamflow forecasting

Estimation of seasonal snowpack, in mountainous regions, is crucial for accurate streamflow prediction. This paper examines the ability of data assimilation (DA) of remotely sensed microwave radiance data to improve snow water equivalent prediction, and ultimately operational streamflow forecasts. Operational streamflow forecasts in the National Weather Service River Forecast Center (NWSRFC) are produced with a coupled SNOW17 (snow model) and SACramento Soil Moisture Accounting (SAC-SMA) model. A comparison of two assimilation techniques, the ensemble Kalman filter (EnKF) and the particle filter (PF), is made using a coupled SNOW17 and the microwave emission model for layered snow pack (MEMLS) model to assimilate microwave radiance data. Microwave radiance data, in the form of brightness temperature (TB), is gathered from the advanced microwave scanning radiometer-earth observing system (AMSR-E) at the 36.5 GHz channel. SWE prediction is validated in a synthetic experiment. The distribution of snowmelt from an experiment with real data is then used to run the SAC-SMA model. Several scenarios on state or joint state-parameter updating with TB data assimilation to SNOW-17 and SAC-SMA models were analyzed, and the results show potential benefit for operational streamflow forecasting.     

Assessing the uncertainties of hydrologic model selection in climate change impact studies

The uncertainties associated with atmosphere‐ocean General Circulation Models (GCMs) and hydrologic models are assessed by means of multi‐modelling and using the statistically downscaled outputs from eight GCM simulations and two emission scenarios. The statistically downscaled atmospheric forcing is used to drive four hydrologic models, three lumped and one distributed, of differing complexity: the Sacramento Soil Moisture Accounting (SAC‐SMA) model, Conceptual HYdrologic MODel (HYMOD), Thornthwaite‐Mather model (TM) and the Precipitation Runoff Modelling System (PRMS). The models are calibrated based on three objective functions to create more plausible models for the study. The hydrologic model simulations are then combined using the Bayesian Model Averaging (BMA) method according to the performance of each models in the observed period, and the total variance of the models. The study is conducted over the rainfall‐dominated Tualatin River Basin (TRB) in Oregon, USA. This study shows that the hydrologic model uncertainty is considerably smaller than GCM uncertainty, except during the dry season, suggesting that the hydrologic model selection‐combination is critical when assessing the hydrologic climate change impact. The implementation of the BMA in analysing the ensemble results is found to be useful in integrating the projected runoff estimations from different models, while enabling to assess the model structural uncertainty. Copyright © 2011 John Wiley & Sons, Ltd.     

Snow water equivalent prediction using Bayesian data assimilation methods

Using the U.S. National Weather Service’s SNOW-17 model, this study compares common sequential data assimilation methods, the ensemble Kalman filter (EnKF), the ensemble square root filter (EnSRF), and four variants of the particle filter (PF), to predict seasonal snow water equivalent (SWE) within a small watershed near Lake Tahoe, California. In addition to SWE estimation, the various data assimilation methods are used to estimate five of the most sensitive parameters of SNOW-17 by allowing them to evolve with the dynamical system. Unlike Kalman filters, particle filters do not require Gaussian assumptions for the posterior distribution of the state variables. However, the likelihood function used to scale particle weights is often assumed to be Gaussian. This study evaluates the use of an empirical cumulative distribution function (ECDF) based on the Kaplan–Meier survival probability method to compute particle weights. These weights are then used in different particle filter resampling schemes. Detailed analyses are conducted for synthetic and real data assimilation and an assessment of the procedures is made. The results suggest that the particle filter, especially the empirical likelihood variant, is superior to the ensemble Kalman filter based methods for predicting model states, as well as model parameters.     

A regional Bayesian hierarchical model for flood frequency analysis

Stochastic Environmental Research and Risk Assessment - In this study, we propose a regional Bayesian hierarchical model for flood frequency analysis. The Bayesian method is an alternative to the...     

Towards improved post‐processing of hydrologic forecast ensembles

Forecast ensembles of hydrological and hydrometeorologial variables are prone to various uncertainties arising from climatology, model structure and parameters, and initial conditions at the forecast date. Post‐processing methods are usually applied to adjust the mean and variance of the ensemble without any knowledge about the uncertainty sources. This study initially addresses the drawbacks of a commonly used statistical technique, quantile mapping (QM), in bias correction of hydrologic forecasts. Then, an auxiliary variable, the failure index (γ), is proposed to estimate the ineffectiveness of the post‐processing method based on the agreement of adjusted forecasts with corresponding observations during an analysis period prior to the forecast date. An alternative post‐processor based on copula functions is then introduced such that marginal distributions of observations and model simulations are combined to create a multivariate joint distribution. A set of 2500 hypothetical forecast ensembles with parametric marginal distributions of simulated and observed variables are post‐processed with both QM and the proposed multivariate post‐processor. Deterministic forecast skills show that the proposed copula‐based post‐processing is more effective than the QM method in improving the forecasts. It is found that the performance of QM is highly correlated with the failure index, unlike the multivariate post‐processor. In probabilistic metrics, the proposed multivariate post‐processor generally outperforms QM. Further evaluation of techniques is conducted for river flow forecast of Sprague River basin in southern Oregon. Results show that the multivariate post‐processor performs better than the QM technique; it reduces the ensemble spread and is a more reliable approach for improving the forecast. Copyright © 2012 John Wiley & Sons, Ltd.     

Analysis of precipitation extremes with the assessment of regional climate models over the Willamette River Basin,USA

An appropriate, rapid and effective response to extreme precipitation and any potential flood disaster is essential. Providing an accurate estimate of future changes to such extreme events due to climate change are crucial for responsible decision making in flood risk management given the predictive uncertainties. The objective of this article is to provide a comparison of dynamically downscaled climate models simulations from multiple model including 12 different combinations of General Circulation Model (GCM)–regional climate model (RCM), which offers an abundance of additional data sets. The three major aspects of this study include the bias correction of RCM scenarios, the application of a newly developed performance metric and the extreme value analysis of future precipitation. The dynamically downscaled data sets reveal a positive overall bias that is removed through quantile mapping bias correction method. The added value index was calculated to evaluate the models' simulations. Results from this metric reveal that not all of the RCMs outperform their host GCMs in terms of correlation skill. Extreme value theory was applied to both historic, 1980–1998, and future, 2038–2069, daily data sets to provide estimates of changes to 2‐ and 25‐year return level precipitation events. The generalized Pareto distribution was used for this purpose. The Willamette River basin was selected as the study region for analysis because of its topographical variability and tendency for significant precipitation. The extreme value analysis results showed significant differences between model runs for both historical and future periods with considerable spatial variability in precipitation extremes. Copyright © 2012 John Wiley & Sons, Ltd.