The global distribution of cropping intensity (CI) is essential to our understanding of agricultural land use management on Earth. Optical remote sensing has revolutionized our ability to map CI over large areas in a repeated and cost-efficient manner. Previous studies have mainly focused on investigating the spatiotemporal patterns of CI ranging from regions to the entire globe with the use of coarse-resolution data, which are inadequate for characterizing farming practices within heterogeneous landscapes. To fill this knowledge gap, in this study, we utilized multiple satellite data to develop a global, spatially continuous CI map dataset at 30 m resolution (GCI30). Accuracy assessments indicated that GCI30 exhibited high agreement with visually interpreted validation samples and in situ observations from the PhenoCam network. We carried out both statistical and spatial comparisons of GCI30 with six existing global CI estimates. Based on GCI30, we estimated that the global average annual CI during 2016–2018 was 1.05, which is close to the mean (1.09) and median (1.07) CI values of the existing six global CI estimates, although the spatial resolution and temporal coverage vary significantly among products. A spatial comparison with two satellite-based land surface phenology products further suggested that GCI30 was not only capable of capturing the overall pattern of global CI but also provided many spatial details. GCI30 indicated that single cropping was the primary agricultural system on Earth, accounting for 81.57 % (12.28×10⁶ km²) of the world's cropland extent. Multiple-cropping systems, on the other hand, were commonly observed in South America and Asia. We found large variations across countries and agroecological zones, reflecting the joint control of natural and anthropogenic drivers on regulating cropping practices. As the first global-coverage, fine-resolution CI product, GCI30 is expected to fill the data gap for promoting sustainable agriculture by depicting worldwide diversity of agricultural land use intensity. The GCI30 dataset is available on Harvard Dataverse: https://doi.org/10.7910/DVN/86M4PO (Zhang et al., 2020).

Fusing multiple change detection results has great potentials in dealing with the spectral variability in multitemporal very high-resolution (VHR) remote sensing images. However, it is difficult to solve the problem of uncertainty, which mainly includes the inaccuracy of each candidate change map and the conflicts between different results. Dempster–Shafer theory (D–S) is an effective method to model uncertainties and combine multiple evidences. Therefore, in this paper, we proposed an urban change detection method for VHR images by fusing multiple change detection methods with D–S evidence theory. Change vector analysis (CVA), iteratively reweighted multivariate alteration detection (IRMAD), and iterative slow feature analysis (ISFA) were utilized to obtain the candidate change maps. The final change detection result is generated by fusing the three evidences with D–S evidence theory and a segmentation object map. The experiment indicates that the proposed method can obtain the best performance in detection rate, false alarm rate, and comprehensive indicators.

The timely and reliable estimation of imperviousness is essential for the scientific understanding of human-Earth interactions. Due to the unique capacity of capturing artificial light luminosity and long-term data records, the Defense Meteorological Satellite Program (DMSP)’s Operational Line-scan System (OLS) nighttime light (NTL) imagery offers an appealing opportunity for continuously characterizing impervious surface area (ISA) at regional and continental scales. Although different levels of success have been achieved, critical challenges still remain in the literature. ISA results generated by DMSP-OLS NTL alone suffer from limitations due to systemic defects of the sensor. Moreover, the majority of developed methodologies seldom consider spatial heterogeneity, which is a key issue in coarse imagery applications. In this study, we proposed a novel method for multi-temporal ISA estimation. This method is based on a linear regression model developed between the sub-pixel ISA fraction and a multi-source index with the integrated use of DMSP-OLS NTL and MODIS NDVI. In contrast with traditional regression analysis, we incorporated spatial information to the regression model for obtaining spatially adaptive coefficients at the per-pixel level. To produce multi-temporal ISA maps using a mono-temporal reference dataset, temporally stable samples were extracted for model training and validation. We tested the proposed method in the Yangtze River Delta and generated annual ISA fraction maps for the decade 2000–2009. According to our assessments, the proposed method exhibited substantial improvements compared with the standard linear regression model and provided a feasible way to monitor large-scale impervious surface dynamics.

It is of great significance to understand the extent and distribution of bamboo for its valuable ecological services and economic benefits. However, it is challenging to map bamboo using remote sensing images over a large area because of the similarity between bamboo and other vegetation types, the availability of clear optical images, huge workload of image processing, and sample collection. In this study, we use the Landsat 8 times series images archive to map bamboo forests in China via the Google Earth engine. Several spectral indices were calculated and used as classification features, including the normalized difference vegetation index (NDVI), the normalized difference moisture index (NDMI) and textural features of the gray-level co-occurrence matrix (GLCM). We found that the bamboo forest covered an area of 709.92 × 10⁴ hectares, with the provinces of Fujian, Jiangxi, and Zhejiang containing the largest area concentrations. The bamboo forest map was accurate and reliable with an average producer’s accuracy of 89.97%, user’s accuracy of 78.45% and kappa coefficient of 0.7789. In addition, bamboo was mainly distributed in forests with an elevation of 300–1200 m above sea level, average annual precipitation of 1200–1500 mm and average day land surface temperature of 19–25 °C. The NDMI is particularly useful in differentiating bamboo from other vegetation because of the clear difference in canopy moisture content, whilst NDVI and elevation are also helpful to improve the bamboo classification accuracy. The bamboo forest map will be helpful for bamboo forest industry planning and could be used for evaluating the ecological service of the bamboo forest.

Impervious surfaces are commonly acknowledged as major components of human settlements. The expansion of impervious surfaces could lead to a series of human−dominated environmental and ecological issues. Tracing impervious surface dynamics at a finer temporal−spatial scale is a critical way to better understand the increasingly human-dominated system of Earth. In this study, we put forward a new scheme to conduct long-term monitoring of impervious−relevant land disturbances using high frequency Landsat archives and the Google Earth Engine (GEE). First, the developed region was identified using a classification-based approach. Then, the GEE-version LandTrendr (Landsat-based detection of Trends in Disturbance and Recovery) was used to detect land disturbances, characterizing the conversion from vegetation to impervious surfaces. Finally, the actual disturbance areas within the developed regions were derived and quantitatively evaluated. A case study was conducted to detect impervious surface dynamics in Nanjing, China, from 1988 to 2018. Results show that our scheme can efficiently monitor impervious surface dynamics at yearly intervals with good accuracy. The overall accuracy (OA) of the classification results for 1988 and 2018 are 95.86% and 94.14%. Based on temporal−spatial accuracy assessments of the final detection result, the temporal accuracy is 90.75%, and the average detection time deviation is −1.28 a. The OA, precision, and recall of the sampling inspection, respectively, are 84.34%, 85.43%, and 96.37%. This scheme provides new insights into capturing the expansion of impervious−relevant land disturbances with high frequency Landsat archives in an efficient way.