Abiotic controls on macroscale variations of humid tropical forest height

Spatial variation of tropical forest tree height is a key indicator of ecological processes associated with forest growth and carbon dynamics. Here we examine the macroscale variations of tree height of humid tropical forests across three continents and quantify the climate and edaphic controls on these variations. Forest tree heights are systematically sampled across global humid tropical forests with more than 2.5 million measurements from Geoscience Laser Altimeter System (GLAS) satellite observations (2004–2008). We used top canopy height (TCH) of GLAS footprints to grid the statistical mean and variance and the 90 percentile height of samples at 0.5 degrees to capture the regional variability of average and large trees globally. We used the spatial regression method (spatial eigenvector mapping-SEVM) to evaluate the contributions of climate, soil and topography in explaining and predicting the regional variations of forest height. Statistical models suggest that climate, soil, topography, and spatial contextual information together can explain more than 60% of the observed forest height variation, while climate and soil jointly explain 30% of the height variations. Soil basics, including physical compositions such as clay and sand contents, chemical properties such as PH values and cation-exchange capacity, as well as biological variables such as the depth of organic matter, all present independent but statistically significant relationships to forest height across three continents. We found significant relations between the precipitation and tree height with shorter trees on the average in areas of higher annual water stress, and large trees occurring in areas with low stress and higher annual precipitation but with significant differences across the continents. Our results confirm other landscape and regional studies by showing that soil fertility, topography and climate may jointly control a significant variation of forest height and influencing patterns of aboveground biomass stocks and dynamics. Other factors such as biotic and disturbance regimes, not included in this study, may have less influence on regional variations but strongly mediate landscape and small-scale forest structure and dynamics.The research was funded by Gabon National Park (ANPN) under the contract of 011-ANPN/2012/SE-LJTW at UCLA. We thank IIASA, FAO, USGS, NASA, Worldclim science teams for making their data available. (011-ANPN/2012/SE-LJTW - Gabon National Park (ANPN) at UCLA

ICESat/GLAS Data as a Measurement Tool for Peatland Topography and Peat Swamp Forest Biomass in Kalimantan, Indonesia

Indonesian peatlands are one of the largest near-surface pools of terrestrial organic carbon. Persistent logging, drainage and recurrent fires lead to huge emission of carbon each year. Since tropical peatlands are highly inaccessible, few measurements on peat depth and forest biomass are available. We assessed the applicability of quality filtered ICESat/GLAS (a spaceborne LiDAR system) data to measure peatland topography as a proxy for peat volume and to estimate peat swamp forest Above Ground Biomass (AGB) in a thoroughly investigated study site in Central Kalimantan, Indonesia. Mean Shuttle Radar Topography Mission (SRTM) elevation was correlated to the corresponding ICESat/GLAS elevation. The best results were obtained from the waveform centroid (R2 = 0.92; n = 4,186). ICESat/GLAS terrain elevation was correlated to three 3D peatland elevation models derived from SRTM data (R2 = 0.90; overall difference = −1.0 m, ±3.2 m; n = 4,045). Based on the correlation of in situ peat swamp forest AGB and airborne LiDAR data (R2 = 0.75, n = 36) an ICESat/GLAS AGB prediction model was developed (R2 = 0.61, n = 35). These results demonstrate that ICESat/GLAS data can be used to measure peat topography and to collect large numbers of forest biomass samples in remote and highly inaccessible peatland forests

Geographic variability in lidar predictions of forest stand structure in the Pacific Northwest

Estimation of the amount of carbon stored in forests is a key challenge for understanding the global carbon cycle, one which remote sensing is expected to help address. However, carbon storage in moderate to high biomass forests is difficult to estimate with conventional optical or radar sensors. Lidar (light detection and ranging) instruments measure the vertical structure of forests and thus hold great promise for remotely sensing the quantity and spatial organization of forest biomass. In this study, we compare the relationships between lidar measured canopy structure and coincident field measurements of forest stand structure at five locations in the Pacific Northwest of the U.S.A. with contrasting composition. Coefficient of determination values (r2) ranged between 41% and 96%. Correlations for two important variables, LAI (81%) and above ground biomass (92%), were noteworthy, as was the fact that neither variable showed an asymptotic response. Of the 17 stand structure variables considered in this study, we were able to develop eight equations that were valid for all sites, including equations for two variables generally considered to be highly important (aboveground biomass and leaf area index). The other six equations that were valid for all sites were either related to height (which is most directly measured by lidar) or diameter at breast height (which should be closely related to height). Four additional equations (a total of 12) were applicable to all sites where either Douglas-fir (Pseudotsuga menziesii), western hemlock (Tsuga heterophylla) or Sitka spruce (Picea sitchensi) were dominant. Stand structure variables in sites dominated by true firs (Abies sp.) or ponderosa pine (Pinus ponderosa) had biases when predicted by these four additional equations. Productivity-related variables describing the edaphic, climatic and topographic environment of the sites where available for every regression, but only two of the 17 equations (maximum diameter at breast height, stem density) incorporated them. Given the wide range of these environmental conditions sampled, we conclude that the prediction of stand structure is largely independent of environmental conditions in this study area. Most studies of lidar remote sensing for predicting stand structure have depended on intensive data collections within a relatively small study area. This study indicates that the relationships between many stand structure indices and lidar measured canopy structure have generality at the regional scale. This finding, if replicated in other regions, would suggest that mapping of stand structure using lidar may be accomplished by distributing field sites extensively over a region, thus reducing the overall inventory effort required

Observing the Forest Canopy with a New Ultra-Violet Compact Airborne Lidar

We have developed a new airborne UV lidar for the forest canopy and deployed it in the Landes forest (France). It is the first one that: (i) operates at 355 nm for emitting energetic pulses of 16 mJ at 20 Hz while fulfilling eye-safety regulations and (ii) is flown onboard an ultra-light airplane for enhanced flight flexibility. Laser footprints at ground level were 2.4 m wide for a flying altitude of 300 m. Three test areas of ∼500 × 500 m2 with Maritime pines of different ages were investigated. We used a threshold method adapted for this lidar to accurately extract from its waveforms detailed forest canopy vertical structure: canopy top, tree crown base and undergrowth heights. Good detection sensitivity enabled the observation of ground returns underneath the trees. Statistical and one-to-one comparisons with ground measurements by field foresters indicated a mean absolute accuracy of ∼1 m. Sensitivity tests on detection threshold showed the importance of signal to noise ratio and footprint size for a proper detection of the canopy vertical structure. This UV-lidar is intended for future innovative applications of simultaneous observation of forest canopy, laser-induced vegetation fluorescence and atmospheric aerosols

Validation of Canopy Height Profile methodology for small-footprint full-waveform airborne LiDAR data in a discontinuous canopy environment

A Canopy Height Profile (CHP) procedure presented in Harding et al. (2001) for large footprint LiDAR data was tested in a closed canopy environment as a way of extracting vertical foliage profiles from LiDAR raw-waveform. In this study, an adaptation of this method to small-footprint data has been shown, tested and validated in an Australian sparse canopy forest at plot- and site-level. Further, the methodology itself has been enhanced by implementing a dataset-adjusted reflectance ratio calculation according to Armston et al. (2013) in the processing chain, and tested against a fixed ratio of 0.5 estimated for the laser wavelength of 1550nm. As a by-product of the methodology, effective leaf area index (LAIe) estimates were derived and compared to hemispherical photography-derived values. To assess the influence of LiDAR aggregation area size on the estimates in a sparse canopy environment, LiDAR CHPs and LAIes were generated by aggregating waveforms to plot- and site-level footprints (plot/site-aggregated) as well as in 5m grids (grid-processed). LiDAR profiles were then compared to leaf biomass field profiles generated based on field tree measurements. The correlation between field and LiDAR profiles was very high, with a mean R2 of 0.75 at plot-level and 0.86 at site-level for 55 plots and the corresponding 11 sites. Gridding had almost no impact on the correlation between LiDAR and field profiles (only marginally improvement), nor did the dataset-adjusted reflectance ratio. However, gridding and the dataset-adjusted reflectance ratio were found to improve the correlation between raw-waveform LiDAR and hemispherical photography LAIe estimates, yielding the highest correlations of 0.61 at plot-level and of 0.83 at site-level. This proved the validity of the approach and superiority of dataset-adjusted reflectance ratio of Armston et al. (2013) over a fixed ratio of 0.5 for LAIe estimation, as well as showed the adequacy of small-footprint LiDAR data for LAIe estimation in discontinuous canopy forests

Estimating Canopy Gap Fraction Using ICESat GLAS within Australian Forest Ecosystems

Spaceborne laser altimetry waveform estimates of canopy Gap Fraction (GF) vary withrespect to discrete return airborne equivalents due to their greater sensitivity to reflectance differencesbetween canopy and ground surfaces resulting from differences in footprint size, energy thresholding,noise characteristics and sampling geometry. Applying scaling factors to either the ground or canopyportions of waveforms has successfully circumvented this issue, but not at large scales. This studydevelops a method to scale spaceborne altimeter waveforms by identifying which remotely-sensedvegetation, terrain and environmental attributes are best suited to predicting scaling factors basedon an independent measure of importance. The most important attributes were identified as: soilphosphorus and nitrogen contents, vegetation height, MODIS vegetation continuous fields productand terrain slope. Unscaled and scaled estimates of GF are compared to corresponding ALS datafor all available data and an optimized subset, where the latter produced most encouraging results(R2 = 0.89, RMSE = 0.10). This methodology shows potential for successfully refining estimates ofGF at large scales and identifies the most suitable attributes for deriving appropriate scaling factors.Large-scale active sensor estimates of GF can establish a baseline from which future monitoringinvestigations can be initiated via upcoming Earth Observation missions

A new concept for estimating the influence of vegetation on throughfall kinetic energy using aerial laser scanning

Soil loss caused by erosion has enormous economic and social impacts. Splash effects of rainfall are an important driver of erosion processes; however, effects of vegetation on splash erosion are still not fully understood. Splash erosion processes under vegetation are investigated by means of throughfall kinetic energy (TKE). Previous studies on TKE utilized a heterogeneous set of plant and canopy parameters to assess vegetation’s influence on erosion by rain splash but remained on individual plant- or plotlevels. In the present study we developed a method for the area-wide estimation of the influence of vegetation on TKE using remote sensing methods. In a literature review we identified key vegetation variables influencing splash erosion and developed a conceptual model to describe the interaction of vegetation and raindrops. Our model considers both amplifying and protecting effect of vegetation layers according to their height above the ground and aggregates them into a new indicator: the Vegetation Splash Factor (VSF). It is based on the proportional contribution of drips per layer, which can be calculated via the vegetation cover profile from airborne LiDAR datasets. In a case study, we calculated the VSF using a LiDAR dataset for La Campana National Park in central Chile. The studied catchment comprises a heterogeneous mosaic of vegetation layer combinations and types and is hence well suited to test the approach.We calculated a VSF map showing the relation between vegetation structure and its expected influence on TKE. Mean VSF was 1.42, indicating amplifying overall effect of vegetation on TKE that was present in 81% of the area. Values below 1 indicating a protective effect were calculated for 19% of the area. For future work, we recommend refining the weighting factor by calibration to local conditions using field-reference data and comparing the VSF with TKE field measurements

Quantifying Dynamics in Tropical Peat Swamp Forest Biomass with Multi- Temporal LiDAR Datasets

Tropical peat swamp forests in Indonesia store huge amounts of carbon and are responsible for enormous carbon emissions every year due to forest degradation and deforestation. These forest areas are in the focus of REDD+ (reducing emissions from deforestation, forest degradation, and the role of conservation, sustainable management of forests and enhancement of forest carbon stocks) projects, which require an accurate monitoring of their carbon stocks or aboveground biomass (AGB). Our study objective was to evaluate multi-temporal LiDAR measurements of a tropical forested peatland area in Central Kalimantan on Borneo. Canopy height and AGB dynamics were quantified with a special focus on unaffected, selective logged and burned forests. More than 11,000 ha were surveyed with airborne LiDAR in 2007 and 2011. In a first step, the comparability of these datasets was examined and canopy height models were created. Novel AGB regression models were developed on the basis of field inventory measurements and LiDAR derived height histograms for 2007 (r(2) = 0.77, n = 79) and 2011 (r(2) = 0.81, n = 53), taking the different point densities into account. Changes in peat swamp forests were identified by analyzing multispectral imagery. Unaffected forests accumulated on average 20 t/ha AGB with a canopy height increase of 2.3 m over the four year time period. Selective logged forests experienced an average AGB loss of 55 t/ha within 30 m and 42 t/ha within 50 m of detected logging trails, although the mean canopy height increased by 0.5 m and 1.0 m, respectively. Burned forests lost 92% of the initial biomass. These results demonstrate the great potential of repetitive airborne LiDAR surveys to precisely quantify even small scale AGB and canopy height dynamics in remote tropical forests, thereby featuring the needs of REDD+

Estimates of forest canopy height and aboveground biomass using ICESat

Exchange of carbon between forests and the atmosphere is a vital component of the global carbon cycle. Satellite laser altimetry has a unique capability for estimating forest canopy height, which has a direct and increasingly well understood relationship to aboveground carbon storage. While the Geoscience Laser Altimeter System (GLAS) onboard the Ice, Cloud and land Elevation Satellite (ICESat) has collected an unparalleled dataset of lidar waveforms over terrestrial targets, processing of ICESat data to estimate forest height is complicated by the pulse broadening associated with large-footprint, waveform-sampling lidar. We combined ICESat waveforms and ancillary topography from the Shuttle Radar Topography Mission to estimate maximum forest height in three ecosystems; tropical broadleaf forests in Brazil, temperate broadleaf forests in Tennessee, and temperate needleleaf forests in Oregon. Final models for each site explained between 59% and 68% of variance in field-measured forest canopy height (RMSE between 4.85 and 12.66 m). In addition, ICESat-derived heights for the Brazilian plots were correlated with field-estimates of aboveground biomass (r(2) = 73%, RMSE = 58.3 Mgha(-1))