Using spatial aggregation of soil multifunctionality maps to support uncertainty-aware planning decisions

To ensure soil preservation, it is essential to incorporate the soil's ability to provide ecosystem services into the spatial planning process. For well-informed planning decisions, stakeholders need spatially explicit information on the state of the soils and the functions they fulfil, with sufficient spatial resolution and quantified uncertainty. It has been shown that Digital Soil Mapping (DSM) products can provide such information. However, in some cases, fine spatial resolution coupled with high levels of uncertainty may lead stakeholders to overlook the inherent uncertainties in the information. Spatial aggregation of DSM products opens up a promising avenue for obtaining maps that are more tailored to the users' scales of decision making while facilitating uncertainty communication. In this perspective, we propose a new spatial aggregation approach relying on spatially constrained agglomerative clustering (AC). The spatial aggregation approach is applied to a 25-m-resolution soil potential multifunctionality index (SPMI) map developed for the coastal plain of the Occitanie Region. This DSM product was increasingly aggregated to obtain SPMI maps of different resolutions displaying two distinct areal metrics: proportions of area above a given threshold of SPMI, and mean SPMI. Each map was evaluated through a set of indicators selected for their potential impact on user decision making: mean spatial resolution, overall predicted uncertainty, quantity of information and mean within-unit variability. The maps were compared with respect to these indicators to other maps obtained with alternative aggregation methods employed in DSM literature (maps aggregated according to some administrative units and QuadMaps). We show that all the tested aggregation methods produced a substantial decrease of the map uncertainty with moderate loss of spatial resolution. However, only AC preserved the fine spatial pattern of the initial DSM product while enabling fine tuning of the uncertainty displayed to end-users. We show that AC can simplify the identification of extensive regions characterized by low uncertainty without losing information regarding soil multifunctionality, thereby facilitating and enhancing the efficiency of planning decisions.