Spectral Invariants Theory Research Articles

Estimation of canopy photon recollision probability from airborne laser scanning

The past two decades have witnessed the widespread application of the spectral invariant theory (p-theory) in modeling the shortwave radiation absorbed or scattered by vegetation. The basic principle of this theory is that the canopy reflectance is solely determined by the optical properties of leaves and the photon recollision probability p. As one of the spectral invariants, p serves as a critical bond between the spectral characteristics of the canopy at any wavelength and the reflectance, transmittance, or absorptance of the vegetation canopy, and is intimately associated with the solution of the radiative transfer equation. Currently, estimation of p-value is predominantly accomplished through canopy structural parameters, such as the spherically averaged silhouette to total area ratio (STAR¯) or Leaf Area Index (LAI). In this work, we provide an approach to estimate canopy photon recollision probability directly from Airborne Laser Scanning (ALS) data. We showed that the recollision probability can be interpreted as the interceptions of virtual rays emitted from sampling points within tree crowns described with turbid media, which avoided detailed reconstruction of leaf structures. The approach was evaluated with both virtual experiments as well as field measurements. The virtual experiments employed the Large-scalE remote Sensing data and image Simulation model (LESS) to simulate virtual ALS scanning data based on three RAdiation transfer Model Intercomparison (RAMI) actual canopies, with each divided into 25 segments, for which the p-values could be obtained through LESS. Our findings showed that p-values can be accurately estimated from ALS point clouds. Specifically, it exhibited great consistency with RMSE% less than 5% for broadleaved forest scenes, 25% for mixed forest scenes, and less than 10% for coniferous forest scenes in virtual experiments, respectively, and less than 25% compared to field measurements. This study displays the potential of ALS as a promising tool for forest structure characterization and short-wave radiation modeling in forest ecosystems.

Physically based illumination correction for sub-centimeter spatial resolution hyperspectral data

Vegetation biophysical- and chemical traits, defined on the basis of leaf area, can be retrieved from their spectral reflectance. Ultra high resolution hyperspectral images, such as ones collected from drones, allows measuring the spectra of individual leaves. The reflectance signal of such data is calibrated with respect to the top-of-canopy (TOC) irradiance, as the local illumination conditions on leaf surfaces are largely unknown and can vary significantly from the TOC conditions. We developed an inversion algorithm that uses the PROSPECT leaf radiative transfer model and the theory of spectral invariants to retrieve the actual leaf reflectance from TOC-calibrated hyperspectral images. Compared with more traditional canopy reflectance models, this approach accounts for the spatial variation in leaf-level irradiance visible in sub-centimeter-resolution images and is computationally more efficient. We used simulated and measured leaf and canopy reflectance data to validate the approach and found the retrieved leaf reflectances to match closely the actual reflectances (relative RMSD was 12% for simulated data on the average and below 10% for measured data). The proposed method provides an efficient approach for illumination correction, enabling reliable, physically based applications for monitoring vegetation biochemical and biophysical properties from ultra-high-resolution spectral imagery.

First validation of Earth Reflector Type Index (p) parameter from DSCOVR EPIC VESDR data product using Terrestrial Ecosystem Research Network observing sites in Australia

The spectral invariants theory predicts that amount of radiation scattered by a canopy in a given wavelength is a simple function of the leaf albedo and a spectrally invariant parameter p. This study provides the first assessment of p parameter obtained from the Earth Reflector Type Index (ERTI) provided in Deep Space Climate Observatory (DSCOVR) Earth Polychromatic Imaging Camera (EPIC) Vegetation Earth System Data Record (VESDR) product. Our results suggest that ERTI parameter can possibly deliver information about photon recollision probability, albeit with a risk of lower accuracy, even for conditions that differ from the original assumptions in the spectral invariants theory (non-reflecting soil, dense canopies). We verified the expected link between ERTI p parameter in the EPIC VESDR product, foliage clumping, and canopy cover with empirical data. Importantly, this was done with study sites in Australia that included types of forests that are still relatively unexplored with spectral invariants theory. Our findings contribute to understanding product uncertainty and spatial representativeness of the in situ measurements from the Terrestrial Ecosystem Research Network (TERN) SuperSites in Australia.

Improving the estimation of canopy structure using spectral invariants: Theoretical basis and validation

Directional area scattering factor (DASF) is a critical canopy structural parameter for vegetation monitoring as it provides an efficient tool for detection of vegetation functional traits and their dynamics. Current standard approach to estimate DASF from canopy bidirectional reflectance factor (BRF) is based on the assumption that in the weakly absorbing 710 to 790 nm spectral interval, leaf scattering does not change much with the concentration of dry matter and thus its variation can be neglected. This results in biased estimates of DASF and consequently leads to uncertainty in DASF-related applications. This study proposes a new approach to account for variations in concentrations of this biochemical constituent, which additionally uses the canopy BRF at 2260 nm. Analysis of the proposed approach based on radiative transfer models suggests significant decrease in relative root mean square error (rRMSE) by nearly a half for one- and three dimensional dense vegetated scenes, and the improvement is insensitive to the reflectivity of background. Robust refinement is also achieved in a larger actual forest scene. When the new approach is adopted to indoor multi-angular hyperspectral measurements on individual tree crowns, the relative error has also reduced, with more improvement seen in needle-leaf than in broadleaf species. Thus, the proposed DASF estimation approach outperforms the current one and can be used more reliably in DASF-related applications, such as retrieval of vegetation biochemical contents, grasping periodic changes of equatorial forests, and other ecological and environmental vegetation monitoring tasks based on spectral invariants theory.

Spectral invariants in ultra-high spatial resolution hyperspectral images

Current operational optical satellite-based Earth observation methods targeting at vegetation are optimised for the high- and medium-resolution satellites with pixels sizes starting at ten meters. This resolution is coarser than the typical size of a vegetation structural element, such as a tree crown, and much coarser that of scattering elements, such as leaves. For this reason, vegetation can be treated as a continuous medium and the variation in the local illumination conditions on individual leaves or tree crowns can be ignored. This does not hold anymore for very and ultra-high resolution imagery, obtained from new satellite systems or unmanned aerial vehicles, where individual tree crowns or even leaves can be discerned. We tested the applicability of the spectral invariant theory to this type of imagery for characterising the local illumination conditions on plant leaves using Monte Carlo ray tracing simulations. The simulations corroborated the direct link between the spectral invariant parameter ρ and the sunlit fraction of visible leaves, earlier alleged for hyperspectral remote sensing data based on mathematical considerations. The approach allowed us to separate the direct beam and multiple scattering irradiance components and provided intuitive interpretations of the recollision probability and canopy scattering coefficient computed for each image pixel.

A TIR forest reflectance and transmittance (FRT) model for directional temperatures with structural and thermal stratification

Land surface temperature (LST) is listed as an essential climate variable (ECV) and supports quantitative estimates of the energy budget while serving as a proxy for measuring the effects of climate change and extreme events. Forested areas are considered a major land unit impacted by temperature rise; therefore, thorough monitoring is mandatory. An accuracy assessment of the LST of forests must consider their directional anisotropy (DA). This latter can be well depicted by thermal infrared (TIR) radiative transfer models, but the problem is complex for forests because many of the shaded areas generate multiscale gradients of temperature. In this paper, we adapted a mature and widely used visible and near-infrared (VNIR) radiative transfer model called forest reflectance and transmittance (FRT) to enhance the characterization of the DA of forest temperature. In the FRT model, the vertical heterogeneity of the forest is quantified by using the discrete elements of multilayer scene components (i.e., the tree crown, trunk, understory vegetation, and soil), thus inferring vertical thermal gradients. The Planck function and spectral-invariant theory are considered to assess the thermal emissions of the scene components and their multiple scattering processes. The FRT model is validated using directional forest brightness temperatures (BT) measured from an unmanned aerial vehicle (UAV) and simulated by using the three-dimensional ray-tracing LESS (large-scale remote sensing data and image simulation framework over heterogeneous 3D scenes) model. The results show that FRT behaves reliably since the root mean square error (RMSE) is lower than 1.0 °C for UAV measurements obtained at 09:20 and 13:10 and with coefficients of determination (R2) larger than 0.74 and 0.56, respectively; these results are better than the simulated results by existing models. Moreover, the comparison with ray-tracing simulations was also deemed satisfactory. According to the analysis, large variations in BT DAs may appear for different forests and seasonal changes staged by structural and thermal stratification, thus indicating the necessity of using the FRT model for complex and dynamic forest canopies.

Simulation-Based Evaluation of the Estimation Methods of Far-Red Solar-Induced Chlorophyll Fluorescence Escape Probability in Discontinuous Forest Canopies

The escape probability of Solar-induced chlorophyll fluorescence (SIF) can be remotely estimated using reflectance measurements based on spectral invariants theory. This can then be used to correct the effects of canopy structure on canopy-leaving SIF. However, the feasibility of these estimation methods is untested in heterogeneous vegetation such as the discontinuous forest canopy layer under evaluation here. In this study, the Discrete Anisotropic Radiative Transfer (DART) model is used to simulate canopy-leaving SIF, canopy total emitted SIF, canopy interceptance, and the fraction of absorbed photosynthetically active radiation (fAPAR) in order to evaluate the estimation methods of SIF escape probability in discontinuous forest canopies. Our simulation results show that the normalized difference vegetation index (NDVI) can be used to partly eliminate the effects of background reflectance on the estimation of SIF escape probability in most cases, but fails to produce accurate estimations if the background is partly or totally covered by vegetation. We also found that SIF escape probabilities estimated at a high solar zenith angle have better estimation accuracy than those estimated at a lower solar zenith angle. Our results show that additional errors will be introduced to the estimation of SIF escape probability with the use of satellite products, especially when the product of leaf area index (LAI) and clumping index (CI) was underestimated. In other results, fAPAR has comparable estimation accuracy of SIF escape probability when compared to canopy interceptance. Additionally, fAPAR for the entire canopy has better estimation accuracy of SIF escape probability than fPAR for leaf only in sparse forest canopies. These results help us to better understand the current estimation results of SIF escape probability based on spectral invariants theory, and to improve its estimation accuracy in discontinuous forest canopies.

Quantifying leaf optical properties with spectral invariants theory

Leaf optical spectra reflect the combination of leaf biochemical, morphological and physiological properties, and play an important role in many ecological and Earth system processes. Radiative transfer models are widely used to simulate leaf spectra by quantifying photon transfer processes of reflection, transmission and absorption within a plant leaf. Recent advances in spectral invariants theory offer a unique and efficient approach for modeling the canopy-scale radiative transfer processes, but remain underexplored for applications at the leaf scale. In this study, we developed a leaf-scale optical property model based on the spectrally invariant properties (leaf-SIP) of plant leaves. Similar to the canopy-scale model, the leaf-SIP model decouples leaf-scale radiative transfer process into two parts: wavelength-dependent contribution from leaf chemical components and wavelength-independent contribution from leaf structures, described by two spectrally invariant parameters (i.e., a photon recollision probability p and a scattering asymmetry parameter q). We implemented the leaf-SIP model by parameterizing p and q with a measurable leaf morphological trait, the leaf mass per area (LMA). We evaluated the performance of the leaf-SIP model with two in situ datasets (i.e., LOPEX and ANGERS) and the widely used PROSPECT leaf optical model. The results show that the leaf spectra simulated by the leaf-SIP model agreed well with in situ datasets and the simulations of the PROSPECT model, with a small root mean squared error (RMSE), bias, and high coefficients of determination (R2) of 0.026, 0.035, 0.95 and 0.037, 0.049, 0.91 for leaf reflectance and leaf transmittance, respectively. Our results also show that the leaf-SIP model can be used with measured leaf spectra to accurately estimate several key leaf functional traits, such as the leaf chlorophyll content, equivalent water thickness, and LMA. The leaf-SIP model provides an efficient and physical way of accurately simulating leaf spectra and retrieving key leaf functional traits from hyperspectral measurements.

A radiative transfer model for solar induced fluorescence using spectral invariants theory

Solar Induced chlorophyll Fluorescence (SIF) shows promise as an approach for estimating gross primary production (GPP) remotely. However, sun-target-sensor geometry and within-canopy absorption of SIF can alter the relationship between measured SIF and GPP, because sensors can only retrieve some unknown fraction of the total emitted SIF. Radiative transfer models that allow for variation in canopy structure and sensor angles are therefore needed to properly interpret SIF measurements. Spectral invariants allow decoupling of the wavelength-independent canopy structure and the wavelength-dependent leaf and soil spectrum in the radiative transfer process. Here we develop a simple analytical Fluorescence Radiative Transfer model based on Escape and Recollision probability (FluorRTER) to investigate the impact of canopy structure and sun-target-sensor geometry on SIF emissions. SIF simulations using the FluorRTER model agreed well the one-dimensional Soil-Canopy Observation of Photochemistry and Energy balance (SCOPE) model and the three-dimensional Fluorescence model with Weighted Photon Spread (FluorWPS) model. The fractional vegetation cover (FVC) and clumping effect have a large influence the SIF emission of 3D discontinuous canopies. For a moderate solar zenith angle (30°) and a clumped canopy (FVC = 0.6), the difference between the directional observed SIF of a 3D discontinuous canopy and a 1D homogeneous canopy was as large as 43.2% and 38.4% for Photosystem I + II fluorescence at 685 nm and at 740 nm, respectively. By bridging the gap between observed SIF and total emitted SIF over 3D heterogeneous vegetation canopies, the FluorRTER model can assist with the angular normalization of SIF measurements and enable the more robust interpretation of how variations in SIF from directional and hemispherical in-situ, airborne and satellite observations relate to leaf and whole-canopy physiological processes.

Reduction of structural impacts and distinction of photosynthetic pathways in a global estimation of GPP from space-borne solar-induced chlorophyll fluorescence

Quantifying global photosynthesis remains a challenge due to a lack of accurate remote sensing proxies. Solar-induced chlorophyll fluorescence (SIF) has been shown to be a good indicator of photosynthetic activity across various spatial scales. However, a global and spatially challenging estimate of terrestrial gross primary production (GPP) based on satellite SIF remains unresolved due to the confounding effects of species-specific physical and physiological traits and external factors, such as canopy structure or photosynthetic pathway (C3 or C4). Here we analyze an ensemble of far-red SIF data from OCO-2 satellite and ground observations at multiple sites, using the spectral invariant theory to reduce the effects of canopy structure and to retrieve a structure-corrected total canopy SIF emission (SIFtotal). We find that the relationships between observed canopy-leaving SIF and ecosystem GPP vary significantly among biomes. In contrast, the relationships between SIFtotal and GPP converge around two unique models, one for C3 and one for C4 plants. We show that the two single empirical models can be used to globally scale satellite SIF observations to terrestrial GPP. We obtain an independent estimate of global terrestrial GPP of 129.56 ± 6.54 PgC/year for the 2015–2017 period, which is consistent with the state-of-the-art data- and process-oriented models. The new GPP product shows improved sensitivity to previously undetected ‘hotspots’ of productivity, being able to resolve the double-peak in GPP due to rotational cropping systems. We suggest that the direct scheme to estimate GPP presented here, which is based on satellite SIF, may open up new possibilities to resolve the dynamics of global terrestrial GPP across space and time.

Fluorescence Correction Vegetation Index (FCVI): A physically based reflectance index to separate physiological and non-physiological information in far-red sun-induced chlorophyll fluorescence

Sun-induced chlorophyll fluorescence (SIF) has been used to track vegetation photosynthetic activity for improving estimation of gross primary productivity (GPP) and detecting plant stress. There are both physical and physiological controls of SIF measured at the surface and retrieved from remote sensing including satellite observations. In order to accurately use SIF for monitoring of plant physiology, the effects of physically-based radiation processes related to leaf and canopy structure, notably photosynthetically active radiation (PAR) absorption and SIF scattering and re-absorption, must be characterized. In this study, we investigate both PAR absorption and SIF scattering processes and find that although it is difficult to quantify their effects individually by using just reflectance, the combined effects of the two processes can be well approximated by a reflectance index. This index, referred to as FCVI (Fluorescence Correction Vegetation Index), is defined as the difference between near-infrared (NIR) and broad-band visible (VIS, 400–700 nm) reflectance acquired under identical sun-canopy-observer geometry of the SIF measurements. The development of the index was based on the physical connection between reflectance and far-red SIF, which was revealed by using the spectral invariant theory. The utility of FCVI to correct far-red SIF for PAR absorption and scattering effects, thus improving the link to photosynthesis, was tested with data from: (i) a field experiment for a growing season; and (ii) a numerical experiment which included a number of scenarios generated by a radiative transfer model. For both the observations and simulations, the FCVI provided a promising estimate of the impact of the physically-based radiation processes on far-red SIF of moderately dense canopies (i.e., FCVI ≥ 0.18). Normalizing the TOC far-red SIF by both the incident PAR (iPAR) and the FCVI provided a good estimate of the far-red fluorescence emission efficiency of the canopies examined. This approach enhances our ability to generalize retrievals for vegetation processes as they change through natural growth phases and seasons. Taken together, far-red SIF and FCVI may enable the assessment of the light partitioning of vegetation canopies, an essential step to facilitate the use of far-red SIF for tracking physiological processes.

From Canopy‐Leaving to Total Canopy Far‐Red Fluorescence Emission for Remote Sensing of Photosynthesis: First Results From TROPOMI

AbstractSolar‐induced chlorophyll fluorescence (SIF) from the TROPOspheric Monitoring Instrument (TROPOMI), which has substantially improved spatial and temporal resolutions, will improve the global estimations of gross primary production (GPP) than previous satellite SIF data. However, the canopy‐leaving SIF observed by sensors (SIFobs) represents only a portion of the total canopy SIF emission (SIFtotal). This portion is sensitive to the canopy structure and observation direction, resulting in uncertainties in GPP estimations using SIFobs. Here we used the spectral invariant theory to derive global soil‐resistant SIFtotal (SIFtotal‐SR) from TROPOMI SIFobs and evaluated the SIFtotal‐SR performance in estimating GPP. We found that the clear differences in SIFobs between needleleaf forest and crops diminished for SIFtotal‐SR. SIFtotal‐SR increased the coefficient of determination (R2) by 0.09 and 0.11 against the flux tower instantaneous and daily GPP, respectively. This derived SIFtotal‐SR can be used to develop more robust GPP models and better constrain carbon cycle models.

Empirical Methods for Remote Sensing of Nitrogen in Drylands May Lead to Unreliable Interpretation of Ecosystem Function

Nitrogen (N) has been linked to different ecosystem processes, and retrieving this important foliar biochemical constituent from remote sensing observations is of widespread interest. Since N is not explicitly represented in physically based radiative transfer models, empirical methods have been used as an alternative. The spectral bands selected during the calibration of empirical methods have been interpreted in the context of light-N interactions and, consequently, ecosystem function. The low amount of leaves on shrubs and their sparse distribution in drylands create an environment, in which the canopy structure and the bright background soil play an important role in the canopy total radiation budget. In this paper, we examine the assumption that removing the impact of canopy structure and soil will result in improved N retrieval using the empirical methods. We report the inconsistencies in the selection of spectral bands among the empirical approaches. Moreover, these methods are sensitive to the scale of the study and band transformations. Using the generalized theory of canopy spectral invariants, we found that at the canopy scale, a combination of canopy structure and soil dominates the total canopy radiation. At the plot scale, soil contributes up to 95% of the total reflectance. Correction for these two confounding factors leads to no correlation between N and vegetation reflectance at both scales. We conclude that while cross-validated predictive models may be statistically achievable, caution should be taken when interpreting their results in the context of N-light interactions and ecosystem function. Our approach can be extended to all terrestrial ecosystems with multiple layers of canopy and understory.

Bayesian inversion of a forest reflectance model using Sentinel-2 and Landsat 8 satellite images

The inversion of reflectance models is a generalizable tool to obtain estimates on forest biophysical parameters, such as leaf area index, with theoretically little information need from a study area, instead relying on the knowledge about physical processes in the forest radiation regime. The use of prior information can greatly improve the reflectance model inversion, however, the literature does not yet provide much information on the selection of priors and their influence on the inversion results. In this study, we used a Bayesian approach to invert the PARAS forest reflectance model and retrieve leaf area index from Sentinel-2 MSI and Landsat 8 OLI multispectral satellite images. The PARAS model is based on the theory of spectral invariants, which describes the influence of wavelength-independent parameters on forest radiative transfer. The Bayesian inversion approach is highly flexible, provides uncertainty quantification, and enables the explicit incorporation of prior knowledge into the inversion process. We found that the choice of prior information is crucial in inverting a forest reflectance model to predict leaf area index. Regularizing and informative priors for leaf area index strongly improved the predictions, relative to an uninformative prior, in that they counteracted the saturation effect of the optical signal occuring at high values for leaf area index. The predictions of leaf area index were more accurate for Landsat 8 than for Sentinel-2, due to potential inconsistencies in the visible bands of Sentinel-2 in our data, and the higher spectral resolution.

Decoupling Canopy Structure and Leaf Biochemistry: Testing the Utility of Directional Area Scattering Factor (DASF)

Biochemical properties retrieved from remote sensing data are crucial sources of information for many applications. However, leaf and canopy scattering processes must be accounted for to reliably estimate information on canopy biochemistry, carbon-cycle processes and energy exchange. A coupled leaf-canopy model based on spectral invariants theory has been proposed, that uses the so-called Directional Area Scattering Factor (DASF) to correct hyperspectral remote sensing data for canopy structural effects. In this study, the reliability of DASF to decouple canopy structure and biochemistry was empirically tested using simulated reflectance spectra modelled using a Monte Carlo Ray Tracing (MCRT) radiative transfer model. This approach allows all canopy and radiative properties to be specified a priori. Simulations were performed under idealised conditions of directional-hemispherical reflectance, isotropic Lambertian leaf reflectance and transmittance and sufficiently dense (high LAI) canopies with black soil where the impact of canopy background is negligible, and also departures from these conditions. It was shown that both DASF and total canopy scattering could be accurately extracted under idealised conditions using information from both the full 400–2500 nm spectral interval and the 710–790 nm interval alone, even given no prior knowledge of leaf optical properties. Departures from these idealised conditions: varying view geometry, bi-directional reflectance, LAI and soil effects, were tested. We demonstrate that total canopy scattering could be retrieved under conditions of varying view geometry and bi-directional reflectance, but LAI and soil effects were shown to reduce the accuracy with which the scattering can be modelled using the DASF approach. We show that canopy architecture, either homogeneous or heterogeneous 3D arrangements of canopy scattering elements, has important influences over DASF and consequently the accuracy of retrieval of total canopy scattering. Finally, although DASF and total canopy scattering could be retrieved to within 2.4% of the modelled total canopy scattering signal given no prior knowledge of leaf optical properties, spectral invariant parameters were not accurately retrieved from the simulated signal. This has important consequences since these parameters are quite widely used in canopy reflectance modelling and have the potential to help derive new, more accurate canopy biophysical information. Understanding and quantifying the limitations of the DASF approach as we have done here, is an important step in allowing the wider use of these methods for decoupling canopy structure and biochemistry.

An Interplay between Photons, Canopy Structure, and Recollision Probability: A Review of the Spectral Invariants Theory of 3D Canopy Radiative Transfer Processes

Earth observations collected by remote sensors provide unique information to our ever-growing knowledge of the terrestrial biosphere. Yet, retrieving information from remote sensing data requires sophisticated processing and demands a better understanding of the underlying physics. This paper reviews research efforts that lead to the developments of the stochastic radiative transfer equation (RTE) and the spectral invariants theory. The former simplifies the characteristics of canopy structures with a pair-correlation function so that the 3D information can be succinctly packed into a 1D equation. The latter indicates that the interactions between photons and canopy elements converge to certain invariant patterns quantifiable by a few wavelength independent parameters, which satisfy the law of energy conservation. By revealing the connections between plant structural characteristics and photon recollision probability, these developments significantly advance our understanding of the transportation of radiation within vegetation canopies. They enable a novel physically-based algorithm to simulate the “hot-spot” phenomenon of canopy bidirectional reflectance while conserving energy, a challenge known to the classic radiative transfer models. Therefore, these theoretical developments have a far-reaching influence in optical remote sensing of the biosphere.

Downscaling of solar-induced chlorophyll fluorescence from canopy level to photosystem level using a random forest model

Solar-induced chlorophyll fluorescence (SIF), an electromagnetic signal that can potentially indicate vegetation photosynthetic activity, can be retrieved from ground-based, airborne and satellite measurements. However, due to the scattering and re-absorption effects inside the leaves and canopy, SIF measured at the canopy level is only a small part of the total SIF emission at the photosystem level. Therefore, a downscaling mechanism of SIF from the canopy level to the photosystem level is important for better understanding the relationship between SIF and the vegetation gross primary production (GPP). In this study, firstly, we analyzed the canopy scattering effects using a simple parameterization model based on the spectral invariant theory. The probability for SIF photons to escape from the canopy was found to be related to the anisotropic spectral reflectance, canopy interception of the upward solar radiation, and leaf absorption. An empirical approach based on a Random Forest (RF) regression algorithm was applied to downscale SIF constrained by the red, red-edge and far-red anisotropic reflectance. The RF was trained using simulations conducted with the Soil Canopy Observation, Photochemistry and Energy fluxes (SCOPE) model. The performance of the SIF downscaling method was evaluated with SCOPE and Discrete Anisotropic Radiative Transfer (DART) model simulations, ground measurements and airborne data. Results show that estimated SIF at the photosystem level matches well with simulated reference data, and the relationship between SIF and photosynthetically active radiation absorbed by chlorophyll is improved by SIF downscaling. This finding in combination with other evaluation criteria suggests the downscaling of canopy SIF as an efficient strategy to normalize species dependent effects of canopy structure and varying solar-view geometries. Based on our results for the SIF-APAR relationship, we expect that such normalization approaches can be helpful to improve estimates of photosynthesis using remote sensing measurements of SIF.

Data synergy between leaf area index and clumping index Earth Observation products using photon recollision probability theory

Clumping index (CI) is a measure of foliage aggregation relative to a random distribution of leaves in space. The CI can help with estimating fractions of sunlit and shaded leaves for a given leaf area index (LAI) value. Both the CI and LAI can be obtained from global Earth Observation data from sensors such as the Moderate Resolution Imaging Spectrometer (MODIS). Here, the synergy between a MODIS-based CI and a MODIS LAI product is examined using the theory of spectral invariants, also referred to as photon recollision probability (‘p-theory’), along with raw LAI-2000/2200 Plant Canopy Analyzer data from 75 sites distributed across a range of plant functional types. The p-theory describes the probability (p-value) that a photon, having intercepted an element in the canopy, will recollide with another canopy element rather than escape the canopy. We show that empirically-based CI maps can be integrated with the MODIS LAI product. Our results indicate that it is feasible to derive approximate p-values for any location solely from Earth Observation data. This approximation is relevant for future applications of the photon recollision probability concept for global and local monitoring of vegetation using Earth Observation data.

Generating Global Products of LAI and FPAR From SNPP-VIIRS Data: Theoretical Background and Implementation

Leaf area index (LAI) and fraction of photosynthetically active radiation (FPAR) absorbed by vegetation have been successfully generated from the Moderate Resolution Imaging Spectroradiometer (MODIS) data since early 2000. As the Visible Infrared Imaging Radiometer Suite (VIIRS) instrument onboard, the Suomi National Polar-orbiting Partnership (SNPP) has inherited the scientific role of MODIS, and the development of a continuous, consistent, and well-characterized VIIRS LAI/FPAR data set is critical to continue the MODIS time series. In this paper, we build the radiative transfer-based VIIRS-specific lookup tables by achieving minimal difference with the MODIS data set and maximal spatial coverage of retrievals from the main algorithm. The theory of spectral invariants provides the configurable physical parameters, i.e., single scattering albedos (SSAs) that are optimized for VIIRS-specific characteristics. The effort finds a set of smaller red-band SSA and larger near-infrared-band SSA for VIIRS compared with the MODIS heritage. The VIIRS LAI/FPAR is evaluated through comparisons with one year of MODIS product in terms of both spatial and temporal patterns. Further validation efforts are still necessary to ensure the product quality. Current results, however, imbue confidence in the VIIRS data set and suggest that the efforts described here meet the goal of achieving the operationally consistent multisensor LAI/FPAR data sets. Moreover, the strategies of parametric adjustment and LAI/FPAR evaluation applied to SNPP-VIIRS can also be employed to the subsequent Joint Polar Satellite System VIIRS or other instruments.

Estimation of leaf area index and its sunlit portion from DSCOVR EPIC data: Theoretical basis

This paper presents the theoretical basis of the algorithm designed for the generation of leaf area index and diurnal course of its sunlit portion from NASA's Earth Polychromatic Imaging Camera (EPIC) onboard NOAA's Deep Space Climate Observatory (DSCOVR). The Look-up-Table (LUT) approach implemented in the MODIS operational LAI/FPAR algorithm is adopted. The LUT, which is the heart of the approach, has been significantly modified. First, its parameterization incorporates the canopy hot spot phenomenon and recent advances in the theory of canopy spectral invariants. This allows more accurate decoupling of the structural and radiometric components of the measured Bidirectional Reflectance Factor (BRF), improves scaling properties of the LUT and consequently simplifies adjustments of the algorithm for data spatial resolution and spectral band compositions. Second, the stochastic radiative transfer equations are used to generate the LUT for all biome types. The equations naturally account for radiative effects of the three-dimensional canopy structure on the BRF and allow for an accurate discrimination between sunlit and shaded leaf areas. Third, the LUT entries are measurable, i.e., they can be independently derived from both below canopy measurements of the transmitted and above canopy measurements of reflected radiation fields. This feature makes possible direct validation of the LUT, facilitates identification of its deficiencies and development of refinements. Analyses of field data on canopy structure and leaf optics collected at 18 sites in the Hyytiälä forest in southern boreal zone in Finland and hyperspectral images acquired by the EO-1 Hyperion sensor support the theoretical basis.