Uncertainty analysis of statistically downscaled temperature and precipitation regimes in Northern Canada

Summary Uncertainty analysis is used to make a quantitative evaluation of the reliability of statistically downscaled climate data representing local climate conditions in the northern coastlines of Canada. In this region, most global climate models (GCMs) have inherent weaknesses to adequately simulate the climate regime due to difficulty in resolving strong land/sea discontinuities or heterogeneous land cover. The performance of the multiple regression-based statistical downscaling model in reproducing the observed daily minimum/maximum temperature, and precipitation for a reference period (1961–1990) is evaluated using climate predictors derived from NCEP reanalysis data and those simulated by two coupled GCMs (the Canadian CGCM2 and the British HadCM3). The Wilcoxon Signed Rank test and bootstrap confidence-interval estimation techniques are used to perform uncertainty analysis on the downscaled meteorological variables. The results show that the NCEP-driven downscaling results mostly reproduced the mean and variability of the observed climate very well. Temperatures are satisfactorily downscaled from HadCM3 predictors while some of the temperatures downscaled from CGCM2 predictors are statistically significantly different from the observed. The uncertainty in precipitation downscaled with CGCM2 predictors is comparable to the ones downscaled from HadCM3. In general, all downscaling results reveal that the regression-based statistical downscaling method driven by accurate GCM predictors is able to reproduce the climate regime over these highly heterogeneous coastline areas of northern Canada. The study also shows the applicability of uncertainty analysis techniques in evaluating the reliability of the downscaled data for climate scenarios development. Authors’ addresses: Dr. Yonas B. Dibike, NSERC Research Fellow, OURANOS Consortium, 550 Sherbrooke Street West, 19^th Floor, Montreal (QC) H3A 1B9, Canada; Philippe Gachon, Adaptation and Impact Research Division (AIRD), Atmospheric Science and Technology Directorate, Environment Canada at Ouranos, Montreal (QC), Canada; André St-Hilaire and Taha B. M. J. Ouarda, Institut National de la Recherche Scientifique Centre Eau, Terre & Environnement (INRS-ETE), University of Québec, 490 Rue de La Couronne, Québec (QC) G1K 9A9, Canada; Van T.-V. Nguyen, Department of Civil Engineering and Applied Mechanics, McGill University, 817 Sherbrooke Street West, Montreal (QC) H3A 2K6, Canada.     

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.     

Precipitation,temperature and wind in Norway: dynamical downscaling of ERA40

A novel downscaling approach of the ERA40 (ECMWF 40-years reanalysis) data set has been taken and results for comparison with observations in Norway are shown. The method applies a nudging technique in a stretched global model, focused in the Norwegian Sea (67°N, 5°W). The effective resolution is three times the one of the ERA40, equivalent to about 30 km grid spacing in the area of focus. Longer waves (<T42) in the downscaled solution are nudged towards the ERA40 solution, and thus the large-scale circulation is similar in the two data sets. The shorter waves are free to evolve, and produce high intensities of winds and precipitation. The comparison to observations incorporate numerous station data points of (1) precipitation (#357), (2) temperature (#98) and (3) wind (#10), and for the period 1961–1990, the downscaled data set shows large improvements over ERA40. The daily precipitation shows considerable reduction in bias (from 50 to 11%), and twofold reduction at the 99.9 percentile (from −59 to −29%). The daily temperature showed a bias reduction of about a degree in most areas, and relative large RMSE reduction (from 7.5 to 5.0°C except winter). The wind comparison showed a slight improvement in bias, and significant improvements in RMSE.     

基于数值模式的月尺度近地层气象要素预报技术研究

利用WRF模式对美国NCEP发布的CFS气候预测业务产品在中国区域内进行动力降尺度预报,可得到预报时效为45天的逐6小时、30 km分辨率基础气象要素预测产品。再利用全国气象站观测资料和3个风电场70 m高度风速、温度观测资料对2015年冬季预测结果进行检验评估和分析,最后通过线性方法对地面要素预测结果和70 m高度风速、温度预测结果进行统计订正。结果表明:（1）2 m温度和相对湿度的全国预报平均绝对误差分别为4.71 ℃和18.81%,在华东、华中和华南地区误差较小;（2）10 m风速预报平均绝对误差为2.42 m/s,在东北、华北和西北地区误差较小;（3）线性订正后,2 m气温、相对湿度和10 m风速的预报绝对误差分别减小1.05 ℃、5.29%和1.47 m/s,并且订正后误差随时间变化更平稳;（4）订正后70 m高度风速和温度的预报绝对误差均减小,风速平均误差减小最大可达1.29 m/s（B塔）,气温平均绝对误差减小最大可达3 ℃（C塔）。研究结果表明,基于CFS产品和WRF模式的、与月尺度风电预报关系密切的气象要素预报性能较好,未来可将该方法尝试于风电场的月尺度功率预测产品研发。      

Future climate change scenarios over Korea using a multi-nested downscaling system: A pilot study

This study examines a scenario of future summer climate change for the Korean peninsula using a multi-nested regional climate system. The global-scale scenario from the ECHAM5, which has a 200 km grid, was downscaled to a 50 km grid over Asia using the National Centers for Environmental Prediction (NCEP) Regional Spectral Model (RSM). This allowed us to obtain large-scale forcing information for a one-way, double-nested Weather and Research Forecasting (WRF) model that consists of a 12 km grid over Korea and a 3 km grid near Seoul. As a pilot study prior to the multi-year simulation work the years 1995 and 2055 were selected for the present and future summers. This RSM-WRF multi-nested downscaling system was evaluated by examining a downscaled climatology in 1995 with the largescale forcing from the NCEP/Department of Energy (DOE) reanalysis. The changes in monsoonal flows over East Asia and the associated precipitation change scenario over Korea are highlighted. It is found that the RSM-WRF system is capable of reproducing large-scale features associated with the East-Asian summer monsoon (EASM) and its associated hydro-climate when it is nested by the NCEP/DOE reanalysis. The ECHAM5-based downscaled climate for the present (1995) summer is found to suffer from a weakening of the low-level jet and sub-tropical high when compared the reanalysis-based climate. Predicted changes in summer monsoon circulations between 1995 and 2055 include a strengthened subtropical high and an intensified mid-level trough. The resulting projected summer precipitation is doubled over much of South Korea, accompanied by a pronounced surface warming with a maximum of about 2 K. It is suggested that downscaling strategy of this study, with its cloud-resolving scale, makes it suitable for providing high-resolution meteorological data with which to derive hydrology or air pollution models.     

农田灌溉过程与土壤参数更新对中亚费尔干纳盆地WRF/Noah模拟精度的提升

通过中亚费尔干纳盆地2007～2011年气候的模拟试验,揭示了新增农田灌溉过程与更新土壤参数对WRF（Weather Research and Forecasting）/Noah模式模拟精度的提升作用。通过对比标准版本与嵌入灌溉过程参数化方案后的WRF/Noah模式的模拟结果,研究发现农业灌溉提升了土壤含水量,导致地表蒸发增强,潜热增加,感热减少,致使近地层大气降温、增湿,这一效应降低了WRF/Noah模拟的暖、干偏差,模拟2 m气温和大气比湿均方根误差分别由6.52°C降低至5.81°C,由1.66 g/kg降低至1.13 g/kg。进而针对WRF默认配置的费尔干纳盆地内土壤数据精度欠佳的问题,再利用国际土壤参比与信息中心（ISRIC）数据（主要是粉砂粘壤土和粉砂壤土）替换了WRF默认的数据（主要是粘土和壤土）,降低了土壤凋萎系数,使得有效土壤水增多,缩小了灌溉需水量的模拟误差,并使得蒸散发进一步增强,潜热增多,感热减少,导致近地层降温、增湿,进一步降低了WRF/Noah模拟的暖、干偏差,模拟温度、湿度的均方根误差分别由5.81°C降低至5.46°C,由1.13 g/kg降低至1.08 g/kg。上述结果表明:充分农业灌溉对陆面过程产生影响,以及采用高精度的土壤数据能够显著提高WRF/Noah模式在中亚费尔干纳盆地的模拟精度。     

Joint variable spatial downscaling

Joint Variable Spatial Downscaling (JVSD), a new statistical technique for downscaling gridded climatic variables, is developed to generate high resolution gridded datasets for regional watershed modeling and assessments. The proposed approach differs from previous statistical downscaling methods in that multiple climatic variables are downscaled simultaneously and consistently to produce realistic climate projections. In the bias correction step, JVSD uses a differencing process to create stationary joint cumulative frequency statistics of the variables being downscaled. The functional relationship between these statistics and those of the historical observation period is subsequently used to remove GCM bias. The original variables are recovered through summation of bias corrected differenced sequences. In the spatial disaggregation step, JVSD uses a historical analogue approach, with historical analogues identified simultaneously for all atmospheric fields and over all areas of the basin under study. Analysis and comparisons are performed for 20th Century Climate in Coupled Models (20C3M), broadly available for most GCMs. The results show that the proposed downscaling method is able to reproduce the sub-grid climatic features as well as their temporal/spatial variability in the historical periods. Comparisons are also performed for precipitation and temperature with other statistical and dynamic downscaling methods over the southeastern US and show that JVSD performs favorably. The downscaled sequences are used to assess the implications of GCM scenarios for the Apalachicola-Chattahoochee-Flint river basin as part of a comprehensive climate change impact assessment.     

Conditional stochastic simulation model for spatial downscaling for assessing the effects of climate change on hydro-meteorological variables

The current study examines the recently proposed “bias correction and stochastic analogues” (BCSA) statistical spatial downscaling technique and attempts to improve it by conditioning coarse resolution data when generating replicates. While the BCSA method reproduces the statistical features of the observed fine data, this existing model does not replicate the observed coarse spatial pattern, and subsequently, the cross-correlation between the observed coarse data and downscaled fine data with the model cannot be preserved. To address the dissimilarity between the BCSA downscaled data and observed fine data, a new statistical spatial downscaling method, “conditional stochastic simulation with bias correction” (BCCS), which employs the conditional multivariate distribution and principal component analysis, is proposed. Gridded observed climate data of mean daily precipitation (mm/day) covering a month at 1/8° for a fine resolution and at 1° for a coarse resolution over Florida for the current and future periods were used to verify and cross-validate the proposed technique. The observed coarse and fine data cover the 50-year period from 1950 to1999, and the future RCP4.5 and RCP8.5 climate scenarios cover the 100-year period from 2000 to 2099. The verification and cross-validation results show that the proposed BCCS downscaling method serves as an effective alternative means of downscaling monthly precipitation levels to assess climate change effects on hydrological variables. The RCP4.5 and RCP8.5 GCM scenarios are successfully downscaled.     

Underestimation of the Warming Trend over the Tibetan Plateau during 1998–2013 by Global Land Data Assimilation Systems and Atmospheric Reanalyses

Accurate surface air temperature (T2m) data are key to investigating eco-hydrological responses to global warming. Because of sparse in-situ observations, T2m datasets from atmospheric reanalysis or multi-source observation-based land data assimilation system (LDAS) are widely used in research over alpine regions such as the Tibetan Plateau (TP). It has been found that the warming rate of T2m over the TP accelerates during the global warming slowdown period of 1998–2013, which raises the question of whether the reanalysis or LDAS datasets can capture the warming feature. By evaluating two global LDASs, five global atmospheric reanalysis datasets, and a high-resolution dynamical downscaling simulation driven by one of the global reanalysis, we demonstrate that the LDASs and reanalysis datasets underestimate the warming trend over the TP by 27%–86% during 1998–2013. This is mainly caused by the underestimations of the increasing trends of surface downward radiation and nighttime total cloud amount over the southern and northern TP, respectively. Although GLDAS2.0, ERA5, and MERRA2 reduce biases of T2m simulation from their previous versions by 12%-94%, they do not show significant improvements in capturing the warming trend. The WRF dynamical downscaling dataset driven by ERA-Interim shows a great improvement, as it corrects the cooling trend in ERA-Interim to an observation-like warming trend over the southern TP. Our results indicate that more efforts are needed to reasonably simulate the warming features over the TP during the global warming slowdown period, and the WRF dynamical downscaling dataset provides more accurate T2m estimations than its driven global reanalysis dataset ERA-Interim for producing LDAS products over the TP.     

WRF模式对中国夏季降水的动力降尺度模拟研究

采用NCEP的FNL再分析资料驱动WRF模式,对中国10 a(2000—2009年)夏季降水进行双重动力降尺度(双重嵌套)模拟,将子、母区域模拟结果和观测进行对比,以检验双重动力降尺度对中国夏季降水模拟的"增值"能力。结果表明:单重动力降尺度(单重嵌套)方法能较好模拟出中国10 a夏季平均降水的空间分布,对季风雨带"北跳"特征模拟较好,但模拟降水具有系统性正偏差。在母区域的强迫下,双重动力降尺度模拟的降水分布与单重动力降尺度相比,没有发生根本性变化。但由于子区域的分辨率要高于母区域,双重动力降尺度比单重动力降尺度能提供更多有价值的降水细节。双重动力降尺度的这种"增值"能力存在地域依赖性,在华南地区和江淮地区,双重动力降尺度模拟出的降水分布、量值和逐日演变都要好于单重动力降尺度。但在华北地区,双重动力降尺度没有表现出明显的"增值"。     

长时间序列多源土壤湿度产品在中国地区的比较分析

研究利用实测数据对遥感产品(ECV CCI)、全球陆面同化系统产品(GLDAS/CLM、NOAH2.7、NOAH3.3、VIC、MOSAIC)、再分析资料(ERA-Interim、NCEP/DOE)等不同来源的8套长时间序列土壤湿度产品进行时空比较。结果表明:上述产品均可以模拟出中国不同区域土壤湿度的空间分布;从各产品平均态与实测数据的偏差、均方根误差、相关系数统计结果来看,NOAH模式和ERA-Interim表现最好,MOSAIC和NCEP/DOE与实测相差较大;从区域尺度上来看,所有产品在东北区域表现最好,在华南、西南南区和西北西区较差;ERA-Interim有效的模拟了土壤湿度的年际变化。     

Field significance of performance measures in the context of regional climate model evaluation. Part 1: temperature

A new approach for rigorous spatial analysis of the downscaling performance of regional climate model (RCM) simulations is introduced. It is based on a multiple comparison of the local tests at the grid cells and is also known as “field” or “global” significance. New performance measures for estimating the added value of downscaled data relative to the large-scale forcing fields are developed. The methodology is exemplarily applied to a standard EURO-CORDEX hindcast simulation with the Weather Research and Forecasting (WRF) model coupled with the land surface model NOAH at 0.11 ^° grid resolution. Monthly temperature climatology for the 1990–2009 period is analysed for Germany for winter and summer in comparison with high-resolution gridded observations from the German Weather Service. The field significance test controls the proportion of falsely rejected local tests in a meaningful way and is robust to spatial dependence. Hence, the spatial patterns of the statistically significant local tests are also meaningful. We interpret them from a process-oriented perspective. In winter and in most regions in summer, the downscaled distributions are statistically indistinguishable from the observed ones. A systematic cold summer bias occurs in deep river valleys due to overestimated elevations, in coastal areas due probably to enhanced sea breeze circulation, and over large lakes due to the interpolation of water temperatures. Urban areas in concave topography forms have a warm summer bias due to the strong heat islands, not reflected in the observations. WRF-NOAH generates appropriate fine-scale features in the monthly temperature field over regions of complex topography, but over spatially homogeneous areas even small biases can lead to significant deteriorations relative to the driving reanalysis. As the added value of global climate model (GCM)-driven simulations cannot be smaller than this perfect-boundary estimate, this work demonstrates in a rigorous manner the clear additional value of dynamical downscaling over global climate simulations. The evaluation methodology has a broad spectrum of applicability as it is distribution-free, robust to spatial dependence, and accounts for time series structure.     

Statistical downscaling of hourly and daily climate scenarios for various meteorological variables in South-central Canada

Summary A regression-based methodology was used to downscale hourly and daily station-scale meteorological variables from outputs of large-scale general circulation models (GCMs). Meteorological variables include air temperature, dew point, and west–east and south–north wind velocities at the surface and three upper atmospheric levels (925, 850, and 500 hPa), as well as mean sea-level air pressure and total cloud cover. Different regression methods were used to construct downscaling transfer functions for different weather variables. Multiple stepwise regression analysis was used for all weather variables, except total cloud cover. Cumulative logit regression was employed for analysis of cloud cover, since cloud cover is an ordered categorical data format. For both regression procedures, to avoid multicollinearity between explanatory variables, principal components analysis was used to convert inter-correlated weather variables into uncorrelated principal components that were used as predictors. The results demonstrated that the downscaling method was able to capture the relationship between the premises and the response; for example, most hourly downscaling transfer functions could explain over 95% of the total variance for several variables (e.g. surface air temperature, dew point, and air pressure). Downscaling transfer functions were validated using a cross-validation scheme, and it was concluded that the functions for all weather variables used in the study are reliable. Performance of the downscaling method was also evaluated by comparing data distributions and extreme weather characteristics of downscaled GCM historical runs and observations during the period 1961–2000. The results showed that data distributions of downscaled GCM historical runs for all weather variables are significantly similar to those of observations. In addition, extreme characteristics of the downscaled meteorological variables (e.g. temperature, dew point, air pressure, and total cloud cover) were examined. Authors’ addresses: Chad Shouquan Cheng, Guilong Li, Qian Li, Atmospheric Science and Applications Unit, Meteorological Service of Canada Branch-Ontario, Environment Canada, 4905 Dufferin Street, Toronto, Ontario, Canada M3H 5T4; Heather Auld, Adaptation and Impacts Research Division, MSC Branch, Environment Canada, Toronto, Canada.     

Validating a regional climate model’s downscaling ability for East Asian summer monsoonal interannual variability

Performance of a regional climate model (RCM), WRF, for downscaling East Asian summer season climate is investigated based on 11-summer integrations associated with different climate conditions with reanalysis data as the lateral boundary conditions. It is found that while the RCM is essentially unable to improve large-scale circulation patterns in the upper troposphere for most years, it is able to simulate better lower-level meridional moisture transport in the East Asian summer monsoon. For precipitation downscaling, the RCM produces more realistic magnitude of the interannual variation in most areas of East Asia than that in the reanalysis. Furthermore, the RCM significantly improves the spatial pattern of summer rainfall over dry inland areas and mountainous areas, such as Mongolia and the Tibetan Plateau. Meanwhile, it reduces the wet bias over southeast China. Over Mongolia, however, the performance of precipitation downscaling strongly depends on the year: the WRF is skillful for normal and wet years, but not for dry years, which suggests that land surface processes play an important role in downscaling ability. Over the dry area of North China, the WRF shows the worst performance. Additional sensitivity experiments testing land effects in downscaling suggest the initial soil moisture condition and representation of land surface processes with different schemes are sources of uncertainty for precipitation downscaling. Correction of initial soil moisture using the climatology dataset from GSWP-2 is a useful approach to robustly reducing wet bias in inland areas as well as to improve spatial distribution of precipitation. Despite the improvement on RCM downscaling, regional analyses reveal that accurate simulation of precipitation over East China, where the precipitation pattern is strongly influenced by the activity of the Meiyu/Baiu rainfall band, is difficult. Since the location of the rainfall band is closely associated with both lower-level meridional moisture transport and upper-level circulation structures, it is necessary to have realistic upper-air circulation patterns in the RCM as well as lower-level moisture transport in order to improve the circulation-associated convective rainfall band in East Asia.     

我国地面降水的分级回归统计降尺度预报研究

利用TIGGE资料中欧洲中期天气预报中心(ECMWF,the European Centre for Medium-Range Weather Forecasts)、日本气象厅(JMA,the Japan Meteorological Agency)、美国国家环境预报中心(NCEP,the National Centers for Environmental Prediction)以及英国气象局(UKMO,the UK Met Office)4个中心1~7 d预报的日降水量集合预报资料,并以中国降水融合产品作为"观测值",对我国地面降水量预报进行统计降尺度处理。采用空间滑动窗口增加中雨和大雨雨量样本,建立分级雨量的回归方程,并与未分级雨量的统计降尺度预报进行对比。结果表明,对于不同模式、不同预报时效以及不同降水量级,统计降尺度的预报技巧改进程度不尽相同。统计降尺度的预报技巧依赖于模式本身的预报效果。相比雨量未分级回归,雨量分级回归的统计降尺度预报与观测值的距平相关系数更高,均方根误差更小,不同量级降水的ETS评分明显提高。对雨量分级回归统计降尺度预报结果进行二次订正,可大大减少小雨的空报。     

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.     

A Quick Report on a Dynamical Downscaling Simulation over China Using the Nested Model

This paper describes a dynamical downscaling simulation over China using the nested model system, which consists of the modified Weather Research and Forecasting Model (WRF) nested with the NCAR Community Atmosphere Model (CAM). Results show that dynamical downscaling is of great value in improving the model simulation of regional climatic characteristics. WRF simulates regional detailed temperature features better than CAM. With the spatial correlation coefficient between the observation and the simulation increasing from 0.54 for CAM to 0.79 for WRF, the improvement in precipitation simulation is more perceptible with WRF. Furthermore, the WRF simulation corrects the spatial bias of the precipitation in the CAM simulation.     

High-resolution precipitation and temperature downscaling for glacier models

The spatial resolution gap between global or regional climate models and the requirements for local impact studies motivates the need for climate downscaling. For impact studies that involve glacier modelling, the sparsity or complete absence of climate monitoring activities within the regions of interest presents a substantial additional challenge. Downscaling methods for this application must be independent of climate observations and cannot rely on tuning to station data. We present new, computationally-efficient methods for downscaling precipitation and temperature to the high spatial resolutions required to force mountain glacier models. Our precipitation downscaling is based on an existing linear theory for orographic precipitation, which we modify for large study regions by including moist air tracking. Temperature is downscaled using an interpolation scheme that reconstructs the vertical temperature structure to estimate surface temperatures from upper air data. Both methods are able to produce output on km to sub-km spatial resolution, yet do not require tuning to station measurements. By comparing our downscaled precipitation (1 km resolution) and temperature (200 m resolution) fields to station measurements in southern British Columbia, we evaluate their performance regionally and through the annual cycle. Precipitation is improved by as much as 30% (median relative error) over the input reanalysis data and temperature is reconstructed with a mean bias of 0.5°C at locations with high vertical relief. Both methods perform best in mountainous terrain, where glaciers tend to be concentrated.     

Coupling statistical and dynamical methods for spatial downscaling of precipitation

The resolution of General Circulation Models (GCMs) is too coarse for climate change impact studies at the catchment or site-specific scales. To overcome this problem, both dynamical and statistical downscaling methods have been developed. Each downscaling method has its advantages and drawbacks, which have been described in great detail in the literature. This paper evaluates the improvement in statistical downscaling (SD) predictive power when using predictors from a Regional Climate Model (RCM) over a GCM for downscaling site-specific precipitation. Our approach uses mixed downscaling, combining both dynamic and statistical methods. Precipitation, a critical element of hydrology studies that is also much more difficult to downscale than temperature, is the only variable evaluated in this study. The SD method selected here uses a stepwise linear regression approach for precipitation quantity and occurrence (similar to the well-known Statistical Downscaling Model (SDSM) and called SDSM-like herein). In addition, a discriminant analysis (DA) was tested to generate precipitation occurrence, and a weather typing approach was used to derive statistical relationships based on weather types, and not only on a seasonal basis as is usually done. The existing data record was separated into a calibration and validation periods. To compare the relative efficiency of the SD approaches, relationships were derived at the same sites using the same predictors at a 300km scale (the National Center for Environmental Prediction (NCEP) reanalysis) and at a 45km scale with data from the limited-area Canadian Regional Climate Model (CRCM) driven by NCEP data at its boundaries. Predictably, using CRCM variables as predictors rather than NCEP data resulted in a much-improved explained variance for precipitation, although it was always less than 50?% overall. For precipitation occurrence, the SDSM-like model slightly overestimated the frequencies of wet and dry periods, while these were well-replicated by the DA-based model. Both the SDSM-like and DA-based models reproduced the percentage of wet days, but the wet and dry statuses for each day were poorly downscaled by both approaches. Overall, precipitation occurrence downscaled by the DA-based model was much better than that predicted by the SDSM-like model. Despite the added complexity, the weather typing approach was not better at downscaling precipitation than approaches without classification. Overall, despite significant improvements in precipitation occurrence prediction by the DA scheme, and even going to finer scales predictors, the SD approach tested here still explained less than 50?% of the total precipitation variance. While going to even smaller scale predictors (10–15?km) might improve results even more, such smaller scales would basically transform the direct outputs of climate models into impact models, thus negating the need for statistical downscaling approaches.