SEARCH · Search NASA
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.
In modeling turbulent combustion, decoupling the chemistry from the turbulence is of great practical significance. In cases in which the equilibrium chemistry model breaks down, laminar flamelet modeling (LFM) is a promising approach to decoupling. Here, the validity of this approach is investigated using direct numerical simulation of a simple chemical reaction in two-dimensional turbulence.
A spectral 3D model of the stratosphere has been used to study the sensitivity of polar ozone with respect to a carbon dioxide increase. The lower stratospheric cooling associated with an imposed CO2 doubling may increase the probability of polar stratospheric cloud (PSC) formation and this affect ozone. The ozone perturbation obtained with the inclusion of a simple parameterization for heterogeneous chemistry on PSCs is compared to that relative to a pure homogeneous chemistry. In both cases the temperature perturbation is determined by a CO2 doubling, while the total chlorine content is kept at the present level. It is shown that the lower temperature may increase the depth and the extension of the ozone hole by extending the area amenable to PSC formation. It may be argued that this effect, coupled with an increasing amount of chlorine, may produce a positive feedback on the ozone destruction.
Throughout the year 2022, the Langley Mobile Ozone Lidar (LMOL) has been frequently collecting ozone and aerosol measurements in the lower troposphere. The lidar is located in the parking lot behind the Atmospheric Sciences building at NASA Langley in Hampton, Virginia, approximately 37.095 N, -76.389 W. The data was collected during a range of several distinct atmospheric conditions, including stratospheric intrusions, surface frontal passages, and long-range transported wildfire smoke plumes. We use this data to evaluate and intercompare the forecast accuracy of two atmospheric chemistry models: the GEOS Composition Forecasting model (GEOS-CF) and the Weather Research and Forecasting model with Chemistry (WRF-Chem). Both models make daily three-dimensional forecasts of many trace gases and species for the contiguous United States. Over a range of ten distinct collection events, each spanning 36 hours to nearly seven days, the forecast models predict tropospheric ozone at NASA Langley with reasonable accuracy. The models best predict the timing and shape of the observed stratospheric intrusions in the middle troposphere but often vary in the magnitude of the ozone mixing ratios. For other atmospheric conditions, there is more variability in the model accuracy. Here, we summarize the accuracy of the models for all events and investigate reasons for differences between the models and the lidar.
Estimates of ground-level ozone concentrations have been improved through data fusion of observations and atmospheric chemistry models. Our previous global ozone estimates for the Global Burden of Disease study corrected for bias uniformly across continents and then corrected near monitoring stations using the Bayesian Maximum Entropy (BME) framework for data fusion. Here, we use the Regionalized Air Quality Model Performance (RAMP) framework to correct model bias over a much larger spatial range than BME can, accounting for the spatial inhomogeneity of bias and nonlinearity as a function of modeled ozone. RAMP bias correction is applied to a composite of 9 global chemistry-climate models, based on the nearest set of monitors. These estimates are then fused with observations using BME, which matches observations at measurement stations, with the influence of observations declining with distance in space and time. We create global ozone maps for each year from 1990 to 2017 at fine spatial resolution. RAMP is shown to create unrealistic discontinuities due to the spatial clustering of ozone monitors, which we overcome by applying a weighting for RAMP based on the number of monitors nearby. Incorporating RAMP before BME has little effect on model performance near stations, but strongly increases R 2 by 0.15 at locations farther from stations, shown through a checkerboard cross-validation. Corrections to estimates differ based on location in space and time, confirming heterogeneity. We quantify the likelihood of exceeding selected ozone levels, finding that parts of the Middle East, India, and China are most likely to exceed 55 parts per billion (ppb) in 2017. About 96% of the global population was exposed to ozone levels above the World Health Organization guideline of 60 µg m −3 (30 ppb) in 2017. Our annual fine-resolution ozone estimates may be useful for several applications including epidemiology and assessments of impacts on health, agriculture, and ecosystems.
Emissions are a central component of atmospheric chemistry models. The Harmonized Emissions Component (HEMCO) is a software component for computing emissions from a user-selected ensemble of emission inventories and algorithms. It allows users to re-grid, combine, overwrite, subset, and scale emissions from different inventories through a configuration file and with no change to the model source code. The configuration file also maps emissions to model species with appropriate units. HEMCO can operate in offline stand-alone mode, but more importantly it provides an online facility for models to compute emissions at runtime. HEMCO complies with the Earth System Modeling Framework (ESMF) for portability across models. We present a new version here, HEMCO 3.0, that features an improved three-layer architecture to facilitate implementation into any atmospheric model and improved capability for calculating emissions at any model resolution including multiscale and unstructured grids. The three-layer architecture of HEMCO 3.0 includes (1) the Data Input Layer that reads the configuration file and accesses the HEMCO library of emission inventories and other environmental data, (2) the HEMCO Core that computes emissions on the user-selected HEMCO grid, and (3) the Model Interface Layer that re-grids (if needed) and serves the data to the atmospheric model and also serves model data to the HEMCO Core for computing emissions dependent on model state (such as from dust or vegetation). The HEMCO Core is common to the implementation in all models, while the Data Input Layer and the Model Interface Layer are adaptable to the model environment. Default versions of the Data Input Layer and Model Interface Layer enable straightforward implementation of HEMCO in any simple model architecture, and options are available to disable features such as re-gridding that may be done by independent couplers in more complex architectures. The HEMCO library of emission inventories and algorithms is continuously enriched through user contributions so that new inventories can be immediately shared across models. HEMCO can also serve as a general data broker for models to process input data not only for emissions but for any gridded environmental datasets. We describe existing implementations of HEMCO 3.0 in (1) the GEOS-Chem “Classic” chemical transport model with shared-memory infrastructure, (2) the high-performance GEOS-Chem (GCHP) model with distributed-memory architecture, (3) the NASA GEOS Earth System Model (GEOS ESM), (4) the Weather Research and Forecasting model with GEOS-Chem (WRF-GC), (5) the Community Earth System Model Version 2 (CESM2), and (6) the NOAA Global Ensemble Forecast System – Aerosols (GEFS-Aerosols), as well as the planned implementation in the NOAA Unified Forecast System (UFS). Implementation of HEMCO in CESM2 contributes to the Multi-Scale Infrastructure for Chemistry and Aerosols (MUSICA) by providing a common emissions infrastructure to support different simulations of atmospheric chemistry across scales.
We have analysed time-slice simulations from 17 global models, participating in the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP), to explore changes in present-day (2000) hydroxyl radical (OH) concentration and methane (CH4) lifetime relative to preindustrial times (1850) and to 1980. A comparison of modeled and observation-derived methane and methyl chloroform lifetimes suggests that the present-day global multi-model mean OH concentration is overestimated by 5 to 10% but is within the range of uncertainties. The models consistently simulate higher OH concentrations in the Northern Hemisphere (NH) compared with the Southern Hemisphere (SH) for the present-day (2000; inter-hemispheric ratios of 1.13 to 1.42), in contrast to observation-based approaches which generally indicate higher OH in the SH although uncertainties are large. Evaluation of simulated carbon monoxide (CO) concentrations, the primary sink for OH, against ground-based and satellite observations suggests low biases in the NH that may contribute to the high north–south OH asymmetry in the models. The models vary widely in their regional distribution of present-day OH concentrations (up to 34%). Despite large regional changes, the multi-model global mean (mass-weighted) OH concentration changes little over the past 150 yr, due to concurrent increases in factors that enhance OH (humidity, tropospheric ozone, nitrogen oxide (NOx) emissions, and UV radiation due to decreases in stratospheric ozone), compensated by increases in OH sinks (methane abundance, carbon monoxide and non-methane volatile organic carbon (NMVOC) emissions). The large inter-model diversity in the sign and magnitude of preindustrial to present-day OH changes (ranging from a decrease of 12.7% to an increase of 14.6%) indicate that uncertainty remains in our understanding of the long-term trends in OH and methane lifetime. We show that this diversity is largely explained by the different ratio of the change in global mean tropospheric CO and NOx burdens (Delta CO/Delta NOx, approximately represents changes in OH sinks versus changes in OH sources) in the models, pointing to a need for better constraints on natural precursor emissions and on the chemical mechanisms in the current generation of chemistry-climate models. For the 1980 to 2000 period, we find that climate warming and a slight increase in mean OH (3.5 +/- 2.2%) leads to a 4.3 +/- 1.9% decrease in the methane lifetime. Analysing sensitivity simulations performed by 10 models, we find that preindustrial to present-day climate change decreased the methane lifetime by about four months, representing a negative feedback on the climate system. Further, we analysed attribution experiments performed by a subset of models relative to 2000 conditions with only one precursor at a time set to 1860 levels. We find that global mean OH increased by 46.4 +/- 12.2% in response to preindustrial to present-day anthropogenic NOx emission increases, and decreased by 17.3 +/-2.3%, 7.6 +/- 1.5%, and 3.1 +/- 3.0% due to methane burden, and anthropogenic CO, and NMVOC emissions increases, respectively.
The GSFC two-dimensional transport and chemistry model has been used for a wide variety of scientific and assessment studies of stratospheric ozone. Transport is a key element in the ozone simulations, and we have recently upgraded our model transport formulation to include much of the information about atmospheric transport processes available from existing data sets. To properly evaluate the model transport, it is desirable to examine the effects of transport and photochemistry separately. Recently, high quality observations of several long lived stratospheric tracers have become available from aircraft, balloon, and satellite measurement systems. This data provides a means to do a detailed model transport evaluation, as has been done in the recent Models and Measurements Intercomparison Project II. In this paper, we will discuss the GSFC 2D model simulations of ozone together with model-data comparisons of long lived tracers such as methane and the age of air transport diagnostic. We will show that the model can reproduce many of the transport-sensitive features observed in the stratosphere, and can compare reasonably well with measurements of both total ozone and long lived tracers simultaneously. We will also discuss the model deficiencies in simulating some of the detailed aspects of the observations.
Measurements from actinic flux spectroradiometers on board the NASA DC-8 during the Atmospheric Tomography (ATom) mission provide an extensive set of statistics on how clouds alter photolysis rates (J values) throughout the remote Pacific and Atlantic Ocean basins. J values control tropospheric ozone and methane abundances, and thus clouds have been included for more than three decades in tropospheric chemistry modeling. ATom made four profiling circumnavigations of the troposphere capturing each of the seasons during 2016–2018. This work examines J values from the Pacific Ocean flights of the first deployment, but publishes the complete Atom-1 data set (29 July to 23 August 2016). We compare the observed J values (every 3 s along flight track) with those calculated by nine global chemistry– climate/transport models (globally gridded, hourly, for a mid-August day). To compare these disparate data sets, we build a commensurate statistical picture of the impact of clouds on J values using the ratio of J -cloudy (standard, sometimes cloudy conditions) to J -clear (artificially cleared of clouds). The range of modeled cloud effects is inconsistently large but they fall into two distinct classes: (1) models with large cloud effects showing mostly enhanced J values aloft and or diminished at the surface and (2) models with small effects having nearly clear-sky J values much of the time. The ATom-1 measurements generally favor large cloud effects but are not precise or robust enough to point out the best cloud-modeling approach. The models here have resolutions of 50–200 km and thus reduce the occurrence of clear sky when averaging over grid cells. In situ measurements also average scattered sunlight over a mixed cloud field, but only out to scales of tens of kilometers. A primary uncertainty remains in the role of clouds in chemistry, in particular, how models average over cloud fields, and how such averages can simulate measurements.
Numerical solutions for the Aeroassist Flight Experiment vehicle were obtained from three methods at a trajectory point corresponding to the maximum aerodynamic heating. The flow regime and vehicle's speed require a viscous model of laminar flow and finite-rate chemical and thermal modeling. The computational domain covers both forebody and base such that the shock layer and near wake flowfield are included. Because of differences in computational grids, methods of solution, and models of rate equations, the results are generally in poor agreement. Temperature and species concentrations are strongly affected by the physical model equations and associated parameters. The chemistry model based on 11 species is found to yield lower translational temperature profile near stagnation than those from a seven-species model. The vibrational temperature varies according to the modeling details. All solutions indicate that strong neutral and/or molecular dissociation and weak ionization take place at the forebody and vibrational freezing is present in the afterbody expansion region where the vibrational temperature is higher than the translational temperature. Some forms of shear layer emanating from the aerobrake skirt coalesce in the region of reversed flow behind the vehicle.
The chemical nonequilibrium flow around a projectile propelled by supersonic combustion inside a tube of premixed, stoichiometric H2-air mixture is examined using a quasi-one-dimensional flow model and a 28-reactions/13 species chemistry model. Two series of calculations have been performed; one at the initial pressure 5 atm and the other 0.01 atm. Three projectile Mach numbers, namely 8, 10 and 5 were used. The minimum cone angle of the projectile forebody needed for triggering ignition can be estimated from the results of the calculations. The numerical results also show that high initial pressure reduces the induction and reaction zones to a single detonation front, and thrust is produced by the high pressure in the combustor and nozzle. Low initial pressure lengthens the reaction zone. In all the cases, the momentum of the flow changes only slightly because the tube in which the projectile travels limits the area for the nozzle flow to expand.
An axisymmetric single engine Computational Fluid Dynamics calculation of the Saturn V/S 1-C vehicle base region and F1 engine plume is described. There were two objectives of this work, the first was to calculate an axisymmetric approximation of the nozzle, plume and base region flow fields of S1-C/F1, relate/scale this to flight data and apply this scaling factor to a NLS/STME axisymmetric calculations from a parallel effort. The second was to assess the differences in F1 and STME plume shear layer development and concentration of combustible gases. This second piece of information was to be input/supporting data for assumptions made in NLS2 base temperature scaling methodology from which the vehicle base thermal environments were being generated. The F1 calculations started at the main combustion chamber faceplate and incorporated the turbine exhaust dump/nozzle film coolant. The plume and base region calculations were made for ten thousand feet and 57 thousand feet altitude at vehicle flight velocity and in stagnant freestream. FDNS was implemented with a 14 species, 28 reaction finite rate chemistry model plus a soot burning model for the RP-1/LOX chemistry. Nozzle and plume flow fields are shown, the plume shear layer constituents are compared to a STME plume. Conclusions are made about the validity and status of the analysis and NLS2 vehicle base thermal environment definition methodology.
We describe the coupling of the Goddard Institute for Space Studies (GISS) general circulation model (GCM) to an online sulfur chemistry model and source models for organic matter and sea-salt that is used to estimate the aerosol indirect effect. The cloud droplet number concentration is diagnosed empirically from field experiment datasets over land and ocean that observe droplet number and all three aerosol types simultaneously; corrections are made for implied variations in cloud turbulence levels. The resulting cloud droplet number is used to calculate variations in droplet effective radius, which in turn allows us to predict aerosol effects on cloud optical thickness and microphysical process rates. We calculate the aerosol indirect effect by differencing the top-of-the-atmosphere net cloud radiative forcing for simulations with present-day vs. pre-industrial emissions. Both the first (radiative) and second (microphysical) indirect effects are explored. We test the sensitivity of our results to cloud parameterization assumptions that control the vertical distribution of cloud occurrence, the autoconversion rate, and the aerosol scavenging rate, each of which feeds back significantly on the model aerosol burden. The global mean aerosol indirect effect for all three aerosol types ranges from -1.55 to -4.36 W/sq m in our simulations. The results are quite sensitive to the pre-industrial background aerosol burden, with low pre-industrial burdens giving strong indirect effects, and to a lesser extent to the anthropogenic aerosol burden, with large burdens giving somewhat larger indirect effects. Because of this dependence on the background aerosol, model diagnostics such as albedo-particle size correlations and column cloud susceptibility, for which satellite validation products are available, are not good predictors of the resulting indirect effect.
We describe the coupling of the Goddard Institute for Space Studies (GISS) general circulation model (GCM) to an online sulfur chemistry model and source models for organic matter and sea-salt that is used to estimate the aerosol indirect effect. The cloud droplet number concentration is diagnosed empirically from field experiment datasets over land and ocean that observe droplet number and all three aerosol types simultaneously; corrections are made for implied variations in cloud turbulence levels. The resulting cloud droplet number is used to calculate variations in droplet effective radius, which in turn allows us to predict aerosol effects on cloud optical thickness and microphysical process rates. We calculate the aerosol indirect effect by differencing the top-of-the-atmosphere net cloud radiative forcing for simulations with present-day vs. pre-industrial emissions. Both the first (radiative) and second (microphysical) indirect effects are explored. We test the sensitivity of our results to cloud parameterization assumptions that control the vertical distribution of cloud occurrence, the autoconversion rate, and the aerosol scavenging rate, each of which feeds back significantly on the model aerosol burden. The global mean aerosol indirect effect for all three aerosol types ranges from -1.55 to -4.36 W m(exp -2) in our simulations. The results are quite sensitive to the pre-industrial background aerosol burden, with low pre-industrial burdens giving strong indirect effects, and to a lesser extent to the anthropogenic aerosol burden, with large burdens giving somewhat larger indirect effects. Because of this dependence on the background aerosol, model diagnostics such as albedo-particle size correlations and column cloud susceptibility, for which satellite validation products are available, are not good predictors of the resulting indirect effect.
AdapChem software enables high efficiency, low computational cost, and enhanced accuracy on computational fluid dynamics (CFD) numerical simulations used for combustion studies. The software dynamically allocates smaller, reduced chemical models instead of the larger, full chemistry models to evolve the calculation while ensuring the same accuracy to be obtained for steady-state CFD reacting flow simulations. The software enables detailed chemical kinetic modeling in combustion CFD simulations. AdapChem adapts the reaction mechanism used in the CFD to the local reaction conditions. Instead of a single, comprehensive reaction mechanism throughout the computation, a dynamic distribution of smaller, reduced models is used to capture accurately the chemical kinetics at a fraction of the cost of the traditional single-mechanism approach.
The Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP) consists of a series of time slice experiments targeting the long-term changes in atmospheric composition between 1850 and 2100, with the goal of documenting composition changes and the associated radiative forcing. In this overview paper, we introduce the ACCMIP activity, the various simulations performed (with a requested set of 14) and the associated model output. The 16 ACCMIP models have a wide range of horizontal and vertical resolutions, vertical extent, chemistry schemes and interaction with radiation and clouds. While anthropogenic and biomass burning emissions were specified for all time slices in the ACCMIP protocol, it is found that the natural emissions are responsible for a significant range across models, mostly in the case of ozone precursors. The analysis of selected present-day climate diagnostics (precipitation, temperature, specific humidity and zonal wind) reveals biases consistent with state-of-the-art climate models. The model-to- model comparison of changes in temperature, specific humidity and zonal wind between 1850 and 2000 and between 2000 and 2100 indicates mostly consistent results. However, models that are clear outliers are different enough from the other models to significantly affect their simulation of atmospheric chemistry.
We use the LaRC Lagrangian Chemistry and Transport Model (LCTM) [Considine et al., 2007; Pierce et al., 2003] to intercompare ACE, Aura, and HALOE observations of long-lived trace species. The LCTM calculates the transport, mixing, and photochemical evolution of an ensemble of parcels that have been initialized from ACE-FTS measurements. Here we focus on late November, 2004 comparisons, due to the previous 3-week period of continuous HALOE observations and MLS v2.2 data on November 29, 2004.
Sulfur dioxide (SO2) is a major atmospheric pollutant with a strong anthropogenic component mostly produced by the combustion of fossil fuel and other industrial activities. As a precursor of sulfate aerosols that affect climate, air quality, and human health, this gas needs to be monitored on a global scale. Global climate and chemistry models including aerosol processes along with their radiative effects are important tools for climate and air quality research. Validation of these models against in-situ and satellite measurements is essential to ascertain the credibility of these models and to guide model improvements. In this study the Goddard Chemistry, Aerosol, Radiation, and Transport (GOCART) module running on-line inside the Goddard Earth Observing System version 5 (GEOS-5) model is used to simulate aerosol and SO2 concentrations. Data taken in November 2010 over Frostburg, Maryland during an SO2 field campaign involving ground instrumentation and aircraft are used to evaluate GEOS-5 simulated SO2 concentrations. Preliminary data analysis indicated the model overestimated surface SO2 concentration, which motivated the examination of mixing processes in the model and the specification of SO2 anthropogenic emission rates. As a result of this analysis, a revision of anthropogenic emission inventories in GEOS-5 was implemented, and the vertical placement of SO2 sources was updated. Results show that these revisions improve the model agreement with observations locally and in regions outside the area of this field campaign. In particular, we use the ground-based measurements collected by the United States Environmental Protection Agency (US EPA) for the year 2010 to evaluate the revised model simulations over North America.
Estimates of ground-level ozone concentrations are necessary to determine the human health burden of ozone. To support the Global Burden of Disease Study, we produce yearly fine resolution global surface ozone estimates from 1990 to 2017 through a data fusion of observations and models. As ozone observations are sparse in many populated regions, we use a novel combination of the M3Fusion and Bayesian Maximum Entropy (BME) methods. With M3Fusion, we create a multi-model composite by bias-correcting and weighting nine global atmospheric chemistry models based on their ability to predict observations (8,834 sites globally)in each region and year. BME is then used to integrate observations, such that estimates match observations at each monitoring site with the observational influence decreasing smoothly across space and time until the output matches the multi-model composite. After estimating at 0.5° resolution using BME, we add fine spatial detail from an additional model, yielding estimates at 0.1° resolution. Observed ozone is predicted more accurately (R2=0.81 at test point, 0.63 at 0.1°,0.62 at 0.5°) than the multi-model mean (R2=0.28 at 0.5°). Global ozone exposure is estimated to be increasing, driven by highly populated regions of Asia and Africa, despite decreases in the United States and Russia.