|
|
|
Among these models, radiative transfer models have been successfully used in the foreward mode to calculate leaf reflectance and transmittance, and in inversion to estimate leaf biophysical properties. Up to the present, these models have been mainly used to estimate chlorophyll and/or water contents as input parameters (Jacquemoud and Baret, 1990; Fukshansky et al., 1991; Yamada and Fujimura, 1991; Martinez v. Remisowsky et al., 1992). The influence of protein, cellulose, lignin, and starch on leaf reflectance has been recently introduced by Conel et all (1993), who proposed a two-stream Kubelka-Munk model, but in fact, the estimation of leaf biochemistry from remote sensing remains an open question. In order to clarify it, a laboratory experiment associating visible / infrared spectra of plant leaves both with physical measurements and biochemical analyses was conducted at the joint Research Centre in Italy during the summer of 1993. Thousands of measurements were collected in a unique data set partially used to upgrade the PROSPECT model (Jacquemoud and Baret, 1990) by including leaf biochemistry. In an earlier article, Fourty et al. (1995) analyzed the optical properties of the dry leaves and decomposed the absorption spectra of dry vegetation into six specific absorption coefficients having the characteristics of protein, cellulose, hemicellulose, lignin, starch, and sugar. However, in inversion, the estimation of the biochemical constituents was poor.
In this article, we develop a general model applicable both to fresh
and dry leaves that include a wide range of internal cellular structures
and biochemical compositions. In the first part we describe the experiment
and the correlation relationships between the biochemicals. The construction
of PROSPECT is only described where improvements were made in the original
model; otherwise the reader is referred to Jacquemoud and Baret (1990).
The validation, that is, the comparison between the measured and estimated
leaf biophysical and biochemical properties, is a key section of this article,
which justifies the hypotheses of the model. Finally we perform a sensitivity
analysis, as a last step before applying the model at the canopy level
in the future.
Among the physical and biological measurements performed in the frame of LOPEX93, the blade thickness, the specific leaf area (1 / SLA = dry weight per unit leaf area), the equivalent water thickness (EWT = water mass per unit leaf area), the photosynthetic pigments (chlorophyll a, b and total carotenoids expressed in gm cm-2), some biochemical components (protein, cellulose, hemicellulose, lignin, and starch), and finally the elementary composition (C, H, 0, N) were available. In near-infrared reflectance spectroscopy (NIRS) studies, the leaf biochemistry is typically given as a percentage of dry weight. We believe that this unit is not suitable for our study because fractions do not represent the amount of matter interacting with light. This assertion needs some explanation. First, consider the composition of plant foliage. Water represents on average 66.4% of the fresh weight (Fig. la). The remaining part is dry matter composed of cellulose, hemicellulose, lignin, protein, starch, and minerals (Fig. lc), All these constituents explain 85.8% of the dry mass of monocotyledons and 67.8% of the dry mass of dicotyledons, The missing matter may be attributed to lipids, soluble sugars, amino acids, and other primary and secondary metabolites not measured in this study (Fourty et al., 1995). It is not surprising that terrestrial plants have such similar chemistries since they share basic metabolic pathways. On the other hand, the basis for differences between the two groups of flowering plants is unclear, although it may reflect their ecological differentiation and the more herbaceous nature of the monocotyledons chosen for this study. Details of the decomposition can be found in Table 1. Although the concentration of the carbon based constituents may vary, their global fractions are remarkably stable in accordance with the very stable concentration of carbon in plant leaves which averages 47 g g-1 of dry matter. That kind of low variance information is not very useful! Consequently, as for water and pigments, we expressed the other concentrations in mass per unit leaf area using the SLA. Figures 1b and 1e illustrate the increased variability expressed when concentration units are in g cm-2. On an area basis, the biochemical variation increases by a factor between 1 and 10.
Several correlative relationships among biochemicals were also established,
including leaf thickness and EWT, protein and SLA or total chlorophyll,
For instance, we showed that 1 / SLA varied inversely with the weight-based
measure of leaf protein, consistent with values in the literature (Field
and Mooney, 1986; Dijkstra, 1989). The strongest relationships were obtained
between total nitrogen (N) and protein, and between total carbon (C) and
cellulose + lignin when expressed in g cm-2 (Fig.
2). This equivalence is very important because it suggests that the
C / N ratio which drives the decomposition rates of forest litter, affecting
nutrient cycling and trace gas fluxes, could be replaced by the indirect
measure of the ratio of cellulose + lignin to protein. Thus, spectral measures
of water, photosynthetic pigments, total C. and total N could provide significant
information related to canopy nutrient status, physiological state, and
allocation of photosynthate to aboveground canopy components.
Modeling absorption processes first implies that the effects of mesophyll
structure are well accounted for by the model. These effects influence
the whole spectrum but are maximum in the NIR where the absorption features
of chlorophyll and water are minimal if not negligible, On fresh leaves,
this low absorptance is materialized by a plateau of constant reflectance
and trmsmittance values at about 10% incident light. This plateau is somewhat
disturbed in artificially dried leaves due to the appearance of brown pigments
or denatured proteins that absorb light shorter than 1100 nm. In the original
version of PROSPECT, the leaf optical properties in the NIR were only driven
by the parameter N2 the number of stacked elementary layers; the absorption
by one of these layers was small and assumed to be constant. Although the
origin of this absorption is uncertain, it cannot be attributed to either
chlorophyll or water, Likewise, no major leaf-soluble cell fractions have
absorptions across this wavelength interval. So we have hypothesized that
the absorption was due to some component in the cell walls. Since the amount
of dry matter per unit area varies from sample to sample, the absorption
was also allowed to vary. In consequences leaf optical properties in the
NIR are now modeled by the N parameter and by the absorption coefficient
k(l) of this elementary layer. To determine
N, we defined three wavelengths corresponding to the maximum reflectance
(lr), the maximum transmittance
(lt), and the minimum absorptance
(la). For fresh leaves, these
three wavelengths are located in the NIR plateau and may be confounded;
for drv leaves, they are shifted towards longer wavelengths (Fig.
3). The structure parameter of each leaf was adjusted at the same time
as the three absorption coefficients by minimizing:
|
|
(1)
|
The wavelength independent mesophyll structure parameter N allows
the inversion of the Stokes equations: using measured reflectance R(l)
and transmittance T(l), the optical
properties of the compact layer (N = 1) are easily calculated for
each leaf, permitting the determination of a spectral absorption coefficient
k0(l). If the assumption
is made that the leaf is a homogeneous mixture of biochemical components,
this coefficient can be written as
|
|
(2)
|
The first strategy consists in predicting the constituent concentrations and, then, in comparing predictions with measured values. This approach presumes that the specific absorption coefficients we known, for example, deduced from optical measurements performed on pure substances. Thus, the spectral specific infrared absorption coefficients for distilled water have been carefully measured by Curcio and Petty (1951), The data for chlorophylls and, to a lesser extent, for accessory pigments (carotenoids, xantophylls) are also available from the literature (Lichtenthaler, 1987). In this case, though, the use of these absorption coefficients presents a problem: The absorption spectrum obtained from extracts of chlorophyll in a solvent does not correspond closely to the in vivo measurement; spectral shifts of the order of 10 nm are observed (Buschmann and Nagel, 1991; Chappelle et al., 1992), which are attributed both to the influence of the solvent and to the fact that the chlorophylls inside leaf tissues are complexed with other pigments and proteins. The complex three-dimensional macromolecular structure of chlorophyll has been associated with various distortions that contribute to wavelength spreading and shifting of the in vivo optical properties of pigments. For different reasons, the spectral signatures of the other biochemicals are also complex: even if some substances we composed of well characterized repeating units (e.g., starch, sugar), molecular weights can vary, while others are families of biochemical substances which cannot be precisely defined or even isolated with the molecular structure intact (e.g., protein, cellulose, lignin). Moreover, these large classes of macromolecules contain many chemical bonds in common (C-H, N-H, C-0, 0-H, etc.) which occur in various proportions, inducing an overlapping variation in absorption features (Barton et al., 1992).
In order to bypass these difficulties, another strategy was adopted: Using the absorption coefficients k0(l) and the measured concentrations, we deduced the specific absorption coefficients of leaf biochemical components ki(l). Different combinations of leaf biochemical composition and leaf water status have been tested. For instance, in Eq. (2) we decomposed the absorption into chlorophyll a + b, water, protein, and structural biochemicals like cellulose, hemicellulose, or lignin. We also investigated coefficients determined on fresh leaves, dry intact leaves, and fresh + dry leaves, It allowed us to address a delicate problem experienced by those using NIRS techniques to estimate leaf biochemistry. Generally, a regression equation is established between the leaf biochemistry and the optical properties of entire blades or dry vegetation powders. For a given component, Jacquemoud et al. (1959) showed that the wavelengths selected by multiple stepwise regression analysis depended on whether the basis for comparisons was reflectance or transmittance values on fresh/dry single leaves, or leaf stacks. This discrepancy contradicts the idea that a specific biochemical should produce a consistent effect due to its absorption of light at specific allowable energy states. Nonetheless, many statistical analyses have resulted in selections of significantly different wavelengths and the need for taxon specific relationships to particular biochemicals. Let us consider the biochemical coefficients derived in this analysis. The coefficients for water, protein and cellulose + lignin over the 800-2500 nm range are shown in Figure 4. One can see that protein and cellulose + lignin coefficients occur at the same wavelength position whatever the leaf water status. Even where water tends to mask the absorption peaks of these constituents, as expected in fresh leaves, it does not fundamentally shift them to other wavelengths in our analysis. This important result allows us to build a very general model suitable for many kind of flowering plant leaves.
We calculated the specific absorption coefficients of leaf biochemicals
for the following three combinations: chlorophyll, water, protein, cellulose
+ hemicellulose, and lignin [C1], chlorophyll, water, protein, and cellulose
+ hemicellulose + lignin [C2], and chlorophyll, water, protein, and cellulose
+ lignin [C3]. Kab was determined on fresh leaves in
the 400-800 nm region, Kw on fresh +dry leaves in the
800-2500 nm region, and the other coefficients on dry leaves in the 800-2500
nm region. To describe these results, lets examine the case for C3: Kab(l)
displays classical features of photosynthetic pigments with spectral shifts
toward longer wavelengths of the principal absorption peaks of chlorophyll
compared to in vitro observations as discussed earlier (Fig.
5a). Kw(l) shows
good agreement with the fundamental constants published for pure liquid
water (Fig. 5b). Results we surprisingly
also very convincing for both protein (Kp) and cellulose
+ lignin (Kc1): except in a few cases, the absorption
peaks are well represented (Figs. 5c and 5d).
Cl and C2 produced very similar results. The specific absorption coefficients
of chlorophyll a+ b and water equal zero, respectively, after and before
800 nm; due to the appearance of brown pigments or denatured protein during
the drying of the leaves, those of the other constituents have been fixed
to the value calculated at 1100 nm. This assumption is reasonable since
the optical properties of an albino leaf devoid of pigments we constant
along the visible and the NIR wavelength region.
The validation was carried out with the same data sets. It consisted
in estimating the model parameters symbolized by the vector q
by minimizing the following merit function:
|
|
(3)
|
In this sensitivity analysis, the leaf biochemicals were assumed to
be physiologically independent, and concentrations were varied independently.
In reality, physiological functions limit the range of variation that can
be tolerated and still maintain foliar integrity (Field et al., 1992).
For instance, protein and chlorophyll are often positively correlated in
leaves (Field and Mooney, 1986). This means that an increase in protein
content given some time for metabolic adjustment-should indirectly induce
change in leaf optical properties in the visible that results in a greener
leaf. Some very high correlations have been established on particular species
under controled conditions, for instance, corn (Ercoli et al., 1993); but
when a large number of species (e.g., more than fifty in this study) are
compared, and when species exhibit a wide range of adaptive growth patterns
(evergreen, deciduous, crops, C3, C4, wildland species), these correlations
are lower. Another artifact may be produced when varying the water content.
As shown by Woolley (1973), a wilting plant leaf tends to lose some water
which is replaced by air spaces. In consequence, two different optical
effects are combined: a decrease in the absorption (Cw)
and an increase in the multiple scattering (N). Simulations of Figure
6e only show a variation of reflectance in the middle infrared due
to a variation of Cw while the effect on N was
not modeled. One could imagine that other examples of intercorrelations
between the model variables could be identified, for instance, between
water and pigments. Because of the absence of information about their correct
relationships, it is very difficult to take them properly into account.
Buschmann, C., and Nagel, E. (1991), Reflection spectra of terrestrial vegetation as influenced by pigment-protein complexes and the internal optics of the leaf tissue, in Proc. Int. Geosci. Remote Sens. Symp., Helsinki, Finland, pp. 1909-1912.
Chappelle, E. W., Kim, M. S., and McMurtrey, J. E., III (1992), Ratio analysis of reflectance spectra (RARS): an algorithm for the remote estimation of the concentrations of chlorophyll A, chlorophyll B, and carotenoids in soybean leaves, Remote Sens. Environ. 39:239-247.
Conel, J. E., van den Bosch, J., and Grove, C. I. (1993), Application of a two-stream radiative transfer model for leaf lignin and cellulose concentrations from spectral reflectance measurements. Parts 1 & 2, in Proc. 4th Annual JPL Airborne Geoscience Workshop. Vol. 1. AVIRIS Workshop (R. 0. Green, Ed.), 25-29 October, Washington, DC, NASA-JPL Publication 93-26, Pasadena, CA, pp. 39-51.
Curcio, J. A., and Petty, C. C. (1951). The near infrared absorption spectrum of liquid water, J. Opt. Soc. Am., 41(5): 302-304.
Demmig-Adams, B., Winter, K., Kruger, A., and Czygan, F.C. (1989), Light stress and photoprotection related to the carotenoid zeaxanthin in higher plants, in Photosynthesis (W. R. Briggs, ed.), Liss, New York, pp. 375-391.
Dijkstra, P. (1989), Cause and effect of differences in specific leaf area, in Causes and Cosequences of Variation in Growth Rate and Productivity of Higher Plants (H. Lambers et a]., Eds.) Academic, The Hague, pp. 125-140.
Ercoli, L., Mwiotti, M., Masoni, A., and Massantini, F. (1993). Relationship between nitrogen and chlorophyll content and spectral properties in maize leaves, Eur. J. Agron. 2(2): 113-117,
Field, C. B., and Mooney, H. A. (1986), The photosynthesis nitrogen relationship in wild plants, in On the Economy of Plant Form and Function (T. J. Givnisb, Ed.), Cambridge University Press, Cambridge, pp. 25-55.
Field, C. B., Chapin, F. S., III, Matson, P. A., and Mooney, H. A. (1992), Responses of terrestrial ecosystems to the changing atmosphere: a resource-based approach, Ann. Rev. Ecol. Syst. 23:201-236.
Fourty, Th., Buet, F., Jacquemoud, S., Schmuck, C., and Verdebout, J. (1996), Optical properties of dry plant leaves with explicate description of their biochemical composition: direct and inverse problems, Remote Sens. Environ. 56: 104-1174
Fukshansky, L., Fukshansky-Kazarinova, N., and Martinez v. Remisowsky, A. (1991), Estimation of optical parameters in a living tissue by solving the inverse problem to the multiflux radiative transfer, Appl. Opt. 30(22):3145-3153.
Hosgood, B., Jacquemoud, S., Andreoli, G., Verdebout, J., Pedrini, G., and Schmuck, G. (1995), Leaf Optical Properties EXperiment 93 (LOPEX93) Report EUR-16095-EN, European Commission, Joint Research Centre, Institute for Remote Sensing Applications, Ispra, Italy
Jacquemoud, S., and Baret, F (1990), PROSPECT: a model of leaf optical properties spectra, Remote Sens. Environ. 34:75-914
Jacquemoud, S., Verdebout, J., Schmuck, G., Andreoli, G., and Hosgood, B. (1995), Investigation of leaf biochemistry by statistics, Remote Sens. Environ. 54:180-188.
Liebtenthaler, H. K. (1987), Chlorophylls and carotenoids: pigments of photosynthetic biomembranes, Methods Enzymol. 148:350-382.
Martinez v. Remisowsky, A., McClendon, J. E., and Fukshansky, L. (1992), Estimation of the optical parameters and light gradients in leaves: multi-flux versus two-flux treatment, Photochem. Photobiol. 55(6):857-865.
Peterson, D. L., and Hubbard, G. S. (1992), Scientific issues and potential remote-sensing requirements for plant biochemical content. J. Imaging Sci. Technol. 36(5):446-456.
Verdebout, J., Jacquemoud, S., and Schmuck, G. (1994), Optia properties of leaves: modelling and experimental studies, in Imaging Spectrometry a Tool for Environmental Observations (J. Hill and J. Megier, Eds.), ECSC, EEC, EAEC, Brussels and Luxembourg, pp. 169-191.
Woolley, J. T. (1973), Change of leaf dimensions and air volume with change in water content, Plant Physiol. 41: 815-816.
Yamada, N., and Fujimura, S. (1991), Nondestructive measurement of chlorophyll pigment content in plant leaves from three-color reflectance and transmittance, Appl. Opt. 30 (27):3964-3973.