SEARCH · Engineering Papers

Results for “species distribution modeling”

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.

SourceResearch connectionFrom yearThrough year

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

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

NASA NTRS1987

Middle atmosphere trace constituents; Proceedings of Workshop 9 of the Twenty-sixth COSPAR Plenary Meeting, Toulouse, France, June 30-July 11, 1986

Papers are presented on a reference model for stratospheric NO2, a proposed reference model for nitric acid, proposed reference models for ozone, and a background stratospheric aerosol reference model. Other topics include a comparison between the observed and calculated distributions of trace species in the middle atmosphere, recent measurements of trace constituents from rocket and balloon probes, NO2 field observations in the stratosphere during the MAP/GLOBUS Campaign in 1983 and 1985, and the vertical distribution of ozone in the equatorial zone. Also considered are precision measurements of stratospheric ozone profiles by rocket-borne optical ozonesondes, vertical profiles of source gases at tropical latitudes, and measurements of stratospheric trace constituent distributions from balloon-borne far infrared observations.

Keating, G. M.↗

NASA NTRS1974

Reacting viscous-shock-layer solutions with multicomponent diffusion and mass injection

Numerical solutions of the viscous-shock-layer equation where the chemistry is treated as being either frozen, equilibrium, or nonequilibrium are presented. Also the effects of the diffusion model, surface catalysis, and mass injection on surface transport and flow parameters are considered. The flow is treated as a mixture of five inert and thermally perfect species. The viscous-shock-layer equations are solved by using an implicit-difference scheme. All calculations are for hyperboloids with included angles of 20 and 45. The flight conditions are those for various altitudes and velocities in the earth's atmosphere. Data are presented to show the effects of the chemical models; diffusion models; surface catalysis; and mass injection of air on heat transfer; skin friction; shock standoff distance; wall pressure distribution; and tangential velocity, temperature, and species profiles. The results show that an equilibrium analysis can substantially overpredict the heat-transfer rates for flow conditions experienced by earth-orbital entry vehicles. Moreover, at such conditions surface catalysis significantly influences heat-transfer and flow-field properties. If a binary rather than a multicomponent diffusion model is assumed, negligible errors in most flow properties result. Quantitative results are presented that show the effect of mass injection on flow properties within and downstream of the injection region.

Moss, J. N.↗

DOE OSTI2024

Integration and validation of some modules for modelling of high-speed chemically reactive flows in two-phase gas-droplet mixtures

Three modules are integrated into the built-in OpenFOAM rhoCentralFoam solver towards accurate and efficient modelling of high-speed chemically reactive flows in two-phase gas-droplet mixtures within the OpenFOAM 10.0 framework. The first module is the mixture-averaged diffusion model. The second module is the built-in OpenFOAM Lagrangian solver coupled with optimised droplet drag coefficient and convective heat transfer coefficient sub-models. The last module is a sparse stiff chemistry solver based on dynamic adaptive hybrid integration (AHI-S). The optimised droplet sub-models are first verified in correct implementation for subsequent simulations in this work. Further, they show good accuracy against experimental and analytical data in the modelling of ammonia droplet acceleration and cooling in the flowing and/or low-temperature air. The accuracy and efficiency gains related to the mixture-averaged diffusion model and the AHI-S chemistry solver are examined by simulating 1-D detonation propagation in ammonia droplet-free/laden ammoniaoxygen mixtures. Numerical results of detonation propagation speed, gaseous temperature, density, and species distributions around the induction zone show good agreement with experimental data and analytical solutions. Compared to the built-in OpenFOAM diffusion model, the mixture-averaged diffusion model provides different numerical predictions of pulsating instabilities in detonation propagation. It shows better accuracy in depicting the detonation structure within the droplet-free section attributed to improved multi-component diffusion modelling. Compared to the built-in OpenFOAM solver EulerImplicit (backward Euler), the AHI-S chemistry solver reduces the computational cost by around 50%. It achieves satisfactory accuracy in calculating detonation propagation speed within the droplet-free section with the optimal efficiency when the safety factor, β, equals 0.5.

42 ENGINEERING↗

NASA NTRS1996

Measurements and Modeling of Soot Formation and Radiation in Microgravity Jet Diffusion Flames

This is a computational and experimental study for soot formation and radiative heat transfer in jet diffusion flames under normal gravity (1-g) and microgravity (0-g) conditions. Instantaneous soot volume fraction maps are measured using a full-field imaging absorption technique developed by the authors. A compact, self-contained drop rig is used for microgravity experiments in the 2.2-second drop tower facility at NASA Lewis Research Center. On modeling, we have coupled flame structure and soot formation models with detailed radiation transfer calculations. Favre-averaged boundary layer equations with a k-e-g turbulence model are used to predict the flow field, and a conserved scalar approach with an assumed Beta-pdf are used to predict gaseous species mole fraction. Scalar transport equations are used to describe soot volume fraction and number density distributions, with formation and oxidation terms modeled by one-step rate equations and thermophoretic effects included. An energy equation is included to couple flame structure and radiation analyses through iterations, neglecting turbulence-radiation interactions. The YIX solution for a finite cylindrical enclosure is used for radiative heat transfer calculations. The spectral absorption coefficient for soot aggregates is calculated from the Rayleigh solution using complex refractive index data from a Drude- Lorentz model. The exponential-wide-band model is used to calculate the spectral absorption coefficient for H20 and C02. It is shown that when compared to results from true spectral integration, the Rosseland mean absorption coefficient can provide reasonably accurate predictions for the type of flames studied. The soot formation model proposed by Moss, Syed, and Stewart seems to produce better fits to experimental data and more physically sound than the simpler model by Khan et al. Predicted soot volume fraction and temperature results agree well with published data for a normal gravity co-flow laminar flames and turbulent jet flames. Predicted soot volume fraction results also agree with our data for 1-g and 0-g laminar jet names as well as 1-g turbulent jet flames.

Ku, Jerry C.↗

NASA NTRS2023

Arctic Tropospheric Ozone: Assessment of Current Knowledge and Model Performance

As the third most important greenhouse gas (GHG) after carbon dioxide (CO2) and methane (CH4), tropospheric ozone (O3) is also an air pollutant causing damage to human health and ecosystems. This study brings together recent research on observations and modeling of tropospheric O3 in the Arctic, a rapidly warming and sensitive environment. At different locations in the Arctic, the observed surface O3 seasonal cycles are quite different. Coastal Arctic locations, for example, have a minimum in the springtime due to O3 depletion events resulting from surface bromine chemistry. In contrast, other Arctic locations have a maximum in the spring. The 12 state-of-the-art models used in this study lack the surface halogen chemistry needed to simulate coastal Arctic surface O3 depletion in the springtime; however, the multi-model median (MMM) has accurate seasonal cycles at non-coastal Arctic locations. There is a large amount of variability among models, which has been previously reported, and we show that there continues to be no convergence among models or improved accuracy in simulating tropospheric O3 and its precursor species. The MMM underestimates Arctic surface O3 by 5 % to 15 % depending on the location. The vertical distribution of tropospheric O3 is studied from recent ozonesonde measurements and the models. The models are highly variable, simulating free-tropospheric O3 within a range of ±50 % depending on the model and the altitude. The MMM performs best, within ±8 % for most locations and seasons. However, nearly all models overestimate O3 near the tropopause (∼300 hPa or ∼8 km), likely due to ongoing issues with underestimating the altitude of the tropopause and excessive downward transport of stratospheric O3 at high latitudes. For example, the MMM is biased high by about 20 % at Eureka. Observed and simulated O3 precursors (CO, NOx, and reservoir PAN) are evaluated throughout the troposphere. Models underestimate wintertime CO everywhere, likely due to a combination of underestimating CO emissions and possibly overestimating OH. Throughout the vertical profile (compared to aircraft measurements), the MMM underestimates both CO and NOx but overestimates PAN. Perhaps as a result of competing deficiencies, the MMM O3 matches the observed O3 reasonably well. Our findings suggest that despite model updates over the last decade, model results are as highly variable as ever and have not increased in accuracy for representing Arctic tropospheric O3.

Arctic↗

NASA NTRS1987

Odd nitrogen during the MAP/GLOBUS 1983 campaign - Theoretical considerations

The paper reviews the chemical processes involved in the formation and destruction of active nitrogen in the stratosphere and indicates that transport processes play an important role in the behavior of these species. Vertical profiles measured during the MAP/GLOBUS 1983 campaign are compared to distributions calculated by a number of one-dimensional models. An intercomparison between six of these models shows that, if some quantitative differences attributed to different model input are noticeable, the same type of behavior is predicted for NO and NO2 by all models. Although the observed concentration of species belonging to the odd nitrogen family should be highly variable, as a result of dynamical disturbances, most data obtained during the campaign appear to be consistent with the chemical scheme currently used in photochemical models.

Brasseur, Guy↗

DOE OSTI2023

Dark Exciton in 2D Hybrid Halide Perovskite Films Revealed by Magneto-Photoluminescence at High Magnetic Field

Here, a comprehensive study of the exciton fine structure (EFS) is presented in 2D-phenethylammonium lead iodide films using magnetic field-induced polarization of photoluminescence (PL) in both Faraday and Voigt configurations at fields up to 25 Tesla. Three exciton bands are identified in the PL spectrum associated with bound, dark, and bright excitons, respectively. Under a high magnetic field in Faraday/Voigt configuration, large field-induced circular/linear polarization is observed in the PL band related to the dark exciton, which is magnetically activated. Furthermore, it is found that the dark exciton has an anomalous field-induced circular polarization, which cannot be explained by the classical Boltzmann distribution of spin-polarized species. These findings are well explained by an effective mass model that includes exchange terms unique to the monoclinic symmetry as a perturbation of the EFS in the approximate tetragonal symmetry. It is also confirmed that the field-induced linear polarization is sensitive to the monoclinic exchange term, whereas the field-induced circular polarization is immune to such term.

14 SOLAR ENERGY↗

NASA NTRS1985

Transport processes in the stratosphere: Model simulations and comparisons with satellite observations

A three dimensional atmospheric model was used to study transport processes and to simulate the distribution of chemically active species in the stratosphere. The present results are part of a long term simulation of the seasonally varying structure of stratospheric trace constituents. The occurrence of a midwinter stratospheric warming and the associated transport of O3 and HNO3 during the simulation are described. Comparison of the simulated distributions of O3 and HNO3 are made with data from the Limb Infrared Monitor of the Stratosphere (LIMS) experiment. In addition, distributions of Ertel's potential vorticity on isentropic surfaces (IPV) are evaluated as a diagnostic for interpreting transport processes. Comparisons are made with IPV distributions inferred from LIMS temperature data.

Grose, W. L.↗

NASA NTRS2004

Remote Sensing of Icy Galilean Moon Surface and Atmospheric Composition Using Low Energy (1 eV-4 keV) Neutral Atom Imaging

We describe a low energy neutral atom imager suitable for composition measurements Europa and other icy Galilean moons in the Jovian magnetosphere. This instrument employs conversion surface technology and is sensitive to either neutrals converted to negative ions, neutrals converted to positive ions and the positive ions themselves depending on the power supply. On a mission such as the Jupiter Icy Moons Orbiter (JIMO), two back-to-back sensors would be flown with separate power supplies fitted to the neutral atom and iodneutral atom sides. This will allow both remote imaging of 1 eV < E < 4 keV neutrals from icy moon surfaces and atmospheres, and in situ measurements of ions at similar energies in the moon ionospheres and Jovian magnetospheric plasma. The instrument provides composition measurements of the neutrals and ions that enter the spectrometer with a mass resolution dependent on the time-of-flight subsystem and capable of resolving molecules. The lower energy neutrals, up to tens of eV, arise from atoms and molecules sputtered off the moon surfaces and out of the moon atmospheres by impacts of more energetic (keV to MeV) ions from the magnetosphere. Direct Simulation Monte Carlo (DSMC) models are used to convert measured neutral abundances to compositional distributions of primary and trace species in the sputtered surfaces and atmospheres. The escaping neutrals can also be detected as ions after photo- or plasma-ionization and pickup. Higher energy, keV neutrals come from charge exchange of magnetospheric ions in the moon atmospheres and provide information on atmospheric structure. At the jovicentric orbits of the icy moons the presence of toroidal gas clouds, as detected at Europa's orbit, provide M e r opportunities to analyze both the composition of neutrals and ions originating from the moon surfaces, and the characteristics of magnetospheric ions interacting with neutral cloud material. Charge exchange of low energy ions near the moons, and directional distributions of the resultant neutrals, allow indirect global mapping of magnetic field structures around the moons. Temporal variation of the magnetic structures can be linked to induced magnetic fields associated with subsurface oceans.

Collier, M. R.↗

NASA NTRS1995

Tracer transport for realistic aircraft emission scenarios calculated using a three-dimensional model

A three-dimensional transport model, which uses winds from a stratospheric data assimilation system, is used to study the transport of supersonic aircraft exhaust in the lower stratosphere. A passive tracer is continuously injected into the transport model. The tracer source distribution is based on realistic scenarios for the daily emission rate of reactive nitrogen species for all forecasted flight routes. Winds are from northern hemisphere winter/spring months for 1979 and 1989; there are minimal differences between the tracer integrations for the 2 years. During the integration, peak tracer mixing ratios in the flight corridors are compared with the zonal mean and found to be greater by a factor of 2 or less. This implies that the zonal mean assumption used in two dimensional models is reasonable during winter and spring. There is a preference for pollutant buildup in the heavily traveled North Pacific and North Atlantic flight corridors. Pollutant concentration in the corridors depends on the position of the Aleutian anticyclone and the northern hemisphere polar vortex edge.

Weaver, Clark J.↗

DOE OSTI2020

Leaf isoprene and monoterpene emission distribution across hyperdominant tree genera in the Amazon basin

Tropical forests are acknowledged to be the largest global source of isoprene (C 5 H 8 ) and monoterpenes (C 10 H 16 ) emissions, with current synthesis studies suggesting few tropical species emit isoprenoids (20-38%) and do so with highly variable emission capacities, including within the same genera. This apparent lack of a clear phylogenetic thread has created difficulties both in linking isoprenoid function with evolution and for the development of accurate biosphere-atmosphere models. In this work, we present a systematic emission study of "hyperdominant" tree species in the Amazon Basin. Across 162 individuals, distributed among 25 botanical families and 113 species, isoprenoid emissions were widespread among both early and late successional species (isoprene: 61.9% of the species; monoterpenes: 15.0%; both isoprene and monoterpenes: 9.7%). The hyperdominant species (69) across the top five most abundant genera, which make up about 50% of all individuals in the Basin, had a similar abundance of isoprenoid emitters (isoprene: 63.8%; monoterpenes: 17.4%; both 11.6%). Among the abundant genera, only Pouteria had a low frequency of isoprene emitting species (15.8% of 19 species). In contrast, Protium, Licania, Inga, and Eschweilera were rich in isoprene emitting species (83.3% of 12 species, 61.1% of 18 species, 100% of 8 species, and 100% of 12 species, respectively). Light response curves of individuals in each of the five genera showed light-dependent, photosynthesis-linked emission rates of isoprene and monoterpenes. Importantly, in every genus, we observed species with light-dependent isoprene emissions together with monoterpenes including β-ocimene. These observations support the emerging view of the evolution of isoprene synthases from β-ocimene synthases. Our results have important implications for understanding isoprenoid function-evolution relationships and the development of more accurate Earth System Models.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Anisotropic angular scattering models of elastic electron-neutral collisions for Monte Carlo plasma simulations

Many laboratory and industrial plasma applications require accurate modeling techniques to understand the interplay between microscopic and macroscopic processes. A prime example of this interplay is how particle and Monte Carlo (MC) simulation codes describe angular scattering of electrons following elastic scattering events. The forward peaked nature of high energy electron elastic scattering is relatively trivial to accurately describe in plasma simulations. However, for lower energy collisions, which produce near isotropic or backward peaked differential cross sections, there is not a strong consensus among the plasma modeling community on how to best describe these angular scattering trends. Here, in this study, we propose a systematic method to approximate the aforementioned non-trivial angular scattering behavior with a formula that can be readily implemented in particle-in-cell (PIC) and/or MC plasma simulation codes. The present approach is specifically applied to fusion relevant atomic hydrogen and helium, as well as for molecular hydrogen, and results are also applicable to the atomic isotopes and homonuclear molecular isotopologues of these species. Comparisons between the present angular distribution function and benchmark scattering data were used to validate the proposed models. In addition, two-term Boltzmann calculations and PIC direct simulation MC simulations revealed that the proposed angular distribution function is accurate, agreeing very well with benchmark convergent close-coupling scattering calculations, and electron transport measurements. These studies confirmed that the present angular distribution function model can be utilized without the need of renormalization to the momentum transfer cross section (as opposed to using the elastic scattering integrated cross section), which has been suggested by several studies in order to correct for deficient angular scattering models, and to agree with transport measurements. Hence, the present anisotropic angular scattering model can be utilized to accurately model the momentum transfer as well as the electron trajectories of elastic collisions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1973

Solutions of the chemical kinetic equations for initially inhomogeneous mixtures.

Following the recent discussions by O'Brien (1971) and Donaldson and Hilst (1972) of the effects of inhomogeneous mixing and turbulent diffusion on simple chemical reaction rates, the present report provides a more extensive analysis of when inhomogeneous mixing has a significant effect on chemical reaction rates. The analysis is then extended to the development of an approximate chemical sub-model which provides much improved predictions of chemical reaction rates over a wide range of inhomogeneities and pathological distributions of the concentrations of the reacting chemical species. In particular, the development of an approximate representation of the third-order correlations of the joint concentration fluctuations permits closure of the chemical sub-model at the level of the second-order moments of these fluctuations and the mean concentrations.

Hilst, G. R.↗

NASA NTRS1994

A new numerical model of the middle atmosphere. 2: Ozone and related species

A new two-dimensional model with detailed photochemistry is presented. The model includes descriptions of planetary wave and gravity wave propagation and dissipation to characterize the wave forcing and associated mixing in the stratosphere and mesosphere. Such a representation allows for explicit calculation of the regions of strong mixing in the middle atmosphere required for accurate simulation of trace gas transport. The new model also includes a detailed description of photochemical processes in the stratosphere and mesosphere. The downward transport of H2, H2O, and NO(y) from the mesosphere to the stratosphere is examined, and it is shown that mesospheric processes can influence the distributions of these chemical species in polar regions. For HNO3 we also find that small concentrations of liquid aerosols above 30 km could play a major role in determining the abundance in polar winter at high latitudes. The model is also used to examine the chemical budget of ozone in the midlatitude stratosphere and to set constraints on the effectiveness of bromine relative to chlorine for ozone loss and the role of the HO2 + BrO reaction. Recent laboratory data used in this modeling study suggest that this process greatly enhances the effectiveness of bromine for ozone destruction, making bromine-catalyzed chemistry second only to HO(x)-catalyzed ozone destruction in the contemporary stratosphere at midlatitudes below about 18 km. The calculated vertical distribution of ozone in the lower stratosphere agrees well with observations, as does the total column ozone during most seasons and latitudes, with the important exception of southern hemisphere winter and spring.

Garcia, Rolando R.↗

DOE OSTI2023

Quantification of hydraulic trait control on plant hydrodynamics and risk of hydraulic failure within a demographic structured vegetation model in a tropical forest (FATES–HYDRO V1.0)

Abstract. Vegetation plays a key role in the global carbon cycle and thus is an important component within Earth system models (ESMs) that project future climate. Many ESMs are adopting methods to resolve plant size and ecosystem disturbance history, using vegetation demographic models. These models make it feasible to conduct more realistic simulation of processes that control vegetation dynamics. Meanwhile, increasing understanding of the processes governing plant water use, and ecosystem responses to drought in particular, has led to the adoption of dynamic plant water transport (i.e., hydrodynamic) schemes within ESMs. However, the extent to which variations in plant hydraulic traits affect both plant water stress and the risk of mortality in trait-diverse tropical forests is understudied. In this study, we report on a sensitivity analysis of an existing hydrodynamic scheme (HYDRO) model that is updated and incorporated into the Functionally Assembled Terrestrial Ecosystem Simulator (FATES) (FATES–HYDRO V1.0). The size- and canopy-structured representation within FATES is able to simulate how plant size and hydraulic traits affect vegetation dynamics and carbon–water fluxes. To better understand this new model system, and its functionality in tropical forest systems in particular, we conducted a global parameter sensitivity analysis at Barro Colorado Island, Panama. We assembled 942 observations of plant hydraulic traits on 306 tropical plant species for stomata, leaves, stems, and roots and determined the best-fit statistical distribution for each trait, which was used in model parameter sampling to assess the parametric sensitivity. We showed that, for simulated leaf water potential and loss of hydraulic conductivity across different plant organs, the four most important traits were associated with xylem conduit taper (buffers increasing hydraulic resistance with tree height), stomatal sensitivity to leaf water potential, maximum stem hydraulic conductivity, and the partitioning of total hydraulic resistance above vs. belowground. Our analysis of individual ensemble members revealed that trees at a high risk of hydraulic failure and potential tree mortality generally have a lower conduit taper, lower maximum xylem conductivity, lower stomatal sensitivity to leaf water potential, and lower resistance to xylem embolism for stem and transporting roots. We expect that our results will provide guidance on future modeling studies using plant hydrodynamic models to predict the forest responses to droughts and future field campaigns that aim to better parameterize plant hydrodynamic models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Microphysical Sensitivity of Superparameterized Precipitation Extremes in the Contiguous United States Due to Feedbacks on Large-Scale Circulation

Superparameterized (SP) global climate models have been shown to better simulate various features of precipitation relative to conventional models, including its diurnal cycle as well as its extremes. While various studies have focused on the effect of differing microphysics parameterizations on precipitation within limited-area cloud-resolving models, we examine here the effect on contiguous U.S. (CONUS) extremes in a global SP model. We vary the number of predicted moments for hydrometeor distributions, the character of the rimed ice species, and the representation of raindrop self-collection and breakup. Using a likelihood ratio test and accounting for the effects of multiple hypothesis testing, we find that there are some regional differences, particularly during spring and summer in the Southwest and the Midwest, in both the current climate and a warmer climate with uniformly increased sea surface temperatures. These differences are most statistically significant and widespread when the number of moments is changed. To determine whether these results are due to (fast) local effects of the different microphysics or the (slower) ensuing feedback on the large-scale atmospheric circulation, we run a series of short, 5-day simulations initialized from reanalysis data. We find that the differences largely disappear in these runs and therefore infer that the different parameterizations impact precipitation extremes indirectly via the large-scale circulation. Finally, we compare the present-day results with hourly rain gauge data and find that SP underestimates extremes relative to observations regardless of which microphysics scheme is used given a fixed model configuration and resolution.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1983

Charged particle distributions in Jupiter's magnetosphere

In situ data from the Pioneer and Voyager spacecraft, supplemented by earth-based observations and theoretical considerations, are used as the basis for the present quantitative, compact model of the 1 eV-several MeV charged particle distribution in the Jovian magnetosphere. The thermal plasma parameters of convection speed, number density, and characteristic energy, are specified as functions of position for electrons and for the ion species H(+), O(+), O(2+), S(+), S(2+), S(3+), and Na(+). Major features of the magnetic field, thermal plasma, and trapped particle distributions, are modeled and results for each plasma region are compared with observed spectra. Comparisons show that the model represents the data to within a factor of 2 + or - 1, except where time variations are significant. Practical applications of the model to spacecraft near Jupiter are given.

Divine, N.↗

DOE OSTI2023

Arctic tropospheric ozone: assessment of current knowledge and model performance

As the third most important greenhouse gas (GHG) after carbon dioxide (CO 2 ) and methane (CH 4 ), tropospheric ozone (O 3 ) is also an air pollutant causing damage to human health and ecosystems. This study brings together recent research on observations and modeling of tropospheric O 3 in the Arctic, a rapidly warming and sensitive environment. At different locations in the Arctic, the observed surface O 3 seasonal cycles are quite different. Coastal Arctic locations, for example, have a minimum in the springtime due to O 3 depletion events resulting from surface bromine chemistry. In contrast, other Arctic locations have a maximum in the spring. The 12 state-of-the-art models used in this study lack the surface halogen chemistry needed to simulate coastal Arctic surface O 3 depletion in the springtime; however, the multi-model median (MMM) has accurate seasonal cycles at non-coastal Arctic locations. There is a large amount of variability among models, which has been previously reported, and we show that there continues to be no convergence among models or improved accuracy in simulating tropospheric O 3 and its precursor species. The MMM underestimates Arctic surface O 3 by 5% to 15% depending on the location. The vertical distribution of tropospheric O 3 is studied from recent ozonesonde measurements and the models. The models are highly variable, simulating free-tropospheric O 3 within a range of ±50% depending on the model and the altitude. The MMM performs best, within ±8% for most locations and seasons. However, nearly all models overestimate O 3 near the tropopause (~300 hPa or ~8 km), likely due to ongoing issues with underestimating the altitude of the tropopause and excessive downward transport of stratospheric O 3 at high latitudes. For example, the MMM is biased high by about 20% at Eureka. Observed and simulated O 3 precursors (CO, NO x , and reservoir PAN) are evaluated throughout the troposphere. Models underestimate wintertime CO everywhere, likely due to a combination of underestimating CO emissions and possibly overestimating OH. Throughout the vertical profile (compared to aircraft measurements), the MMM underestimates both CO and NO x but overestimates PAN. Perhaps as a result of competing deficiencies, the MMM O 3 matches the observed O 3 reasonably well. Our findings suggest that despite model updates over the last decade, model results are as highly variable as ever and have not increased in accuracy for representing Arctic tropospheric O 3 .

54 ENVIRONMENTAL SCIENCES↗