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

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

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))

Characterizing 3D Vegetation Structure from Space: Mission Requirements

Human and natural forces are rapidly modifying the global distribution and structure of terrestrial ecosystems on which all of life depends, altering the global carbon cycle, affecting our climate now and for the foreseeable future, causing steep reductions in species diversity, and endangering Earth s sustainability. To understand changes and trends in terrestrial ecosystems and their functioning as carbon sources and sinks, and to characterize the impact of their changes on climate, habitat and biodiversity, new space assets are urgently needed to produce high spatial resolution global maps of the three-dimensional (3D) structure of vegetation, its biomass above ground, the carbon stored within and the implications for atmospheric green house gas concentrations and climate. These needs were articulated in a 2007 National Research Council (NRC) report (NRC, 2007) recommending a new satellite mission, DESDynI, carrying an L-band Polarized Synthetic Aperture Radar (Pol-SAR) and a multi-beam lidar (Light RAnging And Detection) operating at 1064 nm. The objectives of this paper are to articulate the importance of these new, multi-year, 3D vegetation structure and biomass measurements, to briefly review the feasibility of radar and lidar remote sensing technology to meet these requirements, to define the data products and measurement requirements, and to consider implications of mission durations. The paper addresses these objectives by synthesizing research results and other input from a broad community of terrestrial ecology, carbon cycle, and remote sensing scientists and working groups. We conclude that: (1) current global biomass and 3-D vegetation structure information is unsuitable for both science and management and policy. The only existing global datasets of biomass are approximations based on combining land cover type and representative carbon values, instead of measurements of actual biomass. Current measurement attempts based on radar and multispectral data have low explanatory power outside low biomass areas. There is no current capability for repeatable disturbance and regrowth estimates. (2) The science and policy needs for information on vegetation 3D structure can be successfully addressed by a mission capable of producing (i) a first global inventory of forest biomass with a spatial resolution 1km or finer and unprecedented accuracy (ii) annual global disturbance maps at a spatial resolution of 1 ha with subsequent biomass accumulation rates at resolutions of 1km or finer, and (iii) transects of vertical and horizontal forest structure with 30 m along-transect measurements globally at 25 m spatial resolution, essential for habitat characterization. We also show from the literature that lidar profile samples together with wall-to53 wall L-band quad-pol-SAR imagery and ecosystem dynamics models can work together to satisfy these vegetation 3D structure and biomass measurement requirements. Finally we argue that the technology readiness levels of combined pol-SAR and lidar instruments are adequate for space flight. Remaining to be worked out, are the particulars of a lidar/pol-SAR mission design that is feasible and at a minimum satisfies the information and measurement requirement articulated herein

Interannual variability of carbon uptake of secondary forests in the Brazilian Amazon (2004‐2014)

Tropical secondary forests (SF) play an important role in the global carbon cycle as a major terrestrial carbon sink. Here, we use high-resolution TerraClass data set for tracking land use activities in the Brazilian Amazon from 2004–2014 to detect spatial patterns and carbon sequestration dynamics of secondary forests (SF). By integrating satellite lidar and radar observations, we found the SF area in the Brazilian Amazon increased from approximately 22 Mha (10^6 ha) in 2004 to 28 Mha in 2014. However, the expansion in area was also accompanied by a dynamic land use activity that resulted in about 50% recycling of SF area annually from frequent clearing and abandonment. Consequently, the average age of SF remained less than 10 years (age ~8.2 with one standard deviation of 3.2 spatially) over the period of the study. Estimation of changes of carbon stocks shows that SF accumulates approximately 8.5 Mg ha^−1 year^−1 aboveground biomass during the first 10 years after clearing and abandonment, 4.5 Mg ha^−1 year^−1 for the next 10 years followed by a more gradual increase of 3 Mg ha^−1 year^−1 from 20 to 30 years with much slower rate thereafter. The effective carbon uptake of SF in Brazilian Amazon was negligible (0.06 ± 0.03 PgC year^−1) during this period, but the interannual variability was significantly larger (±0.2 PgC year^−1). If the SF areas were left to grow without further clearing for 100 years, it would absorb about 0.14 PgC year^−1 from the atmosphere, partially compensating the emissions from current rate of deforestation in the Brazilian Amazon.Published versio

Patterns of covariance between forest stand and canopy structure in the Pacific Northwest

In the past decade, lidar (light detection and ranging) has emerged as a powerful tool for remotely sensing forest canopy and stand structure, including the estimation of aboveground biomass and carbon storage. Numerous papers have documented the use of lidar measurements to predict important aspects of forest stand structure, including aboveground biomass. Other papers have documented the ability to transform lidar measurements to approximate common field measures, such as cover, stand height, and vertical distributions of foliage density and light transmittance. However, only a small number of existing works have thoroughly examined relationships between comprehensive assemblages of forest canopy and forest stand structure indices. In this work, canonical correlation analysis of coincident lidar and field datasets in western Oregon and Washington is used to define seven statistically significant pairs of canonical variables, each defining an axis of variation that stand and canopy structure have in common. The first major axis relates mean stand height, and related variables, to aboveground biomass. The second relates canopy cover and volume to leaf area index and stem density. The third relates canopy height variability to mean stem diameter and the basal area of deciduous species. Of the four remaining axes, three are related to contrasts between mature and old-growth stands. Canonical correlation analysis provides a method for ranking the importance of these effects, and for placing both canopy and stand structure indices within the overall covariance structure of the two datasets. In this sense, and for the study area involved, the first three factors (mean height, cover or leaf index area, height variability) represent the same kind of enhancement of lidar data that the tasseled cap indices [Crist, C.P., R.C. Cicone, 1984. A physically-based transformation of thematic mapper data—the TM tasseled cap. IEEE Transactions on Geoscience and Remote Sensing 22, 256–263.] represent for optical remote sensing

Predictions of Tropical Forest Biomass and Biomass Growth Based on Stand Height or Canopy Area Are Improved by Landsat-Scale Phenology across Puerto Rico and the U.S. Virgin Islands

Remotely-sensed estimates of forest biomass are usually based on various measurements of canopy height, area, volume or texture, as derived from LiDAR, radar or fine spatial resolution imagery. These measurements are then calibrated to estimates of stand biomass that are primarily based on tree stem diameters. Although humid tropical forest seasonality can have low amplitudes compared with temperate regions, seasonal variations in growth-related factors like temperature, humidity, rainfall, wind speed and day length affect both tropical forest deciduousness and tree height-diameter relationships. Consequently, seasonal patterns in spectral measures of canopy greenness derived from satellite imagery should be related to tree height-diameter relationships and hence to estimates of forest biomass or biomass growth that are based on stand height or canopy area. In this study, we tested whether satellite image-based measures of tropical forest phenology, as estimated by indices of seasonal patterns in canopy greenness constructed from Landsat satellite images, can explain the variability in forest deciduousness, forest biomass and net biomass growth after already accounting for stand height or canopy area. Models of forest biomass that added phenology variables to structural variables similar to those that can be estimated by LiDAR or very high-spatial resolution imagery, like canopy height and crown area, explained up to 12% more variation in biomass. Adding phenology to structural variables explained up to 25% more variation in Net Biomass Growth (NBG). In all of the models, phenology contributed more as interaction terms than as single-effect terms. In addition, models run on only fully-forested plots performed better than models that included partially-forested plots. For forest NBG, the models produced better results when only those plots with a positive growth, from Inventory Cycle 1 to Inventory Cycle 2, were analyzed, as compared to models that included all plot

Forest carbon densities and uncertainties from Lidar, QuickBird, and field measurements in California

Greenhouse gas inventories and emissions reduction programs require robust methods to quantify carbon sequestration in forests. We compare forest carbon estimates from Light Detection and Ranging (Lidar) data and QuickBird high-resolution satellite images, calibrated and validated by field measurements of individual trees. We conducted the tests at two sites in California: (1) 59 km2 of secondary and old-growth coast redwood (Sequoia sempervirens) forest (Garcia–Mailliard area) and (2) 58 km2 of old-growth Sierra Nevada forest (North Yuba area). Regression of aboveground live tree carbon density, calculated from field measurements, against Lidar height metrics and against QuickBird-derived tree crown diameter generated equations of carbon density as a function of the remote sensing parameters. Employing Monte Carlo methods, we quantified uncertainties of forest carbon estimates from uncertainties in field measurements, remote sensing accuracy, biomass regression equations, and spatial autocorrelation. Validation of QuickBird crown diameters against field measurements of the same trees showed significant correlation (r = 0.82, P \u3c 0.05). Comparison of stand-level Lidar height metrics with field-derived Lorey\u27s mean height showed significant correlation (Garcia–Mailliard r = 0.94, P \u3c 0.0001; North Yuba R = 0.89, P \u3c 0.0001). Field measurements of five aboveground carbon pools (live trees, dead trees, shrubs, coarse woody debris, and litter) yielded aboveground carbon densities (mean ± standard error without Monte Carlo) as high as 320 ± 35 Mg ha− 1 (old-growth coast redwood) and 510 ± 120 Mg ha− 1 (red fir [Abies magnifica] forest), as great or greater than tropical rainforest. Lidar and QuickBird detected aboveground carbon in live trees, 70–97% of the total. Large sample sizes in the Monte Carlo analyses of remote sensing data generated low estimates of uncertainty. Lidar showed lower uncertainty and higher accuracy than QuickBird, due to high correlation of biomass to height and undercounting of trees by the crown detection algorithm. Lidar achieved uncertainties of \u3c 1%, providing estimates of aboveground live tree carbon density (mean ± 95% confidence interval with Monte Carlo) of 82 ± 0.7 Mg ha− 1 in Garcia–Mailliard and 140 ± 0.9 Mg ha− 1 in North Yuba. The method that we tested, combining field measurements, Lidar, and Monte Carlo, can produce robust wall-to-wall spatial data on forest carbon