In this paper, the random forests method and the support vector machine in machine learning are explored and compared to the traditional statistical-based maximum likelihood method with 126 features from Sentinel-2A images. The spectral reflectance of 12 bands, 96 texture parameters, 7 vegetation indices, and 11 phenological parameters are successfully extracted from Sentinel-2A images in 2017. The classification result shows that the optimal combination of 13 features yields overall accuracies of traditional classification and machine learning classification of 88.96% and 98%, respectively. Short-wave infrared information shows a significant effect on distinguishing rice, corn, and soybean. The water vapor band plays a significant role in distinguishing between corn and rice. In the multiclassification problem, the machine learning methods have robustness with the identification accuracy of greater than 95% for each crop type, whereas the traditional classification result shows imbalanced accuracies for different crops.

Vegetation phenology is a key biological indicator for monitoring terrestrial ecosystems and global change, and regions with the most obvious phenological changes in vegetation are primarily located at high latitudes and altitudes. Over the past three decades, investigations of obvious phenological changes in vegetation at middle and high latitudes in the Northern Hemisphere have provided significant contributions to understanding global climate change. In this study, phenological parameters were extracted from the Global Inventory Modeling and Mapping Studies (GIMMS) Normalized Difference Vegetation Index (NDVI3g) to analyze the spatial and temporal characteristics of vegetation phenological changes above 40°N in the Northern Hemisphere from 1982–2013. The results showed that the start of season (SOS) was significantly advanced (−2.2 ± 0.6 days·decade⁻¹, p < 0.05) and that the end of season (EOS) was slightly delayed (0.78 ± 0.6 days·decade⁻¹, p = 0.21) over the entire study area in the initial 21 years (1982–2002). When the time scale was extended to 2013, the change rate of the SOS and EOS was significantly reduced; in addition, the SOS was delayed (3.2 ± 1.7 days·decade⁻¹, p < 0.05), and the EOS was advanced (4.5 ± 0.9 days·decade⁻¹, p < 0.05) over the entire study area in the last 11 years (2003–2013). The trends of advanced SOS and delayed EOS over the past three decades were slower than those over the initial two decades on a hemispheric scale. The change trends showed obvious variability with different vegetation types and were greater for woody plants than for herbaceous plants. For broad-leaved forest, the SOS was significantly advanced (2 ± 0.5 days·decade⁻¹, p < 0.05) and the EOS was significantly delayed (2.7 ± 0.6 days·decade⁻¹, p < 0.05) from 1982–2013. The trend of delayed EOS was greater than that of advanced SOS for different vegetation types. With respect to the spatial distribution of phenological trends in the Northern Hemisphere, the trends of advanced SOS and delayed EOS were strongest in Europe followed by North America, and the trends were least significant in Asia. Coniferous forest, shrub forest, grassland, and the entire study area have the same change trends for the two time periods (1982–2002 and 2003–2013), and the increased rate of the phenology parameters has decelerated over the most recent decade. The length of season (LOS) of broad-leaved forest and mixed forest over the past 32 years shows a strong increased trend, and simultaneously, the SOS and EOS show an advanced trend and a delayed trend, respectively

Northeast China is located at high northern latitudes and is a typical region of relatively high sensitivity to global climate change. Studies of the land surface phenology in Northeast China and its response to climate change are important for understanding global climate change. In this study, the land surface phenology parameters were calculated using the third generation dataset from the Global Inventory Modeling and Mapping Studies (GIMMS 3g) that was collected from 1982 to 2013 were estimated to analyze the variations of the land surface phenology in Northeast China at different scales and to discuss the internal relationships between phenology and climate change. We examined the phonological changes of all ecoregions. The average start of the growing season (SOS) did not exhibit a significant trend throughout the study area; however, the end of the growing season (EOS) was significantly delayed by 4.1 days or 0.13 days/year (p < 0.05) over the past 32 years. The SOS for the Hulunbuir Plain, Greater Khingan Mountains and Lesser Khingan Mountains was earlier, and the SOS for the Sanjing, Songnen and Liaohe Plains was later. In addition, the EOS of the Greater Khingan Mountains, Lesser Khingan Mountains and Changbai Mountains was later than the EOS of the Liaohe Plain. The spring temperature had the greatest impact on the SOS. Precipitation had an insignificant impact on forest SOS and a relatively large impact on grassland SOS. The EOS was affected by both temperature and precipitation. Furthermore, although temperature had a lag effect on the EOS, no significant lag effect was observed for the SOS.

Time series of Normalized Difference Vegetation Index (NDVI) derived from multiple satellite sensors are crucial data to study vegetation dynamics. The Land Long Term Data Record Version 4 (LTDR V4) NDVI dataset was recently released at a 0.05 × 0.05° spatial resolution and daily temporal resolution. In this study, annual NDVI time series that are composited by the LTDR V4 and Moderate Resolution Imaging Spectroradiometer (MODIS) NDVI datasets (MOD13C1) are compared and evaluated for the period from 2001 to 2014 in China. The spatial patterns of the NDVI generally match between the LTDR V4 and MOD13C1 datasets. The transitional zone between high and low NDVI values generally matches the boundary of semi-arid and sub-humid regions. A significant and high coefficient of determination is found between the two datasets according to a pixel-based correlation analysis. The spatially averaged NDVI of LTDR V4 is characterized by a much weaker positive regression slope relative to that of the spatially averaged NDVI of the MOD13C1 dataset because of changes in NOAA AVHRR sensors between 2005 and 2006. The measured NDVI values of LTDR V4 were always higher than that of MOD13C1 in western China due to the relatively lower atmospheric water vapor content in western China, and opposite observation appeared in eastern China. In total, 18.54% of the LTDR V4 NDVI pixels exhibit significant trends, whereas 35.79% of the MOD13C1 NDVI pixels show significant trends. Good agreement is observed between the significant trends of the two datasets in the Northeast Plain, Bohai Economic Rim, Loess Plateau, and Yangtze River Delta. By contrast, the datasets contrasted in northwestern desert regions and southern China. A trend analysis of the regression slope values according to the vegetation type shows good agreement between the LTDR V4 and MOD13C1 datasets. This study demonstrates the spatial and temporal consistencies and discrepancies between the AVHRR LTDR and MODIS MOD13C1 NDVI products in China, which could provide useful information for the choice of NDVI products in subsequent studies of vegetation dynamics.

We improved the CASA model based on differences in the types of land use, the values of the maximum light use efficiency, and the calculation methods of solar radiation. Then, the parameters of the model were examined and recombined into 16 cases. We estimated the net primary productivity (NPP) using the NDVI3g dataset, meteorological data, and vegetation classification data from the Greater Khingan Mountain region, China. We assessed the accuracy and temporal-spatial distribution characteristics of NPP in the Greater Khingan Mountain region from 1982 to 2013. Based on a comparison of the results of the 16 cases, we found that different values of maximum light use efficiency affect the estimation more than differences in the fraction of photosynthetically active radiation (FPAR). However, the FPARmax and the constant Tε₂ values did not show marked effects. Different schemes were used to assess different model combinations. Models using a combination of parameters established by scholars from China and the United States produced different results and had large errors. These ideas are meaningful references for the estimation of NPP in other regions. The results reveal that the annual average NPP in the Greater Khingan Mountain region was 760 g C/m²·a in 1982–2013 and that the inter-annual fluctuations were not dramatic. The NPP estimation results of the 16 cases exhibit an increasing trend. In terms of the spatial distribution of the changes, the model indicated that the values in 75% of this area seldom or never increased. Prominent growth occurred in the areas of Taipingling, Genhe, and the Oroqen Autonomous Banner. Notably, NPP decreased in the southeastern region of the Greater Khingan Mountains, the Hulunbuir Pasture Land, and Holingol.

Grassland, as one of the most important ecosystems on Earth, experiences fires that affect the local ecology, economy and society. Notably, grassland fires occur frequently each year in northeastern China. Fire occurrence is a complex problem with multiple causes, such as natural factors, human activities and land use. This paper investigates the disruptive effects of grassland fire in the northeastern Inner Mongolia Autonomous Region of China. In this study, we relied on thermal anomaly detection from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor to identify fire occurrences, and land use data were acquired by Landsat Thematic Mapper/Enhanced Thematic Mapper (TM/ETM). We discussed the relationship between land use and the spatial distribution of grassland fires. The results showed that the impact of land use on grassland fires was significant. Spatially, approximately 80% of grassland fires were clustered within 10 km of cultivated land, and grassland fires generally occurred in areas of intense human activity. The correlation between the spatial distribution of grassland fires and the land use degree in 2000, 2005 and 2010 was high, with R² values of 0.686, 0.716, 0.633, respectively (p < 0.01). These results highlight the importance of the relationship between land use and grassland fire occurrence in the northeastern Inner Mongolia Autonomous Region. This study provides significance for local fire management and prevention.