Surface albedo plays a controlling role in the surface energy budget, and albedo-induced radiative forcing has a significant impact on climate and environmental change (e.g., global warming, snow and ice melt, soil and vegetation degradation, and urban heat islands (UHIs)). Several existing review papers have summarized the algorithms and products of surface albedo as well as climate feedback at certain surfaces, while an overall understanding of various land types remains insufficient, especially with increasing studies on albedo management methods regarding mitigating global warming in recent years. In this paper, we present a comprehensive literature review on the variance pattern of surface albedo, the subsequent climate impact, and albedo management strategies. The results show that using the more specific term “surface albedo” is recommended instead of “albedo” to avoid confusion with similar terms (e.g., planetary albedo), and spatiotemporal changes in surface albedo can indicate subtle changes in the energy budget, land cover, and even the specific surface structure. In addition, the close relationships between surface albedo change and climate feedback emphasize the important role of albedo in climate simulation and forecasting, and many albedo management strategies (e.g., the use of retroreflective materials (RRMs)) have been demonstrated to be effective for climate mitigation by offsetting CO₂ emissions. In future work, climate effects and management strategies regarding surface albedo at a multitude of spatiotemporal resolutions need to be systematically evaluated to promote its application in climate mitigation, where a life cycle assessment (LCA) method considering both climate benefits and side effects (e.g., thermal comfort) should be followed.

Methods that link different models for investigating the retrieval of canopy biophysical/structural variables have been substantially adopted in the remote sensing community. To retrieve global biophysical parameters from multiangle data, the kernel-driven bidirectional reflectance distribution function (BRDF) model has been widely applied to satellite multiangle observations to model (interpolate/extrapolate) the bidirectional reflectance factor (BRF) in an arbitrary direction of viewing and solar geometries. Such modeled BRFs, as an essential information source, are then input into an inversion procedure that is devised through a large number of simulation analyses from some widely used physical models that can generalize such an inversion relationship between the BRFs (or their simple algebraic composite) and the biophysical/structural parameter. Therefore, evaluation of such a link between physical models and kernel-driven models contributes to the development of such inversion procedures to accurately retrieve vegetation properties, particularly based on the operational global BRDF parameters derived from satellite multiangle observations (e.g., MODIS). In this study, the main objective is to investigate the potential for linking a popular physical model (PROSAIL) with the widely used kernel-driven Ross-Li models. To do this, the BRFs and albedo are generated by the physical PROSAIL in a forward model, and then the simulated BRFs are input into the kernel-driven BRDF model for retrieval of the BRFs and albedo in the same viewing and solar geometries. To further strengthen such an investigation, a variety of field-measured multiangle reflectances have also been used to investigate the potential for linking these two models. For simulated BRFs generated by the PROSAIL model at 659 and 865 nm, the two models are generally comparable to each other, and the resultant root mean square errors (RMSEs) are 0.0092 and 0.0355, respectively, although some discrepancy in the simulated BRFs can be found at large average leaf angle (ALA) values. Unsurprisingly, albedos generated by the method are quite consistent, and 99.98% and 97.99% of the simulated white sky albedo (WSA) has a divergence less than 0.02. For the field measurements, the kernel-driven model presents somewhat better model-observation congruence than the PROSAIL model. The results show that these models have an overall good consistency for both field-measured and model-simulated BRFs. Therefore, there is potential for linking these two models for looking into the retrieval of canopy biophysical/structural variables through a simulation method, particularly from the current archive of the global routine MODIS BRDF parameters that were produced by the kernel-driven BRDF model; however, erectophile vegetation must be further examined.

Forest-canopy height is an important parameter for the estimation of forest biomass and terrestrial carbon flux and climate-change research at regional and global scales. Currently, various methods combining Light Detection and Ranging (LiDAR) data with various auxiliary data, particularly satellite remotely sensed reflectances, have been widely used to produce spatially continuous canopy-height products. However, current methods in use for remote sensing reflectances mainly focus on the nadir view direction, while anisotropic reflectances, which are theoretically more sensitive to the forest canopy height in the multiangle remote sensing field, have rarely been explored. Here, we attempted to examine the potential of using modeled multiangle reflectances at three typical viewing angles (i.e., from the hotspot, darkspot, and nadir directions) to estimate forest-canopy height as auxiliary data sources. First, the sensitivities of the typical angular reflectances as a function of forest canopy height were fully examined using the Extended Fourier Amplitude Sensitivity Test (EFAST) method based on the 4-scale Bidirectional Reflectance Distribution Function (BRDF) model simulations. This indicated that reflectances in the off-nadir viewing directions are generally sensitive to canopy-height variations. Then, the canopy heights were extracted from airborne Laser Vegetation Imaging Sensor (LVIS) data, which were further divided into training and validation data. Moderate Resolution Imaging Spectroradiometer (MODIS) multiangle reflectances at typical viewing angles were calculated from the MODIS BRDF parameter product (MCD43A1, version 6) as partial training-input data, based on a hotspot-adjusted, kernel-driven linear BRDF model. Subsequently, the Random Forest (RF) machine learning model was trained to acquire the relationship between the extracted canopy heights and the corresponding MODIS typical viewing reflectances. The trained model was further applied to estimate the canopy height metrics in the study areas of Howland Forest, Harvard Forest, and Bartlett Forest. Finally, the estimated canopy heights were independently validated by canopy heights extracted from the LVIS data. The results indicate that the canopy heights modeled through this method exhibit generally high accordance with the LVIS-derived canopy heights (R = 0.65−0.67; RMSE = 3.63−5.78). The results suggest that the MODIS multiangle reflectance data at typical observation angles contain important information regarding forest canopy height and can, therefore, be used to estimate forest canopy height for various ecological applications.

The foliage Clumping Index (CI) is an important vegetation structure parameter that allows for the accurate separation of sunlit and shaded leaves in a canopy. The CI and its seasonality are critical for global Leaf Area Index (LAI) estimating and ecological modelling. However, the cover of snow tends to reduce the reflectance anisotropy of the vegetation canopy and thus probably influences CI estimates. In this paper, we investigate the influence of snow on the magnitude and seasonal variation of the CI retrieved from Moderate-resolution Imaging Spectroradiometer (MODIS) Bidirectional Reflectance Distribution Function (BRDF) products based on field-measured CI and statistics from the global MODIS CI product. We find that the backup algorithm can effectively correct abnormally large CI values and obtain more reasonable CI retrievals than the main algorithm without any constraints in snow-covered areas. Validation indicates that the time-series CI product shows the potential in investigating the trajectories of the clumping effect in snow seasons. For evergreen forests, the clumping effect is relatively stable throughout the year; however, for deciduous vegetation types, CI values tend to display significant seasonal variations. This study suggests that the latest version of the global MODIS CI product, in which the backup algorithm is used to process the snow-covered pixels, has improved accuracy for CI retrievals in snow-covered areas and thus is probably more suitable as the input parameter for ecological and meteorological models.

The original kernel-driven bidirectional reflectance distribution function (BRDF) models were developed based on soil-vegetation systems. To further improve the ability of the models to characterize the snow surface scattering properties, a snow kernel was derived from the asymptotic radiative transfer (ART) model and used in the kernel-driven BRDF model framework. However, there is a need to further evaluate the influence of using this snow kernel to improve the original kernel-driven models in snow albedo retrieval applications. The aim of this study is to perform such an evaluation using a variety of snow BRDF data. The RossThick-Roujean (RTR) model is used as a framework for taking in the new snow kernel (hereafter named the RTS model) since the Roujean geometric-optical (GO) kernel captures a neglectable hotspot effect and represents a more prominent dome-shaped BRDF, especially at a small solar zenith angle (SZA). We obtained the following results: (1) The RTR model has difficulties in reconstructing the snow BRDF shape, especially at large SZAs, which tends to underestimate the reflectance in the forward direction and overestimate reflectance in the backward direction for various data sources. In comparison, the RTS model performs very well in fitting snow BRDF data and shows high accuracy for all data. (2) The RTR model retrieved snow albedos at SZAs = 30°–70° are underestimated by 0.71% and 0.69% in the red and near-infrared (NIR) bands, respectively, compared with the simulation results of the bicontinuous photon tracking (bic-PT) model, which serve as “real” values. However, the albedo retrieved by the RTS model is significantly improved and generally agrees well with the simulation results of the bic-PT model, although the improved model still somewhat underestimates the albedo by 0.01% in the red band and overestimates the albedo by 0.05% in the NIR band, respectively, at SZAs = 30°–70°, which may be negligible. (3) The albedo derived by these two models shows a high correlation (R² > 0.9) between the field-measured and Polarization and Directionality of the Earth’s Reflectances (POLDER) data, especially for the black-sky albedo. However, the albedo derived using the RTR model is significantly underestimated compared with the RTS model. The RTR model underestimates the black-sky albedo (white-sky albedo) retrievals by 0.62% (1.51%) and 0.93% (2.08%) in the red and NIR bands, respectively, for the field-measured data. The shortwave black-sky and white-sky albedos derived using the RTR model for the POLDER data are underestimated by 1.43% and 1.54%, respectively, compared with the RTS model. These results indicate that the snow kernel in the kernel-driven BRDF model frame is more accurate in snow albedo retrievals and has the potential for application in the field of the regional and global energy budget.

Clumping index (CI) is a canopy structural variable important for modeling the terrestrial biosphere, but its retrieval from remote sensing data remains one of the least reliable. The majority of regional or global CI products available so far were generated from multiangle optical reflectance data. However, these reflectance-based estimates have well-known limitations, such as the mere use of a linear relationship between the normalized difference hotspot and darkspot (NDHD) and CI, uncertainties in bidirectional reflectance distribution function (BRDF) models used to calculate the NDHD, and coarse spatial resolutions (e.g., hundreds of meters to several kilometers). To remedy these limitations and develop alternative methods for large-scale CI mapping, here we explored the use of spaceborne lidar—the Geoscience Laser Altimeter System (GLAS)—and proposed a semi-physical algorithm to estimate CI at the footprint level. Our algorithm was formulated to leverage the full vertical canopy profile information of the GLAS full-waveform data; it converted raw waveforms to forest canopy gap distributions and gap fractions of random canopies, which was used to estimate CI based on the radiative transfer theory and a revised Beer–Lambert model. We tested our algorithm over two areas in China—the Saihanba National Forest Park and Heilongjiang Province—and assessed its relative accuracies against field-measured CI and MODIS CI products. We found that reliable estimation of CI was possible only for GLAS waveforms with high signal-to-noise ratios (e.g., >65) and at gentle slopes (e.g., <12°). Our GLAS-based CI estimates for high-quality waveforms compared well to field-based CI (i.e., R² = 0.72, RMSE = 0.07, and bias = 0.02), but they showed less correlation to MODIS CI (e.g., R² = 0.26, RMSE = 0.12, and bias = 0.04). The difference highlights the impact of the scale effect in conducting comparisons of products with huge differences resolution. Overall, our analyses represent the first attempt to use spaceborne lidar to retrieve high-resolution forest CI and our algorithm holds promise for mapping CI globally.