Next Article in Journal
Measurement Characteristics of Near-Surface Currents from Ultra-Thin Drifters, Drogued Drifters, and HF Radar
Next Article in Special Issue
Application of a Three-Dimensional Radiative Transfer Model to Retrieve the Species Composition of a Mixed Forest Stand from Canopy Reflected Radiation
Previous Article in Journal
Hyperspectral Image Restoration under Complex Multi-Band Noises
Previous Article in Special Issue
Implications of Whole-Disc DSCOVR EPIC Spectral Observations for Estimating Earth’s Spectral Reflectivity Based on Low-Earth-Orbiting and Geostationary Observations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Leaf Specular Reflection on Canopy Radiative Regime Using an Improved Version of the Stochastic Radiative Transfer Model

1
College of Electrical and Information Engineering, Hunan University, Changsha 410082, China
2
Department of Earth and Environment, Boston University, 685 Commonwealth Avenue, Boston, MA 02215, USA
3
State Key Laboratory of Remote Sensing Science, School of Geography, Beijing Normal University, Beijing 100875, China
4
Guangxi Colleges and Universities Key Laboratory of Unmanned Aerial Vehicle (UAV) Remote Sensing, Guilin University of Aerospace Technology, Guilin 541004, China
5
Changchun Institute of Optics, Fine Mechanics and Physics, Chinese Academy of Sciences, Changchun 130033, China
6
Guangxi Zhuang Autonomous Region Environmental Monitoring Centre, Nanning 530028, China
*
Authors to whom correspondence should be addressed.
Remote Sens. 2018, 10(10), 1632; https://doi.org/10.3390/rs10101632
Submission received: 9 July 2018 / Revised: 27 September 2018 / Accepted: 11 October 2018 / Published: 14 October 2018
(This article belongs to the Special Issue Radiative Transfer Modelling and Applications in Remote Sensing)

Abstract

:
Interpreting remotely-sensed data requires realistic, but simple, models of radiative transfer that occurs within a vegetation canopy. In this paper, an improved version of the stochastic radiative transfer model (SRTM) is proposed by assuming that all photons that have not been specularly reflected enter the leaf interior. The contribution of leaf specular reflection is considered by modifying leaf scattering phase function using Fresnel reflectance. The canopy bidirectional reflectance factor (BRF) estimated from this model is evaluated through comparisons with field-measured maize BRF. The result shows that accounting for leaf specular reflection can provide better performance than that when leaf specular reflection is neglected over a wide range of view zenith angles. The improved version of the SRTM is further adopted to investigate the influence of leaf specular reflection on the canopy radiative regime, with emphases on vertical profiles of mean radiation flux density, canopy absorptance, BRF, and normalized difference vegetation index (NDVI). It is demonstrated that accounting for leaf specular reflection can increase leaf albedo, which consequently increases canopy mean upward/downward mean radiation flux density and canopy nadir BRF and decreases canopy absorptance and canopy nadir NDVI when leaf angles are spherically distributed. The influence is greater for downward/upward radiation flux densities and canopy nadir BRF than that for canopy absorptance and NDVI. The results provide knowledge of leaf specular reflection and canopy radiative regime, and are helpful for forward reflectance simulations and backward inversions. Moreover, polarization measurements are suggested for studies of leaf specular reflection, as leaf specular reflection is closely related to the canopy polarization.

1. Introduction

Solar radiation intercepted by vegetation canopies undergoes various physical and biological processes until it is either absorbed or escapes from the canopy boundaries. Monitoring these processes has been one of the biggest objectives of vegetation remote sensing. Radiative transfer (RT) theory bridges remote sensing observations and these processes that generate signals and, thus, is important for research regarding vegetation remote sensing [1,2,3,4].
Numerous radiative transfer models have been developed based on the RT theory. They can be broadly categorized as 1D [5,6,7], 3D [8,9,10], and stochastic [11,12] models. The 1D model is simple but less realistic, whereas the 3D model is realistic but complicated. Given this, the stochastic radiative transfer model (SRTM) was developed. In the SRTM, 3D canopy structure is considered using canopy pair-correlation function, which is defined as the probability of simultaneously finding phytoelements at two points [11]. It has been shown that the SRTM integrates advantages of 3D and 1D RT models; it is as realistic as the 3D RT model and as simple as the 1D RT model [11]. It has also been successfully adopted in generating leaf area index (LAI) and its sunlit portion from data collected from the EPIC instrument in the Deep Space Climate Observatory (DSCOVR) mission [13].
Like the majority of the widely used RT models (e.g., the SAIL model [14]), the SRTM assumes that leaves are ideal bi-Lambertian scatters [11,12], i.e., the incident radiation scattered by a leaf follows a cosine distribution about the leaf normal. Such assumption can reduce the number of model parameters and, consequently, simplify the photon-vegetation interactions. However, this is not the case in reality. Scattering from a leaf consists of two parts, diffuse and specular [15,16]. The former is primarily generated from photon interactions with the leaf interior. This part follows a near bi-Lambertian distribution [17]. The latter results from radiation speculary reflected at the air-cuticle boundary. Unlike the diffuse part, the specular part exhibits a strong dependence on sun-sensor geometry. It varies with species and can be up to 50% in the visible bands [18]. Individual leaves are, thus, non-Lambertian scatters.
There are certain attempts to incorporate leaf specular reflection into RT modeling. The general idea is to modify the leaf scattering phase function by accounting for leaf specular reflection, i.e., both diffuse and specular parts are considered [17,19,20,21]. By doing so, the specularly reflected radiation was shown to change the amount of radiation registered by the sensor, which is usually normalized and parameterized as the canopy bidirectional reflectance factor (BRF, the ratio of surface leaving radiance and radiance from an ideal Lambertian surface under identical illumination conditions [22]), especially at strongly absorbing bands. Additionally, few efforts have been devoted to fully exploring the influence of leaf specular reflection on canopy radiative regime. Compared to BRF, the canopy radiative regime provides more comprehensive information about leaf specular reflection and its influences. Knowledge of leaf specular reflection and the canopy radiative regime will substantially help to improve forward canopy reflectance simulations, as well as backward inversions.
The goal of this paper is to propose an improved version of the SRTM by accounting for leaf specular reflection and exploring the influence of leaf specular reflection on the canopy radiative regime. In particular, we investigate its influences on vertical profiles of mean radiation flux density, canopy absorptance, BRF, and normalized difference vegetation index (NDVI, the ratio between the difference and sum of the canopy BRF in NIR and red bands).
The paper is organized as follows. A description of the improved version of the SRTM is presented in Section 2. The evaluation of the model is presented in Section 3. The analyses and discussion of the improved SRTM and influence of leaf specular reflection on the canopy radiative regime are given in Section 4 and Section 5, respectively. Finally, Section 6 summarizes the results.

2. Methods

2.1. Stochastic Radiative Transfer Model

The stochastic radiative transfer model provides the most linkage among remotely measured mean intensity, leaf optical property, and canopy 3D structure. For a stochastic canopy medium that is confined to 0 z H , the horizontal mean intensity I ¯ ( z , Ω ) for downward ( μ < 0 ) and upward ( μ > 0 ) directions at depth z can be descripted as the following equations [11,12]:
I ¯ ( z , Ω ) = 1 | μ | 0 z σ ( Ω ) p ( ξ ) U ( ξ , Ω ) d ξ + 1 | μ | 0 z p ( ξ ) S ( ξ , Ω ) d ξ + I ¯ 0 ( Ω ) , μ < 0
I ¯ ( z , Ω ) = 1 | μ | z H σ ( Ω ) p ( ξ ) U ( ξ , Ω ) d ξ + 1 | μ | z H p ( ξ ) S ( ξ , Ω ) d ξ + I ¯ H ( Ω ) , μ > 0
Here, μ = cos θ and θ is the polar angle of direction Ω . σ ( Ω ) is the extinction coefficient. p ( ξ ) is the probability of finding a foliated point at depth ξ . The second moment U ( ξ , Ω ) represents the mean intensity that incidents on the leaf surface at depth ξ along direction Ω . The term S ( ξ , Ω ) denotes the spherical integration of scattering as:
S ( ξ , Ω ) = 4 π σ s ( Ω Ω ) U ( ξ , Ω ) d Ω
where, σ s ( Ω Ω ) is the differential scattering coefficient, which can be represented as u L Γ ( Ω Ω ) / π . u L is the leaf area volume density. Γ ( Ω Ω ) is the area scattering phase function (or Γ -function). I ¯ 0 ( Ω ) and I ¯ H ( Ω ) describe the upper and lower boundary conditions.
Note that in the classical version of the SRTM, Γ ( Ω Ω ) is parametrized under the assumption that leaves are ideal bi-Lambertian scatters, i.e., leaf specular reflection is neglected [11,13]. In this case, the Γ -function can be expressed as:
1 π Γ ( Ω Ω ) = 1 2 π 2 π + g L ( Ω L ) | Ω · Ω L | γ L ( Ω L , Ω Ω ) d Ω L
Here, g L represents the leaf normal distribution function; γ L denotes leaf scattering phase function; Ω L , Ω , and Ω are directions of the leaf normal, incident, and scattered radiation, respectively.
Equations (1)–(4) describe the key idea of the classical version of the SRTM. It should be noted that the most straightforward variable wherein leaf specular reflection matters is the leaf scattering phase function. Our next step is to modify this variable by accounting for leaf specular reflection.

2.2. Improvement of the SRTM by Accounting for Leaf Specular Reflection

Accounting for leaf specular reflection in the SRTM needs parameterization of the scattering process that appears at the leaf scale. Radiation scattered from a leaf consists of two parts, diffuse ( γ L D ) and specular ( γ L S ), as Figure 1a shows. Leaf scattering phase function ( γ L ) can be expressed as:
γ L ( Ω L , Ω Ω ) = γ L D ( Ω L , Ω Ω ) + γ L S ( Ω L , Ω Ω )
Spherical integration of Equation (5) results in leaf albedo ω L = ( 1 s L ) ω L D + s L , where s L is the scattering part for leaf specular reflection, and ω L D is the scattering part for leaf diffuse scattering, given that photons interact with leaf internal constituents.
It follows from Equation (5), one gets the total Γ -function as:
Γ ( Ω Ω ) = Γ D ( Ω Ω ) + Γ S ( Ω Ω )
Here Γ D ( Ω Ω ) and Γ S ( Ω Ω ) are the Γ -functions for the diffuse and specular part, respectively.
The specular area scattering phase function, Γ S ( Ω Ω ) , describes the specularly reflected radiation resulting from photon interactions on the leaf surface. It can be evaluated as [23]:
1 π Γ S ( Ω Ω ) = 1 8 π Κ ( κ , α ) F r ( n , α ) g L ( Ω L )
where Κ is a reduction factor; F r is Fresnel reflectance describing the magnitude of specularly reflected radiation. In this model, Γ S depends on the index of refraction, n, a parameter κ that quantifies the hairy cuticular structures ( 0.1 0.3 ) [23], and the incident angle α ( = acos ( Ω · Ω L ) ) .
The reduction factor proposed by Nilson and Kuusk is incorporated [24], i.e.:
Κ ( κ , α ) = exp ( κ tan | α | )
It has the desirable property as Κ 1 for small incident angles and Κ 0 for oblique incident angles. For small incident angles, the leaf surface roughness is smaller and, therefore, specularly reflected radiation is less reduced (i.e., Κ 1 ). However, for oblique incident angles, the roughness is greater and consequently lowers specularly reflected radiation (i.e., Κ 0 ).
The Fresnel reflectance is given by [15]:
F r = 1 2 [ sin 2 ( α i ) sin 2 ( α + i ) + tan 2 ( α i ) tan 2 ( α + i ) ]
where the refractive angle, i, is determined by the Snell’s law, i.e., sin i = n 1 sin α . F r varies with the incidence-view geometry. Variation of F r as a function of scattering angle is illustrated in Figure 1b. Note that when the incident radiation goes along the leaf normal, i.e., α = i = 0 , Equation (9) becomes meaningless. In this special case, the Fresnel reflectance is defined as:
F r = ( n 1 n + 1 ) 2
The leaf normal distribution function, g L ( Ω L ) , is expressed as [25]:
g L ( Ω L ) = 2 π ( 1 + a   c o s b θ L s i n θ L )
Here, θ L is the polar angle of the leaf normal Ω L . Variables a and b represent planophile (a = 1, b = 2), erectophile (a = −1, b = 2), plagiophile (a = −1, b = 4), extremophile (a = 1, b = 4), and uniform (a = 0) distributions. For the spherical distribution, g L ( Ω L ) = 1 . Note that foliage is assumed to be azimuthally randomly oriented in this paper.
The diffuse area scattering phase function, Γ D , describes the scattered radiation resulting from photon interactions within the leaf interior. The fraction of radiation that enters the leaf interior is ϰ ( α , n , κ ) = 1 Κ F r . It can be described by the following equation:
1 π Γ D ( Ω Ω ) = 1 2 π 2 π + g L ( Ω L ) | Ω · Ω L | ϰ ( α , n , κ ) γ L ( Ω L , Ω Ω ) d Ω L
The leaf specular reflection is thus incorporated into the SRTM via Equations (5)–(12). Note that the improved version of the SRTM suggests that all photons that have not been specularly reflected enter the leaf interior. We will evaluate the model using field-measured data in the following section.

3. Evaluation of the Improved Version of the SRTM

The BRF dataset collected from summer maize at Luancheng Agricultural Ecosystems Experimental Station, Chinese Academy of Sciences, Hebei province, China (37°52′N, 114°39′E) [17] is used to evaluate the improved version of the SRTM. Summer maize represents one of the most typical food and forage crops in Northern China. It was also the subject of several completed sets of field campaigns carried out in July 2002, as part of a “Spectral Knowledge Library of Typical Land Surface Objects in China”. The field dataset includes (1) maize structural parameters, such as average maize height, LAI, average row-width, plant-to-plant interval, etc., and (2) optical properties of leaves, soil background and maize canopy [17]. Characteristics of the maize canopy are presented in Table 1. Note that the dataset has been used for an evaluation of a 3D Scene BRDF model, the radiosity-graphics combined model (RGM) [17]. A detailed description on this station and the dataset is in [17].
Nadir hemispherical reflectance (NHR) and transmittance (NHT) of ~10 pieces of summer maize leaves, ranging from new to mature, were measured on 21 July 2002 under laboratory conditions using a spectrophotometer (SE590, Spectron Engineering, Denver, USA) and an integrating sphere (Analytical Spectral Devices, Inc., Boulder, USA) in the spectral interval from 375 to 1100 nm with a spectral resolution of about 12 nm and a 3 nm sampling interval. Their averages denote the representatives of the maize leaf’s NHR and NHT.
The maize canopy BRF was measured in the near principal plane using the SE590 spectrophotometer at local time 4:00 p.m., on 23 July 2002. The sensor was mounted on the arm of an automatic bracket device (Beijing Normal University, Beijing, China) at a height of 1.88 m, which ensured the maize canopy in the sensor’s field of view ( 15 ° ). The solar zenith angle (SZA) and solar azimuth angles (SAA) were 38.14 ° and 259.18 ° , respectively. For a given view zenith angle (VZA), the BRF measurement could be done in 15 s. As a result, the BRF were collected in the near principal plane in about 6 min. In addition, the ratios of direct to total incident radiation was measured simultaneously. Their values are shown in Table 1. The nadir reflectance of soil background was measured using the SE590 right after the maize canopy BRF measurements.
We selected four bands, blue (436 nm), green (547 nm), red (700 nm) and NIR (804 nm) as representatives for the evaluation. Optical properties of the leaf and soil background are listed in Table 2. The Poisson germ-grain model with an identical ellipsoidal shape is used to parameterize the maize structure. The maize stem density, d, is estimated with the ratio of unity to the product of row-width and plant-to-plant interval, i.e., d = 1 / ( 0.6 × 0.25 ) = 6.67   stem / m 2 . The maximum maize crown radius, r, is assumed to be identical to the plant-to-plant interval, i.e., r = 0.25   m . Given stem density and crown radius, ground cover is calculated as g c = 1 exp ( π d r 2 ) = 0.73 . Leaf angles are assumed to follow the uniform distribution function [17]. The refractive index of the leaf wax, n, and κ are set to 1.5 and 0.3, respectively. The ratio of maize height to radius is adopted in the improved version of the SRTM to simulate the hot spot effect [12].
BRFs are firstly simulated with bi-Lambertian leaves using the classical version of the SRTM and then with non-Lambertian leaves using the improved version of the SRTM. Given that the illumination condition was a mixture of direct solar radiation and diffuse sky radiation (Table 1), each version of the SRTM is executed twice under purely direct solar radiation and purely diffuse sky radiation. The final modeled BRF is a weighted sum of the BRF calculated from both illumination conditions using the ratios of direct to total incident radiation. The modeled BRF is then compared with the field-measured BRF. Their values at blue, green, red, and NIR bands in the principal plane as a function of the view zenith angle is illustrated in Figure 2.
We use the root mean squared error (RMSE) to quantify the proximity between field-measured BRF, B R F m , and modeled BRF, B R F s , i.e.:
R M S E ( B R F m , B R F s ) = 1 M i = 1 M ( B R F s , i B R F m , i ) 2
Here, B R F s , i and B R F m , i represent components of B R F s and B R F m , respectively. M is the number of view directions in the near principal plane. The RMSEs between field-measured and modeled BRF at blue, green, red, and NIR bands are shown in Figure 3.
As shown in Figure 2, the modeled BRF using the improved version of the SRTM fits better with the field-measured data than that using the classical version, indicating that accounting for leaf specular reflection in the SRTM provides more realistic results than when this specular part is neglected. The RMSE is decreased by 16.5%, 24.2%, 35.8%, and 0.15% at blue, green, red, and NIR bands, respectively. One may see the drop of BRF at NIR band around VZA = 25 ° . It is caused by the detector’s or device shelf’s shadow [17]. This effect significantly increases the RMSE at the NIR band, as Figure 3 illustrates. In addition, the heterogeneity of the maize canopy, measurement uncertainty and small environment changes may also attribute to the variations of maize BRF. Figure 2 and Figure 3 thus support the improved version of the SRTM. Our next step is to explore the influence of leaf specular reflection on canopy radiative regime using the improved version of the SRTM.

4. Influence of Leaf Specular Reflection on Canopy Radiative Regime

The modification of the leaf scattering phase function by accounting for leaf specular reflection allows for a more realistic modeling of the canopy radiative regime, which definitely is different from that using the classical version of the SRTM. The goal of this section is to illustrate their differences and show how leaf specular reflection can influence the canopy radiative regime.
In our calculations, the stochastic canopy consists of identical cylindrical vegetation with crown height equal to canopy depth, H , i.e., the stochastic medium is confined to 0 z H . The crown base diameter ( D B ) of the cylindrical-shaped vegetation is assumed to be equal to half of the canopy depth, D B = 0.5 H . Leaves are randomly distributed (i.e., turbid medium) within the crowns. In this case, the probability of finding a vegetated point at depth z is a constant and coincides with ground cover, gc. To specify the relationship between canopy LAI and ground cover, plant LAI, L 0 , which is defined as the amount of leaf area within canopy crowns (i.e., L 0 = u L H ), is introduced. The following equation is thus derived: LAI = L 0 · g c . A detailed description on the above canopy structural parameters is documented in [11].
Leaf angles are assumed to follow the spherical distribution function. The bi-Lambertian model is adopted to describe the leaf diffuse scattering part. Leaf diffuse hemispherical reflectance and transmittance are assumed to be identical and set to 0.07 at red and 0.38 at NIR bands. The leaf specular scattering part is parameterized using Fresnel reflectance with a constant refractive index, n = 1.5 , at both bands. The reduction factor is set to κ = 0.3 . The vegetation canopy is illuminated by a parallel beam of unity intensity with a SZA and SAA of 30 ° and 0 ° , respectively.
The calculation includes two steps. Firstly, the improved version of SRTM is solved for mean intensity I ¯ ( z , Ω ) with leaf specular reflection accounted for in the leaf scattering phase function. Secondly, the first step is repeated using the classical version of SRTM, i.e., leaf specular reflection is neglected. In order to distinguish the mean intensity calculated from the two steps, we use subscripts “s” and “ns”, i.e., I ¯ s ( z , Ω ) and I ¯ n s ( z , Ω ) , to represent the mean intensities calculated from Step 1 and Step 2, respectively. The differences between I ¯ s ( z , Ω ) and I ¯ n s ( z , Ω ) will be utilized as measures to quantify the influence of leaf specular reflection on the canopy radiative regime.

4.1. Vertical Profiles of Mean Radiation Flux Density

In this subsection, we investigate the influence of leaf specular reflection on the canopy radiative regime by comparing the vertical profiles of mean radiation flux densities. The mean downward and upward flux densities from calculations using the improved version of the SRTM, F I ¯ , s ( z ) , and the classical version of the SRTM, F I ¯ , n s ( z ) , are defined as follows [11]:
F J ( z ) = 1 i 0 | μ 0 | 2 π J ( z , Ω ) | μ | d Ω
Here, J ( z , Ω ) denotes either I ¯ s ( z , Ω ) or I ¯ n s ( z , Ω ) . i 0 is the intensity of the incident radiation, which is unity in our calculation. μ 0 and μ are cosine values of the SZA and the polar angle of direction Ω . 2 π and 2 π + represent integrations over downward and upward hemispherical spaces, respectively.
Figure 4 illustrates the vertical profiles of mean downward and upward radiation flux densities ( F I ¯ , s ( z ) and F I ¯ , n s ( z ) , respectively) at red and NIR bands for four types of ground cover conditions. In the former case when leaf specular reflection is accounted for, leaf albedo, ω L , is specified as the sum of leaf specular reflection and bi-Lambertian scattering, i.e., ω L = ( 1 s L ) ω L D + s L , where s L is the leaf specular part and ω L D is the leaf diffuse scattering part, given that photons interact with leaf internal constituents. However, in the latter case when leaf specular reflection is neglected ( s L = 0 ), leaf albedo, ω L = ω L D . The value of ω L is always greater than ω L D because of the difference ω L ω L D = s L ( 1 ω L D ) 0 . Thus, the leaf scattering is stronger in the case when leaf specular reflection is accounted for. In this example, the flux density can either be scattered downward or upward, indicating that accounting for leaf specular radiation increases both mean downward and upward radiation flux densities, i.e., F I ¯ , s ( z ) F I ¯ , n s ( z ) , as Figure 4 illustrates. The difference, F I ¯ , s ( z ) F I ¯ , n s ( z ) , is thus 0 , as Figure 5 shows.
For the mean downward radiation flux density, a smaller ground cover results in less incident radiation being intercepted by the vegetation canopy. Since the intercepted radiation triggers future leaf specular reflection and consequently generates the difference ( F I ¯ , s ( z ) F I ¯ , n s ( z ) ), the smaller ground cover is, the smaller difference will be, as Figure 5a,c illustrates. The intercepted radiation is thus primarily responsible for propagation of leaf surface reflection in the downward direction. It follows that the difference also increases with depth, z. It reaches its maximum at the bottom of the vegetation canopy where the normalized depth, z/H = 1. However, for a given canopy LAI (LAI = 1.5 in this example), the difference can be saturated. In such a case, the leaf area volume density within the crowns is dense and the intercepted radiation saturates with the increase of depth z, as the green rectangles (gc = 0.25 and u L = 6 m 2 / m 3 ) in Figure 5a,c illustrates.
As for the difference in mean upward radiation flux density, it has the opposite tendency, i.e., the difference decreases with normalized depth. The deeper the depth is, the lower probability of the leaf specular reflection being generated. This consequently involves a decrease in the difference, F I ¯ , s ( z ) F I ¯ , n s ( z ) , as one can see from Figure 5b,d. The difference also varies with ground cover, especially at strong absorbing band, i.e., the red band in this example. At the red band, the amount of leaf specular reflection generated from single scattering determines the difference. For a given canopy LAI, such a single scattering is an increasing function with respect to ground cover. However, at the NIR band, leaves absorb little radiation. Photons can undergo even more than 10 interactions with leaves before they are either absorbed or exit the canopy [26]. Multiple scattering is so strong that the difference under different ground cover conditions is minimized, as Figure 5d illustrates.
Overall, values of the differences are below 1% at red and NIR bands. However, the difference at red can account for up to 1.8% (=0.0079/0.4473) and 13% (=0.0040/0.0287) of the mean radiation flux density for the downward and upward directions. As for the NIR, the influences are smaller with up to 1.1% (=0.0062/0.5769) and 0.6% (=0.0012/0.1960) for the downward and upward directions.

4.2. Canopy Absorptance

Canopy absorptance conveys information on the amount of incident radiation that is potentially utilized in the photosynthetic process. Accurate monitoring of canopy absorptance is critical for many areas of research in hydrology, ecology, and climate change. The distribution of incident radiation within the vegetation canopy is mainly determined by canopy structure, which partitions canopy absorption, transmission and reflection. Additionally, leaf specular reflection decreases the amount of radiation that can enter leaf interiors, and consequently impacts canopy absorptance. In this subsection, we show examples of the influence of leaf specular reflection on canopy absorptance.
Figure 6a illustrates mean canopy absorptance for a vegetation canopy bounded from below by a non-reflecting background, which is calculated as 1 F J ( 0 ) F J ( 1 ) . Here, F J ( 0 ) represents canopy reflectance and F J ( 1 ) denotes canopy transmittance, as Equation (14) defines. Let A s and A n s be canopy absorptances calculated with and without accounting for leaf specular reflection, respectively. Neglecting leaf specular reflection underestimates leaf albedo and, therefore, overestimates canopy absorptance when the underestimated leaf albedo is applied to the stochastic radiative transfer simulations, i.e., A s A n s , as can be seen in Figure 6a.
The difference, A n s A s , is given in Figure 6b. As for the 1D canopy medium, the difference firstly increases and then decreases with canopy LAI. This is as expected because when the vegetation canopy is sparse, canopy interceptance will increase with the increase of canopy LAI, which indicates more incident radiation can be specular reflected and escapes the vegetation canopy from upper or lower boundaries. This in turn causes greater difference between A n s and A s . However, when the vegetation canopy is dense enough, canopy interceptance does not increase rapidly with the increase of canopy LAI. Despite there still being more incident radiation which can be intercepted and specularly reflected by the vegetation canopy, a lesser amount of specular reflected radiation can exit the vegetation canopy. These two opposite features cannot compensate for each other and causes a decrease of difference, A n s A s , with the increase of canopy LAI. As for the stochastic canopy medium, we get similar tendencies, as the hollow symbols in Figure 6b illustrate. However, the variation rate is lower than that of a 1D canopy because considering the within and between crown radiation regimes provides a more accurate saturation domain, i.e., the absorption is saturated at a slower rate for the stochastic medium than that of a 1D medium. A detailed relationship between the saturation domain and the stochastic or 1D canopy medium is discussed in [11].

4.3. BRF and NDVI

The reflective property of the vegetated surface is widely interpreted in the remote sensing community. Such a property is usually described by BRF. Figure 7 shows canopy BRF at red and NIR bands in the nadir view direction. The canopy is bounded from below by a reflecting background with a surface albedo equals 0.18 at both bands. As one can see, canopy BRFs vary with ground cover and have opposite tendencies for red and NIR bands. For the red band, canopy BRF decreases with ground cover. At this band, little intercepted solar radiation is scattered by leaves. The contribution of background, therefore, dominates canopy BRF. Such contribution decreases with ground cover and consequently canopy BRF decreases. However, for the NIR band, canopy BRF increases with ground cover. Despite the decreasing contribution of the background, the strong multiple scattering between leaves increases canopy BRF.
As discussed earlier, accounting for leaf specular reflection increases leaf albedo, which indicates that more radiation can be scattered and escape the canopy from the upper boundary. Consequently, the canopy BRF increases, as Figure 7 shows. The difference can account for up to 17% (=(0.035 − 0.030)/0.030) and 2% (=(0.317 − 0.311)/0.311) of the canopy nadir BRF at red and NIR bands when leaf specular reflection is neglected.
Canopy BRF is often transformed to vegetation indices for the remote sensing of canopy nitrogen content, LAI, and fraction of absorbed photosynthetically active radiation (FPAR) [27,28,29,30], etc. Numerous vegetation indices are thus proposed. Among them, NDVI is quite easy to achieve and usually provides satisfactory results regarding such researches. Thus, it is one of the most popular vegetation indices. Let N D V I n s and N D V I s be the NDVI calculated with and without accounting for leaf specular reflection, respectively. Here, we focus on the NDVI at the nadir direction and investigate the influence of leaf specular reflection on nadir NDVI. Figure 8 illustrates N D V I n s , N D V I s and their difference, N D V I n s N D V I s , as a function of canopy LAI for four types of plant LAI. One can see that the results are quite similar to those presented in Figure 6. The difference in the leaf scattering phase function is primarily responsible for this effect. When leaf specular reflection is accounted for, more radiation is scattered and escapes the canopy from upper or lower boundaries, which increases canopy BRF at red and NIR bands (Figure 7) and, consequently, decreases canopy NDVI, as Figure 8a illustrates. It follows from Figure 6b that the difference, N D V I n s N D V I s , firstly increases and then decreases with canopy LAI at difference variation rates for the 1D and stochastic canopy medium, as Figure 8b shows.
To summarize, the classical version of the SRTM treats leaves as ideal bi-Lambertian scatters and does not account for leaf specular reflection and does not discriminate between radiation scattered from leaf surface and leaf interior. For a vegetation canopy with spherically inclined leaves, this results in, on one hand, underestimation of mean downward/upward radiation flux densities (Figure 4 and Figure 5) and canopy nadir BRF (Figure 7). On the other hand, it causes an overestimation of canopy absorptance (Figure 6) and NDVI (Figure 8). Accounting for leaf specular reflection in RT modeling is thus necessary for the remote sensing of vegetation because it provides a more accurate estimate of the canopy radiative regime.

5. Discussion

Solar radiation intercepted by vegetation canopies can either be scattered or absorbed. A passive optical instrument measures the amount of radiation below or/and above the canopy, which is further transferred to many variables that are linked to the vegetation functioning, such as the canopy BRF, vegetation indices, LAI/FPAR, nitrogen content, etc. A fraction of the instrument registered radiation is caused by leaf specular reflection (Figure 1). This fraction of radiation does not enter leaf interiors and exhibits a strong dependence on sun-sensor geometry. We use the Fresnel reflectance (Equation (9)) to parameterize the leaf specular reflection and propose an improved version of the SRTM. In this model, the leaf scattering phase function is modified accordingly (Equations (5), (7) and (12)) by assuming that all photons that have not been specularly reflected enter the leaf interior.
Canopy BRF estimated from the improved model is evaluated through comparisons with field-measured maize BRF. The maize in this paper is planted in rows. The structure of the maize canopy cannot be simplified as a classical 1D medium if one wants to model the maize BRF at a high accuracy. The structure, however, is critical in the BRF modeling because it impacts the distribution of specularly reflected and diffusely scattered radiation. Indeed, without increasing much computation cost, the row structure can be accurately parameterized using a parameter called pair correlation function in the SRTM [11]. Our results in Figure 2 are quite similar to those reported in [17], which were achieved using a 3D RGM. The improved SRTM is, thus, a powerful tool for the remote sensing of vegetation because of its efficiency and accuracy. A detailed description on the pair correlation function and RGM are documented in [11,17], respectively. The improved model can, therefore, be used for investigating the influences of leaf specular reflection on the canopy radiative regime.
Accounting for the leaf specular reflection in the leaf scattering phase function increases leaf albedo, i.e., for a given amount of radiation incident upon a leaf, more radiation can be scattered and less radiation can be consequently absorbed (Figure 4, Figure 5 and Figure 6). The distribution of canopy scattered radiation is determined by the canopy structure (e.g., LAI, leaf angles and ground cover). Despite overall canopy scattered radiation increasing, the reflectance or transmittance may increase or decrease for some directions. In our example, the nadir BRF increases (Figure 7). However, in other cases, it may decrease, as shown in [17]. The variation of NDVI (Figure 8) may also be different as it is a mathematical product of canopy BRFs. Therefore, one should pay special attention to leaf specular reflection when it comes to remote sensing variables that are direction dependent, especially for the strongly absorbing band where the canopy specularly reflected radiation may be comparable or even greater than the canopy diffusely scattered radiation [17,31].
The influence is greater for downward/upward radiation flux densities and canopy nadir BRF than that for canopy absorptance and NDVI. This is mainly attributed to: (1) canopy absorptance being usually much greater than the contribution from leaf specular reflection and (2) NDVI tending to reduce (but not eliminate) the influence of leaf specular reflection. Our results suggest that the difference caused by leaf specular reflection can be up to 13% and 17% of the vertical profiles of radiation flux density and canopy nadir BRF, respectively. This may indicate that leaf specular reflection should be accounted for in the case when the error limit for modeling and measurement of vertical profiles of radiation flux density and canopy BRF is smaller than 13%.
Note that canopy BRF at red and NIR bands and NDVI have been used in operational algorithm for generating global LAI/FPAR products from the MODIS, VIIRIS, and EPIC datasets [13,27,28,29]. Neglecting leaf specular reflection in the photon-canopy interactions causes underestimations of canopy BRF at red and NIR bands and overestimation of canopy NDVI which, in turn, will result in overestimation of canopy LAI/FPAR [28,29,32,33,34].
Besides the above influences, leaf specular reflection has recently been reported to decreases the accuracy of LAI estimates from LAI-2000 field measurements [35] and remote sensing of leaf biochemical constituents [31,36,37]. Determining the influence of leaf specular reflection on the canopy radiative regime is, therefore, critical because the inputs in these studies are generally related to, or determined, by the canopy radiative regime. Moreover, field measurements of vertical profiles of mean flux density are also necessary for a better understanding of the influence of leaf specular reflection, as well as for further validations and applications of the improved SRTM. It should be noted that many empirical and theoretical analyses show that leaf specular reflection is partly polarized, whereas the leaf diffuse scattering is not [38,39,40,41]. Polarization measurement can, thus, be useful for a full understanding of leaf specular reflection and its influences in the remote sensing community [42,43].

6. Conclusions

Leaf specular reflection is a particular feature that should be accounted for in canopy RT modeling. For this reason, an improved version of the SRTM is proposed in this paper. Compared to the classical version of the SRTM, the leaf scattering phase function is modified to include leaf specular reflection by assuming that all photons that have not been specularly reflected enter the leaf interior. The improved SRTM satisfactorily estimates maize BRF over a wide range of zenith view angles. From a practical point of view, it can efficiently and accurately provide estimates of the canopy radiative regime. It is thus desirable if one considers that leaf specular reflection lowers the sensitivity of sensor-registered signals to canopy structure, leaf biochemistry, etc.
Leaf specular reflection influences the canopy radiative regime in vertical profiles of mean radiation flux density, canopy absorptance, BRF, NDVI, etc. Analyses of the improved version of SRTM for a vegetation canopy with spherically inclined leaves shows that accounting for leaf specular reflection increases mean upward/downward radiation flux density and nadir BRF, and it decreases canopy absorptance and nadir NDVI. The influence is greater for vertical profiles of mean radiation flux density and nadir BRF than that for canopy absorptance and NDVI. Ignoring this part of specularly reflected radiation, therefore, may introduce unpredictable errors when interpreting remotely sensed data. Moreover, a small portion of incident radiation is scattered by the dust and hair on leaf surfaces. The dust and hair can somehow influence the canopy radiative regime and their influences need to be addressed in future studies.
It also should be noted that polarization measurement may be adopted for studies of leaf specular reflection. In the future, we would also expect to extend the improved version of the SRTM by accounting for canopy polarization, which can be helpful for processing remote-sensing data collected by polarization sensors, such as the POLDER, RSP, 3MI, and GF-5, etc.

Author Contributions

B.Y. and Y.K. conceived and designed the model; D.X. collected the data; B.Y. and D.X. analyzed the data; H.Z., J.Z. and Y.W. contributed to the analysis section; and B.Y. wrote the paper.

Funding

This work was supported by the Fundamental Research Funds for the Central Universities under Grant No. 531107051063, the Rainbow Plan of Yusense Information Technology and Equipment (Qingdao) Ltd. under Grant No. YCCH20080101, the National Key Research and Development Program of China under Grant No. 2017YFB0503004-4 and the National Natural Science Foundation of China under grant No. 41701391.

Acknowledgments

The author would like to thank Dong Huang from Brookhaven National Laboratory for the fruitful discussion on stochastic radiative transfer model.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, Q.; Liang, S.; Xiao, Z.; Fang, H. Retrieval of leaf area index using temporal, spectral, and angular information from multiple satellite data. Remote Sens. Environ. 2014, 145, 25–37. [Google Scholar] [CrossRef]
  2. Knyazikhin, Y.; Martonchik, J.; Myneni, R.; Diner, D.; Running, S. Synergistic algorithm for estimating vegetation canopy leaf area index and fraction of absorbed photosynthetically active radiation from modis and misr data. J. Geophys. Res. Atmos. 1998, 103, 32257–32275. [Google Scholar] [CrossRef]
  3. Diner, D.J.; Asner, G.P.; Davies, R.; Knyazikhin, Y.; Muller, J.P.; Nolin, A.W.; Pinty, B.; Schaaf, C.B.; Stroeve, J. New directions in earth observing: Scientific applications of multiangle remote sensing. Bull. Am. Meteorol. Soc. 1999, 80, 2209–2228. [Google Scholar] [CrossRef]
  4. Zeng, Y.; Li, J.; Liu, Q.; Huete, A.R.; Yin, G.; Xu, B.; Fan, W.; Zhao, J.; Yan, K.; Mu, X. A radiative transfer model for heterogeneous agro-forestry scenarios. IEEE Trans. Geosci. Remote Sens. 2016, 54, 4613–4628. [Google Scholar] [CrossRef]
  5. Kuusk, A. A two-layer canopy reflectance model. J. Quant. Spectrosc. Radiat. Transf. 2001, 71, 1–9. [Google Scholar] [CrossRef]
  6. Pinty, B.; Lavergne, T.; Dickinson, R.E.; Widlowski, J.L.; Gobron, N.; Verstraete, M.M. Simplifying the interaction of land surfaces with radiation for relating remote sensing products to climate models. J. Geophys. Res. Atmos. 2006, 111, 255–266. [Google Scholar] [CrossRef]
  7. Verhoef, W. Light scattering by leaf layers with application to canopy reflectance modeling: The sail model. Remote Sens. Environ. 1984, 16, 125–141. [Google Scholar] [CrossRef]
  8. Gastellu-Etchegorry, J.P.; Martin, E.; Gascon, F. Dart: A 3d model for simulating satellite images and studying surface radiation budget. Int. J. Remote Sens. 2004, 25, 73–96. [Google Scholar] [CrossRef]
  9. Kuusk, A.; Nilson, T. A directional multispectral forest reflectance model. Remote Sens. Environ. 2000, 72, 244–252. [Google Scholar] [CrossRef]
  10. Myneni, R.B.; Asrar, G.; Hall, F.G. A three-dimensional radiative transfer method for optical remote sensing of vegetated land surfaces. Remote Sens. Environ. 1992, 41, 105–121. [Google Scholar] [CrossRef]
  11. Huang, D.; Knyazikhin, Y.; Wang, W.; Deering, D.W.; Stenberg, P.; Shabanov, N.; Tan, B.; Myneni, R.B. Stochastic transport theory for investigating the three-dimensional canopy structure from space measurements. Remote Sens. Environ. 2008, 112, 35–50. [Google Scholar] [CrossRef]
  12. Shabanov, N.V.; Knyazikhin, Y.; Baret, F.; Myneni, R.B. Stochastic modeling of radiation regime in discontinuous vegetation canopies. Remote Sens. Environ. 2000, 74, 125–144. [Google Scholar] [CrossRef]
  13. Yang, B.; Knyazikhin, Y.; Mõttus, M.; Rautiainen, M.; Stenberg, P.; Yan, L.; Chen, C.; Yan, K.; Choi, S.; Park, T.; et al. Estimation of leaf area index and its sunlit portion from dscovr epic data: Theoretical basis. Remote Sens. Environ. 2017, 198, 69–84. [Google Scholar] [CrossRef] [PubMed]
  14. Verhoef, W.; Jia, L.; Xiao, Q.; Su, Z. Unified optical-thermal four-stream radiative transfer theory for homogeneous vegetation canopies. IEEE Trans. Geosci. Remote Sens. 2007, 45, 1808–1822. [Google Scholar] [CrossRef]
  15. Vanderbilt, V.C.; Grant, L.; Ustin, S.L. Polarization of Light by Vegetation; Springer: Berlin/Heidelberg, Germany, 1991. [Google Scholar]
  16. Grant, L. Diffuse and specular characteristics of leaf reflectance. Remote Sens. Environ. 1987, 22, 309–322. [Google Scholar] [CrossRef]
  17. Xie, D.; Qin, W.; Wang, P.; Shuai, Y.; Zhou, Y.; Zhu, Q. Influences of leaf-specular reflection on canopy BRF characteristics: A case study of real maize canopies with a 3-d scene brdf model. IEEE Trans. Geosci. Remote Sens. 2017, 55, 619–631. [Google Scholar] [CrossRef]
  18. Goel, N.S. Models of vegetation canopy reflectance and their use in estimation of biophysical parameters from reflectance data. Remote Sens. Rev. 1988, 4, 1–212. [Google Scholar] [CrossRef]
  19. Ross, J.; Marshak, A. The influence of leaf orientation and the specular component of leaf reflectance on the canopy bidirectional reflectance. Remote Sens. Environ. 1989, 27, 251–260. [Google Scholar] [CrossRef]
  20. Gastellu-Etchegorry, J.P.; Demarez, V.; Pinel, V.; Zagolski, F. Modeling radiative transfer in heterogeneous 3-d vegetation canopies. Remote Sens. Environ. 1996, 58, 131–156. [Google Scholar] [CrossRef]
  21. Andrieu, B.; Baret, F.; Jacquemoud, S.; Malthus, T.; Steven, M. Evaluation of an improved version of sail model for simulating bidirectional reflectance of sugar beet canopies. Remote Sens. Environ. 1997, 60, 247–257. [Google Scholar] [CrossRef]
  22. Schaepman-Strub, G.; Schaepman, M.; Painter, T.; Dangel, S.; Martonchik, J. Reflectance quantities in optical remote sensing—definitions and case studies. Remote Sens. Environ. 2006, 103, 27–42. [Google Scholar] [CrossRef]
  23. Knyazikhin, Y.M.A. Fundamental equations of radiative transfer in leaf canopies, and iterative methods for their solution. In Photon-Vegetation Interactions; Myneni, R.B., Ross, J., Eds.; Springer: Berlin/Heidelberg, Germany, 1991. [Google Scholar]
  24. Nilson, T.; Kuusk, A. Approximate analytic relationships for the reflectance of agricultural vegetation canopies. Sov. J. Remote Sens. 1984, 4, 814–826. [Google Scholar]
  25. Bunnik, N.J.J. The Multispectral Reflectance of Shortwave Radiation by Agricultural Crops in Relation with Their Morphological and Optical Properties. Doctoral Dissertation, Wageningen University, Wageningen, The Netherlands, 1978. [Google Scholar]
  26. Huang, D.; Knyazikhin, Y.; Dickinson, R.E.; Rautiainen, M.; Stenberg, P.; Disney, M.; Lewis, P.; Cescatti, A.; Tian, Y.; Verhoef, W. Canopy spectral invariants for remote sensing and model applications. Remote Sens. Environ. 2007, 106, 106–122. [Google Scholar] [CrossRef] [Green Version]
  27. Chen, C.; Knyazikhin, Y.; Park, T.; Yan, K.; Lyapustin, A.; Wang, Y.; Yang, B.; Myneni, R. Prototyping of lai and fpar retrievals from modis multi-angle implementation of atmospheric correction (maiac) data. Remote Sens. 2017, 9, 370. [Google Scholar] [CrossRef]
  28. Yan, K.; Park, T.; Chen, C.; Xu, B.; Song, W.; Yang, B.; Zeng, Y.; Liu, Z.; Yan, G.; Knyazikhin, Y. Generating global products of lai and fpar from snpp-viirs data: Theoretical background and implementation. IEEE Trans. Geosci. Remote Sens. 2018, 2119–2137. [Google Scholar] [CrossRef]
  29. Yan, K.; Park, T.; Yan, G.; Chen, C.; Yang, B.; Liu, Z.; Nemani, R.; Knyazikhin, Y.; Myneni, R. Evaluation of modis lai/fpar product collection 6. Part 1: Consistency and improvements. Remote Sens. 2016, 8, 359. [Google Scholar] [CrossRef]
  30. Wang, Z.; Wang, T.; Darvishzadeh, R.; Skidmore, A.; Jones, S.; Suarez, L.; Woodgate, W.; Heiden, U.; Heurich, M.; Hearne, J. Vegetation indices for mapping canopy foliar nitrogen in a mixed temperate forest. Remote Sens. 2016, 8, 491. [Google Scholar] [CrossRef]
  31. Yang, B.; Knyazikhin, Y.; Lin, Y.; Yan, K.; Chen, C.; Park, T.; Choi, S.; Mõttus, M.; Rautiainen, M.; Myneni, R.; et al. Analyses of impact of needle surface properties on estimation of needle absorption spectrum: Case study with coniferous needle and shoot samples. Remote Sens. 2016, 8, 563. [Google Scholar] [CrossRef] [PubMed]
  32. Xu, B.; Park, T.; Yan, K.; Chen, C.; Zeng, Y.; Song, W.; Yin, G.; Li, J.; Liu, Q.; Knyazikhin, Y. Analysis of global lai/fpar products from viirs and modis sensors for spatio-temporal consistency and uncertainty from 2012–2016. Forests 2018, 9, 73. [Google Scholar] [CrossRef]
  33. Yan, K.; Park, T.; Yan, G.; Liu, Z.; Yang, B.; Chen, C.; Nemani, R.; Knyazikhin, Y.; Myneni, R. Evaluation of modis lai/fpar product collection 6. Part 2: Validation and intercomparison. Remote Sens. 2016, 8, 460. [Google Scholar] [CrossRef]
  34. Zeng, Y.; Li, J.; Liu, Q.; Hu, R.; Mu, X.; Fan, W.; Xu, B.; Yin, G.; Wu, S. Extracting leaf area index by sunlit foliage component from downward-looking digital photography under clear-sky conditions. Remote Sens. 2015, 7, 13410–13435. [Google Scholar] [CrossRef]
  35. Kuusk, A. Specular reflection in the signal of lai-2000 plant canopy analyzer. Agric. For. Meteorol. 2016, 221, 242–247. [Google Scholar] [CrossRef]
  36. Knyazikhin, Y.; Schull, M.A.; Stenberg, P.; Mottus, M.; Rautiainen, M.; Yang, Y.; Marshak, A.; Carmona, P.L.; Kaufmann, R.K.; Lewis, P.; et al. Hyperspectral remote sensing of foliar nitrogen content. Proc. Natl. Acad. Sci. USA 2013, 110, E185–E192. [Google Scholar] [CrossRef] [PubMed]
  37. Latorre-Carmona, P.; Knyazikhin, Y.; Alonso, L.; Moreno, J.F.; Pla, F.; Yan, Y. On hyperspectral remote sensing of leaf biophysical constituents: Decoupling vegetation structure and leaf optics using chris–proba data over crops in barrax. IEEE Geosci. Remote Sens. Lett. 2014, 11, 1579–1583. [Google Scholar] [CrossRef]
  38. Diner, D.; Xu, F.; Martonchik, J.; Rheingans, B.; Geier, S.; Jovanovic, V.; Davis, A.; Chipman, R.; McClain, S. Exploration of a polarized surface bidirectional reflectance model using the ground-based multiangle spectropolarimetric imager. Atmosphere 2012, 3, 591–619. [Google Scholar] [CrossRef]
  39. Grant, L.; Daughtry, C.; Vanderbilt, V. Polarized and specular reflectance variation with leaf surface features. Physiol. Plant. 1993, 88, 1–9. [Google Scholar] [CrossRef]
  40. Sun, Z.; Wu, D.; Lv, Y.; Zhao, Y. Bidirectional polarized reflectance factors of vegetation covers: Influence on the BRF models results. IEEE Trans. Geosci. Remote Sens. 2017, 55, 5687–5701. [Google Scholar] [CrossRef]
  41. Suomalainen, J.; Hakala, T.; Puttonen, E.; Peltoniemi, J. Polarised bidirectional reflectance factor measurements from vegetated land surfaces. J. Quant. Spectrosc. Radiat. Transf. 2009, 110, 1044–1056. [Google Scholar] [CrossRef]
  42. Breon, F.M.; Maignan, F. A brdf–bpdf database for the analysis of earth target reflectances. Earth Syst. Sci. Data 2017, 9, 31–45. [Google Scholar] [CrossRef]
  43. Yang, B.; Zhao, H.; Chen, W. Semi-empirical models for polarized reflectance of land surfaces: Intercomparison using space-borne polder measurements. J. Quant. Spectrosc. Radiat. Transf. 2017, 202, 13–20. [Google Scholar] [CrossRef]
Figure 1. (a) Distribution of leaf specularly reflected, diffusely reflected and transmitted intensity. Radiation scattered by a leaf consists of two parts, diffuse and specular. The first part exists in the 4 π spherical space and can be described using the bi-Lambertian model. The specular part appears at the direction of specular reflection and can be expressed in terms of Fresnel reflectance. Ω , Ω L , Ω , and Ω * are the directions of incident radiation, leaf normal, scattered radiation, and specular reflection, respectively. α represents the incident angle, i.e., the polar angle between Ω L and Ω . (b) Variation of Fresnel reflectance as a function of scattering angle.
Figure 1. (a) Distribution of leaf specularly reflected, diffusely reflected and transmitted intensity. Radiation scattered by a leaf consists of two parts, diffuse and specular. The first part exists in the 4 π spherical space and can be described using the bi-Lambertian model. The specular part appears at the direction of specular reflection and can be expressed in terms of Fresnel reflectance. Ω , Ω L , Ω , and Ω * are the directions of incident radiation, leaf normal, scattered radiation, and specular reflection, respectively. α represents the incident angle, i.e., the polar angle between Ω L and Ω . (b) Variation of Fresnel reflectance as a function of scattering angle.
Remotesensing 10 01632 g001
Figure 2. Comparison between field-measured and modeled maize BRF at blue (a), green (b), red (c) and NIR (d) bands. Legend “Simulated, Lam” represents the modeled BRF with bi-Lambertian leaves, i.e., leaf specular reflection is neglected. Legend “Simulated, Spec” denotes that accounting for leaf specular reflection. Legend “measured” means field-measured maize BRF.
Figure 2. Comparison between field-measured and modeled maize BRF at blue (a), green (b), red (c) and NIR (d) bands. Legend “Simulated, Lam” represents the modeled BRF with bi-Lambertian leaves, i.e., leaf specular reflection is neglected. Legend “Simulated, Spec” denotes that accounting for leaf specular reflection. Legend “measured” means field-measured maize BRF.
Remotesensing 10 01632 g002aRemotesensing 10 01632 g002b
Figure 3. The RMSEs between field-measured and modeled BRF at blue, green, red and NIR bands. “Lam” represents the RMSE calculated with bi-Lambertian leaves, i.e., leaf specular reflection is neglected, whereas “Spec” denotes that accounting for leaf specular reflection.
Figure 3. The RMSEs between field-measured and modeled BRF at blue, green, red and NIR bands. “Lam” represents the RMSE calculated with bi-Lambertian leaves, i.e., leaf specular reflection is neglected, whereas “Spec” denotes that accounting for leaf specular reflection.
Remotesensing 10 01632 g003
Figure 4. Vertical profiles of mean downward (a) and upward (b) radiation flux densities over the entire horizontal plane at red band for four types of ground cover conditions. The 1D vegetation canopy (“1D”) represents the case when ground cover is unity, i.e., gc = 1.0. Solid symbols and lines correspond to the mean radiation flux densities when leaf specular reflection is accounted for, whereas hollow symbols and dashed lines denote that when leaves are bi-Lambertian scatters. The canopy is bounded from below by a non-reflecting background. Canopy LAI is constant and equals 1.5. Plant LAI varies with ground cover as 1.5/gc. The vertical axis shows the normalized depth, which is defined as z / H were H is the canopy depth.
Figure 4. Vertical profiles of mean downward (a) and upward (b) radiation flux densities over the entire horizontal plane at red band for four types of ground cover conditions. The 1D vegetation canopy (“1D”) represents the case when ground cover is unity, i.e., gc = 1.0. Solid symbols and lines correspond to the mean radiation flux densities when leaf specular reflection is accounted for, whereas hollow symbols and dashed lines denote that when leaves are bi-Lambertian scatters. The canopy is bounded from below by a non-reflecting background. Canopy LAI is constant and equals 1.5. Plant LAI varies with ground cover as 1.5/gc. The vertical axis shows the normalized depth, which is defined as z / H were H is the canopy depth.
Remotesensing 10 01632 g004
Figure 5. Difference between mean radiation flux densities when the leaf surface reflection is accounted for and when it is neglected (i.e., leaves are bi-Lambertian scatters) over the entire horizontal plane at red (a,b) and near infrared (c,d) bands for four types of ground cover conditions. (a,c) denote the vertical profiles of downward differences, whereas (b,d) represent that of upward differences. Solid lines correspond to the 1D vegetation canopy (gc = 1.0) and symbols represent gc = 0.25, 0.5, and 0.75, respectively. Other parameters are the same as in Figure 4.
Figure 5. Difference between mean radiation flux densities when the leaf surface reflection is accounted for and when it is neglected (i.e., leaves are bi-Lambertian scatters) over the entire horizontal plane at red (a,b) and near infrared (c,d) bands for four types of ground cover conditions. (a,c) denote the vertical profiles of downward differences, whereas (b,d) represent that of upward differences. Solid lines correspond to the 1D vegetation canopy (gc = 1.0) and symbols represent gc = 0.25, 0.5, and 0.75, respectively. Other parameters are the same as in Figure 4.
Remotesensing 10 01632 g005
Figure 6. (a) Mean canopy absorptance at red and NIR bands as a function of canopy LAI. Solid symbols and lines correspond to the mean canopy absorptance when leaf specular reflection is accounted for, whereas hollow symbols and dashed lines denote that when leaves are bi-Lambertian scatters. (b) Difference between mean canopy absorptance when leaf specular reflection is neglected (i.e., leaves are bi-Lambertian scatters) and when it is accounted for at red and NIR bands. The canopy is bounded from below by a non-reflecting background. Plant LAI is constant and equals 7.0. Ground cover varies with canopy LAI as LAI/7.0. Other parameters are the same as in Figure 4.
Figure 6. (a) Mean canopy absorptance at red and NIR bands as a function of canopy LAI. Solid symbols and lines correspond to the mean canopy absorptance when leaf specular reflection is accounted for, whereas hollow symbols and dashed lines denote that when leaves are bi-Lambertian scatters. (b) Difference between mean canopy absorptance when leaf specular reflection is neglected (i.e., leaves are bi-Lambertian scatters) and when it is accounted for at red and NIR bands. The canopy is bounded from below by a non-reflecting background. Plant LAI is constant and equals 7.0. Ground cover varies with canopy LAI as LAI/7.0. Other parameters are the same as in Figure 4.
Remotesensing 10 01632 g006
Figure 7. Bidirectional reflectance factor (BRF) at red and NIR bands in the nadir view direction as a function of ground cover. Solid symbols correspond to the canopy BRF when leaf specular reflection is accounted for, whereas hollow symbols denote that when leaves are bi-Lambertian scatters. Canopy LAI is constant and equals 7.0. Plant LAI varies with ground cover as 7/gc. Albedo of the canopy background is 0.18 at both red and NIR bands. SZA and SAA are 30 ° and 0 ° , respectively.
Figure 7. Bidirectional reflectance factor (BRF) at red and NIR bands in the nadir view direction as a function of ground cover. Solid symbols correspond to the canopy BRF when leaf specular reflection is accounted for, whereas hollow symbols denote that when leaves are bi-Lambertian scatters. Canopy LAI is constant and equals 7.0. Plant LAI varies with ground cover as 7/gc. Albedo of the canopy background is 0.18 at both red and NIR bands. SZA and SAA are 30 ° and 0 ° , respectively.
Remotesensing 10 01632 g007
Figure 8. (a) Normalized difference vegetation index (NDVI) as a function of canopy LAI for four types of plant LAI. The solid symbols and line correspond to NDVI when leaf specular reflection is accounted for, whereas the hollow symbols and dashed line denote that when leaves are bi-Lambertian scatters. (b) Difference of NDVI when leaf specular reflection is neglected and when it is accounted for. The canopy is bounded from below by a reflective background with a surface albedo equals 0.18 at red and NIR bands. Ground cover varies with canopy LAI as LAI/7.0. Other parameters are as in Figure 6.
Figure 8. (a) Normalized difference vegetation index (NDVI) as a function of canopy LAI for four types of plant LAI. The solid symbols and line correspond to NDVI when leaf specular reflection is accounted for, whereas the hollow symbols and dashed line denote that when leaves are bi-Lambertian scatters. (b) Difference of NDVI when leaf specular reflection is neglected and when it is accounted for. The canopy is bounded from below by a reflective background with a surface albedo equals 0.18 at red and NIR bands. Ground cover varies with canopy LAI as LAI/7.0. Other parameters are as in Figure 6.
Remotesensing 10 01632 g008
Table 1. Characteristics of maize canopy and illumination condition of canopy BRF measurements.
Table 1. Characteristics of maize canopy and illumination condition of canopy BRF measurements.
Characteristics of the Maize CanopyRatios of Direct to Total Incident Radiation
Height (m)Row-Width (m)Plant-to-Plant Interval (m)LAI (/)Blue (436) (/)Green (547) (/)Red (700) (/)NIR (804) (/)
1.090.60.252.240.880.920.940.92
Table 2. Optical properties of maize leaf and soil background.
Table 2. Optical properties of maize leaf and soil background.
CategoriesWavelength (nm)Blue (436)Green (547)Red (700)NIR (804)
LeafNHR0.070.150.140.40
NHT0.020.160.140.50
SoilReflectance0.080.150.230.27

Share and Cite

MDPI and ACS Style

Yang, B.; Knyazikhin, Y.; Xie, D.; Zhao, H.; Zhang, J.; Wu, Y. Influence of Leaf Specular Reflection on Canopy Radiative Regime Using an Improved Version of the Stochastic Radiative Transfer Model. Remote Sens. 2018, 10, 1632. https://doi.org/10.3390/rs10101632

AMA Style

Yang B, Knyazikhin Y, Xie D, Zhao H, Zhang J, Wu Y. Influence of Leaf Specular Reflection on Canopy Radiative Regime Using an Improved Version of the Stochastic Radiative Transfer Model. Remote Sensing. 2018; 10(10):1632. https://doi.org/10.3390/rs10101632

Chicago/Turabian Style

Yang, Bin, Yuri Knyazikhin, Donghui Xie, Haimeng Zhao, Junqiang Zhang, and Yi Wu. 2018. "Influence of Leaf Specular Reflection on Canopy Radiative Regime Using an Improved Version of the Stochastic Radiative Transfer Model" Remote Sensing 10, no. 10: 1632. https://doi.org/10.3390/rs10101632

APA Style

Yang, B., Knyazikhin, Y., Xie, D., Zhao, H., Zhang, J., & Wu, Y. (2018). Influence of Leaf Specular Reflection on Canopy Radiative Regime Using an Improved Version of the Stochastic Radiative Transfer Model. Remote Sensing, 10(10), 1632. https://doi.org/10.3390/rs10101632

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop