Spatial Patterns of Precipitation Anomalies for 30-yr Warm Periods in China During the Past 2000 Years

The spatial patterns of precipitation anomalies during five 30-yr warm periods of 691-720, 1231-1260, 1741-1770, 1921-1950, and 1981-2000 were investigated using a dryness/wetness grading dataset covering 48 stations from Chinese historical documents and 22 precipitation proxy series from natural archives. It was found that the North China Plain (approximately 35 -40 N, east of 105 E) was dry in four warm periods within the centennial warm epochs of 600-750, the Medieval Warm Period (about 900-1300) and after 1900. A wet condition prevailed over most of China during 1741-1770, a 30-yr warm peak that occurred during the Little Ice Age (about 1650-1850). The spatial pattern of the precipitation anomaly in 1981-2000 over East China (25 -40 N, east of 105 E) is roughly consistent with that in 1231-1260, but a difference in the precipitation anomaly appeared over the Tibetan Plateau. The spatial patterns of the precipitation anomalies over China varied between all five 30-yr warm periods, which implies that the matching pattern between temperature and precipitation change is multiform, and the precipitation anomaly could be positive or negative when a decadal warm climate occurs in different climate epochs. This result may provide a primary reference for the mechanism detection and climate simulation of the precipitation anomaly of the future warm climate.     

1980-2020年山西省极端降水的研究

基于1980-2020年山西省109个气象观测站点的逐日降水资料,选取10个极端降水指数,采用气候倾向率、相关分析、因子分析、R/S预测方法等方法,对山西省极端降水进行了时空分布的研究,以期为山西省的气候变化、生态环境保护、防灾减灾、气象服务工作提供参考依据,结果表明：(1)从时间尺度来看,1980-2020年期间,山西省极端降水的强度和极值都有明显增加,连续干旱日数和连续湿日日数呈下降趋势,其余均表现出不同程度的增加,其中年总降水量增加幅度最明显;从空间尺度来看,年总降水量、降水强度、降水频率、极值均为从西北向东南逐渐增多,空间差异较明显;从各站点的空间分布来看,北部和中部地区的极端事件增加最显著,北部地区的干旱日数仍以增加趋势为主,连续湿日日数气候倾向率的空间差异较大,中部地区站点显著增加,南北部以减少趋势为主;(2)基于相关分析方法表明各极端降水指数(除干旱日数外)与年总降水量都有很好的相关关系,强降水量和极强降水量对年总降水量的贡献值呈现出增加趋势;采用因子分析方法提取了3个公共因子,方差贡献率累计达到了87％,可以看出极端降水强度和降水量指数在对极端降水方面影响较大;利用R/S分析法可以得到年总降水量、中雨日数、大雨日数、最大5日降水量这几个指数未来呈现弱减少趋势,而干旱日数仍为减少趋势,连续湿日日数为持续弱增加趋势。总体看来,山西省极端降水近年来呈现出增加趋势,在空间分布有明显差异。     

Simulation of extreme precipitation indices in the Yangtze River basin by using statistical downscaling method (SDSM)

In this study, the applicability of the statistical downscaling model (SDSM) in modeling five extreme precipitation indices including R10 (no. of days with precipitation ≥10?mm?day^?1), SDI (simple daily intensity), CDD (maximum number of consecutive dry days), R1d (maximum 1-day precipitation total) and R5d (maximum 5-day precipitation total) in the Yangtze River basin, China was investigated. The investigation mainly includes the calibration and validation of SDSM model on downscaling daily precipitation, the validation of modeling extreme precipitation indices using independent period of the NCEP reanalysis data, and the projection of future regional scenarios of extreme precipitation indices. The results showed that: (1) there existed good relationship between the observed and simulated extreme precipitation indices during validation period of 1991–2000, the amount and the change pattern of extreme precipitation indices could be reasonably simulated by SDSM. (2) Under both scenarios A2 and B2, during the projection period of 2010–2099, the changes of annual mean extreme precipitation indices in the Yangtze River basin would be not obvious in 2020s; while slightly increase in the 2050s; and significant increase in the 2080s as compared to the mean values of the base period. The summer might be the more distinct season with more projected increase of each extreme precipitation indices than in other seasons. And (3) there would be distinctive spatial distribution differences for the change of annual mean extreme precipitation indices in the river basin, but the most of Yangtze River basin would be dominated by the increasing trend.     

Spatiotemporal variability of drought on a short–medium time scale in the Calabria Region (Southern Italy)

The study analyses spatial and temporal patterns of drought in an area with a wide range of precipitation characteristics (the Calabria region in southern Italy) during the period 1921?C2007. The short-time (2, 3 and 6?months) Standardised Precipitation Index (SPI) was estimated to analyse drought especially from the agricultural point of view. Principal component analysis (PCA) was applied to the SPI to assess the spatial variability of drought. During the period of observation, moderate to severe drought occurred at a frequency of almost 13?% in wet seasons (autumn and winter). Almost half of the region was affected by drought in the years 1981?C1990 when the area experienced its most severe drought. Although the spatial patterns of drought estimated by PCA were logical and consistent with precipitation distribution, very complex patterns were observed for all the time scales looked at. The high fragmentation of the maps obtained makes them of limited value, and caution is recommended in classifying the region into homogeneous areas.     

The best fitting distribution of annual maximum rainfall in Peninsular Malaysia based on methods of L-moment and LQ-moment

This paper discusses patterns of annual and monthly precipitation variability at seven weather stations in east central Europe (1851–2007). Precipitation patterns were compared to three simple regional indices of atmospheric circulation, i.e., western circulation, southern circulation and the cyclonicity (C) index and a relationship between precipitation and the North Atlantic Oscillation index was identified. Correlations of the monthly records and multiple regression, using a principal components’ analysis, helped determine the statistical significance of the dependence of precipitation on the circulation indices. The Mann–Kendall test revealed no trend to change in any of the precipitation series, but a certain spatial regularity could be discerned in the phase of the annual periodic component. A common feature of the variability in central European annual precipitation is the dry period identified in the 1980s and the first half of the 1990s. In the northern part of the region, above-average precipitation was noted from the 1960s through to the mid-1970s as a result of the frequent prevalence of depressions. South of the divide, the wettest period was recorded at the turn of 1930s/1940s. After a number of very wet years in the last decade of the twentieth century and the beginning of the twenty-first century, precipitation began to fall at all of the region’s weather stations. The C index is the strongest circulation-linked factor influencing precipitation in central Europe and it accounts for more than 40% of the variance in spatially averaged wintertime precipitation.     

中国长序列资料测站月降水、月均温相关联系的稳定性分析

选取我国长序列资料测站点中月降水与月均温资料较全的6个站点：西安、广州、汉口、沈阳、重庆、上海,以10a为一次成组实验,时域分辨率定为季,计算出各季月降水与月均温的相关系数,然后用相关分析和贝努里概型分析月均温与月降水之间的相关联系,初步得出我国长序列资料测站点春、夏、秋、冬月降水与月均温的相关稳定性联系。结果表明：①我国月降水与月均温在春、夏、秋、冬季均存在相关关系,以夏、秋季节最显著且为负相关,故我国常见干热或湿凉型夏、秋季;②我国华南地区春季月降水与月均温负相关关系显著,故常呈现干热或湿冷型春季;西北地区冬季月降水与月均温负相关关系较为明显,常呈现干热或湿冷型冬季。     

Spatio-temporal patterns of precipitation in Serbia

The monthly precipitation data from 29 synoptic stations for the period 1946–2012 were analyzed using a number of different multivariate statistical analysis methods to investigate the spatial variability and temporal patterns of precipitation across Serbia. R-mode principal component analysis was used to study the spatial variability of the precipitation. Three distinct sub-regions were identified by applying the agglomerative hierarchical cluster analysis to the two component scores: C1 includes the north and the northeast part of Serbia, while C2 includes the western part of Central Serbia and southwestern part of Serbia and C3 includes central, east, south and southeast part of Serbia. The analysis of the identified sub-regions indicated that the monthly and seasonal precipitation in sub-region C2 had the values above average, while C1 and C3 had the precipitation values under average. The analysis of the linear trend of the mean annual precipitation showed an increasing trend for the stations located in Serbia and three sub-regions. From the result of this analysis, one can plan land use, water resources and agricultural production in the region.     

模糊C均值聚类法在新疆年降水气候分区中的应用

利用新疆98个气象站1960-2011年的年降水量资料,采用模糊C均值聚类法,对新疆年降水量进行分区研究;同时利用线性趋势、累积距平、M-K检验、t检验相结合的方法,对新疆年降水量在不同区域上的变化趋势以及突变时间等进行了对比诊断分析,得到了新疆地区年降水量分布的空间特征。结果表明：（1）新疆年降水的分布大致可分为7个区域;（2）根据模糊C均值方法所分区域能够体现出由于地形差异导致的降水分布不均匀的特点,这一结果与其他研究干旱区降水分布差异产生的原因相同。分区结果合理,说明该聚类方法适用于区域气候区划;（3）新疆大部分区域的降水量在20世纪80年代中后期以前偏少,低于多年平均值,自80年代中后期以后才开始偏多;不同区域的降水量依次从70年代初、80年代中后期、90年代初开始增加;整个南疆盆地、阿勒泰地区、准噶尔盆地降水的突变时间较为接近;（4）新疆降水量整体异常表现为降水一致多（或少）、北多南少、西多东少。     

Spatiotemporal patterns of drought at various time scales in Shandong Province of Eastern China

The temporal variations and spatial patterns of drought in Shandong Province of Eastern China were investigated by calculating the standardized precipitation evapotranspiration index (SPEI) at 1-, 3-, 6-, 12-, and 24-month time scales. Monthly precipitation and air temperature time series during the period 1960–2012 were collected at 23 meteorological stations uniformly distributed over the region. The non-parametric Mann-Kendall test was used to explore the temporal trends of precipitation, air temperature, and the SPEI drought index. S-mode principal component analysis (PCA) was applied to identify the spatial patterns of drought. The results showed that an insignificant decreasing trend in annual total precipitation was detected at most stations, a significant increase of annual average air temperature occurred at all the 23 stations, and a significant decreasing trend in the SPEI was mainly detected at the coastal stations for all the time scales. The frequency of occurrence of extreme and severe drought at different time scales generally increased with decades; higher frequency and larger affected area of extreme and severe droughts occurred as the time scale increased, especially for the northwest of Shandong Province and Jiaodong peninsular. The spatial pattern of drought for SPEI-1 contains three regions: eastern Jiaodong Peninsular and northwestern and southern Shandong. As the time scale increased to 3, 6, and 12 months, the order of the three regions was transformed into another as northwestern Shandong, eastern Jiaodong Peninsular, and southern Shandong. For SPEI-24, the location identified by REOF1 was slightly shifted from northwestern Shandong to western Shandong, and REOF2 and REOF3 identified another two weak patterns in the south edge and north edge of Jiaodong Peninsular, respectively. The potential causes of drought and the impact of drought on agriculture in the study area have also been discussed. The temporal variations and spatial patterns of drought obtained in this study provide valuable information for water resources planning and drought disaster prevention and mitigation in Eastern China.     

Spatial and temporal variability of precipitation indices during 1961–2010 in Hunan Province,central south China

Based on daily precipitation records at 75 meteorological stations in Hunan Province, central south China, the spatial and temporal variability of precipitation indices is analyzed during 1961–2010. For precipitation extremes, most of precipitation indices suggest that both the amount and the intensity of extreme precipitation are increasing, especially the mean precipitation amount on a wet day, showing a significant positive trend. Meanwhile, both of the monthly rainfall heterogeneity and the contribution of the days with the greatest rainfall show an upward trend. When it comes to rainfall erosivity, most of this province is characterized by high values of annual rainfall erosivity. Although the directions of trends in annual rainfall erosivity at most stations are upward, only 6 of the 75 stations have significant trends. Furthermore, the spatial and temporal variation of dryness/wetness has been assessed by the standardized precipitation index (SPI). The principal component analysis (PCA) was applied to the SPI series computed on 24-month time scales. The results demonstrated a noticeable spatial variability with three subregions characterized by different trends: a remarkable wet tendency prevails in the central and southern areas, while the northern areas are dominated by a remarkable dry tendency.     

High resolution regional climate model simulations for Germany: Part II—projected climate changes

The projected climate change signals of a five-member high resolution ensemble, based on two global climate models (GCMs: ECHAM5 and CCCma3) and two regional climate models (RCMs: CLM and WRF) are analysed in this paper (Part II of a two part paper). In Part I the performance of the models for the control period are presented. The RCMs use a two nest procedure over Europe and Germany with a final spatial resolution of 7 km to downscale the GCM simulations for the present (1971–2000) and future A1B scenario (2021–2050) time periods. The ensemble was extended by earlier simulations with the RCM REMO (driven by ECHAM5, two realisations) at a slightly coarser resolution. The climate change signals are evaluated and tested for significance for mean values and the seasonal cycles of temperature and precipitation, as well as for the intensity distribution of precipitation and the numbers of dry days and dry periods. All GCMs project a significant warming over Europe on seasonal and annual scales and the projected warming of the GCMs is retained in both nests of the RCMs, however, with added small variations. The mean warming over Germany of all ensemble members for the fine nest is in the range of 0.8 and 1.3 K with an average of 1.1 K. For mean annual precipitation the climate change signal varies in the range of ?2 to 9 % over Germany within the ensemble. Changes in the number of wet days are projected in the range of ±4 % on the annual scale for the future time period. For the probability distribution of precipitation intensity, a decrease of lower intensities and an increase of moderate and higher intensities is projected by most ensemble members. For the mean values, the results indicate that the projected temperature change signal is caused mainly by the GCM and its initial condition (realisation), with little impact from the RCM. For precipitation, in addition, the RCM affects the climate change signal significantly.     

Characterizing the distribution of observed precipitation and runoff over the continental United States

This paper describes the development of a comprehensive geographic database of historical precipitation and runoff measurements for the conterminous U.S. The database is used in a spatial analysis to characterize large scale precipitation and runoff patterns and to assess the utility and limitations of using historical hydro-meteorological data for providing spatially distributed precipitation estimates at regional and continental scales. Long-term annual average precipitation (P) and runoff (Q) surfaces (geographically referenced, digital representations of a continuous spatial distribution) generated from interpolation of point measurements are used in a distributed water balance calculation to check the reliability of precipitation estimates. The resulting input-output values (P- Q) illustrate the deficiency (sparse distribution and low elevation bias) of historical precipitation measurements in the mountainous western U.S. where snowmelt is an important component of the annual runoff. The incorporation of high elevation snow measurements into the precipitation record significantly improves the water balance estimates in some areas and enhances the utility of historical data for providing spatially distributed precipitation estimates in topographically diverse regions. Regions where the use of historical precipitation data may be most limited for precipitation estimation are identified and alternatives to the use of interpolated historical data for precipitation estimation across large heterogenous regions are suggested. The research establishes a database for continental scale studies and provides direction for the successful development of spatially distributed regional scale water balance models.     

CMIP6不同分辨率全球气候模式对中国降水模拟能力评估

基于参与CMIP6高分辨率模式比较计划(HighResMIP)9个模式组的18个全球气候模式模拟数据,通过与CN05.1观测资料的对比,评估了不同分辨率气候模式对中国区域1961—2014年降水特征的模拟能力.结果表明:低、高分辨率模式均能模拟出中国区域多年平均降水的总体空间分布特征,以及降水冬弱夏强的季节变化特征,但...     

Spatial modeling of annual minimum and maximum temperatures in Iceland

Bayesian hierarchical modeling is applied to the analysis of annual minimum and maximum temperatures. In both cases the generalized extreme value distribution is selected as the marginal distribution at each site due to its flexibility and theoretical basis. This distribution has three unknown parameters, a location parameter, a scale parameter and a shape parameter. The location and scale parameters are assumed to vary across sites while the shape parameter is assumed to be the same for all sites. The location and the scale parameters are modeled as two independent Gaussian spatial processes which are governed by the Matérn correlation function. These spatial processes play a central role in spatial predictions at unobserved sites. The location parameter is also modeled with a common time effect which includes a linear trend and a linear site effect with altitude, distance to open sea, latitude and longitude as covariates. This model allows for spatial predictions for any set of unobserved sites, the scale of the grid can be as fine as possible as long as the covariates are observed at each of the unobserved sites. Data on the annual minimum and maximum temperatures in Iceland from 1961 to 2009?at 72 sites are analyzed and used to predict the 2nd percentile of the minimum temperature for the 12-month period from July 2011 to June 2012 and the 98th percentile of maximum temperature for the year 2011 for a large set of unobserved sites across Iceland at which the covariates are available. The spatial predictions reveal that the 2nd percentile of minimum temperature for 2011 to 2012 reaches ?35°C in the central part of Iceland and ranges from ?12 to ?22°C around the coast of Iceland with the lowest coastal temperatures in the Northern part and the Northeastern part. The 98th percentile of maximum temperature in 2011 ranges from 10 to 30°C where the warmest areas are found in the Eastern part, the North-Eastern part and the Southern part of the country. The estimated increase in minimum and maximum temperatures over the years 1961 to 2009 is 0.71 and 0.47°C per decade, respectively, while the average annual temperature increased 0.24°C per decade over the same period.     

Simulation of regional climate change under the IPCC A2 scenario in southeast China

A variable-grid atmospheric general circulation model, LMDZ, with a local zoom over southeast China is used to investigate regional climate changes in terms of both means and extremes. Two time slices of 30?years are chosen to represent, respectively, the end of the 20th century and the middle of the 21st century. The lower-boundary conditions (sea-surface temperature and sea-ice extension) are taken from the outputs of three global coupled climate models: Institut Pierre-Simon Laplace (IPSL), Centre National de Recherches Météorologiques (CNRM) and Geophysical Fluid Dynamics Laboratory (GFDL). Results from a two-way nesting system between LMDZ-global and LMDZ-regional are also presented. The evaluation of simulated temperature and precipitation for the current climate shows that LMDZ reproduces generally well the spatial distribution of mean climate and extreme climate events in southeast China, but the model has systematic cold biases in temperature and tends to overestimate the extreme precipitation. The two-way nesting model can reduce the ??cold bias?? to some extent compared to the one-way nesting model. Results with greenhouse gas forcing from the SRES-A2 emission scenario show that there is a significant increase for mean, daily-maximum and minimum temperature in the entire region, associated with a decrease in the number of frost days and an increase in the heat wave duration. The annual frost days are projected to significantly decrease by 12?C19?days while the heat wave duration to increase by about 7?days. A warming environment gives rise to changes in extreme precipitation events. Except two simulations (LMDZ/GFDL and LMDZ/IPSL2) that project a decrease in maximum 5-day precipitation (R5d) for winter, other precipitation extremes are projected to increase over most of southeast China in all seasons, and among the three global scenarios. The domain-averaged values for annual simple daily intensity index (SDII), R5d and fraction of total rainfall from extreme events (R95t) are projected to increase by 6?C7, 10?C13 and 11?C14%, respectively, relative to their present-day values. However, it is clear that more research will be needed to assess the uncertainties on the projection in future of climate extremes at local scale.     

Abrupt changes in rainfall in the Eastern area of La Pampa Province, Argentina

The eastern area of La Pampa Province, Argentina, lies in a transition zone between the humid temperate climate stretching east and the steppe climate stretching west. The area is thus very sensitive to abrupt changes in rainfall. In order to determine the long-term occurrence of such phenomena, long-term annual precipitation series (1921?C2004) from 17 stations in the study area were analyzed using the Buishand and Pettitt tests. Results showed a sharp increase in annual rainfall at the southern stations in the 1960s and at the northern and central stations in the 1970s. Increased rainfall can be considered one of the reasons for the subsequent expansion in land planted to crops in the region. While a rapid increase in rainfall can be seen as positive, some researchers believe that if an abrupt decrease in rainfall occurred in future and continued for long, the carrying capacity of the environment could be exceeded, leading to decreased production and environmental degradation.     

Climate change in Central America and Mexico: regional climate model validation and climate change projections

Central America has high biodiversity, it harbors high-value ecosystems and it??s important to provide regional climate change information to assist in adaptation and mitigation work in the region. Here we study climate change projections for Central America and Mexico using a regional climate model. The model evaluation shows its success in simulating spatial and temporal variability of temperature and precipitation and also in capturing regional climate features such as the bimodal annual cycle of precipitation and the Caribbean low-level jet. A variety of climate regimes within the model domain are also better identified in the regional model simulation due to improved resolution of topographic features. Although, the model suffers from large precipitation biases, it shows improvements over the coarse-resolution driving model in simulating precipitation amounts. The model shows a dry bias in the wet season and a wet bias in the dry season suggesting that it??s unable to capture the full range of precipitation variability. Projected warming under the A2 scenario is higher in the wet season than that in the dry season with the Yucatan Peninsula experiencing highest warming. A large reduction in precipitation in the wet season is projected for the region, whereas parts of Central America that receive a considerable amount of moisture in the form of orographic precipitation show significant decreases in precipitation in the dry season. Projected climatic changes can have detrimental impacts on biodiversity as they are spatially similar, but far greater in magnitude, than those observed during the El Ni?o events in recent decades that adversely affected species in the region.     

Assessing future climate changes and extreme indicators in east and south Asia using the RegCM4 regional climate model

This paper assesses future climate changes over East and South Asia using a regional climate model (RegCM4) with a 50?km spatial resolution. To evaluate the model performance, RegCM4 is driven with ??perfect boundary forcing?? from the reanalysis data during 1970?C1999 to simulate the present day climate. The model performs well in reproducing not only the mean climate and seasonality but also most of the chosen indicators of climate extremes. Future climate changes are evaluated based on two experiments driven with boundary forcing from the European-Hamburg general climate model (ECHAM5), one for the present (1970?C1999) and one for the SRES A1B future scenario (2070?C2099). The model predicts an annual temperature increase of about 3°?C5° (smaller over the ocean and larger over the land), and an increase of annual precipitation over most of China north of 30°N and a decrease or little change in the rest of China, India and Indochina. For temperature-related extreme indicators in the future, the model predicts a generally longer growing season, more hot days in summer, and less frost days in winter. For precipitation-related extremes, the number of days with more than 10?mm of rainfall is predicted to increase north of 30°N and decrease in the south, and the maximum five-day rainfall amount and daily intensity will increase across the whole model domain. In addition, the maximum number of consecutive dry days is predicted to increase over most of the model domain, south of 40°N. Most of the Yangtze River Basin in China stands out as ??hotspots?? of extreme precipitation changes, with the strongest increases of daily rain intensity, maximum five-day rain amount, and the number of consecutive dry days, suggesting increased risks of both floods and droughts.     

Searching for Added Value in Simulating Climate Extremes with a High-Resolution Regional Climate Model over Western Canada. II: Basin-Scale Results

We evaluate the capacity of a regional climate model to simulate the statistics of extreme events, and also examine the effect of differing horizontal resolution, at the scale of individual hydrological basins in the topographically complex province of British Columbia, Canada. Two climate simulations of western Canada (WCan) were conducted with the Canadian Regional Climate Model (version 4) at 15 (CRCM15) and 45?km (CRCM45) horizontal resolution driven at the lateral boundaries by global reanalysis over the period 1973–1995. The simulations were evaluated with ANUSPLIN, a daily observational gridded surface temperature and precipitation product and with meteorological data recorded at 28 stations within the upper Peace, Nechako, and upper Columbia River basins. In this work, we focus largely on a comparison of the skill of each model configuration in simulating the 90th percentile of daily precipitation (PR90). The companion paper describes the results for a wider range of temperature and precipitation extremes over the entire WCan domain.

Over all three watersheds, both simulations exhibit cold biases compared with observations, with the bias exacerbated at higher resolution. Although both simulations generally display wet biases in median precipitation, CRCM15 features a reduced bias in PR90 in all three basins in summer and throughout the year in the upper Columbia River basin. However, the higher resolution model is inferior to CRCM45 with respect to rarer heavy precipitation events and also displays high spatial variability and lower spatial correlations with ANUSPLIN compared with the coarser resolution model. A reduction in the range of PR90 biases over the upper Columbia basin is noted when the 15?km results are averaged to the 45?km grid. This improvement is partly attributable to the averaging of errors between different elevation data used in the gridded observations and CRCM, but the sensitivity of CRCM15 to resolved topography is also clear from spatial maps of seasonal extremes. At the station scale, modest but systematic reductions in the bias of PR90 relative to ANUSPLIN are again found when the CRCM15 results are averaged to the 45?km grid. Furthermore, the annual cycle of inter-station spatial variance in the upper Columbia River basin is well reproduced by CRCM15 but not by ANUSPLIN or CRCM45. The former result highlights the beneficial effect of spatial averaging of small-scale climate variability, whereas the latter is evidently a demonstration of the added value at high resolution vis-à-vis the improved simulation of precipitation at the resolution limit of the model.     

Bulgarian mountains air temperatures and precipitation—statistical downscaling of global climate models and some projections

Measured air temperature and precipitation data from three high mountainous Bulgarian stations were used along with data from 18 global climate models (GCMs). Air temperature and precipitation outputs of preindustrial control experiment were compared with actually observed values. GCM with the best overall performance is BCCR BCM 2.0 for air temperatures (period 1941?C2009) and CGCM 3.1/T47 for precipitation (period 1947?C2009). Statistical methods were used in this research??nonparametric Spearman correlation, Mann?CWhitney test, multiple linear regression, etc. Projections were made for the following future decades: 2015?C2024, 2045?C2054 and 2075?C2084. The best months, described by multiple linear regression (MLR) model of air temperatures, are November, January, March, and May. The worst described are summer months. There is not any pattern in the relationship between constructed MLR models and measured precipitation. Models that perform the best in different months at the three investigated stations are MIUB ECHO-G, GISS AOM, CGCM 3.1/T63, and CNRM CM3 for air temperatures and GFDL CM 2.1, GISS AOM, and MIUB ECHO-G for precipitation. The fit between statistical models' outputs and values observed at stations is different, better in cold part of the year. There will be mixed future changes of air temperatures at all the three high mountainous stations. An increase of temperatures is expected in April, November, and December. A decrease will happen in February, July, and October. Mean annual temperatures are expected to rise by 0.1?°C (Botev) to 0.2?°C (Musala and Cherni vrah) in the decade 2075?C2084, but mean annual temperatures at the end of the period with measurements (2009) has already exceeded by far projected values. Trends in precipitation are mixed both in spatial and in temporal directions. Observed decrease of precipitation, especially in the warm half of the year, is not described well in MLR models. The same is valid for annual amounts, which are projected to be higher than those measured in the end of instrumental period (2009). This is opposite to observed trends in recent decades, especially at stations Cherni vrah and Botev, where a significant decrease of precipitation amounts has happened.