Rui Xie

Multi-angular remote sensing observation contains crucial information on forest structure parameters. Here, our goal is to examine the ability of multi-angular indices, which are constructed by the typical-angular reflectances in red and NIR bands from MODIS observations, for the retrieval of forest biomass based on the field-measured above-ground biomass (AGB) data. Specifically, we employed the updated version of the MCD43A1 BRDF parameter product as an input for BRDF models to reconstruct the MODIS typical-angular reflectances. Furthermore, we evaluated the effects of different configurations of BRDF models and solar zenith angles (SZA) on forest AGB estimation using our developed multi-angular indices. The semivariogram analysis strategy combined with Landsat ground-surface reflectance data was employed to determine the MODIS pixel heterogeneity; the survey data from field sites of homogeneous pixels was used in our analysis and validation. The results show that our developed multi-angular indices based on a hot-revised BRDF model, under a SZA of 45°, when combined with forest cover information, can account for up to 72% of the variation forest AGB, with an RMSE = 45 Mg/ha. We also found that different kernels for the BRDF models influenced the weight parameters of the biomass inversion equation but did not significantly affect the estimated AGB. In conclusion, our method can enable the better usage of MODIS multi-angular observations for forest AGB estimation.

The clumping index (CI) quantifies the foliage grouping within distinct canopies relative to randomly distributed canopies, which plays an important role in the vegetation radiative regime. Among the vegetation structure parameters, the global CI map can be retrieved by using multiangle remote sensing observations. The bidirectional reflectance distribution function (BRDF)/Albedo product (MCD43) of the Moderate-Resolution Imaging Spectroradiometer (MODIS) is the crucial input data of the global CI product, which provides validated spatiotemporal continuous directional reflectance data. To determine the impacts of updating the MCD43 product on the consistency and performance of global CI products, CIs retrieved from different MCD43 versions (Collection V005/V006, C5/6) were compared on a global scale and validated with field-measured CI data. The results showed that the global and seasonal comparisons of C5 and C6 CI data are generally consistent. Compared to C5 CI data, C6 CI data have improved quality with more main algorithm retrievals and fewer case of missing data. The comparisons over the field measurements indicate that both versions of CI data agree with field-measured CI data in terms of values and seasonal variations, while C6 CI data (R² = 0.89, RMSE = 0.05, bias = 0.02) are closer to field CIs than C5 CI data (R² = 0.80, RMSE = 0.07, bias = 0.03), indicating a higher accuracy for C6 CI data. The monthly CI is recommended for characterizing the overall seasonal patterns of surface CIs.

Canopy structure parameters (e.g., leaf area index (LAI)) are key variables of most climate and ecology models. Currently, satellite-observed reflectances at a few viewing angles are often directly used for vegetation structure parameter retrieval; therefore, the information content of multi-angular observations that are sensitive to canopy structure in theory cannot be sufficiently considered. In this study, we proposed a novel method to retrieve LAI based on modelled multi-angular reflectances at sufficient sun-viewing geometries, by linking the PROSAIL model with a kernel-driven Ross-Li bi-directional reflectance function (BRDF) model using the MODIS BRDF parameter product. First, BRDF sensitivity to the PROSAIL input parameters was investigated to reduce the insensitive parameters. Then, MODIS BRDF parameters were used to model sufficient multi-angular reflectances. By comparing these reference MODIS reflectances with simulated PROSAIL reflectances within the range of the sensitive input parameters in the same geometries, the optimal vegetation parameters were determined by searching the minimum discrepancies between them. In addition, a significantly linear relationship between the average leaf angle (ALA) and the coefficient of the volumetric scattering kernel of the Ross-Li model in the near-infrared band was built, which can narrow the search scope of the ALA and accelerate the retrieval. In the validation, the proposed method attains a higher consistency (root mean square error (RMSE) = 1.13, bias = −0.19, and relative RMSE (RRMSE) = 36.8%) with field-measured LAIs and 30-m LAI maps for crops than that obtained with the MODIS LAI product. The results indicate the vegetation inversion potential of sufficient multi-angular data and the ALA relationship, and this method presents promise for large-scale LAI estimation.

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., 2 = 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.

The vertical foliage profile (VFP) and leaf area index (LAI) are critical descriptors in terrestrial ecosystem modeling. Although light detection and ranging (lidar) observations have been proven to have potential for deriving the VFP and LAI, existing methods depend only on the received waveform information and are sensitive to additional input parameters, such as the ratio of canopy to ground reflectance. In this study, we proposed a new method for retrieving forest VFP and LAI from Ice, Cloud and land Elevation Satellite (ICESat) Geoscience Laser Altimeter System (GLAS) data over two sites similar in their biophysical parameters. Our method utilized the information from not only the interaction between the laser and the forest but also the sensor configuration, which brought the benefit that our method was free from an empirical input parameter. Specifically, we first derived the transmitted energy profile (TEP) through the lidar 1-D radiative transfer model. Then, the obtained TEP was utilized to calculate the vertical gap distribution. Finally, the vertical gap distribution was taken as input to derive the VFP based on the Beer–Lambert law, and the LAI was calculated by integrating the VFP. Extensive validations of our method were carried out based on the discrete anisotropic radiative transfer (DART) simulation data, ground-based measurements, and the Moderate Resolution Imaging Spectroradiometer (MODIS) LAI product. The validation based on the DART simulation data showed that our method could effectively characterize the VFP and LAI under various canopy architecture scenarios, including homogeneous turbid and discrete individual-tree scenes. The ground-based validation also proved the feasibility of our method: the VFP retrieved from the GLAS data showed a similar trend with the foliage distribution density in the GLAS footprints; the GLAS LAI had a high correlation with the field measurements, with a determination coefficient (R²) of 0.79, root mean square error (RMSE) of 0.49, and bias of 0.17. Once the outliers caused by low data quality and large slope were identified and removed, the accuracy was further improved, with R² = 0.85, RMSE = 0.35, and bias = 0.10. However, the MODIS LAI product did not present a good relationship with the GLAS LAI. Relative to the GLAS LAI, the MODIS LAI showed an overestimation in the low and middle ranges of the LAI and a saturation at high values of approximately LAI = 5.5. Overall, this method has the potential to produce continental- and global-scale VFP and LAI datasets from the spaceborne lidar system.

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 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.