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.     

Towards estimation of canopy foliar biomass with spectral reflectance measurements

Canopy foliar biomass, defined as the product of leaf dry matter content and leaf area index, is an important measurement for global biogeochemical cycles. This study explores the potential for retrieving foliar biomass in green canopies using a spectral index, the Normalized Dry Matter Index (NDMI). This narrow-band index is based on absorption at the C-H bond stretch overtone and is correlated with leaf dry matter content in fresh green leaves. PROSPECT and SAIL model simulations suggest that the NDMI at the canopy scale is able to minimize the effects of leaf thickness and leaf water content and to maximize sensitivity to variation in canopy foliar biomass. The simulation outputs were analyzed with an ANOVA, and 87% of the variation in the NDMI is explained by leaf dry matter content. The NDMI was linearly related to foliar biomass (g cm^− 2) from model simulations (R² = 0.97). The NDMI calculated from spectral reflectances for one to four stacked leaves was also correlated with total leaf biomass (R² = 0.59). These results suggest that it may be possible to determine foliar biomass from airborne and satellite-borne imaging spectrometers, such as NASA's HyspIRI mission.     

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.     

A Landsat time series approach to characterize bark beetle and defoliator impacts on tree mortality and surface fuels in conifer forests

Insects are important forest disturbance agents, and mapping their effects on tree mortality and surface fuels represents a critical research challenge. Although various remote sensing approaches have been developed to monitor insect impacts, most studies have focused on single insect agents or single locations and have not related observed changes to ground-based measurements. This study presents a remote sensing framework to (1) characterize spectral trajectories associated with insect activity of varying duration and severity and (2) relate those trajectories to ground-based measurements of tree mortality and surface fuels in the Cascade Range, Oregon, USA. We leverage a Landsat time series change detection algorithm (LandTrendr), annual forest health aerial detection surveys (ADS), and field measurements to investigate two study landscapes broadly applicable to conifer forests and dominant insect agents of western North America. We distributed 38 plots across multiple forest types (ranging from mesic mixed-conifer to xeric lodgepole pine) and insect agents (defoliator [western spruce budworm] and bark beetle [mountain pine beetle]). Insect effects were evident in the Landsat time series as combinations of both short- and long-duration changes in the Normalized Burn Ratio spectral index. Western spruce budworm trajectories appeared to show a consistent temporal evolution of long-duration spectral decline (loss of vegetation) followed by recovery, whereas mountain pine beetle plots exhibited both short- and long-duration spectral declines and variable recovery rates. Although temporally variable, insect-affected stands generally conformed to four spectral trajectories: short-duration decline then recovery, short- then long-duration decline, long-duration decline, long-duration decline then recovery. When comparing remote sensing data with field measurements of insect impacts, we found that spectral changes were related to cover-based estimates (tree basal area mortality [R²_adj = 0.40, F_1,34 = 24.76, P < 0.0001] and down coarse woody detritus [R²_adj = 0.29, F_1,32 = 14.72, P = 0.0006]). In contrast, ADS changes were related to count-based estimates (e.g., ADS mortality from mountain pine beetle positively correlated with ground-based counts [R²_adj = 0.37, F_1,22 = 14.71, P = 0.0009]). Fine woody detritus and forest floor depth were not well correlated with Landsat- or aerial survey-based change metrics. By characterizing several distinct temporal manifestations of insect activity in conifer forests, this study demonstrates the utility of insect mapping methods that capture a wide range of spectral trajectories. This study also confirms the key role that satellite imagery can play in understanding the interactions among insects, fuels, and wildfire.     

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.     

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.     

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

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.     

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.     

Remote chlorophyll-a retrieval in turbid, productive estuaries: Chesapeake Bay case study

Accurate remote assessment of phytoplankton chlorophyll a (chla) concentration is particularly challenging in turbid, productive waters. Recently a conceptual model containing reflectance in three spectral bands in the red and near infra-red range of the spectrum was suggested for retrieving chla concentrations in turbid productive waters; it was calibrated and validated in lakes and reservoirs in Nebraska and Iowa. The objective of this paper is to evaluate the performance of this three band model as well as its special case, the two-band model to estimate chla concentration in Chesapeake Bay, as representative of estuarine Case II waters, and to assess the accuracy of chla retrieval. To evaluate the model performance, dual spectroradiometers were used to measure subsurface spectral radiance reflectance in the visible and near infra-red range of the spectrum. Water samples were collected concurrently and contained widely variable chla (9 to 77.4 mg/m³) and total suspended solids (7-65 mg/L dry wt). Colored dissolved organic matter (CDOM) absorption at 440 nm was 0.20 to 2.50 m^− 1; Secchi disk transparency ranged from 0.28 to 1.5 m. The two- and three-band models were spectrally tuned to select the spectral bands for most accurate chla estimation. Strong linear relationships were established between analytically measured chla and both the three-band model [R^− 1(675)-R^− 1(695)] × R(730) and the two-band model R(720)/R(670), where R(λ) is reflectance at wavelength λ. The three-band model accounted for 81% of variation in chla and allowed estimation of chla with a root mean square error (RMSE) of less than 7.9 mg/m³, whereas the two-band model accounted for 79% of chla variability and RMSE of chla estimation was below 8.4 mg/m³. The three-band model with MERIS spectral bands allows accurate chla estimation with RMSE below 9.1 mg/m³. Two-band model with SeaWiFS bands and MODIS 667 nm and 748 nm bands can estimate chla with RMSE below 11 mg/m³. The findings underlined the rationale behind the conceptual model and demonstrated the robustness of this algorithm for chla retrieval in turbid, productive estuarine waters.     

Enhancing a leaf radiative transfer model to estimate concentrations and in vivo specific absorption coefficients of total carotenoids and chlorophylls a and b from single-needle reflectance and transmittance

The relative concentrations of different pigments within a leaf have significant physiological and spectral consequences. Photosynthesis, light use efficiency, mass and energy exchange, and stress response are dependent on relationships among an ensemble of pigments. This ensemble also determines the visible characteristics of a leaf, which can be measured remotely and used to quantify leaf biochemistry and structure. But current remote sensing approaches are limited in their ability to resolve individual pigments. This paper focuses on the incorporation of three pigments—chlorophyll a, chlorophyll b, and total carotenoids—into the LIBERTY leaf radiative transfer model to better understand relationships between leaf biochemical, biophysical, and spectral properties.Pinus ponderosa and Pinus jeffreyi needles were collected from three sites in the California Sierra Nevada. Hemispheric single-leaf visible reflectance and transmittance and concentrations of chlorophylls a and b and total carotenoids of fresh needles were measured. These data were input to the enhanced LIBERTY model to estimate optical and biochemical properties of pine needles. The enhanced model successfully estimated reflectance (RMSE = 0.0255, BIAS = 0.00477, RMS%E = 16.7%), had variable success estimating transmittance (RMSE = 0.0442, BIAS = 0.0294, RMS%E = 181%), and generated very good estimates of carotenoid concentrations (RMSE = 2.48 µg/cm², BIAS = 0.143 µg/cm², RMS%E = 20.4%), good estimates of chlorophyll a concentrations (RMSE = 10.7 µg/cm², BIAS = − 0.992 µg/cm², RMS%E = 21.1%), and fair estimates of chlorophyll b concentrations (RMSE = 7.49 µg/cm², BIAS = − 2.12 µg/cm², RMS%E = 43.7%). Overall root mean squared errors of reflectance, transmittance, and pigment concentration estimates were lower for the three-pigment model than for the single-pigment model. The algorithm to estimate three in vivo specific absorption coefficients is robust, although estimated values are distorted by inconsistencies in model biophysics. The capacity to invert the model from single-leaf reflectance and transmittance was added to the model so it could be coupled with vegetation canopy models to estimate canopy biochemistry from remotely sensed data.     

Estimating chlorophyll a concentrations from remote-sensing reflectance in optically shallow waters

A multi-spectral classification and quantification technique is developed for estimating chlorophyll a concentrations, Chl, in shallow oceanic waters where light reflected by the bottom can contribute significantly to the above-water remote-sensing reflectance spectra, R_rs(λ). Classification criteria for determining bottom reflectance contributions for shipboard R_rs(λ) data from the west Florida shelf and Bahamian waters (1998-2001; n = 451) were established using the relationship between R_rs(412)/R_rs(670) and the spectral curvature about 555 nm, [R_rs(412) ? R_rs(670)]/R_rs(555)². Chlorophyll concentrations for data classified as “optically deep” and “optically shallow” were derived separately using best-fit cubic polynomial functions developed from the band-ratios R_rs(490)/R_rs(555) and R_rs(412)/R_rs(670), respectively. Concentrations for transitional data were calculated from weighted averages of the two derived values. The root-mean-square error (RMSE_log10) calculated for the entire data set using the new technique was 14% lower than the lowest error derived using the best individual band-ratio. The standard blue-to-green, band-ratio algorithm yields a 26% higher RMSE_log10 than that calculated using the new method. This study demonstrates the potential of quantifying chlorophyll a concentrations more accurately from multi-spectral satellite ocean color data in oceanic regions containing optically shallow waters.     

Assessing the potential of SeaWiFS and MODIS for estimating chlorophyll concentration in turbid productive waters using red and near-infrared bands

Bio-optical algorithms for remote estimation of chlorophyll-a concentration (Chl) in case-1 waters exploit the upwelling radiation in the blue and green spectral regions. In turbid productive waters other constituents, that vary independently of Chl, absorb and scatter light in these spectral regions. As a consequence, the accurate estimation of Chl in turbid productive waters has so far not been feasible from satellite sensors. The main purpose of this study was to evaluate the extent to which near-infrared (NIR) to red reflectance ratios could be applied to the Sea Wide Field-of-View Sensor (SeaWiFS) and the Moderate Imaging Spectrometer (MODIS) to estimate Chl in productive turbid waters. To achieve this objective, remote-sensing reflectance spectra and relevant water constituents were collected in 251 stations over lakes and reservoirs with a wide variability in optical parameters (i.e. 4 ≤ Chl ≤ 240 mg m^− 3; 18 ≤ Secchi disk depth ≤ 308 cm). SeaWiFS and MODIS NIR and red reflectances were simulated by using the in-situ hyperspectral data. The proposed algorithms predicted Chl with a relative random uncertainty of approximately 28% (average bias between − 1% and − 4%). The effects of reflectance uncertainties on the predicted Chl were also analyzed. It was found that, for realistic ranges of R_rs uncertainties, Chl could be estimated with a precision better than 40% and an accuracy better than ± 35%. These findings imply that, provided that an atmospheric correction scheme specific for the red-NIR spectral region is available, the extensive database of SeaWiFS and MODIS images could be used to quantitatively monitor Chl in turbid productive waters.     

Possibility investigation of tree diversity mapping using Landsat ETM+ data in the Hyrcanian forests of Iran

Lack of data often limits understanding and management of biodiversity in forested areas. Remote sensing imagery has considerable potential to aid in the monitoring and prediction of biodiversity across many spatial and temporal scales. In this paper, we explored the possibility of defining relationships between species diversity indices and Landsat ETM+ reflectance values for Hyrcanian forests in Golestan province of Iran. We used the COST model for atmospheric correction of the imagery. Linear regression models were implemented to predict measures of biodiversity (species richness and reciprocal of Simpson indices) using various combinations of Landsat spectral data. Species richness was modeled using the band set ETM5, ETM7, DVI, wetness and variances of ETM1, ETM2 and ETM5 (adjusted R² = 0.59, RMSE = 1.51). Reciprocal of Simpson index was modeled using the band set NDVI, brightness, greenness, variances of ETM2, ETM5 and ETM7 (adjusted R² = 0.459 RMSE = 1.15). The results demonstrated that spectral reflectance from Landsat can be used to effectively model tree species diversity. Predictive map derived from the presented methodology can help evaluate spatial aspects and monitor tree species diversity of the studied forest. The methodology also facilitates the evaluation of forest management and conservation strategies in northern Iran.     

Testing the potential of multi-spectral remote sensing for retrospectively estimating fire severity in African Savannahs

The remote sensing of fire severity is a noted goal in studies of forest and grassland wildfires. Experiments were conducted to discover and evaluate potential relationships between the characteristics of African savannah fires and post-fire surface spectral reflectance in the visible to shortwave infrared spectral region. Nine instrumented experimental fires were conducted in semi-arid woodland savannah of Chobe National Park (Botswana), where fire temperature (T_max) and duration (dt) were recorded using thermocouples positioned at different heights and locations. These variables, along with measures of fireline intensity (FLI), integrated temperature with time (T_sum) and biomass (and carbon/nitrogen) volatilised were compared to post-fire surface spectral reflectance. Statistically significant relationships were observed between (i) the fireline intensity and total nitrogen volatilised (r² = 0.54, n = 36, p < 0.001), (ii) integrated temperature (T_sum−μ) and total biomass combusted (r² = 0.72, n = 32, p < 0.001), and (iii) fire duration as measured at the top-of-grass sward thermocouple (dt_T) and total biomass combusted (r² = 0.74, n = 34, p < 0.001) and total nitrogen volatilised (r² = 0.73, n = 34, p < 0.001). The post-fire surface spectral reflectance was found to be related to dt and T_sum via a quadratic relationship that varied with wavelength. The use of visible and shortwave infrared band ratios produced statistically significant linear relationships with fire duration as measured by the top thermocouple (dt_T) (r² = 0.76, n = 34, p < 0.001) and the mean of T_sum (r² = 0.82, n = 34, p < 0.001). The results identify fire duration as a versatile measure that relates directly to the fire severity, and also illustrate the potential of spectrally-based fire severity measures. However, the results also point to difficulties when applying such spectrally-based techniques to Earth Observation satellite imagery, due to the small-scale variability noted on the ground. Results also indicate the potential for surface spectral reflectance to increase following higher severity fires, due to the laying down of high albedo white mineral ash. Most current techniques for mapping burned area rely on the general assumption that surface albedo decreases following a fire, and so if the image spatial resolution was high enough such methods may fail. Determination of the effect of spatial resolution on a sensor's ability to detect white ash was investigated using a validated optical mixture modelling approach. The most appropriate mixing model to use (linear or non-linear) was assessed using laboratory experiments. A linear mixing model was shown most appropriate, with results suggesting that sensors having spatial resolutions significantly higher than those of Landsat ETM+ will be required if patches of white ash are to be used to provide EO-derived information on the spatial variation of fire severity.     

A simple band ratio technique to quantify the colored dissolved and detrital organic material from ocean color remotely sensed data

Spectral reflectances of the ocean, R, as derived from ocean color remote sensing data at four wavelengths (412, 443, 490, and 555 nm), can be used to form two ratios of spectral reflectance, namely R(412)/R(443), and R(490)/R(555), thereafter denoted R₄₄₃⁴¹² and R₅₅₅⁴⁹⁰. The former is mainly sensitive to the colored dissolved organic material (CDOM), albeit influenced by the algal content as depicted by the chlorophyll concentration, ([Chl]); in contrast, the latter is essentially depending on [Chl], although it is also influenced by CDOM. Therefore the signatures of CDOM and [Chl] which are not truly separable, can nevertheless be identified by considering simultaneously the two ratios. The concomitant variations in these ratios can be established via a bio-optical model developed for Case 1 waters. This model implicitly includes a “mean” relationship between CDOM and [Chl], and thus produces a unique curve relating R₄₄₃⁴¹² to R₅₅₅⁴⁹⁰. Deviations with respect to this mean relationship can be introduced through a factor Φ, with Φ > 1 (excess) or < 1 (deficit), applied to the CDOM-[Chl] ratio. A family of R₄₄₃⁴¹²-R₅₅₅⁴⁹⁰ curves is thus generated, in correspondence with the discrete values given to Φ; this “grid” (or numerically, a 2-D lookup table) allows the Φ-[Chl] couple to be unambiguously derived for any R₄₄₃⁴¹²-R₅₅₅⁴⁹⁰ couple. By applying this straightforward algorithm to actual reflectance ratios derived from ocean color imagery, the relative anomalies in CDOM with respect to its standard (Chl-related) values can be efficiently assessed. Within the global ocean (discarding the coastal zones), the Φ factor is widely varying, between at least ? and 3, and is roughly log-normally distributed around ~ 1 (no anomaly). The spatial distributions of the Φ factor in the whole ocean are strongly featured according to latitude, season, and hydrographic regimes, and these features are regularly reproducible, from year to year (2002-2007). This simple method is also validated against available in situ data, and its results compare favorably, for instance, to those of the GSM (Garver-Siegel-Maritorena) inversion method, in terms of retrieved CDOM concentrations and distribution patterns.     

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.