Search

Understanding the heterogeneity of biomass accumulation in second-growth tropical forests following land use abandonment is important for informing ecosystem carbon models and forest restoration efforts. There is an urgent need for a broad sample of second-growth forests to enhance our knowledge of carbon accumulation in human-dominated landscapes, especially for older forests. Puerto Rico has predominantly second-growth forests, ranging in age from approximately 25 to more than 80 years. We used an island-wide sample of airborne lidar from the NASA Goddard Lidar, Hyperspectral, and Thermal (G-LiHT) Airborne Imager collected on March 2017, forest inventory data, and data on forest age, precipitation, soils, and land use to estimate aboveground biomass stocks in moist and wet, second-growth tropical forests. Biomass accumulation rates in Puerto Rico were lower, on average, than in other Neotropical forests. Median biomass across >16,700 ha of older second-growth forests was 105 Mg ha−1, and sampled biomass rarely surpassed 250 Mg ha−1. Differences in biomass by age were large and persistent across different substrates and land uses, with a plateau in the pattern of island-wide biomass accumulation after about 33 years. A spatial regression model showed that multiple factors were related to biomass accumulation, including time since abandonment, geologic substrate, past land use as coffee or pasture, precipitation, topographic wetness index, and slope. Our findings have important consequences for the total carbon storage and expected climate mitigation benefits of large-scale reforestation efforts, and highlight the value of airborne lidar for quantifying biomass variability in complex tropical landscapes.

Hurricane Maria, a Category 4 storm, snapped and uprooted canopy trees, removed large branches, and defoliated vegetation across Puerto Rico. The magnitude of forest damages and the rates and mechanisms of forest recovery following Maria provide important benchmarks for understanding the ecology of extreme events. We used airborne Lidar data acquired before (2017) and after Maria (2018, 2020) to quantify landscape-scale changes in forest structure along a 439-ha elevational gradient (100–800 m) in the Luquillo Experimental Forest. Damages from Maria were widespread, with 73% of the study area losing ≥ 1 m in canopy height (mean = −7.1 m). Taller forests at lower elevations suffered more damage than shorter forests above 600 m. Yet only 13.5% of the study area had canopy heights ≤ 2 m in 2018, a typical threshold for forest gaps, highlighting the importance of damaged trees and advanced regeneration on post-storm forest structure. Heterogeneous patterns of regrowth and recruitment yielded shorter and more open forests by 2020. Nearly 45% of forests experienced initial height loss > 1 m (2017–2018) followed by rapid height gain > 1 m (2018–2020), whereas 21.6% of forests with initial height losses showed little or no height gain, and 17.8% of forests exhibited no height changes larger than ± 1 m in either period. Canopy layers < 10 m tall accounted for most increases in canopy height and fractional cover between 2018 and 2020, with gains split evenly between height growth and lateral crown expansion by surviving individuals. These findings benchmark rates of gap formation, crown expansion, and canopy closure following hurricane damage and highlight the diversity of ecosystem impacts from heterogeneous spatial patterns and vertical stratification of forest regrowth following a major disturbance event.

Coarse dead wood is an important component of forest carbon stocks, but it is rarely measured in Amazon forests and is typically excluded from regional forest carbon budgets. Our study is based on line intercept sampling for fallen coarse dead wood conducted along 103 transects with a total length of 48 km matched with forest inventory plots where standing coarse dead wood was measured in the footprints of larger areas of airborne lidar acquisitions. We developed models to relate lidar metrics and Landsat time series variables to coarse dead wood stocks for intact, logged, burned, or logged and burned forests. Canopy characteristics such as gap area produced significant individual relations for logged forests. For total fallen plus standing coarse dead wood (hereafter defined as total coarse dead wood), the relative root mean square error for models with only lidar metrics ranged from 33 % in logged forest to up to 36 % in burned forests. The addition of historical information improved model performance slightly for intact forests (31 % against 35 % relative root mean square error), not justifying the use of a number of disturbance events from historical satellite images (Landsat) with airborne lidar data. Lidar-derived estimates of total coarse dead wood compared favorably with independent ground-based sampling for areas up to several hundred hectares. The relations found between total coarse dead wood and variables quantifying forest structure derived from airborne lidar highlight the opportunity to quantify this important but rarely measured component of forest carbon over large areas in tropical forests.