Next Article in Journal
A Framework for Evaluating Land Use and Land Cover Classification Using Convolutional Neural Networks
Next Article in Special Issue
Relation between Changes in Photosynthetic Rate and Changes in Canopy Level Chlorophyll Fluorescence Generated by Light Excitation of Different Led Colours in Various Background Light
Previous Article in Journal
Class-Specific Anchor Based and Context-Guided Multi-Class Object Detection in High Resolution Remote Sensing Imagery with a Convolutional Neural Network
Previous Article in Special Issue
Upscaling Solar-Induced Chlorophyll Fluorescence from an Instantaneous to Daily Scale Gives an Improved Estimation of the Gross Primary Productivity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Diurnal and Seasonal Solar Induced Chlorophyll Fluorescence and Photosynthesis in a Boreal Scots Pine Canopy

1
School of GeoSciences, University of Edinburgh, Alexander Crum Brown Road, Edinburgh EH9 3FF, Scotland, UK
2
Direction de la recherche forestière, Ministère des Forêts, de la Faune et des Parcs, 2700 rue Einstein, Québec, QC G1P 3W8, Canada
3
Optics of Photosynthesis Laboratory, Institute for Atmospheric and Earth System Research Forest Sciences, University of Helsinki, P.O. Box 27, 00014 Helsinki, Finland
4
Image Processing Laboratory (IPL), Parc Científic, Universitat de València, 46980 Paterna, València, Spain
5
Finnish Meteorological Institute, Erik Palmenin Aukio 1, P.O. Box 501, FI-00101 Helsinki, Finland
6
NASA Goddard Space Flight Center, Greenbelt, MD 20740, USA
7
NERC Field Spectroscopy Facility, School of GeoSciences, Grant Institute, West Mains Road, Edinburgh EH9 3JW, UK
8
Institute for Atmospheric and Earth System Research/Physics, Faculty of Science, University of Helsinki, Hyytiäläntie 124, FI-35500 Korkeakoski, Finland
9
Institute for Atmospheric and Earth System Research/Physics, Faculty of Science, University of Helsinki, PO Box 48, FI-00014 Helsinki, Finland
*
Author to whom correspondence should be addressed.
Remote Sens. 2019, 11(3), 273; https://doi.org/10.3390/rs11030273
Submission received: 13 December 2018 / Revised: 17 January 2019 / Accepted: 28 January 2019 / Published: 30 January 2019

Abstract

:
Solar induced chlorophyll fluorescence has been shown to be increasingly an useful proxy for the estimation of gross primary productivity (GPP), at a range of spatial scales. Here, we explore the seasonality in a continuous time series of canopy solar induced fluorescence (hereafter SiF) and its relation to canopy gross primary production (GPP), canopy light use efficiency (LUE), and direct estimates of leaf level photochemical efficiency in an evergreen canopy. SiF was calculated using infilling in two bands from the incoming and reflected radiance using a pair of Ocean Optics USB2000+ spectrometers operated in a dual field of view mode, sampling at a 30 min time step using custom written automated software, from early spring through until autumn in 2011. The optical system was mounted on a tower of 18 m height adjacent to an eddy covariance system, to observe a boreal forest ecosystem dominated by Scots pine. (Pinus sylvestris) A Walz MONITORING-PAM, multi fluorimeter system, was simultaneously mounted within the canopy adjacent to the footprint sampled by the optical system. Following correction of the SiF data for O2 and structural effects, SiF, SiF yield, LUE, the photochemicsl reflectance index (PRI), and the normalized difference vegetation index (NDVI) exhibited a seasonal pattern that followed GPP sampled by the eddy covariance system. Due to the complexities of solar azimuth and zenith angle (SZA) over the season on the SiF signal, correlations between SiF, SiF yield, GPP, and LUE were assessed on SZA <50° and under strictly clear sky conditions. Correlations found, even under these screened scenarios, resulted around ~r2 = 0.3. The diurnal responses of SiF, SiF yield, PAM estimates of effective quantum yield (ΔF/Fm′), and meteorological parameters demonstrated some agreement over the diurnal cycle. The challenges inherent in SiF retrievals in boreal evergreen ecosystems are discussed.

1. Introduction

Photosynthesis, the light-driven conversion by plants of atmospheric carbon dioxide (CO2) into carbohydrates, is one of the Earth’s largest sinks for CO2, thus playing a key role in determining the global carbon (C) balance [1,2]. The spatial and temporal dynamics of photosynthetic CO2 uptake (Gross Primary Productivity, or GPP) can provide key information about the magnitude and variability of the fraction of atmospheric CO2 that is absorbed by the terrestrial surface. Direct quantification of GPP can be carried out via micrometeorological methods, although this approach measures the Net Ecosystem Exchange (NEE), which is a combination of both GPP and ecosystem respiration, and therefore these processes need to be disentangled. NEE measurement, although being routinely made at hundreds of flux sites around the world, falls short in representing the spatial and temporal variability of these processes globally [3]. Current methods used for estimating GPP include: (1) A combination of eddy covariance, remote sensing and gridded satellite climate products, e.g., Reference [4], (2) satellite metrics of greenness and climate variables through modelling approaches [5], and (3) process-based models integrated into Earth system models [6]. However, in all these methods uncertainties propagate through to the final GPP estimates, and for this reason a need for independent estimates via direct observation still remains.
Improving our ability to monitor photosynthetic C uptake globally has become of major interest to climate and earth system modelers and to the remote sensing community, as this has direct implications for the development of climate mitigation scenarios or for determining targets in reduction of fossil fuel emissions (IPCC 2007). Currently, two main approaches exist to estimate photosynthesis directly from remotely sensed measurements. The first is based on measurements of the Photochemical Reflectance Index (PRI), a numerical index that uses the spectral reflectance values measured in a narrow detection band at 531 nm and in a reference band at 570 nm [7]. At short time scales (e.g., seconds, hours), variations in the PRI can be related to the amount of absorbed energy in the leaf that is directed to non-photochemical quenching (NPQ), a process by which plants under stress will dissipate excess absorbed energy as heat [8,9,10]. During that process, xanthophyll pigments in the chloroplast undergo reversible or sustained changes through de-epoxidation of violaxanthin into zeaxanthin, resulting in changes in leaf-level reflectance at 531 nm that can be related to canopy-level photosynthetic down-regulation. At longer time scales (e.g., days, months), variations in the PRI become increasingly related to changes in carotenoids and chlorophyll pigment pools [11,12], which have also been related to canopy photosynthetic light-use efficiency (LUE) [13,14]. The second approach for estimating photosynthesis from remotely sensed measurements is based on solar-induced chlorophyll fluorescence (SiF): The weak and near-instantaneous emission by chlorophyll molecules of redshifted light (i.e., lower energy relative to absorbed light) during photosynthesis [15]. In actively photosynthesizing plants, the amount SiF changes continuously in response to variations in environmental conditions.
While SiF is very small compared to the amount of energy that is reflected, it is possible to indirectly measure SiF through the application of appropriate algorithms such as those based on Fraunhofer Line Discriminator, Spectral fitting, or indirectly by using spectral indices based on the reflectance measured at SiF-related wavelengths [16,17,18,19,20]. Remote sensing of SiF is complex and challenging; apart from fluctuating in response to environmental conditions, the small amount of SiF also depends on the leaf or canopy chlorophyll content and on the fraction of energy that is directed towards NPQ relative to SiF and photochemistry [15,21,22]. As NPQ and SiF are both closely linked to photosynthetic efficiency, an increase/decrease in the amount of energy directed to any of them will lead to a decrease/increase in the energy available for the remaining two photosynthetic pathways. Therefore, if we can accurately estimate SiF and NPQ, it is then possible to have a measure of the level of efficiency of carbon uptake by plants. While both the PRI and SiF-based approaches were successfully applied at the plant and small canopy levels [23,24,25], their application to larger spatial scales using airborne and satellite data remains challenging [26,27,28,29,30,31,32,33,34].
However, recent articles have highlighted, with varying degrees of success, the ability to extract SiF from orbiting platforms (e.g., from the Japanese Greenhouse gases Observing SATellite, GOSAT, and from the Global Ozone Monitoring Experiment-2, GOME-2) for estimating GPP [27,35,36,37]. These findings highlighted that spatial and temporal patterns of satellite retrieved SiF were highly correlated with GPP at the biome and global levels [33], but such correlations weaken at individual higher temporal and spatial resolutions [38] or during extreme climatic events [39]. Field based studies are continuing to emerge [20,25,29,40,41,42,43,44] though long term data sets, covering multiple seasons, are still rare [45,46,47]. Additionally, many of these studies have been restricted to crop or shrub canopies (e.g., References [47,48,49,50]) though one recent study has been in a deciduous forest [25].
In this study, we obtained a high frequency time series of canopy SiF, from early spring, during physiological inactivity, to peak physiological activity in summer, and decrease into autumn, over an evergreen boreal scots pine canopy, alongside measurements of PAM fluorescence, GPP from eddy covariance and environmental variables. We aimed to address the following questions: (1) How does SiF change over the growing season in a boreal evergreen ecosystem experiencing modest changes in APAR? (2) Is the SiF-GPP relationship as good as previously reported for crops and deciduous forests? (3) Is SiF able to capture changes in canopy physiological dynamics (LUE)?

2. Materials and Methods

2.1. Site Description

The study site was located at the Station for Measuring Forest Ecosystem-Atmosphere Relations (SMEAR II) in Hyytiälä, Southern Finland (61.847 N, 24.294 E, elevation 181 m), which has been measuring ecosystem fluxes (e.g., CO2, H2O), climatic and ecological variables continuously since 1995 [51]. The studied canopy was a Scots pine (Pinus sylvestris L.) stand established in 1962, with stem density, mean tree height and mean diameter at breast height (1.3 m) of 755 stems ha−1, 16 m and 18 cm, respectively [52]. Four years prior to our study, all-sided leaf area index (LAI) of pine trees around the eddy covariance tower was estimated at 6.5 [53], which corresponds to a projected LAI of approximately 3. Ground vegetation around SMEAR II is composed mostly of the ericaceous shrubs Calluna vulgaris L., Vaccinium vitis-idaea L. and V. myrtillus L., and of the moss species Dicranum undulatum. The terrain at the site is relatively flat, with a thin soil (50–150 cm) classified as a Typic Haplocryod [51,54]. The mean annual temperature is 4.3°C and the mean annual precipitation is 590 mm [52].

2.2. Spectral Data

The dual-field-of-view optical system which used a pair of Ocean Optics USB2000+ spectrometers operated in dual field of view mode, sampling at a 30 min time step using custom written automated software, and fully described in Reference [55] was installed on the SMEAR II eddy covariance (EC) flux tower during in early April 2011. It operated continuously between April 23 and September 30 (days of year 113 to 273), acquiring simultaneous incoming solar irradiance and canopy-reflected radiance every 15 min.
The canopy-pointing optical fiber, with a FOV of 24.8° (yielding an off nadir instantaneous FOV of approximately 400 m2 on the canopy) was attached to the tower at 7 m above the tree tops, with azimuth and viewing angles of 280° and 70° (relative to nadir), respectively. This configuration allowed us to measure about 400 m2 of mostly sunlit pine crowns when the sun was at 35° degrees or more above the horizon (see Figure 5 in Reference [55]) minimizing the impact of ground cover dynamics on APAR (e.g., snow vs green ground vegetation) and facilitating the investigation of the impact of canopy LUE dynamics on SiF. The system cosine-corrected optical fiber (CC-3, Ocean Optics Inc., Dunedin, FL, USA) was attached to the tip of a vertical pole at the top of the tower, sitting at about 1 m above the canopy-pointing fiber. This resulted in an almost completely unobstructed view of the sky; only a small part of the nearby SMEAR II tall tower was included in the fiber IFOV.
Intensity counts recorded by the spectrometers were converted to radiation fluxes, namely downwelling irradiance and upwelling radiance, using the procedure described in Reference [55]. Calibration coefficients used to calculate fluxes were obtained from an on-site calibration of the optical system prior to the start of the study period. Calculated radiation fluxes were also used to derive a spectral reflectance spectrum for each periodic acquisition by the optical system.
Diurnal variations in cloud cover contribute to making physiological interpretation of SiF and spectral reflectance data more difficult. Apart from the obvious direct effects of clouds decreasing the overall PPFD and subsequently SiF signal, clouds also increase the ratio of diffuse-to-direct radiation, which will change the radiation field within the canopy and the relative contribution of different canopy components to the measured SiF signal. To minimize potential cloud effects on the spectral data, we visually inspected time series of 30 min photosynthetically active photon flux density (PPFD), which was measured for the whole sampling period by a quantum sensor (Li-190SZ, LI-COR Inc., Lincoln, NE, USA) on the nearby SMEAR II radiation tower. We split the dataset into clear and cloudy observations through a visual inspection of the shape of PPFD diurnal patterns: We considered as clear spectral data points only those ones that occurred on a Gaussian response in PPFD over the diurnal course [56]. Finally, to minimize the impact of structural effects on spectral data, we calculated solar zenith angle (SZA) values for each spectral acquisition using the site geographic coordinates and date-time of measurement; these were used as filtering conditions in regression analyses.

2.3. SiF and Spectral Indexes

A number of methods are reported for the extraction of SiF using the FLD method and variants on this approach [57] The approach taken in this paper utilised two band FLD approaches due to the relatively coarse resolution of the spectral data (1 nm resolution), obtaining the exact position of the two shoulders of the absorption feature accurately was not possible, and thus the more popular 3FLD approach was not used. SiF values were extracted from the spectral data using [16] applied to the O2-A absorption region of the electromagnetic spectrum (~755–765 nm). Damm [57] reported that all FLD approaches were similarly affected when coarser resolution data was utilised, and the two wavelength approaches tended to overestimate SiF. Accordingly, because the absolute SiF values reported in the present study may not be free from overestimation, we focus the analysis on the patterns of variation rather than absolute levels.
Wavelength matching between pairs of temporally concomitant spectra, a requirement of the FLD method, was achieved by interpolating down-welling irradiance spectra to the same wavelengths as those of reflected radiance spectra using locally-weighted polynomial regression. Then, for each pair of measurements, wavelengths of local minima between 758 nm and 762 nm were identified in both irradiance and radiance spectra and were used as absorption wavelengths for retrieving energy fluxes. FLD shoulder wavelengths were identified by visual inspection of irradiance and radiance spectra; a constant value at 757.5 nm was used for both spectra and all retrievals.
To account for energy re-absorption by molecular oxygen along the photon path between the sensor and the canopy, we applied the correction proposed in Reference [58]. Atmospheric transmittance values in the O2-A region were calculated for each SiF observation using the HITRAN molecular spectroscopic database [59,60] along with air temperature and atmospheric pressure values measured at the flux tower. All data presented in this paper have therefore been corrected for atmospheric effects.
Last but not least, because the amount of energy re-emitted by plants as fluorescence is directly proportional to the amount of incoming radiation, the functional interpretation of the SiF signal (either in terms of Fluorescence Yield or in terms of photosynthetic energy absorption or APAR) is not straightforward. Signal normalization by incoming PAR or indeed APAR is commonly a challenge because downwelling PAR radiation registered over the canopy does not necessarily reflect the actual PAR experienced by the foliage within the upwelling sensor FOV, due for example to diurnal changes in sun/shade fractions. In an attempt to control for variation in PAR directly within the sensors IFOV, we calculated a relative SiF yield ( S i F   L 650 1 ), which is adimensional, by normalizing SiF observations by the reflected radiance measured in the red region of the spectrum (650 nm). In other words, diurnal and seasonal variations in L650 were used as proxy of average incoming PAR within the scene. The usage of reflected radiance (Lλ) as a proxy of variations in incoming PAR within the IFOV is based on two assumptions (See Equation (1) below): (i) Constancy in canopy reflectance in wavelength λ (ρλ), and (ii) proportionality between irradiance in the PAR region (EPAR) compared to that at wavelength λ (Eλ).
L λ = E λ   ρ λ =   k   E P A R   ρ λ
We preliminarily selected 650nm band, but other wavelengths within PAR spectral domain could be equally useful. Given the common use however of a normalisation of SiF using APAR, we further computed SiF/APAR for purposes of comparison with the above normalisation scheme used in this study. APAR is described and calculated according to Equation (6).

2.4. Leaf-Level Chlorophyll Fluorescence

To understand seasonal changes in SiF values derived from the optical system, leaf-level fluorescence measurements were carried out on needles from three Scots pine trees located near the IFOV of the optical system. These measurements were made using a system of field fluorometers (MONITORING-PAM, Heinz Walz GmbH, Effeltrich, Germany) [61,62]. Fluorometers were installed in the top canopy and South-facing shoots across three different trees. Every 15–30 min, a saturating pulse was applied to each sampled branch and actual and maximal fluorescence yield from light-exposed leaves (Ft and F′m, respectively), PPFD at sample level (PPFDSL, in μmol m−2 s−1) and leaf temperature (°C) were recorded (see Reference [61] for more details). Quantum yield (i.e., operating efficiency) of photosystem II (PSII) photochemistry ( Δ F / F m ), an indicator of the acclimation of photosynthetic reactions to environmental conditions, was calculated for each measurement as:
Δ F F m = Φ PSII = ( F m F t ) F m
Values of Δ F / F m from all three fluorometers were averaged together for each measurement time. To match averaged Δ F / F m values with the times of the SiF observations, we applied a simple weighted average interpolation method:
Φ PSII t = Φ PSII t 1   w t 1 +   Φ PSII t + 1   w t + 1
where Φ PSII t is the quantum yield of photochemistry at time t of a given SiF observation, Φ PSII t 1 and Φ PSII t + 1 are Φ PSII values at the closest point in time before and after time t, respectively, and w t 1 and w t + 1 are the normalized distances from time t and calculated as:
w t 1 = 1 ( t t t 1 ) ( t t + 1 t t 1 )   ;   w t + 1 = 1 ( t t + 1 t ) ( t t + 1 t t 1 )
where t − 1 and t + 1 are the closest observations in time before and after time t, respectively. Finally, time series of Φ PSII t were visually inspected and outliers removed. Φ PSII t will be referred to as Δ F / F m in the remainder of the paper.

2.5. Gross Primary Production and Light-Use Efficiency

Time series of gross primary production (GPP, in μmol m−2 s−1) and light-use efficiency (LUE) were derived from CO2 net ecosystem exchange (NEE) measured using the eddy covariance technique [9]. NEE was calculated on a half-hourly basis, according to [63] from high-frequency measurements of CO2 concentration and vertical wind speed obtained from a closed-path infrared gas analyser (Li-6262, LI-COR, Inc., Lincoln, NE, USA) and a three dimensional ultrasonic anemometer (Solent 1012R, Gill Instruments Ltd., Lymington, UK), respectively. Partitioning into gross ecosystem productivity (GEP) and ecosystem respiration was done according to Reference [53]. In this study, we assumed tower GEP to be equivalent to GPP and only the later term is used in the paper. Half-hourly LUE values (LUE30) were calculated as:
L U E 30 =   G P P 30 P A R 30
where G P P 30 and PAR30 are the 30-min mean values of GPP and photosynthetically active radiation (PAR). We further calculated that absorbed by the canopy (APAR), respectively. APAR30 was calculated as:
A P A R 30 =   P A R 30 × f P A R 30
Reflected shortwave radiation measurements from a pyranometer (TP-3, Reeman) mounted on the radiation tower were used as a proxy for P A R 30 in f P A R 30 calculations:
f P A R 30 = 1 S W r 30   /SWi30, where S W r 30 are 30-min means of reflected shortwave radiation and SWi30 are 30 min means of incoming shortwave radiance matched to the time of P A R 30 observations.

2.6. Environmental Variables

Environmental conditions during the measurement period were assessed using time series of meteorological variables from sensors located on and around the SMEAR II towers: air temperature (°C) at 17 m (PT100 RTD, Omega Engineering, Inc., Norwalk, CT, USA) and precipitation (F12P, Vaisala, Helsinki, Finland). Measurements from NDVI and PRI sensors (SKR 1800 2-channel sensor, Skye Instruments Ltd., Llandrindod Wells, UK) mounted at 31 m on the SMEAR II tall mast were used to assess the seasonal dynamics of the canopy structure and function, as well as in comparisons with SiF, S i F   L 650 1 , G P P 30 , and L U E 30 .

2.7. Statistical Analyses

Variables calculated from spectral measurements were matched with corresponding 30 min averages of fluxes, meteorological and leaf-level fluorescence variables. To achieve this, only spectral observations for the 15th and 45th minute of each hour were retained since they both fall in the middle of an averaging period for flux and meteorological variables (i.e., from start to the 30th minute, and from the 30th to the 60th minute of each hour, respectively). The resulting dataset was used in regression analyses presented in the next section.

3. Results

3.1. Seasonal Patterns of Environmental and Physiological Variables

Daily photosynthetic photon flux density (PPFD), precipitation (including snow), and air temperature throughout the sampling period (DOY 113–273) are shown in Figure 1. From day 113 through to 130, air temperatures remained close to zero degrees, then started to increase until they reach a summer maximum of 29°C around day 160 (June 9), and then decreased towards the autumn. Maximum PPFD values during the summer period were around 1500–1600 µmol m−2 s−1. A week of almost cloud-free days stretched from June 5 (day 156) through to June 11 (days 156–162), coinciding with the warmest days of the summer. Precipitations as rain occurred during each month of the sampling period, but were more abundant in September (164 mm, 19 days with rain ≥ 1 mm) than in any other months of the sampling period, for which monthly precipitations varied between 22 mm (April) and 96 mm (August), and the number of days with rain accumulation of at least 1 mm was between 6 (April) and 11 days (May and June).
Figure 2 shows samples of irradiance and reflected radiance spectra measured by the optical system on a clear day of the sampling period. The absorption feature related to the oxygen-A absorption band around 760 nm can be clearly appreciated in both curves (shaded area). In this region, the amount of energy reflected by the dark Scots pine canopy represents about 20% of the incoming energy, whereas in the visible region of the spectrum (400–700 nm), this reduces to between 3% and 5% only.
The seasonal course of SiF observations (day and night), including all ranges of solar zenith angles, diurnal variation, and without accounting for structural effects, light levels or cloud cover, exhibited a weak and noisy seasonal pattern throughout the sampling period (Figure 3A). After normalization by our proxy of incoming PAR in the IFOV (L650), S i F   L 650 1 presented a more consistent seasonal pattern, with a much smaller range of variation, and that more closely resembled that of G P P 30 . The quantum yield of photochemistry measured by the PAM system also showed a typical seasonal pattern similar to those of G P P 30 and S i F   L 650 1 , (Figure 3D). Note that the top envelope of points in Figure 3D represents night data and therefore corresponds to the widely used parameter Fv/Fm. Similarly, the lower envelope of points represents noon data points. An instrument malfunction resulted in a three week data gap in the Δ F / F m time series between July 26 and August 19 (days of year 207 to 231). Skye sensor PRI and NDVI both showed a slight incline through the measurement period into summer, although the seasonal pattern was flatter for PRI than NDVI. Skye sensor PRI and NDVI both showed a slight incline through the measurement period into summer, although the seasonal pattern was flatter for PRI than NDVI. Importantly, the rapid increase in PRI towards day 130 (mid May) took place concomitantly with an increase in the apparent quantum yield of PSII (Figure 4G) and slight increase in LUE (Figure 4F), whereas S i F   L 650 1 remained stable during that period.
Seasonal patterns similar to those shown in Figure 3A–F, but with much less noise, were obtained by performing a midday average and removing observations from rain and snow days (Figure 4). Observations acquired during these periods are characterized by lower solar irradiance levels and thus, contribute to reducing the signal-to-noise ratio of the spectral data. They are also increasing noise levels in GPP estimates due to the small differences between NEE and ecosystem respiration fluxes around dawn and dusk. While accounting for low irradiance observations removed noise in the SiF time series, it did not improve its overall seasonality (Figure 4B). The seasonality in S i F   L 650 1 improved (Figure 4C) when the SiF signal was normalized by the radiance at 650 nm, but not when compared to the SiF normalized to APAR.

3.2. Relationships between SiF, SiF L650−1, GPP30 and LUE30

To explore the relationships between SiF, SiF L650−1 with GPP30 and LUE30, three steps were undertaken. Firstly all data were included in the analysis; these also included cloudy and clear observations as well as all solar angles (and shown in grey circles in Figure 5A–H). In a second step, data points were filtered to remove all observations that were not acquired on completely clear days (dark circles in Figure 5A,C,E,G). Finally, data points were filtered to also remove SZA values larger than 50°, thus reducing spectral data variations due to changing illuminated/shadowed canopy fractions (Figure 5B,D,F,H). The relationships between SiF, SiF L650−1 with GPP30 and LUE30, when considering data collected at all SZAs showed scattered and non-significant relationships (Figure 5A,C,E,G). There was marginal improvement in the relationships when data were filtered to include only higher SZAs (>50°), as shown in Figure 5B,D,F,H. Linear regression analyses applied to these data showed only weak correlations between the analyzed variables (R2 0.24–29, p < 0.0001) (Table 1).
Interestingly, the relationship between SiF and GPP was better than between SiF L650−1 and GPP, which could be expected since both SIF and GPP are strongly controlled by PAR, whereas SiF L650−1 is readily normalized by PAR via SiF L650−1. In contrast, LUE was slightly better related to SiF L650−1 than to SiF, which was also expected since both LUE and SiF L650−1 undergo normalization by PAR.
We selected one clear day from each month from April through September and explored the relationships between SiF, GPP, PPFD, and SZA over a diurnal time course (Figure 6). GPP typically followed PPFD throughout the day in all days across the months, although its magnitude was smaller in the spring (April and May) than later during the growing season. During the early hours following sunrise, SiF increased following GPP and PPFD and peeked around 6 UTC. After this point, SiF decreased throughout the day, while GPP remained relatively stable or continued to increase until it peaked around solar noon. This pattern in SiF was consistent across clear days, but showed more variations on days when PPFD was less stable (July and September). Maximum SiF values occurred in May and July. Note that the time in Figure 6 is UTC, with solar noon occurring at 10AM.
Using the same clear days used in Figure 6, the diurnal relationships between leaf-level effective quantum yield from the PAM data, SiF L650−1 and SZA were also explored (Figure 7). During clear summer days, SiF L650−1 consistently showed a small peak around 4–5 UTC (Figure 7B–E). It then remained constant until around noon, when it started to increase until reaching a second peak of the same magnitude as that observed in the morning and which occurred from late- to mid-afternoon as the summer progressed. Past that afternoon peak, SiF L650−1 usually dropped rapidly toward the evening. PAM ΔF/Fm′ showed diurnal courses going in opposite direction to those of SiF L650−1. Most months showed a maximum ΔF/Fm′ in the early hours of the morning, which decreased until it reached a minimum value at a time that occurred increasingly later as the season progressed. Generally, ΔF/Fm′ and SiF L650−1 both decreased from the morning until midday to early afternoon, after which they either moved in the same direction (Figure 7B,C) or diverged (Figure 7D,E). However, the time at which each variables reached a minimum daily value and started to increase and move in opposite directions as the season progressed. For SiF L650−1, this time occurred earlier during the day from April to September, while for ΔF/Fm′ it occurred later as the months passed.

4. Discussion

In this study, we deployed a tower-based optical system [55] to collect a continuous time series of canopy SiF over an evergreen canopy alongside measurements of PAM fluorescence, GPP from eddy covariance and environmental variables, NDVI and PRI, with the aim of describing both the diurnal and seasonal patterns of SiF in connection with ecosystem GPP dynamics. We aimed at exploring the following questions: (1) How does SiF change over the growing season in a boreal evergreen ecosystem experiencing modest changes in APAR? (2) Is the SiF-GPP relationship as good as previously reported for crops and deciduous forests? (3) Is SiF able to capture changes in canopy physiological dynamics (LUE)? To the best of our knowledge, this is the first long-term time-series of SiF in an evergreen canopy from field-based instruments.
While a strong SiF-GPP relationship has been repeatedly observed at the ecosystem level using ground based instruments in crop and deciduous canopy types, [18,23,24,25,29,40,41,42,45,47,50,64,65,66,67], as well as at landscape/regional level from satellite SiF retrievals [26,27,31,33,34,68,69,70,71,72,73], the factors that drive these relationships remain controversial. Biome dependent differences in canopy structure have been suggested to affect the slope of the relationship [35,64,68,74], but recent work at the satellite pixel scale is suggestive of a universal relationship [31]. Similarly, the relative role that ecosystem APAR and LUE dynamics exert on linking SiF to GPP is also expected to be biome dependent, but its characterization is still limited. Observations across crops and deciduous forests suggest that APAR dynamics could be the main factor connecting SiF to GPP [25,73,75]. This is logical, since annual GPP dynamics for these vegetation types are dominated by APAR, which in turn controls both the SIF emission and photosynthetic CO2 assimilation. The question remains as to how well SIF and GPP are related for an evergreen ecosystem, where only modest changes in APAR occur.

4.1. Comparison between Evergreen Forest SiF and Other Terrestrial Ecosystems

The seasonal midday signal of SiF, without any structural correction, was generally very noisy and did not exhibit any seasonal pattern (Figure 4B), this is in contrast to previous studies (albeit over crops and deciduous forests) that demonstrated a strongly seasonal pattern in SiF, [18,23,24,25,29,40,41,42,45,47,50,64,65,66,67]. We found only a very modest seasonal increase in NDVI (Figure 4E), as recorded by an independent broadband sensor installed 30 m above the canopy, likely indicative of minor changes in ecosystem APAR. In addition, because our upwelling radiation sensor was purposefully pointing towards the tree crowns, we would expect that APAR changes in the IFOV would be even smaller due to absent contribution from the ground vegetation. Taken together, these results indicate that the variation in absolute SiF of our evergreen pine canopy was dominated by incoming PAR. Normalization of SiF by incoming PAR has previously been carried out using PAR data acquired above the canopy [42], but unfortunately because of the dynamic sun/shade patterns within the IFOV during the course of the day and passing of the seasons, the actual PAR received by the foliage under examination will differ from that recorded above the canopy following complex temporal patterns and partly undermining the normalization. In an attempt to overcome this limitation, we used reflected radiance at 650 nm (L650, see Equation (1)) as a proxy of actual PAR received by the foliage within the IFOV. After normalization by L650, a clear seasonal pattern in SiF L650−1 emerged (Figure 4C), which more closely followed the time course of GPP (Figure 4A). The SiF L650−1 and indeed GPP, increased through the spring transition towards summer then declined into the autumn period. It is important to note that L650 normalizes the signal by PAR, but not APAR, and therefore, increase in SiF L650−1 could still be due to both changes in APAR (e.g., new needle cohort and shoot elongation during June, and subsequent senescence of old needle cohorts during September), as well as leaf-level adjustments if fluorescence quantum yield (related to LUE). It has been demonstrated that normalisation (using various means such as PAR and APAR) can improve the SiF-GPP relationship in some cases [76], but it is recognised that it is not always possible due to unavailability of measurements like APAR. Certainly attention to various normalisation schemes is warranted in order to understand and separate the physiological and non-physiological components.
Our correlation analysis of SiF, SiF L650−1, and GPP30, LUE30 resulted in only weak relationships (Figure 5), which were again generally weaker than the SiF-GPP relationships reported in earlier seasonal and canopy-level studies for crops and deciduous forests or at larger scales using satellite data [35,43,47,50,73,77,78]. This demonstrated the complexity of extracting the physiological signal component from tower based SiF studies in complex forest ecosystems. To date few studies have focused on evergreen ecosystems. Walther [79], using SiF retrieved from the GOME-2 instrument, explored the seasonality in SiF in a high latitude evergreen ecosystem, and reported a strong seasonal relationship between satellite SiF and modelled GPP. Furthermore, Wolfhart [39] studied a short time series (10 days in length) of SiF in an evergreen canopy that had experienced a short but intense heatwave, and found that during a short period of unchanging APAR, SiF was only very weakly related to the change in GPP, a change attributed to the fact APAR was almost unchanged during the analysis period. In our study, statistically significant relationships only emerged following filtering to limit analyses to midday, cloud free days with high SZA. The correction of SiF to SiF L650−1 improved the strength of the relationship with LUE but not with GPP (Table 1). This is logical since both SiF and GPP are strongly controlled by incoming PAR (indeed APAR), in contrast, SiF L650−1 is intrinsically normalized by PAR via L650−1. In fact, because SiF L650−1 is normalized by PAR, it becomes a relative measure of fluorescence yield and therefore should be expected to reflect the seasonal variations in photosynthetic LUE much better than SiF. Similarly long-term SiF (and SiF yield) measurements in a deciduous forest site at Harvard Forest [25] highlighted stronger relationships between SiF, SiF Yield, GPP and LUE, with similarly strong results between SiF and GPP reported in a mixed forest site [64]. Interestingly, the slight increase in apparent quantum yields of photochemistry (Figure 4G) and LUE (Figure 4H) observed here during early May, was registered by the PRI (Figure 4F) but not by SiF L650−1. Although this observation will require further validation and studies with even broader temporal coverage it preliminary points to potential limitations of SiF to track LUE in evergreens.
The nature of the relationship between SiF and GPP, continues to be questioned. Is it solely driven by the dependence on APAR or a combination of APAR and the photosynthetic light and dark-reactions? [66,78,80,81]. Long term and in situ measurement will prove particularly attractive in answering this question. Previous studies have highlighted that heat dissipation (NPQ) is the main driver of variations in fluorescence and photosystem yields, and the correlation between SiF yield and LUE is consistent with field studies and model-based assessments [25,82]. In the present study however, we found only a very weak correlation between SiF and GPP, and SiF L650−1 and LUE for the canopy component of an evergreen Scots pine forest.

4.2. Diurnal Relationships and Sun-Sensor Geometry

The diurnal patterns presented in this study (Figure 6 and Figure 7) are likely to be controlled by both physiological and optical (directional) factors. Although the tower mounted optical system used in this study was oriented in a so called “hot spot” region of the canopy, when the solar angle varies through the day, so too will the fractions of sunlit and shaded foliage in the sensors field of view. The influence of directionality will therefore inevitably generate differing proportion of SiF coming from the scene viewed by the spectrometer. Understanding and decoupling the influence of both structure and physiological response from SiF retrievals is a much needed next step, and one that models have thus far been used to understand [82,83,84]. The passage of clouds may add orders of magnitude difference to the retrieved SiF (Figure 6D) via a significant change in the measured irradiance and subsequent depth of the O2-A band [26]. In the results presented here the peaks of SiF and PPFD do not match, which is likely due to be dominated by the changes in solar geometry. Indeed the hourly fraction of sunlit and shaded leaves within the sensors FOV will be varying, as will the presence of multiple scattering within the canopy. Future work should undoubtedly focus on such shadow fraction changes, along with the impacts on SiF, in this ecosystem.
The influence of sun-sensor geometry on solar induced fluorescence itself is receiving increasing attention [85,86] with evidence from satellite retrieval analysis and indeed modelling studies mounting. The have reported noticeable angular influences each generating differing ranges of SiF values, which is similar to the effect of sun-sensor geometry on canopy reflectance data [65,84]. A bowl shaped response has been reported from backward to forward scattering directions consistently in ground and model based observation [65,84] (though not evergreen canopy specific). While it is not possible to disentangle the impact of sun-sensor geometry explicitly here, its influence remains highly probable.

4.3. Atmospheric Influence and Instrument Resolution

Retrievals of SiF in the oxygen absorption lines, from aircraft platforms hundreds of meters above the surface, have routinely been atmospherically corrected [48,64,87,88]. However, there has been little consensus until recently as to whether this correction would be needed for SiF in a proximal sensing context, i.e., from a flux tower mounted system within 20 m from the canopy. Historically SiF retrievals were formulated to be applied to top of canopy data, with (1) the assumptions of atmospheric path length between target and sensor is short enough to be neglected, and (2) that the solar irradiance is measured at the same height as the target. When atmospheric path length increases, for example with an instrument mounted on a tower or where large solar zenith and/or view angles are common, these assumptions cannot be met. Oxygen absorption is proportional to air pressure [89] and thus at this lowest level of the atmosphere even a few meters difference can result in significant error in retrieved SiF. In this study, oxygen transmittance was computed by making use of the HITRAN molecular spectroscopic database, therefore modelling the oxygen transmittance according to the experimental site configuration, e.g., fixing the optical path between the canopy and the sensor, and taking into account variations in the oxygen transmittance caused by changes in the environmental conditions, mainly temperature and pressure. Accounting for variations in the environmental conditions as part of the SiF retrieval strategy guarantees, especially in experimental sites subjected to abrupt changes of temperature and pressure, an accurate observation of the SiF seasonal patterns, as pointed out in Reference [58].
The complexity of the influences of instrument spectral resolution, signal-to-noise, atmospheric correction, canopy structure, leaf biochemical parameters and directional effects all play a critical role in shaping the reliability of SiF to quantify GPP [58,90,91]. Only recently have intercomparisons in SiF retrievals been carried out across optical instruments. Julitta [90] carried out tests on four spectrometers and investigated their ability to retrieve both red and far-red SiF. The work presented highlighted that an “accurate” far-red SiF could be retrieved from spectrometers with an ultra-fine resolution (less than 1 nm) with the red SiF estimation requiring a significantly higher resolution (less than 0.5 nm). The Ocean Optics (USB-2000) system used in this study had a 1 nm resolution, and would therefore fall at the upper end of the recommended resolution and thus may have hampered the SiF retrieval.

5. Conclusions

By linking high temporal resolution measurements of solar induced fluorescence with gross primary productivity, canopy light use efficiency and pulse amplitude measure of efficiency, we were able to, for the first time, explore the complex nature of understanding both SiF and its use in understanding evergreen canopy physiological properties. The results of study indicated that SiF alone did not prove useful over the seasonal cycle, and only after a correction term was applied to account for variation in illumination conditions, did the findings elucidate a seasonal cycle that mirrored GPP. While only a simple two-band FLD approach was utilized here, credible seasonal and diurnal trends were not dissimilar to published studies, though we highlighted the challenges in SiF retrieval in an evergreen canopy with a low LAI and low solar angles. The diurnal responses when compared to conventional PAM fluorescence highlighted structural and solar influence, with the nature of the diurnal relationship varying throughout the season. Further studies should invariably focus on higher resolution optical data and a longer time series. Only now are published works becoming available that highlight the very minimum in instrument setup and response, but no review of instruments, resolution, and methodology exists at the time of writing.

Author Contributions

C.J.N. was Principal Investigator for the project and carried out all analysis of the presented data. G.D. collected the Ocean Optics data and carried out its post processing. T.W., with input from C.M., built the Ocean Optics optical system and performed system checks before its use in this study. A.P.-C. and J.M.A. collected the PAM data and oversaw its post processing. A.P.-C. further contributed to the discussion of the results and instigated the use of the reflected radiance normalization idea. J.L. oversaw the collection of the eddy covariance data and I.M. carried out its post processing. The manuscript was written by C.J.N. with significant written contributions from G.D., A.P-C., T.W., N.S., E.M., C.M., J.L., I.M., T.V., J.A.

Funding

This work was funded by the UK Natural Environmental Research Council NERC (NE/F01749), Finnish Academy (#293443, #288039), Academy of Finland Center of Excellence (projects No. 272041 and 118780) and ICOS-Finland (project No. 281255).

Acknowledgments

We thank Sergio Cogliati and Petya Campbell for their helpful discussions of the results presented in this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Valentini, R.; Matteucci, G.; Dolman, A.J.; Schulze, E.D.; Rebmann, C.; Moors, E.J.; Granier, A.; Gross, P.; Jensen, N.O.; Pilegaard, K.; et al. Respiration as the main determinant of carbon balance in European forests. Nature 2000, 404, 861–865. [Google Scholar] [CrossRef] [PubMed]
  2. Beer, C.; Reichstein, M.; Tomelleri, E.; Ciais, P.; Jung, M.; Carvalhais, N.; Rodenbeck, C.; Arain, M.A.; Baldocchi, D.; Bonan, G.B.; et al. Terrestrial Gross Carbon Dioxide Uptake: Global Distribution and Covariation with Climate. Science 2010, 329, 834–838. [Google Scholar] [CrossRef] [PubMed]
  3. Chu, H.; Baldocchi, D.D.; John, R.; Wolf, S.; Reichstein, M. Fluxes all of the time? A primer on the temporal representativeness of FLUXNET. J. Geophys. Res. Biogeosci. 2017, 122, 289–307. [Google Scholar] [CrossRef]
  4. Jung, M.; Reichstein, M.; Margolis, H.A.; Cescatti, A.; Richardson, A.D.; Arain, M.A.; Arneth, A.; Bernhofer, C.; Bonal, D.; Chen, J.Q.; et al. Global patterns of land-atmosphere fluxes of carbon dioxide, latent heat, and sensible heat derived from eddy covariance, satellite, and meteorological observations. J. Geophys. Res. Biogeosci. 2011, 116. [Google Scholar] [CrossRef] [Green Version]
  5. Ryu, Y.; Baldocchi, D.D.; Kobayashi, H.; van Ingen, C.; Li, J.; Black, T.A.; Beringer, J.; van Gorsel, E.; Knohl, A.; Law, B.E.; et al. Integration of MODIS land and atmosphere products with a coupled-process model to estimate gross primary productivity and evapotranspiration from 1 km to global scales. Glob. Biogeochem. Cycles 2011, 25. [Google Scholar] [CrossRef] [Green Version]
  6. Schaefer, K.; Schwalm, C.R.; Williams, C.; Arain, M.A.; Barr, A.; Chen, J.M.; Davis, K.J.; Dimitrov, D.; Hilton, T.W.; Hollinger, D.Y.; et al. A model-data comparison of gross primary productivity: Results from the North American Carbon Program site synthesis. J. Geophys. Res. Biogeosci. 2012, 117. [Google Scholar] [CrossRef] [Green Version]
  7. Gamon, J.A.; Penuelas, J.; Field, C.B. A Narrow-Waveband Spectral Index That Tracks Diurnal Changes in Photosynthetic Efficiency. Remote Sens. Environ. 1992, 41, 35–44. [Google Scholar] [CrossRef]
  8. DemmigAdams, B.; Adams, W.W.; Barker, D.H.; Logan, B.A.; Bowling, D.R.; Verhoeven, A.S. Using chlorophyll fluorescence to assess the fraction of absorbed light allocated to thermal dissipation of excess excitation. Physiol. Plant. 1996, 98, 253–264. [Google Scholar] [CrossRef]
  9. Demmig-Adams, B. Survey of thermal energy dissipation and pigment composition in sun and shade leaves. Plant Cell Physiol. 1998, 39, 474–482. [Google Scholar] [CrossRef]
  10. Atherton, J.; Nichol, C.J.; Porcar-Castell, A. Using spectral chlorophyll fluorescence and the photochemical reflectance index to predict physiological dynamics. Remote Sens. Environ. 2016, 176, 17–30. [Google Scholar] [CrossRef]
  11. Gamon, J.A.; Berry, J.A. Facultative and constitutive pigment effects on the Photochemical Reflectance Index (PRI) in sun and shade conifer needles. Isr. J. Plant Sci. 2012, 60, 85–95. [Google Scholar] [CrossRef]
  12. Porcar-Castell, A.; Garcia-Plazaola, J.I.; Nichol, C.J.; Kolari, P.; Olascoaga, B.; Kuusinen, N.; Fernandez-Marin, B.; Pulkkinen, M.; Juurola, E.; Nikinmaa, E. Physiology of the seasonal relationship between the photochemical reflectance index and photosynthetic light use efficiency. Oecologia 2012, 170, 313–323. [Google Scholar] [CrossRef] [PubMed]
  13. Sims, D.A.; Gamon, J.A. Relationships between leaf pigment content and spectral reflectance across a wide range of species, leaf structures and developmental stages. Remote Sens. Environ. 2002, 81, 337–354. [Google Scholar] [CrossRef]
  14. Wong, C.Y.S.; Gamon, J.A. Three causes of variation in the photochemical reflectance index (PRI) in evergreen conifers. New Phytol. 2015, 206, 187–195. [Google Scholar] [CrossRef] [PubMed]
  15. Porcar-Castell, A.; Tyystjarvi, E.; Atherton, J.; van der Tol, C.; Flexas, J.; Pfundel, E.E.; Moreno, J.; Frankenberg, C.; Berry, J.A. Linking chlorophyll a fluorescence to photosynthesis for remote sensing applications: Mechanisms and challenges. J. Exp. Bot. 2014, 65, 4065–4095. [Google Scholar] [CrossRef] [PubMed]
  16. Plascyk, J.A. The MK II Fraunhofer line discriminator (FLD-II) for airborne and orbital remote sensing of solar-stimulated luminescence. Opt. Eng. 1975, 14, 339–346. [Google Scholar]
  17. Plascyk, J.A.; Gabriel, F.C. Fraunhofer Line Discriminator Mkii—Airborne Instrument For Precise And Standardized Ecological Luminescence Measurement. IEEE Trans. Instrum. Meas. 1975, 24, 306–313. [Google Scholar] [CrossRef]
  18. Meroni, M.; Colombo, R. Leaf level detection of solar induced chlorophyll fluorescence by means of a subnanometer resolution spectroradiometer. Remote Sens. Environ. 2006, 103, 438–448. [Google Scholar] [CrossRef]
  19. Alonso, L.; Gomez-Chova, L.; Vila-Frances, J.; Amoros-Lopez, J.; Guanter, L.; Calpe, J.; Moreno, J. Improved Fraunhofer Line Discrimination Method for Vegetation Fluorescence Quantification. IEEE Geosci. Remote Sens. Lett. 2008, 5, 620–624. [Google Scholar] [CrossRef]
  20. Meroni, M.; Rossini, M.; Guanter, L.; Alonso, L.; Rascher, U.; Colombo, R.; Moreno, J. Remote sensing of solar-induced chlorophyll fluorescence: Review of methods and applications. Remote Sens. Environ. 2009, 113, 2037–2051. [Google Scholar] [CrossRef]
  21. Adams, W.W.; Demmigadams, B.; Winter, K.; Schreiber, U. The Ratio of Variable to Maximum Chlorophyll Fluorescence From Photosystem-II, Measured in Leaves at Ambient-Temperature and at 77 k, as an Indicator of The Photon Yield of Photosynthesis. Planta 1990, 180, 166–174. [Google Scholar] [CrossRef] [PubMed]
  22. Logan, B.A.; Adams, W.W.; Demmig-Adams, B. Viewpoint: Avoiding common pitfalls of chlorophyll fluorescence analysis under field conditions. Funct. Plant Boil. 2007, 34, 853–859. [Google Scholar] [CrossRef]
  23. Campbell, P.K.E.; Middleton, E.M.; McMurtrey, J.E.; Corp, L.A.; Chappelle, E.W. Assessment of vegetation stress using reflectance or fluorescence measurements. J. Environ. Qual. 2007, 36, 832–845. [Google Scholar] [CrossRef] [PubMed]
  24. Campbell, P.K.E.; Middleton, E.M.; Corp, L.A.; Kim, M.S. Contribution of chlorophyll fluorescence to the apparent vegetation reflectance. Sci. Total. Environ. 2008, 404, 433–439. [Google Scholar] [CrossRef] [PubMed]
  25. Yang, X.; Tang, J.W.; Mustard, J.F.; Lee, J.E.; Rossini, M.; Joiner, J.; Munger, J.W.; Kornfeld, A.; Richardson, A.D. Solar-induced chlorophyll fluorescence that correlates with canopy photosynthesis on diurnal and seasonal scales in a temperate deciduous forest. Geophys. Res. Lett. 2015, 42, 2977–2987. [Google Scholar] [CrossRef] [Green Version]
  26. Frankenberg, C.; Fisher, J.B.; Worden, J.; Badgley, G.; Saatchi, S.S.; Lee, J.E.; Toon, G.C.; Butz, A.; Jung, M.; Kuze, A.; et al. New global observations of the terrestrial carbon cycle from GOSAT: Patterns of plant fluorescence with gross primary productivity. Geophys. Res. Lett. 2011, 38. [Google Scholar] [CrossRef] [Green Version]
  27. Joiner, J.; Yoshida, Y.; Vasilkov, A.P.; Yoshida, Y.; Corp, L.A.; Middleton, E.M. First observations of global and seasonal terrestrial chlorophyll fluorescence from space. Biogeosciences 2011, 8, 637–651. [Google Scholar] [CrossRef] [Green Version]
  28. Joiner, J.; Yoshida, Y.; Vasilkov, A.P.; Middleton, E.M.; Campbell, P.K.E.; Yoshida, Y.; Kuze, A.; Corp, L.A. Filling-in of near-infrared solar lines by terrestrial fluorescence and other geophysical effects: Simulations and space-based observations from SCIAMACHY and GOSAT. Atmos. Meas. Tech. 2012, 5, 809–829. [Google Scholar] [CrossRef]
  29. Guanter, L.; Rossini, M.; Colombo, R.; Meroni, M.; Frankenberg, C.; Lee, J.E.; Joiner, J. Using field spectroscopy to assess the potential of statistical approaches for the retrieval of sun-induced chlorophyll fluorescence from ground and space. Remote Sens. Environ. 2013, 133, 52–61. [Google Scholar] [CrossRef]
  30. Middleton, E.M.; Rascher, U.; Corp, L.A.; Huemmrich, K.F.; Cook, B.D.; Noormets, A.; Schickling, A.; Pinto, F.; Alonso, L.; Damm, A.; et al. The 2013 FLEX-US Airborne Campaign at the Parker Tract Loblolly Pine Plantation in North Carolina, USA. Remote Sens. 2017, 9, 612. [Google Scholar] [CrossRef]
  31. Sun, Y.; Frankenberg, C.; Wood, J.D.; Schimel, D.S.; Jung, M.; Guanter, L.; Drewry, D.T.; Verma, M.; Porcar-Castell, A.; Griffis, T.J.; et al. OCO-2 advances photosynthesis observation from space via solar-induced chlorophyll fluorescence. Science 2017, 358. [Google Scholar] [CrossRef] [PubMed]
  32. Gentine, P.; Alemohammad, S.H. Reconstructed Solar-Induced Fluorescence: A Machine Learning Vegetation Product Based on MODIS Surface Reflectance to Reproduce GOME-2 Solar-Induced Fluorescence. Geophys. Res. Lett. 2018, 45, 3136–3146. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Li, X.; Xiao, J.F.; He, B.B. Chlorophyll fluorescence observed by OCO-2 is strongly related to gross primary productivity estimated from flux towers in temperate forests. Remote Sens. Environ. 2018, 204, 659–671. [Google Scholar] [CrossRef]
  34. Sun, Y.; Frankenberg, C.; Jung, M.; Joiner, J.; Guanter, L.; Kohler, P.; 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. 2018, 209, 808–823. [Google Scholar] [CrossRef]
  35. Guanter, L.; Frankenberg, C.; Dudhia, A.; Lewis, P.E.; Gomez-Dans, J.; Kuze, A.; Suto, H.; Grainger, R.G. Retrieval and global assessment of terrestrial chlorophyll fluorescence from GOSAT space measurements. Remote Sens. Environ. 2012, 121, 236–251. [Google Scholar] [CrossRef]
  36. Joiner, J.; Guanter, L.; Lindstrot, R.; Voigt, M.; Vasilkov, A.P.; Middleton, E.M.; Huemmrich, K.F.; Yoshida, Y.; 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. 2013, 6, 2803–2823. [Google Scholar] [CrossRef]
  37. Frankenberg, C.; O’Dell, C.; Berry, J.; Guanter, L.; Joiner, J.; Kohler, P.; Pollock, R.; Taylor, T.E. Prospects for chlorophyll fluorescence remote sensing from the Orbiting Carbon Observatory-2. Remote Sens. Environ. 2014, 147, 1–12. [Google Scholar] [CrossRef] [Green Version]
  38. Colombo, R.; Celesti, M.; Bianchi, R.; Campbell, P.K.E.; Cogliati, S.; Cook, B.D.; Corp, L.A.; Damm, A.; Domec, J.C.; Guanter, L.; et al. Variability of sun-induced chlorophyll fluorescence according to stand age-related processes in a managed loblolly pine forest. Glob. Chang. Boil. 2018, 24, 2980–2996. [Google Scholar] [CrossRef]
  39. Wohlfahrt, G.; Gerdel, K.; Migliavacca, M.; Rotenberg, E.; Tatarinov, F.; Muller, J.; Hammerle, A.; Julitta, T.; Spielmann, F.M.; Yakir, D. Sun-induced fluorescence and gross primary productivity during a heat wave. Sci. Rep. 2018, 8, 14169. [Google Scholar] [CrossRef]
  40. Perez-Priego, O.; Zarco-Tejada, P.J.; Miller, J.R.; Sepulcre-Canto, G.; Fereres, E. Detection of water stress in orchard trees with a high-resolution spectrometer through chlorophyll fluorescence in-filling of the O-2-A band. IEEE Trans. Geosci. Remote Sens. 2005, 43, 2860–2869. [Google Scholar] [CrossRef]
  41. Rossini, M.; Meroni, M.; Migliavacca, M.; Manca, G.; Cogliati, S.; Busetto, L.; Picchi, V.; Cescatti, A.; Seufert, G.; Colombo, R. High resolution field spectroscopy measurements for estimating gross ecosystem production in a rice field. Agric. For. Meteorol. 2010, 150, 1283–1296. [Google Scholar] [CrossRef]
  42. Daumard, F.; Goulas, Y.; Champagne, S.; Fournier, A.; Ounis, A.; Olioso, A.; Moya, I. Continuous Monitoring of Canopy Level Sun-Induced Chlorophyll Fluorescence During the Growth of a Sorghum Field. IEEE Trans. Geosci. Remote Sens. 2012, 50, 4292–4300. [Google Scholar] [CrossRef]
  43. Zarco-Tejada, P.J.; Catalina, A.; Gonzalez, M.R.; Martin, P. Relationships between net photosynthesis and steady-state chlorophyll fluorescence retrieved from airborne hyperspectral imagery. Remote Sens. Environ. 2013, 136, 247–258. [Google Scholar] [CrossRef] [Green Version]
  44. Rossini, M.; Nedbal, L.; Guanter, L.; Ac, A.; Alonso, L.; Burkart, A.; Cogliati, S.; Colombo, R.; Damm, A.; Drusch, M.; et al. Red and far red Sun-induced chlorophyll fluorescence as a measure of plant photosynthesis. Geophys. Res. Lett. 2015, 42, 1632–1639. [Google Scholar] [CrossRef] [Green Version]
  45. Damm, A.; Elbers, J.; Erler, A.; Gioli, B.; Hamdi, K.; Hutjes, R.; Kosvancova, M.; Meroni, M.; Miglietta, F.; Moersch, A.; et al. Remote sensing of sun-induced fluorescence to improve modeling of diurnal courses of gross primary production (GPP). Glob. Chang. Boil. 2010, 16, 171–186. [Google Scholar] [CrossRef] [Green Version]
  46. Balzarolo, M.; Anderson, K.; Nichol, C.; Rossini, M.; Vescovo, L.; Arriga, N.; Wohlfahrt, G.; Calvet, J.C.; Carrara, A.; Cerasoli, S.; et al. Ground-Based Optical Measurements at European Flux Sites: A Review of Methods, Instruments and Current Controversies. Sensors 2011, 11, 7954–7981. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Cheng, Y.B.; Middleton, E.M.; Zhang, Q.Y.; Huemmrich, K.F.; Campbell, P.K.E.; Corp, L.A.; Cook, B.D.; Kustas, W.P.; Daughtry, C.S. Integrating Solar Induced Fluorescence and the Photochemical Reflectance Index for Estimating Gross Primary Production in a Cornfield. Remote Sens. 2013, 5, 6857–6879. [Google Scholar] [CrossRef] [Green Version]
  48. Rascher, U.; Agati, G.; Alonso, L.; Cecchi, G.; Champagne, S.; Colombo, R.; Damm, A.; Daumard, F.; de Miguel, E.; Fernandez, G.; et al. CEFLES2: The remote sensing component to quantify photosynthetic efficiency from the leaf to the region by measuring sun-induced fluorescence in the oxygen absorption bands. Biogeosciences 2009, 6, 1181–1198. [Google Scholar] [CrossRef]
  49. Zarco-Tejada, P.J.; Gonzalez-Dugo, V.; Berni, J.A.J. Fluorescence, temperature and narrow-band indices acquired from a UAV platform for water stress detection using a micro-hyperspectral imager and a thermal camera. Remote Sens. Environ. 2012, 117, 322–337. [Google Scholar] [CrossRef] [Green Version]
  50. Goulas, Y.; Fournier, A.; Daumard, F.; Champagne, S.; Ounis, A.; Marloie, O.; Moya, I. Gross Primary Production of a Wheat Canopy Relates Stronger to Far Red Than to Red Solar-Induced Chlorophyll Fluorescence. Remote Sens. 2017, 9, 97. [Google Scholar] [CrossRef]
  51. Hari, P.; Kulmala, M. Station for measuring ecosystem-atmosphere relations (SMEAR II). Boreal Environ. Res. 2005, 10, 315–322. [Google Scholar]
  52. Kulmala, L.; Read, J.; Nojd, P.; Rathgeber, C.B.K.; Cuny, H.E.; Hollmen, J.; Makinen, H. Identifying the main drivers for the production and maturation of Scots pine tracheids along a temperature gradient. Agric. For. Meteorol. 2017, 232, 210–224. [Google Scholar] [CrossRef]
  53. Kolari, P.; Kulmala, L.; Pumpanen, J.; Launiainen, S.; Ilvesniem, H.; Hari, P.; Nikinmaa, E. CO2 exchange and component CO2 fluxes of a boreal Scots pine forest. Boreal Environ. Res. 2009, 14, 761–783. [Google Scholar]
  54. Ilvesniemi, H.; Giesler, R.; van Hees, P.; Magnussson, T.; Melkerud, P.A. General description of the sampling techniques and the sites investigated in the Fennoscandinavian podzolization project. Geoderma 2000, 94, 109–123. [Google Scholar] [CrossRef]
  55. Drolet, G.; Wade, T.; Nichol, C.J.; MacLellan, C.; Levula, J.; Porcar-Castell, A.; Nikinmaa, E.; Vesala, T. A temperature-controlled spectrometer system for continuous and unattended measurements of canopy spectral radiance and reflectance. Int. J. Remote Sens. 2014, 35, 1769–1785. [Google Scholar] [CrossRef]
  56. Kovac, D.; Malenovsky, Z.; Urban, O.; Spunda, V.; Kalina, J.; Ac, A.; Kaplan, V.; Hanus, J. Response of green reflectance continuum removal index to the xanthophyll de-epoxidation cycle in Norway spruce needles. J. Exp. Bot. 2013, 64, 1817–1827. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Damm, A.; Erler, A.; Hillen, W.; Meroni, M.; Schaepman, M.E.; Verhoef, W.; Rascher, U. Modeling the impact of spectral sensor configurations on the FLD retrieval accuracy of sun-induced chlorophyll fluorescence. Remote Sens. Environ. 2011, 115, 1882–1892. [Google Scholar] [CrossRef]
  58. Sabater, N.; Vicent, J.; Alonso, L.; Verrelst, J.; Middleton, E.M.; Porcar-Castell, A.; Moreno, J. Compensation of Oxygen Transmittance Effects for Proximal Sensing Retrieval of Canopy-Leaving Sun-Induced Chlorophyll Fluorescence. Remote Sens. 2018, 10, 1551. [Google Scholar] [CrossRef]
  59. Rothman, L.S.; Gordon, I.E.; Babikov, Y.; Barbe, A.; Benner, D.C.; Bernath, P.F.; Birk, M.; Bizzocchi, L.; Boudon, V.; Brown, L.R.; et al. The HITRAN2012 molecular spectroscopic database. J. Quant. Spectrosc. Radiat. Transf. 2013, 130, 4–50. [Google Scholar] [CrossRef] [Green Version]
  60. Kochanov, R.V.; Gordon, I.E.; Rothman, L.S.; Wcislo, P.; Hill, C.; Wilzewski, J.S. HITRAN Application Programming Interface (HAPI): A comprehensive approach to working with spectroscopic data. J. Quant. Spectrosc. Radiat. Transf. 2016, 177, 15–30. [Google Scholar] [CrossRef]
  61. Porcar-Castell, A.; Pfundel, E.; Korhonen, J.F.J.; Juurola, E. A new monitoring PAM fluorometer (MONI-PAM) to study the short- and long-term acclimation of photosystem II in field conditions. Photosynth. Res. 2008, 96, 173–179. [Google Scholar] [CrossRef] [PubMed]
  62. Porcar-Castell, A. A high-resolution portrait of the annual dynamics of photochemical and non-photochemical quenching in needles of Pinus sylvestris. Physiol. Plant. 2011, 143, 139–153. [Google Scholar] [CrossRef] [PubMed]
  63. Damm, A.; Guanter, L.; Paul-Limoges, E.; van der Tol, C.; Hueni, A.; Buchmann, N.; Eugster, W.; Ammann, C.; Schaepman, M.E. Far-red sun-induced chlorophyll fluorescence shows ecosystem-specific relationships to gross primary production: An assessment based on observational and modeling approaches. Remote Sens. Environ. 2015, 166, 91–105. [Google Scholar] [CrossRef]
  64. Liu, L.Y.; Liu, X.J.; Wang, Z.H.; Zhang, B. Measurement and Analysis of Bidirectional SIF Emissions in Wheat Canopies. IEEE Trans. Geosci. Remote Sens. 2016, 54, 2640–2651. [Google Scholar] [CrossRef]
  65. Miao, G.F.; Guan, K.Y.; Yang, X.; Bernacchi, C.J.; Berry, J.A.; DeLucia, E.H.; Wu, J.; Moore, C.E.; Meacham, K.; Cai, Y.P.; et al. Sun-Induced Chlorophyll Fluorescence, Photosynthesis, and Light Use Efficiency of a Soybean Field from Seasonally Continuous Measurements. J. Geophys. Res. Biogeosci. 2018, 123, 610–623. [Google Scholar] [CrossRef]
  66. Zhang, Y.G.; Guanter, L.; Joiner, J.; Song, L.; Guan, K.Y. Spatially-explicit monitoring of crop photosynthetic capacity through the use of space-based chlorophyll fluorescence data. Remote Sens. Environ. 2018, 210, 362–374. [Google Scholar] [CrossRef]
  67. Guanter, L.; Zhang, Y.G.; Jung, M.; Joiner, J.; Voigt, M.; Berry, J.A.; Frankenberg, C.; Huete, A.R.; Zarco-Tejada, P.; Lee, J.E.; et al. Global and time-resolved monitoring of crop photosynthesis with chlorophyll fluorescence. Proc. Natl. Acad. Sci. USA 2014, 111, E1327–E1333. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Joiner, J.; Yoshida, Y.; Vasilkov, A.; Schaefer, K.; Jung, M.; Guanter, L.; Zhang, Y.; Garrity, S.; Middleton, E.M.; Huemmrich, K.F.; et al. The seasonal cycle of satellite chlorophyll fluorescence observations and its relationship to vegetation phenology and ecosystem atmosphere carbon exchange. Remote Sens. Environ. 2014, 152, 375–391. [Google Scholar] [CrossRef] [Green Version]
  69. Guanter, L.; Kohler, P.; Walther, S.; Zhang, Y.G. Recent advances in global monitoring of terrestrial sun-induced chlorophyll fluorescence. In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium, Beijing, China, 10–15 July 2016; pp. 1714–1716. [Google Scholar] [CrossRef]
  70. Sanders, A.F.J.; Verstraeten, W.W.; Kooreman, M.L.; van Leth, T.C.; Beringer, J.; Joiner, J. Spaceborne Sun-Induced Vegetation Fluorescence Time Series from 2007 to 2015 Evaluated with Australian Flux Tower Measurements. Remote Sens. 2016, 8, 895. [Google Scholar] [CrossRef]
  71. Kohler, P.; Frankenberg, C.; Magney, T.S.; Guanter, L.; Joiner, J.; Landgraf, J. Global Retrievals of Solar-Induced Chlorophyll Fluorescence With TROPOMI: First Results and Intersensor Comparison to OCO-2. Geophys. Res. Lett. 2018, 45, 10456–10463. [Google Scholar] [CrossRef]
  72. Zuromski, L.M.; Bowling, D.R.; Kohler, P.; Frankenberg, C.; Goulden, M.L.; Blanken, P.D.; Lin, J.C. Solar-Induced Fluorescence Detects Interannual Variation in Gross Primary Production of Coniferous Forests in the Western United States. Geophys. Res. Lett. 2018, 45, 7184–7193. [Google Scholar] [CrossRef]
  73. Parazoo, N.C.; Bowman, K.; Fisher, J.B.; Frankenberg, C.; Jones, D.B.A.; Cescatti, A.; Perez-Priego, O.; Wohlfahrt, G.; Montagnani, L. Terrestrial gross primary production inferred from satellite fluorescence and vegetation models. Glob. Chang. Boil. 2014, 20, 3103–3121. [Google Scholar] [CrossRef] [PubMed]
  74. 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.; et al. Forest productivity and water stress in Amazonia: Observations from GOSAT chlorophyll fluorescence. Proc. R. Soc. B-Boil. Sci. 2013, 280. [Google Scholar] [CrossRef] [PubMed]
  75. He, L.M.; Chen, J.M.; Liu, J.; Mo, G.; Joiner, J. Angular normalization of GOME-2 Sun-induced chlorophyll fluorescence observation as a better proxy of vegetation productivity. Geophys. Res. Lett. 2017, 44, 5691–5699. [Google Scholar] [CrossRef]
  76. Zhang, Y.G.; Guanter, L.; Berry, J.A.; van der Tol, C.; Yang, X.; Tang, J.W.; Zhang, F.M. Model-based analysis of the relationship between sun-induced chlorophyll fluorescence and gross primary production for remote sensing applications. Remote Sens. Environ. 2016, 187, 145–155. [Google Scholar] [CrossRef]
  77. Zhang, Y.; Xiao, X.M.; Wolf, S.; Wu, J.; Wu, X.C.; Gioli, B.; Wohlfahrt, G.; Cescatti, A.; van der Tol, C.; Zhou, S.; et al. Spatio-Temporal Convergence of Maximum Daily Light-Use Efficiency Based on Radiation Absorption by Canopy Chlorophyll. Geophys. Res. Lett. 2018, 45, 3508–3519. [Google Scholar] [CrossRef]
  78. Walther, S.; Voigt, M.; Thum, T.; Gonsamo, A.; Zhang, Y.G.; Kohler, P.; Jung, M.; Varlagin, A.; Guanter, L. Satellite chlorophyll fluorescence measurements reveal large-scale decoupling of photosynthesis and greenness dynamics in boreal evergreen forests. Glob. Chang. Boil. 2016, 22, 2979–2996. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Du, S.S.; Liu, L.Y.; Liu, X.J.; Hu, J.C. Response of Canopy Solar-Induced Chlorophyll Fluorescence to the Absorbed Photosynthetically Active Radiation Absorbed by Chlorophyll. Remote Sens. 2017, 9, 911. [Google Scholar] [CrossRef]
  80. Li, X.; Xiao, J.F.; He, B.B. Higher absorbed solar radiation partly offset the negative effects of water stress on the photosynthesis of Amazon forests during the 2015 drought. Environ. Res. Lett. 2018, 13. [Google Scholar] [CrossRef]
  81. Van der Tol, C.; Berry, J.A.; Campbell, P.K.E.; Rascher, U. Models of fluorescence and photosynthesis for interpreting measurements of solar-induced chlorophyll fluorescence. J. Geophys. Res. Biogeosci. 2014, 119, 2312–2327. [Google Scholar] [CrossRef] [Green Version]
  82. Hilker, T.; Coops, N.C.; Wulder, M.A.; Black, T.A.; Guy, R.D. The use of remote sensing in light use efficiency based models of gross primary production: A review of current status and future requirements. Sci. Total. Environ. 2008, 404, 411–423. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Van der Tol, C.; Verhoef, W.; Timmermans, J.; Verhoef, A.; Su, Z. An integrated model of soil-canopy spectral radiances, photosynthesis, fluorescence, temperature and energy balance. Biogeosciences 2009, 6, 3109–3129. [Google Scholar] [CrossRef] [Green Version]
  84. Zhang, Z.Y.; Zhang, Y.G.; Joiner, J.; Migliavacca, M. Angle matters: Bidirectional effects impact the slope of relationship between gross primary productivity and sun-induced chlorophyll fluorescence from Orbiting Carbon Observatory-2 across biomes. Glob. Chang. Boil. 2018, 24, 5017–5020. [Google Scholar] [CrossRef] [PubMed]
  85. Li, X.; Xiao, J.F.; He, B.B.; Arain, M.A.; Beringer, J.; Desai, A.R.; Emmel, C.; Hollinger, D.Y.; Krasnova, A.; Mammarella, I.; et al. Solar-induced chlorophyll fluorescence is strongly correlated with terrestrial photosynthesis for a wide variety of biomes: First global analysis based on OCO-2 and flux tower observations. Glob. Chang. Boil. 2018, 24, 3990–4008. [Google Scholar] [CrossRef] [PubMed]
  86. Daumard, F.; Goulas, Y.; Ounis, A.; Pedros, R.; Moya, I. Measurement and Correction of Atmospheric Effects at Different Altitudes for Remote Sensing of Sun-Induced Fluorescence in Oxygen Absorption Bands. IEEE Trans. Geosci. Remote Sens. 2015, 53, 5180–5196. [Google Scholar] [CrossRef]
  87. Rascher, U.; Alonso, L.; Burkart, A.; Cilia, C.; Cogliati, S.; Colombo, R.; Damm, A.; Drusch, M.; Guanter, L.; Hanus, J.; et al. Sun-induced fluorescence—A new probe of photosynthesis: First maps from the imaging spectrometer HyPlant. Glob. Chang. Boil. 2015, 21, 4673–4684. [Google Scholar] [CrossRef] [PubMed]
  88. Pierluissi, J.H.; Tsai, C.M. Molecular transmittance band model for oxygen in the visible. Appl. Opt. 1986, 25, 2458–2460. [Google Scholar] [CrossRef] [PubMed]
  89. Julitta, T.; Corp, L.A.; Rossini, M.; Burkart, A.; Cogliati, S.; Davies, N.; Hom, M.; Mac Arthur, A.; Middleton, E.M.; Rascher, U.; et al. Comparison of Sun-Induced Chlorophyll Fluorescence Estimates Obtained from Four Portable Field Spectroradiometers. Remote Sens. 2016, 8, 122. [Google Scholar] [CrossRef]
  90. Zhao, F.; Li, R.; Verhoe, W.; Cogliati, S.; Liu, X.J.; Huang, Y.B.; Guo, Y.Q.; Huang, J.X. Reconstruction of the full spectrum of solar-induced chlorophyll fluorescence: Intercomparison study for a novel method. Remote Sens. Environ. 2018, 219, 233–246. [Google Scholar] [CrossRef]
  91. Garzonio, R.; Di Mauro, B.; Colombo, R.; Cogliati, S. Surface Reflectance and Sun-Induced Fluorescence Spectroscopy Measurements Using a Small Hyperspectral UAS. Remote Sens. 2017, 9, 472. [Google Scholar] [CrossRef]
Figure 1. Seasonal course of meteorological conditions measured at the SMEAR II station during the sampling period. Top panel shows 30 min means of photosynthetic photon flux density (PPFD). The time series in the bottom panel shows 30 min means of air temperature measured at 17 m above the ground while vertical bars represent 30 min sums of precipitation.
Figure 1. Seasonal course of meteorological conditions measured at the SMEAR II station during the sampling period. Top panel shows 30 min means of photosynthetic photon flux density (PPFD). The time series in the bottom panel shows 30 min means of air temperature measured at 17 m above the ground while vertical bars represent 30 min sums of precipitation.
Remotesensing 11 00273 g001
Figure 2. Sample irradiance and radiance spectra measured at 15-min intervals by the Ocean Optics spectrometer system used in this study. Each spectra is the average of 25 consecutive scans. The grey area shows the O2-A absorption region around 760 nm used to extract SiF in this study. Note the difference between the ranges of the left and right Y axes, emphasizing the weak reflected signal from the Scots pine canopy.
Figure 2. Sample irradiance and radiance spectra measured at 15-min intervals by the Ocean Optics spectrometer system used in this study. Each spectra is the average of 25 consecutive scans. The grey area shows the O2-A absorption region around 760 nm used to extract SiF in this study. Note the difference between the ranges of the left and right Y axes, emphasizing the weak reflected signal from the Scots pine canopy.
Remotesensing 11 00273 g002
Figure 3. Complete (day and night, all sky conditions) time series of instantaneous observations of (A) SiF and (B) SiF L650−1, both without correction for structural effects and centered on the flux data acquisitions, (C) 30 min means of gross primary productivity (GPP) calculated from the flux tower eddy covariance data, (D) effective quantum yield ( Δ F / F m ) measured by the PAM system, (E) the normalized difference vegetation index (NDVI), and (F) the photochemical reflectance index (PRI), both measured by the Skye sensors mounted on the SMEAR II tall tower. A malfunction of the PAM system resulted in a gap in the Δ F / F m time series between days of the year 207 and 231.
Figure 3. Complete (day and night, all sky conditions) time series of instantaneous observations of (A) SiF and (B) SiF L650−1, both without correction for structural effects and centered on the flux data acquisitions, (C) 30 min means of gross primary productivity (GPP) calculated from the flux tower eddy covariance data, (D) effective quantum yield ( Δ F / F m ) measured by the PAM system, (E) the normalized difference vegetation index (NDVI), and (F) the photochemical reflectance index (PRI), both measured by the Skye sensors mounted on the SMEAR II tall tower. A malfunction of the PAM system resulted in a gap in the Δ F / F m time series between days of the year 207 and 231.
Remotesensing 11 00273 g003
Figure 4. Seasonal patterns of a midday average (10–12 UTC) of: GPP30 from flux tower data (A), SiF (B) and SiF L650−1 (C) from canopy level measurements from the Ocean optics spectrometers, SIF/APAR (D), PRI, (E) NDVI (F) PRI, and ΔF/Fm′ (G) from the PAM system; LUE (H), (I) SiF L650-1 to SiF/APAR. Rain and snow days were removed. Clear and cloudy observations are presented.
Figure 4. Seasonal patterns of a midday average (10–12 UTC) of: GPP30 from flux tower data (A), SiF (B) and SiF L650−1 (C) from canopy level measurements from the Ocean optics spectrometers, SIF/APAR (D), PRI, (E) NDVI (F) PRI, and ΔF/Fm′ (G) from the PAM system; LUE (H), (I) SiF L650-1 to SiF/APAR. Rain and snow days were removed. Clear and cloudy observations are presented.
Remotesensing 11 00273 g004
Figure 5. Relationships between SiF and SiF L650−1, and GPP30 and LUE30 (10–12 UTC). Grey circles include data acquired on clear and cloudy days as well as all solar angles. Black circles represent data acquired during clear days (A,C,E,G), and also when solar zenith angle was smaller than 50° (n = 118 for B,D, 116 for F,H) on clear days (B,D,F,H).
Figure 5. Relationships between SiF and SiF L650−1, and GPP30 and LUE30 (10–12 UTC). Grey circles include data acquired on clear and cloudy days as well as all solar angles. Black circles represent data acquired during clear days (A,C,E,G), and also when solar zenith angle was smaller than 50° (n = 118 for B,D, 116 for F,H) on clear days (B,D,F,H).
Remotesensing 11 00273 g005
Figure 6. Diurnal relationship in SiF (-), GPP (solid circles), PPFD (open circles) and SZA (-) on select clear days throughout the growing season. (A) April 28th, (B) May 10th, (C) June 5th, (D) July 29th, (E) August 27th, and (F) September 5th.
Figure 6. Diurnal relationship in SiF (-), GPP (solid circles), PPFD (open circles) and SZA (-) on select clear days throughout the growing season. (A) April 28th, (B) May 10th, (C) June 5th, (D) July 29th, (E) August 27th, and (F) September 5th.
Remotesensing 11 00273 g006
Figure 7. Diurnal course of SiF L650−1 (black circles), quantum yield of PSII (ΔF/Fm′, black triangles) and solar zenith angle (SZA, blue curve) for the same clear days shown in Figure 6. (A) April 28th, (B) May 10th, (C) June 5th, (D) August 27th, and (E) September 5th. Note that there are no PAM data for July. Solar noon indicated by the grey line.
Figure 7. Diurnal course of SiF L650−1 (black circles), quantum yield of PSII (ΔF/Fm′, black triangles) and solar zenith angle (SZA, blue curve) for the same clear days shown in Figure 6. (A) April 28th, (B) May 10th, (C) June 5th, (D) August 27th, and (E) September 5th. Note that there are no PAM data for July. Solar noon indicated by the grey line.
Remotesensing 11 00273 g007
Table 1. Parameter estimates, R-square values and significance levels from linear regressions for the four lower panels in Figure 5.
Table 1. Parameter estimates, R-square values and significance levels from linear regressions for the four lower panels in Figure 5.
SlopeInterceptR2Std ErrorFdfp-Value
B0.00020.00030.260.00139.76117<0.0001
D0.02510.44680.250.16638.15117<0.0001
F0.21230.0032 0.240.00137.18115<0.0001
H34.62130.40910.290.16344.72115<0.0001

Share and Cite

MDPI and ACS Style

Nichol, C.J.; Drolet, G.; Porcar-Castell, A.; Wade, T.; Sabater, N.; Middleton, E.M.; MacLellan, C.; Levula, J.; Mammarella, I.; Vesala, T.; et al. Diurnal and Seasonal Solar Induced Chlorophyll Fluorescence and Photosynthesis in a Boreal Scots Pine Canopy. Remote Sens. 2019, 11, 273. https://doi.org/10.3390/rs11030273

AMA Style

Nichol CJ, Drolet G, Porcar-Castell A, Wade T, Sabater N, Middleton EM, MacLellan C, Levula J, Mammarella I, Vesala T, et al. Diurnal and Seasonal Solar Induced Chlorophyll Fluorescence and Photosynthesis in a Boreal Scots Pine Canopy. Remote Sensing. 2019; 11(3):273. https://doi.org/10.3390/rs11030273

Chicago/Turabian Style

Nichol, Caroline J., Guillaume Drolet, Albert Porcar-Castell, Tom Wade, Neus Sabater, Elizabeth M. Middleton, Chris MacLellan, Janne Levula, Ivan Mammarella, Timo Vesala, and et al. 2019. "Diurnal and Seasonal Solar Induced Chlorophyll Fluorescence and Photosynthesis in a Boreal Scots Pine Canopy" Remote Sensing 11, no. 3: 273. https://doi.org/10.3390/rs11030273

APA Style

Nichol, C. J., Drolet, G., Porcar-Castell, A., Wade, T., Sabater, N., Middleton, E. M., MacLellan, C., Levula, J., Mammarella, I., Vesala, T., & Atherton, J. (2019). Diurnal and Seasonal Solar Induced Chlorophyll Fluorescence and Photosynthesis in a Boreal Scots Pine Canopy. Remote Sensing, 11(3), 273. https://doi.org/10.3390/rs11030273

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop