Abstract

We explored ability of reflectance vegetation indexes (VIs) related to chlorophyll fluorescence emission (, ) and de-epoxidation state of xanthophyll cycle pigments (PRI, calculated as ) to track changes in the CO2 assimilation rate and Light Use Efficiency (LUE) in montane grassland and Norway spruce forest ecosystems, both at leaf and also canopy level. VIs were measured at two research plots using a ground-based high spatial/spectral resolution imaging spectroscopy technique. No significant relationship between VIs and leaf light-saturated CO2 assimilation () was detected in instantaneous measurements of grassland under steady-state irradiance conditions. Once the temporal dimension and daily irradiance variation were included into the experimental setup, statistically significant changes in VIs related to tested physiological parameters were revealed. ฮ”PRI and ฮ”() of grassland plant leaves under dark-to-full sunlight transition in the scale of minutes were significantly related to (). In the daily course, the variation of VIs measured in one-hour intervals correlated well with the variation of Gross Primary Production (GPP), Net Ecosystem Exchange (NEE), and LUE estimated via the eddy-covariance flux tower. Statistical results were weaker in the case of the grassland ecosystem, with the strongest statistical relation of the index with NEE and GPP.

1. Introduction

The consequences of the ongoing global climate change [1] might turn terrestrial ecosystems over large areas from carbon sink to carbon source [2โ€“4]; therefore, an accurate assessment of the global carbon balance is of increasing importance. Remote sensing (RS) provides the only global and cost effective tool to monitor spatiotemporal properties of the CO2 assimilation in various plant ecosystems [5โ€“9]. One of the traditional approaches to quantify RS information is the transformation of reflectance spectra into vegetation indexes (VIs). VIs are mathematical transformations of reflectance at specifically selected spectral bands that maximize sensitivity to target biophysical variables and minimize confounding environmental factors (e.g., [10]). Since the start of the satellite era-several VIs have been developed for remote quantification of leaf chlorophyll content (Chl(a+b)) [11], leaf area index [12], and/or other biophysical variables [13] that are important for the assessment of the health status and functioning of terrestrial ecosystems. However, most of these derived variables and VIs are rather insensitive to daily dynamical changes in ecosystem carbon fluxes (however, see [14]).

Recent advances in RS provided the possibility to estimate physiological parameters related to the dynamic processes of CO2 fixation by photosynthesis [15โ€“17], namely, to (i) xanthophyll cycle pigments and (ii) chlorophyll fluorescence (Chl-F). Such approach of direct estimation of the light use efficiency (LUE) has a potential to contribute to a significant reduction of existing large uncertainties in the estimation of the carbon balance [16].

The xanthophyll cycle pigments protect the photosynthetic apparatus from overexcitation [18]. Under high light conditions, plants activate a photoprotective pathway via deepoxidative reactions of xanthophyll pigments, resulting in zeaxanthin as the final product (e.g., [19]). This process enables the dissipation of excessive energy in the form of heat [20] and is associated, along with other processes, with the reduction in photochemical efficiency of Photosystem II (PSII) [18]. These rapid and light intensity-dependent changes were firstly detected remotely as changes in absorbance around 515โ€‰nm [21] and later as changes in reflectance around 530โ€‰nm [22]. Subsequently, Gamon et al. [23] proposed Photochemical Reflectance Index (also called Plant Physiological Reflectance Index, PRI, calculated as (R531 โˆ’ R570)/(R531 + R570)), where is a reflectance intensity at subscripted wavelength. This index has been related to changes in the zeaxanthin pigment content [23โ€“25], LUE (defined as the amount of assimilated carbon per unit of incoming light) measured at leaf [25] and later also at canopy [26โ€“30] and even at satellite levels [7, 31, 32]. Several studies linked PRI measured by eddy-covariance flux towers directly to carbon sink capacity in grassland [33] and forest ecosystems [34]. PRI index has been also shown to correlate with the foliar pigment content of Chl(a+b), carotenoids (Carx+c), Car/Chl ratio and Car/Chl ratio [35, 36] and various Chl-F parameters, especially under the drought stress conditions [37โ€“40].

The Chl-F emission is directly linked to the first steps of photosynthesis, to photochemistry, and it is widely used in photosynthesis research. Chl-F is produced in plant photosynthetic tissues after the absorption of light energy as one of the two deexcitation pathways competing with photochemistry and heat loss [41โ€“43]. The coupling between the photochemistry and Chl-F is the strongest in PSII. Thus, the competition between Chl-F and photosynthesis makes the Chl-F an ideal noninvasive reporter of the photosynthetic activity in plant tissue. The Chl-F emission reaches its maximum intensity at 690 and 740โ€‰nm [44]. Also changes in the zeaxanthin concentration are related to Chl-F, particularly to Nonphotochemical Quenching (NPQ) [45].

Chl-F expresses dynamic changes when a previously dark adapted healthy plant is exposed to actinic light [33]. After the illumination Chl-F increases to a high peak value, it then decreases to a steady state level (FS) [46]. Under saturating light, the peak Chl-F equals to the maximum level (FM). The relative Chl-F decrease (RFd) from the dark-to-light adapted state, quantified as RFd = (FM โˆ’ FS)/FS, was shown to be linearly correlated with the light-saturated photosynthetic rate [47].

However, the only sun-induced Chl-F signal that can be tracked via passive RS sensors is โ€œsteady-stateโ€ chlorophyll fluorescence (Chl-FS) [48โ€“50]. Results of several studies suggested that solar-induced Chl-FS [51โ€“55] is linked to CO2 assimilation. Consequently, several VIs, particularly, the reflectance ratios R686/R630, R690/R630, and R740/R800 have been proposed to quantify Chl-FS [37, 50], where the first wavelength corresponds to the maximum intensity of Chl-FS emission, whereas the second wavelength is unaffected by Chl-FS thus serving as a normalization factor. The difference ratio R760.59 โˆ’ R759.5 utilizes the Chl-F infilling effect [56].

The Frauenhofer Line Discriminator principle [57] is another method used in RS for Chl-FS signal extraction from reflectance measurements [34]. Nevertheless, most of the studies using the retrieval of Chl-FS from hyperspectral reflectance were carried out either in laboratory conditions or with nonimaging sensors. Only a limited number of studies attempted to map the spatial distribution of carbon fluxes which were implemented using ground-based [58, 59] airborne [60, 61] or spaceborne [7, 31] spectroscopy sensors.

In this study, we explored the ability of reflectance VIs related to the FS signal, such as R686/R630, R740/R800 and PRI, to estimate the CO2 assimilation at canopy level in grassland and forest ecosystems. First, we evaluated whether the heterogeneity of VIs of a measured grassland plot can reveal a heterogeneity in light-saturated CO2 assimilation rates () under steady state conditions. Second, we followed VI changes under a change of irradiance intensity in the scale of minutes and hours and searched for correlations between VI variations and variations of physiological variables, particularly LUE, gLUE (gross light use efficiency; GPP/APAR), net ecosystem exchange (NEE) and gross primary production (GPP), measured using the eddy-covariance systems at the scale of grassland and forest ecosystems. The aim of this study was to explore potentials and limitations of scaling the carbon cycle-related physiological processes, observable via RS methods, from leaf to canopy level. This way the feasibility of larger scale RS observation of the vegetation processes at ecosystem level can be advanced. To our knowledge, this is the first study using the imaging spectroscopy to investigate plant canopy properties in relation to eddy-covariance data in the diurnal course [17].

2. Material and Methods

2.1. Study Site: Montane Grassland and Forest Ecosystems

The montane grassland and forest research plots are located at the experimental study site Bรญlรฝ Kล™รญลพ (Moravian-Silesian Beskydy Mts., the Czech Republic, 18.54ยฐE, 49.49ยฐN, 855โ€‰m and 908โ€‰m a.s.l., resp.), which is part of the CARBOEUROPE flux tower network (http://www.carboeurope.org/). This site is characterized by cool (annual mean temperature of 5.5ยฐC) and humid (annual mean relative air humidity of about 80%) climate with high annual precipitation (ca. 1000โ€“1400โ€‰mm). The geological bedrock is formed by Mesozoic Godula sandstone (flysch type) with ferric podzols.

The most abundant species occurring at the natural grassland plot (association Nardo-Callunetea, class Nardo-Agrostion tenuis) that were investigated are Festuca rubra agg. (L.), Hieracium sp., Plantago sp., Nardus stricta (L.), and Jacea pseudophrygia (C. A. Meyer). Additional details on the experimental grassland and physiological properties of the dominant plant species are presented in Urban et al. [62].

The investigated forest stand was planted in 1981 with 4-year old seedlings of Picea abies (L.) Karst. (99%) and Abies alba (Mill). (1%) on a slope of 11โ€“16ยฐ with south-southwest exposition over an area of 0.062โ€‰km2 [63]. At the time of the experimental measurements (2005), the stand density was approximately 2 600โ€‰trees/ha (hemisurface leaf area index of about 11โ€‰m2โ€‰mโˆ’2), with mean (ยฑ standard deviation) tree height of 8.5โ€‰m (ยฑ1โ€‰mm) and stem diameter at 1.3โ€‰m of 10.1โ€‰m ยฑ (1โ€‰mm).

2.2. Imaging Spectroscopy Measurements and Data Processing

The nadir viewing canopy reflectance was obtained under clear sky sunny conditions with a visible near-infrared hyperspectral Airborne Imaging Spectrometer for Applications (AISA Eagle, Specim, Oulu, FI). In total 260 spectral bands between 400โ€“940โ€‰nm were collected with a full width at half maximum (FWHM) of 2.2โ€‰nm, and with a ground pixel resolution of 1โ€‰cm. The AISA sensor was mounted to the tower at a height of 20 meters (i.e., approximately 10โ€‰m above the top of the canopy). The acquired area of about 50โ€‰m2 was located within a footprint of an eddy-covariance tower system. In the case of grassland, the AISA sensor was attached to a ladder at the height of 4 meters above the top of the canopy, sensing in a total area of about 4โ€‰m2 (located completely within a footprint of an eddy-covariance system), with a ground pixel resolution of about 0.3โ€‰cm. Figure 1. shows the representative part of the Norway spruce forest, which is homogeneous even-aged monoculture, thus the assumption here is that the smaller footprint of reflectance measurements is possible to be compared with the larger footprint of eddy-covariance data. Subsequent airborne measurement covering the whole footprint area has confirmed this assumption for both forest as well as grassland ecosystem (data not shown).

The acquired AISA images were transformed into radiance values using the laboratory-derived calibration files, and converted into at-sensor reflectance images (Figure 1(a)) by means of the empirical line method, using five near-Lambertian calibration panels of known flat reflectance response. An automatic supervised Maximum Likelihood classification [64] was used to distinguish and to mask shaded pixels. To ensure that only photosynthetically active pixels were used in the subsequent statistical analysis, an appropriate threshold of the green normalized difference VI (green NDVI = (R554 โˆ’ R677)/(R554 + R677)) was applied to separate the green sunlit grass leaves from dry litter. Finally, three narrow-band VIs (PRI, R686/R630, R740/R800) related to photosynthetic processes were derived per pixel from the AISA Eagle reflectance scans (Figure 1(b), Table 1). Selection of these VIs is based on our previous work with AISA Eagle, where detailed diurnal hyperspectral measurements were coupled with Chl-F measurements (see [48] for full description). No other optical VIs were analyzed in this study.

2.3. Leaf Gas Exchange Measurements of Grassland Species

The in situ CO2 assimilation rate of intact leaves of Festuca rubra agg., Hieracium sp., Jacea pseudophrygia, Plantago lanceolata, Tarraxacum spp., and Veronica chamaedrys (L.) was measured at saturating light intensities (Photosynthetic Photons Flux Density (PPFD) โ‰ˆ 1300 ยฑ 10โ€‰ฮผmol (photons) m2 sโˆ’1, [63]) with the open infrared gasometrical analyzer Li-6400 (Li-Cor, USA). All environmental conditions inside the assimilation chamber remained unchanged and referred to the field conditions: leaf temperature between 34โ€“37ยฐC, relative humidity between 35โ€“45%. All leaves, still being attached to the plant, were measured in their natural environment and marked. This way, they were later identified on the AISA images, and all corresponding green pixels were averaged in order to minimize the angular anisotropy of the vegetation bidirectional reflectance distribution function (BRDF).

2.4. Canopy Eddy-Covariance Measurements

Fluxes of CO2 and H2O, as well as latent and sensible heat exchange between the grass stand and atmosphere, were measured using closed-path eddy-covariance systems Edisol (University of Edinburgh), containing the infrared gas analyzer (Li-6262, Li-Cor, Lincoln, USA), sonic anemometer Gill R2 (Gill Instruments, UK), and Edisol software package. In the case of the forest research stand, the InSituFlux (InSituFlux, Sweden), consisting of the infrared gas analyzer (LI-7000 Li-Cor, Lincoln, USA), anemometer (Gill R3, Gill Instruments Ltd., Lymington, UK) and EcoFlux software, was used. The eddy-covariance systems provide time series of half an hour integral measurements of CO2 and H2O fluxes. The footprint of the system covers a circular area with a radius of ca 500โ€‰m [65]. The postprocessing of eddy-covariance data was carried out based on the methodology paper Aubinet et al. [66], accommodating several modifications according to the most recent CarboEurope and FLUXNET protocols. The Quality Control (QC) Software performed the spike removal and quality check of the raw signals. After gap filling the GPP and NEE, values were modeled according to Urban et al. [63]. LUE and gLUE were calculated as follows: LUE = NEE/PPFD and gLUE = GPP/PPFD, where PPFD stands for Photosynthetic Photons Flux Density. For LUE and gLUE calculations we used PPFD values instead of commonly used absorbed PPFD, since the amount of transmitted light through the young dense forest stand is negligible (<2%; data not shown).

2.5. Laboratory Analysis of Foliar Pigments

Plant leaves of grassland (n = 24) and current one-year and two-year old needles of Norway spruce trees (Picea abies (L.) Karst., 100โ€‰mg, n = 6) were sampled after the spectral measurements during the morning, midday, and afternoon (i.e., at 9:00, 12:00 and 16:00 of local time). Foliar pigments of samples, frozen into laboratory in liquid nitrogen for transportation, were extracted in 80% acetone with a small amount of MgCO3. After centrifugation at 480โ€‰g for 3โ€‰min. the supernatant was used for spectrophotometric (UV/VIS 550, Unicam, Leeds, UK) estimation of chlorophyll a and b per area, according to Lichtenthaler [67], and also for High-Performance Liquid Chromatography (HPLC) quantification (TSP Analytical, USA) of individual carotenoids, as described in ล troch et al. [68]. The deepoxidation state of the xanthophyll cycle pigments (DEPS) [69] was calculated using conversion factors published by Farber [70] as DEPS = (Z + A)/(V + A + Z) [71], where Z is zeaxanthin, A is antheraxanthin, and V is the violaxanthin content.

2.6. Setup of Field Experiments

The variation of VIs was measured under three irradiance regimes: (i) instantaneous steady-state measurements under full sun irradiation, (ii) variation of VIs within dark-to-light transition in scale of minutes, and (iii) daily variation of VIs measured in one-hour intervals.

Within the frame of the โ€œsteady-stateโ€ experiment (i) we mapped VIs of selected experimental grassland plots (an area of about 4โ€‰m2) at midday under full sunlight. 24 leaves of 6 plant species were labeled (each leaf was located on a separate plant), and their maximal CO2 assimilation rate () was measured under saturating midday irradiance (PPFD โ‰ˆ 1300 ยฑ 10โ€‰ฮผmol (photons) mโˆ’2โ€‰sโˆ’1). These leaves were identified on the AISA image, captured simultaneously with CO2 assimilation rate measurements. The computed VIs of 10โ€“20 sunlit pixels of each leaf were averaged and correlated to the measurements. The experiment was carried out on a clear-sky day on 18th of July 2007 at 13:00 local time.

In the โ€œshade removalโ€ experiment (ii), run in a time scale of minutes, we tested whether the VIs could estimate the CO2 assimilation rate of labeled grass leaves exposed to saturating sun irradiation () not in the static steady state, but in a dynamic irradiation environment inducing the plant photoprotective reactions and consequent onset of photosynthesis. The experimental plot was divided into two parts (twice 0.5 ร— 4โ€‰m). The first part, containing 13 labeled sample leaves, was covered for 30โ€‰min using a black nontransparent blanket in order to dark-adapt the grass canopy prior to the imaging spectroscopy measurement. Our method of darkening enabled a lateral air convection, which prevented the physiological changes being induced by higher than ambient air temperature. The second control plot was exposed to a natural irradiation regime. Two AISA Eagle images of both plots were obtained at the 5th and 90th seconds after blanket removal. Subsequently, we computed the dynamic change of the VIs (ฮ”VIs) during the dark-to-light transition per leaf sample (average of 10โ€“20 pixels) as subtraction VIDARK โˆ’ VILIGHT, where VIDARK and VILIGHT are values obtained from the scan taken ~5 and ~90 seconds after removal of the blanket, respectively. These ฮ”VIs were statistically related to measurements of the sampled leaves.

Finally, the โ€œdaily courseโ€ experiment (iii), conducted at a time scale of hours, was to correlate the daily variations of the VIsโ€™ mean values of both grassland and forest ecosystems with a daily variation of physiological variables measured in half-hour intervals via eddy-covariance flux systems. Grassland hyperspectral image data were obtained on the September 2, 2005 in one-hour intervals from 9:30 to 16:00 local time. The forest measurements were carried out at one-hour intervals from 13:30 to 16:30 on August 30, 2005 and on August 31, 2005 from 9:30 to 16:30 local time. The daily variations of the selected VIs measured as well as the carbon fluxes were jointly investigated for grassland and forest ecosystems.

2.7. Statistical Analysis

The best-fitting models were computed to evaluate the relationship between the VIs and the physiological variables of carbon assimilation. The amount of variability of a dependent variable explained by the independent variable within the established regression model was expressed by the coefficient of determination (R2). The significance of the relationships between the variables was tested using the analysis of variance (ANOVA) at two significance levels: and , respectively. All statistical tests were performed with Statistica 7.0 software (StatSoft Inc., USA).

3. Results

3.1. โ€œSteady-Stateโ€ Experiment

The relationships between of individual leaf samples and their instantaneous VIs under steady-state saturating irradiance are shown as transparent dots in Figures 2(a)โ€“2(c). No statistically significant relationship was found, neither between PRI nor the other reflectance ratios (R686/R630, R740/R800) and .

3.2. โ€œShade Removalโ€ Experiment

Dark dots in Figure 2 present the relationships between of the 24 leaf sample and the relative change of their VIs (ฮ”VIs) extracted from the AISA images during the dynamic transition from the dark-to-light adapted state of grassland plants. Here, significant exponential relationships () were found between and ฮ”(R686/R630) (R2 = 0.51, y = 0.0484e0.1647x) (Figure 1(a)), and ฮ”PRI (R2 = 0.51, y = 0.015e0.0621x) (Figure 1(b)). ฮ”(R740/R800) (Figure 2(c)) showed no significant correlation with . These results indicate that, after the inclusion of temporal dynamic of the irradiance variability in the experimental setup, the relative change in VIs (PRI and R686/R630) positively correlates with the carbon assimilation rate.

3.3. โ€œDaily Courseโ€ Experiment

Differences in the xanthophyll cycle pigmentsโ€™ daily variation and concentration of chlorophylls between grassland and forest stand are shown in Figure 3. Figure 3(a) displays the daily courses of xanthophyll cycle pigments measured in the leaves of grassland plants and shoots of spruce trees. The diurnal behavior is typical for sunny days, that is, the amount of de-epoxidated xanthophyll pigments (i.e., the amount of Z and A) was highest during midday. The chlorophyll content (Chl(a+b)โ€‰โ€‰per unit leaf area) was relatively during the day in both grassland as well as in forest ecosystems (Figure 3(b)). The forest canopy contained a higher relative amount of xanthophyll pigments in the de-epoxidated state (Figure 3(a)), and also a higher absolute amount of total chlorophylls (Figure 3(b)) compared to the grassland ecosystem.

In this experiment, we have only presented the results on the variation of PRI and R686/R630 due to their ability to follow plant dynamic changes as had been proved in the previous โ€œshade-removalโ€ experiment. The daily courses and also the relationships between the daily variations of PRI and R686/R630, as well as the NEE and LUE measurements of grassland and spruce forest are given in Figures 4 and 5, respectively.

In grassland, NEE gradually decreased during the day (Figure 4a), while LUE showed the typical daily course with lowest values around midday (Figure 4b). A different yet consistent pattern was observed for two consecutive days in the forest ecosystem. NEE slightly increased in the morning (till 11:00). The highest values around midday were followed by a sharp decrease in the afternoon (after 14:00) (Figure 5(a)). The daily pattern of LUE was similar for both days, showing constant values during the day and a sharp decrease after 15:00 (Figure 5(b)). Both NEE and LUE were higher in the forest compared to the grassland ecosystem.

The daily course of the R686/R630 ratio in grassland was quite stable, with the highest values at 10:30 and 11:30 (Figure 4c). In the forest R686/R630 increased in the morning hours, reached the highest values during the midday hours, and then rapidly declined (Figure 5(c)). The daily patterns of PRI variations of grassland and forest measurements were similar; however, the time when PRI reached its daily minimum was shifted. The grassland PRI decreased between 10:00 to 12:30 (Figure 4(d)), reached its lowest value at the time of the highest irradiation intensity, and increased again at 13:30 to its morning values. In the forest ecosystem PRI decreased to the lowest value at 11:00 and then steeply increased, reaching maximum values from 13:30 to 15:30. In the late afternoon (15:30โ€“17:00) it decreased again (Figure 5(d)). In grassland we were not able to record these late afternoon measurements due to the mountainous terrain and the shading of the measuring plot in the late afternoon. When comparing the absolute values, PRI for forest was higher than for grassland, which corresponds with a higher xanthophyll cycle pigment content in the de-epoxidated state of spruce needles and with a higher Chl(a+b) content.

The statistically significant relationships (at ) between the daily variations in R686/R630 and NEE (Table 2) were found in both ecosystems. Moreover, a strong relationship was observed also between R686/R630 and LUE in the forest ecosystem (Table 2). PRI showed a significant relationship (at ) only with LUE in grassland (Table 2). The results of the regression analysis between the daily variations of physiological variables derived from eddy-covariance measurements and four VIs are summarized in Table 2. The ratio R686/R630 showed a significant relationship with a high R2 for both GPP and NEE in grassland, and also in the forest. Moreover, R686/R630 was highly related to gLUE and LUE in the forest (Table 2), and to gLUE in grassland (Table 2). The ratio R740/R800 showed a relatively high dependency on all physiological variables examined for forest but none for grassland ecosystems (Table 2).

4. Discussion

Our main target was to determine if VIs can indicate dynamic physiological processes related to CO2 assimilation rate. Under steady-state conditions, typical for RS data aquisition, we were unable to find statistically significant relationships between VIs and , indicating that a spatial scaling from the leaf up to canopy level cannot be pursued in this fashion. Grace et al. [16] pointed out that VIs related to Chl-F have the potential to provide information on the diurnal, and also on seasonal changes in photosynthesis (see also [46]). Once the temporal component with its dynamic irradiation regime had been included in our dark-to-light transition experimental setup, the variability observed in VIs increased, probably as a result of induced changes in the electron transport rate and the onset of the xanthophyll de-epoxidation processes. Consequently, also the statistically significant positive relationships between ฮ”(R686/R630) and ฮ”PRI, and have confirmed the ability of VIs to follow the light-saturated rate of photosynthesis in the investigated species. The reflectance change at wavelengths around 690โ€‰nm was previously shown to correspond with the Chl-F emission from PSII [49, 72]. Therefore, the relative change observed in ฮ”(R686/R630) is comparable to the trends of the โ€œvitality indexโ€ RFd, significantly correlated with of outdoor grown plant species [47, 73] in earlier studies. Guo and Trotter [35] published that ฮ”PRI for plants with contrasting photosynthetic capacities (expressed in their study as PRI at low light intensity subtracted from PRI at saturating light intensity) was negatively related to . In contrast to this, yet similar to the studies by Demmig-Adams and Adams [74], we observed a significant positive relationship, showing that a larger relative change in Z pigment content is associated with lower in both sunlit and shaded leaves. We also observed that a larger ฮ”PRI was associated with a higher chlorophyll index ratio (R750/R710) [75] (see Supplementary Material available online at doi: 10.1100/2012/705872), which might explain its significant relationship with .

A number of studies [25, 29, 30, 76] investigated the possibility to use PRI derived at canopy level as a proxy for LUE in Monteithโ€™s equation [77] estimating the vegetation GPP. For instance, Rahman et al. [7] successfully used PRI derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite sensor for the estimation of LUE and NPP in boreal forest. However, our results showed that even if the negative effects of shade fraction and nongreen vegetation are excluded, PRI might still not be a good indicator of canopy level photosynthetic processes. Our PRI measurements were significantly related only to LUE of natural grassland and yet PRI was able to explain only 42% of LUE variability. Our result suggests that PRI relationship with LUE should be interpreted with caution, especially in ecosystems with diverse species composition, where species specific relationships might complicate the overall result.

The use of PRI at larger scales originated from leaf level studies showing that the PRI variation of leaves is driven by changes in the de-epoxidation state of xanthophyll cycle pigments in the antennae that reduce the amount of light used in leaf CO2 assimilation, and this way protect the leaf from excessive radiation [78, 79]. A similar trend was observed in our flux tower LUE measurements of grassland, which slightly decreased during midday. In contrast to this, the midday depression in LUE was not visible in the forest data. This might be explained by the fact that coniferous canopy contains a large proportion of shaded compared to sun exposed needles, which might be driving the LUE diurnal course. Apparently, VIs extracted from forest AISA images corresponded mainly to the sunlit shoots on top of the canopy. Our results also stressed a necessity to exclude from VI analysis the confounding effects coming from shade fraction and nongreen vegetation. All these, and other structural factors of forest canopy, can disturb a potential relationship between VIs and the eddy-covariance observations. The lack of correlation at the canopy level might, however, be also weakened by the fact that the relative changes in the pool of xanthophylls represent only a fraction of total carotenoid pigment content [80] and that light induced changes in the leaf also include some other processes, such as changes in pH of photosynthetic membrane [80] and conformational changes in chloroplasts [81]. Porcar-Castell et al. [82] provided further evidence that pool of total carotenoids is closely linked to monoterpene emission capacity, which might further complicate the interpretation of spectral response at 531โ€‰nm. Inoue and Peรฑuelas [83] found out that exponential relationship between LUE and PRI was strongly affected by the soil water sontent. Ollinger [84] recently summarized and highlighted structural factors, which are strongly affecting canopy level relationships between VIs and plant biophysical properties. Apart from plant physiological and structural limitations, canopy level PRI is also affected by viewing and illumination geometry [30]. Hernรกndes-Clemente et al. [85] found that structural sensitivity of PRI can be reduced by selecting the reference wavelength at 512โ€‰nm, instead of 570โ€‰nm. However, the experimental results of our study found a significant relationship between canopy VIs and CO2 assimilation-related plant physiological variables when considering their daily cycle variations. Obtained R2 coefficients between PRI and LUE are lower than the ones reported in the recent meta-analysis [17]. This might be caused by a narrower dynamic range of PRI and LUE measured during our experiment.

The narrow-band reflectance ratios of Chl-F wavelengths (R686/R630, R740/R800) performed better than PRI, and stronger relations were observed for forest than for grassland canopy. This might be partially explained by lower forest water stress, which increased the leaf stomatal conductance resulting in a stronger Chl-F emission. Although most of the studies investigated fluorescence VIs in relation to various Chl-F-related parameters [40, 49, 50], several studies reported a positive relationship of photosynthesis-related variables with reflectance ratios related to Chl-F. The laboratory experiment of Dobrowski et al. [37] resulted in a positive relationship between the reflectance ratio R690/R600 and CO2 assimilation (R2 = 0.36). Furthermore, the reflectance ratio R695/R805 measured in a 5-year-old canopy of pine seedlings was significantly related to [86] during the diurnal course. In general, our results showed higher R2 values of the reflectance ratios to LUE than to gLUE, which can be explained by the respiration component included in the โ€œgrossโ€ measurements. The reflectance ratio R686/R630, tested in this study, was found to be a statistically significant proxy of GPP, NEE, and LUE in grassland and also of gLUE in forest ecosystems. Thus, we can hypothesize that R686/R630 might be considered a potential RS indicator of the ecosystem vegetation productivity variables.

5. Conclusions

The first experimental results conducted at the leaf level revealed no statistically significant relation of โ€œprocess-relatedโ€ narrow-band VIs (PRI; R686/R630; R740/R800) to the maximum CO2 assimilation rate measured instantaneously under the steady-state irradiation conditions. Significant correlations were only found once the relationships between VIs and physiological variables were investigated at a temporal scale (minutes and hours), taking into account the changes in irradiance intensity. When taking into account the 1.5-minute duration of the dark-to-light transition, ฮ”(R686/R630) and ฮ”PRI of natural grassland, we were able to establish significant regressions between VIs and the carbon assimilation measured under saturating light conditions. Finally, in the daily course of several hours, the simple reflectance ratio R686/R630, of the VIs investigated, showed the highest R2 value to the measured eddy-covariance ecosystem variables with better statistical relations found in the case of spruce forest than in grassland canopy. Considering that we observed a strong statistical relation to GPP and NEE, we conclude that the R686/R630 ratio may potentially act as a suitable RS indicator of vegetation productivity at the scale of the whole canopy. In general, the results of this study provide an insight relevant to the further development of leaf-canopy upscaling remote sensing approaches as well as basis for more accurate estimation of carbon assimilation and the productivity in natural canopies. However, considering spatial and temporal limitations of the presented study, many problems are yet to be resolved before a remote sensing application at larger spatial scales.

Acknowledgments

This work is part of the research supported by Grants SP/2d1/70/08 (ForChange), SP/2d1/93/07 (CzechTerra), and 522/09/0468 (GA CR). This paper is a product of the CzechGlobe Centre that is being developed within the OP RDI and cofinanced from EU funds and the State Budget of the Czech Republic (CZ.1.05/1.1.00/02.0073). The experimental site Bรญlรฝ Kล™รญลพ is within the National Infrastructure for Carbon Observationsโ€”CzeCOS/ICOS supported by Ministry of Education CR (LM2010007). We also thank anonymous reviewers for useful comments on the paper.

Supplementary Materials

Supplementary Figure: Relationship between the Photochemical Reflectance Index (PRI, (R532 โ€“ R570)/(R532 + R570)) and chlorophyll index (R750/R710) (a) and between PRI (โ—) and ฮ”PRI (โ—) and R750/R710 (b). Asterisks denotes for statistical significance (**P < 0.001), error bars show ยฑ SD (n = 8โ€“10).

  1. Supplementary Figure