Estimating forest variables from top-of-atmosphere radiance satellite measurements using coupled radiative transfer models

Traditionally, it is necessary to pre-process remote sensing data to obtain top of canopy (TOC) reflectances before applying physically-based model inversion techniques to estimate forest variables. Corrections for atmospheric, adjacency, topography, and surface directional effects are applied sequentially and independently, accumulating errors into the TOC reflectance data, which are then further used in the inversion process. This paper presents a proof of concept for demonstrating the direct use of measured top-of-atmosphere (TOA) radiance data to estimate forest biophysical and biochemical variables, by using a coupled canopy-atmosphere radiative transfer model. Advantages of this approach are that no atmospheric correction is needed and that atmospheric, adjacency, topography, and surface directional effects can be directly and more accurately included in the forward modelling.In the case study, we applied both TOC and TOA approaches to three Norway spruce stands in Eastern Czech Republic. We used the SLC soil-leaf-canopy model and the MODTRAN4 atmosphere model. For the TOA approach, the physical coupling between canopy and atmosphere was performed using a generic method based on the 4-stream radiative transfer theory which enables full use of the directional reflectance components provided by SLC. The method uses three runs of the atmosphere model for Lambertian surfaces, and thus avoids running the atmosphere model for each new simulation. We used local sensitivity analysis and singular value decomposition to determine which variables could be estimated, namely: canopy cover, fraction of bark, needle chlorophyll, and dry matter content. TOC and TOA approaches resulted in different sets of estimates, but had comparable performance. The TOC approach, however, was at its best potential because of the flatness and homogeneity of the area. On the contrary, the capacities of the TOA approach would be better exploited in heterogeneous rugged areas. We conclude that, having similar performance, the TOA approach should be preferred in situations where minimizing the pre-processing is important, such as in data assimilation and multi-sensor studies.     

Multi-angular reflectance properties of a hemiboreal forest: An analysis using CHRIS PROBA data

Forest types differ in their hyperspectral anisotropy patterns mainly due to species-specific geometrical structure, spatial arrangement of canopies and subsequent shadow patterns. This paper examines the multi-angular, hyperspectral reflectance properties of typical hemiboreal forests during summer time using three simultaneous CHRIS PROBA (mode 3) scenes and stand inventory data from the Järvselja Training and Experimental Forestry District in southeastern Estonia. We investigated the magnitude and reasons for the differences in the anisotropy patterns of deciduous and coniferous stands at three backward viewing angles. A forest reflectance model (FRT) was used as a tool to provide a theoretical basis to the discussion, and to estimate the directional contribution of scattering from crowns and ground to total stand reflectance for the two forest types. The FRT model simulated successfully the HDRF (hemispherical–directional reflectance factor) curves of the study stands to match those obtained from the CHRIS image, yet it produced a smaller and less wavelength-dependent angular reflectance effect than was observed in the satellite image. The main results of this study provide new information for separating the spectral contribution of the forest floor (or understory layer) from the tree canopy layer: (1) the red edge domain was identified to have the largest contribution from forest understory, and (2) the more oblique the viewing angle, the smaller the contribution from the understory. In addition, coniferous stands were observed to have a specific angular effect at the red and red edge domain, possibly as a result of the hierarchical structure and arrangement of coniferous canopies.     

Validation of the forest radiative transfer model FRT

The components of the forest radiative transfer model FRT developed at Tartu Observatory, Estonia are compared to the results of measurements of forest canopy downward radiance under forest canopy. Measurements were performed with a hemispheric-view imaging CCD-radiometer which was specially designed for this task. A thorough study of metrological properties of the radiometer was carried out and the respective preprocessing algorithms were created. The angular distribution of view probabilities and NIR radiances for tree crowns, trunks and canopy gaps are measured and compared to the model simulations for two coniferous and a broadleaf forests. Generally, the model reproduced the angular courses of component probabilities and NIR radiances rather well, however, in some cases problems with both the absolute levels of radiance and the angular course arose.     

Atmospheric correction of PROBA/CHRIS data in an urban environment

The Compact High Resolution Imaging Spectrometer (CHRIS) is an imaging spectrometer onboard the European Space Agency (ESA) Project for On-Board Autonomy (PROBA) satellite. However, it has been shown that CHRIS presents some miscalibration trends over the spectral region covered. This paper reports a practical procedure for the atmospheric correction of CHRIS images based on field recalibration in an urban environment. In the first stage, the spectra of surface targets are measured and used to simulate the spectral radiance at the top of the atmosphere (TOA) for each channel and to determine the recalibration coefficients of the CHRIS images. In the second stage, two methods for atmospheric correction are examined: the radiative transfer model (RTM) and the improved dark-object subtraction (IDOS) method. For comparison purposes, the empirical line method (ELM) is also evaluated. The accuracy assessment shows that the RTM with the Moderate Resolution Transmittance (MODTRAN) code provides the most accurate atmospheric correction for the multiangular CHRIS images when using the proposed procedure.     

A dataset for the validation of reflectance models

Three mature stands at the forest test site Järvselja, Estonia were extensively measured for using as a validation dataset for heterogeneous canopy reflectance models. In order to enable the reconstruction of the 3-D architecture of these 100 × 100 m² test plots, individual tree positions and crown dimensions were inventoried. In addition, leaf, needle, stem bark and branch bark visible and near-infrared (VNIR) reflectance spectra, and VNIR reflectance spectra of ground vegetation were measured. This in situ dataset is supported by atmospherically and radiometrically corrected Mode 3 CHRIS reflectance spectra for three view directions, and top of canopy VNIR nadir spectra from airborne measurements. Details of measurements, instruments in use, data processing, and access to data are described in a technical report which is available on-line.     

Coupled soil-leaf-canopy and atmosphere radiative transfer modeling to simulate hyperspectral multi-angular surface reflectance and TOA radiance data

Coupling radiative transfer models for the soil background and vegetation canopy layers is facilitated by means of the four-stream flux interaction concept and use of the adding method. Also the coupling to a state-of-the-art atmospheric radiative transfer model like MODTRAN4 can be established in this way, thus enabling the realistic simulation of top-of-atmosphere radiances detected by space-borne remote sensing instruments. Possible applications of coupled modeling vary from mission design to parameter retrieval and data assimilation. This paper introduces a modified Hapke soil BRDF model, a robust version of the PROSPECT leaf model, and a modernized canopy radiative transfer model called 4SAIL2. The latter is a hybrid two-layer version of SAIL accommodating horizontal and vertical heterogeneities, featuring improved modeling of the hot spot effect and output of canopy absorptances. The integrated model is simply called SLC (soil-leaf-canopy) and has been implemented as a speed-optimized Windows DLL which allows efficient use of computer resources even when simulating massive amounts of hyperspectral multi-angular observations. In this paper various examples of possible model output are shown, including simulated satellite image products. First validation results have been obtained from atmospherically corrected hyperspectral multi-angular CHRIS-PROBA data of the Upper Rhine Valley in Germany.     

航天成像光谱仪CHRIS在内陆水质监测中的应用

欧空局2001年10月22日成功发射的PROBA卫星上搭载了紧密型高分辨率成像光谱仪（CHRIS），它可以提供高光谱分辨率、高空间分辨率和多角度的遥感数据，它代表了新一代的地球观测数据源。CHRIS有5种工作模式，其中模式2是专门为水体研究而设计的，它在400～1 050 nm的可见光至近红外有18个波段，每个波段数据的空间分辨率是17 m。CHRIS数据的高光谱分辨率、高空间分辨率和多时相覆盖的特点为内陆水质监测提供了有利条件。为了验证CHRIS在内陆水质监测中的具体应用，在太湖梅梁湾开展了水面综合试验，在梅梁湾均匀分布的14个水面采样点分别测量了水面光谱和水质参数。利用这些数据，同时结合CHRIS数据的光谱特征，建立了叶绿素浓度反演半经验模型，应用于CHRIS图像反演了太湖梅梁湾的叶绿素浓度分布图，并取得了较好的结果。最后指出CHRIS数据不但在内陆水质监测中具有巨大潜力，而且CHRIS遥感器是今后内陆水质监测卫星遥感器的典范。     

Improving the estimation of leaf area index by using remotely sensed NDVI with BRDF signatures

A new vegetation index, the Normalized Hotspot-signature Vegetation Index (NHVI), is proposed for a better quantitative estimation of leaf area index (LAI) than with the remotely sensed normalized difference vegetation index (NDVI), especially in the boreal forest. To obtain this new index, the Hotspot-Dark-spot index (HDS) (Lacaze et al., 2002) was introduced. HDS is calculated by the difference between the strongest vector (hotspot) and the weakest vector (dark-spot) of bi-directional reflectance, a given tract of vegetation returns in the reflecting solar position, and the geometric structure of the vegetation canopy, which are poorly represented by NDVI alone. The validity of NHVI was statistically tested using two field data sets of multi-angular observations and LAI from the boreal forests of Canada; one set was our own observations, and the other was from the Boreal Ecosystem-Atmosphere Study (BOREAS). The range of linear correspondence of NHVI with LAI is much wider than that of NDVI alone, indicating significant representation of leaf biomass in the canopy geometry captured by HDS. With the technical innovation of multi-angular remote-sensing and kernel-driven models in the future, this index has the potential to provide a more accurate evaluation of regional and global LAIs.     

Estimating crop chlorophyll content with hyperspectral vegetation indices and the hybrid inversion method

A hybrid inversion method was developed to estimate the leaf chlorophyll content (LCC) and canopy chlorophyll content (CCC) of crops. Fifty hyperspectral vegetation indices (VIs), such as the photochemical reflectance index (PRI) and canopy chlorophyll index (CCI), were compared to identify the appropriate VIs for crop LCC and CCC inversion. The hybrid inversion models were then generated from different modelling methods, including the curve-fitting and least squares support vector regression (LS-SVR) and random forest regression (RFR) algorithms, by using simulated Compact High Resolution Imaging Spectrometer (CHRIS) datasets that were generated by a radiative transfer model. Finally, the remote-sensing mapping of a CHRIS image was completed to test the inversion accuracy. The results showed that the remote-sensing mapping of the CHRIS image yielded an accuracy of R² = 0.77 and normalized root mean squared error (NRMSE) = 17.34% for the CCC inversion, and an accuracy of only R² = 0.33 and NRMSE = 26.03% for LCC inversion, which indicates that the remote-sensing technique was more appropriate for obtaining chlorophyll content at the canopy scale (CCC) than at the leaf scale (LCC). The estimated results of various VIs and algorithms suggested that the PRI and CCI were the optimal VIs for LCC and CCC inversion, respectively, and RFR was the optimal method for modelling.     

Top-of-atmosphere flux retrievals from CERES using artificial neural networks

The Clouds and the Earth's Radiant Energy System (CERES) instruments on the Terra spacecraft provide accurate shortwave (SW), longwave (LW) and window (WN) region top-of-atmosphere (TOA) radiance measurements from which TOA radiative flux values are obtained by applying Angular Distribution Models (ADMs). These models are developed empirically as functions of the surface and cloud properties provided by coincident high-resolution imager measurements over CERES field-of-view. However, approximately 5.6% of the CERES/Terra footprints lack sufficient imager information for a reliable scene identification. To avoid any systematic biases in regional mean radiative fluxes, it is important to provide TOA fluxes for these footprints. For this purpose, we apply a feedforward error-backpropagation Artificial Neural Network (ANN) technique to reproduce CERES/Terra ADMs relying only on CERES measurements. All-sky ANN-based angular distribution models are developed for 10 surface types separately for shortwave, longwave and window TOA flux retrievals. To optimize the ANN performance, we use a partially connected first hidden neuron layer and compact training sets with reduced data noise. We demonstrate the performance of the ANN-based ADMs by comparing TOA fluxes inferred from ANN and CERES anisotropic factors. The global annual average bias in ANN-derived fluxes relative to CERES is less than 0.5% for all ANN scene types. The maximum bias occurs over sea ice and permanent snow surfaces. For all surface types, instantaneous ANN-derived TOA fluxes are self-consistent in viewing zenith angle to within 9% for shortwave, 3.5% and 3% longwave daytime and nighttime, respectively.     

Effective tree scattering and opacity at L-band

This paper investigates vegetation effects at L-band by using a first-order radiative transfer (RT) model and truck-based microwave measurements over natural conifer stands to assess the applicability of the τ ? ω (tau–omega) model over trees. The tau–omega model is a zero-order RT solution that accounts for vegetation effects with two vegetation parameters (vegetation opacity and single-scattering albedo), which represent the canopy as a whole. This approach inherently ignores multiple-scattering effects and, therefore, has a limited validity depending on the level of scattering within the canopy. The fact that the scattering from large forest components such as branches and trunks is significant at L-band requires that zero-order vegetation parameters be evaluated (compared) along with their theoretical definitions to provide a better understanding of these parameters in the retrieval algorithms as applied to trees.This paper compares the effective vegetation opacities, computed from multi-angular pine tree brightness temperature data, against the results of two independent approaches that provide theoretical and measured optical depths. These two techniques are based on forward scattering theory and radar corner reflector measurements, respectively. The results indicate that the effective vegetation opacity values are smaller than but of similar magnitude to both radar and theoretical estimates. The effective opacity of the zero-order model is thus set equal to the theoretical opacity and an explicit expression for the effective albedo is then obtained from the zero- and first-order RT model comparison. The resultant albedo is found to have a similar magnitude as the effective albedo value obtained from brightness temperature measurements. However, both are less than half of the single-scattering albedo estimated using the theoretical calculations (0.5?0.6 for tree canopies at L-band). This lower observed effective albedo balances the scattering darkening effect of the large theoretical single-scattering albedo with a first-order multiple-scattering contribution. The retrieved effective albedo is different from theoretical definitions and not the albedo of single forest elements anymore, but it becomes a global parameter, which depends on all the processes taking place within the canopy, including multiple-scattering and canopy ground interaction.     

Estimation of vegetation clumping index using MODIS BRDF data

The foliage clumping index quantifies the degree of the deviation of leaf spatial distribution in the canopy from the random case. It is of comparable importance for ecological models as the leaf area index for quantifying radiation interception and distribution in plant canopies. Previously, an improved angular index named normalized difference between hotspot and darkspot was proposed for retrieving the clumping index using multi-angle remote sensing data. Global maps of clumping index have been derived successfully from multi-angular Polarization and Directionality of Earth Reflectance (POLDER) data at ～6 km resolution. In this article, we investigate whether it is feasible to derive the clumping index at 500 m resolution with the 16-day Moderate Resolution Imaging Spectroradiometer (MODIS) bidirectional reflectance distribution function model parameters product. The results are compared with an assembled set of field measurements from 63 different sites, covering five continents and diverse biomes.     

Remote sensing of ocean color: Assessment of the water-leaving radiance bidirectional effects on the atmospheric diffuse transmittance for SeaWiFS and MODIS intercomparisons

The portion of the radiance exiting the ocean and transmitted to the top of the atmosphere (TOA) in a particular direction depends on the angular distribution of the exiting radiance, not just the radiance exiting in the direction of interest. The diffuse transmittance t relates the water component of the TOA radiance to that exiting the water in the same direction. As such t is a property of the ocean–atmosphere system and not just the atmosphere. Its computation requires not only the properties of the atmosphere but the angular distribution of the exiting radiance as well. The latter is not known until a determination of the water properties can be made (which is the point of measuring the radiance in the first place). Because of this, it has been customary to assume an angular distribution (uniform upward radiance beneath the water surface) in the computation of t, which is referred to as t. However, it is known that replacing t with t can result in an error of several percent in the retrieved water-leaving radiance. Since the error depends on sun-viewing direction, this error could be particularly important when water-leaving radiance from two or more sensors in different orbits are compared. Even given an estimate of the angular distribution of the water-leaving radiance, computation of t using full radiative transfer theory in an image processing environment is not practical. Thus, we developed a first-order correction to t for bidirectional effects in the water-leaving radiance that captures much of the variability of t with viewing direction. The correction computed across a SeaWiFS scan line shows that a t-induced error of as much as 5–6% could occur near the edges of the scan; however, limiting the scan to polar viewing angles (θ) < 60° reduces the error to ~ 1%. Direct application to SeaWIFS and MODIS (AQUA) suggests that the bidirectionally-induced error in the diffuse transmittance will result in an error less than about 1% in the comparison of their normalized water-leaving radiances, as long as θ is less than about 60°. We conclude that, given this constraint, the normalized water-leaving retrievals from these two sensors at a given location can be merged without regard for the bidirectionally-induced error in the diffuse transmittance, as the resulting uncertainty is well below that from other sources. It is important to note that this result is likely to apply to any other polar-orbiting sensor with equatorial crossing times (similar to SeaWiFS and MODIS) between 1030 and 1330 h (local time).     

Detection of foliage conditions and disturbance from multi-angular high spectral resolution remote sensing

Disturbance of forest ecosystems, an important component of the terrestrial carbon cycle, has become a focus of research over recent years, as global warming is about to increase the frequency and severity of natural disturbance events. Remote sensing offers unique opportunities for detection of forest disturbance at multiple scales; however, spatially and temporally continuous mapping of non-stand replacing disturbance remains challenging. First, most high spatial resolution satellite sensors have relatively broad spectral ranges with bandwidths unsuitable for detection of subtle, stress induced, features in canopy reflectance. Second, directional and background reflectance effects, induced by the interactions between the sun-sensor geometry and the observed canopy surface, make up-scaling of empirically derived relationships between changes in spectral reflectance and vegetation conditions difficult. Using an automated tower based spectroradiometer, we analyse the interactions between canopy level reflectance and different stages of disturbance occurring in a mountain pine beetle infested lodgepole pine stand in northern interior British Columbia, Canada, during the 2007 growing season. Directional reflectance effects were modelled using a bidirectional reflectance distribution function (BRDF) acquired from high frequency multi-angular spectral observations. Key wavebands for observing changes in directionally corrected canopy spectra were identified using discriminant analysis and highly significant correlations between canopy reflectance and field measured disturbance levels were found for several broad and narrow waveband vegetation indices (for instance, r²_NDVI = 0.90; r²_CHL3 = 0.85; p < 0.05). Results indicate that multi-angular observations are useful for extraction of disturbance related changes in canopy reflectance, in particular the temporally and spectrally dense data detected changes in chlorophyll content well. This study will help guide and inform future efforts to map forest health conditions at landscape and over increasingly coarse scales.     

A vector radiative transfer model for sea-surface salinity retrieval from space: a non-raining case

ABSTRACT

Sea-surface salinity (SSS) can be measured from space using a microwave sensor. However, achieving the desired accuracy in SSS retrieval is challenging due to the lower sensitivity of the brightness temperature to SSS especially at low sea-surface temperature conditions. The retrieval accuracy can be further degraded due to the atmospheric and sea-surface effects (including emission and reflection), which require more accurate correction methods based on the radiative transfer model. In this article, a vector radiative transfer model (VRTM) was developed based on a matrix operator method that considers the ocean–atmosphere system under non-raining conditions. The results from this model were compared with measurement data provided by the Soil Moisture and Ocean Salinity (SMOS) satellite sensor and the results from two other RT models (RT4 model and a forward model of the European Space Agency, ESA). Statistical evaluation of these results revealed that estimation errors of top of atmosphere (TOA) radiance by the VRTM model was less than 0.3% as compared to the RT4 model results. The difference of the brightness temperatures predicted by the VRTM model and measured by the SMOS was within 1.5 K which was better than the ESA’s forward model predictions. These results suggest that the VRTM is relatively more accurate and has high computational efficiency for simulating the TOA brightness temperature for various scientific research and remote-sensing applications.     

Separating physiologically and directionally induced changes in PRI using BRDF models

Monitoring of photosynthetic efficiency (ε) over space and time is a critical component of climate change research as it is a major determinant of the amount of carbon accumulated by terrestrial ecosystems. While the past decade has seen progress in the remote estimation of ε at the leaf, canopy and stand level using the photochemical reflectance index PRI (based on the normalized difference of reflectance at 531 and 570 nm), little is known about the temporal and spatial requirements for up-scaling PRI to landscape and global levels using satellite observations. One potential way to investigate these requirements is using automated tower-based remote sensing platforms, which observe stand level reflectance with high spatial, temporal, and spectral resolution. Prediction of ε from PRI diurnally or over a full year requires observations of canopy reflectance over multiple view and sun-angles. As a result, these observations are subject to directional reflectance effects which can be interpreted in terms of the bidirectional reflectance distribution function (BRDF) using semi-empirical kernel driven models. These semi-empirical models use a combination of physically based BRDF shapes and empirical observations to standardize multi-angular observations to a common viewing and illumination geometry. Directional reflectance effects are thereby modeled as a linear superposition of mathematical kernels, representing the bi-direction variation in reflectance from isotropic, geometric, and volumetric scattering components of the vegetation canopy. However, because variations in plant physiological conditions can also introduce bidirectional reflectance variations, we introduce an approach to separate bidirectional effects arising purely from plant physiological status from other effects by stratifying PRI observations into categories based on environmental conditions for which the expected physiological variability is low. Within each of these PRI strata, the derived physically based BRDF shapes were used to standardize multi-angular PRI measurements to a common viewing and illumination geometry. The method significantly enhanced the relationship found between PRI and ε (from r² = 0.38 for the directionally uncorrected case to r² = 0.82 for the directionally corrected case) from data measured continuously over the course of 1 year over an evergreen conifer forest using an automated platform. Results show that isotropic PRI scattering is highly correlated to changes in ε, while geometric scattering can be related to canopy level shading. Instrumentation and approaches such as the one demonstrated in this study may be integrated into current efforts aiming at predicting ε at global scales using satellite observations.     

地形因子对森林冠层BRF模拟的影响分析

针对利用计算机模拟模型进行森林冠层BRF模拟时,地形因子——高程、坡度和坡向对森林冠层接收太阳辐射的影响等问题,本文探讨了利用RGM(Radiosity-Graphics combined Model)模型进行森林冠层反射特性的研究,研究了将森林生长ZELIG模型和L-系统结合起来构建三维森林场景,并采用RGM模型模拟森林冠层的二向性反射特性。研究结果表明地形因子对BRF分布影响非常明显。     

Relationships between polarimetric SAR backscatter and forest canopy and sub-canopy biophysical properties

This paper highlights the potential of multiwaveband polarimetric SAR data for the estimation of both canopy (percentage canopy closure) and sub-canopy (stem biomass) biophysical variables of a Sitka spruce forest in upland Wales. Stand stem biomass was estimated using forest survey data on diameter at breast height (DBH) and tree height from 0.01 ha plots. Photographs of the forest canopy were taken using a camera fitted with a wide-angle fisheye lens from a number of locations within a stand. The photographs were later digitized and estimates of stand percentage canopy closure were derived using image processing software. It was found that C-band HV and VV, and L-band HV and VV polarization backscatter were significantly related to stem biomass. There was no sensitivity to percentage canopy closure using single polarization backscatter but highly significant relationships were obtained using ratios of single polarization backscatter and variables derived from the polarization signatures. The strong correlations between C-band backscatter and stem biomass indicated a relationship between the structure of the top crown layer and sub-canopy biomass.     

Spatial analysis of radiometric fractions from high-resolution multispectral imagery for modelling individual tree crown and forest canopy structure and health

Research was conducted in a forest adjacent to an abandoned acid mine tailings site to assess forest structural health using high spatial and spectral resolution digital camera imagery. Conventional approaches to this problem involve the use image spectral information, basic spectral transformations, or occasionally spatial transformations of image brightness. This research introduces fractional textures and semivariance analysis of image fractions. They were integrated with conventional image measures in stepwise multiple regression modelling of forest structure (canopy and crown closure, stem density, tree height, crown size) and health (a visual stress index). The goal was to conduct a relative comparison of the potential of the various image variable types in modelling of forest structure and health. Analysis was conducted for both canopy (crowns and shadows) and individual tree crown sample data sets extracted from 10 nm bandwidth spectral bands at three resolutions (0.25, 0.5, 1.0 m). Spatial transformations (texture, semivariogram range) of image brightness (DN) and image fractions (IF) were consistently the most significant and first entered variables in the best models of the forest parameters. At the canopy-scale, despite a limited number of available plots (6), stable models were produced that demonstrated the potential for spatially transformed variables. Semivariogram range explained 88% of the total variation of 9 of the 18 models and represented 56% of the variables used in all models while texture variables explained 51% of model variance in 8 of the 18 models and represented 40% of the variables used. At the tree crown scale (n=31), 88% of the total variation of six of eight models was explained by texture variables and 6% by semivariogram variables. DN and IF variables that were not spatially transformed contributed little to the models at both scales. They represented 4% and 6%, respectively, of the variables used in all models. Spatial information in image fractions and image brightness has proven to be more significant than spectral information in these analyses. Of the spatial resolutions evaluated, 0.5 m consistently produced similar or better models than those using the 0.25 or 1.0 m resolutions. These results demonstrate the potential for integration of spatial transforms of image fractions and raw brightness in high-resolution modelling of forest structure and health.     

Needle chlorophyll content estimation through model inversion using hyperspectral data from boreal conifer forest canopies

Leaf chlorophyll content in coniferous forest canopies, a measure of stand condition, is the target of studies and models linking leaf reflectance and transmittance and canopy hyperspectral reflectance imagery. The viability of estimation of needle chlorophyll content from airborne hyperspectral optical data through inversion of linked leaf level and canopy level radiative transfer models is discussed in this paper. This study is focused on five sites of Jack Pine (Pinus banksiana Lamb.) in the Algoma Region (Canada), where field, laboratory and airborne data were collected in 1998 and 1999 campaigns. Airborne hyperspectral CASI data of 72 bands in the visible and near-infrared region and 2 m spatial resolution were collected from 20×20 m study sites of Jack Pine in 2 consecutive years. It was found that needle chlorophyll content could be estimated at the leaf level (r²=0.4) by inversion of the PROSPECT leaf model from needle reflectance and transmittance spectra collected with a special needle carrier apparatus coupled to the Li-Cor 1800 integrating sphere. The Jack Pine forest stands used for this study with LAI>2, and the high spatial resolution hyperspectral reflectance collected, allowed the use of the SPRINT canopy reflectance model coupled to PROSPECT for needle chlorophyll content estimation by model inversion. The optical index R750/R710 was used as the merit function in the numerical inversion to minimize the effect of shadows and LAI variation in the mean canopy reflectance from the 20×20 m plots. Estimates of needle pigment content from airborne hyperspectral reflectance using this linked leaf-canopy model inversion methodology showed an r²=0.4 and RMSE=8.1 μg/cm² when targeting sunlit crown pixels in Jack Pine sites with pigment content ranging between 26.8 and 56.8 μg/cm² (1570-3320 μg/g).