Inversion of a forest reflectance model to estimate structural canopy variables from hyperspectral remote sensing data

The potential of canopy reflectance modelling to retrieve simultaneously several structural variables in managed Norway spruce stands was investigated using the “Invertible Forest Reflectance Model”, INFORM. INFORM is an innovative extension of the FLIM model, with crown transparency, infinite crown reflectance and understory reflectance simulated using physically based sub-models (SAILH, LIBERTY and PROSPECT). The INFORM model was inverted with hyperspectral airborne HyMap data using a neural network approach. INFORM based estimates of forest structural variables were produced using site-specific ranges of stand structural variables. A relatively simple three layer feed-forward backpropagation neural network with two input neurons, one neuron in the hidden layer and three output neurons was employed to map leaf area index (LAI), crown coverage and stem density.To identify the optimum 2-band spectral subset to be used in the inversion process, all 2-band combinations of the HyMap dataset were systematically evaluated for model inversion. Field measurements of structural variables from 39 forest stands were used to validate the maps produced from HyMap imagery. Using two HyMap wavebands at 837 nm and 1148 nm the obtained accuracy of the LAI map amounts to an rmse of 0.58 (relative rmse = 18% of mean, R² = 0.73). With HyMap data resampled to Landsat TM spectral bands and using two “optimum” bands at 840 nm and 1650 nm, rmse was 0.66 and relative rmse 21%. In contrast to approaches based on empirical relations between spectral vegetation indices and structural variables, the main advantage of the inversion approach is that it does not require previous calibration.     

Imaging spectroscopy of jarosite cement in the Jurassic Navajo Sandstone

Imaging spectroscopy is a powerful tool for mapping surface mineralogy. Interactions of energy and surface materials on atomic and molecular levels result in specific absorption features that are diagnostic of mineralogy. Hyperspectral airborne sensors such as HyMap have sufficient spectral resolution to identify subtle features over narrow wavelength ranges. An anomalous zone of jarosite cement within the Jurassic Navajo Sandstone in southern Utah was identified with airborne hyperspectral data. To date no other locations of sulfate-rich cements such as this have been documented in the Navajo Sandstone. Here, we use multiple spectral analysis techniques over a micro to macro spatial continuum in order to focus our evaluation on the distribution and relative abundance of jarosite in this area. In this study, imaging spectroscopy is used to help determine the extent of the distribution of jarosite [KFe₃(SO₄)₂(OH)₆], at “Mollie's Nipple” (MN), a significant geomorphic landmark located within the Grand Staircase-Escalante National Monument, in the Jurassic Navajo Sandstone of south central Utah. The extent of the jarosite across the butte is identified in this study by mapping two diagnostic absorption features, referred to as the ∼ 2.26 μm and near-infrared ferric iron absorption features (∼ 0.70-1.15 μm), independently. This analysis shows that there is a single circular zone of abundant jarosite ∼ 1 km in diameter that includes both in-situ and weathered out accumulations. Jarosite formation requires acidic and sulfate-rich fluids, which are unusual for the Navajo Sandstone. Imaging and field spectroscopy provides the spectral resolution needed to map and analyze the mineralogic characteristics of this area; characteristics that may help constrain the conditions under which this atypical butte formed.     

Combining Isomap ordination and imaging spectroscopy to map continuous floristic gradients in a heterogeneous landscape

Mapping requires a meaningful generalization of information. For vegetation maps, classification is frequently used to generalize the species composition of (semi-)natural plant assemblages. As an alternative to classification, ordination methods aim to extract major floristic gradients describing the prevailing compositional variation in a floristic data set as metric variables. This ability has been used previously to derive gradient maps of homogeneous landscapes that show plant species composition in continuous fields. In the present study, gradient mapping was used in a more heterogeneous landscape with intricate environmental gradients and higher variation in vegetation physiognomy. Since established ordination methods may have difficulties to cope with the highly variable plant species composition, we tested the novel method Isometric Feature Mapping (Isomap) against conventional methods (Detrended Correspondence Analysis and Nonmetric Multidimensional Scaling). The resulting floristic gradients were related to hyperspectral imagery (HyMap) using partial least squares regression (PLSR) and subsequently mapped. Prediction uncertainties are provided as additional map. Isomap was able to preserve 74% of the original variation inherent to the floristic data set in a three-dimensional solution. This was considerably more than the established techniques achieved. The PLSR models for the floristic gradients extracted with Isomap showed model fits ranging from R² = 0.59 to R² = 0.73 in calibration and from R² = 0.55 to R² = 0.69 in tenfold cross-validation. The resulting gradient map provides detailed information on compositional vegetation patterns.     

Distribution of hydrothermally altered rocks in the Reko Diq, Pakistan mineralized area based on spectral analysis of ASTER data

The Reko Diq, Pakistan mineralized study area, approximately 10 km in diameter, is underlain by a central zone of hydrothermally altered rocks associated with Cu-Au mineralization. The surrounding country rocks are a variable mixture of unaltered volcanic rocks, fluvial deposits, and eolian quartz sand. Analysis of 15-band Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) data of the study area, aided by laboratory spectral reflectance and spectral emittance measurements of field samples, shows that phyllically altered rocks are laterally extensive, and contain localized areas of argillically altered rocks.In the visible through shortwave-infrared (VNIR + SWIR) phyllically altered rocks are characterized by Al-OH absorption in ASTER band 6 because of molecular vibrations in muscovite, whereas argillically altered rocks have an absorption feature in band 5 resulting from alunite. Propylitically altered rocks form a peripheral zone and are present in scattered exposures within the main altered area. Chlorite and muscovite cause distinctive absorption features at 2.33 and 2.20 μm, respectively, although less intense 2.33 μm absorption is also present in image spectra of country rocks.Important complementary lithologic information was derived by analysis of the spectral emittance data in the 5 thermal-infrared (TIR) bands. Silicified rocks were not distinguished in the 9 VNIR + SWIR bands because of the lack of diagnostic spectral absorption features in quartz in this wavelength region. Quartz-bearing surficial deposits, as well as hydrothermally silicified rocks, were mapped in the TIR bands by using a band 13/band 12 ratio image, which is sensitive to the intensity of the quartz reststrahlen feature. Improved distinction between the quartzose surficial deposits and silicified bedrock was achieved by using matched-filter processing with TIR image spectra for reference.     

Using ICESat's Geoscience Laser Altimeter System (GLAS) to assess large-scale forest disturbance caused by hurricane Katrina

In 2005, hurricane Katrina resulted in a large disturbance to U.S. forests. Recent estimates of damage from hurricane Katrina have relied primarily on optical remote sensing and field data. This paper is the first large-scale study to use satellite-based lidar data to quantify changes in forest structure from that event. GLAS data for the years prior to and following hurricane Katrina were compared to wind speed, forest cover, and damage data to assess the adequacy of sensor sampling, and to estimate changes in Mean Canopy Height (MCH) over all areas that experienced tropical force winds and greater. Statistically significant decreases in MCH post-Katrina were found to increase with wind intensity: Tropical Storm ?MCH = − 0.5 m, Category 1 ?MCH = − 2 m, and Category 2 ?MCH = − 4 m. A strong relationship was also found between changes in non-photosynthetic vegetation (?NPV), a metric previously shown to be related to storm damage, and post-storm MCH. The season of data acquisition was shown to influence calculations of MCH and MCH loss, but did not preclude the detection of major large-scale patterns of damage. Results from this study show promise for using space-borne lidar for large-scale assessments of forest disturbance, and highlight the need for future data on vegetation structure from space.     

Physically based vertical vegetation structure retrieval from ICESat data: Validation using LVIS in White Mountain National Forest, New Hampshire, USA

Vegetation structure retrieval accuracies from spaceborne Geoscience Laser Altimeter System (GLAS) on the Ice, Cloud and land Elevation Satellite (ICESat) data are affected by surface topography, background noise and sensor saturation. This study uses a physical approach to remove surface topography effect from lidar returns to retrieve vegetation height from ICESat/GLAS data over slope terrains. Slope-corrected vegetation heights from ICESat/GLAS data were compared to airborne Laser Vegetation Imaging Sensor (LVIS) (20 m footprint size) and small-footprint lidar data collected in White Mountain National Forest, NH. Impact of slope on LVIS vegetation height estimates was assessed by comparing LVIS height before and after slope correction with small-footprint discrete-return lidar and field data.Slope-corrected GLAS vegetation heights match well with 98 percentile heights from small-footprint lidar (R² = 0.77, RMSE = 2.2 m) and top three LVIS mean (slope-corrected) heights (R² = 0.64, RMSE = 3.7 m). Impact of slope on LVIS heights is small, however, comparison of LVIS heights (without slope correction) with either small footprint lidar or field data indicates that our scheme improves the overall LVIS height accuracy by 0.4-0.7 m in this region. Vegetation height can be overestimated by 3 m over a 15° slope without slope correction. More importantly, both slope-corrected GLAS and LVIS height differences are independent of slope. Our results demonstrate the effectiveness of the physical approach to remove surface topography from large footprint lidar data to improve accuracy of maximum vegetation height estimates.GLAS waveforms were compared to aggregated LVIS waveforms in Bartlett Experimental Forest, NH, to evaluate the impact of background noise and sensor saturation on vegetation structure retrievals from ICESat/GLAS. We found that GLAS waveforms with sensor saturation and low background noise match well with aggregated LVIS waveforms, indicating these waveforms capture vertical vegetation structure well. However, waveforms with large noise often lead to mismatched waveforms with LVIS and underestimation of waveform extent and vegetation height. These results demonstrate the quality of ICESat/GLAS vegetation structure estimates.     

Mapping vegetation in Yellowstone National Park using spectral feature analysis of AVIRIS data

Knowledge of the distribution of vegetation on the landscape can be used to investigate ecosystem functioning. The sizes and movements of animal populations can be linked to resources provided by different plant species. This paper demonstrates the application of imaging spectroscopy to the study of vegetation in Yellowstone National Park (Yellowstone) using spectral feature analysis of data from the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS). AVIRIS data, acquired on August 7, 1996, were calibrated to surface reflectance using a radiative transfer model and field reflectance measurements of a ground calibration site. A spectral library of canopy reflectance signatures was created by averaging pixels of the calibrated AVIRIS data over areas of known forest and nonforest vegetation cover types in Yellowstone. Using continuum removal and least squares fitting algorithms in the US Geological Survey's Tetracorder expert system, the distributions of these vegetation types were determined by comparing the absorption features of vegetation in the spectral library with the spectra from the AVIRIS data. The 0.68 μm chlorophyll absorption feature and leaf water absorption features, centered near 0.98 and 1.20 μm, were analyzed. Nonforest cover types of sagebrush, grasslands, willows, sedges, and other wetland vegetation were mapped in the Lamar Valley of Yellowstone. Conifer cover types of lodgepole pine, whitebark pine, Douglas fir, and mixed Engelmann spruce/subalpine fir forests were spectrally discriminated and their distributions mapped in the AVIRIS images. In the Mount Washburn area of Yellowstone, a comparison of the AVIRIS map of forest cover types to a map derived from air photos resulted in an overall agreement of 74.1% (kappa statistic=0.62).     

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.     

Quantifying the spatial distribution of soil mass wasting processes after the 2008 earthquake in Wenchuan, China: A case study of the Longmenshan area

The beautiful Longmenshan area is one of the main tourist attractions in Sichuan Province, China. The epicenter of a catastrophic earthquake measured at 8.0 Ms (China Seismological Bureau), occurred within this area at Wenchuan (31°01′16″N, 103°22′01″E) at 14:28 May 12, 2008 (Beijing time). The earthquake triggered numerous types of landslide transport and hazards, including soil and debris avalanches, rockfalls, slumps, debris flows, creation of barrier lakes and slope flattenings. This paper examines the landslide hazards in the Longmenshan area caused by the earthquake using remotely sensed images, mainly Beijing-1 Microsatellite data before and after the earthquake, compared to digital elevation maps and slope gradient maps, land use and vegetation cover maps. Areas of erosion and loss of vegetation were compared from pre- and post-earthquake data, from which were calculated changes in vegetated areas, bare slopes, and mass movement during the earthquake. These events occurred over altitudes from 1000 to 4000 m and on slope angles between 25 and 55°. The results show that the total area of erosion and land movement due to the earthquake increased by 86.3 km² (19.2% of the study area). Compared with pre-earthquake, the areas of very low intensity soil erosion and moderate intensity soil erosion were respectively reduced by 3.6 km², 24.3 km² and 30.9 km². On the other hand, the areas of severe and very severe intensity soil erosion were substantially increased by 45.8 km² and 99.2 km². In the post-earthquake stage, the bare areas (vegetation cover < 15%) have increased by 65.8 km². Without vegetation, the denuded earthquake damaged slopes and other high risk sites have become severe erosion problems. Thus, it is essential to continue long-term monitoring of mass wasting in the denuded areas and evaluate potential risk sites for future landslides and debris flows. We anticipate that these results will be helpful in decision making and policy planning for recovery and reconstruction in the earthquake-affected area.     

Estimating forage quantity and quality using aerial hyperspectral imagery for northern mixed-grass prairie

Sustainable rangeland stewardship calls for synoptic estimates of rangeland biomass quantity (kg dry matter ha^− 1) and quality [carbon:nitrogen (C:N) ratio]. These data are needed to support estimates of rangeland crude protein in forage, either by percent (CP_c) or by mass (CP_m). Biomass derived from remote sensing data is often compromised by the presence of both photosynthetically active (PV) and non-photosynthetically active (NPV) vegetation. Here, we explicitly quantify PV and NPV biomass using HyMap hyperspectral imagery. Biomass quality, defined as plant C:N ratio, was also estimated using a previously published algorithm. These independent algorithms for forage quantity and quality (both PV and NPV) were evaluated in two northern mixed-grass prairie ecoregions, one in the Northwestern Glaciated Plains (NGGP) and one in the Northwestern Great Plains (NGP). Total biomass (kg ha^− 1) and C:N ratios were mapped with 18% and 8% relative error, respectively. Outputs from both models were combined to quantify crude protein (kg ha^− 1) on a pasture scale. Results suggest synoptic maps of rangeland vegetation mass (both PV and NPV) and quality may be derived from hyperspectral aerial imagery with greater than 80% accuracy.     

A comparison of multi-spectral, multi-angular, and multi-temporal remote sensing datasets for fractional shrub canopy mapping in Arctic Alaska

Shrub cover appears to be increasing across many areas of the Arctic tundra biome, and increasing shrub cover in the Arctic has the potential to significantly impact global carbon budgets and the global climate system. For most of the Arctic, however, there is no existing baseline inventory of shrub canopy cover, as existing maps of Arctic vegetation provide little information about the density of shrub cover at a moderate spatial resolution across the region. Remotely-sensed fractional shrub canopy maps can provide this necessary baseline inventory of shrub cover. In this study, we compare the accuracy of fractional shrub canopy (> 0.5 m tall) maps derived from multi-spectral, multi-angular, and multi-temporal datasets from Landsat imagery at 30 m spatial resolution, Moderate Resolution Imaging SpectroRadiometer (MODIS) imagery at 250 m and 500 m spatial resolution, and MultiAngle Imaging Spectroradiometer (MISR) imagery at 275 m spatial resolution for a 1067 km² study area in Arctic Alaska. The study area is centered at 69 °N, ranges in elevation from 130 to 770 m, is composed primarily of rolling topography with gentle slopes less than 10°, and is free of glaciers and perennial snow cover. Shrubs > 0.5 m in height cover 2.9% of the study area and are primarily confined to patches associated with specific landscape features. Reference fractional shrub canopy is determined from in situ shrub canopy measurements and a high spatial resolution IKONOS image swath. Regression tree models are constructed to estimate fractional canopy cover at 250 m using different combinations of input data from Landsat, MODIS, and MISR. Results indicate that multi-spectral data provide substantially more accurate estimates of fractional shrub canopy cover than multi-angular or multi-temporal data. Higher spatial resolution datasets also provide more accurate estimates of fractional shrub canopy cover (aggregated to moderate spatial resolutions) than lower spatial resolution datasets, an expected result for a study area where most shrub cover is concentrated in narrow patches associated with rivers, drainages, and slopes. Including the middle infrared bands available from Landsat and MODIS in the regression tree models (in addition to the four standard visible and near-infrared spectral bands) typically results in a slight boost in accuracy. Including the multi-angular red band data available from MISR in the regression tree models, however, typically boosts accuracy more substantially, resulting in moderate resolution fractional shrub canopy estimates approaching the accuracy of estimates derived from the much higher spatial resolution Landsat sensor. Given the poor availability of snow and cloud-free Landsat scenes in many areas of the Arctic and the promising results demonstrated here by the MISR sensor, MISR may be the best choice for large area fractional shrub canopy mapping in the Alaskan Arctic for the period 2000-2009.     

Assessment of the impacts of surface topography, off-nadir pointing and vegetation structure on vegetation lidar waveforms using an extended geometric optical and radiative transfer model

In order to prioritize the measurement requirements and accuracies of the two new lidar missions, a physical model is required for a fundamental understanding of the impact of surface topography, footprint size and off-nadir pointing on vegetation lidar waveforms and vegetation height retrieval. In this study, we extended a well developed Geometric Optical and Radiative Transfer (GORT) vegetation lidar model to take into account for the impacts of surface topography and off-nadir pointing on vegetation lidar waveforms and vegetation height retrieval and applied this extended model to assess the aforementioned impacts on vegetation lidar waveforms and height retrieval.Model simulation shows that surface topography and off-nadir pointing angle stretch waveforms and the stretching effect magnifies with footprint size, slope and off-nadir pointing angle. For an off-nadir pointing laser penetrating vegetation over a slope terrain, the waveform is either stretched or compressed based on the relative angle. The stretching effect also results in a disappearing ground peak return when slope or off-nadir pointing angle is larger than the “critical slope angle”, which is closely related to various vegetation structures and footprint size. Model simulation indicates that waveform shapes are affected by surface topography, off-nadir pointing angle and vegetation structure and it is difficult to remove topography effects from waveform extent based only on the shapes of waveform without knowing any surface topography information.Height error without correction of surface topography and off-nadir pointing angle is the smallest when the laser beams at the toward-slope direction and the largest from the opposite direction. Further simulation reveals within 20° of slope and off-nadir pointing angle, given the canopy height as roughly 25 m and the footprint size as 25 m, the error for vegetation height (RH100) ranges from − 2 m to greater than 12 m, and the error for the height at the medium energy return (RH50) from − 1 m to 4 m. The RH100 error caused by unknown surface topography and without correction of off-nadir pointing effect can be explained by an analytical formula as a function of vegetation height, surface topography, off-nadir pointing angle and footprint size as a first order approximation. RH50 is not much affected by topography, off-nadir pointing and footprint size. This forward model simulation can provide scientific guidance on prioritizing future lidar mission measurement requirements and accuracies.     

Floodplain roughness parameterization using airborne laser scanning and spectral remote sensing

Floodplain roughness parameterization is one of the key elements of hydrodynamic modeling of river flow, which is directly linked to exceedance levels of the embankments of lowland fluvial areas. The present way of roughness mapping is based on manually delineated floodplain vegetation types, schematized as cylindrical elements of which the height (m) and the vertical density (the projected plant area in the direction of the flow per unit volume, m^− 1) have to be assigned using a lookup table. This paper presents a novel method of automated roughness parameterization. It delivers a spatially distributed roughness parameterization in an entire floodplain by fusion of CASI multispectral data with airborne laser scanning (ALS) data. The method consists of three stages: (1) pre-processing of the raw data, (2) image segmentation of the fused data set and classification into the dominant land cover classes (KHAT = 0.78), (3) determination of hydrodynamic roughness characteristics for each land cover class separately. In stage three, a lookup table provides numerical values that enable roughness calculation for the classes water, sand, paved area, meadows and built-up area. For forest and herbaceous vegetation, ALS data enable spatially detailed analysis of vegetation height and density. The hydrodynamic vegetation density of forest is mapped using a calibrated regression model. Herbaceous vegetation cover is further subdivided in single trees and non-woody vegetation. Single trees were delineated using a novel iterative cluster merging method, and their height is predicted (R² = 0.41, rse = 0.84 m). The vegetation density of single trees was determined in an identical way as for forest. Vegetation height and density of non-woody herbaceous vegetation were also determined using calibrated regression models. A 2D hydrodynamic model was applied with the results of this novel method, and compared with a traditional roughness parameterization approach. The modeling results showed that the new method is well able to provide accurate output data. The new method provides a faster, repeatable, and more accurate way of obtaining floodplain roughness, which enables regular updating of river flow models.     

MODIS spectral signals at a flux tower site: Relationships with high-resolution data, and CO2 flux and light use efficiency measurements

In this study, spectral indices were calculated from single date HyMap (3 m; 126 bands), Hyperion (30 m; 242 bands), ASTER (15/30 m; 9 bands), and a time series of MODIS nadir BRDF-adjusted reflectance (NBAR; 1 km, 7 bands) for a study area surrounding the Tumbarumba flux tower site in eastern Australia. The study involved: a) the calculation of a range of physiologically-based vegetation indices from ASTER, HyMap, Hyperion and MOD43B NBAR imagery over the flux tower site; b) comparison across scales between HyMap, Hyperion and MODIS for the normalized difference water index (NDWI) and the Red-Green ratio; c) analysis of relationships between tower-based flux and light use efficiency (LUE) measurements and seasonal and climatic constraints on growth; and d) examination of relationships between fluxes, LUE and time series of NDVI, NDWI and Red-Green ratio. Strong seasonal patterns of variation were observed in NDWI and Red/Green ratio from MODIS NBAR. Correlations between fine (3 and 30 m) and coarse (1 km) scale indices for a small region around the flux tower site were moderately good for Red/Green ratio, but poor for NDWI. Hymap NDWI values for the understorey canopy were much lower than values for the tree canopy. MODIS NDWI was negatively correlated with CO₂ fluxes during warm and cool seasons. The correlation indicated that surface reflectance, affected by a spectrally bright grassland understorey canopy, was decoupled from growth of trees with access to deep soil moisture. The application of physiologically-based indices at earth observation scale requires careful attention to applicability of band configurations, contribution of vegetation components to reflectance signals, mechanistic relationships between biochemical processes and spectral indexes, and incorporation of ancillary information into any analysis.     

Characterizing forest ecological structure using pulse types and heights of airborne laser scanning

Characterizing forest structure is an important part of any comprehensive biodiversity assessment. However, current methods for measuring structural complexity require a laborious process that involves many logistically expensive point based measurements. An automated or semi-automated method would be ideal. In this study, the utility of airborne laser scanning (LiDAR; Light Detection and Ranging) for characterizing the ecological structure of a forest landscape is examined. The innovation of this paper is to use different laser pulse return properties from a full waveform LiDAR to characterize forest ecological structure. First the LiDAR dataset is stratified into four vertical layers: ground, low vegetation (0-1 m from the ground), medium vegetation (1-5 m from the ground) and high vegetation (> 5 m). Subsequently the “Type” of LiDAR return is analysed: Type 1 (singular returns); Type 2 (first of many returns); Type 3 (intermediate returns); and Type 4 (last of many returns). A forest characterization scheme derived from LiDAR point clouds is proposed. A validation of the scheme is then presented using a network of field sites that recorded commonly used metrics of biodiversity. The proposed forest characterization categories allow for quantification of gaps (above bare ground, low vegetation and medium vegetation), canopy cover and its vertical density as well as the presence of various canopy strata (low, medium and high). Regression analysis showed that LiDAR derived variables were good predictors of field recorded variables (R² = 0.82, P < 0.05 between LiDAR derived presence of low vegetation and field derived LAI for low vegetation). The proposed scheme clearly shows the potential of full waveform LiDAR to provide information on the complexity of habitat structure.     

Responses of the reflectance indices PRI and NDVI to experimental warming and drought in European shrublands along a north-south climatic gradient

The aim of this study was to evaluate the use of ground-based canopy reflectance measurements to detect changes in physiology and structure of vegetation in response to experimental warming and drought treatment at six European shrublands located along a North-South climatic gradient. We measured canopy reflectance, effective green leaf area index (green LAIe) and chlorophyll fluorescence of dominant species. The treatment effects on green LAIe varied among sites. We calculated three reflectance indices: photochemical reflectance index PRI [531 nm; 570 nm], normalized difference vegetation index NDVI₆₈₀ [780 nm; 680 nm] using red spectral region, and NDVI₅₇₀ [780 nm; 570 nm] using the same green spectral region as PRI. All three reflectance indices were significantly related to green LAIe and were able to detect changes in shrubland vegetation among treatments. In general warming treatment increased PRI and drought treatment reduced NDVI values. The significant treatment effect on photochemical efficiency of plants detected with PRI could not be detected by fluorescence measurements. However, we found canopy level measured PRI to be very sensitive to soil reflectance properties especially in vegetation areas with low green LAIe. As both soil reflectance and LAI varied between northern and southern sites it is problematic to draw universal conclusions of climate-derived changes in all vegetation types based merely on PRI measurements. We propose that canopy level PRI measurements can be more useful in areas of dense vegetation and dark soils.     

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.     

Comparison between field measurements and airborne visible and infrared mapping spectrometry (AVIRIS and HyMap) of the Ronda peridotite massif (south-west Spain)

The high spectral resolution of Airborne Visible Infrared Imaging Spectrometer (AVIRIS; 224 channels from 400 to 2455?nm) and HyMap (127 channels between 437 and 2485?nm) images is necessary to conduct geological analysis with remote petrological determinations of rock types or soils, or to determine vegetation groups. When airborne images and field spectra are well adjusted between each other, and when the vegetation does not interfere in the analysis, the spectral shape analysis (SSA) method represents an easy treatment to reveal a large amount of geological information. The method presented in the paper takes into account both wells and peaks of spectra resulting from a combination of absorption features and continuum shapes. It was conducted on the Ronda peridotite, in the south-west of Spain, which was imaged by AVIRIS in 1991 and by HyMap in 2000, and which was partially sampled in the field using a GER 3700 spectrometer in 1997, 2000 and 2001.

In this study, the AVIRIS and HyMap images are processed to infer geological features using first photo-interpretation of colour composite images and then using the SSA method. This allows us to distinguish easily the peridotite massif from its surrounding rocks, to identify petrological variations inside the peridotite, including different varieties of serpentines, but also inside the surrounding rocks, particularly the carbonate-bearing ones which comprise marbles and dolomites.     

Comparative evaluation of airborne LiDAR and ship-based multibeam SoNAR bathymetry and intensity for mapping coral reef ecosystems

Large areas of the world's coastal marine environments remain poorly characterized because they have not been mapped with sufficient accuracy and at spatial resolutions high enough to support a wide range of societal needs. Expediting the rate of seafloor mapping requires the collection of multi-use datasets that concurrently address hydrographic charting needs and support decision-making in ecosystem-based management. While active optical and acoustic sensors have previously been compared for the purpose of hydrographic charting, few studies have evaluated the performance and cost effectiveness of these systems for providing benthic habitat maps. Bathymetric and intensity data were collected in shallow water (< 50 m depth) coral reef ecosystems using two conventional remote sensing technologies: (1) airborne Light Detection and Ranging (LiDAR), and (2) ship-based multibeam (MBES) Sound Navigation and Ranging (SoNAR). A comparative assessment using a suite of twelve metrics demonstrated that LiDAR and MBES were equally capable of discriminating seafloor topography (r = > 0.9), although LiDAR depths were found to be consistently shallower than MBES depths. The intensity datasets were not significantly correlated at a broad 4 × 5 km spatial scale (r = − 0.11), but were moderately correlated in flat areas at a fine 4 × 500 m spatial scale (r = 0.51), indicating that the LiDAR intensity algorithm needs to be improved before LiDAR intensity surfaces can be used for habitat mapping. LiDAR cost 6.6% less than MBES and required 40 fewer hours to map the same study area. MBES provided more detail about the seafloor by fully ensonifying high-relief features, by differentiating between fine and coarse sediments and by collecting data with higher spatial resolutions. Surface fractal dimensions and fast Fourier transformations emerged as useful methods for detecting artifacts in the datasets. Overall, LiDAR provided a more cost effective alternative to MBES for mapping and monitoring shallow water coral reef ecosystems (< 50 m depth), although the unique advantages of MBES may make it a more appropriate choice for answering certain ecological or geological questions requiring very high resolution data.     

Modeling fire severity in black spruce stands in the Alaskan boreal forest using spectral and non-spectral geospatial data

Biomass burning in the Alaskan interior is already a major disturbance and source of carbon emissions, and is likely to increase in response to the warming and drying predicted for the future climate. In addition to quantifying changes to the spatial and temporal patterns of burned areas, observing variations in severity is the key to studying the impact of changes to the fire regime on carbon cycling, energy budgets, and post-fire succession. Remote sensing indices of fire severity have not consistently been well-correlated with in situ observations of important severity characteristics in Alaskan black spruce stands, including depth of burning of the surface organic layer. The incorporation of ancillary data such as in situ observations and GIS layers with spectral data from Landsat TM/ETM+ greatly improved efforts to map the reduction of the organic layer in burned black spruce stands. Using a regression tree approach, the R² of the organic layer depth reduction models was 0.60 and 0.55 (p < 0.01) for relative and absolute depth reduction, respectively. All of the independent variables used by the regression tree to estimate burn depth can be obtained independently of field observations. Implementation of a gradient boosting algorithm improved the R² to 0.80 and 0.79 (p < 0.01) for absolute and relative organic layer depth reduction, respectively. Independent variables used in the regression tree model of burn depth included topographic position, remote sensing indices related to soil and vegetation characteristics, timing of the fire event, and meteorological data. Post-fire organic layer depth characteristics are determined for a large (> 200,000 ha) fire to identify areas that are potentially vulnerable to a shift in post-fire succession. This application showed that 12% of this fire event experienced fire severe enough to support a change in post-fire succession. We conclude that non-parametric models and ancillary data are useful in the modeling of the surface organic layer fire depth. Because quantitative differences in post-fire surface characteristics do not directly influence spectral properties, these modeling techniques provide better information than the use of remote sensing data alone.