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.     

Disentangling the relationships between plant pigments and the photochemical reflectance index reveals a new approach for remote estimation of carotenoid content

The Photochemical Reflectance Index (PRI) is used as an indicator of leaf and plant canopy photosynthetic efficiency. However, the photosynthetic efficiency-PRI relationship has been shown to be inconsistent over time, likely due to changes in foliar pigment content.We measured reflectance spectra and biochemical properties from 24 leaves of two deciduous tree species and acquired pigment and reflectance data from the Leaf Optical Properties EXperiment database for an additional nine species. These data were used as inputs for the PROSPECT-5 leaf optical model. We found measurements of PRI to be significantly (p < 0.05) correlated with chlorophyll content, carotenoid content, and the carotenoid/chlorophyll ratio. However, only the PRI-carotenoid/chlorophyll ratio relationship was consistent across all analyses. Two predictive equations were derived from PROSPECT-5 simulations: a curvilinear PRI model (PRI(clm)) predicted the carotenoid/chlorophyll ratio (r² = 0.99), and a linear model using the chlorophyll index (CI(lm)) predicted chlorophyll content (r² = 0.98). Multiplying PRI(clm) with CI(lm) canceled the influence of chlorophyll content on PRI(clm) and thus allowed for prediction of carotenoid content from 11 deciduous tree species (r² = 0.83). Our results confirm that the PRI is significantly influenced by chlorophyll and carotenoid pools and demonstrate a new approach for non-destructive estimation of leaf carotenoid content using the PRI. Because variation in foliar physiological status is known to relate to leaf carotenoid content and the carotenoid/chlorophyll ratio, convolving the PRI with a chlorophyll index is likely to be useful for understanding the photosynthetic performance of deciduous vegetation across a wide range of temporal periods, ranging from daily to seasonal time scales.     

Estimating dry matter content from spectral reflectance for green leaves of different species

Efficient and accurate detection of the temporal dynamics and spatial variations of leaf dry matter content would help monitor key properties and processes in vegetation and the wider ecosystem. However, leaf water content strongly absorbs at shortwave infrared wavelengths, reducing the signal from dry matter. The major objective of this study was to examine relationship between spectral reflectance of fresh leaves and the ratio of leaf dry mass to leaf area, across a wide range of species at the leaf scale. A narrow-band, normalized index combining two distinct wavebands centred at 1649 and 1722 nm achieved the highest overall performance and discriminatory power compared to either single band or first derivatives. The normalized index was evaluated using the PROSPECT (leaf optical properties spectra) simulated reflectance spectra and field measurements from the Leaf Optical Properties Experiment (LOPEX) data set. Both evaluations show that leaf dry matter contents were retrievable with R ² of 0.845 and 0.681 and regression slopes of 0.903 and 0.886. This study suggests that spectral reflectance measurements hold promise for the assessment of dry matter content for green leaves. Further investigation needs to be conducted to evaluate the effectiveness of this normalized index at canopy scales.     

Spectroscopic determination of leaf water content using continuous wavelet analysis

The gravimetric water content (GWC, %), a commonly used measure of leaf water content, describes the ratio of water to dry matter for each individual leaf. To date, the relationship between spectral reflectance and GWC in leaves is poorly understood due to the confounding effects of unpredictably varying water and dry matter ratios on spectral response. Few studies have attempted to estimate GWC from leaf reflectance spectra, particularly for a variety of species. This paper investigates the spectroscopic estimation of leaf GWC using continuous wavelet analysis applied to the reflectance spectra (350-2500 nm) of 265 leaf samples from 47 species observed in tropical forests of Panama. A continuous wavelet transform was performed on each of the reflectance spectra to generate a wavelet power scalogram compiled as a function of wavelength and scale. Linear relationships were built between wavelet power and GWC expressed as a function of dry mass (LWC_D) and fresh mass (LWC_F) in order to identify wavelet features (coefficients) that are most sensitive to changes in GWC. The derived wavelet features were then compared to three established spectral indices used to estimate GWC across a wide range of species.Eight wavelet features observed between 1300 and 2500 nm provided strong correlations with LWC_D, though correlations between spectral indices and leaf GWC were poor. In particular, two features captured amplitude variations in the broad shape of the reflectance spectra and three features captured variations in the shape and depth of dry matter (e.g., protein, lignin, cellulose) absorptions centered near 1730 and 2100 nm. The eight wavelet features used to predict LWC_D and LWC_F were not significantly different; however, predictive models used to determine LWC_D and LWC_F differed. The most accurate estimates of LWC_D and LWC_F obtained from a single wavelet feature showed root mean square errors (RMSEs) of 28.34% (R² = 0.62) and 4.86% (R² = 0.69), respectively. Models using a combination of features resulted in a noticeable improvement predicting LWC_D and LWC_F with RMSEs of 26.04% (R² = 0.71) and 4.34% (R² = 0.75), respectively. These results provide new insights into the role of dry matter absorption features in the shortwave infrared (SWIR) spectral region for the accurate spectral estimation of LWC_D and LWC_F. This emerging spectral analytical approach can be applied to other complex datasets including a broad range of species, and may be adapted to estimate basic leaf biochemical elements such as nitrogen, chlorophyll, cellulose, and lignin.     

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.     

Mapping tree and shrub leaf area indices in an ombrotrophic peatland through multiple endmember spectral unmixing

Leaf area index (LAI) is an important parameter used by most process-oriented ecosystem models. LAI of forest ecosystems has routinely been mapped using spectral vegetation indices (SVI) derived from remote sensing imagery. The application of SVI-based approaches to map LAI in peatlands presents a challenge, mainly due to peatlands characteristic multi-layer canopy comprising shrubs and open, discontinuous tree canopies underlain by a continuous ground cover of different moss species, which reduces the greenness contrast between the canopy and the background.Our goal is to develop a methodology to map tree and shrub LAI in peatlands and similar ecosystems based on multiple endmember spectral mixture analysis (MESMA). This new mapping method is validated using LAI field measurements from a precipitation-fed (ombrotrophic) peatland near Ottawa, Ontario, Canada. We demonstrate first that three commonly applied SVI are not suitable for tree and shrub LAI mapping in ombrotrophic peatlands. Secondly, we demonstrate for a three-endmember model the limitations of traditional linear spectral mixture analysis (SMA) due to the unique and widely varying spectral characteristics of Sphagnum mosses, which are significantly different from vascular plants. Next, by using a geometric-optical radiative transfer model, we determine the nature of the equation describing the empirical relationship between shadow fraction and tree LAI using nonlinear ordinary least square (OLS) regression. We then apply this equation to describe the empirical relationships between shadow and shrub fractions obtained from mixture decomposition with SMA and MESMA, respectively, and tree and shrub LAI, respectively. Less accurate fractions obtained from SMA result in weaker relationships between shadow fraction and tree LAI (R² = 0.61) and shrub fraction and shrub LAI (R² = 0.49) compared to the same relationships based on fractions obtained from MESMA with R² = 0.75 and R² = 0.68, respectively. Cross-validation of tree LAI (R² = 0.74; RMSE = 0.48) and shrub LAI (R² = 0.68; RMSE = 0.42) maps using fractions from MESMA shows the suitability of this approach for mapping tree and shrub LAI in ombrotrophic peatlands. The ability to account for a spectrally varying, unique Sphagnum moss ground cover during mixture decomposition and a two layer canopy is particularly important.     

Multiple shadow fractions in spectral mixture analysis of a cotton canopy

Shadows are being used more frequently to estimate plant canopy biophysical characteristics. Typically, a zero value is assumed or a threshold value is derived from histogram analysis of imagery to determine the shadow endmember (EM). Here, two distinct shadow EMs were measured in situ for use in spectral mixture analysis of a cotton canopy on five dates in 2003. The four EMs used in the analysis were: sunlit green leaf, sunlit dry soil, self-shadowed leaf, shadowed dry soil. This 4-EM model was compared to a 3-EM model where a zero-value shade EM was used for unmixing with the two sunlit EMs. Multiple endmember spectral mixture analysis (MESMA) was used to allow EM composition to vary across each scene. The analysis and EMs were applied to fine-scale hyperspectral image data collected in the wavelength range, 440 to 810 nm. Ground data collected included percent cover, height, SPAD (a measure of leaf greenness), and chlorophyll a concentration. The normalized difference vegetation index (NDVI) was also compared to the unmixing results. Regression analysis showed that NDVI was equal to the 4-EM model for estimation of percent cover (r² = 0.95, RMSE = 6.6) although the NDVI y-intercept was closer to zero. The 4-EM model was best for estimating height (r² = 0.79, RMSE = 0.07 m) and chlorophyll a concentration (r² = 0.46, RMSE = 7.0 μg/cm²). The 3-EM model and NDVI performed poorly when estimating chlorophyll a concentration. Inclusion of two distinct shadow EMs in the model improved relationships to crop biophysical parameters and was better than assuming one, zero-value shade EM. Since MESMA operates at the pixel level and allows variable EM assignment to each pixel, mapping the spatial variability of shadows and other variables of interest is possible, providing a powerful input to canopy and ecosystem models as well as precision farming.     

PROSPECT-4 and 5: Advances in the leaf optical properties model separating photosynthetic pigments

The PROSPECT leaf optical model has, to date, combined the effects of photosynthetic pigments, but a finer discrimination among the key pigments is important for physiological and ecological applications of remote sensing. Here we present a new calibration and validation of PROSPECT that separates plant pigment contributions to the visible spectrum using several comprehensive datasets containing hundreds of leaves collected in a wide range of ecosystem types. These data include leaf biochemical (chlorophyll a, chlorophyll b, carotenoids, water, and dry matter) and optical properties (directional-hemispherical reflectance and transmittance measured from 400 nm to 2450 nm). We first provide distinct in vivo specific absorption coefficients for each biochemical constituent and determine an average refractive index of the leaf interior. Then we invert the model on independent datasets to check the prediction of the biochemical content of intact leaves. The main result of this study is that the new chlorophyll and carotenoid specific absorption coefficients agree well with available in vitro absorption spectra, and that the new refractive index displays interesting spectral features in the visible, in accordance with physical principles. Moreover, we improve the chlorophyll estimation (RMSE = 9 µg/cm²) and obtain very encouraging results with carotenoids (RMSE = 3 µg/cm²). Reconstruction of reflectance and transmittance in the 400-2450 nm wavelength domain using PROSPECT is also excellent, with small errors and low to negligible biases. Improvements are particularly noticeable for leaves with low pigment content.     

Towards red‐edge positions less sensitive to canopy biophysical parameters for leaf chlorophyll estimation using properties optique spectrales des feuilles (PROSPECT) and scattering by arbitrarily inclined leaves (SAILH) simulated data

Several methods for extracting the chlorophyll sensitive red‐edge position (REP) from hyperspectral data are reported in literature. This study is a continuation of a recent paper published as ‘A new technique for extracting the red edge position from hyperspectral data: the linear extrapolation method’. The method was validated experimentally for estimation of foliar nitrogen concentrations of rye, maize and mixed grass/herb. The objective of this study was to test the utility of the linear extrapolation method under different conditions including variable canopy biophysical parameters, solar zenith angle, sensor noise and spectral bandwidth. REPs were extracted from synthetic canopy spectra that were simulated using properties optique spectrales des feuilles (PROSPECT) and scattering by arbitrarily inclined leaves (SAILH) radiative transfer models. REPs extracted by the linear extrapolation method involving wavebands at 680, 694, 724 and 760 nm produced the highest correlation (R ² = 0.75) with leaf chlorophyll content with minimal effects of leaf and canopy biophysical confounders (leaf area index, leaf inclination distribution and leaf dry matter content) compared to traditional techniques including the linear interpolation, inverted Gaussian modelling and polynomial fitting techniques. In addition, the new technique is insensitive to changes in solar zenith angle. However, the advantage of using the linear extrapolation method compared to the various alternative methods diminishes with increasing sensor noise and decreasing spectral resolution. In summary, the linear extrapolation technique confirms its high potential for leaf chlorophyll estimation. The efficacy of the technique under field conditions needs to be established.     

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.     

Assessing structural effects on PRI for stress detection in conifer forests

The retrieval of indicators of vegetation stress from remote sensing imagery is an important issue for the accurate assessment of forest decline. The Photochemical Reflectance Index (PRI) has been demonstrated as a physiological index sensitive to the epoxidation state of the xanthophyll cycle pigments and to photosynthetic efficiency, serving as a proxy for short-term changes in photosynthetic activity, stress condition, and pigment absorption, but highly affected by illumination conditions, viewing geometry and canopy structure. In this study, a diurnal airborne campaign was conducted over Pinus sylvestris and Pinus nigra forest areas with the Airborne Hyperspectral Scanner (AHS) to evaluate the effects of canopy structure on PRI when used as an indicator of stress in a conifer forest. The AHS airborne sensor was flown at two times (8:00 GMT and 12:00 GMT) over forest areas under varying field-measured stress levels, acquiring 2 m spatial resolution imagery in 80 spectral bands in the 0.43-12.5 μm spectral range. Five formulations of PRI (based on R₅₃₁ as a xanthophyll-sensitive spectral band) were calculated using different reference wavelengths, such as PRI₅₇₀ (reference band R_REF = R₅₇₀), and the PRI modifications PRI_m1 (R_REF = R₅₁₂), PRI_m2 (R_REF = R₆₀₀), PRI_m3 (R_REF = R₆₇₀), and PRI_m4 (R_REF = R₅₇₀, R₆₇₀), along with other structural indices such as NDVI, SR, OSAVI, MSAVI and MTVI2. In addition, thermal bands were used for the retrieval of the land surface temperature. A radiative transfer modeling method was conducted using the LIBERTY and INFORM models to assess the structural effects on the PRI formulations proposed, studying the sensitivity of PRI_m indices to detect stress levels while minimizing the effects caused by the conifer architecture. The PRI indices were related to stomatal conductance, xanthophyll epoxidation state (EPS) and crown temperature. The modeling analysis showed that the coefficient of variation (CV) for PRI was 50%, whereas the CV for PRIm1 (band R512 as a reference) was only 20%. Simulation and experimental results demonstrated that PRI_m1 (R_REF = R₅₁₂) was less sensitive than PRI (R_REF = R₅₇₀) to changes in Leaf Area Index (LAI) and tree densities. PRI₅₁₂ was demonstrated to be sensitive to EPS at both leaf (r² = 0.59) and canopy level (r² = 0.40), yielding superior performance than PRI₅₇₀ (r² = 0.21) at the canopy level. In addition, PRI₅₁₂ was significantly related to water stress indicators such as stomatal conductance (Gs; r² = 0.45) and water potential (Ψ; r² = 0.48), yielding better results than PRI₅₇₀ (Gs, r² = 0.21; Ψ, r² = 0.21) due to the structural effects found on the PRI₅₇₀ index at the canopy level.     

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 live fuel moisture content from remotely sensed reflectance

Fuel moisture content (FMC) is used in forest fire danger models to characterise the moisture status of the foliage. FMC expresses the amount of water in a leaf relative to the amount of dry matter and differs from measures of leaf water content which express the amount of water in a leaf relative to its area. FMC is related to both leaf water content and leaf dry matter content, and the relationships between FMC and remotely sensed reflectance will therefore be affected by variation in both leaf biophysical properties. This paper uses spectral reflectance data from the Leaf Optical Properties EXperiment (LOPEX) and modelled data from the Prospect leaf reflectance model to examine the relationships between FMC, leaf equivalent water thickness (EWT) and a range of spectral vegetation indices (VI) designed to estimate leaf and canopy water content. Significant correlations were found between FMC and all of the selected vegetation indices for both modelled and measured data, but statistically stronger relationships were found with leaf EWT; overall, the water index (WI) was found to be most strongly correlated with FMC. The accuracy of FMC estimation was very low when the global range of FMC was examined, but for a restricted range of 0-100%, FMC was estimated with a root-mean-square error (RMSE) of 15% in the model simulations and 51% with the measured data. The paper shows that the estimation of live FMC from remotely sensed vegetation indices is likely to be problematic when there is variability in both leaf water content and leaf dry matter content in the target leaves. Estimating FMC from remotely sensed data at the canopy level is likely to be further complicated by spatial and temporal variations in leaf area index (LAI). Further research is required to assess the potential of canopy reflectance model inversion to estimate live fuel moisture content where a priori information on vegetation properties may be used to constrain the inversion process.     

Water content estimation in vegetation with MODIS reflectance data and model inversion methods

Statistical and radiative-transfer physically based studies have previously demonstrated the relationship between leaf water content and leaf-level reflectance in the near-infrared spectral region. The successful scaling up of such methods to the canopy level requires modeling the effect of canopy structure and viewing geometry on reflectance bands and optical indices used for estimation of water content, such as normalized difference water index (NDWI), simple ratio water index (SRWI) and plant water index (PWI). This study conducts a radiative transfer simulation, linking leaf and canopy models, to study the effects of leaf structure, dry matter content, leaf area index (LAI), and the viewing geometry, on the estimation of leaf equivalent water thickness from canopy-level reflectance. The applicability of radiative transfer model inversion methods to MODIS is studied, investigating its spectral capability for water content estimation. A modeling study is conducted, simulating leaf and canopy MODIS-equivalent synthetic spectra with random input variables to test different inversion assumptions. A field sampling campaign to assess the investigated simulation methods was undertaken for analysis of leaf water content from leaf samples in 10 study sites of chaparral vegetation in California, USA, between March and September 2000. MODIS reflectance data were processed from the same period for equivalent water thickness estimation by model inversion linking the PROSPECT leaf model and SAILH canopy reflectance model. MODIS reflectance data, viewing geometry values, and LAI were used as inputs in the model inversion for estimation of leaf equivalent water thickness, dry matter, and leaf structure. Results showed good correlation between the time series of MODIS-estimated equivalent water thickness and ground measured leaf fuel moisture (LFM) content (r²=0.7), demonstrating that these inversion methods could potentially be used for global monitoring of leaf water content in vegetation.     

Spectral response of benthic diatoms with different sediment backgrounds

Remote sensing of terrestrial vegetation uses a wide range of vegetation indices (VIs) to monitor plant characteristics, but these indices can be very sensitive to canopy background reflectance. This study investigated background influences on VIs applied to intertidal microphytobenthos, using a synthetic spectral library constituted by a spectral combination of three contrasting types of sediment (sand, fine sand, and mud) and reflectance spectra of benthic diatom monospecific cultures obtained in controlled conditions. The spectral database exhibited, for the same biomass range (3-182 mg chlorophyll a m^− 2), marked differences in albedo and spectral contrast linked to sediment variability in water content, grain size, and organic matter content. Several VIs were evaluated, from ratios using visible and near infrared wavelengths, to hyperspectral indices (derivative analysis, continuum removal). Among the ratios, the Normalized Difference Vegetation Index (NDVI) appeared less sensitive to background effects than VIs with soil corrections such as the Perpendicular Vegetation Index (PVI), the Soil-Adjusted Vegetation Index (SAVI), the Modified second Soil-Adjusted Vegetation Index (MSAVI2) or the Transformed Soil-Adjusted Vegetation Index (TSAVI). The lower efficacy of soil-corrected VIs may be explained by the structural differences and optical behavior of soil vs. canopies compared to sediment vs. microphytobenthos biofilms. The background effects were minimized using Modified Gaussian Model indices at 632 nm and 675 nm, and the second derivative at 632 nm, while poor results were obtained with the red-edge inflection point (REIP) and the second derivative at 675 nm. The least sensitive index was the Phytobenthos Index which is very similar to the NDVI, but uses a red wavelength at 632 nm instead of 675 nm, to account for the absorption by chlorophyll c. The modified NDVI₇₀₅, where the 705 nm wavelength replaces the red band, showed moderate background sensitivity. Moreover, the NDVI₇₀₅ and the Phytobenthos Index have the additional relevant property of being less sensitive to the index saturation response with increasing biomass. Unfortunately, these VIs cannot be applied to broad-band multispectral satellite images, and require sensors with a hyperspectral resolution. Nevertheless, this study showed that the background influence was not a limitation to applying the ubiquitous NDVI to map intertidal microphytobenthos using multispectral satellite images.     

Satellite passive microwave remote sensing for monitoring global land surface phenology

Vegetation phenology characterizes seasonal life-cycle events that influence the carbon cycle and land-atmosphere water and energy exchange. We analyzed global phenology cycles over a six year record (2003-2008) using satellite passive microwave remote sensing based Vegetation Optical Depth (VOD) retrievals derived from daily time series brightness temperature (T_b) measurements from the Advanced Microwave Scanning Radiometer on EOS (AMSR-E) and other ancillary data inputs. The VOD parameter derives vegetation canopy attenuation at a given microwave frequency (18.7 GHz) and varies with canopy height, density, structure and water content. An error sensitivity analysis indicates that the retrieval algorithm can resolve the VOD seasonal cycle over a majority of global vegetated land areas. The VOD results corresponded favorably (p < 0.01) with vegetation indices (VIs) and leaf area index (LAI) information from satellite optical-infrared (MODIS) remote sensing, and phenology cycles determined from a simple bioclimatic growing season index (GSI) for over 82% of the global domain. Lower biomass land cover classes (e.g. savannas) show the highest correlations (R = 0.66), with reduced correspondence at higher biomass levels (0.03 < R < 0.51) and higher correlations for homogeneous land cover areas (0.41 < R < 0.83). The VOD results display a unique end-of-season signal relative to VI and LAI series, and may reflect microwave sensitivity to the timing of vegetation biomass depletion (e.g. leaf abscission) and associated changes in canopy water content (e.g. dormancy preparation). The VOD parameter is independent of and synergistic with optical-infrared remote sensing based vegetation metrics, and contributes to a more comprehensive view of land surface phenology.     

Foliar spectral properties following leaf clipping and implications for handling techniques

After leaves are clipped their reflectance properties change over time at variable rates. Spectral change can in part be attributed to the changing water content of the leaf, which affects absorption in the VIS, NIR and the SWIR. Maintaining water volume within samples has been the motivation behind many leaf handling techniques. This study has assessed the time constraints between leaf collection and spectral measurement. Specifically the relationship between leaf water content and foliar spectra (350-2500 nm) was examined over time for five tropical trees (common guava (Psidium guajava), purple guava (Psidium littorale), weeping fig (Ficus benjamina), floss silk (Chorisia speciosa), and coffee (Coffea arabica)). This investigation was carried for leaves wrapped with moist gauze around their petiole (treatment leaves) and leaves with no treatment. Spectral measurements and mass measurements were repeated for each leaf once every hour for the first 12 h, then every 4-6 h for 18 h, followed by one measurement after 12 h, and finally once a day until the control samples became air-dry. Foliar reflectance in the visible spectrum was not immediately responsive to water content changes and did not change until wilting of the leaf was observed. The NIR and SWIR wavelength regions were affected immediately by small changes in water content. Thus, by the time wilting was first observed the NIR and SWIR foliar reflectance differed considerably from corresponding fresh leaf reflectance. No common time limit could be observed for leaf clipping and reflectance measurement. Leaves have a variety of water contents and dehydration rates hence measurement time constraints are dependent on the properties of the leaf or species. Rather than using a time limit it is recommended that leaf handling techniques be based upon managing leaf water content and leaf structure. The results of this study indicate that leaves with petioles wrapped in moist paper towel and placed within plastic bags will maintain leaf reflectance longer than equivalent leaves without treatment; samples tested here lasted a minimum of 7 days. θ and D indices (“angle difference” and “root mean square difference”, respectively) revealed a stronger relationship between leaf water content and spectral shape than between leaf water and raw reflectance magnitude. The ratio of 1187/1096 nm, when compared with θ and D indices and individual reflectance bands, showed the highest coefficient of determination with leaf water content (r² = 0.952).     

Influence of woody elements of a Norway spruce canopy on nadir reflectance simulated by the DART model at very high spatial resolution

A detailed sensitivity analysis investigating the effect of woody elements introduced into the Discrete Anisotropic Radiative Transfer (DART) model on the nadir bidirectional reflectance factor (BRF) for a simulated Norway spruce canopy was performed at a very high spatial resolution (modelling resolution 0.2 m, output pixel size 0.4 m). We used such a high resolution to be able to parameterize DART in an appropriate way and subsequently to gain detailed understanding of the influence of woody elements contributing to the radiative transfer within heterogeneous canopies. Three scenarios were studied by modelling the Norway spruce canopy as being composed of i) leaves, ii) leaves, trunks and first order branches, and finally iii) leaves, trunks, first order branches and small woody twigs simulated using mixed cells (i.e. cells approximated as composition of leaves and/or twigs turbid medium, and large woody constituents). The simulation of each scenario was performed for 10 different canopy closures (CC = 50-95%, in steps of 5%), 25 leaf area index (LAI = 3.0-15.0 m² m^− 2, in steps of 0.5 m² m^− 2), and in four spectral bands (centred at 559, 671, 727, and 783 nm, with a FWHM of 10 nm). The influence of woody elements was evaluated separately for both, sunlit and shaded parts of the simulated forest canopy, respectively. The DART results were verified by quantifying the simulated nadir BRF of each scenario with measured Airborne Imaging Spectroradiometer (AISA) Eagle data (pixel size of 0.4 m). These imaging spectrometer data were acquired over the same Norway spruce stand that was used to parameterise the DART model.The Norway spruce canopy modelled using the DART model consisted of foliage as well as foliage including robust woody constituents (i.e. trunks and branches). All results showed similar nadir BRF for the simulated wavelengths. The incorporation of small woody parts in DART caused the canopy reflectance to decrease about 4% in the near-infrared (NIR), 2% in the red edge (RE) and less than 1% in the green band. The canopy BRF of the red band increased by about 2%. Subsequently, the sensitivity on accounting for woody elements for two spectral vegetation indices, the normalized difference vegetation index (NDVI) and the angular vegetation index (AVI), was evaluated. Finally, we conclude on the importance of including woody elements in radiative transfer based approaches and discuss the applicability of the vegetation indices as well as the physically based inversion approaches to retrieve the forest canopy LAI at very high spatial resolution.