Internal Reflectance Modelling of Hordeum vulgare Leaves During Drying
Spectral reflectance, or indexes that characterize spectral reflectance at concrete wavelength, is commonly used as an indicator of plant stress, or its photosynthetic apparatus status. In this paper, new leaf optical model is presented. Within this paper, experimental determination of surface and internal reflectance of Spring barley leaves and mathematical-physical modelling of internal reflectance were performed. It was proven that a new proposed theoretical model and the experimental spectra of internal reflectance are strongly correlated. It can be concluded that the total reflectance is not a function of epidermis condition, but it testifies about overall functional condition of Spring barley leaves.
Leaf optical properties are associated with physiological stress, that affects these plants. Stress factors can be different: dehydration, freezing, ozone, diseases, herbicides, intraspecies and interspecies competition, insect, lack of mineral nutrition, high salinity, extreme temperatures, and so forth [1–3]. Reflectance spectrum changes are similar for different plant species . Because of this fact, spectral reflectance, or indexes that characterize spectral reflectance at concrete wavelengths, is commonly used such as reliable indicators of plant condition and/or their photosynthetic apparatus status. The study of spectral characteristics has not only scientific, agricultural, and environmental reasons but it is also important from the economical point of view. The spectral characteristics are predominantly used for detection and prevention of potential stress factors [3, 4].
When the electromagnetic radiation impacts air—leaf interface, three competition processes may occur (absorption, reflection, and transmission). Mathematically, it can be described using (1) as a where is electromagnetic flux that impacts the sample, is wavelength, is reflected radiation to the upper halfspace, is angle of incidence, and is leaf surface anisotropy function. is absorbed radiation and is radiation, which passed through the leaf to the lower halfspace. Let us assume that the radiation is not leaking by leaf sides and the leaf fluorescence is not counted into the energy balance described in (1). If the following variables are assumed: reflectance , absorbance , and transmittance , (1) can be rewritten in the form of Total light intensity reflected from the leaf surface is formed by two components, surface reflectance and internal reflectance . Reflectance (surface and internal) is a function of a wavelength and angle of impact if it is measured considering normal . Surface reflectance can be divided into the two components: mirror reflectance and diffuse reflectance . Mirror reflectance is typical for homogeneous objects with perfectly smooth surface. For ideal mirror-like surfaces, it holds that the angle of incidence is equal to the angle of reflection. For surface diffuse reflectance, it is typical that the light is reflected randomly due to rough surface structures . Leaf internal reflectance is always diffuse.
In fact, there are no surfaces that can be characterized as strictly mirror-like or diffuse. Surface reflectance has always two components, mirror-like and diffuse . If the simplification is used, that leaf surface reflectance is mirror-like, Fresnel equations can be used for surface reflectance components calculation. Fresnel equations represent relations among vector amplitudes of incident, reflected, and refracted waves. Quantity of refracted transverse electrical wave is defined as where and are angles of incidence and reflection. Transverse reflected magnetic wave quantity can be defined as
Due to the fact that the sunlight is nonpolarized electromagnetic radiation, neither transverse reflected electric nor magnetic waves are equal zero. Surface reflectance can be expressed through the squares of relevant amplitude coefficients
The light, that is not reflected to the upper halfspace, enters to the leaf, according to the Snell law: where , respectively, and are leaf refractive indexes. With respect to the energy conservation law, light intensity for radiation that passes through the air-leaf interface can be defined as
It is necessary to define leaf surface anisotropy factor for calculation of diffuse surface reflectance. The anisotropy factor is the input parameter for calculation of Henyey-Greenstein probability function . Let us assume that the light beam is reflected to the upper halfspace, then the following can be used:
2. Materials and Methods
Spectrophotometric measurements were realized by double beam scanning spectrophotometer UV-VIS 550 UNICAM within 400–800 nm wavelength intervals. Diffuse reflectance measurements were performed using integral spheroid RSA-UC-40 Labsphere, which was incorporated into the spectrophotometer. Microscopic measurements were realized by light microscope Amplival. All reflectance measurements were performed on the adaxial leaf side. Primary leaf blades of Spring barley (Hordeum vulgare, cv. Akcent) were used for all measurements. The plants were growing in the growth chamber with the following parameters: temperature 21°C, light regime: 16 hours of light, 8 hours of darkness, and the light source was fluorescence tubes with light intensity 90 μmol·m−2·s−1. According to Feekes phenomenological scale, barley leaves can be characterized by growing phase 1,3. Firstly, the reference group that consisted of nonstressed plants (no stress) was measure; then the plants stressed by drought were measured, too (3rd, 6th, and 9th day). Diffuse reflectance spectra were achieved by the following process. At first, diffuse reflectance was measured separately for whole leaf and epidermis. For epidermis reflectance measurements were used strongly absorbing background to minimize light reflected at epidermis-background interface. Internal reflectance spectrum was calculated as a difference between total diffuse reflectance and the whole epidermis reflectance. Pigments content was calculated according to equations proposed by Lichtenthaler  as an average of six samples. To achieve more trustworthy model, rough measurements of cells size and shape were realized by light microscope. Also the size, shape, and number of chloroplasts in the cells were studied. In presented model, cube shape of cells and spherical shape of chloroplasts are assumed. For microscopic measurements, the control groups of nonstressed and barley stressed plants for nine days were used.
3. Results and Discussion
The presented model is based on few presumptions. The cell walls have transmittance and reflectance both of interval ; light during its propagation passes through chloroplasts; in chloroplasts the Lambert-Beer law holds (it is presumed that chlorophyll content in chloroplasts is constant). Chloroplasts are the only one absorbing components in the intracellular medium. Absorption of nonchloroplasts structures is neglected.
In the simplest case, when the light penetrates into the cell, reflects from the wall, and passes chloroplasts and after the next reflection it leaves the cell; internal reflectance can be defined by where represents absorption properties of one or more chloroplasts, is the molar extinction coefficient, is the chlorophyll concentration in chloroplasts, and is the optical density of one chloroplast. A more sophisticated case is represented by multiple scattering in one cell with one chloroplast. Those numbers may be theoretically infinite. For the light that passes through the cell wall and after reflection it passes through the chloroplast; the following relation can be used: Thereafter, the same light beam reflects from the wall and leaves the cell; in this case intensity of internally reflected light is: Using the same procedure, other wall reflections and chloroplast passes are: For cell reflections, the sum of geometrical series is realized: There are some extreme situations that may occur. If the (nontransparent) walls, than it is clear that . The second extreme situation occurs when , but (nonreflecting walls); in this case internal reflectance is also equal zero. The third extreme situation occurs when ; in this case plant tissue does not absorb: If the last term is changed in comparison with the previous case, , especially if converges to zero (strong absorption), then also .
Generally, if the light beam passes through cell walls and reflects on interfaces, passes through cells and in every of those cells there is chloroplasts, relation for calculation is modified to After approximation, when is set to zero, the relation for calculation is simplified to It is nothing else than an equation for a simple case (two passes, two reflections, and chloroplasts).
Within the discussed model, it is considered that light beam passes through the cell (); light beam geometrical length is . In the cell with the volume , there are chloroplasts. Another presumptions are chloroplasts spherical shape and cells cubic shape. Chloroplasts radius is and the diameter is . Chloroplast geometrical cross-section is (main cut):
Chloroplasts concentration in cell is ; said in different way, number of chloroplasts in cell volume unit can be expressed as
Let’s assume light beam surrounding, if the chloroplast center is inside the cylindrical volume: So if the distance between chloroplast center and light beam is smaller than , the light beam passes through the chloroplast. In the other case, when the distance between the light beam and chloroplast center is bigger than , the light beam does not pass through the chloroplast. Number of the chloroplasts with the centers inside the can be expressed as On the light beam length , there are chloroplasts, while the mean distance among chloroplasts is ; if we assume their regular arrangement, can be defined by or equivalently according to based on calculation manner. Equation (26) holds in case when the light beam hits the chloroplast after it has passed length equal to the average chloroplasts distance. On the contrary, (27) holds in the case when the light beam hits the chloroplast without any preceding pass . For simplicity, the first case is mentioned in the following text. If the spherical shape of the chloroplast (with radius and diameter ) is presumed, refractive indexes in the chloroplasts and its surrounding are the same and the light beams have the same optical density (intensity of incident collimated radiation is homogeneous); then average optical length, that is passed by the light beam in the spherical chloroplast, can be determined. Let us divide chloroplasts into the cylindrical sections (with radius ), where . Optical length in cylindrical section is defined by The mean optical length is or alternatively, after substitution, The mean optical length can be expressed as It is clear that average optical length in spherical chloroplast is . According to Poisson probability distribution, on the light trajectory , there are chloroplasts which can be defined as General formula for internal reflectance calculation is expressed by where indexes and define the number of the light that passes through and number of reflections realized on the cell walls. With respect to (11), (33) can be redefined to the following form: after modification of The formula that fundamentally includes all parameters that are necessary for internal reflectance theoretical spectra modelling is achieved: Some of these parameters can be measured independently, the others may be achieved using fitting procedures. The relation expressed by (36) can be simplified to
There are two unknown variables remaining (chlorophylls concentration in chloroplasts and light beam geometrical length); then if one of these parameters is determined properly, (see following text) the last unknown parameter can be calculated by modification of (37). Let us make a substitution according to After which modification of (37) can be expressed via
Formula for chlorophylls concentration calculation (for known light beam geometrical length) is defined as
On the contrary, if the chlorophylls concentration is known, the light beam geometrical length can be calculated after modification: Then can be expressed as
Another procedure, alternative to the method mentioned above, is a presumption that light beam hits chloroplasts, when it passes through the cells which can be defined as where means the sum of all partial probabilities, while it holds that . is the probability that the light beam does not hit any chloroplast during its pass through the leaf tissue, is the probability that the light beam hits exactly one chloroplast during its pass through the leaf tissue, and so on. If all necessary parameters, such as cells and chloroplasts size and shape, cell walls transmittance and reflectance, chlorophylls concentration in chloroplasts, and the distance, which have to be passed by light beam in the cell, are known, the formula for internal reflectance calculation can be expressed by With respect to the fact that the main contributors to the internal reflectance are especially light beams that pass only through a few chloroplasts, (44) can be rewritten in the form of progression: The other progression members are not mentioned anymore due to the fact that their contribution can be neglected. For determination of average distance among chloroplasts, the following procedure can be used, assuming that chloroplasts distribution in cells is homogeneous. The mean chloroplasts distance was calculated using Monte Carlo method (MATLAB 7.1). For every chloroplast, three pseudorandom numbers were generated, for interval . Meanwhile, chloroplasts geometrical centers in Cartesian coordinates system are defined according to where , and are pseudorandom numbers and is size. Due to the fact that within the presented model there are only chloroplasts mentioned with the volume completely located in the cell, the following criterion was used. Chloroplasts must accomplish boundary conditions expressed by where , and are the distances from the origin of coordinates system, and is the chloroplast radius. Another limiting criterion is prohibited overlap of the chloroplasts: where is the distance between two chloroplasts in the cell. The same procedure was repeated until all chloroplasts center positions were successfully simulated. In experiment, it was investigated that . In general, if two chloroplasts are mentioned, namely, chloroplasts and , where , , indexes Cartesian coordinates, then their mutual position can be calculated from
Chloroplasts distance is not defined by the distance of their geometrical centers but by the distance of two nearest mutual points. The mean distance among all chloroplasts present in the cell is ; this is expressed by When all necessary calculations were performed, the following results were achieved. The mean distance among chloroplasts in Spring barley cells, which was not stressed by dryness, was μm. The barley stressed by drought for 9 days is μm. For the other calculations according to spectrum of pigments molar extinction coefficients (Figure 1) was used.
Using optical microscope, it was determined that the number of chloroplasts was the same for both groups (15). There was also difference in chloroplasts sizes. For the control group, the average chloroplasts diameter was 3,6 μm; meanwhile, in case of stress plants it was only 2,1 μm. There were also differences in cells sizes. The barley stressed by drying had approximately 20% smaller cells diameter ( μm for stressed plants and μm for control group).
If the chlorophylls content per surface area is known, it is possible to calculate rough chlorophylls concentration in chloroplasts. If we assume that during the stress influencing period cells and chloroplasts size are changing linearly, but chloroplasts number remains constant, then the chlorophylls concentration can be calculated for each day of stress conditions: where (molL−1) is chlorophylls concentration in chloroplasts for concrete day, (μm) is chloroplasts radius that was not stressed by drought, (μm) is the radius of the chloroplasts that was stressed by dryness, for days, (μm) is the barley cell that was not stressed by drought and (μm) is the size of the cell that was stressed by dryness for days. If the chlorophylls concetration in chloroplasts is known, it is possible to determine the probability of the light beam pass through the chloroplasts for different wavelengths:
Based on the knowledge of molar extinction coefficients for pigments and spectrophotometry measurements of pigments content probability of light beam passes through the chloroplasts for wavelengths 550 nm and 680 nm was calculated; see Figures 2(a) and 2(b).
It is necessary to do some simplifications in case of the internal reflectance theoretical spectra modelling. The first of them is precondition, in which cell walls transmittance and reflectance are not a function of wavelength (at least for concrete day of measurement). Chlorophylls concetration in chloroplasts, pigments molar extinction coefficients, and cell walls transmittance and reflectance are input parameters for (44). In case of cells and chloroplasts size, it is presumed that these parameters are changing linearly during the dryness. In the next step, the mean optical length, number of cell walls passes, and reflections are optimized; see Tables 1 and 2.
Good agreement between the experimental internal reflectance spectrum and proposed theoretical model was achieved; see Figure 3.
Experimentally measured spectra of diffuse reflectance may be quantitatively very different compared to theoretical spectra. It can happen due to the total reflection of some light beams on the cell walls. The cells do not have to be close to each other, so in fact, total reflection may occur when the light passes from the optically more dense environment (palisade cells) to the optically less dense environment (air in the intercell space).
If the experimentally measured spectra of the internal reflectance are compared with the theoretical spectra calculated according to proposed model, good agreement is achieved. This good agreement is achieved both for control group and plants stressed for 9 days by dryness. Theoretical spectra exhibit interesting informations. Firstly, it is useful to divide the reflectance spectrum into three sections according to the wavelength interval (400–500 nm; 500–640 nm and 640–800 nm). Wavelength interval 400–500 nm is described by relative strong chlorophylls and carotenoids absorption; the mean optical length of the light beam is significantly shorter than for the light with wavelength from interval 500–640 nm; this interval may be characterized by a relative weak pigments absorption. The shortest light mean optical length was calculated for the wavelength 680 nm due to the strong Chl a absorption. The rapid increase of the internal reflectance for wavelengths longer than 700 nm is mainly caused by water content decrease during the drought. The changes in the internal reflectance spectrum are not influenced only by cells and chloroplasts number, size, shape, and pigments content. Very important is also total internal reflection contribution. Weight of the total reflection increases with decreasing water content in the intracellular space.
Experimentally determined spectra of the internal reflectance of Spring barley leafs were compared with proposed theoretical model. Good agreement between theoretical model and experimental data was observed. It was revealed that most influencing factors for nonstressed leafs internal reflectance spectra modelling are cells and chloroplasts number, size and shape, pigments content, cell walls transmittance, and reflectance. However, for the Spring barley stressed by drought, water content was critical. Although water absorption can be neglected in visible light spectrum, intercellular spaces fulfilled by air are responsible for the total reflection, especially for the light with wavelength longer than 700 nm. However, internal reflectance is two orders lower than surface reflectance; it has significant information value about plant physiological state.
This paper has been elaborated in the framework of the project opportunity for young researchers, reg. no. CZ.1.07/2.3.00/30.0016, supported by Operational Programme Education for Competitiveness and cofinanced by the European Social Fund and the state budget of the Czech Republic. The authors are grateful to the Ministry of Education of the Czech Republic for the financial support of the project of OP VaVPI “Regional Materials Science and Technology Centre (RMSTC)” no. CZ.1.05/2.1.00/01.0040.
G. A. Carter, “Responses of leaf spectral reflectance to plant stress,” American Journal of Botany, vol. 80, pp. 239–243, 1993.View at: Google Scholar
G. A. Carter and A. K. Knapp, “Leaf optical properties in higher plants: linking spectral chracteristics to stress and chlorophyll concentration,” American Journal of Botany, vol. 88, pp. 667–684, 2001.View at: Google Scholar
G. A. Carter, W. G. Cibula, and R. L. Miller, “Narrow-band reflectance imagery compared with thermal imagery for early detection of plant stress,” Journal of Plant Physiology, vol. 148, no. 5, pp. 515–522, 1996.View at: Google Scholar
G. V. G. Baranoski and J. G. Rokne, “Efficiently simulating scattering of light by leaves,” Visual Computer, vol. 17, no. 8, pp. 491–505, 2001.View at: Google Scholar
K. Koutajoki, “Bidirectional surface scattering distribution function,” Computer Science, vol. 12, pp. 258–2271, 2002.View at: Google Scholar
V. C. Vanderbilt and L. Grant, “Polarization of light by vegetation,” Computer Science, vol. 8, pp. 111–126, 2000.View at: Google Scholar
L. G. Henyey and J. L. Greenstein, “Diffuse radiation in the galaxy,” Astrophysical Journal, vol. 93, pp. 70–83, 1941.View at: Google Scholar