MODIStsp is a new R package allowing automating the creation of raster time series derived from MODIS Land Products. It allows performing several preprocessing steps (e.g. download, mosaicing, reprojection and resize) on MODIS products on a selected time period and area. All processing parameters can be set with a user-friendly GUI, allowing users to select which specific layers of the original MODIS HDF files have to be processed and which Quality Indicators have to be extracted from the aggregated MODIS Quality Assurance layers. Moreover, the tool allows on-the-fly computation of time series of Spectral Indexes (either standard or custom-specified by the user through the GUI) from surface reflectance bands. Outputs are saved as single-band rasters corresponding to each available acquisition date and output layer. Virtual files allowing easy access to the entire time series as a single file using common image processing/GIS software or R scripts can be also created. Non-interactive execution within an R script and stand-alone execution outside an R environment exploiting a previously created Options File are also possible, the latter allowing scheduling execution of MODIStsp to automatically update a time series when a new image is available. The proposed software constitutes a very useful tool for the Remote Sensing community, since it allows performing all the main preprocessing steps required for the creation of time series of MODIS data within a common framework, and without requiring any particular programming skills by its users

Leaf area index (LAI) is a key biophysical variable fundamental in natural vegetation and agricultural land monitoring and modelling studies. This paper is aimed at comparing, validating and discussing different LAI satellite products from operational services and customized solution based on innovative Earth Observation (EO) data such as Landsat-7/8 and Sentinel-2A. The comparison was performed to assess overall quality of LAI estimates for rice, as a fundamental input of different scale (regional to local) operational crop monitoring systems such as the ones developed during the “An Earth obseRvation Model based RicE information Service” (ERMES) project. We adopted a multiscale approach following international recognized protocols of the Committee on Earth Observation Satellites (CEOS) Land Product Validation (LPV) guidelines in different steps: (1) acquisition of representative field sample measurements, (2) validation of decametric satellite product (10–30 m spatial resolution), and (3) exploitation of such data to assess quality of medium-resolution operational products (~1000 m). The study areas were located in the main European rice areas in Spain, Italy and Greece. Field campaigns were conducted during three entire rice seasons (2014, 2015 and 2016—from sowing to full-flowering) to acquire multi-temporal ground LAI measurements and to assess Landsat-7/8 LAI estimates. Results highlighted good correspondence between Landsat-7/8 LAI estimates and ground measurements revealing high correlations (R² ≥ 0.89) and low root mean squared errors (RMSE ≤ 0.75) in all seasons. Landsat-7/8 as well as Sentinel-2A high-resolution LAI retrievals, were compared with satellite LAI products operationally derived from MODIS (MOD15A2), Copernicus PROBA-V (GEOV1), and the recent EUMETSAT Polar System (EPS) LAI product. Good agreement was observed between high- and medium-resolution LAI estimates. In particular, the EPS LAI product was the most correlated product with both Landsat/7-8 and Sentinel-2A estimates, revealing R² ≥ 0.93 and RMSE ≤ 0.53 m²/m². In addition, a comparison exercise of EPS, GEOV1 and MODIS revealed high correlations (R² ≥ 0.90) and RMSE ≤ 0.80 m²/m² in all cases and years. The temporal assessment shows that the three satellite products capture well the seasonality during the crop phenological cycle. Discrepancies are observed mainly in absolute values retrieved for the peak of rice season. This is the first study that provides a quantitative assessment on the quality of available operational LAI product for rice monitoring to both the scientific community and users of agro-monitoring operational services.

The ERMES agromonitoring system for rice cultivations integrates EO data at different resolutions, crop models, and user-provided in situ data in a unified system, which drives two operational downstream services for rice monitoring. The first is aimed at providing information concerning the behavior of the current season at regional/rice district scale, while the second is dedicated to provide farmers with field-scale data useful to support more efficient and environmentally friendly crop practices. In this contribution, we describe the main characteristics of the system, in terms of overall architecture, technological solutions adopted, characteristics of the developed products, and functionalities provided to end users. Peculiarities of the system reside in its ability to cope with the needs of different stakeholders within a common platform, and in a tight integration between EO data processing and information retrieval, crop modeling, in situ data collection, and information dissemination. The ERMES system has been operationally tested in three European rice-producing countries (Italy, Spain, and Greece) during growing seasons 2015 and 2016, providing a great amount of near-real-time information concerning rice crops. Highlights of significant results are provided, with particular focus on real-world applications of ERMES products and services. Although developed with focus on European rice cultivations, solutions implemented in the ERMES system can be, and are already being, adapted to other crops and/or areas of the world, thus making it a valuable testing bed for the development of advanced, integrated agricultural monitoring systems.