Engineering PapersSearch

Engineering topics

Alexei Lyapustin

Publications and source records attributed to Alexei Lyapustin.

At least 37 records · Page 2

Surface Polarized Reflectance Analysis For Aerosol Remote Sensing

We study the Earth surface polarized reflectance using data collected by a space-based lidar. Accurate modelling of the surface reflectance supports retrieval algorithm development for the current and future Earth Science missions. Strong polarization of the laser light from Cloud-Aerosol Transport System (CATS) instrument, operated in 2015-2017, and nighttime measurements yield higher signal-to-noise ratio for polarization compared to the previous analysis of reflected, initially unpolarized, solar light.

Earth surface polarized reflectance

MAIAC Thermal Technique for Smoke Injection Height from MODIS

We present a new algorithm to derive smoke plume height (Ha) using thermal contrast from the rising mixture of aerosol and emitted gases in the MODIS 11m channel. Validation shows good agreement with wind-corrected MISR-MINX values, with about 60% of the MODIS Terra thermal retrievals within 500m of MISR Ha, and 450m lower on average. The bias is expected because the thermal technique represents an effective rather than a top plume height from MISR MINX. Comparison of MODIS Aqua retrievals with CALIOP CALIPSO shows similar statistics, with standard deviation of 458m for the mean plume height and 216m lower on average. Ha is part of the MAIAC MODIS Collection 6 suite of products (MCD19), accessible via the Land Product Distributed Active Archive Center (LP DAAC). Aerosol injection height is reported in the daily MAIAC atmospheric product MCD19A2 along with the cloud mask, column water vapor, aerosol optical depth, AOD uncertainty and aerosol type, at 1km resolution on global Sinusoidal grid. Despite some limitations, the vastly increased coverage from MODIS observations makes it a valuable dataset complementing the established MISR and CALIOP products.

atmosphere

Analysis of A Severe Dust Storm and Its Impact on Air Quality Conditions Using Wrf-Chem Modeling, Satellite Imagery, and Ground Observations

This study presents a comprehensive analysis of an extreme dust event recorded in the Arabian Peninsula and the United Arab Emirates (UAE) between 31 March and 3 April 2015. Simulations of the dust event with the Weather Research and Forecasting model coupled with the Chemistry module (WRF-Chem) were analyzed and verified using MSG-SEVIRI imagery and aerosol optical depth (AOD) from the recent 1-km Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm for MODIS Terra/Aqua. Data from the National Centers for Atmospheric Prediction/National Center for Atmospheric Research (NCEP/NCAR) and the upper-air radiosonde observations were used to understand the synoptic of the event. In addition, the impact of the event on atmospheric and air quality conditions is investigated. The Air Quality Index (AQI) was calculated prior, during, and after the event to assess the degradation of air quality conditions. Simulated temperature, relative humidity, wind speed, and surface radiation were compared to observations at six monitoring stations in the UAE giving R2 values of 0.84, 0.63, 0.60, and 0.84, respectively. From 1 to 2 April 2015, both observations and simulations showed an average drop in temperature from 33 to 26 °C and radiance reduction from about 950 to 520 Wm−2. The AOD modeled by WRF-Chem showed a good correlation with Aerosol Robotic Network (AERONET) measurements in the UAE with R2 of 0.83. The AQI over the UAE reached hazardous levels during the peak of the dust event before rapidly decreasing to moderate–good air quality levels. This work is the first attempt to demonstrate the potential of using WRF-Chem to estimate AQI over the UAE along with two satellite products (MODIS-MAIAC and MSG-SEVIRI) for dust detection and tracking.

dust storm

Matrix Exponential in C/C++ Version of Vector Radiative Transfer Code IPOL

We use only left eigenvectors to evaluate the matrix exponential in the method of discrete ordinates for the vector radiative transfer equation, which neglects circular polarization, in a plane-parallel atmosphere. This is contrary to a common practice of using the right eigenvectors to evaluate the matrix exponential combined with the left eigenvectors to avoid the inversion of the matrix of the right ones. Two numerical tests for Rayleigh and Aerosol scattering confirm our idea. For better explanation of our approach and for independent crosscheck of our results, we distribute an example in C/C++.

Sergey Korkin

An Example of Aerosol Pattern Variability Over Bright Surface Using High Resolution MODIS MAIAC: The Eastern and Western Areas of the Dead Sea and Environs

The extreme rate of evaporation of the Dead Sea (DS) has serious implications for the surrounding area, including atmospheric conditions. This study analyzes the aerosol properties over the western and eastern parts of the DS during the year 2013, using MAIAC (Multi-Angle Implementation of Atmospheric Correction) for MODIS, which retrieves aerosol optical depth (AOD) data at a resolution of 1km. The main goal of the study is to evaluate MAIAC over the study area and determine, for the first time, theprevailing aerosol spatial patterns. First, the MAIAC-derived AOD data was compared with data from three nearby AERONET sites (Nes Ziona - an urban site, and Sede Boker and Masada - two arid sites), and with the conventional Dark Target (DT) and Deep Blue (DB) retrievals for the same days and locations, on a monthly basis throughout 2013. For the urban site, the correlation coefficient (r) for DTDB products showed better performance than MAIAC (r = 0.80, 0.75, and 0.64 respectively) year-round. However, in the arid zones, MAIAC showed better correspondence to AERONET sites than the conventional retrievals (r = 0.58 - 0.60 and 0.48 - 0.50 respectively). We investigated the difference in AOD levels, and its variability, between the Dead Sea coasts on a seasonal basis and calculated monthlyseasonal AOD averages for presenting AOD patterns over arid zones. Thus, we demonstrated that aerosol concentrations show a strong preference for the western coast, particularly during the summer season. This preference, is most likely a result of local anthropogenic emissions combined with the typical seasonal synoptic conditions, the Mediterranean Sea breeze, and the region complex topography. Our results also indicate that a large industrial zone showed higher AOD levels compared to an adjacent reference-site, i.e., 13 percent during the winter season.

AOD

Accuracy of RT Code SORD for Realistic Atmospheric Profiles

We discuss accuracy of our recently developed RT code SORD using 2 benchmark scenarios published by the IPRT group in 2015. These scenarios define atmospheres with a complicate dependence of scattering and absorption properties over height (profile). Equal step, dh=1km, is assumed in the profiles. We developed subroutines that split such atmospheres into layers of the same optical thickness, d tau. We provide full text of the subroutines with comments in Appendix. The d tau is a step for vertical integration in the method of successive orders. Modification of the input profiles from "equal step over h" to "equal step over tau" changes input for RT simulations. This may cause errors at or above the acceptable level of the measurement uncertainty. We show errors of the RT code SORD for both intensity and polarization. In addition to that, using our discrete ordinates RT code IPOL, we discuss one more IPRT scenario, in which changes in height profile indeed cause unacceptable errors. Clear understanding of source and magnitude of these errors is important, e.g. for the AERONET retrieval algorithm. SORD is available from ftp://maiac.gsfc.nasa.gov/pub/skorkin/SORD_IP_16B/ or by email request from the first author.

polarization

Performance of the Dot Product Function in Radiative Transfer Code SORD

The successive orders of scattering radiative transfer (RT) codes frequently call the scalar (dot) product function. In this paper, we study performance of some implementations of the dot product in the RT code SORD using 50 scenarios for light scattering in the atmosphere-surface system. In the dot product function, we use the unrolled loops technique with different unrolling factor. We also considered the intrinsic Fortran functions. We show results for two machines: ifort compiler under Windows, and pgf90 under Linux. Intrinsic DOT_PRODUCT function showed best performance for the ifort. For the pgf90, the dot product implemented with unrolling factor 4 was the fastest. The RT code SORD together with the interface that runs all the mentioned tests are publicly available from ftp:maiac.gsfc.nasa.govpubskorkinSORD_IP_16B (current release) or by email request from the corresponding (first) author.

polarized radiative transfer

Scaling Estimates of Vegetation Structure in Amazonian Tropical Forests Using Multi-Angle MODIS Observations

Detailed knowledge of vegetation structure is required for accurate modelling of terrestrial ecosystems, but direct measurements of the three dimensional distribution of canopy elements, for instance from LiDAR, are not widely available. We investigate the potential for modelling vegetation roughness, a key parameter for climatological models, from directional scattering of visible and near-infrared (NIR) reflectance acquired from NASA's Moderate Resolution Imaging Spectroradiometer (MODIS). We compare our estimates across different tropical forest types to independent measures obtained from: (1) airborne laser scanning (ALS), (2) spaceborne Geoscience Laser Altimeter System (GLAS)/ICESat, and (3) the spaceborne SeaWinds/QSCAT. Our results showed linear correlation between MODIS-derived anisotropy to ALS-derived entropy (r(exp 2)= 0.54, RMSE= 0.11), even in high biomass regions. Significant relationships were also obtained between MODIS-derived anisotropy and GLAS-derived entropy(0.52 less than or equal to r(exp 2) less than or equal to 0.61; p less than 0.05), with similar slopes and offsets found throughout the season, and RMSE between 0.26 and 0.30 (units of entropy). The relationships between the MODIS-derived anisotropy and backscattering measurements (sigma(sup 0)) from SeaWinds/QuikSCAT presented an r(exp 2) of 0.59 and a RMSE of 0.11. We conclude that multi-angular MODIS observations are suitable to extrapolate measures of canopy entropy across different forest types, providing additional estimates of vegetation structure in the Amazon.

MODIS

Spatiotemporal Prediction of Fine Particulate Matter Using High-Resolution Satellite Images in the Southeastern US 2003-2011

Numerous studies have demonstrated that fine particulate matter (PM(sub 2.5), particles smaller than 2.5 micrometers in aerodynamic diameter) is associated with adverse health outcomes. The use of ground monitoring stations of PM(sub 2.5) to assess personal exposure, however, induces measurement error. Land-use regression provides spatially resolved predictions but land-use terms do not vary temporally. Meanwhile, the advent of satellite-retrieved aerosol optical depth (AOD) products have made possible to predict the spatial and temporal patterns of PM(sub 2.5) exposures. In this paper, we used AOD data with other PM(sub 2.5) variables, such as meteorological variables, land-use regression, and spatial smoothing to predict daily concentrations of PM(sub 2.5) at a 1 sq km resolution of the Southeastern United States including the seven states of Georgia, North Carolina, South Carolina, Alabama, Tennessee, Mississippi, and Florida for the years from 2003 to 2011. We divided the study area into three regions and applied separate mixed-effect models to calibrate AOD using ground PM(sub 2.5) measurements and other spatiotemporal predictors. Using 10-fold cross-validation, we obtained out of sample R2 values of 0.77, 0.81, and 0.70 with the square root of the mean squared prediction errors of 2.89, 2.51, and 2.82 cu micrograms for regions 1, 2, and 3, respectively. The slopes of the relationships between predicted PM2.5 and held out measurements were approximately 1 indicating no bias between the observed and modeled PM(sub 2.5) concentrations. Predictions can be used in epidemiological studies investigating the effects of both acute and chronic exposures to PM(sub 2.5). Our model results will also extend the existing studies on PM(sub 2.5) which have mostly focused on urban areas because of the paucity of monitors in rural areas.

PM2.5

IPRT Polarized Radiative Transfer Model Intercomparison Project-Phase A

The polarization state of electromagnetic radiation scattered by atmospheric particles such as aerosols, cloud droplets, or ice crystals contains much more information about the optical and microphysical properties than the total intensity alone. For this reason an increasing number of polarimetric observations are performed from space, from the ground and from aircraft. Polarized radiative transfer models are required to interpret and analyse these measurements and to develop retrieval algorithms exploiting polarimetric observations. In the last years a large number of new codes have been developed, mostly for specific applications. Benchmark results are available for specific cases, but not for more sophisticated scenarios including polarized surface reflection and multi-layer atmospheres. The International Polarized Radiative Transfer (IPRT) working group of the International Radiation Commission (IRC) has initiated a model intercomparison project in order to fill this gap. This paper presents the results of the first phase A of the IPRT project which includes ten test cases, from simple setups with only one layer and Rayleigh scattering to rather sophisticated setups with a cloud embedded in a standard atmosphere above an ocean surface. All scenarios in the first phase A of the intercomparison project are for a one-dimensional plane-parallel model geometry. The commonly established benchmark results are available at the IPRT website

radiative transfer

GeoNEX: A geostationary earth observatory at NASA Earth eXchange: Earth monitoring from operational geostationary satellite systems

The latest generation of geostationary satellites (Himawari 8/9, GOES-16/17, FY-4, GK-2A) carries sensors that closely mimic the spatial and spectral characteristics of widely used polar-orbiting, global monitoring sensors such as MODIS and VIIRS. When combined, data from various currently operating/planned geostationary platforms provide a geo-ring of hyper-temporal (5-10 minutes), multispectral observations at spatial resolutions as high as 500 m. These high frequency observations offer exciting new possibilities for monitoring our planet, including better retrievals of geophysical variables by overcoming cloud cover, enabling studies of diurnally varying phenomena in the atmosphere, land, and the oceans, and support operational decision-making in agriculture, hydrology and disaster management. The NASA Earth Exchange (NEX) team, in collaboration with scientists from JAXA, KARI, NOAA and other international institutions, created the GeoNEX (www.nasa.gov/geonex) pipeline to integrate data from all available geostationary platforms and produce and distribute spatially, temporally, and radiometrically consistent data for the earth science community. We envision various institutions adapting the Geo component (e.g., GeoNOAA, GeoKARI, GeoChiba, GeoJAXA, GeoCMA) and customizing the pipeline and downstream products to serve the local/regional research and applied science communities.

Ramakrishna R Nemani

Update on EPIC MAIAC V2 Algorithm: Atmospheric Correction

Alexei Lyapustin presented a new v2 MAIAC atmospheric correction algorithm for EPIC. The v2 algorithm introduced new projection removing spatial distortions, a more accurate retrieval of spectral surface reflectance using ancillary MAIAC MODIS BRDF dataset, and a new technique for simultaneous retrieval of aerosol optical depth and spectral absorption. AERONET validation shows that v2 significantly improved AOD accuracy compared to v1, while the retrieved single scattering albedo in the Blue shows a very good accuracy comparable to that of AERONET for both biomass burning smoke and mineral dust aerosols.

Alexei Lyapustin

Using the Diurnal Variability in GeoNEX TOA Reflectances for Earth Monitoring

Observations from the third-generation geostationary satellite instruments (GOES 16/17 ABI, Himawari 8/9 AHI, and etc.) have spatial resolution and spectral band configurations comparable to flagship LEO sensors (e.g., MODIS/VIIRS). More importantly, these data are acquired at very high temporal resolution, faithfully recording the variations of the full disk of Earth at every 5-10 minutes. They thus provide unique information about Earth’s atmosphere and surface. In order to explore the unique information content of geostationary data, this study systematically analyzes the diurnal variability in the GeoNEX L1G TOA reflectance products and compares them to simulated results by state-of-the-art radiative transfer codes. Our results show that • The smoothness of the TOA reflectance diurnal cycle provides a convenient and reliable way to identify stable atmospheric conditions and filter out passing clouds/shadows. • The diurnal variability of the blue band (0.47µm) reflectance is regulated mainly by atmospheric optical conditions over a majority of land cover types. As such, the diurnal variability of the blue band data allows us to retrieve AOD without invoking the use of spectral band ratios (SRC) as in previous algorithms. • In comparison, the diurnal variability of the short-wave infrared band (2.2µm) BRFs is mainly regulated by surface reflectance and the sun-target-satellite geometry. This information allows us to test and, if suitable, retrieve surface BRDF parameters. • Spectral band ratios, especially those between the 2.2µm and 0.47µm bands, are not constant but vary by locations and sun-target-satellite geometries. Our analysis clearly demonstrates that the information provided in high-frequent geostationary observations is unique and complementary to LEO sensors. Therefore, a synergy of GEO and LEO (and other) sensors has the great potential to improve existing remote sensing models and algorithms for better Earth monitoring.

Diurnal Variability

Development of the GeoNEX Level 2G Products: Exploiting the Diurnal Variability of TOA Reflectance in Atmospheric Correction

This study develops a new atmospheric correction algorithm to generate the Level 2G products, in particular the gap-filled Surface Reflectance at 10-minute time steps, for the Geostationary-NASA Earth Exchange (GeoNEX) project. The algorithm is based on the MODIS MAIAC (Multi-Angle Implementation of Atmospheric Correction) framework but with significant modifications to exploit angular/temporal information from the diurnal variability of the GeoNEX L1G TOA (Top-of-Atmosphere) reflectance. The algorithm starts by evaluating the roughness/smoothness of the diurnal time series of the TOA reflectance. Because rapid changes in TOA reflectance are generally caused by passing clouds or shadows, rough segments of the time series are automatically filtered out while the smooth segments are further tested for brightness and temperature to identify clear-sky and snow-free observations. Next the algorithm runs the MAIAC RTM (Radiative Transfer Model) to retrieve the Ross-Thick-Li-Sparse (RTLS) BRDF model parameters and the daily-mean atmospheric optical depth (AOD) that allow the RTM to optimally simulate the observed diurnal variability of clear-sky TOA reflectance. Once the initial RTLS parameters are retrieved after the algorithm’s burn-in period, they are used as the prior information to predict the AOD level for the next days, while the subsequent clear-sky observations are used to make necessary adjustments to the RTLS parameters in an continuous fashion. This “prediction-analysis” cycle is then iterated to process the full time series of the L1G data, skipping only total-cloudy days or when surface snow is detected. We tested the algorithm over a list of selected AERONET sites. The retrieved results (the daily mean AOD and the RTLS parameters) reasonably agree with the ground-based measurements. Importantly, the results indicate that the diurnal cycles of surface reflectance are continuous functions of the illumination-view geometry. Thus we can use the retrieved RTLS model to accurately fill in data gaps on partial cloudy days. Also, our algorithm is totally independent from the traditional approaches based on the use of spectral band ratios between the shortwave infrared (e.g., 2.2µm) and the visible (e.g., 0.47µm and 0.64µm) bands. Our results thus demonstrate that the high-frequent diurnal geostationary observations contain unique information that helps us improve atmospheric correction of remote sensing data.

GeoNEX

Joint Retrieval of Surface BRDF from Geostationary and Polar-Orbiting Satellite Sensors

The latest geostationary sensors like GOES 16/17 ABI and Himawari 8/9 AHI provide high frequent observations of the Earth surface with continuously changing solar illumination geometries, which allow us to retrieve the surface Bidirectional Reflectance Distribution Function (BRDF) with dedicated atmospheric correction algorithms like MAIAC (Multi-Angle Implementation of Atmospheric Correction). However, because the viewing geometry of a specific location from the geostationary satellites are fixed, the angular sampling of surface BRDF by GEO (Geostationary Earth Orbit) sensors is far from comprehensive. This study tries to address this issue by exploring a GEO-LEO (Low-Earth-Orbit) synergy, in particular, jointly retrieving surface BRDF parameters with concurrent ABI/AHI and VIIRS top-of-atmosphere (TOA) reflectance for the near-infrared (NIR) band. The NIR band is chosen because the ABI, AHI, and VIIRS instruments have very similar spectral response functions in this band and therefore simplifies the requirements for cross-sensor radiometric calibration. We compile ABI/AHI and VIIRS TOA data with the best quality (cloud-free and low aerosol loading) at chosen AERONET sites. We then run the GeoNEX MAIAC algorithm to retrieve the Ross-Thick-Li-Sparse (RTLS) surface BRDF parameters with or without the AEORNET measured atmospheric aerosol optical depth (AOD) as inputs. The joint retrieval results are considered the best estimate of surface BRDF. We compare the joint BRDF retrievals with the corresponding MAIAC BRDF products, retrieved with ABI/AHI or VIIRS separately, to evaluate their differences. We expect that the jointly retrieved BRDF data are more robust than the standard products and may help us reduce uncertainties in higher-level earth observation satellite products.

Remote Sensing

The Complex Refractive Indices of Mineral Aerosols and Why They Matter

Aerosol refractive indices are fundamental parameters that are generally measured by spec-troscopists with specialized knowledge. We in the Earth science community frequently utilizethese refractive indices because they are essential for computing aerosol radiative effects andretrieving aerosol composition. Unfortunately, there are a wide variety of refractive indiceswith significant differences for some aerosol species (e.g., hematite) and a lack of refractiveindex choices for other aerosols (e.g., clay minerals, goethite), and this hinders our ability toaccurately compute the radiative effect of mineral dust. Additionally, the mineral refractiveindices used in atmospheric science are not necessarily linked to the mineral reflectancesused to identify surface mineralogy; this creates a disconnect between the atmosphere andthe surface that frustrates closure analyses.In this talk, we will present an overview of some refractive indices of radiative importancein aeolian dust (illite, kaolinite, montmorillonite, hematite, goethite). We will discuss howmineral refractive indices are used in aerosol retrievals, and how we can use remote sensingretrievals to narrow the range of viable choices. We will also discuss how we can use pub-lished spectroscopic measurements to extrapolate the refractive indices that are inferred ata handful of visible and near-infrared wavelengths to the longwave regime. Finally, we willdiscuss how working groups like MIRA (Models, In situ, and Remote sensing of Aerosols;https://science.larc.nasa.gov/mira-wg/) and community repositories like TAO (Tables ofAerosol Optics) can improve radiative closure by enhancing interactions between the threedisciplines.1

aerosols