You're using an outdated browser. Please upgrade to a modern browser for the best experience.
Modern Fertilization in Precision Agriculture: Comparison
Please note this is a comparison between Version 2 by Jason Zhu and Version 3 by Jason Zhu.

Over the past decade, the development of new technologies and sensors and sensor minimization has led to the development and improvement of the soil mapping and fertilization process in precision agriculture. Both soil prediction methods based on machine learning and remote sensing data are increasingly utilized in such studies.

  • remote sensing
  • machine learning
  • Fertilization
  • precision agriculture
  • variable rate technology

1. Modern Approach to Fertilization in Precision Agriculture

1.1. Remote Sensing Data

Novel multispectral and hyperspectral cameras on various platforms, from satellites to drones, enable high spatial and spectral remote sensing data that can be used in continuous agricultural field monitoring. Furthermore, with the development of Synthetic Aperture Radar (SAR; e.g., Sentinel-1) and optical satellite missions (e.g., Sentinel-2, PlanetScope), the Earth’s surface can be observed in high spatial resolution almost on a daily basis [1]. Satellite missions, but also new platforms such as unmanned aerial vehicles (UAVs) have led to the development and improvement of the quality of digital terrain models (DTM) and other important data in precision agriculture, e.g., the canopy height model (CHM). The development of remote sensing technologies and data lead to an increased application in precision agriculture. In the past few years, an increasing number of scientific articles in the field of remote sensing data in precision fertilization have been noted (Figure 1). A significant increase in the use of multispectral, hyperspectral and radar data is present, while the data on digital terrain models in precision fertilization are not growing at such a rate. Although the numbers of research with satellite data are significantly higher than the data obtained by UAV, a rapid growth in UAV application in the last five years has occurred. As in previous research [2], the prevalence of multispectral images is still higher than for radar. With the development of radar imaging technologies [3][4][5][6], mainly in soil moisture assessment, a significant increase in the number of studies in this area in the coming years is expected.
Figure 1. The most frequently used remote sensing data in precision fertilization according to the number of scientific papers indexed in the Web of Science Core Collection database.
The above-mentioned remote sensing technology and sensors enable fast and accurate acquisition of the earth surface data. This data can describe the earth’s surface with 3D (e.g., digital elevation models, digital surface model) or in the 2D data type (e.g., spectral bands, composites) acquired by sensors located on various platforms, e.g., satellite, airplane, or UAV. The acquired data can be used to produce various crucial remote sensing data such as [7]:
(1) various indices for an improved description of the earth’s surface (e.g., water, vegetation, soil) based on multispectral images or
(2) various derivatives of the digital elevation models such as slope, curvature, or flow accumulation analysis.
Radar images must be emphasized while they find their application in detecting the amount of moisture in the soil [3][6], which is extremely important in precision agriculture. For example, the most used remote sensing data in precision fertilization for the exact location is shown in Figure 2.
Figure 2. The most widely used remote sensing data in the modern approach to precision fertilization.

1.2. Modern Remote Sensing Methods for Optimal Fertilization in Precision Agriculture

In the recent research [8][9][10][11][12][13], remote sensing methods for assessing precision agriculture and fertilization are divided into two different approaches:
(1) Multivariate regressions based on various remote sensing data and
(2) Machine and deep learning methods for predictions.
For many years, researchers have used multivariate regressions to model the variables in precision agriculture [11][14]. With the development of remote sensing data, researchers have introduced new measurements collected using remote sensing data [11] into the process of multivariate regressions. These data have greatly improved the process of estimating the required variables to create better and more productive crops [15]. Multivariate regressions based on various remote sensing data have greatly improved the estimates of soil parameters for precision fertilization in regards to the conventional spatial interpolation methods.
In the last decade, hardware and computer development have enabled and popularized the application of and development of advanced methods for processing large amounts of spatial data, such as machine learning and deep learning methods [16][17][18][19]. Nowadays, these methods are also often used in remote sensing to quickly and accurately classify large amounts of satellite images to obtain land use and land cover (LULC) maps of the earth’s surface [20][21]. The great importance of the machine and deep learning method lies in the rapid processing of various spatial and attribute data to improve the estimation and prediction of the variable, e.g., estimation of the fertilizer in crops from available satellite images [9][22][23]. Trends of the application of remote sensing methods in precision fertilization based on the research papers in the last decade are shown in Figure 3. Application of the random forest and support vector machine algorithm have shown an increasing trend during the previous five years, while the artificial neural network and decision tree methods applications have stagnated in precision fertilization.
Figure 3. The most frequently used remote sensing methods in precision fertilization according to the number of scientific papers indexed in the Web of Science Core Collection database.
Despite the existing conventional modern remote sensing methods, developing novel hybrid approaches would significantly improve the fertilization process in precision agriculture. Novel hybrid approaches would enable better spatial resolution and more accurate crop and soil conditions assessment, thus improving the fertilization process. The preliminary proposal of a novel hybrid approach is based on remote sensing data such as for vegetation, water and soil spectral indices, a digital elevation model and other remote sensing data (e.g., spectral bands) supported by raster data obtained by conventional spatial interpolation methods. Machine learning methods could quickly and accurately model all the mentioned data based on which of the above can be realized into a real and accurate prediction of the state of crops and soil.

2. A Representative Overview of Modern and Conventional Approaches for Fertilization in Precision Agriculture

To provide a more insightful overview of the efficiency of conventional and modern approaches to precision fertilization, a case study for a common agricultural parcel in Croatia was performed. Phosphorous pentoxide (P2O5) and potassium oxide (K2O) were used for the prediction, representing two of the most important soil properties in agricultural fertilization [24]. A total of 121 samples were used in the agricultural parcel of 4.1 km2, representing a micro-location such as in [25][26]. The descriptive statistics of the input soil sample set is presented in Table 1. All predicted results were calculated in a spatial resolution of 30 m, according to the specifications of Hengl [27]. Soil prediction accuracy was assessed using R2 and RMSE, as the most commonly applied interpolation metrics in similar previous studies. The cross-validation using the leave-one-out technique was used for the accuracy assessment.
Table 1.
The descriptive statistics of the representative soil sample set used for the comparison of conventional and modern approaches to precision fertilization.
Soil

Property
Average

(mg 100 g–1)
Value Range

(mg 100 g–1)
CV SK KT Shapiro–Wilk Test Moran’s I
W p
phosphorous

pentoxide (P2O5)
23.2 8.9–41.0 0.364 0.587 –0.592 0.941 0.0005 0.209
potassium oxide

(K2O)
26.1 17.2–50.5 0.253 1.517 3.092 0.877 < 0.0001 0.124
OK, as the most commonly applied geostatistical interpolation method, along with IDW, its deterministic counterpart, were used for the representation of the conventional prediction methods in fertilization. As the range of the OK interpolation is conditioned by spatial autocorrelation of the input values, some of the most commonly used mathematical models in previous studies were evaluated, as the primary parameter of the OK interpolation. Analogously, the most common power parameters of the IDW were evaluated. Due to the lack of data normality, a logarithmic transformation was performed in the preprocessing to OK interpolation.
Comparative displays of interpolation results produced by conventional interpolation methods for P2O5 and K2O on the representative soil sample set are shown in Figure 4 and Figure 5. The interpolation results for both soil properties indicated a strong dependence of the prediction accuracy on the input parameters, indicating the importance of evaluating multiple methods, as well as their parameters, as noted in [28]. The R2 of the OK ranged from 0.331 to 0.414 for the P2O5 and from 0.082 to 0.120 for the K2O, indicating a proportionally lower accuracy for the input values with lower spatial autocorrelation, which is one of the main constraints of its prediction accuracy [29]. Due to its deterministic nature, IDW was resistant to this property, with its interpolation accuracy ranging from 0.233 to 0.405 for the P2O5 and from 0.234 to 0.374 for the K2O. It generally produced a lower accuracy but with a more balanced approach regarding sensitivity to the input values, as noted in [25]. Besides varying the interpolation accuracy, the resulting value ranges and CV values were severely affected by the selection of the interpolation method and its parameters.
Figure 4. Comparative presentation of interpolation results using common parameters of OK and IDW for P2O5.
Figure 5. Comparative presentation of interpolation results using common parameters of OK and IDW for K2O.
A total of twelve relevant covariates for the P2O5 and K2O prediction used as the basis of the modern prediction approach are presented in Table 2. These were defined with accordance to the specifications of soil mapping by Hengl and MacMillan [30] and which were used in similar soil prediction studies recently [31][32][7]. Six covariates were derived from a digital elevation model and six from Landsat 8 images, fully based on freely and widely available data. These covariates for the area covering the representative soil sample set are visually represented in Figure 1. Four of the most commonly applied machine learning methods in previous studies indexed in the WoSCC were used: random forest (RF), support vector machine (SVM), artificial neural networks (ANN) and decision tree (DT). These methods gained popularity in the modern approach to fertilization recently, allowing the integration of big data, that is highly accurate and with a computationally efficient prediction [33].
Table 2. Covariates used for the modern prediction using the representative soil sample set for precision fertilization.
Data Source Environmental Segment Covariate Reference
digital elevation model

(EU-DEM v1.1)
morphometry slope [34]
aspect
total curvature
convergence index
hydrology flow accumulation [35]
topographic wetness index [36]
multispectral satellite images

(Landsat 8, sensed on 15th September 2021)
vegetation normalized difference vegetation index (NDVI) [37]
enhanced vegetation index (EVI) [38]
normalized green-red vegetation index (NGRDI) [39]
soil normalized difference soil index (NDSI) [40]
brightness index (BI) [41]
moisture normalized difference moisture index (NDMI) [42]
A comparative display of the modern soil prediction approach for fertilization in precision agriculture, along with the most accurate results of the conventional approach are displayed in Figure 6. Besides an improved prediction accuracy and resistance to the particular properties of input sample values, the modern approach included more soil heterogeneity in the result with higher CV values. Previous studies have also noted a superior prediction accuracy of the modern approach compared to conventional methods, especially in the cases of lower spatial autocorrelation indicated by high nugget values [43]. Besides the spectral indices and topographic indicators, which are applicable at both minor- and major-scales [7], climate data and auxiliary soil information are commonly included in the modern approach [32]. These values are generally more suitable for the macro-location studies due to their local homogeneity, as well as the lack of available spatial data at the higher spatial resolution to match those of satellite images and DEMs [44]. Despite the same spatial resolution of the P2O5 and K2O rasters produced by the conventional and modern approaches, modern machine learning methods have resulted in much less smooth areas, retaining specific local information about field conditions, which are a backbone for precision agriculture [45]. Chen et al. [46] noted the improved spatial resolution as one of the key advantages of the modern approach over the conventional interpolation methods, alongside improved prediction accuracy and time- and cost-efficiency.
Figure 6. A comparative display of the modern soil prediction with the most accurate results of the conventional approach on a representative soil sample set.

References

  1. Radocaj, D.; Obhodas, J.; Jurisic, M.; Gasparovic, M. Global Open Data Remote Sensing Satellite Missions for Land Monitoring and Conservation: A Review. Land 2020, 9, 402.
  2. Gasparović, M.; Dobrinić, D. Comparative Assessment of Machine Learning Methods for Urban Vegetation Mapping Using Multitemporal Sentinel-1 Imagery. Remote Sens. 2020, 12, 1952.
  3. El Hajj, M.; Baghdadi, N.; Zribi, M.; Bazzi, H. Synergic Use of Sentinel-1 and Sentinel-2 Images for Operational Soil Moisture Mapping at High Spatial Resolution over Agricultural Areas. Remote Sens. 2017, 9, 1292.
  4. Uddin, K.; Matin, M.A.; Meyer, F.J. Operational Flood Mapping Using Multi-Temporal Sentinel-1 SAR Images: A Case Study from Bangladesh. Remote Sens. 2019, 11, 1581.
  5. Gasparovic, M.; Klobucar, D. Mapping Floods in Lowland Forest Using Sentinel-1 and Sentinel-2 Data and an Object-Based Approach. Forests 2021, 12, 553.
  6. Balenzano, A.; Mattia, F.; Satalino, G.; Lovergine, F.P.; Palmisano, D.; Peng, J.; Marzahn, P.; Wegmuller, U.; Cartus, O.; Dabrowska-Zielinska, K.; et al. Sentinel-1 Soil Moisture at 1 Km Resolution: A Validation Study. Remote Sens. Environ. 2021, 263, 112554.
  7. Zeraatpisheh, M.; Garosi, Y.; Reza Owliaie, H.; Ayoubi, S.; Taghizadeh-Mehrjardi, R.; Scholten, T.; Xu, M. Improving the Spatial Prediction of Soil Organic Carbon Using Environmental Covariates Selection: A Comparison of a Group of Environmental Covariates. CATENA 2022, 208, 105723.
  8. Zha, H.; Miao, Y.; Wang, T.; Li, Y.; Zhang, J.; Sun, W.; Feng, Z.; Kusnierek, K. Improving Unmanned Aerial Vehicle Remote Sensing-Based Rice Nitrogen Nutrition Index Prediction with Machine Learning. Remote Sens. 2020, 12, 215.
  9. Maes, W.H.; Steppe, K. Perspectives for Remote Sensing with Unmanned Aerial Vehicles in Precision Agriculture. Trends Plant Sci. 2019, 24, 152–164.
  10. Chlingaryan, A.; Sukkarieh, S.; Whelan, B. Machine Learning Approaches for Crop Yield Prediction and Nitrogen Status Estimation in Precision Agriculture: A Review. Comput. Electron. Agric. 2018, 151, 61–69.
  11. Erler, A.; Riebe, D.; Beitz, T.; Loehmannsroeben, H.-G.; Gebbers, R. Soil Nutrient Detection for Precision Agriculture Using Handheld Laser-Induced Breakdown Spectroscopy (LIBS) and Multivariate Regression Methods (PLSR, Lasso and GPR). Sensors 2020, 20, 418.
  12. Das, K.; Twarakavi, N.; Khiripet, N.; Chattanrassamee, P.; Kijkullert, C. A Machine Learning Framework for Mapping Soil Nutrients with Multi-Source Data Fusion. In Proceedings of the 2021 IEEE International Geoscience and Remote Sensing Symposium IGARSS, Brussels, Belgium, 11–16 July 2021; pp. 3705–3708.
  13. Qiu, Z.; Ma, F.; Li, Z.; Xu, X.; Ge, H.; Du, C. Estimation of Nitrogen Nutrition Index in Rice from UAV RGB Images Coupled with Machine Learning Algorithms. Comput. Electron. Agric. 2021, 189, 106421.
  14. Selige, T.; Boehner, J.; Schmidhalter, U. High Resolution Topsoil Mapping Using Hyperspectral Image and Field Data in Multivariate Regression Modeling Procedures. Geoderma 2006, 136, 235–244.
  15. Paustian, M.; Theuvsen, L. Adoption of Precision Agriculture Technologies by German Crop Farmers. Precis. Agric. 2017, 18, 701–716.
  16. Zhang, L.; Zhang, L.; Du, B. Deep Learning for Remote Sensing Data A Technical Tutorial on the State of the Art. IEEE Geosci. Remote Sens. Mag. 2016, 4, 22–40.
  17. Lary, D.J.; Alavi, A.H.; Gandomi, A.H.; Walker, A.L. Machine Learning in Geosciences and Remote Sensing. Geosci. Front. 2016, 7, 3–10.
  18. Zhu, X.X.; Tuia, D.; Mou, L.; Xia, G.-S.; Zhang, L.; Xu, F.; Fraundorfer, F. Deep Learning in Remote Sensing: A Comprehensive Review and List of Resources. IEEE Geosci. Remote Sens. Mag. 2017, 5, 8–36.
  19. Maxwell, A.E.; Warner, T.A.; Fang, F. Implementation of Machine-Learning Classification in Remote Sensing: An Applied Review. Int. J. Remote Sens. 2018, 39, 2784–2817.
  20. Chatziantoniou, A.; Petropoulos, G.P.; Psomiadis, E. Co-Orbital Sentinel 1 and 2 for LULC Mapping with Emphasis on Wetlands in a Mediterranean Setting Based on Machine Learning. Remote Sens. 2017, 9, 1259.
  21. Zhang, C.; Sargent, I.; Pan, X.; Li, H.; Gardiner, A.; Hare, J.; Atkinson, P.M. Joint Deep Learning for Land Cover and Land Use Classification. Remote Sens. Environ. 2019, 221, 173–187.
  22. Sibanda, M.; Mutanga, O.; Rouget, M. Examining the Potential of Sentinel-2 MSI Spectral Resolution in Quantifying above Ground Biomass across Different Fertilizer Treatments. ISPRS J. Photogramm. Remote Sens. 2015, 110, 55–65.
  23. Basso, B.; Fiorentino, C.; Cammarano, D.; Schulthess, U. Variable Rate Nitrogen Fertilizer Response in Wheat Using Remote Sensing. Precis. Agric. 2016, 17, 168–182.
  24. Ge, Y.; Zhang, J.; Zhang, L.; Yang, M.; He, J. Long-Term Fertilization Regimes Affect Bacterial Community Structure and Diversity of an Agricultural Soil in Northern China. J. Soils Sediments 2008, 8, 43–50.
  25. Radočaj, D.; Jug, I.; Vukadinović, V.; Jurišić, M.; Gašparović, M. The Effect of Soil Sampling Density and Spatial Autocorrelation on Interpolation Accuracy of Chemical Soil Properties in Arable Cropland. Agronomy 2021, 11, 2430.
  26. Wang, Y.; Xiao, Z.; Aurangzeib, M.; Zhang, X.; Zhang, S. Effects of Freeze-Thaw Cycles on the Spatial Distribution of Soil Total Nitrogen Using a Geographically Weighted Regression Kriging Method. Sci. Total Environ. 2021, 763, 142993.
  27. Hengl, T. Finding the Right Pixel Size. Comput. Geosci. 2006, 32, 1283–1298.
  28. Radočaj, D.; Jurišić, M.; Antonić, O. Determination of Soil C:N Suitability Zones for Organic Farming Using an Unsupervised Classification in Eastern Croatia. Ecol. Indic. 2021, 123, 107382.
  29. Ye, H.; Lu, C.; Lin, Q. Investigation of the Spatial Heterogeneity of Soil Microbial Biomass Carbon and Nitrogen under Long-Term Fertilizations in Fluvo-Aquic Soil. PLoS ONE 2019, 14, e0209635.
  30. Hengl, T.; MacMillan, R.A. Predictive Soil Mapping with R.; OpenGeoHub Foundation: Wageningen, The Netherlands, 2019; ISBN 978-0-359-30635-0.
  31. Hengl, T.; de Jesus, J.M.; Heuvelink, G.B.M.; Gonzalez, M.R.; Kilibarda, M.; Blagotić, A.; Shangguan, W.; Wright, M.N.; Geng, X.; Bauer-Marschallinger, B.; et al. SoilGrids250m: Global Gridded Soil Information Based on Machine Learning. PLoS ONE 2017, 12, e0169748.
  32. Poggio, L.; de Sousa, L.M.; Batjes, N.H.; Heuvelink, G.B.M.; Kempen, B.; Ribeiro, E.; Rossiter, D. SoilGrids 2.0: Producing Soil Information for the Globe with Quantified Spatial Uncertainty. SOIL 2021, 7, 217–240.
  33. Dash, P.K.; Panigrahi, N.; Mishra, A. Identifying Opportunities to Improve Digital Soil Mapping in India: A Systematic Review. Geoderma Reg. 2022, 28, e00478.
  34. Olaya, V. Chapter 6 Basic Land-Surface Parameters. In Developments in Soil Science; Geomorphometry; Hengl, T., Reuter, H.I., Eds.; Elsevier: Amsterdam, The Netherlands, 2009; Volume 33, pp. 141–169.
  35. Freeman, T.G. Calculating Catchment Area with Divergent Flow Based on a Regular Grid. Comput. Geosci. 1991, 17, 413–422.
  36. Conrad, O.; Bechtel, B.; Bock, M.; Dietrich, H.; Fischer, E.; Gerlitz, L.; Wehberg, J.; Wichmann, V.; Boehner, J. System for Automated Geoscientific Analyses (SAGA) v. 2.1.4. Geosci. Model Dev. 2015, 8, 1991–2007.
  37. Rouse, J.; Haas, R.H.; Schell, J.A.; Deering, D.W. Monitoring Vegetation Systems in the Great Plains with ERTS. NASA Spec. Publ. 1974, 351, 309.
  38. Huete, A.; Didan, K.; Miura, T.; Rodriguez, E.P.; Gao, X.; Ferreira, L.G. Overview of the Radiometric and Biophysical Performance of the MODIS Vegetation Indices. Remote Sens. Environ. 2002, 83, 195–213.
  39. Hunt, E.R.; Cavigelli, M.; Daughtry, C.S.T.; Mcmurtrey, J.E.; Walthall, C.L. Evaluation of Digital Photography from Model Aircraft for Remote Sensing of Crop Biomass and Nitrogen Status. Precis. Agric. 2005, 6, 359–378.
  40. Deng, Y.; Wu, C.; Li, M.; Chen, R. RNDSI: A Ratio Normalized Difference Soil Index for Remote Sensing of Urban/Suburban Environments. Int. J. Appl. Earth Obs. Geoinf. 2015, 39, 40–48.
  41. Escadafal, R. Remote Sensing of Arid Soil Surface Color with Landsat Thematic Mapper. Adv. Space Res. 1989, 9, 159–163.
  42. Jin, S.; Sader, S.A. Comparison of Time Series Tasseled Cap Wetness and the Normalized Difference Moisture Index in Detecting Forest Disturbances. Remote Sens. Environ. 2005, 94, 364–372.
  43. Song, X.-D.; Liu, F.; Wu, H.-Y.; Cao, Q.; Zhong, C.; Yang, J.-L.; Li, D.-C.; Zhao, Y.-G.; Zhang, G.-L. Effects of Long-Term K Fertilization on Soil Available Potassium in East China. Catena 2020, 188, 104412.
  44. Radočaj, D.; Jurišić, M.; Gašparović, M.; Plaščak, I.; Antonić, O. Cropland Suitability Assessment Using Satellite-Based Biophysical Vegetation Properties and Machine Learning. Agronomy 2021, 11, 1620.
  45. Vizzari, M.; Santaga, F.; Benincasa, P. Sentinel 2-Based Nitrogen VRT Fertilization in Wheat: Comparison between Traditional and Simple Precision Practices. Agronomy 2019, 9, 278.
  46. Chen, S.; Arrouays, D.; Leatitia Mulder, V.; Poggio, L.; Minasny, B.; Roudier, P.; Libohova, Z.; Lagacherie, P.; Shi, Z.; Hannam, J.; et al. Digital Mapping of GlobalSoilMap Soil Properties at a Broad Scale: A Review. Geoderma 2022, 409, 115567.
More
Academic Video Service