The Effects of Asymmetric Diurnal Warming on Vegetation Growth of the Tibetan Plateau over the Past Three Decades

Temperatures over the past three decades have exhibited an asymmetric warming pattern between night and day throughout the Tibetan Plateau. However, the implications of such diurnally heterogeneous warming on vegetation growth is still poorly understood. In this paper, we evaluate how vegetation growth has responded to daytime and night-time warming at the regional, biome, and pixel scales based on normalized difference vegetation index (NDVI) and meteorological data from 1982 to 2015. We found a persistent increase in the growing seasonal minimum temperature (Tmin) and maximum temperature (Tmax) over the Tibetan Plateau between 1982–2015, whereas the rate of increase of Tmin was 1.7 times that of Tmax. After removing the correlations between Tmin, precipitation, and solar radiation, we found that the partial correlation between Tmax and NDVI was positive in wetter and colder areas and negative in semi-arid and arid regions. In contrast, the partial correlation between Tmin and NDVI was positive in high-cold steppe and meadow steppe and negative in montane steppe or wet forest. We also found diverse responses of vegetation type to daytime and night-time warming across the Tibetan Plateau. Our results provide a demonstration for studying regional responses of vegetation to climate extremes under global climate change.


Introduction
As one of the key environmental factors impacting the spatial and temporal distribution of vegetation ecosystems, temperature changes directly alter the vegetation growth environment.This affects the vegetation growth dynamics and vegetation structure and function [1].The fifth Assessment Report of the United Nations Intergovernmental Panel on Climate Change (IPCC) noted that the global average temperature is 0.85 • C higher than of 1880.Additionally, the last 30 years may be the warmest 30 years in the northern hemisphere in the past 1400 years [2].Meanwhile, the rate of global warming is heterogeneous across spatial scales.For example, the rate of warming for coasts, mountains, and high latitudes is relatively high [3].Global warming also exhibits temporal asymmetry.
For example, the rate of increase of global surface daily minimum temperatures (T min ) is 1.4 times higher than that of the daily maximum temperature (T max ) [3].This asymmetric day-night warming trend will have a major impact on the carbon uptake and carbon consumption of vegetation globally.The most critical stages of photosynthesis occur primarily during the daytime so the process is more sensitive to changes in T max .However, respiration can occur day and night so it is sensitive to both T max and T min [4].While increases in T min enhance both the respiration of vegetation and carbon consumption, increased photosynthesis during the daytime and prolonged growing seasons elevate carbon sequestration.Therefore, it is difficult to predict the response of ecosystem structure, processes, and functions to global warming when only considering daily mean temperature and daytime temperature.Accordingly, it is necessary to further examine the asymmetric effects of daytime and night-time warming on natural ecosystems.
Generally, the method for assessing the vegetation response to climate changes can be grouped into three categories, which include controlled experiments, model simulations, and quantitative remote sensing surveys.Controlled experiments include the transplantation of undisturbed soil vegetation to plots along different altitude gradients, artificial warming experiments, and soil heating experiments [5][6][7].Although controlled experiments can directly reveal the asymmetric effects of warming on vegetation, it is difficult to quantitatively assess such responses at a large scale and across different ecosystems owing to the restriction of manpower and material resources.Model simulation can use quantitative models to simulate the impact of historical warming periods on vegetation and to forecast the influence of future global warming on vegetation [8,9].However, asymmetric diurnal warming is not currently taken into account by many global carbon cycle models, which affects the accuracy of regional and global models.Quantitative remote sensing methods can quantitatively evaluate large-scale responses of vegetation to T max and T min based on time series of remote sensing data [4,10].Additionally, previous studies have shown that the response of vegetation to changes in T max and T min may differ between regions and ecosystem [4,[10][11][12][13][14].Therefore, it is vital to employ a spatially explicit and quantitative tool to assess the effects of daytime and night-time warming to enable better ecosystem management and adaption.
Satellite-based measurements can provide this type of broad perspective for identifying the spatial heterogeneity of land surface dynamics over the past three decades and to discern the spatiotemporal patterns of vegetation responses to asymmetric diurnal warming [4].The Normalized Difference Vegetation Index (NDVI), which is related to green leaf biomass and vigor, has proven to be a robust indicator of vegetation growth [15].The Global Inventory Modeling and Mapping Studies (GIMMS) NDVI dataset, which is the most complete and longest remote sensing dataset, has been used successfully in long-term monitoring studies of vegetation dynamics and climate change impacts on vegetation at regional and global scales [13,16].Therefore, in order to better understand, simulate, and forecast the response of terrestrial ecosystems to global changes, it is essential that research on the responses of terrestrial vegetation to global warming is based on remote sensing time series data, especially when examining areas that are sensitive to climate change.
As the third pole of the Earth, the Tibetan Plateau, has a unique pattern of natural and hydrothermal spatial differentiation that exhibits gradual change in climate from its warm and humid southeastern regions to its cold and dry northwest.The plateau contains rich, diverse, and fragile vegetation ecosystems and is a natural laboratory for climate change and earth science research [17].In recent decades, the warming of the Tibetan Plateau has been accompanied by a significant decrease in the diurnal temperature difference, which has been mainly caused by the increase of the minimum temperature in winter and at night.The rate of increase of the land-surface T min over the past five decades is twice that of T max [18].Warming effects on the Tibetan Plateau have a significant impact on both the dynamics of the vegetation during the growing season as well as a profound impact on the plateau's geographical and ecological patterns.Despite an increasing recognition of the critical importance of Tibetan Plateau vegetation for water resource conservation, for biodiversity protection, for land degradation, and as a carbon source and sink [17], the relationships between vegetation changes and climate warming are not firmly established [19,20].
Additionally, the climate and vegetation across the physico-geographical zone is relatively uniform.Therefore, research on the response of vegetation to climatic elements in different physico-geographical sub-regions can better reveal the characteristics of the regional differences across the plateau.Meanwhile, there are few studies that have compared the different physico-geographical zones to determine the vulnerability of various vegetation types to climate change.Limited efforts have been made to investigate the relative role of daytime and night-time warming in response to different vegetation types.Accordingly, we analyzed how vegetation growth trends have varied across different physico-geographical zones and evaluated the vegetation response to daytime and night-time warming across regions, biomes, and vegetation types.
The objective of this study was to determine the spatiotemporal responses of the terrestrial vegetation on the Tibetan Plateau to daytime and night-time warming.Based on the GIMMS NDVI3g data from 1982 to 2015, we investigated where and to what extent vegetation growth trends on the Tibetan Plateau have been affected by daytime and night-time warming.We also investigated information on the climatic factors, physico-geographical zones, and land cover types.These results could provide insights into the regions of the Tibetan Plateau that are most vulnerable to environmental change and, therefore, support sustainable land management.

Data Sources and Processing
NDVI is a sensitive indicator of vegetation productivity, vegetation biomass, and fractional vegetation coverage.It has been widely used to quantify the processes, status, and trends of vegetation growth [21].This study used the third-generation of NDVI dataset (NDVI3g), which was produced by the GIMMS group based on NOAA/AVHRR series satellites (NOAA-7, 9, 11, 14, 16, and 17).The data at a spatial resolution of 1/12 • in 15-day time intervals during 1982-2015 were provided by the Ecological Forecasting Lab at the NASA Ames Research Center (http://ecocast.arc.nasa.gov/).The dataset has been corrected to reduce noise resulting from sensor degradation, orbital drift, solar zenith angles, volcanic eruptions, sensor shift among the six satellites, and other extraneous factors that have exerted nonlinear and non-stationary effects on the data [19,22].To further mitigate the contamination of the clouds and the atmospheric turbulence, we obtained the monthly NDVI3g dataset using the Maximum Value Composite (MVC) method.Growing seasonal NDVI (April-October) was defined as the average monthly composite NDVI during the growing season, which has been widely used to assess the photosynthetic activities of vegetation [4,10,23].In order to minimize the influence of soil variation on spectral signal in bare and sparsely vegetated zones, we only considered the pixels with an average growing-season NDVI greater than 0.1 [13,24].
Based on the time period of GIMMS NDVI3g data, T max , T min , precipitation, and solar radiation were collected from 135 meteorological stations across and around the Tibetan Plateau (see Figure 1).These data were provided by the China Meteorological Data Sharing Service System (http://cdc.cma.gov.cn).In terms of the geographical coordinates and altitude of each site, these meteorological data except for the solar radiation data were interpolated into a raster dataset at an 8 km resolution by using the thin plate smoothing splines interpolation method, which incorporated topographic effects on spatial climate interpolation [25].As the solar radiation data were obtained from a few stations, it was difficult to determine reliable patterns of solar radiation using spatial interpolation methods based on ground-based solar radiation observations.Therefore, we used the products of three hourly gridded downward shortwave radiations, which improved the accuracy using estimates by both the surface data-based model and the satellite data [26].The data were provided by Cold and Arid Regions Science Data Center at Lanzhou, China (http://westdc.westgis.ac.cn/data).The maximum and minimum air temperatures during the growing season were calculated as the average of the daily maximum and minimum air temperature in the same days over years, respectively.The growing season precipitation and radiation were defined as the accumulated precipitation and radiation for the corresponding days, respectively.
The digital elevation model (DEM) data at a 250-m spatial resolution were provided by the NASA Shuttle Radar Topographic Mission (SRTM) website (http://www.glcf.umd.edu/; Figure 1).The DEM data were resampled at an 8 km spatial resolution to match the GIMMS NDVI3g data by a cubic spline method.Vegetation distribution information was obtained from a digitized 1:1,000,000 vegetation map of China [27].The main vegetation type was reclassified as forest, shrubland, cropland, alpine meadow, alpine steppe, and others type (see Figure 1).The vegetation map was resampled to match the spatial resolution of GIMMS NDVI3g data using a majority filter method.There are 11 physico-geographical zones of the Tibetan Plateau [28] including IB1 or Golog-Nagqu high-cold shrub-meadow zone, IC1 or Southern Qinghai high-cold meadow steppe zone, IC2 or Qangtang high-cold steppe zone, ID1 or Kunlun high-cold desert zone, IIAB1 or Western Sichuan-eastern Tibet montane coniferous forest zone, IIC1 or Southern Tibet montane shrub-steppe zone, IIC2 or Eastern Qinghai-Qilian montane steppe zone, IID1 or Ngari montane desert-steppe and desert zone, IID2 or Qaidam montane desert zone, IID3 or Northern slopes of Kunlun montane desert zone, and OA1 or Southern slopes of Himalaya montane evergreen broad-leaved forest zone.In order to reduce error and the uncertainty of the results, we used the same coordinate systems for the GIMMS NDVI3g data, climate data, DEM data, vegetation map, and physico-geographical zone map of the Tibetan Plateau, i.e., World Geodetic System 1984 (WGS84).

Linear Regression Trends
In order to analyze the changes of T max (and T min ) across the Tibetan Plateau, a linear regression was fitted using Equation (1).The statistical significance of the trend called the slope was evaluated based on the correlation coefficient in the regression equation.
where y denotes the T max (or T min ) for year x and e is the deviation of the data from the straight line defined by the intercept a and slope b.Positive and negative values of the correlation coefficient b represent increases or decreases, respectively, in T max (or T min ).The coefficients a and b were determined using least-squares fitting.

Partial Correlation Analysis
Partial correlation analyses were used to detect the relationship between each dependent variable and each particular independent variable by excluding the confounding effects of other variables [4,10,29,30].The absolute values of each partial correlation coefficient were categorized as greater than or equal to zero and less than or equal to one.When the absolute value of a partial correlation is equal to one, the two variables are completely related.In contrast, when a partial correlation is equal to zero, the two variables are completely unrelated.Since precipitation and solar radiation also influence NDVI [9], variation in precipitation and short solar radiation were taken into account in the partial correlation analysis.
The statistical significance of the partial correlation coefficient between growing season NDVI and T max (or T min ) after controlling for T min (or T max ), precipitation, and short solar radiation were evaluated by the Student's t-test as seen in Equation (2).
where n, q, and r are the total number of years, the number of independent variables, and the partial correlation coefficient between growing-season NDVI and T max (or T min ) after controlling for T min (or T max ), precipitation, and short solar radiation.

Trends of Daytime and Night-Time Warming in Tibetan Plateau
The increasing trend in T min over the past three decades is 1.7 times that in T max over the Tibetan Plateau during the growing season (see Figure A1).From 1982 to 2015, the warming trend of T max and T min has occurred over 94.2% and 99.8% of the total study area, respectively (Figure 2a,b).Additionally, T max and T min show statistically significant warming across 68.1% (p < 0.05) and 95.9% (p < 0.05) of the Tibetan Plateau, respectively.T max in the northwest and south of the Tibetan Plateau showed a statistically significant downward trend and the declining area occupied only 5.8% of the study area.
Additionally, the warming trend and increasing amplitude for T max and T min were not consistent across the Tibetan Plateau. Figure 2c shows the difference between T max and T min (DTR) on the Tibetan Plateau during 1982-2015.The trend in DTR is heterogeneous across the Tibetan Plateau.DTR decreased in 69.8% of the Tibetan Plateau particularly in the Midwest portion of the plateau (see Figure 2c) of which 43.5% was statistically significant (p < 0.05).In contrast, DTR increased over about 30.2% of the total study area of which 10.6% was statistically significant (p < 0.05).

Partial Correlations between NDVI and Diurnal Extreme Temperature on the Tibetan Plateau
When the effects of growing season T min , precipitation, and solar radiation were removed from the partial correlation, the individual effect of growing season T max inter-annual changes on inter-annual NDVI was obtained, which revealed a remarkable spatial pattern (see Figure 3a).Overall, 74.4% of the area exhibited a positive correlation between NDVI and T max and 18.5% of the region showed a significantly positive correlation (p < 0.05; Figure 3a).These significant positive correlations between NDVI and T max were found mainly within the wetter or colder regions including IB1, IC1, IIC2, OA1, and the western portion of IIAB1 (see Figure 3a).In addition, 25.6% of the area exhibited a negative correlation between NDVI and T max with 1.1% of the total region showing a statistically significant negative correlation between these variables (p < 0.05).The negative correlations between NDVI and T max occurred across the semi-arid and arid region including IIC1, IID3, IID2, and the western portion of IC2.
While the effects of growing-season T max , precipitation, and solar radiation were removed in the partial correlation analysis, the individual effect of growing season T min interannual changes on interannual NDVI was also obtained (see Figure 3b).While 46.6% of the region exhibited positive correlations between NDVI and T min during the growing season, 7.8% of the total area exhibited significant positive correlations (p < 0.05; Figure 3b).These statistically significant positive correlations between NDVI and T min were found mainly in the high-cold steppe and meadow steppe zones (i.e., IC1, IIC1, IC2, IID3, and ID1).In contrast, 53.4% of the region showed negative correlations between NDVI and T min during the growing season with 8.3% of the total area exhibiting a significant negative correlation (p < 0.05).The negative correlations between NDVI and T min exhibit a more complex pattern and occur in the montane steppe zone (i.e., IID2) and wet forest (i.e., OA1) regions.

Partial Correlations between NDVI and Asymmetric Diurnal Warming in Different Physico-Geographical Regions
Dynamic changes in vegetation were closely related to changes in T max and T min simultaneously.The unique topography and location of the Tibetan Plateau forms a unique alpine plateau climate.The synergistic effect of T max or T min on the alpine vegetation can be categorized into four types (see Figure 4).The first type includes the partial correlation coefficient between NDVI of the alpine vegetation and T max or T min , which were both positive with 25.7% of the plateau area were in this category.The percentage of the area for each physico-geographical zone meeting the criteria was the highest in IC1 (48.3%), followed by IB1 (35.3%),IC2 (25.1%),IIC2 (22.7%),IIAB1 (17.2%),IIC1 (12.7%), and OA1 (2.8%) in descending order.These were high-cold steppe or meadow zones in which vegetation growth is mainly limited by heat.Therefore, increases in T max and T min both accelerate the growth of vegetation.The second type includes the correlations between the NDVI of the alpine vegetation and T max or T min were both negative.This category occupies only 4.7% of the plateau, which is ranked at the lowest of the four categories by area.The physico-geographical zone that was most comprised of this category of area was IIC1 (9.6%), followed by IC2 (6.7%), IIAB1 (6.6%), OA1 (6.3%), IB1 (2.5%), IIC2 (1.3%), and IC1 (0.6%) in descending order.In other words, there were few areas where the night and day warming have each had effects on the vegetation of the Tibetan Plateau.The third type includes the correlation coefficient between the NDVI of the alpine vegetation and the T max or T min values, which were positive and negative, respectively.A total of 48.7% of the plateau area was in this category, which was the largest among the four categories.The physico-geographical zone that was most comprised of this category of area corresponding to this category was OA1 (90.3%), followed by IIC2 (69.6%),IIAB1 (62.2%),IB1 (50.2%),IC2 (40.3%),IC1 (30.9%), and IIC1 (29.3%) in descending order.In general, daytime warming has had a positive impact on most of the vegetation and nighttime warming enhanced carbon consumption which has had a negative impact on the vegetation NDVI on the Tibetan Plateau.In the last category, the correlation coefficient between the NDVI of the alpine vegetation and the T max and T min values were negative and positive, respectively, with 20.9% of the plateau area in this category.The physico-geographical zone comprised the largest part of this category was IIC1 (48.4%), followed by IC2 (27.9%),IC1 (20.1%),IIAB1 (14.0%),IB1 (12.0%),IIC2 (6.5%), and OA1 (0.6%) in descending order.In these areas, vegetation grew more during the day when the maximum temperature dropped and the nighttime minimum temperature increased.As ID1, IID1, IID2, and IID3 were desert zones, the NDVI in most of these regions was lower than 0.1.Therefore, these four regions were not analyzed in this particular analysis.

Partial Correlations between NDVI and Asymmetric Diurnal Warming among Different Vegetation Types
In order to assess correlations between NDVI and diurnal warming during the growing season for various vegetation types, we analyzed the correlation between NDVI and T max (or T min ) in the growing season after eliminating the effects of T min (or T max ), precipitation, and solar radiation (see Table 1).Partial correlations between growing season NDVI and T max were positively correlated with forest and crop vegetation areas exhibiting high significance (p < 0.01) while shrubland and others type were also significantly correlated (p < 0.05) with T max and alpine meadow appearing marginally significant (p < 0.1).However, partial correlations between the growing season NDVI and T min exhibited obvious divergence among the different vegetation types with all but alpine steppes and other types of vegetation that exhibit negative correlations.Forest (p < 0.01), crop (p < 0.05), and shrublands (p < 0.1) showed statistical significance in these correlations.Generally speaking, daytime warming enhanced vegetation growth while night-time warming had a positive impact on alpine steppe vegetation and the other types.In contrast, night-time warming had an adverse effect on forest, shrublands, alpine meadow, and crop vegetation types.To assess the spatial patterns of the responses of different vegetation types to daytime and night-time warming, we compared percentages of statistically significant (i.e., p < 0.05) pixels among each of the different vegetation types (see Table 2).The partial correlation coefficient between NDVI and T max relative to that between NDVI and T min for all vegetation types was more significant including the higher proportion of crop (41.46%), forest (34.64%), shrublands (23.04%), and alpine meadow (18.68%).In contrast, the percentage of pixels with significant partial correlations between NDVI and T min were generally relatively low for all vegetation types except for alpine steppe (19.86%) and others type (20.49%).The positive partial correlation between NDVI and T max exhibited overwhelming majority, but for partial correlation between NDVI and T min , forest (25.73%), shrublands (10.94%), and crop (39.02%) mainly showed negative correlations.Alpine steppe (14.01%), alpine meadow (5.92%), and others type (12.25%) mainly exhibited positive correlations.Taken together, the response of vegetation to daytime and night-time warming varied across vegetation types (Tables 1 and 2).We found that daytime warming has significantly positive influence on all types of vegetation while night-time warming has significantly negative influences on forest, shrublands, and crop vegetation types and has statistically significant positive influence on alpine steppe, alpine meadow, and other types of vegetation.

Discussion
We found that annual T max and T min of the Tibetan Plateau increased significantly from 1982 to 2015 and the spatial patterns of warming were heterogeneous.The increasing rate of seasonal T min over the past three decades was 1.7 times that of seasonal T max .This result suggested that the asymmetric diurnal warming exited in Tibetan Plateau, which was consistent with previous studies.For example, Liu and Chen (2000) found that the majority of the Tibetan Plateau has experienced statistically significant warming since the mid-1950s [31].Jun et al. (2013) suggested that there has been significant warming in T max and T min over the past three decades, especially for T min [32].Ma and Li (2003) demonstrated that the annual increases in T min over the past three decades have been one to three times that of T max across weather stations [33].
Vegetation NDVI of the Tibetan Plateau was predominantly positively correlated with changes in T max and 18.5% of the region showed significantly positive correlations (p < 0.05).These statistically significant positive correlations between NDVI and T max were found mainly in the colder or wetter regions which were IB1, IC1, IIC2, the west parts of IIAB1, and OA1.Peng et al. (2013) found that growing season NDVI was positively correlated with growing season T max in Tibetan Plateau [4].Liang et al. (2015) found that NDVI showed a significantly positive partial correlation with T max in the wetland vegetation of Nansi Lake [34].In the high-cold humid mountain regions, photosynthetic activity of the vegetation is limited by temperature more than water [13,35].Therefore, the photosynthetic enzyme activity is generally enhanced as T max rises [36].In contrast, in drier steppe regions (IIC1, IID1, IID2, IID3, and the western portion of IC2), interannual T max anomalies exhibit negative partial correlations with NDVI.Peng et al. (2013) also suggested that such negative correlations between NDVI and T max would occur in drier temperate regions [4].This negative correlation between T max and NDVI could be related to increased water stress via warming-induced soil moisture depletion [23,37].Moreover, decreased vegetation NDVI can reduce latent heat and increase sensible heat, which can result in additional warming [38].Previous field experiments also observed that the soil moisture in root zone significantly changed with warming temperature, which negatively affects vegetation photosynthetic activity [39,40].
However, in arid and semi-arid regions, we found that vegetation NDVI exhibited statistically significant positive correlations (p < 0.05) between NDVI and growing season T min while the T max , precipitation, and solar radiation were controlled for.Additionally, in the wettest forest and coldest montane steppe zone, vegetation NDVI presenting statistically significant negative correlations (p < 0.05) between NDVI and growing season T min after the T max , precipitation, and solar radiation were controlled for.These opposite responses of NDVI to T min between the wet and dry regions of the Tibetan Plateau could be attributed to night-time warming enhancing vegetation respiration, which leads to the consumption of more material and energy.As a result, it stimulates photosynthesis of vegetation the next day.Therefore, improving vegetation transpiration is accompanied by helping atmospheric CO 2 enter the leaves for photosynthesis, which enhances vegetation growth as indicated by NDVI.For instance, Wan et al. (2009) suggested that an increase in T min will enhance growth at a temperate, dry grassland site in China [41].However, Peng et al. (2004) demonstrated that an increase in T min will reduce rice yields by 10% per degree Celsius in the Philippines [12].Peng et al.
(2013) also observed these opposite responses of NDVI to T min between the wet and dry regions of the Northern Hemisphere [4].The reasons for this contrasting behavior may be related to the different mechanisms of temperature change across the plateau vegetation.Wang and Zhou (2004) studied the response characteristics of Leymus chinensis steppe in Inner Mongolia to temperature changes by using multi-year observational data from given plots and found that the response characteristics of different plant species to winter T min changes differed [42].Both this study and the previous research [4,12,41] showed that there were differences in the characteristics of diurnal warming which affected vegetation activities, and that the responses of different ecosystem types also differ.Additional detail research that integrates physiological, ecological, and hydrological observations would assist to better understand and explain the unique climate on the Tibetan Plateau.
In addition, we found that vegetation NDVI in different undisturbed natural areas of the Tibetan Plateau have had different responses to diurnal warming.The combination of altitudes, latitudes, and longitudes across the Tibetan Plateau forms hydrothermal patterns and alpine resources and environments, which lead to differences in the response of the same vegetation types to diurnal warming.The phenomenon of compensation is ubiquitous among life forms affected by stress and injuries and it is considered an adaptation of organisms to adverse environments [43].Therefore, the compensatory physiological responses of individual plants may also promote the different responses on broader scales and among the same vegetation to diurnal warming.
The temporal and spatial change in vegetation NDVI was the result of a combination of natural and anthropogenic factors.Certainly, it should be noted that grassland NDVI is also regulated by grassland degradation [44], thawing-freezing processes [45], and declining snow cover [46].In this study, the correlation between diurnal temperature and vegetation NDVI change was analyzed.Anthropogenic factors were not taken into account.

Conclusions
Our results suggested that the Tibetan Plateau has experienced a warming trend and the increasing rate of T min over the past three decades was much higher than T max .The warming trend and magnitude of T min and T max have not been uniform across the Tibetan Plateau.After removing the correlation between T max and T min , the partial correlation between T max and NDVI was significantly positive in wetter and colder regions of the Tibetan Plateau, but negative in the semi-arid southwest of the Tibetan Plateau.In contrast, the partial correlation between T min and NDVI is positive in the high-cold steppe and meadow steppe zones and exhibits a significant negative behavior in the montane steppe zone and wet forest zone.Meanwhile, asymmetric diurnal warming in the Tibetan Plateau has different effects on different types of vegetation.Daytime warming has significantly positive influences on all the vegetation, while night-time warming has mainly significantly negative influences on forest, shrublands, and crop vegetation types and has statistically significant positive influences on alpine steppe, alpine meadow, and others type.
The results of this study suggested that the variation in vegetation responses to diurnal temperature changes is important for understanding the changes in vegetation photosynthetic activity in a warming world.As future climate data shows the potential for asymmetric warming in global land surface climates, our results provide a reference for assessing and predicting vegetation responses to global climate change.Future research is needed to integrate T max and T min using remote sensing-based ecosystem models to reduce uncertainty regarding the terrestrial carbon cycle and corresponding climate feedback effects.

Figure 1 .
Figure 1.Spatial distribution of the altitude, meteorological stations, vegetation types, and the physico-geographical zones across the Tibetan Plateau.

Figure 2 .
Figure 2. Air temperature trends in the Tibetan Plateau from 1982 to 2015.(a) The slope of T max ; (b) the slope of T min ; (c) the slope of DTR; and (d) the significance level of DTR trends.

Figure 3 .
Figure 3. Spatial patterns of the correlations between growing season (April-October) NDVI and corresponding daily maximum (T max ) or minimum temperatures (T min ) in the Tibetan Plateau during 1982-2015.(a) The spatial patterns of partial correlation coefficients between NDVI and T max .The corresponding T min , precipitation, and solar radiation are controlled for in the calculation.(b) The partial correlation coefficients between NDVI and T min during the growing season with T max , precipitation, and solar radiation controlled for.The color-coded values ± 0.449, ±0.349, ±0.295, and ±0.233 correspond to p-values of 0.01, 0.05, 0.10, and 0.20 according to two-tailed Student's t-tests.

Figure 4 .
Figure 4. Spatial patterns of the correlation between NDVI in growing season (April-October) and corresponding daily maximum (T max ) or minimum temperatures (T min ) of the Tibetan Plateau from 1982-2015.Rmax represents the spatial patterns of partial correlation coefficients between NDVI and T max .The corresponding T min , precipitation, and solar radiation were controlled for in the calculation.Rmin shows the partial correlation coefficients between NDVI and T min during the growing season with T max , precipitation, and solar radiation are controlled for.+ denotes positive partial correlations while − denotes negative partial correlations.

Table 1 .
Partial correlation coefficients between NDVI and T max and T min for different vegetation types.