The Google Earth Engine (GEE) has emerged as an essential cloud-based platform for land-cover classification as it provides massive amounts of multi-source satellite data and high-performance computation service. This paper proposed an automatic land-cover classification method using time-series Landsat data on the GEE cloud-based platform. The Moderate Resolution Imaging Spectroradiometer (MODIS) land-cover products (MCD12Q1.006) with the International Geosphere–Biosphere Program (IGBP) classification scheme were used to provide accurate training samples using the rules of pixel filtering and spectral filtering, which resulted in an overall accuracy (OA) of 99.2%. Two types of spectral–temporal features (percentile composited features and median composited monthly features) generated from all available Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper Plus (ETM+) data from the year 2010 ± 1 were used as input features to a Random Forest (RF) classifier for land-cover classification. The results showed that the monthly features outperformed the percentile features, giving an average OA of 80% against 77%. In addition, the monthly features composited using the median outperformed those composited using the maximum Normalized Difference Vegetation Index (NDVI) with an average OA of 80% against 78%. Therefore, the proposed method is able to generate accurate land-cover mapping automatically based on the GEE cloud-based platform, which is promising for regional and global land-cover mapping.

Fine resolution land cover information is a vital foundation of Earth science. In this paper, a novel SPECLib-based operational method is presented for the classification of multi-temporal Landsat imagery using reflectance spectra from the spatial-temporal spectral library (SPECLib) for 30 m land-cover mapping for the whole of China. Firstly, using the European Space Agency (ESA) Climate Change Initiative Global Land Cover (CCI_LC) product and the MODIS Version 6 Nadir bidirectional reflectance distribution function adjusted reflectance (NBAR) product (MCD43A4), a global SPECLib with a spatial resolution of 158.85 km (equivalent to 1.43° at the equator) and a temporal resolution of eight days was developed in the sinusoidal projection. Then, the Landsat datacube covering the whole of China was developed using all available observations of Landsat OLI imagery in 2015. Thirdly, the multi-temporal random forest method based on SPECLib was presented to produce an annual land-cover map with 22 land-cover types using the Landsat datacube. Finally, the annual China land-cover map was validated by two different validation systems using approximately 11,000 interpretation points. The mapping results achieved the overall accuracy of 71.3% and 80.7% and the kappa coefficient of 0.664 and 0.757 for the level-2 validation system (19 land-cover types) and the level-1 validation system (nine land-cover types), respectively. Therefore, the case study in China indicates that the proposed SPECLib method is an operational and accurate method for regional/global fine land-cover mapping at a spatial resolution of 30 m.

Land-cover plays an important role in the Earth’s energy balance, the hydrological cycle, and the carbon cycle. Therefore, it is important to evaluate the current global land-cover (GLC) products and to understand the differences between these products so that they can be used effectively in different applications. In this study, three 30-m GLC products, namely GlobeLand30-2010, GLC_FCS30-2015, and FROM_GLC30-2015, were evaluated in terms of areal consistency and spatial consistency using the Land Use/Cover Area frame statistical Survey (LUCAS) reference dataset over the European Union (EU). Given the limitations of the traditional confusion matrix used in accuracy assessment, we adjusted the confusion matrices from sample counts by accounting for the class proportions of the map and reported the standard errors of the descriptive accuracy measures in the accuracy assessment. The results revealed the following. (1) The overall accuracy of the GlobeLand30-2010 product was the highest at 88.90 ± 0.68%; this was followed by GLC_FCS30-2015 (84.33 ± 0.80%) and FROM_GLC2015 (65.31 ± 1.0%). (2) The consistency between the GLC_FCS30-2015 and GlobeLand30-2010 is higher than the consistency between other products, with an area correlation coefficient of 0.930 and a proportion of consistent pixels of 52.41%, respectively. (3) Across the area of the EU, the dominant land-cover types such as forest and cropland are the most consistent across the three products, whereas the spatial consistency for bare land, grassland, shrubland, and wetland is relatively low. (4) The proportion of pixels for which the consistency is low accounts for less than 16.17% of pixels, whereas the proportion of pixels for which the consistency is high accounts for about 39.12%. The disagreement between these products primarily occurs in transitional zones with mixed land cover types or in mountain areas. Overall, the GlobeLand30 and GLC-FCS30 products were found to be the most consistent and to have good classification accuracy in the EU, with the disagreement between the three 30-m GLC products mainly occurring in heterogeneous regions.

Over the past decades, benefiting from the development of computing capacity and the free access to Landsat and Sentinel imagery, several fine-resolution global land cover (GLC) products (with a resolution of 10 m or 30 m) have been developed (GlobeLand30, FROM-GLC30, GLC_FCS30, FROM-GLC10, European Space Agency (ESA) WorldCover and ESRI Land Cover). However, there is still a lack of consistency analysis or comprehensive accuracy assessment using a common validation dataset for these GLC products. In this study, a novel stratified random sampling GLC validation dataset (SRS_Val) containing 79,112 validation samples was developed using a visual interpretation method, significantly increasing the number of samples of heterogeneous regions and rare land-cover types. Then, we quantitatively assessed the accuracy of these six GLC products using the developed SRS_Val dataset at global and regional scales. The results reveal that ESA WorldCover achieved the highest overall accuracy (of 70.54% ± 9%) among the global 10 m land cover products, followed by FROM-GLC10 (68.95% ± 8%) and ESRI Land Cover (58.90% ± 7%) and that GLC_FCS30 had the best overall accuracy (of 72.55% ± 9%) among the global 30 m land cover datasets, followed by GlobeLand30 (69.96% ± 9%) and FROM-GLC30 (66.30% ± 8%). The mapping accuracy of the GLC products decreased significantly with the increased heterogeneity of landscapes, and all GLC products had poor mapping accuracies in countries with heterogeneous landscapes, such as some countries in Central and Southern Africa. Finally, we investigated the consistency of six GLC products from the perspective of area distributions and spatial patterns. It was found that the area consistencies between the five GLC products (except ESRI Land Cover) were greater than 85% and that the six GLC products showed large discrepancies in area consistency for grassland, shrubland, wetlands and bare land. In terms of spatial patterns, the totally inconsistent pixel proportions of the 10 m and 30 m GLC products were 23.58% and 14.12%, respectively, and these inconsistent pixels were mainly distributed in transition zones, complex terrains regions, heterogeneous landscapes, or mixed land-cover types. Therefore, the SRS_Val dataset well supports the quantitative evaluation of fine-resolution GLC products, and the assessment results provide users with quantitative metrics to select GLC products suitable for their needs.

Water color is an important parameter in water quality assessment. However, the existing water color investigations have mostly focused on the lakes with areas greater 1 km ² . In order to improve the understanding of the color of water bodies in China, a cloud-free composite image of China for the summer of 2015 was generated using time-series of Landsat-8 imagery and the best-available-pixel (BAP) compositing algorithm. Then, the first Forel–Ule index (FUI) water color product with a resolution of 30 m was produced for China using the generated BAP composite and the Google Earth Engine computing platform. Finally, the first national-scale assessment of the FUI of natural lakes with an area >0.01 km ² ( N = 60026) was conducted based on the generated FUI product. The generated FUI product was shown to have a high degree of consistency with in situ water surface reflectance-derived FUI (R ² = 0.90, P < 0.001). Also, it had a high degree of consistency with the in situ Secchi depth (R ² = 0.90, P < 0.001) and trophic level index (R ² = 0.62, P < 0.001) datasets. In addition, we found that the most prevalent lake colors in China were yellow (about 49%) and green (about 41%). Besides, the proportion of small lakes (areas < 1 km ² ) found to be yellow was much larger than for large lakes (area ≥ 1 km ² ) (50% against 28%). Our results will provide important information that can be used for preserving and restoring inland water resources.

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.