Engineering PapersSearch

Engineering topics

Nickolay Krotkov

Publications and source records attributed to Nickolay Krotkov.

At least 19 records

Inverse Modeling of the Initial Stage of the 1991 Pinatubo Volcanic Cloud Accounting for Radiative Feedback of Volcanic Ash

The way volcanic clouds evolve is very sensitive to the initial spatial 3D distributions of volcanic materials, which are often unknown. In this study, we conducted inverse modeling of the Mt. Pinatubo cloud using TOMS 2D mapping of Aerosol Index and SO 2 loading during the first three post-eruption days to estimate the time-dependent emissions profiles and initial 3D spatial distributions of volcanic ash and SO 2 . We account for aerosol radiative feedback and dynamic lofting of volcanic ash in the inversion calculations for the first time. This resulted in a lower ash injection height (by 1.5 km for ash) than without ash radiative feedback. The Pinatubo eruption ejected ≈77% of fine ash at 12 to 23 km, ≈65% of SO 2 at 18–25 km. In contrast with previous studies, which suggested that all volcanic materials were emitted above the tropopause, a significant fraction of SO 2 (5.1 of 15.5 Mt) and fine ash (37.2 of 66.5 Mt) were ejected in the troposphere, where SO 2 quickly oxidized into sulfate aerosol that is short-lived in the troposphere. This explains the early presence of sulfate aerosols in the plume and why the models can reproduce the observed volcanic aerosols' optical depth (AOD), assuming lower-than-observed SO 2 emission in the stratosphere. Despite the quicker than in observations build-up of sulfate AOD, in a month after the eruption, the evolution of the Pinatubo AOD simulated using the obtained ash and SO 2 initial distributions converges with the available stratospheric aerosol and gas experiment (SAGE) observations.

Pinatubo volcanic ash and SO2 emissions

Use of Machine Learning and Principal Component Analysis to Retrieve Nitrogen Dioxide (NO 2 ) With Hyperspectral Imagers and Reduce Noise in Spectral Fitting

Nitrogen dioxide (NO 2 ) is an important trace-gas pollutant and climate agent whose presence also leads to spectral interference in ocean color retrievals. NO 2 column densities have been retrieved with satellite UV–Vis spectrometers such as the Ozone Monitoring Instrument (OMI) and the Tropospheric Monitoring Instrument (TROPOMI) that typically have spectral resolutions of the order of 0.5 nm or better and spatial footprints as small as 3.6 km × 5.6 km. These NO 2 observations are used to estimate emissions, monitor pollution trends, and study effects on human health. Here, we investigate whether it is possible to retrieve NO 2 amounts with lower-spectral-resolution hyperspectral imagers such as the Ocean Color Instrument (OCI) that will fly on the Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) satellite set for launch in early 2024. OCI will have a spectral resolution of 5 nm and a spatial resolution of ∼ 1 km with global coverage in 1–2 d. At this spectral resolution, small-scale spectral structure from NO 2 absorption is still present. We use real spectra from the OMI to simulate OCI spectra that are in turn used to estimate NO 2 slant column densities (SCDs) with an artificial neural network (NN) trained on target OMI retrievals. While we obtain good results with no noise added to the OCI simulated spectra, we find that the expected instrumental noise substantially degrades the OCI NO 2 retrievals. Nevertheless, the NO 2 information from OCI may be of value for ocean color retrievals. OCI retrievals can also be temporally averaged over timescales of the order of months to reduce noise and provide higher-spatial-resolution maps that may be useful for downscaling lower-spatial-resolution data provided by instruments such as OMI and TROPOMI; this downscaling could potentially enable higher-resolution emissions estimates and be useful for other applications. In addition, we show that NNs that use coefficients of leading modes of a principal component analysis of radiance spectra as inputs appear to enable noise reduction in NO 2 retrievals. Once trained, NNs can also substantially speed up NO 2 spectral fitting algorithms as applied to OMI, TROPOMI, and similar instruments that are flying or will soon fly in geostationary orbit.

NO2

Version 2 of the global catalogue of large anthropogenic and volcanic SO 2 sources and emissions derived from satellite measurements

Sulfur dioxide (SO 2 ) measurements from the Ozone Monitoring Instrument (OMI), Ozone Mapping and Profiler Suite (OMPS), and TROPOspheric Monitoring Instrument (TROPOMI) satellite spectrometers were used to update and extend the previously developed global catalogue of large SO 2 emission sources. This version 2 of the global catalogue covers the period of 2005–2021 and includes a total of 759 continuously emitting point sources releasing from about 10 kt yr −1 to more than 4000 kt yr −1 of SO 2 , that have been identified and grouped by country and primary source origin: volcanoes (106 sources); power plants (477); smelters (74); and sources related to the oil and gas industry (102). There are several major improvements compared to the original catalogue: it combines emissions estimates from three satellite instruments instead of just OMI, uses a new version 2 of the OMI and OMPS SO 2 dataset, and updated consistent site-specific air mass factors (AMFs) are used to calculate SO 2 vertical column densities (VCDs). The newest TROPOMI SO 2 data processed with the Covariance-Based Retrieval Algorithm (COBRA), used in the catalogue, can detect sources with emissions as low as 8 kt yr −1 (in 2018–2021) compared to the 30 kt yr −1 limit for OMI. In general, there is an overall agreement within ±12 % in total emissions estimated from the three satellite instruments for large regions. For individual emission sources, the spread is larger: the annual emissions estimated from OMI and TROPOMI agree within ±13 % in 50 % of cases and within ±28 % in 90 % of cases. The version 2 catalogue emissions were calculated as a weighted average of emission estimates from the three satellite instruments using an inverse-variance weighting method. OMI, OMPS, and TROPOMI data contribute 7 %, 5 %, and 88 % to the average, respectively, for small (<30 kt yr −1 ) sources and 33 %, 20 %, and 47 %, respectively, for large (>300 kt yr −1 ) sources. The catalogue data show an approximate 50 % decline in global SO 2 emissions between 2005 and 2021, although emissions were relatively stable during the last 3 years. The version 2 of the global catalogue has been posted at the NASA global SO 2 monitoring website (https://doi.org/10.5067/MEASURES/SO2/DATA406, Fioletov et al., 2022).

volcanic

Using Machine Learning for Timely Estimates of Ocean Color Information From Hyperspectral Satellite Measurements in the Presence of Clouds, Aerosols, and Sunglint

Retrievals of ocean color from space are important for better understanding of the ocean ecosystem but can be limited under conditions such as clouds, aerosols, and sunglint. Many ocean color algorithms use a few selected spectral bands to perform an atmospheric correction and then derive the upwelling radiance from the ocean. The limitations in the atmospheric correction under certain conditions lead to many gaps in daily spatial coverage of ocean color retrievals. To address these limitations, we introduce a new approach that uses machine learning to estimate ocean color from top of atmosphere radiances or reflectance measurements. In this approach, a principal component analysis is used to decompose the hyperspectral measurements into spectral features that describe the scattering and absorption of the atmosphere and the underlying surface. The coefficients of the principal components are then used to train a neural network to predict ocean color properties derived from the MODIS atmospheric correction algorithm. This machine learning approach is independent of a priori information and does not rely on any radiative transfer modeling. We apply the approach to two hyperspectral UV/VIS instruments, the ozone monitoring instrument (OMI) and the TROPOspheric Monitoring Instrument (TROPOMI), using measurements from 320–500 nm to show that it can be used to reproduce ocean color properties in less-than-ideal conditions. This machine learning approach complements the current atmospheric correction ocean color retrievals by filling in the gaps resulting from cloud, aerosol, and sunglint contamination. This method can be applied to the future hyperspectral Ocean Color Instrument (OCI), which will be onboard NASA’s Plankton, Aerosol Cloud, ocean Ecosystem (PACE) ocean color satellite set to launch in 2024.

Ocean color

Tracking aerosols and SO2 clouds from the Raikoke eruption: 3D view from satellite observations

The 21 June 2019 eruption of the Raikoke volcano (Kuril Islands, Russia; 48° N, 153° E) produced significant amounts of volcanic aerosols (sulfate and ash) and sulfur dioxide (SO2) gas that penetrated into the lower stratosphere. The dispersed SO2 and sulfate aerosols in the stratosphere were still detectable by multiple satellite sensors for many months after the eruption. For this study of SO2 and aerosol clouds we use data obtained from two of the Ozone Mapping and Profiler Suite sensors on the Suomi National Polar-orbiting Partnership satellite: total column SO2 from the Nadir Mapper and aerosol extinction profiles from the Limb Profiler as well as other satellite data sets. We evaluated the limb viewing geometry effect (the “arch effect”) in the retrieval of the LP standard aerosol extinction product at 674 nm. It was shown that the amount of SO2 decreases with a characteristic period of 8–18 d and the peak of stratospheric aerosol optical depth recorded at a wavelength of 674 nm lags the initial peak of SO2 mass by 1.5 months. Using satellite observations and a trajectory model, we examined the dynamics of an unusual atmospheric feature that was observed, a stratospheric coherent circular cloud of SO2 and aerosol from 18 July to 22 September 2019.

Eruption of the Raikoke volcano

A Sulfur Dioxide Covariance-Based Retrieval Algorithm (COBRA): Application to TROPOMI Reveals New Emission Sources

Sensitive and accurate detection of sulfur dioxide (SO 2 ) from space is important for monitoring and estimating global sulfur emissions. Inspired by detection methods applied in the thermal infrared, we present here a new scheme to retrieve SO 2 columns from satellite observations of ultraviolet back-scattered radiances. The retrieval is based on a measurement error covariance matrix to fully represent the SO 2 -free radiance variability, so that the SO 2 slant column density is the only retrieved parameter of the algorithm. We demonstrate this approach, named COBRA, on measurements from the TROPOspheric Monitoring Instrument (TROPOMI) aboard the Sentinel-5 Precursor (S-5P) satellite. We show that the method reduces significantly both the noise and biases present in the current TROPOMI operational DOAS SO 2 retrievals. The performance of this technique is also benchmarked against that of the Principal Component Algorithm (PCA) approach. We find that the quality of the data is similar and even slightly better with the proposed COBRA approach. The ability of the algorithm to retrieve SO 2 accurately is also further supported by comparison with ground-based observations. We illustrate the great sensitivity of the method with a high-resolution global SO 2 map, considering two and a half years of TROPOMI data. In addition to the known sources, we detect many new SO 2 emission hotspots worldwide. For the largest sources, we use the COBRA data to estimate SO 2 emission rates. Results are comparable to other recently published TROPOMI-based SO 2 emissions estimates, but the associated uncertainties are significantly lower than with the operational data. Next, for a limited number of weak sources, we demonstrate the potential of our data for quantifying SO 2 emissions with a detection limit of about 8 kt yr -1 , a factor of 4 better than the emissions derived from the Ozone Monitoring Instrument (OMI). We anticipate that the systematic use of our TROPOMI COBRA SO 2 column data set at a global scale will allow identifying and quantifying missing sources, and help improving SO 2 emission inventories.

SO2

Inconsistencies in Sulphur Dioxide Emissions from the Canadian Oil Sands and Potential Implications

Satellite-derived and reported sulfur dioxide (SO2) emissions from the Canadian oil sands are shown to have been consistent up to 2013. Post-2013, these sources of emissions data diverged, with reported emissions dropping by a factor of two, while satellite-derived emissions for the region remained relatively constant, with the discrepancy (satellite-derived emissions minus reported emissions) peaking at 50 kt(SO2) yr−1 around 2016. The 2013–2014 period corresponds to when new flue-gas desulfurization units came on-line. Previous work has established a high level of consistency between at-stack SO2 emissions observations and satellite estimates, and surface monitoring network SO2 concentrations over the same multi-year period show similar trends as the satellite data, with a slight increase in concentrations post-2013. No clear explanation for this discrepancy currently exists. The implications of the discrepancy towards estimated total sulfur deposition to downwind ecosystems were estimated relative to 2013 emissions levels, with the satellite-derived values leaving the area of regional critical load exceedances of aquatic ecosystems largely unchanged from 2013 values, 335 000 km2, and reported values potentially decreasing this area to 185,000 km2.

Sulphur dioxide (SO2) emissions

Global Distribution and 14-Year Changes in Erythemal Irradiance, UV Atmospheric Transmission, and Total Column Ozone 2005 – 2018 Estimated from OMI and EPIC Observations

Satellite data from the Ozone Measuring Instrument (OMI) and Earth Polychromatic Imaging Camera (EPIC) are used to study changes and global distribution of UV erythemal irradiance in mW/m(exp 2) E(ζ,φ,z,t) and the dimensionless UV index (E/25 mWm(exp 2)) over the Earth’s surface as a function of latitude ζ, longitude φ, altitude z, and time t. OMI time series data starting in January 2005 to December 2018 are used to estimate 14-year changes in erythemal irradiance, total column ozone TCO3, scene reflectivity LER, and aerosol absorption τA at specific land sites in the Northern and Southern Hemispheres using publicly available OMI data. Estimates of changes in atmospheric transmission ΔCτ(ζ,φ,z,t) derived from OMI measured cloud and haze reflectivity show an increase between 60°S to 45°S almost no average 14-year change from 55°S to 30°N but show an increase from 35°N to 55°N. This implies increased solar insolation at higher latitudes that suggests positive feedback for global warming. Changes in aerosol absorption of UV irradiance ΔCA(ζ,φ,z,t) derived from changes in τA are mostly in the equatorial region and at higher southern latitudes resulting in correlated changes in E. TCO3 derived from OMI data has increased at a rate of 2% per decade for the latitudes between 60°S to 10°N changing to a decrease of 1% per decade between 40°N to 60°N. The result is an average decrease in E(ζ,φ,z,t) at a rate of 2% per decade in the Southern Hemisphere and an increase between 40°N to 60°N. For some specific sites (latitudes from 55°S to 45°N) there has been little or no change in E(ζ,φ,z,t) for the period 2005 – 2018. Nearly half the sites show strong correlation effects of both short- and long-term cloud and absorbing aerosol change as well as anticorrelation with total column ozone change ΔTCO3. Synoptic EPIC data from the sunlit Earth are used to derive ozone and reflectivity needed for global images of the distribution of E(ζ,φ,z,t) from sunrise to sunset centered on the Americas, Europe-Africa, and Asia. EPIC data are used to show the latitudinal distribution of E(ζ,φ,z,t) from the equator to 75° for specific longitudes. Dangerously high amounts of erythemal irradiance (12 < UV index < 18) are found for many low latitude and high-altitude sites (e.g., San Pedro, Chile (2.45 km), La Paz, Bolivia (3.78 km). Lower UV indices at some equatorial or high-altitude sites (e.g., Quito, Ecuador) are moderated by the presence of persistent cloud effects. High UVI levels (UVI > 6) are also found at most mid-latitude sites during the summer months.

Ozone Measuring Instrument (OMI) and Earth Polychr

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

Lightning NOx Production in the Tropics as Determined Using OMI NO2 Retrievals and WWLLN Stroke Data

Nitrogen oxide (NOx) production by lightning in the tropics is estimated using tropospheric NOx amounts (LNOx*) over deep convective grid boxes derived from OMI (Ozone Monitoring Instrument) nitrogen dioxide (NO2) slant columns and detection efficiency adjusted WWLLN (World Wide Lightning Location Network) flashes. The lightning NOx production efficiency (LNOx PE) in the tropics is determined for the austral and boreal summers of 2007 to 2011 by regressing regional mean daily values of LNOx* for individual seasons against daily flash totals during flash windows prior to the OMI overpass. LNOx PE is determined to be approximately two times larger over marine locations than over continental locations possibly because marine flashes are more energetic. Overall, the mean LNOx PE for the tropics is calculated to be 170 ± 100 mol per flash with values over the tropical Pacific (low flash rate region) being largest. The main contributors to uncertainties in PE are uncertainties in WWLLN flash detection efficiency, upper tropospheric NOx lifetime in the near field of convection, and air mass factor biases.

Dale J. Allen

Midlatitude Lightning NOx Production Efficiency Inferred From OMI and WWLLN Data

Oxides of nitrogen are critical trace gases in the troposphere and are precursors for nitrate aerosol and ozone, which is an important pollutant and greenhouse gas. Lightning is the major source of NOx(NO + NO2) in the mid-to upper troposphere. We estimate the production efficiency (PE) of lightning NOx (LNOx) using satellite data from the Ozone Monitoring Instrument (OMI) and the ground-based WorldWide Lightning Location Network (WWLLN) in three northern midlatitude, primarily continental regions that include much of North America, Europe and East Asia. Data were obtained over 5 boreal summers, 2007 –2011 and comprise the largest number of midlatitude convective events to date for estimating the LNOx PE with satellite NO2and ground-based lightning measurements. In contrast to some previous studies, the algorithm assumes no minimum flash-rate threshold and estimates freshly produced LNOxby subtracting a background of aged NOx estimated from the OMI dataset itself. We infer an average value of 180 ± 100 moles LNOx produced per lightning flash. We also show evidence of a dependence of PE on lightning flash rate and find an approximate empirical power function relating moles LNOxto flashes. PE decreases by an order of magnitude for a 2-order of magnitude increase in flash rate. This phenomenon has not been reported in previous satellite LNOxstudies but is consistent with ground-based observations suggesting an inverse relationship between flash rate and size.

LNOx

High‐Resolution Mapping of Nitrogen Dioxide With TROPOMI: First Results and Validation Over the Canadian Oil Sands

TROPOspheric Monitoring Instrument (TROPOMI), on-board the Sentinel-5 Precursersatellite, is a nadir-viewing spectrometer measuring reflected sunlight in the ultraviolet, visible,near-infrared, and shortwave infrared. From these spectra several important air quality and climate-relatedatmospheric constituents are retrieved, including nitrogen dioxide (NO2) at unprecedented spatialresolution from a satellite platform. We present the first retrievals of TROPOMI NO2over the Canadian OilSands, contrasting them with observations from the Ozone Monitoring Instrument satellite instrument,and demonstrate TROPOMI's ability to resolve individual plumes and highlight its potential for derivingemissions from individual mining facilities. Further, the first TROPOMI NO2validation is presented,consisting of aircraft and surface in situ NO2observations, and ground-based remote-sensingmeasurements between March and May 2018. Our comparisons show that the TROPOMI NO2verticalcolumn densities are highly correlated with the aircraft and surface in situ NO2observations, and theground-based remote-sensing measurements with a low bias (15–30 %); this bias can be reduced byimproved air mass factors. Plain Language Summary: Nitrogen dioxide (NO2) is a pollutant that is linked to respiratoryhealth issues and has negative environmental impacts such as soil and water acidification. Near thesurface the most significant sources of NO2are fossil fuel combustion and biomass burning. With a recentlylaunched satellite instrument (TROPOspheric Monitoring Instrument [TROPOMI]), NO2can be measuredwith an unprecedented combination of accuracy, spatial coverage, and resolution. This work presents thefirst TROPOMI NO2measurements near the Canadian Oil Sands and shows that these measurements havean outstanding ability to detect NO2on a very high horizontal resolution that is unprecedented for satelliteNO2observations. Further, these satellite measurements are in excellent agreement with aircraft andground-based measurements.

TROPOspheric Monitoring Instrument (TROPOMI)

The NASA Applied Sciences Program: Volcanic Ash Observations and Applications

Since 2000, the NASA Applied Sciences Program has been actively transitioning observations and research to operations. Particular success has been achieved in developing applications for NASA Earth Observing Satellite (EOS) sensors, integrated observing systems, and operational models for volcanic ash detection, characterization, and transport. While other observations and applications remain to be developed for the current constellation of NASA EOS sensors and integrated with observing and forecast systems, future requirements and capabilities for volcanic ash observations and applications are also being developed. Many of these are based on technologies currently being tested on NASA aircraft, Unmanned Aerial Systems (UAS) and balloons. All of these efforts and the potential advances that will be realized by integrating them are shared in this presentation.

John J Murray

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

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

TROPOMI

Decoupling the Effects of Anthropogenic Emission Reductions from the Meteorology and Natural Emissions in TROPOMI NO 2 Retrievals During the 2020 COVID-19 Lockdowns

Satellite measurements during the COVID-19 lockdowns that began in 2020 revealed unprecedented reductions in NO 2 tropospheric vertical column densities (VCD). These reductions have largely been attributed to reduced anthropogenic emissions associated with abrupt decreases in road traffic and other power consuming business activities. Although decreased emissions tended to be the main contributor to the observed NO 2 VCD reduction, meteorological variability also played a role. Whereas the observed VCD changes were predominantly negative in places where public health policies were strictly enforced, meteorology had both positive and negative effects over short time intervals. Here, we present results from a global study of the NO 2 reductions aimed at disentangling the meteorological and natural emission variability from the anthropogenic emissions over the world’s most populated megacities. For this study, NASA’s TROPOMI NO 2 algorithm was used in conjunction with the Global Modeling Initiative (GMI) chemical transport model to separate the contributions due emissions and meteorology. A priori NO 2 profiles were generated from two GMI simulations performed for 2020 at a resolution of 0.25° longitude x 0.25° latitude. The first simulation used updated, COVID-impacted NOX emissions based on Forster et al. (2020), while the second simulation used the 2019 emissions with the 2020 meteorology, referred to here as 2020BAU. The 2020BAU data set allowed for the decoupling of the emission component from the meteorology. When compared to the same period in 2019, NO 2 column amounts during the lockdowns were reduced in 35 out of 36 cities. While reduced emissions contributed most to the observed total change in NO 2 during the lockdowns, the effects of meteorology were significant, ranging between 40% (Chennai) and 15% (Beijing). In China, an increase in NO 2 levels due to meteorology were observed in five out of the seven cities considered in the study. Use of different a-priori NO 2 profiles from the two simulations in our TROPOMI retrievals allowed for a determination of retrieval errors that ranged from -1.7% to -11.0%. We used the quality assurance flag > 0.75 to select the highest quality scenes in the study period. Using the GMI, we estimated the sampling biases to be in the range -11% to 10% of the total change for the cities in our study.

Brad Fisher