Creating accurate digital maps of the agrochemical properties of soils on a field scale with a limited data set is a problem that slows down the introduction of precision farming. The use of machine learning methods based on the use of direct and indirect predictors of spatial changes in the agrochemical properties of soils is promising. Spectral indicators of open soil based on remote sensing data, as well as soil properties, were used to create digital maps of available forms of nitrogen, phosphorus, and potassium. It was shown that machine learning methods based on support vectors (SVMr) and random forest (RF) using spectral reflectance data are similarly accurate at spatial prediction. An acceptable prediction was obtained for available nitrogen and available potassium; the variability of available phosphorus was modeled less accurately. The coefficient of determination (R2) of the best model for nitrogen is R2SVMr = 0.90 (Landsat 8 OLI) and R2SVMr = 0.79 (Sentinel 2), for potassium—R2SVMr = 0.82 (Landsat 8 OLI) and R2SVMr = 0.77 (Sentinel 2), for phosphorus—R2SVMr = 0.68 (Landsat 8 OLI), R2SVMr = 0.64 (Sentinel 2). The models based on remote sensing data were refined when soil organic matter (SOC) and fractions of texture (Silt, Clay) were included as predictors. The SVMr models were the most accurate. For Landsat 8 OLI, the SVMr model has a R2 value: nitrogen—R2 = 0.95, potassium—R2 = 0.89 and phosphorus—R2 = 0.65. Based on Sentinel 2, nitrogen—R2 = 0.92, potassium—R2 = 0.88, phosphorus—R2 = 0.72. The spatial prediction of nitrogen content is influenced by SOC, potassium—by SOC and texture, phosphorus—by texture. The validation of the final models was carried out on an independent sample on soils from a chernozem zone. For nitrogen based on Landsat 8 OLI R2 = 0.88, for potassium R2 = 0.65, and for phosphorus R2 = 0.31. Based on Sentinel 2, for nitrogen R2 = 0.85, for potassium R2 = 0.62, and for phosphorus R2 = 0.71. The inclusion of SOC and texture in remote sensing-based machine learning models makes it possible to improve the spatial prediction of nitrogen, phosphorus and potassium availability of soils in chernozem zones and can potentially be widely used to create digital agrochemical maps on the scale of a single field.
Assessing the performance of machine learning algorithms for soil salinity mapping in Google Earth Engine Platform using Sentinel-2A and Landsat-8 OLI data
