Chlorophyll fluorescence remote sensing from space in scattering atmospheres : implications for its retrieval and interferences with atmospheric CO 2 retrievals

Chlorophyll fluorescence remote sensing from space in scattering atmospheres: implications for its retrieval and interferences with atmospheric CO2 retrievals C. Frankenberg, C. O’Dell, L. Guanter, and J. McDuffie Jet Propulsion Laboratory, California Institute of Technology, Pasadena, USA Colorado State University, Fort Collins, CO, USA Atmospheric, Oceanic and Planetary Physics, University of Oxford, UK


Introduction
Both the Greenhouse Gases Observing Satellite (GOSAT) (Hamazaki et al., 2005;Kuze et al., 2009) and the Orbiting Carbon Observatory (OCO) (Crisp et al., 2004) aim at retrieving atmospheric greenhouse gas column averaged mixing ratios (denoted Figures X CO 2 or X CH 4 ) with high enough accuracy and precision to improve our current understanding of land-atmosphere fluxes on regional scales.OCO, targeting CO 2 exclusively, suffered a launch failure in early 2009, but an instrument copy is being built and the launch of OCO-2 is planned for 2014.
A typical full-physics algorithm for the retrieval of X CO 2 (B ösch et al., 2006;Butz et al., 2009;O'Dell et al., 2011) concurrently employs three spectral bands, centered around 0.76 µm (O 2 A-band), 1.61 µm (weak CO 2 band) and 2.06 µm (strong CO 2 band).By using this multi-channel approach, X CO 2 , surface albedos as well as aerosol properties can be retrieved concurrently.The aerosol properties deserve special interest in this study as they are mostly constrained by the strong absorption features of atmospheric O 2 at 0.76 µm as well as strong CO 2 features at 2.06 µm.Since atmospheric scattering can enhance or shorten the light-path of recorded photons, it will affect the X CO 2 retrievals if retrieved scattering properties are biased.Frankenberg et al. (2011a) found that chlorophyll fluorescence (F s ) cannot be unambiguously distinguished from the effect of scattering on the depth and shape of atmospheric O 2 absorption features in the 0.76 µm range.However, the potential propagated impact on X CO 2 retrievals was not explicitly quantified.The goals of this paper are twofold: first, to quantify and mitigate the potential impact of chlorophyll fluorescence on X CO 2 retrievals, and second, to quantify how accurately TOA (top-of-atmosphere) fluorescence can be retrieved and how it relates to fluorescence at the surface level.While the first goal is of a technical nature and relates to the CO 2 remote sensing community, the second one has implications regarding the applicability of remotely sensed chlorophyll fluorescence, especially given that real retrievals are now being successfully performed using GOSAT data (Joiner et al., 2011;Frankenberg et al., 2011b;Guanter et al., 2012).
The paper is structured as follows: Section 2 gives an overview of the chlorophyll fluorescence signal and its application in plant physiology and vegetation remote sensing.Section 3 describes the implementation of the fluorescence emission in a radiative transfer model for simulations as well as in a a full-physics retrieval algorithm.Section 4 Figures

Back Close
Full shows simulation results as to how fluorescence impacts scattering and X CO 2 retrievals.Section 5 focuses on the retrieval of the chlorophyll fluorescence signal itself and discusses how F s retrieved at TOA relates to the emission from the plant at surface level in the presence of scattering.It also provides an overview on how future satellites with similar measurement principles can improve upon current and planned ones.The audience for this manuscript may be rather diverse, with interest either in X CO 2 or fluorescence retrievals.Thus, we tried to keep Sect.4 (focused on X CO 2 biases) and Sect. 5 (focused on chlorophyll fluorescence itself) as self-contained as possible in case the general reader wants to skip sections.

Chlorophyll fluorescence
During photosynthesis, visible solar energy absorbed by chlorophyll can either be used for carbon fixation, be dissipated into heat, or be re-emitted via fluorescence at longer wavelengths in the 660-800 nm window.This so-called solar-induced chlorophyll fluorescence (Krause and Weis, 1991;Baker, 2008, and references therein) thus offers a very direct measure of photosynthetic activity.At the surface level, the fluorescence emission (F s ) adds a small offset of up to about 10 W m −2 sr −1 µm −1 (Entcheva Campbell et al., 2008) in the two emission peaks around 690 and 740 nm, encompassing the strongly saturated O 2 A-band around 765 nm.The broad-band emission leads to a filling-in of solar and telluric absorption features (Plascyk and Gabriel, 1975), and current retrieval techniques at the ground and airborne scales mostly focus on using the O 2 absorption features for its retrieval.See Meroni et al. (2009) for a literature overview of common retrieval techniques for fluorescence signals used to study the photosynthetic apparatus (e.g.Krause andWeis, 1984, 1991;Flexas et al., 2002;Moya et al., 2004;Rascher et al., 2009) and stress related changes in gross primary production (GPP) (Damm et al., 2010;Daumard et al., 2010).
In Frankenberg et al. (2011a) we have shown that the use of O 2 lines in fluorescence retrievals cannot be readily applied if only the O 2 A-band is measured.However, Introduction

Conclusions References
Tables Figures

Back Close
Full Fraunhofer lines (absorption features in the solar atmosphere) can be used to disentangle the fluorescence emission from scattering effects, a technique used in the quantification of chlorophyll emission from GOSAT data (Joiner et al., 2011;Frankenberg et al., 2011b).
To put the analysis in the following sections into perspective, Fig. 1 shows a 2-yr average of chlorophyll fluorescence in terms of % of continuum level radiance.Regional averages can reach 1.5 % with sharp gradients occurring at transitions to barren surfaces.Systematic biases caused by the neglect of fluorescence in space-based retrievals of X CO 2 will thus resemble this spatial structure, which, in turn, resembles GPP.

Fluorescence implementation in simulations and retrievals
From a radiative transfer perspective, chlorophyll fluorescence can be considered as a spectrally smooth source term at the surface.As will be explained in the following, we implemented the chlorophyll fluorescence term directly in the radiative transfer solver in forward model simulations.However, we ignore scattering (but not absorption) effects on the fluorescence source term in the implementation in the retrieval scheme, greatly simplifying its parameterization.

Accurate implementation in a radiative transfer scheme for simulation purposes
In order to provide realistic simulations of the O 2 A-band at TOA impacted by fluorescence, we implemented the hitherto neglected chlorophyll emission into the satellite simulator used for GOSAT and OCO-2 studies (O'Dell et al., 2011).The simulator models spectra at the TOA using realistic and diverse aerosol and cloud distributions, surface and instrument properties, coupled with an accurate radiative transfer model (O'Brien et al., 2009).The main purpose of the simulator is to quantify systematic Introduction

Conclusions References
Tables Figures

Back Close
Full errors in retrievec X CO 2 , as the retrieval approach cannot capture the full variability of atmospheric and surface properties, most importantly clouds and aerosols.
For the sake of simplicity, we model fluorescence as an isotropic, unpolarized source term at the surface.The spectral shape of the surface emission, given for instance in Guanter et al. (2010), is modeled as two co-added Gaussians: where F s,755 is the fluorescence intensity at 755 nm.The parameters in Eq. ( 1) are given in Table 1.This model has the benefit that it enables future studies using other wavelength regions and that the spectral shape in the A-band is fairly realistic.This allows for the assessment of errors due to the assumption of linearity used in the retrieval described in Sect.3.2.
The fluorescence emission from the surface has to pass through the atmosphere to be detected by the satellite at TOA.It will undergo absorption by gases, Rayleigh scattering, as well as possible extinction by clouds and aerosols.These effects are accurately modeled using a 24-stream adding-doubling radiative transfer model (see e.g.Liou, 2002).In this radiative transfer approach, global scalar reflectance and transmittance matrices are derived for the entire atmosphere using doubling for individual layers, and adding in order to combine all the atmospheric layers together.The surface emission at TOA is modeled as: where T a (R a ) is the atmospheric transmittance (reflectance) matrix from surface to space, R g is the surface reflectance matrix, F surf is the isotropic fluorescence emission vector at the surface, and F TOA is the corresponding vector at TOA. Polarization effects of the atmosphere on the unpolarized fluorescence emission are assumed to be small.Introduction

Conclusions References
Tables Figures

Back Close
Full As for the scattered solar irradiance, the method of "low-streams interpolation" (O'Dell, 2010) is used to accelerate the radiative transfer.
The spatial distribution of F s in the simulator is based upon 0.5 • × 0.5 • degree GPP distribution from Beer et al. (2010) and the empirical linear relationship between F s and GPP derived in Frankenberg et al. (2011b).For the following analysis, we simulated three polar orbits of a GOSAT-like instrument in nadir viewing geometry.This is a subset of the orbits used in a recent simulation-based study of satellite-based X CO 2 retrievals (O'Dell et al., 2011).

Simplified implementation for retrieval purposes
A typical X CO 2 full-physics retrieval algorithm uses a radiative transfer model to simulate radiances at TOA as a function of atmospheric parameters such as surface albedo as well as absorption and scattering profiles (including scattering phase functions).
The Jacobians of the simulated radiances are computed with respect to retrieval parameters such as the CO 2 profile, surface albedo as well as aerosol optical depth and profile.Given a true set of radiances measured by a space-borne spectrometer, the retrieval parameters can be inverted using a non-linear least squares approach minimizing the differences between modeled and measured radiances.The most common setup makes concurrent use of three bands in a combined retrieval: the O 2 A-band at 760 nm as well as weak and strong CO 2 bands at 1.6 and 2 µm, respectively.Chlorophyll fluorescence, if neglected, will impact the radiances in one of the retrieval windows, viz. the O 2 A-band.We now describe modifications to the standard ACOS/OCO-2 X CO 2 retrieval algorithm to include chlorophyll fluorescence.This algorithm has been described at length in O' Dell et al. (2011) and references therein.Briefly, this optimal estimation-based algorithm (Rodgers, 2000) attempts to minimize differences between observed and modeled spectra in the three spectral bands described above, moderated by a side constraint involving prior knowledge.The algorithm fits simultaneously for a profile of CO 2 , surface pressure, a linearly-varying surface albedo in each of the three fitted Introduction

Conclusions References
Tables Figures

Back Close
Full bands, spectral wavelength offset in each band, and profiles of four scattering particle types with fixed microphysical and optical properties.These include a water cloud type, an ice cloud type, and two different aerosol types.The idea is that by mixing together these different types appropriately, virtually any type of scattering situation can be modeled.When applied to actual GOSAT data, it was found that an additional state vector parameter representing a zero-level offset in the O 2 A-band was also needed to correct for an instrument calibration problem (Crisp et al., 2012;Butz et al., 2011).
We implemented the chlorophyll emission as two additional state vector elements in the ACOS retrieval algorithm, as follows.We simplify the simulation of the chlorophyll emission inside the retrieval forward model by decoupling it entirely from the radiative transfer mode.Instead of calculating the full scattering effect on the chlorophyll emission observed at TOA, we solely model the spectral shape as well as the absorption features by O 2 .The fluorescence is then treated as a simple additive term to the radiance calculated by the radiative transfer code before convolution with the instrument line shape is performed.The fluorescence term at the surface (using nm as λ units) in the retrieval is: where the two state vector elements F s,755 and s represent the chlorophyll emission at 755 nm and its spectral slope, respectively.Compared to a simpler polynomial representation of F s (λ) surf , this parameterization allows a scaling of a predefined spectral shape with some freedom for adjusting the slope itself.This model explicitly assumes a linear spectral shape over the relatively narrow O 2 A-band, an assumption that is not made in the F s simulation (Sect.3.1).At TOA, the modeled fluorescence ignoring scattering reads where µ represents the cosine of the viewing zenith angle and τ O 2 (λ) the vertical optical thickness of O 2 .F s (λ) TOA is simply added to the radiative transfer scheme implemented Introduction

Conclusions References
Tables Figures

Back Close
Full in the ACOS/OCO-2 full-physics retrieval algorithm, largely facilitating the calculation of Jacobians with respect to F s,755 and s.

Interference errors between scattering properties, X CO 2 and chlorophyll fluorescence
The impact of chlorophyll fluorescence on the retrieval of X CO 2 is rather complex, as it is propagated through its interference with retrieved scattering properties and, potentially, surface pressure.The impact can vary not only scene-by-scene but also between different retrieval algorithms, depending on how scattering and surface pressure are treated.However, there always is a propagation of biases from the O 2 A-band into X CO 2 , unless the retrieval is set up to entirely decouple the O 2 A-band from the CO 2 retrieval, which would forfeit the initial purpose of the O 2 A-band in CO 2 retrievals.In this study, we focus on the ACOS/OCO-2 3-band retrieval algorithm described above.
To gain a deeper understanding of the retrieval problem, we perform simulations and retrievals of varying complexity, listed in Table 2.We first focus on the case where when fluorescence is entirely ignored in the X CO 2 retrieval step.To our knowledge, this is currently the standard approach in all retrievals that utilize the O 2 A-band (e.g.Schneising et al., 2008;Butz et al., 2011;Parker et al., 2011;O'Dell et al., 2011;Yoshida et al., 2011), and is thus highly relevant.We then analyze the results of a retrieval with fluorescence included as part of the state vector according to Sect.3.2.To quantify systematic biases, we run all retrievals twice: once using simulations with no fluorescence added and once with fluorescence included.In the following, all biases caused by fluorescence are reported as the difference between the two runs.Both runs are performed using noise-less simulations while a realistic GOSAT-like noise is assumed in the measurement error covariance matrix for the retrieval.Introduction

Conclusions References
Tables Figures

Back Close
Full

Case I: ignoring fluorescence
As stated previously, most retrieval algorithms of CO 2 which employ the O 2 A-band currently entirely ignore chlorophyll fluorescence.They will thus be subject to systematic biases depending on their specific retrieval setup as well as instrument characteristics, most importantly spectral resolution.Here and in the following, we focus on simulations performed for GOSAT spectral resolution.Contrary to SCIAMACHY, for instance, Fraunhofer lines can be better resolved and will thus reduce interference errors.

Aerosol-free atmosphere
In the absence of aerosol scattering in the atmosphere, chlorophyll fluorescence nevertheless impacts the retrieval of X CO 2 , as the reduction in fractional depth of O 2 absorption lines can resemble a change in surface pressure.Case #1a represents this most simple case, viz. a completely aerosol free atmosphere assumed in both simulations and retrievals.The only free variables in the retrieval related to the O 2 A-band are thus spectral albedo, spectral shift, temperature offset and surface pressure (given a 4 hPa 1−σ prior uncertainty).As F s fills in the O 2 lines, surface pressure is consequently underestimated with increasing F s .
Figure 3 shows the errors propagating into CO 2 via the negative surface pressure bias with increasing fluorescence.At about 1 % relative F s , the overestimation in X CO 2 amounts to 1 ppm.Changes in retrieved surface pressure are about 2.5 hPa, slightly higher than current uncertainties in surface pressure from meteorological models (1-Introduction

Conclusions References
Tables Figures

Back Close
Full

Aerosol-laden atmosphere
In the presence of scattering in the atmosphere (both in simulation and retrieval), the impact on X CO 2 is more complex as various fitting parameters for atmospheric aerosols and/or clouds have to be included in the state vector.These parameters have been shown to interfere with fluorescence as well (Frankenberg et al., 2011a).In a scattering atmosphere, post-processing has to be performed, discarding non-converging and bad retrievals according to O'Dell et al. ( 2011), with a retrieved AOD threshold of < 0.2. Figure 4 shows results from case #2a, where both aerosols as well as surface pressure are free retrieval parameters.As expected, the results are more mixed, with higher scatter in the X CO 2 biases.The overall fluorescence-induced bias in X CO 2 is somewhat lower than in case #1a, mainly because aerosols can now partially account for the fluorescence term.The propagation of aerosols biases is, however, not as direct as in the case of surface pressure, thus explaining the smaller X CO 2 as well surface pressure changes.Intuitively, one would expect that retrieved aerosols increase with increasing fluorescence as higher aerosol backscatter typically reduces the absorption depth of O 2 lines.However, at high surface albedos, the opposite can also be true as aerosols can as also lengthen the light path.This is probably the reason why changes in retrieved aerosol properties are not particularly systematic.Overall, we see an insignificant decrease in retrieved total aerosol optical depth but a significant redistribution of aerosols to higher atmospheric layers, which more efficiently shorten the atmospheric light-path.
To investigate the isolated impact of aerosols, we performed case #3a, where surface pressure is assumed known, an assumption made by some CO 2 retrievals (e.g., Butz et al., 2011).Figure 5 shows these results where naturally only aerosols as well as X CO 2 change systematically.As the error propagation through surface pressure is eliminated, the bias induced in X CO 2 is lower, though still apparent with an average slope of 0.55 ppm per W m −2 sr −1 µm −1 .However, in both case #2a and case #3a, substantial Introduction

Conclusions References
Tables Figures

Back Close
Full random scatter in retrieved X CO 2 and AOD is induced as well.The impact on aerosols is naturally larger when the surface pressure is fixed.
Looking at cases #1a-3a, one can conclude that the neglect of fluorescence in our particular retrieval (and instrument characteristics) can result in systematic biases of ≈ 0.5−1 ppm per % fluorescence contribution to the continuum radiance.Biases in retrieved CO 2 are thus within ≈ 1.5 ppm but vary frequently in the spatial domain (see Fig. 1).While these biases are still below 1 % of the CO 2 column amount, they exhibit spatial and temporal correlations with GPP, which is why they must not be neglected.

Case II: fitting fluorescence
As we have seen in the previous section, fluorescence can cause systematic biases in retrieved X CO 2 , surface pressure and aerosols.As interferences shown in Frankenberg et al. (2011a) point to problems in fitting all variables simultaneously, we investigate the impact of including fluorescence as a state vector element, as described in Sect.3.2.This is done from the perspective of reducing systematic biases in X CO 2 , rather than to retrieve fluorescence accurately.Several studies have already fitted fluorescence using GOSAT data (Frankenberg et al., 2011a,b;Joiner et al., 2011) but these have so far focussed solely on Fraunhofer lines, avoiding telluric O 2 absorption features.These are typically used in ground-based studies and suggested for the FLEX mission but are, as will be corroborated in the following section, problematic if observed from TOA.
In the pure Rayleigh case, where aerosols are present neither in the simulations nor in the retrieval, the fit can virtually eliminate the systematic error previously existent in X CO 2 .As shown in Fig. 6, the X CO 2 bias changes from 1.07 to 0.07 ppm % −1 compared to case #1a.Also fluorescence is fitted quite accurately, with the caveat that these are noise-free simulations.However, the retrieval is not able to reproduce the true fluorescence 1 : 1, potentially because of the neglect of Rayleigh scattering in the retrieval forward modeling of fluorescence.Results, especially in terms of scatter unrelated to noise, become more complex when aerosols are explicitly included.Figure 7 shows results for case #2b, where Introduction

Conclusions References
Tables Figures

Back Close
Full aerosols are present in the simulations, and aerosols as well as surface pressure are fitted in the retrieval.The systematic bias in X CO 2 as well as surface pressure is not significant and seems to be dominated by few outliers.However, there is high scatter in the retrieved fluorescence parameter (noting that this is an accuracy, not a precision error as simulations are noise-free).There is also some scatter induced in X CO 2 but most deviations are below 1 ppm and appear unsystematic.Perfect knowledge of surface pressure (#3b, Fig. 8) slightly improves the situation, with a slope of close to 1 in retrieved vs. true fluorescence.However, the scatter in retrieved F s in case #3b is still very high, but systematic biases in X CO 2 as a function of signal level are small.The inclusion of fluorescence as fitting parameters caused some outliers in retrieved X CO 2 at low true F s values though.
In summary, the results confirm what has been alluded to in Frankenberg et al. ( 2011a): (a) there is a strong but hard to generalize interference between fluorescence, surface pressure and aerosol properties; (b) if ignored, systematic biases in X CO 2 on the order of 1 ppm can occur.Inclusion of fluorescence as a state vector element partially eliminates the bias but is not ideal if chlorophyll fluorescence itself is the target quantity.Using the potentially added information contained in the O 2 absorption structures seems to do more harm than good for the fluorescence retrieval for interferences are introduced.If the retrieval from space would be purely based on the O 2 A-band (as is done frequently in ground-based studies), ancillary information on aerosols from other sensors thus won't necessarily be accurate enough to provide a useful constraint, as the impact on the O 2 lines is a complex function of AOD, aerosol type, height and surface albedo.The fact that we didn't see very systematic changes in retrieved aerosol parameters when ignoring fluorescence supports this conclusion.Similar concerns would be valid for airborne studies, if aerosol scattering between the surface and flight altitude exists.A combination of O 2 A and B-bands, such as envisioned in the FLEX mission concept, would be preferred if spectral resolution is not high enough to exploit Fraunhofer line features (Guanter et al., 2010;Sanghavi et al., 2011).Introduction

Conclusions References
Tables Figures

Back Close
Full  2011b) discovered a zero-level offset in GOSAT O 2 A-band spectra, which was subsequently confirmed by an independent analysis of the JAXA GOSAT team.This adds another dimension to the problem for this particular case for a zero level offset in FTS systems causes a rather flat offset to the true radiance at all wavelengths.For Fraunhofer lines, this effect is thus virtually indistinguishable from a chlorophyll fluorescence signal as both simply add an additive term to the spectrum.Within O 2 lines, however, the effect is different from fluorescence for it adds a spectrally flat offset as opposed to fluorescence, where the spectrally flat emission at the surface level is partially reabsorbed by O 2 absorption along the path to the detector.Hence, both signals have an identical impact in Fraunhofer lines and the continuum but differ in O 2 absorption features.
For previous dedicated fluorescence retrievals using Fraunhofer lines, the zero-level offset contribution was parameterized as a function of average signal level (representative for the zero-path difference readout in FTS systems) and later subtracted from the retrieved offset in order to isolate the fluorescence term (Frankenberg et al., 2011b;Guanter et al., 2012).Here, we investigate the impact of fluorescence on the fullphysics retrieval in view of the apparent zero-level offset.There are two questions to be answered: (a) What is the potential bias induced by fluorescence on current operational retrievals fitting the zero-level offset but not fluorescence and (b) will a simultaneous retrieval of fluorescence and zero-level offset be stable and also minimize biases in X CO 2 ?For this purpose, we ran two additional tests, where we (a) mimicked the current ACOS operational retrieval by fitting aerosols, surface pressure and zero-level offset but not fluorescence and (b) fitted fluorescence in addition.
Figure 9 shows case (a) where we fit for zero-level offset but ignore fluorescence.
As expected, fitting the zero-level offset somewhat diminishes the bias in X CO 2 , with a slope reduced from 0.86 to 0.58 ppm % −1 .However, the bias is still positive and not Introduction

Conclusions References
Tables Figures

Back Close
Full eliminated as in the case of fitting fluorescence.It can be concluded that a remaining small bias exists in most current X CO 2 retrievals.
Figure 10 shows case b) where we fit for both fluorescence and zero-level offset.The retrieval appears as stable as when zero-level offset is not fitted and also the bias in X CO 2 is diminished.It can be concluded that current X CO 2 retrievals can include both fluorescence and zero-level offset in their retrieval, thereby eliminating the bias through fluorescence even in the presence of systematic instrument errors.The main reason why the fit can partially disentangle fluorescence from zero-level offset is their strong difference in the impact on deeply saturated O 2 lines (where zero-level offset has the highest fractional impact but fluorescence is almost entirely reabsorbed in the atmosphere).

Chlorophyll fluorescence as primary retrieval target
In the previous section, we largely focussed on the impact of fluorescence on X CO 2 retrievals.Here, we focus on the implications of atmospheric scattering by clouds and aerosols on the fluorescence retrieval itself.The most appropriate way to deal with fluorescence is to avoid the interferences altogether by allowing for a processing step using Fraunhofer lines only, as in Frankenberg et al. (2011a).In the following, this method is used for fluorescence retrievals as the primary target.

Implications of scattering on retrieved chlorophyll fluorescence at top-ofatmosphere
We have shown previously (Frankenberg et al., 2011a) that pure retrieval interference errors between fluorescence and scattering properties can be eliminated in the retrieval if only Fraunhofer lines are used.The retrieved fluorescence signal thus accurately represents the true signal at TOA (see Fig. 11).Introduction

Conclusions References
Tables Figures

Back Close
Full However, scattering impacts the propagation of the fluorescence emission from the surface to TOA, hence TOA fluorescence doesn't necessarily represent fluorescence at the plant level.Given an extinction optical thickness of τ, which includes gas absorption, Rayleigh scattering and extinction by clouds and aerosols, the direct beam fluorescence would be reduced by e −τ/µ .However, the fluorescence signal is partially conserved in the diffuse radiation field.This essentially smears out the effective FOV representative for the F s signal at TOA, which should be negligible for footprint sizes on the order of multiple kilometers.
To quantify the potential reduction in fluorescence at TOA, we computed the fraction of fluorescence at TOA vs. surface.Figure 12 depicts this fraction using the simulator runs, which employed a nadir observation angle (µ = 1).Please note that this only holds outside of atmospheric absorption bands, i.e.only if the Fraunhofer line method is applied.At low optical depths, the ratio is basically 1.Most importantly, the ratio is still well above 0.8 for τ < 1 and mostly above 0.6 for τ = 2−8.The overall shape of the simulated fractions mostly resembles a e −x•τ/µ curve with x ∼ 0.05.As long as severely absorbing aerosols are absent, a large fraction of the surface fluorescence source is thus conserved in the diffuse radiation field and still contributes to the TOA radiance.We therefore conclude that chlorophyll fluorescence retrievals based on the Fraunhofer line method are hardly susceptible to atmospheric aerosols.In the future, chlorophyll fluorescence remote sensing may even allow for the quantification of photosynthesis in the presence of clouds, a parameter hitherto entirely inaccessible from space and impossible to achieve using reflectance-based retrievals.
The impact of atmospheric extinction between the surface and the TOA on F s retrieval has been further analyzed with real GOSAT F s retrievals and AOD data from the AErosol RObotic NETwork (AERONET) (Holben et al., 1998).Temporal series of F s with and without the compensation of atmospheric transmittance have been produced.The aim of this exercise is to evaluate the error in monthly averages of F s due to not accounting for atmospheric extinction between the TOA and the top-of-canopy (TOC) levels.TOA F s retrievals between June 2009 and March 2011 have been generated Introduction

Conclusions References
Tables Figures

Back Close
Full with the method described in Guanter et al. (2012).AERONET level 2.0 data (cloudscreened and quality assured) are used.Since only an estimate of what the impact of realistic AOD levels can have on F s , one year of AOD values have been replicated to cover the entire 22 months of the F s series.Temporal series of GOSAT-based F s retrievals have been extracted from 4 • ×4 • boxes centered at different AERONET stations.The atmospheric transmittance associated to each of the single F s retrievals in the temporal series has been estimated with a lookup table providing total atmospheric transmittance (Rayleigh plus aerosol, both direct and diffuse radiation fields) as a function of the AOD at 550 nm from AERONET observations.A rural aerosol model has been assumed for these simulations.Monthly averages are then calculated from both the TOA and the atmospherically-corrected F s series.
Results for three different areas are shown in Fig. 13.The top panel corresponds to the Ilorin station in Central Africa ( 8• N, 4 • E), which is the site for which the highest AOD values were found.High AOD events in this area are associated to wildfires happening over the dry season.It can be stated that the impact of those high aerosol loadings in the F s temporal series is relatively small, with differences in the F s monthly averages being normally below 0.1 mW m −2 sr −1 nm −1 .There seems to be a very low seasonality in the difference between the TOA and the TOC F s , whereas there are strong seasonal signals in both F s and AOD.This is due to the fact that atmospheric extinction is a multiplicative error source on F s , which makes the seasonality of The GOSAT chlorophyll fluorescence retrievals, for instance, typically exhibit a 1−σ single sounding precision error of 0.5 % of the continuum level radiance (on the order of ±0.5 W m −2 sr −1 µm −1 at 755 nm in absolute flux terms for standard scenarios).Precision errors for OCO-2 are expected to be of similar magnitude, potentially somewhat smaller (see Frankenberg et al., 2011a).
A hypothetical instrument covering the larger spectral range from 749-760 nm, for instance, could still measure the O 2 A-band but provide a largely improved chlorophyll fluorescence retrieval, which would not only help improve the fluorescence for its own sake but also result in an improved scattering properties retrieval as the chlorophyll fluorescence term is better constrained.We repeated the sensitivity analysis in Frankenberg et al. (2011a) but with an extended spectral range and taking the spectral shape of the fluorescence into account (all errors are normalized to F s at 755).
Figure 15 shows expected instrument performance using the extended fitting range and 0.05 nm spectral resolution (full width at half maximum).At SNR = 1000, the single measurement noise error is only 0.1 % of the continuum radiance, 5 times lower than for

Conclusions References
Tables Figures

Back Close
Full a typical GOSAT sounding (with lower SNR and spectral range but higher spectral resolution).This implies that 25 GOSAT retrievals would be needed to achieve the same standard error (or about 6 for an SNR of 500).At this level of precision, single soundings will become useful individually, as even relatively small changes in fluorescence could be distinguished.Hence, a more dedicated chlorophyll fluorescence measurement purely based on Fraunhofer (not oxygen, like FLEX) line retrievals is feasible with adequate instrumentation because the precision error is not intrinsically high but can be largely reduced by a proper (and realistic) choice of wavelength ranges, spectral resolution and SNR.

Conclusions
We have quantified the impact of chlorophyll fluorescence on current multi-band X CO 2 retrievals employing the O 2 A-band.Neglect of fluorescence leads to biases in retrieved aerosol height and/or surface pressure, which propagate into X CO 2 .Typical errors of a variety of retrieval setups range from a positive 0.5-1 ppm bias if fluorescence constitutes 1 % of the continuum level radiance.Even though this bias is well below 1 %, its will be spatially correlated with gross primary production and thus must not be neglected.We show that the bias can be eliminated if fluorescence is implemented as a state vector element in the full-physics retrieval algorithm.However, interference of fluorescence with scattering state vector elements leads to a sub-optimal fit of the fluorescence signal itself, though this may be due to deficiencies in our chosen implementation.We show that a proposed fluorescence fit using Fraunhofer lines only (Frankenberg et al., 2011a) results in a much more stable and accurate fit of the chlorophyll fluorescence signal.It is further shown that chlorophyll fluorescence retrieved in this way is only moderately susceptible to the impact of atmospheric scattering, with most of the signal retained even in the presence of moderate clouds (optical depths < 5).It can be concluded that chlorophyll fluorescence not only offers the most direct proxy for gross primary production (as shown in Frankenberg et al., 2011b) but that its retrieval

Conclusions References
Tables Figures

Back Close
Full  Full  Full    (Frankenberg et al., 2011a).Points are for true F s at TOA vs. retrieved F s (black) and true F s at the surface vs. retrieved F s (red).Small deviations from the 1.0 slope can be explained by small differences in the spectral shift and squeeze as well as the assumed constancy of F s (Frankenberg et al., 2011a).Error-bars indicate a typical 1−σ precision error for a GOSAT specific retrieval.Dark gray indicates the single measurement noise and light gray the standard error in an average of 10 soundings.Introduction

Conclusions References
Tables Figures

Back Close
Full Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Fig.11).However, current instruments were not designed to retrieve fluorescence, thus the fitting ranges inFrankenberg et al. (2011a);Joiner et al. (2011);Frankenberg et al. (2011b) are chosen out of necessity, limited by the spectral coverage of the respective instrument.Relatively small changes to the GOSAT and OCO-2 instrument design, however, could largely improve retrieval precision: (I) the fluorescence signal strongly increases towards shorter wavelength, peaking at about 740 nm Guanter et al. (see 2010); (II) there are numerous isolated and strong Fraunhofer lines in the range between 740-760 nm.The GOSAT chlorophyll fluorescence retrievals, for instance, typically exhibit a 1−σ single sounding precision error of 0.5 % of the con- Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |is, unlike traditional optical parameters, much more robust and less affected by fractional cloud-cover and seasonal variations in aerosol optical depth.The only drawback of current satellites such as GOSAT or in the future OCO-2 is the high single measurement noise.However, we show that simple modifications to the instrument design, in particular the extension of the spectral range to shorter wavelengths, can largely reduce measurement noise, paving the way towards future, more dedicated fluorescence missionsDiscussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Fig. 1. 2-yr average of retrieved chlorophyll fluorescence at 757 nm, retrieved and corrected according to Frankenberg et al. (2011b).

Fig. 2 .Fig. 3 .Fig. 4 .Fig. 5 .Fig. 7 .
Fig. 2. Simulated noise-free spectrum with GOSAT spectral resolution.Solar zenith angle = 26.35• ; surface albedo ≈ 0.2; aerosol + cloud optical depth = 0.46.The fluorescence spectrum at the TOA is modeled once in a pure Rayleigh atmosphere and once including aerosol and cloud scattering of the emission emanating from the surface.Positions of Fraunhofer lines most suitable for fluorescence retrievals in the GOSAT spectral range are indicated by green lines.

Fig. 11 .
Fig.11.Retrievals of F s in a scattering atmosphere using the Fraunhofer line method(Frankenberg et al., 2011a).Points are for true F s at TOA vs. retrieved F s (black) and true F s at the surface vs. retrieved F s (red).Small deviations from the 1.0 slope can be explained by small differences in the spectral shift and squeeze as well as the assumed constancy of F s(Frankenberg et al., 2011a).Error-bars indicate a typical 1−σ precision error for a GOSAT specific retrieval.Dark gray indicates the single measurement noise and light gray the standard error in an average of 10 soundings.

Fig. 12 .
Fig. 12. Fraction of the surface fluorescence emission observed at 755 nm at the top-ofatmosphere (TOA) as a function of total optical depth.The dots indicate simulated measurements using realistic aerosols and cloud scenarios.The red line indicate theoretical curves assuming a simple exp(−x • AOD) relationship where x takes values of 1, 0.5, 0.2, 0.1, 0.05 (curves left to right).

Fig. 15 .
Fig. 15.Single measurement precision error for an imaginary spectrometer with 0.05 nm FWHM spectral resolution, varying signal-to-noise ratio (SNR) and spectra fitting windows for fluorescence.

Table 1 .
Parameters of the fluorescence spectral shape model used in simulations.

Table 2 .
Simulation and retrieval setups used in the sensitivity studies.Baseline runs are the benchmark for business as usual (fluorescence neither simulated nor retrieved), runs #xa have fluorescence in the simulation but ignore them in the retrieval and runs #xb attempt to fit fluorescence in the retrieval as well.