Abstract. The construction of terraces is a key soil conservation
practice on agricultural land in China providing multiple valuable
ecosystem services. Accurate spatial information on terraces is needed for
both management and research. In this study, the first 30 m resolution
terracing map of the entire territory of China is produced by a supervised
pixel-based classification using multisource and multi-temporal data based
on the Google Earth Engine (GEE) platform. We extracted time-series spectral
features and topographic features from Landsat 8 images and the Shuttle
Radar Topography Mission digital elevation model (SRTM DEM) data,
classifying cropland area (cultivated land of Globeland30) into terraced and
non-terraced types through a random forest classifier. The overall accuracy
and kappa coefficient were evaluated by 10 875 test samples and achieved
values of 94 % and 0.72, respectively. For terrace class, the producer's
accuracy (PA) was 79.945 %, and the user's accuracy (UA) was 71.149 %.
The classification performed best in the Loess Plateau and southwestern
China, where terraces are most numerous. Some northeastern, eastern-central,
and southern areas had relatively high uncertainty. Typical errors in the
mapping results are from the sloping cropland (non-terrace cropland with a slope
of ≥ 5∘), low-slope terraces, and non-crop vegetation.
Terraces are widely distributed in China, and the total terraced area was
estimated to be 53.55 Mha (i.e., 26.43 % of China's cropland area) by pixel
counting (PC) method and 58.46 ± 2.99 Mha (i.e., 28.85 % ± 1.48 % of China's cropland area) by error-matrix-based model-assisted
estimation (EM) method. Elevation and slope were identified as the main
features in the terrace/non-terrace classification, and multi-temporal
spectral features (such as percentiles of NDVI, TIRS2, and BSI) were also
essential. Terraces are more challenging to identify than other land use
types because of the intra-class feature heterogeneity, interclass feature
similarity, and fragmented patches, which should be the focus of future
research. Our terrace mapping algorithm can be used to map large-scale
terraces in other regions globally, and our terrace map will serve as a
landmark for studies on multiple ecosystem service assessments including
erosion control, carbon sequestration, and biodiversity conservation. The
China terrace map is available to the public at
https://doi.org/10.5281/zenodo.3895585 (Cao et al., 2020).