Forests play a critical role in stabilizing Earth’s climate. Establishing protected areas (PAs) represents one approach to forest conservation, but PAs were rarely created to mitigate climate change. The global impact of PAs on the carbon cycle has not previously been quantified due to a lack of accurate global-scale carbon stock maps. Here we used ~412 million lidar samples from NASA’s GEDI mission to estimate a total PA aboveground carbon (C) stock of 61.43 Gt (+/− 0.31), 26% of all mapped terrestrial woody C. Of this total, 9.65 + /− 0.88 Gt of additional carbon was attributed to PA status. These higher C stocks are primarily from avoided emissions from deforestation and degradation in PAs compared to unprotected forests. This total is roughly equivalent to one year of annual global fossil fuel emissions. These results underscore the importance of conservation of high biomass forests for avoiding carbon emissions and preserving future sequestration.

Plant biodiversity supports life on Earth and provides a range of important ecosystem services, but is under severe pressure by global change. Structural diversity plays a crucial role for carbon, water and energy cycles and animal habitats. However, it is very difficult to map and monitor over large areas, limiting our ability to assess the status of biodiversity and predict change. NASA’s Global Ecosystem Dynamics Investigation (GEDI) provides a new opportunity to measure 3D plant canopy structure of the world’s temperate, Mediterranean and tropical ecosystems, but its potential to map structural diversity is not yet tested. Here, we use wall-to-wall airborne laser scanning (ALS) to simulate GEDI data (GEDI _sim ) over 7380 km² in the southern Sierra Nevada mountains in California and evaluate how well GEDI’s sampling scheme captures patterns of structural diversity. We evaluate functional richness and functional beta diversity in a biodiversity hot spot. GEDI _sim performed well for trait retrievals (r² = 0.68) and functional richness mapping (r² = 0.75) compared to ALS retrievals, despite lower correlations in complex terrain with steep slopes. Functional richness patterns were strongly associated with soil organic carbon stocks and density as well as variables related to water availability and could be appropriately mapped by GEDI _sim with and without cloud cover. Functional beta diversity was more strongly related to local changes in topography and more challenging to map, especially with decreasing sampling density. The reduced number of GEDI _sim shots when simulating cloud cover lead to a strong overestimation of beta diversity and a reduction of r² from 0.64 to 0.40 compared to ALS. The ability to map functional richness has been demonstrated with potential application at continental scales that could be transformative for our understanding of large-scale patterns of plant canopy structure, diversity and potential links to animal diversity, movement and habitats.

Forest structural complexity is a key element of ecosystem functioning, impacting light environments, nutrient cycling, biodiversity, and habitat quality. Addressing the need for a comprehensive global assessment of actual forest structural complexity, we derive a near-global map of 3D canopy complexity using data from the GEDI spaceborne lidar mission. These data show that tropical forests harbor most of the high complexity observations, while less than 20% of temperate forests reached median levels of tropical complexity. Structural complexity in tropical forests is more strongly related to canopy attributes from lower and middle waveform layers, whereas in temperate forests upper and middle layers are more influential. Globally, forests exhibit robust scaling relationships between complexity and canopy height, but these vary geographically and by biome. Our results offer insights into the spatial distribution of forest structural complexity and emphasize the importance of considering biome-specific and fine-scale variations for ecological research and management applications. The GEDI Waveform Structural Complexity Index data product, derived from our analyses, provides researchers and conservationists with a single, easily interpretable metric by combining various aspects of canopy structure.

Canopy structure, the vertical distribution of canopy material, is an important element of forest ecosystem dynamics and habitat preference. Although vertical stratification, or “canopy layering,” is a basic characterization of canopy structure for research and forest management, it is difficult to quantify at landscape scales. In this paper we describe canopy structure and develop methodologies to map forest vertical stratification in a mixed temperate forest using full-waveform lidar. Two definitions—one categorical and one continuous—are used to map canopy layering over Hubbard Brook Experimental Forest, New Hampshire with lidar data collected in 2009 by NASA’s Laser Vegetation Imaging Sensor (LVIS). The two resulting canopy layering datasets describe variation of canopy layering throughout the forest and show that layering varies with terrain elevation and canopy height. This information should provide increased understanding of vertical structure variability and aid habitat characterization and other forest management activities.

Long-term trends in photosynthetic capacity measured with the satellite-derived Normalized Difference Vegetation Index (NDVI) are usually associated with climate change. Human impacts on the global land surface are typically not accounted for. Here, we provide the first global analysis quantifying the effect of the earth’s human footprint on NDVI trends. Globally, more than 20% of the variability in NDVI trends was explained by anthropogenic factors such as land use, nitrogen fertilization, and irrigation. Intensely used land classes, such as villages, showed the greatest rates of increase in NDVI, more than twice than those of forests. These findings reveal that factors beyond climate influence global long-term trends in NDVI and suggest that global climate change models and analyses of primary productivity should incorporate land use effects.