Variability in Surface BRDF at Different Spatial Scales (30 m-500 m) Over a Mixed Agricultural Landscape as Retrieved from Airborne and Satellite Spectral Measurements

Over the past decade, the role of multiangle remote sensing has been central to the development of algorithms for the retrieval of global land surface properties including models of the bidirectional reflectance distribution function (BRDF), albedo, land cover/dynamics, burned area extent, as well as other key surface biophysical quantities represented by the anisotropic reflectance characteristics of vegetation. In this study, a new retrieval strategy for fine-to-moderate resolution multiangle observations was developed, based on the operational sequence used to retrieve the Moderate Resolution Imaging Spectroradiometer (MODIS) Collection 5 reflectance and BRDF/albedo products. The algorithm makes use of a semiempirical kernel-driven bidirectional reflectance model to provide estimates of intrinsic albedo (i.e., directional-hemispherical reflectance and bihemispherical reflectance), model parameters describing the BRDF, and extensive quality assurance information. The new retrieval strategy was applied to NASA's Cloud Absorption Radiometer (CAR) data acquired during the 2007 Cloud and Land Surface Interaction Campaign (CLASIC) over the well-instrumented Atmospheric Radiation Measurement Program (ARM) Southern Great Plains (SGP) Cloud and Radiation Testbed (CART) site in Oklahoma, USA. For the case analyzed, we obtained approx.1.6 million individual surface bidirectional reflectance factor (BRF) retrievals, from nadir to 75 off-nadir, and at spatial resolutions ranging from 3 m - 500 m. This unique dataset was used to examine the interaction of the spatial and angular characteristics of a mixed agricultural landscape; and provided the basis for detailed assessments of: (1) the use of a priori knowledge in kernel-driven BRDF model inversions; (2) the interaction between surface reflectance anisotropy and instrument spatial resolution; and (3) the uncertain ties that arise when sub-pixel differences in the BRDF are aggregated to a moderate resolution satellite pixel. Results offer empirical evidence concerning the influence of scale and spatial heterogeneity in kernel-driven BRDF models; providing potential new insights into the behavior and characteristics of different surface radiative properties related to land/use cover change and vegetation structure

Parameterizing anisotropic reflectance of snow surfaces from airborne digital camera observations in Antarctica

The surface reflection of solar radiation comprises an important boundary condition for solar radiative transfer simulations. In polar regions above snow surfaces, the surface reflection is particularly anisotropic due to low Sun elevations and the highly anisotropic scattering phase function of the snow crystals. The characterization of this surface reflection anisotropy is essential for satellite remote sensing over both the Arctic and Antarctica. To quantify the angular snow reflection properties, the hemispherical-directional reflectance factor (HDRF) of snow surfaces was derived from airborne measurements in Antarctica during austral summer in 2013/14. For this purpose, a digital 180∘ fish-eye camera (green channel, 490–585 nm wavelength band) was used. The HDRF was measured for different surface roughness conditions, optical-equivalent snow grain sizes, and solar zenith angles. The airborne observations covered an area of around 1000 km × 1000 km in the vicinity of Kohnen Station (75∘0′ S, 0∘4′ E) at the outer part of the East Antarctic Plateau. The observations include regions with higher (coastal areas) and lower (inner Antarctica) precipitation amounts and frequencies. The digital camera provided upward, angular-dependent radiance measurements from the lower hemisphere. The comparison of the measured HDRF derived for smooth and rough snow surfaces (sastrugi) showed significant differences, which are superimposed on the diurnal cycle. By inverting a semi-empirical kernel-driven bidirectional reflectance distribution function (BRDF) model, the measured HDRF of snow surfaces was parameterized as a function of solar zenith angle, surface roughness, and optical-equivalent snow grain size. This allows a direct comparison of the HDRF measurements with the BRDF derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite product MCD43. For the analyzed cases, MODIS observations (545–565 nm wavelength band) generally underestimated the anisotropy of the surface reflection. The largest deviations were found for the volumetric model weight fvol (average underestimation by a factor of 10). These deviations are likely linked to short-term changes in snow properties

Efficient hyper-parameter determination for regularised linear BRDF parameter retrieval

Linear kernel driven models of the surface Bidirectional Reflectance Distribution Function (BRDF) are valuable tools for exploiting Earth observation data acquired at different sun–sensor geometries. Here we present a method that efficiently determines linear BRDF model weights using Tikhonov smoothing where the smoothing parameter λ is determined via a Generalized Singular Value Decomposition with the root mean square error prescribed depending on the MODIS band. We applied this method to twenty-six different deciduous broadleaf sites across an entire year using the MODIS Terra and Aqua reflectance data products. Kernel weights and white sky albedo derived from this GSVD method were generally consistent with those provided by the MCD43 data products. The GSVD derived results had less sample variability compared to the MCD43 data products, attributable to the assumed smoothness between kernel weights in the Tikhonov smoothing method. The GSVD technique consistently outperforms MCD43 in the reconstruction of observed MODIS reflectance data, of which retrievals from this method will do a better job of estimating albedo and normalizing data to specified geometries

Effect of the aerosol type selection for the retrieval of shortwave ground net radiation: case study using landsat 8 data

This paper discusses the aerosol radiative effects involved in the accuracy of shortwave net radiation, Rn.sw, with sw 2 (400–900) nm, retrieved by the Operational Land Imager (OLI), the new generation sensor of the Landsat mission. Net radiation is a key parameter for the energy exchange between the land and atmosphere; thus, Rn.sw retrieval from space is under investigation by exploiting the increased spatial resolution of the visible and near-infrared OLI data. We adopted the latest version of the Second Simulation of a Satellite Signal in the Solar Spectrum (6SV) atmospheric radiative transfer model implemented in the atmospheric correction algorithm (OLI Atmospherically-Corrected Reflectance Imagery (OLI@CRI)) developed specifically for OLI data. The values of Rn.sw were obtained by varying the microphysical properties of the aerosol during the OLI@CRI retrieval of both the OLI surface reflectance, roli pxl , and the incoming solar irradiance at the surface. The analysis of the aerosol effects on the Rn.sw was carried out on a spectrally-homogeneous desert area located in the southwestern Nile Delta. The results reveal that the Rn.sw available for energy exchange between the land and atmosphere reduces the accuracy (NRMSE ' 14%) when the local aerosol microphysical properties are not considered during the processing of space data. Consequently, these findings suggest that the aerosol type should be considered for variables retrieved by satellite observations concerning the energy exchange in the natural ecosystems, such as Photosynthetically-Active Radiation (PAR). This will also improve the accuracy of land monitoring and of solar energy for power generation when space data are used

Intercomparison of Surface Albedo Retrievals from MISR, MODIS, CGLS Using Tower and Upscaled Tower Measurements

Surface albedo is of crucial interest in land–climate interaction studies, since it is a key parameter that affects the Earth’s radiation budget. The temporal and spatial variation of surface albedo can be retrieved from conventional satellite observations after a series of processes, including atmospheric correction to surface spectral bi-directional reflectance factor (BRF), bi-directional reflectance distribution function (BRDF) modelling using these BRFs, and, where required, narrow-to-broadband albedo conversions. This processing chain introduces errors that can be accumulated and then affect the accuracy of the retrieved albedo products. In this study, the albedo products derived from the multi-angle imaging spectroradiometer (MISR), moderate resolution imaging spectroradiometer (MODIS) and the Copernicus Global Land Service (CGLS), based on the VEGETATION and now the PROBA-V sensors, are compared with albedometer and upscaled in situ measurements from 19 tower sites from the FLUXNET network, surface radiation budget network (SURFRAD) and Baseline Surface Radiation Network (BSRN) networks. The MISR sensor onboard the Terra satellite has 9 cameras at different view angles, which allows a near-simultaneous retrieval of surface albedo. Using a 16-day retrieval algorithm, the MODIS generates the daily albedo products (MCD43A) at a 500-m resolution. The CGLS albedo products are derived from the VEGETATION and PROBA-V, and updated every 10 days using a weighted 30-day window. We describe a newly developed method to derive the two types of albedo, which are directional hemispherical reflectance (DHR) and bi-hemispherical reflectance (BHR), directly from three tower-measured variables of shortwave radiation: downwelling, upwelling and diffuse shortwave radiation. In the validation process, the MISR, MODIS and CGLS-derived albedos (DHR and BHR) are first compared with tower measured albedos, using pixel-to-point analysis, between 2012 to 2016. The tower measured point albedos are then upscaled to coarse-resolution albedos, based on atmospherically corrected BRFs from high-resolution Earth observation (HR-EO) data, alongside MODIS BRDF climatology from a larger area. Then a pixel-to-pixel comparison is performed between DHR and BHR retrieved from coarse-resolution satellite observations and DHR and BHR upscaled from accurate tower measurements. The experimental results are presented on exploring the parameter space associated with land cover type, heterogeneous vs. homogeneous and instantaneous vs. time composite retrievals of surface albedo

Assessing change in the Earth's land surface albedo with moderate resolution satellite imagery

Land surface albedo describes the proportion of incident solar radiant flux that is reflected from the Earth's surface and therefore is a crucial parameter in modeling and monitoring attempts to capture the current climate, hydrological, and biogeochemical cycles and predict future scenarios. Due to the temporal variability and spatial heterogeneity of land surface albedo, remote sensing offers the only realistic method of monitoring albedo on a global scale. While the distribution of bright, highly reflective surfaces (clouds, snow, deserts) govern the vast majority of the fluctuation, variations in the intrinsic surface albedo due to natural and human disturbances such as urban development, fire, pests, harvesting, grazing, flooding, and erosion, as well as the natural seasonal rhythm of vegetation phenology, play a significant role as well. The development of times series of global snow-free and cloud-free albedo from remotely sensed observations over the past decade and a half offers a unique opportunity to monitor and assess the impact of these alterations to the Earth's land surface. By utilizing multiple satellite records from the MODerate-resolution Imaging Spectroradiometer (MODIS), the Multi-angle Imaging Spectroradiometer (MISR) and the Visible Infrared Imaging Radiometer Suite (VIIRS) instruments, and developing innovative spectral conversion coefficients and temporal gap-filling strategies, it has been possible to utilize the strengths of the various sensors to improve the spatial and temporal coverage of global land surface albedo retrievals. The availability of these products is particularly important in tropical regions where cloud cover obscures the forest for significant periods. In the Amazon, field ecologists have noted that some areas of the forest ecosystem respond rapidly with foliage growth at the beginning of the dry season, when sunlight can finally penetrate fully to the surface and have suggested this phenomenon can continue until reductions in water availability (particularly in times of drought) impact the growth cycle. While it has been difficult to capture this variability from individual optical satellite sensors, the temporally gap-filled albedo products developed during this research are used in a case study to monitor the Amazon during the dry season and identify the extent of these regions of foliage growth

A General Method to Normalize Landsat Reflectance Data to Nadir BRDF Adjusted Reflectance

The Landsat satellites have been providing spectacular imagery of the Earth\u27s surface for over 40 years. However, they acquire images at view angles ±7.5° from nadir that cause small directional effects in the surface reflectance. There are also variations with solar zenith angle over the year that can cause apparent change in reflectance even if the surface properties remain constant. When Landsat data from adjoining paths, or from long time series are used, a model of the surface anisotropy is required to adjust all Landsat observations to a uniform nadir view (primarily for visual consistency, vegetation monitoring, or detection of subtle surface changes). Here a generalized approach is developed to provide consistent view angle corrections across the Landsat archive. While this approach is not applicable for generation of Landsat surface albedo, which requires a full characterization of the surface bidirectional reflectance distribution function (BRDF), or for correction to a constant solar illumination angle across a wide range of sun angles, it provides Landsat nadir BRDF-adjusted reflectance (NBAR) for a range of terrestrial monitoring applications. The Landsat NBAR is derived as the product of the observed Landsat reflectance and the ratio of the reflectances modeled using MODIS BRDF spectral model parameters for the observed Landsat and for a nadir view and fixed solar zenith geometry. In this study, a total of 567 conterminous United States (CONUS) January and July 2010 Landsat 5 Thematic Mapper (TM) and Landsat 7 Enhanced Thematic Mapper (ETM+) images that have swath edge overlapping paths sensed in alternating backscatter and forward scattering orientations were used. The average difference between Landsat 5 TM and Landsat 7 ETM+ surface reflectance in the forward and backward scatter directions at the overlapping Landsat scan edges was quantified. The CONUS July view zenith BRDF effects were about 0.02 in the Landsat visible bands, and about 0.03, 0.05 and 0.06, in the 2.1 μm, 1.6 μm and nearinfrared bands respectively. Comparisons of Landsat 5 TM and Landsat 7 ETM+ NBAR derived using MODIS BRDF spectral model parameters defined with respect to different spatial and temporal scales, and defined with respect to different land cover types, were undertaken. The results suggest that, because the BRDF shapes of different terrestrial surfaces are sufficiently similar over the narrow 15° Landsat field of view, a fixed set of MODIS BRDF spectral model parameters may be adequate for Landsat NBAR derivation with little sensitivity to the land cover type, condition, or surface disturbance. A fixed set of BRDF spectral model parameters, derived from a global year of highest quality snow-free MODIS BRDF product values, are provided so users may implement the described Landsat NBAR generation method

Intercomparison of MODIS Albedo Retrievals and In Situ Measurements Across the Global FLUXNET Network

Surface albedo is a key parameter in the Earth's energy balance since it affects the amount of solar radiation directly absorbed at the planet surface. Its variability in time and space can be globally retrieved through the use of remote sensing products. To evaluate and improve the quality of satellite retrievals, careful intercomparisons with in situ measurements of surface albedo are crucial. For this purpose we compared MODIS albedo retrievals with surface measurements taken at 53 FLUXNET sites that met strict conditions of land cover homogeneity. A good agreement between mean yearly values of satellite retrievals and in situ measurements was found (R(exp 2)= 0.82). The mismatch is correlated to the spatial heterogeneity of surface albedo, stressing the relevance of land cover homogeneity when comparing point to pixel data. When the seasonal patterns of MODIS albedo is considered for different plant functional types, the match with surface observation is extremely good at all forest sites. On the contrary, in non-forest sites satellite retrievals underestimate in situ measurements across the seasonal cycle. The mismatch observed at grasslands and croplands sites is likely due to the extreme fragmentation of these landscapes, as confirmed by geostatistical attributes derived from high resolution scenes

Sensitivity of estimated total canopy SIF emission to remotely sensed LAI and BRDF products

Remote sensing of solar-induced chlorophyll fluorescence (SIF) provides new possibilities to estimate terrestrial gross primary production (GPP). To mitigate the angular and canopy structural effects on original SIF observed by sensors (SIFobs), it is recommended to derive total canopy SIF emission (SIFtotal) of leaves within a canopy using canopy interception (i0) and reflectance of vegetation (RV). However, the effects of the uncertainties in i0 and RV on the estimation of SIFtotal have not been well understood. Here, we evaluated such effects on the estimation of GPP using the Soil-Canopy-Observation of Photosynthesis and the Energy balance (SCOPE) model. The SCOPE simulations showed that the R2 between GPP and SIFtotal was clearly higher than that between GPP and SIFobs and the differences in R2 (ΔR2) tend to decrease with the increasing levels of uncertainties in i0 and RV. The resultant ΔR2 decreased to zero when the uncertainty level in i0 and RV was ~30% for red band SIF (RSIF, 683 nm) and ~20% for far-red band SIF (FRSIF, 740 nm). In addition, as compared to the TROPOspheric Monitoring Instrument (TROPOMI) SIFobs at both red and far-red bands, SIFtotal derived using any combination of i0 (from MCD15, VNP15, and CGLS LAI products) and RV (from MCD34, MCD19, and VNP43 BRDF products) showed comparable improvements in estimating GPP. With this study, we suggest a way to advance our understanding in the estimation of a more physiological relevant SIF datasets (SIFtotal) using current satellite products