Search NASASearch

Engineering topics

Robert Spurr

Publications and source records attributed to Robert Spurr.

Revised and extended benchmark results for Rayleigh scattering of sunlight in spherical atmospheres

While most of traditional Earth-atmosphere satellite remote sensing relies on radiative transfer in the plane parallel geometry, effects of sphericity are important at high sun and view zenith angles. Broad understanding of these effects is limited and, contrary to the plane-parallel case, finding accurate numerical results to test spherical RT codes is not easy. This paper aims to partially fill in this gap. Using the full-spherical RT code MYSTIC (Monte Carlo), and the plane-parallel RT code VLIDORT (discrete ordinates) corrected for atmospheric sphericity in the single and multiple scattering, we reproduced with better accuracy and extended the benchmark results by Adams & Kattawar [1978].

spherical atmospheres

GMAO Simulated Data

Updates on GMAO simulated observations for PACE OCI.

Patricia Castellanos

TROPOMI Geometry-dependent Lambertian-Equivalent surface Reflectivity (GLER) Product for Improved Trace-Gas Retrieval

Accurate information about the reflectivity of the Earth's surface is required for most satellite retrievals of atmospheric composition, and this information is generally taken from monthly surface reflectivity climatology that neglects angular dependence. Previously we introduced Geometry-dependent Lambertian-equivalent surface reflectivity (GLER) which captures solar and satellite viewing angle dependence as well as daily and seasonal changes. GLER is calculated from simulations of Rayleigh-only top-of-atmosphere (TOA) radiances over non-Lambertian surfaces. We use NASA's Moderate Resolution Imaging Spectroradiometer (MODIS) bidirectional reflectance distribution function (BRDF) product over land and the wind-dependent Cox–Munk wave-facet slope distribution including water-leaving radiance over water to accounts for surface BRDF. We have developed global GLER product, previously for the Ozone Monitoring Instrument (OMI) and recently for Sentinel-5 Precursor (S5P) TROPOspheric Monitoring Instrument (TROPOMI) with several new improvements and updates. We have implemented the near real time daily V006 MODIS MCD43C1 BRDF data and gap-filled with a daily BRDF coefficient climatology created from 2002-2017 V006 MCD43GF data. The NASA’s Global Modelling Initiative hourly 0.25 x 0.25 deg Replay simulations are used for more accurate determination of pixel specific terrain pressure. To improve detection of seasonal snow/ice scenes, we use the 4-km snow cover product from the Interactive Multi-sensor Snow and Ice Mapping System (IMS). Finally, we use an improved version of the vector linearized discrete ordinate radiative transfer (VLIDORT) for update of the top-of-atmosphere (TOA) radiance look-up-tables (LUTs). We demonstrate how the use of GLER is beneficial to TROPOMI’s high spatial resolution (up to 3.5 km x 3.5 km) measurements to monitor atmospheric trace gas pollutants down to the sub-city scale.

TROPOMI

Development of An Improved BRDF Hotspot Model and its Use in VLIDORT to Study the Impact of Atmospheric Scattering on Hotspot Directional Signatures in the Atmosphere

The term “hotspot” refers to the sharp increase of reflectance occurring when incident (solar) and reflected (viewing) directions almost coincide in the backscatter direction. The accurate simulation of hotspot directional signatures is important for many remote sensing applications. The RossThick-LiSparse-Reciprocal (RTLSR) Bidirectional Reflectance Distribution Function (BRDF) model is widely used in radiative transfer simulations, and the hotspot model mostly used is from Maignan- Bréon but it typically requires large values of numerical quadrature and Fourier expansion terms in order to represent the hotspot accurately. To improve its use in atmospheric radiative transfer (RT) model simulations, in this paper we have developed a modified version based on the Maignan-Bréon’s hotspot BRDF model that converge much faster numerically, making it more practical for use in RT models that require Fourier expansion of BRDF to simulate the top-of-atmosphere (TOA) hotspot signatures. Using the vector linearized discrete ordinate radiative transfer model (VLIDORT), we found that reasonable TOA hotspot accuracy can be obtained with just 23 Fourier terms for clear atmospheres, and 63 Fourier terms for atmospheres with aerosol scattering. One advantage of this modified model is that the new hotspot model agrees very well with the original RossThick model away the hotspot region, making it is very convenient to use in the condition with and without hotspot in applications. This model can calculate the amplitude of hot spot accurately, and has been added in the most recent version of VLIDORT. However, there are some difference of this modified model with the original model for scattering angle close the hot spot, and it may not be appropriate for those who need an exact representation of the hot spot angular signature close to hot spot.

Xiaozhen (Shawn) Xiong

A New BRDF Hotspot Model and its Use in Improving Radiative Transfer Model Efficiency and Application to VLIDORT-based PCRTM Model

The term “hotspot” refers to the sharp increase of reflectance occurring when incident (solar) and reflected (viewing) directions almost coincide in the backscatter direction. The accurate simulation of hotspot directional signatures is important for many remote sensing applications. The RossThick-LiSparse-Reciprocal (RTLSR) Bidirectional Reflectance Distribution Function (BRDF) model is widely used in radiative transfer simulations, and the hotspot model mostly used is from Maignan- Bréon but it typically requires large values of numerical quadrature and Fourier expansion terms in order to represent the hotspot accurately. To improve its use in atmospheric radiative transfer (RT) model simulations, we have developed a modified version based on the Maignan-Bréon’s hotspot BRDF model that converge much faster numerically, making it more practical for use in RT models that require Fourier expansion of BRDF to simulate the top-of-atmosphere (TOA) hotspot signatures. On basis of the vector linearized discrete ordinate radiative transfer model (VLIDORT), we have built a line-by-line based simulation system that can simulate the TOA radiance for atmospheres with 26 gases, different aerosols, dust and clouds, and surface with different BRDF. Using simulated results under diverse atmospheric conditions, we have developed VLIDORT-PCRTM model.

Xiaozhen Xiong

Data Simulation Using VLIDORT and PCRTM and its Use for Radiometric Calibration of Hyperspectral UV-NIR-SWIR Sensors

Using the vector linearized discrete ordinate radiative transfer model (VLIDORT), a multiple-scattering multi-layer discrete ordinate scattering codes, coupled with a fast principal component (PC) radiative transfer model (PCRTM), we can simulate the radiation fields for hyperspectral UV-NIR-SWIR sensors from 250nm to 2500nm with a very high speed. Since the accurate gas absorptions, molecules and aerosol/cloud particle scatterings, and the surface reflectance have been considered in these simulations, these model data can be used for satellite or airborne sensor design, characterization, or end-to-end system level uncertainty assessment. Comparison of the modelled data with measured results has been a general practice for operational satellite sensor calibrations, such as the Ozone Mapping and Profiler Suite (OMPS) onboard on SNPP and NOAA-20. This presentation will give more details of this simulation system and models. Some examples of the generation of simulations data for CLARREO PathFinder (CPF) and the Surface Biology and Geology (SBG), as well as some examples to use these data for CPF and OMPS radiometric calibration will be presented.

Xiaozhen Xiong

Developing Fast and Accurate Radiative Transfer Models to Meet the Needs of Modern Satellite Remote Sensing Applications

Modern hyperspectral satellite remote sensors provide highly accurate measurements the Earth’s Top-of-Atmosphere (TOA) radiance, reflectance, or polarized spectra with hundreds to thousands of spectral channels and with millions of observations per day. The large data volume and high spectral dimensionality of the data pose challenges for retrieval algorithms. To process the satellite Level-1 data (e.g. calibrated TOA spectra) into Level-2 products (e.g. atmospheric and surface properties) using physical-based retrieval algorithms, accurate and fast Radiative Transfer Models (RTMs) are needed. RTMs are usually the limiting factor in determining the speed of a level-2 algorithm. For example, more than one million Line-by-Line (LBL) radiative transfer (RT) calculations are needed in order to properly capture the spectral contributions of important atmospheric molecules for an IR hyperspectral sensor with a spectral coverage from 3.5 m to 15 m or a solar hyperspectral sensor with spectral coverage from 0.25 m to 2.5 m. In this presentation, we will discuss advantages and disadvantages of different ways (e.g. correlated k and effective transmittance) to accelerate the speed of a fast RTM. We finally describe a Principal Component-based Radiative Transfer Model (PCRTM), which can calculate TOA radiance or reflectance spectra from 50 cm-1 to 40,000 cm-1 (200 m to 0.25 m). It has demonstrated very good accuracy relative to reference LBL RTMs and saves orders of magnitude in computational time. The PCRTM has been used in many satellite remote sensing applications. Examples include forward modeling in Level-2 and Level-3 retrieval algorithms, high fidelity satellite instrument simulators and instrument performance trade studies, spectral and radiometric accuracy characterizations of satellite Level-1 data, tools for inter-satellite calibrations, tools for satellite RTM lookup table generations, and tools for generating physically based training datasets for Artificial Intelligence (AI) algorithms.

climate data record

Using Polarized PCRTM-Solar Model to Validate EMIT Surface Reflectance L2A Product

Using the vector linearized discrete ordinate radiative transfer model (VLIDORT), a multiple-scattering multi-layer scattering code, coupled with a fast principal component (PC) radiative transfer model (PCRTM), we can simulate hyperspectral radiation fields from visible to shortwave infrared (VSWIR), i.e., 250 nm to 2500 nm. Since only limited accurate line-by-line calculation using VLIDORT is needed, this coupled PCRTM-VLIDORT approach is very efficient for simulating hyperspectral VSWIR spectra. These simulations incorporate accurately defined optical properties (gas absorptions, molecular and aerosol/cloud particle scattering) and detailed surface reflectance inputs. In this study we used the current L2A products of the Earth Surface Mineral Dust Source Investigation (EMIT) imaging spectrometer as the inputs to PCRTM-VLIDORT to simulate EMIT spectra. Comparing the simulated spectra with EMIT measurements provides an efficient method to validate the current EMIT product.

Xiaozhen Xiong