Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transport models”

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 811 records · Page 45

Evaluation of several relevant fractionation processes as possible explanation for radioxenon isotopic activity ratios in samples taken near underground nuclear explosions in shafts and tunnels

Gas samples taken from two historic underground nuclear tests done in 1989 at the Nevada National Security Site (NNSS), formerly the Nevada Test Site (NTS), were examined to determine how xenon isotopes fractionate because of early-time cavity processes, transport through the rock, or dispersal through tunnels. Xenon isotopes are currently being used to distinguish civilian sources of xenon in the atmosphere from sources associated with underground nuclear explosions (UNEs). The two nuclear tests included (1) BARNWELL, a test conducted in a vertical shaft approximately 600 m below ground surface at Pahute Mesa, and (2) DISKO ELM, a horizontal line-of-sight test done in P-tunnel approximately 261 m below the surface of Aqueduct Mesa. Numerical flow and transport models developed for the two sites had mixed success when attempting to match the observed xenon isotope ratios. At the BARNWELL site, the simulated xenon isotope ratios were consistent with measurements from the chimney and ground surface and appeared to have been affected primarily by fractionation during subsurface transport. At the DISKO ELM site, samples taken from two elevations in the chimney failed to show the degree of fractionation predicted by the models during transport and did not show evidence for significant fractionation due to early-time condensation of refractory xenon-precursor radionuclides into the melt glass. Gas samples taken from the adjacent tunnels in the days following the test showed mixed evidence for early-time separation of xenon isotopes from their iodine precursors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Towards an A-train Aerosol Assimilation System: Assimilation of MODIS Aerosol Optical Thickness Retrievals into a Global Aerosol Transport and Radiation Model

To understand the role of atmospheric aerosols in Earth's climate system we ideally require a global, three-dimensional, and temporally continuous system of measurements. The view of aerosols provided by a single polar-orbiting satellite instrument (e.g., MODIS), while providing the greatest spatial coverage of any measurement system available, can offer only a single relevant observation per day of a particular location. Furthermore, the quantitative interpretation of those observations is subject to considerable uncertainty because of we do not know the aerosol optical properties or vertical profiles well enough. Ground-based and in situ measurements can provide missing detail about aerosol vertical and optical properties, but lack the spatial (and sometimes temporal) coverage desired. Chemical transport models treating the aerosol lifecycle complement this picture by providing three-dimensional, global, and continuous aerosol fields, but suffer from poor constraints on their initial conditions (e.g., the aerosol source, sink, and transformation mechanisms). Here we present a global, three- dimensional aerosol transport model and assimilation system. The model assimilates the MODIS aerosol optical thickness product to constrain its aerosol distributions and burdens. We provide an analysis of the model aerosol distributions as constrained by ground-based sun photometer measurements and aerosol products from other satellite measurements.

Colarco, P.↗

The right conditions for high-precision dynamic temperature and heat capacity measurement via pyrometry and conductivity

The pursuit of accurate bulk temperature T under extreme conditions has been a long-standing goal of the high pressure science community, complicated by a lack of data to inform models. To reach these extremely high-pressure, high-temperature (high P − T) conditions, a combination of dynamic and heated static experiments (e.g., diamond or gem anvil cel experiments) are used. For example, in a diamond anvil cell (DAC) experiment, a sample placed in the DAC is first pressurized. Following pressurization, the sample T is increased either by heating the entire DAC (usually using resistive heating, and limited to ∼1000K) or by applying intense laser power to the sample surfaces. In a dynamic experiment, the process of pressurizing the sample also heats it. In the case of shock physics experiments, such heating is substantial, easily reaching thousands of Kelvin; in our work we have seen T ∼17000K. Most methods of measuring temperature at ambient are not compatible with experiments under these high-pressure, high-temperature conditions: thermocouples break, melt, or have conductivity properties that differ from ambient where they are calibrated; thermometers would melt; both are too slow. As a result most methods are based on non-contact techniques such as x-ray diffraction broadening, neutron scattering, or optical methods. Of these, optical methods using the visible and near-infrared region of the spectrum are the most commonly used as the sources and detectors are readily available. In the case of optical methods the optical depth, and therefore the measurement location, is limited to the surface. When a window or anvil material is used, heat flows from the sample into the window/anvil. Likewise, if the sample undergoes a change in thermodynamic state, such as expansion upon release, different T may be expected. As a result, the surface or apparent temperature T app measurement will differ from the bulk or interior temperature that is desired. This surface measurement must be related to the bulk measurement using thermal transport models and material models. While it is tempting to conclude that one should just use x-ray methods that directly probe the interior, even these methods have been shown to depend on thermal transport and material models. Regardless of the method used to create the high P − T condition, therefore, we must understand the role of thermal transport and material models upon our interpretation of the T measurement, as well as the errors and uncertainties associated with the choice of models used in the analysis. This is a substantial area of research and this paper is by no means a complete survey of the relevant sources of uncertainty. For example, we have yet to begin to address alternate transport models in a detailed manner (e.g., Tan-Ahrens), or the many models that use additional layers to approximate melting, turbulence, or epitaxial phenomena). Likewise, we have not explored the impact upon uncertainty of thermal models that use temperature-dependent thermal transport coefficients, or the wide range of material models that can be applied. Instead, this paper focuses on using one simple model, the Urtiew-Grover model, to understand the sources of error in T measurement so that we may identify how best to focus future research efforts to return the best improvements and avoid working on over-optimizing a single type of measurement. To this end, we work through some of the best and worst case scenarios for T measurement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Isotopic Dependence of GCR Fluence behind Shielding

In this paper we consider the effects of the isotopic composition of the primary galactic cosmic rays (GCR), nuclear fragmentation cross-sections, and isotopic-grid on the solution to transport models used for shielding studies. Satellite measurements are used to describe the isotopic composition of the GCR. For the nuclear interaction data-base and transport solution, we use the quantum multiple-scattering theory of nuclear fragmentation (QMSFRG) and high-charge and energy (HZETRN) transport code, respectively. The QMSFRG model is shown to accurately describe existing fragmentation data including proper description of the odd-even effects as function of the iso-spin dependence on the projectile nucleus. The principle finding of this study is that large errors (+/-100%) will occur in the mass-fluence spectra when comparing transport models that use a complete isotopic-grid (approx.170 ions) to ones that use a reduced isotopic-grid, for example the 59 ion-grid used in the HZETRN code in the past, however less significant errors (<+/-20%) occur in the elemental-fluence spectra. Because a complete isotopic-grid is readily handled on small computer workstations and is needed for several applications studying GCR propagation and scattering, it is recommended that they be used for future GCR studies.

Cucinotta, Francis A.↗

Model simulations of Plum Island Ecosystems LTER low marsh site using ELM-PFLOTRAN

Model simulations using the E3SM Land Model (ELM) coupled to the PFLOTRAN reactive transport model via the Alquimia interface. The simulations were conducted for a tidal salt marsh at the Plum Island Ecosystems LTER near Rowley, Massachusetts, USA. Model simulations were forced using site-specific tidal cycles and salinity, and the simulations used a biogeochemical reaction network including aerobic decomposition, sulfate reduction, iron reduction, and methanogenesis. Model outputs include simulated carbon stocks, carbon dioxide and methane fluxes, and porewater concentrations of key solutes related to sulfur, iron, and carbon cycling. The model simulations included a saline simulation (with tidal sulfate inputs), a fresh simulation (with low salinity and low sulfate inputs), and a saline simulation with lower vegetation productivity to represent the effect of salinity on vegetation. These simulations were conducted to demonstrate that a new model framework incorporating subsurface redox and biogeochemical interactions into a land surface model could reproduce measured surface greenhouse gas fluxes and biogeochemical dynamics in tidal marsh ecosystems, and to test whether including redox interactions in a land surface model would allow the model to resolve contrasts in biogeochemical cycling and greenhouse gas production between saline and freshwater wetlands.The data package includes gzipped tar archives (which can be expanded using standard tar and gzip utilities) of model outputs from three model configurations: saline subsurface and reduced vegetation productivity related to salinity; saline subsurface with vegetation productivity not reduced; and freshwater. Also included are code for the modified E3SM model, Alquimia interface, and PFLOTRAN reactive transport simulator in gzipped tar format; plain text parameter and configuration files; python code files for visualizing model output and defining model configurations; and model output, tide and salinity forcing, and configuration files in netCDF format. See the README.md file in the data package for a detailed description of all files contained in the package. All files are in netCDF (.nc), gzipped tar archive (.tar.gz or .tgz), or text (all other files).Updated: May 13, 2024. Model output, E3SM code, PFLOTRAN input files, and python codes for visualizing results were updated to reflect changes made for the manuscript revision. The updated archive reflects the code and model output from the final accepted manuscript. Changes included updated reaction parameters reflecting improved parameterization and additional comparisons with field measurements. E3SM code changes included better support for multiple grid cells and improved flow and transport parameterization.

54 ENVIRONMENTAL SCIENCES↗

Interpreting Space-Based Trends in Carbon Monoxide with Multiple Models

We use a series of chemical transport model and chemistry climate model simulations to investigate the observed negative trends in MOPITT CO over several regions of the world, and to examine the consistency of timedependent emission inventories with observations. We find that simulations driven by the MACCity inventory, used for the Chemistry Climate Modeling Initiative (CCMI), reproduce the negative trends in the CO column observed by MOPITT for 2000-2010 over the eastern United States and Europe. However, the simulations have positive trends over eastern China, in contrast to the negative trends observed by MOPITT. The model bias in CO, after applying MOPITT averaging kernels, contributes to the model-observation discrepancy in the trend over eastern China. This demonstrates that biases in a model's average concentrations can influence the interpretation of the temporal trend compared to satellite observations. The total ozone column plays a role in determining the simulated tropospheric CO trends. A large positive anomaly in the simulated total ozone column in 2010 leads to a negative anomaly in OH and hence a positive anomaly in CO, contributing to the positive trend in simulated CO. These results demonstrate that accurately simulating variability in the ozone column is important for simulating and interpreting trends in CO.

carbon monoxide↗

Stratospheric Water Vapor and the Asian Monsoon: An Adjoint Model Investigation

A new adjoint model of the Goddard Parameterized Chemistry and Transport Model is used to investigate the role that the Asian monsoon plays in transporting water to the stratosphere. The adjoint model provides a unique perspective compared to non-diffusive and non-mixing Lagrangian trajectory analysis. The quantity of water vapor transported from the monsoon and the pathways into the stratosphere are examined. The emphasis is on the amount of water originating from the monsoon that contributes to the tropical tape recorder signal. The cross-tropopause flux of water from the monsoon to the midlatitude lower stratosphere will also be discussed.

Olsen, Mark A.↗

Attribution of Trends and Variability in Surface Ozone over the United States

Concentrations of tropospheric ozone, a greenhouse gas and air pollutant, are impacted by changes in precursor emissions as well meteorology and influx from the stratosphere. Observations show a decreasing trend in summertime surface ozone at rural stations in the eastern United States, while some western stations show increasing trends, particularly in springtime. We use the Global Modeling Initiative (GMI) global chemical transport model to investigate the roles of precursor emission changes, meteorological variability, and stratosphere-troposphere exchange (STE) in explaining observed trends in surface ozone from rural sites in the United States from 1991-2010. The model's interannual variability shows significant correlations with observations from many of the surface sites. We also compare the simulated ozone to ozonesonde data for several locations with sufficiently long records. We compare a simulation with time-dependent precursor emissions, including emission reductions over the United States and Europe and increases over Asia, to a simulation with fixed emissions to quantify the impact of changing emissions on the surface trends. The simulation with varying emissions reproduces much of the east-west difference in summertime ozone over the U.S., although it generally underestimates the negative trend in the East. In contrast, the fixed-emission simulation shows increasing ozone at both eastern and western sites. We will discuss possible causes of this behavior, including long-range transport and STE.

ozone↗

A three-dimensional model simulation of atmospheric nitrous oxide

The NASA Langley 3D GCM chemical transport model is used to investigate the distribution of atmospheric N2O up to 60 km altitude. The transport characteristics of the model is evaluated without the complications of a detailed chemical formulation for all of the relevant stratospheric minor constituents. Interpretation of the yearly average zonal mean N2O distribution in terms of transport by the yearly averaged meridional circulation and stratospheric photochemical loss indicates large regions in the Northern Hemisphere stratosphere where dynamical mixing apparently plays a large role in maintaining the N2O distribution. In these regions, slopes of the N2O mixing ratio isopleths are maintained by competition between advection by the meridional circulation acting to steepen and dynamical mixing acting to flatten the slopes.

Turner, R. E.↗

Estimating the Contribution of the Nickel to Protium Loading of Full-Length Getters

The full-length getter (FLG) is a critical component of the tritium producing burnable absorber rod (TPBAR), designed to capture tritium produced within the lithium aluminate pellets. However, due to the chemical similarity between protium (¹H) and tritium (³H), the FLG readily absorbs protium, reducing its capacity to absorb tritium and increasing the risk of tritium permeation into the reactor coolant. This study evaluates protium produced from neutron irradiation of nickel plating on the FLG through 5?Ni(n,p) reactions, quantifies its contribution to the total protium observed, and informs models of tritium and hydrogen transport within TPBARs. During irradiation, neutron capture by 58Ni results in the formation of 5?Ni, which undergoes neutron bombardment to produce 4He via 5?Ni(n,a) reactions and protium via 5?Ni(n,p) reactions. The measured helium content post-irradiation provides insight into the neutron capture processes within the getter. The 4He measured within the getter may be useful in determining hydrogen produced by the nickel plating on FLG because 59Ni also produces protium in a 59Ni(n, p) reaction. The contribution of 1H from the nickel in FLG contributed less than 1% of the measured H2 gas in PIE, ranging from 0.005% to 0.614%. These findings refine current knowledge of the protium-tritium interplay in TPBARs and support the development of improved transport models for tritium and hydrogen, ultimately aiding in the optimization of TPBAR design and reactor operations.

Arbova, Dana L.↗

Key advancements toward eliminating the “drive deficit” in ICF hohlraum simulations

This paper, extended from a Letter [Chen et al., Phys. Rev. E 110, L013201 (2024)] published in Phys. Rev. E in 2024, presents detailed results from a set of experiments studying the x-ray drive produced in laser-heated ignition-scale hohlraums at the National Ignition Facility. The results reveal that deficiencies in non-local thermodynamic equilibrium atomics physics calculations of gold emissivity lead to the overestimation of the “M-band” radiation and are responsible for discrepancies in drive predictions vs measurements in these hohlraums. These experiments were performed on ViewFactor targets with the same gas fill and laser pulse shapes as currently used in ongoing ignition experiments. The common hohlraum model with a magnetohydrodynamics heat transport model and a 10%–20% reduction to the M-band absorption and emission opacity can bring the simulations into agreement with several independent measurements, including the absolute radiation flux, the radiation flux in the 2–4 keV spectral range, and spectroscopic measurements of both K-shell dopant and L-shell gold. This new model is applicable to ICF and high-energy density experiments that use gold hohlraums as drivers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hard jet substructure in a multistage approach

We present predictions and postdictions for a wide variety of hard jet-substructure observables using a multistage model within the framework. The details of the multistage model and the various parameter choices are described in []. A novel feature of this model is the presence of two stages of jet modification: a high-virtuality phase [modeled using the modular all twist transverse-scattering elastic-drag and radiation model ()], where modified coherence effects diminish medium-induced radiation, and a lower virtuality phase [modeled using the linear Boltzmann transport model ()], where parton splits are fully resolved by the medium as they endure multiple scattering induced energy loss. Energy-loss calculations are carried out on event-by-event viscous fluid dynamic backgrounds constrained by experimental data. The uniform and consistent descriptions of multiple experimental observables demonstrate the essential role of modified coherence effects and the multistage modeling of jet evolution. Using the best choice of parameters from [], and with no further tuning, we present calculations for the medium modified jet fragmentation function, the groomed jet momentum fraction z g and angular separation r g distributions, as well as the nuclear modification factor of groomed jets. These calculations provide accurate descriptions of published data from experiments at the Large Hadron Collider. Furthermore, we provide predictions from the multistage model for future measurements at the BNL Relativistic Heavy Ion Collider. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Deuterium Permeation Through a Self-Supported Palladium-Silver Membrane in Helium Gas Mixtures

Palladium membranes and membrane reactors can separate and purify tritium from impurities in the plasma exhaust processing system for the fusion energy fuel cycle. Membranes can also act as a continuous separation method to remove tritium from helium streams in the breeder blanket tritium extraction system, such as from the purge gas of solid breeders. To investigate the potential of these membranes for these applications, we performed a deuterium permeation campaign with a self-supported palladium-silver (Pd-25Ag wt%) tube of 0.15 m length, 76 µ m wall thickness, and 3.0 mm inner diameter. A gas mixture of 3.95% D 2 and a balance of He was delivered to the inside of the tube and D permeated radially outwards through the membrane into a vacuum chamber. Permeation experiments were conducted at 300 ° C, 350 ° C, 400° C, and 450°C, from 100 to 1000 sccm total flow rate, and with total pressures of 90, 150, 190, and 250 kPa. Further, deuterium permeation was consistently lower than predicted from diffusion-limited permeation models, thus we developed a transport model that included gas-phase mass transfer and surface reactions to model experimental results. The dissociation constant was optimized to fit the developed model to experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fabrication of model ultrafiltration membranes with uniform, high aspect ratio pores

Here, in this manuscript, we report the facile fabrication of large-area model membranes with highly uniform and high aspect ratio pores with diameters <20 nm. These membranes are useful for fundamental investigations of separation by size exclusion in the ultrafiltration regime, where species to be separated from solution have dimensions of 1–100 nm. Such investigations require membranes with narrow pores and high aspect ratios such that the Hagen–Poiseuille equation is followed, enabling well-known models such as the hindered transport model to be evaluated and other affecting factors to be ignored. We demonstrate that the sub-20 nm pores in the membrane are of sufficiently high aspect ratio such that water flux through the membrane is consistent with the Hagen–Poiseuille equation. The fabrication relies on self-assembling block copolymers to form uniform, densely packed patterns with sub-20 nm resolution, sequential infiltration synthesis to convert the block copolymer in situ into a mask with adequate contrast to etch pores with an aspect ratio >5, and low-resolution photolithography to transfer the pattern over a large area into a silicon nitride membrane. Model membranes with narrow pore-size distribution fabricated in this way provide the means to investigate parameters that impact size-selective ultrafiltration separations such as the relationships between solute or particle size and pore size, their distributions, and rejection profiles, and, therefore, test the validity or limits of separation models.

Chen, Wen↗

Model Data Archive Associated with Manuscript "Fire-altered Carbon Pools Create Disturbance Memory in Stream Dissolved Organic Carbon"

This data package supports the publication “Fire-altered Carbon Pools Create Disturbance Memory in Stream Dissolved Organic Carbon” by Li et al. (2026). The package contains processed model inputs, configuration files, restart files, simulation outputs, scripts, and visualization products used to evaluate post-fire dissolved organic carbon (DOC) dynamics in the Naches River Watershed, Washington, USA, following the 2021 Schneider Springs Fire. The modeling workflow couples ELM-BGC, the biogeochemistry-enabled Energy Exascale Earth System Model Land Model; ATS, the Advanced Terrestrial Simulator for integrated surface-subsurface hydrology; and PFLOTRAN, a reactive transport model for multicomponent aqueous geochemistry. Together, these models simulate how wildfire-induced changes in vegetation, litter, coarse woody debris, and soil organic matter influence DOC production, transport, and reaction from burned hillslopes to stream networks. The archive includes preprocessed meteorological, geospatial, hydrologic, and biogeochemical forcing data; ELM-BGC-derived DOC source terms; ATS mesh files; PFLOTRAN reactive-transport inputs; model configuration files; spin-up and transient restart files; watershed-scale diagnostic outputs; stream concentration time series; and figures or visualization files used to inspect and reproduce key results. File types include Hierarchical Data Format 5 (HDF5) files for gridded forcing and model-coupling data, model input and configuration files for ELM-BGC, ATS, and PFLOTRAN, restart and simulation-output files generated by the modeling workflow, tabular or time-series diagnostic outputs, scripts for post-processing and figure generation, and image or visualization products associated with the manuscript. Use of the package depends on the intended task. Re-running the simulations requires the relevant modeling software, including ELM-BGC, ATS, and PFLOTRAN as ATS's geochemical engine. Inspecting outputs and reproducing figures requires Python with scientific plotting libraries such as Matplotlib, and three-dimensional model outputs may be viewed with ParaView. Geographic information system files or maps may be inspected with ArcGIS Pro or comparable GIS software. The data package is intended to enable traceability, reuse, and partial reproduction of the coupled land-to-watershed hydro-biogeochemical modeling workflow used to test how wildfire disturbance affects terrestrial carbon pools and downstream DOC dynamics.

ATS↗

Multiscale Development and Validation of the Stainless Steel Alloy Corrosion (SStAC) Tool for High Temperature Engine Materials

The goal of this project was to create a simulation tool to assist in part design, reducing costs by eliminating the need for conservative material selection, and enabling alloy optimization to improve corrosion resistance. Thus, it helps decrease costs for existing engines and assists in the material and engine design for the engines of the future. To meet our stated goal, we developed the Stainless-Steel Alloy Corrosion (SStAC) tool, which: (1) Models corrosion of engine valves fabricated from 21-2N valve steel in an engine environment at temperatures up to 800 °C. (2) Runs in 1D for a fast estimate of the corrosion rate, and in 2D or 3D for more detailed simulations that represent any part geometry and predict the precise location and rate of corrosion and its impact on the mechanical and thermal behavior. (3) Was implemented using the open-source Multiphysics Object Oriented Simulation Environment (MOOSE). (4) Couples the corrosion model with mechanics and thermal transport models. The SStAC tool can predict the corrosion rate of 21-2N engine valves with less than 10% error.

36 MATERIALS SCIENCE↗

A practical concept for powered or tethered weight-lifting LTA vehicles

A concept for a multi-hull weightlifting airship is presented. The concept is based upon experience in the design and handling of gas-filled balloons for commercial purposes, it was first tested in April, 1972. In the flight test, two barrage balloons were joined side-by-side, with an intermediate frame, and launched in captive flight. The success of this flight test led to plans for a development program calling for a powered, piloted prototype, a follow-on 40 ton model, and a 400 ton transport model. All of these airships utilize a tetrehedric three-line tethering method for loading and unloading phases of flight, which bypasses many of the difficulties inherent in the handling of a conventional airship near the ground. Both initial and operating costs per ton of lift capability are significantly less for the subject design than for either helicopters or airships of conventional mono-hull design.

Balleyguier, M. A.↗

Evaluation of Space-Based Constraints on Global Nitrogen Oxide Emissions with Regional Aircraft Measurements over and Downwind of Eastern North America

We retrieve tropospheric nitrogen dioxide (NO 2) columns for May 2004 to April 2005 from the SCIAMACHY satellite instrument to derive top-down emissions of nitrogen oxides (NO(x) = NO + NO2) via inverse modeling with a global chemical transport model (GEOS-Chem). Simulated NO 2 vertical profiles used in the retrieval are evaluated with airborne measurements over and downwind of North America (ICARTT); a northern midlatitude lightning source of 1.6 Tg N/yr minimizes bias in the retrieval. Retrieved NO2 columns are validated (r2 = 0.60, slope = 0.82) with coincident airborne in situ measurements. The top-down emissions are combined with a priori information from a bottom-up emission inventory with error weighting to achieve an improved a posteriori estimate of the global distribution of surface NOx emissions. Our a posteriori NOx emission inventory for land surface NOx emissions (46.1 Tg N/yr) is 22% larger than the GEIA-based a priori bottom-up inventory for 1998, a difference that reflects rising anthropogenic emissions, especially from East Asia A posteriori NOx emissions for East Asia (9.8 Tg N/yr) exceed those from other continents. The a posteriori inventory improves the GEOS-Chem simulation of NOx, peroxyacetylnitrate, and nitric acid with respect to airborne in situ measurements over and downwind of New York City. The a posteriori is 7% larger than the EDGAR 3.2FT2000 global inventory, 3% larger than the NEI99 inventory for the United States, and 68% larger than a regional inventory for 2000 for eastern Asia. SCIAMACHY NO2 columns over the North Atlantic show a weak plume from lightning NO(x).

Martin, Randall V.↗