Engineering topics

Alexei Lyapustin

Publications and source records attributed to Alexei Lyapustin.

At least 19 records

NASA NTRS2023

Calibration of the SNPP and NOAA 20 VIIRS Sensors for Continuity of the MODIS Climate Data Records

Accurate long-term sensor calibration and periodic re-processing to ensure consistency and continuity of atmospheric, land and ocean geophysical retrievals from space within the mission period and across different missions is a major requirement of climate data records. In this work, we applied the Multi-Angle Implementation of Atmospheric Correction (MAIAC)-based vicarious calibration technique over Libya-4 desert site to perform calibration analysis of Visible Infrared Imaging Radiometer Suite (VIIRS) on Suomi National Polar-orbiting Partnership (SNPP) and NOAA-20 satellites. For both VIIRS sensors we characterized residual linear calibration trends and cross-calibrated both sensors to MODerate resolution Imaging Spectroradiometer (MODIS) Aqua regarded as a calibration standard. The relative spectral response (RSR) differences were accounted for using the German Aerospace Center (DLR) Earth Sensing Imaging Spectrometer (DESIS) hyperspectral surface reflectance data. Our results agree with independent vicarious calibration results of both the MODIS/VIIRS Characterization Support Team as well as the CERES Imager and Geostationary Calibration Group within estimated uncertainty of 1–2%. Analysis of MAIAC geophysical products with the new calibration shows a high level of agreement of MAIAC aerosol, surface reflectance and NDVI records between MODIS and VIIRS. Excluding high aerosol optical depth (AOD), all three sensors agree in AOD with mean difference (MD) less than 0.01 and residual mean squared difference rmsd ∼ 0.04. Spectral geometrically normalized surface reflectance agrees within rmsd of 0.003–0.005 in the visible and 0.01–0.012 at longer wavelengths. The residual surface reflectance differences are fully explained by differences in spectral filter functions. Finally, difference in NDVI is characterized by rmsd ∼ 0.02 and MD less than 0.003 for NDVI based on VIIRS imagery bands I1/I2 and less than 0.01 for NDVI based on VIIRS radiometric bands M5/M7. In practical sense, these numbers indicate consistency and continuity in MAIAC records ensuring the smooth transition from MODIS to VIIRS.

MAIAC↗

NASA NTRS2023

Assessment of the Impact of Discontinuity in Satellite Instruments and Retrievals on Global PM2.5

Exposure to fine particulate matter (PM2.5) is the leading environmental risk factor for mortality globally. Satellite-derived estimates of surface PM2.5 developed from a combination of satellites, simulations, and ground monitor data are relied upon for health impact studies. The ability to develop satellite-derived PM2.5 estimates requires the continued availability of aerosol optical depth (AOD) sources. This work examines the impact of the addition or loss of satellite AOD data sources on global PM2.5 estimation and the impact of continuing the longterm record with AOD from the Visible Infrared Imaging Radiometer Suite (VIIRS) on the Suomi-National Polar orbiting Partnership (S-NPP) satellite after the loss of the MODIS (MODerate resolution Imaging Spectroradiometer) and MISR (Multi-angle Imaging Spectroradiometer) instruments on board the Terra and Aqua satellites. We find that the addition of VIIRS S-NPP AOD products to geophysical PM2.5 estimates from satellites and simulations causes regional differences that correspond to differences in the VIIRS and MODIS Deep Blue AOD algorithms and sampling. Changes in long-term trends and timeseries due to the addition or loss of AOD data sources are generally within their uncertainties. Statistical fusion with ground monitor data partially corrects for changes due to sampling differences when introducing the VIIRS AOD products, but uncertainty remains over desert regions where ground monitor coverage is sparse. This work provides promise for the sustained development of global satellite-derived PM2.5 estimates, despite discontinuities in instruments and retrieval methods.

Remote sensing↗

NASA NTRS2023

Radiative interaction of atmosphere and surface: write up with elements of code

In passive satellite remote sensing of the Earth, separation of the path radiance (atmosphere-only contribution) from the surface reflection remains a “significant challenge”. Recent literature names it among the gaps in radiative transfer (RT) topics that “require continued research in the near future”. The challenge comes from multiple reflections (bouncing) between the atmosphere and surface – radiative interaction. In this paper we use a known RT technique, the matrix-operator method (MOM), and a new modification of the monochromatic vector RT (vRT) code IPOL (Intensity and POLarization) to simulate the interaction of a plane-parallel atmosphere and a few widely used surface reflection models. Following the idea of the Green’s function method, IPOL no longer takes the surface model parameters on input. Instead, it provides the path radiance, and the atmospheric reflection and transmission matrices as output. Despite many RT codes use the MOM formalism, this output does not seem common. The surface reflection matrix is computed externally. Therefore, this paper extends the Green’s function atmospheric correction technique to the case of polarized light. Aiming clarity rather than performance, we explain in Python the structure of the surface matrices for the isotropic (Lambertian), directional unpolarized, and polarized ocean reflection models. We then combine these surface matrices and the precomputed IPOL output to get numerically accurate signal at the top of atmosphere (TOA) and test it vs. published benchmarks. Then, for each benchmark scenario we show how to get the surface from the TOA signal, i.e. perform the RT-based atmospheric correction.

radiative transfer↗

NASA NTRS2022

A Novel Atmospheric Correction Algorithm to Exploit the Diurnal Variability in Hypertemporal Geostationary Observations

This study developed a new atmospheric correction algorithm, GeoNEX-AC, that is independent from the traditional use of spectral band ratios but dedicated to exploiting information from the diurnal variability in the hypertemporal geostationary observations. The algorithm starts by evaluating smooth segments of the diurnal time series of the top-of-atmosphere (TOA) reflectance to identify clear-sky and snow-free observations. It then attempts to retrieve the Ross-Thick–Li-Sparse (RTLS) surface bi-directional reflectance distribution function (BRDF) parameters and the daily mean atmospheric optical depth (AOD) with an atmospheric radiative transfer model (RTM) to optimally simulate the observed diurnal variability in the clear-sky TOA reflectance. Once the initial RTLS parameters are retrieved after the algorithm’s burn-in period, they serve as the prior information to estimate the AOD levels for the following days and update the surface BRDF information with the new clear-sky observations. This process is iterated through the full time span of the observations, skipping only totally cloudy days or when surface snow is detected. We tested the algorithm over various Aerosol Robotic Network (AERONET) sites and the retrieved results well agree with the ground-based measurements. This study demonstrates that the high-frequency diurnal geostationary observations contain unique information that can help to address the atmospheric correction problem from new directions.

atmospheric correction↗

NASA NTRS2022

Inferring iron-oxide species content in atmospheric mineral dust from DSCOVR EPIC observations

The iron-oxide content of dust in the atmosphere and most notably its apportionment between hematite (α-Fe2O3) and goethite (α-FeOOH) are key determinants in quantifying dust’s light absorption, its top of atmosphere UV radiances used for dust monitoring, and ultimately shortwave dust direct radiative effects (DRE). Hematite and goethite column mass concentrations and iron-oxide mass fractions of total dust mass concentration were retrieved from the DeepSpace Climate Observatory (DSCOVR) Earth Polychromatic Imaging Camera (EPIC) measurements in the ultraviolet–visible (UV–Vis) channels. The retrievals were performed for dust-identified aerosol plumes over land using aerosol optical depth (AOD) and spectral imaginary refractive index provided by the Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm over six continental regions (North America, North Africa, West Asia, Central Asia, EastAsia, and Australia). The dust particles are represented as an internal mixture of non-absorbing host and absorbing hematite and goethite. We use the Maxwell–Garnett effective medium approximation with carefully selected complex refractive indices of hematite and goethite that produce mass fractions of iron oxides species consistent with in situ values found in the literature to derive the hematite and goethite volumetric/mass concentrations from MAIAC EPIC products. We compared the retrieved hematite and goethite concentrations with in situ dust aerosol mineralogical content measurements, as well as with published data. Our data display variations within the published range of hematite, goethite, and iron-oxide mass fractions for pure mineral dust cases. A specific analysis is presented for 15 sites over the main dust source regions. Sites in the central Sahara, Sahel, and Middle East exhibit a greater temporal variability of iron oxides relative to other sites. Niger site(13.52°N, 2.63°E) is dominated by goethite over Harmattan season with median of ~2 weight percentage (wt.%) of iron-oxide. Saudi Arabia site (27.49°N, 41.98°E) over Middle East also exhibited surge of goethite content with the beginning of Shamal season. The Sahel dust is richer in iron-oxide than Saharan and northern China dust except in Summer. The Bodélé Depression area shows a distinctively lower iron-oxide concentration (~1 wt.%) throughout the year. Finally, we show that EPIC data allow to constrain the hematite refractive index. Specifically, we select 5 out of 13 different number of hematite refractive indices widely variable in published laboratory studies by constraining the iron-oxide mass ratio to the known measured values. Provided climatology of hematite and goethite mass fractions across main dust regions of the Earth will be useful for dust shortwave DRE studies and climate modeling.

Sujung Go↗

NASA NTRS2021

Evaluation and Intercomparison of Wildfire Smoke Forecasts from Multiple Modeling Systems for the 2019 Williams Flats Fire

Wildfire smoke is one of the most significant concerns of human and environmental health, associated with its substantial impacts on air quality, weather, and climate. However, biomass burning emissions and smoke remain among the largest sources of uncertainties in air quality forecasts. In this study, we evaluate the smoke emissions and plume forecasts from 12 state-of-the-art air quality forecasting systems during the Williams Flats fire in Washington State, US, August 2019, which was intensively observed during the Fire Influence on Regional to Global Environments and Air Quality (FIREX-AQ) field campaign. Model forecasts with lead times within 1 d are intercompared under the same framework based on observations from multiple platforms to reveal their performance regarding fire emissions, aerosol optical depth (AOD), surface PM2.5, plume injection, and surface PM2.5 to AOD ratio. The comparison of smoke organic carbon (OC) emissions suggests a large range of daily totals among the models, with a factor of 20 to 50. Limited representations of the diurnal patterns and day-to-day variations of emissions highlight the need to incorporate new methodologies to predict the temporal evolution and reduce uncertainty of smoke emission estimates. The evaluation of smoke AOD (sAOD) forecasts suggests overall underpredictions in both the magnitude and smoke plume area for nearly all models, although the high-resolution models have a better representation of the fine-scale structures of smoke plumes. The models driven by fire radiative power (FRP)-based fire emissions or assimilating satellite AOD data generally outperform the others. Additionally, limitations of the persistence assumption used when predicting smoke emissions are revealed by substantial underpredictions of sAOD on 8 August 2019, mainly over the transported smoke plumes, owing to the underestimated emissions on 7 August. In contrast, the surface smoke PM2.5 (sPM2.5) forecasts show both positive and negative overall biases for these models, with most members presenting more considerable diurnal variations of sPM2.5. Overpredictions of sPM2.5 are found for the models driven by FRP-based emissions during nighttime, suggesting the necessity to improve vertical emission allocation within and above the planetary boundary layer (PBL). Smoke injection heights are further evaluated using the NASA Langley Research Center's Differential Absorption High Spectral Resolution Lidar (DIAL-HSRL) data collected during the flight observations. As the fire became stronger over 3–8 August, the plume height became deeper, with a day-to-day range of about 2–9 km a.g.l. However, narrower ranges are found for all models, with a tendency of overpredicting the plume heights for the shallower injection transects and underpredicting for the days showing deeper injections. The misrepresented plume injection heights lead to inaccurate vertical plume allocations along the transects corresponding to transported smoke that is 1 d old. Discrepancies in model performance for surface PM2.5 and AOD are further suggested by the evaluation of their ratio, which cannot be compensated for by solely adjusting the smoke emissions but are more attributable to model representations of plume injections, besides other possible factors including the evolution of PBL depths and aerosol optical property assumptions. By consolidating multiple forecast systems, these results provide strategic insight on pathways to improve smoke forecasts.

AOD↗

NASA NTRS2021

Evaluation of Novel NASA MODIS and VIIRS Aerosol Products and Assessment of Smoke Height Boundary Layer Ratio During Extreme Smoke Events in the Western U.S.

We analyze new aerosol products from NASA satellite retrievals over the western U.S. during August 2013, with special attention to locally generated wildfire smoke and downwind plume structures. Aerosol optical depth (AOD) at 550 nm from MODIS (Terra and Aqua Collections 6 and 6.1) and VIIRS Deep Blue (DB) and MODIS (Terra and Aqua) Multi‐Angle Implementation of Atmospheric Correction (MAIAC) retrievals are evaluated against ground‐based AErosol RObotic NETwork (AERONET) observations. We find a significant improvement in correlation with AERONET and other metrics in the latest DB AOD (MODIS C6.1 r^2=0.75, VIIRS r^2=0.79) compared to MODIS C6 (r2=0.62). In general, MAIAC (r^2=0.84) and DB (MODIS C6.1 and VIIRS) present similar statistical evaluation metrics for the western U.S. and are useful tools to characterize aerosol loading associated with wildfire smoke. We also evaluate three novel NASA MODIS plume injection height (PIH) products, one from MAIAC and two from the Aerosol Single scattering albedo and layer Height Estimation (ASHE) (MODIS and VIIRS) algorithm. Both Terra and Aqua MAIAC PIHs statistically agree with ground‐based and satellite lidar observations near the fire source, as do ASHE, although the latter is sensitive to assumptions about aerosol absorption properties. We introduce a first‐order approximation Smoke Height Boundary Layer Ratio (SHBLR) to qualitatively distinguish between aerosol pollution within the planetary boundary layer and the free troposphere. We summarize the scope, limitations, and suggestions for scientific applications of surface level aerosol concentrations specific to wildfire emissions and smoke plumes using these novel NASA MODIS and VIIRS aerosol products.

Plume injection height↗

NASA NTRS2021

Retrievals of Aerosol Optical Depth and Spectral Absorption from DSCOVR EPIC

A new algorithm is described for joint retrievals of the aerosol optical depth and spectral absorption from EPIC observations in the UV—Vis spectral range. The retrievals are illustrated on examples of the wildfire smoke events over North America, and dust storms over greater Sahara region in 2018. An initial evaluation of single scattering albedo (SSA) at 443 nm over these regions shows a good agreement with AERONET data, generally within the uncertainty of AERONET SSA of ± 0.03. A particularly good agreement is achieved for dust with R∼0.62, rmse∼0.02, negligible bias, and 85% points within the expected error. This new capability is part of version 2 MAIAC EPIC algorithm. The v2 algorithm has recently completed reprocessing of the EPIC record covering the period of 2015–2020.

EPIC↗

NASA NTRS2021

The AOD Sensitivity Comparison between MODIS Multi-Angle Implementation of Atmospheric Correction (MAIAC) and Standard MODIS Surface Reflectance

This study compares the Multi-Angle Implementation of Atmospheric Correction (MAIAC) and standard MODIS surface reflectance (SR) products (MOD09)for various aerosol optical depth (AOD) levels. Data from one MODIS tile in northern China for the whole year 2018 is compared. The results reveal that MAIAC SR has good stability for the full range of AOD for which SR is reported, whereas MOD09 SR shows increasing high bias with AOD increases. The cross-comparison of spectral characteristics between MAIAC SR and MOD09 SR shows a systematic MAIAC-MOD09 difference increasing from NIR to Blue. This pattern is consistent with the bias caused by Lambertian assumption used in MOD09 SR algorithm.

MAIAC↗

NASA NTRS2020

An AeroCom–AeroSat study: intercomparison of satellite AOD datasets for aerosol model evaluation

To better understand and characterize current uncertainties in the important observational constraint of climate models of aerosol optical depth (AOD), we evaluate and intercompare 14 satellite products, representing nine different retrieval algorithm families using observations from five different sensors on six different platforms. The satellite products (super-observations consisting of 1°×1° daily aggregated retrievals drawn from the years 2006, 2008 and 2010) are evaluated with AErosol RObotic NETwork (AERONET) and Maritime Aerosol Network (MAN) data. Results show that different products exhibit different regionally varying biases (both under- and overestimates) that may reach ±50 %, although a typical bias would be 15 %–25 % (depending on the product). In addition to these biases, the products exhibit random errors that can be 1.6 to 3 times as large. Most products show similar performance, although there are a few exceptions with either larger biases or larger random errors. The intercomparison of satellite products extends this analysis and provides spatial context to it. In particular, we show that aggregated satellite AOD agrees much better than the spatial coverage (often driven by cloud masks) within the 1°×1° grid cells. Up to ∼50 % of the difference between satellite AOD is attributed to cloud contamination. The diversity in AOD products shows clear spatial patterns and varies from 10 % (parts of the ocean) to 100 % (central Asia and Australia). More importantly, we show that the diversity may be used as an indication of AOD uncertainty, at least for the better performing products. This provides modellers with a global map of expected AOD uncertainty in satellite products, allows assessment of products away from AERONET sites, can provide guidance for future AERONET locations and offers suggestions for product improvements. We account for statistical and sampling noise in our analyses. Sampling noise, variations due to the evaluation of different subsets of the data, causes important changes in error metrics. The consequences of this noise term for product evaluation are discussed.

Aerosol↗

NASA NTRS2020

Spatial validation reveals poor predictive performance of large-scale ecological mapping models

Mapping aboveground forest biomass is central for assessing the global carbon balance. However, current large-scale maps show strong disparities, despite good validation statistics of their underlying models. Here, we attribute this contradiction to a flaw in the validation methods, which ignore spatial autocorrelation (SAC) in data, leading to overoptimistic assessment of model predictive power. To illustrate this issue, we reproduce the approach of large-scale mapping studies using a massive forest inventory dataset of 11.8 million trees in central Africa to train and validate a random forest model based on multispectral and environmental variables. A standard nonspatial validation method suggests that the model predicts more than half of the forest biomass variation, while spatial validation methods accounting for SAC reveal quasi-null predictive power. This study underscores how a common practice in big data mapping studies shows an apparent high predictive power, even when predictors have poor relationships with the ecological variable of interest, thus possibly leading to erroneous maps and interpretations.

biomass↗

NASA NTRS2020

Advancing Methodologies for Applying Machine Learning and Evaluating Spatiotemporal Models of Fine Particulate Matter (PM 2.5 ) Using Satellite Data Over Large Regions

Reconstructing the distribution of fine particulate matter (PM 2.5 ) in space and time, even far from ground monitoring sites, is an important exposure science contribution to epidemiologic analyses of PM 2.5 health impacts. Flexible statistical methods for prediction have demonstrated the integration of satellite observations with other predictors, yet these algorithms are susceptible to overfitting the spatiotemporal structure of the training datasets. We present a new approach for predicting PM 2.5 using machine-learning methods and evaluating prediction models for the goal of making predictions where they were not previously available. We apply extreme gradient boosting (XGBoost) modeling to predict daily PM 2.5 on a 1 x 1 km 2 resolution for a 13 state region in the Northeastern USA for the years 2000–2015 using satellite-derived aerosol optical depth and implement a recursive feature selection to develop a parsimonious model. We demonstrate excellent predictions of withheld observations but also contrast an RMSE of 3.11 μg/m 3 in our spatial cross-validation withholding nearby sites versus an overfit RMSE of 2.10 μg/m 3 using a more conventional random ten-fold splitting of the dataset. As the field of exposure science moves forward with the use of advanced machine-learning approaches for spatiotemporal modeling of air pollutants, our results show the importance of addressing data leakage in training, overfitting to spatiotemporal structure, and the impact of the predominance of ground monitoring sites in dense urban sub-networks on model evaluation. The strengths of our resultant modeling approach for exposure in epidemiologic studies of PM 2.5 include improved efficiency, parsimony, and interpretability with robust validation while still accommodating complex spatiotemporal relationships.

air pollution↗

NASA NTRS2020

The AERONET Version 3 Aerosol Retrieval Algorithm, Associated Uncertainties and Comparisons to Version 2

The Aerosol Robotic Network (AERONET) Version 3 (V3) aerosol retrieval algorithm is described, which is based on the Version 2 (V2) algorithm with numerous updates. Comparisons of V3 aerosol retrievals to those of V2 are presented, along with a new approach to estimate uncertainties in many of the retrieved aerosol parameters. Changes in the V3 aerosol retrieval algorithm include (1) a new polarized radiative transfer code (RTC), which replaced the scalar RTC of V2, (2) detailed characterization of gas absorption by adding NO2 and H2O to specify total gas absorption in the atmospheric column, specification of vertical profiles of all the atmospheric species, (3) new bidirectional reflectance distribution function (BRDF) parameters for land sites adopted from the MODIS BRDF/Albedo product, (4) a new version of the extraterrestrial solar flux spectrum, and (5) a new temperature correction procedure of both direct Sun and sky radiance measurements. The potential effect of each change in V3 on single scattering albedo (SSA) retrievals was analyzed. The operational almucantar retrievals of V2 versus V3 were compared for four AERONET sites: GSFC, Mezaira, Mongu, and Kanpur. Analysis showed very good agreement in retrieved parameters of the size distributions. Comparisons of SSA retrievals for dust aerosols (Mezaira) showed a good agreement in 440 nm SSA, while for longer wavelengths V3 SSAs are systematically higher than those of V2, with the largest mean difference at 675 nm due to cumulative effects of both extraterrestrial solar flux and BRDF changes. For non-dust aerosols, the largest SSA deviation is at 675 nm due to differences in extraterrestrial solar flux spectrums used in each version. Further, the SSA 675 nm mean differences are very different for weakly (GSFC) and strongly (Mongu) absorbing aerosols, which is explained by the lower sensitivity to a bias in aerosol scattering optical depth by less absorbing aerosols. A new hybrid (HYB) sky radiance measurement scan is introduced and discussed. The HYB combines features of scans in two different planes to maximize the range of scattering angles and achieve scan symmetry, thereby allowing for cloud screening and spatial averaging, which is an advantage over the principal plane scan that lacks robust symmetry. We show that due to an extended range of scattering angles, HYB SSA retrievals for dust aerosols exhibit smaller variability with solar zenith angles (SZAs) than those of almucantar (ALM), which allows extension of HYB SSA retrievals to SZAs less than 50∘ to as small as 25∘. The comparison of SSA retrievals from closely time-matched HYB and ALM scans in the 50 to 75∘ SZA range showed good agreement with the differences below ∼0.005. We also present an approach to estimate retrieval uncertainties which utilizes the variability in retrieved parameters generated by perturbing both measurements and auxiliary input parameters as a proxy for retrieval uncertainty. The perturbations in measurements and auxiliary inputs are assumed as estimated biases in aerosol optical depth (AOD), radiometric calibration of sky radiances combined with solar spectral irradiance, and surface reflectance. For each set of Level 2 Sun/sky radiometer observations, 27 inputs corresponding to 27 combinations of biases were produced and separately inverted to generate the following statistics of the inversion results: average, standard deviation, minimum and maximum values. From these statistics, standard deviation (labeled U27) is used as a proxy for estimated uncertainty, and a lookup table (LUT) approach was implemented to reduce the computational time. The U27 climatological LUT was generated from the entire AERONET almucantar (1993–2018) and hybrid (2014–2018) scan databases by binning U27s in AOD (440 nm), Angström exponent (AE, 440–870 nm), and SSA (440, 675, 870, 1020 nm). Using this LUT approach, the uncertainty estimates U27 for each individual V3 Level 2 retrieval can be obtained by interpolation using the corresponding measured and inverted combination of AOD, AE, and SSA.

Aliaksandr Siniuk↗

NASA NTRS2020

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↗

NASA NTRS2020

Global Estimates and Long-Term Trends of Fine Particulate Matter Concentrations (1998-2018)

Exposure to outdoor fine particulate matter (PM2.5) is a leading risk factor for mortality. We develop global estimates of annual PM2.5 concentrations and trends for 1998–2018 using advances in satellite observations, chemical transport modeling, and ground-based monitoring. Aerosol optical depths (AODs) from advanced satellite products including finer resolution, increased global coverage, and improved long-term stability are combined and related to surface PM2.5 concentrations using geophysical relationships between surface PM2.5 and AOD simulated by the GEOS-Chem chemical transport model with updated algorithms. The resultant annual mean geophysical PM2.5 estimates are highly consistent with globally distributed ground monitors (R2 = 0.81; slope = 0.90). Geographically weighted regression is applied to the geophysical PM2.5 estimates to predict and account for the residual bias with PM2.5 monitors, yielding even higher cross validated agreement (R2 = 0.90–0.92; slope = 0.90–0.97) with ground monitors and improved agreement compared to all earlier global estimates. The consistent long-term satellite AOD and simulation enable trend assessment over a 21 year period, identifying significant trends for eastern North America (−0.28 ± 0.03 μg/m3/yr), Europe (−0.15 ± 0.03 μg/m3/yr), India (1.13 ± 0.15 μg/m3/yr), and globally (0.04 ± 0.02 μg/m3/yr). The positive trend (2.44 ± 0.44 μg/m3/yr) for India over 2005–2013 and the negative trend (−3.37 ± 0.38 μg/m3/yr) for China over 2011–2018 are remarkable, with implications for the health of billions of people.

Melanie S. Hammer↗

NASA NTRS2020

Evaluating impacts of snow, surface water, soil and vegetation on empirical vegetation and snow indices for the Utqiaġvik tundra ecosystem in Alaska with the LVS3 model

Satellite observations for the Arctic and boreal region may contain information of vegetation, soil, snow, snowmelt, and/or other surface water bodies. We investigated the impacts of vegetation, soil, snow and surface water on empirical vegetation/snow indices on a tundra ecosystem area located around Utqiaġvik (formerly Barrow) of Alaska with the Moderate Resolution Imaging Spectrometer (MODIS) images in 2001–2014. Empirical vegetation indices, such as normalized difference vegetation index (NDVI), enhanced vegetation index (EVI), the index of near infrared of vegetation (NIRv), and modified EVI (EVI2), have been used to monitor vegetation. Normalized difference snow index (NDSI) has been widely applied to monitor snow. The vegetation cover fraction (VGCF), the soil cover fraction (SOILCF), the snow cover fraction (SNOWCF), the surface water body cover fraction (WaterBodyCF), the fractional absorption of photosynthetically active radiation (PAR) by vegetation chlorophyll (fAPARchl), the fractional absorption of PAR by non-chlorophyll components of the vegetation (fAPARnon-chl), and the fractional absorption of PAR by the entire canopy (fAPARcanopy) are retrieved with the MODIS images and a coupled Leaf-Vegetation-Soil-Snow-Surface water body radiative transfer model, LVS3. The vegetation indices (NDVI, EVI, EVI2 and NIRv) differ from VGCF, fAPARchl, fAPARnon-chl, and fAPARcanopy. In addition to vegetation, we find that soil, snow and surface water also have impacts on vegetation indices NDVI, EVI (EVI2), and NIRv. Presence of snow makes lower the observed values of NDVI, EVI2 and NIRv. After snowmelt is gone, the vegetation indices (NDVI, EVI, EVI2 and NIRv) linearly decrease with SOILCF and WaterBodyCF, and WaterBodyCF has stronger impacts on these vegetation indices than SOILCF. The relationship between EVI and snow is complicated. NDSI non-linearly increases with SNOWCF, but linearly increases with sum of SNOWCF and WaterBodyCF (sum = 0.5893 × NDSI +0.4342, R2 = 0.976). NDSI linearly decreases with VGCF, and the relationship between NDSI and SOILCF is complex. Retrievals of VGCF, fAPARchl, fAPARnon-chl and fAPARcanopy with the LVS3 model provide alternatives for vegetation monitoring and ecological modeling.

Qingyuan Zhang↗

NASA NTRS2019

First Provisional Land Surface Reflectance Product from Geostationary Satellite Himawari-8 AHI

A provisional surface reflectance (SR) product from the Advanced Himawari Imager (AHI) on-board the new generation geostationary satellite (Himawari-8) covering the period between July 2015 and December 2018 is made available to the scientific community. The Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm is used in conjunction with time series Himawari-8 AHI observations to generate 1-km gridded and tiled land SR every 10 minutes during day time. This Himawari-8 AHI SR product includes retrieved atmospheric properties (e.g., aerosol optical depth at 0.47μm and 0.51μm), spectral surface reflectance (AHI bands 1–6), parameters of the RTLS BRDF model, and quality assurance flags. Product evaluation shows that Himawari-8 AHI data on average yielded 35% more cloud-free, valid pixels in a single day when compared to available data from the low earth orbit (LEO) satellites Terra/Aqua with MODIS sensor. Comparisons of Himawari-8 AHI SR against corresponding MODIS SR products (MCD19A1) over a variety of land cover types with the similar viewing geometry show high consistency between them, with correlation coefficients (r) being 0.94 and 0.99 for red and NIR bands, respectively. The high-frequency geostationary data are expected to facilitate studies of ecosystems on daily to diurnal time scales, complementing observations from networks such as the FLUXNET.

Himawari-8 AHI↗

NASA NTRS2019

Connecting Crop Productivity, Residue Fires, and Air Quality over Northern India

Northwestern India is known as the “breadbasket” of the country producing two-thirds of food grains, with wheat and rice as the principal crops grown under the crop rotation system. Agricultural data from India indicates a 25% increase in the post-monsoon rice crop production in Punjab during 2002–2016. NASA’s A-train satellite sensors detect a consistent increase in the vegetation index (net 21%) and post-harvest agricultural fire activity (net ~60%) leading to nearly 43% increase in aerosol loading over the populous Indo-Gangetic Plain in northern India. The ground-level particulate matter (PM2.5) downwind over New Delhi shows a concurrent uptrend of net 60%. The effectiveness of a robust satellite-based relationship between vegetation index—a proxy for crop amounts, and post-harvest fires—a precursor of extreme air pollution events, has been further demonstrated in predicting the seasonal agricultural burning. An efficient crop residue management system is critically needed towards eliminating open field burning to mitigate episodic hazardous air quality over northern India.

Hiren Jethva↗