Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numeric profiling”

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 1,567 records · Page 87

Distributed Disdrometer and Rain Gauge Measurement Infrastructure Developed for GPM Ground Validation

Global Precipitation Mission (GPM)retrieval algorithm validation requires datasets characterizing the 4-D structure, variability, and correlation properties of hydrometeor particle size distributions (PSD) and accumulations over satellite fields of view (FOV;<10 km). Collection of this data provides a means to assess retrieval errors related to beam filling and algorithm PSD assumptions. Hence, GPM Ground Validation is developing a deployable network of precipitation gauges and disdrometers to provide fine-scale measurements of PSD and precipitation accumulation variability. These observations will be combined with dual-frequency, polarimetric, and profiling radar data in a bootstrapping fashion to extend validated PSD measurements to a large coverage domain. Accordingly, a total of 24 Parsivel disdrometers(PD), 5 3rd-generation 2D Video Disdrometers (2DVD), 70 tipping bucket rain gauges (TBRG),9 weighing gauges, 7 Hot-Plate precipitation sensors (HP), and 3 Micro Rain Radars (MRR) have been procured. In liquid precipitation the suite of TBRG, PD and 2DVD instruments will quantify a broad spectrum of rain rate and PSD variability at sub-kilometer scales. In the envisioned network configuration 5 2DVDs will act as reference points for 16 collocated PD and TBRG measurements. We find that PD measurements provide similar measures of the rain PSD as observed with collocated 2DVDs (e.g., D0, Nw) for rain rates less than 15 mm/hr. For heavier rain rates we will rely on 2DVDs for PSD information. For snowfall we will combine point-redundant observations of SWER distributed over three or more locations within a FOV. Each location will contain at least one fenced weighing gauge, one HP, two PDs, and a 2DVD. MRRs will also be located at each site to extend the measurement to the column. By collecting SWER measurements using different instrument types that employ different measurement techniques our objective is to separate measurement uncertainty from natural variability in SWER and PSD. As demonstrated using C3VP polarimetric radar, gauge, and 2DVD/PD datasets these measurements can be combined to bootstrap an area wide SWER estimate via constrained modification of density-diameter and radar reflectivity-snowfall relationships. These data will be combined with snowpack, airborne microphysics, radar, radiometer, and tropospheric sounding data to refine GPM snowfall retrievals. The gauge and disdrometer instruments are being developed to operate autonomously when necessary using solar power and wireless communications. These systems will be deployed in numerous field campaigns through 2016. Planned deployment of these systems include field campaigns in Finland (2010), Oklahoma (2011), Canada (2012) and North Carolina (2013). GPM will also deploy 20 pairs of TBRGs within a 25 km2 region along the Virginia coast under NASA NPOL radar coverage in order to quantify errors in point-area rainfall measurements.

Petersen, Walter A.↗

Evaluation of the Terminal Area Precision Scheduling and Spacing System for Performance-Based Navigation Arrivals

The growth of global demand for air transportation has put increasing strain on the nation's air traffic management system. To relieve this strain, the International Civil Aviation Organization has urged all nations to adopt Performance-Based Navigation (PBN), which can help to reduce air traffic congestion, decrease aviation fuel consumption, and protect the environment. NASA has developed a Terminal Area Precision Scheduling and Spacing (TAPSS) system that can support increased use of PBN during periods of high traffic, while supporting fuel-efficient, continuous descent approaches. In the original development of this system, arrival aircraft are assigned fuel-efficient Area Navigation (RNAV) Standard Terminal Arrival Routes before their initial descent from cruise, with routing defined to a specific runway. The system also determines precise schedules for these aircraft that facilitate continuous descent through the assigned routes. To meet these schedules, controllers are given a set of advisory tools to precisely control aircraft. The TAPSS system has been evaluated in a series of human-in-the-loop (HITL) air traffic simulations during 2010 and 2011. Results indicated increased airport arrival throughput up to 10 over current operations, and maintained fuel-efficient aircraft decent profiles from the initial descent to landing with reduced controller workload. This paper focuses on results from a joint NASA and FAA HITL simulation conducted in 2012. Due to the FAA rollout of the advance terminal area PBN procedures at mid-sized airports first, the TAPSS system was modified to manage arrival aircraft as they entered Terminal Radar Approach Control (TRACON). Dallas-Love Field airport (DAL) was selected by the FAA as a representative mid-sized airport within a constrained TRACON airspace due to the close proximity of a major airport, in this case Dallas-Ft Worth International Airport, one of the busiest in the world. To address this constraint, RNAV routes and Required Navigation Performance with the particular capability known as Radius-to-Fix (RNP-RF) approaches to a short final were used. The purpose of this simulation was to get feedback on how current operations could benefit with the TAPSS system and also to evaluate the efficacy of the advisory tools to support the broader use of PBN in the US National Airspace System. For this NASA-FAA joint experiment, an Air Traffic Control laboratory at NASA Ames was arranged to simulate arrivals into DAL in Instrument Meteorological Conditions utilizing parallel dependent approaches, with two feeder positions that handed off traffic to one final position. Four FAA controllers participated, alternately covering these three positions. All participants were Full-Performance Level terminal controllers and members of the National Air Traffic Controllers Association. During the simulation, PBN arrival operations were compared and contrasted in three conditions. They were the Baseline, where none of the TAPSS systems TRACON controller decision support advisories were provided, the Limited Advisories, reflecting the existing but dormant capabilities of the current terminal automation equipment with providing a subset of the TAPSS systems advisories; numerical delay, landing sequence, and runway assignment information, and the Full Advisories, with providing the following in addition to the ones in the Limited condition; trajectory slot markers, timelines of estimated times of arrivals and sched

Performance based navigation↗

Impact of MODIS High-Resolution Sea-Surface Temperatures on WRF Forecasts at NWS Miami, FL

Over the past few years,studies at the Short-term Prediction Research and Transition (SPoRT) Center have suggested that the use of Moderate Resolution Imaging Spectroradiometer (MODIS) composite sea-surface temperature (SST) products in regional weather forecast models can have a significant positive impact on short-term numerical weather prediction in coastal regions. The recent paper by LaCasse et al. (2007, Monthly Weather Review) highlights lower atmospheric differences in regional numerical simulations over the Florida offshore waters using 2-km SST composites derived from the MODIS instrument aboard the polar-orbiting Aqua and Terra Earth Observing System satellites. To help quantify the value of this impact on NWS Weather Forecast Offices (WFOs), the SPoRT Center and the NWS WFO at Miami, FL (MIA) are collaborating on a project to investigate the impact of using the high-resolution MODIS SST fields within the Weather Research and Forecasting (WRF) prediction system. The scientific hypothesis being tested is: More accurate specification of the lower-boundary forcing within WRF will result in improved land/sea fluxes and hence, more accurate evolution of coastal mesoscale circulations and the associated sensible weather elements. The NWS MIA is currently running the WRF system in real-time to support daily forecast operations, using the National Centers for Environmental Prediction Nonhydrostatic Mesoscale Model dynamical core within the NWS Science and Training Resource Center's Environmental Modeling System (EMS) software; The EMS is a standalone modeling system capable of downloading the necessary daily datasets, and initializing, running and displaying WRF forecasts in the NWS Advanced Weather Interactive Processing System (AWIPS) with little intervention required by forecasters. Twenty-seven hour forecasts are run daily with start times of 0300,0900, 1500, and 2100 UTC on a domain with 4-km grid spacing covering the southern half of Florida and the far western portions of the Bahamas, the Florida Keys, the Straights of Florida, and adjacent waters of the Gulf of Mexico and Atlantic Ocean. Each model run is initialized using the Local Analysis and Prediction System (LAPS) analyses available in AWIPS, invoking the diabatic. "hot-start" capability. In this WRF model "hot-start", the LAPS-analyzed cloud and precipitation features are converted into model microphysics fields with enhanced vertical velocity profiles, effectively reducing the model spin-up time required to predict precipitation systems. The SSTs are initialized with the NCEP Real-Time Global (RTG) analyses at l/12 degree resolution (approx. 9 km); however, the RTG product does not exhibit fine-scale details consistent with its grid resolution. SPoRT is conducting parallel WRF EMS runs identical to the operational runs at NWS MIA in every respect except for the use of MODIS SST composites in place of the RTG product as the initial and boundary conditions over water. The MODIS SST composites for initializing the SPoRT WRF runs are generated on a 2-km grid four times daily at 0400, 0700, 1600, and 1900 UTC, based on the times of the overhead passes of the Aqua and Terra satellites. The incorporation of the MODIS SST composites into the SPoRTWRF runs is staggered such that the 0400UTC composite initializes the 0900 UTC WRF, the 0700 UTC composite initializes the 1500 UTC WRF, the 1600 UTC composite initializes the 2100 UTC WRF, and the 1900 UTC composite initializes the 0300 UTC WRF. A comparison of the SPoRT and Miami forecasts is underway in 2007, and includes quantitative verification of near-surface temperature, dewpoint, and wind forecasts at surface observation locations. In addition, particular days of interest are being analyzed to determine the impact of the MODIS SST data on the development and evolution of predicted sea/land-breeze circulations, clouds, and precipitation. This paper will present verification results comparing the NWS MIA forecasts the SPoRT experimental WRF forecasts, and highlight any substantial differences noted in the predicted mesoscale phenomena.

Case, Jonathan L.↗

Mid-Infrared Study of Samples from Multiple Stones from the Sutters Mill Meteorite

The Sutter's Mill meteorite fell in N. California on April 22, 2012 and numerous pieces have been recovered and studied. We present Fourier transform infrared (FTIR) spectra of fragments from several stones of the meteorite. Methods and analysis: Infrared spectra of the samples were recorded with a Nicolet iN10 MX FTIR microscope in the mid-IR range (4000-675/cm; spectral resolution 4/cm). All samples were deposited on a clean glass slide, crushed with a stainless steel roller tool, and placed directly on the focal plane of the microscope. IR spectra were collected by averaging 128 scans. Results: Preliminary IR spectra of the non-fusion crust samples show mineral compositions that are dominated by phyllosilicates, carbonates, or mixtures of both [2]. The carbonates display a dominant, broad band centered at 1433/cm, with additional bands at 2515/cm, 1797/cm, 882/cm, and 715/cm). Features associated with phyllosilicates include a symmetric Si-O stretching mode band centered at 1011/cm and several O-H stretching mode bands. The O-H shows up in two forms (1) a broadband centered at 3415/cm that is probably largely due to adsorbed H2O and (2) a much weaker, narrower feature centered near 3680/cm due to structural -OH. Features observed in the 2985-2855/cm range suggest the presence of aliphatic -CH3 and -CH2- groups. The relative intensities of the bands in this range are somewhat unusual. Typically, the asymmetric aliphatic CH stretching bands are stronger than the symmetric stretching bands, but in this case the reverse is true. This unusual pattern is well matched by the aliphatic features seen in the spectrum of a terrestrial calcite (CaCO3) standard. This observation, and the fact that the strength of the carbonate and aliphatic bands seem to correlate, suggest the organics are associated with the carbonates. Conclusions: IR spectra of samples from the Sutter's Mill meteorite show absorption features associated with carbonates, phyllosilicates, and organics. Both the unusual band profile of the aliphatic C-H stretching mode bands and their correlation with the strength of the carbonate bands suggests the organics and the carbonates are associated in some manner.

Sutter's mill↗

Dynamics and Structure of Dusty Reacting Flows: Inert Particles in Strained, Laminar, Premixed Flames

A detailed numerical study was conducted on the dynamics and thermal response of inert, spherical particles in strained, laminar, premixed hydrogen/air flames. The modeling included the solution of the steady conservation equations for both the gas and particle phases along and around the stagnation streamline of an opposed-jet configuration, and the use of detailed descriptions of chemical kinetics and molecular transport, For the gas phase, the equations of mass, momentum, energy, and species are considered, while for the particle phase, the model is based on conservation equations of the particle momentum balance in the axial and radial direction, the particle number density, and the particle thermal energy equation. The particle momentum equation includes the forces as induced by drag, thermophoresis, and gravity. The particle thermal energy equation includes the convective/conductive heat exchange between the two phases, as well as radiation emission and absorption by the particle. A one-point continuation method is also included in the code that allows for the description of turning points, typical of ignition and extinction behavior. As expected, results showed that the particle velocity can be substantially different than the gas phase velocity, especially in the presence of large temperature gradients and large strain rates. Large particles were also found to cross the gas stagnation plane, stagnate, and eventually reverse as a result of the opposing gas phase velocity. It was also shown that the particle number density varies substantially throughout the flowfield, as a result of the straining of the flow and the thermal expansion. Finally, for increased values of the particle number density, substantial flame cooling to extinction states and modification of the gas phase fluid mechanics were observed. As also expected, the effect of gravity was shown to be important for low convective velocities and heavy particles. Under such conditions, simulations indicate that the magnitude and direction of the gravitational force can substantially affect the profiles of the particle velocity, number density, mass flux, and temperature.

Egolfopoulos, Fokion N.↗

Dynamics and Structure of Dusty Reacting Flows: Inert Particles in Strained, Laminar, Premixed Flames

A detailed numerical study was conducted on the dynamics and thermal response of inert spherical particles in strained, laminar, premixed hydrogen/air flames. The modeling included the solution of the steady conservation equations for both the gas and particle phases along and around the stagnation streamline of an opposed-jet configuration, and the use of detailed descriptions of chemical kinetics and molecular transport. For the gas phase, the equations of mass, momentum, energy, and species are considered, while for the particle phase, the model is based on conservation equations of the particle momentum balance in the axial and radial direction, the particle number density, and the particle thermal energy equation. The particle momentum equation includes the forces as induced by drag, thermophoresis, and gravity. The particle thermal energy equation includes the convective/conductive heat exchange between the two phases, as well as radiation emission and absorption by the particle. A one-point continuation method is also included in the code that allows for the description of turning points, typical of ignition and extinction behavior. As expected, results showed that the particle velocity can be substantially different than the gas phase velocity, especially in the presence of large temperature gradients and large strain rates. Large particles were also found to cross the gas stagnation plane, stagnate, and eventually reverse as a result of the opposing gas phase velocity. It was also shown that the particle number density varies substantially throughout the flowfield, as a result of the straining of the flow and the thermal expansion. Finally, for increased values of the particle number density, substantial flame cooling to extinction states and modification of the gas phase fluid mechanics were observed. As also expected, the effect of gravity was shown to be important for low convective velocities and heavy particles. Under such conditions, simulations indicate that the magnitude and direction of the gravitational force can substantially affect the profiles of the particle velocity, number density, mass flux, and temperature.

Egolfopoulos, Fokion N.↗

Meso-Beta scale numerical simulation studies of terrain-induced jet streak mass/momentum perturbations

The Control-B simulation experiment was designed to improve on the Control-A experiment performed with the GMASS model. This experiment addressed several inadequacies with the first smooth terrain numerical simulation by including: (1) increased nested-grid resolution to better define the simulated gravity waves, (2) increased horizontal diffusion to remove outflow boundary condition noise, and (3) the use of reanalyzed rawinsonde data and surface observations in the initial state to increase the definition of the observed jet streak as well as other low-level features. A smoothed-terrain dry simulation with the nested-grid GMASS model has revealed many important aspects of the processes which resulted in the generation of gravity waves in the region and time when and where they were observed. However, the vertical structure, number, and characteristics of the waves are still quite different from observed waves as diagnosed thus necessitating future improved simulations. However, this control simulation has produced substantial insight into processes which occur on many spatial scales over a 30 hour time period thus allowing one to draw promising inferences as to the mechanisms for the complex process which occurred in nature during the CCOPE case study. The theoretical aspects of the project have focus on understanding the nature of the ageostrophic circulations which are produced in idealized models of the atmosphere in which the troposphere is modeled in one of two ways. The first model assumes that the lower atmosphere can be represented as a single layer of homogeneous fluid whose upper surface is free to exhibit vertical displacement. Two-dimensional internal convergence (divergence) occurring during the adjustment to an asymptotic equilibrium state from an ageostrophic initial state whose momentum structure is representative of a midlatitude localized zonal wind anomaly will cause the free upper surface of the homogeneous atmosphere to rise (fall), and therefore the response can be viewed as being physically three-dimensional. The second model assumes that the troposphere can be represented by an unbounded continuously stratified Boussinesq fluid of constant Brunt-Vaisala frequency N, where the vertical gradient of the basic state potential temperature profile allows for the existence of vertically propagating internal inertia-gravity waves.

Lin, Yuh-Lang↗

Free-Energy Calculations. A Mathematical Perspective

Ion channels are pore-forming assemblies of transmembrane proteins that mediate and regulate ion transport through cell walls. They are ubiquitous to all life forms. In humans and other higher organisms they play the central role in conducting nerve impulses. They are also essential to cardiac processes, muscle contraction and epithelial transport. Ion channels from lower organisms can act as toxins or antimicrobial agents, and in a number of cases are involved in infectious diseases. Because of their important and diverse biological functions they are frequent targets of drug action. Also, simple natural or synthetic channels find numerous applications in biotechnology. For these reasons, studies of ion channels are at the forefront of biophysics, structural biology and cellular biology. In the last decade, the increased availability of X-ray structures has greatly advanced our understanding of ion channels. However, their mechanism of action remains elusive. This is because, in order to assist controlled ion transport, ion channels are dynamic by nature, but X-ray crystallography captures the channel in a single, sometimes non-native state. To explain how ion channels work, X-ray structures have to be supplemented with dynamic information. In principle, molecular dynamics (MD) simulations can aid in providing this information, as this is precisely what MD has been designed to do. However, MD simulations suffer from their own problems, such as inability to access sufficiently long time scales or limited accuracy of force fields. To assess the reliability of MD simulations it is only natural to turn to the main function of channels - conducting ions - and compare calculated ionic conductance with electrophysiological data, mainly single channel recordings, obtained under similar conditions. If this comparison is satisfactory it would greatly increase our confidence that both the structures and our computational methodologies are sufficiently accurate. Channel conductance, defined as the ratio of ionic current through the channel to applied voltage, can be calculated in MD simulations by way of applying an external electric field to the system and counting the number of ions that traverse the channel per unit time. If the current is small, a voltage significantly higher than the experimental one needs to be applied to collect sufficient statistics of ion crossing events. Then, the calculated conductance has to be extrapolated to the experimental voltage using procedures of unknown accuracy. Instead, we propose an alternative approach that applies if ion transport through channels can be described with sufficient accuracy by the one-dimensional diffusion equation in the potential given by the free energy profile and applied voltage. Then, it is possible to test the assumptions of the equation, recover the full voltage/current dependence, determine the reliability of the calculated conductance and reconstruct the underlying (equilibrium) free energy profile, all from MD simulations at a single voltage. We will present the underlying theory, model calculations that test this theory and simulations on ion conductance through a channel that has been extensively studied experimentally. To our knowledge this is the first case in which the complete, experimentally measured dependence of the current on applied voltage has been reconstructed from MD simulations.

free energy↗

4-D Cloud Water Content Fields Derived from Operational Satellite Data

In order to improve operational safety and efficiency, the transportation industry, including aviation, has an urgent need for accurate diagnoses and predictions of clouds and associated weather conditions. Adverse weather accounts for 70% of all air traffic delays within the U.S. National Airspace System. The Federal Aviation Administration has determined that as much as two thirds of weather-related delays are potentially avoidable with better weather information and roughly 20% of all aviation accidents are weather related. Thus, it is recognized that an important factor in meeting the goals of the Next Generation Transportation System (NexGen) vision is the improved integration of weather information. The concept of a 4-D weather cube is being developed to address that need by integrating observed and forecasted weather information into a shared 4-D database, providing an integrated and nationally consistent weather picture for a variety of users and to support operational decision support systems. Weather analyses and forecasts derived using Numerical Weather Prediction (NWP) models are a critical tool that forecasters rely on for guidance and also an important element in current and future decision support systems. For example, the Rapid Update Cycle (RUC) and the recently implemented Rapid Refresh (RR) Weather Research and Forecast (WRF) models provide high frequency forecasts and are key elements of the FAA Aviation Weather Research Program. Because clouds play a crucial role in the dynamics and thermodynamics of the atmosphere, they must be adequately accounted for in NWP models. The RUC, for example, cycles at full resolution five cloud microphysical species (cloud water, cloud ice, rain, snow, and graupel) and has the capability of updating these fields from observations. In order to improve the models initial state and subsequent forecasts, cloud top altitude (or temperature, T(sub c)) derived from operational satellite data, surface observations of cloud base altitude, radar reflectivity, and lightning data are used to help build and remove clouds in the models assimilation system. Despite this advance and the many recent advances made in our understanding of cloud physical processes and radiative effects, many problems remain in adequately representing clouds in models. While the assimilation of cloud top information derived from operational satellite data has merit, other information is available that has not yet been exploited. For example, the vertically integrated cloud water content (CWC) or cloud water path (CWP) and cloud geometric thickness (delta Z) are standard products being derived routinely from operational satellite data. These and other cloud products have been validated under a variety of conditions. Since the uncertainties have generally been found to be less than those found in model analyses and forecasts, the satellite products should be suitable for data assimilation, provided an appropriate strategy can be developed that links the satellite-derived cloud parameters with cloud parameters specified in the model. In this paper, we briefly outline such a strategy and describe a methodology to retrieve cloud water content profiles from operational satellite data. Initial results and future plans are presented. It is expected that the direct assimilation of this new product will provide the most accurate depiction of the vertical distribution of cloud water ever produced at the high spatial and temporal resolution needed for short term weather analyses and forecasts.

Smith, William L., Jr.↗

A Meteorological Model's Dependence on Radiation Update Frequency

Numerical weather models are used to simulate circulations in the atmosphere including clouds and precipitation by applying a set of mathematical equations over a three-dimensional grid. The grid is composed of discrete points at which the meteorological variables are defined. As computing power continues to rise these models are being used at finer grid spacing, but they must still cover a wide range of scales. Some of the physics that must be accounted for in the model cannot be explicitly resolved, and their effects, therefore, must be estimated or "parameterized". Some of these parameterizations are computationally expensive. To alleviate the problem, they are not always updated at the time resolution of the model with the assumption being that the impact will be small. In this study, a coupled land-atmosphere model is used to assess the impact of less frequent updates of the computationally expensive radiation physics for a case on June 6, 2002, that occurred during a field experiment over the central plains known as International H20 Project (IHOP). The model was tested using both the original conditions, which were dry, and with modified conditions wherein moisture was added to the lower part of the atmosphere to produce clouds and precipitation (i.e., a wet case). For each of the conditions (i.e., dry and wet), four set of experiments were conducted wherein the model was run for a period of 24 hours and the radiation fields (including both incoming solar and outgoing longwave) were updated every 1, 3, 10, and 100 time steps. Statistical tests indicated that average quantities of surface variables for both the dry and wet cases were the same for the various update frequencies. However, spatially the results could be quite different especially in the wet case after it began to rain. The near-surface wind field was found to be different most of the time even for the dry case. In the wet case, rain intensities and average vertical profiles of heating associated with cloudy areas were found to differ for the various radiation update frequencies. The latter implies that the mean state of the model could be different as a result of not updating the radiation fields every time step and has important implications for longer term climate studies

Eastman, Joseph L.↗

Differences in Water Vapor Radiative Transfer Among 1D Models Can Significantly Affect the Inner Edge of the Habitable Zone

An accurate estimate of the inner edge of the habitable zone is critical for determining which exoplanets are potentially habitable and for designing future telescopes to observe them. Here, we explore differences in estimating the inner edge among seven one-dimensional radiative transfer models: two line-by-line codes (SMART and LBLRTM) as well as five band codes (CAM3, CAM4_Wolf, LMDG, SBDART, and AM2) that are currently being used in global climate models. We compare radiative fluxes and spectra in clear-sky conditions around G and M stars, with fixed moist adiabatic profiles for surface temperatures from 250 to 360 K. We find that divergences among the models arise mainly from large uncertainties in water vapor absorption in the window region (10 μm) and in the region between 0.2 and 1.5 μm. Differences in outgoing longwave radiation increase with surface temperature and reach 10–20Wm(exp −2); differences in shortwave reach up to 60 Wm(exp −2), especially at the surface and in the troposphere, and are larger for an M-dwarf spectrum than a solar spectrum. Differences between the two line-by-line models are significant, although smaller than among the band models. Our results imply that the uncertainty in estimating the insolation threshold of the inner edge (the runaway greenhouse limit) due only to clear-sky radiative transfer is ≈10% of modern Earth’s solar constant (i.e., ≈34 Wm(exp −2) in global mean) among band models and ≈3% between the two line-by-line models. These comparisons show that future work is needed that focuses on improving water vapor absorption coefficients in both shortwave and longwave, as well as on increasing the resolution of stellar spectra in broadband models.

Numerical methods↗

Global Reference Atmospheric Model (GRAM) Advancements and Additions

Introduction:The Global Reference Atmospheric Model (GRAM) is one of the most widely used engineering models of the atmosphere. GRAM development and maintenance has been led by NASA Marshall Space Flight Center (MSFC). The NASA Science Mission Directorate (SMD) has provided funding support to upgrade the GRAMs since Fiscal Year 2018. NASA Langley Research Center has been working with MSFC on the upgrades.This presentation will provide details regarding the upgrades that have been made to the existing GRAMs, the development of new GRAMs, as well as the ongoing objectives, tasks, and milestones re-lated to the GRAM upgrades funded by NASA SMD. GRAM: The GRAMs are engineering-oriented atmospheric models that estimate mean values and statistical variations of the atmospheric properties for numerous planetary destinations.They provide mean values and variability for any point in the atmosphere as well as seasonal, geographic, and altitude variations. GRAM outputs include atmospheric density, temperature, pressure, winds, and chemical composition along a user-defined path.Theyhave been widely used by the engineering community because of their ability to create realistic dispersions. GRAMs have been integrated into high fidelity flight dynamic simulations of launch, entry, descent and landing (EDL), aerobraking and aerocapture. MSFC has been developing and updating GRAMs since 1974; GRAMs are currently available for Earth, Mars, Venus, Neptune, and Titan. GRAM Upgrade Status: Code Moderization.The planetary GRAMs are being rearchitected from Fortran to a common object-ori-entedC++ framework called the GRAM Suite. This new architecture creates a common GRAM library of data models and utilities. The first C++ releases of the existing planetary GRAMs (Mars, Venus, Neptune, and Titan-GRAM) in the GRAM Suite are straight conversions from the latest Fortran version. Model Upgrades.The focus of the model upgrade task is to improve the atmosphere models in the existing GRAMs and to establish a foundation for developing GRAMs for additional destinations. The GRAM ephemerishas been upgraded to the NASA Navigation and Ancillary Information Facility (NAIF) SPICE toolkit (version N0066). The calculation of the speed of sound has also been improved in the GRAMs. The GRAM team has received updated Mars General Circulation Model (MGCM) datasets from NASA Ames Research Center. Mars Global Ionosphere-Thermosphere Model (M-GITM) data is being obtained toreplace the Mars Thermospheric General Circulation Model (MTGCM) data in legacy Mars-GRAM. M-GITM and updated MGCM data will be incorporated into a future GRAM Suite release.Twoprojects that will improve the atmospheric model data in the GRAMs have been funded by the GRAM team since Fiscal Year 2020. Sanjay Limaye and Patrick Fry at the University of Wisconsin are reanalyzing the Venus Express radio occultation observations and analyzing the Akatsukiradio occultation observations.This will lead tothe calculation of number density, temperature, and pressure profiles for the 40-90 km altitude range. Kunio Sayanagi, Justin Garland, and Ryan McCabeat Hampton University are developing empirical global models for Venus,Jupiter, Saturn, Uranus, Neptune, and Titan that incorporates the latest data available for each of these planetary destinations. Upgraded GRAM Releases. GRAM Suite Version 1.0 was released in May 2020 and contains the rearchitected Neptune-GRAM, including the common GRAM framework and planet–specific code. GRAM Suite Version 1.1 was released in September 2020 and addedthe rearchitected Titan-GRAM to the GRAM Suite. A User Guide and Programmer’s Manualarereleased with all GRAMs. The rearchitected Mars and Venus-GRAMs will be released in upcomingversions of the GRAM Suite. New GRAM Releases. New GRAMs have been de-veloped for Uranus and Jupiter. Uranus-GRAM is based on an individual profile generated by Gary Allen (ARC) from Voyager 2 occultation data and will be released in GRAM Suite Version 1.2. Jupiter-GRAM is based on individual profile produced from Al Seiff’s Ju-piter model[1].Jupiter-GRAM will be released in GRAM Suite Version 1.3. Saturn-GRAM is currently under development and will be released in a future version of the GRAM Suite. Conclusions: GRAMs are vital and frequently used toolsets. Releases of the GRAM Suite, upgrades of the existing planetary GRAMs, and development of new planetary GRAMs are ongoing. NASA SMD funding has been essential to addressing current limitations and accomplishing GRAM developmental goals. Continua-tion ofSMD fundingwill ensure the development, up-grades, and maintenance of the GRAMs. References: [1] Seiff, A., et al. (1998) JGR, 103, 22,857-22,889. Acknowledgments: The authors gratefully acknowledge support from the NASA SMD

atmospheric models↗

Lipid Biomarkers for Methanogens in Hypersaline Cyanobacterial Mats for Guerrero Negro, Baja California Sur

Analyses of sediments from the vicinity of active methane seeps have uncovered a particular suite of lipid biomarker patterns that characterize methane consuming archaea and their syntrophic, sulfate reducing partners. These isoprenoid biomarkers, largely identified by their anomalously light carbon isotopic signatures, have been a topic of intense research activity and are recorded in numerous methane-rich environments from Holocene to Cenozoic. This phenomenon has implications for depleted kerogens at 2.7 Ga on early Earth (Hinrichs 2002). In contrast, the lipid biosignatures of methane producing archaea are not readily identified through distinct isotopic labels and have received comparably little attention in analyses of archaea in environmental samples. Indeed, environmental analyses generally detect only free archaeal lipids, not the intact, polar molecules found in the membrane of living organisms. As part of the Ames NAI, the 'Early Microbial Ecosystem Research Group' (EMERG) is working to understand microbial processes in the hypersaline cyanobacterial mats growing in the salt evaporation ponds of the Exportadora de Sal at Guerrero Negro, Baja California Sur, Mexico. The aim of this study was to develop methods by which we could identify the organisms responsible for methane generation in this environment. While the ester-bound fatty acids, hopanoids and wax esters provide a means to identify most of the bacterial components of these mats, the archaea which Ere evidently present through genomic assays and the fact of intense methane production (Hoehler et al. 200l), have not been identified through their corresponding lipid signatures. Archaeal core lipids present a number of analytical challenges. The core lipids of methanogens comprise C20, C40 and sometimes C25 isoprenoid chains, linked through ether bonds to glycerol. As well as archaeal (C20), sn-2- and sn-3-hydroxyarchaeol are associated particularly with methylotrophic methanogens. Recently, we have also identified a dihydroxyarchaeol in a hyperthermophilic methanogen (Summons et al. 2002). Additional structural diversity is encoded into the polar head groups that are attached to the glycerol ether cores. The C20 core lipids are readily analyzed by GC-MS as their volatile trimethylsilyl derivatives while compounds with intact polar head groups can only be detected using LC-MS approaches. Our approach was to utilize the alternative of an ether cleavage reagent (BBr3 vs. HI) and a hydride reducing agent to convert all ether lipids to hydrocarbon in order to provide a vertical profile of quantitative information that might be matched to methane fluxes. We have found that while conventional acid hydrolysis and HI treatment will destroy hydroxyarchaeols, molecular information remains intact through use of BBr3 for ether cleavage. This method revealed the presence of traces of biphytane and various ether alkyls associated with some sulfate reducing bacteria within the mat structure. An interesting, and potentially valuable, byproduct of the method utilizing HI was the identification of abundant homohopanoids after superhydride reduction. Evidently present as sulfur-bound diagenetic products these hopanoids are likely cyanobacterial biomarkers in the early stages of diagenetic preservation.

Jahnke, Linda L.↗

North African savanna fires and atmospheric carbon dioxide

The effect of north African savanna fires on atmospheric CO2 is investigated using a tracer transport model. The model uses winds from operational numerical weather prediction analyses and provides CO2 concentrations as a function of space and time. After a spin-up period of several years, biomass-burning sources are added, and model experiments are run for an additional year, utilizing various estimates of CO2 sources. The various model experiments show that biomass burning in the north African savannas significantly affects CO2 concentrations in South America. The effect is more pronounced during the period from January through March, when biomass burning in South America is almost nonexistent. During this period, atmospheric CO2 concentrations in parts of South America typically may increase by 0.5 to 0.75 ppm at 970 mbar, the average pressure of the lowest model layer. These figures are above the probable uncertainty level, as model runs with biomass-burning sources estimated from independent studies using distinct data sets and techniques indicate. From May through September, when severe biomass burning occurs in South America, the effect of north African savanna fires over South America has become generally small at 970 mbar, but north of the equator it may be of the same magnitude or larger than the effect of South American fires. The CO2 concentration increase in the extreme northern and southern portions of South America, however, is mostly due to southern African fires, whose effect may be 2-3 times larger than the effect of South American fires at 970 mbar. Even in the central part of the continent, where local biomass-burning emissions are maximum, southern African fires contribute to at least 15% of the CO2 concentration increase at 970 mbar. At higher levels in the atmosphere, less CO2 emitted by north African savanna fires reaches South America, and at 100 mbar no significant amount of CO2 is transported across the Atlantic Ocean. The vertical structure of the CO2 concentration increase due to biomass burning differs substantially, depending on whether sources are local or remote. A prominent maximum of CO2 concentration increase in the lower layers characterizes the effect of local sources, whereas a more homogeneous profile of CO2 concentration increase characterizes the effect of remote sources. The results demonstrate the strong remote effects of African biomoass burning which, owing to the general circulation of the atmosphere, are felt as far away as South America.

Iacobellis, Sam F.↗

GEOS Atmospheric Model: Challenges at Exascale

The Goddard Earth Observing System (GEOS) model at NASA's Global Modeling and Assimilation Office (GMAO) is used to simulate the multi-scale variability of the Earth's weather and climate, and is used primarily to assimilate conventional and satellite-based observations for weather forecasting and reanalysis. In addition, assimilations coupled to an ocean model are used for longer-term forecasting (e.g., El Nino) on seasonal to interannual times-scales. The GMAO's research activities, including system development, focus on numerous time and space scales, as detailed on the GMAO website, where they are tabbed under five major themes: Weather Analysis and Prediction; Seasonal-Decadal Analysis and Prediction; Reanalysis; Global Mesoscale Modeling, and Observing System Science. A brief description of the GEOS systems can also be found at the GMAO website. GEOS executes as a collection of earth system components connected through the Earth System Modeling Framework (ESMF). The ESMF layer is supplemented with the MAPL (Modeling, Analysis, and Prediction Layer) software toolkit developed at the GMAO, which facilitates the organization of the computational components into a hierarchical architecture. GEOS systems run in parallel using a horizontal decomposition of the Earth's sphere into processing elements (PEs). Communication between PEs is primarily through a message passing framework, using the message passing interface (MPI), and through explicit use of node-level shared memory access via the SHMEM (Symmetric Hierarchical Memory access) protocol. Production GEOS weather prediction systems currently run at 12.5-kilometer horizontal resolution with 72 vertical levels decomposed into PEs associated with 5,400 MPI processes. Research GEOS systems run at resolutions as fine as 1.5 kilometers globally using as many as 30,000 MPI processes. Looking forward, these systems can be expected to see a 2 times increase in horizontal resolution every two to three years, as well as less frequent increases in vertical resolution. Coupling these resolution changes with increases in complexity, the computational demands on the GEOS production and research systems should easily increase 100-fold over the next five years. Currently, our 12.5 kilometer weather prediction system narrowly meets the time-to-solution demands of a near-real-time production system. Work is now in progress to take advantage of a hybrid MPI-OpenMP parallelism strategy, in an attempt to achieve a modest two-fold speed-up to accommodate an immediate demand due to increased scientific complexity and an increase in vertical resolution. Pursuing demands that require a 10- to 100-fold increases or more, however, would require a detailed exploration of the computational profile of GEOS, as well as targeted solutions using more advanced high-performance computing technologies. Increased computing demands of 100-fold will be required within five years based on anticipated changes in the GEOS production systems, increases of 1000-fold can be anticipated over the next ten years.

ESMF↗

Magellan stereo images and Venusian geology

Areas of Venus imaged by Magellan radar with multiple viewing conditions provide unique data that will contribute to the solution of venusian geologic problems and provide a basis for quantitative comparison of venusian landforms with those on other planetary bodies. Three sets of images with different viewing conditions have been acquired: (1) left-looking with variable incidence angles (cycle 1 profile), (2) right-looking with nearly constant incidence angles (cycle 2 profile), and (3) left-looking with variable incidence angles that are almost always smaller than those in (1) (cycle 3 profiles). The unique data provided by paired images of the same scene with different incidence angles arises from image displacements caused by the relief of individual landforms at scales comparable to the ground-range and azimuth resolutions of the images. There are two aspects of the data: (1) Stereopsis achieved by simultaneous viewing of paired left-looking images of the same scene permits three-dimensional perception and interpretation of the morphologies of landforms at resolutions much finer than the altimetry footprints. (2) Measurements of differences of image displacements (parallax) on paired images with known imaging geometries provide quantitative estimates of the relief and shapes of landforms. The potential scientific contributions of the data can be grouped into two interrelated classes: (A) geologic mapping, analysis, and interpretation and (B) topical studies that involve topographic measurements. Stereopsis, without quantitative measurements, enhances geologic mapping, analysis, and interpretation of the rock units of Venus to a degree that cannot be overestimated. In geologic mapping, assemblages of landforms, assessments of backscatter and variations in backscatter, and fine-scale topography are used to define and characterize geologic map units that represent laterally continuous deposits or rock units. Stereopsis adds the important dimension of local relief for characterization of geologic units at a scale that is not possible with Magellan altimetry or products derived from the altimetry. Relative ages of the geologic units are determined using the well-known principles of superposition and intersection. Here, the perception of relief is invaluable because superposition relations among the geological units are more readily and clearly established. The recognition of folds, faults, and fault systems, regardless of their orientations, is facilitated with stereopsis so that sequences of deformation of the geologic units can be determined and structural analyses vastly improved. Shapes of landforms are readily perceived so that they can be properly interpreted. The end result of the mapping, analyses, and interpretations is a geologic history of Venus that includes the sequences of formation and deformation of various geologic units. Measurements of relief at the finest scale possible are necessary for numerous topical studies. Standard altimetry will provide the necessary information on the relief of most large landforms, but it tends to underestimate the relief of small landforms and distorts their shapes. Although special processing of the altimeter echoes improves the estimates of the relief and shapes of some landforms, there are uncertainties in the interpretations of the echoes. Examples of topical studies requiring measurements of relief are given.

Moore, H. J.↗

Low-shear modeled microgravity: a global environmental regulatory signal affecting bacterial gene expression, physiology, and pathogenesis

Bacteria inhabit an impressive variety of ecological niches and must adapt constantly to changing environmental conditions. While numerous environmental signals have been examined for their effect on bacteria, the effects of mechanical forces such as shear stress and gravity have only been investigated to a limited extent. However, several important studies have demonstrated a key role for the environmental signals of low shear and/or microgravity in the regulation of bacterial gene expression, physiology, and pathogenesis [Chem. Rec. 1 (2001) 333; Appl. Microbiol. Biotechnol. 54 (2000) 33; Appl. Environ. Microbiol. 63 (1997) 4090; J. Ind. Microbiol. 18 (1997) 22; Curr. Microbiol. 34(4) (1997) 199; Appl. Microbiol. Biotechnol. 56(3-4) (2001) 384; Infect Immun. 68(6) (2000) 3147; Cell 109(7) (2002) 913; Appl. Environ. Microbiol. 68(11) (2002) 5408; Proc. Natl. Acad. Sci. U. S. A. 99(21) (2002) 13807]. The response of bacteria to these environmental signals, which are similar to those encountered during prokaryotic life cycles, may provide insight into bacterial adaptations to physiologically relevant conditions. This review focuses on the current and potential future research trends aimed at understanding the effect of the mechanical forces of low shear and microgravity analogues on different bacterial parameters. In addition, this review also discusses the use of microgravity technology to generate physiologically relevant human tissue models for research in bacterial pathogenesis.

NASA Discipline Cell Biology↗

Quantitative Species Measurements in Microgravity Combustion Flames using Near-Infrared Diode Lasers

Understanding the physical phenomena controlling the ignition and spread of flames in microgravity has importance for space safety as well as for characterizing dynamical and chemical combustion processes which are normally masked by buoyancy and other gravity-related effects. Unfortunately, combustion is highly complicated by fluid mechanical and chemical kinetic processes, requiring the use of numerical modeling to compare with carefully designed experiments. More sophisticated diagnostic methods are needed to provide the kind of quantitative data necessary to characterize the properties of microgravity combustion as well as provide accurate feedback to improve the predictive capabilities of the models. Diode lasers are a natural choice for use under the severe conditions of low gravity experiments. Reliable, simple solid state operation at low power satisfies the operational restrictions imposed by drop towers, aircraft and space-based studies. Modulation wavelength absorption spectroscopy (WMS) provides a means to make highly sensitive and quantitative measurements of local gas concentration and, in certain cases, temperature. With near-infrared diode lasers, detection of virtually all major combustion species with extremely rapid response time is possible in an inexpensive package. Advancements in near-infrared diode laser fabrication technology and concurrent development of optical fibers for these lasers led to their use in drop towers. Since near-infrared absorption line strengths for overtone and combination vibrational transitions are weaker than the mid-infrared fundamental bands, WMS techniques are applied to increase detection sensitivity and allow measurement of the major combustion gases. In the first microgravity species measurement, Silver et al. mounted a fiber-coupled laser at the top of the NASA 2.2-sec drop tower and piped the light through a single-mode fiber to the drop rig. A fiber splitter divided the light into eight channels that directed the laser beam across a methane or propane diffusion jet flame. The light beams were recaptured by a set of gradient index lenses, coupled back into separate fiber optic lines, and transmitted back to detectors and electronics in the instrument package. In these experiments a 6-mm od fiber cable (containing the nine optical fibers) fell with the drop rig. Using separate detection and demodulation channels, spatial and temporal (up to 20 Hz) maps of water vapor and methane concentrations were obtained at differing heights in the flames. While this apparatus was useful from a demonstration standpoint, several drawbacks needed attention before useful scientific measurements could be obtained. First, eight lines of sight are somewhat insufficient for detailing the spatial profiles of the gas. Second, multiple detection channels operating in parallel are both expensive and present a challenge for accurate calibration. As a result, a newer scanning system was developed in our first contract under this program. The primary characteristic of this system is that it contains a single detection channel and achieves "continuous" spatial resolution by scanning the laser beam across the flame region, then directing this beam onto a single detector. Thus spatial measurements are converted to a temporal series of data. The true spatial resolution is limited only by the beam diameter and width of the sweep. In these experiments the beam is focused to about 1-mm diameter and scans across a region up to 4-cm wide.

Silver, Joel A.↗