Integrating Environmental Variables into Geostatistical Interpolation: Enhancing Soil Mapping for the MEDALUS Model in Montenegro

Stefan Miletić (), Jelena Beloica and Predrag Miljković
Additional contact information
Stefan Miletić: Department of Ecological Engineering for Soil and Water Resources Protection, Faculty of Forestry, University of Belgrade, 11000 Belgrade, Serbia
Jelena Beloica: Department of Ecological Engineering for Soil and Water Resources Protection, Faculty of Forestry, University of Belgrade, 11000 Belgrade, Serbia
Predrag Miljković: Department of Ecological Engineering for Soil and Water Resources Protection, Faculty of Forestry, University of Belgrade, 11000 Belgrade, Serbia

Land, 2025, vol. 14, issue 4, 1-23

Abstract: Geostatistical methods are important in analyzing natural resources providing input data for complex mathematical models that address environmental processes and their spatial distribution. Ten interpolation methods and one empirical-based classification grounded in empirical knowledge, with a total of 929 soil samples, were used to create the most accurate spatial prediction maps for clay, sand, humus, and soil depth in Montenegro. These analyses serve as a preparatory phase and prioritize the practical application of the obtained results for the implementation and improvement of the MEDALUS model. This model, used to assess sensitivity to land degradation, effectively integrates into broader current and future research. The study emphasizes the importance of incorporating auxiliary variables, such as topography, climate, and vegetation data, enhancing explanatory power and accuracy in delineating the environmental characteristics, ensuring better adaptability to the studied area. The results were validated by the coefficient of determination (R 2 ) and root mean square error (RMSE). For the clay, EBKRP (empirical Bayesian kriging regression prediction) achieved R 2 = 0.35 and RMSE = 6.95%, for the sand, it achieved R 2 = 0.34 and RMSE = 17.38%, for the humus, it achieved R 2 = 0.50 and RMSE = 3.80%, and for the soil depth, it achieved R 2 = 0.76 and RMSE = 5.36 cm. These results indicate that EBKRP is the optimal method for accurately mapping soil characteristics in future research in Montenegro.

Keywords: soil mapping; environmental covariates; GIS; auxiliary variables; Montenegro (search for similar items in EconPapers)
JEL-codes: Q15 Q2 Q24 Q28 Q5 R14 R52 (search for similar items in EconPapers)
Date: 2025
References: View references in EconPapers View complete reference list from CitEc
Citations:

Downloads: (external link)
https://www.mdpi.com/2073-445X/14/4/702/pdf (application/pdf)
https://www.mdpi.com/2073-445X/14/4/702/ (text/html)

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:gam:jlands:v:14:y:2025:i:4:p:702-:d:1620579

Access Statistics for this article

Land is currently edited by Ms. Carol Ma

More articles in Land from MDPI
Bibliographic data for series maintained by MDPI Indexing Manager ().