A simple parameterization of canopy reflectance using photon recollision probability

Information on the fractions of incident radiation reflected, transmitted and absorbed by a plant canopy is crucial in remote sensing of vegetation and modeling of canopy microclimate. Photon recollision probability p allows to calculate easily the spectral behavior of canopy scattering, i.e. the sum of canopy reflectance and transmittance. However, to divide the scattered radiation into reflected and transmitted fluxes, additional models are needed. In this paper, we present a simple formula to estimate the fraction of radiation scattered upwards by a canopy. The new method is semi-empirical, makes use of the concept of photon recollision probability, and is derived from an analysis of modeling results. Although a physical interpretation is given for the single additional parameter needed in the formula, the scattering asymmetry parameter q, the method is not strictly based on the radiative transfer equation. Our results indicate that the method is accurate for low to moderate leaf area index (LAI) values, and provides a reasonable approximation even at LAI = 8. In addition, we present a method to compute p using numerical radiative transfer models.     

Simple analytical formula for calculating average photon recollision probability in vegetation canopies

The concept of recollision probability originates from the theory of canopy spectral invariants (‘p-theory’) but is a simplification that involves several heuristic assumptions. Nonetheless, the concept has been shown to be a useful tool for incorporating the effects of 3D structure on canopy absorptive and reflective properties in forest reflectance models. A method is presented by which an average value of the canopy recollision probability (pˆ) can be calculated from canopy gap fraction data, which are provided for example by the LAI-2000 plant canopy analyzer or can be extracted from fisheye photographs. The new method was used to calculate the average recollision probabilities (pˆ values) of uniform leaf and shoot canopies, showing good agreement with results from Monte Carlo simulations. Strengths of the method presented here are the explicitly formulated relationship between recollision probability and canopy structure, and its direct applicability in canopy RT studies.     

Photon recollision probability in discrete crown canopies

The photon recollision probability in vegetation canopies, defined as the probability that a photon, after having interacted with a canopy element, will interact again, is a useful tool in remote sensing and ecological applications, enabling to link canopy optical properties at different wavelength and to estimate radiation absorption. In this work, a method is presented to estimate the photon recollision probability for horizontally homogeneous leaf canopies with arbitrary leaf angle distribution as well as for discrete crown canopies. The estimation is based on analytical approximation of the first-order recollision probability. Using the analytical solution of the two-stream equations of radiative transfer and Monte Carlo modeling, the first-order photon recollision probability is shown to slightly underestimate the mean recollision probability. Also, an approximation formula for the mean recollision probability in a horizontally homogeneous canopy is presented as a function of leaf area index. The method to calculate photon recollision probability in discrete crown canopies requires only the knowledge of total and between-crown canopy transmittance and is thus independent of the geometric-optical model used.     

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.     

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.     

Utility of an image-based canopy reflectance modeling tool for remote estimation of LAI and leaf chlorophyll content at the field scale

This paper presents a physically-based approach for estimating critical variables describing land surface vegetation canopies, relying on remotely sensed data that can be acquired from operational satellite sensors. The REGularized canopy reFLECtance (REGFLEC) modeling tool couples leaf optics (PROSPECT), canopy reflectance (ACRM), and atmospheric radiative transfer (6SV1) model components, facilitating the direct use of at-sensor radiances in green, red and near-infrared wavelengths for the inverse retrieval of leaf chlorophyll content (C_ab) and total one-sided leaf area per unit ground area (LAI). The inversion of the canopy reflectance model is constrained by assuming limited variability of leaf structure, vegetation clumping, and leaf inclination angle within a given crop field and by exploiting the added radiometric information content of pixels belonging to the same field. A look-up-table with a suite of pre-computed spectral reflectance relationships, each a function of canopy characteristics, soil background effects and external conditions, is accessed for fast pixel-wise biophysical parameter retrievals. Using 1 m resolution aircraft and 10 m resolution SPOT-5 imagery, REGFLEC effectuated robust biophysical parameter retrievals for a corn field characterized by a wide range in leaf chlorophyll levels and intermixed green and senescent leaf material. Validation against in-situ observations yielded relative root-mean-square deviations (RMSD) on the order of 10% for the 1 m resolution LAI (RMSD = 0.25) and C_ab (RMSD = 4.4 μg cm^− 2) estimates, due in part to an efficient correction for background influences. LAI and C_ab retrieval accuracies at the SPOT 10 m resolution were characterized by relative RMSDs of 13% (0.3) and 17% (7.1 μg cm^− 2), respectively, and the overall intra-field pattern in LAI and C_ab was well established at this resolution. The developed method has utility in agricultural fields characterized by widely varying distributions of model variables and holds promise as a valuable operational tool for precision crop management. Work is currently in progress to extend REGFLEC to regional scales.     

Coupling forest canopy and understory reflectance in the Arctic latitudes of Finland

The Arctic region is predicted to experience considerable climatic and environmental changes as the global atmospheric CO₂ increases. Growing awareness of the role of tundra and taiga ecosystems and their transition zone in the climate change process has resulted in a recent increase in remote sensing studies focusing on the Arctic latitudes. Remote sensing of biophysical properties of the canopy layer in the forested part of the region is often, however, challenged by the dominating role of the understory in the spectral signal. In this paper, we examine the influence of understory vegetation on forest reflectance in the Arctic region of Finland during no-snow conditions. The study is based on SPOT HRVIR images, field goniospectrometry, 300 ground reference plots and a physically-based forest reflectance model (PARAS). The results indicate that lichen-dominated forest site types can be distinguished from sites dominated by dwarf shrubs. The paper also contains results from applying an analytical method for calculating photon recollision probability from canopy transmittance data for forest stands, and then using it to simulate the reflectance of the same stands.     

NDVI responses to the forest canopy and floor from spring to summer observed by airborne spectrometer in eastern Siberia

A popular method of satellite-based monitoring of the photosynthetic potential of vegetation is to calculate the normalised difference vegetation index (NDVI) from measurements of the red (RED) and near-infrared (NIR) bands. Enormous amounts of vegetation information have been obtained over continental to global areas based on NDVI derived from NOAA-AVHRR, Terra/Aqua-MODIS, and SPOT-VEGETATION satellite observations. In eastern Siberia, where sparse boreal forests are dominant, the lack of landscape-scale canopy-reflectance observations impedes interpretation of how NDVI seasonality is controlled by the forest canopy and floor status. We discuss the NDVI of the canopy and floor separately based on airborne spectral reflectance measurements and simultaneous airborne land surface images acquired around Yakutsk, Siberia, using a hedgehopping aircraft from spring to summer 2000. The aerial land surface images (4402 scenes) were visually classified into four types according to the forest condition: no-green canopy and snow floor (Type 1), green canopy and snow floor (Type 2), no-green canopy and no-snow floor (Type 3), and green canopy and no-snow floor (Type 4). The spectral reflectance from 350 to 1200 nm was then calculated for these four types. Type 1 had almost no difference in reflectance between the RED and NIR bands, and the resultant NDVI was slightly negative (− 0.03). Although Type 2 showed a significant difference between the two bands because of canopy greenness, the resultant NDVI was rather low (0.17) because of high reflection from the snow cover on the floor. In Type 3, the significant difference between the two bands was mainly caused by the greenness of the floor, and the NDVI was relatively high (0.45). The NDVI for Type 4 was the highest (0.75) among the four types. The contributions of reflectance from the forest canopy and floor to the total reflectance were tested with a forest radiative transfer model. The reflectance difference between NIR and RED bands (NIR − RED) of Type 4 (15.6%) was approximately double the differences of Type 2 (7.0%) and of Type 3 (7.9%), suggesting half-and-half contributions of forest canopy greenness and floor greenness to the total greenness. The result also suggested that the satellite-derived NDVI in the larch forest around Yakutsk reaches 85% of the maximum NDVI owing to the forest floor greenness, and only the other 15% of the increase in NDVI should be attributed to the canopy foliation. These results quantitatively reveal that the NDVI depends considerably on forest floor greenness and snow cover in addition to canopy greenness in the case of relatively sparse forest in Siberia.     

Assessment of the MODIS LAI product for Australian ecosystems

The leaf area index (LAI) product from the Moderate Resolution Imaging Spectroradiometer (MODIS) is important for monitoring and modelling global change and terrestrial dynamics at many scales. The algorithm relies on spectral reflectances and a six biome land cover classification. Evaluation of the specific behaviour and performance of the product for regions of the globe such as Australia are needed to assist with product refinement and validation. We made an assessment of Collection 4 of the MODIS LAI product using four approaches: (a) assessment against a continental scale Structural Classification of Australian Vegetation (SCAV); (b) assessment against a continental scale land use classification (LUC); (c) assessment against historical field-based measurement of LAI collected prior to the Terra Mission; and (d) direct comparison of MODIS LAI with coincident field measurements of LAI, mostly from hemispherical photography. The MODIS LAI product produced a wide variety of geographically and structurally specific temporal response profiles between different classes and even for sub-groups within classes of the SCAV. Historical and concurrent field measurements indicated that MODIS LAI was giving reasonable estimates for LAI for most cover types and land use types, but that major overestimation of LAI occurs in some eastern Australian open forests and woodlands. The six biome structural land cover classification showed some significant deviations in class allocation compared to the SCAV particularly where grasslands are allocated to shrubland, savanna woodlands are allocated to shrubland, savanna and broadleaf forest, and open forests are allocated to savanna and broadleaf forest. The land cover and LAI products could benefit from some additional examination of Australian data addressing the structural representation of Eucalypt canopies in the “space of canopy realisation” for savanna and broadleaf forest classes.     

Mapping leaf chlorophyll and leaf area index using inverse and forward canopy reflectance modeling and SPOT reflectance data

Reflectance data in the green, red and near-infrared wavelength region were acquired by the SPOT high resolution visible and geometric imaging instruments for an agricultural area in Denmark (56°N, 9°E) for the purpose of estimating leaf chlorophyll content (C_ab) and green leaf area index (LAI). SPOT reflectance observations were atmospherically corrected using aerosol data from MODIS and profiles of air temperature, humidity and ozone from the Atmospheric Infrared Sounder (AIRS), and used as input for the inversion of a canopy reflectance model. Computationally efficient inversion schemes were developed for the retrieval of soil and land cover-specific parameters which were used to build multiple species and site dependent formulations relating the two biophysical properties of interest to vegetation indices or single spectral band reflectances. Subsequently, the family of model generated relationships, each a function of soil background and canopy characteristics, was employed for a fast pixel-wise mapping of C_ab and LAI.The biophysical parameter retrieval scheme is completely automated and image-based and solves for the soil background reflectance signal, leaf mesophyll structure, specific dry matter content, Markov clumping characteristics, C_ab and LAI without utilizing calibration measurements.Despite the high vulnerability of near-infrared reflectances (ρ_nir) to variations in background properties, an efficient correction for background influences and a strong sensitivity of ρ_nir to LAI, caused LAI-ρ_nir relationships to be very useful and preferable over LAI-NDVI relationships for LAI prediction when LAI > 2. Reflectances in the green waveband (ρ_green) were chosen for producing maps of C_ab.The application of LAI-NDVI, LAI-ρ_nir and C_ab-ρ_green relationships provided reliable quantitative estimates of C_ab and LAI for agricultural crops characterized by contrasting architectures and leaf biochemical constituents with overall root mean square deviations between estimates and in-situ measurements of 0.74 for LAI and 5.0 μg cm^− 2 for C_ab.The results of this study illustrate the non-uniqueness of spectral reflectance relationships and the potential of physically-based inverse and forward canopy reflectance modeling techniques for a reasonably fast and accurate retrieval of key biophysical parameters at regional scales.     

Estimating aboveground forest biomass from canopy reflectance model inversion in mountainous terrain

The amount and spatial distribution of aboveground forest biomass (AGB) are required inputs to forest carbon budgets and ecosystem productivity models. Satellite remote sensing offers distinct advantages for large area and multi-temporal applications, however, conventional empirical methods for estimating forest canopy structure and AGB can be difficult in areas of high relief and variable terrain. This paper introduces a new method for obtaining AGB from forest structure estimates using a physically-based canopy reflectance (CR) model inversion approach. A geometric-optical CR model was run in multiple forward mode (MFM) using SPOT-5 imagery to derive forest structure and biomass at Kananaskis, Alberta in the Canadian Rocky Mountains. The approach first estimates tree crown dimensions and stem density for satellite image pixels which are then related to tree biomass and AGB using a crown spheroid surface area approach. MFM estimates of AGB were evaluated for 36 deciduous (trembling aspen) and conifer (lodgepole pine) field validation sites and compared against spectral mixture analysis (SMA) and normalised difference vegetation index (NDVI) biomass predictions from atmospherically and topographically corrected (SCS+C) imagery. MFM provided the lowest error for all validation plots of 31.7 tonnes/hectare (t/ha) versus SMA (32.6 t/ha error) and NDVI (34.7 t/ha) as well as for conifer plots (MFM: 23.0 t/ha; SMA 27.9 t/ha; NDVI 29.7 t/ha) but had higher error than SMA and NDVI for deciduous plots (by 4.5 t/ha and 2.1 t/ha, respectively). The MFM approach was considerably more stable over the full range of biomass values (67 to 243 t/ha) measured in the field. Field plots with biomass > 1 standard deviation from the field mean (over 30% of plots) had biomass estimation errors of 37.9 t/ha using MFM compared with 65.5 t/ha and 67.5 t/ha error from SMA and NDVI, respectively. In addition to providing more accurate overall results and greater stability over the range of biomass values, the MFM approach also provides a suite of other biophysical structural outputs such as density, crown dimensions, LAI, height and sub-pixel scale fractions. Its explicit physical-basis and minimal ground data requirements are also more appropriate for larger area, multi-scene, multi-date applications with variable scene geometry and in high relief terrain. MFM thus warrants consideration for applications in mountainous and other, less complex terrain for purposes such as forest inventory updates, ecological modeling and terrestrial biomass and carbon monitoring studies.     

Use of coupled canopy structure dynamic and radiative transfer models to estimate biophysical canopy characteristics

Leaf area index (LAI) is a key variable for the understanding of several eco-physiological processes within a vegetation canopy. The LAI could thus provide vital information for the management of the environment and agricultural practices when estimated continuously over time and space thanks to remote sensing sensors.This study proposed a method to estimate LAI spatial and temporal variation based on multi-temporal remote sensing observations processed using a simple semi-mechanistic canopy structure dynamic model (CSDM) coupled with a radiative transfer model (RTM). The CSDM described the temporal evolution of the LAI as function of the accumulated daily air temperature as measured from classical ground meteorological stations.The retrieval performances were evaluated for two different data sets: first, a data set simulated by the RTM but taking into account realistic measurement conditions and uncertainties resulting from different error sources; second, an experimental data set acquired over maize crops the Blue Earth City area (USA) in 1998. Results showed that the proposed approach improved significantly the retrieval performances for LAI mainly by smoothing the residual errors associated to each individual observation. In addition it provides a way to describe in a continuous manner the LAI time course from a limited number of observations during the growth cycle.     

Evaluation of the MODIS LAI product using independent lidar-derived LAI: A case study in mixed conifer forest

This study presents an alternative assessment of the MODIS LAI product for a 58,000 ha evergreen needleleaf forest located in the western Rocky Mountain range in northern Idaho by using lidar data to model (R² = 0.86, RMSE = 0.76) and map LAI at higher resolution across a large number of MODIS pixels in their entirety. Moderate resolution (30 m) lidar-based LAI estimates were aggregated to the resolution of the 1-km MODIS LAI product and compared to temporally-coincident MODIS retrievals. Differences in the MODIS and lidar-derived values of LAI were grouped and analyzed by several different factors, including MODIS retrieval algorithm, sun/sensor geometry, and sub-pixel heterogeneity in both vegetation and terrain characteristics. Of particular interest is the disparity in the results when MODIS LAI was analyzed according to algorithm retrieval class. We observed relatively good agreement between lidar-derived and MODIS LAI values for pixels retrieved with the main RT algorithm without saturation for LAI LAI ≤ 4. Moreover, for the entire range of LAI values, considerable overestimation of LAI (relative to lidar-derived LAI) occurred when either the main RT with saturation or back-up algorithm retrievals were used to populate the composite product regardless of sub-pixel vegetation structural complexity or sun/sensor geometry. These results are significant because algorithm retrievals based on the main radiative transfer algorithm with or without saturation are characterized as suitable for validation and subsequent ecosystem modeling, yet the magnitude of difference appears to be specific to retrieval quality class and vegetation structural characteristics.     

A robust new method for the remote estimation of LAI in montane and boreal forests

Two promising techniques for estimating Leaf Area Index (LAI) using remote sensing are Linear Spectral Mixture Analysis (LSMA) and Modification of Spectral Vegetation Indices (MSVI). The Normalized Distance Method (ND), which uses principles employed by the LSMA and MSVI techniques, is introduced in this study. These three methods are applied to a region of montane forest in Kananaskis Country, Alberta, Canada, in order to estimate LAI. In situ measurements of LAI in 10 deciduous and 10 coniferous plots, and a SPOT‐4 image taken at the height of the growing season, provided test data that produced relationships for LAI in pure stands of either coniferous or deciduous vegetation using each of the three methods. All methods exhibited varying degrees of performance and demonstrated significant dependence on vegetation type. The ND method produced relationships with coefficients of determination (R ²) of 0.86 and 0.65 for coniferous and deciduous vegetation, respectively; the MSVI method (when using the adjusted Normalized Difference Vegetation Index) produced relationships with R ² values of 0.79 and 0.59 for coniferous and deciduous vegetation, respectively; and the LSMA technique produced relationships with R ² values of 0.83 and 0.0 for coniferous and deciduous vegetation, respectively.     

Lidar-based mapping of leaf area index and its use for validating GLOBCARBON satellite LAI product in a temperate forest of the southern USA

Lidar provides enhanced abilities to remotely map leaf area index (LAI) with improved accuracies. We aim to further explore the capability of discrete-return lidar for estimating LAI over a pine-dominated forest in East Texas, with a secondary goal to compare the lidar-derived LAI map and the GLOBCARBON moderate-resolution satellite LAI product. Specific problems we addressed include (1) evaluating the effects of analysts and algorithms on in-situ LAI estimates from hemispherical photographs (hemiphoto), (2) examining the effectiveness of various lidar metrics, including laser penetration, canopy height and foliage density metrics, to predict LAI, (3) assessing the utility of integrating Quickbird multispectral imagery with lidar for improving the LAI estimate accuracy, and (4) developing a scheme to co-register the lidar and satellite LAI maps and evaluating the consistency between them. Results show that the use of different analysts or algorithms in analyzing hemiphotos caused an average uncertainty of 0.35 in in-situ LAI, and that several laser penetration metrics in logarithm models were more effective than other lidar metrics, with the best one explaining 84% of the variation in the in-situ LAI (RMSE = 0.29 LAI). The selection of plot size and height threshold in calculating laser penetration metrics greatly affected the effectiveness of these metrics. The combined use of NDVI and lidar metrics did not significantly improve estimation over the use of lidar alone. We also found that mis-registration could induce a large artificial discrepancy into the pixelwise comparison between the coarse-resolution satellite and fine-resolution lidar-derived LAI maps. By compensating for a systematic sub-pixel shift error, the correlation between two maps increased from 0.08 to 0.85 for pines (n = 24 pixels). However, the absolute differences between the two LAI maps still remained large due to the inaccuracy in accounting for clumping effects. Overall, our findings imply that lidar offers a superior tool for mapping LAI at local to regional scales as compared to optical remote sensing, accuracies of lidar-estimate LAI are affected not only by the choice of models but also by the absolute accuracy of in-situ reference LAI used for model calibration, and lidar-derived LAI maps can serve as reliable references for validating moderate-resolution satellite LAI products over large areas.     

Retrieval of the canopy leaf area index in the BOREAS flux tower sites using linear spectral mixture analysis

This paper reports on the use of linear spectral mixture analysis for the retrieval of canopy leaf area index (LAI) in three flux tower sites in the Boreal Ecosystem-Atmosphere Study (BOREAS) southern study area: Old Black Spruce, Old Jack Pine, and Young Jack Pine (SOBS, SOJP, and SYJP). The data used were obtained by the Compact Airborne Spectrographic Imager (CASI) with a spatial resolution of 2 m in the winter of 1994. The convex geometry method was used to select the endmembers: sunlit crown, sunlit snow, and shadow. Along transects for these flux tower sites, the fraction of sunlit snow was found to have a higher correlation with the field-measured canopy LAI than the fraction of sunlit crown or the fraction of shadow. An empirical equation was obtained to describe the relation between canopy LAI and the fraction of sunlit snow. There is a strong correlation between the estimated LAI and the field-measured LAI along transects (with R² values of 0.54, 0.71, and 0.60 obtained for the SOBS, SYJP, and SOJP sites, respectively). The estimated LAI for the whole tower site is consistent with that obtained by the inversion of a canopy model in our previous study where values of 0.94, 0.92, and 0.63 were obtained for R² for the SOBS, SYJP and SOJP sites, respectively.The CASI 2-m summer data over the SOBS site was also employed to investigate the possibility of deriving canopy LAI from the summer data using linear mixture analysis. At a spatial resolution of 10 m, the correlation between the field-measured LAI and the estimated LAI along transects is small at R² less than 0.3, while R² increases to 0.6 at a spatial resolution of 30 m. The difficulty in canopy LAI retrieval from the summer data at a spatial resolution of 10 m is likely due to the variation of the understory reflectance across the scene, although spatial misregistration of the CASI data used may also be a possible contributing factor.     

Spatial non-stationarity and scale-dependency of prediction accuracy in the remote estimation of LAI over a tropical rainforest in Sulawesi, Indonesia

The technique of Geographically Weighted Regression (GWR) was used for estimation of Leaf Area Index (LAI) from remote sensing-based multi-spectral vegetation indices (VI) such as Normalized Difference Vegetation Index (NDVI), the mid-infrared corrected Normalized Difference Vegetation Index (NDVIc), Simple Ratio (SR), Soil-Adjusted Vegetation Index (SAVI) and Reduced Simple Ratio (RSR) in a region of equatorial rainforest in Central Sulavesi, Indonesia. The linear regressions between NDVI, NDVIc, SR, SAVI and RSR as explanatory variables and ground measurements of LAI at 166 plots as a dependent variable were produced using common modelling approach — Ordinary Least Squares (OLS) regression fitted to all data points, as well as GWR. Accuracy and precision statistics indicate that the GWR method made significantly better predictions of LAI in all simulations than OLS did. The relationships between LAI and the explanatory variables were found to be significantly spatially variable and scale-dependent. GWR has the potential to reveal local patterns in the spatial distribution of parameter estimates, it demonstrated sensitivity of the model's accuracy and performance to scale variation. The GWR approach enables finding the most appropriate scale for data analysis. This scale was different for each VI. The results suggest that spatial non-stationarity and scale-dependency in the relationship between LAI and remote sensing data has important implications for estimations of LAI based on empirical transfer functions.     

Upscaling ground observations of vegetation water content, canopy height, and leaf area index during SMEX02 using aircraft and Landsat imagery

Microwave-based remote sensing algorithms for mapping soil moisture are sensitive to water contained in surface vegetation at moderate levels of canopy cover. Correction schemes require spatially distributed estimates of vegetation water content at scales comparable to that of the microwave sensor footprint (10¹ to 10⁴ m). This study compares the relative utility of high-resolution (1.5 m) aircraft and coarser-resolution (30 m) Landsat imagery in upscaling an extensive set of ground-based measurements of canopy biophysical properties collected during the Soil Moisture Experiment of 2002 (SMEX02) within the Walnut Creek Watershed. The upscaling was accomplished using expolinear relationships developed between spectral vegetation indices and measurements of leaf area index, canopy height, and vegetation water content. Of the various indices examined, a Normalized Difference Water Index (NDWI), derived from near- and shortwave-infrared reflectances, was found to be least susceptible to saturation at high levels of leaf area index. With the aircraft data set, which did not include a short-wave infrared water absorption band, the Optimized Soil Adjusted Vegetation Index (OSAVI) yielded best correlations with observations and highest saturation levels. At the observation scale (10 m), LAI was retrieved from both NDWI and OSAVI imagery with an accuracy of 0.6, vegetation water content at 0.7 kg m⁻², and canopy height to within 0.2 m. Both indices were used to estimate field-scale mean canopy properties and variability for each of the intensive soil-moisture-sampling sites within the watershed study area. Results regarding scale invariance over the SMEX02 study area in transformations from band reflectance and vegetation indices to canopy biophysical properties are also presented.     

Airborne discrete-return LIDAR data in the estimation of vertical canopy cover, angular canopy closure and leaf area index

Remote sensing of forest canopy cover has been widely studied recently, but little attention has been paid to the quality of field validation data. Ecological literature has two different coverage metrics. Vertical canopy cover (VCC) is the vertical projection of tree crowns ignoring within-crown gaps. Angular canopy closure (ACC) is the proportion of covered sky at some angular range around the zenith, and can be measured with a field-of-view instrument, such as a camera. We compared field-measured VCC and ACC at 15° and 75° from the zenith to different LiDAR (Light Detection and Ranging) metrics, using several LiDAR data sets and comprehensive field data. The VCC was estimated to a high precision using a simple proportion of canopy points in first-return data. Confining to a maximum 15° scan zenith angle, the absolute root mean squared error (RMSE) was 3.7-7.0%, with an overestimation of 3.1-4.6%. We showed that grid-based methods are capable of reducing the inherent overestimation of VCC. The low scan angles and low power settings that are typically applied in topographic LiDARs are not suitable for ACC estimation as they measure in wrong geometry and cannot easily detect small within-crown gaps. However, ACC at 0-15° zenith angles could be estimated from LiDAR data with sufficient precision, using also the last returns (RMSE 8.1-11.3%, bias -6.1-+4.6%). The dependency of LiDAR metrics and ACC at 0-75° zenith angles was nonlinear and was modeled from laser pulse proportions with nonlinear regression with a best-case standard error of 4.1%. We also estimated leaf area index from the LiDAR metrics with linear regression with a standard error of 0.38. The results show that correlations between airborne laser metrics and different canopy field characteristics are very high if the field measurements are done with equivalent accuracy.     

A method to account for shoot scale clumping in coniferous canopy reflectance models

The three-dimensional structure of a coniferous shoot gives rise to multiple scattering of light between the needles of the shoot, causing the shoot spectral reflectance to differ from that of a flat leaf. Forest reflectance models based on the radiative transfer equation handle shoot level clumping by correcting the radiation attenuation coefficient with a clumping index. The clumping index causes a reduction in the interception of radiation by the canopy at a fixed leaf area index (LAI). In this study, we show how within-shoot multiple scattering is related to shoot scale clumping and derive a similar, but wavelength dependent, correction to the scattering coefficient. The results provide a method for integrating shoot structure into current radiative transfer equation based forest reflectance models. The method was applied to explore the effect of shoot scale clumping on canopy spectral reflectance using simple model canopies with a homogeneous higher level structure. The clumping of needles into shoots caused a wavelength dependent reduction in canopy reflectance, as compared to that of a leaf canopy with similar interception. This is proposed to be one reason why coniferous and broad-leaved canopies occupy different regions in the spectral space and exhibit different dependency of spectral vegetation indices on LAI.