Impact of understory vegetation on forest canopy reflectance and remotely sensed LAI estimates

Forest leaf area index (LAI), is an important variable in carbon balance models. However, understory vegetation is a recognized problem that limits the accuracy of satellite-estimated forest LAI. A canopy reflectance model was used to investigate the impact of the understory vegetation on LAI estimated from reflectance values estimated from satellite sensor data. Reflectance spectra were produced by the model using detailed field data as input, i.e. forest LAI, tree structural parameters, and the composition, distribution and reflectance of the forest floor. Common deciduous and coniferous forest types in southern Sweden were investigated. A negative linear relationship (r² = 0.6) was observed between field estimated LAI and the degree of understory vegetation, and the results indicated better agreement when coniferous and deciduous stands were analysed separately. The simulated spectra verified that the impact of the understory on the reflected signal from the top of the canopy is important; the reflectance values varying by up to ± 18% in the red and up to ± 10% in the near infra-red region of the spectra due to the understory. In order to predict the variation in LAI due to the understory vegetation, model inversions were performed where the input spectra were changed between the minimum, average and maximum reflectance values obtained from the forward runs. The resulting variation in LAI was found to be 1.6 units on average. The LAI of the understory could be predicted indirectly from simple stand data on forest characteristics, i.e. from allometric estimates, as an initial step in the process of estimating LAI. It is suggested here that compensation for the effect of the understory would improve the accuracy in the estimates of canopy LAI considerably.     

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.     

Application of photon recollision probability in coniferous canopy reflectance simulations

A new semi-physical forest reflectance model, PARAS, is presented in the paper. PARAS is a simple parameterization model for taking into account the effect of within-shoot scattering on coniferous canopy reflectance. Multiple scattering at the small scale represented by a shoot is a conifer-specific characteristic which causes the spectral signature of coniferous forests to differ from that of broadleaved forests. This has for long led to problems in remote sensing of canopy structural variables in coniferous dominated regions. The PARAS model uses a relationship between photon recollision probability and leaf area index (LAI) for simulating forest reflectance. The recollision probability is a measurable, wavelength independent variable which is defined as the probability with which a photon scattered in the canopy interacts with a phytoelement again. In this study, we present application results using PARAS in simulating reflectance of coniferous forests for approximately 800 Scots pine and Norway spruce dominated stands. The results of this study clearly indicate that a major improvement in simulating canopy reflectance in near-infrared (NIR) is achieved by simply accounting for the within-shoot scattering. In other words, the low NIR reflectance observed in coniferous areas is mainly due to within-shoot scattering. In the red wavelength the effect of within-shoot scattering was not pronounced due to the high level of needle absorption in the red range. To conclude the paper, further application possibilities of the presented parameterization model are discussed.     

Evaluate Quantitatively Effects of Understory on Leaf Area Index (LAI) Estimation Combining Active and Passive Remote Sensing Data

Natural forests have the vertical three\|dimensional structure of canopy and understory vegetation (shrubs,grasslands,and bare soil).Accurate and quantitative separation of understory vegetation is of great scientific significance and practicality on improving the precision of inversion of forest canopy leaf area index.value.Due to the limitations of traditional passive optical remote sensing data on directly acquiring 3D information,this study intends to combine active and passive ALS and HyperMap data with the Washington Botanic Garden as the key research area.On the basis of individual tree segmentation,the vertical stratification of the forest (forest canopy and undergrowth vegetation layer) is achieved.On this basis,the forest canopy laser point cloud data was used to remove the understory information from the optical image data.By comparing the results of the forest effective leaf area index obtained from aerial optical images and ground measurements,it was found that：(1) forest canopy density has a significant impact on the penetration of ALS data;(2) removal of understory information can effectively improve the forest crown accuracy of LAIe estimated.The correlation between Normalized Difference Vegetation Index (NDVI) and ground surface measured effective leaf area index increased from 0.087 to 0.591.In addition,the optical remote sensing image based on the removal of understory vegetation information was compared with the Simple Ratio vegetation index (SR) (the correlation increased from 0.209 to 0.559) and the simplified simple Ratio vegetation index (RSR) (the correlation increased from 0.147 to 0.358).The NDVI was most sensitive to changes in canopy leaf area index (correlation increased by 0.5).The method of quantitatively separating understory vegetation with the combined active and passive remote sensing data proposed in this study can effectively improve the accuracy of inversion of forest canopy leaf area index,and provide a solid foundation for quantitative and accurate estimate of forest biophysical parameters and study of carbon and water cycle processes.     

Effects of woody elements on simulated canopy reflectance: Implications for forest chlorophyll content retrieval

An important bio-indicator of actual plant health status, the foliar content of chlorophyll a and b (Cab), can be estimated using imaging spectroscopy. For forest canopies, however, the relationship between the spectral response and leaf chemistry is confounded by factors such as background (e.g. understory), canopy structure, and the presence of non-photosynthetic vegetation (NPV, e.g. woody elements)—particularly the appreciable amounts of standing and fallen dead wood found in older forests. We present a sensitivity analysis for the estimation of chlorophyll content in woody coniferous canopies using radiative transfer modeling, and use the modeled top-of-canopy reflectance data to analyze the contribution of woody elements, leaf area index (LAI), and crown cover (CC) to the retrieval of foliar Cab content. The radiative transfer model used comprises two linked submodels: one at leaf level (PROSPECT) and one at canopy level (FLIGHT). This generated bidirectional reflectance data according to the band settings of the Compact High Resolution Imaging Spectrometer (CHRIS) from which chlorophyll indices were calculated. Most of the chlorophyll indices outperformed single wavelengths in predicting Cab content at canopy level, with best results obtained by the Maccioni index ([R₇₈₀ − R₇₁₀] / [R₇₈₀ − R₆₈₀]). We demonstrate the performance of this index with respect to structural information on three distinct coniferous forest types (young, early mature and old-growth stands). The modeling results suggest that the spectral variation due to variation in canopy chlorophyll content is best captured for stands with medium dense canopies. However, the strength of the up-scaled Cab signal weakens with increasing crown NPV scattering elements, especially when crown cover exceeds 30%. LAI exerts the least perturbations. We conclude that the spectral influence of woody elements is an important variable that should be considered in radiative transfer approaches when retrieving foliar pigment estimates in heterogeneous stands, particularly if the stands are partly defoliated or long-lived.     

Seasonal reflectance dynamics of common understory types in a northern European boreal forest

The influence of the seasonal cycle of boreal forest understory has been noticed in global remote sensing of vegetation, especially in remote sensing of biophysical properties (e.g. leaf area index) of the tree-layer in a forest. A general problem in the validation of operationally produced global biophysical vegetation products is the lack of ground reference data on the seasonal variability of different land surface types. Currently, little is known about the spectral properties of the understory layers of boreal forests, and even less is known about the seasonal dynamics of the spectra. In this paper, we report seasonal trajectories of understory reflectance spectra measured in a European boreal forest. Four study sites representing different forest fertility site types were selected from central Finland. The understory composition was recorded and its spectra measured with an ASD FieldSpec Hand-Held UV/VNIR Spectroradiometer ten times during the growing period (from May to September) in 2010. Our results show that the spectral differences between and within understory types are the largest at the peak of the growing season in early July whereas in the beginning and end of the growing season (i.e. early May and late September, respectively) the differences between the understory types are marginal. In general, the fertile sites had the brightest NIR spectra throughout the growing season whereas infertile types appeared darker in NIR. Our results also indicated that a mismatch in the seasonal development of understory and tree layers does not occur in boreal forests: the understory and tree layer vegetation develop at a similar pace in the spring (i.e. there are no or only few spring ephemerals present), and the forests with the strongest seasonal dynamics in tree canopy structure (LAI) have also the strongest dynamics in understory spectra.     

Estimating biomass for boreal forests using ASTER satellite data combined with standwise forest inventory data

In the present study, the suitability of optical ASTER satellite data (with 9 spectral bands) for estimating the biomass of boreal forest stands in mineral soils was tested. The remote sensing data were analysed and tested together with standwise forest inventory data. Stand volume estimates were converted to aboveground tree biomass using biomass expansion factors, and the aboveground biomass of understory vegetation was predicted according to the stand age. Non-linear regression analysis and neural networks were applied to develop models for predicting biomass according to standwise ASTER reflectance. All ASTER bands appeared to be sensitive to tree biomass, in particular the green band 1. The relative estimation errors (RMSE_r) of the total aboveground biomass of the forest stands were 44.7% and 41.0% using multiple regression analysis and neural networks, respectively. Although the estimation errors remained large, the predictions were relatively accurate in comparison to previous studies. Furthermore, the predictions obtained here were significantly close to the municipality-level mean values provided by the National Forest Inventory of Finland.     

The effect of crown shape on the reflectance of coniferous stands

The Kuusk-Nilson forest reflectance model was used to study the effect of crown shape on the reflectance of Scots pine and Norway spruce stands. In the first part of the study, we examined spruce and pine stands with an age range of 20-100 years and compared their simulated hemispherical-directional reflectance factors (HDRFs) at nadir in red (661 nm), NIR (838 nm) and MIR (1677 nm) when crowns were modeled as ellipsoids or cones. In all the cases, when a stand was modeled with conical crowns, it had a smaller reflectance factor than the same stand with ellipsoidal crowns.To analyze the sensitivity of HDRF on crown shape, in the second part of the study we simulated the angular distributions of HDRF of two pine stands with different leaf area index (LAI) and canopy closure values at 661 nm assuming four different crown shapes (cone, cylinder, ellipsoid, and cylinder bottom, cone top) and separated the components forming the HDRF. Considerable difference in the HDRF between the four crown shapes was observed: The larger the crown volume, the higher the canopy reflectance at similar LAI and canopy closure. A comparison of the two stands revealed that in denser stands (with a higher canopy closure) single scattering from tree crowns was responsible for the difference in HDRF between the different crown shapes, whereas in stands with a smaller canopy closure the single scattering from ground dominated the HDRF. Finally, the role of crown shape for the retrieval of LAI by inversion from remotely sensed data is discussed.     

3D modelling of forest canopy structure for remote sensing simulations in the optical and microwave domains

A detailed 3D structural model of a conifer forest canopy was developed in order to simulate the reflectance (optical) and backscatter (microwave) signals measured remotely. We show it is feasible to model forest canopy scattering using detailed 3D models of tree structure including the location and orientation of individual needles. An existing structural growth model of Scots pine (Pinus sylvestris L.), Treegrow, was modified to simulate observed growth stages of a Scots pine canopy from age 5 to 50 years. The 3D tree models showed close structural agreement with in situ measurements. Needles were added to the structural models according to observed phyllotaxy (distribution). Individual trees were used to generate model canopies, which in turn were used to drive optical and microwave models of canopy scattering. Simulated canopy radiometric response was compared with airborne hyperspectral reflectance data (HyMAP) and airborne synthetic aperture RADAR (ASAR) backscatter data. Model simulations agreed well in general with observations, particularly at optical wavelengths where model simulations of low and high density canopy stands were shown to bracket observations. Relatively small sensitivity of observed reflectance to canopy age was captured reasonably well by the simulations. The choice of needle shape and phyllotaxy was shown to have a significant impact on multiple scattering behaviour at the branch scale. In the microwave domain, simulated backscatter values agreed reasonably well with observations at L-band, less so at X-band. L-band simulated backscatter significantly underestimated observed backscatter at younger canopy ages, probably as a result of inappropriate modelling of soil/understory. It is demonstrated that a combined structural and radiometric modelling approach provides a flexible and powerful method for simulating the remotely sensed signal of a forest canopy in the optical and microwave domains. This is particularly useful for exploring the impact of canopy structure on the resulting signal and also for combined retrievals of forest structural parameters from optical and microwave data.     

Mapping and monitoring conifer mortality using remote sensing in the Lake Tahoe Basin

A prolonged drought in the western United States has resulted in alarming levels of mortality in conifer forests. Satellite remote sensing holds the potential for mapping and monitoring the effects of such environmental changes over large geographic areas in a timely manner. Results from the application of a forest canopy reflectance model using multitemporal Landsat TM imagery and field measurements, indicate conifer mortality can be effectively mapped and inventoried. The test area for this project is the Lake Tahoe Basin Management Unit in the Sierra Nevada of California. The Landsat TM images are from the summers of 1988 and 1991. The Li-Strahler canopy model estimates several forest stand parameters, including tree size and canopy cover for each conifer stand, from reflectance values in satellite imagery. The difference in cover estimates between the dates forms the basis for stratifying stands into mortality classes, which are used as both themes in a map and the basis of the field sampling design. Field measurements from 61 stands collected in the summer of 1992 indicate 15 % of the original timber volume in the true fir zone died between 1988 and 1992. The resulting low standard error of 11 % for this estimate indicates the utility of these mortality classes for detecting areas of high mortality. Also, the patterns in the estimated mean timber volume loss for each class follow the expected trends. The results of this project are immediately useful for fire hazard management, by providing both estimates of the degree of overall mortality and maps showing its location. They also indicate current remote sensing technology may be useful for monitoring the changes in vegetation that are expected to result from climate change.     

Satellite estimation of tropical secondary forest above-ground biomass: Data from Brazil and Bolivia

This paper reports on a test of the ability to estimate above-ground biomass of tropical secondary forest from canopy spectral reflectance using satellite optical data. Landsat Thematic Mapper data were acquired concurrent with field surveys conducted in secondary forest fallows near Manaus, Brazil and Santa Cruz de la Sierra, Bolivia. Measurements of age and above-ground live biomass were made in 34 regrowth stands. Satellite data were converted to surface reflectances and compared with regrowth stand age, biomass and structural variables. Among the Brazilian stands, significant relationships were observed between middle-infrared reflectance and stand age, height, volume and biomass. The canopy reflectance-biomass relationship saturated at around 15.0 kg m^-2, or over 15 years of age (r > 0.80, p < 0.01). In the Bolivian study area, no significant relationship between canopy spectral reflectance and biomass was observed. These contrasting results are probably caused by a low Sun angle during the satellite measurements from Bolivia. However, regrowth structural and general compositional differences between the two study areas could explain the lack of a significant relationship in Bolivia. The results demonstrate a current potential for biomass estimation of secondary forests with satellite optical data in some, but not all, tropical regions. A discussion of the potential for regional extrapolation of spectral relationships and future satellite imagery is included.     

Remote sensing of upland conifer plantations using Landsat TM data: A case study from Galloway,south-west Scotland

The relationship between reflectance and forest structure is investigated for 52 forest stands in Galloway. Height and basal area are most strongly correlated with reflectance in TM bands 3, 5 and 7 (R² >= 0.77, P < 0.01). The strength of these relationships decreases substantially after canopy closure has occurred. Stem density cannot be predicted from canopy reflectance. Tree age can only be reliably predicted for areas where growth is uniform, which is rare in upland environments. Analysis of TM imagery provides a cost effective method of mapping certain canopy parameters provided that limited detail is required for older mature plantations with a closed canopy.     

Leaf chlorophyll content retrieval from airborne hyperspectral remote sensing imagery

Hyperspectral remote sensing has great potential for accurate retrieval of forest biochemical parameters. In this paper, a hyperspectral remote sensing algorithm is developed to retrieve total leaf chlorophyll content for both open spruce and closed forests, and tested for open forest canopies. Ten black spruce (Picea mariana (Mill.)) stands near Sudbury, Ontario, Canada, were selected as study sites, where extensive field and laboratory measurements were carried out to collect forest structural parameters, needle and forest background optical properties, and needle biophysical parameters and biochemical contents chlorophyll a and b. Airborne hyperspectral remote sensing imagery was acquired, within one week of ground measurements, by the Compact Airborne Spectrographic Imager (CASI) in a hyperspectral mode, with 72 bands and half bandwidth 4.25-4.36 nm in the visible and near-infrared region and a 2 m spatial resolution. The geometrical-optical model 4-Scale and the modified leaf optical model PROSPECT were combined to estimate leaf chlorophyll content from the CASI imagery. Forest canopy reflectance was first estimated with the measured leaf reflectance and transmittance spectra, forest background reflectance, CASI acquisition parameters, and a set of stand parameters as inputs to 4-Scale. The estimated canopy reflectance agrees well with the CASI measured reflectance in the chlorophyll absorption sensitive regions, with discrepancies of 0.06%-1.07% and 0.36%-1.63%, respectively, in the average reflectances of the red and red-edge region. A look-up-table approach was developed to provide the probabilities of viewing the sunlit foliage and background, and to determine a spectral multiple scattering factor as functions of leaf area index, view zenith angle, and solar zenith angle. With the look-up tables, the 4-Scale model was inverted to estimate leaf reflectance spectra from hyperspectral remote sensing imagery. Good agreements were obtained between the inverted and measured leaf reflectance spectra across the visible and near-infrared region, with R² = 0.89 to R² = 0.97 and discrepancies of 0.02%-3.63% and 0.24%-7.88% in the average red and red-edge reflectances, respectively. Leaf chlorophyll content was estimated from the retrieved leaf reflectance spectra using the modified PROSPECT inversion model, with R² = 0.47, RMSE = 4.34 μg/cm², and jackknifed RMSE of 5.69 μg/cm² for needle chlorophyll content ranging from 24.9 μg/cm² to 37.6 μg/cm². The estimates were also assessed at leaf and canopy scales using chlorophyll spectral indices TCARI/OSAVI and MTCI. An empirical relationship of simple ratio derived from the CASI imagery to the ground-measured leaf area index was developed (R² = 0.88) to map leaf area index. Canopy chlorophyll content per unit ground surface area was then estimated, based on the spatial distributions of leaf chlorophyll content per unit leaf area and the leaf area index.     

Quantifying empirical relations between planted species mixtures and canopy reflectance with PROTEST

Mapping plant species composition of mixed vegetation stands with remote sensing is a complicated task. Uncertainties may arise from similar spectral signatures of different plant species as well as from variable influences of prevailing plant states (e.g., growth stages, vigor, or stress levels). Despite these uncertainties, empirical approaches may often be able to take up the challenge. However, their performance is likely to be affected by the temporal variability of empirical relations between reflectance and plant species composition. To assess some aspects of this temporal variability, we performed a greenhouse study. Three mixed stands of grassland species were planted with defined spatial variation in species proportions. The canopy reflectance of these mixed stands was measured with a field spectrometer over a period of three months. Confounding external influences on plant states apart from maturation were minimized.The suitability of canopy reflectance and derivative reflectance to draw conclusions on differences in qualitative species mixtures between the stands was tested with a classification approach (Spectral Angle Mapper, SAM). Procrustean randomization test (PROTEST), which is to our knowledge new to the field of remote sensing, was applied in combination with Isometric Feature Mapping to quantify the spectral variation caused by within-stand spatial variation in species proportions. Model fits in both analyses increased with progressing plant development; further, utilization of derivative reflectance improved the model fits. Regardless of the within-stand variation, SAM enabled a successful discrimination of the three stands with an average overall accuracy of 85% (reflectance) and 92% (derivative reflectance). In PROTEST analysis, spatial variation in reflectance was successfully related to within-stand variation in species proportions. However, observed influences of variable growth stages and health states on these relations were considerable. The temporal variation of these relations (r = 0.27-0.73 for reflectance and 0.48-0.73 for derivative reflectance) was quantified for the first time under controlled conditions.     

Simple reflectance indices track heat and water stress-induced changes in steady-state chlorophyll fluorescence at the canopy scale

Non-invasive remote sensing techniques for monitoring plant stress and photosynthetic status have received much attention. The majority of published vegetation indices are not sensitive to rapid changes in plant photosynthetic status brought on by common environmental stressors such as diurnal fluxes in irradiance and heat. This is due to the fact that most vegetation indices have no direct link to photosynthetic functioning beyond their sensitivity to canopy structure and pigment concentration changes. In contrast, this study makes progress on a more direct link between passive reflectance measurements and plant physiological status through an understanding of photochemical quenching (qP) and non-photochemical quenching processes. This is accomplished through the characterization of steady-state fluorescence (Fs) and its influence on apparent reflectance in the red-edge spectral region. A series of experiments were conducted under controlled environmental conditions, linking passive reflectance measurements of a grapevine canopy (Vitis vinifera L. cv. Cabernet Sauvignon) to leaf level estimates of CO₂ assimilation (A), stomatal conductance (g), qP, and Fs. Plant stress was induced by imposing a diurnal heat stress and recovery event and by withholding water from the plant canopy over the course of the experiment. We outlined evidence for a link between Fs and photosynthetic status, identified the Fs signal in passive remote sensing reflectance data, and related reflectance-derived estimates of Fs to plant photosynthetic status. These results provide evidence that simple reflectance indices calculated in the red-edge spectral region can track temperature and water-induced changes in Fs and, consequently, provide a rapid assessment of plant stress that is directly linked to plant physiological processes.     

A multi-scale approach to mapping effective Leaf Area Index in Boreal Picea mariana stands using high spatial resolution CASI imagery

This study evaluates the performance of a simple geometric-optics reflectance model, used in combination with multi-spectral clustering, to map spatial patterns of effective Leaf Area Index ( L e ) within Boreal Picea mariana stands. Two metre resolution Compact Airborne Spectrographic Imager (CASI) images, acquired during the winter to minimize variability in understory reflectance, are used to map L e over BOREAS northern and southern Old Black Spruce tower flux sites. A combined multi-spectral clustering and ray-tracing approach is used to map open areas in each site at 2 m scale. A modified version of the Forest Light Interaction Model (FLIM) is then applied over canopy areas using 30 m scale red and near-infrared reflectance values derived from CASI images. Comparison of the combined FLIM and clustering approach (FLIMCLUS) with surface L e measurements in areas with overstory cover indicate a R 2 of 0.67 for the SSA-OBS site and 0.16 for the NSA-OBS site. The poor NSAOBS performance may be due to the low observed range of L e along the transect selected since additional measurements in a sparse canopy area closely match FLIM estimates. The relative standard error at both sites is under 10% and is close to the 5% precision error in surface L e estimates. Comparison of FLIMCLUS with FLIM and the simple ratio indicate substantial differences over open areas where the latter methods map-zero L e values. Further validation over other study sites, including surface data mapping edges between canopy and open areas, is proposed. The FLIM-CLUS L e maps may be useful for testing scale dependent assumptions within remote sensing algorithms and ecosystem flux models applied to the study sites and similar Picea mariana stands.     

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).     

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.     

Forest canopy recovery from the 1938 hurricane and subsequent salvage damage measured with airborne LiDAR

The structure of a forest canopy often reflects its disturbance history. Such signatures of past disturbances or legacies can influence how the ecosystem functions across broad spatio-temporal scales. The 1938 hurricane and ensuing salvage operations which swept through New England represent the most recent large, infrequent disturbance (LID) in this region. Though devastating (downing ∼ 70% of the timber at Harvard Forest), the disturbance was not indiscriminate; it left behind a heterogeneous landscape comprised of different levels of canopy damage. We analyzed large-footprint LiDAR, from the Prospect Hill tract at Harvard Forest in central Massachusetts, to assess whether damage to the forest structure from the hurricane and subsequent timber extraction could be discerned after ∼ 65 years. Differences in LiDAR-derived measures of canopy height and vertical diversity were a function of the degree of damage from the 1938 hurricane and the predominant tree species which is, in part, a function of land use history. Higher levels of damage corresponded to slightly shorter canopies with a less even vertical distribution of return from the ground to the top. In addition, differences in canopy topography as revealed by spatial autocorrelation of canopy top heights were found among the damage classes. Less disturbed stands were characterized by lower levels of local autocorrelation for canopy height and higher levels of vertical diversity of LiDAR returns. These differences in canopy structure reveal that the forest tract has not completely recovered from the 1938 LID and salvage regime, which may have implications on arboreal and understory habitat and other ecosystem functions.     

Direct retrieval of the shape of leaf spectral albedo from multiangular hyperspectral Earth observation data

Physically-based retrieval of vegetation canopy properties from remote sensing data presumes a knowledge of the spectral albedo of the basic scattering unit, i.e. leaf. In this paper, we present a novel method to directly retrieve the spectral dependence of leaf single-scattering albedo of a closed broadleaf forest canopy from multiangular hyperspectral satellite imagery. The new algorithm is based on separating the reflected signal into a linear (first-order) and non-linear (diffuse) reflectance component. A limitation of the proposed algorithm is that the leaf single-scattering albedo ω(λ) is retrieved with an accuracy of a structural parameter (called a₀) which, in turn, depends on canopy bidirectional gap probability, ratio of leaf reflectance to transmittance, and distribution of leaf normals. The structural parameter (a₀) was found to depend on tree-level structural parameters, such as tree height and volume of a single crown, but not the amount of leaf area.