# Temporal Trends and Spatial Patterns of Temperature and Its Extremes over the Beijing-Tianjin Sand Source Region (1960–2014), China

Authors: Wei Wei; Baitian Wang; Kebin Zhang; Zhongjie Shi; Genbatu Ge; Xiaohui Yang
Journal: Advances in Meteorology
 (2018)

---

## Abstract

In order to examine temperature changes and extremes in the Beijing-Tianjin Sand Source Region (BTSSR), ten extreme temperature indices were selected, categorized, and calculated spanning the period 1960–2014, and the spatiotemporal variability and trends of temperature and extremes on multitimescales in the BTSSR were investigated using the Mann-Kendall (M-K) test, Sen’s slope estimator, and linear regression. Results show that mean temperatures have increased and extreme temperature events have become more frequent. Annual temperature has recorded a significant increasing trend over the BTSSR, in which 51 stations exhibited significant increasing trends (p<0.05); winter temperature recorded the most significant increasing trend in the northwest subregion. All extreme temperature indices showed warming trends at most stations; a higher warming slope in extreme temperature mainly occurred along the northeast border and northwest border and in the central-southern mountain area. As extreme low temperature events decrease, vegetation damage due to freezing temperatures will reduce and low cold-tolerant plants may expand their distribution range northward to revegetate barren areas in the BTSSR. However, in water-limited areas of the BTSSR, increasing temperatures in the growing season may exacerbate stress associated with plants relying on precipitation due to higher temperatures combining with decreasing precipitation.


---

## Body

## 1. Introduction

As highlighted by the IPCC, climate change is the most important environmental problem and the greatest challenge facing mankind [1]. From 1951 to 2012, global average surface temperatures are reported to have risen by 0.72 degrees [1]; as temperatures have increased, extreme climatic events, which can be sudden and destructive, have also increased in frequently, and their occurrence directly or indirectly has resulted in prominent social problems and natural disasters [2], especially in ecologically fragile areas. Recent global temperature extremes [3] indicate that the number of cold days, cold nights, and frost days has decreased while the number of warm days and warm nights has increased [3], these changes having obvious regional characteristics [4–7].Global regional-scale temperature changes are currently well documented. For example, Rehman and Al-Hadhrami [8] showed that the annual average maximum temperature on the west coast of Saudi Arabia has increased more than the average minimum temperature, while Khomsi et al. [9] recorded that hot and very hot events during the summer in Morocco (Safi and Marrakech) showed statistically significant decreasing trends. The number of very cold nights recorded annually and during the winter has significantly decreased across the Basilicata region (southern Italy) over the period 1951–2010, while the frequency of cold days has shown a weak increase [10]. Extreme temperature indices have indicated steady warming trends for arid and semiarid areas of southeastern Kenya [11], while days with 1-day extreme minimum or maximum apparent temperature have exhibited an increasing trend across many stations in the United States [12]. Temporal and spatial distribution of extreme temperatures across China shows that hot extremes indicate an upward trend, while cold extremes show downward trends; between 1961 and 2011, more frequent extreme temperature events were recorded in the southern area of the major grain producing region in China [13] and some trends of extreme temperature indices in the Songhua River Basin were stronger than those recorded in the Yangtze River Basin, southwestern China, and the Tibetan Plateau. Zhong et al. [14] also found the number of warm nights, warm days, and summer days would be more significant at higher latitudes. Sun et al. [15], investigating extreme temperature events in the Loess Plateau, found that the trend magnitudes of cold extremes were greater than those of warm extremes and that the growing season increased and diurnal temperature ranges declined in the region.Surface temperature changes have been mainly considered to be due to anthropogenic activities, especially due to changes in concentrations of greenhouse gases (GHG) [16, 17]. In areas surrounding Beijing and Tianjin, a large amount of land was reclaimed in the 1960s which destroyed the surface vegetation and caused mobile and semimobile sand dunes to expand [18]. Since the implementation of the reform and development policy (1979), environmental degradation has continued to intensify in these areas due to over grazing and excessive reclamation of grasslands [19, 20]; as a result, the Beijing, Tianjin, and Hebei areas experienced huge dust storms during March to April in the year 2000. To improve the environment of these areas, the Chinese government established the sand source control project (phaseI) in the Beijing-Tianjin Sand Source Region (BTSSR) in the year 2000 [21, 22]. Following the implementation of this project, more than half of the land in this region has experienced an increase in vegetation productivity during 2000–2010 [23]. Studies on the characteristics of climate change and the response of vegetation cover to weather changes in this region have found that NDVI has recorded a slight increasing trend in the growing season [24, 25]. In order to consolidate the results of phase I of the Blown-Sand Control Project, and to further reduce sandstorm hazards and construct the northern ecological barrier, phaseII has been implemented, spanning the period 2013–2022. However, with continued increasing temperatures, decreasing precipitation, and enlarging of the drought area, it is very difficult to consolidate the achievements of phase I and to implement and improve the construction of phaseII at the same time [23, 26]. Evaluation of extreme climate changes in the BTSSR can identify factors causing regional environmental and social problems and provide reference for the construction of phaseIIto cope with extreme climatic events and to improve disaster prevention and mitigation work in this area.The objectives of this study, therefore, are (1) to undertake trend analysis of temperature and extremes in the BTSSR; (2) to investigate the spatial patterns of trends of temperature and extremes; (3) and to discuss the possible impact temperature and extremes have on vegetation.


## 2. Study Area

The study area, located at the scope of the phaseII project in the BTSSR, includes six provinces with 138 counties. This area includes the Mu Us Sandland, situated between Yulin city (Shaanxi province) and Erdos city (Inner Mongolia), and the Kubuqi Desert, situated in the northern Ordos Plateau. The region covers an area of 71.05 × 104 km2 (approximately 36° ~46°N, 107° ~119°E; Figure 1). The area of land desertification covers an area of 22.69 × 104 km2, exceeding one-third of the total area of phaseII [22].Figure 1
Location of the study area and meteorological stations. Greek capital numbers indicate the subregion division: I is the southwest subregion (SR-I), II is the northwest subregion (SR-II), III is the south-central subregion (SR-III), IV is the east-central subregion (SR-IV), and V is the northeast subregion (SR-V). Source: [27].The climate in the study region is complex, including temperate semihumid, semiarid, arid, and extremely arid climate types [27]. Annual precipitation across the area varies from 105 to 743 mm, with the majority of precipitation occurring from June to September. Average temperature is 6.2°C with significant regional differences, such as a maximum temperature of 12.7°C in Beijing and a minimum temperature of −2.6°C in Arxan, Inner Mongolia. The average sunshine duration is 7.9 h/d, and the regional mean wind speed is 4.5 m/s. The average growing season length is 206 days. In general, the climate of the BTSSR is dry, windy, and has a low temperature, with climatic characteristics varying significantly across the region.


## 3. Data and Methods

### 3.1. Data Source, Quality Control, and Homogeneity Test

Daily temperature data from the BTSSR was provided by the China Meteorological Data Network (http://data.cma.cn/). When more than 10% of daily temperature data was missing, or more than three months contained more than 20% missing days, annual temperature records were deemed as missing. To ensure data integrity and consistency, meteorological stations lacking more than 75% of a year’s temperature records were omitted; missing days were not superseded by estimated daily temperature values [28, 29]. A total of 53 stations evenly distributed across the BTSSR were selected for this study; data from these stations spanned the period from 1960 to 2014 (Figure 1).Quality control of the data was undertaken using the RClimDex package developed by Zhang and Yang [30], a control package that has been previously applied to test climate data [31]. The RClimDex package enabled errors in the data to be eliminated, including errors in manual keying, daily maximum temperatures lower than daily minimum temperatures, and any identified outliers. Daily temperature outliers were identified by manually examining visual data graphs and histograms; suspicious outliers were identified using statistical tests, local knowledge, and comparisons with adjacent days or the same day at neighboring stations. Data that was clearly flawed were adjusted or removed.Due to relocation of meteorological stations, instrument changes, or observing procedures, step changes were identified at specific stations which resulted in observed meteorological data to not be fully comparable. To eliminate problems associated with step changes, daily temperature values were used to calculate the monthly temperatures; log transformations were subsequently used to test data homogeneity using the RHtest V4 [32]. Any step changes identified using RHtest V4 were analyzed using metadata to evaluate the changes. Results from these tests showed that all data from the 53 meteorological stations were acceptable for use in this investigation.After quality control and homogeneity tests, monthly temperature values were used to calculate mean annual temperatures, the standard climatological seasons (spring: March–May; summer: June–August; autumn: September–November; Winter: December–February) and the growing season (GS: April–October).


### 3.2. Calculation of Mean Temperature and Indices

To identify spatial changes in average temperature and extremes, the study region was divided into five subregions: the southeast subregion (SR-I), the northwest subregion (SR-II), the south-central subregion (SR-III), the east-central subregion (SR-IV), and the northeast subregion (SR-V), as per the investigation of Lu and Wu [27] (Figure 1). Temporally, two time cut points were determined as the starting year of reform and opening policy (ROP) (1979) and the starting year of the Blown-Sand Control Project Phase I in the BTSSR (hereafter referred to as project phase I) (2000). To identify the impact of land use/land cover on temperature and extremes, the study period (S: 1960–2014) was divided into three time stages: 1960–1978 (S1) (prior to ROP); 1979–1999 (S2) (ROP to before project phase I); and 2000–2014 (S3) (Project phaseI).Ten core indices for extreme temperature from ETCCDMI, indices which have been widely applied to evaluate extreme temperature shifts [33–40], were selected and classified into three categories: cold-related indices, warm-related indices, and variability indices (Table 1). All calculations of these indices were undertaken using RClimDex 1.0 software [30].Table 1
Definitions of temperature and extreme indices used in this study.
Type Index Descriptive name Definition Units Mean temperature TA Annual temperature Annual mean temperature °C T s p r i n g Spring temperature Mean temperature in spring °C T s u m m e r Summer temperature Mean temperature in summer °C T a u t u m n Autumn temperature Mean temperature in autumn °C T w i n t e r Winter temperature Mean temperature in winter °C TG Growing season temperature Mean temperature in growing season °C  Cold-related indices TN10p Cold nights Percentage of days when TN < 10th percentile Days TX10p Cold days Percentage of days when TX < 10th percentile Days CSDI Cold spell duration indicator Annual count of days with at least 6 consecutive days when TN < 10th percentile Days TNn Min Tmin Monthly minimum value of daily minimum temperature °C  Warm-related indices TN90p Warm nights Percentage of days when TN > 90th percentile Days TX90p Warm days Percentage of days when TX > 90th percentile Days WSDI Warm spell duration indicator Annual count of days with at least 6 consecutive days when TX > 90th percentile Days TXx Max Tmax Monthly maximum value of daily maximum temperature °C  Variability indices DTR Diurnal temperature range Monthly mean difference between TX and TN °C GSL Growing season length Annual (1st Jan to 31st Dec in NH, 1st July to 30th June in SH) count between first span of at least 6 days with TG > 5°C and first span after July 1 (January 1 in SH) of 6 days with TG < 5°C Days


### 3.3. Trend Analysis of Mean Temperature and Indices

Trends in mean temperature and indices were analyzed for the three time stages and for the duration of the study period. Trends for individual station data and regional series were estimated; the significance of the trends was determined using the Mann-Kendall test [41, 42]. This test does not assume that data are normally distributed, and it robustly responds to the effects of outliers in the series; this test has been previously used to test trends in hydrological and meteorological data [26, 43, 44].The Mann-Kendall test, enabling data with missing values and values below a detection limit to be analyzed, is one of the most important and commonly used statistical methods to detect trends as nonparametric tests in hydroclimatic time series data. The null hypothesis thatZc is not statistically significant or has no significant trend will be accepted if -Z1-α/2≤Zc≤Z1-α/2, where ±Z1-α/2 are the standard normal deviates and α is the significance level for the test [45].Kendall’s statisticS assessing the monotonic trend was calculated as(1)S=∑i=1n-1∑k=i+1nsgnxk-xi,varS=nn-12n+5-∑i=1meiei-12ei+518,where xk and xi are sequential data values representing the annual values in years k and i, respectively; n is the time length of the dataset; sgn(θ)=-1 if θ<0, sgn(θ)=0 if θ=0, and sgn(θ)=1 if θ>0; ei is the number of ties of the ith value and m is the number of tied values. The Mann-Kendall test statistic Zc was calculated as(2)Zc=S-1varS,ifS>0,Zc=0,ifS=0,Zc=S+1varS,ifS<0,where Zc is a standard normal variable. A positive Zc indicates an upward trend while a negative Zc indicates a downward trend [46]. A two-sided test at the significance level of 5% was undertaken in this study, and Sen’s slope estimator was used to estimate the true slope of an existing trend (the change per year) [47]. All of the trend calculations for these indices were undertaken using MAKESENS 1.0 [48]. Prior to performing the Mann-Kendall test, the trend-free, prewhitening method of Yue et al. [49] was used to limit the effect of serial correlations on the Mann-Kendall test [50, 51].


### 3.4. Spatial Interpolation of Mean Temperature and Extremes Trends

Mean temperature and indices trends were spatially interpolated using the spline method. This is characterized by fitting a smooth and continuous surface with the observation points, and it does not require a preliminary estimate for the structure of a temporal variance and statistical hypothesis [52]. In comparison with other methods, previous studies, for example, those using Kriging [53, 54] and those analyzing climatic and meteorological change [55–58], have shown the spline method to be more advantageous around topographic features. Spatial interpolation with the spline method was completed using ArcGIS software (Version 10.2).


## 3.1. Data Source, Quality Control, and Homogeneity Test

Daily temperature data from the BTSSR was provided by the China Meteorological Data Network (http://data.cma.cn/). When more than 10% of daily temperature data was missing, or more than three months contained more than 20% missing days, annual temperature records were deemed as missing. To ensure data integrity and consistency, meteorological stations lacking more than 75% of a year’s temperature records were omitted; missing days were not superseded by estimated daily temperature values [28, 29]. A total of 53 stations evenly distributed across the BTSSR were selected for this study; data from these stations spanned the period from 1960 to 2014 (Figure 1).Quality control of the data was undertaken using the RClimDex package developed by Zhang and Yang [30], a control package that has been previously applied to test climate data [31]. The RClimDex package enabled errors in the data to be eliminated, including errors in manual keying, daily maximum temperatures lower than daily minimum temperatures, and any identified outliers. Daily temperature outliers were identified by manually examining visual data graphs and histograms; suspicious outliers were identified using statistical tests, local knowledge, and comparisons with adjacent days or the same day at neighboring stations. Data that was clearly flawed were adjusted or removed.Due to relocation of meteorological stations, instrument changes, or observing procedures, step changes were identified at specific stations which resulted in observed meteorological data to not be fully comparable. To eliminate problems associated with step changes, daily temperature values were used to calculate the monthly temperatures; log transformations were subsequently used to test data homogeneity using the RHtest V4 [32]. Any step changes identified using RHtest V4 were analyzed using metadata to evaluate the changes. Results from these tests showed that all data from the 53 meteorological stations were acceptable for use in this investigation.After quality control and homogeneity tests, monthly temperature values were used to calculate mean annual temperatures, the standard climatological seasons (spring: March–May; summer: June–August; autumn: September–November; Winter: December–February) and the growing season (GS: April–October).


## 3.2. Calculation of Mean Temperature and Indices

To identify spatial changes in average temperature and extremes, the study region was divided into five subregions: the southeast subregion (SR-I), the northwest subregion (SR-II), the south-central subregion (SR-III), the east-central subregion (SR-IV), and the northeast subregion (SR-V), as per the investigation of Lu and Wu [27] (Figure 1). Temporally, two time cut points were determined as the starting year of reform and opening policy (ROP) (1979) and the starting year of the Blown-Sand Control Project Phase I in the BTSSR (hereafter referred to as project phase I) (2000). To identify the impact of land use/land cover on temperature and extremes, the study period (S: 1960–2014) was divided into three time stages: 1960–1978 (S1) (prior to ROP); 1979–1999 (S2) (ROP to before project phase I); and 2000–2014 (S3) (Project phaseI).Ten core indices for extreme temperature from ETCCDMI, indices which have been widely applied to evaluate extreme temperature shifts [33–40], were selected and classified into three categories: cold-related indices, warm-related indices, and variability indices (Table 1). All calculations of these indices were undertaken using RClimDex 1.0 software [30].Table 1
Definitions of temperature and extreme indices used in this study.
Type Index Descriptive name Definition Units Mean temperature TA Annual temperature Annual mean temperature °C T s p r i n g Spring temperature Mean temperature in spring °C T s u m m e r Summer temperature Mean temperature in summer °C T a u t u m n Autumn temperature Mean temperature in autumn °C T w i n t e r Winter temperature Mean temperature in winter °C TG Growing season temperature Mean temperature in growing season °C  Cold-related indices TN10p Cold nights Percentage of days when TN < 10th percentile Days TX10p Cold days Percentage of days when TX < 10th percentile Days CSDI Cold spell duration indicator Annual count of days with at least 6 consecutive days when TN < 10th percentile Days TNn Min Tmin Monthly minimum value of daily minimum temperature °C  Warm-related indices TN90p Warm nights Percentage of days when TN > 90th percentile Days TX90p Warm days Percentage of days when TX > 90th percentile Days WSDI Warm spell duration indicator Annual count of days with at least 6 consecutive days when TX > 90th percentile Days TXx Max Tmax Monthly maximum value of daily maximum temperature °C  Variability indices DTR Diurnal temperature range Monthly mean difference between TX and TN °C GSL Growing season length Annual (1st Jan to 31st Dec in NH, 1st July to 30th June in SH) count between first span of at least 6 days with TG > 5°C and first span after July 1 (January 1 in SH) of 6 days with TG < 5°C Days


## 3.3. Trend Analysis of Mean Temperature and Indices

Trends in mean temperature and indices were analyzed for the three time stages and for the duration of the study period. Trends for individual station data and regional series were estimated; the significance of the trends was determined using the Mann-Kendall test [41, 42]. This test does not assume that data are normally distributed, and it robustly responds to the effects of outliers in the series; this test has been previously used to test trends in hydrological and meteorological data [26, 43, 44].The Mann-Kendall test, enabling data with missing values and values below a detection limit to be analyzed, is one of the most important and commonly used statistical methods to detect trends as nonparametric tests in hydroclimatic time series data. The null hypothesis thatZc is not statistically significant or has no significant trend will be accepted if -Z1-α/2≤Zc≤Z1-α/2, where ±Z1-α/2 are the standard normal deviates and α is the significance level for the test [45].Kendall’s statisticS assessing the monotonic trend was calculated as(1)S=∑i=1n-1∑k=i+1nsgnxk-xi,varS=nn-12n+5-∑i=1meiei-12ei+518,where xk and xi are sequential data values representing the annual values in years k and i, respectively; n is the time length of the dataset; sgn(θ)=-1 if θ<0, sgn(θ)=0 if θ=0, and sgn(θ)=1 if θ>0; ei is the number of ties of the ith value and m is the number of tied values. The Mann-Kendall test statistic Zc was calculated as(2)Zc=S-1varS,ifS>0,Zc=0,ifS=0,Zc=S+1varS,ifS<0,where Zc is a standard normal variable. A positive Zc indicates an upward trend while a negative Zc indicates a downward trend [46]. A two-sided test at the significance level of 5% was undertaken in this study, and Sen’s slope estimator was used to estimate the true slope of an existing trend (the change per year) [47]. All of the trend calculations for these indices were undertaken using MAKESENS 1.0 [48]. Prior to performing the Mann-Kendall test, the trend-free, prewhitening method of Yue et al. [49] was used to limit the effect of serial correlations on the Mann-Kendall test [50, 51].


## 3.4. Spatial Interpolation of Mean Temperature and Extremes Trends

Mean temperature and indices trends were spatially interpolated using the spline method. This is characterized by fitting a smooth and continuous surface with the observation points, and it does not require a preliminary estimate for the structure of a temporal variance and statistical hypothesis [52]. In comparison with other methods, previous studies, for example, those using Kriging [53, 54] and those analyzing climatic and meteorological change [55–58], have shown the spline method to be more advantageous around topographic features. Spatial interpolation with the spline method was completed using ArcGIS software (Version 10.2).


## 4. Results and Discussion

### 4.1. Mean Temperature Spatiotemporal Trends

#### 4.1.1. Annual Temperature

MK test results across the study period for the BTSSR show that annual temperature (TA) significantly increased (p<0.001). The TA slope was 0.34°C/decade (Figures 2(a) and 3(a), Table 2), a result which coincides with TA change trends identified in the China territory [59, 60], in regions such as the Yellow River basin, the Yangtze River Basin [61], arid and semiarid regions in northwest China [62], and those in the monsoon area of eastern China [63]. TA results for the different stages showed some differences to those recorded for the entire study period. TA recorded a slightly downward trend in S1, a significantly upward trend in S2 (with a slope of 0.77°C/decade), and an insignificant downward trend in S3 (Figure 3(a), Table 2), results which may be due to alternating surface albedo due to large areas of desertification in S2 and rapid revegetation activities in S3 [23, 37, 64–70].Table 2
Trend slopes for average temperature and extreme temperature indices in the BTSSR during the different stages.
Index S1 S2 S3 S Increasing (%) Decreasing (%) Units T A −0.01 0.77∗∗ −0.14 0.34∗∗∗ 52 96.23 % 1.89% 1 0.00 % 1.89 % °C/decade T s p r i n g −0.01 0.76∗ 0.00 0.38∗∗∗ 52 98.11 % 0.00% 1 0.00 % 1.89 % °C/decade T s u m m e r −0.24 0.35 −0.32 0.24∗∗∗ 52 84.91 % 13.21% 1 0.00 % 1.89 % °C/decade T a u t u m n −0.03 0.64 0.29 0.30∗∗∗ 51 83.02 % 13.21% 2 3.77 % 0.00 % °C/decade T w i n t e r 0.64 1.28∗∗ −1.10 0.45∗∗∗ 53 79.25 % 20.75% 0 0.00 % 0.00 % °C/decade T G −0.02 0.47∗ −0.15 0.27∗∗∗ 52 94.34 % 3.77% 1 1.89 % 0.00 % °C/decade TN10p −0.98 −3.40∗∗ −0.35 −2.10∗∗∗ 3 0.00 % 5.66% 50 81.13 % 13.21 % days/decade TX10p 1.38 −2.19∗ −0.73 −0.96∗ 0 0.00 % 0.00% 53 56.60 % 43.40 % days/decade CSDI −0.22 −0.73∗ −0.81∗ −0.80∗∗ 4 0.00 % 7.55% 49 49.06 % 43.40 % days/decade TNn −0.17 1.61∗ 0.50 0.47∗∗ 48 52.83 % 37.74% 5 1.89 % 7.55 % °C/decade TN90p −0.17 3.49∗ −0.12 1.74∗∗∗ 51 88.68 % 7.55% 2 3.77 % 0.00 % days/decade TX90p −1.06 2.77 −1.37 1.13∗∗ 52 54.72 % 43.40% 1 0.00 % 0.00 % days/decade WSDI −0.27 3.52∗∗ −1.02 1.11∗∗ 53 79.25 % 20.75% 0 0.00 % 0.00 % days/decade TXx −0.23 0.36 −0.95 0.16 48 18.87 % 71.70% 5 0.00 % 9.43 % °C/decade DTR −0.43∗ −0.11 −0.12 −0.17∗∗∗ 10 13.21 % 5.66% 43 66.04 % 15.09 % °C/decade GSL 1.33 4.41 0.62 3.20∗∗∗ 52 79.25 % 18.87% 1 0.00 % 1.89 % days/decade S1 is 1960–1978, S2 is 1979–2000, and S3 is 2000–2014; ∗∗∗ is 0.001 level of significance, ∗∗ is 0.01 level of significance, and ∗ is 0.05 level of significance; percentage of stations with a significant (p<0.05) trend (bold), and percentage of stations with an insignificant (p>0.05) trend (italics).Figure 2
Anomaly of regional average temperature curves over the BTSSR, 1960–2014. Theblackline represents the linear trends for (a) annual, (b) spring, (c) summer, (d) autumn, (e) winter, and (f) growing season. Theredline represents the 5-year moving average.
(a) (b) (c) (d) (e) (f)Figure 3
Box-Plot of Slopes for regional average temperature over the BTSSR in different stages. (a) Annual, (b) spring, (c) summer, (d) autumn, (e) winter, and (f) growing season. S1 is 1960–1978, S2 is 1979–2000, S3 is 2000–2014, and S is 1960–2014. ∗∗∗ is 0.001 level of significance, ∗∗ is 0.01 level of significance, and ∗ is 0.05 level of significance.
(a) (b) (c) (d) (e) (f)The variations inTA trends for the different stages in the five subregions are shown in Figure 4(a) and Table 3. Results show that TA trends in these subregions significantly increased from 1960 to 2014 (p<0.001); the highest TA slope was 0.41°C/decade in the northwest subregion (SR-II) while the east-central subregion (SR-IV) recorded the lowest slope (0.30°C/decade) (Figures 3(a) and 4(a), Table 3). Results for all five subregions recorded a slight fluctuation during S1 and a significant increasing trend in S2, especially in the northwest subregion (SR-II; 0.92°C/decade), a result which may be attributed to serious steppe degradation in this area [71, 72]. During the last decade, TA recorded slightly negative trends, except for the northwest subregion (SR-II).Table 3
Trend slopes (°C/decade) of mean temperature and extreme temperature indices in the five subregions during the different stages.
Stage S1 S2 S3 S Index I II III IV V I II III IV V I II III IV V I II III IV V Units  T A 0.00 0.08 −0.02 −0.02 −0.08 0.74∗ 0.92∗∗ 0.77∗∗ 0.68∗∗ 0.85∗∗ −0.09 0.08 −0.22 −0.18 −0.29 0.31∗∗∗ 0.41∗∗∗ 0.36∗∗∗ 0.30∗∗∗ 0.37∗∗∗ °C/decade T s p r i n g −0.01 −0.18 0.05 0.03 −0.11 0.85∗ 0.82 0.73∗ 0.56 0.91 0.36 0.51 −0.04 −0.58 −0.27 0.36∗∗∗ 0.38∗∗∗ 0.40∗∗∗ 0.34∗∗ 0.41∗∗∗ °C/decade T s u m m e r −0.09 −0.35 −0.42 −0.20 0.04 0.61∗ 0.49 0.38 0.08 −0.01 −0.40 −0.12 −0.36 −0.09 −0.75 0.19∗∗ 0.31∗∗∗ 0.25∗∗∗ 0.22∗∗∗ 0.31∗∗∗ °C/decade T a u t u m n −0.10 0.03 −0.16 0.15 0.13 0.64 0.71 0.62 0.56 0.75 0.43 0.33 0.28 0.19 0.08 0.24∗∗∗ 0.40∗∗∗ 0.29∗∗∗ 0.29∗∗∗ 0.38∗∗∗ °C/decade T w i n t e r 0.54 1.17 0.77 0.42 0.25 1.07∗ 1.37∗ 1.34∗∗ 1.28∗∗ 1.50∗ −1.13 −0.94 −1.09 −1.12 −1.21 0.45∗∗∗ 0.46∗∗ 0.50∗∗∗ 0.41∗∗ 0.37∗ °C/decade T G 0.11 −0.12 −0.07 −0.11 0.08 0.53∗∗ 0.58∗∗ 0.45 0.34 0.49∗ −0.04 0.02 −0.15 −0.28 −0.46 0.23∗∗∗ 0.34∗∗∗ 0.27∗∗∗ 0.26∗∗∗ 0.35∗∗∗ °C/decade TN10p −0.93 −1.62 −0.23 −1.55 −2.07 −3.32 −3.04∗ −4.72∗∗ −2.65∗∗ −2.64∗∗ −0.20 −0.64 0.39 −1.04 −0.53 −1.72∗∗∗ −2.30∗∗∗ −2.30∗∗∗ −2.14∗∗∗ −2.28∗∗ days/decade TX10p 1.10 1.13 1.74 1.29 1.52 −2.38 −2.59∗ −2.20∗ −1.80 −1.95∗ −1.33 −1.20 −0.71 −0.07 0.09 −0.97∗∗ −0.86∗ −1.13∗∗ −0.82∗ −0.81 days/decade CSDI 0.74 −2.64 0.29 −0.35 −0.92 −0.45 −0.66 −0.89∗ −0.83 −0.92∗ −0.61 −0.53 0.18 −2.54∗∗ −1.14∗∗ −0.72∗ −1.38∗∗ −0.54∗∗ −0.81∗∗ −1.04∗∗ days/decade TNn −0.07 0.98 0.09 0.05 0.19 1.03 1.26 1.78∗∗ 2.10∗∗ 2.05∗∗ 0.81 0.90 0.52 0.17 −0.20 0.25 0.55∗∗ 0.60∗∗∗ 0.49∗∗ 0.49∗ °C/decade TN90p −0.33 0.20 −0.58 0.11 0.48 3.96∗ 4.01∗ 3.34∗ 2.96 3.09∗ 1.17 1.96 −1.31 −1.18 −0.48 1.60∗∗∗ 2.07∗∗∗ 1.80∗∗∗ 1.63∗∗∗ 1.70∗∗∗ days/decade TX90p −0.55 −1.44∗ −1.44 −1.07 −0.71 3.29 3.69 2.69 1.98 2.03 −0.56 0.45 −3.16 −1.53 −2.13 1.31 1.07∗ 1.17∗∗∗ 0.87∗ 1.13∗∗ days/decade WSDI 0.58∗∗ −1.57∗∗ −0.32 −0.28 −0.46 5.42∗∗ 3.29∗∗ 4.31∗∗ 0.57∗ 2.45∗∗ 1.20 2.16 4.31 −2.76∗∗ −2.92 1.34∗∗ 1.00∗∗ 1.40∗∗ 0.57∗∗ 0.88∗∗ days/decade TXx 0.33 −0.03 −0.78 −0.42 0.11 0.54 1.02 0.33 −0.16 0.22 −1.22 −0.13 −1.08 −1.12 −0.54 0.11 0.23 0.23 0.05 0.24 °C/decade DTR −0.31 −0.60∗∗ −0.33 −0.54∗∗ −0.58∗∗ −0.04 −0.06 −0.18 −0.10 −0.21 −0.09 −0.01 −0.07 −0.29 −0.14 −0.09 −0.23∗∗∗ −0.17∗∗∗ −0.20∗∗∗ −0.21∗∗∗ °C/decade GSL −4.26 5.92 3.56 1.11 3.90 3.39 5.17 3.03 5.64 7.95 0.14 3.52 −0.84 0.35 3.23 2.92∗∗ 3.42∗∗∗ 3.40∗∗∗ 3.08∗∗∗ 3.28∗∗∗ days/decade S1 is 1960–1978, S2 is 1979–2000, and S3 is 2000–2014; ∗∗∗ is 0.001 level of significance, ∗∗ is 0.01 level of significance, and ∗ is 0.05 level of significance.Figure 4
Box-Plot of Slopes for regional average temperature over the five subregions in different stages. (a) Annual, (b) spring, (c) summer, (d) autumn, (e) winter, and (f) growing season. S1 is 1960–1978, S2 is 1979–2000, S3 is 2000–2014, and S is 1960–2014. I is the southwest subregion (SR-I), II is the northwest subregion (SR-II), III is the south-central subregion (SR-III), IV is the east-central subregion (SR-IV), and V is the northeast subregion (SR-V). ∗∗∗ is 0.001 level of significance, ∗∗ is 0.01 level of significance, and ∗ is 0.05 level of significance.
(a) (b) (c) (d) (e) (f)Results for the spatial distribution ofTA trends (Figure 5(a)) show that most areas of the BTSSR had positive slopes, with the largest occurring in the central-southern mountain area; however, negative slopes were identified in the southeastern region in the Luanhe River Basin. Based on data from 51 stations (96.23% of all stations), the spatial distribution of TA showed a significant increase in the BTSSR (p<0.01) (Table 2), corresponding to the average slope of 0.35°C/decade. A strong increasing trend was observed at the Wutaishan station and a weak negative trend was found at the Chengde station, both meteorological stations being situated in the south-central subregion (SR-III).Figure 5
Spatial patterns of the trend slopes of mean temperature (1960–2014) over the BTSSR and anomaly of subregional indices curves. Theblackline represents linear trends for (a) annual, (b) spring, (c) summer, (d) autumn, (e) winter, and (f) growing season. Theredline represents the 5-year moving average.
(a) (b)  (c) (d)  (e) (f)In general, an increasingTA trend in the BTSSR can accelerate plant growth by altering plant photosynthesis processes [73, 74] and lengthening the growing season [75, 76]. However, increasing TA trends can also increase regional evapotranspiration [26] which, in conjunction with a decreasing trend in precipitation in the BTSSR [77], will increase the risk of drought [78]. The tradeoff of change trends of these meteorological factors is currently unclear.


#### 4.1.2. Seasonal Temperature

MK test results indicate that the seasonal temperature of the BTSSR significantly increased (p<0.05) during the study period; the lowest (0.24°C/decade) and highest (0.45°C/decade) rate of temperature rise were recorded in summer and winter, respectively (Figures 2(b)–2(e) and 3(b)–3(e), Table 2). This result was similar to that recorded by Wang et al. [79] for seasonal temperature changes across the arid region of northwestern China. Results for the different stages of the BTSSR showed differing changes. Results during S1 recorded seasonal temperature to have a slightly negative trend, except for winter temperature (Twinter). Seasonal temperature in S2 increased, with Twinter recording a significant increase (having a slope of 1.28°C/decade) and Tsummer having an insignificant increase (with a slope of 0.35°C/decade), results which were in agreement with seasonal temperature changes across northeastern China [80]. Seasonal temperatures during S3 recorded a mixture of slightly negative and positive trends (Table 2).Results for seasonal temperature trends across the five subregions (Figures4(b)–4(e), Table 3) recorded a significant increase (p<0.05) during the study period, especially for Twinter in SR-III (having a slope of 0.50°C/decade) (Figures 5(b)–5(e)). During S1, insignificant increasing trends for Twinter were recorded in all subregions, while other seasonal temperature series showed a mixture of slightly negative and positive trends across the five subregions. During S2, Twinter showed a significant increasing trend in all subregions, especially in SR-V (with a slope of 1.50°C/decade), and other seasonal temperature series showed insignificant positive trends for most subregions. Insignificant decreasing trends were recorded in the five subregions during the winter months and summer months, while spring temperatures showed a mixture of slightly negative and positive trends, with more negative than positive trends being recorded, and slightly increasing trends were found in all five subregions during the autumn months over the last decade (Figures 4(b)–4(e), Table 3).The spatial distribution of seasonal temperature trends (Figures5(b)–5(e)) shows that a positive temperature slope was identified for almost the entire BTSSR in the four seasons; the southeastern area of SR-III and the eastern area of SR-I were exceptions, these areas recording negative Tspring, Tsummer, and Tautumn slopes and a negative Tautum slope, respectively (Figures 5(b)–5(e)). Most stations exhibited increasing trends, with the lowest winter slope of 0.01°C/decade (Hequ station) and the highest winter slope of 1.51°C/decade (Wutaishan station) being located in SR-I and SR-III, respectively. In spring, data from 98.11% of stations recorded a significant increasing trend (Figure 5(b), Table 2). In summer and autumn, the number of stations recording significant positive trends was similar, the trends mainly occurring in SR-I, SR-III, and SR-IV. Similarly, the number of stations recording negative trends was few (Figures 5(c) and 5(d), Table 2). In the winter, 42 stations recorded significant positive trends (79.25% of all stations), these being mainly located in SR-I, SR-III, and SR-IV, with no station recording a negative trend (Figure 5(e), Table 2). This result implies that the southern and east-central areas of the BTSSR experienced a warmer climate each year.The winter warming trends identified in our study across the BTSSR are similar to those previously identified in Europe [81], North Asia and North America [82], subarctic [83], China [84, 85], and other areas across the globe [86]. Winter warming will undoubtedly be beneficial in alleviating vegetation damage due to freezing temperatures, to both native and planted vegetation, and it will increase the possibility of using low cold-tolerant plants to revegetate the BTSSR.


#### 4.1.3. Growing Season Temperature

MK test results for growing season temperature (TG) from 1960 to 2014 across the BTSSR recorded a significant increase (p<0.001), with a slope of 0.27°C/decade (Figures 2(f) and 3(f), Table 2). Results during the different stages showed TG to record a slight decrease in S1, a significant increase in S2 (with a slope of 0.47°C/decade), and an insignificant decline in S3 (Figure 3(f), Table 2). The overall trend of TG was similar with that of TA.Results for the variations inTG trends for the five subregions in the different stages recorded a significant increase (p<0.001; Figure 4(f), Table 3) from 1960 to 2014, especially in SR-V, with a slope of 0.35°C/decade (Figures 5(f) and 4(f), Table 3). As temperatures increased, the date of initiation of the growing season advanced and earlier growth resulted in vegetation covering the land earlier, thus reducing the occurrence of dust storms [87]. However, earlier vegetation growth could increase water stress in the reproductive growth period of water-limited vegetation in the BTSSR [24]. TG trends for all five subregions recorded slight fluctuations during S1. During S2, a slope of 0.58°C/decade was recorded, indicating a significant increasing trend, especially in SR-II. This result was consistent with previous findings showing TG to have significantly increased across most regions of China since 1980 [88]. During the last decade (S3), TG showed slightly negative trends, except in SR-II (Figure 4(f), Table 3).Results for theTG slope showed a similar spatial pattern with those for TA (Figure 5(f)). A larger positive slope was recorded in the central-southern mountain area and the northeast region, and a negative slope occurred in the southeastern area of SR-III. Based on data from 50 stations (94.34% of all stations; Table 2), the spatial distribution of TG showed a significant increase over the BTSSR (p<0.05) (Figure 5(f)), corresponding to an average slope of 0.29°C/decade. A strong increasing trend was observed at the Wutaishan station and a significant negative trend was only found at the Chengde station, both stations being located in SR-III (Figure 5(f)).Our results showed that a trend of increasingTG can enhance plant photosynthetic activities and further productivity without causing water stress [89, 90]. However, in water-limited areas of the BTSSR, increasing temperature in the growing season may exacerbate rain-fed plant stress due to high temperatures combined with decreasing precipitation [77, 91], especially for vegetation planted (shrubs or trees) as part of the degraded land restoration project in the BTSSR.


### 4.2. Spatiotemporal Trends of Indices

#### 4.2.1. Cold-Related Indices

Results for temporal variation of cold nights (TN10p), cold days (TX10p), and cold spell duration indicator (CSDI), representing the frequency and duration of extreme low temperature for the BTSSR, are shown in Figures6(a)–6(c). In general, from 1960 to 2014, the fluctuations of these indices continuously significantly decreased (p<0.05), especially TN10p (having a slope of −2.10 days/decade; Figure 7(a), Table 2). Changes observed for TN10p are more significant than those of TX10p across the BTSSR (Figures 7(a) and 7(b), Table 2); a reduction in TN10p will result in a reduction in frost damage during the night [92, 93]. The results of our study were in accordance with those over the arid region of northwestern China [79] and in Iran [94]. Results for the different stages showed that, during S1, the frequency and duration of extreme low temperature recorded a slightly negative trend, except for TX10p. During S2, the frequency and duration of extreme low temperatures significantly decreased, especially TN10p (with a slope of −3.40 days/decade), and results during S3 showed CSDI to significantly decrease (a slope of −0.81 days/decade), while TN10p and TX10p insignificantly decreased (Figures 7(a)–7(c), Table 2). Figure 6(d) shows the temporal variation of Min Tmin (TNn) across the BTSSR. From 1960 to 2014, the fluctuations of TNn continuously significantly increased, having a slope of 0.47°C/decade (p<0.01). During the different stages, TNn showed a slightly positive trend during S1, a significant increase in S2 (with a slope of 1.61°C/decade), and an insignificant increased in S3 (Figure 7(d), Table 2).Figure 6
Anomaly of extreme temperature indices curves across the BTSSR (1960–2014). Theblack line represents the linear trends for (a) TN10p, (b) TX10p, (c) CSDI, (d) TNn, (e) TN90p, (f) TX90p, (g) WSDI, (h) TXx, (i) DTR, and (j) GSL. Theredline represents the 5-year moving average.
(a) (b) (c) (d) (e) (f) (g) (h) (i) (j)Figure 7
Box-Plot of Slopes for extreme temperature indices across the BTSSR in different stages. (a) TN10p, (b) TX10p, (c) CSDI, (d) TNn, (e) TN90p, (f) TX90p, (g) WSDI, (h) TXx, (i) DTR, and (j) GSL. S1 is 1960–1978, S2 is 1979–2000, S3 is 2000–2014, and S is 1960–2014. ∗∗∗ is 0.001 level of significance, ∗∗ is 0.01 level of significance, and ∗ is 0.05 level of significance.
(a) (b) (c) (d) (e) (f) (g) (h) (i) (j)Variations in the trends for the extreme cold indices at the different stages in the five subregions (Figures8(a)–8(d), Table 3) show that their frequency and duration for the majority of the subregions significantly decreased from 1960 to 2014 (p<0.05). TN10p results in SR-III recording a slope of −2.30 days/decade (Figures 8(a) and 9(a), Table 3) and, for the majority of subregions, TNn significantly increased during the study period (p<0.05), notably in SR-III (with a slope of 0.60°C/decade, Figures 8(d) and 9(d), Table 3). During the period S1, all subregions showed a mixture of negative and positive trends, except for TN10p and TX10p. These fluctuating trends stabilized during S2 as a significant warming trend occurred, especially for TN10p in SR-III (with a slope of −4.72 days/decade), and for TNn in SR-IV (with a slope of 2.10°C/decade), these results being in accordance with observations over China during 1979–1999 [95]. During S3, TN10p, TX10p, and CSDI showed negative trends for most subregions while TNn showed a positive trend for most subregions, especially in SR-II (with a slope of −0.90°C/decade).Figure 8
Box-Plot of Slopes for extreme temperature indices across the five subregions in different stages. (a) TN10p, (b) TX10p, (c) CSDI, (d) TNn, (e) TN90p, (f) TX90p, (g) WSDI, (h) TXx, (i) DTR, and (j) GSL. S1 is 1960–1978, S2 is 1979–2000, S3 is 2000–2014, S is 1960–2014. ∗∗∗ is 0.001 level of significance, ∗∗ is 0.01 level of significance, and ∗ is 0.05 level of significance.
(a) (b) (c) (d) (e) (f) (g) (h) (i) (j)Figure 9
Spatial patterns of the trend slopes of cold-related indices across the BTSSR and anomaly of subregional indices curves, 1960–2014. Theblack line represents the linear trends for (a) TN10p, (b) TX10p, (c) CSDI, and (d) TNn. Theredline represents the 5-year moving average.
(a) (b)  (c) (d)The slope distribution of TN10p showed a spatial disparity across the BTSSR (Figure9(a)). A positive slope was identified in the eastern area of SR-I, the northern area of SR-II, and the southeastern area of SR-III, this being the area around the Chengde station; a negative slope prevailed across the rest of the region where higher negative slopes occurred along the northeast border and northwest border and in the south-central part. Results for the spatial distribution of TN10p, based on data from 43 stations (81.13% of all stations), showed a significant decrease across the BTSSR (p<0.05) (Table 2), corresponding to an average slope of −5.05 days/decade. A strong decreasing trend was observed at the Wutaishan station in SR-III and a weak positive trend was identified at the Sonid Left station (SR-II; Figure 9(a)).A negative slope for TX10p was identified across the entire region, with the largest slope occurring in the central-southern mountain area and in the east Siramulen River Basin (Figure9(b)). Data from 56.60% of stations recorded a significant decreasing trend (p<0.05) with no stations recording an increasing trend (Table 2). The Wutaishan station in SR-III recorded the strongest decreasing trend in the BTSSR (Figure 9(b)).Results for CSDI also recorded a negative slope across the majority of the region, an exception being in the eastern parts of SR-I and SR-III where positive slopes occurred (Figure9(c)). Data from 49.06% of stations recorded a significant decreasing trend, concentrated in the northern part of the study area (Figure 9(c), Table 2); a strong decreasing trend was observed at the Linhe station (SR-I) and a weak positive trend was identified at the Yuanping station (SR-III; Figure 9(c)).The slopes for TNn and TN10p recorded similar patterns (Figure9(d)), with larger positive slopes being recorded along the northeast-northwest border and in the south-central part, and negative slopes in the eastern and southern areas of SR-I, the northern area of SR-II, and the eastern area of SR-III. Based on data from 28 stations (52.83% of all stations; Table 2), the spatial distribution of TNn showed a significant increase across the BTSSR (p<0.05), corresponding to an average slope of 0.75°C/decade. A strong increasing trend was observed at the Wutaishan station (SR-III) and a significant negative trend was identified at the Hequ station (SR-I; Figure 9(d)).Results for the study period across the BTSSR showed a general trend of warming for the cold-related indices, trends which were also identified in northern Mongolia [96], South Africa [97], northeast China [98], and globally [99]. As extreme low temperature events decreased in the BTSSR, vegetation damage due to freezing conditions was mitigated, thus resulting in an expansion of range northwards of low, cold-tolerant plants which can result in the enrichment of plant species in this area.


#### 4.2.2. Warm-Related Indices

The frequency and duration of extreme high temperature events across the BTSSR are shown by the results for the temporal variation of warm nights (TN90p), warm days (TX90p), and the warm spell duration indicator (WSDI) (Figures6(e)–6(g)). In general, these indices fluctuated continuously, having a significant increasing trend during the study period (p<0.01). The results for TN90p (with a slope of 1.74 days/decade; Figures 7(e)–7(g), Table 2) were noticeable, these being in accordance with previous findings across China [100]. Results for the different stages showed that the frequency and duration of extreme high temperature recorded an insignificant negative trend during S1. Results during S2 only recorded significant increases for TN90p and WSDI; TN90p had the most noticeable increase with a slope of 3.49 days/decade; during S3, the frequency and duration of extreme high temperature decreased (Figures 7(e)–7(g), Table 2). The temporal variation of Max Tmax (TXx) across the BTSSR (Figure 6(h)) from 1960 to 2014 recorded fluctuations to insignificantly increase, having a slope of 0.16°C/decade. TXx results during the different stages recorded a slightly negative trend in S1, an insignificant increase in S2, and an insignificant decrease in S3 (Figure 7(h), Table 2).Variations in extreme warm indices in the five subregions (Figures8(e)–8(h), Table 3) showed that frequency and duration of extreme high temperature recorded significant increasing trends (1960–2014) for the majority of the subregions, especially TN90p (Figures 10(a)–10(c)). The change of TN90p is more pronounced than that of TX90p, a finding which is in accordance with previous studies [79, 94]. The temporal variation of TXx for the five subregions (Figure 10(d)) shows an increasing trend for all subregions, especially SR-V (with a slope of 0.24°C/decade). As daily maximum temperatures increase, transpiration and plant water stress will be exacerbated which may lead to an increase in water demand for vegetation production [101]. All five subregions showed that there was a combination of negative and positive trends during S1, with more negative trends being identified. During S2, an increasing trend occurred for the majority of the subregions, with more insignificant than significant trends, especially for TN90p in SR-III (with a slope of 3.34 days/decade). These findings are also in accordance with changes of temperature extremes in China as a whole [95]. During the last decade (S3), the five subregions showed a mixture of negative and positive trends, with more negative trends recorded than positive (Figures 8(e)–8(h), Table 3).Figure 10
Spatial patterns of the trend slopes of warm-related indices across the BTSSR and anomaly of subregional indices curves, 1960–2014. Theblack line represents the linear trends for (a) TN90p, (b) TX90p, (c) WSDI, and (d) TXx. Theredline represents the 5-year moving average.
(a) (b)  (c) (d)Results across the BTSSR showed that positive slopes in TN90p occurred across the majority of the region, with the largest trend occurring in the central-southern mountain area. However, negative slopes were identified in the Luanhe River Basin of the southeastern region and the east area of the southwest subregion (SR-I) (Figure10(a)). Based on data from 47 stations (88.68% of all stations) (Table 2), the spatial distribution of TN90p showed a significant increase for the BTSSR (p<0.05), corresponding to an average slope of 4.58 days/decade. Significant decreasing trends (at the 99% confidence level) were identified at the Hequ and Chengde stations (Figure 10(a)), accounting for 3.77% of the total stations, and a strong increasing trend was observed at the Wutaishan station (SR-III).A positive slope for TX90p was identified across almost the entire BTSSR, the only exception being the southeastern area of SR-I where a negative slope was recorded. The larger positive slopes were recorded in the central-southern mountain area, and the number of stations recording positive trends was much larger than those recording negative trends (Figure10(b)). Twenty-nine stations recorded significant increases (54.72% of the total stations), with the Wutaishan station in SR-III recording the strongest increasing trend, and only one station recorded a decrease (Table 2; Figure 10(b)).WSDI results across the region recorded positive slopes, with the largest being in the central-southern mountain area (Figure10(c)). Positive trends were significant at 42 stations (79.25%), these being evenly distributed across the region; there were no negative trends recorded for this variable (Figure 10(c), Table 2). The strongest increasing trend was observed at the Wutaishan station in SR-III (Figure 10(c)).TXx slope distribution results across the BTSSR showed a spatial disparity. Negative slopes were identified in the eastern and southern areas of SR-I, in southeastern areas of SR-IV, and in the northern area of SR-V. Positive slopes were present in the other areas of the study region, with the highest positive slopes occurring in the central-southern mountain area (Figure10(d)). The number of stations recording positive trends was much greater than those recording negative trends; 10 stations recorded significant increases (18.87%), with no significant decreases being recorded at any station (Table 2). The strongest increasing trend was recorded at the Wutaishan station (SR-III; Figure 10(d)).Differences in results between the minimum nighttime temperature (TNn) and the maximum daytime temperature (TXx) (Figures7(d) and 7(h), Table 2) indicated that temperatures across the BTSSR showed a general warming trend, a result which is consistent with global increases [102]. The rate of warming of TXx was lower than that of TNn over all of the subregions (Figures 8(d) and 8(h), Table 3), an occurrence which is not conducive to daytime vegetation growth; this warming will have the opposite effect on nighttime vegetation growth by affecting carbon absorption and consumption [103, 104]. Moreover, the increase of maximum daily temperature may inhibit the growth of grass across the BTSSR; an increase of the minimum temperature may promote grassland growth via the compensatory effect of increasing nocturnal respiration [105], as well as reducing desert vegetation productivity by accelerating plant respiration [106], increasing the rate of nutrient metabolism, and shortening the duration of grain filling [107, 108].


#### 4.2.3. Variability Indices

The temporal variation of the diurnal temperature range (DTR) across the BTSSR showed DTR fluctuations to continuously significantly decrease (p<0.001), with a slope of −0.17°C/decade (Figures 6(i) and 7(i), Table 2). During the different stages, DTR showed a significant negative trend during S1 (with a slope of −0.43°C/decade; p<0.05), and insignificant decreases during S2 and S3 (Figure 7(i), Table 2). DTR trend variations for the different stages for the majority of the five subregions (Figure 8(i), Table 3) showed a significant decreasing trend across the study period, especially in SR-II which had a slope of −0.23°C/decade (Figures 8(i) and 11(a), Table 3). All five subregions showed negative trends during S1; SR-II recorded a more significant trend, having a slope of −0.60°C/decade. During S2, insignificant decreasing trends were recorded for all subregions, and during the last decade (S3) the five subregions showed negative trends (Figure 8(i), Table 3). Higher DTR is propitious to the accumulation of biomass [93, 109], therefore, in the BTSSR, as DTR decreased the quality of vegetation may decline.Figure 11
Spatial patterns of the trend slopes of variability indices over the BTSSR and anomaly of subregional indices curves, 1960–2014. Theblack line represents the linear trends for (a) DTR, (b) GSL. Theredline represents 5-year moving average.
(a) (b)The slope distribution of DTR shows a spatial disparity in the BTSSR (Figure11(a)). A positive slope was recorded in the north, east, and southwestern areas of the region, and a negative slope prevailed in the other areas where higher negative slopes occurred along the northeast-northwest border and the central-southern mountain area. Based on data from 35 stations (66.04% of all stations), the spatial distribution of DTR showed a significant decrease across the BTSSR (p<0.05) (Table 2), corresponding to an average slope of −0.29°C/decade. Seven stations (13.21% of total stations) had a significant increasing trend at the 95% confidence level. A strong decreasing trend was observed at the Weixian station (SR-III) and a strong increasing trend was observed at the Hequ station (SR-I; Figure 11(a)).The temporal variation of the growing season length (GSL) across the BTSSR during the study period continuously increased significantly by 3.20 days/decade from 1960 to 2014 (p<0.001) (Figures 6(j) and 7(j), Table 2). Results for the three stages showed that GSL across the BTSSR showed an insignificant positive trend during S1, S2, and S3 (Figure 7(j), Table 2). The variations in GSL trends for the different stages in the five subregions showed significant increasing trends for the majority of the subregions during the study period, especially for GSL results in SR-II (with a slope of 3.42 days/decade; Table 3, Figures 8(j) and 11(b)). As the length of the growing season increased, the potential number of harvests and seasonal yields for perennial forage crops may be promoted at high latitudes [110]. All five subregions showed a mixture of negative and positive trends, with more positive trends recorded during S1. During S2, insignificant increasing trends occurred for all subregions and, during the last decade (S3), the five subregions recorded a mixture of negative and positive trends, with more positive than negative trends recorded (Figure 8(j), Table 3).The slope distribution of GSL shows a spatial disparity across the BTSSR (Figure11(b)). Negative slopes were recorded in the eastern areas of SR-III and positive slopes were recorded for all of the other areas in the region; the highest positive slopes occurred in the central-southern mountain area. The number of stations recording positive trends was much greater than those recording negative trends; 42 stations recorded significant increases (79.25%) and no station recorded a significant decrease (Table 2). A strong increasing trend was observed at the Wutaishan station (SR-III) and a negative trend was only found at the Chengde station (SR-III; Figure 11(b)).


## 4.1. Mean Temperature Spatiotemporal Trends

### 4.1.1. Annual Temperature

MK test results across the study period for the BTSSR show that annual temperature (TA) significantly increased (p<0.001). The TA slope was 0.34°C/decade (Figures 2(a) and 3(a), Table 2), a result which coincides with TA change trends identified in the China territory [59, 60], in regions such as the Yellow River basin, the Yangtze River Basin [61], arid and semiarid regions in northwest China [62], and those in the monsoon area of eastern China [63]. TA results for the different stages showed some differences to those recorded for the entire study period. TA recorded a slightly downward trend in S1, a significantly upward trend in S2 (with a slope of 0.77°C/decade), and an insignificant downward trend in S3 (Figure 3(a), Table 2), results which may be due to alternating surface albedo due to large areas of desertification in S2 and rapid revegetation activities in S3 [23, 37, 64–70].Table 2
Trend slopes for average temperature and extreme temperature indices in the BTSSR during the different stages.
Index S1 S2 S3 S Increasing (%) Decreasing (%) Units T A −0.01 0.77∗∗ −0.14 0.34∗∗∗ 52 96.23 % 1.89% 1 0.00 % 1.89 % °C/decade T s p r i n g −0.01 0.76∗ 0.00 0.38∗∗∗ 52 98.11 % 0.00% 1 0.00 % 1.89 % °C/decade T s u m m e r −0.24 0.35 −0.32 0.24∗∗∗ 52 84.91 % 13.21% 1 0.00 % 1.89 % °C/decade T a u t u m n −0.03 0.64 0.29 0.30∗∗∗ 51 83.02 % 13.21% 2 3.77 % 0.00 % °C/decade T w i n t e r 0.64 1.28∗∗ −1.10 0.45∗∗∗ 53 79.25 % 20.75% 0 0.00 % 0.00 % °C/decade T G −0.02 0.47∗ −0.15 0.27∗∗∗ 52 94.34 % 3.77% 1 1.89 % 0.00 % °C/decade TN10p −0.98 −3.40∗∗ −0.35 −2.10∗∗∗ 3 0.00 % 5.66% 50 81.13 % 13.21 % days/decade TX10p 1.38 −2.19∗ −0.73 −0.96∗ 0 0.00 % 0.00% 53 56.60 % 43.40 % days/decade CSDI −0.22 −0.73∗ −0.81∗ −0.80∗∗ 4 0.00 % 7.55% 49 49.06 % 43.40 % days/decade TNn −0.17 1.61∗ 0.50 0.47∗∗ 48 52.83 % 37.74% 5 1.89 % 7.55 % °C/decade TN90p −0.17 3.49∗ −0.12 1.74∗∗∗ 51 88.68 % 7.55% 2 3.77 % 0.00 % days/decade TX90p −1.06 2.77 −1.37 1.13∗∗ 52 54.72 % 43.40% 1 0.00 % 0.00 % days/decade WSDI −0.27 3.52∗∗ −1.02 1.11∗∗ 53 79.25 % 20.75% 0 0.00 % 0.00 % days/decade TXx −0.23 0.36 −0.95 0.16 48 18.87 % 71.70% 5 0.00 % 9.43 % °C/decade DTR −0.43∗ −0.11 −0.12 −0.17∗∗∗ 10 13.21 % 5.66% 43 66.04 % 15.09 % °C/decade GSL 1.33 4.41 0.62 3.20∗∗∗ 52 79.25 % 18.87% 1 0.00 % 1.89 % days/decade S1 is 1960–1978, S2 is 1979–2000, and S3 is 2000–2014; ∗∗∗ is 0.001 level of significance, ∗∗ is 0.01 level of significance, and ∗ is 0.05 level of significance; percentage of stations with a significant (p<0.05) trend (bold), and percentage of stations with an insignificant (p>0.05) trend (italics).Figure 2
Anomaly of regional average temperature curves over the BTSSR, 1960–2014. Theblackline represents the linear trends for (a) annual, (b) spring, (c) summer, (d) autumn, (e) winter, and (f) growing season. Theredline represents the 5-year moving average.
(a) (b) (c) (d) (e) (f)Figure 3
Box-Plot of Slopes for regional average temperature over the BTSSR in different stages. (a) Annual, (b) spring, (c) summer, (d) autumn, (e) winter, and (f) growing season. S1 is 1960–1978, S2 is 1979–2000, S3 is 2000–2014, and S is 1960–2014. ∗∗∗ is 0.001 level of significance, ∗∗ is 0.01 level of significance, and ∗ is 0.05 level of significance.
(a) (b) (c) (d) (e) (f)The variations inTA trends for the different stages in the five subregions are shown in Figure 4(a) and Table 3. Results show that TA trends in these subregions significantly increased from 1960 to 2014 (p<0.001); the highest TA slope was 0.41°C/decade in the northwest subregion (SR-II) while the east-central subregion (SR-IV) recorded the lowest slope (0.30°C/decade) (Figures 3(a) and 4(a), Table 3). Results for all five subregions recorded a slight fluctuation during S1 and a significant increasing trend in S2, especially in the northwest subregion (SR-II; 0.92°C/decade), a result which may be attributed to serious steppe degradation in this area [71, 72]. During the last decade, TA recorded slightly negative trends, except for the northwest subregion (SR-II).Table 3
Trend slopes (°C/decade) of mean temperature and extreme temperature indices in the five subregions during the different stages.
Stage S1 S2 S3 S Index I II III IV V I II III IV V I II III IV V I II III IV V Units  T A 0.00 0.08 −0.02 −0.02 −0.08 0.74∗ 0.92∗∗ 0.77∗∗ 0.68∗∗ 0.85∗∗ −0.09 0.08 −0.22 −0.18 −0.29 0.31∗∗∗ 0.41∗∗∗ 0.36∗∗∗ 0.30∗∗∗ 0.37∗∗∗ °C/decade T s p r i n g −0.01 −0.18 0.05 0.03 −0.11 0.85∗ 0.82 0.73∗ 0.56 0.91 0.36 0.51 −0.04 −0.58 −0.27 0.36∗∗∗ 0.38∗∗∗ 0.40∗∗∗ 0.34∗∗ 0.41∗∗∗ °C/decade T s u m m e r −0.09 −0.35 −0.42 −0.20 0.04 0.61∗ 0.49 0.38 0.08 −0.01 −0.40 −0.12 −0.36 −0.09 −0.75 0.19∗∗ 0.31∗∗∗ 0.25∗∗∗ 0.22∗∗∗ 0.31∗∗∗ °C/decade T a u t u m n −0.10 0.03 −0.16 0.15 0.13 0.64 0.71 0.62 0.56 0.75 0.43 0.33 0.28 0.19 0.08 0.24∗∗∗ 0.40∗∗∗ 0.29∗∗∗ 0.29∗∗∗ 0.38∗∗∗ °C/decade T w i n t e r 0.54 1.17 0.77 0.42 0.25 1.07∗ 1.37∗ 1.34∗∗ 1.28∗∗ 1.50∗ −1.13 −0.94 −1.09 −1.12 −1.21 0.45∗∗∗ 0.46∗∗ 0.50∗∗∗ 0.41∗∗ 0.37∗ °C/decade T G 0.11 −0.12 −0.07 −0.11 0.08 0.53∗∗ 0.58∗∗ 0.45 0.34 0.49∗ −0.04 0.02 −0.15 −0.28 −0.46 0.23∗∗∗ 0.34∗∗∗ 0.27∗∗∗ 0.26∗∗∗ 0.35∗∗∗ °C/decade TN10p −0.93 −1.62 −0.23 −1.55 −2.07 −3.32 −3.04∗ −4.72∗∗ −2.65∗∗ −2.64∗∗ −0.20 −0.64 0.39 −1.04 −0.53 −1.72∗∗∗ −2.30∗∗∗ −2.30∗∗∗ −2.14∗∗∗ −2.28∗∗ days/decade TX10p 1.10 1.13 1.74 1.29 1.52 −2.38 −2.59∗ −2.20∗ −1.80 −1.95∗ −1.33 −1.20 −0.71 −0.07 0.09 −0.97∗∗ −0.86∗ −1.13∗∗ −0.82∗ −0.81 days/decade CSDI 0.74 −2.64 0.29 −0.35 −0.92 −0.45 −0.66 −0.89∗ −0.83 −0.92∗ −0.61 −0.53 0.18 −2.54∗∗ −1.14∗∗ −0.72∗ −1.38∗∗ −0.54∗∗ −0.81∗∗ −1.04∗∗ days/decade TNn −0.07 0.98 0.09 0.05 0.19 1.03 1.26 1.78∗∗ 2.10∗∗ 2.05∗∗ 0.81 0.90 0.52 0.17 −0.20 0.25 0.55∗∗ 0.60∗∗∗ 0.49∗∗ 0.49∗ °C/decade TN90p −0.33 0.20 −0.58 0.11 0.48 3.96∗ 4.01∗ 3.34∗ 2.96 3.09∗ 1.17 1.96 −1.31 −1.18 −0.48 1.60∗∗∗ 2.07∗∗∗ 1.80∗∗∗ 1.63∗∗∗ 1.70∗∗∗ days/decade TX90p −0.55 −1.44∗ −1.44 −1.07 −0.71 3.29 3.69 2.69 1.98 2.03 −0.56 0.45 −3.16 −1.53 −2.13 1.31 1.07∗ 1.17∗∗∗ 0.87∗ 1.13∗∗ days/decade WSDI 0.58∗∗ −1.57∗∗ −0.32 −0.28 −0.46 5.42∗∗ 3.29∗∗ 4.31∗∗ 0.57∗ 2.45∗∗ 1.20 2.16 4.31 −2.76∗∗ −2.92 1.34∗∗ 1.00∗∗ 1.40∗∗ 0.57∗∗ 0.88∗∗ days/decade TXx 0.33 −0.03 −0.78 −0.42 0.11 0.54 1.02 0.33 −0.16 0.22 −1.22 −0.13 −1.08 −1.12 −0.54 0.11 0.23 0.23 0.05 0.24 °C/decade DTR −0.31 −0.60∗∗ −0.33 −0.54∗∗ −0.58∗∗ −0.04 −0.06 −0.18 −0.10 −0.21 −0.09 −0.01 −0.07 −0.29 −0.14 −0.09 −0.23∗∗∗ −0.17∗∗∗ −0.20∗∗∗ −0.21∗∗∗ °C/decade GSL −4.26 5.92 3.56 1.11 3.90 3.39 5.17 3.03 5.64 7.95 0.14 3.52 −0.84 0.35 3.23 2.92∗∗ 3.42∗∗∗ 3.40∗∗∗ 3.08∗∗∗ 3.28∗∗∗ days/decade S1 is 1960–1978, S2 is 1979–2000, and S3 is 2000–2014; ∗∗∗ is 0.001 level of significance, ∗∗ is 0.01 level of significance, and ∗ is 0.05 level of significance.Figure 4
Box-Plot of Slopes for regional average temperature over the five subregions in different stages. (a) Annual, (b) spring, (c) summer, (d) autumn, (e) winter, and (f) growing season. S1 is 1960–1978, S2 is 1979–2000, S3 is 2000–2014, and S is 1960–2014. I is the southwest subregion (SR-I), II is the northwest subregion (SR-II), III is the south-central subregion (SR-III), IV is the east-central subregion (SR-IV), and V is the northeast subregion (SR-V). ∗∗∗ is 0.001 level of significance, ∗∗ is 0.01 level of significance, and ∗ is 0.05 level of significance.
(a) (b) (c) (d) (e) (f)Results for the spatial distribution ofTA trends (Figure 5(a)) show that most areas of the BTSSR had positive slopes, with the largest occurring in the central-southern mountain area; however, negative slopes were identified in the southeastern region in the Luanhe River Basin. Based on data from 51 stations (96.23% of all stations), the spatial distribution of TA showed a significant increase in the BTSSR (p<0.01) (Table 2), corresponding to the average slope of 0.35°C/decade. A strong increasing trend was observed at the Wutaishan station and a weak negative trend was found at the Chengde station, both meteorological stations being situated in the south-central subregion (SR-III).Figure 5
Spatial patterns of the trend slopes of mean temperature (1960–2014) over the BTSSR and anomaly of subregional indices curves. Theblackline represents linear trends for (a) annual, (b) spring, (c) summer, (d) autumn, (e) winter, and (f) growing season. Theredline represents the 5-year moving average.
(a) (b)  (c) (d)  (e) (f)In general, an increasingTA trend in the BTSSR can accelerate plant growth by altering plant photosynthesis processes [73, 74] and lengthening the growing season [75, 76]. However, increasing TA trends can also increase regional evapotranspiration [26] which, in conjunction with a decreasing trend in precipitation in the BTSSR [77], will increase the risk of drought [78]. The tradeoff of change trends of these meteorological factors is currently unclear.


### 4.1.2. Seasonal Temperature

MK test results indicate that the seasonal temperature of the BTSSR significantly increased (p<0.05) during the study period; the lowest (0.24°C/decade) and highest (0.45°C/decade) rate of temperature rise were recorded in summer and winter, respectively (Figures 2(b)–2(e) and 3(b)–3(e), Table 2). This result was similar to that recorded by Wang et al. [79] for seasonal temperature changes across the arid region of northwestern China. Results for the different stages of the BTSSR showed differing changes. Results during S1 recorded seasonal temperature to have a slightly negative trend, except for winter temperature (Twinter). Seasonal temperature in S2 increased, with Twinter recording a significant increase (having a slope of 1.28°C/decade) and Tsummer having an insignificant increase (with a slope of 0.35°C/decade), results which were in agreement with seasonal temperature changes across northeastern China [80]. Seasonal temperatures during S3 recorded a mixture of slightly negative and positive trends (Table 2).Results for seasonal temperature trends across the five subregions (Figures4(b)–4(e), Table 3) recorded a significant increase (p<0.05) during the study period, especially for Twinter in SR-III (having a slope of 0.50°C/decade) (Figures 5(b)–5(e)). During S1, insignificant increasing trends for Twinter were recorded in all subregions, while other seasonal temperature series showed a mixture of slightly negative and positive trends across the five subregions. During S2, Twinter showed a significant increasing trend in all subregions, especially in SR-V (with a slope of 1.50°C/decade), and other seasonal temperature series showed insignificant positive trends for most subregions. Insignificant decreasing trends were recorded in the five subregions during the winter months and summer months, while spring temperatures showed a mixture of slightly negative and positive trends, with more negative than positive trends being recorded, and slightly increasing trends were found in all five subregions during the autumn months over the last decade (Figures 4(b)–4(e), Table 3).The spatial distribution of seasonal temperature trends (Figures5(b)–5(e)) shows that a positive temperature slope was identified for almost the entire BTSSR in the four seasons; the southeastern area of SR-III and the eastern area of SR-I were exceptions, these areas recording negative Tspring, Tsummer, and Tautumn slopes and a negative Tautum slope, respectively (Figures 5(b)–5(e)). Most stations exhibited increasing trends, with the lowest winter slope of 0.01°C/decade (Hequ station) and the highest winter slope of 1.51°C/decade (Wutaishan station) being located in SR-I and SR-III, respectively. In spring, data from 98.11% of stations recorded a significant increasing trend (Figure 5(b), Table 2). In summer and autumn, the number of stations recording significant positive trends was similar, the trends mainly occurring in SR-I, SR-III, and SR-IV. Similarly, the number of stations recording negative trends was few (Figures 5(c) and 5(d), Table 2). In the winter, 42 stations recorded significant positive trends (79.25% of all stations), these being mainly located in SR-I, SR-III, and SR-IV, with no station recording a negative trend (Figure 5(e), Table 2). This result implies that the southern and east-central areas of the BTSSR experienced a warmer climate each year.The winter warming trends identified in our study across the BTSSR are similar to those previously identified in Europe [81], North Asia and North America [82], subarctic [83], China [84, 85], and other areas across the globe [86]. Winter warming will undoubtedly be beneficial in alleviating vegetation damage due to freezing temperatures, to both native and planted vegetation, and it will increase the possibility of using low cold-tolerant plants to revegetate the BTSSR.


### 4.1.3. Growing Season Temperature

MK test results for growing season temperature (TG) from 1960 to 2014 across the BTSSR recorded a significant increase (p<0.001), with a slope of 0.27°C/decade (Figures 2(f) and 3(f), Table 2). Results during the different stages showed TG to record a slight decrease in S1, a significant increase in S2 (with a slope of 0.47°C/decade), and an insignificant decline in S3 (Figure 3(f), Table 2). The overall trend of TG was similar with that of TA.Results for the variations inTG trends for the five subregions in the different stages recorded a significant increase (p<0.001; Figure 4(f), Table 3) from 1960 to 2014, especially in SR-V, with a slope of 0.35°C/decade (Figures 5(f) and 4(f), Table 3). As temperatures increased, the date of initiation of the growing season advanced and earlier growth resulted in vegetation covering the land earlier, thus reducing the occurrence of dust storms [87]. However, earlier vegetation growth could increase water stress in the reproductive growth period of water-limited vegetation in the BTSSR [24]. TG trends for all five subregions recorded slight fluctuations during S1. During S2, a slope of 0.58°C/decade was recorded, indicating a significant increasing trend, especially in SR-II. This result was consistent with previous findings showing TG to have significantly increased across most regions of China since 1980 [88]. During the last decade (S3), TG showed slightly negative trends, except in SR-II (Figure 4(f), Table 3).Results for theTG slope showed a similar spatial pattern with those for TA (Figure 5(f)). A larger positive slope was recorded in the central-southern mountain area and the northeast region, and a negative slope occurred in the southeastern area of SR-III. Based on data from 50 stations (94.34% of all stations; Table 2), the spatial distribution of TG showed a significant increase over the BTSSR (p<0.05) (Figure 5(f)), corresponding to an average slope of 0.29°C/decade. A strong increasing trend was observed at the Wutaishan station and a significant negative trend was only found at the Chengde station, both stations being located in SR-III (Figure 5(f)).Our results showed that a trend of increasingTG can enhance plant photosynthetic activities and further productivity without causing water stress [89, 90]. However, in water-limited areas of the BTSSR, increasing temperature in the growing season may exacerbate rain-fed plant stress due to high temperatures combined with decreasing precipitation [77, 91], especially for vegetation planted (shrubs or trees) as part of the degraded land restoration project in the BTSSR.


## 4.1.1. Annual Temperature

MK test results across the study period for the BTSSR show that annual temperature (TA) significantly increased (p<0.001). The TA slope was 0.34°C/decade (Figures 2(a) and 3(a), Table 2), a result which coincides with TA change trends identified in the China territory [59, 60], in regions such as the Yellow River basin, the Yangtze River Basin [61], arid and semiarid regions in northwest China [62], and those in the monsoon area of eastern China [63]. TA results for the different stages showed some differences to those recorded for the entire study period. TA recorded a slightly downward trend in S1, a significantly upward trend in S2 (with a slope of 0.77°C/decade), and an insignificant downward trend in S3 (Figure 3(a), Table 2), results which may be due to alternating surface albedo due to large areas of desertification in S2 and rapid revegetation activities in S3 [23, 37, 64–70].Table 2
Trend slopes for average temperature and extreme temperature indices in the BTSSR during the different stages.
Index S1 S2 S3 S Increasing (%) Decreasing (%) Units T A −0.01 0.77∗∗ −0.14 0.34∗∗∗ 52 96.23 % 1.89% 1 0.00 % 1.89 % °C/decade T s p r i n g −0.01 0.76∗ 0.00 0.38∗∗∗ 52 98.11 % 0.00% 1 0.00 % 1.89 % °C/decade T s u m m e r −0.24 0.35 −0.32 0.24∗∗∗ 52 84.91 % 13.21% 1 0.00 % 1.89 % °C/decade T a u t u m n −0.03 0.64 0.29 0.30∗∗∗ 51 83.02 % 13.21% 2 3.77 % 0.00 % °C/decade T w i n t e r 0.64 1.28∗∗ −1.10 0.45∗∗∗ 53 79.25 % 20.75% 0 0.00 % 0.00 % °C/decade T G −0.02 0.47∗ −0.15 0.27∗∗∗ 52 94.34 % 3.77% 1 1.89 % 0.00 % °C/decade TN10p −0.98 −3.40∗∗ −0.35 −2.10∗∗∗ 3 0.00 % 5.66% 50 81.13 % 13.21 % days/decade TX10p 1.38 −2.19∗ −0.73 −0.96∗ 0 0.00 % 0.00% 53 56.60 % 43.40 % days/decade CSDI −0.22 −0.73∗ −0.81∗ −0.80∗∗ 4 0.00 % 7.55% 49 49.06 % 43.40 % days/decade TNn −0.17 1.61∗ 0.50 0.47∗∗ 48 52.83 % 37.74% 5 1.89 % 7.55 % °C/decade TN90p −0.17 3.49∗ −0.12 1.74∗∗∗ 51 88.68 % 7.55% 2 3.77 % 0.00 % days/decade TX90p −1.06 2.77 −1.37 1.13∗∗ 52 54.72 % 43.40% 1 0.00 % 0.00 % days/decade WSDI −0.27 3.52∗∗ −1.02 1.11∗∗ 53 79.25 % 20.75% 0 0.00 % 0.00 % days/decade TXx −0.23 0.36 −0.95 0.16 48 18.87 % 71.70% 5 0.00 % 9.43 % °C/decade DTR −0.43∗ −0.11 −0.12 −0.17∗∗∗ 10 13.21 % 5.66% 43 66.04 % 15.09 % °C/decade GSL 1.33 4.41 0.62 3.20∗∗∗ 52 79.25 % 18.87% 1 0.00 % 1.89 % days/decade S1 is 1960–1978, S2 is 1979–2000, and S3 is 2000–2014; ∗∗∗ is 0.001 level of significance, ∗∗ is 0.01 level of significance, and ∗ is 0.05 level of significance; percentage of stations with a significant (p<0.05) trend (bold), and percentage of stations with an insignificant (p>0.05) trend (italics).Figure 2
Anomaly of regional average temperature curves over the BTSSR, 1960–2014. Theblackline represents the linear trends for (a) annual, (b) spring, (c) summer, (d) autumn, (e) winter, and (f) growing season. Theredline represents the 5-year moving average.
(a) (b) (c) (d) (e) (f)Figure 3
Box-Plot of Slopes for regional average temperature over the BTSSR in different stages. (a) Annual, (b) spring, (c) summer, (d) autumn, (e) winter, and (f) growing season. S1 is 1960–1978, S2 is 1979–2000, S3 is 2000–2014, and S is 1960–2014. ∗∗∗ is 0.001 level of significance, ∗∗ is 0.01 level of significance, and ∗ is 0.05 level of significance.
(a) (b) (c) (d) (e) (f)The variations inTA trends for the different stages in the five subregions are shown in Figure 4(a) and Table 3. Results show that TA trends in these subregions significantly increased from 1960 to 2014 (p<0.001); the highest TA slope was 0.41°C/decade in the northwest subregion (SR-II) while the east-central subregion (SR-IV) recorded the lowest slope (0.30°C/decade) (Figures 3(a) and 4(a), Table 3). Results for all five subregions recorded a slight fluctuation during S1 and a significant increasing trend in S2, especially in the northwest subregion (SR-II; 0.92°C/decade), a result which may be attributed to serious steppe degradation in this area [71, 72]. During the last decade, TA recorded slightly negative trends, except for the northwest subregion (SR-II).Table 3
Trend slopes (°C/decade) of mean temperature and extreme temperature indices in the five subregions during the different stages.
Stage S1 S2 S3 S Index I II III IV V I II III IV V I II III IV V I II III IV V Units  T A 0.00 0.08 −0.02 −0.02 −0.08 0.74∗ 0.92∗∗ 0.77∗∗ 0.68∗∗ 0.85∗∗ −0.09 0.08 −0.22 −0.18 −0.29 0.31∗∗∗ 0.41∗∗∗ 0.36∗∗∗ 0.30∗∗∗ 0.37∗∗∗ °C/decade T s p r i n g −0.01 −0.18 0.05 0.03 −0.11 0.85∗ 0.82 0.73∗ 0.56 0.91 0.36 0.51 −0.04 −0.58 −0.27 0.36∗∗∗ 0.38∗∗∗ 0.40∗∗∗ 0.34∗∗ 0.41∗∗∗ °C/decade T s u m m e r −0.09 −0.35 −0.42 −0.20 0.04 0.61∗ 0.49 0.38 0.08 −0.01 −0.40 −0.12 −0.36 −0.09 −0.75 0.19∗∗ 0.31∗∗∗ 0.25∗∗∗ 0.22∗∗∗ 0.31∗∗∗ °C/decade T a u t u m n −0.10 0.03 −0.16 0.15 0.13 0.64 0.71 0.62 0.56 0.75 0.43 0.33 0.28 0.19 0.08 0.24∗∗∗ 0.40∗∗∗ 0.29∗∗∗ 0.29∗∗∗ 0.38∗∗∗ °C/decade T w i n t e r 0.54 1.17 0.77 0.42 0.25 1.07∗ 1.37∗ 1.34∗∗ 1.28∗∗ 1.50∗ −1.13 −0.94 −1.09 −1.12 −1.21 0.45∗∗∗ 0.46∗∗ 0.50∗∗∗ 0.41∗∗ 0.37∗ °C/decade T G 0.11 −0.12 −0.07 −0.11 0.08 0.53∗∗ 0.58∗∗ 0.45 0.34 0.49∗ −0.04 0.02 −0.15 −0.28 −0.46 0.23∗∗∗ 0.34∗∗∗ 0.27∗∗∗ 0.26∗∗∗ 0.35∗∗∗ °C/decade TN10p −0.93 −1.62 −0.23 −1.55 −2.07 −3.32 −3.04∗ −4.72∗∗ −2.65∗∗ −2.64∗∗ −0.20 −0.64 0.39 −1.04 −0.53 −1.72∗∗∗ −2.30∗∗∗ −2.30∗∗∗ −2.14∗∗∗ −2.28∗∗ days/decade TX10p 1.10 1.13 1.74 1.29 1.52 −2.38 −2.59∗ −2.20∗ −1.80 −1.95∗ −1.33 −1.20 −0.71 −0.07 0.09 −0.97∗∗ −0.86∗ −1.13∗∗ −0.82∗ −0.81 days/decade CSDI 0.74 −2.64 0.29 −0.35 −0.92 −0.45 −0.66 −0.89∗ −0.83 −0.92∗ −0.61 −0.53 0.18 −2.54∗∗ −1.14∗∗ −0.72∗ −1.38∗∗ −0.54∗∗ −0.81∗∗ −1.04∗∗ days/decade TNn −0.07 0.98 0.09 0.05 0.19 1.03 1.26 1.78∗∗ 2.10∗∗ 2.05∗∗ 0.81 0.90 0.52 0.17 −0.20 0.25 0.55∗∗ 0.60∗∗∗ 0.49∗∗ 0.49∗ °C/decade TN90p −0.33 0.20 −0.58 0.11 0.48 3.96∗ 4.01∗ 3.34∗ 2.96 3.09∗ 1.17 1.96 −1.31 −1.18 −0.48 1.60∗∗∗ 2.07∗∗∗ 1.80∗∗∗ 1.63∗∗∗ 1.70∗∗∗ days/decade TX90p −0.55 −1.44∗ −1.44 −1.07 −0.71 3.29 3.69 2.69 1.98 2.03 −0.56 0.45 −3.16 −1.53 −2.13 1.31 1.07∗ 1.17∗∗∗ 0.87∗ 1.13∗∗ days/decade WSDI 0.58∗∗ −1.57∗∗ −0.32 −0.28 −0.46 5.42∗∗ 3.29∗∗ 4.31∗∗ 0.57∗ 2.45∗∗ 1.20 2.16 4.31 −2.76∗∗ −2.92 1.34∗∗ 1.00∗∗ 1.40∗∗ 0.57∗∗ 0.88∗∗ days/decade TXx 0.33 −0.03 −0.78 −0.42 0.11 0.54 1.02 0.33 −0.16 0.22 −1.22 −0.13 −1.08 −1.12 −0.54 0.11 0.23 0.23 0.05 0.24 °C/decade DTR −0.31 −0.60∗∗ −0.33 −0.54∗∗ −0.58∗∗ −0.04 −0.06 −0.18 −0.10 −0.21 −0.09 −0.01 −0.07 −0.29 −0.14 −0.09 −0.23∗∗∗ −0.17∗∗∗ −0.20∗∗∗ −0.21∗∗∗ °C/decade GSL −4.26 5.92 3.56 1.11 3.90 3.39 5.17 3.03 5.64 7.95 0.14 3.52 −0.84 0.35 3.23 2.92∗∗ 3.42∗∗∗ 3.40∗∗∗ 3.08∗∗∗ 3.28∗∗∗ days/decade S1 is 1960–1978, S2 is 1979–2000, and S3 is 2000–2014; ∗∗∗ is 0.001 level of significance, ∗∗ is 0.01 level of significance, and ∗ is 0.05 level of significance.Figure 4
Box-Plot of Slopes for regional average temperature over the five subregions in different stages. (a) Annual, (b) spring, (c) summer, (d) autumn, (e) winter, and (f) growing season. S1 is 1960–1978, S2 is 1979–2000, S3 is 2000–2014, and S is 1960–2014. I is the southwest subregion (SR-I), II is the northwest subregion (SR-II), III is the south-central subregion (SR-III), IV is the east-central subregion (SR-IV), and V is the northeast subregion (SR-V). ∗∗∗ is 0.001 level of significance, ∗∗ is 0.01 level of significance, and ∗ is 0.05 level of significance.
(a) (b) (c) (d) (e) (f)Results for the spatial distribution ofTA trends (Figure 5(a)) show that most areas of the BTSSR had positive slopes, with the largest occurring in the central-southern mountain area; however, negative slopes were identified in the southeastern region in the Luanhe River Basin. Based on data from 51 stations (96.23% of all stations), the spatial distribution of TA showed a significant increase in the BTSSR (p<0.01) (Table 2), corresponding to the average slope of 0.35°C/decade. A strong increasing trend was observed at the Wutaishan station and a weak negative trend was found at the Chengde station, both meteorological stations being situated in the south-central subregion (SR-III).Figure 5
Spatial patterns of the trend slopes of mean temperature (1960–2014) over the BTSSR and anomaly of subregional indices curves. Theblackline represents linear trends for (a) annual, (b) spring, (c) summer, (d) autumn, (e) winter, and (f) growing season. Theredline represents the 5-year moving average.
(a) (b)  (c) (d)  (e) (f)In general, an increasingTA trend in the BTSSR can accelerate plant growth by altering plant photosynthesis processes [73, 74] and lengthening the growing season [75, 76]. However, increasing TA trends can also increase regional evapotranspiration [26] which, in conjunction with a decreasing trend in precipitation in the BTSSR [77], will increase the risk of drought [78]. The tradeoff of change trends of these meteorological factors is currently unclear.


## 4.1.2. Seasonal Temperature

MK test results indicate that the seasonal temperature of the BTSSR significantly increased (p<0.05) during the study period; the lowest (0.24°C/decade) and highest (0.45°C/decade) rate of temperature rise were recorded in summer and winter, respectively (Figures 2(b)–2(e) and 3(b)–3(e), Table 2). This result was similar to that recorded by Wang et al. [79] for seasonal temperature changes across the arid region of northwestern China. Results for the different stages of the BTSSR showed differing changes. Results during S1 recorded seasonal temperature to have a slightly negative trend, except for winter temperature (Twinter). Seasonal temperature in S2 increased, with Twinter recording a significant increase (having a slope of 1.28°C/decade) and Tsummer having an insignificant increase (with a slope of 0.35°C/decade), results which were in agreement with seasonal temperature changes across northeastern China [80]. Seasonal temperatures during S3 recorded a mixture of slightly negative and positive trends (Table 2).Results for seasonal temperature trends across the five subregions (Figures4(b)–4(e), Table 3) recorded a significant increase (p<0.05) during the study period, especially for Twinter in SR-III (having a slope of 0.50°C/decade) (Figures 5(b)–5(e)). During S1, insignificant increasing trends for Twinter were recorded in all subregions, while other seasonal temperature series showed a mixture of slightly negative and positive trends across the five subregions. During S2, Twinter showed a significant increasing trend in all subregions, especially in SR-V (with a slope of 1.50°C/decade), and other seasonal temperature series showed insignificant positive trends for most subregions. Insignificant decreasing trends were recorded in the five subregions during the winter months and summer months, while spring temperatures showed a mixture of slightly negative and positive trends, with more negative than positive trends being recorded, and slightly increasing trends were found in all five subregions during the autumn months over the last decade (Figures 4(b)–4(e), Table 3).The spatial distribution of seasonal temperature trends (Figures5(b)–5(e)) shows that a positive temperature slope was identified for almost the entire BTSSR in the four seasons; the southeastern area of SR-III and the eastern area of SR-I were exceptions, these areas recording negative Tspring, Tsummer, and Tautumn slopes and a negative Tautum slope, respectively (Figures 5(b)–5(e)). Most stations exhibited increasing trends, with the lowest winter slope of 0.01°C/decade (Hequ station) and the highest winter slope of 1.51°C/decade (Wutaishan station) being located in SR-I and SR-III, respectively. In spring, data from 98.11% of stations recorded a significant increasing trend (Figure 5(b), Table 2). In summer and autumn, the number of stations recording significant positive trends was similar, the trends mainly occurring in SR-I, SR-III, and SR-IV. Similarly, the number of stations recording negative trends was few (Figures 5(c) and 5(d), Table 2). In the winter, 42 stations recorded significant positive trends (79.25% of all stations), these being mainly located in SR-I, SR-III, and SR-IV, with no station recording a negative trend (Figure 5(e), Table 2). This result implies that the southern and east-central areas of the BTSSR experienced a warmer climate each year.The winter warming trends identified in our study across the BTSSR are similar to those previously identified in Europe [81], North Asia and North America [82], subarctic [83], China [84, 85], and other areas across the globe [86]. Winter warming will undoubtedly be beneficial in alleviating vegetation damage due to freezing temperatures, to both native and planted vegetation, and it will increase the possibility of using low cold-tolerant plants to revegetate the BTSSR.


## 4.1.3. Growing Season Temperature

MK test results for growing season temperature (TG) from 1960 to 2014 across the BTSSR recorded a significant increase (p<0.001), with a slope of 0.27°C/decade (Figures 2(f) and 3(f), Table 2). Results during the different stages showed TG to record a slight decrease in S1, a significant increase in S2 (with a slope of 0.47°C/decade), and an insignificant decline in S3 (Figure 3(f), Table 2). The overall trend of TG was similar with that of TA.Results for the variations inTG trends for the five subregions in the different stages recorded a significant increase (p<0.001; Figure 4(f), Table 3) from 1960 to 2014, especially in SR-V, with a slope of 0.35°C/decade (Figures 5(f) and 4(f), Table 3). As temperatures increased, the date of initiation of the growing season advanced and earlier growth resulted in vegetation covering the land earlier, thus reducing the occurrence of dust storms [87]. However, earlier vegetation growth could increase water stress in the reproductive growth period of water-limited vegetation in the BTSSR [24]. TG trends for all five subregions recorded slight fluctuations during S1. During S2, a slope of 0.58°C/decade was recorded, indicating a significant increasing trend, especially in SR-II. This result was consistent with previous findings showing TG to have significantly increased across most regions of China since 1980 [88]. During the last decade (S3), TG showed slightly negative trends, except in SR-II (Figure 4(f), Table 3).Results for theTG slope showed a similar spatial pattern with those for TA (Figure 5(f)). A larger positive slope was recorded in the central-southern mountain area and the northeast region, and a negative slope occurred in the southeastern area of SR-III. Based on data from 50 stations (94.34% of all stations; Table 2), the spatial distribution of TG showed a significant increase over the BTSSR (p<0.05) (Figure 5(f)), corresponding to an average slope of 0.29°C/decade. A strong increasing trend was observed at the Wutaishan station and a significant negative trend was only found at the Chengde station, both stations being located in SR-III (Figure 5(f)).Our results showed that a trend of increasingTG can enhance plant photosynthetic activities and further productivity without causing water stress [89, 90]. However, in water-limited areas of the BTSSR, increasing temperature in the growing season may exacerbate rain-fed plant stress due to high temperatures combined with decreasing precipitation [77, 91], especially for vegetation planted (shrubs or trees) as part of the degraded land restoration project in the BTSSR.


## 4.2. Spatiotemporal Trends of Indices

### 4.2.1. Cold-Related Indices

Results for temporal variation of cold nights (TN10p), cold days (TX10p), and cold spell duration indicator (CSDI), representing the frequency and duration of extreme low temperature for the BTSSR, are shown in Figures6(a)–6(c). In general, from 1960 to 2014, the fluctuations of these indices continuously significantly decreased (p<0.05), especially TN10p (having a slope of −2.10 days/decade; Figure 7(a), Table 2). Changes observed for TN10p are more significant than those of TX10p across the BTSSR (Figures 7(a) and 7(b), Table 2); a reduction in TN10p will result in a reduction in frost damage during the night [92, 93]. The results of our study were in accordance with those over the arid region of northwestern China [79] and in Iran [94]. Results for the different stages showed that, during S1, the frequency and duration of extreme low temperature recorded a slightly negative trend, except for TX10p. During S2, the frequency and duration of extreme low temperatures significantly decreased, especially TN10p (with a slope of −3.40 days/decade), and results during S3 showed CSDI to significantly decrease (a slope of −0.81 days/decade), while TN10p and TX10p insignificantly decreased (Figures 7(a)–7(c), Table 2). Figure 6(d) shows the temporal variation of Min Tmin (TNn) across the BTSSR. From 1960 to 2014, the fluctuations of TNn continuously significantly increased, having a slope of 0.47°C/decade (p<0.01). During the different stages, TNn showed a slightly positive trend during S1, a significant increase in S2 (with a slope of 1.61°C/decade), and an insignificant increased in S3 (Figure 7(d), Table 2).Figure 6
Anomaly of extreme temperature indices curves across the BTSSR (1960–2014). Theblack line represents the linear trends for (a) TN10p, (b) TX10p, (c) CSDI, (d) TNn, (e) TN90p, (f) TX90p, (g) WSDI, (h) TXx, (i) DTR, and (j) GSL. Theredline represents the 5-year moving average.
(a) (b) (c) (d) (e) (f) (g) (h) (i) (j)Figure 7
Box-Plot of Slopes for extreme temperature indices across the BTSSR in different stages. (a) TN10p, (b) TX10p, (c) CSDI, (d) TNn, (e) TN90p, (f) TX90p, (g) WSDI, (h) TXx, (i) DTR, and (j) GSL. S1 is 1960–1978, S2 is 1979–2000, S3 is 2000–2014, and S is 1960–2014. ∗∗∗ is 0.001 level of significance, ∗∗ is 0.01 level of significance, and ∗ is 0.05 level of significance.
(a) (b) (c) (d) (e) (f) (g) (h) (i) (j)Variations in the trends for the extreme cold indices at the different stages in the five subregions (Figures8(a)–8(d), Table 3) show that their frequency and duration for the majority of the subregions significantly decreased from 1960 to 2014 (p<0.05). TN10p results in SR-III recording a slope of −2.30 days/decade (Figures 8(a) and 9(a), Table 3) and, for the majority of subregions, TNn significantly increased during the study period (p<0.05), notably in SR-III (with a slope of 0.60°C/decade, Figures 8(d) and 9(d), Table 3). During the period S1, all subregions showed a mixture of negative and positive trends, except for TN10p and TX10p. These fluctuating trends stabilized during S2 as a significant warming trend occurred, especially for TN10p in SR-III (with a slope of −4.72 days/decade), and for TNn in SR-IV (with a slope of 2.10°C/decade), these results being in accordance with observations over China during 1979–1999 [95]. During S3, TN10p, TX10p, and CSDI showed negative trends for most subregions while TNn showed a positive trend for most subregions, especially in SR-II (with a slope of −0.90°C/decade).Figure 8
Box-Plot of Slopes for extreme temperature indices across the five subregions in different stages. (a) TN10p, (b) TX10p, (c) CSDI, (d) TNn, (e) TN90p, (f) TX90p, (g) WSDI, (h) TXx, (i) DTR, and (j) GSL. S1 is 1960–1978, S2 is 1979–2000, S3 is 2000–2014, S is 1960–2014. ∗∗∗ is 0.001 level of significance, ∗∗ is 0.01 level of significance, and ∗ is 0.05 level of significance.
(a) (b) (c) (d) (e) (f) (g) (h) (i) (j)Figure 9
Spatial patterns of the trend slopes of cold-related indices across the BTSSR and anomaly of subregional indices curves, 1960–2014. Theblack line represents the linear trends for (a) TN10p, (b) TX10p, (c) CSDI, and (d) TNn. Theredline represents the 5-year moving average.
(a) (b)  (c) (d)The slope distribution of TN10p showed a spatial disparity across the BTSSR (Figure9(a)). A positive slope was identified in the eastern area of SR-I, the northern area of SR-II, and the southeastern area of SR-III, this being the area around the Chengde station; a negative slope prevailed across the rest of the region where higher negative slopes occurred along the northeast border and northwest border and in the south-central part. Results for the spatial distribution of TN10p, based on data from 43 stations (81.13% of all stations), showed a significant decrease across the BTSSR (p<0.05) (Table 2), corresponding to an average slope of −5.05 days/decade. A strong decreasing trend was observed at the Wutaishan station in SR-III and a weak positive trend was identified at the Sonid Left station (SR-II; Figure 9(a)).A negative slope for TX10p was identified across the entire region, with the largest slope occurring in the central-southern mountain area and in the east Siramulen River Basin (Figure9(b)). Data from 56.60% of stations recorded a significant decreasing trend (p<0.05) with no stations recording an increasing trend (Table 2). The Wutaishan station in SR-III recorded the strongest decreasing trend in the BTSSR (Figure 9(b)).Results for CSDI also recorded a negative slope across the majority of the region, an exception being in the eastern parts of SR-I and SR-III where positive slopes occurred (Figure9(c)). Data from 49.06% of stations recorded a significant decreasing trend, concentrated in the northern part of the study area (Figure 9(c), Table 2); a strong decreasing trend was observed at the Linhe station (SR-I) and a weak positive trend was identified at the Yuanping station (SR-III; Figure 9(c)).The slopes for TNn and TN10p recorded similar patterns (Figure9(d)), with larger positive slopes being recorded along the northeast-northwest border and in the south-central part, and negative slopes in the eastern and southern areas of SR-I, the northern area of SR-II, and the eastern area of SR-III. Based on data from 28 stations (52.83% of all stations; Table 2), the spatial distribution of TNn showed a significant increase across the BTSSR (p<0.05), corresponding to an average slope of 0.75°C/decade. A strong increasing trend was observed at the Wutaishan station (SR-III) and a significant negative trend was identified at the Hequ station (SR-I; Figure 9(d)).Results for the study period across the BTSSR showed a general trend of warming for the cold-related indices, trends which were also identified in northern Mongolia [96], South Africa [97], northeast China [98], and globally [99]. As extreme low temperature events decreased in the BTSSR, vegetation damage due to freezing conditions was mitigated, thus resulting in an expansion of range northwards of low, cold-tolerant plants which can result in the enrichment of plant species in this area.


### 4.2.2. Warm-Related Indices

The frequency and duration of extreme high temperature events across the BTSSR are shown by the results for the temporal variation of warm nights (TN90p), warm days (TX90p), and the warm spell duration indicator (WSDI) (Figures6(e)–6(g)). In general, these indices fluctuated continuously, having a significant increasing trend during the study period (p<0.01). The results for TN90p (with a slope of 1.74 days/decade; Figures 7(e)–7(g), Table 2) were noticeable, these being in accordance with previous findings across China [100]. Results for the different stages showed that the frequency and duration of extreme high temperature recorded an insignificant negative trend during S1. Results during S2 only recorded significant increases for TN90p and WSDI; TN90p had the most noticeable increase with a slope of 3.49 days/decade; during S3, the frequency and duration of extreme high temperature decreased (Figures 7(e)–7(g), Table 2). The temporal variation of Max Tmax (TXx) across the BTSSR (Figure 6(h)) from 1960 to 2014 recorded fluctuations to insignificantly increase, having a slope of 0.16°C/decade. TXx results during the different stages recorded a slightly negative trend in S1, an insignificant increase in S2, and an insignificant decrease in S3 (Figure 7(h), Table 2).Variations in extreme warm indices in the five subregions (Figures8(e)–8(h), Table 3) showed that frequency and duration of extreme high temperature recorded significant increasing trends (1960–2014) for the majority of the subregions, especially TN90p (Figures 10(a)–10(c)). The change of TN90p is more pronounced than that of TX90p, a finding which is in accordance with previous studies [79, 94]. The temporal variation of TXx for the five subregions (Figure 10(d)) shows an increasing trend for all subregions, especially SR-V (with a slope of 0.24°C/decade). As daily maximum temperatures increase, transpiration and plant water stress will be exacerbated which may lead to an increase in water demand for vegetation production [101]. All five subregions showed that there was a combination of negative and positive trends during S1, with more negative trends being identified. During S2, an increasing trend occurred for the majority of the subregions, with more insignificant than significant trends, especially for TN90p in SR-III (with a slope of 3.34 days/decade). These findings are also in accordance with changes of temperature extremes in China as a whole [95]. During the last decade (S3), the five subregions showed a mixture of negative and positive trends, with more negative trends recorded than positive (Figures 8(e)–8(h), Table 3).Figure 10
Spatial patterns of the trend slopes of warm-related indices across the BTSSR and anomaly of subregional indices curves, 1960–2014. Theblack line represents the linear trends for (a) TN90p, (b) TX90p, (c) WSDI, and (d) TXx. Theredline represents the 5-year moving average.
(a) (b)  (c) (d)Results across the BTSSR showed that positive slopes in TN90p occurred across the majority of the region, with the largest trend occurring in the central-southern mountain area. However, negative slopes were identified in the Luanhe River Basin of the southeastern region and the east area of the southwest subregion (SR-I) (Figure10(a)). Based on data from 47 stations (88.68% of all stations) (Table 2), the spatial distribution of TN90p showed a significant increase for the BTSSR (p<0.05), corresponding to an average slope of 4.58 days/decade. Significant decreasing trends (at the 99% confidence level) were identified at the Hequ and Chengde stations (Figure 10(a)), accounting for 3.77% of the total stations, and a strong increasing trend was observed at the Wutaishan station (SR-III).A positive slope for TX90p was identified across almost the entire BTSSR, the only exception being the southeastern area of SR-I where a negative slope was recorded. The larger positive slopes were recorded in the central-southern mountain area, and the number of stations recording positive trends was much larger than those recording negative trends (Figure10(b)). Twenty-nine stations recorded significant increases (54.72% of the total stations), with the Wutaishan station in SR-III recording the strongest increasing trend, and only one station recorded a decrease (Table 2; Figure 10(b)).WSDI results across the region recorded positive slopes, with the largest being in the central-southern mountain area (Figure10(c)). Positive trends were significant at 42 stations (79.25%), these being evenly distributed across the region; there were no negative trends recorded for this variable (Figure 10(c), Table 2). The strongest increasing trend was observed at the Wutaishan station in SR-III (Figure 10(c)).TXx slope distribution results across the BTSSR showed a spatial disparity. Negative slopes were identified in the eastern and southern areas of SR-I, in southeastern areas of SR-IV, and in the northern area of SR-V. Positive slopes were present in the other areas of the study region, with the highest positive slopes occurring in the central-southern mountain area (Figure10(d)). The number of stations recording positive trends was much greater than those recording negative trends; 10 stations recorded significant increases (18.87%), with no significant decreases being recorded at any station (Table 2). The strongest increasing trend was recorded at the Wutaishan station (SR-III; Figure 10(d)).Differences in results between the minimum nighttime temperature (TNn) and the maximum daytime temperature (TXx) (Figures7(d) and 7(h), Table 2) indicated that temperatures across the BTSSR showed a general warming trend, a result which is consistent with global increases [102]. The rate of warming of TXx was lower than that of TNn over all of the subregions (Figures 8(d) and 8(h), Table 3), an occurrence which is not conducive to daytime vegetation growth; this warming will have the opposite effect on nighttime vegetation growth by affecting carbon absorption and consumption [103, 104]. Moreover, the increase of maximum daily temperature may inhibit the growth of grass across the BTSSR; an increase of the minimum temperature may promote grassland growth via the compensatory effect of increasing nocturnal respiration [105], as well as reducing desert vegetation productivity by accelerating plant respiration [106], increasing the rate of nutrient metabolism, and shortening the duration of grain filling [107, 108].


### 4.2.3. Variability Indices

The temporal variation of the diurnal temperature range (DTR) across the BTSSR showed DTR fluctuations to continuously significantly decrease (p<0.001), with a slope of −0.17°C/decade (Figures 6(i) and 7(i), Table 2). During the different stages, DTR showed a significant negative trend during S1 (with a slope of −0.43°C/decade; p<0.05), and insignificant decreases during S2 and S3 (Figure 7(i), Table 2). DTR trend variations for the different stages for the majority of the five subregions (Figure 8(i), Table 3) showed a significant decreasing trend across the study period, especially in SR-II which had a slope of −0.23°C/decade (Figures 8(i) and 11(a), Table 3). All five subregions showed negative trends during S1; SR-II recorded a more significant trend, having a slope of −0.60°C/decade. During S2, insignificant decreasing trends were recorded for all subregions, and during the last decade (S3) the five subregions showed negative trends (Figure 8(i), Table 3). Higher DTR is propitious to the accumulation of biomass [93, 109], therefore, in the BTSSR, as DTR decreased the quality of vegetation may decline.Figure 11
Spatial patterns of the trend slopes of variability indices over the BTSSR and anomaly of subregional indices curves, 1960–2014. Theblack line represents the linear trends for (a) DTR, (b) GSL. Theredline represents 5-year moving average.
(a) (b)The slope distribution of DTR shows a spatial disparity in the BTSSR (Figure11(a)). A positive slope was recorded in the north, east, and southwestern areas of the region, and a negative slope prevailed in the other areas where higher negative slopes occurred along the northeast-northwest border and the central-southern mountain area. Based on data from 35 stations (66.04% of all stations), the spatial distribution of DTR showed a significant decrease across the BTSSR (p<0.05) (Table 2), corresponding to an average slope of −0.29°C/decade. Seven stations (13.21% of total stations) had a significant increasing trend at the 95% confidence level. A strong decreasing trend was observed at the Weixian station (SR-III) and a strong increasing trend was observed at the Hequ station (SR-I; Figure 11(a)).The temporal variation of the growing season length (GSL) across the BTSSR during the study period continuously increased significantly by 3.20 days/decade from 1960 to 2014 (p<0.001) (Figures 6(j) and 7(j), Table 2). Results for the three stages showed that GSL across the BTSSR showed an insignificant positive trend during S1, S2, and S3 (Figure 7(j), Table 2). The variations in GSL trends for the different stages in the five subregions showed significant increasing trends for the majority of the subregions during the study period, especially for GSL results in SR-II (with a slope of 3.42 days/decade; Table 3, Figures 8(j) and 11(b)). As the length of the growing season increased, the potential number of harvests and seasonal yields for perennial forage crops may be promoted at high latitudes [110]. All five subregions showed a mixture of negative and positive trends, with more positive trends recorded during S1. During S2, insignificant increasing trends occurred for all subregions and, during the last decade (S3), the five subregions recorded a mixture of negative and positive trends, with more positive than negative trends recorded (Figure 8(j), Table 3).The slope distribution of GSL shows a spatial disparity across the BTSSR (Figure11(b)). Negative slopes were recorded in the eastern areas of SR-III and positive slopes were recorded for all of the other areas in the region; the highest positive slopes occurred in the central-southern mountain area. The number of stations recording positive trends was much greater than those recording negative trends; 42 stations recorded significant increases (79.25%) and no station recorded a significant decrease (Table 2). A strong increasing trend was observed at the Wutaishan station (SR-III) and a negative trend was only found at the Chengde station (SR-III; Figure 11(b)).


## 4.2.1. Cold-Related Indices

Results for temporal variation of cold nights (TN10p), cold days (TX10p), and cold spell duration indicator (CSDI), representing the frequency and duration of extreme low temperature for the BTSSR, are shown in Figures6(a)–6(c). In general, from 1960 to 2014, the fluctuations of these indices continuously significantly decreased (p<0.05), especially TN10p (having a slope of −2.10 days/decade; Figure 7(a), Table 2). Changes observed for TN10p are more significant than those of TX10p across the BTSSR (Figures 7(a) and 7(b), Table 2); a reduction in TN10p will result in a reduction in frost damage during the night [92, 93]. The results of our study were in accordance with those over the arid region of northwestern China [79] and in Iran [94]. Results for the different stages showed that, during S1, the frequency and duration of extreme low temperature recorded a slightly negative trend, except for TX10p. During S2, the frequency and duration of extreme low temperatures significantly decreased, especially TN10p (with a slope of −3.40 days/decade), and results during S3 showed CSDI to significantly decrease (a slope of −0.81 days/decade), while TN10p and TX10p insignificantly decreased (Figures 7(a)–7(c), Table 2). Figure 6(d) shows the temporal variation of Min Tmin (TNn) across the BTSSR. From 1960 to 2014, the fluctuations of TNn continuously significantly increased, having a slope of 0.47°C/decade (p<0.01). During the different stages, TNn showed a slightly positive trend during S1, a significant increase in S2 (with a slope of 1.61°C/decade), and an insignificant increased in S3 (Figure 7(d), Table 2).Figure 6
Anomaly of extreme temperature indices curves across the BTSSR (1960–2014). Theblack line represents the linear trends for (a) TN10p, (b) TX10p, (c) CSDI, (d) TNn, (e) TN90p, (f) TX90p, (g) WSDI, (h) TXx, (i) DTR, and (j) GSL. Theredline represents the 5-year moving average.
(a) (b) (c) (d) (e) (f) (g) (h) (i) (j)Figure 7
Box-Plot of Slopes for extreme temperature indices across the BTSSR in different stages. (a) TN10p, (b) TX10p, (c) CSDI, (d) TNn, (e) TN90p, (f) TX90p, (g) WSDI, (h) TXx, (i) DTR, and (j) GSL. S1 is 1960–1978, S2 is 1979–2000, S3 is 2000–2014, and S is 1960–2014. ∗∗∗ is 0.001 level of significance, ∗∗ is 0.01 level of significance, and ∗ is 0.05 level of significance.
(a) (b) (c) (d) (e) (f) (g) (h) (i) (j)Variations in the trends for the extreme cold indices at the different stages in the five subregions (Figures8(a)–8(d), Table 3) show that their frequency and duration for the majority of the subregions significantly decreased from 1960 to 2014 (p<0.05). TN10p results in SR-III recording a slope of −2.30 days/decade (Figures 8(a) and 9(a), Table 3) and, for the majority of subregions, TNn significantly increased during the study period (p<0.05), notably in SR-III (with a slope of 0.60°C/decade, Figures 8(d) and 9(d), Table 3). During the period S1, all subregions showed a mixture of negative and positive trends, except for TN10p and TX10p. These fluctuating trends stabilized during S2 as a significant warming trend occurred, especially for TN10p in SR-III (with a slope of −4.72 days/decade), and for TNn in SR-IV (with a slope of 2.10°C/decade), these results being in accordance with observations over China during 1979–1999 [95]. During S3, TN10p, TX10p, and CSDI showed negative trends for most subregions while TNn showed a positive trend for most subregions, especially in SR-II (with a slope of −0.90°C/decade).Figure 8
Box-Plot of Slopes for extreme temperature indices across the five subregions in different stages. (a) TN10p, (b) TX10p, (c) CSDI, (d) TNn, (e) TN90p, (f) TX90p, (g) WSDI, (h) TXx, (i) DTR, and (j) GSL. S1 is 1960–1978, S2 is 1979–2000, S3 is 2000–2014, S is 1960–2014. ∗∗∗ is 0.001 level of significance, ∗∗ is 0.01 level of significance, and ∗ is 0.05 level of significance.
(a) (b) (c) (d) (e) (f) (g) (h) (i) (j)Figure 9
Spatial patterns of the trend slopes of cold-related indices across the BTSSR and anomaly of subregional indices curves, 1960–2014. Theblack line represents the linear trends for (a) TN10p, (b) TX10p, (c) CSDI, and (d) TNn. Theredline represents the 5-year moving average.
(a) (b)  (c) (d)The slope distribution of TN10p showed a spatial disparity across the BTSSR (Figure9(a)). A positive slope was identified in the eastern area of SR-I, the northern area of SR-II, and the southeastern area of SR-III, this being the area around the Chengde station; a negative slope prevailed across the rest of the region where higher negative slopes occurred along the northeast border and northwest border and in the south-central part. Results for the spatial distribution of TN10p, based on data from 43 stations (81.13% of all stations), showed a significant decrease across the BTSSR (p<0.05) (Table 2), corresponding to an average slope of −5.05 days/decade. A strong decreasing trend was observed at the Wutaishan station in SR-III and a weak positive trend was identified at the Sonid Left station (SR-II; Figure 9(a)).A negative slope for TX10p was identified across the entire region, with the largest slope occurring in the central-southern mountain area and in the east Siramulen River Basin (Figure9(b)). Data from 56.60% of stations recorded a significant decreasing trend (p<0.05) with no stations recording an increasing trend (Table 2). The Wutaishan station in SR-III recorded the strongest decreasing trend in the BTSSR (Figure 9(b)).Results for CSDI also recorded a negative slope across the majority of the region, an exception being in the eastern parts of SR-I and SR-III where positive slopes occurred (Figure9(c)). Data from 49.06% of stations recorded a significant decreasing trend, concentrated in the northern part of the study area (Figure 9(c), Table 2); a strong decreasing trend was observed at the Linhe station (SR-I) and a weak positive trend was identified at the Yuanping station (SR-III; Figure 9(c)).The slopes for TNn and TN10p recorded similar patterns (Figure9(d)), with larger positive slopes being recorded along the northeast-northwest border and in the south-central part, and negative slopes in the eastern and southern areas of SR-I, the northern area of SR-II, and the eastern area of SR-III. Based on data from 28 stations (52.83% of all stations; Table 2), the spatial distribution of TNn showed a significant increase across the BTSSR (p<0.05), corresponding to an average slope of 0.75°C/decade. A strong increasing trend was observed at the Wutaishan station (SR-III) and a significant negative trend was identified at the Hequ station (SR-I; Figure 9(d)).Results for the study period across the BTSSR showed a general trend of warming for the cold-related indices, trends which were also identified in northern Mongolia [96], South Africa [97], northeast China [98], and globally [99]. As extreme low temperature events decreased in the BTSSR, vegetation damage due to freezing conditions was mitigated, thus resulting in an expansion of range northwards of low, cold-tolerant plants which can result in the enrichment of plant species in this area.


## 4.2.2. Warm-Related Indices

The frequency and duration of extreme high temperature events across the BTSSR are shown by the results for the temporal variation of warm nights (TN90p), warm days (TX90p), and the warm spell duration indicator (WSDI) (Figures6(e)–6(g)). In general, these indices fluctuated continuously, having a significant increasing trend during the study period (p<0.01). The results for TN90p (with a slope of 1.74 days/decade; Figures 7(e)–7(g), Table 2) were noticeable, these being in accordance with previous findings across China [100]. Results for the different stages showed that the frequency and duration of extreme high temperature recorded an insignificant negative trend during S1. Results during S2 only recorded significant increases for TN90p and WSDI; TN90p had the most noticeable increase with a slope of 3.49 days/decade; during S3, the frequency and duration of extreme high temperature decreased (Figures 7(e)–7(g), Table 2). The temporal variation of Max Tmax (TXx) across the BTSSR (Figure 6(h)) from 1960 to 2014 recorded fluctuations to insignificantly increase, having a slope of 0.16°C/decade. TXx results during the different stages recorded a slightly negative trend in S1, an insignificant increase in S2, and an insignificant decrease in S3 (Figure 7(h), Table 2).Variations in extreme warm indices in the five subregions (Figures8(e)–8(h), Table 3) showed that frequency and duration of extreme high temperature recorded significant increasing trends (1960–2014) for the majority of the subregions, especially TN90p (Figures 10(a)–10(c)). The change of TN90p is more pronounced than that of TX90p, a finding which is in accordance with previous studies [79, 94]. The temporal variation of TXx for the five subregions (Figure 10(d)) shows an increasing trend for all subregions, especially SR-V (with a slope of 0.24°C/decade). As daily maximum temperatures increase, transpiration and plant water stress will be exacerbated which may lead to an increase in water demand for vegetation production [101]. All five subregions showed that there was a combination of negative and positive trends during S1, with more negative trends being identified. During S2, an increasing trend occurred for the majority of the subregions, with more insignificant than significant trends, especially for TN90p in SR-III (with a slope of 3.34 days/decade). These findings are also in accordance with changes of temperature extremes in China as a whole [95]. During the last decade (S3), the five subregions showed a mixture of negative and positive trends, with more negative trends recorded than positive (Figures 8(e)–8(h), Table 3).Figure 10
Spatial patterns of the trend slopes of warm-related indices across the BTSSR and anomaly of subregional indices curves, 1960–2014. Theblack line represents the linear trends for (a) TN90p, (b) TX90p, (c) WSDI, and (d) TXx. Theredline represents the 5-year moving average.
(a) (b)  (c) (d)Results across the BTSSR showed that positive slopes in TN90p occurred across the majority of the region, with the largest trend occurring in the central-southern mountain area. However, negative slopes were identified in the Luanhe River Basin of the southeastern region and the east area of the southwest subregion (SR-I) (Figure10(a)). Based on data from 47 stations (88.68% of all stations) (Table 2), the spatial distribution of TN90p showed a significant increase for the BTSSR (p<0.05), corresponding to an average slope of 4.58 days/decade. Significant decreasing trends (at the 99% confidence level) were identified at the Hequ and Chengde stations (Figure 10(a)), accounting for 3.77% of the total stations, and a strong increasing trend was observed at the Wutaishan station (SR-III).A positive slope for TX90p was identified across almost the entire BTSSR, the only exception being the southeastern area of SR-I where a negative slope was recorded. The larger positive slopes were recorded in the central-southern mountain area, and the number of stations recording positive trends was much larger than those recording negative trends (Figure10(b)). Twenty-nine stations recorded significant increases (54.72% of the total stations), with the Wutaishan station in SR-III recording the strongest increasing trend, and only one station recorded a decrease (Table 2; Figure 10(b)).WSDI results across the region recorded positive slopes, with the largest being in the central-southern mountain area (Figure10(c)). Positive trends were significant at 42 stations (79.25%), these being evenly distributed across the region; there were no negative trends recorded for this variable (Figure 10(c), Table 2). The strongest increasing trend was observed at the Wutaishan station in SR-III (Figure 10(c)).TXx slope distribution results across the BTSSR showed a spatial disparity. Negative slopes were identified in the eastern and southern areas of SR-I, in southeastern areas of SR-IV, and in the northern area of SR-V. Positive slopes were present in the other areas of the study region, with the highest positive slopes occurring in the central-southern mountain area (Figure10(d)). The number of stations recording positive trends was much greater than those recording negative trends; 10 stations recorded significant increases (18.87%), with no significant decreases being recorded at any station (Table 2). The strongest increasing trend was recorded at the Wutaishan station (SR-III; Figure 10(d)).Differences in results between the minimum nighttime temperature (TNn) and the maximum daytime temperature (TXx) (Figures7(d) and 7(h), Table 2) indicated that temperatures across the BTSSR showed a general warming trend, a result which is consistent with global increases [102]. The rate of warming of TXx was lower than that of TNn over all of the subregions (Figures 8(d) and 8(h), Table 3), an occurrence which is not conducive to daytime vegetation growth; this warming will have the opposite effect on nighttime vegetation growth by affecting carbon absorption and consumption [103, 104]. Moreover, the increase of maximum daily temperature may inhibit the growth of grass across the BTSSR; an increase of the minimum temperature may promote grassland growth via the compensatory effect of increasing nocturnal respiration [105], as well as reducing desert vegetation productivity by accelerating plant respiration [106], increasing the rate of nutrient metabolism, and shortening the duration of grain filling [107, 108].


## 4.2.3. Variability Indices

The temporal variation of the diurnal temperature range (DTR) across the BTSSR showed DTR fluctuations to continuously significantly decrease (p<0.001), with a slope of −0.17°C/decade (Figures 6(i) and 7(i), Table 2). During the different stages, DTR showed a significant negative trend during S1 (with a slope of −0.43°C/decade; p<0.05), and insignificant decreases during S2 and S3 (Figure 7(i), Table 2). DTR trend variations for the different stages for the majority of the five subregions (Figure 8(i), Table 3) showed a significant decreasing trend across the study period, especially in SR-II which had a slope of −0.23°C/decade (Figures 8(i) and 11(a), Table 3). All five subregions showed negative trends during S1; SR-II recorded a more significant trend, having a slope of −0.60°C/decade. During S2, insignificant decreasing trends were recorded for all subregions, and during the last decade (S3) the five subregions showed negative trends (Figure 8(i), Table 3). Higher DTR is propitious to the accumulation of biomass [93, 109], therefore, in the BTSSR, as DTR decreased the quality of vegetation may decline.Figure 11
Spatial patterns of the trend slopes of variability indices over the BTSSR and anomaly of subregional indices curves, 1960–2014. Theblack line represents the linear trends for (a) DTR, (b) GSL. Theredline represents 5-year moving average.
(a) (b)The slope distribution of DTR shows a spatial disparity in the BTSSR (Figure11(a)). A positive slope was recorded in the north, east, and southwestern areas of the region, and a negative slope prevailed in the other areas where higher negative slopes occurred along the northeast-northwest border and the central-southern mountain area. Based on data from 35 stations (66.04% of all stations), the spatial distribution of DTR showed a significant decrease across the BTSSR (p<0.05) (Table 2), corresponding to an average slope of −0.29°C/decade. Seven stations (13.21% of total stations) had a significant increasing trend at the 95% confidence level. A strong decreasing trend was observed at the Weixian station (SR-III) and a strong increasing trend was observed at the Hequ station (SR-I; Figure 11(a)).The temporal variation of the growing season length (GSL) across the BTSSR during the study period continuously increased significantly by 3.20 days/decade from 1960 to 2014 (p<0.001) (Figures 6(j) and 7(j), Table 2). Results for the three stages showed that GSL across the BTSSR showed an insignificant positive trend during S1, S2, and S3 (Figure 7(j), Table 2). The variations in GSL trends for the different stages in the five subregions showed significant increasing trends for the majority of the subregions during the study period, especially for GSL results in SR-II (with a slope of 3.42 days/decade; Table 3, Figures 8(j) and 11(b)). As the length of the growing season increased, the potential number of harvests and seasonal yields for perennial forage crops may be promoted at high latitudes [110]. All five subregions showed a mixture of negative and positive trends, with more positive trends recorded during S1. During S2, insignificant increasing trends occurred for all subregions and, during the last decade (S3), the five subregions recorded a mixture of negative and positive trends, with more positive than negative trends recorded (Figure 8(j), Table 3).The slope distribution of GSL shows a spatial disparity across the BTSSR (Figure11(b)). Negative slopes were recorded in the eastern areas of SR-III and positive slopes were recorded for all of the other areas in the region; the highest positive slopes occurred in the central-southern mountain area. The number of stations recording positive trends was much greater than those recording negative trends; 42 stations recorded significant increases (79.25%) and no station recorded a significant decrease (Table 2). A strong increasing trend was observed at the Wutaishan station (SR-III) and a negative trend was only found at the Chengde station (SR-III; Figure 11(b)).


## 5. Conclusions

Based on daily temperature data from 53 meteorological stations across the Beijing-Tianjin Sand Source Region (BTSSR), spanning 1960–2014, annual and seasonal temperature temporal and spatial trends were analyzed, as well as extreme temperature, using the Mann-Kendall test, Sen’s slope estimator, and linear regression.Results from 1960 to 2014 showed that a long-term warming trend is evident in this region which is affecting weather conditions. Across the whole area, annual temperature recorded a significant increase of 0.34°C/decade, with the trend test for annual temperature showing that the majority of meteorological stations (96.23% of all stations) exhibited an increasing trend at the 95% confidence level. The highest upward trends were recorded during the winter for the five subregions, with winter temperature in SR-III exhibiting the most significant increasing trend (with a value of 0.50°C/decade).Significant upward trends were also recorded in warm-related extreme temperature indices (TN90p, TX90p, and WSDI), a cold extreme index (TNn) and growing season length (GSL) during the study period, while significant downward trends were identified in cold-related extreme temperature indices (TN10p, TX10p, and CSDI) and for the diurnal temperature range (DTR). In addition, cold nights (TN10p) in SR-III and warm nights (TN90p) in SR-II exhibited higher warming trends with values of −2.30 and 2.07 days/decade, respectively. A higher positive slope for extreme temperature indices was recorded in the west, south-central, and northeastern regions; higher negative slopes occurred in the northwest, south-central, and eastern regions. A higher warming slope was also recorded in extreme temperature along the northeast border, northwest border, and in the central-south mountain part. The trend test in extreme temperature showed that most stations (>79.25% of all stations) exhibited a warming trend.A warming trend for temperature and extremes can result in accelerated plant growth, alleviate plant damage due to freezing, and increase the possibility of using low cold-tolerant plants to revegetate in the BTSSR. As the length of the growing season increases, the potential number of harvests and seasonal yields for perennial forage crops may also be promoted at high latitudes. However, in water-limited areas of the BTSSR, increasing temperatures during the growing season may exacerbate rain-fed plant stress due to higher temperatures and a decrease in precipitation, especially vegetation (shrubs or trees) used for restoration of degraded areas. The increase of the minimum temperature may also reduce productivity of desert vegetation by accelerating plant respiration across the BTSSR.The accurate estimation of changes of temperature and extremes across the BTSSR may be useful for designing and operating appropriate irrigation systems, as well as revegetation strategies, to protect vulnerable arid and semiarid ecosystems against degradation. These findings may be very important for the implementation of future phases of the regeneration project in this area.



---

*Source: 5473105-2018-07-10.xml*