This paper describes the selected algorithm for the ESA climate change initiative vegetation parameters project. Multi- and hyper-spectral, multi-angular, or multi-sensor top-of-canopy reflectance data call for an efficient generic retrieval system which can improve the consistent retrieval of standard canopy parameters as albedo, Leaf Area Index (LAI), Fraction of Absorbed Photosynthetically Active Radiation (fAPAR) and their uncertainties, and exploit the information to retrieve additional parameters (e.g. leaf pigments). We present a retrieval system for canopy and sub-canopy parameters (OptiSAIL), which is based on a model comprising SAIL (canopy reflectance), PROSPECT-D (leaf properties), TARTES (snow properties), a soil model (soil reflectance anisotropy, moisture effect), and a cloud contamination model. The inversion is gradient based and uses codes created by Automatic Differentiation. The full per pixel covariance-matrix of the retrieved parameters is computed. For this demonstration, single observation data from the Sentinel-3 SY_2_SYN (synergy) product is used. The results are compared with the MODIS 4-day LAI/fAPAR product and PhenoCam site photography. OptiSAIL produces generally consistent and credible results, at least matching the quality of the technically quite different MODIS product. The system is computationally efficient with a rate of 150 pixel s(-1) (7 ms per pixel) for a single thread on a current desktop CPU using observations on 26 bands. Not all of the model parameters are well determined in all situations. Significant correlations between the parameters are found, which can change sign and magnitude over time. OptiSAIL appears to meet the design goals and puts real-time processing with this kind of system into reach.

The Granger Causality (GC) statistical test explores the causal relationships between different time series variables. By employing the GC method, the underlying causal links between environmental drivers and global vegetation properties can be untangled, which opens possibilities to forecast the increasing strain on ecosystems by droughts, global warming, and climate change. This study aimed to quantify the spatial distribution of four distinct satellite vegetation products' (VPs) sensitivities to four environmental land variables (ELVs) at the global scale given the GC method. The GC analysis assessed the spatially explicit response of the VPs: (i) the fraction of absorbed photosynthetically active radiation (FAPAR), (ii) the leaf area index (LAI), (iii) solar-induced fluorescence (SIF), and, finally, (iv) the normalized difference vegetation index (NDVI) to the ELVs. These ELVs can be categorized as water availability assessing root zone soil moisture (SM) and accumulated precipitation (P), as well as, energy availability considering the effect of air temperature (T) and solar shortwave (R) radiation. The results indicate SM and P are key drivers, particularly causing changes in the LAI. SM alone accounts for 43%, while P accounts for 41%, of the explicitly caused areas over arid biomes. SM further significantly influences the LAI at northern latitudes, covering 44% of cold and 50% of polar biome areas. These areas exhibit a predominant response to R, which is a possible trigger for snowmelt, showing more than 40% caused by both cold and polar biomes for all VPs. Finally, T's causality is evenly distributed amongst all biomes with fractional covers between similar to 10 and 20%. By using the GC method, the analysis presents a novel way to monitor the planet's ecosystem, based on solely two years as input data, with four VPs acquired by the synergy of Sentinel-3 (S3) and 5P (S5P) satellite data streams. The findings indicated unique, biome-specific responses of vegetation to distinct environmental drivers.