the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Global near real-time 500 m 10-day FPAR dataset from MODIS and VIIRS for operational agricultural monitoring and crop yield forecasting
Abstract. Climate change and extreme weather events pose challenges to food security, emphasizing the need for reliable and timely monitoring of crop and rangeland conditions. For this purpose, long-term consistent Earth Observation datasets on vegetation conditions are typically used in early warning and crop yield forecast systems. However, the near-real-time (NRT) production of high quality datasets and the need to guarantee long-term records present various challenges. To address these, we present a NRT global dataset of Fraction of Photosynthetically Active Radiation (FPAR) at 500 m resolution, optimized for agricultural applications. Our dataset combines MODIS-FPAR (Collection 6.1) and VIIRS-FPAR (Collection 2) data, ensuring continuity from 2000 to well beyond 2030. We applied a robust filtering approach based on the Whittaker smoother to produce reliable FPAR estimates in NRT, accounting for sparse and irregular spaced observations due to cloud cover. The dataset is composed of two 10-day filtered timeseries: 1) MODIS-FPAR for 2000 to 2023, being the reference dataset, and 2) intercalibrated VIIRS-FPAR for 2018 onward. While several methods can effectively smooth and gap-fill FPAR data (i.e., using observations before and after the estimation date), our method is designed for optimal filtering in NRT (i.e., using only prior observations). Our approach yields six successive estimates of the same FPAR data point with increasing quality: a inital estimate immediately after the 10-day reference period, four subsequent estimates every 10 days using new observations, and a final consolidated estimate 90 days later. The implemented filtering ingests the available FPAR observations and their original quality assessment (QA) layers. To avoid unrealistic extrapolation when observations are sparse, we impose constraints, season and location specific, to FPAR estimates. We then intercalibrated the VIIRS-FPAR with the MODIS-FPAR filtered timeseries, using a mean difference correction approach, to ensure consistency between both series. This paper describes the filtering and intercalibration method used, the quality assessment of resulting timeseries, and details the obtained products and the corresponding QA layers. The NRT FPAR dataset is publicly available through the Joint Research Centre Data Catalogue, https://data.jrc.ec.europa.eu/dataset/1aac79d8-0d68-4f1c-a40f-b6e362264e50 (Seguini et al., 2025).
- Preprint
(13916 KB) - Metadata XML
- BibTeX
- EndNote
Status: final response (author comments only)
-
RC1: 'Comment on essd-2025-287', Hongliang Fang, 24 Jun 2025
-
RC2: 'Reply on RC1', Hongliang Fang, 25 Jun 2025
For L73-75, there are a few FPAR products existing (doi: 10.1109/TGRS.2018.2818929). For example:
1. GEOV2-FPAR product
Verger, A., Baret , F., Weiss, M., 2020. Algorithm Theoretical Basis Document - GEOV2/AVHRR: Leaf Area Index (LAI), Fraction of Absorbed Photosynthetically Active Radiation (FAPAR) and Fraction of green Vegetation Cover (FCOVER) from LTDR AVHRR (Issue I2.50), Report for THEIA-SP-44-0207-CREAF, https://www.theia-land.fr/wp-content/uploads/2022/03/THEIA-SP-44-0207-CREAF_I2.50.pdf, p. 51.
2. FPAR3g from Boston University (http://www.mdpi.com/2072-4292/5/2/927)
3. JRC-TIP FAPAR (https://doi.org/10.3390/rs11243055)
Citation: https://doi.org/10.5194/essd-2025-287-RC2 -
AC2: 'Reply on RC2', Lorenzo Seguini, 05 Sep 2025
RC2: For L73-75, there are a few FPAR products existing (doi: 10.1109/TGRS.2018.2818929). For example:
- GEOV2-FPAR product
Verger, A., Baret , F., Weiss, M., 2020. Algorithm Theoretical Basis Document - GEOV2/AVHRR: Leaf Area Index (LAI), Fraction of Absorbed Photosynthetically Active Radiation (FAPAR) and Fraction of green Vegetation Cover (FCOVER) from LTDR AVHRR (Issue I2.50), Report for THEIA-SP-44-0207-CREAF, https://www.theia-land.fr/wp-content/uploads/2022/03/THEIA-SP-44-0207-CREAF_I2.50.pdf, p. 51.
- FPAR3g from Boston University (http://www.mdpi.com/2072-4292/5/2/927)
- JRC-TIP FAPAR (https://doi.org/10.3390/rs11243055)
Authors: the reviewer is correct that several valuable FPAR timeseries have been published in recent years. Nonetheless, the mentioned sources did not respond to the operational needs of near-real-time (NRT) monitoring systems.
The QA4ECV black-sky FAPAR long-term record (referred to in RC1) aims to provide a long timeseries of intercalibrated NOAA-AVHHR FPAR data between 1982 and 2006 for retrospective studies. The product did not aim to serve operational systems with continuous NRT data flow. Additionally, the ground resolution of the product (0.05° x 0.05°, roughly 5km at the equator) is considered sub-optimal for crop monitoring.
Regarding the other three suggestions provided:
- GEOV2-FPAR products are operationally delivered within the Copernicus Land Monitoring Service (CLMS). It consists of two timeseries. One has a spatial resolution of 1 km (based on SPOT/VEGETATION and Proba-V) and covers the time period 1998-2020. The second has a 300 m spatial resolution (based on Proba-V and Sentinel 3), it starts in 2014 and is updated in NRT. However, the latter offers only 12 years of data, which is too short to be used for robust statistical analysis (e.g. anomalies, crop yield forecast). No continuous long term (i.e. since 2000) and NRT product is currently available from CLMS. A possible combination of the three sources (SPOT/VEGETATION, Proba-V and Sentinel-3) would be difficult as requires the harmonization of FPAR data produced by different algorithms and different spatial resolutions.
- GIMMS FPAR3g offers a timeseries from 1981 to 2011 while the update product (GIMMS FPAR4g, Zhao et al., 2024) released in 2024, covers from 1982 to 2022. The products offer low temporal (15-days) and spatial (1/12 degree, roughly 8km) resolution and no NRT production is available. These elements frame both products as unsuitable for operational agricultural monitoring purposes.
- JRC-TIP FAPAR: the reference proposed for this product is the same as QA4ECV black-sky FAPAR long-term records. We acknowledge that TIP (Two-stream Inversion Package) FPAR is a well-known and used FPAR retrieval algorithm, but no operational timeseries has been released so far.
Given these considerations, we have changed the introduction citing the products from GIMMS FPAR4g, CLMS and the timeseries cited by RC3 (HiQ-LAI). After line 84 we inserted the following paragraph:
“Several products exist that offer high-quality timeseries of biophysical variables, such as HiQ-LAI (Yan et al., 2025) or GIMMS FPAR4g (Zhao et al., 2024), but they usually lack some of the mandatory features needed by operational agricultural monitoring systems (i.e., long-term record and NRT availability). Typically, there is no guaranteed NRT data delivery into the future, nor are datasets filtered to reduce atmospheric influences. An exception is the Copernicus Land Monitoring Service that provides continuous, NRT, and filtered timeseries of biophysical variables from Proba-V and Sentinel-3 satellites (https://land.copernicus.eu/en/products/vegetation/fraction-of-absorbed-photosynthetically-active-radiation-v1-0-300m). Nevertheless, the timeseries length offered is too short (around 12 years) for robust statistical analysis needed for anomaly computation or crop yield forecasting. In this study we fill this gap by proposing a new dataset meeting the requirements of continuous, NRT, and filtered biophysical timeseries for more than 20 years.”
Citation: https://doi.org/10.5194/essd-2025-287-AC2
-
AC2: 'Reply on RC2', Lorenzo Seguini, 05 Sep 2025
-
AC1: 'Reply on RC1', Lorenzo Seguini, 05 Sep 2025
RC1:the work presents a NRT global dataset of FPAR generated from existing MODIS and VIIRS datasets. The new dataset combines MODIS-FPAR for 2000 to 2023, and the intercalibrated VIIRS-FPAR for 2018 onward. Detailed procedures to smooth and filter the original datasets were introduced, together with fine quality assessment and intercalibration studies. The dataset is openly available and should be valuable for crop monitoring studies. I recommend its publication in ESSD.
Authors: We thank Hongliang Fang for reading our paper and recommending it for publication.
RC1: L73-75 I think JRC once generated a FPAR product from AVHRR. Please check it. https://www.mdpi.com/2072-4292/11/24/3055
Authors: We agree that JRC already generated a timeseries of FPAR data within the Quality Assurance for Essential Climate Variables (QA4ECV) project, but such timeseries was never released operationally (i.e., providing continuos and near-real-time data). Please see our reply to RC2 for a more reasoned analysis of the FPAR timeseries already available in the literature and our proposed modification to the manuscript.
RC1: L146. For the MCD12Q1 biome map, I think a better reference would be Friedl and Sulla-Menash (2019), as that cited for Fig. B1. Indeed, the MODIS LAI C6.1 (Myneni, 2020) has used the MCD12Q1 biome map.
Authors: We agree with the reviewer and have changed our reference as below.
“From the remaining, non-watered pixels, we selected all vegetated pixels accordingly to MCD12Q1 biome map (Friedl and Sulla-Menashe, 2019) which cover crops, shrubs, savanna, and forests. Since the main scope of our dataset is agricultural monitoring, we used two further masks (cropland and rangeland) from the JRC-ASAP system (Fritz et al., 2024) to analyse our results.”
Citation: https://doi.org/10.5194/essd-2025-287-AC1
-
RC2: 'Reply on RC1', Hongliang Fang, 25 Jun 2025
-
RC3: 'Comment on essd-2025-287', Kai Yan, 23 Jul 2025
This paper presents a comprehensive methodology for generating a global near real-time (NRT) FPAR dataset at 500m resolution with 10-day temporal resolution, combining MODIS and VIIRS satellite data for agricultural monitoring and crop yield forecasting applications. The authors employ a Whittaker smoother-based filtering approach with quality-based weighting and constraint mechanisms to handle sparse observations due to cloud cover. The dataset provides multiple consolidation stages (C0 to CF) representing increasing data quality over time, and includes intercalibrated VIIRS-FPAR data to ensure continuity beyond the MODIS mission lifetime. The work addresses an important operational need for consistent, gap-filled vegetation monitoring data for early warning systems and crop yield forecasting. Nevertheless, there are still improvements for this manuscript. My detailed comments are as follows:
Major Comments
- This work presents a valuable operational implementation of established techniques, specifically for agricultural monitoring. I'm curious about the specific adaptations made to the Whittaker smoother approach compared to previous implementations. Could the authors provide more detail on what aspects of their constraint mechanisms or gap-filling procedures are novel or represent improvements over existing methods? This would help readers better understand the technical contributions.
- The quality layers provide valuable information about processing conditions, but I'm curious about how measurement uncertainties from the original MODIS and VIIRS products affect the final filtered results. This information would be particularly valuable for users needing to understand confidence levels in their applications.
- The evaluation metrics are appropriate. I wonder if seasonal or regional performance variations might provide additional insights into the method's behavior under different conditions.
Minor Comments
- The paper is generally well-written with clear explanations of the methodology. There are a few minor issues (e.g., "reangeland" should be "rangeland" in several places) that could be addressed during revision.
- More justification for the λ=3000 parameter selection would be helpful
- I note that "hindacsting" in line 215 should be corrected to "hindcasting" and "poral resolution" in line 90 to "temporal resolution". Please check for similar cases.
Citation: https://doi.org/10.5194/essd-2025-287-RC3 - AC3: 'Reply on RC3', Lorenzo Seguini, 05 Sep 2025
Data sets
A global near real-time filtered 500m 10-day FPAR dataset from MODIS and VIIRS instruments, suited for operational agricultural monitoring and crop yield forecasting Lorenzo Seguini, Anja Klisch, Michele Meroni, Clement Atzberger, Anton Vrieling, Giacinto Manfron, Felix Rembold https://doi.org/10.2905/1aac79d8-0d68-4f1c-a40f-b6e362264e50
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
847 | 138 | 28 | 1,013 | 19 | 42 |
- HTML: 847
- PDF: 138
- XML: 28
- Total: 1,013
- BibTeX: 19
- EndNote: 42
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1
The work presents a NRT global dataset of FPAR generated from existing MODIS and VIIRS datasets. The new dataset combines MODIS-FPAR for 2000 to 2023, and the intercalibrated VIIRS-FPAR for 2018 onward. Detailed procedures to smooth and filter the original datasets were introduced, together with fine quality assessment and intercalibration studies. The dataset is openly available and should be valuable for crop monitoring studies. I recommend its publication in ESSD.
L73-75 I think JRC once generated a FPAR product from AVHRR. Please check it
https://www.mdpi.com/2072-4292/11/24/3055
L146. For the MCD12Q1 biome map, I think a better reference would be Friedl and Sulla-Menash (2019), as that cited for Fig. B1. Indeed, the MODIS LAI C6.1 (Myneni, 2020) has used the MCD12Q1 biome map.