Responses of the Carbon Storage and Sequestration Potential of Forest Vegetation to Temperature Increases in Yunnan Province , SW China

The distribution of forest vegetation and forest carbon sequestration potential are significantly influenced by climate change. In this study, a map of the current distribution of vegetation in Yunnan Province was compiled based on data from remote sensing imagery from the Advanced Land Observing Satellite (ALOS) from 2008 to 2011. A classification and regression tree (CART) model was used to predict the potential distribution of the main forest vegetation types in Yunnan Province and estimate the changes in carbon storage and carbon sequestration potential (CSP) in response to increasing temperature. The results show that the current total forest area in Yunnan Province is 1.86 × 107 ha and that forest covers 48.63% of the area. As the temperature increases, the area of forest distribution first increases and then decreases, and it decreases by 11% when the temperature increases from 1.5 to 2 ◦C. The mean carbon density of the seven types of forest vegetation in Yunnan Province is 84.69 Mg/ha. The total carbon storage of the current forest vegetation in Yunnan Province is 871.14 TgC, and the CSP is 1100.61 TgC. The largest CSP (1114.82 TgC) occurs when the temperature increases by 0.5 ◦C. Incremental warming of 2 ◦C will sharply decrease the forest CSP, especially in those regions with mature coniferous forest vegetation. Semi-humid evergreen broad-leaved forests were highly sensitive to temperature changes, and the CSP of these forests will decrease with increasing temperature. Warm-hot coniferous forests have the greatest CSP in all simulation scenarios except the scenario of a 2 ◦C temperature increase. These results indicate that temperature increases can influence the CSP in Yunnan Province, and the largest impact emerged in the 2 ◦C increase scenario.


Introduction
Forest vegetation is a major terrestrial ecosystem and can provide low-cost options to mitigate climate change.However, climate change, especially temperature change [1], has significantly influenced the distribution [2], structure and ecology [3] of forest vegetation.Understanding the distribution of forest vegetation and assessing the present and future carbon balance of forest ecosystems are of scientific and policy interest.
Studies have explored the carbon storage and sequestration potential of forest vegetation in China with different models.Some used mature forests to estimate the carbon sequestration potential (CSP) of forest vegetation as a reference for the forests in China [4].Some utilized logistic equations to estimate the biomass carbon storage at a species level in China [5,6].Some used the data from the forest resource inventory to predict the forest vegetation carbon storage from 2005 to 2050 with a stage-classified matrix model [7].However, these approaches cannot be used to predict future whole-forest CSP, especially under climate change.In addition, the prediction results at the national scale are not suitable for prediction at smaller scales.Therefore, other studies have been performed to focus on the regional scale in China.There had researchers used the TRIPLEX model to simulate the forest growth and carbon dynamics of the boreal and temperate forest in Northeast China [8]; however, this model is limited for modeling vegetation succession and future dynamics.Some studies have explored the carbon storage and carbon sequestration in Yunnan Province [9][10][11].However, those studies have assumed that all forests are mature, and they have not considered the succession and distribution changes of forest vegetation.In addition, climate change causes high uncertainty, which is a key factor that regulates carbon sequestration in a forest.
Bioclimatic classification schemes are commonly used to predict changes in forest vegetation distribution caused by climate change [12].The most popular method is the Holdridge Scheme [13], but this scheme does not provide a biological interpretation of the influences of plant life forms and attributes [14].The dynamic vegetation model [15] considers the succession and physiological responses of forest vegetation to climate change but requires abundant data and many physiological parameters [16].General linear models (GLM), general addition models (GAM) and classification and regression (CART) [17] models have been extensively used to assess climate-vegetation relationships.Researchers compared these three types of predictive models for Fagus crenata forests (Fagus crenata Blume) in Japan and concluded that when accounting for the interactions of predictor variables, CART could explain the relationships between climate and forest vegetation distribution with greater accuracy than GLM and GAM [18].CART models are a practical technology that can be used to explore the relationships between environmental variables and vegetation [19], and they have been used to predict plant habitat distributions [20].
The Paris climate agreement aims to hold the global average temperature increase to well below 2 • C and pursues efforts to limit the temperature increases to 1.5 • C above preindustrial levels [21].Some previous studies have suggested that the global temperature will reach the 2 • C increase between 2026 and 2060 [22].Increasing temperatures affect forests by changing their geographical distributions [23] and exhibit dominant control over the natural distribution of forest ecosystems [24].Temperature increases will change the regeneration capacities [25] of species and therefore alter the area of suitable habitat of forest vegetation.Changes in the forest distribution area entail changes in the forest biomass and hence changes in the CSP.Predicting how forest distributions respond to ongoing and anticipated global warming is a challenge with great ecological relevance [26].
Forest vegetation is characterized by large biomass carbon stocks that are vulnerable to both biological and non-biological factors [27].Temperate and high-latitude forests have been shown to be carbon sinks; however, the carbon balance of subtropical forests has been less well studied [28].Achieving an understanding of the current and potential future role of forest sequestration is required by international negotiations and is very important for both managed and unmanaged forests [29].However, to date, few studies have investigated the potential impacts of further temperature increases on the distribution and carbon sequestration of natural vegetation.
Yunnan Province is one of the most climatically and biologically diverse areas in the world and has been noted to be sensitive to climate changes [2].The landform, climate, ecosystem, and species diversity of the province are higher than in any comparably sized region in China.In this study, the space-for-time method based on climax theory was adopted, which hypothesizes that forest vegetation ultimately reaches a climax status by succession [30].A spatial modeling approach based on a statistically derived bioclimatic factor is used to predict and understand the vegetation distribution and CSP of projected temperature increases within Yunnan Province.Our objectives are as follows: (1) determine the changes in forest vegetation distribution with temperature increases in Yunnan Province and (2) evaluate the responses of carbon storage and sequestration potential in forest vegetation in Yunnan Province to temperature increases.

Study Area
Yunnan Province is located in Southwest China between 21 • 08 32"-29 • 15 08" N and 97 • 31 39"-106 • 11 47" E. It is situated at the meeting point of three geographic regions: the eastern Asia monsoon region, the Tibetan Plateau region and the tropical monsoon region of southern Asia and Indo-China.The province has topographic complexity, with a large altitudinal range of 76 to 6740 m.Yunnan Province spans only 8 • in latitude but exhibits all of the climate zones and land ecosystem types of China [31].The climate is generally mild with a long growing period.There are cold winters at the higher elevations in the northwestern mountain regions, moderate temperatures in the middle plateau region, and tropical, hot and humid conditions at the lower elevations and valley bottoms in the southern region [32].Yunnan's rich ecosystems comprise over 30 ecosystem types according to the Chinese classification, which span from the lower tropical valleys and basins in the southern part to the barren high peaks and deep valleys in the northwestern part [33].The ecosystems of Yunnan Province are sensitive to environmental change and are less stable than those in temperate zones [34].

Climate and Forest Vegetation Data
We selected six bioclimatic factors from the WorldClim dataset: annual mean temperature (TMA, • C), mean temperature of the warmest quarter (TMS, • C), minimum temperature of the coldest month (TMW, • C), annual precipitation (PRA, mm), precipitation of the warmest quarter (PRS, mm) and precipitation of the coldest quarter (PRW, mm).The resolution of all datasets is 1 × 1 km.TMS is a measure of the effective heat required for plant growth.TMW is a measure of extreme cold, which controls the altitudinal and northern range limit of evergreen broad-leaved forests [35].PRS and PRW are measures of water supply during the growing and winter seasons, respectively.By performing paired samples t-tests, we confirmed that there were no significant differences between the WorldClim dataset and the dataset derived from 134 climate observation stations in Yunnan [36].
In this study, we investigated different climate change scenarios of air temperature due to anthropogenic influences based on the current climate background.Precipitation has not changed significantly over the past several decades in southwestern China [37].As the temperature increases in the future, the change in precipitation is not sure [38].Therefore, we considered the influence of temperature increases only.Here, we considered temperature (TMA, TMS and TMW) increases from 0 to 2 • C at 0.5 • C intervals with no changes in precipitation in each scenario.
We set up a vegetation database by combining field investigations with remote sensing identification.We drew a vegetation map based on remote sensing images from the Advanced Land Observing Satellite (ALOS) from 2008 to 2011 (Table 1), in the Albers equal-area conic projection and Krasovsky datum.We used the 1:50,000 topographic map for geometric correction of remote sensing images.We were building the 1 × 1 km grid to resample vegetation data which was drawn from ALOS remote sensing by ArcGIS 10.2 (ESRI Inc., Redlands, CA, USA).The longitude and latitude records of the cells were taken from this map, which yielded the same spatial resolution as the climate data.A 10 × 10 km grid and global positioning system (GPS) point layer were built for accuracy tests.The Beijing 54 coordinates were adopted in all layers and subsequent processes.The 1 × 1 km spatial resolution was used to combine these data.The accuracy of the two methods was tested using the 10 × 10 km geographic grid point and the GPS point that was recorded in the field investigation and not used for interpretation.When assessing the accuracy by using the 10 × 10 km grid, we evaluated whether the altitude and image characteristics of the point matched the vegetation type.When testing the accuracy by using the GPS point, we determined whether the vegetation type was consistent with that recorded within a 0.25 km buffer area of the point.The validation results indicated that the vegetation map was highly accurate [39,40].Thus, the vegetation map was used for the geographic distribution of vegetation data in Yunnan Province.In this study, we chose 7 dominant vegetation types, which are shown in Table 2.

Description of the CART Model
A CART model is built using the binary recursive partitioning method, and the results are a simple binary tree structure.Each branch of the binary tree represents a test result.CART models are a practical technology that can be used to explore the relationships between environmental variables and vegetation [41] and have been widely used to predict the distribution of plant habitat [20].Minimum Gini coefficient values are used to evaluate the attributes of the model.A smaller Gini value indicates a higher "purity" of the samples and better divided results.The CART algorithm builds a tree before pruning, and the accuracy is often higher than that of the multi-tree algorithms because the binary tree cannot easily generate data fragments.Therefore, the CART model uses binary recursive partitioning and utilizes Boolean tests at the branch nodes.If a condition is met, then the sample is divided into the left branch; otherwise, the sample is divided down the right branch, and a binary decision tree is eventually formed.
We obtained the spatial vegetation data by encoding the forest vegetation and 1 × 1 km grid samples.We established the CART model by using the "Tree" program in R3.3.3 software (University of California, Berkeley, CA, USA).All climate-vegetation sampling data were divided into two data sets, one for validation and another for testing.The model was built using 70% of the data that was randomly selected; the rest of the data were used for pruning.Overfitting may occur during the modeling process; therefore, we used cross-validation rules to eliminate the least important tree using the "prune tree" function on the final model.
The area under the curve (AUC) derived from receiver operating characteristic (ROC) analysis was calculated to validate the performance of the CART model.The AUC values were interpreted for model accuracy using the following standards: 0.90-1.00:excellent, 0.8-0.90:good and 0.7-0.80:fair [42].We divided the habitats into three categories: non-habitats, marginal habitats and suitable habitats, and the occurrence probability thresholds were obtained from the ROC analysis [43].The areas where the predicted probability of occurrence was less than a low occurrence probability (0.01) were defined as non-habitats; the areas where the probability was equal to or greater than 0.01 but smaller than the optimal threshold were defined as marginal habitats; and the areas where the probability was greater than the optimal threshold were defined as suitable habitats.In this study, we selected the suitable habitats as potential habitats.It has been proved that climate factors can be used to build CART models [36].

Calculation of CSP
Using the Forest Identity concept [44], the area of forest vegetation (S, ha), the mean aboveground biomass (W a , Mg/ha) (Table 3) and the forest vegetation biomass (W, TgC) can be linked using Equation (1).The forest carbon storage was based on the biomass in the study area and was calculated according to the different forest carbon sequestration rates.We used a conversion coefficient of 0.5 [52] (the carbon content of per gram of dry matter) to calculate the forest carbon storage.We obtained the forest carbon storage by multiplying the carbon density by the forest vegetation distribution area.Suitable habitat and actual carbon storage were calculated by the following formula: where W c is forest carbon storage (TgC), W is the vegetation biomass of the forest (TgC), and 0.5 is the conversion coefficient.We simulated the habitat area suitable for forest vegetation by using the CART model.We obtained the carbon sequestration potential (CSP) by subtracting the carbon storage of suitable habitat from the actual carbon storage: where W p is the CSP (TgC), W sui is the carbon storage of suitable habitat (TgC), and W act is the actual carbon storage (TgC).

Prediction Accuracy and Contribution of Climate Variables
The AUC values are greater than 0.8 and less than 0.9 for all simulation scenarios.Deviance-weighted scores (DWS) were applied to evaluate the contributions of each predictor variable to the model (Table 4).TMW was an overwhelmingly potent factor among the six climate variables, indicating that extremely cold temperatures in a year play a decisive role in the broad-scale distribution of the forest vegetation in Yunnan Province.TMS, PRS and PRW also affected the distribution to a small extent.TMA and PRA showed no contribution; thus, these variables are not shown in Table 4.

Distribution Area of Current Forest Vegetation
The map of the vegetation in Yunnan Province (Figure 1) shows that the main forest vegetation area in Yunnan Province was 1.86 × 10 7 ha and that the forest coverage was 48.63%.The TCC and CTC types had small distribution areas and coverage.They were mainly distributed in northwestern Yunnan.The areas of TCC and CTC were 7.10 × 10 5 ha and 1.35 × 10 6 ha, respectively, and accounted for 1.85% and 3.52% of Yunnan's land area, respectively.WHC was mainly distributed in Pu'er District and had an area of 1.94 × 10 6 ha.MHEB had the smallest area (3.79 × 10 5 ha) among the forest vegetation types and was mainly distributed in mountainous areas with higher elevations in Yunnan; MHEB accounted for 0.99% of Yunnan's total land area.MEB was mainly distributed in southern Yunnan and had an area of 3.16 × 10 6 ha.The SEB and WTC forests were intermixed and had the largest distributions, being spread throughout the central region of Yunnan.The areas of SEB and WTC were 1.47 × 10 6 ha and 9.62 × 10 6 ha, respectively, and accounted for 3.84% and 25.11% of the land area, respectively.

Potential Forest Vegetation Distribution
Under the T0.0 scenario, the forest vegetation ultimately reaches a climax status by succession.The entire forest area in Yunnan Province increased by 130.42% compared with the current forest vegetation distribution area.The vegetation distribution areas of the 7 selected forest types all increased under this scenario, especially WHC, which increased by 572.30%.As the temperature increased, the distribution area of all forest vegetation types (except WTC) decreased to varying extents (Table 5).The distribution area of WTC increased by approximately 30% in all of the simulation scenarios.When the temperature increased from 0.5 °C to 1.0 °C , the distribution area of MHEB decreased significantly.When the temperature increased from 1.5 °C to 2.0 °C , the distribution area of WHC decreased significantly.SEB was the forest type most sensitive to temperature, with its distribution area decreasing by more than half with a temperature increase of 0.5 °C (Figure 2).In general, the mean forest distribution area first increased and then decreased with increasing temperature.When the temperature increased by 2 °C , the mean forest distribution area in Yunnan decreased by 11% (Table 5).

Potential Forest Vegetation Distribution
Under the T0.0 scenario, the forest vegetation ultimately reaches a climax status by succession.The entire forest area in Yunnan Province increased by 130.42% compared with the current forest vegetation distribution area.The vegetation distribution areas of the 7 selected forest types all increased under this scenario, especially WHC, which increased by 572.30%.As the temperature increased, the distribution area of all forest vegetation types (except WTC) decreased to varying extents (Table 5).The distribution area of WTC increased by approximately 30% in all of the simulation scenarios.When the temperature increased from 0.5 • C to 1.0 • C, the distribution area of MHEB decreased significantly.When the temperature increased from 1.5 • C to 2.0 • C, the distribution area of WHC decreased significantly.SEB was the forest type most sensitive to temperature, with its distribution area decreasing by more than half with a temperature increase of 0.5 • C (Figure 2).In general, the mean forest distribution area first increased and then decreased with increasing temperature.When the temperature increased by 2 • C, the mean forest distribution area in Yunnan decreased by 11% (Table 5).

Potential Forest Vegetation Distribution
Under the T0.0 scenario, the forest vegetation ultimately reaches a climax status by succession.The entire forest area in Yunnan Province increased by 130.42% compared with the current forest vegetation distribution area.The vegetation distribution areas of the 7 selected forest types all increased under this scenario, especially WHC, which increased by 572.30%.As the temperature increased, the distribution area of all forest vegetation types (except WTC) decreased to varying extents (Table 5).The distribution area of WTC increased by approximately 30% in all of the simulation scenarios.When the temperature increased from 0.5 °C to 1.0 °C , the distribution area of MHEB decreased significantly.When the temperature increased from 1.5 °C to 2.0 °C , the distribution area of WHC decreased significantly.SEB was the forest type most sensitive to temperature, with its distribution area decreasing by more than half with a temperature increase of 0.5 °C (Figure 2).In general, the mean forest distribution area first increased and then decreased with increasing temperature.When the temperature increased by 2 °C , the mean forest distribution area in Yunnan decreased by 11% (Table 5).

Carbon Storage and CSP
The total current actual carbon storage of the seven types of forest vegetation in Yunnan Province was 871.14 TgC.Among the forest types, MEB and WTC had the highest levels of actual carbon storage due to their large distribution areas, with values of 205.42 TgC and 172.72 TgC, respectively.Although MHEB has the highest mean biomass, it had lower actual carbon storage, with a value of 75.93 TgC.TCC had a high mean biomass, but because of its small distribution area and coverage, it had the lowest value (56.40 TgC) of actual carbon storage among the forest types.In contrast, WTC had a low mean biomass but a wide distribution, resulting in a high value of actual carbon storage.SEB, WHC and CTC had intermediate levels of actual carbon storage, which were 99.91 TgC, 137.78 TgC and 122.98 TgC, respectively.The total actual carbon storage in coniferous forests was 489.88 TgC, accounting for 56.23% of the total.
Overall, the CSP of the seven types of forest vegetation in Yunnan Province increased and then decreased as the temperature increased.In particular, the CSP decreased sharply from 1059.19 TgC to 647.24 TgC when the temperature increased by 2 °C .On the whole, the forest vegetation exhibited the largest CSP (1114.82TgC) when the temperature increased by 0.5 °C .The different forest vegetation types in Yunnan Province showed varying increases in carbon sink ability as the temperature increased except SEB, which exhibited a negative CSP value.Among the forest types, WHC and WTC exhibited the largest CSP values when the temperature increased.The results show

Carbon Storage and CSP
The total current actual carbon storage of the seven types of forest vegetation in Yunnan Province was 871.14 TgC.Among the forest types, MEB and WTC had the highest levels of actual carbon storage due to their large distribution areas, with values of 205.42 TgC and 172.72 TgC, respectively.Although MHEB has the highest mean biomass, it had lower actual carbon storage, with a value of 75.93 TgC.TCC had a high mean biomass, but because of its small distribution area and coverage, it had the lowest value (56.40 TgC) of actual carbon storage among the forest types.In contrast, WTC had a low mean biomass but a wide distribution, resulting in a high value of actual carbon storage.SEB, WHC and CTC had intermediate levels of actual carbon storage, which were 99.91 TgC, 137.78 TgC and 122.98 TgC, respectively.The total actual carbon storage in coniferous forests was 489.88 TgC, accounting for 56.23% of the total.
Overall, the CSP of the seven types of forest vegetation in Yunnan Province increased and then decreased as the temperature increased.In particular, the CSP decreased sharply from 1059.19 TgC to 647.24 TgC when the temperature increased by 2 • C. On the whole, the forest vegetation exhibited the largest CSP (1114.82TgC) when the temperature increased by 0.5 • C. The different forest vegetation types in Yunnan Province showed varying increases in carbon sink ability as the temperature increased except SEB, which exhibited a negative CSP value.Among the forest types, WHC and WTC exhibited the largest CSP values when the temperature increased.The results show that incremental warming of 2 • C will sharply decrease forest carbon sequestration in Yunnan Province.Much of the observed decrease was due to WHC, which showed a decrease in CSP from 734.52 TgC to 234.96 TgC (Table 6).

Model Accuracy Test
The AUC values of all models were greater than 0.8 and less than 0.9 in all simulation scenarios (Table 5).The high accuracy of the model suggested that the distribution of vegetation at the province scale can be explained by climatic variables.The accuracy of the carbon density estimate is very important and should be discussed.Although future vegetation distributions were predicted by the CART model under a rising temperature scenario, the present carbon density values were used to calculate the future carbon storage in Yunnan Province in this study.However, this method is not entirely accurate because carbon density in the future will change under future climate conditions [53,54].Carbon sequestration depends on both vegetation composition and the climate [53].Although the past and future vegetation compositions may be similar to the present composition, they may have different carbon densities under different future climate conditions [55].Hence, the application of the modern carbon density database [56] might underestimate or overestimate past and future terrestrial carbon storage values.Thus, changes in carbon density under future climate conditions could have a significant impact on the estimates of carbon storage in Yunnan Province.In this sense, new approaches, such as dynamic global vegetation models (DGVMs) that couple biogeography and biogeochemistry models are encouraging [57] because they consider the effects of changes in both the climate and atmospheric CO 2 concentrations on both vegetation redistribution and carbon density [53].

The Status of Forest Vegetation Carbon Storage in Yunnan
The carbon density of the forest vegetation in Yunnan Province was 84.69 Mg/ha, which was twice the average carbon density in China (41.32 Mg/ha) [58], similar to the worldwide average carbon density (86.00 Mg/ha) [59] and greater than the subtropical forest carbon density (66-77 Mg/ha) [60].These results can be interpreted as follows.First, the forest cover in Yunnan Province is high and mostly consists of mature coniferous forest vegetation [4].Furthermore, the International Biosphere Project (IBP) may have overestimated worldwide forest carbon storage because of an insufficient number of sampling sites [61].
The vegetation carbon storage in the forests of Yunnan Province was approximately 871.14 TgC, which was 14.9% of the total forest vegetation carbon storage in China (5.85 PgC;Fang et al., 2007) and 0.24% of the worldwide forest vegetation carbon storage (360 PgC; Pan et al., 2011).The main reason for this difference was that the other studies used the national forest inventory data to estimate carbon storage, which led to the underestimation of the forest vegetation distribution area.In addition, differences in vegetation classification schemes and estimation methods result in different carbon storage estimates.The forest vegetation in Yunnan Province is a large carbon sink in China and plays an important role in the world.

The Effects of Temperature Increases on Forest Vegetation Carbon Storage and CSP
The successful simulation of the distribution of forest vegetation under different climate conditions achieved here is partly attributable to the accurate simulation of the relative distribution of forest carbon.Previous studies have indicated that temperature is the main factor that affects forest vegetation carbon storage in China, the Midwestern United States [62], Russia [63], Canada, and the Netherlands [64], among other places.The effect of temperature differs among places and vegetation types.Temperature increases have been found to increase forest carbon storage in colder and wetter ecoregions [65] but reduce the rates of net above-ground biomass increases in the Amazon rainforest [66] and the growth rates of mature rainforests [67].Rising temperatures affect forest carbon storage in two ways.First, rising temperatures alter vegetation types [68] and vegetation boundaries [69] and can transform coniferous forest regions into broad-leaved forest regions [68].Second, rising temperatures increase species diversity [70], which reduces the sensitivity to temperature and influences CO 2 and energy exchange [71].A temperature increase of 4 • C will increase the absorption and saturation of CO 2 in forest vegetation [72].Forests respond positively to the influences of rising temperatures [73,74], and temperature increases have positive effects on forest growth and wood production in the short-to-medium term [75].In turn, stand age has an influence on forest carbon storage.With increasing stand age, the carbon storage of forest increases quickly, then reaches a maximum value [76].At last, the carbon storage of forest decreases to a relatively stable level because of the limit of hydraulic resistance [77], and the growth of wood becomes extremely slow or is almost not change [78].Future studies can consider further the effect of stand age.
In Yunnan Province, the forest distribution area changed by 5.74%, 3.30%, 5.31% and −11.10% when the temperature increased by 0.5, 1.0, 1.5 and 2.0 • C, respectively.The total amount of carbon sequestration in the worldwide forest is approximately 360 PgC [32], and the CSP of mature forest constitutes 49.32 to 58.37% of total forest CSP [78].Among forest types, subtropical forests have higher CSP [4] but are more sensitive to temperature changes.The carbon sequestration rate of coniferous forests is more impacted by climate change than is that of deciduous forests [79] in subtropical regions.Much of the coniferous forest vegetation in Yunnan Province is mature vegetation [4] and is easily affected by rising temperatures (Table 5).Temperature increases from 0 • C to 1.5 • C did not have severe impacts on the CSP of the forests in Yunnan Province.However, a temperature increase of 2 • C resulted in sharp decreases in the CSP of coniferous forest vegetation (Table 6).This result indicated that the CSP of mature forest vegetation was more easily influenced by temperature increases than broad-leaved forest vegetation.Mature coniferous forests are mainly distributed in the Hengduan Mountains of northwestern Yunnan Province.In this region, the mean annual temperature over the last two decades has increased at a rate of 0.6 • C/10 year, and the vegetation distribution has changed in the past [80].In conclusion, rising temperature has impacted the CSP in Yunnan Province by affecting the mature forest vegetation, especially coniferous forest vegetation.In this paper, we hypothesize only that the forest is balanced with current climate conditions, and the forest will develop to the climax.At last, the places where are uitable for forest growth will full of this kind of forest vegetaiton in turn.With the development of forest vegetation, the CO 2 will increase and have fertilization effects [81], revealing underestimates of forest carbon storage in our projections.However, nitrogen limitation is a driving factor of the forest carbon storage responses to elevated CO 2 [82].It suppresses the positive vegetation response to enhanced CO 2 fertilization [83], and this same limitation effect has been observed in a grassland ecosystem [84].Other problems are anthropogenic effects, because the need for cropland will lead to over-estimates of the carbon sink and coupling effects of temperature with other climate factors.In this study, the influence of temperature change was taken into consideration, and there was a gap between it and the actual effect.Future studies can consider the comprehensive effect.

The CSP of Warm-Hot Coniferous Forest
The warm-hot coniferous forest had the highest CSP in Yunnan Province under each simulation scenario except the temperature increase of 2 • C (Table 6).These results are similar to those obtained in India, that noted that the total forest biomass and carbon pool of Pinus kesiya forest were greater than those of other pine forests studied in other regions of the world [85].In the present study, the forest distribution area was found to change by −5.83%, −4.89%, −5.83% and −59.76% when the temperature increased by 0.5, 1.0, 1.5 and 2.0 • C, respectively.The main species composition was P. kesiya which can occupy a variety of habitats and is mainly distributed in Southwestern Yunnan Province [86].In the distribution area, the annual precipitation ranges from 1000 to 1500 mm, the relative humidity reaches 80%, and the elevation ranges from 600 to 1950 m [87].P. kesiya trees have rapid growth [88] and strong natural regeneration ability.The average carbon sequestration rate of P. kesiya is 12.7 Mg C ha −1 yr −1 in the Philippines [88].If P. kesiya and Pinus yunnanensis share the same habitat, the two species exhibit converging development [87].The emission-reduction effect of the P. kesiya afforestation project in Yunnan Province has been pronounced [89].

Conclusions
Classification and regression tree models are capable of modeling the distribution of forest vegetation, especially suitable habitat and the carbon sequestration potential under changing climatic conditions.Based on the simulated results, several conclusions can be drawn.The carbon sequestration potential of coniferous forests is more strongly influenced by temperature increases than the carbon sequestration potential of deciduous forests.Temperature increases can influence the carbon storage and carbon sequestration potential in Yunnan Province.Warm-hot coniferous forests have a large current carbon storage and carbon sequestration potential, and this vegetation type is especially sensitive to temperature increases of 2 • C compared with other forest vegetation types.However, warm-temperate coniferous forests have the largest distribution area, and the carbon sequestration potential will increase when temperatures increase by 2 • C. The carbon sequestration of semi-humid evergreen broad-leaved forests is most sensitive to temperature increases, and it will decrease when temperatures increase by 0.5 • C. Overall, the forest vegetation in Yunnan Province has high carbon density, carbon storage and carbon sequestration potential.Temperature increases of 2 • C will sharply decrease the carbon sequestration potential of the forests in Yunnan Province.

Figure 1 .
Figure 1.The distribution of forest vegetation in Yunnan Province in the 2010s.

Figure 1 .
Figure 1.The distribution of forest vegetation in Yunnan Province in the 2010s.

Forests 2018, 9 , 15 Figure 1 .
Figure 1.The distribution of forest vegetation in Yunnan Province in the 2010s.

Table 1 .
Waveband parameters of the ALOS remote sensing images.

Table 2 .
Vegetation types and representative tree species of the study.

Table 3 .
The mean forest vegetation biomass of seven forest types in Yunnan Province.

Table 4 .
DWS and the percentage of each climatic variable and AUC of each simulation scenario.

Table 5 .
The area and rate of change in forest distribution area under different temperature scenarios.

Table 5 .
The area and rate of change in forest distribution area under different temperature scenarios.Units: Area (×10 5 ha); Rate of change (%).