Random forest regression kriging modeling for soil organic carbon density estimation using multi-source environmental data in central Vietnamese forests

• Published in: Modeling earth systems and environment 2024-09, Vol.10 (6), p.7137-7158
• Main Authors: Ho, Viet Hoang, Morita, Hidenori, Bachofer, Felix, Ho, Thanh Ha
• Format: Article
• Language: English
• Subjects: Autocorrelation; Carbon; Carbon cycle; Chemistry and Earth Sciences; Climate change; Climate change mitigation; Climatic data; Computer Science; Data sources; Density; Digital Elevation Models; Digital imaging; Earth and Environmental Science; Earth Sciences; Earth System Sciences; Ecosystems; Environment; Environmental management; Forest ecosystems; Forest management; Forest soils; Imagery; Math. Appl. in Environmental Science; Mathematical Applications in the Physical Sciences; Modelling; Organic carbon; Organic soils; Original Article; Phased arrays; Physics; Predictions; Radar arrays; Root-mean-square errors; SAR (radar); Satellite imagery; Satellite observation; Soil; Soil analysis; Spatial data; Spatial distribution; Statistics for Engineering; Sustainability management; Sustainable forestry; Synthetic aperture radar; Terrestrial ecosystems
• ISSN: 2363-6203; 2363-6211
• DOI: 10.1007/s40808-024-02158-1

Forest soil organic carbon plays a vital role in the terrestrial carbon cycle. Accurately analyzing the spatial distribution of soil organic carbon density (SOCD) is therefore necessary for sustainable forest management and climate change mitigation. Previous studies explored the potential of random forest (RF) in modeling forest SOCD using various environmental data sources. However, how forest SOCD prediction would be affected by using random forest regression kriging (RFRK), which integrates the predictive power of RF in generating deterministic trends and the capability of the ordinary kriging (OK) in handling spatial autocorrelation structure of residuals, based on the environmental data sources and their combinations remains elusive and deserves further exploration. For this purpose, 104 soil samples were collected at a depth of 30 cm in forest ecosystems of Central Vietnam, and 33 environmental covariates were derived from Sentinel-2 (S2) imagery, Advanced Land Observing Satellite-2 Phased Array L-band Synthetic Aperture Radar-2 (AL2) imagery, digital elevation model (DEM), and climatic data. Using a leave-one-out cross-validation procedure to evaluate and compare the model performances, four metrics, including coefficient of determination (R 2 ), mean absolute error (MAE), root mean square error (RMSE), and relative improvement (RI), were calculated. The results showed that enhanced RFRK performance for forest SOCD estimation was found with the inclusion of additional environmental data sources, with RFRK based on all data sources achieving a high accuracy (R 2  = 0.78, MAE = 8.28 t·ha −1 , and RMSE = 10.54 t·ha −1 ). The comparison of the RF and RFRK models exhibited that additionally interpolated residuals by OK were more accurate than only considering the influences of predictor covariates. The relative improvement of the RFRK models over the RF models in forest SOCD estimation was notable, with RI R 2 ranging from 8.20 to 65.00%, RI MAE ranging from 8.18 to 21.07%, and RI RMSE ranging from 6.76 to 18.18%. The result from our case study emphasizes the robustness of RFRK using S2, AL2, DEM, and climatic data in accurately predicting forest SOCD.