Mapping the maximum peat thickness of cultivated organic soils in the southwest plain of Montreal

Large organic deposits in the southwestern plain of Montreal have been converted to agricultural land for vegetable pro-duction. In addition to the variable depth of the organic deposits, these soils commonly have an impermeable coprogenous layer between the peat and the underlying mineral substratum. Estimations of the depth and thickness of these materials are critical for soil management. Therefore, five drained and cultivated peatlands were studied to estimate their maximum peat thickness (MPT)--a potential key soil property that can help identify management zones for their conservation. MPT can be defined as the depth to the mineral layer (DML) minus the coprogenous layer thickness (CLT). The objective of this study was to estimate DML, CLT, and MPT at a regional scale using environmental covariates derived from remote sensing. Three machine-learning models (Cubist, Random Forest, and k-Nearest Neighbor) were compared to produce maps of DML and CLT, which were combined to generate MPT at a spatial resolution of 10 m. The Cubist model performed the best for predicting both features of interest, yielding Lin's concordance correlation coefficients of 0.43 and 0.07 for DML and CLT, respectively, using a spatial cross-validation procedure. Interpretation of the drivers of CLT was limited by the poor predictive power of the final model. More precise data on MPT are needed to support soil conservation practices, and more CLT field observations are required to obtain a higher prediction accuracy. Nonetheless, digital soil mapping using open-access geospatial data shows promise for understanding and managing cultivated peatlands.

Automated summary

This study used remote sensing data and machine learning models to estimate the soil thickness of cultivated peatlands in the southwestern plain of Montreal. The results showed that the Cubist model performed the best in predicting the depth and thickness. However, the interpretation of the drivers for coprogenous layer thickness was limited by the predictive power of the model.