Detection of grassland mowing frequency using time series of vegetation indices from Sentinel-2 imagery

Management intensity deeply influences meadow structure and functioning, therefore affecting grassland ecosystem services. Conservation and management measures, including European Common Agricultural Policy subsidies, should therefore be based on updated and publicly available data about management intensity. The mowing frequency is a crucial trait to describe meadows management intensity, but the potential of using vegetation indices from Sentinel-2 imagery for its retrieval has not been fully exploited. In this work we developed on the Google Earth Engine platform a four-phases algorithm to identify mowing frequency, including i) vegetation index time-series computing, ii) smoothing and resampling, iii) mowing detection, and iv) majority analysis. Mowing frequency during 2020 of 240 ha of grassland fields in the Italian Alps was used for algorithm optimization and evaluation. Six vegetation indexes (EVI, GVMI, MTCI, NDII, NDVI, RENDVI 783.740 ) were tested as input to the proposed algorithm. The Normalized Difference Infrared Index (NDII) showed the best performance, resulting in mean absolute error of 0.07 and 93% overall accuracy on average at the four sites used for optimization, at pixel resolution. A slightly lower accuracy (mean absolute error = 0.10, overall accuracy = 90%) was obtained aggregating the maps to management parcels. The algorithm showed a good generalization ability, with a similar performance between global and local optimization and an average mean absolute error of 0.12 and an overall accuracy of 89% on average on the sites not used for parameters optimization. The lowest accuracies occurred in intensively managed grasslands surveyed by one satellite orbit only. This study demonstrates the suitability of the proposed algorithm to monitor very fragmented grasslands in complex mountain ecosystems. Google Earth Engine was used to develop the model and will enable researchers, agencies and practitioners to easily and quickly apply the code to map grassland mowing frequency for extensive grasslands protection and conservation, for mowing event verification, or for forage system characterization.