林分平均高度卫星遥感新进展

林分平均高度是生态系统模型重要的输入参数之一,与生物量估算与碳循环研究高度相关。通过系统回顾林分平均高度研究的发展历史和最新进展,总结了不同传感器林分平均高度(SAR树高与LiDAR冠层高度)研究的主要模型和方法,通过单一传感器技术进行林分平均高度研究的内在特征的不同,分析了多传感器的区域林分平均高度联合反演方法的优劣性,并从科学发展趋势和社会需求出发,认识目前存在的主要问题与难点及未来面临的挑战和机遇,为今后更好地进行森林垂直结构和全球碳循环研究提供思路和借鉴。     

Extracting LiDAR indices to characterise multilayered forest structure using mixture distribution functions

Discrete Light Detection and Ranging (LiDAR) data is used to stratify a multilayered eucalyptus forest and characterise the structure of the vertical profile. We present a methodology that may prove useful for a very broad range of forest management applications, particularly for timber inventory evaluation and forest growth modelling. In this study, we use LiDAR data to stratify a multilayered eucalyptus forest and characterise the structure of specific vegetation layers for forest hydrology research, as vegetation dynamics influence a catchment's streamflow yield. A forest stand's crown height, density, depth, and closure, influence aerodynamic properties of the forest structure and the amount of transpiring leaf area, which in turn determine evapotranspiration rates. We present a methodology that produces canopy profile indices of understorey and overstorey vegetation using mixture models with a wide range of theoretical distribution functions. Mixture models provide a mechanism to summarise complex canopy attributes into a short list of parameters that can be empirically analysed against stand characteristics.Few studies have explored theoretical distribution functions to represent the vertical profile of vegetation structure in LiDAR data. All prior studies have focused on a Weibull distribution function, which is unimodal. In a complex native forest ecosystem, the form of the distribution of LiDAR points may be highly variable between forest types and age classes. We compared 44 probability distributions within a two component mixture model to determine the most suitable bimodal distributions for representing LiDAR density estimates of Mountain Ash forests in south-eastern Australia. An elimination procedure identified eleven candidate distributions for representing the eucalyptus component of the mixture model.We demonstrate the methodology on a sample of plots to predict overstorey stand volumes and basal area, and understorey basal area of 18-, 37-, and 70-year old Mountain Ash forest with variable density classes. The 70-year old forest has been subjected to a range of treatments including: thinning of the eucalyptus layer with two distinct retention rates, removal of the understorey, and clear felling of patches that have 37 year old regenerating forest. We demonstrate that the methodology has clear potential, as observed versus predicted values of eucalyptus basal area and stand volume were highly correlated, with bootstrap based r² ranging from 0.61 to 0.89 and 0.67 to 0.88 respectively. Non-eucalyptus basal area r² ranged from 0.5 to 0.91.     

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 snags and understory shrubs for a LiDAR-based assessment of wildlife habitat suitability

The lack of maps depicting forest three-dimensional structure, particularly as pertaining to snags and understory shrub species distribution, is a major limitation for managing wildlife habitat in forests. Developing new techniques to remotely map snags and understory shrubs is therefore an important need. To address this, we first evaluated the use of LiDAR data for mapping the presence/absence of understory shrub species and different snag diameter classes important for birds (i.e. ≥ 15 cm, ≥ 25 cm and ≥ 30 cm) in a 30,000 ha mixed-conifer forest in Northern Idaho (USA). We used forest inventory plots, LiDAR-derived metrics, and the Random Forest algorithm to achieve classification accuracies of 83% for the understory shrubs and 86% to 88% for the different snag diameter classes. Second, we evaluated the use of LiDAR data for mapping wildlife habitat suitability using four avian species (one flycatcher and three woodpeckers) as case studies. For this, we integrated LiDAR-derived products of forest structure with available models of habitat suitability to derive a variety of species-habitat associations (and therefore habitat suitability patterns) across the study area. We found that the value of LiDAR resided in the ability to quantify 1) ecological variables that are known to influence the distribution of understory vegetation and snags, such as canopy cover, topography, and forest succession, and 2) direct structural metrics that indicate or suggest the presence of shrubs and snags, such as the percent of vegetation returns in the lower strata of the canopy (for the shrubs) and the vertical heterogeneity of the forest canopy (for the snags). When applied to wildlife habitat assessment, these new LiDAR-based maps refined habitat predictions in ways not previously attainable using other remote sensing technologies. This study highlights new value of LiDAR in characterizing key forest structure components important for wildlife, and warrants further applications to other forested environments and wildlife species.     

Assessment of forest structure with airborne LiDAR and the effects of platform altitude

Airborne scanning LiDAR is a spatial technology increasingly used for forestry and environmental applications. However, the accuracy and coverage of LiDAR observations is highly dependent on both the extrinsic specifications of the LiDAR survey as well as the intrinsic effects such as the underlying forest structure. Extrinsic parameters which are set as part of the LiDAR survey include platform altitude, scan angle (half max. angle off nadir), and beam cross sectional diameter at the reflecting surface (referred to as footprint size). In this paper we investigate the effect of a number of these extrinsic parameters, including three different platform altitudes (1000, 2000, and 3000 m), two scan angles at 1000 m (10° and 15° half max. angle off nadir), and three footprint sizes (0.2, 0.4, and 0.6 m). The comparison was undertaken in eucalypt forests at three sites, varying in vegetation structure and topography within the Wedding Bells State Forest, Coffs Harbour, Australia. Results at the plot scale (40 × 90 m areas) indicate that tree heights computed from the 1000 m LiDAR data set (10° half max. angle off nadir) are well correlated with maximum plot heights (difference < 3 m) and field measured canopy volume (r² > 0.75, p < 0.001). Using normalised canopy height profiles (CHP) derived for sites, from data recorded at each altitude, we observed no significant difference between the relative distribution of LiDAR returns, indicating that platform altitude and footprint size have not had a major influence on CHP estimation. Interestingly, comparisons of first and last returns for individual pulses at increasing altitudes identified progressively fewer discrete first/last pulse combinations with more than 70% of pulses recorded as a single return at the highest altitude (3000 m). A possible hypothesis is that greater platform altitude and footprint size reduces the intensity of laser beam incident on a given surface area thus decreasing the probability of recording a last return above the noise threshold. Furthermore, tree scale analysis found a positive relationship between platform altitude and the underestimation of crown area and crown volume. The implications of this work for forest management are: (i) platform altitudes as high as 3000 m can be used to quantify the vertical distribution of phyto-elements, (ii) higher platform altitudes record a lower proportion of first/last return combinations that will further reduce the number of points available for forest structural assessment and development of digital elevation models, and (iii) for discrete LiDAR data, increasing platform altitude will record a lower frequency of returns per crown, resulting in larger underestimates of individual tree crown area and volume if standard algorithms are applied.     

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.     

Monitoring forest changes in the southwestern United States using multitemporal Landsat data

Landsat time series data sets were acquired for the Santa Fe National Forest in New Mexico. This area includes the San Pedro Parks Wilderness area, which was designated as an official wilderness in 1964. Eight autumnal Landsat Thematic Mapper (TM) scenes acquired from 1988 to 2006 were analyzed to determine whether significant changes have occurred throughout the region during the past 18 years and, if so, to assess whether the changes are long-term and gradual or short-term and abrupt. It was found that, starting in about 1995, many of the conifer stands within the Wilderness area showed consistently gradual and marked increases in the Shortwave Infrared/Near Infrared Index. These trends generally imply decreases in canopy greenness or increases in mortality. Other high-elevation conifer forests located outside of the Wilderness area showed similar spectral trends, indicating that changes are potentially widespread. The spatial patterns of forest damage as inferred from the image analyses were very similar to the general patterns of insect defoliation damage mapped via aerial sketch mapping by the United States Department of Agriculture Forest Service Forest Health Monitoring Program. A field visit indicated that zones of spectral change are associated with high levels of forest damage and mortality, likely caused by a combination of insects and drought. The study demonstrates the effectiveness of using historical Landsat data for providing objective and consistent long-term assessments of the gradual ecosystem changes that are occurring within the western United States.     

Development of a simulation model to predict LiDAR interception in forested environments

Airborne scanning LiDAR systems are used to predict a range of forest attributes. However, the accuracy with which this can be achieved is highly dependent on the sensor configuration and the structural characteristics of the forest examined. As a result, there is a need to understand laser light interactions with forest canopies so that LiDAR sensor configurations can be optimised to assess particular forest types. Such optimisation will not only ensure the targeted forest attributes can be accurately and consistently quantified, but may also minimise the cost of data acquisition and indicate when a survey configuration will not deliver information needs.In this paper, we detail the development and application of a model to simulate laser interactions within forested environments. The developed model, known as the LiDAR Interception and Tree Environment (LITE) model, utilises a range of structural configurations to simulate trees with variable heights, crown dimensions and foliage clumping. We developed and validated the LITE model using field data obtained from three forested sites covering a range of structural classes. Model simulations were then compared to coincident airborne LiDAR data collected over the same sites. Results indicate that the LITE model can be used to produce comparable estimates of maximum height of trees within plots (differences < 2.42 m), mean heights of first return data (differences < 2.27 m), and canopy height percentiles (r² = 0.94, p < 0.001) when compared to airborne LiDAR data. In addition, the distribution of airborne LiDAR hits through the canopy profile was closely matched by model predictions across the range of sites. Importantly, this demonstrates that the structural differences between forest stands can be characterised by LITE. Models that are capable of interpreting the response of small-footprint LiDAR waveforms can facilitate algorithm development, the generation of corrections for actual LiDAR data, and the optimisation of sensor configurations for differing forest types, benefiting a range of experimental and commercial LiDAR applications. As a result, we also performed a scenario analysis to demonstrate how differences in forest structure, terrain, and sensor configuration can influence the interception of LiDAR beams.     

Retrieval of canopy height using moderate-resolution imaging spectroradiometer (MODIS) data

In this study we use the 500 m Moderate Resolution Imaging Spectroradiometer (MODIS) Bidirectional Reflectance Distribution Function (BRDF) product to develop multivariate linear regression models that estimate canopy heights over study sites at Howland Forest, Maine, Harvard Forest, Massachusetts and La Selva Forest, Costa Rica using (1) directional escape probabilities that are spectrally independent and (2) the directional spectral reflectances used to derive the directional escape probabilities. These measures of canopy architecture are compared with canopy height information retrieved from the airborne Laser Vegetation Imaging Sensor (LVIS). Both the escape probability and the directional reflectance approaches achieve good results, with correlation coefficients in the range 0.54-0.82, although escape probability results are usually slightly better. This suggests that MODIS 500 m BRDF data can be used to extrapolate canopy heights observed by widely-spaced satellite LIDAR swaths to larger areas, thus providing wide-area coverage of canopy height.     

Fusion of LiDAR and imagery for estimating forest canopy fuels

Due to increased fuel loading as a result of fire suppression, land managers in the American west are in need of precise information about the fuels they manage, including canopy fuels. Canopy fuel metrics such as canopy height (CH), canopy base height (CBH), canopy bulk density (CBD) and available canopy fuel (ACF) are specific inputs for wildfire behavior models such as FARSITE and emission models such as FOFEM. With finer spatial resolution data, accurate quantification of these metrics with detailed spatial heterogeneity can be accomplished. Light Detection and Ranging (LiDAR) and color near-infrared imagery are active and passive systems, respectively, that have been utilized for measuring a range of forest structure characteristics at high resolution. The objective of this research was to determine which remote sensing dataset can estimate canopy fuels more accurately and whether a fusion of these datasets produces more accurate estimates. Regression models were developed for ponderosa pine (Pinus ponderosa) stand representative of eastern Washington State using field data collected in the Ahtanum State Forest and metrics derived from LiDAR and imagery. Strong relationships were found with LiDAR alone and LiDAR was found to increase canopy fuel accuracy compared to imagery. Fusing LiDAR with imagery and/or LiDAR intensity led to small increases in estimation accuracy over LiDAR alone. By improving the ability to estimate canopy fuels at higher spatial resolutions, spatially explicit fuel layers can be created and used in wildfire behavior and smoke emission models leading to more accurate estimations of crown fire risk and smoke related emissions.     

Model-assisted regional forest biomass estimation using LiDAR and InSAR as auxiliary data: A case study from a boreal forest area

There is a need for accurate inventory methods that produce relevant and timely information on the forest resources and carbon stocks for forest management planning and for implementation of national strategies under the United Nations Collaborative Program on Reduced Emissions from Deforestation and Forest Degradation in Developing Countries (REDD). Such methods should produce information that is consistent across various geographical scales. Airborne scanning Light Detection and Ranging (LiDAR) is among the most promising remote sensing technologies for estimation of forest resource information such as timber volume and biomass, while acquisition of three dimensional data with Interferometric Synthetic Aperture Radar (InSAR) from space is seen as a relevant option for inventory in the tropics because of its ability to “see through the clouds” and its potential for frequent updates at low costs. Based on a stratified probability sample of 201 field survey plots collected in a 960 km² boreal forest area in Norway, we demonstrate how total above-ground biomass (AGB) can be estimated at three distinct geographical levels in such a way that the estimates at a smaller level always sum up to the estimate at a larger level. The three levels are (1) a district (the entire study area), (2) a village, local community or estate level, and (3) a stand or patch level. The LiDAR and InSAR data were treated as auxiliary information in the estimation. At the two largest geographical levels model-assisted estimators were employed. A model-based estimation was conducted at the smallest level. Estimates of AGB and corresponding error estimates based on (1) the field sample survey were compared with estimates obtained by using (2) LiDAR and (3) InSAR data as auxiliary information. For the entire study area, the estimates of AGB were 116.0, 101.2, and 111.3 Mg ha⁻¹, respectively. Corresponding standard error estimates were 3.7, 1.6, and 3.2 Mg ha⁻¹. At the smallest geographical level (stand) an independent validation on 35 large field plots was carried out. RMSE values of 17.1-17.3 Mg ha⁻¹ and 42.6-53.2 Mg ha⁻¹ were found for LiDAR and InSAR, respectively. A time lag of six years between acquisition of InSAR data and field inventory has introduced some errors. Significant differences between estimates and reference values were found, illustrating the risk of using pure model-based methods in the estimation when there is a lack of fit in the models. We conclude that the examined remote sensing techniques can provide biomass estimates with smaller estimated errors than a field-based sample survey. The improvement can be highly significant, especially for LiDAR.     

Forest canopy height from the Multiangle Imaging SpectroRadiometer (MISR) assessed with high resolution discrete return lidar

In this study retrievals of forest canopy height were obtained through adjustment of a simple geometric-optical (GO) model against red band surface bidirectional reflectance estimates from NASA's Multiangle Imaging SpectroRadiometer (MISR), mapped to a 250 m grid. The soil-understory background contribution was partly isolated prior to inversion using regression relationships with the isotropic, geometric, and volume scattering kernel weights of a Li-Ross kernel-driven bidirectional reflectance distribution function (BRDF) model. The height retrievals were assessed using discrete return lidar data acquired over sites in Colorado as part of the Cold Land Processes Experiment (CLPX) and used with fractional crown cover retrievals to obtain aboveground woody biomass estimates. For all model runs with reasonable backgrounds and initial b/r (vertical to horizontal crown radii) values < 2.0, root mean square error (RMSE) distributions were centered between 2.5 and 3.7 m while R² distributions were centered between 0.4 and 0.7. The MISR/GO aboveground biomass estimates predicted via regression on fractional cover and mean canopy height for the CLPX sites showed good agreement with U.S. Forest Service Interior West map data (adjusted R² = 0.84). The implication is that multiangle sensors such as MISR can provide spatially contiguous retrievals of forest canopy height, cover, and aboveground woody biomass that are potentially useful in mapping distributions of aboveground carbon stocks, tracking disturbance, and in initializing, constraining, and validating ecosystem models. This is important because the MISR record is spatially comprehensive and extends back to the year 2000 and the launch of the NASA Earth Observing System (EOS) Terra satellite; it might thus provide a ~ 10-year baseline record that would enhance exploitation of data from the NASA Deformation, Ecosystem Structure and Dynamics of Ice (DESDynI) mission, as well as furthering realization of synergies with active instruments.     

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.     

Delineation of secondary succession mechanisms for tropical dry forests using LiDAR

Research suggests that secondary forest structure and composition subsequent to human disturbance could be related to wind or vertebrate dependent seed dispersal mechanisms. This paper examines the capability of a waveform LiDAR system (LVIS: LiDAR Vegetation Imaging System) to detect the location of wind-dispersed or vertebrate-dispersed forest patches in a mosaic of a secondary Tropical Dry Forest in the Santa Rosa National Park, Guanacaste, Costa Rica. Canopy-height estimations derived from LVIS data were used to identify and locate two types of dispersal-dependent arrangements of canopy heights: flat-topped (windborne seed dispersed) or dome-shaped (vertebrate-dispersed) fragments. Following identification, a Ripley's K indicator was used to compare the locations for each confirmed canopy type to those reported in existing maps. Results indicated that dome-shaped patches tended to be arranged in clusters at the middle of former pastures (far from older forests), and that flat-topped patches tended to be located adjacent to older forests, forming extensions from the old growth patches and at the downwind side of the slope. Such characteristics agree with the theoretical observation of the spatial configuration of vertebrate- and wind-dispersed fragments reported in the study site. Our findings demonstrate that processes controlling forest regeneration in the TDF can be successfully identified through LiDAR data.     

Simulated impact of sample plot size and co-registration error on the accuracy and uncertainty of LiDAR-derived estimates of forest stand biomass

Regression has been widely applied in Light Detection And Ranging (LiDAR) remote sensing to spatially extend predictions of total aboveground biomass (TAGB) and other biophysical properties over large forested areas. Sample (field) plot size has long been considered a key sampling design parameter and focal point for optimization in forest surveys, because of its impact on sampling effort and the estimation accuracy of forest inventory attributes. In this study, we demonstrate how plot size and co-registration error interact to influence the estimation of LiDAR canopy height and density metrics, regression model coefficients, and the prediction accuracy of least-squares estimators of TAGB. We made use of simulated forest canopies and synthetic LiDAR point clouds, so that we could maintain strict control over the spatial scale and complexity of forest scenes, as well as the magnitude and type of planimetric error inherent in ground-reference and LiDAR datasets. Our results showed that predictions of TAGB improved markedly as plot size increased from 314 (10 m radius) to 1964 m² (25 m radius). The co-registration error (spatial overlap) between ground-reference and LiDAR samples negatively impacted the estimation of LiDAR metrics, regression model fit, and the prediction accuracy of TAGB. We found that larger plots maintained a higher degree of spatial overlap between ground-reference and LiDAR datasets for any given GPS error, and were therefore more resilient to the ill effects of co-registration error compared to small plots. The impact of co-registration error was more pronounced in tall, spatially heterogeneous stands than short, homogeneous stands. We identify and briefly discuss three possible ways that LiDAR data could be used to optimize plot size, sample selection, and the deployment of GPS resources in forest biomass surveys.     

基于机载激光雷达数据的森林结构参数反演

机载激光雷达(Light Detection And Ranging,LiDAR)技术对植被空间结构和地形的探测能力较强,在植被参数定量测量和反演方面具有显著优势。首先利用野外调查并结合高分辨率Geoeye-1影像数据,对黑河上游天涝池流域植被类型进行分类,提取研究区森林分布,然后结合0.5m×0.5m机载激光雷达(LiDAR)数据对森林结构参数(树高、冠幅、胸径和叶面积指数)进行反演,最后利用实际观测数据对反演结果进行验证。结果表明:机载激光雷达数据能够精确地反演森林结构参数,树高、冠幅、胸径和叶面积指数的实测值与估测值决定系数分别为0.98、0.84、0.57和0.73。本研究获得流域森林覆盖区域高精度树冠高度和叶面积指数空间分布图,同时分析了冠层高度和叶面积指数随高度的变化。本研究的结果为该流域分布式生态水文模型提供了重要的输入参数。     

Quantifying forest above ground carbon content using LiDAR remote sensing

The UNFCCC and interest in the source of the missing terrestrial carbon sink are prompting research and development into methods for carbon accounting in forest ecosystems. Here we present a canopy height quantile-based approach for quantifying above ground carbon content (AGCC) in a temperate deciduous woodland, by means of a discrete-return, small-footprint airborne LiDAR. Fieldwork was conducted in Monks Wood National Nature Reserve UK to estimate the AGCC of five stands from forest mensuration and allometric relations. In parallel, a digital canopy height model (DCHM) and a digital terrain model (DTM) were derived from elevation measurements obtained by means of an Optech Airborne Laser Terrain Mapper 1210. A quantile-based approach was adopted to select a representative statistic of height distributions per plot. A forestry yield model was selected as a basis to estimate stemwood volume per plot from these heights metrics. Agreement of r=0.74 at the plot level was achieved between ground-based AGCC estimates and those derived from the DCHM. Using a 20×20 m grids superposed to the DCHM, the AGCC was estimated at the stand level and at the woodland level. At the stand level, the agreement between the plot data upscaled in proportion to area and the LiDAR estimates was r=0.85. At the woodland level, LiDAR estimates were nearly 24% lower than those from the upscaled plot data. This suggests that field-based approaches alone may not be adequate for carbon accounting in heterogeneous forests. Conversely, the LiDAR 20×20 m grid approach has an enhanced capability of monitoring the natural variability of AGCC across the woodland.     

Strengths and limitations of assessing forest density and spatial configuration with aerial LiDAR

Changes in the structural state of forests of the semi-arid U.S.A., such as an increase in tree density, are widely believed to be leading to an ecological crisis, but accurate methods of quantifying forest density and configuration are lacking at landscape scales. An individual tree canopy (ITC) method based on aerial LiDAR has been developed to assess forest structure by estimating the density and spatial configuration of trees in four different height classes. The method has been tested against field measured forest inventory data from two geographically distinct forests with independent LiDAR acquisitions. The results show two distinct patterns: accurate, unbiased density estimates for trees taller than 20 m, and underestimation of density in trees less than 20 m tall. The underestimation of smaller trees is suggested to be a limitation of LiDAR remote sensing. Ecological applications of the method are demonstrated through landscape metrics analysis of density and configuration rasters.     

Relationship of a Landsat cumulative disturbance index to canopy nitrogen and forest structure

Spatially-explicit knowledge of the timing, frequency, and intensity of forest disturbances is essential for forest management, yet little is known about how disturbances such as forest harvests and insect outbreaks might accumulate in their effects over time. Capturing the many forest harvest and insect defoliation events occurring over twenty-five years, we transformed a series of Landsat images into cumulative disturbance maps covering Green Ridge State Forest (GRSF) and Savage River State Forest (SRSF) in western Maryland. These maps summed yearly ΔDI images, which were defined as the change in a yearly tasseled cap disturbance index (DI), relative to a synthetic reference condition map created by finding the minimum DI value for all years. Intensive field-plot surveys and AVIRIS imagery collected during the summer of 2009 provided measurements of forest structure and canopy nitrogen. With these data, we found that while the most recent year's ΔDI had little relation, increases in the cumulative DI were related to decreased field-measured current canopy cover (R² = 0.66 at GRSF, 0.67 at SRSF and 0.34 combined) and watershed-averaged AVIRIS canopy N (R² = 0.40 at GRSF, 0.57 at SRSF and 0.54 combined). The latter relationship was obscured at the field-plot level of analysis, suggesting that fine scale studies will also need to account for other drivers (e.g. species composition) of variability in canopy N. Nevertheless, our study demonstrates that Landsat time series data can be synthesized into cumulative metrics incorporating multiple disturbance types, which help explain important cumulative disturbance-mediated changes in ecosystem functioning.     

Exploring LiDAR-RaDAR synergy—predicting aboveground biomass in a southwestern ponderosa pine forest using LiDAR, SAR and InSAR

Scanning Light Detecting and Ranging (LiDAR), Synthetic Aperture Radar (SAR) and Interferometric SAR (InSAR) were analyzed to determine (1) which of the three sensor systems most accurately predicted forest biomass, and (2) if LiDAR and SAR/InSAR data sets, jointly considered, produced more accurate, precise results relative to those same data sets considered separately. LiDAR ranging measurements, VHF-SAR cross-sectional returns, and X- and P-band cross-sectional returns and interferometric ranges were regressed with ground-estimated (from dbh) forest biomass in ponderosa pine forests in the southwestern United States. All models were cross-validated. Results indicated that the average canopy height measured by the scanning LiDAR produced the best predictive equation. The simple linear LiDAR equation explained 83% of the biomass variability (n = 52 plots) with a cross-validated root mean square error of 26.0 t/ha. Additional LiDAR metrics were not significant to the model. The GeoSAR P-band (λ = 86 cm) cross-sectional return and the GeoSAR/InSAR canopy height (X-P) captured 30% of the forest biomass variation with an average predictive error of 52.5 t/ha. A second RaDAR-FOPEN collected VHF (λ ∼ 7.8 m) and cross-polarized P-band (λ = 88 cm) cross-sectional returns, none of which proved useful for forest biomass estimation (cross-validated R² = 0.09, RMSE = 63.7 t/ha). Joint consideration of LiDAR and RaDAR measurements produced a statistically significant, albeit small improvement in biomass estimation precision. The cross-validated R² increased from 83% to 84% and the prediction error decreased from 26.0 t/ha to 24.9 t/ha when the GeoSAR X-P interferometric height is considered along with the average LiDAR canopy height. Inclusion of a third LiDAR metric, the 60th decile height, further increased the R² to 85% and decreased the RMSE to 24.1 t/ha. On this 11 km² ponderosa pine study area, LiDAR data proved most useful for predicting forest biomass. RaDAR ranging measurements did not improve the LiDAR estimates.