Bidirectional Reflectance Distribution Function Effects Research Articles

Using Geostationary Satellite Observations to Improve the Monitoring of Vegetation Phenology

Geostationary satellite data enable frequent observations of the Earth’s surface, facilitating the rapid monitoring of land covers and changes. However, optical signals over vegetation, represented by the vegetation index (VI), exhibit an anisotropic effect due to the diurnal variation in the solar angle during data acquisition by geostationary satellites. This effect, typically characterized by the bi-directional reflectance distribution function (BRDF), can introduce uncertainties in vegetation monitoring and the estimation of phenological transition dates (PTDs). To address this, we investigated the diurnal variation in the normalized difference vegetation index (NDVI) with solar angles obtained from geostationary satellites since the image had fixed observation angles. By establishing a temporal conversion relationship between instantaneous NDVI and daily NDVI at the local solar noon (LSNVI), we successfully converted NDVIs obtained at any time during the day to LSNVI, increasing cloud-free observations of NDVI by 34%. Using different statistics of the time series vegetation index, including LSNVI, daily averaged NDVI (DAVI), and angular corrected NDVI (ACVI), we extracted PTD at five typical sites in China. The results showed a difference of up to 41.5 days in PTD estimation, with the highest accuracy achieved using LSNVI. The use of the proposed conversion approach, utilizing time series LSNVI, reduced the root mean square error (RMSE) of PTD estimation by 9 days compared with the use of actual LSNVI. In conclusion, this study highlights the importance of eliminating BRDF effects in geostationary satellite observations and demonstrates that the proposed angular normalization method can enhance the accuracy of time series NDVI in vegetation monitoring.

Radiometric Improvement of Spectral Indices Using Multispectral Lightweight Sensors Onboard UAVs

Close-range remote sensing techniques employing multispectral sensors on unoccupied aerial vehicles (UAVs) offer both advantages and drawbacks in comparison to traditional remote sensing using satellite-mounted sensors. Close-range remote sensing techniques have been increasingly used in the field of precision agriculture. Planning the flight, including optimal flight altitudes, can enhance both geometric and temporal resolution, facilitating on-demand flights and the selection of the most suitable time of day for various applications. However, the main drawbacks stem from the lower quality of the sensors being used compared to satellites. Close-range sensors can capture spectral responses of plants from multiple viewpoints, mitigating satellite remote sensing challenges, such as atmospheric interference, while intensifying issues such as bidirectional reflectance distribution function (BRDF) effects due to diverse observation angles and morphological variances associated with flight altitude. This paper introduces a methodology for achieving high-quality vegetation indices under varied observation conditions, enhancing reflectance by selectively utilizing well-geometry vegetation pixels, while considering factors such as hotspot, occultation, and BRDF effects. A non-parametric ANOVA analysis demonstrates significant statistical differences between the proposed methodology and the commercial photogrammetric software AgiSoft Metashape, in a case study of a vineyard in Fuente-Alamo (Albacete, Spain). The BRDF model is expected to substantially improve vegetation index calculations in comparison to the methodologies used in satellite remote sensing and those used in close-range remote sensing.

The influence of BRDF effects and representativeness of training data on tree species classification using multi-flightline airborne hyperspectral imagery

Airborne hyperspectral imagery provides detailed spatial and spectral information that can be used for tree species classification, but its application in practical forest inventories over large areas is limited. This study investigates the impact of (i) bidirectional reflectance distribution function (BRDF) effects on multi-flightline airborne imagery and (ii) the selection strategy for training data on tree species classification accuracy over large forested areas. Hyperspectral imagery was acquired using the Chinese Academy of Forestry’s LiDAR, CCD, and Hyperspectral system (CAF-LiCHy) over the Mengjiagang Forest Farm (MFF, ∼155 km2) in northeastern China, which includes common keystone coniferous species in northern China’s temperate-boreal ecozones. A kernel-driven BRDF model based on mNDVI704.50 stratification was applied to airborne hyperspectral data acquired from multi-flightlines. We developed three training data selection strategies, namely #1 sampling over the entire area, #2 localized sampling over the area for each flight date, and #3 gradually decreasing sampling in size. Additionally, two validation assessment strategies were respectively employed: independent data over the entire area and independent data in areas with inconsistent classification labeling between using BRDF corrected and uncorrected hyperspectral imagery. These strategies were implemented to perform a comprehensive cross-comparison of tree species classification accuracy using BRDF corrected and uncorrected airborne hyperspectral imagery. The study found that the representativeness of the training data (in terms of solar-viewing geometry distribution and size) has a stronger impact on the classification accuracy of tree species than the BRDF effects in airborne hyperspectral imagery. By utilizing BRDF corrected hyperspectral imagery with a majority of representative training data, the mean overall accuracy of tree species classification improves by ≤ 2.9 % compared to uncorrected imagery. Notably, BRDF correction has improved classification accuracy by 11.5 % in inconsistently labeled areas, mainly concentrated in the spatial overlap regions of adjacent flightlines. These findings provide valuable insights for improving tree species classification accuracy in large forested areas using multi-flightline airborne hyperspectral imagery.

Radiometric Cross-Calibration of Wide-Field-of-View Cameras Based on Gaofen-1/6 Satellite Synergistic Observations Using Landsat-8 Operational Land Imager Images: A Solution for Off-Nadir Wide-Field-of-View Associated Problems

The Gaofen-1 satellite is equipped with four wide-field-of-view (WFV) instruments, enabling an impressive spatial resolution of 16 m and a combined swath exceeding 800 km. These WFV images have shown their valuable applications across diverse fields. However, achieving accurate radiometric calibration is an essential prerequisite for establishing reliable connections between satellite signals and biophysical, as well as biochemical, parameters. However, observations with large viewing angles (>20°) pose new challenges due to the bidirectional reflectance distribution function (BRDF) effect having a pronounced impact on the accuracy of cross-radiation calibrations, especially for the off-nadir WFV1 and WFV4 cameras. To overcome this challenge, a novel approach was introduced utilizing the combined observations from the Gaofen-1 and Gaofen-6 satellites, with Landsat-8 OLI serving as a reference sensor. The key advantage of this synergistic observation strategy is the ability to obtain a greater number of image pairs that closely resemble Landsat-8 OLI reference images in terms of geometry and observation dates. This increased availability of matching images ensures a more representative dataset of the observation geometry, enabling the derived calibration coefficients to be applicable across various sun–target–sensor geometries. Then, the geometry angles and bidirectional reflectance information were put into a Particle Swarm Optimization (PSO) algorithm incorporating radiative transfer modeling. This PSO-based approach formulates cross-calibration as an optimization problem, eliminating the reliance on complex BRDF models and satellite-based BRDF products that can be affected by cloud contamination. Extensive validation experiments involving satellite data and in situ measurements demonstrated an average uncertainty of less than eight percent for the proposed cross-radiation calibration scheme. Comparisons of top-of-atmosphere (TOA) results calibrated using our proposed scheme, the previous traditional radiative transfer modeling using MODIS BRDF products for BRDF correction (RTM-BRDF) method, and official coefficients reveal the superior accuracy of our method. The proposed scheme achieves a 36.99% decrease in root mean square error (RMSE) and a 38.13% increase in mean absolute error (MAE) compared to official coefficients. Moreover, it achieves comparable accuracy to the RTM-BRDF method while eliminating the need for MODIS BRDF products, with a decrease in RMSE exceeding 14% for the off-nadir WFV1 and WFV4 cameras. The results substantiate the efficacy of the proposed scheme in enhancing cross-calibration accuracy by improving image match-up selection, efficiently removing BRDF effects, and expanding applicability to diverse observation geometries.

An unmixing-based BRDF correction in spectral remote sensing data

While extensive research addressed the Bidirectional Reflectance Distribution Function (BRDF) effect and approaches to correct it, too few works have considered the influence of spectral mixture on the correction results.This work studies the BRDF effect in spectral data and presents an approach to correct its undesired impact, considering the likely presence of mixed pixels. We propose an unmixing-based semiempirical model, incorporating endmembers’ (EMs) fractions within the data correction.To evaluate the performance of the proposed methodology, we conducted experiments with laboratory and aerial hyperspectral image data. The outcomes from all experiments reveal vital insights into the influence of the spectral mixture on the BRDF correction accuracy. Nevertheless, the results clearly show that the accuracy of the corrected reflectance significantly improved using the proposed model for reducing the BRDF effect, regardless of the pixel’s microtopography arrangement.Most importantly, a quantitative assessment of the difference in reflectance values within the overlapping region between two aerial images shows that the unmixing-based model outperforms the commonly used one.While the traditional method reduces the influence of the BRDF in the corrected data by a factor of two, with an average Mean Absolute Error (MAE) of≈2.5%, the proposed approach reduces it by a factor of four with an average MAE≈1.2%.

Temporal Variation in Surface Bidirectional Reflectance of the Railroad Valley Vicarious Calibration Test Site in Nevada

Spectral reflectance-based vicarious calibration (VicCal) requires accurate characterization of the bidirectional reflectance distribution function (BRDF) of the ground-based target. Railroad Valley (RRV) Playa, Nevada, has been used as a VicCal test site since 1995 as it is large, appears stable over time, and has a reasonably consistent surface. This study presents the results of a diurnal measurement cycle that closely replicated illumination geometries for Earth Observing (EO) satellites over a year. By measuring the rate of change of the BRDF with respect to time, we recorded the range of BRDF effects while holding the surface constant with respect to moisture and surface condition variation. The rate of spectral reflectance change increased rapidly with view angle in the backscatter direction, reaching rates of change that are 2.3 and 10.5 times greater in the backscatter than in the forward scatter direction for view angles of 20° and 40°, respectively. The results show that larger off-nadir viewing angles in the backscatter direction are particularly sensitive to changes in solar/view geometries.

MULTIVIEW ANALYSIS OF MIXED PIXELS IN THE FRACTION AND REFLECTANCE DOMAINS FOR UNDERSTANDING SUB-PIXEL TOPOGRAPHIC STRUCTURE

Abstract. The spectral mixture analysis (SMA) plays a vital role in spectral data analysis and extraction of subpixel information. However, this technique provides only quantitative information regarding the materials’ abundance fractions within the pixel. On the other hand, the Bidirectional Reflectance Distribution Function (BRDF) indicates that sub-pixel topography affects the surface’s directional reflection to a large extent. Unfortunately, despite the high importance of the BRDF effect and the SMA in remote sensing, only very few research works addressed their mutual influence. Thus, in this work, we propose a study that addresses this mutual influence and suggests an approach for extracting sub-pixel topographic information from mixed pixels. For this purpose, we conducted two multiview imaging experiments under controlled conditions using artificial mixed surfaces. Each surface type is made of two materials and has a varying structural pattern. Then we measured the BiConical Reflectance Factor (BCRF) of each surface from various viewing zenith angles. Next, we applied spectral unmixing to estimate the abundance fraction of three endmembers (EMs) in each surface’s pattern. Finally, we tested the relationship between the sup-pixel topography and the fraction variation vs. the multiple imaging directions. The first experiment results showed that multiview spectral measurements allow the separability between surfaces combining the same materials’ composition but with different sub-pixel structural arrangements.Moreover, such separability is more accurate in the fraction space than in reflectance space. Besides, and most importantly, the second experiment revealed exciting outcomes regarding the relationship between the sub-pixel topographic feature and the variation of the EM fraction vs. the imaging viewing direction. Specifically, we showed a high correlation between the EMs’ fractions and the height of a repetitive element within the sub-pixel topography with a determination coefficient that reaches 0.89.

IMAGE-BASED BRDF MEASUREMENT

Abstract. One of the most challenging effects of remote sensing is landcover materials' Bidirectional Reflectance Distribution Function (BRDF). A wide range of approaches and measuring methods address the BRDF in various studies. However, there is a requirement for an accurate measurement setup and costly special equipment. Furthermore, the measurements and calculations are applied to model the BRDF for a single point on the object's surface. Considering these limitations, we propose a new modular framework and methodology for measuring, modeling, and analyzing the BRDF without the need for unique instruments. Instead, we suggest acquiring multiple overlapping images in a simple and time-saving way, sampling the desired object's Region Of Interest (ROI) in one image and automatically tracking it in the other images. Experimental results using laboratory data acquired under controlled conditions clearly show the advantages of our framework in retrieving the camera positions, tracking ROIs in the different images, and accurately measuring the BRDF of various land-cover types. Moreover, we observed the variability of the obtained measurements before and after applying the kernel-driven approach to minimize the BRDF effect. The results show that the applied correction reduces this variability significantly, indicating the high accuracy of measuring the directional reflectance using the proposed approach.

Semi-Empirical Models for the Bidirectional Water-Leaving Radiance: An Analysis of a Turbid Inland Lake

When a precise quantitative analysis of satellite measurements over bodies of water is required, the bidirectional effects of water-leaving radiance must be considered. The bidirectional reflectance distribution function (BRDF) is used to estimate the directional dependency of the radiance. Previous research on BRDF has focused on oceanic waters; few studies on turbid inland waters have been conducted. In this article, using multi-angle MISR measurements, five semi-empirical BRDF models (MAG2002, Lee2004, Park-Ruddick 2005, Woerd-Pasterkamp2008, and Lee2011) were quantitatively compared, and an adaptive algorithm was proposed over a typical turbid inland lake, Taihu Lake, China. Our results reveal the following: 1) the Woerd-Pasterkamp2008 and Lee2011 models provide the best fits with correlation coefficients greater than 0.8; 2) when prior modeling parameters were used, the Lee2011 model was still the most accurate with RMSEs less than 1.1%, while the accuracy of the Woerd-Pasterkamp2008 model varied; and 3) the use of an adaptive algorithm including an empirical rule based on the ratio ofbb/aimproved the accuracy. The results provide a theoretical basis for BRDF models and BRDF effects over inland Case II waters. They also providea prioriknowledge for future studies on water constituents and the quantitative inversion of atmospheric parameters.

A Semi-Empirical Anisotropy Correction Model for UAS-Based Multispectral Images of Bare Soil

The recent developments in the performance and miniaturization of uncrewed aircraft systems (UAS) and multispectral imaging sensors provide new tools for the assessment of the spatial and temporal variability of soil properties at sub-meter resolution and at relatively low costs, in comparison to traditional chemical analysis. The accuracy of multispectral data is nevertheless influenced by the anisotropic behaviour of natural surfaces, framed in the general theory of the bidirectional reflectance distribution function (BRDF). Accounting for BRDF effects in multispectral data is paramount before formulating any scientific interpretation. This study presents a semi-empirical spectral normalization methodology for UAS-based multispectral imaging datasets of bare soils to account for the effects of the BRDF, based on the application of an anisotropy factor (ANIF). A dataset of images from 15 flights over bare soil fields in the Belgian loam belt was used to calibrate a model relating the ANIF to a wide range of illumination geometry conditions by using only two angles: relative sensor-pixel-sun zenith and relative sensor-pixel-sun azimuth. The employment of ANIF-corrected images for multispectral orthomosaic generation with photogrammetric software provided spectral maps free of anisotropic-related artefacts in most cases, as assessed by several ad hoc indexes, and was also tested on an independent validation set. Most notably, the standard deviation in the measured reflectance of the same georeferenced point by different pictures decreased from 0.032 to 0.023 (p < 0.05) in the calibration dataset and from 0.037 to 0.030 in the validation dataset. The validation dataset, however, showed the presence of some systematic errors, the causes of which require further investigation.

FlexBRDF: A Flexible BRDF Correction for Grouped Processing of Airborne Imaging Spectroscopy Flightlines.

Bidirectional reflectance distribution function (BRDF) effects are a persistent issue for the analysis of vegetation in airborne imaging spectroscopy data, especially when mosaicking results from adjacent flightlines. With the advent of large airborne imaging efforts from NASA and the U.S. National Ecological Observatory Network (NEON), there is increasing need for methods that are flexible and automatable across images with diverse land cover. Flexible bidirectional reflectance distribution function (FlexBRDF) is built upon the widely used kernel method, with additional features including stratified random sampling across flightline groups, dynamic land cover stratification by normalized difference vegetation index (NDVI), interpolation of correction coefficients across NDVI bins, and the use of a reference solar zenith angle. We demonstrate FlexBRDF using nine long (150–400 km) airborne visible/infrared imaging spectrometer (AVIRIS)‐Classic flightlines collected on 22 May 2013 over Southern California, where diverse land cover and a wide range of solar illumination yield significant BRDF effects. We further test the approach on additional AVIRIS‐Classic data from California, AVIRIS‐Next Generation data from the Arctic and India, and NEON imagery from Wisconsin. Comparison of overlapping areas of flightlines show that models built from multiple flightlines performed better than those built for single images (root mean square error improved up to 2.3% and mean absolute deviation 2.5%). Standardization to a common solar zenith angle among a flightline group improved performance, and interpolation across bins minimized between‐bin boundaries. While BRDF corrections for individual sites suffice for local studies, FlexBRDF is an open source option that is compatible with bulk processing of large airborne data sets covering diverse land cover needed for calibration/validation of forthcoming spaceborne imaging spectroscopy missions.

Airborne Hyperspectral Data Acquisition and Processing in the Arctic: A Pilot Study Using the Hyspex Imaging Spectrometer for Wetland Mapping

A pilot study for mapping the Arctic wetlands was conducted in the Yukon Flats National Wildlife Refuge (Refuge), Alaska. It included commissioning the HySpex VNIR-1800 and the HySpex SWIR-384 imaging spectrometers in a single-engine Found Bush Hawk aircraft, planning the flight times, direction, and speed to minimize the strong bidirectional reflectance distribution function (BRDF) effects present at high latitudes and establishing improved data processing workflows for the high-latitude environments. Hyperspectral images were acquired on two clear-sky days in early September, 2018, over three pilot study areas that together represented a wide variety of vegetation and wetland environments. Steps to further minimize BRDF effects and achieve a higher geometric accuracy were added to adapt and improve the Hyspex data processing workflow, developed by the German Aerospace Center (DLR), for high-latitude environments. One-meter spatial resolution hyperspectral images, that included a subset of only 120 selected spectral bands, were used for wetland mapping. A six-category legend was established based on previous U.S. Geological Survey (USGS) and U.S. Fish and Wildlife Service (USFWS) information and maps, and three different classification methods—hybrid classification, spectral angle mapper, and maximum likelihood—were used at two selected sites. The best classification performance occurred when using the maximum likelihood classifier with an averaged Kappa index of 0.95; followed by the spectral angle mapper (SAM) classifier with a Kappa index of 0.62; and, lastly, by the hybrid classifier showing lower performance with a Kappa index of 0.51. Recommendations for improvements of future work include the concurrent acquisition of LiDAR or RGB photo-derived digital surface models as well as detailed spectra collection for Alaska wetland cover to improve classification efforts.

A Dual-Channel Aerosol Optical Depth Retrieval Algorithm Incorporating the BRDF Effect from AVHRR over Eastern Asia

A NOAA/AVHRR dual-channel method over land is proposed to simultaneously retrieve aerosol optical depth (AOD) at 0.55 μm, and surface reflectance at 0.63 and 0.85 μm. Compared with previous well-established one-channel retrieval algorithms, this algorithm takes advantage of the surface reflectance ratio between 0.63 and 0.85 μm in an attempt to account for the effect induced by the surface bidirectional reflectance distribution function (BRDF). This effect cannot be negligible due to the orbit drift and phasing running of NOAA satellites, both of which potentially cause large solar angular variation. Meanwhile, the observation posture change of AVHRR would cause large sensor angular variation in time series measurements. The used surface reflectance ratio based on dual channels at 0.63 and 0.85 μm is found more reasonable to be assumed as unchanged during a certain period of time, compared to the traditional ratio when addressing the BRDF issue. AOD retrievals have been carried out over Eastern Asia. Validation against aerosol robotic network (AERONET) measurements shows that up to 83% of AOD validation collocations are within error lines (±0.15 ± 0.15τ, τ is AOD) with an R of 0.88 and an root mean square error (RMSE) of 0.15. The dual-channel algorithm taking into account the surface BRDF effect is proved outperforming the conventional 0.63 μm-channel method. It indicates that our algorithm has the potential to be applied to the retrieval of long series of AOD, especially to the AOD retrieval of the sensors which lack a shortwave infrared channel required in the MODerate resolution Imaging Spectroradiometer (MODIS) dark target AOD algorithm.

Ozone Monitoring Instrument (OMI) Aura nitrogen dioxide standard product version 4.0 with improved surface and cloud treatments

Abstract. We present a new and improved version (V4.0) of the NASA standard nitrogen dioxide (NO2) product from the Ozone Monitoring Instrument (OMI) on the Aura satellite. This version incorporates the most salient improvements for OMI NO2 products suggested by expert users and enhances the NO2 data quality in several ways through improvements to the air mass factors (AMFs) used in the retrieval algorithm. The algorithm is based on the geometry-dependent surface Lambertian equivalent reflectivity (GLER) operational product that is available on an OMI pixel basis. GLER is calculated using the vector linearized discrete ordinate radiative transfer (VLIDORT) model, which uses as input high-resolution bidirectional reflectance distribution function (BRDF) information from NASA's Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) instruments over land and the wind-dependent Cox–Munk wave-facet slope distribution over water, the latter with a contribution from the water-leaving radiance. The GLER combined with consistently retrieved oxygen dimer (O2–O2) absorption-based effective cloud fraction (ECF) and optical centroid pressure (OCP) provide improved information to the new NO2 AMF calculations. The new AMFs increase the retrieved tropospheric NO2 by up to 50 % in highly polluted areas; these differences arise from both cloud and surface BRDF effects as well as biases between the new MODIS-based and previously used OMI-based climatological surface reflectance data sets. We quantitatively evaluate the new NO2 product using independent observations from ground-based and airborne instruments. The new V4.0 data and relevant explanatory documentation are publicly available from the NASA Goddard Earth Sciences Data and Information Services Center (https://disc.gsfc.nasa.gov/datasets/OMNO2_V003/summary/, last access: 8 November 2020), and we encourage their use over previous versions of OMI NO2 products.

Effects of clouds on the UV Absorbing Aerosol Index from TROPOMI

Abstract. The ultraviolet (UV) Absorbing Aerosol Index (AAI) is widely used as an indicator for the presence of absorbing aerosols in the atmosphere. Here we consider the TROPOMI AAI based on the 340 nm/380 nm wavelength pair. We investigate the effects of clouds on the AAI observed at small and large scales. The large-scale effects are studied using an aggregate of TROPOMI measurements over an area mostly devoid of absorbing aerosols (Pacific Ocean). The study reveals that several structural features can be distinguished in the AAI, such as the cloud bow, viewing zenith angle dependence, sunglint, and a previously unexplained increase in AAI values at extreme viewing and solar geometries. We explain these features in terms of the bidirectional reflectance distribution function (BRDF) of the scene in combination with the different ratios of diffuse and direct illumination of the surface at 340 and 380 nm. To reduce the dependency on the BRDF and homogenize the AAI distribution across the orbit, we present three different AAI retrieval models: the traditional Lambertian scene model (LSM), a Lambertian cloud model (LCM), and a scattering cloud model (SCM). We perform a model study to assess the propagation of errors in auxiliary databases used in the cloud models. The three models are then applied to the same low-aerosol region. Results show that using the LCM and SCM gives on average a higher AAI than the LSM. Additionally, a more homogeneous distribution is retrieved across the orbit. At the small scale, related to the high spatial resolution of TROPOMI, strong local increases and decreases in AAI are observed in the presence of clouds. The BRDF effect presented here is a first step – more research is needed to explain the small-scale cloud effects on the AAI.

Complete and accurate data correction for seamless mosaicking of airborne hyperspectral images: A case study at a mining site in Inner Mongolia, China

Airborne hyperspectral remote sensing is an important application in the ecological monitoring of the environment in mining areas, and accurate preprocessing of the original images is the key to quantitative information retrieval. The original image data need radiation correction to acquire surface reflectance data. Due to the impact of the field angle, incidental radiance, and the bidirectional reflectance distribution function (BRDF), there can be a brightness gradient between adjacent strips, which leads to radiance difference and obvious chromatic aberration of the mosaicked images. We propose a novel data correction method for seamless mosaicking of airborne hyperspectral images. Firstly, visible and near-infrared (VNIR) and shortwave infrared (SWIR) sensors are calibrated in the laboratory, and the radiation calibration model of the sensor is established by an integrating-sphere system. A correction function is then established by combining the BRDF effect and the radiation attenuation coefficients. We also normalize the exposure time, sun altitude angle, and sensor altitude angle according to the flight strip. The results showed that this method is able to eliminate the signal distortion, allowing the seamless mosaicking of 37 strip images which were taken in different date and conditions in the study area. After the atmospheric correction of the imagery was completed, the accuracy of the preprocessing results was evaluated by field-measured ASD spectroradiometer data. The coefficient of determination R2 of the results for the reflectance was greater than 0.9. The experiments show that the proposed method has a good performance in radiation accuracy, and can provide high-quality hyperspectral data for the follow-up application of the ecological monitoring of a mining area.

A Kernel-Driven BRDF Approach to Correct Airborne Hyperspectral Imagery over Forested Areas with Rugged Topography

Airborne hyper-spectral imaging has been proven to be an efficient means to provide new insights for the retrieval of biophysical variables. However, quantitative estimates of unbiased information derived from airborne hyperspectral measurements primarily require a correction of the anisotropic scattering properties of the land surface depicted by the bidirectional reflectance distribution function (BRDF). Hitherto, angular BRDF correction methods rarely combined viewing-illumination geometry and topographic information to achieve a comprehensive understanding and quantification of the BRDF effects. This is in particular the case for forested areas, frequently underlaid by rugged topography. This paper describes a method to correct the BRDF effects of airborne hyperspectral imagery over forested areas overlying rugged topography, referred in the reminder of the paper as rugged topography-BRDF (RT-BRDF) correction. The local viewing and illumination geometry are calculated for each pixel based on the characteristics of the airborne scanner and the local topography, and these two variables are used to adapt the Ross-Thick-Maignan and Li-Transit-Reciprocal kernels in the case of rugged topography. The new BRDF model is fitted to the anisotropy of multi-line airborne hyperspectral data. The number of pixels is set at 35,000 in this study, based on a stratified random sampling method to ensure a comprehensive coverage of the viewing and illumination angles and to minimize the fitting error of the BRDF model for all bands. Based on multi-line airborne hyperspectral data acquired with the Chinese Academy of Forestry’s LiDAR, CCD, and Hyperspectral system (CAF-LiCHy) in the Pu’er region (China), the results applying the RT-BRDF correction are compared with results from current empirical (C, and sun-canopy-sensor (SCS) adds C (SCS+C)) and semi-physical (SCS) topographic correction methods. Both quantitative assessment and visual inspection indicate that RT-BRDF, C, and SCS+C correction methods all reduce the topographic effects. However, the RT-BRDF method appears more efficient in reducing the variability in reflectance of overlapping areas in multiple flight-lines, with the advantage of reducing the BRDF effects caused by the combination of wide field of view (FOV) airborne scanner, rugged topography, and varying solar illumination angle over long flight time. Specifically, the average decrease in coefficient of variation (CV) is 3% and 3.5% for coniferous forest and broadleaved forest, respectively. This improvement is particularly marked in the near infrared (NIR) region (i.e., >750 nm). This finding opens new possible applications of airborne hyperspectral surveys over large areas.

Assessment of bidirectional reflectance effects on desert and forest for radiometric cross-calibration of satellite sensors

Cross-calibration is a wildly used approach to radiometrically calibrate satellite sensors to meet uncertainty requirements. However, its accuracy is greatly influenced by bidirectional reflectance distribution function (BRDF) effects. To understand such influences, we analyzed the long-term BRDF characteristics of two commonly used cross-calibration targets (‘bright’ desert and ‘dark’ forest) using long-term (2002 – 2016) Moderate Resolution Imaging Spectroradiometer (MODIS) BRDF products. BRDF characteristics derived by MODIS observations of both land cover types had 365-day cycles. Furthermore, both solar and sensor zenith angles showed more influence on BRDF than relative azimuth angles, and BRDF variations were approximately 10 times sensitive when solar or sensor zenith angles surpassed the boundary of 75°, which is therefore not recommended for selecting image pairs between target and reference satellite sensors for cross-calibration. Due to cloud cover, three prevailing interpolation methods (nearest, linear, and spline interpolation) were also assessed in their capacity to fill temporally missing BRDF products. Linear interpolation was preferred when valid MODIS scenes were more than two in five days, while nearest interpolation outperformed the other two methods for no more than two valid MODIS scenes in five days. In addition, BRDF look up tables (LUTs) were established for radiometrically cross-calibrating satellite sensors when BRDF products were unavailable. The LUTs were further validated with Landsat-8/Operational Land Imager (OLI) top of atmosphere (TOA) radiance simulated by the Second Simulation of the Satellite Signal in the Solar Spectrum model (6S) with BRDF effects corrected using the BRDF LUTs and MODIS BRDF products for typical geometries at two sites. High consistency was achieved with mean biases less than 0.3%, which suggests that BRDF LUTs could be used as alternatives when BRDF products are unavailable. These results provide guidelines to improve the accuracy of cross-calibration including image-matchup selecting criteria, efficient BRDF effect removal, and selection of potential calibration sites.

Radiometric Top-of-Atmosphere Reflectance Consistency Assessment for Landsat 8/OLI, Sentinel-2/MSI, PROBA-V, and DEIMOS-1 over Libya-4 and RadCalNet Calibration Sites

There is a clear trend toward the use of higher spatial resolution satellite sensors. Due to the low revisit time of these sensors and frequent cloud coverage, many applications require data from different sensors to be combined in order to have more frequent observations. This raises concerns regarding data interoperability and consistency. The initial pre-requisite is that there are no radiometric differences in top-of-atmosphere (TOA) observations. This paper aims to quantitatively assess differences in the TOA signal provided by PROBA-V, Sentinel-2A and Sentinel-2B, Landsat-8, and Deimos-1 by using observations over both the Libya-4 desert calibration site and the RadCalNet sites. The results obtained over the Libya-4 site indicate that for all sensors investigated, the inter-sensor deviations are negligible, i.e., within ±2% for comparable spectral bands, with the exception of the Deimos-1 Green band. Clear BRDF (bi-directional reflectance distribution function) effects were observed over the RadCalNet sites, thereby preventing consistent conclusions on inter-sensor deviations from being made. In order to fully explore the potential of the RadCalNet sites, it is recommended that BRDF characterizations be additionally incorporated into the RadCalNet simulations and made publicly available through the distribution portal.

Evaluating the Performance of Satellite-Derived Vegetation Indices for Estimating Gross Primary Productivity Using FLUXNET Observations across the Globe

Satellite-derived vegetation indices (VIs) have been widely used to approximate or estimate gross primary productivity (GPP). However, it remains unclear how the VI-GPP relationship varies with indices, biomes, timescales, and the bidirectional reflectance distribution function (BRDF) effect. We examined the relationship between VIs and GPP for 121 FLUXNET sites across the globe and assessed how the VI-GPP relationship varied among a variety of biomes at both monthly and annual timescales. We used three widely-used VIs: normalized difference vegetation index (NDVI), enhanced vegetation index (EVI), and 2-band EVI (EVI2) as well as a new VI - NIRV and used surface reflectance both with and without BRDF correction from the moderate resolution imaging spectroradiometer (MODIS) to calculate these indices. The resulting traditional (NDVI, EVI, EVI2, and NIRV) and BRDF-corrected (NDVIBRDF, EVIBRDF, EVI2BRDF, and NIRV, BRDF) VIs were used to examine the VI-GPP relationship. At the monthly scale, all VIs were moderate or strong predictors of GPP, and the BRDF correction improved their performance. EVI2BRDF and NIRV, BRDF had similar performance in capturing the variations in tower GPP as did the MODIS GPP product. The VIs explained lower variance in tower GPP at the annual scale than at the monthly scale. The BRDF-correction of surface reflectance did not improve the VI-GPP relationship at the annual scale. The VIs had similar capability in capturing the interannual variability in tower GPP as MODIS GPP. VIs were influenced by temperature and water stresses and were more sensitive to temperature stress than to water stress. VIs in combination with environmental factors could improve the prediction of GPP than VIs alone. Our findings can help us better understand how the VI-GPP relationship varies among indices, biomes, and timescales and how the BRDF effect influences the VI-GPP relationship.