Abstract

Based on time series of Moderate Resolution Imaging Spectroradiometer (MODIS) Enhanced Vegetation Index (EVI) data (2000–2009), we extracted forest phenological variables in Northeast China using a threshold-based method, which included the start of the growing season (SOS), end of the growing season (EOS), and length of the growing season (LOS). The spatial variation of phenological trends was analyzed using the linear regression method. In Northeast China, SOS was delayed at the rate of <1.5 days per year. The delay trend of EOS was well distributed in the entire region with almost the same rates. LOS increased slightly. The analysis of the relationship between forest phenology and meteorological variations shows that SOS was mainly affected by spring temperature, whereas SOS had a negative relationship with precipitation in the warm-temperate deciduous broadleaf forest region. The EOS in temperate steppe region was affected by temperature and precipitation in August, whereas the others were significantly affected by temperature. Because of the increased temperature in spring, the LOS of the temperate steppe region and temperate mixed forest region increased, and the LOS was positively correlated with the mean temperature of summer in the cool-temperate needleleaf forest region.

1. Introduction

Phenology has proven to be a sensitive and integrative indicator of climate variability and vegetation growth responses to climate change [1, 2]. The understanding of phenology brings significant insight into both climate and vegetation interactions and their impacts on different spatial and temporal scales [3]. There is growing evidence that the ecological equilibrium has been altered due to global climate change, resulting in changes in vegetation cover over time and in space. Phenology monitoring can serve as an efficient way to understand the interactions between vegetation and climate, and repeated observations from satellite sensors provide the mechanism to move from plant-specific to regional scale studies of phenology. Satellite-derived vegetation indices are commonly used as indicators of vegetation phenology [410]. Justice et al. [4] used Normalized Difference Vegetation Index (NDVI) to qualitatively assess the global phenology of numerous land cover types. Satellites were later used to interpret phenology as an indicator of land cover changes in South America [5] and to detect phenological dynamics in shrublands [6]. White et al. [7] integrated the basic concepts of traditional meteorologically based phenology modeling with intensive satellite phenology observations and produced biome-specific ecosystem phenology models. Duchemin et al. [8] developed a method that consists of a fit of NDVI predicted by line segment to advanced very high resolution radiometer (AVHRR) NDVI time series to monitor two key stages, budburst and senescence, in the phenological cycle of deciduous forests. Zhang et al. [9, 10] used a series of piecewise logistic functions fit to MODIS Vegetation Index (VI) data to monitor four key transition phases of vegetation dynamics at annual temporal scales. Much related research has demonstrated that environmental drivers such as climate, topography, and soil properties affect vegetation dynamics at different spatial and temporal scales, ranging from instant to long-term and from local to regional scales [1118].

Some literatures have described Chinese vegetation phenology research using remotely sensed data. Chen et al. studied the relationship between plant phenology and satellite sensor derived measures of greenness in Eastern China based on field phenological data and NOAA AVHRR data [1923]. Zhang et al. found that green-up dates in the Tibetan Plateau have continuously advanced from 1982 to 2011 [24]. Guo et al. analyzed Global Inventory Modeling and Mapping Studies (GIMMS) NDVI time series between 1982 and 2003 and found that the start of the growing season of vegetation in Northeast China was significantly influenced by spring temperature [25]. Mao et al. demonstrated that precipitation was a major factor in determining the characteristics of phenology in permafrost regions [26]. Most of these previous studies were limited by spatial resolution (from 1 km to 8 km) and temporal resolution. Since the launch of Terra satellite in late 1999, its MODIS sensor on board provided daily coverage with 250 m spatial resolution.

The objective of this study is to analyze the spatial pattern of key forest phenological variables and to explore the relationship between phenology and meteorological variables in Northeast China. In this study, the 500 m, 8-day composite product from the Terra satellite of Earth Observing System was used to calculate Enhanced Vegetation Index (EVI) for 2000–2009. Spatial distribution maps of vegetation phenological variables were established based on 10-year EVI data, and phenological variables were then analyzed at the regional scale. The relationships between phenology and meteorological variables were analyzed after coupling with meteorological data.

2. Data and Methodology

2.1. Study Area

The study area is in Northeast China (115°0–135°52′E and 38°72′–53°55′N), including Jilin, Heilongjiang, and Liaoning Provinces and eastern Inner Mongolia (Figure 1). The climates in Northeast China are warm-temperate, temperate, and cool-temperate. According to the vegetation regionalization map of 2001, Northeast China is categorized into four vegetation zones, that is, a cool-temperate needleleaf forest region, temperate needleleaf and broadleaf mixed-forest region, warm-temperate deciduous broadleaf forest region, and temperate steppe region (Figure 2) (Editorial Board of Vegetation Map of China, 2001). Northeast China has abundant tree species and a variety of forest types [27]. The forests are widely distributed over mountainous terrain (e.g., Daxinganling and Xiaoxinganling Ranges) and show large variation in species composition across latitudinal domains, elevation gradients, and moisture gradients. At low elevations, the dominant tree species include Korean pine (Pinus koraiensis), basswood (Tilia amurensis), oak (Quercus mongolica), painted maple (Acer mono), and ash (Fraxinus mandshurica). At high elevations, the major tree species include spruce (Picea jezoensis var. microsperma), fir (Abies nephrolepis), Mongolia pine (Pinus sylvestris var. mongolica), aspen (Populus davidiana), and birch (Betula platyphylla). Deciduous coniferous larch (Larix gmelinii) forests are widely distributed in the most northern part of Northeast China.

2.2. Materials
2.2.1. MODIS EVI Data

The MODIS Land Science Team provides 8-day composite products for users. The 8-day Land Surface Reflectance product (MOD09A1) from 2000 to 2009 was used in this study. The data were downloaded from the EROS Data Center, US Geological Survey (https://lpdaac.usgs.gov/products/modis_products_table/mod09a1). The MOD09A1 product contains seven spectral bands with a spatial resolution of 500 m. The Enhanced Vegetation Index (EVI) was calculated using three-band reflectance data.

EVI directly adjusts the reflectance in the red band as a function of the reflectance in the blue band, accounting for residual atmospheric contamination (e.g., aerosols), variable soil, and canopy background reflectance [28]. EVI is linearly correlated with the leaf area index and has a higher sensitivity than that of NDVI in areas of high biomass [29]. The equation for EVI is as follows: where , , and ; , , and are the reflectance of the blue, red, and near infrared bands, respectively.

2.2.2. Land Use Data

To study forest phenology individually, the forest zones of Northeast China were extracted from the entire region based on the 1 : 100,000 Land Use Map of China of 2000 from the Data Center for Resources and Environment Sciences, Chinese Academy of Sciences [30].

2.2.3. Climate Data

Air temperature data and rainfall data were acquired from the Chinese Meteorological Data Sharing Service System (http://cdc.cma.gov.cn/), including daily average air temperature and precipitation at 343 meteorological stations within the study area from 2000 to 2009. In this study, four representative stations located in Northeast China were used: Xinlin Station (N51., E124.3) in the cool-temperate needleleaf forest region, Boketu Station (N48.7, E121.9) in the temperate steppe region, Donggang Station (N42., E127.5) in the Temperate needleleaf and broadleaf mixed-forest region, and Yingkou Station (N40.6, E122.1) in the warm-temperate deciduous broadleaf forest region.

2.3. Methods

As shown in Figure 3, the double logistic fitting method was used to reduce noise in the time series MODIS EVI data. The threshold-based method was then used to extract forest phenological variables. The least square method was used to analyze the trend of 10 years of phenological variables. Lastly, the change trend was linked to the meteorological variables.

2.3.1. Double Logistic Fitting Method

The double logistic (D_L) function [31, 32] can be used as a basis function as follows: where and determine the position of the left and right inflection points of the curve, respectively, and and determine the rate of change at the left and right inflection points, respectively. The TIMESAT software [33] was used to fit the asymmetric Gaussians (AG), D_L, and Savizky-Glolay (SG) models. The AG and D_L functions are superior to the SG with regard to noise reduction [34]. The D_L model was used to fit the time series data in this study. After the data preprocessing of the Double Logistic fitting, the output data had litter noise, and the coefficients of fitted logistic functions were saved for each pixel.

2.3.2. Forest Phenological Variables

Considering the regional characteristics of the forests in Northeast China and the basic rationale of the above methods, the threshold-based method was used to identify forest phenology.

EVI percentage thresholds of 0.2 and 0.27 were used to justify the start of the growing season (SOS) and the end of the growing season (EOS), respectively, according to previous studies [3537]. For a series of EVI in a given year, we detected as the maximum EVI and as the minimum EVI in the first half of the year. The SOS EVI value for a given pixel was calculated using the formula Then, the day of EVIstart was determined as the SOS. The EOS was defined using the same method as follows: where is the minimum EVI in the last half of the year.

The processing was performed pixel by pixel to determine EVIstart and EVIend. SOS and EOS were then determined using the double logistic fitting results, which can be determined for a specific day with the above EVI thresholds. By using the percentage thresholds, the absolute EVI value could be adapted to each pixel, which makes the phenology more reasonable.

2.3.3. Trend Analysis and Its Linkage to Meteorological Variations

The linear regression method was used to analyze the trends of phenological variables. The least square method is a commonly used method in the analysis of vegetation growth [38] as follows: where is the number of the years, is the SOS of the -year, and is the slope of the trend line. If , the SOS had a delaying trend during 2000–2009; otherwise, the SOS had an advancing trend. SPSS 16.0 was used for the regression analysis for each pixel.

The interannual variations of phenology variables might be associated with the interannual variations of meteorological variables at regional scales. The connection between the changes in forest phenology (SOS, EOS, and LOS) and meteorological variations was analyzed. The effects of climate change on forest phenology trends were analyzed based on the phenological changes in the four major forest regions.

3. Results

3.1. Spatial Pattern of Forest Phenology

The 10-year annual EVI was calculated for each pixel, and the spatial pattern of forest phenology was generated based on the above-defined thresholds. The spatial distribution of the 10-year average phenological variables (SOS, EOS, and LOS) for Northeast China from 2000 to 2009 is shown in Figure 4. Table 1 shows the maximum, minimum, and average values of the phenological variables in the study area.

Figure 4(a) shows the SOS in Northeast China during DOY (day of the year) 100–140. SOS began earlier in the south, and greening occurred gradually toward the north. This is consistent with the period of tree leaf unfolding in spring. Figure 4(b) shows that the EOS ranges from the DOY 265 to 300, with the end of growth arriving later in lower latitudes, which corresponds to the period of defoliation in autumn. In the cool-temperate needleleaf forest region, there were some patches with an earlier SOS and later EOS, which were related to the evergreen needleleaf species of Mongolian pine. In general, the LOS mainly ranged from DOY 130 to 200 along latitudes from north to south in Northeast China, as shown in Figure 4(c).

The spatial patterns of forest phenology in this study are consistent with previous research results [35, 39] and published field-observed phenological data [1921]. Yu et al. used the threshold-based method to determine the SOS and EOS of Northeast China and found the SOS from DOY 100 to 150, the EOS from DOY 260 to 290, and the LOS mainly from 140–180 DOY [35]. Guo et al. found a DOY of 118–135 for SOS and 252–263 for EOS in Northeast China [39]. Comparison between these results showed that the results derived from the three phenological variables are reasonable in the present study.

3.2. Spatial Change of Phenology Trends

The spatial distribution of trends in the variables SOS and EOS for Northeast China over 10 years from 2000 to 2009 is provided in Figure 5. To evaluate annual trends in the phenological variables, these variables were analyzed at the vegetation zone scale. Figure 5(a) shows that the SOS was delayed with rates of 0~1.5 days per year for most pixels in the region. Although both advanced and delayed EOS were observed in the vegetation regionalization, the delay trend was more pervasive, with almost the same rates (Figure 5(b)). As a result, the LOS increased slightly (Figure 5(c)).

Table 2 demonstrates the linear trends of the area-averaged phenological variables (SOS, EOS, and LOS) in different vegetation zones of Northeast China. The average SOS advanced by 0.06 days per year in the warm-temperate deciduous broadleaf forest region, whereas SOS was delayed in the other three vegetation zones. The average EOS was delayed in all of the vegetation zones. The increased delay in SOS in the cool-temperate needleleaf forest zone caused a decreased LOS (Table 2). The average LOS of the other vegetation zones was increased, and the most positive trend was approximately 0.58 days per year in the warm-temperate deciduous broadleaf forest zone.

3.3. Effects of Climate Change on Forest Phenology Trends
3.3.1. Climate Change

According to the meteorological data, we built the significant trends in meteorological parameters (mean temperature and mean precipitation) for Northeast China over 10 years from 2000 to 2009 (Table 3).

The positive linear trend of temperature in the entire study area showed increasing temperature. The largest increase in temperature was observed in the temperate steppe region, at 0.07°C per year. The precipitation in the warm-temperate deciduous broadleaf forest region decreased at the rate of 0.03 mm per year. The average precipitation in the other three vegetation zones increased by the average rate of 0.03 mm per year.

3.3.2. Phenology Trends near Meteorological Stations

Using a window of 20 × 20 pixels centered at the selected meteorological stations, we calculated the average phenology metrics for each area. Figures 6, 7, and 8 show the trends of the average forest phenology for each meteorological station. These trends near the meteorological stations were consistent with each vegetation region.

3.3.3. Correlation Analysis between Phenology and Meteorological Variations

To analyze the effects of climate change on forest phenology trends, a correlation analysis between meteorological variations and phenology was used.

During 2000–2009, a series of temporal windows of February, March, April, May, spring (from March to May), the winter of the previous year (from December to the next February), and the entire previous year were used to calculate the mean meteorological variables (temperature and precipitation). These meteorological statistics were used as independent variables because they can represent the regional climate environment when SOS occurred in the study area.

Table 4 shows the correlation between meteorological variations and SOS (Table 4). The SOS of the cool-temperate needleleaf forest region had a relatively high negative relationship with the mean temperature of May (, ). In the temperate steppe region, SOS had a negative relationship with the mean spring temperature (, ), especially with the mean temperature of April (, ). For the total needleleaf forest and temperate steppe region, no significant correlation was found between SOS and mean precipitation in all seasons. Thus, the SOS of the cool-temperate needleleaf forest and temperate steppe region was strongly affected by temperature. The SOS of the temperate needleleaf and broadleaf mixed-forest region showed a negative relationship correlation with the mean temperature of spring (, ) and March (, ), whereas it was positively correlated with the mean precipitation of the previous year (, ). In addition, in the warm-temperate deciduous broadleaf forest region, SOS had a strongly negative relationship with the mean temperature of March (, ), and the mean precipitation of March had a negative relationship with the SOS (, ). The SOS of the cool-temperate needleleaf forest and temperate steppe region was separately affected by the mean temperature of May and April.

Similar to the SOS correlation analysis, a series of temporal windows of July–November, summer (from June to August), and autumn (from September to November) were used to calculate the mean meteorological variables (temperature and precipitation) for the EOS correlation analysis. Table 5 shows that the EOS of the cool-temperate needleleaf forest region was significantly positively correlated with the mean temperature of autumn (, ) and the mean temperature of November (, ). In the temperate needleleaf and broadleaf mixed-forest region, the EOS was positively correlated with the mean temperature of October (, ), and the EOS of the deciduous broadleaf forest region had a strongly positive correlation with the mean temperature of autumn (, ). For the above-mentioned three regions, EOS was weakly correlated with the mean precipitation in all seasons. Table 5 shows that EOS was affected by temperature more strongly than precipitation. The temperature and precipitation of August both affected the EOS in the temperate steppe region.

To carry out the LOS correlation analysis, the mean meteorological variables (temperature and precipitation) of spring (from March to May), summer (from June to August), autumn (from September to November), and the entire year were used as the independent variables to calculate the correlation between climate and phenology (Table 6). The results showed that the LOS of the cool-temperate needleleaf forest region was positively correlated with the temperature of summer (, ). In the temperate steppe region and temperate mixed-forest region, the LOS was positively correlated with the temperature of spring (, and , , resp.). For these three regions, LOS was weakly correlated with the mean precipitation in all seasons. In the deciduous broadleaf forest region, the LOS had a positive correlation with the mean temperature of the year (, ), yet it had a negative correlation with the precipitation (, ). Table 6 shows that LOS was positively correlated with temperature throughout the study area.

4. Conclusion and Discussion

In this study, we derived forest phenological variables (start of the growing season, end of the growing season, and length of the growing season) from the MODIS EVI time series data by the threshold-based method. In Northeast China, the average of SOS and EOS are 114 DOY and 278 DOY, respectively. Moreover, the forest phenological variables were found to be related to the distribution of forest types.

Based on the phenological variables extracted from EVI, we built spatial patterns of three forest phenological variables and the linear trends using the linear regression method. In Northeast China, the SOS was delayed, with the rates of 0~1.5 days per year. Although both advanced and delayed EOS were observed, the delay trend was more pervasive, with almost the same rates. As a result, the LOS increased slightly.

The analysis of the relationship between phenology and climate showed that the SOS (start of season) of each region was mainly affected by the spring temperature. Guo et al. came out with the same result using GIMMS NDVI time series between 1982 and 2003, which found that SOS of vegetation in Northeast China was significantly influenced by spring temperature [25]. Only the mean precipitation of March in the warm-temperate deciduous broadleaf forest region had a negative relationship with SOS (, ). Except for the EOS (end of the growing season) of the temperate steppe region, which was affected by the temperature and precipitation of August, the EOS was significantly affected by temperature in the other study areas. The climate of different seasons had different influences in each area. Because of the increased temperature in spring, the LOS of the temperate steppe region and temperate mixed-forest region increased. The LOS of the cool-temperate needleleaf forest region was positively correlated with the temperature of summer (, ), and the LOS of the deciduous broadleaf forest region was affected by both temperature and precipitation. Furthermore, the LOS was positively correlated with the mean temperature of the year (, ) and negatively correlated with precipitation (, ).

The average start of the growing season in the deciduous broadleaf forest region was advanced, but delayed in the needleleaf forest region, temperate steppe region, and temperate needleleaf and broadleaf mixed forest region. Due to the stronger delay in SOS compared to EOS, the LOS was decreased in the needleleaf forest region, a result that is not consistent with the result derived from the NDVI that needleleaf forest has advanced by rate of 2.5 days per year [39]. The discrepancy is most likely due to the different data and derivation method. The availability of methods and data specially developed for extracting phenological characteristics from remotely sensed data has simplified the data processing and made it more efficient. Vegetation phenology dynamics and its responses to meteorological variations can be described with such research. This study is helpful for phenology-linked climate change research and for implementing climate-informed monitoring in the context of adaptive management [40].

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This study was funded by the National Science Foundation of China (no. 41001279), Youth Science Funds of LREIS, CAS and Strategic Priority Research Program—Climate Change: Carbon Budget and Related Issues of the Chinese Academy of Sciences (nos. XDA05050102 and XDA05050602).