Mapping and uncertainty of predictions based on multiple primary variables from joint co-simulation with Landsat TM image and polynomial regression

In the management of natural resources, multiple variables correlated with each other usually need to be mapped jointly. However, joint mapping and spatial uncertainty analyses are very difficult mainly because of interactions among variables and imperfection of existing methods. There is abundant evidence that considering interactions among variables and spatial information from neighbors can result in improved maps. This study presents a remote sensing-aided method for that purpose. The method is based on the integration of joint sequential co-simulation with Landsat TM image for mapping and polynomial regression for spatial uncertainty analysis. The method was applied to a case study in which ground cover (GC), canopy cover (CC), and vegetation height (VH) were jointly mapped to derive a map of the vegetation cover factor for predicting soil loss. The variance contributions from the variables, their interactions, and the spatial information from neighbors leading to uncertainty of predicted vegetation cover factor were assessed. The results showed that in addition to unbiased maps, this method reproduced the spatial variability of the variables and the spatial correlation among them, and successfully quantified the effect of variation from all the components on the prediction of the vegetation cover factor.