Landsat-8 and Sentinel-2 Based Prediction of Forest Plantation C Stock Using Spatially Varying Coefficient Bayesian Hierarchical Models

Abstract

This study sought to establish the performance of Spatially Varying Coefficient (SVC) Bayesian Hierarchical models using Landsat-8, and Sentinel-2 derived auxiliary information in predicting plantation forest carbon (C) stock in the eastern highlands of Zimbabwe. The development and implementation of Zimbabwe’s land reform program undertaken in the year 2000 and the subsequent redistribution and resizing of large-scale land holdings are hypothesized to have created heterogeneity in aboveground forest biomass in plantation ecosystems. The Bayesian hierarchical framework, accommodating residual spatial dependence and non-stationarity of model predictors, was evaluated. Firstly, SVC models utilizing Normalized Difference Vegetation Index (NDVI), Soil Adjusted Vegetation Index (SAVI), and Enhanced Vegetation Index (EVI), derived from Landsat-8 and Sentinel-2 data and 191 sampled C stock observations were constructed. The SVC models built for each of the two multispectral remote sensing data sets were assessed based on the goodness of fit criterion as well as the predictive performance using a 10-fold cross-validation technique. The introduction of spatial random effects in the form of Landsat-8 and Sentinel-2 derived covariates to the model intercept improved the model fit and predictive performance where residual spatial dependence was dominant. For the Landsat-8 C stock predictive model, the RMSPE for the non-spatial, Spatially Varying Intercept (SVI) and SVC models were 8 MgCha−1, 7.77 MgCha−1, and 6.42 MgCha−1 whilst it was 7.85 MgCha−1, 7.69 MgCha−1 and 6.23 MgCha−1 for the Sentinel-2 C stock predictive models, respectively. Overall, the Sentinel-2-based SVC model was preferred for predicting C stock in plantation forest ecosystems as its model provided marginally tighter credible intervals, [1.17–1.60] MgCha−1 when compared to the Landsat-8 based SVC model with 95% credible intervals of [1.13–1.62] Mg Cha−1. The built SVC models provided an understanding of the performance of the multispectral remote sensing derived predictors for modeling C stock and thus provided an essential foundation required for updating the current carbon forest plantation databases.