Evaluation of soil fertility using combination of Landsat 8 and Sentinel‑2 data in agricultural lands

Today, remote sensing is widely used to estimate soil properties. Because it is an easy and accessible way to estimate soil properties that are difficult to estimate in the field. Based on this, to evaluate the soil fertility (SF), soil sampling was performed irregularly from the surface depth of 0–30 cm in 216 points, 11 soil properties were measured, and the soil fertility index (SFI) was calculated by soil properties. Simultaneously, we combined satellite images of Landsat 8 and Sentinel-2 using the Gram-Schmidt algorithm. Finally, multiple linear regression SFI was calculated using satellite data, as well as the spatial distribution of SFI was obtained in very low, low, moderate, high, and very high classes. Our findings showed that the combination of Landsat 8 and Sentinel-2 data using the Gram-Schmidt algorithm has a higher correlation with SFI than when these data are individually. Therefore, combined Landsat 8 and Sentinel 2 data were used for SFI modeling. Using model selection procedure indices (including Cp, AIC, and ρc criteria), the visible range bands, notably blue (r = 0.65), green (r = 0.63), and red (r = 0.61), provide the best model for estimating SFI (R2 = 0.43, Cp = 3.34, AIC = -277.4, and ρc = 0.44). Therefore, these bands were used to estimate the SFI index. Also, the spatial distribution of the SIF index showed that the most significant area was related to the low class, and the lowest area belonged to the high and very high fertility classes. According to these results, it can be concluded that using the combination of Landsat 8 and Sentinel 2 bands to estimate soil fertility index in agricultural lands can increase the accuracy of soil fertility estimation.