Mapping of continuous floristic gradients in grasslands using hyperspectral imagery

Transitions between plant species assemblages are often continuous with the form of the transition dependent on the ‘slope’ of environmental gradients and on the style of self-organization in vegetation. Image segmentation can present misleading or even erroneous results if applied to continuous spatial changes in vegetation. Even methods that allow for multiple-class memberships of pixels presuppose the existence of ideal types of species assemblages that constitute mixtures—an assumption that does not fit the case of continua where any section of a gradient is as ‘pure’ as any other section like in modulations of grassland species composition.Thus, we attempted to spatially model floristic gradients in Bavarian meadows by extrapolating axes of an unconstrained ordination of species data. The models were based on high-resolution hyperspectral airborne imagery. We further modelled the distribution of plant functional response types (Ellenberg indicator values) and the cover values of selected species. The models were made with partial least squares (PLS) regression analyses. The realistic utility of the regression models was evaluated by full leave-one-out cross-validation.The modelled floristic gradients showed a considerable agreement with ground-based observations of floristic gradients (R²=0.71 and 0.66 for the first two axes of ordination). Apart from mapping the most important continuous floristic differences, we mapped gradients in the appearance of plant functional response groups as represented by averaged Ellenberg indicator values for soil pH (R²=0.76), water supply (R²=0.66) and nutrient supply (R²=0.75), while models for the cover of single species were weak.Compared to many other vegetation attributes, plant species composition is difficult to detect with remote sensing techniques. This is partly caused by a lack of compatibility between methods of vegetation ecology and remote sensing. We believe that the present study has the potential to increase compatibility as neither spectral nor vegetation information gets lost by a classifying step.     

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

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

Quantitative analysis of salt-affected soil reflectance spectra: A comparison of two adaptive methods (PLSR and ANN)

In this paper the possibility of predicting salt concentrations in soils from measured reflectance spectra is studied using partial least squares regression (PLSR) and artificial neural network (ANN). Performance of these two adaptive methods has been compared in order to examine linear and non-linear relationship between soil reflectance and salt concentration.Experiment-, field- and image-scale data sets were prepared consisting of soil EC measurements (dependent variable) and their corresponding reflectance spectra (independent variables). For each data set, PLSR and ANN predictive models of soil salinity were developed based on soil reflectance data. The predictive accuracies of PLSR and ANN models were assessed against independent validation data sets not included in the calibration or training phase.The results of PLSR analyses suggest that an accurate to good prediction of EC can be made based on models developed from experiment-scale data (R² > 0.81 and RPD (ratio of prediction to deviation) > 2.1) for soil samples salinized by bischofite and epsomite minerals. For field-scale data sets, the PLSR predictive models provided approximate quantitative EC estimations (R² = 0.8 and RPD = 2.2) for grids 1 and 6 and poor estimations for grids 2, 3, 4 and 5. The salinity predictions from image-scale data sets by PLSR models were very reliable to good (R² between 0.86 and 0.94 and RPD values between 2.6 and 4.1) except for sub-image 2 (R² = 0.61 and RPD = 1.2).The ANN models from experiment-scale data set revealed similar network performances for training, validation and test data sets indicating a good network generalization for samples salinized by bischofite and epsomite minerals. The RPD and the R² between reference measurements and ANN outputs of theses models suggest an accurate to good prediction of soil salinity (R² > 0.92 and RPD > 2.3). For the field-scale data set, prediction accuracy is relatively poor (0.69 > R² > 0.42). The ANN predictive models estimating soil salinity from image-scale data sets indicate a good prediction (R² > 0.86 and RPD > 2.5) except for sub-image 2 (R² = 0.6 and RPD = 1.2).The results of this study show that both methods have a great potential for estimating and mapping soil salinity. Performance indexes from both methods suggest large similarity between the two approaches with PLSR advantages. This indicates that the relation between soil salinity and soil reflectance can be approximated by a linear function.     

Mapping burn severity and burning efficiency in California using simulation models and Landsat imagery

Uncertainties in burning efficiency (BE) estimates can lead to large errors in fire emission quantification (from 23% to 46%). One of the main causes of these errors is the spatial variability of fuel consumption within burned areas. This paper studies whether burn severity (BS) maps can be used to improve BE assessment. A burn severity map of two large fires in California was obtained by inverting a simulation model constrained by post-fire observations from Landsat TM imagery. Model output values of BS were validated against field measurements, obtaining a high correlation (R² = 0.85) and low errors (Root Mean Square Error, RMSE = 0.14) throughout a wide range of BS levels. The BS map obtained was then used to adjust BE reference values per vegetation type found in the area before the fire. The adjusted burning efficiency (BEadj) was compared to the burned biomass, which was estimated by subtracting vegetation indices from pre- and post-fire images. Results showed a high correlation for conifers (R² = 0.75) and hardwoods (R² = 0.73), and a moderate correlation (R² ∼ 0.5) for shrubs and grasslands. In general, for all vegetation types BEadj performed better (R² = 0.4-0.75) than literature-based BE (R² < 0.0001). This study demonstrates: (i) the consistency of the simulation model inversion for BS estimation in temperate ecosystems, and (ii) the improvement of BE estimation when the spatial variability of the combustion was quantified in terms of BS.     

Composition versus physiognomy of vegetation as predictors of bird assemblages: The role of lidar

Whether diversity and composition of avian communities is determined primarily by responses of species to the floristic composition or to the structural characteristics of habitats has been an ongoing debate, at least since the publication of MacArthur and MacArthur (1961). This debate, however, has been hampered by two problems: 1) it is notoriously time consuming to measure the physiognomy of habitat, particularly in forests, and 2) rigorous statistical methods to predict the composition of bird assemblages from assemblages of plants have not been available. Here we use airborne laser scanning (lidar) to measure the habitat (vegetation) structure of a montane forest across large spatial extents with a very fine grain. Furthermore, we use predictive co-correspondence and canonical correspondence analyses to predict the composition of bird communities from the composition and structure of another community (i.e. plants). By using these new techniques, we show that the physiognomy of the vegetation is a significantly more powerful predictor of the composition of bird assemblages than plant species composition in the field and as well in the shrub/tree layer, both on a level of p < 0.001. Our results demonstrate that ecologists should consider remote sensing as a tool to improve the understanding of the variation of bird assemblages in space and time. Particularly in complex habitats, such as forests, lidar is a valuable and comparatively inexpensive tool to characterize the structure of the canopy even across large and rough terrain.     

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.     

Mapping an invasive plant, Phragmites australis, in coastal wetlands using the EO-1 Hyperion hyperspectral sensor

Mapping tools are needed to document the location and extent of Phragmites australis, a tall grass that invades coastal marshes throughout North America, displacing native plant species and degrading wetland habitat. Mapping Phragmites is particularly challenging in the freshwater Great Lakes coastal wetlands due to dynamic lake levels and vegetation diversity. We tested the applicability of Hyperion hyperspectral satellite imagery for mapping Phragmites in wetlands of the west coast of Green Bay in Wisconsin, U.S.A. A reference spectrum created using Hyperion data from several pure Phragmites stands within the image was used with a Spectral Correlation Mapper (SCM) algorithm to create a raster map with values ranging from 0 to 1, where 0 represented the greatest similarity between the reference spectrum and the image spectrum and 1 the least similarity. The final two-class thematic classification predicted monodominant Phragmites covering 3.4% of the study area. Most of this was concentrated in long linear features parallel to the Green Bay shoreline, particularly in areas that had been under water only six years earlier when lake levels were 66 cm higher. An error matrix using spring 2005 field validation points (n = 129) showed good overall accuracy—81.4%. The small size and linear arrangement of Phragmites stands was less than optimal relative to the sensor resolution, and Hyperion's 30 m resolution captured few if any pure pixels. Contemporary Phragmites maps prepared with Hyperion imagery would provide wetland managers with a tool that they currently lack, which could aid attempts to stem the spread of this invasive species.     

Modeling forest songbird species richness using LiDAR-derived measures of forest structure

Conservation of biodiversity requires information at many spatial scales in order to detect and preserve habitat for many species, often simultaneously. Vegetation structure information is particularly important for avian habitat models and has largely been unavailable for large areas at the desired resolution. Airborne LiDAR, with its combination of relatively broad coverage and fine resolution provides existing new opportunities to map vegetation structure and hence avian habitat. Our goal was to model the richness of forest songbirds using forest structure information obtained from LiDAR data. In deciduous forests of southern Wisconsin, USA, we used discrete-return airborne LiDAR to derive forest structure metrics related to the height and density of vegetation returns, as well as composite variables that captured major forest structural elements. We conducted point counts to determine total forest songbird richness and the richness of foraging, nesting, and forest edge-related habitat guilds. A suite of 35 LiDAR variables were used to model bird species richness using best-subsets regression and we used hierarchical partitioning analysis to quantify the explanatory power of each variable in the multivariate models. Songbird species richness was correlated most strongly with LiDAR variables related to canopy and midstory height and midstory density (R² = 0.204, p < 0.001). Richness of species that nest in the midstory was best explained by canopy height variables (R² = 0.197, p < 0.001). Species that forage on the ground responded to mean canopy height and the height of the lower canopy (R² = 0.149, p < 0.005) while aerial foragers had higher richness where the canopy was tall and dense and the midstory more sparse (R² = 0.216, p < 0.001). Richness of edge-preferring species was greater where there were fewer vegetation returns but higher density in the understory (R² = 0.153, p < 0.005). Forest interior specialists responded positively to a tall canopy, developed midstory, and a higher proportion of vegetation returns (R² = 0.195, p < 0.001). LiDAR forest structure metrics explained between 15 and 20% of the variability in richness within deciduous forest songbird communities. This variability was associated with vertical structure alone and shows how LiDAR can provide a source of complementary predictive data that can be incorporated in models of wildlife habitat associations across broad geographical extents.     

Mapping of riparian zone attributes using discrete return LiDAR, QuickBird and SPOT-5 imagery: Assessing accuracy and costs

Riparian zones in Australia are exposed to increasing pressures because of disturbance from agricultural and urban expansion, weed invasion, and overgrazing. Accurate and cost-effective mapping of riparian environments is important for assessing riparian zone functions associated with water quality, biodiversity, and wildlife habitats. The objective of this research was to compare the accuracy and costs of mapping riparian zone attributes from image data acquired by three different sensor types, i.e. Light Detection and Ranging (LiDAR) (0.5-2.4 m pixels), and multi-spectral QuickBird (2.4 m pixels) and SPOT-5 (10 m pixels). These attributes included streambed width, riparian zone width, plant projective cover, longitudinal continuity, vegetation overhang, and bank stability. The riparian zone attributes were mapped for a study area along Mimosa Creek in the Fitzroy Catchment, Central Queensland, Australia. Object-based image and regression analyses were used for mapping the riparian zone attributes. The validation of the LiDAR, QuickBird, and SPOT-5 derived maps of streambed width (R = 0.99, 0.71, and 0.44 respectively) and riparian zone width (R = 0.91, 0.87, and 0.74 respectively) against field derived measurements produced the highest accuracies for the LiDAR data and the lowest using the SPOT-5 image data. Cross-validation estimates of misclassification produced a root mean square error of 1.06, 1.35 and 1.51 from an ordinal scale from 0 to 4 of the bank stability score for the LiDAR, QuickBird and SPOT-5 image data, respectively. The validation and empirical modelling showed high correlations for all datasets for mapping plant projective cover (R > 0.93). The SPOT-5 image data were unsuitable for assessment of riparian zone attributes at the spatial scale of Mimosa Creek and associated riparian zones. Cost estimates of image and field data acquisition and processing of the LiDAR, QuickBird, and SPOT-5 image data showed that discrete return LiDAR can be used for costs lower than those for QuickBird image data over large spatial extents (e.g. 26,000 km of streams). With the higher level of vegetation structural and landform information, mapping accuracies, geometric precision, and lower overall costs at large spatial extents, LiDAR data are a feasible means for assessment of riparian zone attributes.     

Monitoring elevation variations in leaf phenology of deciduous broadleaf forests from SPOT/VEGETATION time-series

In mountain forest ecosystems where elevation gradients are prominent, temperature gradient-based phenological variability can be high. However, there are few studies that assess the capability of remote sensing observations to monitor ecosystem phenology along elevation gradients, despite their relevance under climate change. We investigated the potential of medium resolution remotely sensed data to monitor the elevation variations in the seasonal dynamics of a temperate deciduous broadleaf forested ecosystem. Further, we explored the impact of elevation on the onset of spring leafing. This study was based on the analysis of multi-annual time-series of VEGETATION data acquired over the French Pyrenees Mountain Region (FPMR), in conjunction with simultaneous ground-based observations of leaf phenology made for two dominant tree species in the region (oak and beech). The seasonal variations in the perpendicular vegetation index (PVI) were analyzed during a five-year period (2002 to 2006). The five years of data were averaged into a one sole year in order to fill the numerous large spatio-temporal gaps due to cloud and snow presence - frequent in mountains - without altering the temporal resolution. Since a VEGETATION pixel (1 km²) includes several types of land cover, the broadleaf forest-specific seasonal dynamics of PVI was reconstructed pixel-by-pixel using a temporal unmixing method based on a non-parametric statistical approach. The spatial pattern of the seasonal response of PVI was clearly consistent with the relief. Nevertheless the elevational or geographic range of tree species, which differ in their phenology sensitivity to temperature, also has a significant impact on this pattern. The reduction in the growing season length with elevation was clearly observable from the delay in the increase of PVI in spring and from the advance of its decrease in the fall. The elevation variations in leaf flushing timing were estimated from the temporal change in PVI in spring over the study area. They were found to be consistent with those measured in situ (R2 > 0.95). It was deduced that, over FPMR, the mean delay of leaf flushing timing for every 100 m increase in elevation was estimated be approximately 2.3 days. The expected estimation error of satellite-based leaf unfolding date for a given elevation was approximately 2 days. This accuracy can be considered as satisfactory since it would allow us to detect changes in leafing timing of deciduous broadleaf forests with a magnitude equivalent to that due to an elevation variation of 100 m (2.3 days on average), or in other words, to that caused by a variation in the mean annual air temperature of 0.5 °C. Although averaging the VEGETATION data over five years led to a loss of interannual information, it was found to be a robust approach to characterise the elevation variations in spring leafing and its long-term trends.     

Remote sensing of evapotranspiration and carbon uptake at Harvard Forest

A land surface vegetation index, defined as the difference of microwave land surface emissivity at 19 and 37 GHz, was calculated for a heavily forested area in north central Massachusetts. The microwave emissivity difference vegetation index (EDVI) was estimated from satellite SSM/I measurements at the defined wavelengths and used to estimate land surface turbulent fluxes. Narrowband visible and infrared measurements and broadband solar radiation observations were used in the EDVI retrievals and turbulent flux estimations. The EDVI values represent physical properties of crown vegetation such as vegetation water content of crown canopies. The collocated land surface turbulent and radiative fluxes were empirically linked together by the EDVI values. The EDVI values are statistically sensitive to evapotranspiration fractions (EF) with a correlation coefficient (R) greater than 0.79 under all-sky conditions. For clear skies, EDVI estimates exhibit a stronger relationship with EF than normalized difference vegetation index (NDVI). Furthermore, the products of EDVI and input energy (solar and photosynthetically active radiation) are statistically significantly correlated to evapotranspiration (R = 0.95) and CO₂ uptake flux (R = 0.74), respectively.     

Estimating forest canopy height and terrain relief from GLAS waveform metrics

Quantifying aboveground biomass in forest ecosystems is required for carbon stock estimation, aspects of forest management, and further developing a capacity for monitoring carbon stocks over time. Airborne Light Detection And Ranging (LiDAR) systems, of all remote sensing technologies, have been demonstrated to yield the most accurate estimates of aboveground biomass for forested areas over a wide range of biomass values. However, these systems are limited by considerations including large data volumes and high costs. Within the constraints imposed by the nature of the satellite mission, the GeoScience Laser Altimeter System (GLAS) aboard ICESat has provided data conferring information regarding forest vertical structure for large areas at a low end user cost. GLAS data have been demonstrated to accurately estimate forest height and aboveground biomass especially well in topographically smooth areas with homogeneous forested conditions. However in areas with dense forests, high relief, or heterogeneous vegetation cover, GLAS waveforms are more complex and difficult to consistently characterize. We use airborne discrete return LiDAR data to simulate GLAS waveforms and to subsequently deconstruct coregistered GLAS waveforms into vegetation and ground returns. A series of waveform metrics was calculated and compared to topography and vegetation information gleaned from the airborne data. A model to estimate maximum relief directly from waveform metrics was developed with an R² of 0.76 (n = 110), and used for the classification of the maximum relief of the areas sensed by GLAS. Discriminant analysis was also conducted as an alternative classification technique. A model was also developed estimating forest canopy height from waveform metrics for all of the data (R² = 0.81, n = 110) and for the three separate relief classes; maximum relief 0-7 m (R² = 0.83, n = 44), maximum relief 7-15 m (R² = 0.88, n = 41) and maximum relief > 15 m (R² = 0.75, n = 25). The moderate relief class model yielded better predictions of forest height than the low relief class model which is attributed to the increasing variability of waveform metrics with terrain relief. The moderate relief class model also yielded better predictions than the high relief class model because of the mixing of vegetation and terrain signals in waveforms from high relief footprints. This research demonstrates that terrain can be accurately modeled directly from GLAS waveforms enabling the inclusion of terrain relief, on a waveform specific basis, as supplemental model input to improve estimates of canopy height.     

Estimating northern peatland CO2 exchange from MODIS time series data

Studies using satellite sensor-derived data as input to models for CO₂ exchange show promising results for closed forest stands. There is a need for extending this approach to other land cover types, in order to carry out large-scale monitoring of CO₂ exchange. In this study, three years of eddy covariance data from two peatlands in Sweden were averaged for 16-day composite periods and related to data from the Moderate Resolution Imaging Spectroradiometer (MODIS) and modeled photosynthetic photon flux density (PPFD). Noise in the time series of MODIS 250 m vegetation indices was reduced by using double logistic curve fits. Smoothed normalized difference vegetation index (NDVI) showed saturation during summertime, and the enhanced vegetation index (EVI) generally gave better results in explaining gross primary productivity (GPP). The strong linear relationships found between GPP and the product of EVI and modeled PPFD (R² = 0.85 and 0.76) were only slightly stronger than for the product of EVI and MODIS daytime 1 km land surface temperature (LST) (R² = 0.84 and 0.71). One probable reason for these results is that several controls on GPP were related to both modeled PPFD and daytime LST. Since ecosystem respiration (ER) was largely explained by diurnal LST in exponential relationships (R² = 0.89 and 0.83), net ecosystem exchange (NEE) was directly related to diurnal LST in combination with the product of EVI and modeled PPFD in multiple exponential regressions (R² = 0.81 and 0.73). Even though the R² values were somewhat weaker for NEE, compared to GPP and ER, the RMSE values were much lower than if NEE would have been estimated as the sum of GPP and ER. The overall conclusion of this study is that regression models driven by satellite sensor-derived data and modeled PPFD can be used to estimate CO₂ fluxes in peatlands.     

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.     

Capturing tree crown formation through implicit surface reconstruction using airborne lidar data

Forest structure data derived from lidar is being used in forest science and management for inventory analysis, biomass estimation, and wildlife habitat analysis. Regression analysis dominated previous approaches to the derivation of tree stem and crown parameters from lidar. The regression model for tree parameters is locally applied based on vertical lidar point density, the tree species involved, and stand structure in the specific research area. The results of this approach, therefore, are location-specific, limiting its applicability to other areas. For a more widely applicable approach to derive tree parameters, we developed an innovative method called ‘wrapped surface reconstruction’ that employs radial basis functions and an isosurface. Utilizing computer graphics, we capture the exact shape of an irregular tree crown of various tree species based on the lidar point cloud and visualize their exact crown formation in three-dimensional space. To validate the tree parameters given by our wrapped surface approach, survey-grade equipment (a total station) was used to measure the crown shape. Four vantage points were established for each of 55 trees to capture whole-tree crown profiles georeferenced with post-processed differential GPS points. The observed tree profiles were linearly interpolated to estimate crown volume. These fieldwork-generated profiles were compared with the wrapped surface to assess goodness of fit. For coniferous trees, the following tree crown parameters derived by the wrapped surface method were highly correlated (p < 0.05) with the total station-derived measurements: tree height (R² = 0.95), crown width (R² = 0.80), live crown base (R² = 0.92), height of the lowest branch (R² = 0.72), and crown volume (R² = 0.84). For deciduous trees, wrapped surface-derived parameters of tree height (R² = 0.96), crown width (R² = 0.75), live crown base (R² = 0.53), height of the lowest branch (R² = 0.51), and crown volume (R² = 0.89) were correlated with the total station-derived measurements. The wrapped surface technique is less susceptible to errors in estimation of tree parameters because of exact interpolation using the radial basis functions. The effect of diminished energy return causes the low correlation for lowest branches in deciduous trees (R² = 0.51), even though leaf-off lidar data was used. The wrapped surface provides fast and automated detection of micro-scale tree parameters for specific applications in areas such as tree physiology, fire modeling, and forest inventory.     

Monitoring herbaceous biomass and water content with SPOT VEGETATION time-series to improve fire risk assessment in savanna ecosystems

This paper evaluated the capacity of SPOT VEGETATION time-series to monitor the vegetation biomass and water content in order to improve fire risk assessment in the savanna ecosystem of Kruger National Park in South Africa. First, the single date and integrated vegetation index approach, which quantify the amount of herbaceous biomass at the end of the rain season, were evaluated using in situ biomass data. It was shown that the integral of the Ratio Vegetation Index (iRVI) during the rain season was the most suitable index to estimate herbaceous biomass (R² = 0.69). Next, the performance of single, greenness, and accumulated remotely sensed fire risk indices, related to vegetation water content, were evaluated using fire activity data. The Accumulated Relative Normalised Difference Vegetation Index Decrement (ARND) performed the best when estimating fire risk (c-index = 0.76). Finally, results confirmed that the assessment of fire risk was improved by combination of both the vegetation biomass (iRVI) and vegetation water content (ARND) related indices (c-index = 0.80). The monitoring of vegetation biomass and water content with SPOT VEGETATION time-series provided a more suitable tool for fire management and suppression compared to satellite-based fire risk assessment methods, only related to vegetation water content.     

Remote sensing of structural complexity indices for habitat and species distribution modeling

Spatial distribution models are increasingly used in ecological studies, but are limited by the poor accuracy of remote sensing (RS) for mapping microhabitat (< 0.1 ha) features. Mapping accuracy can be improved by combining advanced RS image-processing techniques with microhabitat data expressed as a structural complexity index (SCI). To test this idea, we used principal components analysis (PCA) and an additive SCI method developed for forest ecology (calculated by re-scaling and summing representative structural variables) to summarize 13 microhabitat-scale (0.04 ha) vegetation structure attributes describing the rare mountain bongo antelope's (Tragelaphus eurycerus isaaci) habitat in Kenya's Aberdare mountains. Microhabitat data were collected in 127 plots: 37 related to bongo habitat use, 90 from 1 km-spaced grid points representing overall habitat availability and bongo non-presence. We then assessed each SCI's effectiveness for discerning microhabitat variability and bongo habitat selection, using Wilcoxon Rank Sum tests for differences in mean SCI scores among plots divided into 4 vegetation classes, and the Area Under the Curve (AUC) of receiver operating characteristics from logistic regressions. We also examined the accuracy of predicted SCI scores resulting from regression models based on variables derived from a) ASTER imagery processed with spectral mixture and texture analysis, b) an SRTM DEM and c) rainfall data, using the 90 grid plots for model training and the bongo plots as an independent test dataset. Of the five SCIs derived, two performed best: the PCA-derived Canopy Structure Index (CSI) and an additive index summarizing 8 structural variables (AI₈). CSI and AI₈ showed significant differences between 5 of 6 vegetation class pairs, strong abilities to distinguish bongo-selected from available habitat (AUCs = 0.71 (CSI); 0.70 (AI₈)), and predicted scores 60-110% more accurate than reported by other studies using RS to quantify individual microhabitat structural attributes (CSI model R² = 0.51, RMSE = 0.19 (training) and 0.21 (test); AI₈ model R² = 0.46, RMSE = 0.17 (training) and 0.19 (test)). Repeating the Wilcoxon tests and logistic regressions with RS-predicted SCI values showed that AI₈ most effectively preserved the patterns found with the observed SCIs. These results demonstrate that SCIs effectively characterize microhabitat structure and selection, and boost microhabitat mapping accuracy when combined with enhanced RS image-processing techniques. This approach can improve distribution models and broaden their applicability, makes RS more relevant to applied ecology, and shows that processing field data to be more compatible with RS can improve RS-based habitat mapping accuracy.     

Mapping biomass and stress in the Sierra Nevada using lidar and hyperspectral data fusion

In this paper, we explored fusion of structural metrics from the Laser Vegetation Imaging Sensor (LVIS) and spectral characteristics from the Airborne Visible Infrared Imaging Spectrometer (AVIRIS) for biomass estimation in the Sierra Nevada. In addition, we combined the two sensors to map species-specific biomass and stress at landscape scale. Multiple endmember spectral mixture analysis (MESMA) was used to classify vegetation from AVIRIS images and obtain sub-pixel fractions of green vegetation, non-photosynthetic vegetation, soil, and shade. LVIS metrics, AVIRIS spectral indices, and MESMA fractions were compared with field measures of biomass using linear and stepwise regressions at stand (1 ha) level. AVIRIS metrics such as water band indices and shade fractions showed strong correlation with LVIS canopy height (r² = 0.69, RMSE = 5.2 m) and explained around 60% variability in biomass. LVIS variables were found to be consistently good predictors of total and species specific biomass (r² = 0.77, RMSE = 70.12 Mg/ha). Prediction by LVIS after species stratification of field data reduced errors by 12% (r² = 0.84, RMSE = 58.78 Mg/ha) over using LVIS metrics alone. Species-specific biomass maps and associated errors created from fusion were different from those produced without fusion, particularly for hardwoods and pines, although mean biomass differences between the two techniques were not statistically significant. A combined analysis of spatial maps from LVIS and AVIRIS showed increased water and chlorophyll stress in several high biomass stands in the study area. This study provides further evidence that lidar is better suited for biomass estimation, per se, while the best use of hyperspectral data may be to refine biomass predictions through a priori species stratification, while also providing information on canopy state, such as stress. Together, the two sensors have many potential applications in carbon dynamics, ecological and habitat studies.     

MODIS tree cover validation for the circumpolar taiga-tundra transition zone

A validation of the 2005 500 m MODIS vegetation continuous fields (VCF) tree cover product in the circumpolar taiga-tundra ecotone was performed using high resolution Quickbird imagery. Assessing the VCF's performance near the northern limits of the boreal forest can help quantify the accuracy of the product within this vegetation transition area. The circumpolar region was divided into 7 longitudinal zones and validation sites were selected in areas of varying tree cover where Quickbird imagery is available in Google Earth. Each site was linked to the corresponding VCF pixel and overlaid with a regular dot grid within the VCF pixel's boundary to estimate percent tree crown cover in the area. Percent tree crown cover was estimated using Quickbird imagery for 396 sites throughout the circumpolar region and related to the VCF's estimates of canopy cover for 2000-2005. Regression results of VCF inter-annual comparisons (2000-2005) and VCF-Quickbird image-interpreted estimates indicate that: (1) Pixel-level, inter-annual comparisons of VCF estimates of percent canopy cover were linearly related (mean R² = 0.77) and exhibited an average root mean square error (RMSE) of 10.1% and an average root mean square difference (RMSD) of 7.3%. (2) A comparison of image-interpreted percent tree crown cover estimates based on dot counts on Quickbird color images by two different interpreters were more variable (R² = 0.73, RMSE = 14.8%, RMSD = 18.7%) than VCF inter-annual comparisons. (3) Across the circumpolar boreal region, 2005 VCF-Quickbird comparisons were linearly related, with an R² = 0.57, a RMSE = 13.4% and a RMSD = 21.3%, with a tendency to over-estimate areas of low percent tree cover and anomalous VCF results in Scandinavia. The relationship of the VCF estimates and ground reference indicate to potential users that the VCF's tree cover values for individual pixels, particularly those below 20% tree cover, may not be precise enough to monitor 500 m pixel-level tree cover in the taiga-tundra transition zone.     

Estimation of diurnal air temperature using MSG SEVIRI data in West Africa

Spatially distributed air temperature data with high temporal resolution are desired for several modeling applications. By exploiting the thermal split window channels in combination with the red and near infrared channels of the geostationary MSG SEVIRI sensor, multiple daily air temperature estimates can be achieved using the contextual temperature-vegetation index method. Air temperature was estimated for 436 image acquisitions during the 2005 rainy season over West Africa and evaluated against in situ data from a field test site in Dahra, Northern Senegal. The methodology was adjusted using data from the test site resulting in RMSE = 2.55 K, MBE = − 0.30 K and R² = 0.63 for the estimated versus observed air temperatures. A spatial validation of the method using 12 synoptic weather stations from Senegal and Mali within the Senegal River basin resulted in overall values of RMSE = 2.96 K, MBE = − 1.11 K and R² = 0.68. The daytime temperature curve is interpolated using a sine function based on the multiple daily air temperature estimates from the SEVIRI data. These estimates (covering the 8:00-20:00 UCT time window) were in good agreement with observed values with RMSE = 2.99 K, MBE = − 0.70 K and R² = 0.64. The temperature-vegetation index method was applied as a moving window technique to produce distributed maps of air temperature with 15 min intervals and 3 km spatial resolution for application in a distributed hydrological model.