A spatial temporal downscaling approach to development of IDF relations for Perth airport region in the context of climate change

ABSTRACT

Downscaling of climate projections is the most adapted method to assess the impacts of climate change at regional and local scales. This study utilized both spatial and temporal downscaling approaches to develop intensity–duration–frequency (IDF) relations for sub-daily rainfall extremes in the Perth airport area. A multiple regression-based statistical downscaling model tool was used for spatial downscaling of daily rainfall using general circulation models (GCMs) (Hadley Centre’s GCM and Canadian Global Climate Model) climate variables. A simple scaling regime was identified for 30 minutes to 24 hours duration of observed annual maximum (AM) rainfall. Then, statistical properties of sub-daily AM rainfall were estimated by scaling an invariant model based on the generalized extreme value distribution. RMSE, Nash-Sutcliffe efficiency coefficient and percentage bias values were estimated to check the accuracy of downscaled sub-daily rainfall. This proved the capability of the proposed approach in developing a linkage between large-scale GCM daily variables and extreme sub-daily rainfall events at a given location. Finally IDF curves were developed for future periods, which show similar extreme rainfall decreasing trends for the 2020s, 2050s and 2080s for both GCMs.

Editor M.C. Acreman; Associate editor S. Kanae     

Assessing regression‐based statistical approaches for downscaling precipitation over North America

This paper assesses linear regression‐based methods in downscaling daily precipitation from the general circulation model (GCM) scale to a regional climate model (RCM) scale (45‐ and 15‐km grids) and down to a station scale across North America. Traditional downscaling experiments (linking reanalysis/dynamical model predictors to station precipitation) as well as nontraditional experiments such as predicting dynamic model precipitation from larger‐scale dynamic model predictors or downscaling dynamic model precipitation from predictors at the same scale are conducted. The latter experiments were performed to address predictability limit and scale issues. The results showed that the downscaling of daily precipitation occurrence was rarely successful at all scales, although results did constantly improve with the increased resolution of climate models. The explained variances for downscaled precipitation amounts at the station scales were low, and they became progressively better when using predictors from a higher‐resolution climate model, thus showing a clear advantage in using predictors from RCMs driven by reanalysis at its boundaries, instead of directly using reanalysis data. The low percentage of explained variances resulted in considerable underestimation of daily precipitation mean and standard deviation. Although downscaling GCM precipitation from GCM predictors (or RCM precipitation from RCM predictors) cannot really be considered downscaling, as there is no change in scale, the exercise yields interesting information as to the limit in predictive ability at the station scale. This was especially clear at the GCM scale, where the inability of downscaling GCM precipitation from GCM predictors demonstrates that GCM precipitation‐generating processes are largely at the subgrid scale (especially so for convective events), thus indicating that downscaling precipitation at the station scale from GCM scale is unlikely to be successful. Although results became better at the RCM scale, the results indicate that, overall, regression‐based approaches did not perform well in downscaling precipitation over North America. Copyright © 2013 John Wiley & Sons, Ltd.     

Extraction of information content from stochastic disaggregation and bias corrected downscaled precipitation variables for crop simulation

We applied a simple statistical downscaling procedure for transforming daily global climate model (GCM) rainfall to the scale of an agricultural experimental station in Katumani, Kenya. The transformation made was two-fold. First, we corrected the rainfall frequency bias of the climate model by truncating its daily rainfall cumulative distribution into the station’s distribution based on a prescribed observed wet-day threshold. Then, we corrected the climate model rainfall intensity bias by mapping its truncated rainfall distribution into the station’s truncated distribution. Further improvements were made to the bias corrected GCM rainfall by linking it with a stochastic disaggregation scheme to correct the time structure problem inherent with daily GCM rainfall. Results of the simple and hybridized GCM downscaled precipitation variables (total, probability of occurrence, intensity and dry spell length) were linked with a crop model for a more objective evaluation of their performance using a non-linear measure based on mutual information based on entropy. This study is useful for the identification of both suitable downscaling technique as well as the effective precipitation variables for forecasting crop yields using GCM’s outputs which can be useful for addressing food security problems beforehand in critical basins around the world.     

Estimation of catchment yield and associated uncertainties due to climate change in a mountainous catchment in Australia

This paper examines the impacts of climate change on future water yield with associated uncertainties in a mountainous catchment in Australia using a multi‐model approach based on four global climate models (GCMs), 200 realisations (50 realisations from each GCM) of downscaled rainfalls, 2 hydrological models and 6 sets of model parameters. The ensemble projections by the GCMs showed that the mean annual rainfall is likely to reduce in the future decades by 2–5% in comparison with the current climate (1987–2012). The results of ensemble runoff projections indicated that the mean annual runoff would reduce in future decades by 35%. However, considerable uncertainty in the runoff estimates was found as the ensemble results project changes of the 5^th (dry scenario) and 95^th (wet scenario) percentiles by ?73% to +27%, ?73% to +12%, ?77% to +21% and ?80% to +24% in the decades of 2021–2030, 2031–2040, 2061–2070 and 2071–2080, respectively. Results of uncertainty estimation demonstrated that the choice of GCMs dominates overall uncertainty. Realisation uncertainty (arising from repetitive simulations for a given time step during downscaling of the GCM data to catchment scale) of the downscaled rainfall data was also found to be remarkably high. Uncertainty linked to the choice of hydrological models was found to be quite small in comparison with the GCM and realisation uncertainty. The hydrological model parameter uncertainty was found to be lowest among the sources of uncertainties considered in this study. Copyright © 2015 John Wiley & Sons, Ltd.     

Regional precipitation and temperature scenarios for climate change

Abstract

To investigate the consequences of climate change on the water budget in small catchments, it is necessary to know the change of local precipitation and temperature. General Circulation Models (GCM) cannot provide regional climate parameters yet, because of their coarse resolution and imprecise modelling of precipitation. Therefore downscaling of precipitation and temperature has to be carried out from the GCM grids to a small scale of a few square kilometres. Daily rainfall and temperature are modelled as processes conditioned on atmospheric circulation. Rainfall is linked to the circulation patterns (CPs) using conditional probabilities and conditional rainfall amount distribution. Both temperature and precipitation are downscaled to several locations simultaneously taking into account the CP dependent spatial correlation. Temperature is modelled using a simple autoregressive approach, conditioned on atmospheric circulation and local areal precipitation. The model uses the classification scheme of the German Weather Service and a fuzzy rule-based classification. It was applied in the Aller catchment for validation using observed rainfall and temperature, and observed classified geopotential pressure heights. GCM scenarios of the ECHAM model were used to make climate change predictions (using classified GCM geopotential heights); simulated values agree fairly well with historical data. Results for different GCM scenarios are shown.     

Propagation of uncertainty from rainfall to runoff: A case study with a stochastic rainfall generator

We explore the impact of uncertainties in the spatial–temporal distribution of rainfall on the prediction of peak discharge in a typical mountain basin. To this end, we use a stochastic generator previously developed for rainfall downscaling, and we estimate the basin response by adopting a semi-distributed hydrological model. The results of the analysis provide information on the minimum rainfall resolution needed for operational flood forecasting, and confirm the sensitivity of peak discharge estimates to errors in the determination of the power spectrum of the precipitation field.     

Hydrologic drought prediction under climate change: Uncertainty modeling with Dempster–Shafer and Bayesian approaches

Representation and quantification of uncertainty in climate change impact studies are a difficult task. Several sources of uncertainty arise in studies of hydrologic impacts of climate change, such as those due to choice of general circulation models (GCMs), scenarios and downscaling methods. Recently, much work has focused on uncertainty quantification and modeling in regional climate change impacts. In this paper, an uncertainty modeling framework is evaluated, which uses a generalized uncertainty measure to combine GCM, scenario and downscaling uncertainties. The Dempster–Shafer (D–S) evidence theory is used for representing and combining uncertainty from various sources. A significant advantage of the D–S framework over the traditional probabilistic approach is that it allows for the allocation of a probability mass to sets or intervals, and can hence handle both aleatory or stochastic uncertainty, and epistemic or subjective uncertainty. This paper shows how the D–S theory can be used to represent beliefs in some hypotheses such as hydrologic drought or wet conditions, describe uncertainty and ignorance in the system, and give a quantitative measurement of belief and plausibility in results. The D–S approach has been used in this work for information synthesis using various evidence combination rules having different conflict modeling approaches. A case study is presented for hydrologic drought prediction using downscaled streamflow in the Mahanadi River at Hirakud in Orissa, India. Projections of n most likely monsoon streamflow sequences are obtained from a conditional random field (CRF) downscaling model, using an ensemble of three GCMs for three scenarios, which are converted to monsoon standardized streamflow index (SSFI-4) series. This range is used to specify the basic probability assignment (bpa) for a Dempster–Shafer structure, which represents uncertainty associated with each of the SSFI-4 classifications. These uncertainties are then combined across GCMs and scenarios using various evidence combination rules given by the D–S theory. A Bayesian approach is also presented for this case study, which models the uncertainty in projected frequencies of SSFI-4 classifications by deriving a posterior distribution for the frequency of each classification, using an ensemble of GCMs and scenarios. Results from the D–S and Bayesian approaches are compared, and relative merits of each approach are discussed. Both approaches show an increasing probability of extreme, severe and moderate droughts and decreasing probability of normal and wet conditions in Orissa as a result of climate change.     

Distributed modelling of future changes in hydrological processes of Spencer Creek watershed

The hydrologic impact of climate change has been largely assessed using mostly conceptual hydrologic models. This study investigates the use of distributed hydrologic model for the assessment of the climate change impact for the Spencer Creek watershed in Southern Ontario (Canada). A coupled MIKE SHE/MIKE 11 hydrologic model is developed to represent the complex hydrologic conditions in the Spencer Creek watershed, and later to simulate climate change impact using Canadian global climate model (CGCM 3·1) simulations. Owing to the coarse resolution of GCM data (daily GCM outputs), statistical downscaling techniques are used to generate higher resolution data (daily precipitation and temperature series). The modelling results show that the coupled model captured the snow storage well and also provided good simulation of evapotranspiration (ET) and groundwater recharge. The simulated streamflows are consistent with the observed flows at different sites within the catchment. Using a conservative climate change scenario, the downscaled GCM scenarios predicted an approximately 14–17% increase in the annual mean precipitation and 2–3 °C increase in annual mean maximum and minimum temperatures for the 2050s (i.e., 2046–2065). When the downscaled GCM scenarios were used in the coupled model, the model predicted a 1–5% annual decrease in snow storage for 2050s, approximately 1–10% increase in annual ET, and a 0·5–6% decrease in the annual groundwater recharge. These results are consistent with the downscaled temperature results. For future streamflows, the coupled model indicated an approximately 10–25% increase in annual streamflows for all sites, which is consistent with the predicted changes in precipitation. Overall, it is shown that distributed hydrologic modelling can provide useful information not only about future changes in streamflow but also changes in other key hydrologic processes such as snow storage, ET, and groundwater recharge, which can be particularly important depending on the climatic region of concern. The study results indicate that the coupled MIKE SHE/MIKE 11 hydrologic model could be a particularly useful tool for understanding the integrated effect of climate change in complex catchment scale hydrology. Copyright © 2010 John Wiley & Sons, Ltd.     

Impact of the spatial variability of daily precipitation on hydrological projections: A comparison of GCM‐ and RCM‐driven cases in the Han River basin,Korea

In this study, we investigate the impact of the spatial variability of daily precipitation on hydrological projections based on a comparative assessment of streamflow simulations driven by a global climate model (GCM) and two regional climate models (RCMs). A total of 12 different climate input datasets, that is, the raw and bias‐corrected GCM and raw and bias‐corrected two RCMs for the reference and future periods, are fed to a semidistributed hydrological model to assess whether the bias correction using quantile mapping and dynamical downscaling using RCMs can improve streamflow simulation in the Han River basin, Korea. A statistical analysis of the daily precipitation demonstrates that the precipitation simulated by the GCM fails to capture the large variability of the observed daily precipitation, in which the spatial autocorrelation decreases sharply within a relatively short distance. However, the spatial variability of precipitation simulated by the two RCMs shows better agreement with the observations. After applying bias correction to the raw GCM and raw RCMs outputs, only a slight change is observed in the spatial variability, whereas an improvement is observed in the precipitation intensity. Intensified precipitation but with the same spatial variability of the raw output from the bias‐corrected GCM does not improve the heterogeneous runoff distributions, which in turn regulate unrealistically high peak downstream streamflow. GCM‐simulated precipitation with a large bias correction that is necessary to compensate for the poor performance in present climate simulation appears to distort streamflow patterns in the future projection, which leads to misleading projections of climate change impacts on hydrological extremes.     

Regional impacts of climate change on irrigation water demands

This paper presents an approach to model the expected impacts of climate change on irrigation water demand in a reservoir command area. A statistical downscaling model and an evapotranspiration model are used with a general circulation model (GCM) output to predict the anticipated change in the monthly irrigation water requirement of a crop. Specifically, we quantify the likely changes in irrigation water demands at a location in the command area, as a response to the projected changes in precipitation and evapotranspiration at that location. Statistical downscaling with a canonical correlation analysis is carried out to develop the future scenarios of meteorological variables (rainfall, relative humidity (RH), wind speed (U₂), radiation, maximum (Tmax) and minimum (Tmin) temperatures) starting with simulations provided by a GCM for a specified emission scenario. The medium resolution Model for Interdisciplinary Research on Climate GCM is used with the A1B scenario, to assess the likely changes in irrigation demands for paddy, sugarcane, permanent garden and semidry crops over the command area of Bhadra reservoir, India. Results from the downscaling model suggest that the monthly rainfall is likely to increase in the reservoir command area. RH, Tmax and Tmin are also projected to increase with small changes in U₂. Consequently, the reference evapotranspiration, modeled by the Penman–Monteith equation, is predicted to increase. The irrigation requirements are assessed on monthly scale at nine selected locations encompassing the Bhadra reservoir command area. The irrigation requirements are projected to increase, in most cases, suggesting that the effect of projected increase in rainfall on the irrigation demands is offset by the effect due to projected increase/change in other meteorological variables (viz., Tmax and Tmin, solar radiation, RH and U₂). The irrigation demand assessment study carried out at a river basin will be useful for future irrigation management systems. Copyright © 2012 John Wiley & Sons, Ltd.     

Hydrological modelling using data from monthly GCMs in a regional catchment

This study demonstrates the use of spatially downscaled, monthly general circulation model (GCM) rainfall and temperature data to drive the established HyMOD hydrological model to evaluate the prospective effects of climate change on the fluvial run‐off of the River Derwent basin in the UK. The evaluation results of this monthly hydrological model using readily available, monthly GCM data are consistent with studies on nearby catchments employing high‐temporal resolution data, indicating that useful hydro‐climatic planning studies may be possible using standard datasets and modest computational resources. HyMOD was calibrated against 5 km² gridded UK Climate Projections dataset data and then driven using monthly spatially interpolated (~5 km²) outputs from Hadley Centre Coupled Model, version 3 and the Canadian Centre for Climate Modelling and Analysis for Intergovernmental Panel on Climate Change Special Report on Emissions Scenarios (IPCC‐SRES) A2a and B2a covering the 2020s, 2050s and 2080s. Results for both GCMs project a decrease in annual run‐off in both GCM models and scenarios with higher values in the summer/autumn months, whereas an increase in the later winter months. Both Hadley Centre Coupled Model, version 3 and the Canadian Centre for Climate Modelling and Analysis show higher ranges of uncertainty during the winter season with higher values of run‐off associated with December in all three simulation periods and two scenarios. A seasonal comparison of run‐off simulations shows that both GCMs give similar results in summer and autumn, whereas disparities due to GCM uncertainties are more conspicuous in winter and spring. In this study, both the GCMs under A2a scenario have demonstrated the high possibility of time shift in monthly average peak run‐offs in the Derwent River by 2080s in comparison with the early 21st century. Copyright © 2013 John Wiley & Sons, Ltd.     

A comparison of three methods for downscaling daily precipitation in the Punjab region

Many downscaling techniques have been developed in the past few years for projection of station‐scale hydrological variables from large‐scale atmospheric variables simulated by general circulation models (GCMs) to assess the hydrological impacts of climate change. This article compares the performances of three downscaling methods, viz. conditional random field (CRF), K‐nearest neighbour (KNN) and support vector machine (SVM) methods in downscaling precipitation in the Punjab region of India, belonging to the monsoon regime. The CRF model is a recently developed method for downscaling hydrological variables in a probabilistic framework, while the SVM model is a popular machine learning tool useful in terms of its ability to generalize and capture nonlinear relationships between predictors and predictand. The KNN model is an analogue‐type method that queries days similar to a given feature vector from the training data and classifies future days by random sampling from a weighted set of K closest training examples. The models are applied for downscaling monsoon (June to September) daily precipitation at six locations in Punjab. Model performances with respect to reproduction of various statistics such as dry and wet spell length distributions, daily rainfall distribution, and intersite correlations are examined. It is found that the CRF and KNN models perform slightly better than the SVM model in reproducing most daily rainfall statistics. These models are then used to project future precipitation at the six locations. Output from the Canadian global climate model (CGCM3) GCM for three scenarios, viz. A1B, A2, and B1 is used for projection of future precipitation. The projections show a change in probability density functions of daily rainfall amount and changes in the wet and dry spell distributions of daily precipitation. Copyright © 2011 John Wiley & Sons, Ltd.     

Challenges in Studying and Modeling Heat and Moisture Exchange in Soil–Vegetation/Snow Cover–Surface Air Layer Systems

The current state of the art in the investigation and modeling of soil–vegetation/snow cover–surface air layer system (SVAS) is briefly reviewed. This system plays a decisive role in the formation of heat and moisture exchange between the land surface and atmosphere. The SWAP model, developed by the authors of this paper, is used to illustrate the potentialities of SVAS models in reproducing the dynamics of components of water and heat balances of the land surface under different natural conditions and at different spatial and temporal scales. The challenges and perspectives of further development of SVAS models are analyzed.     

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.     

Dynamical downscaling of regional climate: A review of methods and limitations

The traditional dynamical downscaling (TDD) method employs continuous integration of regional climate models (RCM) with the general circulation model (GCM) providing the initial and lateral boundary conditions. Dynamical downscaling simulations are constrained by physical principles and can generate a full set of climate information, providing one of the important approaches to projecting fine spatial-scale future climate information. However, the systematic biases of climate models often degrade the TDD simulations and hinder the application of dynamical downscaling in the climate-change related studies. New methods developed over past decades improve the performance of dynamical downscaling simulations. These methods can be divided into four groups: the TDD method, the pseudo global warming method, dynamical downscaling with GCM bias corrections, and dynamical downscaling with both GCM and RCM bias corrections. These dynamical downscaling methods are reviewed and compared in this paper. The merits and limitations of each dynamical downscaling method are also discussed. In addition, the challenges and potential directions in progressing dynamical downscaling methods are stated.     

An improved statistical analogue downscaling procedure for seasonal precipitation forecast

Seasonal forecasting can be highly valuable for water resources management. Hydrological models (either lumped conceptual rainfall-runoff models or physically based distributed models) can be used to simulate streamflows and update catchment conditions (e.g. soil moisture status) using rainfall records and other catchment data. However, in order to use any hydrological model for skillful seasonal forecasting, rainfall forecast at relevant spatial and/or temporal scales is required. Together with downscaling, general circulation models are probably the only tools for making such seasonal predictions. The Predictive Ocean Atmosphere Model for Australia (POAMA) is a state-of-the-art seasonal climate forecast system developed by the Australian Bureau of Meteorology. Based on the preliminary assessment on the performance of existing statistical downscaling methods used in Australia, this paper is devoted to develop an analogue downscaling method by modifying the Euclidian distance in the selection of similar weather pattern. Such a modification consists of multivariate Box–Cox transformation and then standardization to make the resulted POAMA and observed climate pattern more similar. For the predictors used in Timbal and Fernadez (CAWCR Technical Report No. 004, 2008), we also considered whether the POAMA precipitation provides useful information in the analogue method. Using the high quality station data in the Murray Darling Basin of Australia, we found that the modified analogue method has potential to improve the seasonal precipitation forecast using POAMA outputs. Finally, we found that in the analogue method, the precipitation from POAMA should not be used in the calculation of similarity. The findings would then help to improve the seasonal forecast of streamflows in Australia.     

Effects of soil cracking on the seismic response of soil–structure systems

A numerical study on the influence that cracks and discontinuities (closed cracks) can have on the seismic response of a hypothetical soil–structure system is presented and discussed. A 2-D finite-difference model of the soil was developed, considering a bilinear failure surface using a Mohr–Coulomb model. The cracks are simulated with interface elements. The soil stiffness is used to characterize the contact force that is generated when the crack closes. For the cases studied herein, it was considered that the crack does not propagate during the dynamic event. Both cases, open and closed cracks, are considered. The nonlinear behavior was accounted for approximately using equivalent linear properties calibrated against several 1-D wave propagation analyses of selected soil columns with variable depth to account for changes in depth to bed rock. Free field boundaries were used at the edges of the 2-D finite-difference model to allow for energy dissipation of the reflected waves. The effect of cracking on the seismic response was evaluated by comparing the results of site response analysis with and without crack, for several lengths and orientations. The changes in the response obtained for a single crack and a family of cracks were also evaluated. Finally, the impact that a crack may have on the structural response of nearby structures was investigated by solving the seismic-soil–structure interaction of two structures, one flexible and one rigid to bracket the response. From the results of this investigation, insight was gained regarding the effect that discontinuities may have both on the seismic response of soil deposits and on nearby soil–structure systems.     

Correcting bias in rainfall inputs to a semi-distributed hydrological model using downstream flow simulation errors

The focus of this study is on bias correcting semi-distributed rainfall inputs into a hydrological model applied in the Okavango River basin in southern Africa, where there are very few local observations and heavy reliance is placed on global rainfall datasets. While the hydrological model, before rainfall bias correction, is able to represent the broad characteristics of the sub-basin streamflow responses, as demonstrated by good agreement between observed and simulated flow duration curves, there are many years where the annual volumes are over- or underestimated. The long records of observed flow at downstream stations are successfully used to bias correct the rainfall inputs to the upstream sub-basins using an analysis of their individual contributions to downstream flow and their annual rainfall–runoff response ratios. The results show improved simulations for the relatively shorter observation periods at the upstream gauging stations.     

Scaling characteristics of spatial patterns of soil moisture from distributed modelling

Characterizing the spatial dynamics of soil moisture fields is a key issue in hydrology, offering an avenue to improve our understanding of complex land surface–atmosphere interactions. In this paper, the statistical structure of soil moisture patterns is examined using modelled soil moisture obtained from the North American Land Data Assimilation System (NLDAS) at 0.125° resolution. The study focuses on the vertically averaged soil moisture in the top 10 cm and 100 cm layers. The two variables display a weak dependence for lower values of surface soil moisture, with the strength of the relationship increasing with the water content of the top layer. In both cases, the variance of the soil moisture follows a power law decay as a function of the averaging area. The superficial layer shows a lower degree of spatial organization and higher temporal variability, which is reflected in rapid changes in time of the slope of the scaling functions of the soil moisture variance. Conversely, the soil moisture in the top 100 cm has lower variability in time and larger spatial correlation. The scaling of these patterns was found to be controlled by the changes in the soil water content. Results have implications for the downscaling of soil moisture to prevent model bias.