Variability and Climate Change Trend in Vegetation Phenology of Recent Decades in the Greater Khingan Mountain Area , Northeastern China

Vegetation phenology has been used in studies as an indicator of an ecosystem’s responses to climate change. Satellite remote sensing techniques can capture changes in vegetation greenness, which can be used to estimate vegetation phenology. In this study, a long-term vegetation phenology study of the Greater Khingan Mountain area in Northeastern China was performed by using the Global Inventory Modeling and Mapping Studies (GIMMS) normalized difference vegetation index version 3 (NDVI3g) dataset from the years 1982–2012. After reconstructing the NDVI time series, the start date of the growing season (SOS), the end date of the growing season (EOS) and the length of the growing season (LOS) were extracted using a dynamic threshold method. The response of the variation in phenology with climatic factors was also analyzed. The results showed that the phenology in the study area changed significantly in the three decades between 1982 and 2012, including a 12.1-day increase in the entire region’s average LOS, a 3.3-day advance in the SOS and an 8.8-day delay in the EOS. However, differences existed between the steppe, forest and agricultural regions, with the LOSs of the steppe region, forest region and OPEN ACCESS Remote Sens. 2015, 7 11915 agricultural region increasing by 4.40 days, 10.42 days and 1.71 days, respectively, and a later EOS seemed to more strongly affect the extension of the growing season. Additionally, temperature and precipitation were closely correlated with the phenology variations. This study provides a useful understanding of the recent change in phenology and its variability in this high-latitude study area, and this study also details the responses of several ecosystems to climate change.


Introduction
Terrestrial vegetation occupies 75% of the Earth's land surface [1] and plays an important role in the interaction between the biosphere and the atmosphere [2].Vegetation phenology is an ecosystem function that responds to terrestrial energy exchange, the hydrological cycle and carbon uptake [3].As a result, variability in vegetation phenology is frequently used as an indicator in the study of the dynamic responses of terrestrial ecosystems to climate change [1].
Variability in vegetation phenology has been widely documented on regional and global scales.Recent literature indicates that the Tibetan Plateau [2], Africa [3], Europe [4] and northern high latitudes [5] have all experienced advances in the start date of the growing season (SOS) and increases in the length of the growing season (LOS).However, the magnitudes of these changes vary with location, species, and the investigation period.Extended vegetation growing seasons have mostly been attributed to various climatic factors, particularly global warming induced by the increase in the atmospheric concentration of greenhouse gases [6] or water availability [7].By examining the long-term changes in the growing seasons of temperate vegetation in the Northern Hemisphere for the period of 1982-2008, Jeong et al. [1] found that the SOS advanced by 5.2 days and the end date of the growing season (EOS) became delayed by 4.3 days in the early period (1982)(1983)(1984)(1985)(1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999), whereas the SOS advanced by only 0.2 days and the EOS became delayed by 2.3 days in the later period (2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008).This difference was attributed to the differing effects of warming temperatures [1].Zhang et al. [2] investigated the SOS trends of alpine vegetation on the Tibetan Plateau between 1982 and 2011 and found a continuous SOS advance of 1.04 days/yr, which was consistent with observed warming in spring and winter.By examining the recent trends in plant phenology in the temperate zone of China between 1982 and 1999, Piao et al. [8] found that the growing season duration increased by 1.16 days/yr with the green-up of vegetation advancing by 0.79 days/yr and the dormancy becoming delayed by 0.37 days/yr.These changes were primarily attributed to increases in the mean temperature.A 1 °C warming caused a 7.5-day advance in the onset of green-up and a 3.8-day delay in the vegetation dormancy.Jeganathan et al. [9] estimated the phenology trends of northern high-latitude terrestrial vegetation at latitudes above 45°N and found an advance of 0.58 days/yr in the SOS and a delay of 0.64 days/yr in the EOS.The maximum observed changes were a 1.07 day/yr advance in the SOS of needle-leaved deciduous vegetation and a 1.06 day/yr delay in EOS for broad-leaved deciduous vegetation in the 55-65°N zone, in which the vegetation cover is relatively homogeneous.Thus, monitoring changes in the vegetation phenology can aid our understanding of future ecosystem dynamics and quantify the effects of climate change on terrestrial ecosystems.
Satellite data and field-based surveys are the two main approaches used to study phenology changes [1,10].However, field-based approaches using station-based phenological records may introduce uncertainties due to limited in situ observations and local perspectives at the species level [1,10].With broad spatial coverage and high temporal resolution, the time series data derived from satellite images have been widely used to study vegetation phenology at the landscape, regional, and global scales.Algorithms developed to track the vegetation phenology include the Normalized Difference Vegetation Index (NDVI) [11] time series based on data collected by satellite sensors that include the Advanced Very High Resolution Radiometer (AVHRR), SPOT VEGETATION, and Moderate Resolution Imaging Spectroradiometer (MODIS) [12,13].The newest version of the global GIMMS (global inventory modeling and mapping studies) NDVI (normalized difference vegetation index) dataset version 3, termed NDVI3g, has been released, and it extended the satellite record to more than 30 years from 1981-2013 with data sensed by AVHRR onboard NOAA [14,15].The data have been shown to represent real responses of vegetation to climate variability [16], suggesting a positive potential for the long-term monitoring of vegetation phenology [17].
The Hulunber region is located in a climatic and topographic transition zone between the humid forest landscape in the Greater Khingan Range and the semiarid steppe in the western plateau of Inner Mongolia.The transitional nature of the physical conditions makes this region sensitive to ecological changes, such as climate change and human land use.Research on the phenology variability in this region is important for understanding the effects of climate change on the regional ecosystem in the future.The primary objective of this study is to evaluate the long-term phenological variability in the growing season in Hulunber from 1982 to 2012 by using the GIMMS NDVI3g dataset and climate data to characterize the interactions between climate change and vegetation phenology.

Study Region
Hulunber (between 47°05′-53°20′N and 115°31′-126°04′E) is located in northeast Inner Mongolia in northeastern China, and covers a total area of approximately 2.53 × 10 5 km 2 .Geographically, the region is part of the Mongolian Plateau, and the northeast-southwest-oriented Greater Khingan Range occupies the middle of the Hulunber region.The region can be characterized as a transitional ecotone between the forest ecosystem of the eastern Great Khingan Mountains and plateau steppe of western Inner Mongolia, as well as small areas of farming in the southeast (Figure 1).The climate in this region is governed by a typical continental climate regime with short summers, long winters, and a short frost-free period of only 100-110 days.The Great Khingan Mountains divide the regional climate between a temperate monsoon climate (in the eastern portion) and a temperate continental climate (in the western portion).The annual mean temperature range is −4-4 °C , with temperate ranges of 1.3-2.4°C in the farming region in the eastern Great Khingan Mountains, −2.0-4.3 °C in the northern forest region, and 0.4-3.0°C in the western steppe region.The annual precipitation varies widely from 195 mm in the west to 450 mm of the east and primarily falls between July and August.The forest region receives 240-400 mm, the area west of the Great Khingan Mountains receives 200-240 mm, and the area east of the Great Khingan Mountain receives 400-450 mm.The main vegetation types include coniferous forest, broad-leaved forest, shrub, steppe, meadow, and cultivated plants [18].

GIMMS NDVI3g Data
The NDVI has long been used to monitor terrestrial vegetation dynamics because it is directly related to the photosynthetic capacity and energy absorption of plant canopies [19,20].In this study, the latest GIMMS NDVI3g product was used.The data (1982-2012) was acquired from the Ecological Forecasting Lab at the NASA Ames Research Center [21] with a consistent time series with a 1/12° (approximately 8 km) spatial and a bi-monthly temporal resolution [14,22].The product is derived from the AVHRR satellite series 7, 9, 11, 14, 16, 17, 18 and 19 and has been corrected for radiation, geometry, clouds and volcanic aerosols to improve the data quality: It has an improved calibration and improved cloud masking compared to the older versions [14,23].NDVI3g is the first data set of its kind that is appropriate for long-term studies of land surface trends in vegetation photosynthetic capacity, seasonality or phenology, and climate-vegetation couplings, and it aimed to improve data quality in the high latitudes [24] and is superior in describing the changes in vegetation activity in the Northern Hemisphere [25].

Climate Data and Auxiliary Data
A climate dataset of monthly mean temperature and monthly precipitation from 18 meteorological stations (black points in Figure 1) across eastern Inner Mongolia and western Heilongjiang province was provided by the China Meteorological Data Sharing Service System [26] for the period of 1982-2012.These monthly climate data were then interpolated and rasterized to a spatial resolution of 8 km in an ArcGIS environment using the Kriging method.Global Land Survey Digital Elevation Model (GLSDEM) data with a resolution of 90 m for the city of Hulunber, published by the Global Land Cover Facility at the University of Maryland in 2008 [27], were obtained to extract topographic information.The vegetation types were obtained from a digitized 1:1,000,000 vegetation map of the People's Republic of China acquired from the Data Center for Resources and Environmental Sciences, Chinese Academy of Sciences (RESDC) [28].

NDVI Time Series Pre-Processing
Before extracting the vegetation phenology from the NDVI time-series, an analysis and reconstruction procedure was needed to remove residual atmospherically contaminated and unfavorable sun-sensor-surface viewed pixels in the dataset.Several methods have been developed for this purpose [29][30][31][32][33][34].The HANTS (Harmonic Analysis of NDVI Time Series) algorithm [32,35] is a harmonic analysis based series reconstruction method, and it has been widely used to process time series NDVI images.This method has been proven globally applicable and appealing for phenological studies [36][37][38][39].The HANTS algorithm uses Fourier analyses and iterative curve fitting to remove pronounced outliers in a time series that may result from atmospheric conditions or snow cover and then to construct a smooth curve.Thus, the HANTS algorithm was adopted to reconstruct NDVI time series in this study.To correctly reconstruct the time series, the parameters of the HANTS algorithm need to be set carefully.After reviewing various parameter settings used in previous applications [35,[38][39][40], the detailed parameters were set in this study and are shown in Table 1. Figure 2 shows the HANTS reconstruction of the NDVI time series, which produced a smoothly varying curve.
Table 1.HANTS parameter settings applied in this study.

Determination of Phenological Parameters
Two phenological parameters, the SOS and the EOS were extracted from the NDVI time series using the dynamic threshold method developed by Burgan and Hartford [41].This method has been widely accepted by ecologists and phenologists for use in phenology studies [5,42,43].In the dynamic threshold method, the greenness of the vegetation is indexed by transforming the NDVI data into a NDVIratio (with values between 0 and 1) between the NDVI value at a given time [NDVI (t)] and the minimum NDVI value (NDVImin) for the time span of interest, normalized by the total range of NDVI values during this time span (NDVImax −NDVImin) (Equation ( 1)), the process was implemented per year and per pixel.The SOS is defined as the day of the year (DOY) when the NDVI value exceeds a local threshold NDVIratio, and the EOS is defined as the DOY when the NDVI value drops below the local threshold NDVIratio.Considering the effects of snow cover on the onset vegetation green-up [44][45][46], the threshold was set to one-fifth of the difference between NDVImax and NDVImin (Equation ( 2)), which was determined from previous literature that applied in the same threshold to different ecosystems [44,[47][48][49][50].The LOS is calculated as the difference between the start and end of the growing season.To minimize the influence of small patches in the phenology extraction, we initially clustered pixels into homogenous clusters using a window size of 5; consequently, the window-averaged NDVI values were used instead of the center pixel values.

Analysis
Three vegetation groups were created to investigate the temporal variation in vegetation phenology: (a) the steppe region; (b) the forest region; and (c) the agricultural region.The divisions were based on the vegetation types classified in the vegetation map.The steppe region includes meadow steppe, typical steppe, desert steppe, etc.The forest region includes coniferous forest, broad-leaved forest, deciduous forest, etc.In this study, the phenological parameters were analyzed within the different vegetation types.Temporal trends in the vegetation phenology from 1982-2012 were discussed in Section 4.2, which determined using the linear regression between the annual phenological parameters and the year.The response of vegetation phenology to climatic variables was assessed based on the temporal trends in the monthly temperature and precipitation data for the month prior to the phenological parameters.Additionally, the variability in the phenological parameters with time between different vegetation types was also evaluated.

Spatial Patterns in SOS, EOS and LOS
Figure 3 displays the distributions of the SOS, EOS and LOS in the whole Hulunber region averaged from the 30 years during 1982-2012.In Figure 3a, the majority of SOS dates were widely distributed between the 90th Julian day and the 150th Julian day.The SOS of the forest region (approximate 105th Julian day (±5); 10 April-20 April), was later than that of the agricultural region (approximate 95th Julian day (±10); 25 March-15 April).The steppe region had the latest SOS, which occurred on the 115th Julian day (±15) (10 April-10 May).The EOS showed the opposite trend relative to the SOS and is shown in Figure 3b.The EOS dates ranged from the 245th Julian day to the 305th Julian day.The steppe region had the earliest EOS, which was approximately the 259th Julian day (±15) (28 August-30 September), and the EOS of the forest region (the 278th Julian day (±10); 26 September-15 October) was slightly earlier than the agricultural region (279th Julian day (±10); 27 September-16 October).The LOS depended on both the SOS and EOS, and regions with earlier mean SOS dates featured longer LOS values (Figure 3c).The maximum value was observed in Honghuaerji Forest Park, which is the largest sandy Scotch pine forest in Asia.The LOS of the agricultural region (184 ± 15 days) was longer than that of the forest (172 ± 10 days) or steppe regions (144 ± 30 days).

Temporal Changes in SOS, EOS, LOS
To evaluate the temporal SOS, EOS and LOS trends among the three regions in Hulunber, Figure 4 illustrates the interannual variations of the SOS, EOS and LOS in Hulunber and in its three regions during the period between 1982 and 2012.The SOS trended downward, whereas the LOS trended upward, both of which suggested an earlier spring leaf-out and a longer vegetation growing season.The SOS for the entire study area was 3.30 days earlier in the 2010 than in the 1982.The SOS of the steppe area and farming region advanced by 7.90 days and 2.80 days, respectively, whereas the forested area was delayed by 0.82 days.Additionally, the EOS was delayed by 8.81 days on average during the analysis period.The EOS of the steppe region, forested area and agricultural region advanced by 12.30 days, 9.60 days and 4.51 days, respectively.The opposite long-term trends of the SOS and the EOS led to an increase in the LOS (Figure 4c), with the average LOS of the Hulunber region increasing by an average of 12.11 days between the beginning and end of the analyzed period.The LOS of the steppe region, forest region and agricultural region increased by 4.4 days, 10.42 days and 1.71 days, respectively.The increase in the length of the growing season in the steppe region during 1982-2012 was caused by the earlier SOS and later EOS, and that of the forest region was associated with the later EOS.The lengthened LOS in the agricultural region was caused by both the earlier SOS and the later EOS.Thus, the later EOS was the most significant factor in extending the growing season.

Variation Ratio in the SOS, EOS, and LOS
To determine the phenology variation ratio, the SOS, EOS and LOS values were divided by the time period of 30 years.Figure 5 shows the spatial distributions of the variation ratios of the SOS, EOS and LOS in Hulunber from 1982-2012.The variations in the SOS were negative in 52% of the whole study region, indicating that half of the region experienced an earlier SOS.The changes in the steppe region were greater than those in the forest and agricultural regions.The SOS date advanced in 79% of the steppe region, 32% of the forest region and 61% of the agricultural region.The advanced SOS distributed around Hulun Lake may have been caused by changes in land use or by the interannual variability in precipitation.The changes in the SOS of the steppe region and agricultural region were mainly in the range of −0.5-0 days/yr.Overall, 91% of the entire region exhibited later autumn leaf loss, and 84% of the steppe region exhibited a later EOS.In the other two regions, 95% of the forest region and 87% of the agricultural region exhibited later autumn leaf loss.The LOS increased in 74% of the study region.The LOS decreased in 36% of the forest region and 24% of the agricultural region but increased in 25% of the steppe region by more than 1 day/yr.

Relationships between Phenology and Climate
To explore the variability in the three phenological parameters (SOS, EOS, and LOS) in relation to climatic factors, we compared the average SOS, EOS, and LOS values in the steppe region, forest region and agricultural region, with 14 climatic factors (annual precipitation, annual temperature, seasonal precipitation, seasonal temperature, seasonal cumulative temperature above 10 °C, and seasonal cumulative temperature above 0 °C.The seasons were defined as spring (April-May), summer (June-August), and autumn (September-October)).For the steppe region we found that the correlation coefficients (CCs) between the SOS and the spring precipitation, annual temperature, and spring cumulative temperature above 0 °C are stronger than between the SOS and the other factors (Table 2).Among these three factors, the correlation between the annual temperature and the SOS is the strongest (Figure 6a, CC = −0.513,P = 0.003), suggesting that a warmer year would result in the advance of the SOS.A spring cumulative temperature above 0 °C also affects the steppe green-up.Additionally, a significant correlation between the summer precipitation and the EOS was found (Figure 6b, CC=0.382,P=0.034), and the correlation between the summer temperature and the EOS is also notable (CC = −0.399,P = 0.032).These observations indicate that greater summer precipitation and lower summer temperatures are correlated with later steppe EOS dates.Interestingly, the annual precipitation and the summer precipitation are the most influential factors for the LOS.In the forest region, a significant correlation between the SOS and the spring precipitation (Table 3; Figure 6c, CC = 0.448, P = 0.011) was found, and a negative correlation exists between the SOS and the spring mean temperature (CC = −0.388,P = 0.031).These correlations confirmed that the precipitation and temperature in the spring are the main factors influencing forest green-up and that less rainfall and warmer springs lead to earlier SOS dates.The correlation coefficient between the EOS and the summer temperature is the strongest (Figure 6d, CC = 0.577, P = 0.001), and the correlations between the EOS and the summer precipitation and summer cumulative temperature above 10 °C are also significant.In the forest region, the summer temperature was the main factor driving the delay in the EOS, and the summer precipitation also contributed to the variations in the EOS.For the LOS, the annual temperature was the most influential factor, and precipitation appeared to have only a slight effect on the LOS variation in the forest region.In the agricultural region (Table 4), no significant correlations existed between the phenological parameters and climatic factors.However, several weak relationships are observed, such as between the spring temperature and the SOS, the summer precipitation and the EOS, and the summer cumulative temperature above 10 °C and the EOS.These observations indicate that precipitation was not the driving factor of the variability in the phenology.The temperature may more significantly influence the variations in the phenology.We believe that the phenology of the agricultural region in Hulunber is primarily related to human disturbances, and it might have been a result of ecological projects by the Chinese government, such as Grain for Green, the Natural Forest Conversion Program and the Beijing-Tianjin Sand Source Control Program.Our results showed that in the steppe region, both the temperature and the precipitation led to a lengthened LOS.However, in the forest region, longer durations of vegetation growth were more closely correlated with the temperature than the precipitation.The rate of temperature increase in Hulunber was 0.025 °C /yr, with spring temperatures increasing by 0.016 °C /yr and autumn temperatures increasing by 0.051 °C /yr.Our results also suggest that warmer spring temperatures trigger earlier vegetation green-up dates at a rate of 6.96 days/°C , and warmer autumn temperatures cause later EOS dates at a rate of 5.76 days/°C .

Comparison with Previous Studies
Previous studies using the NDVI time series dataset from AVHRR or MODIS show long-lasting LOS trends in the Inner Mongolia steppe [51], China [10], Africa [3] and Europe [22,52].Zhou et al. [53] analyzed NDVI data from 1981-1999 and found that north of 40°N in Eurasia, the vegetation growing season was extended nearly 18 days, with spring beginning a week earlier and autumn ending 10 days later.Miao et al. [54] analyzed the phenology in the Mongolian Plateau through an inter-comparison of global vegetation datasets: they found that the temporal average SOS showed no significant "earlier spring" onset during 1982-2012 in the middle and northern Mongolian Plateau.Zhang et al. [2] merged GIMMS-based SOS from 1982 to 2000 with SPOT-VGT-based SOS from 2001 to 2011 in the Tibetan Plateau, which showed an advancing SOS trend at a rate of 1.04 d/y from 1982 to 2011.He et al. [55] analyzed the vegetation phenological variation and the relationship with the climatic factors using MODIS NDVI data of Zhejiang province.The results showed that the average LOS was approximately 222 days, SOS displayed a non-significant advancing trend, and the EOS was significantly delayed.The correlation analysis of phenological metrics and climatic factors indicated that SOS and EOS were positively related with mean temperature, whereas precipitation and the humidity index were negatively correlated, though not significant.Myneni et al. [5] and Cayan et al. [56] analyzed the relationship between temperature and LOS, and their results showed that warmer springs in 1980s enabled earlier onsets of plant growth and longer growing seasons in the Northern Hemisphere mid-and high-latitudes.Those studies indicated that LOS was extended in recent years,a change that was more closely correlated with temperature than with precipitation [57,58].Jeong et al. [58] also reported that the correlation of the extended growing season with autumn temperatures was greater than with spring temperature.In our results, the vegetation growing season was extended 12.11 days, with 3.30 days from a slightly early spring and 8.81 days attributed to a late autumn.The correlation analysis of phenological metrics and climatic factors indicated that both temperature and precipitation led to a lengthened LOS in the steppe area and forest region, and warmer spring temperatures trigger earlier vegetation green-up dates at a rate of 6.96 days/°C , and warmer autumn temperatures cause later EOS dates at a rate of 5.76 days/°C .
Most of the changes in the growing seasons during 1982-2012 were generally consistent with previous research.However, many studies have shown that both earlier spring leafing and later autumn senescence were factors [10,59] that led to the extended growing season.Our research indicated that the later EOS date was the most significant factor in extending the growing season during 1982-2012 in Hulunber, thus longer-lasting vegetation growth in the research years can be attributed to extended leaf senescence.

Study Limitations
Uncertainty in our results may result from the reconstructed time series processed by the HANTS algorithm.Several models have been developed in the past to smooth and reconstruct time series vegetation index data, and differences between the models and fine tuning of model parameters lead to potential differences, uncertainty and bias between the results amongst users [29,60].Thus, a comparative analysis of the differences between these methods in Hulunber must be further performed.
This research must also be expanded to a deeper discussion of the response of phenological variation to climate change.In this study, we only used simple relationships between the phenological and climatic parameters, such as mean temperature and precipitation, to examine changes in LOS related to climate change.Therefore, they reflect only spatial correlations.However, to draw conclusions about phenological changes under a changing climate conditions, it would be much more valuable to analyze temporal correlations.This study also does not consider the effects of latitude [61], solar radiation changes, land use and land cover changes [62], rising atmospheric CO2 concentrations, water controls [7] or the unexpected role of winter precipitation [63], among other factors, and these factors will be addressed in follow-up studies.

Conclusions
In this study, the latest GIMMS NDVI3g dataset which has high quality data for the high latitudes [24], was used in Hulunber, China, from 1982 to 2012.The start date of the growing season (SOS), the end date of the growing season (EOS) and the length of growing season (LOS) were derived from this time series using a dynamic threshold method to describe the distribution and variability of vegetation phenology in the study area, and three ecosystem regions were classified in the study region: the steppe region, forest region, and agricultural region.The response of the variations in phenology to climate change in the three eco-regions was also discussed.The results indicated that Hulunber experienced obvious phenological variation over the past three decades, including a 12.1-day increase in the LOS, a 3.3-day advance of SOS and an 8.8-day delay of EOS.However, the steppe, forest and agricultural regions presented different variation patterns, with the LOS of the steppe, forest and agricultural regions increased by 4.4 days, 10.42 days and 1.71 days, respectively.The later EOS was found to be the most significant factor in extending the growing season.Furthermore, temperature and precipitation were closely correlated with the variations in the phenology.The correlation analysis of the phenological metrics and climatic factors indicated that both temperature and precipitation led to a lengthened LOS in the steppe area and forest region, and warmer spring temperatures trigger earlier vegetation green-up dates at a rate of 6.96 days/°C , and warmer autumn temperatures cause later EOS dates at a rate of 5.76 days/°C .This study provides a useful understanding of the recent changes in phenology and its variability in this high-latitude study area, and details the responses of several ecosystems to climate change.

Figure 1 .
Figure 1.Spatial patterns in the digital elevation model (DEM) of the Hulunber region.

Figure 2 .
Figure 2. Illustration of NDVI (normalized difference vegetation index) time series reconstruction using the HANTS algorithm.

Figure 4 .
Figure 4. Interannual phenology variability in the growing season (April-October) for the three main regions and the entire Hulunber area from 1982 to 2012.(a) SOS; (b) EOS; (c) LOS.

Figure 6 .
Figure 6.Relationships between phenological parameters and climatic variables for two regions of the study area.(a) SOS of the steppe region; (b) EOS of the steppe region; (c) SOS of the forest region; (d) EOS of the forest region.

Table 2 .
Correlation coefficients between phenological parameters and climatic factors in the steppe region.

Table 3 .
Correlation coefficients between phenological parameters and climatic factors in the forest region.

Table 4 .
Correlation coefficients between phenological parameters and climatic factors in the agricultural region.