Next Article in Journal
Amorphous In–Ga–Zn–O Powder with High Gas Selectivity towards Wide Range Concentration of C2H5OH
Previous Article in Journal
Human Thalamic Somatosensory Nucleus (Ventral Caudal, Vc) as a Locus for Stimulation by INPUTS from Tactile, Noxious and Thermal Sensors on an Active Prosthesis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Effect of Leaf Stacking on Leaf Reflectance and Vegetation Indices Measured by Contact Probe during the Season

by
Eva Neuwirthová
1,
Zuzana Lhotáková
1 and
Jana Albrechtová
1,2,*
1
Department of Experimental Plant Biology, Faculty of Science, Charles University in Prague, Vinicna 5, 128 44 Prague 2, Czech Republic
2
Institute of Botany, Academy of Sciences, Zámek 1, 252 43 Pruhonice, Czech Republic
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(6), 1202; https://doi.org/10.3390/s17061202
Submission received: 13 February 2017 / Revised: 18 May 2017 / Accepted: 18 May 2017 / Published: 24 May 2017
(This article belongs to the Section Remote Sensors)

Abstract

:
The aims of the study were: (i) to compare leaf reflectance in visible (VIS) (400–700 nm), near-infrared (NIR) (740–1140 nm) and short-wave infrared (SWIR) (2000–2400 nm) spectral ranges measured monthly by a contact probe on a single leaf and a stack of five leaves (measurement setup (MS)) of two broadleaved tree species during the vegetative season; and (ii) to test if and how selected vegetation indices differ under these two MS. In VIS, the pigment-related spectral region, the effect of MS on reflectance was negligible. The major influence of MS on reflectance was detected in NIR (up to 25%), the structure-related spectral range; and weaker effect in SWIR, the water-related spectral range. Vegetation indices involving VIS wavelengths were independent of MS while indices combining wavelengths from both VIS and NIR were MS-affected throughout the season. The effect of leaf stacking contributed to weakening the correlation between the leaf chlorophyll content and selected vegetation indices due to a higher leaf mass per area of the leaf sample. The majority of MS-affected indices were better correlated with chlorophyll content in both species in comparison with MS-unaffected indices. Therefore, in terms of monitoring leaf chlorophyll content using the contact probe reflectance measurement, these MS-affected indices should be used with caution, as discussed in the paper. If the vegetation indices are used for assessment of plant physiological status in various times of the vegetative season, then it is essential to take into consideration their possible changes induced by the particular contact probe measurement setup regarding the leaf stacking.

1. Introduction

Leaf optical properties (reflectance, absorbance and transmittance) are defined as a ratio of incoming light, which is reflected, absorbed or transmitted by a leaf. Leaf optical properties, particularly reflectance, are exploited by remote sensing techniques to gain information about photosynthetic pigment contents, water content or other biochemical compounds in vegetation [1,2,3]. In laboratory conditions, optical properties at the leaf level are measured by a spectroradiometer usually equipped with an integrating sphere or with a contact probe. Although some spectroradiometers operate in different spectral ranges, the most common portable spectroradiometers for a laboratory and field measurement of optical properties (e.g., FieldSpec 4 and ASD) operate from visible (VIS) to infrared wavelengths of electromagnetic spectrum (350–2500 nm). A non-contact measurement with an optical sensor mounted in nadir above a sample has been used recently less frequently. Based on a measurement setup we gain information about reflectance and transmittance when using the integrating sphere [4,5,6,7,8] or only about reflectance, when using the contact probe [9,10]. Measurements with a contact probe have an advantage of fast processing comparing to time-consuming measurement with an integration sphere.
For laboratory measurement of leaf optical properties using a contact probe it is possible to use harvested planar leaves separately as individual, single leaves [9] or in a stack of several leaves [11]. Contact reflectance measurement in conifers can be measured on needles attached on twigs [12] or only needles cut off twigs [13]. Under the non-contact measurement setup with an optical sensor mounted less than one meter in nadir of the whole plant sample [9], usually harvested planar leaves in a stack [14] or whole seedlings are used [15]. Apart of the measurement on a single leaf in all mentioned setups, the leaves in the field of view of an optical sensor overlay at least partly. The leaf overlapping introduces additional variability in leaf spectra due to higher LAI (leaf area index) [16] of the sample and accumulation of its mass per an area unit. Furthermore, we suppose a different mode of scattering inside a single leaf in comparison to the mode of scattering at the interface of individual leaves in a stack what may influence the resulting sample reflectance. Spectral data acquired at a leaf level by a spectroradiometer then serve as an input in radiative transfer (RT) models at a leaf level (e.g., PROSPECT [17]) and can be coupled with RT models at a canopy level (e.g., PROSAIL [18] or DART [19]) and provide better understanding of airborne remote sensing data in large-scale studies [20]. From this point of view, the effect of leaf stacking and overlapping on its reflectance is important if the spectra from the leaf level would be scaled up to the canopy level.
Vegetation indices are used for assessment of plant physiological status. Indices are derived from reflectance in absorption maxima of particular biochemical leaf compounds [5,9,14] and correlate with the chlorophyll content [4,5,21,22], carotenoid and anthocyanin contents [23] or a ratio of carotenoids to chlorophyll [24] and water content [9]. The most frequently used vegetation indices are calculated as: (1) simple ratios of reflectance in individual wavelengths (in nm), Rλ (reflectance in λ wavelengths), e.g., BGI (Blue green pigment index) = R450/R550, G (Greenness index) = R554/R677 [25]; or (2) normalized differences of reflectances in individual wavelengths, e.g., PSNDc (Pigment specific normalized difference) = (R800 − R470)/(R800 + R470) [24]. Simple ratio pigment-related indices are usually computed from wavelengths from VIS region where the photosynthetically active radiation is absorbed [26,27]. Therefore, these indices are sensitive to changes in photosynthetic pigments, particularly in the chlorophyll content [28]. For example, absorption properties of chlorophyll in 700 nm and 550 nm [23] were included in CRI 550, or CRI 700 (Carotenoid concentration indices) [22], which are commonly used for chlorophyll content estimates.
Normalized difference indices for chlorophyll estimate include red edge wavelengths (680–750 nm), where the reflectance sharply increases [29], and wavelengths from near infrared (NIR) region (740–1140 nm) [10]. Particularly, reflectance in NIR is insensitive to changes in the chlorophyll content, however, sensitive to internal leaf structure and water content [6]. Thus, indices including wavelengths from NIR enable an efficient pigment content estimate, e.g., OSAVI [30], Gitelson2, Datt [31], PSSRa, and PSNDb [24], while these indices do not take into account changes in leaf anatomical properties. On the other hand, accumulation of leaf structural compounds contributes to reflectance changes in infrared wavelengths (NIR, SWIR). In spectroscopic studies the leaf content of dry matter and structural compounds is usually expressed as leaf mass per area (LMA) [32]. Vegetation indices can also serve for empirical LMA retrieval [33]. However, the physical approaches are more common, provided that LMA is one of the main parameters for leaf level radiative transfer models, such as PROSPECT [17]. Reflectance in particular regions in NIR is further influenced by leaf water status [9] that is frequently expressed as an equivalent of water thickness (EWT (water content per unit of leaf area), e.g., [34]). The leaf reflectance in NIR also decreases with leaf senescence [4] and depends on the leaf orientation (adaxial or abaxial side of a leaf) towards a sensor [35]. Regarding the leaf orientation, VIS is more effectively absorbed if incoming from the adaxial leaf side than from the abaxial one. The reason is that the chlorophyll is unevenly distributed within a leaf with a dorsiventral structure when more chlorophyll is below the adaxial epidermis in the dense palisade parenchyma [35]. On the other hand, the situation in NIR is different—more radiation is transmitted by the leaf and less is reflected when the leaf is illuminated from the abaxial side than from the adaxial one. This effect is probably caused by the different density of mesophyll cells and intercellular spaces in palisade and spongy parenchyma of a dorsiventral leaf, when dense palisade parenchyma is adjacent to the adaxial side while a loose spongy parenchyma is facing to the abaxial side.
Leaf optical properties change during the season depending on leaf phenology, accumulation of pigments, development of anatomical structure and leaf senescence at the end of the season. However, leaf optical properties and the shape of the spectral curve can change in a different manner in particular spectral ranges—VIS, NIR and SWIR [10] and thus, they can affect differently absolute values of vegetation indices. The reflectances in VIS and NIR decreased up to 5% and 30% as it was shown for Vitis vinifera L. [4] and Picea abies L. Karst [13], respectively, depending on the leaf/needle senescence. While in SWIR, the reflectance increased up to 20% in Populus ssp. [9] and was connected with water deficit [26]. The seasonal difference in the reflectance curve in VIS range influences the values of vegetation indices, which correlate with physiological parameters (e.g., photochemical reflectance index PRI = (R570 − R530)/(R570 + R530) correlating with the light use efficiency) [36]. For example, the normalized difference vegetation index (NDVI) in Quercus rubra L. and Betula papyrifera Marsh. [37] and PRI in selected tropical evergreen species [36] reached their maximum values in July and then values decreased at the end of the season. The strength of the relationship between vegetation indices and the real chlorophyll content may vary during the season and could depend on a leaf type, e.g., planar or needle-like leaves [5]. The relationship between PRI and light use efficiency in Pinus strobus L. decoupled in early spring due to chlorophyll and carotenoids recovery after winter and then it followed a good fit for a linear function [38]. However, there has not been done yet an extensive study comparing numerous vegetation indices based on leaf reflectance acquired by a contact probe during a vegetation season.
Measuring leaf reflectance in a different setup (a single leaf or more leaves in a stack) usually affects the reflectance curve as we can find several records throughout the literature (e.g., [11,16,26,39]). In the study of Blackburn [16], the leaf stacking was used in purpose to increase the range of LAI and chlorophyll content. Since the majority of studies including contact probe measurements are using only one measurement setup (either a single leaf or a leaf stack) we focused on comparison of both setups on two temperate broadleaved tree species with a similar dorsiventral leaf structure growing in the same field conditions. Moreover, there are several gaps in knowledge on dynamics of different spectral indices during the vegetation season when measuring the leaf reflectance with the contact probe. Considering that a reflectance in particular spectral ranges (VIS, NIR and SWIR) changes with a different extent during the vegetation season, we assume that this effect on values of spectral indices is not negligible. Therefore, filling this gap in knowledge will contribute to interpretation of spectral data acquired by the contact probe at a leaf level in purpose to assess temporal changes in the physiological status of vegetation or leaf phenology.
Our aim was to compare the leaf spectral reflectance measured by the contact probe on a single leaf and a stack of five leaves of two broadleaved tree species in VIS, NIR and SWIR spectral ranges during the vegetative season and to test if and how selected vegetation indices differ under these two measurement setups. We asked following research questions:
(1)
How does the measurement setup (a single leaf or a leaf stack) influence a reflectance curve in selected spectral ranges (VIS, NIR and SWIR)?
(2)
Does the difference between the reflectances measured on a single leaf and a leaf stack in the above selected spectral ranges differ during the season?
(3)
Are the water- and pigment-related indices affected by the contact probe measurement setup (a single leaf or a leaf stack) during the season?

2. Materials and Methods

2.1. Study Site

The study was conducted at non-reclaimed post mining sites in the Sokolov coal mining district (50°14′21″ N, 12°39′24″ E) in the western part of the Czech Republic, Central Europe. The study area was located at an altitude of 500–600 m a.s.l., with mean annual precipitation of 650 mm and mean annual temperature of 6.8 °C [40]. The study site was covered by spontaneous woody vegetation dominated by European aspen (Populus tremula L.) and goat willow (Salix caprea L.). Due to the higher trichome density on the abaxial leaf side and lower specific leaf area, S. caprea leaves have a better xeromorphic adaptation than P. tremula leaves. However, both leaves exhibit a similar dorsiventral structure, i.e., mesophyll is differentiated into palisade and spongy parenchyma. These two species of broadleaved deciduous trees were sampled during the vegetative season from May to October 2014; particular sampling dates were: 26 May, 17 June, 8 July, 12 August, 9 September and 1 October. The trees were approximately 20–25 years old, and 5–7 m tall. Fully developed leaves were taken from the crown transition zone between sun-exposed and shaded parts of a tree crown. Fresh leaves from 10 trees of each species were clipped using telescopic pruning poles, put into the zip lock plastic bags, cooled and immediately transported to the laboratory. The spectra were measured within 24 h after the collecting. The leaf samples for assessment of chlorophyll and water contents were put into plastic vials, cooled and further processed as described below.

2.2. Leaf Spectra Measurements

Spectral measurements were conducted using an ASD FieldSpec 4 Wide-Res portable spectroradiometer (ASD Inc., Boulder, CO, USA). Each spectral measurement was preceded by a dark current measurement and a white reference measurement using a 99% Spectralon panel. The radiance spectra of the leaf reflectance between 350 and 2500 nm were normalized against white reference to produce relative reflectance spectra.
The leaves were placed on a spectrally black surface to minimize the background spectral noise or radiation transmitted through them. A fiber optic contact probe (ASD Plant Probe; ASD Inc., USA) was placed on a leaf surface, which was only illuminated by a constant light source inside the contact probe. The scan average on the spectroradiometer was set to 15 to avoid overheating the scanned foliage [9]. Five leaves per tree were used for spectral analyses. First, spectral reflectance of each of those five leaves was measured separately. Then, the five leaves were arranged into a stack by their adaxial surfaces facing to the top and the reflectance was measured again. Five independent measurements were taken on different parts of each sample. For each leaf or a leaf stack a median reflectance value was calculated and then the curve was used for further data analysis.

2.3. Assessment of Biophysical Leaf Traits

After reflectance measurements, fresh leaves were scanned and the leaf area determined by ImageJ freeware. Then the leaves were dried and weighed for assessment of LMA (leaf mass per area unit). The leaf water content was determined as the percentage of water in the fresh leaves (the fraction of weight decrease after drying at 80 °C for 48 h until constant weight) and EWT (equivalent water thickness) was calculated according to [41]. The total chlorophyll content was determined in dimethylformamide extract using a spectrophotometric assessment based on the equations reported in [42].

2.4. Spectral Ranges and Vegetation Indices

According to Cavender-Bares [10] we divided full leaf spectra to four spectral ranges: VIS (400–700 nm), NIR (740–1140 nm), SWIR1 (1500–1800 nm) and SWIR2 (2000–2400 nm). In purpose of quantification the effect of the MS in those spectral ranges, the difference between the reflectances measured on a leaf stack of 5 leaves and a single leaf (i.e., 1 leaf), (R5L subtraction R1L = ∆R5L–1L), was calculated for each wavelength. Thus, the value of ∆R5L–1L expresses the additional variability in leaf spectra introduced by increasing LAI, LMA and the internal leaf surface area as a potential scattering interface. Stacking of the leaves in a laboratory can be regarded as approximation of the contact measuring setup to the measurement at the canopy level.
Vegetation indices were computed in the statistic program R-3.2.3 (the open source version of Statistica) from the reflectance measured in both measurement setups (a single leaf, a leaf stack). In total, 74 pigment-related indices and 7 water-related indices were determined. For all tested indices and their references see Table 1. For demonstration of temporal dynamics of the selected indices, their values were transformed into standardized z-scores to enable their comparability and independence of their physical dimensions. Z-scores generally express how far from the mean the particular value is in terms of the standard deviations.

2.5. Statistical Analyses

All statistical analyses were performed with NCSS 9 software (NCSS, LCC, Kaysville, UT, USA). Descriptive statistics for the reflectance of a single leaf and a stack of leaves were calculated separately for VIS, NIR, SWIR1 and SWIR2 spectral ranges. The effect of the sampling date on the value of difference ∆R5L–1L and vegetation indices was tested by one-way ANOVA and Tukey–Kramer (in case of a normal distribution) or Kruskal–Wallis (if not normal distribution) tests. All differences were classified as significant if p < 0.05.

3. Results and Discussion

Reflectance curves acquired by a contact probe measurement on a single leaf and a leaf stack and the difference between them (∆R5L–1L) averaged from all measurements during the season are presented for P. tremula and S. caprea in Figure 1 and Figure 2, respectively. The results are very close to each other for both species, which is in accordance with the following facts: (1) The leaves of both tree species had a similar anatomical structure of a typically dorsiventral leaf, which also confirmed their similar values of LMA (Figure 3A,B); (2) Both tree species grew in close vicinity at one site, thus having the same environmental conditions. This is in accordance with the seasonal pattern of the EWT (Figure 3C,D), which was the same for both species; depending probably on precipitation; (3) Both species share the life strategy of pioneer species and similar phenology that is obvious from the chlorophyll content and LMA seasonal dynamics (Figure 3A,B and Figure 4A,B).

3.1. Effect of the Measurement Setup on the Reflectance Curve in Selected Spectral Ranges

In order to capture the effect of the MS on reflectance curves, we designed for the purpose of the present study a metric corresponding to the difference between the reflectances measured on a leaf stack of 5 leaves and a single leaf, i.e., reflectance difference (∆R5L–1L). The highest mean value of the reflectance difference ∆R5L–1L for both species was up to 30% in NIR (Figure 5B and Figure 6B) similarly as in the study by Blackburn [16]. We suppose that a higher ∆R5L–1L was caused by an increased content of structure components in stacks of leaves (e.g., increase in LAI and/or LMA) as reported in other studies [26,82,83]. Relatively small reflectance difference ∆R5L–1L in SWIR 1 (up to 10%) and SWIR2 (up to 5%) (Figure 5C,D and Figure 6C,D) could be explained by water-related absorption depending on EWT [34], which increases with each additional leaf layer [11]. The reflectance in VIS was up to 15% in both measurement setups (Figure 5A and Figure 6A). A similar trend was observed in a stack of 12 leaves in the study by [39]. The value of ∆R5L–1L in VIS was close to zero what confirmed that structure and function of a planar leaf is very well adapted to absorb light in the VIS range [84], regardless if the light comes directly down on a leaf or is transmitted from an upper leaf layer. The radiation reflected from the adaxial leaf side from leaves in a stack is again reflected from the spongy parenchyma of upper leaves, which is helpful for achievement of maximum utilization of light in a canopy [85].
The descriptive statistics for the reflectance difference ∆∆R5L-1L in selected spectral ranges (VIS, NIR, SWIR1 and SWIR 2), throughout the growing season from all samples, are presented in Table 2. Both P. tremula and S. caprea kept the similar values of ∆R5L–1L in the identical spectral ranges what can be explained by the same growing conditions and living strategies of both tree species. The similarity between optical properties of P. tremula and S. caprea is also based on the field-observed biophysical leaf traits—total chlorophyll content, leaf mass per area and equivalent water thickness (Figure 3 and Figure 4), which were almost identical in both species. Moreover, the temporal trends of all studied leaf traits followed the same pattern for both species.
For both tree species, the maximum and the minimum values of ∆R5L–1L were detected in NIR and VIS, respectively. The variance of ∆R5L–1L in NIR was up to two orders higher than in VIS. Measuring the leaf reflectance under a different measurement setup (a single leaf or more leaves in a stack) usually affects the reflectance curve [26]. The reflectance of a leaf stack in NIR was higher about 20–30% than the reflectance of a single leaf as it was shown in the study by [11] when the stack was composed of four or six planar leaves (species not indicated). We confirmed the same effect of the measurement setup in our study using one and five leaves of P. tremula and S. caprea.

3.2. Seasonal Dynamics of the Difference between Reflectances of a Leaf Stack and a Single Leaf

The reflectance difference ∆R5L–1L during the season in four selected spectral ranges (VIS, NIR, SWIR1 and SWIR2), is displayed in Figure 5A–D and Figure 6A–D for P. tremula and S. caprea, respectively. Differences in values of ∆R5L–1L among measurements conducted in different months during the vegetative season were significant, but not in all cases. Although values of ∆R5L–1L in VIS were lower than 1% in both species, the highest difference between values of ∆R5L–1L was detected between 500 and 600 nm in October for both species (Figure 5A and Figure 6A). We suppose that it was most likely due to the progress of the leaf senescence accompanied by the chlorophyll decay shown for both species (Figure 4). In both species, the effect of the measurement setup was more pronounced at the end of the season, comparing to the beginning of the season, showing the highest values of ∆R5L–1L in VIS in October. The ∆R5L–1L values in NIR were considerably higher than the values in other spectral ranges (VIS, SWIR1 and SWIR2) during the whole season. This can be explained that the reflectance in NIR depends on sample structural parameters (LMA, LAI) and leaf sample thickness, which mainly differed between the two used measurement setups, a single leaf and a stack of leaves. During the season, ∆R5L–1L in NIR declined in both P. tremula (Figure 5B) and S. caprea (Figure 6B). On one hand, we explain the decrease of ∆R5L–1L value by declining NIR reflectance in both a single leaf and a stack of leaves (data not shown) in the consequence of leaf aging, similarly as was observed by [4] in V. vinifera and by [8] in Alnus incana L. Moench and Betula pendula L. On the other hand, we expected that the accumulation of the cell wall structural compounds documented by the trend of LMA increase (Figure 3A,B) would lead to a higher NIR reflectance. However, we can attribute the decrease of both, the single leaf reflectance and also a value of ∆R5L–1L during the season, not only to absorption of NIR radiation by the structural compound themselves, but also to changing internal mesophyll structure (e.g., cell wall thickening). Modulating the size of the internal leaf surface may lead to a lower light scattering inside the leaf and, thus, a lower NIR reflectance.
The decrease of ∆R5L–1L value in SWIR1 (Figure 5C and Figure 6C) and SWIR2 (Figure 5D and Figure 6D) kept the same pattern during the season for both species. From May to July the decrease in ∆R5L–1L was gradual, however, from August to October the change in ∆R5L–1L was insignificant and rather inconsistent. The seasonal changes in ∆R5L–1L in SWIR could be partly explained by accumulation of structural compounds (LMA) similarly as in NIR, however, the effect was weaker. We attribute the fluctuation of value of ∆R5L–1L in SWIR between August and October to irregular changes of EWT in both species (Figure 3C,D).

3.3. The Effect of the Measurement Setup on Vegetation Indices during the Season

Subsequent paragraphs discuss differences between the water- and pigment-related indices computed from the reflectance acquired by both contact probe measurement setups (a single leaf, a leaf stack) during the vegetative season. Significant differences between indices calculated from both measurement setups during the season are expressed in Table 1. When the difference between vegetation indices obtained from two contact probe measurement setups was calculated, eight indices (BGI, CRI 550, CRI 700, MCARI, McM94, PRIm1, TCARI, and VI[700]) out of 85 indices tested were not affected by the measurement setup throughout the growing season in both tree species. All eight of these indices were derived from the reflectance in wavelengths up to 700 nm, i.e., in VIS. Five (BGI, MCARI, McM94, TCARI and VI[700]) out of those vegetation indices classified as unaffected by the measurement setup were chlorophyll-related and followed the seasonal course of chlorophyll leaf content in both species (see example of BGI in Figure 4 and Figure 7). The temporal variation of vegetation indices corresponded better with the leaf chlorophyll content in P. tremula in comparison to S. caprea. This difference between species is also documented by higher R2 values for linear regression between chlorophyll content and indices (Table 1) for P. tremula. In these five MS-unaffected chlorophyll indices, the R2 values for linear regression with the chlorophyll content were always higher for indices computed from a reflectance of a single leaf (Table 1).
On the other hand, indices including wavelengths from both VIS and NIR spectral regions (e.g., DD, Macc01, mND705, MNDVIRe, RMCARI, and RMACARIROSAVI) always showed a significant dependence on the measurement setup (a single leaf, a leaf stack). These indices also correlated significantly with the chlorophyll content (Table 1) and exhibited the good fit to a chlorophyll seasonal dynamics (see example of Datt in Figure 4). Here, the effect of the contact probe measurement setup was obvious (Figure 4A,B): stacking the five leaves significantly shifted the absolute values of the index comparing to measurement on a single leaf. We assume that the fact, if a particular index is affected by the contact probe measurement setup, is determined by the extent of contribution of lower leaves in the stack to resulting sample reflectance. Regarding this, Blackburn [16] showed that the lower leaves in the stack have a minor or negligible effect on overall reflectance, particularly in VIS. Therefore, the indices unaffected by the contact probe measurement setup were rather determined by the chlorophyll content in the uppermost leaf.
The z-score standardization of vegetation indices enabled to compare their seasonal course with the chlorophyll content dynamics irrespective of their different absolute values (Figure 7). However, we could still recognize different behaviour of MS-unaffected and MS-affected indices. The effect of the contact probe measurement setup on the latter ones was detectable particularly in the beginning of the season (May and June) and decreased towards the end of the season (Figure 7C,D and the p-values slightly increased since August, Table 1). In the case of processing vegetation indices and spectral data acquired under different contact probe measurement setups (regarding leaf stacking), we recommend to take into account this dependence on measurement setup particularly in the beginning of the season, even if the data standardization was applied.
The results presented above imply that vegetation indices, which involve wavelengths from VIS only, were almost unaffected by a contact probe measurement setup of a single leaf or a leaf stack, in comparison to indices operating with reflectance wavelengths from more than one spectral range. The dependence of vegetation indices on measurement setup is given by a larger, more contrasting difference in ∆R5L–1L between spectral ranges used for index calculation (Figure 1). Vegetation indices including wavelengths from two spectral ranges with different value of ∆R5L–1L were more dependent on a measurement setup in the beginning of the vegetation season than at the end of it (Figure 7C,D). The dynamic change of ∆R5L–1L during the season may affect the correlation strength of these vegetation indices to leaf biophysical traits used for the estimation of a tree physiological status if the reflectance is measured at a leaf level. However, many of MS-affected indices, such as those presented in (Figure 7C,D) better correlated with chlorophyll content in both species in comparison with MS-unaffected indices (Table 1). Therefore, in terms of monitoring leaf biophysical traits using the spectroscopic approach, NIR reflectance should not be excluded and these MS-affected indices should be used with the above described caution.
Currently, there are emerging studies on spectra modeling or reverse radiative transfer modeling that combine leaf optical properties and biophysical traits from independent experiments or open access spectral libraries: LeuvenV and LOPEX datasets [86]; ANGERS, LOPEX and PANAMA datasets [87]; or various independent datasets (see Table 1 in [88]). These studies also test and implement vegetation indices derived from spectra provided in these databases [87]. We are convinced that it is always important to use the data acquired by the most comparable measurement setups. Concerning this, in our previous study of Potůčková [89], we emphasized that in case of processing spectra acquired by two different devices, a linear model transforming spectra between those devices is necessary. On the example of several vegetation indices we have demonstrated, that in most cases the measurement setup should be taken into account, even if the identical instruments (the spectroradiometer and the contact probe) and data standardization are used.

3.4. Leaf Stack as an Intermediate Step towards the Canopy Reflectance

Measuring reflectance on more leaves together in a stack in a present study had a similar effect on the leaf reflectance as changing a structural parameter N in the radiative transfer model PROSPECT [17]. Simulation of the leaf reflectance in PROSPECT includes structural parameter N, corresponding to a number of elementary homogeneous layers in a leaf [87]. Increasing of N parameter results in the increase of the reflectance values systematically in all spectral ranges [17]. Nevertheless, the effect of changes in leaf structural parameters on quantitative values of the reflectance in the VIS range was negligible in comparison to changes induced by pigment contents [83].
In some cases, we can consider the reflectance of a leaf stack as an intermediate step towards the canopy reflectance. Regarding that the change between the reflectance curve of a single leaf and a stack of leaves in our study across all studied spectral ranges (Figure 1 and Figure 2) was very similar to the effect of the increasing LAI in the canopy reflectance simulated by the radiative transfer model PROSAIL [90]. Although we suggest that stacking the leaves can partly approximate the contact probe reflectance measurement to the canopy level, we are aware that the reflectance derived from a leaf stack is not equal to the canopy-scale reflectance. By piling several layers of leaves we technically increased chlorophyll content, leaf area index and leaf mass per area unit (LAI and LMA). In addition, the effect of the internal leaf structure such as both volume and surface of intercellular spaces is enhanced in a leaf stack. These effects of a leaf stacking altogether can influence the correlation between leaf biophysical traits and leaf optical properties including vegetation indices, which we demonstrated with systematically lower R2 values for vegetation indices acquired from the reflectance of the leaf stack (Table 1). Moreover, it is important to mention that the contribution of leaf traits to the resulting leaf stack reflectance may depend on a leaf position within the stack [16]. On the other hand, the real canopy-scale reflectance is affected by additional factors, such as leaf clumping, leaf angle distribution, presence of non-photosynthetic structures (branches and twigs) and soil/understory reflectance background as summarized, e.g., by Homolová et al. [32]. In case of airborne remote sensing, the effect of atmosphere should be considered and atmospheric corrections should be made [91]. Currently, the appropriate upscaling of the reflectance from the single leaf level to the canopy level is usually performed by the coupling of radiative transfer models at the leaf (e.g., PROSPECT) and canopy levels (e.g., SAIL [92] or DART [19]).

5. Conclusions

The major influence of the contact probe measurement setup (one leaf or a stack of five leaves) on a leaf spectral curve was detected in NIR, the structure-related spectral range and a weaker effect was detected in SWIR1 and SWIR2, water-related spectral ranges. In VIS, the pigment-related spectral region, the effect of measurement setup was almost absent. Our results confirmed that in NIR, the reflectance of a leaf stack is about 25% higher than the reflectance of a single leaf. In a consequence of good adaptation of a planar leaf to absorb VIS radiation incident to its adaxial side (including the light transmitted from upper leaf layers), there was no influence connected with measurement setup in VIS. We concluded that the main effect on difference between the two measurement setups in NIR was induced by the change in sample structural parameters originating from additional leaves in the stack. On the contrary, in VIS, the chlorophyll content of leaves placed lower in the stack contributed less to the reflectance of the sample not affecting it significantly.
The considerable influence of a measurement setup was detected at the beginning of the season, according to the highest values of ∆R5L–1L in May in all selected spectral ranges except VIS. The maximum values of ∆R5L–1L in VIS observed in October were most likely due to the leaf senescence and chlorophyll content decrease in both species. The values of ∆R5L–1L in NIR and SWIR were higher at the beginning of the season and diminished towards the end of the season. We attribute the decrease of ∆R5L–1L values during the season to the accumulation of the cell wall structural compounds documented by the trend of LMA increase. We suppose that the cell wall thickening modulates the size of the internal leaf surface and may lead to lower scattering inside the leaf and, thus, lowering the NIR reflectance. We assume that the seasonal changes in values of ∆R5L–1L in SWIR could be partly explained by the fluctuation of leaf equivalent water thickness and the accumulation of structural compounds, similar to in NIR, however, the effect was weaker.
The strong effect of stacking up the leaves while measuring reflectance by a contact probe should be regarded in studies where the reflectance is measured on a stack of leaves. We clearly showed that vegetation indices involving wavelengths up to 700 nm were independent of the measurement setup during the vegetation season. On the contrary, indices combining wavelengths from both VIS and NIR were affected by the measurement setup throughout the season. We assume that the sensitivity of a particular index to the contact probe measurement setup is given by the contribution of lower leaves in the stack to resulting sample reflectance. This effect of leaf stacking contributed to weakening the correlation between the leaf chlorophyll content and vegetation indices what then in practice may complicate chlorophyll retrieval from contact probe measurements from a leaf stack. It also implies that some chlorophyll-related indices unaffected by the contact probe measurement setup could be used in meta-studies even when combining various setups of the contact probe measurements at the leaf level. On the other hand, the effect of the contact probe measurement setup on chlorophyll indices combining wavelengths from both VIS and NIR (indices which are MS-affected) was present even after the z-score standardization, particularly in the beginning of the season. Therefore, we recommend taking into account this dependence on measurement setup, regardless the data standardization. Many MS-affected indices were better correlated with chlorophyll content in both species in comparison with MS-unaffected indices. Therefore, in terms of monitoring leaf chlorophyll content using the spectroscopic approach, these MS-affected indices should be used with the above-described caution.
Even if the vegetation indices are used for assessment of plant physiological status in various times of a vegetative season, it is still essential to take into consideration their possible changes induced by the particular contact probe measurement setup regarding the leaf stacking.

Acknowledgments

The experimental part of the study was performed in the framework of the project of the Czech Science Foundation (GAČR P504/12/1288); the data processing and manuscript writing under the project of the Ministry of Education, Youth and Sports of the Czech Republic (NPUI LO1417). We thank Jan Frouz for support in the field sampling and the GAČR project leadership. We thank Drahomíra Bartáková, Monika Kovářová, and Miroslav Barták for the help in the leaf sampling. We thank two anonymous reviewers for giving us constructive criticism and valuable suggestions how to improve the manuscript.

Author Contributions

Eva Neuwirthová did the literature review, participated in the field sampling, conducted measurements with ASD CP, processed and analyzed all dataset of optical properties and vegetation indices, performed the statistic tests and contributed to overall manuscript writing. Zuzana Lhotáková contributed to the design of the study, participated on measurements of optical properties of plant material and participated in manuscript writing and interpretation of the results. Jana Albrechtová contributed to the design of the study, interpretation of the results and final manuscript writing and editing. All authors contributed to writing and editing of the paper and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Campbell, P.K.E.; Rock, B.N.; Martin, M.E.; Neefus, C.D.; Irons, J.R.; Middleton, E.M.; Albrechtova, J. Detection of initial damage in Norway spruce canopies using hyperspectral airborne data. Int. J. Remote Sens. 2004, 25, 5557–5584. [Google Scholar] [CrossRef]
  2. Mišurec, J.; Kopačková, V.; Lhotáková, Z.; Hanuš, J.; Weyermann, J.; Entcheva-Campbell, P.; Albrechtová, J. Utilization of hyperspectral image optical indices to assess the Norway spruce forest health status. J. Appl. Remote Sens. 2012, 6, 063545. [Google Scholar]
  3. Kopačková, V.; Lhotáková, Z.; Oulehle, F.; Albrechtová, J. Assessing forest health via linking the geochemical properties of a soil profile with the biochemical parameters of vegetation. Int. J. Environ. Sci. Technol. 2015, 12, 1987–2002. [Google Scholar] [CrossRef]
  4. Rapaport, T.; Hochberg, U.; Rachmilevitch, S.; Karnieli, A. The Effect of Differential Growth Rates across Plants on Spectral Predictions of Physiological Parameters. PLoS ONE 2014, 9, e88930. [Google Scholar] [CrossRef] [PubMed]
  5. Croft, H.; Chen, J.M.; Zhang, Y. The applicability of empirical vegetation indices for determining leaf chlorophyll content over different leaf and canopy structures. Ecol. Complex. 2014, 17, 119–130. [Google Scholar] [CrossRef]
  6. Gates, D.M.; Keegan, H.J.; Schleter, J.C.; Weidner, V.R. Spectral properties of plants. Appl. Opt. 1965, 4, 11–20. [Google Scholar] [CrossRef]
  7. Daughtry, C.S.T.; Walthall, C.L.; Kim, M.S.; De Colstoun, E.B.; McMurtrey, J.E. Estimating corn leaf chlorophyll concentration from leaf and canopy reflectance. Remote Sens. Environ. 2000, 74, 229–239. [Google Scholar] [CrossRef]
  8. Mottus, M.; Sulev, M.; Hallik, L. Seasonal Course of the Spectral Properties of Alder and Birch Leaves. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2014, 7, 2496–2505. [Google Scholar] [CrossRef]
  9. Eitel, J.U.H.; Gessler, P.E.; Smith, A.M.S.; Robberecht, R. Suitability of existing and novel spectral indices to remotely detect water stress in Populus spp. For. Ecol. Manag. 2006, 229, 170–182. [Google Scholar] [CrossRef]
  10. Cavender-Bares, J.; Meireles, J.E.; Couture, J.J.; Kaproth, M.A.; Kingdon, C.C.; Singh, A.; Serbin, S.P.; Center, A.; Zuniga, E.; Pilz, G.; et al. Associations of Leaf Spectra with Genetic and Phylogenetic Variation in Oaks: Prospects for Remote Detection of Biodiversity. Remote Sens. 2016, 8, 221. [Google Scholar] [CrossRef]
  11. Gao, B. NDWI—A normalized difference water index for remote sensing of vegetation liquid water from space. Remote Sens. Environ. 1996, 58, 257–266. [Google Scholar] [CrossRef]
  12. Lhotakova, Z.; Brodsky, L.; Kupkova, L.; Kopackova, V.; Potuckova, M.; Misurec, J.; Klement, A.; Kovarova, M.; Albrechtova, J. Detection of multiple stresses in Scots pine growing at post-mining sites using visible to near-infrared spectroscopy. Environ. Sci. Process. Impacts 2013, 15, 2004–2015. [Google Scholar]
  13. Einzmann, K.; Ng, W.-T.; Immitzer, M.; Pinnel, N.; Atzberger, C. Method Analysis for Collecting and Processing in-situ Hyperspectral Needle Reflectance Data for Monitoring Norway Spruce; Methodenanalyse zur Erfassung und Prozessierung hyperspektraler in-situ Nadelreflexionsdaten zum Monitoring von Fichten. Photogramm. Fernerkund. Geoinform. 2014, 2014, 423–434. [Google Scholar] [CrossRef]
  14. Lu, X.; Peng, H. Predicting Cherry Leaf Chlorophyll Concentrations Based on Foliar Reflectance Spectra Variables. J. Indian Soc. Remote Sens. 2015, 43, 109–120. [Google Scholar] [CrossRef]
  15. 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]
  16. Blackburn, G.A. Relationships between Spectral Reflectance and Pigment Concentrations in Stacks of Deciduous Broadleaves. Remote Sens. Environ. 1999, 70, 224–237. [Google Scholar] [CrossRef]
  17. Jacquemoud, S.; Baret, F. PROSPECT: A model of leaf optical properties spectra. Remote Sens. Environ. 1990, 34, 75–91. [Google Scholar] [CrossRef]
  18. Dian, Y.; Le, Y.; Fang, S.; Xu, Y.; Yao, C.; Liu, G. Influence of Spectral Bandwidth and Position on Chlorophyll Content Retrieval at Leaf and Canopy Levels. J. Indian Soc. Remote Sens. 2016, 44, 583–593. [Google Scholar] [CrossRef]
  19. Gastellu-Etchegorry, J.P.; Martin, E.; Gascon, F. DART: A 3D model for simulating satellite images and studying surface radiation budget. Int. J. Remote Sens. 2004, 25, 73–96. [Google Scholar] [CrossRef]
  20. Hernández-Clemente, R.; Navarro-Cerrillo, R.M.; Zarco-Tejada, P.J. Carotenoid content estimation in a heterogeneous conifer forest using narrow-band indices and PROSPECT + DART simulations. Remote Sens. Environ. 2012, 127, 298–315. [Google Scholar] [CrossRef]
  21. Huang, N.; Niu, Z.; Zhan, Y.; Xu, S.; Tappert, M.C.; Wu, C.; Huang, W.; Gao, S.; Hou, X.; Cai, D. Relationships between soil respiration and photosynthesis-related spectral vegetation indices in two cropland ecosystems. Agric. For. Meteorol. 2012, 160, 80–89. [Google Scholar] [CrossRef]
  22. Yi, Q.; Jiapaer, G.; Chen, J.; Bao, A.; Wang, F. Different units of measurement of carotenoids estimation in cotton using hyperspectral indices and partial least square regression. ISPRS J. Photogramm. Remote Sens. 2014, 91, 72–84. [Google Scholar] [CrossRef]
  23. Liu, N.; Lin, Z.-F.; Van Devender, A.; Lin, G.-Z.; Peng, C.-L.; Pan, X.-P.; Chen, S.-W.; Gu, Q. Spectral reflectance indices and pigment functions during leaf ontogenesis in six subtropical landscape plants. Plant Growth Regul. 2009, 58, 73–84. [Google Scholar] [CrossRef]
  24. Blackburn, G.A. Quantifying chlorophylls and caroteniods at leaf and canopy scales: An evaluation of some hyperspectral approaches. Remote Sens. Environ. 1998, 66, 273–285. [Google Scholar] [CrossRef]
  25. Zarco-Tejada, P.J.; Berjón, A.; López-Lozano, R.; Miller, J.R.; Martín, P.; Cachorro, V.; González, M.R.; de Frutos, A. Assessing vineyard condition with hyperspectral indices: Leaf and canopy reflectance simulation in a row-structured discontinuous canopy. Remote Sens. Environ. 2005, 99, 271–287. [Google Scholar] [CrossRef]
  26. Kokaly, R.F.; Asner, G.P.; Ollinger, S.V.; Martin, M.E.; Wessman, C.A. Characterizing canopy biochemistry from imaging spectroscopy and its application to ecosystem studies. Remote Sens. Environ. 2009, 113, S78–S91. [Google Scholar] [CrossRef]
  27. Ustin, S.L.; Gitelson, A.A.; Jacquemoud, S.; Schaepman, M.; Asner, G.P.; Gamon, J.A.; Zarco-Tejada, P. Retrieval of foliar information about plant pigment systems from high resolution spectroscopy. Remote Sens. Environ. 2009, 113, S67–S77. [Google Scholar] [CrossRef]
  28. Castro, K.L.; Sanchez-Azofeifa, G.A. Changes in spectral properties, chlorophyll content and internal mesophyll structure of senescing Populus balsamifera and Populus tremuloides leaves. Sensors 2008, 8, 51–69. [Google Scholar] [CrossRef] [PubMed]
  29. Gitelson, A.A.; Merzlyak, M.N.; Lichtenthaler, H.K. Detection of Red Edge Position and Chlorophyll Content by Reflectance Measurements near 700 nm. J. Plant Physiol. 1996, 148, 501–508. [Google Scholar] [CrossRef]
  30. Rondeaux, G.; Steven, M.; Baret, F. Optimization of soil-adjusted vegetation indices. Remote Sens. Environ. 1996, 55, 95–107. [Google Scholar] [CrossRef]
  31. Main, R.; Cho, M.A.; Mathieu, R.; O’Kennedy, M.M.; Ramoelo, A.; Koch, S. An investigation into robust spectral indices for leaf chlorophyll estimation. ISPRS J. Photogramm. Remote Sens. 2011, 66, 751–761. [Google Scholar] [CrossRef]
  32. Homolová, L.; Malenovský, Z.; Clevers, J.G.P.W.; García-Santos, G.; Schaepman, M.E. Review of optical-based remote sensing for plant trait mapping. Ecol. Complex. 2013, 15, 1–16. [Google Scholar] [CrossRef]
  33. Le Maire, G.; François, C.; Soudani, K.; Berveiller, D.; Pontailler, J.-Y.; Bréda, N.; Genet, H.; Davi, H.; Dufrêne, E. Calibration and validation of hyperspectral indices for the estimation of broadleaved forest leaf chlorophyll content, leaf mass per area, leaf area index and leaf canopy biomass. Remote Sens. Environ. 2008, 112, 3846–3864. [Google Scholar] [CrossRef]
  34. González-Fernández, A.B.; Rodríguez-Pérez, J.R.; Marcelo, V.; Valenciano, J.B. Using field spectrometry and a plant probe accessory to determine leaf water content in commercial vineyards. Agric. Water Manag. 2015, 156, 43–50. [Google Scholar] [CrossRef]
  35. Baránková, B.; Lazár, D.; Nauš, J. Analysis of the effect of chloroplast arrangement on optical properties of green tobacco leaves. Remote Sens. Environ. 2016, 174, 181–196. [Google Scholar] [CrossRef]
  36. Nakaji, T.; Kosugi, Y.; Takanashi, S.; Niiyama, K.; Noguchi, S.; Tani, M.; Oguma, H.; Nik, A.R.; Kassim, A.R. Estimation of light-use efficiency through a combinational use of the photochemical reflectance index and vapor pressure deficit in an evergreen tropical rainforest at Pasoh, Peninsular Malaysia. Remote Sens. Environ. 2014, 150, 82–92. [Google Scholar] [CrossRef]
  37. Dillen, S.Y.; de Beeck, M.O.; Hufkens, K.; Buonanduci, M.; Phillips, N.G. Seasonal patterns of foliar reflectance in relation to photosynthetic capacity and color index in two co-occurring tree species, Quercus rubra and Betula papyrifera. Agric. For. Meteorol. 2012, 160, 60–68. [Google Scholar] [CrossRef]
  38. Fréchette, E.; Wong, C.Y.S.; Junker, L.V.; Chang, C.Y.-Y.; Ensminger, I. Zeaxanthin-independent energy quenching and alternative electron sinks cause a decoupling of the relationship between the photochemical reflectance index (PRI) and photosynthesis in an evergreen conifer during spring. J. Exp. Bot. 2015, 66, 7309–7323. [Google Scholar] [CrossRef] [PubMed]
  39. Bartlett, M.K.; Ollinger, S.V.; Hollinger, D.Y.; Wicklein, H.F.; Richardson, A.D. Canopy-scale relationships between foliar nitrogen and albedo are not observed in leaf reflectance and transmittance within temperate deciduous tree species. Botany 2011, 89, 491–497. [Google Scholar] [CrossRef]
  40. Frouz, J.; Keplin, B.; Pižl, V.; Tajovský, K.; Starý, J.; Lukešová, A.; Nováková, A.; Balı́k, V.; Háněl, L.; Materna, J.; et al. Soil biota and upper soil layer development in two contrasting post-mining chronosequences. Ecol. Eng. 2001, 17, 275–284. [Google Scholar] [CrossRef]
  41. Mobasheri, M.R.; Fatemi, S.B. Leaf Equivalent Water Thickness assessment using reflectance at optimum wavelengths. Theor. Exp. Plant Physiol. 2013, 25, 196–202. [Google Scholar] [CrossRef]
  42. Wellburn, A.R. The Spectral Determination of Chlorophylls a and b, as well as Total Carotenoids, Using Various Solvents with Spectrophotometers of Different Resolution. J. Plant Physiol. 1994, 144, 307–313. [Google Scholar] [CrossRef]
  43. Liu, J.G.; Moore, J.M. Hue image RGB colour composition. A simple technique to suppress shadow and enhance spectral signature. Int. J. Remote Sens. 1990, 11, 1521–1530. [Google Scholar] [CrossRef]
  44. Carter, G.A. Ratios of leaf reflectances in narrow wavebands as indicators of plant stress. Int. J. Remote Sens. 1994, 15, 697–703. [Google Scholar] [CrossRef]
  45. Gitelson, A.A.; Gritz, Y.; Merzlyak, M.N. Relationships between leaf chlorophyll content and spectral reflectance and algorithms for non-destructive chlorophyll assessment in higher plant leaves. J. Plant Physiol. 2003, 160, 271–282. [Google Scholar] [CrossRef] [PubMed]
  46. Zarco-Tejada, P.J.; Miller, J.R.; Noland, T.L.; Mohammed, G.H.; Sampson, P.H. Scaling-up and model inversion methods with narrowband optical indices for chlorophyll content estimation in closed forest canopies with hyperspectral data. IEEE Trans. Geosci. Remote Sens. 2001, 39, 1491–1507. [Google Scholar] [CrossRef]
  47. Datt, B. A New Reflectance Index for Remote Sensing of Chlorophyll Content in Higher Plants: Tests using Eucalyptus Leaves. J. Plant Physiol. 1999, 154, 30–36. [Google Scholar] [CrossRef]
  48. Le Maire, G.; François, C.; Dufrêne, E. Towards universal broad leaf chlorophyll indices using PROSPECT simulated database and hyperspectral reflectance measurements. Remote Sens. Environ. 2004, 89, 1–28. [Google Scholar] [CrossRef]
  49. Jordan, C.F. Derivation of Leaf-Area Index from Quality of Light on the Forest Floor. Ecology 1969, 50, 663–666. [Google Scholar] [CrossRef]
  50. Gitelson, A.; Merzlyak, M.N. Spectral Reflectance Changes Associated with Autumn Senescence of Aesculus hippocastanum L. and Acer platanoides L. Leaves. Spectral Features and Relation to Chlorophyll Estimation. J. Plant Physiol. 1994, 143, 286–292. [Google Scholar] [CrossRef]
  51. Smith, R.C.G.; Adams, J.; Stephens, D.J.; Hick, P.T.; Smith, R.C.G.; Adams, J.; Stephens, D.J.; Hick, P.T. Forecasting wheat yield in a Mediterranean-type environment from the NOAA satellite. Crop Pasture Sci. 1995, 46, 113–125. [Google Scholar] [CrossRef]
  52. Maccioni, A.; Agati, G.; Mazzinghi, P. New vegetation indices for remote measurement of chlorophylls based on leaf directional reflectance spectra. J. Photochem. Photobiol. B 2001, 61, 52–61. [Google Scholar] [CrossRef]
  53. Haboudane, D.; Miller, J.R.; Pattey, E.; Zarco-Tejada, P.J.; Strachan, I.B. Hyperspectral vegetation indices and novel algorithms for predicting green LAI of crop canopies: Modeling and validation in the context of precision agriculture. Remote Sens. Environ. 2004, 90, 337–352. [Google Scholar] [CrossRef]
  54. Wu, C.; Niu, Z.; Tang, Q.; Huang, W. Estimating chlorophyll content from hyperspectral vegetation indices: Modeling and validation. Agric. For. Meteorol. 2008, 148, 1230–1241. [Google Scholar] [CrossRef]
  55. McMurtrey, J.E.; Chappelle, E.W.; Kim, M.S.; Meisinger, J.J.; Corp, L.A. Distinguishing nitrogen fertilization levels in field corn (Zea mays L.) with actively induced fluorescence and passive reflectance measurements. Remote Sens. Environ. 1994, 47, 36–44. [Google Scholar] [CrossRef]
  56. Mutanga, O.; Skidmore, A.K. Narrow band vegetation indices overcome the saturation problem in biomass estimation. Int. J. Remote Sens. 2004, 25, 3999–4014. [Google Scholar] [CrossRef]
  57. Qi, J.; Chehbouni, A.; Huete, A.R.; Kerr, Y.H.; Sorooshian, S. A modified soil adjusted vegetation index. Remote Sens. Environ. 1994, 48, 119–126. [Google Scholar] [CrossRef]
  58. Hunt, E.R.; Rock, B.N. Detection of changes in leaf water content using Near- and Middle-Infrared reflectances. Remote Sens. Environ. 1989, 30, 43–54. [Google Scholar]
  59. Chen, J.M. Evaluation of Vegetation Indices and a Modified Simple Ratio for Boreal Applications. Can. J. Remote Sens. 1996, 22, 229–242. [Google Scholar] [CrossRef]
  60. Dash, J.; Curran, P.J. The MERIS terrestrial chlorophyll index. Int. J. Remote Sens. 2004, 25, 5403–5413. [Google Scholar] [CrossRef]
  61. Serrano, L.; Gamon, J.A.; Berry, J. Estimation of leaf area with an integrating sphere. Tree Physiol. 1997, 17, 571–576. [Google Scholar] [CrossRef] [PubMed]
  62. Rouse, J.W.; Haas, R.H.; Schell, J.A.; Deering, D.W.; Harlan, J.C. Monitoring the Vernal Advancement and Retrogradation (Greenwave Effect) of Natural Vegetation; NASA/GSFC Type III Final Report; NASA: Greenbelt, MD, USA, 1974; p. 371. [Google Scholar]
  63. Stimson, H.C.; Breshears, D.D.; Ustin, S.L.; Kefauver, S.C. Spectral sensing of foliar water conditions in two co-occurring conifer species: Pinus edulis and Juniperus monosperma. Remote Sens. Environ. 2005, 96, 108–118. [Google Scholar] [CrossRef]
  64. Chen, D.; Huang, J.; Jackson, T.J. Vegetation water content estimation for corn and soybeans using spectral indices derived from MODIS near- and short-wave infrared bands. Remote Sens. Environ. 2005, 98, 225–236. [Google Scholar] [CrossRef]
  65. Wang, L.; Qu, J.J. NMDI: A normalized multi-band drought index for monitoring soil and vegetation moisture with satellite remote sensing. Geophys. Res. Lett. 2007, 34, L20405. [Google Scholar] [CrossRef]
  66. Bravo, C.; Moshou, D.; Oberti, R.; West, J.; McCartney, A.; Bodria, L.; Ramon, H. Foliar Disease Detection in the Field Using Optical Sensor Fusion. J. Sci. Res. Dev. 2004, 6, 1–14. [Google Scholar]
  67. Gamon, J.A.; Peñuelas, 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]
  68. Hernández-Clemente, R.; Navarro-Cerrillo, R.M.; Suárez, L.; Morales, F.; Zarco-Tejada, P.J. Assessing structural effects on PRI for stress detection in conifer forests. Remote Sens. Environ. 2011, 115, 2360–2375. [Google Scholar] [CrossRef]
  69. Garrity, S.R.; Eitel, J.U.H.; Vierling, L.A. Disentangling the relationships between plant pigments and the photochemical reflectance index reveals a new approach for remote estimation of carotenoid content. Remote Sens. Environ. 2011, 115, 628–635. [Google Scholar] [CrossRef]
  70. Roujean, J.-L.; Breon, F.-M. Estimating PAR absorbed by vegetation from bidirectional reflectance measurements. Remote Sens. Environ. 1995, 51, 375–384. [Google Scholar] [CrossRef]
  71. Guyot, G.; Baret, F. Utilisation de la Haute Resolution Spectrale pour Suivre L’etat des Couverts Vegetaux. In Proceedings of the Spectral Signature of Objects in Remote Sensing, Aussois, France, 18–22 January 1988; p. 279. [Google Scholar]
  72. Curran, P.J.; Windham, W.R.; Gholz, H.L. Exploring the relationship between reflectance red edge and chlorophyll concentration in slash pine leaves. Tree Physiol. 1995, 15, 203–206. [Google Scholar] [CrossRef] [PubMed]
  73. Peñuelas, J.; Filella, I.; Lloret, P.; Muñoz, F.; Vilajeliu, M. Reflectance assessment of mite effects on apple trees. Int. J. Remote Sens. 1995, 16, 2727–2733. [Google Scholar] [CrossRef]
  74. Zarco-Tejada, P.J.; Miller, J.R. Land cover mapping at BOREAS using red edge spectral parameters from CASI imagery. J. Geophys. Res. 1999, 104, 27921–27933. [Google Scholar] [CrossRef]
  75. Gitelson, A.A.; Merzlyak, M.N. Remote estimation of chlorophyll content in higher plant leaves. Int. J. Remote Sens. 1997, 18, 2691–2697. [Google Scholar] [CrossRef]
  76. Zarco-Tejada, P.; Rueda, C.; Ustin, S. Water content estimation in vegetation with MODIS reflectance data and model inversion methods. Remote Sens. Environ. 2003, 85, 109–124. [Google Scholar] [CrossRef]
  77. Haboudane, D.; Miller, J.R.; Tremblay, N.; Zarco-Tejada, P.J.; Dextraze, L. Integrated narrow-band vegetation indices for prediction of crop chlorophyll content for application to precision agriculture. Remote Sens. Environ. 2002, 81, 416–426. [Google Scholar] [CrossRef]
  78. Broge, N.H.; Leblanc, E. Comparing prediction power and stability of broadband and hyperspectral vegetation indices for estimation of green leaf area index and canopy chlorophyll density. Remote Sens. Environ. 2000, 76, 156–172. [Google Scholar] [CrossRef]
  79. Gitelson, A.A.; Kaufman, Y.J.; Stark, R.; Rundquist, D. Novel algorithms for remote estimation of vegetation fraction. Remote Sens. Environ. 2002, 80, 76–87. [Google Scholar] [CrossRef]
  80. Vogelmann, J.E.; Rock, B.N.; Moss, D.M. Red edge spectral measurements from sugar maple leaves. Int. J. Remote Sens. 1993, 14, 1563–1575. [Google Scholar] [CrossRef]
  81. Peñuelas, J.; Filella, I.; Biel, C.; Serrano, L.; Savé, R. The reflectance at the 950–970 nm region as an indicator of plant water status. Int. J. Remote Sens. 1993, 14, 1887–1905. [Google Scholar] [CrossRef]
  82. Buschmann, C.; Lenk, S.; Lichtenthaler, H.K. Reflectance spectra and images of green leaves with different tissue structure and chlorophyll content. Isr. J. Plant Sci. 2012, 60, 49–64. [Google Scholar] [CrossRef]
  83. Slaton, M.R.; Hunt, E.R.; Smith, W.K. Estimating near-infrared leaf reflectance from leaf structural characteristics. Am. J. Bot. 2001, 88, 278–284. [Google Scholar] [CrossRef] [PubMed]
  84. Smith, W.K.; Vogelmann, T.C.; DeLucia, E.H.; Bell, D.T.; Shepherd, K.A. Leaf Form and Photosynthesis. BioScience 1997, 47, 785–793. [Google Scholar] [CrossRef]
  85. Cordón, G.B.; Lagorio, M.G. Optical properties of the adaxial and abaxial faces of leaves. Chlorophyll fluorescence, absorption and scattering coefficients. Photochem. Photobiol. Sci. 2007, 6, 873–882. [Google Scholar] [CrossRef] [PubMed]
  86. Stuckens, J.; Verstraeten, W.W.; Delalieux, S.; Swennen, R.; Coppin, P. A dorsiventral leaf radiative transfer model: Development, validation and improved model inversion techniques. Remote Sens. Environ. 2009, 113, 2560–2573. [Google Scholar] [CrossRef]
  87. Cheng, T.; Rivard, B.; Sánchez-Azofeifa, A.G.; Féret, J.-B.; Jacquemoud, S.; Ustin, S.L. Deriving leaf mass per area (LMA) from foliar reflectance across a variety of plant species using continuous wavelet analysis. ISPRS J. Photogramm. Remote Sens. 2014, 87, 28–38. [Google Scholar] [CrossRef]
  88. Féret, J.-B.; François, C.; Gitelson, A.; Asner, G.P.; Barry, K.M.; Panigada, C.; Richardson, A.D.; Jacquemoud, S. Optimizing spectral indices and chemometric analysis of leaf chemical properties using radiative transfer modeling. Remote Sens. Environ. 2011, 115, 2742–2750. [Google Scholar] [CrossRef]
  89. Potůčková, M.; Červená, L.; Kupková, L.; Lhotáková, Z.; Lukeš, P.; Hanuš, J.; Novotný, J.; Albrechtová, J. Comparison of Reflectance Measurements Acquired with a Contact Probe and an Integration Sphere: Implications for the Spectral Properties of Vegetation at a Leaf Level. Sensors 2016, 16, 1801. [Google Scholar] [CrossRef] [PubMed]
  90. Jacquemoud, S.; Verhoef, W.; Baret, F.; Bacour, C.; Zarco-Tejada, P.J.; Asner, G.P.; François, C.; Ustin, S.L. PROSPECT + SAIL models: A review of use for vegetation characterization. Remote Sens. Environ. 2009, 113, S56–S66. [Google Scholar] [CrossRef]
  91. Guanter, L.; Estellés, V.; Moreno, J. Spectral calibration and atmospheric correction of ultra-fine spectral and spatial resolution remote sensing data. Application to CASI-1500 data. Remote Sens. Environ. 2007, 109, 54–65. [Google Scholar] [CrossRef]
  92. Verhoef, W. Light scattering by leaf layers with application to canopy reflectance modeling: The SAIL model. Remote Sens. Environ. 1984, 16, 125–141. [Google Scholar] [CrossRef]
Figure 1. Averaged reflectance curves measured by a contact probe for a single leaf of P. tremula, a leaf stack and a difference (∆R5L–1L) between the reflectance measured on a leaf stack (5L) and a single leaf (1L). The mean reflectance (%) during the six months in P. tremula (n = 10 trees).
Figure 1. Averaged reflectance curves measured by a contact probe for a single leaf of P. tremula, a leaf stack and a difference (∆R5L–1L) between the reflectance measured on a leaf stack (5L) and a single leaf (1L). The mean reflectance (%) during the six months in P. tremula (n = 10 trees).
Sensors 17 01202 g001
Figure 2. Averaged reflectance curves measured by a contact probe for a single leaf of S. caprea, a leaf stack and a difference (∆R5L–1L) between the reflectance measured on a leaf stack (5L) and a single leaf (1L). The mean reflectance (%) during the six months in S. caprea (n = 10 trees).
Figure 2. Averaged reflectance curves measured by a contact probe for a single leaf of S. caprea, a leaf stack and a difference (∆R5L–1L) between the reflectance measured on a leaf stack (5L) and a single leaf (1L). The mean reflectance (%) during the six months in S. caprea (n = 10 trees).
Sensors 17 01202 g002
Figure 3. Seasonal dynamics of leaf traits of: P. tremula (A,C); and S. caprea (B,D). LMA, leaf mass per area (A,B); and EWT, equivalent water thickness (C,D). Different letters above the columns correspond to a significant difference, One-way ANOVA, Tukey–Kramer test. Significant differences if p < 0.05.
Figure 3. Seasonal dynamics of leaf traits of: P. tremula (A,C); and S. caprea (B,D). LMA, leaf mass per area (A,B); and EWT, equivalent water thickness (C,D). Different letters above the columns correspond to a significant difference, One-way ANOVA, Tukey–Kramer test. Significant differences if p < 0.05.
Sensors 17 01202 g003
Figure 4. Chlorophyll a + b content of P. tremula (A,C) and S. caprea (B,D) during the vegetative season (bars, left axis) and a seasonal course of selected chlorophyll indices (lines; right axis). Datt is presented as an index affected by the measurement setup, BGI as an index unaffected by the measurement setup. Left graphs (A,B) show absolute values of indices computed from the contact probe reflectance of a single leaf (1L) and a leaf stack (5L). Right graphs (C,D) show identical indices after z-score standardization. Chlorophyll content: different letters above the columns correspond to a significant difference, One-way ANOVA, Tukey–Kramer test. Significant differences if p < 0.05.
Figure 4. Chlorophyll a + b content of P. tremula (A,C) and S. caprea (B,D) during the vegetative season (bars, left axis) and a seasonal course of selected chlorophyll indices (lines; right axis). Datt is presented as an index affected by the measurement setup, BGI as an index unaffected by the measurement setup. Left graphs (A,B) show absolute values of indices computed from the contact probe reflectance of a single leaf (1L) and a leaf stack (5L). Right graphs (C,D) show identical indices after z-score standardization. Chlorophyll content: different letters above the columns correspond to a significant difference, One-way ANOVA, Tukey–Kramer test. Significant differences if p < 0.05.
Sensors 17 01202 g004
Figure 5. The difference (∆R5L–1L) between the reflectances measured on a leaf stack (5L) and a single leaf (1L) during the vegetative season, represented by curves for 6 months (May–October) for selected spectral ranges in P. tremula. The mean value of ∆R5L–1L for each month from May to October within: (A) VIS (400–700 nm); (B) NIR (740–1140 nm); (C) SWIR1 (1500–1800 nm); and (D) SWIR2 (2000–2400 nm). Different capital letters by months´ names correspond to a significant difference in ∆R5L–1L, computed from the mean reflectance values in the particular spectral range and the month of sampling. One-way ANOVA, Tukey–Kramer test. Significant differences if p < 0.05.
Figure 5. The difference (∆R5L–1L) between the reflectances measured on a leaf stack (5L) and a single leaf (1L) during the vegetative season, represented by curves for 6 months (May–October) for selected spectral ranges in P. tremula. The mean value of ∆R5L–1L for each month from May to October within: (A) VIS (400–700 nm); (B) NIR (740–1140 nm); (C) SWIR1 (1500–1800 nm); and (D) SWIR2 (2000–2400 nm). Different capital letters by months´ names correspond to a significant difference in ∆R5L–1L, computed from the mean reflectance values in the particular spectral range and the month of sampling. One-way ANOVA, Tukey–Kramer test. Significant differences if p < 0.05.
Sensors 17 01202 g005
Figure 6. The difference (∆R5L–1L) between the reflectances measured on a leaf stack (5L) and a single leaf (1L) during the vegetative season, represented by curves for 6 months (May–October) for selected spectral ranges in S. caprea. The mean value ∆R5L–1L for each month from May to October within: (A) VIS (400–700 nm), (B) NIR (740–1140 nm) and (C) SWIR1 (1500–1800 nm) and (D) SWIR2 (2000–2400 nm). Different capital letters by months´ names correspond to a significant difference in ∆R5L–1L, computed from the mean reflectance values in the particular spectral range and the month of sampling. One-way ANOVA, Tukey–Kramer test. Significant differences if p < 0.05.
Figure 6. The difference (∆R5L–1L) between the reflectances measured on a leaf stack (5L) and a single leaf (1L) during the vegetative season, represented by curves for 6 months (May–October) for selected spectral ranges in S. caprea. The mean value ∆R5L–1L for each month from May to October within: (A) VIS (400–700 nm), (B) NIR (740–1140 nm) and (C) SWIR1 (1500–1800 nm) and (D) SWIR2 (2000–2400 nm). Different capital letters by months´ names correspond to a significant difference in ∆R5L–1L, computed from the mean reflectance values in the particular spectral range and the month of sampling. One-way ANOVA, Tukey–Kramer test. Significant differences if p < 0.05.
Sensors 17 01202 g006
Figure 7. Chlorophyll a + b content of P. tremula (A,C) and S. caprea (B,D) during the season (bars; left axis) and a seasonal course of selected chlorophyll indices (lines; right axis) after z-score standardization. Left graphs (A,B) present indices unaffected by the measurement setup. Right graphs (C,D) present indices affected by the measurement setup. Graphs show absolute values after z-score standardization of indices computed from contact probe reflectance of a single leaf (1L-full lines) and a leaf stack (5L-dashed lines).
Figure 7. Chlorophyll a + b content of P. tremula (A,C) and S. caprea (B,D) during the season (bars; left axis) and a seasonal course of selected chlorophyll indices (lines; right axis) after z-score standardization. Left graphs (A,B) present indices unaffected by the measurement setup. Right graphs (C,D) present indices affected by the measurement setup. Graphs show absolute values after z-score standardization of indices computed from contact probe reflectance of a single leaf (1L-full lines) and a leaf stack (5L-dashed lines).
Sensors 17 01202 g007
Table 1. The first and the second columns list a name of an index and a formula of the index, respectively. The column R2-1L, 5L—shows the coefficient of determination for linear regression between the vegetation index and the chlorophyll content for a single leaf and a stack of leaves during the season. Underlined values of R2 correspond to a significant correlation. Indices defined as carotenoid- and water-related were excluded from the regression analysis (an empty lines). Months (May–October) show significance of the difference between indices computed from the reflectances of a single leaf and a leaf stack for each month separately tested by two-sample t-test, considered as significant if p < 0.05. According to the p-value, the level of significance is classified to three groups: not significant if p ≥ 0.05 coded in blue; significant if p = (0.05, 0.0000001) coded in red; very strongly significant if p < 0.0000001 coded in dark red. The majority of indices are listed in [5,22,24,31]
Table 1. The first and the second columns list a name of an index and a formula of the index, respectively. The column R2-1L, 5L—shows the coefficient of determination for linear regression between the vegetation index and the chlorophyll content for a single leaf and a stack of leaves during the season. Underlined values of R2 correspond to a significant correlation. Indices defined as carotenoid- and water-related were excluded from the regression analysis (an empty lines). Months (May–October) show significance of the difference between indices computed from the reflectances of a single leaf and a leaf stack for each month separately tested by two-sample t-test, considered as significant if p < 0.05. According to the p-value, the level of significance is classified to three groups: not significant if p ≥ 0.05 coded in blue; significant if p = (0.05, 0.0000001) coded in red; very strongly significant if p < 0.0000001 coded in dark red. The majority of indices are listed in [5,22,24,31]
Populus tremulaSalix caprea
IndexFormulaR2MayJuneJulAugSepOctR2MayJuneJulAugSepOct
1L5L1L5L
BGI=R450/R5500.420.33 0.200.16 [25]
BI=(R800+R670+R550)/sqrt(3)0.010.10 0.000.00 [43]
Carter4=R710/R760 0.520.42 0.400.33 [44]
CI=(R800-R550)/R8000.420.34 0.310.30 [43]
CRI550=(1/R515)-(1/R550) [45]
CRI700=(1/R515)-(1/R700) [45]
CTR=R695/R7600.400.29 0.360.22 [44]
CUR=(R675·R690)/R68330.020.01 0.000.00 [46]
Datt=(R850-R710)/(R850-R680) 0.510.44 0.350.35 [47]
Datt2=R850/R710 0.500.41 0.370.38 [47]
DD=(R749-R720)-(R701-R672) 0.530.42 0.390.32 [48]
DVI=R800-R6700.310.21 0.090.06 [49]
G=R554/R6770.210.15 0.120.10 [25]
Gitelson2=(R750-R800/R695-R740)-1 0.390.24 0.410.31 [45]
GM94a=R750/R7000.470.32 0.430.35 [50]
gNDVI780=(R780-R550)/(R780+R550)0.420.34 0.320.31 [51]
GRg=(R800/R550)-10.400.29 0.300.34 [45]
Macc01=(R780-R710)/(R780-R680)0.520.45 0.380.34 [52]
MCARI=((R700-R670)-0.2·(R700-R550))·(R700/R670)0.430.33 0.300.27 [7]
MCARI1=1.2·(2.5·(R800-R670)-1.3·(R800-R550))0.020.09 0.000.01 [53]
MCARI2=((R750-R705)-0.2·(R750-R550))·(R750/R705) 0.540.36 0.440.31 [54]
MCARI2/OSAV2 0.550.36 0.430.28 [54]
MCARI/OSAVI 0.460.35 0.310.28 [7]
McM94=R700/R6700.370.28 0.280.23 [55]
MND=(R750-R445)/(R750+R705-2·R445)0.530.42 0.420.34 [15]
mND705=(R750-R705)/(R750+R705-2·R445) 0.530.42 0.420.34 [15]
MNDVI1=(R755-R745)/(R755+R745)0.430.44 0.270.30 [56]
MNDVI8=(R755-R730)/(R755+R730)0.520.46 0.360.35 [56]
MNDVIre=(R750-R705)/(R750+R705-R445)0.520.41 0.420.34 [15]
MSAVI=0.5·(2·R800+1-sqrt((2·R800+1)2-8·(R800-R670)))0.280.22 0.060.06 [57]
MSI=R1600/R820 [58]
MSR=((R800-R670)-1)/sqrt((R800/R670)+1)0.080.09 0.110.06 [59]
MTCI=(R754-R709)/(R709-R681) 0.510.44 0.390.38 [60]
N705=(R705-R675)/(R750-R670)0.510.42 0.390.33 [1]
N715=(R715-R675)/(R750-R670)0.520.46 0.370.36 [1]
N725=(R725-R675)/(R750-R670)0.500.48 0.370.36 [1]
NDII2=(R820-R1650)/(R820+R1650) [61]
NDVI1=(R800-R670)/(R800+R670)0.090.10 0.110.06 [20]
NDVI800680=(R800-R680)/(R800+R680)0.100.10 0.110.06 [62]
NDVIre=(R750-R705)/(R750+R705)0.510.39 0.420.33 [50]
NDWI=(R858-R1240)/(R858+R1240) [63]
NDWI2130=(R858-R2130)/(R858+R2130) [64]
NMDI=(R860-(R1640-R2130))/(R860+(R1640-R2130)) [65]
NPCI=(R680-R430)/(R680+R430)0.250.16 0.130.08 [66]
OSAVI=(1+0.16)·(R800-R670)/(R800+R670+0.16)0.250.17 0.130.08 [30]
OSAVI2=(1+0.16)·(R750-R705)/(R750+R705+0.16)0.520.39 0.420.32 [54]
PRI=(R570-R530)/(R570+R530) 0.210.15 0.150.18 [67]
PRI·CI-H=(R680-R500)/R750 [68]
PRI·CI-Y=(R570-R530)/(R570+R530)·(R760/R700-1) [69]
PRIm1=(R515-R530)/(R515+R530)0.250.22 0.000.00 [68]
PSNDb=(R800-R635)/(R800+R635) 0.300.23 0.350.19 [24]
PSNDc=(R800-R470)/(R800+R470) [24]
PSSRa=R800/R680 0.080.05 0.130.05 [24]
PSSRb=R800/R635 0.270.16 0.400.24 [24]
PSSRc=R800/R470 [24]
RDVI=(R800-R670)/(sqrt(R800+R670))0.180.11 0.140.07 [70]
REIP=(700+40·((Rre-R700)/(R740-R700)))/1000.510.43 0.340.27 [71]
REP=700+40·((((R670+R780)/2)-R700)/(R740-R700))0.510.42 0.340.27 [72]
REP-Li=700+40·((R670+R780/2)/(R740-R700))0.510.42 0.340.27 [71]
RMSR=((R750/R705)-1)/sqrt((R750/R705)+1)0.500.38 0.420.35 [54]
RNIRCRI550=((1/R515-1/R550)·R770) [45]
RNIRCRI700=((1/R515-1/R700)·R770) [45]
Rre=(R670+R780)/20.250.22 0.050.05 [71]
SIPI=(R800-R445)/(R800-R680) 0.320.23 0.170.09 [73]
SIPI680=(R800-R455)/(R800-R680)0.320.24 0.120.08 [73]
SIPI705=(R800-R455)/(R800-R705)0.500.41 0.350.26 [73]
SR=R800/R6700.070.06 0.130.06 [49]
SR2=R750/R710 0.510.41 0.410.37 [74]
SR3=R750/R5500.400.29 0.330.34 [75]
SRPI=R430/R680 0.250.15 0.120.08 [73]
SRWI=R858/R1240 [76]
TCARI=3·((R700-R670)-0.2·(R700-R550)·(R700/R670))0.380.30 0.220.22 [77]
TCARI/OSAVI 0.420.33 0.270.25 [77]
TCARI2=3·((R750-R705)-0.2·(R750-R550)·(R750/R705))0.120.29 0.000.26 [54]
TCARI2/OSAVI2 0.410.35 0.280.35 [54]
TVI=0.5·(120·(R750-R550)-200·(R670-R550))0.080.09 0.000.01 [78]
VI[700]=(R700-R670)/(R700+R670)0.380.30 0.280.21 [79]
Vogelmann=R740/R720 0.510.44 0.400.37 [80]
Vogelmann2=(R734-R747)/(R715+R726) 0.520.45 0.380.36 [80]
WI=R900/R970 [81]
Table 2. Descriptive statistics of the reflectance difference (∆R5L–1L) between the reflectances measured on a leaf stack (5L) and a single leaf (1L) in % within VIS (400–700 nm), NIR (740–1140 nm), SWIR1 (1500–1800 nm) and SWIR2 (2000–2400 nm) regions for all data averaged from six months for P. tremula and S. caprea.
Table 2. Descriptive statistics of the reflectance difference (∆R5L–1L) between the reflectances measured on a leaf stack (5L) and a single leaf (1L) in % within VIS (400–700 nm), NIR (740–1140 nm), SWIR1 (1500–1800 nm) and SWIR2 (2000–2400 nm) regions for all data averaged from six months for P. tremula and S. caprea.
P. tremulaVISNIRSWIR1SWIR2
Variance0.0%0.2%0.1%0.0%
Max3.2%43.3%14.0%4.8%
Min−3.9%10.0%0.7%−0.7%
Mean0.0%28.7%7.2%1.3%
SD0.6%4.2%2.4%1.1%
S. capreaVISNIRSWIR1SWIR2
Variance0.0%0.2%0.1%0.0%
Max10.1%37.4%14.0%6.0%
Min−3.0%6.3%1.1%−0.9%
Mean0.1%26.3%7.2%1.6%
SD0.9%5.0%2.4%1.2%

Share and Cite

MDPI and ACS Style

Neuwirthová, E.; Lhotáková, Z.; Albrechtová, J. The Effect of Leaf Stacking on Leaf Reflectance and Vegetation Indices Measured by Contact Probe during the Season. Sensors 2017, 17, 1202. https://doi.org/10.3390/s17061202

AMA Style

Neuwirthová E, Lhotáková Z, Albrechtová J. The Effect of Leaf Stacking on Leaf Reflectance and Vegetation Indices Measured by Contact Probe during the Season. Sensors. 2017; 17(6):1202. https://doi.org/10.3390/s17061202

Chicago/Turabian Style

Neuwirthová, Eva, Zuzana Lhotáková, and Jana Albrechtová. 2017. "The Effect of Leaf Stacking on Leaf Reflectance and Vegetation Indices Measured by Contact Probe during the Season" Sensors 17, no. 6: 1202. https://doi.org/10.3390/s17061202

APA Style

Neuwirthová, E., Lhotáková, Z., & Albrechtová, J. (2017). The Effect of Leaf Stacking on Leaf Reflectance and Vegetation Indices Measured by Contact Probe during the Season. Sensors, 17(6), 1202. https://doi.org/10.3390/s17061202

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