Although China is the largest producer of rice, accounting for about 25% of global production, there are no high-resolution maps of paddy rice covering the entire country. Using time-weighted dynamic time warping (TWDTW), this study developed a pixel- and phenology-based method to identify planting areas of double-season paddy rice in China, by comparing temporal variations of synthetic aperture radar (SAR) signals of unknown pixels to those of known double-season paddy rice fields. We conducted a comprehensive evaluation of the method’s performance at pixel and regional scales. Based on 145,210 field surveyed samples from 2018 to 2020, the producer’s and user’s accuracy are 88.49% and 87.02%, respectively. Compared to county-level statistical data from 2016 to 2019, the relative mean absolute errors are 34.11%. This study produced distribution maps of double-season rice at 10 m spatial resolution from 2016 to 2020 over nine provinces in South China, which account for more than 99% of the planting areas of double-season paddy rice of China. The maps are expected to contribute to timely monitoring and evaluating rice growth and yield.

Sugarcane is an important sugar and biofuel crop with high socio-economic importance, and its planted area has increased rapidly in recent years. China is the world’s third or fourth sugarcane producer. However, to our knowledge, no study has investigated the mapping of sugarcane cultivation areas across entire China. In this study, we developed a phenology-based method to identify sugarcane plantations in China at 30-m spatial resolution from 2016–2020 using the time-series of Landsat and Sentinel-1/2 images derived from Google Earth Engine (GEE) platform. The method worked by comparing the phenological similarity in normalized difference vegetation index (NDVI) series between unknown pixels and sugarcane samples. The phenological similarity was assessed using the time-weighted dynamic time warping method (TWDTW), which has less sensitivity to training samples than machine learning methods and therefore can be easily applied to large areas with limited samples. More importantly, our method introduced multiple and moving time standard phenological curves of sugarcane to the TWDTW by fully considering the variable crop life-cycle of sugarcane, particularly its long harvest season spanning from December to March of the following year. Validations showed the method performed well in 2019, with overall accuracies of 93.47% and 92.74% for surface reflectance (SR) and top of atmosphere reflectance (TOA) data, respectively. The sugarcane maps agreed well with the agricultural statistical areas from 2016–2020. The mapping accuracies using TOA data were comparable to SR data in 2019–2020, but outperformed SR data in 2016–2018 when SR data had lower availability on GEE. The sugarcane maps produced in this study can be used to monitor growing conditions and production of sugarcane and, therefore, can benefit sugarcane management, sustainable sugarcane production, and national food security.

Winter cereals, including wheat, rye, barley, and triticale, are important food crops, and it is crucial to identify the distribution of winter cereals for monitoring crop growth and predicting yield. The production and plating area of winter cereals in Europe both contribute 12.57% to the total global cereal production and plating area in 2020. However, the distribution maps of winter cereals with high spatial resolution are scarce in Europe. Here, we first used synthetic aperture radar (SAR) data from Sentinel-1 A/B, in the Interferometric Wide (IW) swath mode, to distinguish rapeseed and winter cereals; we then used a time-weighted dynamic time warping (TWDTW) method to discriminate winter cereals from other crops by comparing the similarity of seasonal changes in the Normalized Difference Vegetation Index (NDVI) from Landsat and Sentinel-2 images. We generated winter cereal maps for 2016–2020 that cover 32 European countries with 30 m spatial resolution. Validation using field samples obtained from the Google Earth Engine (GEE) platform show that the producer’s and user’s accuracies are 91% ± 7.8% and 89% ± 10.3%, respectively, averaged over 32 countries in Europe. The winter cereal map agrees well with agricultural census data for planted winter cereal areas at municipal and country levels, with the averaged coefficient of determination R² as 0.77 ± 0.15 for 2016–2019. In addition, our method can identify the distribution of winter cereals two months before harvest, with an overall accuracy of 88.4%, indicating that TWDTW is an effective method for timely crop growth monitoring and identification at the continent level. The winter cereal maps in Europe are available via an open-data repository.

Accurate simulations of the spatial and temporal changes in vegetation gross primary production (GPP) play an important role in ecological studies. Previous studies highlighted large uncertainties in GPP datasets based on satellite data with coarse spatial resolutions (>500 m), and implied the need to produce high-spatial-resolution datasets. However, estimating fine spatial resolution GPP is time-consuming and requires an enormous amount of computing storage space. In this study, based on the Eddy Covariance-Light Use Efficiency (EC-LUE) model, we used Google Earth Engine (GEE) to develop a web application (EC-LUE APP) to generate 30-m-spatial-resolution GPP estimates within a region of interest. We examined the accuracy of the GPP estimates produced by the APP and compared them with observed GPP at 193 global eddy covariance sites. The results showed the good performance of the EC-LUE APP in reproducing the spatial and temporal variations in the GPP. The fine-spatial-resolution GPP product (GPP_L) explained 64% of the GPP variations and had fewer uncertainties (root mean square error = 2.34 g C m⁻² d⁻¹) and bias (−0.09 g C m⁻² d⁻¹) than the coarse-spatial-resolution GPP products. In particular, the GPP_L significantly improved the GPP estimations for cropland and dryland ecosystems. With this APP, users can easily obtain 30-m-spatial-resolution GPP at any given location and for any given year since 1984.

Restricted by the design of satellite sensors, the existing satellite-based normalized difference vegetation index (NDVI) cannot simultaneously have a high temporal resolution and spatial resolution, which substantially limits its applications. In recent years, several spatiotemporal fusion models have been developed to produce vegetation index datasets with both high spatial and temporal resolutions, but large uncertainties remain. This study proposes a spatiotemporal fusion model [i.e., Integrating ENvironmental VarIable spatiotemporal fusion (InENVI) model] based on a machine-learning method by incorporating environmental variables to reconstruct NDVI data. Over 14 study areas covering various vegetation types globally, the InENVI method was validated for reproducing spatiotemporal variations in NDVI. On average, the determining coefficients ( $R^{2}$ ) of the reconstructed NDVI compared with satellite-based NDVI observations were above 0.90, reflecting the spatiotemporal variations over all study sites. In addition, we compared the performance of the InENVI model with seven other fusion models over two cropland areas with high vegetation heterogeneity. The results showed that the newly developed InENVI method had the best performance, and the reconstruction error of the InENVI method decreased about 23.68%–59.63% on average over two study areas compared to the other seven methods. Our analyses also highlighted that the integration of environmental variables into spatiotemporal fusion is necessary to improve reconstruction accuracy. The InENVI model provides an alternative approach for reconstructing NDVI datasets with both high spatial and temporal resolutions over large areas.

Gross primary production (GPP) determines the amounts of carbon and energy that enter terrestrial ecosystems. However, the tremendous uncertainty of the GPP still hinders the reliability of GPP estimates and therefore understanding of the global carbon cycle. In this study, using observations from global eddy covariance (EC) flux towers, we appraised the performance of 24 widely used GPP models and the quality of major spatial data layers that drive the models. Results show that global GPP products generated by the 24 models varied greatly in means (from 92.7 to 178.9 Pg C yr⁻¹) and trends (from −0.25 to 0.84 Pg C yr⁻¹). Model structure differences (i.e., light use efficiency models, machine learning models, and process-based biophysical models) are an important aspect contributing to the large uncertainty. In addition, various biases in currently available spatial datasets have found (e.g., only 57% of the observed variation in photosynthetically active radiation at the flux tower locations was explained by the spatial dataset), which not only affect GPP simulation but more importantly hinder the simulation and understanding of the earth system. Moving forward, research into the efficacy of model structures and precision of input data may be more important for global GPP estimation.