Accurate crop planting area information is of significance for understanding regional food security and agricultural development planning. While increasing numbers of medium resolution satellite imagery and improved classification algorithms have been used for crop mapping, limited efforts have been made in feature selection, despite its vital impacts on crop classification. Furthermore, different crop types have their unique spectral and phenology characteristics; however, the different features of individual crop types have not been well understood and considered in previous studies of crop mapping. Here, we examined an optimized strategy to integrate specific features of individual crop types for mapping an improved crop type layer in the Sanjiang Plain, a new food bowl in China, by using all Sentinel-2 time series images in 2018. First, an automatic spectro-temporal feature selection (ASTFS) method was used to obtain optimal features for individual crops (rice, corn, and soybean), including sorting all features by the global separability indices for each crop and removing redundant features by accuracy changes when adding new features. Second, the ASTFS-based optimized feature sets for individual crops were used to produce three crop probability maps with the Random Forest classifier. Third, the probability maps were then composited into the final crop layer by considering the probability of each crop at every pixel. The resultant crop layer showed an improved accuracy (overall accuracy = 93.94%, Kappa coefficient = 0.92) than the other classifications without such a feature optimizing process. Our results indicate the potential of the ASTFS method for improving regional crop mapping.

Unprecedented human-induced land cover changes happened in China after the Reform and Opening-up in 1978, matching with the era of Landsat satellite series. However, it is still unknown whether Landsat data can effectively support retrospective analysis of land cover changes in China over the past four decades. Here, for the first time, we conduct a systematic investigation on the availability of Landsat data in China, targeting its application for retrospective and continuous monitoring of land cover changes. The latter is significant to assess impact of land cover changes, and consequences of past land policy and management interventions. The total and valid observations (excluding clouds, cloud shadows, and terrain shadows) from Landsat 5/7/8 from 1984 to 2017 were quantified at pixel scale, based on the cloud computing platform Google Earth Engine (GEE). The results show higher intensity of Landsat observation in the northern part of China as compared to the southern part. The study provides an overall picture of Landsat observations suitable for satellite-based annual land cover monitoring over the entire country. We uncover that two sub-regions of China (i.e., Northeast China-Inner Mongolia-Northwest China, and North China Plain) have sufficient valid observations for retrospective analysis of land cover over 30 years (1987–2017) at an annual interval; whereas the Middle-Lower Yangtze Plain (MLYP) and Xinjiang (XJ) have sufficient observations for annual analyses for the periods 1989–2017 and 2004–2017, respectively. Retrospective analysis of land cover is possible only at a two-year time interval in South China (SC) for the years 1988–2017, Xinjiang (XJ) for the period 1992–2003, and the Tibetan Plateau (TP) during 2004–2017. For the latter geographic regions, land cover dynamics can be analyzed only at a three-year interval prior to 2004. Our retrospective analysis suggest that Landsat-based analysis of land cover dynamics at an annual interval for the whole country is not feasible; instead, national monitoring at two- or three-year intervals could be achievable. This study provides a preliminary assessment of data availability, targeting future continuous land cover monitoring in China; and the code is released to the public to facilitate similar data inventory in other regions of the world.

China launched the Three-North Shelterbelt Forest Program (TNSFP) in 1978 in northern China to combat desertification and dust storms, but it is still controversial in ecologically fragile arid and semi-arid areas, which is partly due to the uncertainties of monitoring of the spatial-temporal changes of forest distribution. In this study, we aim to provide an overall retrospect of the forest changes (i.e., forest gain and forest loss) in northern China during 2007–2017, and to analyze the forest changes in different precipitation zones. We first generated annual forest maps at 30 m spatial resolution during 2007–2017 in northern China through integrating Landsat and PALSAR/PALSAR-2 data. We found the PALSAR/Landsat-based forest maps outperform other four existing products (i.e., JAXA F/NF, FROM-GLC, GlobeLand30, and NLCD-China) from either PALSAR or Landsat data, with a higher overall accuracy 96% ± 1%. The spatial-temporal analyses of forests showed a substantial forest expansion from 316,898 ± 34,537 km² in 2007 to 384,568 ± 35,855 km² in 2017 in the central and eastern areas. We found a higher forest loss rate (i.e., 35%) in the precipitation zones with the annual mean precipitation less than 400 mm (i.e., the arid and semi-arid areas) comparing to that (i.e., 25%) in the zones with more than 400 mm (i.e., the humid areas), which suggests the lower surviving rate in the drylands. This study provides satellite-based evidence for the forest changes in different precipitation zones, and suggests that the likely impacts of precipitation on afforestation effectiveness should be considered in future implementation of ecological restoration projects like TNSFP.

Information about urban land use is important for urban planning and sustainable development. The emergence of geospatial big data (GBD), increased the availability of remotely sensed (RS) data and the development of new methods for data integration to provide new opportunities for mapping types of urban land use. However, the modes of RS and GBD integration are diverse due to the differences in data, study areas, classifiers, etc. In this context, this study aims to summarize the main methods of data integration and evaluate them via a case study of urban land use mapping in Hangzhou, China. We first categorized the RS and GBD integration methods into decision-level integration (DI) and feature-level integration (FI) and analyzed their main differences by reviewing the existing literature. The two methods were then applied for mapping urban land use types in Hangzhou city, based on urban parcels derived from the OpenStreetMap (OSM) road network, 10 m Sentinel-2A images, and points of interest (POI). The corresponding classification results were validated quantitatively and qualitatively using the same testing dataset. Finally, we illustrated the advantages and disadvantages of both approaches via bibliographic evidence and quantitative analysis. The results showed that: (1) The visual comparison indicates a generally better performance of DI-based classification than FI-based classification; (2) DI-based urban land use mapping is easy to implement, while FI-based land use mapping enables the mixture of features; (3) DI-based and FI-based methods can be used together to improve urban land use mapping, as they have different performances when classifying different types of land use. This study provides an improved understanding of urban land use mapping in terms of the RS and GBD integration strategy.

The phenology-based approach has proven effective for paddy rice mapping due to the unique flooding and transplanting features of rice during the early growing season. However, the method may be greatly affected if no valid observations are available during the flooding and rice transplanting phase. Here, we compare the effects of data availability of different sensors in the critical phenology phase, thereby supporting paddy rice mapping based on phenology-based approaches. Importantly, our study further analyzed the effects of the spatial pattern of the valid observations related to certain factors (i.e., sideslips, clouds, and temporal window lengths of flooding and rice transplanting), which supply the applicable area of the phenology-based approach indications. We first determined the flooding and rice transplanting phase using in situ observational data from agrometeorological stations and remote sensing data, then evaluated the effects of data availability in this phase of 2020 in China using all Landsat-7 and 8 and Sentinel-2 data. The results show that on the country level, the number of average valid observations during the flooding and rice transplanting phase was more than ten for the integration of Landsat and Sentinel images. On the sub-country level, the number of average valid observations was high in the cold temperate zone (17.4 observations), while it was relatively lower in southern China (6.4 observations), especially in Yunnan–Guizhou Plateau, which only had three valid observations on average. Based on the multicollinearity test, the three factors are significantly correlated with the absence of valid observations: (R² = 0.481) and Std.Coef. (Std. Err.) are 0.306 (0.094), −0.453 (0.003) and −0.547 (0.019), respectively. Overall, these results highlight the substantial spatial heterogeneity of valid observations in China, confirming the reliability of the integration of Landsat-7 and 8 and Sentinel-2 imagery for paddy rice mapping based on phenology-based approaches. This can pave the way for a national-scale effort of rice mapping in China while further indicating potential omission errors in certain cloud-prone regions without sufficient optical observation data, i.e., the Sichuan Basin.

The Tibetan Plateau (TP), known as “The Roof of World”, has expansive alpine grasslands and is a hotspot for climate change studies. However, cropland expansion and increasing anthropogenic activities have been poorly documented, let alone the effects of agricultural activities on food security and environmental change in the TP. The existing cropland mapping products do not depict the spatiotemporal characteristics of the TP due to low accuracies and inconsistent cropland distribution, which is affected by complicated topography and impedes our understanding of cropland expansion and its associated environmental impacts. One of the biggest challenges of cropland mapping in the TP is the diverse crop phenology across a wide range of elevations. To decrease the classification errors due to elevational differences in crop phenology, we developed two pixel- and phenology-based algorithms to map croplands using Landsat imagery and the Google Earth Engine platform along the Brahmaputra River and its two tributaries (BRTT) in the Tibet Autonomous Region, also known as the granary of TP, in 2015–2019. Our first phenology-based cropland mapping algorithm (PCM1) used different thresholds of land surface water index (LSWI) by considering varied crop phenology along different elevations. The second algorithm (PCM2) further offsets the phenological discrepancy along elevational gradients by considering the length and peak of the growing season. We found that PCM2 had a higher accuracy with fewer images compared with PCM1. The number of images for PCM2 was 279 less than PCM1, and the Matthews correlation coefficient for PCM2 was 0.036 higher than PCM1. We also found that the cropland area in BRTT was estimated to be 1979 ± 52 km² in the late 2010s. Croplands were mainly distributed in the BRTT basins with elevations of 3800–4000 m asl. Our phenology-based methods were effective for mapping croplands in mountainous areas. The spatially explicit information on cropland area and distribution in the TP aid future research into the effects of cropland expansion on food security and environmental change in the TP.