Research Article

Analysis of Factors Influencing Spatial Distribution of Soil Erosion under Diverse Subwatershed Based on Geospatial Perspective: A Case Study at Citarum Watershed, West Java, Indonesia

Table 1

List of required data inputs for the InVEST-SDR and Geodetector models.

Data descriptionSourceProcessingSpatial resolution

(a) InVEST-SDR model
Digital elevation modelShuttle Radar Topography Mission (SRTM) Global (USGS)30 m × 30 m
Land use/land cover (LULC) dataLandsat/8 OLI imageries with acquisition date in 2000 (https:/www.usgs.gov/ and Google Earth Engine)Supervised classification30 m × 30 m
Precipitation dataMeteorological, Climatological, and Geophysical AgencyNumeric or tabular data, along with location coordinates, are processed using a spline interpolation method30 m × 30 m
Isoerosivity map (R factor)
Soil dataIndonesian Soil Research Institute (the soil map type of 1 : 100,000)Vector file30 m × 30 m
Soil erodibility map (K factor)
Watershed boundary dataCitarum-Ciliwung Watershed Management CenterVector file
Data biophysical table (PUSLE and CUSLE)FAO.org(.csv)30 m × 30 m

(b) Geodetector model
(1) Physical data
(1.1) Topographic factors
Slope (variable name: XA1)Shuttle Radar Topography Mission (SRTM) Global (USGS)Calculation from DEM (X2)30 m × 30 m
Digital elevation model (DEM) (variable name: XA2)Shuttle Radar Topography Mission (SRTM) Global (USGS)30 m × 30 m
(1.2) Climate factor
Annual average temperature (TEM) (variable name: XA3)Meteorological, Climatological, and Geophysical AgencyNumeric or tabular data, along with location coordinates, are processed using a spline interpolation method30 m × 30 m
Annual average precipitation, PRE (variable name: XA4)Meteorological, Climatological, and Geophysical AgencyNumeric or tabular data, along with location coordinates, are processed using a spline interpolation method30 m × 30 m
(1.3) Vegetation factor
Net primary production (NPP) (variable name: XA5)MODIS-MOD17A3HGF V6.1 product in 2000, 2010, and 2020 (https:/www.usgs.gov/ and Google Earth Engine)The sum of all 8-day net photosynthesis is the difference between the gross primary productivity and maintenance respiration. Resampling from 500 m × 500 m30 m × 30 m
Fractional vegetation cover (FVC) (variable name: XA6)MODIS-MOD13Q1 with from NDVI value in 2000, 2010, and 2020 (https:/www.usgs.gov/ and Google Earth Engine)FVC = ((NDVI − 0.2)/0.3) 100. Resampling from 250 m × 250 m30 m × 30 m
(2) Human activities
Income per capita (variable name: XA7)West Java Central Statistics Agency 2010 and 2020Numerical/table data30 m × 30 m
Population density (variable name: XA8)West Java Central Statistics Agency 2010 and 2020Numerical/table data30 m × 30 m

To get more detailed erosion class data, the watershed boundaries with subdistricts are used.