Many single-land-use simulation models are available to simulate and predict Land Use and Land Cover Change (LUCC). However, few studies have used multiple models to simulate LUCC in the same region. The paper utilizes the CA-Markov model, Land Change Modeler (LCM), and Patch-generating Land Use Simulation model (PLUS) with natural and social driving factors to simulate the LUCC on the Western Sichuan Plateau, using Kappa coefficient, overall accuracy (OA), and Figure of Merit (FoM) to verify the accuracy of the model, and selects a suitable model to predict the LUCC and landscape pattern in the study area from 2020 to 2070. The results are as follows: (1) The LCM has the highest simulation effect, and its Kappa coefficient, OA, and FoM are higher than the other two models. (2) The area of land types other than grassland and wetland will increase from 2020 to 2070. Among them, the grassland area will decrease, but is still most prominent land category in this region. The proportion of wetland areas remains unchanged. The fragmentation degree of forest (F), grassland (GL), shrubland (SL), water bodies (WBs), bare areas (BAs), and permanent ice and snow (PIS) decreases, and the distribution shows a trend of aggregation. The dominance of F and C decreases but still dominates in the landscape. The overall landscape aggregation increased and complexity decreased, and each landscape type’s diversity, evenness, and richness increased, presenting as a more reasonable development. Using multiple models to simulate the LUCC in the same region, and choosing the most suitable local land model is of great significance to scientifically manage and effectively allocate the land resources in the field.

Assessing vegetation phenology is very important for better understanding the impact of climate change on the ecosystem, and many vegetation index datasets from different remote sensors have been used to quantify vegetation phenology from a regional to global perspective. This study mainly analyzes the similarities and differences in phenology derived from GIMMS NDVI3g and MODIS NDVI datasets across different biomes throughout temperate China. We applied three commonly used methods to extract the start and end of the growing season (SOS and EOS) from two datasets between 2000 and 2015, and analyzed the spatio-temporal characteristics and trends of key phenological parameters between these two datasets in temperate China. Results showed that the multi-year mean GIMMS NDVI was higher than MODIS NDVI throughout most of temperate China, and the consistencies between GIMMS NDVI and MODIS NDVI for all biomes in the senescence phase were better than those in the green-up phase. NDVI differences between GIMMS and MODIS resulted in some distinctions between phenology derived from the two datasets. The results of SOS and EOS for three methods also showed wide discrepancies in spatial patterns, especially in SOS. For different biomes, differences of SOS in forests were obviously less than that in shrublands, grasslands-IM, grasslands-QT and meadows, whereas the differences of EOS in forests were relatively greater than that in SOS. Moreover, large differences of phenological trends were found between GIMMS and MODIS datasets from 2000 to 2015 in entire region and different biomes, and it is particularly noteworthy that both SOS and EOS showed a low proportion of the identical significant trends. The results suggested NDVI datasets obtained from GIMMS and MODIS sensors could induce the differences of the inversion of vegetation phenology in some degree due to the differences of instrumental characteristics between these two sensors. These findings highlighted that inter-calibrate datasets derived from different satellite sensors for some biomes (e.g., grasslands) should be needed when analyzing land surface phenology and their trends, and also provided baseline information for choosing different NDVI datasets in subsequent studies on vegetation patterns and dynamics.