Downscaling daily precipitation over the Yellow River source region in China: a comparison of three statistical downscaling methods

Three statistical downscaling methods are compared with regard to their ability to downscale summer (June–September) daily precipitation at a network of 14 stations over the Yellow River source region from the NCEP/NCAR reanalysis data with the aim of constructing high-resolution regional precipitation scenarios for impact studies. The methods used are the Statistical Downscaling Model (SDSM), the Generalized LInear Model for daily CLIMate (GLIMCLIM), and the non-homogeneous Hidden Markov Model (NHMM). The methods are compared in terms of several statistics including spatial dependence, wet- and dry spell length distributions and inter-annual variability. In comparison with other two models, NHMM shows better performance in reproducing the spatial correlation structure, inter-annual variability and magnitude of the observed precipitation. However, it shows difficulty in reproducing observed wet- and dry spell length distributions at some stations. SDSM and GLIMCLIM showed better performance in reproducing the temporal dependence than NHMM. These models are also applied to derive future scenarios for six precipitation indices for the period 2046–2065 using the predictors from two global climate models (GCMs; CGCM3 and ECHAM5) under the IPCC SRES A2, A1B and B1scenarios. There is a strong consensus among two GCMs, three downscaling methods and three emission scenarios in the precipitation change signal. Under the future climate scenarios considered, all parts of the study region would experience increases in rainfall totals and extremes that are statistically significant at most stations. The magnitude of the projected changes is more intense for the SDSM than for other two models, which indicates that climate projection based on results from only one downscaling method should be interpreted with caution. The increase in the magnitude of rainfall totals and extremes is also accompanied by an increase in their inter-annual variability.     

Uncertainty analysis of statistical downscaling models using Hadley Centre Coupled Model

Regression-based statistical downscaling is a method broadly used to resolve the coarse spatial resolution of general circulation models. Nevertheless, the assessment of uncertainties linked with climatic variables is essential to climate impact studies. This study presents a procedure to characterize the uncertainty in regression-based statistical downscaling of daily precipitation and temperature over a highly vulnerable area (semiarid catchment) in the west of Iran, based on two downscaling models: a statistical downscaling model (SDSM) and an artificial neural network (ANN) model. Biases in mean, variance, and wet/dry spells are estimated for downscaled data using vigorous statistical tests for 30 years of observed and downscaled daily precipitation and temperature data taken from the National Center for Environmental Prediction reanalysis predictors for the years of 1961 to 1990. In the case of daily temperature, uncertainty is estimated by comparing monthly mean and variance of downscaled and observed daily data at a 95 % confidence level. In daily precipitation, downscaling uncertainties were evaluated from comparing monthly mean dry and wet spell lengths and their confidence intervals, cumulative frequency distributions of monthly mean of daily precipitation, and the distributions of monthly wet and dry days for observed and modeled daily precipitation. Results showed that uncertainty in downscaled precipitation is high, but simulation of daily temperature can reproduce extreme events accurately. Finally, this study shows that the SDSM is the most proficient model at reproducing various statistical characteristics of observed data at a 95 % confidence level, while the ANN model is the least capable in this respect. This study attempts to test uncertainties of regression-based statistical downscaling techniques in a semiarid area and therefore contributes to an improvement of the quality of predictions of climate change impact assessment in regions of this type.     

Comparison of regional climate model and statistical downscaling simulations of different winter precipitation change scenarios over Romania

Summary Regional climate model and statistical downscaling procedures are used to generate winter precipitation changes over Romania for the period 2071–2100 (compared to 1961–1990), under the IPCC A2 and B2 emission scenarios. For this purpose, the ICTP regional climate model RegCM is nested within the Hadley Centre global atmospheric model HadAM3H. The statistical downscaling method is based on the use of canonical correlation analysis (CCA) to construct climate change scenarios for winter precipitation over Romania from two predictors, sea level pressure and specific humidity (either used individually or together). A technique to select the most skillful model separately for each station is proposed to optimise the statistical downscaling signal. Climate fields from the A2 and B2 scenario simulations with the HadAM3H and RegCM models are used as input to the statistical downscaling model. First, the capability of the climate models to reproduce the observed link between winter precipitation over Romania and atmospheric circulation at the European scale is analysed, showing that the RegCM is more accurate than HadAM3H in the simulation of Romanian precipitation variability and its connection with large-scale circulations. Both models overestimate winter precipitation in the eastern regions of Romania due to an overestimation of the intensity and frequency of cyclonic systems over Europe. Climate changes derived directly from the RegCM and HadAM3H show an increase of precipitation during the 2071–2100 period compared to 1961–1990, especially over northwest and northeast Romania. Similar climate change patterns are obtained through the statistical downscaling method when the technique of optimum model selected separately for each station is used. This adds confidence to the simulated climate change signal over this region. The uncertainty of results is higher for the eastern and southeastern regions of Romania due to the lower HadAM3H and RegCM performance in simulating winter precipitation variability there as well as the reduced skill of the statistical downscaling model.     

Performance comparison of three predictor selection methods for statistical downscaling of daily precipitation

Predictor selection is a critical factor affecting the statistical downscaling of daily precipitation. This study provides a general comparison between uncertainties in downscaled results from three commonly used predictor selection methods (correlation analysis, partial correlation analysis, and stepwise regression analysis). Uncertainty is analyzed by comparing statistical indices, including the mean, variance, and the distribution of monthly mean daily precipitation, wet spell length, and the number of wet days. The downscaled results are produced by the artificial neural network (ANN) statistical downscaling model and 50 years (1961–2010) of observed daily precipitation together with reanalysis predictors. Although results show little difference between downscaling methods, stepwise regression analysis is generally the best method for selecting predictors for the ANN statistical downscaling model of daily precipitation, followed by partial correlation analysis and then correlation analysis.     

Comparison of statistically downscaled precipitation in terms of future climate indices and daily variability for southern Ontario and Quebec,Canada

Given the coarse resolution of global climate models, downscaling techniques are often needed to generate finer scale projections of variables affected by local-scale processes such as precipitation. However, classical statistical downscaling experiments for future climate rely on the time-invariance assumption as one cannot know the true change in the variable of interest, nor validate the models with data not yet observed. Our experimental setup involves using the Canadian regional climate model (CRCM) outputs as pseudo-observations to estimate model performance in the context of future climate projections by replacing historical and future observations with model simulations from the CRCM, nested within the domain of the Canadian global climate model (CGCM). In particular, we evaluated statistically downscaled daily precipitation time series in terms of the Peirce skill score, mean absolute errors, and climate indices. Specifically, we used a variety of linear and nonlinear methods such as artificial neural networks (ANN), decision trees and ensembles, multiple linear regression, and k-nearest neighbors to generate present and future daily precipitation occurrences and amounts. We obtained the predictors from the CGCM 3.1 20C3M (1971–2000) and A2 (2041–2070) simulations, and precipitation outputs from the CRCM 4.2 (forced with the CGCM 3.1 boundary conditions) as predictands. Overall, ANN models and tree ensembles outscored the linear models and simple nonlinear models in terms of precipitation occurrences, without performance deteriorating in future climate. In contrast, for the precipitation amounts and related climate indices, the performance of downscaling models deteriorated in future climate.     

Climate variability and trends in downscaled high-resolution simulations and projections over Metropolitan France

In order to fulfill the society demand for climate information at the spatial scale allowing impact studies, long-term high-resolution climate simulations are produced, over an area covering metropolitan France. One of the major goals of this article is to investigate whether such simulations appropriately simulate the spatial and temporal variability of the current climate, using two simulation chains. These start from the global IPSL-CM4 climate model, using two regional models (LMDz and MM5) at moderate resolution (15–20 km), followed with a statistical downscaling method in order to reach a target resolution of 8 km. The statistical downscaling technique includes a non-parametric method that corrects the distribution by using high-resolution analyses over France. First the uncorrected simulations are evaluated against a set of high-resolution analyses, with a focus on temperature and precipitation. Uncorrected downscaled temperatures suffer from a cold bias that is present in the global model as well. Precipitations biases have a season- and model-dependent behavior. Dynamical models overestimate rainfall but with different patterns and amplitude, but both have underestimations in the South-Eastern area (Cevennes mountains) in winter. A variance decomposition shows that uncorrected simulations fairly well capture observed variances from inter-annual to high-frequency intra-seasonal time scales. After correction, distributions match with analyses by construction, but it is shown that spatial coherence, persistence properties of warm, cold and dry episodes also match to a certain extent. Another aim of the article is to describe the changes for future climate obtained using these simulations under Scenario A1B. Results are presented on the changes between current and mid-term future (2021–2050) averages and variability over France. Interestingly, even though the same global climate model is used at the boundaries, regional climate change responses from the two models significantly differ.     

The influence of synoptic airflow on UK daily precipitation extremes. Part I: Observed spatio-temporal relationships

We study the influence of synoptic scale atmospheric circulation on extreme daily precipitation across the United Kingdom, using observed time series from 689 rain gauges. To this end we employ a statistical model, that uses airflow strength, direction and vorticity as predictors for the generalised extreme value distribution of monthly precipitation maxima. The inferred relationships are connected with the dominant westerly flow, the orography, and the moisture supply from surrounding seas. We aggregated the results for individual rain gauges to regional scales to investigate the temporal variability of extreme precipitation. Airflow explains a significant fraction of the variability on subannual to decadal time scales. A large fraction of the especially heavy winter precipitation during the 1980s and 1990s in north Scotland can be attributed to a prevailing positive phase of the North Atlantic Oscillation. Our statistical model can be used for statistical downscaling and to validate regional climate model output.     

Downscaling daily precipitation time series using a combined circulation- and regression-based approach

The aim of this paper is to introduce a new conditional statistical model for generating daily precipitation time series. The generated daily precipitation can thus be used for climate change impact studies, e.g., crop production, rainfall–runoff, and other water-related processes. It is a stochastic model that links local rainfall events to a continuous atmospheric predictor, moisture flux, in addition to classified atmospheric circulation patterns. The coupled moisture flux is proved to be capable of capturing continuous property of climate system and providing extra information to determine rainfall probability and rainfall amount. The application was made to simultaneously downscale daily precipitation at multiple sites within the Rhine River basin. The results show that the model can well reproduce statistical properties of daily precipitation time series. Especially for extreme rainfall events, the model is thought to better reflect rainfall variability compared to the pure CP-based downscaling approach.     

Precipitation and temperature space–time variability and extremes in the Mediterranean region: evaluation of dynamical and statistical downscaling methods

This study evaluates how statistical and dynamical downscaling models as well as combined approach perform in retrieving the space–time variability of near-surface temperature and rainfall, as well as their extremes, over the whole Mediterranean region. The dynamical downscaling model used in this study is the Weather Research and Forecasting (WRF) model with varying land-surface models and resolutions (20 and 50 km) and the statistical tool is the Cumulative Distribution Function-transform (CDF-t). To achieve a spatially resolved downscaling over the Mediterranean basin, the European Climate Assessment and Dataset (ECA&D) gridded dataset is used for calibration and evaluation of the downscaling models. In the frame of HyMeX and MED-CORDEX international programs, the downscaling is performed on ERA-I reanalysis over the 1989–2008 period. The results show that despite local calibration, CDF-t produces more accurate spatial variability of near-surface temperature and rainfall with respect to ECA&D than WRF which solves the three-dimensional equation of conservation. This first suggests that at 20–50 km resolutions, these three-dimensional processes only weakly contribute to the local value of temperature and precipitation with respect to local one-dimensional processes. Calibration of CDF-t at each individual grid point is thus sufficient to reproduce accurately the spatial pattern. A second explanation is the use of gridded data such as ECA&D which smoothes in part the horizontal variability after data interpolation and damps the added value of dynamical downscaling. This explains partly the absence of added-value of the 2-stage downscaling approach which combines statistical and dynamical downscaling models. The temporal variability of statistically downscaled temperature and rainfall is finally strongly driven by the temporal variability of its forcing (here ERA-Interim or WRF simulations). CDF-t is thus efficient as a bias correction tool but does not show any added-value regarding the time variability of the downscaled field. Finally, the quality of the reference observation dataset is a key issue. Comparison of CDF-t calibrated with ECA&D dataset and WRF simulations to local measurements from weather stations not assimilated in ECA&D, shows that the temporal variability of the downscaled data with respect to the local observations is closer to the local measurements than to ECA&D data. This highlights the strong added-value of dynamical downscaling which improves the temporal variability of the atmospheric dynamics with regard to the driving model. This article highlights the benefits and inconveniences emerging from the use of both downscaling techniques for climate research. Our goal is to contribute to the discussion on the use of downscaling tools to assess the impact of climate change on regional scales.     

自组织映射神经网络(SOM)降尺度方法对江淮流域逐日降水量的模拟评估

利用1961~2002年ERA-40逐日再分析资料和江淮流域56个台站逐日观测降水量资料,引入基于自组织映射神经网络(Self-Organizing Maps,简称SOM)的统计降尺度方法,对江淮流域夏季(6~8月)逐日降水量进行统计建模与验证,以考察SOM对中国东部季风降水和极端降水的统计降尺度模拟能力。结果表明,SOM通过建立主要天气型与局地降水的条件转换关系,能够再现与观测一致的日降水量概率分布特征,所有台站基于概率分布函数的Brier评分(Brier Score)均近似为0,显著性评分(Significance Score)全部在0.8以上;模拟的多年平均降水日数、中雨日数、夏季总降水量、日降水强度、极端降水阈值和极端降水贡献率区域平均的偏差都低于11%;并且能够在一定程度上模拟出江淮流域夏季降水的时间变率。进一步将SOM降尺度模型应用到BCCCSM1.1(m)模式当前气候情景下,评估其对耦合模式模拟结果的改善能力。发现降尺度显著改善了模式对极端降水模拟偏弱的缺陷,对不同降水指数的模拟较BCC-CSM1.1(m)模式显著提高,降尺度后所有台站6个降水指数的相对误差百分率基本在20%以内,偏差比降尺度前减小了40%~60%;降尺度后6个降水指数气候场的空间相关系数提高到0.9,相对标准差均接近1.0,并且均方根误差在0.5以下。表明SOM降尺度方法显著提高日降水概率分布,特别是概率分布曲线尾部特征的模拟能力,极大改善了模式对极端降水场的模拟能力,为提高未来预估能力提供了基础。     

Regional climate model projections for the State of Washington

Global climate models do not have sufficient spatial resolution to represent the atmospheric and land surface processes that determine the unique regional climate of the State of Washington. Regional climate models explicitly simulate the interactions between the large-scale weather patterns simulated by a global model and the local terrain. We have performed two 100-year regional climate simulations using the Weather Research and Forecasting (WRF) model developed at the National Center for Atmospheric Research (NCAR). One simulation is forced by the NCAR Community Climate System Model version 3 (CCSM3) and the second is forced by a simulation of the Max Plank Institute, Hamburg, global model (ECHAM5). The mesoscale simulations produce regional changes in snow cover, cloudiness, and circulation patterns associated with interactions between the large-scale climate change and the regional topography and land-water contrasts. These changes substantially alter the temperature and precipitation trends over the region relative to the global model result or statistical downscaling. To illustrate this effect, we analyze the changes from the current climate (1970–1999) to the mid twenty-first century (2030–2059). Changes in seasonal-mean temperature, precipitation, and snowpack are presented. Several climatological indices of extreme daily weather are also presented: precipitation intensity, fraction of precipitation occurring in extreme daily events, heat wave frequency, growing season length, and frequency of warm nights. Despite somewhat different changes in seasonal precipitation and temperature from the two regional simulations, consistent results for changes in snowpack and extreme precipitation are found in both simulations.     

基于统计降尺度方法的夏季降水精细化预报

利用1999—2009年安徽省淮河以南地区60个县市站夏季逐日降水资料和安庆市探空站逐日资料,研究了中低层不同风向配置下局地降水与大尺度降水场之间的关系,以3种不同预报对象及相应的预报因子分别采用神经网络和线性回归方法设计6种预报模型对观测资料进行逼近和优化,从而实现空间降尺度.分析对比6种预报模型46站逐日降水量的拟合和预报效果,结果表明:采取相同的预报对象及预报因子的BP神经网络模型在拟合和预报效果上均好于线性回归模型,可见夏季降水场之间以非线性相关为主;神经网络模型预报结果同常用的Cressman插值预报相比,能很好地反映出降水的基本分布及局地特征;预报对象为单站降水序列的神经网络模型在以平原、河流为主要地形的区域预报效果较好,预报对象为REOF主成分的神经网络模型则在山地和丘陵地形区域预报效果较好.     

Application of physical scaling towards downscaling climate model precipitation data

Physical scaling (SP) method downscales climate model data to local or regional scales taking into consideration physical characteristics of the area under analysis. In this study, multiple SP method based models are tested for their effectiveness towards downscaling North American regional reanalysis (NARR) daily precipitation data. Model performance is compared with two state-of-the-art downscaling methods: statistical downscaling model (SDSM) and generalized linear modeling (GLM). The downscaled precipitation is evaluated with reference to recorded precipitation at 57 gauging stations located within the study region. The spatial and temporal robustness of the downscaling methods is evaluated using seven precipitation based indices. Results indicate that SP method-based models perform best in downscaling precipitation followed by GLM, followed by the SDSM model. Best performing models are thereafter used to downscale future precipitations made by three global circulation models (GCMs) following two emission scenarios: representative concentration pathway (RCP) 2.6 and RCP 8.5 over the twenty-first century. The downscaled future precipitation projections indicate an increase in mean and maximum precipitation intensity as well as a decrease in the total number of dry days. Further an increase in the frequency of short (1-day), moderately long (2–4 day), and long (more than 5-day) precipitation events is projected.     

A comparison of statistical downscaling and climate change factor methods: impacts on low flows in the River Thames, United Kingdom

Strategic-scale assessments of climate change impacts are often undertaken using the change factor (CF) methodology whereby future changes in climate projected by General Circulation Models (GCMs) are applied to a baseline climatology. Alternatively, statistical downscaling (SD) methods apply climate variables from GCMs to statistical transfer functions to estimate point-scale meteorological series. This paper explores the relative merits of the CF and SD methods using a case study of low flows in the River Thames under baseline (1961–1990) and climate change conditions (centred on the 2020s, 2050s and 2080s). Archived model outputs for the UK Climate Impacts Programme (UKCIP02) scenarios are used to generate daily precipitation and potential evaporation (PE) for two climate change scenarios via the CF and SD methods. Both signal substantial reductions in summer precipitation accompanied by increased PE throughout the year, leading to reduced flows in the Thames in late summer and autumn. However, changes in flow associated with the SD scenarios are generally more conservative and complex than that arising from CFs. These departures are explained in terms of the different treatment of multidecadal natural variability, temporal structuring of daily climate variables and large-scale forcing of local precipitation and PE by the two downscaling methods.     

中国夏季降水的组合统计降尺度模型预测研究

现阶段的动力气候模式尚不能满足东亚区域气候预测的实际需求,这就需要动力和统计相结合的方法,将动力模式中具有较高预测技巧的大尺度环流信息应用到降水等气象要素的统计预测模型当中,以改善后者预测效果。本文中所介绍的组合统计降尺度模型,可将动力气候模式预测的大尺度环流变量和前期观测的外强迫信号作为预测因子来预测中国夏季降水异常。交叉检验结果显示,组合统计降尺度预测模型的距平相关系数较原始模式结果有较大提高。在实时夏季降水预测中,2013～2018年平均的预测技巧相对较高,趋势异常综合检验（PS）评分平均为71.5分,特别是2015～2018年平均的PS评分预测技巧达到72.7分,总体上高于业务模式原始预测和业务发布预测的技巧。该组合统计降尺度模型预测性能稳定,为我国季节预测业务提供了一种有效参考。     

Drought episodes over Greece as simulated by dynamical and statistical downscaling approaches

Drought over the Greek region is characterized by a strong seasonal cycle and large spatial variability. Dry spells longer than 10 consecutive days mainly characterize the duration and the intensity of Greek drought. Moreover, an increasing trend of the frequency of drought episodes has been observed, especially during the last 20 years of the 20th century. Moreover, the most recent regional circulation models (RCMs) present discrepancies compared to observed precipitation, while they are able to reproduce the main patterns of atmospheric circulation. In this study, both a statistical and a dynamical downscaling approach are used to quantify drought episodes over Greece by simulating the Standardized Precipitation Index (SPI) for different time steps (3, 6, and 12 months). A statistical downscaling technique based on artificial neural network is employed for the estimation of SPI over Greece, while this drought index is also estimated using the RCM precipitation for the time period of 1961–1990. Overall, it was found that the drought characteristics (intensity, duration, and spatial extent) were well reproduced by the regional climate models for long term drought indices (SPI12) while ANN simulations are better for the short-term drought indices (SPI3).

     

Statistical downscaling of daily mean temperature, pan evaporation and precipitation for climate change scenarios in Haihe River, China

A statistical downscaling method (SDSM) was evaluated by simultaneously downscaling air temperature, evaporation, and precipitation in Haihe River basin, China. The data used for evaluation were large-scale atmospheric data encompassing daily NCEP/NCAR reanalysis data and the daily mean climate model results for scenarios A2 and B2 of the HadCM3 model. Selected as climate variables for downscaling were measured daily mean air temperature, pan evaporation, and precipitation data (1961–2000) from 11 weather stations in the Haihe River basin. The results obtained from SDSM showed that: (1) the pattern of change in and numerical values of the climate variables can be reasonably simulated, with the coefficients of determination between observed and downscaled mean temperature, pan evaporation, and precipitation being 99%, 93%, and 73%, respectively; (2) systematic errors existed in simulating extreme events, but the results were acceptable for practical applications; and (3) the mean air temperature would increase by about 0.7°C during 2011~2040; the total annual precipitation would decrease by about 7% in A2 scenario but increase by about 4% in B2 scenario; and there were no apparent changes in pan evaporation. It was concluded that in the next 30 years, climate would be warmer and drier, extreme events could be more intense, and autumn might be the most distinct season among all the changes.     

Statistical downscaling of daily climate variables for climate change impact assessment over New South Wales, Australia

This paper outlines a new statistical downscaling method based on a stochastic weather generator. The monthly climate projections from global climate models (GCMs) are first downscaled to specific sites using an inverse distance-weighted interpolation method. A bias correction procedure is then applied to the monthly GCM values of each site. Daily climate projections for the site are generated by using a stochastic weather generator, WGEN. For downscaling WGEN parameters, historical climate data from 1889 to 2008 are sorted, in an ascending order, into 6 climate groups. The WGEN parameters are downscaled based on the linear and non-linear relationships derived from the 6 groups of historical climates and future GCM projections. The overall averaged confidence intervals for these significant linear relationships between parameters and climate variables are 0.08 and 0.11 (the range of these parameters are up to a value of 1.0) at the observed mean and maximum values of climate variables, revealing a high confidence in extrapolating parameters for downscaling future climate. An evaluation procedure is set up to ensure that the downscaled daily sequences are consistent with monthly GCM output in terms of monthly means or totals. The performance of this model is evaluated through the comparison between the distributions of measured and downscaled climate data. Kruskall-Wallis rank (K-W) and Siegel-Tukey rank sum dispersion (S-T) tests are used. The results show that the method can reproduce the climate statistics at annual, monthly and daily time scales for both training and validation periods. The method is applied to 1062 sites across New South Wales (NSW) for 9 GCMs and three IPCC SRES emission scenarios, B1, A1B and A2, for the period of 1900–2099. Projected climate changes by 7 GCMs are also analyzed for the A2 emission scenario based on the downscaling results.     

两种统计降尺度方法在秦岭山地的适用性

基于ASD（automated statistical downscaling）统计降尺度模型提供的多元线性回归和岭回归两种统计降尺度方法,采用RCP4.5（representative concentration pathways 4.5）和RCP8.5情景下全球气候模式MPI-ESM-LR输出的预报因子数据、NCEP/NCAR再分析数据和秦岭山地周边10个气象站观测数据,评估两种统计降尺度方法在秦岭山地的适用性及预估秦岭山地未来3个时期（2006-2040年、2041-2070年和2071-2100年）的平均气温和降水。结果表明:率定期和验证期内,两种统计降尺度方法均可以较好地模拟研究区域的平均气温和降水的变化特征,且多元线性回归的模拟效果优于岭回归。在未来气候情景下,两种统计降尺度方法预估的研究区域平均气温均呈明显上升趋势,气温增幅随辐射强迫增加而增大。降水方面,21世纪未来3个时期降水均呈不明显减少趋势,但季节分配发生变化。综合考虑两种统计降尺度方法在秦岭山地对平均气温和降水的模拟效果和情景预估结果,认为多元线性回归降尺度方法更适用于秦岭山地气候变化的降尺度预估研究。     

Statistical downscaling of regional daily precipitation over southeast Australia based on self-organizing maps

This paper presents a novel statistical downscaling method based on a non-linear classification technique known as self-organizing maps (SOMs) and has therefore been named SOM-SD. The relationship between large-scale atmospheric circulation and local-scale surface variable was constructed in a relatively simple and transparent manner. For a specific atmospheric state, an ensemble of possible values was generated for the predictand following the Monte Carlo method. Such a stochastic simulation is essential to explore the uncertainties of climate change in the future through a series of random re-sampling experiments. The novel downscaling method was evaluated by downscaling daily precipitation over Southeast Australia. The large-scale predictors were extracted from the daily NCAR/NCEP reanalysis data, while the predictand was high-resolution gridded daily observed precipitation (1958?C2008) from the Australian Bureau of Meteorology. The results showed that the method works reasonably well across a variety of climatic zones in the study area. Overall, there was no particular zone that stands out as a climatic entity where the downscaling skill in reproducing all statistical indices was consistently lower or higher across seasons than the other zones. The method displayed a high skill in reproducing not only the climatologic statistical properties of the observed precipitation, but also the characteristics of the extreme precipitation events. Furthermore, the model was able to reproduce, to a certain extent, the inter-annual variability of precipitation characteristics.