Soil Organic Carbon Mapping from Remote Sensing: The Effect of Crop Residues

• Published in: Remote sensing (Basel, Switzerland) Vol. 12; no. 12; p. 1913
• Main Authors: Dvorakova, Klara, Shi, Pu, Limbourg, Quentin, van Wesemael, Bas
• Format: Journal Article
• Language: English
• Published: Basel MDPI AG 01-06-2020
• Subjects: Agricultural land; Carbon; Cellulose; cellulose absorption index; Corn; Crop residues; Environmental monitoring; hyperspectral data; Hyperspectral imaging; Landsat; Landsat satellites; Lignin; Lignocellulose; Mapping; Mathematical analysis; Normalized difference vegetative index; Organic carbon; Organic soils; partial least squares regression; Pixels; Prediction models; Regression models; Remote sensing; Residues; Root-mean-square errors; Satellite imagery; Satellites; Short wave radiation; Soil conditions; Soil moisture; soil organic carbon mapping; Soil properties; Soil surfaces; Soils; Spectrum analysis; Thresholds; Vegetation; Vegetation index
• ISSN: 2072-4292; 2072-4292
• DOI: 10.3390/rs12121913

Since the onset of agriculture, soils have lost their organic carbon to such an extent that the soil functions of many croplands are threatened. Hence, there is a strong demand for mapping and monitoring critical soil properties and in particular soil organic carbon (SOC). Pilot studies have demonstrated the potential for remote sensing techniques for SOC mapping in croplands. It has, however, been shown that the assessment of SOC may be hampered by the condition of the soil surface. While growing vegetation can be readily detected by means of the well-known Normalized Difference Vegetation Index (NDVI), the distinction between bare soil and crop residues is expressed in the shortwave infrared region (SWIR), which is only covered by two broad bands in Landsat or Sentinel-2 imagery. Here we tested the effect of thresholds for the Cellulose Absorption Index (CAI), on the performance of SOC prediction models for cropland soils. Airborne Prism Experiment (APEX) hyperspectral images covering an area of 240 km2 in the Belgian Loam Belt were used together with a local soil dataset. We used the partial least square regression (PLSR) model to estimate the SOC content based on 104 georeferenced calibration samples (NDVI < 0.26), firstly without setting a CAI threshold, and obtained a satisfactory result (coefficient of determination (R2) = 0.49, Ratio of Performance to Deviation (RPD) = 1.4 and Root Mean Square Error (RMSE) = 2.13 g kgC−1 for cross-validation). However, a cross comparison of the estimated SOC values to grid-based measurements of SOC content within three fields revealed a systematic overestimation for fields with high residue cover. We then tested different CAI thresholds in order to mask pixels with high residue cover. The best model was obtained for a CAI threshold of 0.75 (R2 = 0.59, RPD = 1.5 and RMSE = 1.75 g kgC−1 for cross-validation). These results reveal that the purity of the pixels needs to be assessed aforehand in order to produce reliable SOC maps. The Normalized Burn Ratio (NBR2) index based on the SWIR bands of the MSI Sentinel 2 sensor extracted from images collected nine days before the APEX flight campaign correlates well with the CAI index of the APEX imagery. However, the NBR2 index calculated from Sentinel 2 images under moist conditions is poorly correlated with residue cover. This can be explained by the sensitivity of the NBR2 index to both soil moisture and residues.