Journal topic
Earth Syst. Sci. Data, 12, 1101–1116, 2020
https://doi.org/10.5194/essd-12-1101-2020
Earth Syst. Sci. Data, 12, 1101–1116, 2020
https://doi.org/10.5194/essd-12-1101-2020

Data description paper 13 May 2020

Data description paper | 13 May 2020

# A spatially downscaled sun-induced fluorescence global product for enhanced monitoring of vegetation productivity

A spatially downscaled sun-induced fluorescence global product for enhanced monitoring of vegetation productivity
Gregory Duveiller1, Federico Filipponi1, Sophia Walther2, Philipp Köhler3, Christian Frankenberg3,4, Luis Guanter5, and Alessandro Cescatti1 Gregory Duveiller et al.
• 1European Commission Joint Research Centre, Ispra, Italy
• 2Max Planck Institute for Biogeochemistry, Jena, Germany
• 3Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, CA, USA
• 4Jet Propulsion Laboratory, California Institute of Technology, CA, USA
• 5Centro de Tecnologías Físicas, Universitat Politècnica de València, Valencia, Spain

Correspondence: Gregory Duveiller (gregory.duveiller@ec.europa.eu)

Abstract

Sun-induced chlorophyll fluorescence (SIF) retrieved from satellite spectrometers can be a highly valuable proxy for photosynthesis. The SIF signal is very small and notoriously difficult to measure, requiring sub-nanometre spectral-resolution measurements, which to date are only available from atmospheric spectrometers sampling at low spatial resolution. For example, the widely used SIF dataset derived from the GOME-2 mission is typically provided in 0.5 composites. This paper presents a new SIF dataset based on GOME-2 satellite observations with an enhanced spatial resolution of 0.05 and an 8 d time step covering the period 2007–2018. It leverages on a proven methodology that relies on using a light-use efficiency (LUE) modelling approach to establish a semi-empirical relationship between SIF and various explanatory variables derived from remote sensing at higher spatial resolution. An optimal set of explanatory variables is selected based on an independent validation with OCO-2 SIF observations, which are only sparsely available but have a high accuracy and spatial resolution. After bias correction, the resulting downscaled SIF data show high spatio-temporal agreement with the first SIF retrievals from the new TROPOMI mission, opening the path towards establishing a surrogate archive for this promising new dataset. We foresee this new SIF dataset becoming a valuable asset for Earth system science in general and for monitoring vegetation productivity in particular. The dataset is available at https://doi.org/10.2905/21935FFC-B797-4BEE-94DA-8FEC85B3F9E1 .

1 Introduction

Mapping and monitoring the spatial and temporal patterns of terrestrial gross primary productivity (GPP) through the use of satellite remote sensing are of paramount interest for vegetation, ecosystem and climate science. While the rate of terrestrial photosynthesis cannot be directly measured from space, recent research has demonstrated that sun-induced chlorophyll fluorescence (SIF) retrieved from satellite spectrometers can be a highly valuable proxy . SIF is generally positively correlated with leaf photochemistry during specific light conditions that are common across the globe and should thus serve as valid proxy for GPP, even though the mechanistic link between the two is complex . The origin of this signal is the fluorescence of chlorophyll a, consisting of a re-emission of absorbed photons at lower-energy wavelengths (from 650 to 850 nm, with peaks at approximately 690 and 740 nm). This physical mechanism allows photosynthetic organisms to adjust instantaneously to rapid perturbations in environmental conditions of light, temperature and water availability, before the heat dissipation mechanism of non-photochemical quenching can be triggered .

Chlorophyll a fluorescence has been extensively studied in laboratories, from the sub-cellular scale up to the leaf scale, for many decades (Baker2008), but only recently has it been possible to retrieve it from space-borne remote-sensing platforms at the global scale. The majority of leaf-scale fluorescence research focused on induction kinetics, e.g. the temporal course of fluorescence (induced by saturating light pulses under various conditions) during the period of light acclimation of dark-adapted leaves. On the contrary, passive remote sensing from satellite platforms allows only measuring the light-adapted steady-state fluorescence (or SIF). However, this has opened the possibility of exploring the carbon cycle dynamics in various terrestrial ecosystems, ranging from tropical forests to the northern tundra , passing by agricultural landscapes . It has been used to analyse the response of plants to water availability and drought but also as an alternative measure to typical remote-sensing greenness indices such as the NDVI (normalised difference vegetation index) and EVI (enhanced vegetation index) for vegetation phenology .

Despite the abundance of recent studies, currently there is no operational satellite instrument specifically dedicated to the measurement of SIF. The SIF signal is notoriously difficult to measure from space, as it represents only 1 % to 5 % of the total reflected radiation in the near-infrared spectrum that is detected by a remote-sensing instrument . Measuring it requires both a high spectral resolution and high signal-to-noise ratio, which generally comes at the expense of spatial detail. The first mission specifically designed for this purpose is the European Spatial Agency's Earth Explorer – FLEX (FLuorescence EXplorer), due to be launched in 2023. All datasets that are currently available come from missions originally conceived for measuring atmospheric trace gas concentrations. These can be separated into those providing highly precise soundings but with sparse and thus spatially discontinuous samples (e.g. GOSAT, TanSat and OCO-2) and those which do provide a spatially continuous coverage but with less accuracy due to a reduced spectral resolution (e.g. SCIAMACHY, GOME-2 and, more recently, TROPOMI). After the initial discovery that such instruments could be serendipitously used to retrieve SIF from space , different retrieval methods have been developed to adjust to the various sensors , each with specific properties such as different spectral fitting windows and cloud filtering procedures, resulting in various distinct SIF datasets.

For a SIF dataset to become truly useful to the Earth system science community as a proxy for GPP, the following properties should be optimised: (1) the temporal archive should be as long as possible, (2) the revisit frequency should be high (ideally at daily or even sub-daily scale), (3) the geographic extent should be global (ideally gap-free), and (4) the spatial resolution should be high enough to relate to distinct land cover elements or plant functional types (PFTs). No current individual SIF dataset adequately satisfies these specifications. Retrievals based on GOME-2 are perhaps at present those closest to the mark, thanks to their spatially continuous sampling design combined with their temporal archive, which has been running since 2007. However, their spatial resolution is 0.5 (approximately 50 km), which is too low for many applications in Earth science. The situation changed in 2018 with the operational arrival of TROPOMI, from which SIF can be retrieved from a ground pixel footprint of around 7 km×3.5 km at nadir , but the shallow temporal depth of this data record will preclude its use for many applications for various years to come if no compatible archive is established.

To remedy the lack of spatial detail, several studies have proposed enhancing the spatial resolution of currently available SIF data. The first of such studies leverages on the concept of light-use efficiency (LUE) used in GPP modelling to constrain the spatial re-allocation of GOME-2 SIF values within a ${\mathrm{0.5}}^{\circ }×{\mathrm{0.5}}^{\circ }$ grid cell. This model assumes that SIF can be estimated as a function of greenness, as described by the NDVI, which is then down-regulated based on water availability and temperature, characterised respectively by evapotranspiration (ET) and land surface temperature (LST). NDVI, ET and LST are all satellite remote-sensing variables that are available at the fine scale (e.g. $\le {\mathrm{0.05}}^{\circ }$), can be aggregated to the SIF resolution to establish the relationship over a local spatio-temporal window and can then be used to predict SIF at higher resolution. Since then, other studies have adopted entirely data-driven approaches relying on machine-learning techniques to either reconstruct SIF based on high-spatial-resolution reflectances from another satellite or to gap-fill spatially sparse OCO-2 data . Compared to these entirely empirical approaches, the method has the particularity that it remains data-driven, as the original SIF signal is preserved at each step, yet the downscaling remains physiologically constrained by the assumptions of a semi-empirical process-based model grounded on theory.

The objective of this work is to present an improved and updated downscaled SIF dataset based on the methodology. Besides an extension of the archive until the end of 2018, the new dataset has a higher temporal frequency and is constructed from both updated explanatory variables and different SIF retrievals. The optimal model configuration is selected based on a comparison with validation data composed of high-spatial-resolution and high-spectral-resolution OCO-2 SIF observations. Finally, the resulting downscaled SIF data are compared to the new TROPOMI retrievals in view of constituting an archive for this promising SIF data stream.

2 Material and methods

## 2.2 SIF data from GOME-2, OCO-2 and TROPOMI

The first source of SIF data is from the GOME-2 instrument aboard the MetOp-A satellite. The SIF retrievals are acquired at around 09:30 local time (at the Equator) and at a spectral wavelength around 740 nm. The spatial footprint of these measurements was approximately 40 km×80 km before 2013 and 40 km×40 km after that date as a consequence of the additional monitoring capacities due to the arrival of MetOp-B. Two retrievals from GOME-2 are used in this paper. The first was developed by and is referred to as JJ in this study. Version 25 of this retrieval was used in the original paper, which was superseded by version 27, used here. The second retrieval was proposed by and is referred to here as PK. For both JJ and PK datasets, the individual retrievals are filtered to keep only those with solar zenith angles below 70, local solar time between 08:00 and 14:00, and effective cloud cover fraction below or equal to 0.5. Then they are gridded into 0.5 cells by taking the mean value over a period of 16 d, in line with the MODIS reflectance products. Since the GOME-2 acquisition time is early in the morning, values need to be multiplied by a daily correction factor to make them comparable with estimates from other sources. proposed a simple approach to convert the instantaneous SIF to a daily average, which accounts for variations in overpass time, length of day and solar zenith angle. For the JJ dataset, such a daily correction factor is already provided in the dataset. For the PK dataset it is applied using the implementation used in as follows:

$\begin{array}{}\text{(1)}& \begin{array}{rl}\mathrm{SIF}=& \mathrm{SIF}\left({t}_{\mathrm{m}}\right)\cdot \frac{\mathrm{1}}{\mathrm{cos}\left(\mathit{\theta }\left({t}_{\mathrm{m}}\right)\right)}\\ & \cdot \underset{t={t}_{\mathrm{m}}-\mathrm{12}\phantom{\rule{0.125em}{0ex}}\mathrm{h}}{\overset{t={t}_{\mathrm{m}}+\mathrm{12}\phantom{\rule{0.125em}{0ex}}\mathrm{h}}{\int }}\mathrm{cos}\left(\mathit{\theta }\left(t\right)\right)\cdot H\left(\mathrm{cos}\left(\mathit{\theta }\left(t\right)\right)\right)\mathrm{d}t,\end{array}\end{array}$

where θ(tm) is the solar zenith angle at the time of measurement tm and the integral is computed numerically in 10 min time steps (dt), with the heavy-side step function H zeroing out negative values of cos (θ). This daily correction factor is applied to the SIF retrievals prior to compositing and gridding.

The second source of SIF data is the OCO-2 platform . These retrievals are made from soundings over footprints of 1.3 km×2.25 km at nadir, which together create a 10 km wide stripe with a revisit time of 16 d. Retrievals are made at two wavelengths, 757 and 771 nm, at 13:00 local time (at the Equator). The SIF values at both wavelengths are first corrected daily, and then, to render them comparable to the GOME-2 SIF at 740 nm, they are combined together using the following formula:

$\begin{array}{}\text{(2)}& {\mathrm{SIF}}_{\mathrm{740}\phantom{\rule{0.125em}{0ex}}\mathrm{nm}}=\mathrm{1.56}\cdot \left({\mathrm{SIF}}_{\mathrm{757}\phantom{\rule{0.125em}{0ex}}\mathrm{nm}}+\mathrm{1.8}\cdot {\mathrm{SIF}}_{\mathrm{771}\phantom{\rule{0.125em}{0ex}}\mathrm{nm}}\right)/\mathrm{2}.\end{array}$

The scaling factors were determined based on a reference SIF emission shape derived from leaf-level measurements conducted by . Although ideally this relationship should be established at the canopy level, this would prove to be a daunting task given the large canopy-scale variability encountered across the globe. Furthermore, the canopy effects are expected to be more prominent at smaller wavelengths than 740 nm, so we deemed that the leaf-level relationship should be sufficient for our purpose. All individual observations between 2015 and 2017 are gridded into a common 0.05 grid to match the MODIS grid for explanatory variables (and which becomes the grid of the final downscaled SIF product). This results in a validation dataset that contains more than 138 million records distributed across the globe. Each record is further attached to ancillary information of the major climate zone group in which it falls (tropical, dry, temperate, continental or polar) based on the revisited Köppen–Geiger classification and the dominant vegetation type derived from the European Space Agency's Climate Change Initiative land cover maps (ESA2017). Due to their superior spectral and spatial resolution, the SIF retrievals from OCO-2 are considered here to be a reference.

The third source of SIF data is from TROPOMI, the single instrument aboard the Sentinel-5 Precursor satellite . These data have been available daily since early 2018, with a footprint of 7 km×3.5 km at nadir. The daily corrected retrieval results with the default filtering as described in were aggregated to the 0.05 grid using the 16 d compositing scheme used for GOME-2. Similar to the other SIF datasets, the filtering consists of thresholds for the fit quality, clouds, and extremely low and high radiance levels and high solar zenith angles.

Table 1Boundary conditions used to initialise the calibration procedure for every local optimisation dependent on the explanatory variables used as a proxy for the vegetation V and for the water W and thermal T stress.

Figure 1Exploring the agreement, absolute bias and correlation between the OCO-2 validation dataset and downscaled GOME-2 SIF based on different retrievals (PK for , or JJ for ) and using different explanatory variables (NDVI, EVI or NIRv for V). Not all combinations were calculated. Yellow circles indicate the use of ET as a proxy for hydric stress (W), while triangles represent the use of NDWI instead. Empty triangles correspond to a subset in which the proxy for thermal stress (T) is changed from MYD to MOD (this was only done for the PK dataset). Each metric is calculated based on all available samples within the period 2015–2017.

## 2.3 Parameterisation of the downscaling methodology

The downscaling procedure follows three steps: (1) an aggregation of the explanatory variables to the low spatial resolution, (2) a calibration of the downscaling function over a local spatio-temporal window of low-spatial-resolution data and (3) the application of the calibrated function to the original explanatory variables at high spatial resolution to result in a downscaled SIF estimation. The downscaling function always takes the following form:

$\begin{array}{}\text{(3)}& \mathrm{SIF}=f\left(V\right)\cdot f\left(W\right)\cdot f\left(T\right),\end{array}$

where f(V) is a function of a vegetation index V, which is down-regulated by two other functions, f(W) and f(T), both outputting a value between 0 and 1 based on indicators W and T, respectively representing water and thermal stress. A quadratic, a sigmoid and a Gaussian function are respectively used to model f(V), f(W) and f(T), resulting in the following expanded expression:

$\begin{array}{}\text{(4)}& \begin{array}{rl}\mathrm{SIF}=& {b}_{\mathrm{2}}{V}^{{b}_{\mathrm{1}}}\cdot \frac{\mathrm{1}}{\mathrm{1}+\mathrm{exp}\left({b}_{\mathrm{3}}\left({b}_{\mathrm{4}}-W\right)\right)}\\ & \cdot \mathrm{exp}\left(-\mathrm{0.5}{\left[\frac{T+{b}_{\mathrm{5}}}{{b}_{\mathrm{6}}}\right]}^{\mathrm{2}}\right).\end{array}\end{array}$

The values of these bi parameters can have some physiological interpretation. b1 determines the degree of linearity of the general relationship between the vegetation index V and SIF, which is allowed to range from having some saturating effect (b1=0.5) to a mild exponential effect (b1=1.5). b2 is a simple multiplicative scaling factor, as the units of V are not the same as those of SIF. The function to down-regulate f(V) based on water stress is characterised by b3, which determines the steepness of the sigmoid function f(W), and b4, which represents the value of W at which f(W)=0.5, thus defining the point on the W scale at which the hydric stress reduces f(V) by half. The thermic down-regulating function f(T) depends on its position, determined by the parameter b5, which can be interpreted as an optimal temperature (in degrees kelvin) for photosynthesis to occur over the region. Finally, b6 defines the width of the Gaussian function f(T), meaning higher values represent lower sensitivity to temperature variations. The bi parameters are estimated in the calibration phase for each time step using an adaptable spatial moving window always composed of 40 observations. These are selected within the limits of a box with an 11×11 GOME-2 pixel window, and they include the central pixel and the 39 available observations that are nearest to it. The calibration of the bi parameters is done over this selection using the quasi-Newton L-BFGS-B optimisation algorithm  implemented in the core R package stats, which allows the setting of a lower and an upper boundary for each parameter.

Figure 2Inter-comparison of the performance of different downscaled products disaggregated per major vegetation type and climate zone. The downscaled products are based on either the PK or JJ retrievals combined with either the old set of explanatory variables (NDVI, ET and MYD) or the new one (NIRv, NDWI and MYD). Each graph plots the agreement of a given downscaled product with OCO-2 over the period 2015–2017 against the agreement of another product with OCO-2 for the same period. The overall bias (B), correlation (r) and agreement index (L) are reported for at the corresponding corner of each graph, with numbers in bold highlighting the better values in each pairwise comparison.

The difference with the original procedure described in lies in using two different SIF retrievals, SIFJJ and SIFPK, and different explanatory variables for V, W and T. NDVI, EVI and NIRv are three alternative spectral indices explored for V; ET and NDWI are used for W; and morning and afternoon LST, labelled MOD and MYD, are used for T. Various combinations of variables are tested, but giving priority to those that are suspected to provide more relevant information, as the computational cost of each run was not negligible. For all tests we used the initial conditions exposed in Table 1 for a period coinciding with that of the OCO-2 validation dataset, i.e. 2015–2017. The individual values of each downscaled dataset are matched in space and time with those of the OCO-2 dataset and grouped per climate zone and dominant vegetation type.

## 2.4 Quantifying agreement

To quantify the agreement between different sources of SIF data, the λ index of agreement proposed by is used in addition to regular metrics of correlation and bias. This metric quantifies the degree of agreement between two sets of values, x and y, considering both the bias between them and their level of correlation, all within a single number ranging from 0 to 1. It has an added advantage of being symmetric, unlike a typical measurement of goodness of fit such as the coefficient of determination, and yet can still be interpreted as a familiar correlation coefficient when there is no bias. The calculation of this index is as follows:

$\begin{array}{}\text{(5)}& \mathit{\lambda }=\mathrm{1}-\frac{{n}^{-\mathrm{1}}{\sum }_{i=\mathrm{1}}^{n}\left({x}_{i}-{y}_{i}{\right)}^{\mathrm{2}}}{{\mathit{\sigma }}_{x}^{\mathrm{2}}+{\mathit{\sigma }}_{y}^{\mathrm{2}}+\left({\mathit{\mu }}_{x}-{\mathit{\mu }}_{y}{\right)}^{\mathrm{2}}+\mathit{\kappa }},\end{array}$

where μ and σ represent the mean and standard deviation, respectively. The numerator of the fraction in Eq. (5) is the mean squared deviation between x and y, while the denominator quantifies the maximum deviation this set of points could take. The term κ represents the covariance between x and y. Including it in the denominator ensures that λ does not take negative values when x and y are anti-correlated, but it also unnecessarily inflates λ non-linearly. Therefore, κ is set to zero when the correlation between x and y is positive and otherwise takes the value of $\mathit{\kappa }=\mathrm{2}{n}^{-\mathrm{1}}|{\sum }_{i=\mathrm{1}}^{n}\left({x}_{i}-{\mathit{\mu }}_{x}\right)\left({y}_{i}-{\mathit{\mu }}_{y}\right)|$.

We also use a variant of λ that quantifies only the unsystematic contribution to the agreement, i.e. after the systematic bias between x and y is removed. To calculate this λu, the numerator in Eq. (5) needs to relate only to the unsystematic component of the deviations instead of the total deviations. The mean squared unsystematic deviation is calculated based on the orthogonal distances h from the principal plane between x and y, resulting in the following expression for λu:

$\begin{array}{}\text{(6)}& {\mathit{\lambda }}_{\mathrm{u}}=\mathrm{1}-\frac{{n}^{-\mathrm{1}}{\sum }_{i=\mathrm{1}}^{n}{h}_{i}^{\mathrm{2}}}{{\mathit{\sigma }}_{x}^{\mathrm{2}}+{\mathit{\sigma }}_{y}^{\mathrm{2}}+\left({\mathit{\mu }}_{x}-{\mathit{\mu }}_{y}{\right)}^{\mathrm{2}}+\mathit{\kappa }}.\end{array}$

To characterise the principal plane between x and y and obtain the distances h, we use an eigendecomposition of the covariance matrix containing vectors x and y, as described in detail in Sect. 5 of the Supplement in . This also provides the slope and intercept of the lines between x and y in a symmetric way; i.e. the values remain unchanged when x and y are interchanged. This slope and intercept are used in our study to quantify and correct the bias between downscaled GOME-2 SIF values and TROPOMI retrievals.

3 Results

## 3.1 Benchmarking the downscaled datasets

The OCO-2 validation dataset is used to identify the optimal combination of explanatory variables and input dataset to produce the enhanced GOME-2 downscaled dataset. Since generating these datasets is computationally expensive, not all combinations were calculated. We instead explore how replacing a single variable at a time affects the results, starting from the initial configuration from the paper: the JJ SIF retrieval downscaled with NDVI, ET and MYD. The results are summarised in Fig. 1.

The first variable to substitute is ET with NDWI. NDWI is a product that is at a lower level than ET, requiring far fewer assumptions and thus rendering the downscaling independent from other sources of information such as the climate re-analysis data and the eddy-covariance towers used in calibrating ET. Substituting ET for NDWI results in a marked reduction in agreement for the JJ dataset due to a loss in correlation and an increase in bias. For the PK dataset, the drop in correlation is smaller and is accompanied with a reduction in the bias, resulting in an increase in agreement.

Figure 3General benchmarking of the performance of GOME-2 SIF downscaling versus OCO-2 validation measures for the period 2015–2017. The old downscaling method refers to that based on the original explanatory variables used in , i.e. NDVI, ET and afternoon LST, while the new method refers to the best explanatory variables identified in this study, i.e. NIRv, NDWI and afternoon LST. PK refers to the GOME-2 SIF retrieval proposed by , while JJ refers to that by . The agreement refers to the index of agreement λ proposed by .

The second variable change that is explored is the overpass time of the LST data serving as a proxy for the thermal stress. Replacing the afternoon overpass (MYD) used initially with the morning one (MOD) could be relevant, as it is closer to the overpass time of the GOME-2 instrument. On the other hand, data from the morning MODIS instrument aboard the Terra platform may be of lower quality due to sensor degradation . As we did not expect this change to make a strong difference, and to prioritise computational resources, this analysis is done only on PK data. It shows an increase in both the correlation (marginally) and the bias (more substantially), and as a result this reduces the overall agreement. We do not expect the results to be much different for the JJ dataset.

The third variable to exchange is the vegetation index from NDVI to EVI and NIRv. For both JJ and PK retrievals there is a progressive increase in correlation and agreement, with the best results stemming from the use of NIRv. Replacing NDVI with EVI or NIRv even reduces slightly the bias for JJ, whilst it marginally increases for PK. For JJ, the loss of agreement due to the replacement of ET by NDWI appears to be considerably mitigated by the use of NIRv instead of NDVI.

In view of these results, a single common configuration of explanatory variables consisting of NIRv, NDWI and MYD is selected for downscaling both JJ and PK retrievals over the longer time period from 2007 to 2018. Both of these downscaled products are made available together as outputs of this study and will be referred to as a single downscaled SIF dataset containing two separate products .

Figure 4Selection of spatio-temporal subsets of the newly downscaled PK SIF product. Each box covers a region of ${\mathrm{8}}^{\circ }×{\mathrm{8}}^{\circ }$, covering parts of (A) the European Alps and their surroundings, (B) Bolivia, (C) the African Great Lakes, (D) the Indus Valley and the Indo-Gangetic Plain, and (E) the US Corn Belt (the locations of these areas are shown in Fig. 5). The corresponding time is mentioned above each image. Although the temporal frequency of the dataset is 8 d, only an image every 24 d is presented here to accentuate the seasonal dynamics. Notice also that the five time series are not synchronised. Grey areas indicate lack of data.

To delve further into the details, Fig. 2 compares how both JJ and PK retrievals fare with respect to OCO-2 when downscaled either with the old configuration of explanatory variables or the new one, but disaggregating the results for different vegetation types within distinct climate zones. Each graph plots the index of agreement of a given downscaled product with OCO-2 against the index of agreement of the other downscaled product with OCO-2 so that for each point, if it falls on a given side of the 1:1 line, the product of that side has higher agreement with the references. Figure 2a shows that, with the original explanatory variables of NDVI, ET and MYD, the JJ retrieval agrees more with the OCO-2 dataset than the PK retrieval, as the latter is generally affected by a bias. However, some vegetation types in the tropics, such as croplands, grasslands and evergreen broadleaf forests, do have higher agreement with OCO-2 in the PK product, and this is accentuated with the new downscaling methodology (Fig. 2c and d). The new methodology only marginally improves the JJ product with respect to the original downscaled SIF data (Fig. 2b). Figure 2d provides an inter-comparison of the two downscaled products made available here, potentially guiding users to prefer one or the other depending on the type of climate or the dominant vegetation cover of their area of interest.

To resume the outcome of this benchmarking with respect to the OCO-2 validation dataset, Fig. 3 illustrates how both SIF retrievals fare with the new set of downscaling explanatory variables compared both to the old downscaling method and to SIF at the original spatial resolution. For both retrievals, the improvement in agreement between GOME-2 and OCO-2 is due more to the actual spatial downscaling than to the choice of the downscaling variables. The original PK retrieval shows higher agreement than the original JJ retrieval due to higher correlation despite a larger bias, but the downscaling procedure improves both to very similar levels of agreement and correlation.

Figure 5Maximum value of SIF over the entire new PK downscaled product. Values are in megawatts per square metre per steradian per nanometre (mW m−2 sr−1 nm−1), the spatial resolution is 0.05 and the temporal range spans from early 2007 to the end of 2018. The white boxes represent the zones illustrated in Fig. 4.

Figure 6Agreement between the downscaled PK SIF product and TROPOMI SIF retrievals for the period from 16 April to 23 December 2018. Panel (a) shows the total agreement using the λ metric, based on all deviations irrespective of whether these are systematic or non-systematic deviations. Panel (b) shows the ratio between systematic and total deviations. Panel (c) shows the agreement based only on the unsystematic component, λu, in which any disagreement due to a systematic bias is removed.

## 3.2 Exploration of the data

As both PK and JJ downscaled products have broadly similar spatio-temporal patterns, for the sake of brevity and simplicity we are going to focus the rest of this section on only one of the two: PK. The dataset spans from 21 January 2007 to 31 December 2018, with an 8 d revisiting frequency. To provide an overview of the data, a selection of spatio-temporal subsets are displayed in Fig. 4. These represent chronosequences of various areas around the globe, with a sub-sampled revisit frequency to accentuate the seasonal dynamics of the signal. The rise and fall of productivity in seasonal vegetation are clearly visible in all cases, ranging from the surroundings of the European Alps (A) to the east of the Andean Cordillera (B) and even areas around the African Great Lakes (C). The impact of agriculture is particularly evident in the irrigated areas of the Indo-Gangetic plains (D) and within the US Corn Belt (E). Figure 4 also provides a glimpse of the data gaps that can occur over mountainous areas and lakes and during rainy periods, which can occur either due to missing coarse-scale SIF input data or missing fine-scale values of explanatory variables.

To resume the SIF patterns across the globe, Fig. 5 provides a map of the maximum SIF value encountered over the entire dataset. This provides a largely gap-free representation of the maximal productivity of terrestrial ecosystems during the 2007–2018 period. The highest values are encountered in intensive agricultural areas of South America, North America and north-eastern China.

## 3.3 Inter-comparison with TROPOMI

The PK downscaled product overlaps with TROPOMI during the period from 18 March to 23 December 2018. These two instruments sound vegetation in potentially different physiological conditions due to the particular overpass times of their respective satellite platforms (morning for GOME-2 and midday for TROPOMI). While this discrepancy may be mitigated by the daily correction factor and the light-use efficiency downscaling procedure, it may also be compounded by other factors, and a proper inter-comparison is warranted before attempting to establish an archive for TROPOMI based on the downscaled GOME-2 values. For all common records over this period the temporal agreement is thus quantified and mapped in Fig. 6a. The λ metric takes on high values in the northern temperate zones, particularly over cultivated areas such as the US Corn Belt and north-eastern China, where high maximum SIF values are observed in Fig. 5. Relatively high agreement ($\mathrm{0.5}\le \mathit{\lambda }\le \mathrm{0.75}$) is found in many highly seasonal areas, such as the Sahel, while areas with no seasonality such as deserts and tropical forests have very low agreement (λ≤0.25). This suggests that a systematic bias is largely responsible for the disagreement, which can be confirmed by mapping the fraction of systematic deviations over total deviations in Fig. 6b. This fraction is elevated everywhere, and particularly over South America, where the South Atlantic Anomaly of the magnetic field is known to deteriorate the quality of the GOME-2 SIF retrievals (e.g. ). When the systematic component is ignored using the unsystematic index of agreement λu, the map in Fig. 6c reveals that the agreement between both data streams is high. Parts with lower agreement remain in tropical forests and deserts, but even these generally have λu≥0.75.

To further illustrate the compatibility of our new downscaled GOME-2 SIF product with TROPOMI SIF retrievals, Fig. 7 shows latitudinal profiles of median SIF for different time slices. The systematic bias between both products is evident and relatively consistent. When the bias is corrected using the slope and intercept values obtained at the pixel level based on the common overlapping time series between the downscaled GOME-2 PK SIF and TROPOMI SIF, the resulting median latitudinal profiles closely match each other (reducing the root-mean-square deviations from 0.188 to 0.0328 mW m−2 sr−1 nm−1).

Figure 7Comparison of latitudinal profiles of TROPOMI SIF, downscaled GOME-2 PK SIF, and the same downscaled GOME-2 SIF but bias-corrected based on the slope and intercept obtained at pixel level over the overlapping time series.

4 Discussion

The newly downscaled SIF dataset presented here should be of general interest for Earth system science and more particularly for those studying vegetation dynamics over terrestrial ecosystems. We foresee it becoming a valuable asset to better calibrate dynamic global vegetation models (DGVMs), such as those used in TRENDY and which serve as a baseline to establish the yearly global carbon budget . This SIF dataset could also be used with data-driven approaches, such as those behind widely used products such as FLUXCOM and GLEAM , to provide enhanced datasets of the important variables such as GPP and ET.

The comparison of the downscaled PK product with TROPOMI SIF retrievals for the current overlap period shows promising results towards creating a surrogate archive for TROPOMI extending back to early 2007. Since the unmixing will never be perfect, this will never fully replace actual high-spatial-resolution SIF retrievals, but to a certain extent it should be a valuable proxy. To determine the degree to which it can be used, as a bench line for detecting changes for instance, more investigations will be necessary once longer time series of synchronous data are available. For the time being, the discrepancies between both appear to consist of a systematic bias that may originate from a combination of various reasons, ranging from differences in the retrieval approach to sun-surface-sensor geometries, acquisition times and post-filtering. Correcting this bias empirically results in large agreement between the two data sources. The spatialised pixel-wise coefficients (slope and intercept) to rescale the GOME-2 downscaled PK SIF retrievals onto the TROPOMI values are also provided along with this dataset. Having been extracted using an eigendecomposition as described in , the coefficients for this linear bias correction are reference-agnostic, i.e. neither dataset source (GOME-2 or TROPOMI) was explicitly chosen as a reference, and are thus symmetric.

The resulting downscaled products still contain gaps. As mentioned before, gaps can occur due to missing data of either the input SIF or the explanatory variables. Another reason for data gaps occurs over islands or peninsulas, which are areas were there is an insufficient number of neighbouring grid cells to perform the downscaling operation. These gaps could be filled in various ways with different levels of complexity, ranging from statistically smoothing the time series to coupling them with a model such as SCOPE  in a data assimilation system . To allow users to have the maximum level of flexibility, we have chosen here not to do any gap-filling beyond what is already being done using the spatial weighted smoothing included in the original downscaling methodology.

Revisiting the downscaling approach from created an opportunity to explore the use of distinct retrievals and different input variables. It emerges that the factor contributing most to the improvement is the downscaling procedure itself rather than the choice of the retrieval or that of the explanatory variables. This highlights the benefit and rationale of seeking to have information coming from a finer spatial support, which is more adequate for characterising the spatial fragmentation of terrestrial ecosystems. The JJ retrieval, which is known to be noisier and have a smaller bias than PK, benefited particularly from the downscaling procedure, probably due to the embedded spatial smoothing step. Regarding the change in explanatory variables, the more important change comes from substituting NDVI by NIRv, which has indeed been shown to be highly correlated to SIF . That improvement partly enabled us to tolerate the replacement of ET with the reflectance-based index NDWI, rendering the output fully independent from eddy-covariance flux towers and thus from ground-based GPP measurements. The change in the timing of the LST estimation appears marginal based on the analysis using the PK dataset, and it was thus decided to keep the afternoon overpass, as it comes from a younger satellite that may be available for longer. However, the effect of LST timing on the JJ retrieval could further be investigated in the future.

The downscaling function based on light-use efficiency theory used here could further be fine-tuned to further increase the performance. In this study and for the present version of our enhanced downscaled SIF dataset, we made a deliberate choice to maintain the overall structure of this function in order to stay compatible with the original work. A possible refinement could come by including information on the incoming photosynthetically active radiation (PAR) at the surface, if this can be obtained at the high spatial resolution necessary for downscaling (e.g. Ryu et al.2018). Having an input for PAR could help remove the implicit assumption that cloud cover has to be similar over the spatio-temporal supports of both the input downscaling explanatory variables and the retrieved GOME-2 SIF. Instead, cloud cover could be expected not to be equal due to various reasons such as differences in overpassing times, spatial resolution and compositing schemes. Another improvement should come from actually considering the fraction of escaping SIF photons relative to the total SIF emitted by the whole canopy, otherwise known as the escape ratio . This would imply using some knowledge of the canopy structure of the vegetation type in question, which may come from estimations of the clumping index based on multi-angular satellite observations . Such improvements should be explored and considered in future developments of the present dataset.

The present approach to downscaling SIF has some distinctive characteristics with respect to those proposed in other recent studies . First, the downscaling is done within some physiological constrains imposed by light-use efficiency theory rather than using a purely empirical machine-learning approach. This should ground the downscaled values within limits of plausibility and further allow for a certain degree of extrapolation. However, the present approach remains data-driven, as the model only disaggregates the SIF signal in space, but does not alter its mean value at a given time and location. The downscaling is also done within a regionalised context using local moving windows in both space (40 nearest pixels) and time (16 d). This ensures a calibration that is tailored to local conditions and further normalises explanatory variables such as ET and NDWI that need to be used in relative terms. By forcing the calibration of parameters to be done locally and independently at every time step, we further ensure that the original information contained in the SIF signal is maintained in the downscaled products. This means that the downscaled SIF is not a mere rescaling of the information provided by the high-spatial-resolution reflectance bands; this also means that any trend present in the input SIF data due to sensor degradation (e.g. Zhang et al.2018a) will necessarily be present in the respective downscaled SIF dataset. Finally, the approach also uses a particular weighted-average smoothing using an ensemble of 3×3 sets of calibrated parameters (for details, see Duveiller and Cescatti2016) that removes tiling artefacts and partially fills gaps where original SIF retrievals are deemed to be too noisy. Despite these differences, a full inter-comparison between downscaled or reconstructed SIF datasets and a benchmarking with a common and independent reference (such as GPP from flux towers) should be considered to guide future algorithmic developments and consolidate our capacity to estimate SIF and GPP from space.

5 Code availability

Code associated with this study, including the script to downscale SIF and the scripts to generate the figures in this paper, is available here: https://doi.org/10.5281/zenodo.3753521 .

6 Data availability

The dataset of daily corrected downscaled SIF described in this document has been labelled Version 2.0 to differentiate it from Version 1.0, which is the original dataset described in . These are all available in the following repository:

– Downscaled-GOME2-SIF, European Commission, Joint Research Centre (JRC) (dataset), https://doi.org/10.2905/21935FFC-B797-4BEE-94DA-8FEC85B3F9E1, PID – http://data.europa.eu/89h/21935ffc-b797-4bee-94da-8fec85b3f9e1 (last access: 11 May 2020).

The files of the Version 2.0 dataset that are described here are grouped by year in distinct NetCDF files for each of the two GOME-2 retrievals (JJ and PK). The product is distributed in an equirectangular projection with a pixel size of 0.05. The temporal coverage spans from 2007 to 2018. The temporal sampling of the product is 8 d. However, every record is based on SIF input data retrieved over a 16 d moving window. This results in a certain amount of temporal auto-correlation, as the 16 d window moves every 8 d, leaving an overlap of 8 d in each successive record. The day reported in the NetCDF file corresponds to the ninth day of the 16 d retrieval period to match the MODIS convention used in the MCD43C4 product. Along with the dataset, we also provide in a separate file the slope and intercept values at the pixel level to allow users to rescale the downscaled PK GOME-2 values to TROPOMI estimates.

Supporting data to this paper are available in this extra repository: https://doi.org/10.5281/zenodo.3753024 . This includes data to generate the figures in this paper as well as OCO-2 and TROPOMI gridded datasets that had to be generated specifically for this study.

7 Conclusions

This paper presents a new daily corrected SIF dataset with a spatial resolution of 0.05 at 8 d time steps for the period 2007–2018 based on two different retrievals of GOME-2 satellite observations. Validation with OCO-2, an independent instrument capable of estimating SIF at higher resolution with a very sparse sampling scheme, has served to identify an adequate combination of explanatory variables to reach this result. A comparison with SIF from the new TROPOMI mission indicates that this downscaled SIF could serve as an archive after a pixel-wise bias correction. As such, we foresee this new SIF dataset being a valuable asset for Earth system science in general and for monitoring vegetation productivity in particular.

This work also holds promise beyond the simple adjusting of GOME-2 towards TROPOMI. First, it could serve to produce a prior SIF dataset that could be used to optimise the SIF retrievals from the future FLEX mission. The downscaling framework could further be used to downscale TROPOMI retrievals to the spatial resolution of FLEX (∼300 m) by using explanatory variables from instruments aboard Sentinel-3. FLEX retrievals could even be downscaled to decametric spatial resolution by leveraging on combinations of multi-spectral and thermic instruments including Sentinel-2, Landsat-8 and the potential ESA candidate mission called the High Spatio-Temporal Resolution Land Surface Temperature Mission (LSTM). Finally, the framework and current dataset could also be adapted towards exploring geostationary satellite data that will be able to provide optimised sub-daily information of plant status.

Author contributions
Author contributions.

GD, FF and AC designed the study. SW, PK and CF prepared and provided input data and guidance on how to use it. FF and GD processed the data. GD analysed the data and made the figures. GD wrote the text, with contributions from all the authors.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Disclaimer
Disclaimer.

The data are provided as is, with no warranties.

Acknowledgements
Acknowledgements.

We thank Joanna Joiner for valuable discussions regarding her GOME-2 SIF retrieval and how to use it in the context of this study. G Most of the heavy calculations were done using JEODPP . The figures where realised using the ggplot2 R package with ancillary spatial vectors from Natural Earth.

Review statement
Review statement.

This paper was edited by Yuyu Zhou and reviewed by Min Chen and two anonymous referees.

References

Badgley, G., Field, C. B., and Berry, J. A.: Canopy near-infrared reflectance and terrestrial photosynthesis, Science Advances, 3, e1602244, https://doi.org/10.1126/sciadv.1602244, 2017. a, b

Baker, N. R.: Chlorophyll fluorescence: a probe of photosynthesis in vivo., Annu. Rev. Plant Biol., 59, 89–113, https://doi.org/10.1146/annurev.arplant.59.032607.092759, 2008. a

Byrd, R. H., Lu, P., Nocedal, J., and Zhu, C.: A Limited Memory Algorithm for Bound Constrained Optimization, SIAM J. Sci. Comput., 16, 1190–1208, https://doi.org/10.1137/0916069, 1995. a

Duveiller, G. and Cescatti, A.: Spatially downscaling sun-induced chlorophyll fluorescence leads to an improved temporal correlation with gross primary productivity, Remote Sens. Environ., 182, 72–89, https://doi.org/10.1016/j.rse.2016.04.027, 2016. a, b, c, d, e, f, g, h, i, j, k, l, m

Duveiller, G. and Filipponi, F.: GregDuveiller/sif-downscaling-essd: code associated with the paper Duveiller et al. 2020 ESSD (Version v1.0), Zenodo, https://doi.org/10.5281/zenodo.3753521, 2020. a

Duveiller, G., Fasbender, D., and Meroni, M.: Revisiting the concept of a symmetric index of agreement for continuous datasets, Sci. Rep., 6, 19401, https://doi.org/10.1038/srep19401, 2016. a, b, c, d

Duveiller, G., Filipponi, F., Walther, S., Kohler, P., Frankenberg, C., Guanter, L., and Cescatti, A.: Downscaled-GOME2-SIF, European Commission, Joint Research Centre (JRC) [Dataset], https://doi.org/10.2905/21935FFC-B797-4BEE-94DA-8FEC85B3F9E1, 2019. a, b, c

Duveiller, G., Filipponi, F., Köhler, P., and Frankenberg, C.: sif-downscaling-essd-supporting-data, Zenodo, https://doi.org/10.5281/zenodo.3753024, 2020. a

ESA: Land Cover CCI Product User Guide Version 2, Tech. rep., available at: http://maps.elie.ucl.ac.be/CCI/viewer/download/ESACCI-LC-Ph2-PUGv2_2.0.pdf (last access: 11 May 2020), 2017. a

Frankenberg, C., Fisher, J. B., Worden, J., Badgley, G., Saatchi, S. S., Lee, J.-E., Toon, G. C., Butz, A., Jung, M., Kuze, A., and Yokota, T.: New global observations of the terrestrial carbon cycle from GOSAT: Patterns of plant fluorescence with gross primary productivity, Geophys. Res. Lett., 38, L17706, https://doi.org/10.1029/2011GL048738, 2011. a, b, c

Gao, B.-C.: NDWI–A normalized difference water index for remote sensing of vegetation liquid water from space, Remote Sens. Environ., 58, 257–266, 1996. a

Gentine, P. and Alemohammad, S. H.: RSIF (Reconstructed Solar Induced Fluorescence): a machine‐learning vegetation product based on MODIS surface reflectance to reproduce GOME‐2 solar induced fluorescence, Geophys. Res. Lett., 45, 3136–3146, https://doi.org/10.1002/2017GL076294, 2018. a, b

Gorelick, N., Hancher, M., Dixon, M., Ilyushchenko, S., Thau, D., and Moore, R.: Google Earth Engine: Planetary-scale geospatial analysis for everyone, Remote Sens. Environ., 202, 18–27, https://doi.org/10.1016/j.rse.2017.06.031, 2017. a

Guan, K., Berry, J. A., Zhang, Y., Joiner, J., Guanter, L., Badgley, G., and Lobell, D. B.: Improving the monitoring of crop productivity using spaceborne solar-induced fluorescence, Glob. Change Biol., 22, 716–726, https://doi.org/10.1111/gcb.13136, 2015. a

Guanter, L., Frankenberg, C., Dudhia, A., Lewis, P. E., Gómez-Dans, J., Kuze, A., Suto, H., and Grainger, R. G.: Retrieval and global assessment of terrestrial chlorophyll fluorescence from GOSAT space measurements, Remote Sens. Environ., 121, 236–251, https://doi.org/10.1016/j.rse.2012.02.006, 2012. a

Guanter, L., Zhang, Y., Jung, M., Joiner, J., Voigt, M., Berry, J. A., Frankenberg, C., Huete, A. R., Zarco-Tejada, P., Lee, J.-E., Moran, M. S., Ponce-Campos, G., Beer, C., Camps-Valls, G., Buchmann, N., Gianelle, D., Klumpp, K., Cescatti, A., Baker, J. M., and Griffis, T. J.: Global and time-resolved monitoring of crop photosynthesis with chlorophyll fluorescence, P. Natl. Acad. Sci. USA, 111, E1327–E1333, https://doi.org/10.1073/pnas.1320008111, 2014. a

Huete, A., Didan, K., Miura, T., Rodriguez, E. P., Gao, X., and Ferreira, L. G.: Overview of the radiometric and biophysical performance of the MODIS vegetation indices, Remote Sens. Environ., 83, 195–213, 2002. a

Jiao, Z., Dong, Y., Schaaf, C. B., Chen, J. M., Román, M., Wang, Z., Zhang, H., Ding, A., Erb, A., Hill, M. J., Zhang, X., and Strahler, A.: An algorithm for the retrieval of the clumping index (CI) from the MODIS BRDF product using an adjusted version of the kernel-driven BRDF model, Remote Sens. Environ., 209, 594–611, https://doi.org/10.1016/J.RSE.2018.02.041, 2018. a

Joiner, J., Yoshida, Y., Vasilkov, A. P., Yoshida, Y., Corp, L. A., and Middleton, E. M.: First observations of global and seasonal terrestrial chlorophyll fluorescence from space, Biogeosciences, 8, 637–651, https://doi.org/10.5194/bg-8-637-2011, 2011. a

Joiner, J., Guanter, L., Lindstrot, R., Voigt, M., Vasilkov, A. P., Middleton, E. M., Huemmrich, K. F., Yoshida, Y., and Frankenberg, C.: Global monitoring of terrestrial chlorophyll fluorescence from moderate-spectral-resolution near-infrared satellite measurements: methodology, simulations, and application to GOME-2, Atmos. Meas. Tech., 6, 2803–2823, https://doi.org/10.5194/amt-6-2803-2013, 2013. a, b, c, d, e

Joiner, J., Yoshida, Y., Vasilkov, A., Schaefer, K., Jung, M., Guanter, L., Zhang, Y., Garrity, S., Middleton, E., Huemmrich, K., Gu, L., and Belelli Marchesini, L.: The seasonal cycle of satellite chlorophyll fluorescence observations and its relationship to vegetation phenology and ecosystem atmosphere carbon exchange, Remote Sens. Environ., 152, 375–391, https://doi.org/10.1016/j.rse.2014.06.022, 2014. a

Jung, M., Koirala, S., Weber, U., Ichii, K., Gans, F., Gustau-Camps-Valls, Papale, D., Schwalm, C., Tramontana, G., and Reichstein, M.: The FLUXCOM ensemble of global land-atmosphere energy fluxes, http://arxiv.org/abs/1812.04951, 2018. a

Köhler, P., Guanter, L., and Joiner, J.: A linear method for the retrieval of sun-induced chlorophyll fluorescence from GOME-2 and SCIAMACHY data, Atmos. Meas. Tech., 8, 2589–2608, https://doi.org/10.5194/amt-8-2589-2015, 2015. a, b, c, d, e, f

Köhler, P., Frankenberg, C., Magney, T. S., Guanter, L., Joiner, J., and Landgraf, J.: Global Retrievals of Solar-Induced Chlorophyll Fluorescence With TROPOMI: First Results and Intersensor Comparison to OCO-2, Geophys. Res. Lett., 45, 456–10, https://doi.org/10.1029/2018GL079031, 2018. a, b, c, d

Kottek, M., Grieser, J., Beck, C., Rudolf, B., and Rubel, F.: World Map of the Köppen-Geiger climate classification updated, Meteorol. Z., 15, 259–263, https://doi.org/10.1127/0941-2948/2006/0130, 2006. a

Lee, J. E., Frankenberg, C., van der Tol, C., Berry, J. A., Guanter, L., Boyce, C. K., Fisher, J. B., Morrow, E., Worden, J. R., Asefi, S., Badgley, G., and Saatchi, S.: Forest productivity and water stress in Amazonia: observations from GOSAT chlorophyll fluorescence, Proc. Roy. Soc. B-Biol. Sci., 280, 20130171, https://doi.org/10.1098/rspb.2013.0171, 2013. a

Le Quéré, C., Andrew, R. M., Friedlingstein, P., Sitch, S., Hauck, J., Pongratz, J., Pickers, P. A., Korsbakken, J. I., Peters, G. P., Canadell, J. G., Arneth, A., Arora, V. K., Barbero, L., Bastos, A., Bopp, L., Chevallier, F., Chini, L. P., Ciais, P., Doney, S. C., Gkritzalis, T., Goll, D. S., Harris, I., Haverd, V., Hoffman, F. M., Hoppema, M., Houghton, R. A., Hurtt, G., Ilyina, T., Jain, A. K., Johannessen, T., Jones, C. D., Kato, E., Keeling, R. F., Goldewijk, K. K., Landschützer, P., Lefèvre, N., Lienert, S., Liu, Z., Lombardozzi, D., Metzl, N., Munro, D. R., Nabel, J. E. M. S., Nakaoka, S., Neill, C., Olsen, A., Ono, T., Patra, P., Peregon, A., Peters, W., Peylin, P., Pfeil, B., Pierrot, D., Poulter, B., Rehder, G., Resplandy, L., Robertson, E., Rocher, M., Rödenbeck, C., Schuster, U., Schwinger, J., Séférian, R., Skjelvan, I., Steinhoff, T., Sutton, A., Tans, P. P., Tian, H., Tilbrook, B., Tubiello, F. N., van der Laan-Luijkx, I. T., van der Werf, G. R., Viovy, N., Walker, A. P., Wiltshire, A. J., Wright, R., Zaehle, S., and Zheng, B.: Global Carbon Budget 2018, Earth Syst. Sci. Data, 10, 2141–2194, https://doi.org/10.5194/essd-10-2141-2018, 2018. a

Lewis, P., Gómez-Dans, J., Kaminski, T., Settle, J., Quaife, T., Gobron, N., Styles, J., and Berger, M.: An Earth Observation Land Data Assimilation System (EO-LDAS), Remote Sens. Environ., 120, 219–235, https://doi.org/10.1016/j.rse.2011.12.027, 2012. a

Li, X. and Xiao, J.: A Global, 0.05-Degree Product of Solar-Induced Chlorophyll Fluorescence Derived from OCO-2, MODIS, and Reanalysis Data, Remote Sensing, 11, 517, https://doi.org/10.3390/rs11050517, 2019. a

Magney, T. S., Frankenberg, C., Köhler, P., North, G., Davis, T. S., Dold, C., Dutta, D., Fisher, J. B., Grossmann, K., Harrington, A., Hatfield, J., Stutz, J., Sun, Y., and Porcar‐Castell, A.: Disentangling changes in the spectral shape of chlorophyll fluorescence: Implications for remote sensing of photosynthesis, J. Geophys. Res.-Biogeo., 124, 1491–1507, https://doi.org/10.1029/2019JG005029, 2019. a

Mattiuzzi, M. and Detsch, F.: MODIS: Acquisition and Processing of MODIS Products, available at:https://cran.r-project.org/package=MODIS (last access: 11 May 2020), 2018. a

Maxwell, K. and Johnson, G.: Chlorophyll fluorescence – a practical guide, J. Exp. Bot., 51, 659–668, https://doi.org/10.1093/jexbot/51.345.659, 2000. a

Meroni, M., Rossini, M., Guanter, L., Alonso, L., Rascher, U., Colombo, R., and Moreno, J.: Remote sensing of solar-induced chlorophyll fluorescence: Review of methods and applications, Remote Sens. Environ., 113, 2037–2051, https://doi.org/10.1016/j.rse.2009.05.003, 2009. a

Miralles, D. G., Holmes, T. R. H., De Jeu, R. A. M., Gash, J. H., Meesters, A. G. C. A., and Dolman, A. J.: Global land-surface evaporation estimated from satellite-based observations, Hydrol. Earth Syst. Sci., 15, 453–469, https://doi.org/10.5194/hess-15-453-2011, 2011. a

Monteith, J. L.: Climate and efficiency of crop production in Britain, Philosophical Transactions of the Royal Society B: Biological Sciences, 281, 277–294, 1977. a

Mu, Q., Zhao, M., and Running, S. W.: Improvements to a MODIS global terrestrial evapotranspiration algorithm, Remote Sens. Environ., 115, 1781–1800, https://doi.org/10.1016/j.rse.2011.02.019, 2011. a

Parazoo, N. C., Bowman, K., Frankenberg, C., Lee, J. E., Fisher, J. B., Worden, J., Jones, D. B., Berry, J., Collatz, G. J., Baker, I. T., Jung, M., Liu, J., Osterman, G., O'Dell, C., Sparks, A., Butz, A., Guerlet, S., Yoshida, Y., Chen, H., and Gerbig, C.: Interpreting seasonal changes in the carbon balance of southern Amazonia using measurements of XCO2 and chlorophyll fluorescence from GOSAT, Geophys. Res. Lett., 40, 2829–2833, https://doi.org/10.1002/grl.50452, 2013. a

Parazoo, N. C., Barnes, E., Worden, J., Harper, A. B., Bowman, K. B., Frankenberg, C., Wolf, S., Litvak, M., and Keenan, T. F.: Influence of ENSO and the NAO on terrestrial carbon uptake in the Texas-northern Mexico region, Global Biogeochem. Cy., 29, 1247–1265, https://doi.org/10.1002/2015GB005125, 2015. a

Porcar-Castell, A., Tyystjärvi, E., Atherton, J., van der Tol, C., Flexas, J., Pfündel, E. E., Moreno, J., Frankenberg, C., and Berry, J. A.: Linking chlorophyll a fluorescence to photosynthesis for remote sensing applications: mechanisms and challenges, J. Exp. Bot., 65, 4065–95, https://doi.org/10.1093/jxb/eru191, 2014. a

Running, S., Mu, Q., and Zhao, M.: MOD16A2 MODIS/Terra Net Evapotranspiration 8-Day L4 Global 500m SIN Grid V006 [Data set], https://doi.org/10.5067/MODIS/MOD16A2.006, 2017. a

Running, S. W., Nemani, R. R., Heinsch, F. A., Zhao, M., Reeves, M., and Hashimoto, H.: A Continuous Satellite-Derived Measure of Global Terrestrial Primary Production, BioScience, 54, 547–560, https://doi.org/10.1641/0006-3568(2004)054[0547:ACSMOG]2.0.CO;2, 2004. a

Ryu, Y., Jiang, C., Kobayashi, H., and Detto, M.: MODIS-derived global land products of shortwave radiation and diffuse and total photosynthetically active radiation at 5 km resolution from 2000, Remote Sens. Environ., 204, 812–825, https://doi.org/10.1016/j.rse.2017.09.021, 2018. a

Ryu, Y., Berry, J. A., and Baldocchi, D. D.: What is global photosynthesis? History, uncertainties and opportunities, Remote Sens. Environ., 223, 95–114, https://doi.org/10.1016/j.rse.2019.01.016, 2019. a

Sayer, A. M., Hsu, N. C., Bettenhausen, C., Jeong, M. J., and Meister, G.: Effect of MODIS terra radiometric calibration improvements on Collection 6 Deep blue aerosol products: Validation and terra/aqua consistency, J. Geophys. Res., 120, 157–12, https://doi.org/10.1002/2015JD023878, 2015. a

Schaaf, C. and Wang, Z.: MCD43C4 MODIS/Terra+Aqua BRDF/Albedo Nadir BRDF-Adjusted Ref Daily L3 Global 0.05Deg CMG V006 [Data set], https://doi.org/10.5067/MODIS/MCD43C4.006, 2015. a

Schaaf, C. B., Gao, F., Strahler, A. H., Lucht, W., Li, X., Tsang, T., Strugnell, N. C., Zhang, X., Jin, Y., Muller, J.-P., Lewis, P., Barnsley, M., Hobson, P., Disney, M., Roberts, G., Dunderdale, M., Doll, C., d'Entremont, R. P., Hu, B., Liang, S., Privette, J. L., and Roy, D.: First operational BRDF, albedo nadir reflectance products from MODIS, Remote Sens. Environ., 83, 135–148, https://doi.org/10.1016/S0034-4257(02)00091-3, 2002. a

Sitch, S., Friedlingstein, P., Gruber, N., Jones, S. D., Murray-Tortarolo, G., Ahlström, A., Doney, S. C., Graven, H., Heinze, C., Huntingford, C., Levis, S., Levy, P. E., Lomas, M., Poulter, B., Viovy, N., Zaehle, S., Zeng, N., Arneth, A., Bonan, G., Bopp, L., Canadell, J. G., Chevallier, F., Ciais, P., Ellis, R., Gloor, M., Peylin, P., Piao, S. L., Le Quéré, C., Smith, B., Zhu, Z., and Myneni, R.: Recent trends and drivers of regional sources and sinks of carbon dioxide, Biogeosciences, 12, 653–679, https://doi.org/10.5194/bg-12-653-2015, 2015. a

Soille, P., Burger, A., De Marchi, D., Kempeneers, P., Rodriguez, D., Syrris, V., and Vasilev, V.: A versatile data-intensive computing platform for information retrieval from big geospatial data, Future Gene. Comp. Sy., 81, 30–40, https://doi.org/10.1016/J.FUTURE.2017.11.007, 2018. a

Sun, Y., Gu, L., Dickinson, R. E., Norby, R. J., Pallardy, S. G., and Hoffman, F. M.: Impact of mesophyll diffusion on estimated global land CO2 fertilization, P. Natl. Acad. Sci. USA, 111, 15774–15779, https://doi.org/10.1073/pnas.1418075111, 2014. a

Sun, Y., Frankenberg, C., Jung, M., Joiner, J., Guanter, L., Köhler, P., and Magney, T.: Overview of Solar-Induced chlorophyll Fluorescence (SIF) from the Orbiting Carbon Observatory-2: Retrieval, cross-mission comparison, and global monitoring for GPP, Remote Sens. Environ., 209, 808–823, https://doi.org/10.1016/j.rse.2018.02.016, 2018. a, b

van der Tol, C., Verhoef, W., Timmermans, J., Verhoef, A., and Su, Z.: An integrated model of soil-canopy spectral radiances, photosynthesis, fluorescence, temperature and energy balance, Biogeosciences, 6, 3109–3129, https://doi.org/10.5194/bg-6-3109-2009, 2009. a

Veefkind, J., Aben, I., McMullan, K., Förster, H., de Vries, J., Otter, G., Claas, J., Eskes, H., de Haan, J., Kleipool, Q., van Weele, M., Hasekamp, O., Hoogeveen, R., Landgraf, J., Snel, R., Tol, P., Ingmann, P., Voors, R., Kruizinga, B., Vink, R., Visser, H., and Levelt, P.: TROPOMI on the ESA Sentinel-5 Precursor: A GMES mission for global observations of the atmospheric composition for climate, air quality and ozone layer applications, Remote Sens. Environ., 120, 70–83, https://doi.org/10.1016/j.rse.2011.09.027, 2012. a

Walther, S., Voigt, M., Thum, T., Gonsamo, A., Zhang, Y., Köhler, P., Jung, M., Varlagin, A., and Guanter, L.: Satellite chlorophyll fluorescence measurements reveal large-scale decoupling of photosynthesis and greenness dynamics in boreal evergreen forests, Glob. Change Biol., 22, 2979–2996, https://doi.org/10.1111/gcb.13200, 2016. a

Walther, S., Guanter, L., Heim, B., Jung, M., Duveiller, G., Wolanin, A., and Sachs, T.: Assessing the dynamics of vegetation productivity in circumpolar regions with different satellite indicators of greenness and photosynthesis, Biogeosciences, 15, 6221–6256, https://doi.org/10.5194/bg-15-6221-2018, 2018. a

Walther, S., Duveiller, G., Jung, M., Guanter, L., Cescatti, A., and Camps-Valls, G.: Satellite Observations of the Contrasting Response of Trees and Grasses to Variations in Water Availability, Geophys. Res. Lett., 46, 1429–1440, https://doi.org/10.1029/2018GL080535, 2019. a

Wan, Z.: New refinements and validation of the MODIS Land-Surface Temperature/Emissivity products, Remote Sens. Environ., 112, 59–74, https://doi.org/10.1016/j.rse.2006.06.026, 2008. a

Wan, Z., Hook, S., and Hulley, G.: MOD11C2 MODIS/Terra Land Surface Temperature/Emissivity 8-Day L3 Global 0.05Deg CMG V006 [Data set], NASA EOSDIS LP DAAC, https://doi.org/10.5067/MODIS/MOD11C2.006, 2015a. a

Wan, Z., Hook, S., and Hulley, G.: MYD11C2 MODIS/Aqua Land Surface Temperature/Emissivity 8-Day L3 Global 0.05Deg CMG V006 [Data set], https://doi.org/10.5067/MODIS/MYD11C2.006, 2015b. a

Yoshida, Y., Joiner, J., Tucker, C., Berry, J. A., Lee, J.-E., Walker, G., Reichle, R., Koster, R., Lyapustin, A., and Wang, Y.: The 2010 Russian drought impact on satellite measurements of solar-induced chlorophyll fluorescence: Insights from modeling and comparisons with parameters derived from satellite reflectances, Remote Sens. Environ., 166, 163–177, https://doi.org/10.1016/j.rse.2015.06.008, 2015.  a

Yu, L., Wen, J., Chang, C. Y., Frankenberg, C., and Sun, Y.: High Resolution Global Contiguous SIF of OCO-2, Geophys. Res. Lett., 46, 1449–1458, https://doi.org/10.1029/2018GL081109, 2018. a, b

Zhang, Y., Joiner, J., Gentine, P., and Zhou, S.: Reduced solar-induced chlorophyll fluorescence from GOME-2 during Amazon drought caused by dataset artifacts, Glob. Change Biol., 24, 2229–2230, https://doi.org/10.1111/gcb.14134, 2018a. a

Zhang, Y., Joiner, J., Alemohammad, S. H., Zhou, S., and Gentine, P.: A global spatially contiguous solar-induced fluorescence (CSIF) dataset using neural networks, Biogeosciences, 15, 5779–5800, https://doi.org/10.5194/bg-15-5779-2018, 2018. a, b