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

A Method of Retrieving BRDF from Surface Reflected Radiance Using Decoupling of Atmospheric Radiative Transfer and Surface Reflection

BRDF defines anisotropy of the surface reflection. It is required to specify the boundary condition for radiative transfer (RT) modeling used in aerosol retrievals, cloud retrievals, atmospheric modeling and other applications. Ground based measurements of reflected radiance draw increasing attention as a source of information about anisotropy of surface reflection. Derivation of BRDF from surface radiance requires atmospheric correction. This study develops a new method of retrieving BRDF on its whole domain making it immediately suitable for further atmospheric RT modeling applications. The method is based on the integral equation relating surface reflected radiance, BRDF and solutions of two auxiliary atmosphere-only RT problems. The method requires kernel-based BRDF. The weights of the kernels are obtained with a quickly converging iterative procedure. RT modeling has to be done only one time before the start of iterative process

Iterative Discrete Ordinates Solution of the Equation for the Surface-Reflected Radiance

This paper presents a new method of numerical solution of the integral equation for the radiance reflected from an anisotropic surface. The equation relates the radiance at the surface level with BRDF and solutions of the standard radiative transfer problems for a slab with no reflection on its surfaces. It is also shown that the kernel of the equation satisfies the condition of the existence of a unique solution and the convergence of the successive approximations to that solution. The developed method features two basic steps: discretization on a 2D quadrature, and solving the resulting system of algebraic equations with successive over-relaxation method based on the Gauss-Seidel iterative process. Presented numerical examples show good coincidence between the surface-reflected radiance obtained with DISORT and the proposed method. Analysis of contributions of the direct and diffuse (but not yet reflected) parts of the downward radiance to the total solution is performed. Together, they represent a very good initial guess for the iterative process. This fact ensures fast convergence. The numerical evidence is given that the fastest convergence occurs with the relaxation parameter of 1 (no relaxation). An integral equation for BRDF is derived as inversion of the original equation. The potential of this new equation for BRDF retrievals is analyzed. The approach is found not viable as the BRDF equation appears to be an ill-posed problem, and it requires knowledge the surface-reflected radiance on the entire domain of both Sun and viewing zenith angles

Passive remote sensing of tropospheric aerosol and atmospheric correction for the aerosol effect

The launch of ADEOS in August 1996 with POLDER, TOMS, and OCTS instruments on board and the future launch of EOS-AM 1 in mid-1998 with MODIS and MISR instruments on board start a new era in remote sensing of aerosol as part of a new remote sensing of the whole Earth system (see a list of the acronyms in the Notation section of the paper). These platforms will be followed by other international platforms with unique aerosol sensing capability, some still in this century (e.g., ENVISAT in 1999). These international spaceborne multispectral, multiangular, and polarization measurements, combined for the first time with international automatic, routine monitoring of aerosol from the ground, are expected to form a quantum leap in our ability to observe the highly variable global aerosol. This new capability is contrasted with present single-channel techniques for AVHRR, Meteosat, and GOES that although poorly calibrated and poorly characterized already generated important aerosol global maps and regional transport assessments. The new data will improve significantly atmospheric corrections for the aerosol effect on remote sensing of the oceans and be used to generate first real-time atmospheric corrections over the land. This special issue summarizes the science behind this change in remote sensing, and the sensitivity studies and applications of the new algorithms to data from present satellite and aircraft instruments. Background information and a summary of a critical discussion that took place in a workshop devoted to this topic is given in this introductory paper. In the discussion it was concluded that the anticipated remote sensing of aerosol simultaneously from several space platforms with different observation strategies, together with continuous validations around the world, is expected to be of significant importance to test remote sensing approaches to characterize the complex and highly variable aerosol field. So far, we have only partial understanding of the information content and accuracy of the radiative transfer inversion of aerosol information from the satellite data, due to lack of sufficient theoretical analysis and applications to proper field data. This limitation will make the anticipated new data even more interesting and challenging. A main concern is the present inadequate ability to sense aerosol absorption, from space or from the ground. Absorption is a critical parameter for climate studies and atmospheric corrections. Over oceans, main concerns are the effects of white caps and dust on the correction scheme. Future improvement in aerosol retrieval and atmospheric corrections will require better climatology of the aerosol properties and understanding of the effects of mixed composition and shape of the particles. The main ingredient missing in the planned remote sensing of aerosol are spaceborne and ground-based lidar observations of the aerosol profiles

Statistically Optimized Inversion Algorithm for Enhanced Retrieval of Aerosol Properties from Spectral Multi-Angle Polarimetric Satellite Observations

The proposed development is an attempt to enhance aerosol retrieval by emphasizing statistical optimization in inversion of advanced satellite observations. This optimization concept improves retrieval accuracy relying on the knowledge of measurement error distribution. Efficient application of such optimization requires pronounced data redundancy (excess of the measurements number over number of unknowns) that is not common in satellite observations. The POLDER imager on board the PARASOL microsatellite registers spectral polarimetric characteristics of the reflected atmospheric radiation at up to 16 viewing directions over each observed pixel. The completeness of such observations is notably higher than for most currently operating passive satellite aerosol sensors. This provides an opportunity for profound utilization of statistical optimization principles in satellite data inversion. The proposed retrieval scheme is designed as statistically optimized multi-variable fitting of all available angular observations obtained by the POLDER sensor in the window spectral channels where absorption by gas is minimal. The total number of such observations by PARASOL always exceeds a hundred over each pixel and the statistical optimization concept promises to be efficient even if the algorithm retrieves several tens of aerosol parameters. Based on this idea, the proposed algorithm uses a large number of unknowns and is aimed at retrieval of extended set of parameters affecting measured radiation

Land Surface Temperature Measurements form EOS MODIS Data

We have developed a physics-based land-surface temperature (LST) algorithm for simultaneously retrieving surface band-averaged emissivities and temperatures from day/night pairs of MODIS (Moderate Resolution Imaging Spectroradiometer) data in seven thermal infrared bands. The set of 14 nonlinear equations in the algorithm is solved with the statistical regression method and the least-squares fit method. This new LST algorithm was tested with simulated MODIS data for 80 sets of band-averaged emissivities calculated from published spectral data of terrestrial materials in wide ranges of atmospheric and surface temperature conditions. Comprehensive sensitivity and error analysis has been made to evaluate the performance of the new LST algorithm and its dependence on variations in surface emissivity and temperature, upon atmospheric conditions, as well as the noise-equivalent temperature difference (NE(Delta)T) and calibration accuracy specifications of the MODIS instrument. In cases with a systematic calibration error of 0.5%, the standard deviations of errors in retrieved surface daytime and nighttime temperatures fall between 0.4-0.5 K over a wide range of surface temperatures for mid-latitude summer conditions. The standard deviations of errors in retrieved emissivities in bands 31 and 32 (in the 10-12.5 micrometer IR spectral window region) are 0.009, and the maximum error in retrieved LST values falls between 2-3 K. Several issues related to the day/night LST algorithm (uncertainties in the day/night registration and in surface emissivity changes caused by dew occurrence, and the cloud cover) have been investigated. The LST algorithms have been validated with MODIS Airborne Simulator (MAS) dada and ground-based measurement data in two field campaigns conducted in Railroad Valley playa, NV in 1995 and 1996. The MODIS LST version 1 software has been delivered

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

Optical Remote Sensing Of Snow On Sea Ice: Ground Measurements, Satellite Data Analysis, And Radiative Transfer Modeling

Thesis (Ph.D.) University of Alaska Fairbanks, 2002The successful launch of the Terra satellite on December 18, 1999 opened a new era of earth observation from space. This thesis is motivated by the need for validation and promotion of the use of snow and sea ice products derived from MODIS, one of the main sensors aboard the Terra and Aqua satellites. Three cruises were made in the Southern Ocean, in the Ross, Amundsen and Bellingshausen seas. Measurements of all-wave albedo, spectral albedo, BRDF, snow surface temperature, snow grain size, and snow stratification etc. were carried out on pack ice floes and landfast ice. In situ measurements were also carried out concurrently with MODIS. The effect of snow physical parameters on the radiative quantities such as all-wave albedo, spectral albedo and bidirectional reflectance are studied using statistical techniques and radiative transfer modeling, including single scattering and multiple scattering. The whole thesis consists of six major parts. The first part (chapter 1) is a review of the present research work on the optical remote sensing of snow. The second part (chapter 2) describes the instrumentation and data-collection of ground measurements of all-wave albedo, spectral albedo and bidirectional reflectance distribution function (BRDF) of snow and sea ice in the visible-near-infrared (VNIR) domain in Western Antarctica. The third part (chapter 3) contains a detailed multivariate correlation and regression analysis of the measured radiative quantities with snow physical parameters such as snow density, surface temperature, single and composite grain size and number density. The fourth part (chapter 4) describes the validation of MODIS satellite data acquired concurrently with the ground measurements. The radiances collected by the MODIS sensor are converted to ground snow surface reflectances by removing the atmospheric effect using a radiative transfer algorithm (6S). Ground measured reflectance is corrected for ice concentration at the subpixel level so that the in situ and space-borne measured reflectance data are comparable. The fifth part (chapter 5) investigates the single scattering properties (extinction optical depth, single albedo, and the phase function or asymmetry factor) of snow grains (single or composite), which were calculated using the geometrical optical method. A computer code, GOMsnow, is developed and is tested against benchmark results obtained from an exact Mie scattering code (MIE0) and a Monte Carlo code. The sixth part (chapter 6) describes radiative transfer modeling of spectral albedo using a multi-layer snow model with a multiple scattering algorithm (DISORT). The effect of snow stratification on the spectral albedo is explored. The vertical heterogeneity of the snow grain-size and snow mass density is investigated. It is found that optical remote sensing of snow physical parameters from satellite measurements should take the vertical variation of snow physical parameters into account. The albedo of near-infrared bands is more sensitive to the grain-size at the very top snow layer (<5cm), while the albedo of the visible bands is sensitive to the grain-size of a much thicker snow layer. Snow parameters (grain-size, for instance) retrieved with near-infrared channels only represent the very top snow layer (most probably 1--3 cm). Multi-band measurements from visible to near-infrared have the potential to retrieve the vertical profile of snow parameters up to a snow depth limited by the maximum penetration depth of blue light

Cloud Detection And Trace Gas Retrieval From The Next Generation Satellite Remote Sensing Instruments

Thesis (Ph.D.) University of Alaska Fairbanks, 2005The objective of this thesis is to develop a cloud detection algorithm suitable for the National Polar Orbiting Environmental Satellite System (NPOESS) Visible Infrared Imaging Radiometer Suite (VIIRS) and methods for atmospheric trace gas retrieval for future satellite remote sensing instruments. The development of this VIIRS cloud mask required a flowdown process of different sensor models in which a variety of sensor effects were simulated and evaluated. This included cloud simulations and cloud test development to investigate possible sensor effects, and a comprehensive flowdown analysis of the algorithm was conducted. In addition, a technique for total column water vapor retrieval using shadows was developed with the goal of enhancing water vapor retrievals under hazy atmospheric conditions. This is a new technique that relies on radiance differences between clear and shadowed surfaces, combined with ratios between water vapor absorbing and window regions. A novel method for retrieving methane amounts over water bodies, including lakes, rivers, and oceans, under conditions of sun glint has also been developed. The theoretical basis for the water vapor as well as the methane retrieval techniques is derived and simulated using a radiative transfer model