A new instrument for passive remote sensing: 2. Measurement of leaf and canopy reflectance changes at 531 nm and their relationship with photosynthesis and chlorophyll fluorescence

A previously described passive remote sensing fluorimeter (see companion paper) was modified to detect changes in the reflectance of vegetation. The utility of this remote sensing technique to measure the Physiological Reflectance Index (PRI) is shown at both leaf level under laboratory conditions and at the canopy level in the field. PRI, defined as the relative changes in reflectance at 531 nm with respect to those at 570 nm (PRI=R531−R570/R531+R570), is related to xanthophyll-related, dynamic changes of non-photochemical quenching of chlorophyll fluorescence. The robustness of this relationship by simultaneous remote sensing of PRI and chlorophyll fluorescence is strengthened. At the leaf level, the existence of two kinetically distinct components of PRI is shown. A fast (within seconds) component that is partly attributed to ΔpH induced chloroplast shrinkage, and a slow (within minutes), main component that is related to xanthophyll de-epoxidation, as demonstrated by its disappearance in the presence of DTT. Overall, PRI correlated better with non-photochemical quenching of chlorophyll fluorescence (NPQ) than with any other measured parameter, including the photochemical efficiency of PSII. Finally, at the canopy level and under field conditions, it is shown that PRI can be a useful tool for remote sensing of water stress in grapevines.     

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.     

Remote sensing of photosynthetic light-use efficiency across two forested biomes: Spatial scaling

Eddy covariance (EC) measurements have greatly advanced our knowledge of carbon exchange in terrestrial ecosystems. However, appropriate techniques are required to upscale these spatially discrete findings globally. Satellite remote sensing provides unique opportunities in this respect, but remote sensing of the photosynthetic light-use efficiency (ε), one of the key components of Gross Primary Production, is challenging. Some progress has been made in recent years using the photochemical reflectance index, a narrow waveband index centered at 531 and 570 nm. The high sensitivity of this index to various extraneous effects such as canopy structure, and the view observer geometry has so far prevented its use at landscape and global scales. One critical aspect of upscaling PRI is the development of generic algorithms to account for structural differences in vegetation. Building on previous work, this study compares the differences in the PRI: ? relationship between a coastal Douglas-fir forest located on Vancouver Island, British Columbia, and a mature Aspen stand located in central Saskatchewan, Canada. Using continuous, tower-based observations acquired from an automated multi-angular spectro-radiometer (AMSPEC II) installed at each site, we demonstrate that PRI can be used to measure ? throughout the vegetation season at the DF-49 stand (r² = 0.91, p < 0.00) as well as the deciduous site (r² = 0.88, p < 0.00). It is further shown that this PRI signal can be also observed from space at both sites using daily observations from the Moderate Resolution Imaging Spectro-radiometer (MODIS) and a multi-angular implementation of atmospheric correction (MAIAC) (r² = 0.54 DF-49; r² = 0.63 SOA; p < 0.00). By implementing a simple hillshade model derived from airborne light detection and ranging (LiDAR) to approximate canopy shadow fractions (α_s), it is further demonstrated that the differences observed in the relationship between PRI and ε at DF-49 and SOA can be attributed largely to differences in α_s. The findings of this study suggest that algorithms used to separate physiological from extraneous effects in PRI reflectance may be more broadly applicable and portable across these two climatically and structurally different biome types, when the differences in canopy structure are known.     

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.     

The impact of soil reflectance on the quantification of the green vegetation fraction from NDVI

The green vegetation fraction (Fg) is an important climate and hydrologic model parameter. A common method to calculate Fg is to create a simple linear mixing model between two NDVI endmembers: bare soil NDVI (NDVI_o) and full vegetation NDVI (NDVI_∞). Usually it is assumed that NDVI_o is close to zero (NDVI_o ∼ 0.05) and is generally chosen from the lowest observed NDVI values. However, the mean soil NDVI computed from 2906 samples is much larger (NDVI = 0.2) and is highly variable (standard deviation = 0.1). We show that the underestimation of NDVI_o yields overestimations of Fg. The largest errors occur in grassland and shrubland areas. Using parameters for NDVI_o and NDVI_∞ derived from global scenes yields overestimations of Fg that are larger than 0.2 for the majority of U.S. land cover types when pixel NDVI values are 0.2 < NDVI_pixel < 0.4. When using conterminous U.S. scenes to derive NDVI_o and NDVI_∞, the overestimation is less (0.10-0.17 for 0.2 < NDVI_pixel < 0.4). As a result, parts of the conterminous U.S. are affected at different times of the year depending on the local seasonal NDVI cycle. We propose using global databases of NDVI_o along with information on historical NDVI_pixel values to compute a statistically most-likely estimate of Fg. Using in situ measurements made at the Sevilleta LTER, we show that this approach yields better estimates of Fg than using global invariant NDVI_o values estimated from whole scenes. At the two studied sites, the Fg estimate was adjusted by 52% at the grassland and 86% at the shrubland. More significant advances will require information on spatial distribution of soil reflectance.     

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.     

Separating physiologically and directionally induced changes in PRI using BRDF models

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

An assessment of photosynthetic light use efficiency from space: Modeling the atmospheric and directional impacts on PRI reflectance

Estimation of photosynthetic light use efficiency (ε) from satellite observations is an important component of climate change research. The photochemical reflectance index, a narrow waveband index based on the reflectance at 531 and 570 nm, allows sampling of the photosynthetic activity of leaves; upscaling of these measurements to landscape and global scales, however, remains challenging. Only a few studies have used spaceborne observations of PRI so far, and research has largely focused on the MODIS sensor. Its daily global coverage and the capacity to detect a narrow reflectance band at 531 nm make it the best available choice for sensing ε from space. Previous results however, have identified a number of key issues with MODIS-based observations of PRI. First, the differences between the footprint of eddy covariance (EC) measurements and the MODIS footprint, which is determined by the sensor's observation geometry make a direct comparison between both data sources challenging and second, the PRI reflectance bands are affected by atmospheric scattering effects confounding the existing physiological signal. In this study we introduce a new approach for upscaling EC based ε measurements to MODIS. First, EC-measured ε values were “translated” into a tower-level optical PRI signal using AMSPEC, an automated multi-angular, tower-based spectroradiometer instrument. AMSPEC enabled us to adjust tower-measured PRI values to the individual viewing geometry of each MODIS overpass. Second, MODIS data were atmospherically corrected using a Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm, which uses a time series approach and an image-based rather than pixel-based processing for simultaneous retrievals of atmospheric aerosol and surface bidirectional reflectance (BRDF). Using this approach, we found a strong relationship between tower-based and spaceborne reflectance measurements (r² = 0.74, p < 0.01) throughout the vegetation period of 2006. Swath (non-gridded) observations yielded stronger correlations than gridded data (r² = 0.58, p < 0.01) both of which included forward and backscatter observations. Spaceborne PRI values were strongly related to canopy shadow fractions and varied with different levels of ε. We conclude that MAIAC-corrected MODIS observations were able to track the site-level physiological changes from space throughout the observation period.     

PHOTOSYNSAT, photosynthesis from space: Theoretical foundations of a satellite concept and validation from tower and spaceborne data

We develop herein the theoretical foundations for a new satellite concept, utilizing multi-angle, along track spectral measurements to infer photosynthesis and gross primary production, at the landscape level over time. We validate the theory using both tower and space-borne sensors. The concept, originated in Hall et al. (2008), and Hilker et al. (2008a) and is based on two principles: (1) The first derivative of the photochemical reflectance index (PRI) with respect to shadow fraction viewed by the sensor ∂PRI/∂α_s, is proportional to light-use efficiency ε. (2) This behavior can be shown both theoretically and empirically to be independent of vegetation structure and optical properties. These two principles provide the basis for a robust photosynthesis algorithm that can be applied consistently both spatially and temporally. We develop the general theoretical concept using a canopy reflectance model that incorporates a dependence of leaf reflectance on illumination strength, permitting the leaf reflectance at 531 nm to depend on the intensity of photosynthetic down-regulation. Using this model we are able to show that using PRI alone to infer ε is confounded by the shadow fraction viewed by a sensor, the PRI value in a non-down-regulated physiological state, and the sunlit canopy reflectance. We are able to demonstrate that these difficulties are mitigated by using ∂PRI/∂α_s—not PRI—as the primary measure of canopy level ε. We demonstrate our concept using tower and satellite data acquired over three years, in two distinct biomes and vegetation types to show that PRI/∂α_s and ε are related by a single function. Building on these ideas we propose the development of a new satellite concept that can utilize a spatially and temporally robust algorithm to map photosynthesis at landscape scales and its temporal variation.     

Normalized difference spectral indices for estimating photosynthetic efficiency and capacity at a canopy scale derived from hyperspectral and CO2 flux measurements in rice

We explored simple and useful spectral indices for estimating photosynthetic variables (radiation use efficiency and photosynthetic capacity) at a canopy scale based on seasonal measurements of hyperspectral reflectance, ecosystem CO₂ flux, and plant and micrometeorological variables. An experimental study was conducted over the simple and homogenous ecosystem of an irrigated rice field. Photosynthetically active radiation absorbed by the canopy (APAR), the canopy absorptivity of APAR (fAPAR), net ecosystem exchange of CO₂ (NEE_CO2) gross primary productivity (GPP), photosynthetic capacity at the saturating APAR (P_max), and three parameters of radiation use efficiency (ε_N: NEE_CO2/APAR; ε_G: GPP/APAR; φ: quantum efficiency) were derived from the data set. Based on the statistical analysis of relationships between these ecophysiological variables and reflectance indicators such as normalized difference spectral indices (NDSI[i,j]) using all combinations of two wavelengths (i and j nm), we found several new indices that would were more effective than conventional spectral indices such as photochemical reflectance index (PRI) and normalized difference vegetation index (NDVI = NDSI[near-infrared, red]). ε_G was correlated well with NDSI[710, 410], NDSI[710, 520], and NDSI[530, 550] derived from nadir measurements. φ was best correlated with NDSI[450, 1330]. NDSI[550, 410] and NDSI[720, 420] had a consistent linear relationships with fAPAR throughout the growing season, whereas conventional indices such as NDVI showed very different relationships before and after heading. Off-nadir measurements were more closely related to the efficiency parameters than nadir measurements. Our results provide useful insights for assessing plant productivity and ecosystem CO₂ exchange, using a wide range of available spectral data as well as useful information for designing future sensors for ecosystem observations.     

Remote sensing of sunlight-induced chlorophyll fluorescence and reflectance of Scots pine in the boreal forest during spring recovery

A measurement campaign to assess the feasibility of remote sensing of sunlight-induced chlorophyll fluorescence (ChlF) from a coniferous canopy was conducted in a boreal forest study site (Finland). A Passive Multi-wavelength Fluorescence Detector (PMFD) sensor, developed in the LURE laboratory, was used to obtain simultaneous measurements of ChlF in the oxygen absorption bands, at 687 and 760 nm, and a reflectance index, the PRI (Physiological Reflectance Index), for a month during spring recovery. When these data were compared with active fluorescence measurements performed on needles they revealed the same trend. During sunny days fluorescence and reflectance signals were found to be strongly influenced by shadows associated with the canopy structure. Moreover, chlorophyll fluorescence variations induced by rapid light changes (due to transient cloud shadows) were found to respond more quickly and with larger amplitude under summer conditions compared to those obtained under cold acclimation conditions. In addition, ChlF at 760 nm was observed to increase with the chlorophyll content. During this campaign, the CO₂ assimilation was measured at the forest canopy level and was found remarkably well correlated with the PRI index.     

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.     

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.     

Estimation of leaf and canopy water content in poplar plantations by means of hyperspectral indices and inverse modeling

This study investigates the applicability of empirical and radiative transfer models to estimate water content at leaf and landscape level. The main goal is to evaluate and compare the accuracy of these two approaches for estimating leaf water content by means of laboratory reflectance/transmittance measurements and for mapping leaf and canopy water content by using airborne Multispectral Infrared and Visible Imaging Spectrometer (MIVIS) data acquired over intensive poplar plantations (Ticino, Italy).At leaf level, we tested the performance of different spectral indices to estimate leaf equivalent water thickness (EWT) and leaf gravimetric water content (GWC) by using inverse ordinary least squares (OLS) regression, and reduced major axis (RMA) regression. The analysis showed that leaf reflectance is related to changes in EWT rather than GWC, with best results obtained by using RMA regression by exploiting the spectral index related to the continuum removed area of the 1200 nm water absorption feature with an explained variance of 61% and prediction error of 6.6%. Moreover, we inverted the PROSPECT leaf radiative transfer model to estimate leaf EWT and GWC and compared the results with those obtained by means of empirical models. The inversion of this model showed that leaf EWT can be successfully estimated with no prior information with mean relative errors of 14% and determination coefficient of 0.65. Inversion of the PROSPECT model showed some difficulties in the simultaneous estimation of leaf EWT and dry matter content, which led to large errors in GWC estimation.At landscape level with MIVIS data, we tested the performance of different spectral indices to estimate canopy water per unit ground area (EWT_canopy). We found a relative error of 20% using a continuum removed spectral index around 1200 nm. Furthermore, we used a model simulation to evaluate the possibility of applying empirical models based on appositely developed MIVIS double ratios to estimate mean leaf EWT at landscape level (). It is shown that combined indices (double ratios) yielded significant results in estimating leaf EWT at landscape level by using MIVIS data (with errors around 2.6%), indicating their potential in reducing the effects of LAI on the recorded signal. The accuracy of the empirical estimation of EWT_canopy and was finally compared with that obtained from inversion of the PROSPECT + SAILH canopy reflectance model to evaluate the potential of both methods for practical applications. A relative error of 27% was found for EWT_canopy and an overestimation of leaf with relative errors around 19%.Results arising from this remote sensing application support the robustness of hyperspectral regression indices for estimating water content at both leaf and landscape level, with lower relative errors compared to those obtained from inversion of leaf and 1D canopy radiative transfer models.     

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.     

Remote sensing of tundra gross ecosystem productivity and light use efficiency under varying temperature and moisture conditions

Satellite observations have shown greening trends in tundra in response to climate change, suggesting increases in productivity. To better understand the ability of remote sensing to detect climate impacts on tundra vegetation productivity, we applied a photosynthetic light use efficiency model to simulated climate change treatments of tundra vegetation. We examined changes in the Normalized Difference Vegetation Index (NDVI) and photosynthetic light use efficiency (ε) in experimental warming and moisture treatments designed to simulate climate change in northern Alaska. Plots were warmed either passively, using Open Top Chambers, or actively using electric heaters in the soil. In one set of plots water table depth was actively altered, while other plots were established in locations that were naturally wet or dry. Over two growing seasons, plot-level carbon flux and spectral reflectance measurements were collected, and the results were used to derive a light use efficiency model that could explore the effects of moisture and temperature treatments using remote sensing.Warming increased values of canopy greenness (NDVI) relative to control plots, this effect being more pronounced in wet plots than in dry plots. Light use efficiency (LUE), the relationship between absorbed photosynthetically active radiation (PAR) and gross ecosystem production (GEP), was consistent across warming treatments, growing season, subsequent years, and sites. However, LUE was affected by vegetation type, which varied with moisture; plots in naturally dry locations showed reduced light use efficiency relative to moist plots. Additionally moss exhibited reduced LUE relative to vascular plants. Understory moss production, not accounted for by the usual definition of the fraction of absorbed PAR (f_APAR), was found to be a significant part of total GEP, particularly in areas with low vascular plant cover. These results support the use of light use efficiency models driven by spectral reflectance for estimating GEP in tundra vegetation, provided effects of vegetation functional type (e.g. mosses versus vascular plants) and microtopography are considered.     

Reflectance assessment of seasonal and annual changes in biomass and CO2 uptake of a Mediterranean shrubland submitted to experimental warming and drought

We aimed to evaluate how the remote sensing vegetation indices NDVI and PRI responded to seasonal and annual changes in an early successional stage Mediterranean coastal shrubland canopy that was submitted to experimental warming and drought simulating predicted climate change for the next decades. These conditions were obtained by using a new non-intrusive methodological approach that increases the temperature and prolongs the drought period by using roofs that automatically cover the vegetation after the sunset or when it rains. On average, warming increased air temperature by 0.7 °C and soil temperature by 1.6 °C, and the drought treatment reduced soil moisture by 22%. We measured spectral reflectance at the canopy level and at the individual plant level seasonally during 4 years. Shrubland NDVI tracked the community development and activity. In control and warming treatments, NDVI increased with the years while it did not change in the drought treatment. There was a good relationship between NDVI and both community and individual plant biomass. NDVI also decreased in summer seasons when some species dry or decolour. The NDVI of E. multiflora plant individuals was lower in autumn and winter than in the other seasons, likely because of flowering. Shrubland PRI decreased only in winter, similarly to the PRI of the most dominant species, G. alypum. At this community scale, NDVI was better related than PRI to photosynthetic activity, probably because photosynthetic fluxes followed canopy seasonal greening in this complex canopy, which includes brevideciduous, annual and evergreen species and variable morphologies and canopy coverage. PRI followed the seasonal variations in photosynthetic rates in E. multiflora and detected the decreased photosynthetic rates of drought treatment. However, PRI did not track the photosynthetic rates of G. alypum plants which have lower LAIs than E. multiflora. In this community, which is in its early successional stages, NDVI was able to track biomass, and indirectly, CO₂ uptake changes, likely because LAI values did not saturate NDVI. Thus, NDVI appears as a valid tool for remote tracking of this community development. PRI was less adequate for photosynthetic assessment of this community especially for its lower LAI canopies. PRI usefulness was also species-dependent and could also be affected by flowering. These results will help to improve the interpretation of remote sensing information on the structure and physiological status of these Mediterranean shrublands, and to gain better insight on ecological and environmental controls on their ecosystem carbon dioxide exchange. They also show the possibility of assessing the impacts of climate change on shrubland communities.     

Utility of spectral vegetation index for estimation of gross CO2 flux under varied sky conditions

To estimate the gross CO₂ flux (F_CO2) of deciduous coniferous forest from canopy spectral reflectance, we introduced spectral vegetation indices (VIs) into a light use efficiency (LUE) model of mature Japanese larch (Larix kaempferi) forest. We measured the eddy covariance CO₂ flux and spectral reflectance of larch canopy at half-hourly intervals during one growing season, and investigated the relationships between the parameters of the LUE model (FAPAR, ?) and 3 types of VIs (NDVI, PRI, EVI) in both clear sky and cloudy conditions.FAPAR (fraction of absorbed photosynthetically active radiation) had a positive linear relationship with both NDVI (normalized difference vegetation index) and EVI (enhanced vegetation index), and the sky condition had little effect on the relationships. The relative RMSE (root mean square error) of the APAR (absorbed photosynthetically active radiation) based on the incoming PAR and estimated FAPAR from a linear function of NDVI was less than 10.5%, irrespective of sky condition.Half-hourly values of ? (conversion efficiency of absorbed energy) showed both seasonal variation related to leaf phenology and short-term variation related to light intensity due to varied sun position and sky condition. Both EVI and PRI (photochemical reflectance index) were significantly correlated with ?. EVI showed a positive linear relationship with ? as a result of their similar seasonal variation. However, since EVI did not detect short-term variation of ?, their relationship differed among sky conditions. On the other hand, although PRI could trace the short-term variation of ? in green needles, the relationship became non-linear due to drastic reduction of PRI in the senescent needles.EVI/(PRI/PRI_min), a combined index based on a 6-day moving minimum value of PRI (PRI_min), showed a linear relationship with half-hourly values of ? throughout the seasons irrespective of sky condition. This index allow us to estimate ? in all sky conditions with a smaller error (rRMSE = 35.2%) than using EVI or PRI alone (38.7%-48.7%). Consequently, this combined index-derived ? and NDVI-based FAPAR gave a low estimation error of F_CO2 (rRMSE = 36.4%, RMSE = 8.3 μmol m^− 2 s^− 1). Although there are still various issues to resolve, including adaptive limit and combination of vegetation index type, we conclude that the combination of PRI and EVI increased the accuracy of estimation of CO₂ uptake in deciduous forest even though sky conditions varied.     

Parallel adjustments in vegetation greenness and ecosystem CO2 exchange in response to drought in a Southern California chaparral ecosystem

Some form of the light use efficiency (LUE) model is used in most models of ecosystem carbon exchange based on remote sensing. The strong relationship between the normalized difference vegetation index (NDVI) and light absorbed by green vegetation make models based on LUE attractive in the remote sensing context. However, estimation of LUE has proven problematic since it varies with vegetation type and environmental conditions. Here we propose that LUE may in fact be correlated with vegetation greenness (measured either as NDVI at constant solar elevation angle, or a red edge chlorophyll index), making separate estimates of LUE unnecessary, at least for some vegetation types. To test this, we installed an automated tram system for measurement of spectral reflectance in the footprint of an eddy covariance flux system in the Southern California chaparral. This allowed us to match the spatial and temporal scales of the reflectance and flux measurements and thus to make direct comparisons over time scales ranging from minutes to years. The 3-year period of this study included both “normal” precipitation years and an extreme drought in 2002. In this sparse chaparral vegetation, diurnal and seasonal changes in solar angle resulted in large variation in NDVI independent of the actual quantity of green vegetation. In fact, one would come to entirely different conclusions about seasonal changes in vegetation greenness depending on whether NDVI at noon or NDVI at constant solar elevation angle were used. Although chaparral vegetation is generally considered “evergreen”, we found that the majority of the shrubs were actually semi-deciduous, leading to large seasonal changes in NDVI at constant solar elevation angle. LUE was correlated with both greenness indices at the seasonal timescale across all years. In contrast, the relationship between LUE and PRI was inconsistent. PRI was well correlated with LUE during the “normal” years but this relationship changed dramatically during the extreme drought. Contrary to expectations, none of the spectral reflectance indices showed consistent relationships with CO₂ flux or LUE over the diurnal time-course, possibly because of confounding effects of sun angle and stand structure on reflectance. These results suggest that greenness indices can be used to directly estimate CO₂ exchange at weekly timescales in this chaparral ecosystem, even in the face of changes in LUE. Greenness indices are unlikely to be as good predictors of CO₂ exchange in dense evergreen vegetation as they were in the sparse, semi-deciduous chaparral. However, since relatively few ecosystems are entirely evergreen at large spatial scales or over long time spans due to disturbance, these relationships need to be examined across a wider range of vegetation types.     

Mapping carbon and water vapor fluxes in a chaparral ecosystem using vegetation indices derived from AVIRIS

Using simple models derived from spectral reflectance, we mapped the patterns of ecosystem CO₂ and water fluxes in a semi-arid site in southern California during a period of extreme disturbance, marked by drought and fire. Employing a combination of low (∼ 2 km) and high (∼ 16 km) altitude images from the hyperspectral Airborne Visible Infrared Imaging Spectrometer (AVIRIS), acquired between April 2002 and September 2003, and ground data collected from an automated tram system, several vegetation indices were calculated for Sky Oaks field station, a FLUXNET and SpecNet site located in northern San Diego County (CA, USA). Based on the relationships observed between the fluxes measured by the eddy covariance tower and the vegetation indices, net CO₂ and water vapor flux maps were derived for the region around the flux tower. Despite differences in the scale of the images (from ∼ 2 m to 16 m pixel size) as well as marked differences in environmental conditions (drought in 2002, recovery in early 2003, and fire in mid 2003), net CO₂ and water flux modeled from AVIRIS-derived reflectance indices (NDVI, PRI and WBI) effectively tracked changes in tower fluxes across both drought and fire, and readily revealed spatial variation in fluxes within this landscape. After an initial period of net carbon uptake, drought and fire caused the ecosystem to lose carbon to the atmosphere during most of the study period. Our study shows the power of integrating optical and flux data in LUE models to better understand factors driving surface-atmosphere carbon and water vapor flux cycles, one of the main goals of SpecNet.