Figures
Abstract
Leaves of various ages and positions in a plant's canopy can present distinct physiological, morphological and anatomical characteristics, leading to complexities in selecting a single leaf for spectral representation of an entire plant. A fortiori, as growth rates between canopies differ, spectral-based comparisons across multiple plants – often based on leaves' position but not age – becomes an even more challenging mission. This study explores the effect of differential growth rates on the reflectance variability between leaves of different canopies, and its implication on physiological predictions made by widely-used spectral indices. Two distinct irrigation treatments were applied for one month, in order to trigger the formation of different growth rates between two groups of grapevines. Throughout the experiment, the plants were physiologically and morphologically monitored, while leaves from every part of their canopies were spectrally and histologically sampled. As the control vines were constantly developing new leaves, the water deficit plants were experiencing growth inhibition, resulting in leaves of different age at similar nodal position across the treatments. This modification of the age-position correlation was characterized by a near infrared reflectance difference between younger and older leaves, which was found to be exponentially correlated (R2 = 0.98) to the age-dependent area of intercellular air spaces within the spongy parenchyma. Overall, the foliage of the control plant became more spectrally variable, creating complications for intra- and inter-treatment leaf-based comparisons. Of the derived indices, the Structure-Insensitive Pigment Index (SIPI) was found indifferent to the age-position effect, allowing the treatments to be compared at any nodal position, while a Normalized Difference Vegetation Index (NDVI)-based stomatal conductance prediction was substantially affected by differential growth rates. As various biotic and abiotic factors may form distinctions in growth, future precision agriculture studies should consider its spectral effect on physiological predictions.
Citation: Rapaport T, Hochberg U, Rachmilevitch S, Karnieli A (2014) The Effect of Differential Growth Rates across Plants on Spectral Predictions of Physiological Parameters. PLoS ONE 9(2): e88930. https://doi.org/10.1371/journal.pone.0088930
Editor: Manuel Reigosa, University of Vigo, Spain
Received: August 26, 2013; Accepted: January 13, 2014; Published: February 11, 2014
Copyright: © 2014 Rapaport et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: The reported project was funded by the Israeli Ministry of Agriculture (grant 857-0614-11). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Over the last few decades, improvements in the spatial and spectral capabilities of remote sensors have significantly promoted the precise monitoring of crop growth and development, from the field scale to the satellite level [1]–[4]. These technological advantages can provide a great assistance to farmers, provided that the variability within a plant's canopy – originating from spatial and phenological differences [5] – is understood. Physiologically, leaves of various ages and positions in a single plant were previously documented to have distinct levels of water content [6]–[8], water use efficiency [9], [10], stomatal conductance [11], [12], nitrogen content and allocation [13]–[15], photosynthesis [16]–[19], chlorophyll content [19]–[21], assimilation rates [12], [22], [23], and cellular structures and processes [24]. These within-canopy variations highlight the difficulty of representing an entire plant and, more so, comparing between different plants using a single leaf – normally chosen on the basis of the ‘age’ and ‘position’ concepts, which many studies refer to as equals (e.g., [25]–[28]). Indeed, when plants grow at a uniform rate, a constant age-position correlation can be assumed across them [22], meaning that leaves of a similar position are plausibly of a similar age. However, differences in leaf longevity or in growth rates, caused by genotypic differences [29] and biotic or abiotic stresses [30], can significantly modify this correlation [22], [31]. This hypothesis was already tested four decades ago by Wilson and Cooper [32], who concluded that “Differential rates of leaf appearance among plants may also affect comparisons”, and is especially notable for deciduous and annual plants that regrow their canopy on a yearly-basis. Dwyer and Stewart [33] reaffirmed this conclusion by showing that maximal photosynthetic capacity in maize correlates with leaf age but not with position. Generally, it seems that leaves can gain mature characteristics during growth inhibition, regardless of their position [34]. The specific control mechanism for this phenomenon, which is likely genetically- and environmentally dependent, is yet to be elucidated.
From a spectral perspective, Shibata [35] was one of the first to address the age-dependent variance within a canopy, when he correlated optical absorbance with chlorophyll formation in Phaseolus vulgaris leaves of different ages. Since then, pigment concentrations have been related to visible and infrared differences between young and mature leaves of Hedera helix [36], Carica papaya [20], Pinus taeda [37], [38], and of many other plant species (e.g., [39]–[43]). Spectral variability was also linked to changes in age-dependent histological variables within the leaf, such as the thickness of its tissues [44]–[46] and the air cavity distribution in its parenchyma layers [47]–[49]. The latter feature has been mainly linked to near infrared (NIR) reflectance, as the difference between the refraction index of cell walls and that of intercellular air voids – spaces that generally spread with leaf maturity – is likely the cause for scattering in this portion of the electromagnetic spectrum [50]–[54]. To date, only few remote sensing studies seemed to emphasize the importance of leaf age variability within a canopy in the context of spectral comparisons between plants. The most relevant and recent exception is the work of Liu et al. [41], who argued that “Few reports had simultaneously analyzed the relationship among the spectral reflectance indices, pigment concentrations, different plant species and developmental stages”. Their conclusion that “Some indices show greater promise as estimators of pigment concentrations than others” supports the ideas of Stone et al. [21], who stated that “Leaf-age composition of crowns needs to be taken into account when applying reflectance-based indices”. Another related example is the work of Isaacson et al. [55], who tried to normalize the phenological effect for common spectral indices using Landsat and MODIS data.
The current study was set to explore the effect of differential growth rates on the reflectance variability of leaves across different plants, and its impact on spectral-based physiological predictions. Specifically, the objectives were (1) to model the variance between individual leaves of canopies in the face of differential growth rates; (2) to relate this spectral variance among leaves of different age and position to the anatomical variance that causes it; and (3) to test the susceptibility of several widely-used reflectance indices to this spectral-anatomical effect, and their ability to predict physiological parameters in light of differential growth rates. To accomplish these objectives, two groups of ten grapevines were subjected to distinct water treatments, so that different growth rates could be established between them. During their development, the plants were physiologically and morphologically monitored in order to associate their water status to growth regulation. Additionally, leaves from every part of their canopies were spectrally and histologically sampled, so that reflectance measurements can be correlated with leaves' strata thickness and air cavity distribution. It was hypothesized that leaves of the same nodal position across the two irrigation groups could be gradually distinguished by their NIR reflectance values, due to anatomical differences, consequently limiting the ability of some spectral indices to physiologically compare between the stress treatments.
Methodology
Plant Materials and Growth Conditions
Twenty Cabernet Sauvignon one-year-old plants (Vitis vinifera L. cv.), grafted on Richter-110 rootstocks, were grown under controlled conditions in 10 liter pots filled with 9 liters of ‘RAM8’ media (Tuff Merom Golan, Merom Golan, Israel). The soil was enriched with a controlled-release fertilizer, and fresh water was applied on a daily basis through drip-irrigation in a rate of 2 L h−1. Day and night temperatures within the greenhouse were maintained between 28±1.5°C and 18±1.5°C, respectively. Plants were pruned to allow the growth of a single branch, while flower buds were removed upon emergence to sustain the vines in a vegetative state.
Experimental Design
Until the beginning of the experiment, all twenty plants were irrigated daily to field capacity. The experiment began once the plants had reached a size of 12–16 leaves and lasted for 30 days. While ten vines were kept on the same irrigation regime and maintained as controls (C), the other ten plants were subjected to a steady level of water deficit (WD). The stress treatment was applied according to the suggestions of Pou et al. [56] and is briefly described as follows. Irrigation was halted from day 1 until a volumetric water content (VWC) range of 5–10% was achieved in the pots. The plants were then kept at this stress level by daily replenishments of the amounts of water consumed, which were determined by weighing the pots every evening and subtracting the results from their corresponding target weights. All pots were arranged in the greenhouse according to a randomized complete block design, and were put sufficiently far from one another and from the greenhouse edges to minimize shading effects.
Differences in the stress level of the vine groups were physiologically assessed by measuring selected gas exchange parameters. These measurements were taken on the 1st day of experiment (DOE) before the stress was applied, on the 7th DOE in order to allow the WD pots to reach the desired VWC, and on the 29th DOE that was the last measuring day. Changes in vegetative development were evaluated by counting the number of nodes and measuring foliage area on the 1st, 15th, and 29th DOE. In order to keep track of leaves' age, they were marked upon reaching a size of 1 cm2. Spectral measurements of the vines' leaves were taken on the 1st, 7th, 15th, and 29th DOE. Tissue samples were collected on the 29th day for complementary histological measurements. All types of measurements (Table 1) were conducted between 10:00–13:00.
Physiological and Morphological Measurements
Stomatal conductance (gs) and net CO2 assimilation rates (AN) of the youngest, fully-matured leaf of each plant were collected using a portable Li-6400 Infrared Gas Analyzer (Li-Cor Biosciences Inc., NE, USA). The open photosynthesis system was equipped with an external carbon dioxide source, in order to maintain the leaf chamber at a concentration of 400 µmol mol−1. Temperature and relative humidity range were 28°C and 30–55%, respectively, and the photosynthetic active radiation (PAR) rate was 1000 µmol photons m−2 s−1. Gas exchange measurements were taken inside the greenhouse.
Foliage area was determined using a regression model that correlated the surface dimensions of a leaf with its actual area. After logging the width and length of 30 leaves, their pictures were then taken against a white background. The images were then made binary using The GIMP 2.6.11 software (http://www.gimp.org), so that the real areas could be derived from the sum of all relevant pixels. Eventually, a least squares linear regression model was found to best predict the actual leaf area from width and length (Image area = 1.72× Squared area +12.1; R2 = 0.95; p<0.01).
Spectral Measurements
At every time-point, an ASD Pro-FR Field Spectroradiometer (Analytical Spectral Devices Inc., CO, USA), equipped with an 1800-12S Integrating Sphere (Li-Cor Biosciences Inc., NE, USA), was used to log the signature of the oldest leaf in every fourth node of each vine, starting from the terminal bud. The system was set to measure mean radiance units (each measurement was averaged out of 20 scans) in the 350–2500 nm range, at 1 nm spectral resolution. Normalized values of leaf reflectance (ρ) were derived from complementary energy flux measurements of white-reference standards (BaSO4 tablets) and dark current, as suggested in the Integrating Sphere's manual. The reflectance data was used to exert several narrow-band, water stress-related indices: Normalized Difference Vegetation Index (NDVI; [57]; (ρ800 nm−ρ670 nm)/(ρ800 nm +ρ670 nm)), Water Index (WI; [58]; (ρ900 nm/ρ970 nm)), Normalized Difference Infrared Index (NDII; [59]; (ρ820 nm−ρ1650 nm)/(ρ820 nm +ρ1650 nm)), Moisture Stress Index (MSI; [59]; (ρ1600 nm/ρ820 nm)) and Structure Insensitive Pigment Index (SIPI; [60]; (ρ800 nm−ρ445 nm)/(ρ800 nm−ρ680 nm)).
Anatomical Measurements
Samples were collected from every 4th leaf in a plant – over five repetitions per treatment – starting from the terminal bud and ending at the twentieth leaf from the top. Cross sections (transversal cuts) were processed according to the protocol of Rewald et al. [61], which was slightly modified. Photomicrographs of the cross sections were then taken with an Axio Imager A1 Light Microscope (Carl Zeiss Microscopy Inc., TH, Germany) and ×100 magnification images were produced using its AxioVision 4.6.3 software. The GIMP software was then used for measuring the thickness of all leaf tissues, i.e., both cuticle, epidermal, and parenchyma layers, and in obtaining the percentage of intercellular air voids' area from the entire adaxial and abaxial mesophyll surface. The latter analysis was performed by making each section image binary, summing all void pixels in each mesophyll layer, and normalizing each area result to the total cell region of the respective parenchyma.
Statistical Analyses
All statistical calculations were performed using the SigmaPlot 12.0 software (Systat Software Inc., IL, USA). Differences between mean values were assessed using one-way Analysis of Variation (ANOVA) test, which was followed by Tukey's Honest Significance Difference (HSD) post-hoc test whenever necessary. Prior to those tests, assumptions of residuals' normality and homoscedasticity were checked and met using Shapiro-Wilk and Bartlett's tests, respectively. Least squares regressions were used in order to describe the continuous relationship between various variables. Results of all comparison tests and correlation analyses were considered significant at p<0.05.
Results
Physiological and Morphological Analyses
The water range for the WD pots was achieved between the 5th and 6th DOE, and was maintained at an average of 8% until the end of the experiment (Figure 1A). This water level was significantly different from that of the control pots, which were continuously refilled to a VWC of 33% by day 30. Significant physiological distinctions between the irrigation groups were noticed on the 7th DOE, as the stomatal conductance (Figure 1B) and net assimilation rates (Figure 1C) of the WD vines were lower than those of the control vines by 0.23 mol H2O m−2 s−1 (a 75% reduction) and 5.35 µmol CO2 m−2 s−1 (a 37% reduction), respectively. Those gas exchange differences remained similar until the end of the experiment.
Excluding day 1, values of the two groups were significantly different for every measuring date. Vertical bars represent means ± standard deviations (n = 10).
While the average control plant began the trial 2 nodes shorter than the average WD plant, the former concluded the experiment with 11 nodes more than the latter (Figure 2A). In terms of total leaf area, although both treatments were found similar at day 1 with 1730 cm2, the WD vines' foliage was 1860 cm2 smaller (a 40% reduction) than that of the control vines by the 29th DOE (Figure 2B). Growth inhibition was apparent in the water deficit treatment from the 3rd week onward, as its morphological values of the 15th DOE were insignificantly different from those of day 29.
Vertical bars represent means ± standard deviations (n = 10). Lower-case letters indicate significance of differences.
Spectral Analyses
By the end of the first week of the experiment, the reflectance signatures of leaves in the control vine (Figure 3A) were still similar to those of the WD plant (Figure 3B). During the second week, from day 7 to day 14, the 4th leaf of the control treatment has moved to the 8th nodal position, presenting a 7% NIR reflectance increase (measured by taking the average of 750–1250 nm; Figure 3C). On the other hand, the 4th leaf of the WD vine maintained its position, displaying only a minor reflectance transition (Figure 3D). As the foliage of the control plant continued to expand by day 29 (Figure 2), the NIR values of its younger leaves, down until the 12th node, were significantly lower than those of its older leaves (downward from the 16th node; Figure 3E). In the WD treatment, however, no new leaves emerged, while the existing ones maintained their position and their reflectance spectrum (Figure 3F). Eventually, by the end of the experiment, significant differences in NIR reflectance were found between control and stressed leaves of the same position down until the 16th node. At that time, the average control and WD plants could have been spectrally-represented solely by leaves in their 24 and 16 lower nodes, respectively, as their reflectance values were similar and the majority in each canopy (Figure 4).
The reflectance signatures presented are of the 7th, 14th, and 29th days of measuring for the control (A, C, and E, respectively) and water deficit (B, D, and F, respectively) groups. Each reflectance signature was averaged out of 10 leaf samples.
Significant differences in the means of near infrared reflectance (750–1250 nm) were found between control and stressed leaves of the same nodal position down until the 16th node. The reflectance values in the upper 15 and 11 nodes of the control and water deficit groups, respectively, were significantly different from those in the rest of the nodes, and thus were not representative for most of the plant. Vertical bars represent means ± standard deviations (n = 10). Lower-case letters indicate significance of differences.
The effect of growth rates on selected spectral indices was assessed on the 29th DOE (Table 2). In the water deficit treatment, characterized by low growth rates, most indices presented uniformity in values for all leaves measured – regardless of their position. This was a direct result of the similar reflectance signatures of leaves in the majority of nodes. In contrast, most control indices displayed value inconsistencies at higher leaf positions (4th-12th nodes). The greatest variation was observed in the NDVI and WI profiles, wherein leaves at the top 12 nodes (between 41–55% of the foliage) showed significantly different values as compared with leaves in the rest of the canopy. The least susceptible indices, second only to the indifferent SIPI, were the NDII and MSI that were insensitive to value changes at the bottom 32 nodes (between 59–73% of the foliage).
To assess the extent of possible deviations caused by spectral variability, the NDVI values of the 29th DOE (Table 2) were used. Throughout the experiment, the index values of twenty plants – calculated from the same youngest, fully-matured leaves that were used in the gas exchange measurements – were regressed against their stomatal conductance values, resulting in a statistically significant linear correlation (R2 = 0.68; p<0.05; Figure 5). For each treatment of the last measuring day, the regression was then used to compare the average gs value of all measured leaves in the canopy to the average value predicted by only the leaves in the representative portion of the foliage (according to the NDVI value distribution shown in Table 2). Regarding the control vine, it was found that the exclusion of the 8 uppermost nodal positions, which were about 41% of the total foliage area, has increased the gs rate from 0.050 to 0.292 mol H2O m−2 s−1. A smaller increase in stomatal conductance, from 0.031 to 0.057 mol H2O m−2 s−1, was observed when the 4 uppermost nodal positions were discarded from the WD vine, because – in comparison to the control canopy – the excluded leaves were only about 32% of the total foliage area and their NDVI values were of closer resemblance to those of older leaves in the canopy.
The diamonds represent values of the youngest, fully-matured leaves of both treatments from all time points, and were used to create a linear regression. The filled square and circle are the 29th day's mean values of leaves from the whole control and water deficit canopies, respectively. The unfilled square and circle are the 29th day's mean values from only the representative portion of the control and water deficit canopies (Table 2), respectively.
Anatomical Analyses
Thickness measurements of all leaf strata, including both cuticle, epidermal, and mesophyll layers, revealed insignificant differences between most nodes in both treatments (Table 3). Leaf position along the vine branches was also found to poorly affect the intercellular void areas within the palisade parenchyma (Table 4), as little evidence of air cavity expansion was found between young leaves of both the control and water stress treatments (Figures 6A and 6B, respectively) and their older counterparts (Figures 6C and 6D, respectively). However, leaf maturation in both treatments did result in a prominent visual increase in air cavities' area within the spongy mesophyll, backed up by significant differences between most nodal positions (Table 4). A statistically significant, strong exponential correlation (R2 = 0.98; p<0.05) was also found between the area of air spaces within the abaxial parenchyma of all leaves and NIR reflectance (Figure 7).
The scale bars of both control and water deficit's 4th leaves (A and B, respectively) and 16th leaves (C and D, respectively) equal 50 µm.
Vertical and horizontal bars represent means ± standard deviations (n = 10 for the near infrared reflectance; n = 5 for the spongy mesophyll void area).
Discussion
In order to investigate the spectral effects of differential growth rates on physiological predictions, Cabernet Sauvignon grapevines were subjected to two distinct irrigation regimes. Throughout the experiment, the stress level that was applied to the water deficit (WD) treatment group (Figure 1) can be regarded as moderate, and its occurrence under field conditions is not rare [62]–[64]. The observed, consequent down-regulation of stomatal conductance (gs) in the WD vines was likely the cause for their steeply-reduced CO2 assimilation rates (AN; Figure 1) [17], [65]–[68]. The AN restrictions, in turn, were presumably responsible for the inhibition of vegetative growth in the water stressed group (Figures 2) [69], [70]. Since the expansion of the control vines' foliage was not limited by water shortage, a distinct correlation between leaves' age and position was gradually forming across the treatment groups. From a spectral perspective, the breaking of the age-position relation was mainly expressed through a growing difference in the NIR signatures of parallel leaves of a higher position (Figure 3). Namely, while a leaf of the WD group maintained its position along the branch but still continued to age and increase in reflectivity, its control counterpart was being constantly replaced by a younger leaf with significantly lower NIR reflectance values. In this regard, it should be noted that a differential effect of the water stress factor itself on the structure of young and mature leaves in the same plant cannot be ruled out [5], [34], [41], [71]. However, the results of previous studies, in which a similar spectral pattern of leaf maturation was shown in light of other stressing agents and conditions (e.g., [5], [20], [47], [52], [72], [73]), seem to support our proposed mechanism and suggest that the age variability mainly influenced the differences in reflectance. Therefore, it is argued that the canopy of the fully-growing plant became, over time, younger and more spectrally-variable than that of the stressed plant (Figure 4), creating complications for the selection of representative leaves for intra- and inter-treatment comparisons.
The aging of grapevine leaves is known to be accompanied by an increase in their mesophyll layers' intercellular air spaces [74]–[76]. As was previously suggested for various Vitis vinifera cultivars (e.g., [8], [77]), the dependence of NIR reflectance on leaf maturity is likely a derivative of compaction changes within the adaxial and abaxial parenchyma, as scattering in this spectrum generally increases with the increase in air void distribution in both mesophyll layers. In contrast to earlier studies that correlated NIR reflectance to the number of air voids in both mesophyll layers, to the total area of cell wall-air cavity refraction zones, or to the thickness of the leaves' strata (e.g., [39], [45]–[48], [53]), the results of this study indicate that NIR reflectance in a Cabernet Sauvignon leaf is mainly a function of the percentage of air spaces' area from the spongy mesophyll's area (Figures 6 & 7 and Table 4) and not of the thickness of any strata (Table 3). These findings are in agreement with those of Gausman et al. [78], Slaton et al. [79], and Sims and Gamon [42], who found relatively weak correlations between NIR reflectance and the thickness of leaves' layers of other plant species, and also with the conclusions of Nobel et al. [80], Slaton et al. [79], Castro and Sanchez-Azofeifa [81], and Ollinger [82], which emphasized the importance of the correlation between NIR reflectance and the normalized area of the intercellular air spaces within the spongy mesophyll. The spectral significance of the age-dependent variability in intercellular void areas among leaves is manifested by its contribution on NIR reflectance. For example, in their study on spectroscopic monitoring of osmotic treatments, Gausman et al. [51] suggested that leaf maturation and internal structure changes have negatively affected their spectral representations. Similarly, in their study on predicting leaf water content through reflectance measurements, Thomas et al. [83] suggested that the poor correlation that was found was presumably the cause of leaf structural variations during development. Thus, overlooking the spectral differences between young and old leaves might prove to be a significant source of bias in reflectance measurements [41].
Growth regulation is also likely to affect many common vegetative spectral indices, especially those that are derived from NIR reflectance, absorbance, or transmittance data to monitor plant physiology. For instance, the NDVI and WI, which were established as indicative parameters for plant status and are widely used (e.g., [8], [39], [41], [57], [58], [84]–[87]), were found to be significantly affected by leaf age and, therefore, in high growth rates, might produce values that are unrepresentative for most leaves in the canopy (Table 2). The NDII, MSI, and SIPI, however, showed little to no sensitivity to aging implications and, thus, have a clear advantage as indicators under differential growth rates. The index sensitivity results are in accordance with those found in a similar, leaf age-directed study by Liu et al. [41], wherein the SIPI was also found to be the most insensitive to leaf ontogenesis, while the NDVI was found to be affected by it. The impact of within-foliage spectral variability was assessed in the current study by regressing values of gs against values of NDVI – two variables that have been correlated in previous Vitis vinifera studies (e.g., [57], [87]). Calculating the weighted-average NDVI value of all leaves in the control plant's canopy has resulted in a gs value of 0.050 mol H2O m−2 s−1, which was 0.242 mol H2O m−2 s−1 lower than the weighted-average NDVI prediction of only the representative leaves (around 45% of the canopy; Figure 5). Although a smaller difference in gs was apparent in the case of the water deficit group, due to a higher ratio of representative leaves, the inclusion of all leaves in the canopy has resulted in a weighted-average NDVI that nearly zeroed the gs of otherwise 0.057 mol H2O m−2 s−1. Similarly, when assessing potential physiological deviations in WI-based regressions of other studies, it seems that the inclusion of younger leaves in the predictions of Rodriguez-Pérez et al. [8] will result in a 30–40% higher water content as percent of total fresh mass (WCt), and their inclusion in the predictions of Sims and Gamon [84] will result in a 0.5 kg m−2 higher water content of thin tissues. These distinctions strongly imply that the use of spectral indices without considering the inherent variability within plants' canopies might lead to either an under- or over-estimation of health status parameters [41]. Furthermore, special attention should be given to inter-treatment comparative images taken by satellite or airborne remote sensors of very high spatial resolution, which are generally limited to the upper part of the foliage and suffer from background effects, complex canopy geometry issues, and disadvantageous viewing angles [49], [88], [89].
As was demonstrated in this study, growth regulation can break the age-position correlation across stress treatments, leading to incomparable physiological values. Although the implications of vegetative growth inhibition were discussed here in light of water deficit alone, the findings of this research are not limited to any other biotic or abiotic factor, as plant development can be affected by diseases [90], shoot trimming [91], pests [92], defoliation [93], [94], light deficiency [95], and many other stressing agents. Moreover, since growth rate differences may occur between various cultivars of the same species (e.g., [96], [97]), and since spectral studies may include combinations of deciduous and non-deciduous plants, growth regulation becomes an even more common and complex problem altogether. Ergo, before attempting to compare leaves of different plants and treatments, future high spatial resolution-remote sensing works should take into account the specific spectral variability within the canopies of their choice, with respect to the specific type of stress induced. This idea may prove to be more difficult to implement than the protocol of the current study, as age- and stress-induced changes in leaf structure, pigment concentrations, and biochemical components need to be spectrally separated. However, combining high spectral resolution spectroscopy with advanced radiative models (e.g., [98], [99]) or advanced statistical tools [5] may lead to a successful discrimination between the impacts of maturation and stress, thus allowing a better understanding of differential growth rates' spectral effects.
Summary and Conclusions
In this study, the spectral-physiological implications of overlooking differential growth rates across stress treatments were discussed. As was shown, the variance in the growth rates has gradually established a new and distinct correlation between leaf age and position for each of the treatments, depicted by NIR reflectance differences among leaves of the same nodal degree. The spectral difference between young leaves of both treatment groups was anatomically attributed to the difference between their spongy mesophylls' intercellular void areas, which was affected by leaf maturation and controlled NIR scattering. Differences in the sensitivity level of derived reflectance indices to the age-position effect have demonstrated that some of them may have more potential than others under terms of differential growth rates, and are likely to predict less biased physiological values, especially for annual and deciduous species. Therefore, the methodological flaw that was described in this study – manifested by many forms of stress on the one hand and rarely discussed in the remote sensing literature on the other hand – should be taken into account in future precision agriculture studies.
Acknowledgments
We would like to express our deep appreciation to Mrs. Tanya Gendler, Mr. Alexander Goldberg, and Mrs. Liron Summerfield for greatly assisting in data retrieval and analysis. We would also like to thank four anonymous reviewers for their feedback and constructive comments that helped greatly improving this manuscript.
Author Contributions
Conceived and designed the experiments: TR UH SR AK. Performed the experiments: TR UH. Analyzed the data: TR UH. Contributed reagents/materials/analysis tools: TR UH SR AK. Wrote the paper: TR UH SR AK.
References
- 1. Moran MS, Inoue Y, Barnes EM (1997) Opportunities and limitations for image-based remote sensing in precision crop management. Remote Sensing of Environment 61: 319–346.
- 2. Hall A, Lamb DW, Holzapfel B, Louis J (2002) Optical remote sensing applications in viticulture – a review. Australian Journal of Grape and Wine Research 8: 36–47.
- 3. Seelan SK, Laguette S, Casady GM, Seielstad GA (2003) Remote sensing applications for precision agriculture: A learning community approach. Remote Sensing of Environment 88: 157–169.
- 4. Mulla DJ (2013) Twenty five years of remote sensing in precision agriculture: Key advances and remaining knowledge gaps. Biosystems Engineering 114: 358–371.
- 5.
Barry K, Corkrey R, Stone C, Mohammed C (2009) Characterizing Eucalypt leaf phenology and stress with spectral analysis. In: Jones S, Reinke K, editors. Innovations in Remote Sensing and Photogrammetry. Berlin: Springer-Verlag. pp. 193–210.
- 6. McCaig T, Romagosa I (1991) Water status measurements of excised wheat leaves: Position and age effects. Crop Science 31: 1538–1588.
- 7. Rabas A, Martin C (2003) Movement of water from old to young leaves in three species of succulents. Annals of Botany 92: 529–536.
- 8. Rodriguez-Pérez JR, Riaño D, Carlisle E, Ustin S, Smart DR (2007) Evaluation of hyperspectral reflectance indexes to detect grapevine water status in vineyards. American Journal of Enology and Viticulture 58: 302–317.
- 9. Lin ZF, Ehleringer JR (1982) Effects of leaf age on photosynthesis and water use efficiency of papaya. Photosynthetica 16: 514–519.
- 10. Medrano H, Pou A, Tomás M, Martorell S, Gulias J, et al. (2012) Average daily light interception determines leaf water use efficiency among different canopy locations in grapevine. Agricultural Water Management 114: 4–10.
- 11. Field C, Mooney HA (1983) Leaf age and seasonal effects on light, water, and nitrogen use efficiency in a California shrub. Oecologia 56: 348–355.
- 12. Snider JL, Choinski JSJ, Wise RR (2009) Juvenile Rhus glabra leaves have higher temperatures and lower gas exchange rates than mature leaves when compared in the field during periods of high irradiance. Journal of Plant Physiology 166: 686–696.
- 13. Connor DJ, Sadras VO, Hall AJ (1995) Canopy nitrogen distribution and the photosynthetic performance of sunflower crops during grain filling – a quantitative analysis. Oecologia 101: 274–281.
- 14. Wilson KB, Baldocchi DD, Hanson PJ (1999) Spatial and seasonal variability of photosynthetic parameters and their relationship to leaf nitrogen in a deciduous forest. Tree Physiology 20: 565–578.
- 15. Prieto JA, Louarn G, Perez-Peña J, Ojeda H, Simonneau T, et al. (2012) A leaf gas exchange model that accounts for intra-canopy variability by considering leaf nitrogen content and local acclimation to radiation in grapevine (Vitis vinifera L.). Plant, Cell & Environment 35: 1313–1328.
- 16. Poni S, Intrieri C, Silvestroni O (1994) Interaction of leaf age, fruiting, and exogeneous cytokinins in Sangiovese grapevines under non-irrigated conditions. American Journal of Enology and Viticulture 45: 71–78.
- 17. Escalona JM, Flexas J, Medrano H (1999) Stomatal and non-stomatal limitations of photosynthesis under water stress in field-grown grapevines. Australian Journal of Plant Physiology 26: 421–433.
- 18. Kennedy RA, Johnson D (1981) Changes in photosynthetic characteristics during leaf development in apple. Photosynthesis Research 2: 213–223.
- 19. Choinski JSJ, Ralph P, Eamus D (2003) Changes in photosynthesis during leaf expansion in Corymbia gummifera. Australian Journal of Botany 51: 111–118.
- 20. Lin ZF, Ehleringer JR (1982) Changes in spectral properties of leaves as related to chlorophyll content and age of papaya. Photosynthetica 16: 520–525.
- 21. Stone C, Chisholm L, McDonald S (2005) Effects of leaf age and psyllid damage on the spectral reflectance properties of Eucalyptus saligna foliage. Australian Journal of Botany 53: 45–54.
- 22. Kitajima K, Mulkey SS, Samaniego M, Wright SJ (2002) Decline of photosynthetic capacity with leaf age and position in two tropical pioneer tree species. American Journal of Botany 89: 1925–1932.
- 23. Escalona JM, Flexas J, Bota J, Medrano H (2003) Distribution of leaf photosynthesis and transpiration within grapevine canopies under different drought conditions. Vitis 42: 57–64.
- 24. Barna B, Györgyi B (1992) Resistence of young versus old tobacco leaves to necrotrophs, fusaric acis, cell wall-degrading enzymes and autolysis of membrane lipids. Physiological and Molecular Plant Pathology 40: 247–257.
- 25. Datt B (1999) A new reflectance index for remote sensing of chlorophyll content in higher plants: Tests using Eucalyptus leaves. Journal of Plant Physiology 154: 30–36.
- 26. El-Shikha DM, Waller P, Hunsaker D, Clarke T, Barnes E (2007) Ground-based remote sensing for assessing water and nitrogen status of broccoli. Agricultural Water Management 92: 183–193.
- 27. Seelig HD, Hoehn A, Stodieck LS, Klaus DM, Adams WW, et al. (2008) Relations of remote sensing leaf water indices to leaf water thickness in cowpea, bean, and sugarbeet plants. Remote Sensing of Environment 112: 445–455.
- 28. Eitel JUH, Vierling LA, Long DS, Hunt ER (2011) Early season remote sensing of wheat nitrogen status using a green scanning laser. Agricultural and Forest Meteorology 151: 1338–1345.
- 29. Poorter H, Remkes C (1990) Leaf area ratio and net assimilation rate of 24 wild species differing in relative growth rate. Oecologia 83: 553–559.
- 30. Hsiao TC, Acevedo E, Fereres E, Henderson D (1976) Water stress, growth, and osmotic adjustment. Philosophical Transactions of the Royal Society of London B Biological Sciences 273: 479–500.
- 31. Ackerly DD, Bazzaz FA (1995) Leaf dynamics, self-shading and carbon gain in seedlings of a tropical pioneer tree. Oecologia 101: 289–298.
- 32. Wilson D, Cooper J (1969) Apparent photosynthesis and leaf characters in relation to leaf position and age, among contrasting Lolium genotypes. New Phytologist 68: 645–655.
- 33. Dwyer L, Stewart D (1986) Effect of leaf age and position on net photosynthetic rates in maize (Zea mays L.). Agricultural and Forest Meteorology 37: 29–46.
- 34. Schurr U, Heckenberger U, Herdel K, Walter A, Feil R (2000) Leaf development in Ricinus communis during drought stress: Dynamics of growth processes, of cellular structure and of sink-source transition. Journal of Experimental Botany 51: 1515–1529.
- 35. Shibata K (1957) Spectroscopic studies on chlorophyll formation in intact leaves. Journal of Biochemistry 44: 147–173.
- 36. Gates DM, Keegan HJ, Schleter JC, Weidner VR (1965) Spectral properties of plants. Applied Optics 4: 11–20.
- 37. Carter GA, Paliwal K, Pathre U, Green T, Mitchell R, et al. (1989) Effect of competition and leaf age on visible and infrared reflectance in pine foliage. Plant, Cell & Environment 12: 309–315.
- 38. Carter GA, Knapp AK (2001) Leaf optical properties in higher plants: Linking spectral characteristics to stress and chlorophyll concentration. American Journal of Botany 88: 677–684.
- 39. Gamon JA, Surfus JS (1999) Assessing leaf pigment content and activity with a reflectometer. New Phytologist 143: 105–117.
- 40.
Sabat SC, Grover A, Mohanty P (1989) Selective alternations in photosynthetic pigment characteristics and photoelectron transport during senescence of wheat leaves. In: Singhal GS, Barber J, Dilley RA, Govindjee, Haselkorn R et al., editors. Photosynthesis: Molecular biology and bioenergetics. New Delhi, India: Springer-Verlag, Narosa Publishing House. pp. 343–351.
- 41. Liu N, Lin ZF, van Devender A, Lin GZ, Peng CL, et al. (2009) Spectral reflectance indices and pigment functions during leaf ontogenesis in six subtropical landscape plants. Plant Growth Regulation 58: 73–84.
- 42. Sims DA, Gamon JA (2002) Relationships between leaf pigment content and spectral reflectance across a wide range of species, leaf structures and developmental stages. Remote Sensing of Environment 81: 337–354.
- 43. Serrano L (2008) Effects of leaf structure on reflectance estimates of chlorophyll content. International Journal of Remote Sensing 29: 5265–5274.
- 44. Gausman HW, Allen WA (1973) Optical parameters of leaves of 30 plant species. Plant Physiology 52: 57–62.
- 45. Knapp AK, Carter GA (1998) Variability in leaf optical properties among 26 species from a broad range of habitats. American Journal of Botany 85: 940–946.
- 46. Ourcival JM, Joffre R, Rambal S (1999) Exploring the relationships between reflectance and anatomical and biochemical properties in Quercus ilex leaves. New Phytologist 143: 351–364.
- 47. Gausman HW, Allen WA, Cardenas R, Richardson AJ (1971) Effect of nodal position on absorption and scattering coefficients and infinite reflectance of cotton leaves, Gossypium hirsutum L. Agronomy Journal 63: 87–91.
- 48. Grant L (1987) Diffuse and specular characteristics of leaf reflectance. Remote Sensing of Environment 22: 309–322.
- 49. Meggio F, Zarco-Tejada PJ, Miller JR, Martín P, Gonzalez MR, et al. (2008) Row orientation and viewing geometry effects on row-structured vine crops for chlorophyll content estimation. Canadian Journal of Remote Sensing 34: 220–234.
- 50. Allen WA, Gausman HW, Richardson AJ, Thomas JW (1969) Interaction of isotropic light with a compact leaf. Journal of the Optical Society of America 59: 1376–1379.
- 51. Gausman HW, Allen WA, Cardenas R (1969) Reflectance of cotton leaves and their structure. Remote Sensing of Environment 1: 19–22.
- 52. Gausman HW, Allen WA, Cardenas R, Richardson AJ (1970) Relation of light reflectance to histological and physical evaluations of cotton leaf maturity. Applied Optics 9: 545–552.
- 53. Knipling EB (1970) Physical and physiological basis for the reflectance of visible and near-infrared radiation from vegetation. Remote Sensing of Environment 1: 155–159.
- 54. Gausman HW, Allen WA, Cardenas R, Escobar DE (1974) Refractive index of plant cell walls. Applied Optics 13: 109–111.
- 55. Isaacson BN, Serbin SP, Townsend PA (2012) Detection of relative differences in phenology of forest species using Landsat and MODIS. Landscape Ecology 27: 529–543.
- 56. Pou A, Medrano H, Tomàs M, Martorell S, Ribas-Carbó M, et al. (2012) Anisohydric behaviour in grapevines results in better performance under moderate water stress and recovery than isohydric behaviour. Plant and Soil 359: 335–349.
- 57. Zarco-Tejada PJ, González-Dugo V, Williams LE, Suárez L, Berni JAJ, et al. (2013) A PRI-based water stress index combining structural and chlorophyll effects: Assessment using diurnal narrow-band airborne imagery and the CWSI thermal index. Remote Sensing of Environment 138: 38–50.
- 58. Peñuelas J, Filella I, Biel C, Serrano L, Savé R (1993) The reflectance at the 950–970 nm region as an indicator of plant water status. International Journal of Remote Sensing 14: 1887–1905.
- 59. Hunt ER, Rock BN (1989) Detection of Changes in leaf water content using near- and middle-infrared reflectance. Remote Sensing of Environment 30: 43–54.
- 60. Peñuelas J, Inoue Y (1999) Reflectance indices indicative of changes in water and pigment contents of peanut and wheat leaves. Photosynthetica 36: 355–360.
- 61. Rewald B, Raveh E, Gendler T, Ephrath JE, Rachmilevitch S (2012) Phenotypic plasticity and water flux rates of citrus root orders under salinity. Journal of Experimental Botany 63: 2717–2727.
- 62. Flexas J, Escalona JM, Evain S, Gulias J, Moya I, et al. (2002) Steady-state chlorophyll fluorescence (FS) measurements as a tool to follow variations of net CO2 assimilation and stomatal conductance during water-stress in C3 plants. Physiologia Plantarum 114: 231–240.
- 63. Medrano H, Escalona JM, Bota J, Gulias J, Flexas J (2002) Regulation of photosynthesis of C3 plants in response to progressive drought: Stomatal conductance as a reference parameter. Annals of Botany 89: 895–905.
- 64. Cifre J, Bota J, Escalona JM, Medrano H, Flexas J (2005) Physiological tools for irrigation scheduling in grapevine (Vitis vinifera L.): An open gate to improve water-use efficiency? Agriculture, Ecosystems & Environment 106: 159–170.
- 65. Flexas J, Bota J, Escalona JM, Sampol B, Medrano H (2002) Effects of drought on photosyntheis in grapevines under field conditions: An evaluation of stomatal and mesophyll limitations. Functional Plant Biology 29: 461–471.
- 66. Flexas J, Baron M, Bota J, Ducruet J-M, Galle A, et al. (2009) Photosynthesis limitations during water stress acclimation and recovery in the drought-adapted Vitis hybrid Richter-110 (V. berlandieri x V. rupestris). Journal of Experimental Botany 60: 2361–2377.
- 67. Santesteban LG, Miranda C, Royo JB (2009) Effect of water deficit and rewatering on leaf gas exchange and transpiration decline of excised leaves of four grapevine (Vitis vinifera L.) cultivars. Scientia Horticulturae 121: 434–439.
- 68. Hochberg U, Degu A, Fait A, Rachmilevitch S (2013) Near isohydric grapevine cultivar displays higher photosynthetic efficiency and photorespiration rates under drought stress as compared with near anisohydric grapevine cultivar. Physiologia Plantarum 147: 443–452.
- 69. Flexas J, Medrano H (2002) Drought-inhibition of photosynthesis in C3 plants: Stomatal and non-stomatal limitations revisited. Annals of Botany 89: 183–189.
- 70. Padgett-Johnson M, Williams LE, Walker MA (2003) Vine growth relations, gas exchange, and vegetative growth of seventeen Vitis species grown under irrigated and nonirrigated conditions in California. Journal of the American Society of Horticultural Science 128: 269–276.
- 71. Aldakheel YY, Danson FM (1997) Spectral reflectance of dehydrating leaves: Measurements and modelling. International Journal of Remote Sensing 18: 3683–3690.
- 72. Gausman HW, Allen WA, Cardenas R, Richardson AJ (1973) Reflectance discrimination of cotton and corn at four growth stages. Agronomy Journal 65: 194–198.
- 73. Sinclair TR, Hoffer RM, Schreiber MM (1971) Reflectance and internal structure of leaves from several crops during a growing season. Agronomy Journal 63: 864–868.
- 74. Kriedemann PE (1968) Photosynthesis in vine leaves as a function of light intensity, temperature and leaf age. Vitis 7: 213–220.
- 75. Kriedemann PE, Kliewer WM, Harris JM (1970) Leaf age and photosynthesis in Vitis vinifera L. Vitis. 9: 97–104.
- 76. Hunter JJ, Visser JH (1988) The effect of partial defoliation, leaf position and developmental stage of the vine on the photosynthetic activity of Vitis vinifera L. cv. Cabernet Sauvignon. South African Journal of Enology and Viticulture 9: 9–15.
- 77. Qin JL, Rundquist D, Gitelson AA, Steele M, Harkins C (2010) A non–linear model for measuring grapevine leaf thickness by means of red-edge/near-infrared spectral reflectance. Acta Ecologica Sinica 30: 297–303.
- 78.
Gausman HW, Wigand CL, Escobar DE, Rodriguez RR, Richardson AJ (1973) The leaf mesophyll of twenty crops, their light spectra, and optical and geometrical parameters. Weslaco, TX: USA: U.S Agricultural Research Service, Soil and Water Conservation Research Devision. 88 p.
- 79. Slaton MR, Hunt ER, Smith WK (2001) Estimating near-infrared leaf reflectance from leaf structural characteristics. American Journal of Botany 88: 278–284.
- 80.
Nobel PS, Zaragoza LJ, Smith WK (1975) Relation between mesophyll surface area, photosynthetic rate and illumination level during development for leaves of Plecranthus parviflorus Henckel. Plant Physiology 55.
- 81. Castro KL, Sanchez-Azofeifa GA (2008) Changes in spectral properties, chlorophyll content and internal mesophyll structure of senescing Populus balsamifera and Populus tremuloides leaves. Sensors 8: 51–69.
- 82. Ollinger SV (2011) Sources of variability in canopy reflectance and the convergent properties of plants. New Phytologist 189: 375–394.
- 83. Thomas JR, Namken LN, Oerther GF, Brown RG (1971) Estimating leaf water content by reflectance measurements. Agronomy Journal 63: 845–847.
- 84. Sims DA, Gamon JA (2003) Estimation of vegetation water content and photosynthetic tissue area from spectral reflectance: A comparison of indices based on liquid water and chlorophyll absorption features. Remote Sensing of Environment 84: 526–537.
- 85. Peñuelas J, Filella I, Llusiá J, Siscart D, Piñol J (1998) Comparative field study of spring and summer leaf gas exchange and photobiology of the mediterranean trees Quercus ilex and Phillyrea latifolia. Journal of Experimental Botany 49: 229–238.
- 86. Rahimzadeh-Bajgiran P, Munehiro M, Omasa K (2012) Relationships between the photochemical reflectance index (PRI) and chlorophyll fluorescence parameters and plant pigment indices at different leaf growth stages. Photosynthesis Research 113: 261–271.
- 87. Baluja J, Diago MP, Balda P, Zorer R, Meggio F, et al. (2012) Assessment of vineyard water status variability by thermal and multispectral imagery using an unmanned aerial vehicle (UAV). Irrigation Science 30: 511–522.
- 88. Zarco-Tejada P, Berjón A, López-Lozano R, Miller J, Martín P, et al. (2005) Assessing vineyard condition with hyperspectral indices: Leaf and canopy reflectance simulation in a row-structured discontinuous canopy. Remote Sensing of Environment 99: 271–287.
- 89. Colombo R, Meroni M, Marchesi A, Busetto L, Rossini M, et al. (2008) Estimation of leaf and canopy water content in popular plantations by means of hyperspectral indices and inverse modeling. Remote Sensing of Environment 112: 1820–1834.
- 90. Zhu Y, Du B, Qian J, Zou B, Hua J (2013) Disease resistance gene-induced growth inhibition is enhanced by rcd1 independent of defense activation in Arabidopsis. Plant Physiology 161: 2005–2013.
- 91. Poni S, Giachino E (2000) Growth, photosynthesis and cropping of potted grapevines (Vitis vinifera L. cv. Cabernet Sauvignon) in relation to shoot trimming. Australian Journal of Grape and Wine Research 6: 216–226.
- 92. Lenz MS, Isaacs R, Flore JA, Howell GS (2009) Vegetative growth responses of Pinot gris (Vitis vinifera L.) grapevines to infestation by Potato Leaf hoppers (Empoasca fabae Harris). American Journal of Enology and Viticulture 60: 130–137.
- 93. Hunter JJ, Visser JH (1990) The effect of partial defoliation on growth characteristics of Vitis vinifera L. cv. Cabernet Sauvignon. II. Reproductive Growth. South African Journal of Enology and Viticulture 11: 26–32.
- 94. Poni S, Bernizzoni F, Civardi S (2008) The effect of early leaf removal on whole-canopy gas exchange and vine performance of Vitis vinifera L. ‘Sangiovese’. Vitis 47: 1–6.
- 95. Carvalho L, Santos P, Amancio S (2002) Effect of light intensity and CO2 concentration on growth and the acquisition of in vivo characteristics during acclimatization of grapevine regenerated in vitro. Vitis 41: 1–6.
- 96. During H, Loveys BR (1982) Diurnal changes in water relations and abscisic acid in field-grown Vitis vinifera cultivars. Vitis 21: 223–232.
- 97. Pérez-Leroux HAJ, Long SP (1994) Growth analysis of contrasting cultivars of Zea mays L. at different rates of nitrogen supply. Annals of Botany 73: 507–513.
- 98. Dawson TP, Curran PJ, Plummer SE (1998) LIBERTY - modeling the effects of leaf biochemical concentration on reflectance spectra. Remote Sensing of Environment 65: 50–60.
- 99. Jacquemoud S, Baret F (1990) PROSPECT: A model of leaf optical properties spectra. Remote Sensing of Environment 34: 75–91.