Evaluation of Global Fraction of Absorbed Photosynthetically Active Radiation (FAPAR) Products at 500 m Spatial Resolution

The fraction of absorbed photosynthetically active radiation (FAPAR) is a key biophysical variable directly associated with the photosynthetic activity of plants. Several global FAPAR products with different spatial resolutions have been generated from remote sensing data, and much work has focused on validating them. However, those studies have primarily evaluated global FAPAR products at a spatial resolution of 1 km or more, whereas few studies have evaluated the global 500 m resolution FAPAR product distributed in recent years. Furthermore, there are a few FAPAR products, including black-sky, white-sky and blue-sky FAPAR datasets, and almost no studies have evaluated these products. In this article, three global FAPAR products at 500 m resolution, namely MODIS (only black-sky FAPAR), MUSES and EBR (black-sky, white-sky and blue-sky FAPAR) were compared to evaluate their temporal and spatial discrepancies and direct validation was conducted to compare these FAPAR products with the FAPAR values derived from the high-resolution reference maps from the Validation of Land European Remote Sensing Instrument (VALERI) and Implementing Multi-Scale Agricultural Indicators Exploiting Sentinels (IMAGINES) projects. The results showed that the MUSES FAPAR product exhibited the best spatial integrity, whereas the MODIS and EBR FAPAR products had many missing pixels in the equatorial rainforest regions and at high latitudes in the Northern Hemisphere. The MODIS, MUSES and EBR FAPAR products were generally consistent in their spatial patterns. However, a relatively large discrepancy among these FAPAR products was present in the equatorial rainforest regions and the middle and high latitude regions where the main vegetation type was forest. The differences between the black-sky and white-sky FAPAR datasets at the global scale were significant. In January, the MUSES and EBR black-sky FAPAR values were larger than their white-sky FAPAR values in the region north of 30° north latitude but they were smaller than their white-sky FAPAR values in the region south of 30° north latitude. In July, the MUSES and EBR black-sky FAPAR values were lower than their white-sky FAPAR values in the region north of 30° south latitude and they were larger than their white-sky FAPAR values in the region south of 30° south latitude. The temporal profiles of the MUSES FAPAR product were continuous and smooth, whereas those of the MODIS and EBR FAPAR products showed many fluctuations, particularly during the growing seasons. Direct validation indicated that the MUSES FAPAR product had the best accuracy (R2 = 0.6932, RMSE = 0.1495) compared to the MODIS FAPAR product (R2 = 0.6202, RMSE = 0.1710) and the EBR FAPAR product (R2 = 0.5746, RMSE = 0.1912).

Cross-ECV consistency at global scale: LAI and FAPAR changes

A framework is proposed for assessing the physical consistency between two terrestrial Essential Climate Variables (ECVs) products retrieved from Earth Observation at global scale. The methodology assessed the level of agreement between the temporal variations of Leaf Area Index (LAI) and Fraction of Absorbed Photosynthetically Active Radiation (FAPAR). The simultaneous changes were classified according to their sign, magnitude and level of confidence, whereby the respective products uncertainties were taken into consideration. A set of proposed agreement metrics were used to identify temporal and spatial biases of non-coherency, non-significance, sensitivity and the overall level of agreement of the temporal changes between two ECVs. We applied the methodology using the Joint Research Center (JRC) Two-stream Inversion Package (TIP) products at 1 km, those provided by the Copernicus Global Land Service (CGLS) based on the SPOT/VGT and Proba-V at 1 km, and the MODIS MCD15A3 at 500 m. In addition, the same analysis was applied with aggregated products at a larger scale over Southern Africa. We found that the CGLS LAI and FAPAR products lacked consistency in their spatial and temporal changes and were severely affected by trends. The MCD15A3 products were characterized by the highest number of non-coherent changes between the two ECVs but temporal inconsistencies were mainly located over the eastern hemisphere. The JRC-TIP products were highly consistent. The results showed the advantages of physically-based retrieval algorithms, in both JRC-TIP and MODIS products, and indicated also that, except for MODIS over forests, aggregated products using an uncertainty-based weighted average led to higher agreement between the ECVs changes.

Fiducial Reference Measurements for Vegetation Bio-Geophysical Variables: An End-to-End Uncertainty Evaluation Framework

With a wide range of satellite-derived vegetation bio-geophysical products now available to users, validation efforts are required to assess their accuracy and fitness for purpose. Substantial progress in the validation of such products has been made over the last two decades, but quantification of the uncertainties associated with in situ reference measurements is rarely performed, and the incorporation of uncertainties within upscaling procedures is cursory at best. Since current validation practices assume that reference data represent the truth, our ability to reliably demonstrate compliance with product uncertainty requirements through conformity testing is limited. The Fiducial Reference Measurements for Vegetation (FRM4VEG) project, initiated by the European Space Agency, is aiming to address this challenge by applying metrological principles to vegetation and surface reflectance product validation. Following FRM principles, and in accordance with the International Standards Organisation’s (ISO) Guide to the Expression of Uncertainty in Measurement (GUM), for the first time, we describe an end-to-end uncertainty evaluation framework for reference data of two key vegetation bio-geophysical variables: the fraction of absorbed photosynthetically active radiation (FAPAR) and canopy chlorophyll content (CCC). The process involves quantifying the uncertainties associated with individual in situ reference measurements and incorporating these uncertainties within the upscaling procedure (as well as those associated with the high-spatial-resolution imagery used for upscaling). The framework was demonstrated in two field campaigns covering agricultural crops (Las Tiesas–Barrax, Spain) and deciduous broadleaf forest (Wytham Woods, UK). Providing high-spatial-resolution reference maps with per-pixel uncertainty estimates, the framework is applicable to a range of other bio-geophysical variables including leaf area index (LAI), the fraction of vegetation cover (FCOVER), and canopy water content (CWC). The proposed procedures will facilitate conformity testing of moderate spatial resolution vegetation bio-geophysical products in future validation exercises.

Parametric Models to Characterize the Phenology of the Lowveld Savanna at Skukuza, South Africa

Mathematical models, such as the logistic curve, have been extensively used to model the temporal evolution of biological processes, though other similarly shaped functions could be (and sometimes have been) used for this purpose. Most previous studies focused on agricultural regions in the Northern Hemisphere and were based on the Normalized Difference Vegetation Index (NDVI). This paper compares the capacity of four parametric double S-shaped models (Gaussian, Hyperbolic Tangent, Logistic, and Sine) to represent the seasonal phenology of an unmanaged, protected savanna biome in South Africa’s Lowveld, using the Fraction of Absorbed Photosynthetically Active Radiation (FAPAR) generated by the Multi-angle Imaging SpectroRadiometer-High Resolution (MISR-HR) processing system on the basis of data originally collected by National Aeronautics and Space Administration (NASA)’s Multi-angle Imaging SpectroRadiometer (MISR) instrument since 24 February 2000. FAPAR time series are automatically split into successive vegetative seasons, and the models are inverted against those irregularly spaced data to provide a description of the seasonal fluctuations despite the presence of noise and missing values. The performance of these models is assessed by quantifying their ability to account for the variability of remote sensing data and to evaluate the Gross Primary Productivity (GPP) of vegetation, as well as by evaluating their numerical efficiency. Simulated results retrieved from remote sensing are compared to GPP estimates derived from field measurements acquired at Skukuza’s flux tower in the Kruger National Park, which has also been operational since 2000. Preliminary results indicate that (1) all four models considered can be adjusted to fit an FAPAR time series when the temporal distribution of the data is sufficiently dense in both the growing and the senescence phases of the vegetative season, (2) the Gaussian and especially the Sine models are more sensitive than the Hyperbolic Tangent and Logistic to the temporal distribution of FAPAR values during the vegetative season, and, in particular, to the presence of long temporal gaps in the observational data, and (3) the performance of these models to simulate the phenology of plants is generally quite sensitive to the presence of unexpectedly low FAPAR values during the peak period of activity and to the presence of long gaps in the observational data. Consequently, efforts to screen out outliers and to minimize those gaps, especially during the rainy season (vegetation’s growth phase), would go a long way to improve the capacity of the models to adequately account for the evolution of the canopy cover and to better assess the relation between FAPAR and GPP.

Upscaling from Instantaneous to Daily Fraction of Absorbed Photosynthetically Active Radiation (FAPAR) for Satellite Products

The fraction of absorbed photosynthetically active radiation (FAPAR) is an essential climate variable (ECV) widely used for various ecological and climate models. However, all the current FAPAR satellite products correspond to instantaneous FAPAR values acquired at the satellite transit time only, which cannot represent the variations in photosynthetic processes over the diurnal period. Most studies have directly used the instantaneous FAPAR as a reasonable approximation of the daily integrated value. However, clearly, FAPAR varies a lot according to the weather conditions and amount of incoming radiation. In this paper, a temporal upscaling method based on the cosine of the solar zenith angle (SZA) at local noon ( c o s ( S Z A n o o n ) ) is proposed for converting instantaneous FAPAR to daily integrated FAPAR. First, the diurnal variations in FAPAR were investigated using PROSAIL (a model of Leaf Optical Properties Spectra (PROSPECT) integrating a canopy radiative transfer model (Scattering from Arbitrarily Inclined Leaves, SAIL)) simulations with different leaf area index (LAI) values corresponding to different latitudes. It was found that the instantaneous black sky FAPAR at 09:30 AM provided a good approximation for the daily integrated black sky FAPAR; this gave the highest correlation (R2 = 0.995) and lowest Root Mean Square Error (RMSE = 0.013) among the instantaneous black sky FAPAR values observed at different times. Secondly, the difference between the instantaneous black sky FAPAR values acquired at different times and the daily integrated black sky FAPAR was analyzed; this could be accurately modelled using the cosine value of solar zenith angle at local noon ( c o s ( S Z A n o o n ) ) for a given vegetation scene. Therefore, a temporal upscaling method for typical satellite products was proposed using a cos(SZA)-based upscaling model. Finally, the proposed cos(SZA)-based upscaling model was validated using both the PROSAIL simulated data and the field measurements. The validated results indicated that the upscaled daily black sky FAPAR was highly consistent with the daily integrated black sky FAPAR, giving very high mean R2 values (0.998, 0.972), low RMSEs (0.007, 0.014), and low rMAEs (0.596%, 1.378%) for the simulations and the field measurements, respectively. Consequently, the cos(SZA)-based method performs well for upscaling the instantaneous black sky FAPAR to its daily value, which is a simple but extremely important approach for satellite remote sensing applications related to FAPAR.

Can We Use the QA4ECV Black-sky Fraction of Absorbed Photosynthetically Active Radiation (FAPAR) using AVHRR Surface Reflectance to Assess Terrestrial Global Change?

NOAA platforms provide the longest period of terrestrial observation since the 1980s. The progress in calibration, atmospheric corrections and physically based land retrieval offers the opportunity to reprocess these data for extending terrestrial product time series. Within the Quality Assurance for Essential Climate Variables (QA4ECV) project, the black-sky Joint Research Centre (JRC)-fraction of absorbed photosynthetically active radiation (FAPAR) algorithm was developed for the AVHRR sensors on-board NOAA-07 to -16 using the Land Surface Reflectance Climate Data Record. The retrieval algorithm was based on the radiative transfer theory, and uncertainties were included in the products. We proposed a time and spatial composite for providing both 10-day and monthly products at 0.05º × 0.05º. Quality control and validation were achieved through benchmarking against third-party products, including Sea-Viewing Wide Field-of-View Sensor (SeaWiFS) datasets produced with the same retrieval algorithm. Past ground-based measurements, providing a proxy of FAPAR, showed good agreement of seasonality values over short homogeneous canopies and mixed vegetation. The average difference between SeaWiFS and QA4ECV monthly products over 2002–2005 is about 0.075 with a standard deviation of 0.091. We proposed a monthly linear bias correction that reduced these statistics to 0.02 and 0.001. The complete harmonized long-term time series was then used to address its fitness for the purpose of analysis of global terrestrial change.

Retrieval of the Fraction of Radiation Absorbed by Photosynthetic Components (FAPARgreen) for Forest using a Triple-Source Leaf-Wood-Soil Layer Approach

The fraction of absorbed photosynthetically active radiation (FAPAR) is generally divided into the fraction of radiation absorbed by the photosynthetic components (FAPARgreen) and the fraction of radiation absorbed by the non-photosynthetic components (FAPARwoody) of the vegetation. However, most global FAPAR datasets do not take account of the woody components when considering the canopy radiation transfer. The objective of this study was to develop a generic algorithm for partitioning FAPARcanopy into FAPARgreen and FAPARwoody based on a triple-source leaf-wood-soil layer (TriLay) approach. The LargE-Scale remote sensing data and image simulation framework (LESS) model was used to validate the TriLay approach. The results showed that the TriLay FAPARgreen had higher retrieval accuracy, as well as a significantly lower bias (R2 = 0.937, Root Mean Square Error (RMSE) = 0.064, and bias = −6.02% for black-sky conditions; R2 = 0.997, RMSE = 0.025 and bias = −4.04% for white-sky conditions) compared to the traditional linear method (R2 = 0.979, RMSE = 0.114, and bias = −18.04% for black-sky conditions; R2 = 0.996, RMSE = 0.106 and bias = −16.93% for white-sky conditions). For FAPAR that did not take account of woody components (FAPARnoWAI), the corresponding results were R2 = 0.920, RMSE = 0.071, and bias = −7.14% for black-sky conditions, and R2 = 0.999, RMSE = 0.043, and bias = −6.41% for white-sky conditions. Finally, the dynamic FAPARgreen, FAPARwoody, FAPARcanopy and FAPARnoWAI products for a North America region were generated at a resolution of 500 m for every eight days in 2017. A comparison of the results for FAPARgreen against those for FAPARnoWAI and FAPARcanopy showed that the discrepancy between FAPARgreen and other FAPAR products for forest vegetation types could not be ignored. For deciduous needleleaf forest, in particular, the black-sky FAPARgreen was found to contribute only about 23.86% and 35.75% of FAPARcanopy at the beginning and end of the year (from January to March and October to December, JFM and OND), and 75.02% at the peak growth stage (from July to September, JAS); the black-sky FAPARnoWAI was found to be overestimated by 38.30% and 28.46% during the early (JFM) and late (OND) part of the year, respectively. Therefore, the TriLay approach performed well in separating FAPARgreen from FAPARcanopy, which is of great importance for a better understanding of the energy exchange within the canopy.

Can We Use Satellite-Based FAPAR to Detect Drought?

Drought in Australia has widespread impacts on agriculture and ecosystems. Satellite-based Fraction of Absorbed Photosynthetically Active Radiation (FAPAR) has great potential to monitor and assess drought impacts on vegetation greenness and health. Various FAPAR products based on satellite observations have been generated and made available to the public. However, differences remain among these datasets due to different retrieval methodologies and assumptions. The Quality Assurance for Essential Climate Variables (QA4ECV) project recently developed a quality assurance framework to provide understandable and traceable quality information for Essential Climate Variables (ECVs). The QA4ECV FAPAR is one of these ECVs. The aim of this study is to investigate the capability of QA4ECV FAPAR for drought monitoring in Australia. Through spatial and temporal comparison and correlation analysis with widely used Moderate Resolution Imaging Spectroradiometer (MODIS), Satellite Pour l’Observation de la Terre (SPOT)/PROBA-V FAPAR generated by Copernicus Global Land Service (CGLS), and the Standardized Precipitation Evapotranspiration Index (SPEI) drought index, as well as the European Space Agency’s Climate Change Initiative (ESA CCI) soil moisture, the study shows that the QA4ECV FAPAR can support agricultural drought monitoring and assessment in Australia. The traceable and reliable uncertainties associated with the QA4ECV FAPAR provide valuable information for applications that use the QA4ECV FAPAR dataset in the future.

Global White-Sky and Black-Sky FAPAR Retrieval Using the Energy Balance Residual Method: Algorithm and Validation

The fraction of absorbed photosynthetically active radiation by vegetation (FAPAR) is a key variable in describing the light absorption ability of the vegetation canopy. Most global FAPAR products, such as MCD15A2H and GEOV1, correspond to FAPAR under black-sky conditions at the satellite overpass time only. In this paper, we aim to produce both the global white-sky and black-sky FAPAR products based on the moderate resolution imaging spectroradiometer (MODIS) visible (VIS) albedo, leaf area index (LAI), and clumping index (CI) products. Firstly, a non-linear spectral mixture model (NSM) was designed to retrieve the soil visible (VIS) albedo. The global soil VIS albedo and its dynamics were successfully mapped at a resolution of 500 m using the MCD43A3 VIS albedo product and the MCD15A2H LAI product. Secondly, a method based on the energy balance residual (EBR) principle was presented to retrieve the white-sky and black-sky FAPAR using the MODIS broadband VIS albedo (white-sky and black-sky) product (MCD43A3), the LAI product (MCD15A2H) and CI products. Finally, the two EBR FAPAR products were compared with the MCD15A2H and Geoland2/BioPar version 1 (GEOV1) black-sky FAPAR products. A comparison of the results indicates that these FAPAR products show similar spatial and seasonal patterns. Direct validation using FAPAR observations from the Validation of Land European Remote sensing Instrument (VALERI) project demonstrates that the EBR black-sky FAPAR product was more accurate and had a lower bias (R2 = 0.917, RMSE = 0.088, and bias = −2.8 %) than MCD15A2H (R2 = 0.901, RMSE = 0.096, and bias = 7.6 % ) and GEOV1 (R2 = 0.868, RMSE = 0.105, and bias = 6.1%).

Mapping Climatological Bare Soil Albedos over the Contiguous United States Using MODIS Data

Surface bare soil albedo is an important variable in climate modeling studies and satellite-based retrievals of land-surface properties. In this study, we used multiyear 500 m albedo products from the Moderate Resolution Imaging Spectroradiometer (MODIS) to derive the bare soil albedo for seven spectral bands and three broadbands over the contiguous United States (CONUS). The soil line based on red and green spectral signatures derived from MODIS data was used as the basis to detect and extract bare soil albedo. A comparison against bare soil albedo derived from 30 m Landsat data has been made, showing that the MODIS bare soil albedo had a bias of 0.003 and a root-mean-square-error (RMSE) of 0.036. We found that the bare soil albedo was negatively correlated with soil moisture from the Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E), with a relatively stable exponential relationship reflecting the darkening effect that moisture has on most soils. However, quantification of the relationship between bare soil albedo and soil moisture still needs to be improved through simultaneous and instantaneous measurements at a finer spatial resolution. Statistics of the multiyear climatological bare soil albedos calculated using soil types and the International Geosphere-Biosphere Programme (IGBP) land cover types suggest that: Land cover type is a better indicator for determining the magnitude of bare soil albedos for the vegetated areas, as the vegetation density is correlated with soil moisture; and soil type is a better indicator for determining the slope of soil lines over sparsely vegetated areas, as it contains information of the soil texture, roughness, and composition. The generated bare soil albedo can be applied to improve the parameterization of surface energy budget in climate and remote sensing models as well as the retrieval accuracy of some satellite products.