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

An Approach to Retrieve BRDF from Satellite and Airborne Measurements of Surface-Reflected Radiance Based on Decoupling of Atmospheric Radiative Transfer and Surface Reflection

Bi-directional Reflection Distribution Function (BRDF) defines anisotropy of the surface reflection. It is required to specify the boundary condition for radiative transfer (RT) modeling. Measurements of reflected radiance by satellite- and air-borne sensors provide information about anisotropy of surface reflection. Atmospheric correction needs to be performed to derive BRDF from the reflected radiance. Common approach for BRDF retrievals consists of the use of kernel-based BRDF and RT modeling that needs to be done anew at every step of the iterative process. The kernels weights are obtained by minimization of the difference between measured and modeled radiance. This study develops a new method of retrieving kernel-based BRDF that requires RT calculations to be done only once. The method employs the exact analytical expression of radiance at any atmospheric level through the solutions of two auxiliary atmosphere-only RT problems and the surface-reflected radiance at the surface level. The latter is related to BRDF and solutions of the auxiliary RT problems by a Fredholm integral equation of the second kind. The approach requires to perform RT calculations one time before the iterations. It can use observations taken at different atmospheric conditions assuming that surface conditions remain unchanged during the time span of observations. The algorithm accurately catches zero weights of the kernels that may be a concern if the number of kernels is greater than 3 in current mainstream approaches. The study presents numerical tests of the BRDF retrieval algorithm for various surface and atmospheric conditions

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

Earth observations from DSCOVR EPIC instrument

The National Oceanic and Atmospheric Administration (NOAA) Deep Space Climate Observatory (DSCOVR) spacecraft was launched on 11 February 2015 and in June 2015 achieved its orbit at the first Lagrange point (L1), 1.5 million km from Earth toward the sun. There are two National Aeronautics and Space Administration (NASA) Earth-observing instruments on board: the Earth Polychromatic Imaging Camera (EPIC) and the National Institute of Standards and Technology Advanced Radiometer (NISTAR). The purpose of this paper is to describe various capabilities of the DSCOVR EPIC instrument. EPIC views the entire sunlit Earth from sunrise to sunset at the backscattering direction (scattering angles between 168.5° and 175.5°) with 10 narrowband filters: 317, 325, 340, 388, 443, 552, 680, 688, 764, and 779 nm. We discuss a number of preprocessing steps necessary for EPIC calibration including the geolocation algorithm and the radiometric calibration for each wavelength channel in terms of EPIC counts per second for conversion to reflectance units. The principal EPIC products are total ozone (O3) amount, scene reflectivity, erythemal irradiance, ultraviolet (UV) aerosol properties, sulfur dioxide (SO2) for volcanic eruptions, surface spectral reflectance, vegetation properties, and cloud products including cloud height. Finally, we describe the observation of horizontally oriented ice crystals in clouds and the unexpected use of the O2 B-band absorption for vegetation properties.The NASA GSFC DSCOVR project is funded by NASA Earth Science Division. We gratefully acknowledge the work by S. Taylor and B. Fisher for help with the SO2 retrievals and Marshall Sutton, Carl Hostetter, and the EPIC NISTAR project for help with EPIC data. We also would like to thank the EPIC Cloud Algorithm team, especially Dr. Gala Wind, for the contribution to the EPIC cloud products. (NASA Earth Science Division)Accepted manuscrip

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

ESTIMATING LAND SURFACE ALBEDO FROM SATELLITE DATA

Land surface albedo, defined as the ratio of the surface reflected incoming and outgoing solar radiation, is one of the key geophysical variables controlling the surface radiation budget. Surface shortwave albedo is widely used to drive climate and hydrological models. During the last several decades, remotely sensed surface albedo products have been generated through satellite-acquired data. However, some problems exist in those products due to instrument measurement inaccuracies and the failure of current retrieving procedures, which have limited their applications. More significantly, it has been reported that some albedo products from different satellite sensors do not agree with each other and some even show the opposite long term trend regionally and globally. The emergence of some advanced sensors newly launched or planned in the near future will provide better capabilities for estimating land surface albedo with fine resolution spatially and/or temporally. Traditional methods for estimating the surface shortwave albedo from satellite data include three steps: first, the satellite observations are converted to surface directional reflectance using the atmospheric correction algorithms; second, the surface bidirectional reflectance distribution function (BRDF) models are inverted through the fitting of the surface reflectance composites; finally, the shortwave albedo is calculated from the BRDF through the angular and spectral integration. However, some problems exist in these algorithms, including: 1) "dark-object" based atmospheric correction methods which make it difficult to estimate albedo accurately over non-vegetated or sparsely vegetated area; 2) the long-time composite albedo products cannot satisfy the needs of weather forecasting or land surface modeling when rapid changes such as snow fall/melt, forest fire/clear-cut and crop harvesting occur; 3) the diurnal albedo signature cannot be estimated in the current algorithms due to the Lambertian approximation in some of the atmospheric correction algorithms; 4) prior knowledge has not been effectively incorporated in the current algorithms; and 5) current observation accumulation methods make it difficult to obtain sufficient observations when persistent clouds exist within the accumulation window. To address those issues and to improve the satellite surface albedo estimations, a method using an atmospheric radiative transfer procedure with surface bidirectional reflectance modeling will be applied to simultaneously retrieve land surface albedo and instantaneous aerosol optical depth (AOD). This study consists of three major components. The first focuses on the atmospheric radiative transfer procedure with surface reflectance modeling. Instead of executing atmospheric correction first and then fitting surface reflectance in the previous satellite albedo retrieving procedure, the atmospheric properties (e.g., AOD) and surface properties (e.g., BRDF) are estimated simultaneously to reduce the uncertainties produced in separating the entire radiative transfer process. Data from the Moderate Resolution Imaging Spectroradiometer (MODIS) onboard Terra and Aqua are used to evaluate the performance of this albedo estimation algorithm. Good agreement is reached between the albedo estimates from the proposed algorithm and other validation datasets. The second part is to assess the effectiveness of the proposed algorithm, analyze the error sources, and further apply the algorithm on geostationary satellite - the Spinning Enhanced Visible and InfraRed Imager (SEVIRI) onboard Meteosat Second Generation (MSG). Extensive validations on surface albedo estimations from MSG/SEVIRI observations are conducted based on the comparison with ground measurements and other satellite products. Diurnal changes and day-to-day changes in surface albedo are accurately captured by the proposed algorithm. The third part of this study is to develop a spatially and temporally complete, continuous, and consistent albedo maps through a data fusion method. Since the prior information (or climatology) of albedo/BRDF plays a vital role in controlling the retrieving accuracy in the optimization method, currently available multiple land surface albedo products will be integrated using the Multi-resolution Tree (MRT) models to mitigate problems such as data gaps, systematic bias or low information-noise ratio due to instrument failure, persistent clouds from the viewing direction and algorithm limitations. The major original contributions of this study are as follows: 1) this is the first algorithm for the simultaneous estimations of surface albedo/reflectance and instantaneous AOD by using the atmospheric radiative transfer with surface BRDF modeling for both polar-orbiting and geostationary satellite data; 2) a radiative transfer with surface BRDF models is used to derive surface albedo and directional reflectance from MODIS and SEVIRI observations respectively; 3) extensive validations are made on the comparison between the albedo and AOD retrievals, and the satellite products from other sensors; 4) the slightly modified algorithm has been adopted to be the operational algorithm of Advanced Baseline Imager (ABI) in the future Geostationary Operational Environmental Satellite-R Series (GOES-R) program for estimating land surface albedo; 5) a framework of using MRT is designed to integrate multiple satellite albedo products at different spatial scales to build the spatially and temporally complete, continuous, and consistent albedo maps as the prior knowledge in the retrieving procedure

First operational BRDF, albedo nadir reflectance products from MODIS

With the launch of NASA’s Terra satellite and the MODerate Resolution Imaging Spectroradiometer (MODIS), operational Bidirectional Reflectance Distribution Function (BRDF) and albedo products are now being made available to the scientific community. The MODIS BRDF/Albedo algorithm makes use of a semiempirical kernel-driven bidirectional reflectance model and multidate, multispectral data to provide global 1-km gridded and tiled products of the land surface every 16 days. These products include directional hemispherical albedo (black-sky albedo), bihemispherical albedo (white-sky albedo), Nadir BRDF-Adjusted surface Reflectances (NBAR), model parameters describing the BRDF, and extensive quality assurance information. The algorithm has been consistently producing albedo and NBAR for the public since July 2000. Initial evaluations indicate a stable BRDF/Albedo Product, where, for example, the spatial and temporal progression of phenological characteristics is easily detected in the NBAR and albedo results. These early beta and provisional products auger well for the routine production of stable MODIS-derived BRDF parameters, nadir reflectances, and albedos for use by the global observation and modeling communities

Simulation of multiangular remote sensing products using small satellite formations

To completely capture the multiangular reflectance of an opaque surface, one must estimate the bidirectional reflectance distribution function (BRDF), which seeks to represent variations in surface reflectance as a function of measurement and illumination angles at any time instant. The gap in angular sampling abilities of existing single satellites in Earth observation missions can be complemented by small satellites in formation flight. The formation would have intercalibrated spectrometer payloads making reflectance measurements, at many zenith and azimuthal angles simultaneously. We use a systems engineering tool coupled with a science evaluation tool to demonstrate the performance impact and mission feasibility. Formation designs are generated and compared to each other and multisensor single spacecraft, in terms of estimation error of BRDF and its dependent products such as albedo, light use efficiency (LUE), and normalized difference vegetation index (NDVI). Performance is benchmarked with respect to data from previous airborne campaigns (NASA's Cloud Absorption Radiometer), and tower measurements (AMSPEC II), and assuming known BRDF models. Simulations show that a formation of six small satellites produces lesser average error (21.82%) than larger single spacecraft (23.2%), purely in terms of angular sampling benefits. The average monolithic albedo error of 3.6% is outperformed by a formation of three satellites (1.86%), when arranged optimally and by a formation of seven to eight satellites when arranged in any way. An eight-satellite formation reduces albedo errors to 0.67% and LUE errors from 89.77% (monolithic) to 78.69%. The average NDVI for an eight satellite, nominally maintained formation is better than the monolithic 0.038