Mapping spatial variability of foliar nitrogen in coffee (Coffea arabica L.) plantations with multispectral Sentinel-2 MSI data

Nitrogen (N) is the most limiting factor to coffee development and productivity. Therefore, development of rapid, spatially explicit and temporal remote sensing-based approaches to determine spatial variability of coffee foliar N are imperative for increasing yields, reducing production costs and mitigating environmental impacts associated with excessive N applications. This study sought to assess the value of Sentinel-2 MSI spectral bands and vegetation indices in empirical estimation of coffee foliar N content at landscape level. Results showed that coffee foliar N is related to Sentinel-2 MSI B4 (R² = 0.32), B6 (R² = 0.49), B7 (R² = 0.42), B8 (R² = 0.57) and B12 (R² = 0.24) bands. Vegetation indices were more related to coffee foliar N as shown by the Inverted Red-Edge Chlorophyll Index – IRECI (R² = 0.66), Relative Normalized Difference Index – RNDVI (R² = 0.48), CIRE1 (R² = 0.28), and Normalized Difference Infrared Index – NDII (R² = 0.37). These variables were also identified by the random forest variable optimisation as the most valuable in coffee foliar N prediction. Modelling coffee foliar N using vegetation indices produced better accuracy (R² = 0.71 with RMSE = 0.27 for all and R² = 0.73 with RMSE = 0.25 for optimized variables), compared to using spectral bands (R² = 0.57 with RMSE = 0.32 for all and R² = 0.58 with RMSE = 0.32 for optimized variables). Combining optimized bands and vegetation indices produced the best results in coffee foliar N modelling (R² = 0.78, RMSE = 0.23). All the three best performing models (all vegetation indices, optimized vegetation indices and combining optimal bands and optimal vegetation indices) established that 15.2 ha (4.7%) of the total area under investigation had low foliar N levels (<2.5%). This study demonstrates the value of Sentinel-2 MSI data, particularly vegetation indices in modelling coffee foliar N at landscape scale.