Search NASA⌕ Search

SEARCH · Search NASA

Results for “SO2”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Evaluation of SO2, SO4(2−) and an Updated SO2 Dry Deposition Parameterization in the United Kingdom Earth System Model

In this study we evaluate simulated surface SO2 and sulphate (SO4(2-)) concentrations from the United Kingdom Earth System Model (UKESM1) against observations from ground based measurement networks in the USA and Europe for the period 1987 to 2014. We find that UKESM1 captures the historical trend for decreasing concentrations of atmospheric SO2 and SO4(2-) in both Europe and the USA over the period 1987 to 2014. However, in the polluted regions of the eastern USA and Europe, UKESM1 over-predicts surface SO2 concentrations by a factor of 3, while under-predicting surface SO4(2-) concentrations by 25-35%. In the cleaner western USA, the model over-predicts both surface SO2 and SO4(2-) concentrations by a factor of 12 and 1.5 respectively. We find that UKESM1’s bias in surface SO2 and SO4(2-) concentrations is variable according to region and season. We also evaluate UKESM1 against total column SO2 from the Ozone Monitoring Instrument (OMI) using an updated data product. This comparison provides information about the model’s global performance, finding that UKESM1 over predicts total column SO2 over much of the globe, including the large source regions of India, China, the USA and Europe as well as over outflow regions. Finally, we assess the impact of a more realistic treatment of the model’s SO2 dry deposition parameterization. This change increases SO2 dry deposition to the land and ocean surfaces, thus reducing the atmospheric loading of SO2 and SO(sup 2-)(sub 4). In comparison with the ground-based and satellite observations, we find that the modified parameterization reduces the models over prediction of surface SO2 concentrations and total column SO2. Relative to the ground-based observations the simulated surface SO4(2-) concentrations are also reduced, while the simulated SO2 dry deposition fluxes increase.

SO2↗

Version 2 Ozone Monitoring Instrument SO2 Product (OMSO2 V2): New Anthropogenic SO2 Vertical Column Density Dataset

The Ozone Monitoring Instrument (OMI) has been providing global observations of SO2 pollution since 2004. Here we introduce the new anthropogenic SO2 vertical column density (VCD) dataset in the version 2 OMI SO2 product (OMSO2 V2). As with the previous version (OMSO2 V1.3), the new dataset is generated with an algorithm based on principal component analysis of OMI radiances, but features several updates. The most important among those is the use of expanded lookup tables and model a priori profiles to estimate SO2 Jacobians for individual OMI pixels, in order to better characterize pixel-to-pixel variations in SO2 sensitivity, including over snow and ice. Additionally, new data screening and spectral fitting schemes have been implemented to improve the quality of the spectral fit. As compared with the planetary boundary layer SO2 dataset in OMSO2 V1.3, the new dataset has substantially better data quality, especially over areas that are relatively clean or affected by the south Atlantic anomaly. The updated retrievals over snow/ice yield more realistic seasonal changes in SO2 at high latitudes and offer enhanced sensitivity to sources during wintertime. An error analysis has been conducted to assess uncertainties in SO2 VCDs from both the spectral fit and Jacobian calculations. The uncertainties from spectral fitting are reflected in SO2 slant column densities (SCDs) and largely depend on the signal-to-noise ratio of the measured radiances, as implied by the generally smaller SCD uncertainties over clouds or for smaller solar zenith angles. The SCD uncertainties for individual pixels are estimated to be~0.15-0.3 DU (Dobson Units) between ~40°S and ~40°N and to be~0.2-0.5 DU at higher latitudes. The uncertainties from the Jacobians are approximately ~50-100% over polluted areas, and primarily attributed to errors in SO2 a priori profiles and cloud pressures, as well as the lack of explicit treatment for aerosols. Finally, the daily mean and median SCDs over the presumably SO2-free equatorial East Pacific have increased by only~0.0035 DU and ~0.003 DU respectively over the entire 15-year OMI record; while the standard deviation of SCDs has grown by only~0.02 DU or ~10%. Such remarkable long-term stability makes the new dataset particularly suitable for detecting regional changes in SO2 pollution.

OMI, SO2, Remote Sensing↗

SO2 Emission Estimates Using OMI SO2 Retrievals for 2005-2017

SO2 column densities from OMI provide important information on emission trends and missing sources, but there are discrepancies between different retrieval products. We employ three OMI SO2 retrieval products (NASA standard (SP), NASA prototype, and BIRA) to study the magnitude and trend of SO2 emissions. SO2 column densities from these retrievals are most consistent when viewing angles and solar zenith angles are small, suggesting more robust emission estimates in summer and at low latitudes. We then apply a hybrid 4DVar/mass balance emission inversion to derive monthly SO2 emissions from the NASA SP and BIRA products. Compared to HTAPv2 emissions in 2010, both posterior emission estimates are lower in US, India and Southeast China, but show different changes of emissions in North China Plain. The discrepancies between monthly NASA and BIRA posterior emissions in 2010 are less than or equal to 17% in China and 34% in India. SO2 emissions increase from 2005 to 2016 by 35% (NASA) 48% (BIRA) in India, but decrease in China by 23% (NASA) 33% (BIRA) since 2008. Compared to insitu measurements, the posterior GEOSChem surface SO2 concentrations have reduced NMB in China, the US, and India but not in South Korea in 2010. BIRA posteriors have better consistency with the annual growth rate of surface SO2 measurement in China and spatial variability of SO2 concentration in China, South Korea and India, whereas NASA SP posteriors have better seasonality. These evaluations demonstrate the capability to recover SO2 emissions using OMI observations.

Data assimilation↗

High-Resolution Mapping of SO2 Using Airborne Observations from the GeoTASO Instrument During the KORUS-AQ Field Study: PCA-Based Vertical Column Retrievals

The Geostationary Trace gas and Aerosol Sensor Optimization (GeoTASO) instrument is an airborne hyperspectral spectrometer measuring backscattered solar radiation in the ultraviolet (290–400 nm) and visible (415–695 nm) wavelength regions. This paper presents high-resolution sulfur dioxide (SO2) maps over the Korean Peninsula, produced by SO2 retrievals from GeoTASO measurements during the Korea–United States Air Quality Field Study (KORUS-AQ) from May to June 2016. The highly sensitive GeoTASO instrument with a spatial resolution of ~250 m × 250 m can detect point emission sources of SO2 within its fields of view, even without merging multiple overlapping observations. To retrieve SO2 vertical columns from the GeoTASO measurements, we apply an algorithm based on principal component analysis (PCA), which is effective in suppressing noise and biases in SO2 retrievals. The retrievals successfully capture SO2 plumes and various point sources such as power plants, a petrochemical complex, and a steel mill, located in South Chungcheong Province, some of which are not detected by a ground-based in situ measurement network. Spatial distributions of SO2 from GeoTASO observations in source areas are consistent with those from the Stack Tele-Monitoring System reports and airborne in situ SO2 measurements. Comparisons of SO2 retrievals from GeoTASO and existing satellite sensors demonstrate the significance of high-resolution SO2 observations, by indicating that GeoTASO detects small SO2 emission sources that are not precisely resolved by single overpasses of satellites. To assess future geostationary SO2 observations having a higher spatial resolution, we upscale the GeoTASO SO2 retrievals to a spatial resolution of the Geostationary Environment Monitoring Spectrometer (GEMS). Since the upscaled GeoTASO retrievals also detect SO2 plumes clearly, we expect from GEMS to identify even small SO2 emission sources over Asia.

Chong, Heesung↗

Ozone Monitoring Instrument Observations of Interannual Increases in SO2 Emissions from Indian Coal-fired Power Plants During 2005-2012

Due to the rapid growth of electricity demand and the absence of regulations, sulfur dioxide (SO2) emissions from coal-fired power plants in India have increased notably in the past decade. In this study, we present the first interannual comparison of SO2 emissions and the satellite SO2 observations from the Ozone Monitoring Instrument (OMI) for Indian coal-fired power plants during the OMI era of 2005-2012. A detailed unit-based inventory is developed for the Indian coal-fired power sector, and results show that its SO2 emissions increased dramatically by 71 percent during 2005-2012. Using the oversampling technique, yearly high-resolution OMI maps for the whole domain of India are created, and they reveal a continuous increase in SO2 columns over India. Power plant regions with annual SO2 emissions greater than 50 Gg year-1 produce statistically significant OMI signals, and a high correlation (R equals 0.93) is found between SO2 emissions and OMI-observed SO2 burdens. Contrary to the decreasing trend of national mean SO2 concentrations reported by the Indian Government, both the total OMI-observed SO2 and average SO2 concentrations in coal-fired power plant regions increased by greater than 60 percent during 2005-2012, implying the air quality monitoring network needs to be optimized to reflect the true SO2 situation in India.

SO2↗

SO2 Emissions and Lifetimes: Estimates from Inverse Modeling Using In Situ and Global, Space-Based (SCIAMACHY and OMI) Observations

Top-down constraints on global sulfur dioxide (SO2) emissions are inferred through inverse modeling using SO2 column observations from two satellite instruments (SCIAMACHY and OMI). We first evaluated the S02 column observations with surface SO2 measurements by applying local scaling factors from a global chemical transport model (GEOS-Chem) to SO2 columns retrieved from the satellite instruments. The resulting annual mean surface SO2 mixing ratios for 2006 exhibit a significant spatial correlation (r=0.86, slope=0.91 for SCIAMACHY and r=0.80, slope = 0.79 for OMI) with coincident in situ measurements from monitoring networks throughout the United States and Canada. We evaluate the GEOS-Chem simulation of the SO2 lifetime with that inferred from in situ measurements to verity the applicability of GEOS-Chem for inversion of SO2 columns to emissions. The seasonal mean SO2 lifetime calculated with the GEOS-Chem model over the eastern United States is 13 h in summer and 48 h in winter, compared to lifetimes inferred from in situ measurements of 19 +/- 7 h in summer and 58 +/- 20 h in winter. We apply SO2 columns from SCIAMACHY and OMI to derive a top-down anthropogenic SO2 emission inventory over land by using the local GEOS-Chem relationship between SO2 columns and emissions. There is little seasonal variation in the top-down emissions (<15%) over most major industrial regions providing some confidence in the method. Our global estimate for annual land surface anthropogenic SO2 emissions (52.4 Tg S/yr from SCIAMACHY and 49.9 Tg S / yr from OMI) closely agrees with the bottom-up emissions (54.6 Tg S/yr) in the GEOS-Chem model and exhibits consistency in global distributions with the bottom-up emissions (r = 0.78 for SCIAMACHY, and r = 0.77 for OMI). However, there are significant regional differences.

Lee, Chulkyu↗

Dispersion and Lifetime of the SO2 Cloud from the August 2008 Kasatochi Eruption

Hemispherical dispersion of the SO2 cloud from the August 2008 Kasatochi eruption is analyzed using satellite data from the Ozone Monitoring Instrument (OMI) and the Goddard Trajectory Model (GTM). The operational OMI retrievals underestimate the total SO2 mass by 20-30% on 8-11 August, as compared with more accurate offline Extended Iterative Spectral Fit (EISF) retrievals, but the error decreases with time due to plume dispersion and a drop in peak SO2 column densities. The GTM runs were initialized with and compared to the operational OMI SO2 data during early plume dispersion to constrain SO2 plume heights and eruption times. The most probable SO2 heights during initial dispersion are estimated to be 10-12 km, in agreement with direct height retrievals using EISF algorithm and IR measurements. Using these height constraints a forward GTM run was initialized on 11 August to compare with the month-long Kasatochi SO2 cloud dispersion patterns. Predicted volcanic cloud locations generally agree with OMI observations, although some discrepancies were observed. Operational OMI SO2 burdens were refined using GTM-predicted mass-weighted probability density height distributions. The total refined SO2 mass was integrated over the Northern Hemisphere to place empirical constraints on the SO2 chemical decay rate. The resulting lower limit of the Kasatochi SO2 e-folding time is approx.8-9 days. Extrapolation of the exponential decay back in time yields an initial erupted SO2 mass of approx.2.2 Tg on 8 August, twice as much as the measured mass on that day.

Krotkov, N. A.↗

Modeling of 2008 Kasatochi Volcanic Sulfate Direct Radiative Forcing: Assimilation of OMI SO2 Plume Height Data and Comparison with MODIS and CALIOP Observations

Volcanic SO2 column amount and injection height retrieved from the Ozone Monitoring Instrument (OMI) with the Extended Iterative Spectral Fitting (EISF) technique are used to initialize a global chemistry transport model (GEOS-Chem) to simulate the atmospheric transport and lifecycle of volcanic SO2 and sulfate aerosol from the 2008 Kasatochi eruption, and to subsequently estimate the direct shortwave, top-of-the-atmosphere radiative forcing of the volcanic sulfate aerosol. Analysis shows that the integrated use of OMI SO2 plume height in GEOS-Chem yields: (a) good agreement of the temporal evolution of 3-D volcanic sulfate distributions between model simulations and satellite observations from the Moderate Resolution Imaging Spectroradiometer (MODIS) and Cloud-Aerosol Lidar with Orthogonal Polarisation (CALIOP), and (b) an e-folding time for volcanic SO2 that is consistent with OMI measurements, reflecting SO2 oxidation in the upper troposphere and stratosphere is reliably represented in the model. However, a consistent (approx. 25 %) low bias is found in the GEOS-Chem simulated SO2 burden, and is likely due to a high (approx.20 %) bias of cloud liquid water amount (as compared to the MODIS cloud product) and the resultant stronger SO2 oxidation in the GEOS meteorological data during the first week after eruption when part of SO2 underwent aqueous-phase oxidation in clouds. Radiative transfer calculations show that the forcing by Kasatochi volcanic sulfate aerosol becomes negligible 6 months after the eruption, but its global average over the first month is -1.3W/sq m, with the majority of the forcing-influenced region located north of 20degN, and with daily peak values up to -2W/sq m on days 16-17. Sensitivity experiments show that every 2 km decrease of SO2 injection height in the GEOS-Chem simulations will result in a approx.25% decrease in volcanic sulfate forcing; similar sensitivity but opposite sign also holds for a 0.03 m increase of geometric radius of the volcanic aerosol particles. Both sensitivities highlight the need to characterize the SO2 plume height and aerosol particle size from space. While more research efforts are warranted, this study is among the first to assimilate both satellite-based SO2 plume height and amount into a chemical transport model for an improved simulation of volcanic SO2 and sulfate transport.

Wang, J.↗

New-Generation NASA Aura Ozone Monitoring Instrument (OMI) Volcanic SO2 Dataset: Algorithm Description, Initial Results, and Continuation with the Suomi-NPP Ozone Mapping and Profiler Suite (OMPS)

Since the fall of 2004, the Ozone Monitoring Instrument (OMI) has been providing global monitoring of volcanic SO2 emissions, helping to understand their climate impacts and to mitigate aviation hazards. Here we introduce a new-generation OMI volcanic SO2 dataset based on a principal component analysis (PCA) retrieval technique. To reduce retrieval noise and artifacts as seen in the current operational linear fit (LF) algorithm, the new algorithm, OMSO2VOLCANO, uses characteristic features extracted directly from OMI radiances in the spectral fitting, thereby helping to minimize interferences from various geophysical processes (e.g., O3 absorption) and measurement details (e.g., wavelength shift). To solve the problem of low bias for large SO2 total columns in the LF product, the OMSO2VOLCANO algorithm employs a table lookup approach to estimate SO2 Jacobians (i.e., the instrument sensitivity to a perturbation in the SO2 column amount) and iteratively adjusts the spectral fitting window to exclude shorter wavelengths where the SO2 absorption signals are saturated. To first order, the effects of clouds and aerosols are accounted for using a simple Lambertian equivalent reflectivity approach. As with the LF algorithm, OMSO2VOLCANO provides total column retrievals based on a set of predefined SO2 profiles from the lower troposphere to the lower stratosphere, including a new profile peaked at 13 km for plumes in the upper troposphere. Examples given in this study indicate that the new dataset shows significant improvement over the LF product, with at least 50% reduction in retrieval noise over the remote Pacific. For large eruptions such as Kasatochi in 2008 (approximately 1700 kt total SO2/ and Sierra Negra in 2005 (greater than 1100DU maximum SO2), OMSO2VOLCANO generally agrees well with other algorithms that also utilize the full spectral content of satellite measurements, while the LF algorithm tends to underestimate SO2. We also demonstrate that, despite the coarser spatial and spectral resolution of the Suomi National Polar-orbiting Partnership (Suomi-NPP) Ozone Mapping and Profiler Suite (OMPS) instrument, application of the new PCA algorithm to OMPS data produces highly consistent retrievals between OMI and OMPS. The new PCA algorithm is therefore capable of continuing the volcanic SO2 data record well into the future using current and future hyperspectral UV satellite instruments.

OMI↗

Using Machine Learning to Estimate Surface-Level SO2 Concentrations from Satellite-Based Measurements

Sulfur dioxide (SO2) is a criteria air pollutant due to its contributions to aerosol formation, rainfall acidification, and harm to human health. The placement of air quality monitoring sites is typically biased towards urban areas, leaving large areas with very limited monitoring data. The Ozone Monitoring Instrument (OMI) has been used to provide estimates of SO2 vertical column densities (VCDs) globally at spatial resolution of 10s of kms once per day. OMI SO2 VCDs have been previously used to estimate surface SO2 concentrations using chemical transport model (CTM) simulations. The CTMs use estimated emissions and assimilated meteorological data, and simulate the chemical and physical processes that determine the vertical profile of SO2, which can be used to derive a ratio between the surface concentrations and VCDs. These models are complex, computationally expensive, and have large uncertainties in the simulated surface-to-VCD ratio due to biases in emissions and relatively coarse resolution. Machine learning techniques are comparatively easier to use, much less computationally expensive to use after training, and can produce more accurate estimations of surface concentrations than the CTM-based method. The interpretation of machine learning models often poses challenges, and in some cases, non-physical variables unrelated to SO2 are used as predictors. In this work, we create an artificial neural network (ANN) to relate OMI retrievals and archived GEOS-FP boundary layer heights to surface SO2 concentrations from the ChinaHighAirPollutants ChinaHighSO2 dataset (CHAP; Wei et al., 2023) on a seasonal average timescale from 2013-2018. Our model only utilizes five variables that are directly relevant to the satellite retrieval, lifetime, and spatial distribution of SO2. The model was trained on 16 seasons (four of each) with independent validation (one of each season) and testing datasets (one of each season) to avoid overfitting. Our ANN generates surface SO2 concentrations that are sensitive (slope = 0.51) and consistent (r = 0.74) with the CHAP data, but are underpredicted by an average of 1.2 ppbv with a mean absolute error of 2.2 ppbv. These results are better than recent studies utilizing the CTM method. To our knowledge, this is the best performing machine learning model that only uses physical variables to predict surface SO2. Our work demonstrates that a carefully constructed, simple ML model can accurately estimate surface-based SO2 concentrations from satellite VCD measurements, and this technique has future promise to expend to newer, higher resolution satellites and other air pollutants.

SO2, air quality, OMI, machine learning↗

The atmospheric SO2 budget for Pinatubo derived from NOAA-11 SBUV/2 spectral data

Spectral scan data from the NOAA-11 SBUV/2 instrument were used to derive SO2 for three days following the eruption of Mt. Pinatubo - June 19, July 1, and July 17, 1991. Band structure between 300 and 310 nm observed in the backscattered albedo uniquely identifies the presence of SO2. Band ratios are used to infer SO2 amounts to better accuracy (10-20 percent) and sensitivity (about 0.5 milli-atm-cm of SO2) than the TOMS retrieval, but with relatively poor spatial coverage because the measurement is nadir only. Only 7 scans showed detectable SO2 on June 19 when the cloud was still very localized. On July 1 there were 29 scans between 35N and 12S with SO2, with the highest concentration detected over the Atlantic, and on July 17 SO2 was detected in 30 scans around the world, but in decreased concentration. Estimates of the total SO2 budget made after the cloud had spread sufficiently for the sparse SBUV/2 sampling to be adequate indicated that there were 8.4 million metric tons (MMT) of SO2 in the stratosphere on July 1, 1991, and 4.1 MMT remaining on July 17. This corresponds to an e-folding time of about 24 days for the conversion of SO2 to aerosol, and is consistent with an initial injection into the stratosphere of 12-15 MMT of SO2.

Mcpeters, Richard D.↗

Next-Generation Aura/OMI NO2 and SO2 Products

The measurement of both SO2 and NO2 gases are recognized as an essential component of atmospheric composition missions. We describe current capabilities and limitations of the operational Aura/OMI NO2 and SO2 data that have been used by a large number of researchers. Analyses of the data and validation studies have brought to light a number of areas in which these products can be expanded and improved. Major improvements for new NASA standard (SP) NO2 product include more accurate tropospheric and stratospheric column amounts, along with much improved error estimates and diagnostics. Our approach uses a monthly NO2 climatology based on the NASA Global Modeling Initiative (GMI) chemistry-transport model and takes advantage of OMI data from cloudy scenes to find clean areas where the contribution from the trap NO2 column is relatively small. We then use a new filtering, interpolation and smoothing techniques for separating the stratospheric and tropospheric components of NO2, minimizing the influence of a priori information. The new algorithm greatly improves the structure of stratospheric features relative to the original SP. For the next-generation OMI SO2 product we plan to implement operationally the offline iterative spectral fitting (ISF) algorithm and re-process the OMI Level-2 SO2 dataset using a priori SO2 and aerosol profiles, clouds, and surface reflectivity appropriate for observation conditions. This will improve the ability to detect and quantify weak tropospheric SO2 loadings. The new algorithm is validated using aircraft in-situ data during field campaigns in China (2005 and 2008) and in Maryland (Frostburg, 2010 and DISCOVER-AQ in July 2011). The height of the SO2 plumes will also be estimated for high SO2 loading cases (e.g., volcanic eruptions). The same SO2 algorithm will be applied to the data from OMPS sensor to be launched on NPP satellite later this year. The next-generation NO2 and SO2 products will provide critical information (e.g., averaging kernels) for evaluation of chemistry-transport models, for data assimilation, and to impose top-down constraints on the SO2 and NO2 emission sources.

Krotkov, Nickolay↗

SO2 over Central China: Measurements, Numerical Simulations and the Tropospheric Sulfur Budget

SO2 in central China was measured in situ from an aircraft and remotely using the Ozone Monitoring Instrument (OMI) from the Aura satellite; results were used to develop a numerical tool for evaluating the tropospheric sulfur budget - sources, sinks, transformation and transport. In April 2008, measured ambient SO2 concentrations decreased from approx.7 ppbv near the surface to approx. 1 ppbv at 1800 m altitude (an effective scale height of approx.800 m), but distinct SO2 plumes were observed between 1800 and 4500 m, the aircraft's ceiling. These free tropospheric plumes play a major role in the export of SO2 and in the accuracy of OMI retrievals. The mean SO2 column contents from aircraft measurements (0.73 DU, Dobson Units) and operational OMI SO2 products (0.63+/-0.26 DU) were close. The OMI retrievals were well correlated with in situ measurements (r = 0.84), but showed low bias (slope = 0.54). A new OMI retrieval algorithm was tested and showed improved agreement and bias (r = 0.87, slope = 0.86). The Community Multiscale Air Quality (CMAQ) model was used to simulate sulfur chemistry, exhibiting reasonable agreement (r = 0.62, slope = 1.33) with in situ SO2 columns. The mean CMAQ SO2 loading over central and eastern China was 54 kT, approx.30% more than the estimate from OMI SO2 products, 42 kT. These numerical simulations, constrained by observations, indicate that ",50% (35 to 61 %) of the anthropogenic sulfur emissions were transported downwind, and the overall lifetime of tropospheric SO2 was 38+/-7 h.

He, Hao↗

Response of SO2 and Particulate Air Pollution to Local and Regional Emission Controls: A Case Study in Maryland

This paper addresses the questions of what effect local regulations can have on pollutants with different lifetimes and how surface observations and remotely sensed data can be used to determine the impacts. We investigated the decadal trends of tropospheric sulfur dioxide (SO2) and aerosol pollution over Maryland and its surrounding states, using surface, aircraft, and satellite measurements. Aircraft measurements indicated fewer isolated SO2 plumes observed in summers, a 40 decrease of column SO2, and a 20 decrease of atmospheric optical depth (AOD) over Maryland after the implementation of local regulations on sulfur emissions from power plants (90 reduction from 2010). Surface observations of SO2 and particulate matter (PM) concentrations in Maryland show similar trends. OMI SO2 and MODIS AOD observations were used to investigate the column contents of air pollutants over the eastern U.S.; these indicate decreasing trends in column SO2 (60 decrease) and AOD (20 decrease). The decrease of upwind SO2 emissions also reduced aerosol loadings over the downwind Atlantic Ocean near the coast by 20, while indiscernible changes of the SO2 column were observed. A step change of SO2 emissions in Maryland starting in 20092010 had an immediate and profound benefit in terms of local surface SO2 concentrations but a modest impact on aerosol pollution, indicating that short-lived pollutants are effectively controlled locally, while long-lived pollutants require regional measures.

satellite measurements↗

The Role of the SO2 Radiative Effect in Sustaining the Volcanic Winter and Soothing the Toba Impact on Climate

Volcanic eruptions are an important climate driver. The impact of Pinatubo-sized eruptions has been observed and is well constrained. The magnitude and duration of volcanic winter effects after super-eruptions such as Toba remain disputed due to disagreement between the strong cooling predicted by models, and much milder climate perturbations according to the paleo data. Here we present a re-evaluated climate impact of a Toba-sized super-eruption based on up-to-date GISS ModelE simulations. In this study, we account for all known primary mechanisms that govern the evolution of the volcanic plume and their nonlinear interactions. The SO2 radiative effects are evaluated for the first time in coupled climate simulations with the interactive atmospheric chemistry module. We found that SO2 effects on photochemistry, dynamics, and radiative forcing are especially prominent. Due to strong absorption in ultraviolet, SO2 feedback on photochemistry partially offsets the limiting effect associated with aerosol microphysical processes. SO2 greenhouse warming soothes the radiative cooling exerted by sulfate aerosols. SO2 absorption in the shortwave and longwave causes radiative heating and lofting of the volcanic plume, and boosts the efficiency of SO2 impact on photochemistry. Our analysis shows that SO2 lifetime and magnitude of effects scale up and increase with the amount of emitted material. For a Pinatubo-sized eruption, SO2 feedbacks on chemistry and dynamics are relevant only during the initial stage of the volcanic plume evolution, while local SO2 concentrations are high. For a Toba-sized eruption, SO2 effects are as important as sulfate aerosols, and produce a less extreme volcanic winter.

Sergey Osipov↗

A New Machine Learning Based Analysis for Improving Satellite Retrieved Atmospheric Composition Data: OMI SO2 as an Example

Despite recent progress, satellite retrievals of anthropogenic SO2 still suffer from relatively low signal-tonoise ratios. In this study, we demonstrate a new machine learning data analysis method to improve the quality of satellite SO2 products. In the absence of large ground-truth datasets for SO2, we start from SO2 slant column densities (SCDs) retrieved from the Ozone Monitoring Instrument (OMI) using a data-driven, physically based algorithm and calculate the ratio between the SCD and the root mean square (rms) of the fitting residuals for each pixel. To build the training data, we select presumably clean pixels with small SCD / rms ratios (SRRs) and set their target SCDs to zero. For polluted pixels with relatively large SRRs, we set the target to the original retrieved SCDs. We then train neural networks (NNs) to reproduce the target SCDs using predictors including SRRs for individual pixels, solar zenith, viewing zenith and phase angles, scene reflectivity, and O3 column amounts, as well as the monthly mean SRRs. For data analysis, we employ two NNs: (1) one trained daily to produce analyzed SO2 SCDs for polluted pixels each day and (2) the other trained once every month to produce analyzed SCDs for less polluted pixels for the entire month. Test results for 2005 show that our method can significantly reduce noise and artifacts over background regions. Over polluted areas, the monthly mean NN-analyzed and original SCDs generally agree to within ±15 %, indicating that our method can retain SO2 signals in the original retrievals except for large volcanic eruptions. This is further confirmed by running both the NN-analyzed and original SCDs through a topdown emission algorithm to estimate the annual SO2 emissions for ∼ 500 anthropogenic sources, with the two datasets yielding similar results. We also explore two alternative approaches to the NN-based analysis method. In one, we employ a simple linear interpolation model to analyze the original SCD retrievals. In the other, we develop a PCA–NN algorithm that uses OMI measured radiances, transformed and dimension-reduced with a principal component analysis (PCA) technique, as inputs to NNs for SO2 SCD retrievals. While the linear model and the PCA–NN algorithm can reduce retrieval noise, they both underestimate SO2 over polluted areas. Overall, the results presented here demonstrate that our new data analysis method can significantly improve the quality of existing OMI SO2 retrievals. The method can potentially be adapted for other sensors and/or species and enhance the value of satellite data in air quality research and applications.

Can Li↗

Calculating the Height of Volcanic Cloud SO2 With a Lagrangian Trajectory Tool: Raikoke (2019) and Pinatubo (1991) Cases

We have developed a new data analysis tool to reconstruct the altitude of SO2 clouds ejected by avolcanic eruption. Starting with 2D fields of satellite observed SO2 column density, known volcanolocation, and reanalysis wind fields, the Goddard Lagrangian trajectory tool allows us to estimatethe altitude and local concentration of SO2 at the time of observation. We demonstrate this tool forthe June 21, 2019 Mt. Raikoke eruption (8h, 30h and 78 hours after eruption). We use SO2 data fromthe Ozone Mapping and Profiler Suite/Nadir Mapper (OMPS/NM) onboard the NASA-NOAASuomi National Polar Partnership satellite and the Tropospheric Monitoring Instrument(TROPOMI) onboard the European Copernicus Sentinel 5 precursor satellite and wind fields fromthe NASA Goddard Earth Observing System (GEOS) model. We obtained a wide distribution ofSO2 altitudes between 1 and 19 kilometers in different parts of the Raikoke SO2 clouds, with thehighest SO2 concentration between 11 and 16 km, in good agreement with data from independentSO2 layer height retrievals from TROPOMI and Infrared Atmospheric Sounding Interferometer(IASI). We then applied this method to the June 15, 1991 Mt. Pinatubo eruption using SO2 columnmeasurements from the NASA Total Ozone Mapping Spectrometer (TOMS) on June 16, 17, 18 andusing wind fields from the National Centers for Environmental Prediction Reanalysis version 2. Weobtained a wide distribution of SO2 heights from 6 km up to 30 km in the northern part of the cloud.We use our new trajectory tool to reconstruct volcanic SO2/aerosol emission sources as functions oftime and altitude, providing inputs to the GEOS model aerosol and chemistry schemes in order tosimulate the dispersion of volcanic SO2 and aerosol clouds from these events.

SO2↗

Continuing Global SO2 Data Record from OMI and SNPP/OMPS to JPSS-1/NOAA-20/OMPS

Since 2004, the Ozone Monitoring Instrument (OMI) aboard NASA's Earth Observing System (EOS) Aura spacecraft has been providing global observations that help to constrain the sources, transport, and environmental impacts of anthropogonic and volcanic SO2. The OMI SO2 data record is now being continued with the NASA/NOAA Suomi National Polar-orbiting Partnership (SNPP)/Ozone Mapping and Profiler Suite (OMPS) launched in 2011. Both OMI and SNPP/OMPS SO2 products are produced with the Goddard principal component analysis (PCA) based spectral fitting algorithm. This data-driven technique inherently accounts for various instrumental factors and geophysical interferences, leading to high-quality, consistent SO2 retrievals between OMI and SNPP/OMPS, despite coarser spectral (~0.5 nm vs. ~1 nm) and spatial (13  24 km2 vs. 50  50 km2 at nadir) resolution for the latter. In this presentation, we describe our effort to continue the long-term SO2 climate data record using measurements from the Joint Polar Satellite System (JPSS)-1/NOAA-20 (N20)/OMPS. Launched in 2017, the N20/OMPS is a follow-on for SNPP/OMPS but features a spatial resolution (17  13 km2) that is comparable with OMI. We will discuss our progress implementing the PCA SO2 algorithm with N20/OMPS, especially algorithmic improvements to further reduce retrieval noise and bias for large volcanic eruptions. We will present examples for both continuously emitting sources (e.g., power plants in India and oil/gas fields in the Middle East) and volcanic eruptions (e.g., Raikoke in 2019). We will also compare N20/OMPS SO2 retrievals with OMI and SNPP/OMPS, as well as other instruments such as the ESA Copernicus Sentinel-5 Precursor (S5P)/TROPOspheric Monitoring Instrument (TROPOMI). To assess the ability of N20/OMPS to monitor and quantify SO2 sources, we will run the level 2 retrievals through a top-down emission algorithm to estimate the SO2 emission strengths for a number of point sources. Finally, we will outline our plan for further algorithm refinement and public data release.

SO2↗