Crop phenology is an important parameter for crop growth monitoring, yield prediction, and growth simulation. The dynamic threshold method is widely used to retrieve vegetation phenology from remotely sensed vegetation index time series. However, crop growth is not only driven by natural conditions, but also modified through field management activities. Complicated planting patterns, such as multiple cropping, makes the vegetation index dynamics less symmetrical. These impacts are not considered in current approaches for crop phenology retrieval based on the dynamic threshold method. Thus, this paper aimed to (1) investigate the optimal thresholds for retrieving the start of the season (SOS) and the end of the season (EOS) of different crops, and (2) compare the performances of the Normalized Difference Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI) in retrieving crop phenology with a modified version of the dynamic threshold method. The reference data included SOS and EOS ground observations of three major crop types in 2015 and 2016, which includes rice, wheat, and maize. Results show that (1) the modification of the original method ensures a 100% retrieval rate, which was not guaranteed using the original method. The modified dynamic threshold method is more suitable to retrieve crop SOS/EOS because it considers the asymmetry of crop vegetation index time series. (2) It is inappropriate to retrieve SOS and EOS with the same threshold for all crops, and the commonly used 20% or 50% thresholds are not the optimal thresholds for all crops. (3) For single and late rice, the accuracies of the SOS estimations based on EVI are generally higher compared to those based on NDVI. However, for spring maize and summer maize, results based on NDVI give higher accuracies. In terms of EOS, for early rice and summer maize, estimates based on EVI result in higher accuracies, but, for late rice and winter wheat, results based on NDVI are closer to the ground records.

Heading and flowering are two key phenological stages in the growth process of winter wheat. It is of great significance for agricultural management and scientific research to accurately monitor and forecast the heading and flowering dates of winter wheat. However, the monitoring accuracy of existing methods based on remote sensing needs to be improved, and these methods cannot realize forecasting in advance. This study proposed an accumulated temperature method (ATM) for monitoring and forecasting the heading and flowering dates of winter wheat from the perspective of thermal requirements for crop growth. The ATM method consists of three key procedures: (1) extracting the green-up date of winter wheat as the starting point of temperature accumulation with the dynamic threshold method from remotely sensed vegetation index (VI) time-series data, (2) calculating the accumulated temperature and determining the thermal requirements from the green-up date to the heading date or the flowering date based on phenology observation samples, and (3) combining the satellite-derived green-up date, daily temperature data, and thermal requirements to monitor and forecast the heading date and flowering date of winter wheat. When applying the ATM method to winter wheat in the North China Plain during 2017–2019, the root mean square error (RMSE) for the estimated heading date was between 4.76 and 6.13 d and the RMSE for the estimated flowering date was between 5.30 and 6.41 d. By contrast, the RMSE for the heading and flowering dates estimated by the widely used maximum vegetation index method was approximately 10 d. Furthermore, the forecasting accuracy of the ATM method was also high, and the RMSE was approximately 6 d. In summary, the proposed ATM method can be used to accurately monitor and forecast the heading and flowering dates of winter wheat in large spatial scales and it performs better than the existing maximum vegetation index method.

Process-based crop models are extensively used to assess the impacts of climate change, environmental variations, and management practices on crop yields. However, parameters sourced from the literature are often not universally applicable, necessitating calibration to enhance the model performance. The in-situ observed crop yield (hereafter referred to as “observed crop yield”) is commonly used for model calibration. Satellite-based data, such as evapotranspiration (ET), offers additional insights into plant growth and holds significant potential for enhancing calibration efforts. The LPJ-GUESS (Lund-Potsdam-Jena General Ecosystem Simulator) model has been applied extensively to simulate crop yields across various scales, but it has not been calibrated for the regional-scale yield simulation. This study aims to enhance the performance of LPJ-GUESS in simulating crop yields in southern Sweden through calibration using observed crop yield and satellite-based ET data. Results demonstrated that calibrating with observed crop yield substantially improved simulation accuracy for both spring barley and winter wheat, reducing the normalized root mean square error (nRMSE) from 50.3% to 12.8% and from 15.5% to 12.2%, respectively. Sensitivity analysis identified four key parameters influencing yield simulations: minimum C:N ratio (CNmin), N demand reduction by leaves (Ndred), and the retranslocation of nitrogen and carbon (Nret and Cret). Calibration using the Penman-Monteith-Leuning Version 2 (PML-V2) ET product moderately enhanced yield simulation accuracy, particularly for winter wheat, achieving an nRMSE of 14.2%, demonstrating its potential as an alternative when especially when the long-term and continuous observed crop yield is not available. The calibrated LPJ-GUESS model effectively simulated crop yield for both crop types under drought and normal conditions, highlighting its robustness across varying environmental scenarios.

Land cover (LC) determines the energy exchange, water and carbon cycle between Earth's spheres. Accurate LC information is a fundamental parameter for the environment and climate studies. Considering that the LC in China has been altered dramatically with the economic development in the past few decades, sequential and fine-scale LC monitoring is in urgent need. However, currently, fine-resolution annual LC dataset produced by the observational images is generally unavailable for China due to the lack of sufficient training samples and computational capabilities. To deal with this issue, we produced the first Landsat-derived annual China Land Cover Dataset (CLCD) on the Google Earth Engine (GEE) platform, which contains 30 m annual LC and its dynamics of China from 1990 to 2019. We first collected the training samples by combining stable samples extracted from China’s Land-Use/Cover Datasets (CLUD), and visually-interpreted samples from satellite time-series data, Google Earth and Google Map. Using 335,709 Landsat images on the GEE, several temporal metrics were constructed and fed to the random forest classifier to obtain classification results. We then proposed a post-processing method incorporating spatial-temporal filtering and logical reasoning to further improve the spatial-temporal consistency of CLCD. Finally, the overall accuracy of CLCD reached 79.31 % based on 5,463 visually-interpreted samples. A further assessment based on 5,131 third-party test samples showed that the overall accuracy of CLCD outperforms that of MCD12Q1, ESACCI_LC, FROM_GLC, and GlobaLand30. Besides, we intercompared the CLCD with several Landsat-derived thematic products, which exhibited good consistencies with the Global Forest Change, the Global Surface Water, and three impervious surface products. Based on the CLCD, the trends and patterns of China’s LC changes during 1985 and 2019 were revealed, such as expansion of impervious surface (+148.71 %) and water (+18.39 %), decrease of cropland (−4.85 %) and grassland (−3.29 %), increase of forest (+4.34 %). In general, CLCD reflected the rapid urbanization and a series of ecological projects (e.g., Gain for Green) in China and revealed the anthropogenic implications on LC under the condition of climate change, signifying its potential application in the global change research. The CLCD dataset introduced in this article is freely available at http://doi.org/10.5281/zenodo.4417810 (Yang and Huang, 2021).