Climate data records of vegetation variables from geostationary SEVIRI / MSG data : products , algorithms and applications

The uncertainties due to the BRDF parameters correspond with the diagonal elements of the uncertainty covariance matrix (C00, C11 and C22 fields of the AL-Ci-CK product) obtained for the three model parameters (SAF/LAND/MF/PUM_AL/1.4). The absolute uncertainty of FVC takes into account the uncertainty due to the propagation of the input errors (σk0 ) and model selection (σmodel): σFVC = √σk0 2 + σmodel 2 . (1)


Introduction
Observations and monitoring are essential components of Global Framework for Climate Services initiative [1,2] that contribute significantly to meet the needs of climate research services.The role that satellites have played in observing the variability and change of the Earth system has increased significantly in the last decades.Thanks to the Earth observation (EO) satellite systems, a large number of variables related to the atmosphere, oceanic and terrestrial domains are accessible [3].
The Satellite Application Facility (SAF) for Land Surface Analysis (LSA) is a leading center for development and operational retrieval of land surface variables from European Organization for the Exploitation of Meteorological Satellites (EUMETSAT), specifically aimed to understand and quantify terrestrial processes and land-atmosphere interactions.Land surface models' assessment and improvement of high quality and robust EO products are among the target applications of LSA SAF products.Particularly, the geostationary Spinning Enhanced Visible and Infrared Imager (SEVIRI) on board the Meteosat Second Generation (MSG) platform provides a very high temporal resolution (15 minutes) required for resolving the diurnal cycles of radiation, surface albedo and surface temperature [4].The MSG program covers a series of four identical satellites MSG 1-4 (Meteosat 8-11) launched in 2002,2005,2012, and 2015 respectively.SEVIRI operates in 12 spectral channels [5], although only channels 1 (VIS0.6), 2 (VIS0.8)and 3 (NIR1.6),centered at about 0.635 µm (red), 0.81 µm (NIR) and 1.64 µm (SWIR), respectively situate in the optical domain.
The LSA SAF provides also variables for the characterization of terrestrial ecosystems and their role in the energy balance of Earth, such as land surface fluxes and three important biophysical variables related with the amount, structure and state of vegetation: Leaf Area Index (LAI), Fraction of Absorbed Photosynthetically Active Radiation (FAPAR) and Fractional Vegetation Cover (FVC) [6].These variables are widely used in many land-surface vegetation, climate, and crop production models [7][8][9][10][11][12][13]. FVC represents the fraction of green vegetation covering a unit area of horizontal soil, corresponding to the gap fraction in the nadir direction.LAI is a measure of the amount of live foliage present in the canopy per unit ground surface.The FAPAR is the fraction of photosynthetically active radiation (400-700 nm) absorbed by the green parts of the canopy, and therefore constitutes an indicator of the presence, health and productivity of live vegetation.
Nowadays, consistent long-term data records of these variables with well-characterized uncertainty are currently required by the scientific community to model terrestrial ecosystems and energy cycles at regional and global scales [14].Models related with land surface and climate usually require Near-Real-Time (NRT) datasets of these variables at high temporal frequencies, and with regional to global coverage [15].Both, LAI and FAPAR are recognized as essential climate variables in the terrestrial domain by Global Climate Observing System (GCOS) [16].
Several methods have been proposed for estimating biophysical parameters from large scale optical sensors, some of which have been implemented in operational processing lines from POLarization and Directionality of Earth Reflectances (POLDER) [17], Moderate Resolution Imaging Spectroradiometer (MODIS) [18], MISR (Multi-angle Imaging SpectroRadiometer) [19], MEdium Resolution Imaging Spectrometer (MERIS) [20,21], Sea-viewing Wide Field-of-view Sensor (SEAWIFS) [22], VEGETATION [23,24] and Advanced Very High Resolution Radiometer (AVHRR) [25].These methods can be roughly classified in two main categories: (i) methods based on optimized vegetation indices [26], and (ii) methods based on the inversion of physical models.Both approaches present advantages and drawbacks.It is referred to [27] for a detail intercomparison on the spread of currently available physical models.Although methods that invert physical models can potentially be applied to varying surface types, they require specifying a large number of model parameters, which may lead to unstable solutions [28].Most inversion methods use one-dimensional (turbid) models (e.g., PROSAIL [29]), thereby assuming surface homogeneity within an image pixel.They are thus limited to address the sub-pixel mixing of land cover types, which is common at coarse resolutions satellites.The process of identifying the sub-pixel proportions of the constituent components is called spectral mixture model (SMM).SMM approaches are methods especially adequate for global studies, since the spatial variability within pixel is high [30,31].
In this paper, we describe the methodology currently used in LSA SAF to operationally generate the suite of SEVIRI/MSG (FVC, LAI and FAPAR) vegetation products.The high frequency of acquisition provided by the SEVIRI instrument assures the availability of cloud-free data for appropriately monitoring the vegetation dynamics in large regions.When one is choosing a retrieval method, one needs to consider the limited spatial and spectral resolution of SEVIRI instrument (i.e., only 3 optical channels and a spatial resolution of 3 km at the sub-satellite point that declines rapidly at higher latitudes).Furthermore, the challenges that the MSG geostationary satellite poses as compared to polar-orbiting satellite sensors due to diurnal variations in the solar zenith angle need to be considered.The algorithm to retrieve FVC relies on a novel stochastic SMM method, which addresses the large variability of soils and vegetation types using statistical (sum of Gaussian) distributions.The LSA SAF operational scheme to retrieve LAI and FAPAR uses statistical relationships general enough for global applications, proposed by Roujean and Lacaze [16] and Roujean and Bréon [32], respectively.
LSA SAF has been providing in NRT vegetation products (FVC, LAI and FAPAR) at two different time resolutions (daily and 10-days) since 2008 over the geostationary Meteosat disk, covering Europe, Africa, the Middle East and parts of South America.The suite of SEVIRI/MSG vegetation products is freely disseminated to users through LSA SAF website (https://landsaf.ipma.pt).During the different operational phases, these vegetation products have been routinely assessed and some improvements have been made in the upgraded versions.LSA SAF has recently reprocessed the entire MSG archive with the latest version of the several retrieval algorithms in order to obtain a continuous and homogeneous Climate Data Records (CDR) of vegetation products.More than fifteen years (January 2004-present) of homogeneous consistent estimates of FVC, LAI and FAPAR fields are available.The 10-days (FVC, LAI and FAPAR) CDRs are part of the suite of EUMETSAT climate products, which include Fundamental CDRs, Thematic CDRs and Data for operational climate monitoring and are also available from EUMETSAT web page (https://navigator.eumetsat.int/).
The paper presents the following structure.Section 2 is dedicated to describe the algorithms to retrieve FVC, LAI and FAPAR.The suite of LSA SAF SEVIRI/MSG vegetation products are described in Section 3, including expert knowledge and a quality assessment of the consistency between the products.Section 4 provides further insight about the added value and introduces examples of its application in numerous environmental applications.

SEVIRI/MSG
The algorithm for retrieving biophysical variables from SEVIRI/MSG data uses as input normalized spectral bidirectional reflectance factor (BRDF) at a fixed angular configuration (i.e., vertical illumination and observation angles of 0 • ) in order to reduce angular effects.This hotspot geometry leads to a minimum contribution of the shadow proportion and a precise estimation of FVC, coinciding with the complement to unity of the gap fraction at nadir direction [16].Estimating the FVC with increased sun zenith angles would lead to an overestimation of FVC.The main drawback of this geometry is the significant contribution of illuminated soil background, which constitutes a confounding factor for the retrieval of FVC. Figure 1 outlines the main ingredients for deriving the SEVIRI/MSG products including the stochastic SMM approach for retrieving FVC, a pragmatic method for LAI retrieval from FVC, and a method based on an optimized vegetation index for FAPAR determination.

FVC Algorithm
Let r be the spectrum of each mixture pixel, i.e., a column vector (r 1 ,r 2 ,...,r n ), where n is the total number of bands.When multiple scattering can be reasonably disregarded, the spectral reflectance of each pixel can be approximated by a linear mixture of pure spectra E (the so-called endmembers) weighted by their corresponding fractional abundances: where E [n × c] is the matrix of endmembers, f is a vector with the c unknown proportions in the mixture, and ε is the residual vector.The mixing equation is accompanied by two constraints: (1) the normalization constraint says proportions should sum up to one, and (2) the positivity constraint says that no endmember can make a negative contribution.The least-square principle establishes that the unknown parameters are those that minimize the Mahalanobis distance between the pixel r and point E f: where V(r) denotes the error matrix of the observations r.The coarse spatial resolution of SEVIRI and the availability of only three spectral bands poses a significant challenge for endmember selection in traditional (deterministic) SMM.The FVC algorithm relies on a novel stochastic SMM which uses a statistical representation of the endmembers to accommodate their variability at a global scale.Each pixel is described by a multiple combination of two pure classes, the target (vegetation) and the background (soil).Both vegetation and soil classes are not treated as deterministic (i.e., using fixed endmembers) but they follow a statistical (multi-modal) distribution, which attempts to capture the variability of soils and vegetation components at a global scale.The main steps of the algorithm are now described (see Figure 2).The BRDF algorithm applies a semi-empirical reflectance model in order to invert SEVIRI/MSG top-of-canopy reflectance factor values into three parameters (k 0 , k 1 , k 2 ).The input of the FVC algorithm is atmospherically corrected BRDF k 0 parameters in the MSG channels 1, 2 and 3, while FAPAR uses all three BRDF parameters (k 0 , k 1 , k 2 ) in order to correct surface's reflectance anisotropy and minimize the effect of soil reflectance (see Figure 1).The algorithm to correct atmospheric effects and generate BRDF model parameters is fully described in [33], and in the ATBD of the Albedo MSG product [34].

FVC Algorithm
Let r be the spectrum of each mixture pixel, i.e., a column vector (r 1 ,r 2 ,...,r n ), where n is the total number of bands.When multiple scattering can be reasonably disregarded, the spectral reflectance of each pixel can be approximated by a linear mixture of pure spectra E (the so-called endmembers) weighted by their corresponding fractional abundances: where E [n × c] is the matrix of endmembers, f is a vector with the c unknown proportions in the mixture, and ε is the residual vector.The mixing equation is accompanied by two constraints: (1) the normalization constraint says proportions should sum up to one, and (2) the positivity constraint says that no endmember can make a negative contribution.The least-square principle establishes that the unknown parameters are those that minimize the Mahalanobis distance between the pixel r and point E f: where V(r) denotes the error matrix of the observations r.The coarse spatial resolution of SEVIRI and the availability of only three spectral bands poses a significant challenge for endmember selection in traditional (deterministic) SMM.The FVC algorithm relies on a novel stochastic SMM which uses a statistical representation of the endmembers to accommodate their variability at a global scale.Each pixel is described by a multiple combination of two pure classes, the target (vegetation) and the background (soil).Both vegetation and soil classes are not treated as deterministic (i.e., using fixed endmembers) but they follow a statistical (multi-modal) distribution, which attempts to capture the variability of soils and vegetation components at a global scale.The main steps of the algorithm are now described (see Figure 2).Step 1: Selection of a training database Given the large diversity of vegetation types and underlying soil background that can be found in the SEVIRI disk, a large training database of spectral signatures needs to be generated.To more reliably identify the pure components, two composite k0 images were generated from all high quality observations (cloud-and snow-free) over a one-year period: a vegetated k0 image corresponding to the peak of season and a devegetated k0 image corresponding to the minimum canopy closure.The samples were chosen to be homogeneous over areas higher than SEVIRI spatial resolution, and the classification is based on the 1-km Global Land Cover 2000 (GLC2000) [35].Training areas of nonvegetated class include desert areas, sparsely vegetated and shrublands.Vegetation areas were identified in crop and forest classes.Samples were further verified using purity methods to filter out possible outliers.In particular, pixels having less than 95% of soil/vegetation were excluded.
Step 2: Gaussian mixture model At the SEVIRI resolution, the soils and vegetation components usually present a high sub-pixel spatial variability.In order to address this variability, a parametric mixture model weighted sum of Gaussian distributions (or clusters) is adopted [36].This probabilistic model is applied separately for the target (vegetation) and the background (soil) classes: where  is class conditional distribution of the k-th Gaussian, with mean  and covariance  , is the probability that an observation belongs to the k-th component ( 0 ; ∑  = 1 ), and Gs and Gv are the number of Gaussian components for soil and vegetation, respectively.The algorithm uses the Expectation-Maximization (E-M) approach [37] to estimate the means  and covariances  of the individual Gaussian components.This simple iterative approach increases the log likelihood of the observed data and usually converges if the data conform reasonably well to the mixture model.The k-means algorithm is used to initialize the E-M parameters ( ,  ).Multiple initializations are made to avoid numerical problems local maxima.We assume ellipsoidal unconstrained component covariance matrices and use the Bayesian Information Criterion (BIC) [38,39] to determine an appropriate number of Gaussian components.The value of BIC is the Step 1: Selection of a training database Given the large diversity of vegetation types and underlying soil background that can be found in the SEVIRI disk, a large training database of spectral signatures needs to be generated.To more reliably identify the pure components, two composite k 0 images were generated from all high quality observations (cloud-and snow-free) over a one-year period: a vegetated k 0 image corresponding to the peak of season and a devegetated k 0 image corresponding to the minimum canopy closure.The samples were chosen to be homogeneous over areas higher than SEVIRI spatial resolution, and the classification is based on the 1-km Global Land Cover 2000 (GLC2000) [35].Training areas of non-vegetated class include desert areas, sparsely vegetated and shrublands.Vegetation areas were identified in crop and forest classes.Samples were further verified using purity methods to filter out possible outliers.In particular, pixels having less than 95% of soil/vegetation were excluded.
Step 2: Gaussian mixture model At the SEVIRI resolution, the soils and vegetation components usually present a high sub-pixel spatial variability.In order to address this variability, a parametric mixture model weighted sum of Gaussian distributions (or clusters) is adopted [36].This probabilistic model is applied separately for the target (vegetation) and the background (soil) classes: where φ k is class conditional distribution of the k-th Gaussian, with mean µ k and covariance Σ k , τ k is the probability that an observation belongs to the k-th component (τ k ≥ 0 ; and G s and G v are the number of Gaussian components for soil and vegetation, respectively.The algorithm uses the Expectation-Maximization (E-M) approach [37] to estimate the means µ k and covariances Σ k of the individual Gaussian components.This simple iterative approach increases the log likelihood of the observed data and usually converges if the data conform reasonably well to the mixture model.The k-means algorithm is used to initialize the E-M parameters (µ k , Σ k ) .Multiple initializations are made to avoid numerical problems local maxima.We assume ellipsoidal unconstrained component covariance matrices and use the Bayesian Information Criterion (BIC) [38,39] to determine an appropriate number of Gaussian components.The value of BIC is the maximized log likelihood with a penalty for the number of parameters.The larger the BIC score, the stronger the evidence for the model.The determination of the number of Gaussian components takes into account not only the BIC score but also the requirement that the distributions must provide a faithful representation of the data.A typical number of 6-7 Gaussians for soil and 3-5 for vegetation has been used to represent the variability of the different SEVIRI geographical areas.
Step 3: Model selection Each SEVIRI pixel is usually a mixture of several backgrounds and/or vegetation components.The stochastic SMM computes all possible models by taking all possible sets of models {M 1 , . . .
. At the SEVIRI resolution, different sets of soil and vegetation may yield to very similar mixtures.Let p(M k |r) be the posterior probability or likelihood of model M k given pixel data r, and π(M K ) the a priori probability of having the model M k at a particular pixel.Based on the Bayes theorem most often used in classification problems, the posterior probability assigned to a model M K is proportional to its likelihood times its prior probability: Although a priori probabilities π(M k ) are often unknown and could be assumed to be equal (uninformative prior), they provide a means to inject prior information in the algorithm and, therefore, reduce the misidentification of the models by considering rules on the basis of ancillary data such as land cover classifications.The likelihood of M k for each individual pixel, p(r|M k ), can be determined as follows: where Γ(x, x , r) is an "accept-reject" function that is the unity when the line joining x (vegetation) and x' (soil) intercepts the region in the feature space centered at the mixture r, and is zero otherwise.We considered around r an envelope volume V(r) given by the radiometric uncertainties attached to the input, i.e., the covariance structure of the k 0 product.The algorithm solves this integral using Monte Carlo methods, i.e., from endmember spectra drawn from multivariate normal distribution.It basically consists in counting the number of model samples from each possible mixture combination that can give rise to the mixture pixel r, as illustrated in Figure 3 considering a simplified bidimensional problem with 3 vegetation clusters (V 1 , V 2 , V 3 ) and 3 soil clusters (S 1 , S 2 , S 3 ).The size and orientation in the elliptical probability density contours are determined by the covariance of clusters Σ k .It is obvious that (S 3 , V 1 ) is the most likely model in the mixture.However, other different model combinations, such as (S 2 , V 1 ) and (S 3 , V 2 ), have non negligible probabilities and may yield to the same mixture spectrum.
drawn from multivariate normal distribution.It basically consists in counting the number of model samples from each possible mixture combination that can give rise to the mixture pixel r, as illustrated in Figure 3    Rather than using a single date, multi-temporal SEVIRI observations {t 1 , t 2 , . . .t M } over a seasonal period can be used to more reliably determine the likelihood of model M K in a mixture pixel r, i.e.,: where r(t i ) represent the input reflectance at date t i .The LSA SAF algorithm uses M=2 dates, associated to the maximum and minimum canopy closure of each SEVIRI pixel, making use of the vegetated/devegetated k 0 images generated in step 1: The simultaneous use of devegetated k 0 t devegetated and vegetated k 0 t vegetated spectral information for extreme situations of the annual growing cycle such as harvested crops and peak of season, are well suited to determine the likelihood of the soil and vegetation components, respectively.Figure 4 allows to better understand the concept.Symbols correspond to pure soil and vegetation samples identified over the Southern Africa SEVIRI geographical area, projected onto the bidimensional k 0 space of SEVIRI channels.Bare areas are predominantly found in sparsely vegetated and open shrublands (GLC2000 classes 19, 14 and 12) whereas purely vegetated areas are mostly found in close forest classes, herbaceous and croplands (GLC2000 classes 1, 2, 13 and 16).In this example, the best suited number of Gaussian components was 3 for vegetation and 6 for soil.associated to the maximum and minimum canopy closure of each SEVIRI pixel, making use of the vegetated/devegetated k0 images generated in step 1: The improved multitemporal identification of vegetation and underlying soil background is illustrated for a SEVIRI pixel occupied by open broadleaved deciduous species.Blue symbols correspond to the migration of the pixel signature during year 2015 (numbers from 1 to 36 correspond to decades).The leaf fall is concentrated during the winter season, reaching devegetated period (filled pentagram in Figure 4) at decade 8, when the spectral signature approaches to the soil line.During the development phase, NIR values increase while red and MIR values decrease, and the pixel departs from the soil line, eventually reaching a high canopy closure (filled square in Figure 4a).Combining the probabilities in both extreme periods (Equation ( 9)), the most likely models for the vegetation and soil are obtained, i.e., vegetation classes V 1 and V 2 with the darkest soil type S 1 .
Step 4. Estimation of FVC The retrieval of vegetation parameters is an ill-posed problem.A feature selection/extraction step, may be useful to transform the original space r ∈ •3 [(k 0 ) red , (k 0 ) NIR , (k 0 ) SWIR ] onto a feature space that makes the retrieval model more sensitive to changes in vegetation and minimize the sensitivity to soil background and noise in the inputs.Vegetation retrieval is specially hampered by the influence of undesired soil variability at SWIR wavelengths (1.6 µm), which is particularly large for the soils in Africa (see Figure 4b).This high variability has shown to be a cause for overestimation of FVC in semi-arid regions over dark soils (further insights will be given in Figure 5b).Latest version projects r onto a new features space w ∈ •5 , increasing the relative influence of channels 1 and 2 with respect to channel 3 ([(k 0 ) red , (k 0 ) red , (k 0 ) NIR , (k 0 ) NIR , (k 0 ) SWIR ]).This strategy has served us to reduce possible biases due to SWIR background reflectance variability.
The results correspond to mean predictions of FVC given a realized value of two SEVIRI channels (considering all possible values of the third channel).We can identify the triangle envelope formed by the mixed pixels which is the physical basis for SMM approaches.The soil/vegetation components are spectrally distinct enough, particularly in the red-NIR space (Figure 5a), where non-vegetation surfaces distribute primarily along the so-called "soil line".Dense canopies situate in the top vertex, although showing a significant spectral variability.The LSA SAF algorithm provides a fair generalization of information from the training data generated in step 1, suggesting that solutions are well constrained.The isolines present realistic values that are consistent with other literature studies [46], and EPS satellite products [24].The smooth FVC variations with gradual change in reflectance suggests that solutions are stable, without discontinuities that lead to large retrieval errors and are usually found in SMM problems with a large number of endmembers.The algorithm mitigates thus the effects of noise in the input data.
The algorithm uses also SWIR information, which conveys useful information about the vegetation water content and a good sensitivity to FVC (Figure 5b).However, the large variability of SWIR reflectance of soils in the African continent has shown to be confounding factor.In particular, inclusion of SWIR bands causes an overestimation trend for FVC in sparse areas over dark soils in Southern Africa, with values up to 0.15 (see Figure 5b).An effective way to reduce this overestimation is underweighting the SWIR band (step 4 of the algorithm).

LAI Algorithm
The LSA SAF algorithm uses a pragmatic solution to the radiative transfer problem, which assumes a tractable physical model for interception of solar direct irradiance by leaf canopies [43].By assuming that leaves are flat with bi-Lambertian properties, the total transmittance, which represents the fraction of incident radiation above the canopy which reaches the soil background level can be expressed as follows [47]: In addition, a standardization transform is applied, which transforms the data to a set of variations about the mean value with a mean value of zero and a standard deviation of one: where r is the standardized vector associated to the pixel vector r, with mean and µ r standard deviation σ r .Using the standardised endmembers, Êi (i = 1, . . .,c), the unmixing is formulated as follows: where fi is the proportion of such endmember in the standardised coordinates, and ε is the residual vector expressed in standardised units.Although the solution is similar to the conventional SMM, the sum-to-one condition is now expressed as fi Through this standardization, SMM is less sensitive to the brightness variability within each vegetation-soil component [40].Standardized signatures are performed for all possible models and image spectra as previous step before computing the fractions of soil and vegetation.In order to compute the fractions, the closed formulation of the method of Lagrange is adopted [36], which provides a fast and unbiased solution.Finally, FVC is estimated as a linear combination of single-model estimates: where FVC(M k ) is unmixed abundance of vegetation in the M k model.The contribution of each model is thus weighted by its a posteriori probability p(M K |r).Although other similar approaches have been developed to address the issue of the spatial variability of endmembers [40][41][42][43][44][45], the LSA SAF approach decomposes the soil/vegetation into a number of subclasses using non-uniform probabilities for the different models.These probabilities p(M K |r) vary from pixel to pixel, depending on its likelihood and were pre-computed beforehand in order to increase the robustness of the algorithm and speed up its computational efficiency.Figure 5 allows understanding better the physical basis of the algorithm and its performance.The results correspond to mean predictions of FVC given a realized value of two SEVIRI channels (considering all possible values of the third channel).We can identify the triangle envelope formed by the mixed pixels which is the physical basis for SMM approaches.The soil/vegetation components are spectrally distinct enough, particularly in the red-NIR space (Figure 5a), where non-vegetation surfaces distribute primarily along the so-called "soil line".Dense canopies situate in the top vertex, although showing a significant spectral variability.The LSA SAF algorithm provides a fair generalization of information from the training data generated in step 1, suggesting that solutions are well constrained.The isolines present realistic values that are consistent with other literature studies [46], and EPS satellite products [24].The smooth FVC variations with gradual change in reflectance suggests that solutions are stable, without discontinuities that lead to large retrieval errors and are usually found in SMM problems with a large number of endmembers.The algorithm mitigates thus the effects of noise in the input data.
The algorithm uses also SWIR information, which conveys useful information about the vegetation water content and a good sensitivity to FVC (Figure 5b).However, the large variability of SWIR reflectance of soils in the African continent has shown to be confounding factor.In particular, inclusion of SWIR bands causes an overestimation trend for FVC in sparse areas over dark soils in Southern Africa, with values up to 0.15 (see Figure 5b).An effective way to reduce this overestimation is underweighting the SWIR band (step 4 of the algorithm).

LAI Algorithm
The LSA SAF algorithm uses a pragmatic solution to the radiative transfer problem, which assumes a tractable physical model for interception of solar direct irradiance by leaf canopies [43].By assuming that leaves are flat with bi-Lambertian properties, the total transmittance, which represents the fraction of incident radiation above the canopy which reaches the soil background level can be expressed as follows [47]: where µ s = cosθ s , being θ s the solar zenith angle, G(θ s ) is the average extinction function [48], T(θ s ) and b is the backscattered parameter, which can be roughly assumed to be equal to 0.945 for all vegetation types [49].The fraction of solar radiation intercepted by the vegetation (FIPAR), which coincides with FVC when the sun and the observer are both at zenith, i.e., FVC = FIPAR(θ s = 0) is expressed as: The above formulation assumes a random distribution of foliage elements, disregarding the increased probability of light penetration in clumped canopies.A clumping index Ω [50], which accounts for the degree of the deviation of foliage distribution from the random case, can be introduced in Equation (15) to correct possible underestimation of LAI [20]: It can be assumed for simplicity a value of 0.5 for the leaf projection factor G(θ s ) considering spherical orientation of the foliage.In order to avoid maximum LAI values in fully vegetated areas exceeding a value about 7, a coefficient a 0 in the range (1.04-1.07) is introduced in Equation ( 16): A land cover-dependent clumping index is considered for each of the GLC2000 classes, which corresponds to the maximum values calculated from global POLDER multiangular data [51].The values range from 0.68 for evergreen forest to 0.83-0.85for herbaceous, shrub and cultivated areas (see Table S1 in Supplementary Material), leading to a conservative first-order correction of the clumping effect.

FAPAR Algorithm
The algorithm used for retrieving daily integrated green FAPAR is not novel but it relies on the method proposed by [32], a statistical relationship general enough for global applications based on simulations using the homogeneous SAIL (Scattering by Arbitrary Inclined Leaves) model [52].The algorithm does not use any prior knowledge on the land cover.Although the SAIL model is widespread in the remote sensing community for the estimation of vegetation biophysical variables, limitations of homogeneous models should be more important in heterogeneous (e.g., savannas and open shrublands) or in canopies showing a complex architecture (e.g., boreal forest).More than 5000 soil/vegetation combinations varying leaf inclination distribution (LIFD), LAI, leaf transmittance, leaf reflectance and soil spectral albedo were considered.The diffuse fraction of incoming radiation was held constant and equal to 0.2, which represents clear sky conditions.For each scenario, red and NIR reflectances with variations of sun and view angles were derived from the SAIL model along with daily-integrated FAPAR, as computed by integration of the instantaneous FAPAR over the day: where t and t' are the time for sunrise and sunset.The FAPAR was integrated over solar angles corresponding to a target located at 45 • N latitude and at the equinox.An optimal geometry based on the criteria of linearity and minimum dispersion between NDVI and daily-integrated FAPAR was found in the solar principal plane (θ s =45 • , θ v =60 • , φ=0 • ).The soil contribution to the canopy reflectance was reduced at this oblique view using a vegetation index, called RDVI (Renormalized Difference Vegetation Index), defined as follows: where DVI [16] is the difference vegetation index, defined as follows: Finally, the RDVI-FAPAR relationship in the optimal geometry is given as follows: where (RDVI) opt refers to the RDVI computed in the optimal geometry.The reflectance in the optimal geometry for each spectral channel is estimated as follows [53]:

Products Uncertainty Estimation
A per-pixel estimate of the uncertainty for each prediction, namely σ FVC , σ LAI and σ FAPAR , is provided, which propagates uncertainty of inputs and other variables involved in the calculations (see Equations in Supplementary Material).

The SEVIRI/MSG Vegetation Products
The LSA SAF generates and disseminates a suite of vegetation products derived from SEVIRI/MSG BRDF composited data for the whole Meteosat disk over land surfaces ideally free of snow and ice cover using the described algorithm at two different time resolutions (see Figure 6), daily (MDLAI, MDFVC, MDFAPAR), and 10-days (MTLAI, MTFVC, MTFAPAR) [54].The characteristics of SEVIRI based products provided by the LSA SAF are described in Table 1.Daily products use as input the MDAL BRDF parameters computed using a daily rolling compositing approach with a characteristic 5-day compositing period, whereas the 10-days product use as input the MTAL BRDF product based on a composite over a 30-day period.We refer for further details about the recursive temporal composition to Geiger et al. [33].Using the latest version of the daily (v3.1) and 10-days (v1.2) algorithms, released in January 2016, the LSA SAF system generates and disseminated the products in Near-Real-Time (NRT) with a time lag of about six hours.NRT users are also encouraged to use EUMETCast (https://navigator.eumetsat.int),i.e., the primary distribution means for EUMETSAT image data and derived products.
field.The error dataset is a pixel-wise estimate of the observation uncertainty, which is indicative of the quality of the retrieval error (see details of calculation in Supplementary Material).These retrieval errors have shown to be within the range of the typical differences found between satellite products.For examples the uncertainty ranges typically between 0.05 and 0.10 for FVC, between 0.5 and 1.0 for LAI and between 0.05 and 0.20 for FAPAR.The FAPAR product presents a lower quality because it uses as input the three BRDF parameters (k0, k1, k2), with k2 presenting large uncertainties and noisy profiles on a short time scale, mainly in Western Africa.The algorithms to retrieve LSA SAF vegetation require ideally free of snow and ice cover observations.Since snow events in previous days, though influencing the signal, may be unidentified BRDF product, an empirical criterion was used to identify and mask out pixels with residual, by combining the SEVIRI channels 1 (0.6 μm) and 3 (1.6 μm):  ( ) −  ( ) > 0   ( ) >  , ( ) + 0.06  The algorithms were reprocessed to generate the full archive of CDRs for 10-days vegetation products for FVC (MTFVC-R), LAI (MTLAI-R) and FAPAR (MTFAPAR-R) over the period 2004-2015, which complement the datasets available in the LSA SAF website since 2016 onwards.Both NRT products and homogeneous CDR 10-days vegetation products are available from the website (https: //landsaf.ipma.pt/en/products/vegetation/).A similar dataset of reprocessed daily products has been produced as internal products for the period 2004-2015.Although the daily products are not foreseen to be distributed as CDR, they may be made available upon request.
The algorithm is firstly applied separately at four SEVIRI geographic regions: Europe (Euro), Northern Africa (NAfr), Southern Africa (SAfr), and South America (SAme) (see location in Figure 6).However, since the latest version (v3.1) onwards, products are disseminated for the whole Meteosat disk.The projection and spatial resolution correspond to the characteristics of Level 1.5 (standard geolocation) MSG/SEVIRI instrument data.Information on geo-location and data distribution is available at the LSA SAF web-site: http://landsaf.ipma.pt.
Each product is delivered in the Hierarchical Data Format version 5 (HDF5), which contains three separate datasets: a biophysical product, its respective error estimate and a Quality Flag (QF) field.The error dataset is a pixel-wise estimate of the observation uncertainty, which is indicative of the quality of the retrieval error (see details of calculation in Supplementary Material).These retrieval errors have shown to be within the range of the typical differences found between satellite products.For examples the uncertainty ranges typically between 0.05 and 0.10 for FVC, between 0.5 and 1.0 for LAI and between 0.05 and 0.20 for FAPAR.The FAPAR product presents a lower quality because it uses as input the three BRDF parameters (k 0 , k 1 , k 2 ), with k 2 presenting large uncertainties and noisy profiles on a short time scale, mainly in Western Africa.
The algorithms to retrieve LSA SAF vegetation require ideally free of snow and ice cover observations.Since snow events in previous days, though influencing the signal, may be unidentified BRDF product, an empirical criterion was used to identify and mask out pixels with residual, by combining the SEVIRI channels 1 (0.6 µm) and 3 (1.6 µm): where k 0,max (λ 1 ) and k 0,min (λ 3 ) refer to the values of the k 0 devegetated image in channel 1 and 3, respectively.The QF records relevant information at the pixel level about the quality and processing status of retrievals (e.g., unrealistic input ranges, traces of inland water, traces of snow and retrieval values out of maximum physical range).The QF codes and several empirical thresholds used to blind problematic areas are shown in Tables S2 and S3 in Supplementary Material, respectively, as in the Product User Manual (available from http://landsaf.ipma.pt).
Figure 7 shows quality information of the SEVIRI/MSG FVC, LAI and FAPAR products, using the mean error of the products along the year 2014.Green color refers to consolidate regions with optimal quality.Reliable areas with medium quality are shown in cyan color.Orange color refers to regions with low quality.Finally, red color refers to generally unusable areas, presenting large view zenith angles and frequent snow cover (e.g., Europe during wintertime).FAPAR is generally unusable in areas with persistent cloud occurrence (e.g., western Africa).The seasonal variations in the quality and coverage of the daily MSG FVC product along 2014 are depicted in Figure 8. Optimal or medium quality retrievals are obtained for north and south Africa continental zones with a negligible percentage of poor quality or missing values.An acceptable performance is also found over Europe from April through September, with a significant quality reduction during late autumn and winter.The worst performance is observed for south America regions, generally due to decreased accuracy and a larger percentage of missing values.Both daily MSG LAI and FAPAR products show similar performance stages with slightly lower rates of high quality values.The seasonal variations in the quality and coverage of the daily MSG FVC product along 2014 are depicted in Figure 8. Optimal or medium quality retrievals are obtained for north and south Africa continental zones with a negligible percentage of poor quality or missing values.An acceptable performance is also found over Europe from April through September, with a significant quality reduction during late autumn and winter.The worst performance is observed for south America regions, generally due to decreased accuracy and a larger percentage of missing values.Both daily MSG LAI and FAPAR products show similar performance stages with slightly lower rates of high quality values.A detailed assessment of the MSG CDR products is provided in a validation report [55], showing good spatial and temporal consistency as compared to currently available validated vegetation products and no presence of artefacts.Noteworthy good inter-annual precision and stability of the long-term time series was found.Based on the validation results, performed over a representative network of sites over global conditions, CDR MSG vegetation products reached the operational status, what means that the products fulfilled all defined requirements and are suitable for distribution to users.Positive results were found for all the considered criteria, evaluated by comparisons with similar products such as Copernicus Global Land Monitoring Service GEOV1 based on SPOT/VGT data and NASA MODIS/TERRA products.Accuracy assessment against ground references provided by the Committee on Earth Observation Satellite (CEOS) through its Working Group in Calibration and Validation (WGCV) showed improved accuracy as compared to similar products and up to 84% of validation samples within target requirements [55].

Internal Consistency between the LSA SAF Products
Temporal profiles of the 10-days MSG vegetation products were verified through a comparison with the daily vegetation products.A very good agreement was found for the three products (see examples in Figure 9).Although one main drawback of the daily FAPAR product with regard to other operational products was the shaky temporal profiles on a daily basis, the 10-days (MTFAPAR) improves notably the smoothness of daily FAPAR.A detailed assessment of the MSG CDR products is provided in a validation report [55], showing good spatial and temporal consistency as compared to currently available validated vegetation products and no presence of artefacts.Noteworthy good inter-annual precision and stability of the long-term time series was found.Based on the validation results, performed over a representative network of sites over global conditions, CDR MSG vegetation products reached the operational status, what means that the products fulfilled all defined requirements and are suitable for distribution to users.Positive results were found for all the considered criteria, evaluated by comparisons with similar products such as Copernicus Global Land Monitoring Service GEOV1 based on SPOT/VGT data and NASA MODIS/TERRA products.Accuracy assessment against ground references provided by the Committee on Earth Observation Satellite (CEOS) through its Working Group in Calibration and Validation (WGCV) showed improved accuracy as compared to similar products and up to 84% of validation samples within target requirements [55].

Internal Consistency between the LSA SAF Products
Temporal profiles of the 10-days MSG vegetation products were verified through a comparison with the daily vegetation products.A very good agreement was found for the three products (see examples in Figure 9).Although one main drawback of the daily FAPAR product with regard to other operational products was the shaky temporal profiles on a daily basis, the 10-days (MTFAPAR) improves notably the smoothness of daily FAPAR.The aim of this section is also to analyse the consistency among the suite of LSA SAF vegetation products (i.e., FAPAR and LAI) and albedo fields.The energy absorbed by the ground below vegetation (FGROUND) can be estimated by combining the FAPAR and an independent LSA SAF field from the albedo SEVIRI product, the bi-hemispherical reflectance integrated over the photosynthetically active spectral region (BHRPAR) [33].From energy conservation, the fraction of PAR absorbed by the ground beneath the canopy is given by [56]: Following previous works [56], FGROUND can be expressed as the downward PAR flux density, F, times the soil absorptance, 1−α, where α is the reflectance of the canopy background and F can be expressed via the Beer's law [57]: where r* is the probability that photon entering through the lower canopy boundary will be reflected back by the vegetated layer [58,59].Figure 10 shows the relationship between MSG LAI and FGROUND for two SEVIRI regions, Europe (Euro), and southern Africa (SAfr), as derived from BHRPAR and FAPAR SEVIRI/MSG products using Equation (24).The distribution can be well approximated by the exponential function, FGROUND = Ae • for different periods, which indicates a strong consistency between LAI and FGROUND at the resolution of SEVIRI products.The aim of this section is also to analyse the consistency among the suite of LSA SAF vegetation products (i.e., FAPAR and LAI) and albedo fields.The energy absorbed by the ground below vegetation (FGROUND) can be estimated by combining the FAPAR and an independent LSA SAF field from the albedo SEVIRI product, the bi-hemispherical reflectance integrated over the photosynthetically active spectral region (BHRPAR) [33].From energy conservation, the fraction of PAR absorbed by the ground beneath the canopy is given by [56]: Following previous works [56], FGROUND can be expressed as the downward PAR flux density, F, times the soil absorptance, 1−α, where α is the reflectance of the canopy background and F can be expressed via the Beer's law [57]: where r* is the probability that photon entering through the lower canopy boundary will be reflected back by the vegetated layer [58,59].Figure 10 shows the relationship between MSG LAI and FGROUND for two SEVIRI regions, Europe (Euro), and southern Africa (SAfr), as derived from BHRPAR and FAPAR SEVIRI/MSG products using Equation (24).The distribution can be well approximated by the exponential function, FGROUND = Ae −B•LAI for different periods, which indicates a strong consistency between LAI and FGROUND at the resolution of SEVIRI products.

Potential applications of SEVIRI vegetation products
Different studies have assessed the quality of SEVIRI vegetation products [60][61][62].MSG/SEVIRI LAI and FAPAR products were included in a detailed scientific analysis for quality information provision within the Evaluation and Quality Control (EQC) of climate data products derived from satellite and in situ observations to be catalogued within the Copernicus Climate Change Service (C3S) Climate Data Store (CDS) [63].The products have been extensively used as input for environmental applications [64,65].Particularly, the MSG LAI product found applicability to examine the feedback of land surface in response to vegetation changes [66], and in a hydrological model for flood early warning systems [67].Gessner et al. [61] reported that the MSG LAI product is well suited to capture critical periods of growing cycle during the rainy season in West Africa as compared to existing global LAI products derived from polar orbit instruments.The use of MSG FVC and LAI products improved evapotranspiration estimates with respect to ECOCLIMAP-I in a land surface model [68].Daily MSG LAI products were used for automatic derivation of detailed phenological information at the African continent [69].The MSG LAI was fused with VGT/SPOT LAI to fill frequent gaps in VGT product during the rainy season in West Africa [70].Further case studies and product applications of the SEVIRI/MSG vegetation products are provided in the LSA SAF web page (https://landsaf.ipma.pt/).
The SEVIRI/MSG vegetation products are routinely used as input in the operational line of several LSA SAF products.Daily MSG FVC and LAI are required inputs to update information about the state of vegetation for the derivation of the LSA SAF evapotranspiration and (latent and sensible) heat flux products [71].The MSG FVC information on the pixel fraction of vegetation cover is

Potential Applications of SEVIRI Vegetation Products
Different studies have assessed the quality of SEVIRI vegetation products [60][61][62].MSG/SEVIRI LAI and FAPAR products were included in a detailed scientific analysis for quality information provision within the Evaluation and Quality Control (EQC) of climate data products derived from satellite and in situ observations to be catalogued within the Copernicus Climate Change Service (C3S) Climate Data Store (CDS) [63].The products have been extensively used as input for environmental applications [64,65].Particularly, the MSG LAI product found applicability to examine the feedback of land surface in response to vegetation changes [66], and in a hydrological model for flood early warning systems [67].Gessner et al. [61] reported that the MSG LAI product is well suited to capture critical periods of growing cycle during the rainy season in West Africa as compared to existing global LAI products derived from polar orbit instruments.The use of MSG FVC and LAI products improved evapotranspiration estimates with respect to ECOCLIMAP-I in a land surface model [68].Daily MSG LAI products were used for automatic derivation of detailed phenological information at the African continent [69].The MSG LAI was fused with VGT/SPOT LAI to fill frequent gaps in VGT product during the rainy season in West Africa [70].Further case studies and product applications of the SEVIRI/MSG vegetation products are provided in the LSA SAF web page (https://landsaf.ipma.pt/).
The SEVIRI/MSG vegetation products are routinely used as input in the operational line of several LSA SAF products.Daily MSG FVC and LAI are required inputs to update information about the state of vegetation for the derivation of the LSA SAF evapotranspiration and (latent and sensible) heat flux products [71].The MSG FVC information on the pixel fraction of vegetation cover is required to operationally retrieve thermal Land Surface Emissivity [72].The daily MSG FAPAR is also included as input on the Monteith's LUE concept used for the computation of the 10-days Gross Primary Production (MGPP) product recently developed in LSA SAF [73,74].The remaining of this section provides further insights about potential use of the LSA SAF vegetation products, considering three illustrative examples of land processes applications.

Application 1: Monitoring of Seasonal Cycle and Phenology
The enhanced data quality of geoestationary SEVIRI products can be demonstrated by comparing time series of FVC and polar-orbiting satellite sensors such as MODIS and SPOT/VGT across the cloud-prone equatorial regions, such as Amazonian, West Africa and Central Africa.As illustrated in a representative site in the western part of Equatorial Africa (Figure 11), the quality of polar-orbiting satellite products is severely hampered by cloud occurrence.The presence of missing values for long periods and unrealistic short-term variations reduces the ability of polar orbit products to capture the seasonal cycle.
Remote Sens. 2019, 11, x FOR PEER REVIEW 18 of 26 required to operationally retrieve thermal Land Surface Emissivity [72].The daily MSG FAPAR is also included as input on the Monteith's LUE concept used for the computation of the 10-days Gross Primary Production (MGPP) product recently developed in LSA SAF [73][74].The remaining of this section provides further insights about potential use of the LSA SAF vegetation products, considering three illustrative examples of land processes applications.

Application 1: Monitoring of Seasonal Cycle and Phenology
The enhanced data quality of geoestationary SEVIRI products can be demonstrated by comparing time series of FVC and polar-orbiting satellite sensors such as MODIS and SPOT/VGT across the cloud-prone equatorial regions, such as Amazonian, West Africa and Central Africa.As illustrated in a representative site in the western part of Equatorial Africa (Figure 11), the quality of polar-orbiting satellite products is severely hampered by cloud occurrence.The presence of missing values for long periods and unrealistic short-term variations reduces the ability of polar orbit products to capture the seasonal cycle.Conv

Application 2: Interrelation between Vegetation and Rainfall
Accumulated precipitation data derived from 10-days rainfall estimates (RFE 2.0) at a spatial resolution of 0.1° were considered to explore the performance of the SEVIRI FVC product over the Africa continent.The Rainfall Estimator (RFE) is produced by NOAA Climate Prediction Center (NOAA/CPC) merging information from polar-orbiting microwave measurements, geostationary infrared satellite estimates and gridded rainfall gauge measurements [76,77].As illustrated in Figure 12, the SEVIRI FVC is an indicator of spatial and temporal rainfall variability, showing a strong temporal correlation with rainfall measurements (Pearson coefficient of determination r ranging from 0.7 to 0.9) in most of Africa.Conversely, temporal profiles of MSG FVC are consistent from year to year showing realistic smooth variations, even during periods with high cloud coverage during the growing season.These findings are in line with previous studies [61,75].This is mainly because the frequent sampling of its geostationary platform allows to accumulate enough cloud-free observations along a single day.

Application 2: Interrelation between Vegetation and Rainfall
Accumulated precipitation data derived from 10-days rainfall estimates (RFE 2.0) at a spatial resolution of 0.1 • were considered to explore the performance of the SEVIRI FVC product over the Africa continent.The Rainfall Estimator (RFE) is produced by NOAA Climate Prediction Center (NOAA/CPC) merging information from polar-orbiting microwave measurements, geostationary infrared satellite estimates and gridded rainfall gauge measurements [76,77].As illustrated in Figure 12, the SEVIRI FVC is an indicator of spatial and temporal rainfall variability, showing a strong temporal correlation with rainfall measurements (Pearson coefficient of determination r ranging from 0.7 to 0.9) in most of Africa.Figure 13 shows an example of the relationships between FVC and 3-month accumulated rainfall for an herbaceous biome located in southern Africa.A clear interrelation is found between interannual variations of both magnitudes, proving the FVC product to be a proper indicator to detect sensitive areas.In particular, the rainfall deficit occurred during the rainy seasons of 2005, 2013, 2015 and 2016 produced a severe impact on the FVC during the growing season.Although not all regions present similar response to rainfall precipitation, NRT monitoring of FVC anomalies can be used to diagnose in below-normal rainfall values in sensitive areas, and inputted in early warning system of water stress initiation and severe drought events.Figure 13 shows an example of the relationships between FVC and 3-month accumulated rainfall for an herbaceous biome located in southern Africa.A clear interrelation is found between inter-annual variations of both magnitudes, proving the FVC product to be a proper indicator to detect sensitive areas.In particular, the rainfall deficit occurred during the rainy seasons of 2005, 2013, 2015 and 2016 produced a severe impact on the FVC during the growing season.Although not all regions present similar response to rainfall precipitation, NRT monitoring of FVC anomalies can be used to diagnose in below-normal rainfall values in sensitive areas, and inputted in early warning system of water stress initiation and severe drought events.A trend analysis over Africa to report long-term variations due to subtle changes over the period 2004-2017 has been performed.A multiresolution analysis (MRA) wavelet transform (WT) has been applied to the 10-days MSG FVC product [78].Figure 14 shows the derived normalized trend of inter-annual changes over Africa for the period 2004-2017.In general, no trend is observed over central Africa (light blue and green colors) with some distinguishable negative trend areas (dark blue colors) localized at eastern (Horn of Africa, region 1), southern Mozambique (region 2) and western Namibia (region 5).The impact on vegetation is particularly evident in these regions since they have been affected by a precipitation deficit due to El Niño-Southern Oscillation (ENSO) events [79][80][81].Conversely, positive trend values (yellow colors) are mostly located in the grassland biomes of the Sahel zone (region 3) and a few patches in east of Namibia and west of Bostwana (region 4).The Sahel results agree with the widely reported re-greening behavior in this area for the last two decades [82,83] whereas the behavior of region 4 is possibly explained by exceptionally wet years in southwestern Africa since circa 2005.

Summary and Conclusions
This paper describes the algorithm currently used in the LSA SAF operational system to produce in NRT daily and 10-days basis global FVC, LAI and FAPAR biophysical products from SEVIRI/MSG data.A novel algorithm has been proposed to retrieve FVC, which is also the base for retrieving LAI using a semi-empirical method.It relies on a stochastic spectral mixture model which characterizes the large variability of soil background and target vegetation using statistical (Gaussian mixture) models with a finite number of classes, and determines a posteriori probabilities of the background/soil classes making use of the pixel multitemporal trajectory.The FAPAR algorithm uses statistical relationships with BRDF in an optimal angular geometry, general enough for global applications.This study evidences a good consistency among the suite of SEVIRI/MSG (daily and 10days) vegetation products, as well as among FAPAR, LAI and BHRPAR, enabling a realistic partitioning of incoming solar radiation between the canopy and the ground below the canopy.

Summary and Conclusions
This paper describes the algorithm currently used in the LSA SAF operational system to produce in NRT daily and 10-days basis global FVC, LAI and FAPAR biophysical products from SEVIRI/MSG data.A novel algorithm has been proposed to retrieve FVC, which is also the base for retrieving LAI using a semi-empirical method.It relies on a stochastic spectral mixture model which characterizes the large variability of soil background and target vegetation using statistical (Gaussian mixture) models with a finite number of classes, and determines a posteriori probabilities of the background/soil classes making use of the pixel multitemporal trajectory.The FAPAR algorithm uses statistical relationships with BRDF in an optimal angular geometry, general enough for global applications.This study evidences a good consistency among the suite of SEVIRI/MSG (daily and 10-days) vegetation products, as well as among FAPAR, LAI and BHRPAR, enabling a realistic partitioning of incoming solar radiation between the canopy and the ground below the canopy.
Main advantages and relevant information derived from this study are summarized below: 1.
NRT daily (MDFVC, MDLAI, MDFAPAR) and 10-days (MTFVC, MTLAI, MTFAPAR) products are generated and disseminated from LSA SAF since January 2004 over the geostationary Meteosat disk offering almost fifteen years of an alternative dataset to the user community.2.
The 10-days (MTFVC-R, MTLAI-R and MTFAPAR-R) CDRs are provided as a suite of EUMETSAT climate products data records estimated consistently along the years using the latest versions of the whole processing chain algorithms.The 10-days products could be suitable for a community of users that requires observations representative of a 30-day period with at frequency of 10 days (e.g., numerical weather and climate models, and flood forecasting systems).

3.
The daily SEVIRI/MSG timeliness of the distribution of the observations and its smaller compositing period avoids possible shifts regarding the actual state of the vegetation (e.g., for an early estimate of key phenological parameters and seasonal production).4.
The absence of gaps and the high temporal frequency and continuity of the products over Africa offer major potentials for NRT monitoring of land cover dynamics for applications that require frequent observations such as agriculture, and food management.

5.
The SEVIRI/MSG vegetation products have demonstrated its suitability to accurately resolving long term changes in large regions, allowing improving the understanding of interactions between land surface and climate.
The MSG programme is expected to provide observations for the next years to ensure long-term operations.In order to improve the coverage of the LSA SAF service, an ongoing work is the clone of MSG algorithms to produce similar LSA SAF products based on the Indian Ocean Data Coverage (IODC) observations.A priority work in LSA SAF is the adaptation of the algorithms to take advantage of enhanced capabilities in terms of spatial and spectral resolution of the future EUMETSAT sensors.The algorithm will be adapted to the enhanced characteristics of the Flexible Combined Imager (FCI)/Meteosat Third Generation satellites (MTG), allowing a more accurate characterization of the soil and vegetation signatures and improved screening of snow pixels.The FCI based products will ensure the continuity of the service during the lifetime of the MTG program.

Supplementary Materials:
The following are available online at http://www.mdpi.com/2072-4292/11/18/2103/s1, Figure S1: Theoretical FAPAR error as a function of input k 0 and k 2 errors for a given k 1 error of 0.01.Two different cases have been considered: Low FAPAR values (a) and high FAPAR values (b)., Table S1: Cover-dependent clumping index values for LAI algorithm based on the GLC2000 land cover classification based on values obtained in [51], Table S2: VEGA products QF information.The default missing value for the product fields is −10.The associated error estimate fields for unprocessed pixels take different negative values, depending on the identified problem (default missing value = −10).Main identified problems in the VEGA products and empirical thresholds used to blind problematic areas.Note that although the missing value for the product fields is unique (−10), associated error estimate fields for unprocessed pixels take different negative values, depending on the identified problem (default missing value = −10).

Figure 1 .
Figure 1.Flow chart of the algorithm for FVC, LAI and FAPAR determination.

Figure 1 .
Figure 1.Flow chart of the algorithm for FVC, LAI and FAPAR determination.

Figure 2 .
Figure 2. Flow chart of the stochastic SMM algorithm for FVC determination.

Figure 2 .
Figure 2. Flow chart of the stochastic SMM algorithm for FVC determination.
considering a simplified bidimensional problem with 3 vegetation clusters (V1, V2, V3) and 3 soil clusters (S1, S2, S3).The size and orientation in the elliptical probability density contours are determined by the covariance of clusters  .It is obvious that (S3, V1) is the most likely model in the mixture.However, other different model combinations, such as (S2, V1) and (S3, V2), have non negligible probabilities and may yield to the same mixture spectrum.

Figure 3 .
Figure 3. Illustration of the model selection in a bidimensional feature space.Dashed lines join random spectra drawn from soil class S3 and vegetation classes V1, V2 and V3.

Figure 3 .
Figure 3. Illustration of the model selection in a bidimensional feature space.Dashed lines join random spectra drawn from soil class S 3 and vegetation classes V 1 , V 2 and V 3 .

Figure 4 .
Figure 4. Illustration of the probabilistic mixing model concept over in the k 0 space of SEVIRI channels:(a) channel 1 and 2; (b) channel 2 and 3. Elliptical probability density contours associated with clusters (from V 1 to V 3 for vegetation and from S 1 to S 6 for soil).BDF and BEF correspond to broadleaved deciduous and broadleaved evergreen forest, respectively.BIC s and BIC v refer to the values of BIC for the soil and vegetation mixture models, respectively.

Figure 5 .
Figure 5. Projection of the SEVIRI/FVC product corresponding to the 15th of June 2015 onto the k0 feature space of (a) channels 1 and 2; (b) channels 2 and 3. Red circles correspond to soil and vegetation pure pixels for the Southern Africa (SAfr) region.

Figure 5 .
Figure 5. Projection of the SEVIRI/FVC product corresponding to the 15th of June 2015 onto the k 0 feature space of (a) channels 1 and 2; (b) channels 2 and 3. Red circles correspond to soil and vegetation pure pixels for the Southern Africa (SAfr) region.

Figure 6 .
Figure 6.MSG Daily LAI (top), FVC (middle) and FAPAR (bottom) LSA SAF product composition corresponding to the 17th of April 2014 products (left panels) and their respective error estimates (right panels).Location of the four LSA SAF geographical areas is also provided.Rectangles in the LAI field are associated to Euro, NAfr, SAfr and Same SEVIRI geographical regions for illustrating purposes.

Figure 6 .
Figure 6.MSG Daily LAI (top), FVC (middle) and FAPAR (bottom) LSA SAF product composition corresponding to the 17th of April 2014 products (left panels) and their respective error estimates (right panels).Location of the four LSA SAF geographical areas is also provided.Rectangles in the LAI field are associated to Euro, NAfr, SAfr and Same SEVIRI geographical regions for illustrating purposes.

Figure 8 .
Figure 8. Monthly fraction of valid pixels for daily MSG FVC product during the entire 2014 year over the four SEVIRI geographical areas.

Figure 8 .
Figure 8. Monthly fraction of valid pixels for daily MSG FVC product during the entire 2014 year over the four SEVIRI geographical areas.

Figure 9 .
Figure 9.Time series of LAI, FAPAR and FVC daily and 10-days products at a Deciduous Broadleaf Forest site.

Figure 9 .
Figure 9.Time series of LAI, FAPAR and FVC daily and 10-days products at a Deciduous Broadleaf Forest site.

Figure 10 .
Figure 10.Joint probability density plots between FGROUND and MSG daily LAI at several different periods of the year.Top Figures correspond to the Euro SEVIRI zone, whereas bottom Figures correspond to the SAfr SEVIRI region.

Figure 10 .
Figure 10.Joint probability density plots between FGROUND and MSG daily LAI at several different periods of the year.Top Figures correspond to the Euro SEVIRI zone, whereas bottom Figures correspond to the SAfr SEVIRI region.

Figure 11 .
Figure 11.Comparison between LSA SAF products and equivalent MODIS and SPOT/VGT over a three-year period representative example at a site in Gabon.

Figure 11 .
Figure 11.Comparison between LSA SAF products and equivalent MODIS and SPOT/VGT over a three-year period representative example at a site in Gabon.

Figure 12 .
Figure 12.Pearson correlation coefficient between time series of FVC and 3-month accumulated precipitation for the 2004-2017 period.

Figure 12 .
Figure 12.Pearson correlation coefficient between time series of FVC and 3-month accumulated precipitation for the 2004-2017 period.

4. 3 .
Application 3: The Detection of Inter-Annual Vegetation Trends over the Period 2004-2017 A trend analysis over Africa to report long-term variations due to subtle changes over the period 2004-2017 has been performed.A multiresolution analysis (MRA) wavelet transform (WT) has been

26 Figure 14 .
Figure 14.Trend of inter-annual changes over Africa derived by applying the MRA-WT method to 10-days FVC time series (2004-2017).

Figure 14 .
Figure 14.Trend of inter-annual changes over Africa derived by applying the MRA-WT method to 10-days FVC time series (2004-2017).

Table 1 .
Product Requirements for MSG vegetation products, in terms resolution and accuracy.