RossThick-LiSparseReciprocal Research Articles

An improved BRDF hotspot model and its use in VLIDORT for studying the impact of atmospheric scattering on hotspot directional signatures in the atmosphere

Abstract. The term “hotspot” refers to the sharp increase in the reflectance occurring when incident (solar) and reflected (viewing) directions almost coincide in the backscatter direction. The accurate simulation of hotspot directional signatures is important for many remote sensing applications. The RossThick–LiSparse–Reciprocal (RTLSR) bidirectional reflectance distribution function (BRDF) model is widely used in radiative transfer simulations, and the hotspot model mostly used is from Maignan–Bréon, but it typically requires large values of numerical quadrature and Fourier expansion terms in order to represent the hotspot accurately for its use coupled with atmospheric radiative transfer modeling (RTM). In this paper, we have developed a modified version based on the Maignan–Bréon's hotspot BRDF model that converges much faster numerically, making it more practical for use in the RTMs that require Fourier expansion of BRDF to simulate the top-of-atmosphere (TOA) hotspot signatures, such as in the RTM models using the Doubling–Adding or discrete ordinate method. Using the vector linearized discrete ordinate radiative transfer model (VLIDORT), we found that reasonable TOA–hotspot accuracy can be obtained with just 23 Fourier terms for clear atmosphere and 63 Fourier terms for atmosphere with aerosol scattering. In order to study the impact of molecular and aerosol scattering on hotspot signatures, we carried out a number of hotspot signature simulations with VLIDORT. We confirmed that (1) atmospheric molecule scattering and the existence of aerosol tend to smooth out the hotspot signature at the TOA and that (2) the hotspot signature at the TOA in the near-infrared is larger than in the visible, and its impact by surface reflectance is more significant. As the hotspot amplitude at the TOA with aerosol scattering included is smaller than that with molecular scattering only, the amplitude of hotspot signature at the surface is likely underestimated in the previous analysis based on the POLDER measurements, where the atmospheric correction was based on a single-scatter Rayleigh-only calculation. This modified model can calculate the amplitude of the hotspot accurately, and, as it agrees very well with the original RossThick model away from the hotspot region, this model can be simply used in conditions with and without hotspots. However, there are some differences in this modified model compared to the original Maignan–Bréon model for the scattering angles close to the hotspot point; thus, it may not be appropriate for those who need an exact representation of the hotspot angular signature.

Analyzing GOES-R ABI BRDF-adjusted EVI2 time series by comparing with VIIRS observations over the CONUS

The Advanced Baseline Imager (ABI) instruments onboard the Geostationary Operational Environmental Satellite-R series (GOES-R) provide new opportunities to expand their applications in global terrestrial ecology and environmental sciences. It offers a 5–15 min temporal resolution and 0.5–1 km spatial resolution that greatly benefit those environmental applications where spatial resolution can be compromised in favor of more frequent imagery. However, ABI captures reflected radiation at varying sun positions throughout the day, causing significant diurnal variations. The varying off-nadir view across ABI full disk presents further challenges for generating comparable and consistent vegetation observations, such as the two-band Enhanced Vegetation Index (EVI2). To address these challenges, three Bidirectional Reflectance Distribution Function (BRDF) models (RossThick-LiSparse-Reciprocal (RTLSR), Walthall, and Tian models) are examined to generate high-quality BRDF-adjusted EVI2 based on GeoNEX GOES-16 ABI surface reflectance. The ABI BRDF retrievals are then systematically compared to the corresponding VIIRS (Visible Infrared Imaging Radiometer Suite) Nadir BRDF-Adjusted Reflectance (NBAR) products (VNP43 V001) retrieved from RTLSR model. The viewing geometry effects are further analyzed over various land cover types using time series analysis, Pearson Correlation Coefficients (PCC), and Relative Difference (RD) between ABI and VIIRS high-quality time series inversions. The results indicate that the RTLSR model fits well the diurnal ABI observed reflectance and EVI2 values (even for the hot-spot phenomena), and outperforms the Walthall and Tian models for generating ABI high-quality BRDF retrievals with consistently low Root Mean Square Error (RMSE<0.02) and Weight of Determination (WoD < 0.5) across the CONUS. The annual proportion of ABI high-quality BRDF retrievals with a 3-day moving window using the RTLSR model is spatially and statistically comparable to that of VIIRS NBAR with a 16-day retrieval period. Compared to the time series of high-quality VIIRS NBAR EVI2, the BRDF adjusted ABI EVI2 time series still show considerable impacts of view zenith angles over evergreen forest, shrubland, and savanna pixels. However, such impacts are not significant over spatially homogeneous pixels such as cropland and deciduous forests. These results emphasize the need for BRDF correction when considering the impact of off-nadir observations from geostationary satellites on vegetation monitoring. They also demonstrate the advantages of geostationary satellites in capturing rapid changes in greenness trajectories.

A Practical Approach to Improve the MODIS MCD43A Products in Snow-Covered Areas

The MODerate Resolution Imaging Spectroradiometer (MODIS) MCD43A products have been extensively applied in the remote sensing field, but recent researchers have demonstrated that these products still had the potential to be further improved by using the latest development of the kernel-driven model [RossThick-LiSparseReciprocal-Snow (RTLSRS)] in snow-covered areas, since the MCD43A product algorithm [RossThick-LiSparseReciprocal (RTLSR)] needed to be improved for the accurate simulation of snow bidirectional reflectance distribution function (BRDF) signatures. In this paper, we proposed a practical approach to improve the MCD43A products, which used the Polarization and Directionality of the Earth's Reflectance (POLDER) observations and random forest algorithm to establish the relationship between the BRDF parameters (MCD43A1) estimated by the RTLSR and RTLSRS models. We applied this relationship to correct the MCD43A1 product and retrieve the corresponding albedo (MCD43A3) and nadir reflectance (MCD43A4). The results obtained highlight several aspects: (a) The proposed approach can perform well in correcting BRDF parameters [root mean square error (RMSE) = ~0.04]. (b) The corrected BRDF parameters were then used to retrieve snow albedo, which matched up quite well with the results of the RTLSRS model. (c) Finally, the snow albedo retrieved by the proposed approach was assessed using ground-based albedo observations. Results indicated that the retrieved snow albedo showed a higher accuracy as compared to the station measurements (RMSE = 0.055, bias = 0.005), which was better than the results of the MODIS albedo product (RMSE = 0.064, bias = −0.018), especially at large angles. These results demonstrated that this proposed approach presented the potential to further improve the MCD43A products in snow-covered areas.

Using HJ-1 CCD and MODIS Fusion Data to Invert HJ-1 NBAR for Time Series Analysis, a Case Study in the Mountain Valley of North China

HJ-1 charge-coupled device (CCD) data with high temporal and medium spatial resolution are widely used in environmental and disaster monitoring in China. However, due to bad weather, it is difficult to obtain sufficient time-continuous HJ-1 CCD data for environmental monitoring. In this study, the mountain valley with farmland and forestland in North China is selected as the experimental area, and HJ-1 CCD and moderate resolution imaging spectroradiometer (MODIS) data are used in the case study. An improved method of fusing data and inversing surface reflectivity is presented to obtain the HJ-1 inversion network-based application resolution (NBAR) data using linear matching of the Ross Thick-Li Sparse Reciprocal (RTLSR) model, and then predicted reflectivity using the seasonal autoregressive integrated moving average (SARIMA) model. The fusion data have advantages of high spatial and temporal resolution, as well as meeting the requirements of high quality and quantity of small-scale regional data. This case study provides a feasibility method for the HJ-1 satellites to produce the secondary products for small-scale remote sensing ground surface research. It also provides a reference for dynamic information acquisition and application of small satellite data, contributing to the improvement in RS estimation of surface environment variables.

Evaluation of Orbital Drift Effect on Proba-V Surface Reflectances Time Series

Multi-temporal consistency of space-borne observations is an essential requirement for studying inter-annual changes and trends of satellite-derived biophysical products. The Proba-V mission, launched in 2013, was designed to ensure the continuity of the SPOT-VEGETATION long-term data record of global daily observations for land applications. The suitability of Proba-V to provide a temporally consistent data record is, however, potentially jeopardized by the orbital drift effect, which is known to induce spurious trends in time series. The aim of this paper is therefore to evaluate, for the first time, the orbital drift effect on Proba-V surface reflectance time series at 1 km resolution. In order to reliably identify such an effect, a two-fold approach is adopted. A simulation study is first defined to predict the temporal anomalies induced by the drifting illumination conditions. The numerical simulations are used as a benchmark to predict the impact of the drift for a range of sun-viewing angles. Real observations are then analyzed over a large set of land sites, globally spread and spanning a wide range of surface and environmental conditions. The surface anisotropy is characterized using the Ross-Thick Li-Sparse Reciprocal (RTLSR) Bidirectional Reflectance Distribution Function (BRDF) model. Both the simulation and the analysis of real observations consistently show that the orbital drift induces distinct and opposite trends in the two sides of the sensor across-track swath. Particularly, a positive drift is estimated in backward and a negative one in the forward scattering direction. When observations from all angular conditions are retained, these opposite trends largely compensate, with no remaining statistically significant drifts in time series of surface reflectances or Normalized Difference Vegetation Index (NDVI). As such, the Proba-V archive at 1 km resolution can be reliably used for inter-annual vegetation studies.

Extending a Linear Kernel-Driven BRDF Model to Realistically Simulate Reflectance Anisotropy Over Rugged Terrain

Bidirectional reflectance distribution function (BRDF) models are used to correct surface bidirectional effects and estimate land surface albedo. Many operational BRDF/albedo algorithms adopt a Roujean linear kernel-driven BRDF (RLKB) model because of its simple form and good performance in fitting multidirectional surface reflectance values. However, this model does not explicitly consider topographic effects, resulting in errors when applied over rugged terrain. To address this issue, we proposed a hybrid algorithm suitable for both flat and rugged terrain, called topographical kernel-driven (Topo-KD). First, we constructed a linear kernel-driven BRDF model considering terrain (LKB_T) which describes the topographic effects with a mountain radiative transfer (MRT) model. Then, the Topo-KD algorithm adaptively selects the most suitable model (RLKB or LKB_T) according to the terrain conditions and fitting residuals. The performances of Topo-KD and RLKB using the RossThick–LiSparseReciprocal (RTLSR) kernel are compared using simulated data sets and moderate-resolution imaging spectroradiometer (MODIS) observations. The results show that the BRDF of the pixel is affected by topography. But the RTLSR model does not specifically account for it, resulting in larger biases over rugged terrain than the Topo-KD algorithm in both the red and near-infrared (NIR) bands. The experiment using MODIS data sets demonstrates that the Topo-KD algorithm reduces fitting residuals in the red and NIR bands by 21.5% and 27.4% compared with the RTLSR model. These results indicate that the Topo-KD algorithm can be a better choice for retrieving land surface parameters and describing the radiative transfer process in mountainous areas.

Assessment of Improved Ross–Li BRDF Models Emphasizing Albedo Estimates at Large Solar Angles Using POLDER Data

Surface albedo is closely related to the Earth’s energy budget and is usually estimated by integrating remotely sensed bidirectional reflectance distribution function (BRDF) data based on the widely used Ross–Li kernel-driven models. However, for large solar zenith angles (i.e., SZAs > 70°), albedo estimation using the operational algorithm of the Moderate Resolution Imaging Spectroradiometer (MODIS), i.e., RossThick-LiSparseReciprocal (RTLSR), is not recommended because it is reported to somewhat underestimate the black-sky albedo (BSA) at large SZAs based on ground albedo measurements. Recently, various combinations of the Ross–Li BRDF models with improved capabilities have been developed, and the assessments of these models based on worldwide satellite BRDF data with good spatial sampling, particularly at the large view and solar angles, will be important to improve an understanding of their performance in estimating intrinsic albedos. Following previous studies, the objective of this study is to further assess a series of hotspot-corrected Ross–Li models by demonstrating their ability to fit the POLarization and Directionality of the Earth’s Reflectances (POLDER) data sets and estimate albedo, especially at large SZAs, based on selected concurrent POLDER and MODIS data. The hotspot-corrected RTLSR model obtained by combining the RossThickChen and LiSparseReciprocalChen kernels (RTLSR_C) shows the best fitting ability, with a high cumulative frequency of small root-mean-square errors (RMSEs), thus confirming previous conclusions. Model differences mainly appear in albedo estimates, especially BSA estimates at large SZAs. The BSAs estimated by other models are significantly different from the RTLSR_C estimates in the near-infrared (NIR) and red bands as the SZA increases to approximately 60° and 70°, respectively. In this case, RossThinChen-LiSparseReciprocalChen (RTNLSR_C) yields higher BSA estimates than those of RTLSR_C. Comparisons of the MODIS and POLDER albedos estimated with Ross–Li models show that models with the RossThinChen kernel yield higher BSA estimates than those of the RTLSR_C model as the SZA increases. The results indicate that the retrieved albedo is likely to be more accurate with appropriately selected kernels for BRDF models at large SZAs, providing guidance for selecting suitable combinations of multiple kernels.

Reconstruction of Spatiotemporally Continuous MODIS-Band Reflectance in East and South Asia from 2012 to 2015

To reconstruct Moderate Resolution Imaging Spectroradiometer (MODIS) band reflectance with optimal spatiotemporal continuity, three bidirectional reflectance distribution function (BRDF) models—the Ross-Thick-Li-Sparse Reciprocal (RTLSR) model, Gao model, and adjusted BF model—were used to retrieve MODIS-band reflectance for cloudy MODIS pixels according to different inversion conditions with a proposed filling algorithm. Then, a spatiotemporally continuous MODIS-band reflectance dataset for most of Asia with more than 98% spatiotemporal coverage was reconstructed from 2012 to 2015. The validation highlighted an evident improvement in filling cloudy MODIS observations; a reasonable spatial distribution, such as in South Asia and Southeast Asia; and acceptable precision for the filled MODIS pixels, with the root mean square error percentage (RMSE%) at 9.7–9.8% and 12–16% for the Gao and adjusted BF models, respectively. In the course of reconstructing the spatiotemporal continuous MODIS-band reflectance, the differences among the three models were discussed further. For a 16-day period with a stable and unchanged land surface, the RTLSR model, as a basic model, accurately derived land surface reflectance (no more than 10% RMSE% for MCD43C1 V006 band 1) and outperformed the other two models. When the inversion period is sufficiently long (e.g., 108 days, 188 days, 268 days, or a full year), the Gao/adjusted BF model provides better precision than the RTLSR model by considering the normalized difference vegetation index (NDVI) and soil moisture/NDVI as intermediate variables used to adjust the BRDF parameters in real time. The Gao model is optimal when the inversion period is sufficiently long. Based on combining the RTLSR model and Gao/adjusted BF model, we proposed a filling algorithm to derive a dataset of MODIS-band reflectance with optimal spatiotemporal continuity.

Sun-Angle Effects on Remote-Sensing Phenology Observed and Modelled Using Himawari-8

Satellite remote sensing of vegetation at regional to global scales is undertaken at considerable variations in solar zenith angle (SZA) across space and time, yet the extent to which these SZA variations matter for the retrieval of phenology remains largely unknown. Here we examined the effect of seasonal and spatial variations in SZA on retrieving vegetation phenology from time series of the Normalized Difference Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI) across a study area in southeastern Australia encompassing forest, woodland, and grassland sites. The vegetation indices (VI) data span two years and are from the Advanced Himawari Imager (AHI), which is onboard the Japanese Himawari-8 geostationary satellite. The semi-empirical RossThick-LiSparse-Reciprocal (RTLSR) bidirectional reflectance distribution function (BRDF) model was inverted for each spectral band on a daily basis using 10-minute reflectances acquired by H-8 AHI at different sun-view geometries for each site. The inverted RTLSR model was then used to forward calculate surface reflectance at three constant SZAs (20°, 40°, 60°) and one seasonally varying SZA (local solar noon), all normalised to nadir view. Time series of NDVI and EVI adjusted to different SZAs at nadir view were then computed, from which phenological metrics such as start and end of growing season were retrieved. Results showed that NDVI sensitivity to SZA was on average nearly five times greater than EVI sensitivity. VI sensitivity to SZA also varied among sites (biome types) and phenological stages, with NDVI sensitivity being higher during the minimum greenness period than during the peak greenness period. Seasonal SZA variations altered the temporal profiles of both NDVI and EVI, with more pronounced differences in magnitude among NDVI time series normalised to different SZAs. When using VI time series that allowed SZA to vary at local solar noon, the uncertainties in estimating start, peak, end, and length of growing season introduced by local solar noon varying SZA VI time series, were 7.5, 3.7, 6.5, and 11.3 days for NDVI, and 10.4, 11.9, 6.5, and 8.4 days for EVI respectively, compared to VI time series normalised to a constant SZA. Furthermore, the stronger SZA dependency of NDVI compared with EVI, resulted in up to two times higher uncertainty in estimating annual integrated VI, a commonly used remote-sensing proxy for vegetation productivity. Since commonly used satellite products are not generally normalised to a constant sun-angle across space and time, future studies to assess the sun-angle effects on satellite applications in agriculture, ecology, environment, and carbon science are urgently needed. Measurements taken by new-generation geostationary (GEO) satellites offer an important opportunity to refine this assessment at finer temporal scales. In addition, studies are needed to evaluate the suitability of different BRDF models for normalising sun-angle across a broad spectrum of vegetation structure, phenological stages and geographic locations. Only through continuous investigations on how sun-angle variations affect spatiotemporal vegetation dynamics and what is the best strategy to deal with it, can we achieve a more quantitative remote sensing of true signals of vegetation change across the entire globe and through time.

Improving Kernel-Driven BRDF Model for Capturing Vegetation Canopy Reflectance With Large Leaf Inclinations

Semiempirical, kernel-driven linear bidirectional reflectance distribution function (BRDF) models are widely used to characterize vegetation reflectance anisotropy and provide land surface bidirectional reflectance factor (BRF) products at the regional and global scales. However, these models usually imply an assumption of spherical leaf inclination. The effects of such an ideal assumption on simulating surface BRF remain inadequately quantified. In this article, we first evaluated the effects of leaf inclination on the most commonly used kernel-driven RossThick-LiSparse-Reciprocal (RTLSR) model by using the reflectance benchmark simulated by the mature PROSAIL (PROSAIL+SAIL) radiative transfer model. Subsequently, we improved the RTLSR model into a four-parameter version (RTLSRV4p) with a new volumetric scattering kernel derived from the assumption of vertical leaf inclination. Finally, the proposed RTLSRV4p model was validated by PROSAIL canopy BRF simulations, in situ canopy BRF measurements, and wide-angle infrared dual-mode line/area Aarray scanner (WIDAS) airborne observations. Validation results demonstrate that RTLSRV4p improves vegetation reflectance characterization for large leaf inclinations compared to the original RTLSR model, especially for the near-infrared (NIR) spectral domain. When validated against the simulated canopy BRFs, the mean root-mean-square error (RMSE), mean absolute percentage error (MAPE), bias, and coefficient of determination ( R 2) were improved from 0.0810, 31.63%, 0.0651, 0.6578 to 0.0453, 13.38%, 0.0326, 0.8734. Using the in situ BRF measurement, the fitted RMSE, MAPE, bias, and R 2 were improved from 0.0917, 14.31%, 0.0728, and 0.5776 to 0.0226, 3.35%, 0.0166, and 0.9744. These validation metrics were improved from 0.0423, 10.85%, 0.0347, and 0.6598 to 0.0258, 5.85%, 0.0181, and 0.8732 when compared against the WIDAS observations. The RTLSRV4p model also shows good performance in calculating the reflectance of vegetation. These extensive validations suggest that RTLSRV4p is promising for capturing vegetation reflectance under large leaf inclinations.

An assessment of the performance of two snow kernels in characterizing snow scattering properties

ABSTRACTThe kernel-driven RossThick-LiSparseReciprocal (RTLSR) bidirectional reflectance distribution function (BRDF) model has been widely used in the quantitative remote sensing community. However, the performance of this model is challenged when modelling the optical scattering properties of pure snow surfaces. Recently, two snow kernels have been developed to improve the snow anisotropic reflectance in the kernel-driven RTLSR model framework. However, the performances of these two snow kernels must be assessed to identify their potential applications. Therefore, we assess the performances of these two kernels using various BRDF data sources. Our findings demonstrate their differences in several aspects. (1) These two kernels differ in characterizing the variability in BRDF shape as a function of the solar zenith angle (SZA). As the SZA increases, the shape of snow kernel derived by the asymptotic radiative transfer (ART) model (hereinafter named the ART method) changes from a dome shape to bowl shape, which agrees well with the simulation data of the bicontinuous photon tracking (bic-PT) model. The shape of the snow kernel proposed by Qu et al. (hereinafter named the Qu method) based on the Rahman-Pinty-Verstraete (RPV) model maintains a bowl shape for all SZAs. These differences in the kernel performances affect their abilities to fit snow BRDF data with different SZAs. (2) The corrected RTLSR models, with their respective snow kernels, are generally able to model the forward-scattering properties of snow surfaces compared with the original RTLSR model. However, the ART method performs better in capturing the BRDF variations in snow surfaces than the Qu method. This assessment provides an improved understanding of the performance of these two snow kernels and, thus, suggests further applications for the ART snow kernel in the kernel-driven BRDF model framework in the near future.

Derivation of Kernel Functions for Kernel-Driven Reflectance Model Over Sloping Terrain

The importance of the bidirectional reflectance distribution function (BRDF) has been well documented in quantitative remote sensing. The semiempirical, kernel-driven BRDF model is widely used to generate operational BRDF/albedo products due to its simplicity and accuracy. However, the effect of topography is rarely coupled with a kernel-based BRDF model. In this paper, a new kernel-driven reflectance model for sloping terrain (KDST) was developed based on the framework of the RossThick-LiSparse-Reciprocal (RTLSR) model. The slope, aspect, geotropic nature of the tree crown governed by gravity, component spectra contrasts, and diffuse irradiance were considered in the KDST model. The performance of KDST was evaluated by 3-D discrete anisotropic radiative transfer (DART) simulations, in situ measurements, and HJ-1A/B constellation charge-coupled device satellite observations. Using DART simulations, KDST reduces the maximum biases of the RTLSR from 15.6% and 29.7% to 7.1% and 4.8% for the surface bidirectional reflectance factor and hemispherical-direction reflectance factor, respectively. Compared with in situ measurements, KDST improves the reflectance simulation accuracy for the red and near-infrared bands from 0.0172 (18.65%) and 0.022 (4.77%) to 0.0060 (6.51%) and 0.0043 (1.05%), respectively. The preliminary comparison results indicate that KDST is promising for reflectance simulation over rugged terrain.

Development of a snow kernel to better model the anisotropic reflectance of pure snow in a kernel-driven BRDF model framework

The linear kernel-driven RossThick-LiSparseReciprocal (RTLSR) bidirectional reflectance distribution function (BRDF) model was originally developed from the simplified scenarios of continuous and discrete vegetation canopies, and has been widely used to fit multiangle observations of vegetation-soil systems of the land surface in many fields. Although this model was not developed explicitly for snow surfaces, it can capture the geometric-optical effect caused by the shadowing of rugged or undulating snow surfaces. However, in this study, this model has been further developed to better characterize the scattering properties of snow surface, which can also exhibit strongly forward-scattering behavior. This study proposes a new snow kernel to characterize the reflectance anisotropy of pure snow based on the asymptotic radiative transfer (ART) model that assumes snow can be modeled as a semi-infinite, plane-parallel, weakly absorbing light scattering layer. This new snow kernel adopts a correction term with a free parameter α to correct the analytic form of the ART model that has been reported to underestimate observed snow reflectance in the forward-scattering direction in the principal plane (PP), particularly in cases of a large viewing zenith angle (>60°). This snow kernel has now been implemented in the kernel-driven RTLSR BRDF model framework in conjunction with two additional kernels (i.e., the volumetric scattering kernel and geometric-optical scattering kernel) and is validated using observed and simulated multiangle data from three data sources. Pure snow targets were selected from the extensive archive of the Polarization and Directionality of the Earth's Reflectance (POLDER) BRDF data. Antarctic snow field measurements, which were taken from the top of a 32-m-tall tower at Dome C Station and include 6336 spectral bidirectional reflectance factors (BRFs), were also utilized. Finally, a set of simulated BRFs, generated by a hybrid scattering snow model that combines the geometric optics with vector radiative transfer theory, were used to further assess the proposed method. We first retrieve the value of the free parameter α for a comprehensive analysis using single multiangle snow data with a sufficient BRDF sampling. Then, we determine the optimally fixed value of the α parameter as prior information for potential users. The new snow kernel method is shown to be quite accurate, presenting a high correlation coefficient (R2 = ~0.9) and a negligible bias between the modeled BRFs and the various snow BRDF validation data. The finding demonstrates that this snow kernel provides an improved potential compared to that of the original kernel-driven model framework for a pure snow surface in many applications, particularly those involving the global water cycle and radiation budget, where snow cover plays an important role.

Characterization of Vegetation Physiognomic Types Using Bidirectional Reflectance Data

This paper presents an assessment of the bidirectional reflectance features for the classification and characterization of vegetation physiognomic types at a national scale. The bidirectional reflectance data at multiple illumination and viewing geometries were generated by simulating the Moderate Resolution Imaging Spectroradiometer (MODIS) Bidirectional Reflectance Distribution Function (BRDF) model parameters with Ross-Thick Li-Sparse-Reciprocal (RT-LSR) kernel weights. This research dealt with the classification and characterization of six vegetation physiognomic types—evergreen coniferous forest, evergreen broadleaf forest, deciduous coniferous forest, deciduous broadleaf forest, shrubs, and herbaceous—which are distributed all over the country. The supervised classification approach was used by employing four machine learning classifiers—k-Nearest Neighbors (KNN), Random Forests (RF), Support Vector Machines (SVM), and Multilayer Perceptron Neural Networks (NN)—with the support of ground truth data. The confusion matrix, overall accuracy, and kappa coefficient were calculated through a 10-fold cross-validation approach, and were also used as the metrics for quantitative evaluation. Among the classifiers tested, the accuracy metrics did not vary much with the classifiers; however, the Random Forests (RF; Overall accuracy = 0.76, Kappa coefficient = 0.72) and Support Vector Machines (SVM; Overall accuracy = 0.76, Kappa coefficient = 0.71) classifiers performed slightly better than other classifiers. The bidirectional reflectance spectra did not only vary with the vegetation physiognomic types, it also showed a pronounced difference between the backward and forward scattering directions. Thus, the bidirectional reflectance data provides additional features for improving the classification and characterization of vegetation physiognomic types at the broad scale.

An algorithm for the retrieval of the clumping index (CI) from the MODIS BRDF product using an adjusted version of the kernel-driven BRDF model

The clumping index (CI) characterizes the grouping of foliage relative to a random spatial distribution of leaves and is an important structural parameter for plant canopies that can influence canopy radiation regimes. Consequently, the CI is very useful for ecological and meteorological models. One method used to retrieve the CIs of plant canopies is to construct a linear relationship between the CI and the normalized difference between hotspot and dark spot (NDHD) angular index. This method requires a particularly accurate reconstruction of hotspot signatures, which are difficult to measure. In this study, we propose a framework to retrieve CIs from Moderate Resolution Imaging Spectroradiometer (MODIS) bidirectional reflectance distribution function (BRDF) parameters, which are generally based on linear CI-NDHD equations. The main algorithm is designed to retrieve CIs in the closed interval [0.33, 1.00]. This range is derived from the CI-NDHD equations and is thus called as the physical range here, although a modified lower boundary can be implemented in the future if necessary. If CIs are outside of this range, then a backup algorithm is designed to reprocess these so-called outlier CIs. The hotspot-adjusted version of the RossThick-LiSparseReciprocal (RTLSR) model (i.e., the RTCLSR model) is employed to reconstruct the hotspot signatures for the MODIS BRDF parameters. This method simplifies the hotspot reconstruction by using two hotspot parameters that are not distinctly scale-dependent particularly in the context of an inhomogeneous coarse spatial resolution. To evaluate this algorithm framework, we collect dozens of global field-measured CIs and calculate their determination coefficient (R2), root mean square error (RMSE) and bias relative to MODIS CIs derived using both the main algorithm and the backup algorithm. Our results show that this framework can derive MODIS CIs with a high accuracy (i.e., R2 = 0.80 (0.72), RMSE = 0.07 (0.12), bias = −0.03 (−0.10)) using the main (backup) algorithms and that it shows promise for various ecological applications, especially in combination with the leaf area index (LAI).

A method for improving hotspot directional signatures in BRDF models used for MODIS

The semi-empirical, kernel-driven, linear RossThick-LiSparseReciprocal (RTLSR) Bidirectional Reflectance Distribution Function (BRDF) model is used to generate the routine MODIS BRDF/Albedo product due to its global applicability and the underlying physics. A challenge of this model in regard to surface reflectance anisotropy effects comes from its underestimation of the directional reflectance signatures near the Sun illumination direction; also known as the hotspot effect. In this study, a method has been developed for improving the ability of the RTLSR model to simulate the magnitude and width of the hotspot effect. The method corrects the volumetric scattering component of the RTLSR model using an exponential approximation of a physical hotspot kernel, which recreates the hotspot magnitude and width using two free parameters (C1 and C2, respectively). The approach allows one to reconstruct, with reasonable accuracy, the hotspot effect by adjusting or using the prior values of these two hotspot variables. Our results demonstrate that: (1) significant improvements in capturing hotspot effect can be made to this method by using the inverted hotspot parameters; (2) the reciprocal nature allow this method to be more adaptive for simulating the hotspot height and width with high accuracy, especially in cases where hotspot signatures are available; and (3) while the new approach is consistent with the heritage RTLSR model inversion used to estimate intrinsic narrowband and broadband albedos, it presents some differences for vegetation clumping index (CI) retrievals. With the hotspot-related model parameters determined a priori, this method offers improved performance for various ecological remote sensing applications; including the estimation of canopy structure parameters.

Influence of solar zenith angle on the enhanced vegetation index of a Guyanese rainforest

In this study, the effect of solar zenith angle () on enhanced vegetation index (EVI) of a Guyanese tropical rainforest was studied. For this sub-crown resolution, hyperspectral data have been collected with an unmanned aerial vehicle (UAV) at five different solar zenith angles in a 1-day period. The hyperspectral data were used to simulate Moderate Resolution Imaging Spectroradiometer (MODIS) spectral bands and generate EVI. The linear trend of EVI with solar zenith angle at nadir viewing conditions was found to be –0.00285 (). The direction of this trend was in agreement with earlier studies, but with a differing magnitude. Analysis of EVI images with sub-crown resolution pointed to strong influence of canopy shadows on EVI, which is supported by other studies. Additionally, the EVI–solar zenith angle trend was investigated in the semi-empirical RossThick-LiSparse-Reciprocal (RTLSR) model implemented in the MODIS MCD43 product suite. A database of model parameters has been created and the EVI–solar zenith angle trend was modelled with each set of parameters. The linear approximated trend was found to be –0.00219 on average, only slightly weaker compared to the trend derived from the UAV. Further analysis of the relationship between the single RTLSR model parameters and the EVI–solar zenith angle trend showed that the RTLSR produces the trend for the right reason, namely canopy shadowing expressed by the near-infrared geometric kernel. In total, this study delivers further evidence that EVI is dependent on solar zenith angle and this effect is mediated through EVI’s sensitivity to within-canopy shadows.

Evaluation of the MODIS Albedo product over a heterogeneous agricultural area

In this article, the Moderate Resolution Imaging Spectroradiometer (MODIS) Bidirectional Reflectance Distribution Function (BRDF)/Albedo product (MCD43) is evaluated over a heterogeneous agricultural area in the framework of the Earth Observation: Optical Data Calibration and Information Extraction (EODIX) project campaign, which was developed in Barrax (Spain) in June 2011. In this method, two models, the RossThick-LiSparse-Reciprocal (RTLSR) (which corresponds to the MODIS BRDF algorithm) and the RossThick-Maignan-LiSparse-Reciprocal (RTLSR-HS), were tested over airborne data by processing high-resolution images acquired with the Airborne Hyperspectral Scanner (AHS) sensor. During the campaign, airborne images were retrieved with different view zenith angles along the principal and orthogonal planes. Comparing the results of applying the models to the airborne data with ground measurements, we obtained a root mean square error (RMSE) of 0.018 with both RTLSR and RTLSR-HS models. The evaluation of the MODIS BRDF/Albedo product (MCD43) was performed by comparing satellite images with AHS estimations. The results reported an RMSE of 0.04 with both models. Additionally, taking advantage of a homogeneous barley pixel, we compared in situ albedo data to satellite albedo data. In this case, the MODIS albedo estimation was (0.210 ± 0.003), while the in situ measurement was (0.204 ± 0.003). This result shows good agreement in regard to a homogeneous pixel.

Analysis of BRDF and Albedo Retrieved by Kernel-Driven Models Using Field Measurements

The algorithm for Model Bidirectional Reflectance Anisotropies of the Land Surface (AMBRALS) includes a series of kernel-driven bidirectional reflectance distribution function (BRDF) models. Among these models, the RossThick-LiSparse-Reciprocal (RTLSR) model has been selected as the operational MODIS BRDF/Albedo algorithm. With the newly developed LiTransit kernel and Ross kernel with the hotspot effect by Breon, there is a need to further evaluate models' potentials in retrieving land surface properties for user community. In this paper, 16 kernel-driven models in reciprocal form have been analyzed using 70 measurements mainly from ground campaigns of the First International Satellite Land Surface Climatology Project (ISLSCP) Field Experiment (FIFE) and the Boreal Ecosystem Atmosphere Study (BOREAS). Results show that all models under investigation generally fit the ground measurements with reasonable accuracy—although the “sparse” geometric-optical (GO) kernels appear to present better behaviors compared with the “dense” GO kernels—in combining with volumetric shapes. Estimated albedos from these models show high correlations, while the black-sky albedos (BSAs) present higher correlations in small solar zenith angles than in large ones, indicating that the difference in these kernels may arise from large solar geometry. A further investigation shows that the hotspot effects for these models display a significant discrepancy. Although the additional hotspot factor for Ross kernels significantly improves the hotspot effects, LiTransit kernel can well characterize the hotspot effects in combination with Ross kernels without the hotspot factor. Such an assessment may be helpful in understanding the models' potentials in retrieving vegetation structural characteristics.

Improving MODIS land cover classification by combining MODIS spectral and angular signatures in a Canadian boreal forest

This study explores the use of reflectance anisotropy as described by the Bidirectional Reflectance Distribution Function (BRDF) as an additional source of information to improve land surface classification accuracies in a Canadian boreal forest region through the use of a decision tree classifier (C4.5). This effort primarily uses a daily rolling version of the operational algorithm developed for Direct Broadcast to generate 500 m 16-day daily rolling data sets in the study region. Descriptive statistic and statistically rigorous techniques are used to assess classification accuracies based on confusion matrices and a 10-fold cross-validation method. The results show that the inclusion of additional 7-band model anisotropic parameter group (volumetric (VOL) plus geometric (GEO)) with spectral feature group (nadir BRDF-adjusted reflectance (NABR) plus Enhanced Vegetation Index (EVI)) is most useful in classification, increasing overall accuracies by 5.68%. The most improvements of per-class accuracies are seen for Wetland shrub class with user's and producer's accuracies increasing by up to 17.7% and 11.3%, respectively. Increases on the order of 5% to 15% are encountered for the classes of Wetland herb, Wetland tree, Coniferous dense, and Coniferous open with no detriments to other candidate classes. The inclusion of the 2-band BRDF shape indicator group in the classification is, however, not as useful as inclusion of the 7-band model anisotropic parameter group in improving the classification accuracies. A further investigation of the classification accuracies regarding reflectance anisotropy for the sampling pixels within each class shows that land cover types that are dominated by geometric-optical scattering type or a mixture of scattering types are relatively difficult to be classified with spectral feature group alone, and the inclusion of additional BRDF features can significantly improve classification accuracies for these land cover types. However, despite their use as ancillary data, this study also confirms that the spectral feature group provided with NBAR and EVI captures the major information content regarding land cover types, exceeding the information content contained in the model anisotropic parameter group provided with the 7-band VOL and GEO parameters of RossThick-LiSparse-Reciprocal (RTLSR) models.