Both leaf inclination angle distribution (LAD) and leaf area index (LAI) dominate optical remote sensing signals. The G-function, which is a function of LAD and remote sensing geometry, is often set to 0.5 in the LAI retrieval of coniferous canopies even though this assumption is only valid for spherical LAD. Large uncertainties are thus introduced. However, because numerous tiny leaves grow on conifers, it is nearly impossible to quantitatively evaluate such uncertainties in LAI retrieval. In this study, we proposed a method to characterize the possible change of G-function of coniferous canopies as well as its effect on LAI retrieval. Specifically, a Multi-Directional Imager (MDI) was developed to capture stereo images of the branches, and the needles were reconstructed. The accuracy of the inclination angles calculated from the reconstructed needles was high. Moreover, we analyzed whether a spherical distribution is a valid assumption for coniferous canopies by calculating the possible range of the G-function from the measured LADs of branches of Larch and Spruce and the true G-functions of other species from some existing inventory data and three-dimensional (3D) tree models. Results show that the constant G assumption introduces large errors in LAI retrieval, which could be as large as 53% in the zenithal viewing direction used by spaceborne LiDAR. As a result, accurate LAD estimation is recommended. In the absence of such data, our results show that a viewing zenith angle between 45 and 65 degrees is a good choice, at which the errors of LAI retrieval caused by the spherical assumption will be less than 10% for coniferous canopies.

Clumping effect denotes the nonrandomness of foliage. It deviates from the random distribution assumption of Beer's law which is usually applied to leaf area index (LAI) retrieval from large-footprint full-waveform light detection and ranging (LiDAR). Some studies correct for large gaps-induced between-crown clumping, yet ignore the within-crown clumping. The error of LAI caused by these clumping effects and the influence of the forest structure parameters on them have not been quantitatively studied. This study quantified the between-crown, within-crown, and total clumping indices through a theoretical derivation, clarifying the mechanism of clumping; we used airborne LiDAR point clouds data in 11 290 footprints (diameter = 25 m) to estimate these indices in real forests. We found that: 1) the underestimation of LAI caused by directly applying Beer's law could be up to 93%, and it decreases with fractional crown coverage but increases with crown length and leaf area density; 2) the method of correcting between-crown clumping improves LAI retrieval for cylindrical canopies effectively; however, 3) considerable underestimation (up to 58%) exists if we neglect the within-crown clumping for other canopies, which has not been realized before; and 4) both the between-crown and the within-crown clumping can be the dominant contributor, and the within-crown clumping was greater than the between-crown clumping in 47% of the studied footprints. In the two physically based LAI retrieval methods, Beer's law has been commonly used due to its simplicity. Pathways to improve future LAI retrieval would be instrument improvement to capture the between-crown gaps and method study to correct the within-crown clumping further.

The demand for leaf area index (LAI) retrieval from spaceborne full-waveform LiDAR increases due to its direct sampling of the three-dimensional forest structure at a near-global scale. However, the nonrandomness (i.e., clumping effect) of canopy composition limits the reliability of LAI derived from two common methods. They either assume a homogeneous scene in the footprint or just correct for the large gaps-induced between-crown clumping. The clumping in the crown is still an unaddressed issue. We proposed a method to compensate occlusion (i.e., lower canopy layers are occluded by the upper canopy in the process of LiDAR measurement), through which the vertical canopy profile can be resolved from the waveform. Further, we developed a method of deriving relative path length distribution that can reflect the heterogeneity of the canopy from the occlusion-corrected waveform. In addition to correcting the between-crown clumping, we corrected the within-crown clumping further using the derived relative path length distribution, based on path length distribution (PATH) theory. We used simulated waveform data with known LAI and GLAS data with corresponding field-measured LAI to test the performance of our and the other two common LAI retrieval methods. Results show that the errors of our approach are the lowest (with an error generally below 10% and the maximum error below 20%, compared with up to 69% and 47% for the other two methods), and it is relatively stable in various scenes. This study demonstrated the potential of improving LAI retrieval through full utilization of full-waveform data.

Sloping terrain distorts the sun-target–sensor geometry, resulting in biases of the optical reflectance measured by remote sensors relative to flat situations. Performing topographic correction (TC) is, therefore, deemed mandatory to foster the full exploitation of satellite images worldwide to support various applications in mountainous regions. Various TC models have already been proposed and developed, while most of them were previously evaluated at local or regional scales using a few images with various evaluation criteria. Therefore, a systematic and comprehensive assessment has yet to be done on these TC models in the global mountainous regions. In the present study, 10 523 Landsat 8 Operational Land Imager (OLI) images, filtered by land cover types and seasons and sampled across global mountain regions, were corrected using ten popular TC models—statistical empirical (SE), b correction, variable empirical coefficient algorithm (VECA), C-correction (CC), sun-canopy–sensor (SCS), DymondShepherd (DS), SCS + C, path length correction (PLC), Minnaert, and Minnaert + SCS—with a unified evaluation criterion on the Google Earth Engine (GEE) platform. The outcomes are that: 1) global TC effects on Landsat 8 OLI images generally increase with sun zenith angles (SZAs) and latitudes; 2) six models (SE, b correction, CC, VECA, Minnaert, and Minnaert + SCS) show good adaptability among the ten models for the global mountainous placing a disregard to land cover types and seasons; and 3) considering permanent snow and ice, needle-leaved forests in winter, and null values might appear in b correction, SE is deemed to be with the most global adaptability. This study pioneers an evaluation of fashionable TC models concerning mountainous regions worldwide and will be useful for applying TC to Landsat images for the benefit of making global TC products in the future and a fair intercomparison of OLI surface reflectance measured in various mountainous areas of the globe.

Phenological changes in forests directly influence the spatiotemporal dynamics of carbon fixation and the carbon and water cycles in terrestrial ecosystems. Previous studies have shown that variations in biological factors (e.g., canopy height, leaf area, water use efficiency) can increase uncertainty in forest phenology, and these variations are closely linked to tree species and forest age-class structure. However, the interaction mechanisms through which tree species and forest age-class structure influence phenological changes remain insufficiently explored. In this study, phenological changes and their interactions and response mechanisms to different dominant tree species and forest age-class structures were analyzed via Sentinel-2 normalized difference vegetation index (NDVI) time series data from 2020 and 2021 across 480 typical deciduous forest plots in northeastern China. The results were as follows: (1) There were significant differences (p < 0.05) in the intra-annual phenological responses of temperate deciduous landscapes to the interaction between tree species and forest age-class structure. (2) The indirect effect of forest age-class structure through tree species on phenology exceeded the indirect effect of tree species through forest age-class structure, with a difference of 30.77%–35.09%. (3) When the dominant tree species and forest age-class structure were not distinguished, phenological differences in temperate forests ranged from 3 to 41 days and 2 to 23 days, respectively. This study highlights the differential impacts of key biological factors and their interactions on regional forest phenology, offering valuable insights into how these factors influence forest landscapes and providing a theoretical basis for improving forest management strategies.

Information about the composition and distribution of dominant tree species is crucial for sustainable forest management. A global ecosystem dynamics investigation (GEDI) offers unique advantages in detecting the vertical spatial structure of forest stands, which may improve the common issues of spectral similarity and saturation in traditional spectral-based tree species mapping. However, the effects of its application have not been explored. This study examines temperate and subtropical forests in eastern China, which are dominated by deciduous and evergreen species, respectively. We employed GEDI-derived forest vertical structure characterization (FVSC) to complement Sentinel-2 spectral features for dominant tree species mapping. The results indicate that FVSC significantly improved the mapping accuracy for 11 tree species in both temperate and subtropical forest regions across seasons, with greater benefits observed for broadleaf species than for coniferous species. During the main phenological stages of spring, summer, and autumn, the accuracy of tree species mapping in the temperate region improved by 6.99-9.85%. The key contributing factors were the cumulative plant area index (PAI) from the ground to the canopy top and cumulative vegetation coverage (COVER) from 5 m to the canopy top. In the subtropical region, the accuracy improvement ranged from 7.75-9.5%, with the highest contributions from the plant gap probability (Pgap_theta) and cumulative COVER from 5 m to the canopy top. These findings demonstrate that FVSC can effectively support spectral feature data in the detailed mapping of dominant tree species at regional scales. Moreover, the method shows good stability and applicability across seasons and climatic regions.