Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model species”

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 721 records · Page 40

Interpreting Methanol v(sub 2)-Band Emission in Comets Using Empirical Fluorescence g-Factors

For many years we have been developing the ability, through high-resolution spectroscopy targeting ro-vibrational emission in the approximately 3 - 5 micrometer region, to quantify a suite of (approximately 10) parent volatiles in comets using quantum mechanical fluorescence models. Our efforts are ongoing and our latest includes methanol (CH3OH). This is unique among traditionally targeted species in having lacked sufficiently robust models for its symmetric (v(sub 3) band) and asymmetric (v(sub 2) and v(sub 9) bands) C-H3 stretching modes, required to provide accurate predicted intensities for individual spectral lines and hence rotational temperatures and production rates. This has provided the driver for undertaking a detailed empirical study of line intensities, and has led to substantial progress regarding our ability to interpret CH3OH in comets. The present study concentrates on the spectral region from approximately 2970 - 3010 per centimeter (3.367 - 3.322 micrometer), which is dominated by emission in the (v(sub 7) band of C2H6 and the v(sub 2) band of CH3OH, with minor contributions from CH3OH (v(sub 9) band), CH4 (v(sub 3)), and OH prompt emissions (v(sub 1) and v(sub 2)- v(sub 1)). Based on laboratory jet-cooled spectra (at a rotational temperature near 20 K)[1], we incorporated approximately 100 lines of the CH3OH v(sub 2) band, having known frequencies and lower state rotational energies, into our model. Line intensities were determined through comparison with several comets we observed with NIRSPEC at Keck 2, after removal of continuum and additional molecular emissions and correcting for atmospheric extinction. In addition to the above spectral region, NIRSPEC allows simultaneous sampling of the CH3OH v(sub 3) band (centered at 2844 per centimeter, or 3.516 micrometers and several hot bands of H2O in the approximately 2.85 - 2.9 micrometer region, at a nominal spectral resolving power of approximately 25,000 [2]. Empirical g-factors for v(sub 2) lines were based on the production rate as determined from the v(sub 3) Q-branch intensity; application to comets spanning a range of rotational temperatures (approximately 50 - 90 K) will be reported. This work represents an extension of that presented for comet 21P/Giacobini-Zinner at the 2010 Division for Planetary Sciences meeting [3]. Our empirical study also allows for quantifying CH3OH in comets using IR spectrometers for which the v(sub 3) and v(sub 2) bands are not sampled simultaneously, for example CSHELL/NASA IRTF or CRIRES/VLT.

DiSanti, Michael↗

Single-nuclei transcriptome analysis of the shoot apex vascular system differentiation in Populus

ABSTRACT Differentiation of stem cells in the plant apex gives rise to aerial tissues and organs. Presently, we lack a lineage map of the shoot apex cells in woody perennials – a crucial gap considering their role in determining primary and secondary growth. Here, we used single-nuclei RNA-sequencing to determine cell type-specific transcriptomes of the Populus vegetative shoot apex. We identified highly heterogeneous cell populations clustered into seven broad groups represented by 18 transcriptionally distinct cell clusters. Next, we established the developmental trajectories of the epidermis, leaf mesophyll and vascular tissue. Motivated by the high similarities between Populus and Arabidopsis cell population in the vegetative apex, we applied a pipeline for interspecific single-cell gene expression data integration. We contrasted the developmental trajectories of primary phloem and xylem formation in both species, establishing the first comparison of vascular development between a model annual herbaceous and a woody perennial plant species. Our results offer a valuable resource for investigating the principles underlying cell division and differentiation conserved between herbaceous and perennial species while also allowing us to examine species-specific differences at single-cell resolution.

59 BASIC BIOLOGICAL SCIENCES↗

A robust and efficient generalized real cascade model

This paper describes a new separation cascade model that allows for arbitrary cascade connectivity and includes the effects of mixing, material losses and generalized separation element performance. The model also accommodates mixtures with an arbitrary number of component species. The mathematical and numerical framework of this model includes both time-dependent and steady state cascade simulation. Robust numerical approximations of the continuous mathematical model are accomplished using the Sundials suite of nonlinear algebraic and differential-algebraic equation solvers. Verification of the model is accomplished using manufactured solutions while validation is demonstrated by comparison to published experimental results and benchmarks.

36 MATERIALS SCIENCE↗

Response of the global climate to changes in atmospheric chemical composition due to fossil fuel burning

Tropospheric ozone and methane might increase in the future as the result of increasing anthropogenic emissions of CO, NOx and CH4 due to fossil fuel burning. Since O3 and CH4 are both greenhouse gases, increases in their concentrations could augment global warming due to larger future amounts of atmospheric CO2. To test this possible climatic impact, a zonal energy-balance climate model has been combined with a vertically-averaged tropospheric chemical model. The latter model includes all relevant chemical reactions which affect species derived from H2O, O2, CH4 and NOx. The climate model correspondingly incorporates changes in the infrared heating of the surface-troposphere system resulting from chemically induced changes in tropospheric ozone and methane. This coupled climate-chemical model indicates that global climate is sensitive to changes in emissions of CO, NOx and CH4, and that future increases in these emissions could enhance global warming due to increasing atmospheric CO2.

Cess, R. D.↗

Investigation of the significance of excited species in wave-heated dense non-ideal microplasmas using particle-in-cell Monte Carlo-collision modeling

A computational model for the interaction of a high-pressure microplasma with an electromagnetic wave is presented. A one-dimensional particle-in-cell Monte Carlo collision model is used to investigate the plasma non-ideality effects in a second-stage laser-heated xenon plasma with a comprehensive chemistry mechanism, including excited species. A xenon microplasma at a temperature of 300 K and pressure of 10 bar becomes non-ideal after the laser heating stage with lower ionization and electronic excitation barriers, which makes the presence of excited species important to the plasma generation process. For these conditions, two-step ionization via excited species becomes a major contributor to plasma ionization. We find that full ionization of the plasma is obtained in ~2 ps when excited species are included in the chemistry model, as opposed to ~2.5 ps for a chemistry mechanism that does not consider their effect. With excited species, the ionization mode transitions from direct ionization via electron-neutral collisions to direct and stepwise ionization as the plasma generation progresses.

36 MATERIALS SCIENCE↗

CORPSE model with litter decomposition parameters derived from the LIDET dataset

This is a version of the CORPSE model (Carbon, Organisms, Rhizosphere and Protection in the Soil Environment, Sulman et al. 2014) that uses litter decomposition parameters derived from a modified Monte Carlo simulation using the LIDET litter decomposition dataset (Long-term Intersite Decomposition Experiment Team, Harmon 2013). The code also includes the Baseline parameters, and the eight other best parameter sets identified in a modified Monte Carlo simulation. Related publication:Juice, S.M., Ridgeway, J.R., Hartman, M.D., Parton, W.J., Berardi, D.M., Sulman, B.N., Allen, K.E., & Brzostek, E.R. Reparameterizing litter decomposition using a simplified Monte Carlo method improves litter decay simulated by a microbial model and alters bioenergy soil carbon estimates. Description of files:The folder "Input Files" contains one folder for each LIDET site with data necessary to run the model. Note that "(site)" in the filenames below indicates where the LIDET site code appears (see Table 1 for site codes). Data streams include: CORPSE_full_spinup_litter.csv, CORPSE_full_spinup_rhizo.csv, CORPSE_full_spinup_bulk.csv, litterbag_init_100g_6spp.csv: initial C and N (kg C or N/m2) pool values for each soil layer, the litterbag_init_100_6spp.csv file is for the litterbag layer and is the same file for all sites. All initial C and N files have the same columns (Column - Description - Units) uFastC - Unprotected fast decomposing carbon - kg carbon/m2 uSlowC - Unprotected slow decomposing carbon - kg carbon/m2 uNecroC - Unprotected necromass carbon - kg carbon/m2 pFastC - Protected fast decomposing carbon - kg carbon/m2 pSlowC - Protected slow decomposing carbon - kg carbon/m2 pNecroC - Protected necromass carbon - kg carbon/m2 livingMicrobeC - Carbon in living microbial biomass - kg carbon/m2 uFastN - Unprotected fast decomposing nitrogen - kg nitrogen/m2 uSlowN - Unprotected slow decomposing nitrogen - kg nitrogen/m2 uNecroN - Unprotected necromass nitrogen - kg nitrogen/m2 pFastN - Protected fast decomposing nitrogen - kg nitrogen/m2 pSlowN - Protected slow decomposing nitrogen - kg nitrogen/m2 pNecroN - Protected necromass nitrogen - kg nitrogen/m2 inorganicN - Inorganic nitrogen - kg nitrogen/m2 CO2 - Carbon in carbon dioxide - kg carbon/m2 livingMicrobeN - Nitrogen in living microbial biomass - kg nitrogen/m2 soilT (site) DOY274start.csv: Average daily soil temperature (oC) interpolated from previously calculated monthly values used in DayCent LIDET simulations (Bonan et al., 2013). soilT (site) DOY274start.csv: Average daily soil volumetric water content (VWC) scalar interpolated from previously calculated monthly values used in DayCent LIDET simulations (Bonan et al., 2013). litter production.csv: Average daily litter production values for each site, data sources listed in Table S3 of related publication. litter (site) CN.csv: C:N ratio for each species from LIDET dataset (Table 2, Harmon 2013). (site).csv: Table indicating number of observations for each species decomposed at each site. Instructions: Save the model code ("CORPSE_LIDET.R") and "Input Files" folder in the same folder. Also make a folder for the model output (e.g., "results_Baseline") in the same folder. Set the working directory (setwd) in the model code to the folder with the files saved in step #1. Select the parameter set to use for the litter and litterbag compartments, comment out all other parameter sets. Run code. Output will be saved in the folder made in step 1. Output destination can be changed as necessary in code section called "Running the model." Table 1 LIDET sites and site codes used in model files. Site Code - Site AND - H.J. Andrews Experimental Forest BNZ - Bonanza Creek Experimental Forest BSF - Blodgett Research Forest CDR - Cedar Creek Natural History Area CPR - Central Plains Experimental Range HBR - Hubbard Brook Experimental Forest HFR - Harvard Forest JUN - Juneau KBS - Kellogg Biological Station KNZ - Konza Prairie Research Natural Area NWT - Niwot Ridge/Green Lakes Valley OLY - Olympic National Park OLY Conifer forest SEV - Sevilleta National Wildlife Refuge SMR - Santa Margarita Ecological Reserve UFL - University of Florida VCR - Virginia Coast Reserve Table 2 LIDET species and species codes used in model files (6 common species). Species - Species Code Sugar maple (Acer saccharum) - ACSA Drypetes (Drypetes glauca) - DRGL Red pine (Pinus resinosa) - PIRE Chestnut oak (Quercus prinus) - QUPR Western redcedar (Thuja plicata) - THPL Wheat (Triticum aestivum) - TRAE References:Bonan, G. B., Hartman, M. D., Parton, W. J., & Wieder, W. R. (2013). Evaluating litter decomposition in earth system models with long-term litterbag experiments: an example using the Community Land Model version 4 (CLM4). Global Change Biology, 19(3), 957-974. https://doi.org/https://doi.org/10.1111/gcb.12031 Harmon, M. (2013). LTER Intersite Fine Litter Decomposition Experiment (LIDET), 1990 to 2002. Long-Term Ecological Research. Forest Science Data Bank, Corvallis, OR. [Data set]. Accessed http://andlter.forestry.oregonstate.edu/data/abstract.aspx?dbcode=TD023. https://doi.org/10.6073/pasta/f35f56bea52d78b6a1ecf1952b4889c5. Sulman, B. N., Phillips, R. P., Oishi, A. C., Shevliakova, E., & Pacala, S. W. (2014). Microbe-driven turnover offsets mineral-mediated storage of soil carbon under elevated CO2. Nature Climate Change, 4, 1099 - 1102. https://doi.org/10.1038/nclimate2436

Juice, Stephanie↗

Polycyclic aromatic hydrocarbons and the unidentified infrared emission bands - Auto exhaust along the Milky Way

The unidentified infrared emission features (UIR bands) are attributed to a collection of partially hydrogenated, positively charged polycyclic aromatic hydrocarbons (PAHs). This assignment is based on a spectroscopic analysis of the UIR bands. Comparison of the observed interstellar 6.2 and 7.7-micron bands with the laboratory measured Raman spectrum of a collection of carbon-based particulates (auto exhaust) shows a very good agreement, supporting this identification. The infrared emission is due to relaxation from highly vibrationally and electronically excited states. The excitation is probably caused by UV photon absorption. The infrared fluorescence of one particular, highly vibrationally excited PAH (chrysene) is modeled. In this analysis the species is treated as a molecule rather than bulk material and the non-thermodynamic equilibrium nature of the emission is fully taken into account. From a comparison of the observed ratio of the 3.3 to 11.3-micron UIR bands with the model calculations, the average number of carbon atoms per molecule is estimated to be about 20. The abundance of interstellar PAHs is calculated to be about 2 x 10 to the -7th with respect to hydrogen.

Allamandola, L. J.↗

Laboratory measurements and methane photochemistry modeling

Methane is photolyzed by the solar UV in the stratosphere of Saturn. Subsequent photochemistry leads to the production of acetylene (C2H2) and diacetylene (C4H2). These species are produced where it is relatively warm (T is greater than or equal to 140 K), but the tropopause temperature of Saturn (approximately 80 K) is low enough that these two species may freeze out to their respective ices. Numerical models which include both photochemistry and condensation loss make predictions about the mixing ratios of these species and haze production rates. These models are dependent upon knowing reaction pathways and their associated kinetic reaction rate constants and vapor pressures. How uncertainties in the chemistry and improvements in the vapor pressures affect model predictions for Saturn are discussed.

Romani, P. N.↗

Hydraulic redistribution supplies a major water subsidy and improves water status of understory species in a longleaf pine ecosystem

Hydraulic redistribution (HR) is a common phenomenon in water-limited ecosystems; however, it remains unclear how the volume of water transported via HR compares to other components of the hydrologic budget and how HR influences water availability for understory plant communities. In this study, we investigate the absolute and relative magnitude of HR on a forest water budget and identify potential impacts of this water subsidy to understory plant communities. We scaled tree-level estimates of transpiration and HR of three common tree species naturally occurring in a longleaf pine woodland with plot-level measurements of basal area to determine their magnitude at the stand scale. We trenched plots containing understory vegetation but devoid of mature trees and their connected roots to exclude HR subsidies to understory plant species. We analysed soil water isotopes and assessed leaf water potential (Ψ L ) in trenched and control plots to determine if HR results in mixing of water among soil strata and improves understory plant moisture status. Water inputs from HR were equivalent to >30% of total rainfall for the site during the observation period and ~40% of total tree water uptake, depending on species. A stable isotope mixing model confirmed that soil water within HR-exposed plots was more similar to groundwater, whereas soil water within trenched plots was more similar to precipitation. Exclusion of HR via trenching decreased soil moisture and pre-dawn Ψ L for all understory species. These three lines of evidence suggest that HR from overstory trees redistributes a sizable portion of water from deeper to shallower soil profiles and that this water subsidy enhances understory plant water status.

54 ENVIRONMENTAL SCIENCES↗

Stabilized bases for high-order, interpolation semi-Lagrangian, element-based tracer transport

In a computational fluid model of the atmosphere, the advective transport of trace species, or tracers, can be computationally expensive. For efficiency, models often use semi-Lagrangian advection methods. High-order interpolation semi-Lagrangian (ISL) methods, in particular, can be extremely efficient, if the problem of property preservation specific to them can be addressed. Atmosphere models often use geometrically and logically nonuniform grids for efficiency and, as a result, element-based discretizations. Such grids and discretizations make stability a particular problem for ISL methods. Generally, high-order, element-based ISL methods that use the natural polynomial interpolant associated with a nodal finite-element discretization are unstable. Here, we derive new bases having order of accuracy up to nine, with positive nodal weights, that stabilize the element-based ISL method. We use these bases to construct the linear advection operator in the property-preserving Interpolation Semi-Lagrangian Element-based Transport (Islet) method. Then we discuss key software implementation details. Finally, we show performance results for the Energy Exascale Earth System Model's atmosphere dynamical core, comparing the original and new transport methods. These simulations used up to 27,600 Graphical Processing Units (GPU) on the Oak Ridge Leadership Computing Facility's Summit supercomputer.

97 MATHEMATICS AND COMPUTING↗

Real-time capable modeling of ICRF heating on NSTX and WEST via machine learning approaches

Abstract A real-time capable core Ion Cyclotron Range of Frequencies (ICRF) heating model on NSTX and WEST is developed. The model is based on two nonlinear regression algorithms, the random forest ensemble of decision trees and the multilayer perceptron neural network. The algorithms are trained on TORIC ICRF spectrum solver simulations of the expected flat-top operation scenarios in NSTX and WEST assuming Maxwellian plasmas. The surrogate models are shown to successfully capture the multi-species core ICRF power absorption predicted by the original model for the high harmonic fast wave and the ion cyclotron minority heating schemes while reducing the computational time by six orders of magnitude. Although these models can be expanded, the achieved regression scoring, computational efficiency and increased model robustness suggest these strategies can be implemented into integrated modeling frameworks for real-time control applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical Computation of Flame Spread over a Thin Solid in Forced Concurrent Flow with Gas-phase Radiation

Excerpts from a paper describing the numerical examination of concurrent-flow flame spread over a thin solid in purely forced flow with gas-phase radiation are presented. The computational model solves the two-dimensional, elliptic, steady, and laminar conservation equations for mass, momentum, energy, and chemical species. Gas-phase combustion is modeled via a one-step, second order finite rate Arrhenius reaction. Gas-phase radiation considering gray non-scattering medium is solved by a S-N discrete ordinates method. A simplified solid phase treatment assumes a zeroth order pyrolysis relation and includes radiative interaction between the surface and the gas phase.

Jiang, Ching-Biau↗

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↗

Overview of NASA GRCs Green Propellant Infusion Mission Thruster Testing and Plume Diagnostics

The Green Propellant Infusion Mission (GPIM) is sponsored by NASA's Space Technology Mission Directorate (STMD) Technology Demonstration Mission (TDM) office. The goal of GPIM is to advance the technology readiness level of a green propulsion system, specifically, one using the monopropellant, AF-M315E, by demonstrating ground handling, spacecraft processing, and on-orbit operations. One of the risks identified for GPIM is potential contamination of sensitive spacecraft surfaces from the effluents in the plumes of AF-M315E thrusters. NASA Glenn Research Center (GRC) is conducting activities to characterize the effects of AF-M315E plume impingement and deposition. GRC has established individual plume models of the 22-N and 1-N thrusters that will be used on the GPIM spacecraft. The models describe the pressure, temperature, density, Mach number, and species concentration of the AF-M315E thruster exhaust plumes. The models are being used to assess the impingement effects of the AF-M315E thrusters on the GPIM spacecraft. The model simulations will be correlated with plume measurement data from Laboratory and Engineering Model 22-N, AF-M315E thrusters. The thrusters will be tested in a small rocket, altitude facility at NASA GRC. The GRC thruster testing will be conducted at duty cycles representatives of the planned GPIM maneuvers. A suite of laser-based diagnostics, including Raman spectroscopy, Rayleigh spectroscopy, Schlieren imaging, and physical probes will be used to acquire plume measurements of AFM315E thrusters. Plume data will include temperature, velocity, relative density, and species concentration. The plume measurement data will be compared to the corresponding simulations of the plume model. The GRC effort will establish a data set of AF-M315E plume measurements and a plume model that can be used for future AF-M315E applications.

green Propellants↗

Photoionization of seeded combustion products as a method of enhancing the efficiency of magnetohydrodynamic power generators

Here, in this study, we performed an experimental and computational investigation into the feasibility of utilizing photoionization to enhance the electrical conductivity of seeded oxy-fuel combustion products and improve the performance of magnetohydrodynamic (MHD) power generators. We applied a variety of optical and microwave diagnostics to study the ionization and recombination processes of potassium excited by an excimer laser in a high-velocity oxy-fuel free jet. Computational fluid dynamic (CFD) simulations were performed to model the thermophysical properties and species densities of the free jet. The CFD results were validated with position-dependent potassium concentration measurements. Electron recombination exponential lifetimes were measured through time-resolved microwave transmission. The experimental electron lifetimes were compared with lifetimes calculated from CFD-predicted species densities and literature recombination rates. It was determined that K + or O 2 are the most likely recombination partners for photoionized electrons. Time-resolved fluorescence measurements provided evidence of an ionization pathway involving a two-photon ionization of KOH . Finally, a zero-dimensional chemical kinetic model was developed to assess the fundamental viability of inducing a non-equilibrium electron population to provide a net energy return in combustion-driven MHD power generators. We determined that a high energy return is feasible for targeting electrode boundary layers with ultraviolet photoionization. We also found that photoionization could potentially lower the required temperature of the bulk gas flow.

20 FOSSIL-FUELED POWER PLANTS↗

The Influences of Airmass Histories on Radical Species During POLARIS

The POLARIS mission focused on understanding the processes associated with the decrease of polar stratospheric ozone from spring to fall at high latitudes. This decrease is linked primarily to in situ photochemical destruction by reactive nitrogen species, NO and NO2, which also control other catalytic loss cycles. Steady state models have been used to test photochemistry and radical behavior but are not always adequate in simulating radical species observations. In some cases, air mass history can be important and trajectory models give an improved simulation of the radical species. Trajectory chemistry models, however, still consistently underestimate NO and NO2 abundances compared to measurements along the ER-2 flight track. The Goddard chemistry on trajectory model has been used to test updated rate constants for NO2 + OH, NO2 + O and OH + HNO3, key reactions that affect NO and NO2 abundances. We present comparisons between the modified Goddard chemistry on trajectory model, the JPL steady state model and observations from selected flights.

Pierson, James M.↗

GEOS Constituent Data Assimilation Beyond Aura MLS: Assimilating NASA SAGE III/ISS Profiles of Stratospheric Water Vapor and Ozone

Ozone and water vapor in the lower stratosphere are important trace gases for atmospheric chemistry and radiative budget. The Stratospheric Aerosol and Gas Experiment (SAGE) missions have been crucial in monitoring the stratospheric ozone loss and the subsequent recovery as well as the trends in water vapor linked to surface temperature trends. The SAGE III instrument aboard the International Space Station (ISS) continues the SAGE mission record, with high vertical resolution profiles of ozone and water vapor available since mid 2017. The NASA GEOS Earth system model has the new capability to assimilate multi-constituents from ground and space-based instruments using the GEOS Constituent Data Assimilation System (CoDAS). The recently released MERRA-2 Stratospheric Composition Reanalysis with Aura MLS (M2-SCREAM) assimilates version 4.2 MLS ozone, water vapor and other chemically-reactive species with the NASA GEOS model coupled to a stratospheric-only chemistry mechanism and transport constrained to the MERRA-2 reanalysis. While the number of solar occultation observations a day from SAGE III/ISS is about 1% of the total number of profiles observed globally by MLS, the chemical timescales of ozone and water vapor in the lower stratosphere are long enough that the SAGE III/ISS data may provide a useful constraint on the assimilated product. Using the same GEOS CoDAS configuration as M2-SCREAM, we will present a series of experiments to investigate if ozone and water vapor trends are consistent with the assimilation of SAGE observations with and without Aura retrievals, and to determine if the assimilation of SAGE observations produces a steady product for trend analysis, especially as the end of the Aura mission nears. In our experiments, assimilating only SAGE III/ISS water vapor profiles results in water vapor fields more consistent with experiments that assimilate MLS v5; however, in the polar regions SAGE III/ISS observations are not available and the values are unconstrained. We are encouraged by the positive benefit assimilating the less frequent SAGE III/ISS observations has on stratospheric composition. Sensitivity experiments such as these will allow us to assess the added value of SAGE data for continued monitoring of the stratospheric composition for climate and ozone recovery assessments.

SAGE↗

GEOS Constituent Data Assimilation Beyond Aura MLS: Assimilating NASA SAGE III/ISS Profiles of Stratospheric Water Vapor

Ozone and water vapor in the lower stratosphere are important trace gases for atmospheric chemistry and radiative budget. The Stratospheric Aerosol and Gas Experiment (SAGE) missions have been crucial in monitoring the stratospheric ozone loss and the subsequent recovery as well as the trends in water vapor linked to surface temperature trends. The SAGE III instrument aboard the International Space Station (ISS) continues the SAGE mission record, with high vertical resolution profiles of ozone and water vapor available since mid 2017. The NASA GEOS Earth system model has the new capability to assimilate multi-constituents from ground and space-based instruments using the GEOS Constituent Data Assimilation System (CoDAS). The recently released MERRA-2 Stratospheric Composition Reanalysis with Aura MLS (M2-SCREAM) assimilates version 4.2 MLS ozone, water vapor and other chemically-reactive species with the NASA GEOS model coupled to a stratospheric-only chemistry mechanism and transport constrained to the MERRA-2 reanalysis. While the number of solar occultation observations a day from SAGE III/ISS is about 1% of the total number of profiles observed globally by MLS, the chemical timescales of ozone and water vapor in the lower stratosphere are long enough that the SAGE III/ISS data may provide a useful constraint on the assimilated product. Using the same GEOS CoDAS configuration as M2-SCREAM, we will present a series of experiments to investigate if water vapor trends are consistent with the assimilation of SAGE observations with and without Aura retrievals, and to determine if the assimilation of SAGE observations produces a steady product for trend analysis, especially as the end of the Aura mission nears. In our experiments, assimilating only SAGE III/ISS water vapor profiles results in water vapor fields more consistent with experiments that assimilate MLS v5; however, in the polar regions SAGE III/ISS observations are not available and the modelled values are unconstrained. We are encouraged by the positive benefit assimilating the less frequent SAGE III/ISS observations has on stratospheric composition. Sensitivity experiments such as these will allow us to assess the added value of SAGE data for continued monitoring of the stratospheric composition for climate and ozone recovery assessments.

K Emma Knowland↗