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.     

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.     

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.     

The use of airborne lidar to assess avian species diversity, density, and occurrence in a pine/aspen forest

Vegetation structure is an important factor that influences wildlife-habitat selection, reproduction, and survival. However, field-based measurements of vegetation structure can be time consuming, costly, and difficult to undertake in areas that are remote and/or contain rough terrain. Light detection and ranging (lidar) is an active remote sensing technology that can quantify three-dimensional vegetation structure over large areas and thus holds promise for examining wildlife-habitat relationships. We used discrete-return airborne lidar data acquired over the Black Hills Experimental Forest in South Dakota, USA in combination with field-collected vegetation and bird data to assess the utility of lidar data in quantifying vegetation structural characteristics that relate to avian diversity, density, and occurrence. Indices of foliage height diversity calculated from lidar data were positively and significantly correlated with indices of bird species diversity, with the highest correlations observed when foliage height diversity categories contained proportionally more foliage layers near the forest floor (< 5 m). In addition, lidar-derived indices of vegetation volume were significantly correlated with bird density. Using lidar-derived vegetation height data in combination with multispectral IKONOS data, we delineated five general habitat types within the study area according to the presence of prominent vegetation layers at lower levels of the forest and predominant tree type (deciduous or conifer). Habitat type delineations were tested by examining the occurrence and relative density of two bird species common to the study area that prefer lower level vegetation for foraging and nesting. Dark-eyed Juncos were significantly associated with the 0.5–2.0 m high vegetation layer in pine-dominated stands, and Warbling Vireos were significantly associated with this same layer in aspen-dominated stands. These results demonstrate that discrete-return lidar can be an effective tool to remotely quantify vegetation structural attributes important to birds, and may be enhanced when used in combination with spectral data.     

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.     

A pseudo-waveform technique to assess forest structure using discrete lidar data

The use of airborne laser scanning systems (lidar) to describe forest structure has increased dramatically since height profiling experiments nearly 30 years ago. The analyses in most studies employ a suite of frequency-based metrics calculated from the lidar height data, which are systematically eliminated from a full model using stepwise multiple linear regression. The resulting models often include highly correlated predictors with little physical justification for model formulation. We propose a method to aggregate discrete lidar height and intensity measurements into larger footprints to create “pseudo-waves”. Specifically, the returns are first sorted into height bins, sliced into narrow discrete elements, and finally smoothed using a spline function. The resulting “pseudo-waves” have many of the same characteristics of traditional waveform lidar data. We compared our method to a traditional frequency-based method to estimate tree height, canopy structure, stem density, and stand biomass in coniferous and deciduous stands in northern Wisconsin (USA). We found that the pseudo-wave approach had strong correlations for nearly all tree measurements including height (cross validated adjusted R² (R²_cv) = 0.82, RMSE_cv = 2.09 m), mean stem diameter (R²_cv = 0.64, RMSE_cv = 6.15 cm), total aboveground biomass (R²_cv = 0.74, RMSE_cv = 74.03 kg ha^− 1), and canopy coverage (R²_cv = 0.79, RMSE_cv = 5%). Moreover, the type of wave (derived from height and intensity or from height alone) had little effect on model formulation and fit. When wave-based and frequency-based models were compared, fit and mean square error were comparable, leading us to conclude that the pseudo-wave approach is a viable alternative because it has 1) an increased breadth of available metrics; 2) the potential to establish new meaningful metrics that capture unique patterns within the waves; 3) the ability to explain metric selection based on the physical structure of forests; and 4) lower correlation among independent variables.     

On a new coordinate system for improved discrimination of vegetation and burned areas using MIR/NIR information

A transformation defined on the near-infrared (NIR) and middle-infrared (MIR) space is presented that allows deriving a new coordinate system appropriate for vegetation and burned area discrimination. The transformation is based on the difference between MIR and NIR in conjunction with the distance from a convergence point in the MIR/NIR space, representative of a totally burnt surface. One of the derived coordinates presents a small scatter for pixels associated to vegetated surfaces (strict scale) whereas the other one covers a wide range of values (large scale) that suggest its use as a proxy of water content of vegetation. The strict scale character of the first coordinate together with the large scale character of the second one make the coordinate system especially adequate to discriminate vegetated surfaces and rank them according to the water content, from green and dry to burned vegetation. The performance of the new coordinate system is then assessed against than traditional ratio or modified ratio indices (namely the Vegetation Index, the Burned Area Index and the Global Environmental Monitoring Index, modified to the MIR/NIR space) and it is shown that the new coordinate system provides better information than traditional indices, opening interesting perspectives for burned area discrimination and other applications like drought monitoring.     

Analysis of teleconnections between AVHRR-based sea surface temperature and vegetation productivity in the semi-arid Sahel

Vegetation productivity across the Sahel is known to be affected by a variety of global sea surface temperature (SST) patterns. Often climate indices are used to relate Sahelian vegetation variability to large-scale ocean-atmosphere phenomena. However, previous research findings reporting on the Sahelian vegetation response to climate indices have been inconsistent and contradictory, which could partly be caused by the variations in spatial extent/definitions of climate indices and size of the region studied. The aim of this study was to analyze the linkage between climate indices, pixel-wise spatio-temporal patterns of global sea surface temperature and the Sahelian vegetation dynamics for 1982-2007. We stratified the Sahel into five subregions to account for the longitudinal variability in rainfall. We found significant correlations between climate indices and the Normalized Difference Vegetation Index (NDVI) in the Sahel, however with different magnitudes in terms of strength for the western, central and eastern Sahel. Also the correlations based on NDVI and global SST anomalies revealed the same East-West gradient, with a stronger association for the western than the eastern Sahel. Warmer than average SSTs throughout the Mediterranean basin seem to be associated with enhanced greenness over the central Sahel whereas colder than average SSTs in the Pacific and warmer than average SSTs in the eastern Atlantic were related to increased greenness in the most western Sahel. Accordingly, we achieved high correlations for SSTs of oceanic basins which are geographically associated to the climate indices yet by far not always these patterns were coherent. The detected SST-NDVI patterns could provide the basis to develop new means for improved forecasts in particular of the western Sahelian vegetation productivity.     

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.     

Laser remote sensing of canopy habitat heterogeneity as a predictor of bird species richness in an eastern temperate forest, USA

Habitat heterogeneity has long been recognized as a fundamental variable indicative of species diversity, in terms of both richness and abundance. Satellite remote sensing data sets can be useful for quantifying habitat heterogeneity across a range of spatial scales. Past remote sensing analyses of species diversity have largely been limited to correlative studies based on the use of vegetation indices or derived land cover maps. A relatively new form of laser remote sensing (lidar) provides another means to acquire information on habitat heterogeneity. Here we examine the efficacy of lidar metrics of canopy structural diversity as predictors of bird species richness in the temperate forests of Maryland, USA. Canopy height, topography and the vertical distribution of canopy elements were derived from lidar imagery of the Patuxent National Wildlife Refuge and compared to bird survey data collected at referenced grid locations. The canopy vertical distribution information was consistently found to be the strongest predictor of species richness, and this was predicted best when stratified into guilds dominated by forest, scrub, suburban and wetland species. Similar lidar variables were selected as primary predictors across guilds. Generalized linear and additive models, as well as binary hierarchical regression trees produced similar results. The lidar metrics were also consistently better predictors than traditional remotely sensed variables such as canopy cover, indicating that lidar provides a valuable resource for biodiversity research applications.     

遥感和地理信息系统技术在湿地研究中的应用

湿地研究的关键在于定量化获取和分析湿地信息,RS和GIS技术为湿地研究提供了新的方法和技术支持。从湿地资源调查、湿地识别和湿地分类、湿地动态变化监测、湿地制图、湿地与环境因素定量模型、湿地信息系统6个方面,对国内外学者利用RS和GIS进行湿地研究的工作进展予以回顾。在此基础上概要说明RS和GIS技术在湿地研究中存在的问题及发展前景,指出RS和GIS技术在湿地研究中的应用从拓宽研究领域和完善技术体系两方面向更高、更深层次发展。     

利用四叉树解决激光表演系统的布点优化问题

为了绘制一幅图像,激光表演器需要的时间依赖于在该图像中线段绘制的顺序和方向。通过重新排列线段方向和顺序,可以大幅降低绘制时间。由于不同等级的激光表演器所标称的每秒绘制点数有所不同,通过降低绘制时间。可以令不同的激光表演器均可有效、完整地展现图像内容。     

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.     

Spatial and temporal dynamics of river channel migration and vegetation in central Amazonian white-water floodplains by remote-sensing techniques

We investigated spatial and temporal migration of the Solimões, the Japurá, and the Aranapu River channels in western Brazilian Amazonia with Landsat TM imagery over a 21-year period. Additionally, we classified and monitored how channel migrations affect the distribution of pioneer vegetation and old-growth forest. The cloud-free study area was 153,032 ha — open water plus 3 km inland on each margin. The channel migration rates, expressed as percent dislocation of the open water body of the river year^− 1, were lowest in the Japurá River (1.2%), and highest in the Aranapu channel (2.5%), the point bars at river confluence being the most affected landforms subject to geomorphic changes. Annual rates of lateral erosion and accretion of vegetated land along the three rivers were well-balanced. They averaged 0.79 and 0.83% of the cloud-free channel area over the 21 years. The Solimões River was more dynamic than the Japurá River, which can be traced to higher water discharge and sediment load. During the 21 years, the area covered by pioneer vegetation increased by 5.8% of the study area, while late-succession areas decreased by a similar amount (5.5%). According to local biomass estimates of the different vegetation types, these values suggest that C-releases by alluvial erosion would be much higher than C-sequestration caused by the creation of areas suitable for colonization by pioneer vegetation at our study site.     

微粒群算法在求解方程组中的应用

研究了利用微粒群算法求解线性、非线性方程组解的问题。对于线性、非线性方程组可以在指定的搜索区间内获得方程组的实数解。最后在计算机上进行了实验,证明了方法的正确性。     

遥感地质线性体与化探信息数字复合处理提取金矿控矿构造信息

从遥感图像获取的地质线性体与土壤岩石化探异常数据处理来进行定量分析提取具有控矿构造意义的遥感地质线性体,并从遥感地质环形体影像与化探金异常数据的关系,间接地分析环形体在成矿、控矿中的作用,为地面地质工作者提供找矿构造数值分析依据。研究结果与过去做过的大量地面成矿、控矿地质构造分析结果非常吻合。