Characteristics and Scenarios Projection of Climate Change on the Tibetan Plateau

The Fourth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC AR4) presents twenty-two global climate models (GCMs). In this paper, we evaluate the ability of 22 GCMs to reproduce temperature and precipitation over the Tibetan Plateau by comparing with ground observations for 1961~1900. The results suggest that all the GCMs underestimate surface air temperature and most models overestimate precipitation in most regions on the Tibetan Plateau. Only a few models (each 5 models for precipitation and temperature) appear roughly consistent with the observations in annual temperature and precipitation variations. Comparatively, GFCM21 and CGMR are able to better reproduce the observed annual temperature and precipitation variability over the Tibetan Plateau. Although the scenarios predicted by the GCMs vary greatly, all the models predict consistently increasing trends in temperature and precipitation in most regions in the Tibetan Plateau in the next 90 years. The results suggest that the temperature and precipitation will both increase in all three periods under different scenarios, with scenario A1 increasing the most and scenario A1B increasing the least.


Introduction
In the context of global warming, climate specialists all over the world show great attention to the projection for future local climate change. Currently, global climate models are considered as an important tool for understanding attributions of past climate change and predicting the future [1][2][3][4][5]. Therefore, it is very important to assess the ability of those models to reproduce the observed climatological features, which will directly affect the "reproduction" of current decadal climate changes, and verify to some extent the credibility for future climate change projections.
IPCC (the Intergovernmental Panel on Climate Change) Data Distribution Center (IPCC-DDC for short) released the Fourth Assessment Report with climate scenarios simulations and projections with 24 climate models. Some climate models predicted climate changes under 9 emission scenarios, among which 3 scenarios are identified as the most important ones for future climate change: high scenario SRES A2 (Special Report on Emission Scenarios A2); middle scenario SRES A1B; and low scenario SRES B1 [1,2]. It has become a hot research issue to assess global and regional simulation ability and predict climate change tendency for different emission scenarios with a single model or multiple models [3,4]. Phillips and Gleckr [6] evaluated the ability of the 20 IPCC-AR4 (the Fourth Assessment Report of the Intergovernmental Panel on Climate Change) models to reproduce global land annual mean precipitation; the result shows that new global climate models have better simulations of global land precipitation than the previous version. Based on the atmospheric circulation features of European climate change, AP van Ulden and van Oldenborgh [7] assessed the simulation ability of global climate models over Europe by calculating spatial correlation between observed and simulated values; the result shows that there are 8 models which have well-projected European atmospheric circulation changes. Johnson and Sharma [8] assessed the credibility of global climate models in time and space using the "Variable  Convergence Score (VCS)" method. Sun et al. [9] evaluated the ability of new global climate models to reproduce frequency, intensity, and other indicators of precipitation, which indicated that precipitation was still too frequent in the latest models. Chinese scientists have been working on future climate change assessment in China. Zhou and Yu [10] assessed the ability of the 19 IPCC AR4 models to reproduce precipitation climatology in China; the result shows that new GCMs simulated much better than the previous with regard to China's climate; most models could represent the average of surface temperature, but the simulation of warming ratio in recent years still needs to be improved. With the 22 IPCC AR4 models, Xu et al. [11] evaluated ability of various models to reproduce China's climate changes in the 20th century; and then pointed out that the temperature simulation of each model was better than precipitation, while that multimodel ensembles mean simulated much better than a single model. Jiang et al. [12,13] analyzed the ability of the IPCC-AR4 models to reproduce precipitation, temperature, and extreme precipitation index and project future climate changes under different emission scenarios in China. In addition, a lot of research has been done by scholars on climate change assessments and projections of GCMs [14][15][16][17].
Messerli and Ives [18] pointed out that mountainous and highland areas were vulnerable and sensitive to climate changes. The climate environments ecosystems in these regions respond most rapidly and significantly to global climate changes [19]. Climate warming on the plateau in past decades has been suggested by meteorological observation and ice core records [20]. The climatic warming trend seems more evident on the plateau than globally during the period from 1950 to 1993 [21]. Therefore, trend analysis of climate changes in high-altitude regions is very important for research on climate changes tendency. Being at an average elevation about 4000 m, the Tibetan Plateau is known as the roof of the world, which covers over 200 million km 2 with steep terrain. Besides it has far-reaching effects of global climate changes [22]. Although the Tibetan Plateau has a significant impact on atmospheric circulation and hydrological cycle, researches on the Tibetan Plateau climate changes have been few because of lack of data.

Data Sets and Methodology
With monthly average temperature and precipitation data of reference period (20 C3 M) of 22 IPCC-AR4 models under contemporary climate condition, this paper compares these data with observed values from 96 meteorological stations on the Tibetan Plateau (Table 2) and assesses generally the ability of GCMs to reproduce precipitation and temperature climatology on the Tibetan Plateau and climate projections under different emission scenarios, in order to make deeper understanding of China's climate changes in the context of sustaining global warming and provide basis for assessing impacts of climate changes on hydrology and water resources. Table 1 lists the basic information about 22 GCMs, as well as the border and number of grids of the Tibetan Plateau. More detailed information about GCMs can be found at http://www-pcmdi.llnl.gov/ipcc/about ipcc.php. Take GFCM21 pattern from American Geophysical Fluid Dynamics Laboratory, for example. Distribution of grids on the Tibetan Plateau is shown in Figure 1; small black dots indicate 96 meteorological stations.
Currently, the method of multimodel ensembles mean [10,12,23] is commonly used for analyzing global climate models data. Compared with a single mode, multi-model ensembles mean reduces errors and reproduces-a more realistic situation. However, with calculating the arithmetic average of multi-model data, characteristics of each model are disappearing in this method, which is not conducive to distinguish simulation features from each model. In order to assess the ability of each model to reproduce temperature and precipitation on the Tibetan Plateau, based on original grid data of monthly temperature and precipitation from GCMs, we have calculated regional monthly average of reference period for each model, analyzed differences between simulated and observed values in time and spatial scales, and taken relative error, absolute error, correlation coefficient, and deterministic coefficient as performance indeices of GCMs [24].  With regard to comparison between multiyear monthly mean temperature in reference period of various models and the observed data (Figure 2), it shows that the simulated values of all modes are smaller than the observed; however, the high ones (monthly mean temperature from June to August) are floating around observation. We selected five patterns with the best simulation ability ( Figure 3): GGMR, GFCM21, HADCM3, HADGEM, and MRCGCM. Figure 3 shows that climate changes with five models are consistent with the fact, which is also confirmed by Table 1. As seen from Table 1, these five models own a high value of uncertainty coefficient and the best simulation ability in annual mean temperature.

Temperature. As seen in
Based on grid output of each model, we have interpolated seasonal mean temperature for each station with Delta-DCSI downscaling method [25]. Traditionally, spring is from March to May; summer is from June to August; fall is from September to November; winter is from December to February of next year. Figure 4 shows the observed values of temperature on the Tibetan Plateau in reference period, which indicates that the average temperature of the Tibetan Plateau is between −5.6 and 14.5 ∘ C, with significant seasonal and spatial differences: high temperature in the southeast, decreasing from south to north and from east to west. Take GFCM21 as an example; we have analyzed differences between simulated and observed values of temperature in reference period (shown in Figure 5). Figure 5 shows that simulated temperature of GFCM21 pattern, which has a deviation of −18.1 ∘ C ∼ 7.2 ∘ C from the observed, is much lower than the observed in most part of the region except for edge of the northern and southern areas. Annual mean temperature has a deviation of −14.0 ∘ C ∼ 5.4 ∘ C. The maximum deviation has met −18.1 ∘ C which occurred in the western region in fall. The deviation was decreasing from west to east.

Precipitation.
The same method has been applied to assess the ability of each model to reproduce precipitation. Annual mean precipitation on the Tibetan Plateau was 4 The Scientific World Journal  492.5 mm in reference period which was lower than simulated values of most models ( Table 2). Different models had great differences in simulation ability while more than half of models had bad simulation performance. The maximum relative error is up to 155.8% with CNCM3, while the minimum is 6.0% with IPCM4. Some models have good correlation coefficient between simulated and observed values. There are more differences in precipitation simulation ability than in temperature of each model.
With regard to comparison between simulated and observed values of precipitation (Figure 6), most models have not well simulated annual changes of precipitation as the simulated values are higher than the observed. As most precipitation occurred in the period from June to September, only half of the models had higher simulated values than the observed. In general, CGMR, five patterns as CSMK3, GFCM20, GFCM21, and HADGEM, have well simulated annual precipitation trends (Figure 7), especially CSMK3, whose simulated monthly precipitation is very similar to the observation. Figure 8 shows observed values of precipitation on the Tibetan Plateau in reference period which indicates that precipitation concentrated in spring and summer and there was more precipitation in southeast. Annual mean precipitation was decreasing from south to north and from east to west. As seen from differences of simulated and observed precipitations of GFCM21 shown in Figure 9, the maximum deviation occurred in summer which equaled 13.7 mm/day, followed by spring and summer. With a deviation of −1.4 ∼ 2.0 mm/day, the precipitation simulation in winter is the best. For the whole region, deviation of average precipitation varied from −0.8 to 5.7 mm/day with a good simulation in the center of the region.

Projected Changes in Temperature and Precipitation
Prediction from these climate models of A1B, A2, and B1 scenarios has almost covered all possible climate changes caused The Scientific World Journal 5  by growing of greenhouse gas emissions on the Tibetan Plateau. Table 3

Conclusions and Discussion
In order to make further climate change projections under A1B, A2, and B1 emission scenarios on the Tibetan Plateau, temperature and precipitation simulation abilities of GCMs have been evaluated which is based on the differences between simulated and observed of reference period with 22 models from IPCC AR4. Some interesting conclusions can be presented and discussed as follows.
(1) 22 climate models have a certain capability to simulation temperature and precipitation on the Tibetan Plateau. The correlation coefficient of temperature of all the models (except INCM3 mode) is above 0.96, but there are still great differences in simulation performance of each model, while only GGMR, GFCM21, HADCM3, HADGEM, and MRCGCM patterns have relatively well simulated climate changes with an annual climate trend similar to the fact. Simulated precipitation of most models is higher than the observed values while the regional simulated values of some models are lower than the observed. However there are more differences between models in precipitation than in temperature. Five models namely, CGMR, CSMK3, GFCM20, GFCM21, and HADGEM have better simulated the precipitation on the Tibetan Plateau which indicates that simulation of most models need to be further improved.
(2) With a general assessment of the simulation ability of temperature and precipitation, it is obvious that GFCM21 and CGMR patterns can basically reproduce climate change on the Tibetan Plateau. Take GFCM21 pattern as an example; simulated temperature in summer is more similar to the observed. At the same time, there are significant spatial differences, while there is less deviation between simulated values of temperature and precipitation and the observed with GFCM21 model in the east part of the Tibetan Plateau.
(3) Under A1B, A2, and B1 emissions scenarios, there is great uncertainty in simulation of temperature and precipitation of the 21st century with climate models on the Tibetan Plateau. Temperature and precipitation of GCMs will all increase with different ratios. Precipitation under A2 scenario will increase most significantly while there is the least increase under A1B scenario. Based on evaluation of simulation performance with 22 climate models from IPCC AR4, there are great differences in the simulated values while only a few models have well simulated temperature and precipitation on the Tibetan Plateau. It is obvious that simulation with GCMs still needs to be further improved. Because of the low resolution of GCMs, parameters of physical process of climate models are also needed to be further improved, so there are great uncertainties in the simulation of temperature and precipitation. In the future, we could think of downscaling methods to predict climate changes in order to find proper method for China and provide more reliable prediction results with multimodel ensembles method. view of many uncertainties of climate models, it is of great importance to investigate climate changes in China with GCMs to reduce those uncertainties.