Canopy layers form essential structural components, affecting stand productivity and wildlife habitats. Airborne laser scanning (ALS) provides horizontal and vertical information on canopy structure simultaneously. Existing approaches to assess canopy layering often require prior information about stand characteristics or rely on pre-defined height thresholds. We developed a multi-scale method using ALS data with point densities >10 pts/m² to determine the number and vertical extent of canopy layers (canopy_layer, canopy_length), seasonal variations in the topmost canopy layer (canopy_type), as well as small-scale heterogeneities in the canopy (canopy_{heterogeneity}). We first tested and developed the method on a small forest patch (800 ha) and afterwards tested transferability and robustness of the method on a larger patch (180,000 ha). We validated the approach using an extensive set of ground data, achieving overall accuracies >77% for canopy_type and canopy_{heterogeneity}, and >62% for canopy_layer and canopy_length. We conclude that our method provides a robust characterization of canopy layering supporting automated canopy structure monitoring.

We outline an approach combining airborne laser scanning (ALS) and imaging spectroscopy (IS) to quantify and assess patterns of tree density (TD) and forest productivity (FP) in a protected heterogeneous alpine forest in the Swiss National Park (SNP). We use ALS data and a local maxima (LM) approach to predict TD, as well as IS data (Airborne Prism Experiment—APEX) and an empirical model to estimate FP. We investigate the dependency of TD and FP on site related factors, in particular on surface exposition and elevation. Based on reference data (i.e., 1598 trees measured in 35 field plots), we observed an underestimation of ALS-based TD estimates of 40%. Our results suggest a limited sensitivity of the ALS approach to small trees as well as a dependency of TD estimates on canopy heterogeneity, structure, and species composition. We found a weak to moderate relationship between surface elevation and TD (R² = 0.18–0.69) and a less pronounced trend with FP (R² = 0.0–0.56), suggesting that both variables depend on gradients of resource availability. Further to the limitations faced in the sensitivity of the applied approaches, we conclude that the combined application of ALS and IS data was convenient for estimating tree density and mapping FP in north-facing forested areas, however, the accuracy was lower in south-facing forested areas covered with multi-stemmed trees.