Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “trajectory inference”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

Chapter 13. Atmospheric Dynamics and Meteorology

Titan, after Venus, is the second example in the solar system of an atmosphere with a global cyclostrophic circulation, but in this case a circulation that has a strong seasonal modulation in the middle atmosphere. Direct measurement of Titan's winds, particularly observations tracking the Huygens probe at 10 deg S, indicate that the zonal winds are mostly in the sense of the satellite's rotation. They generally increase with altitude and become cyclostrophic near 35 km above the surface. An exception to this is a sharp minimum centered near 75 km, where the wind velocity decreases to nearly zero. Zonal winds derived from temperatures retrieved from Cassini orbiter measurements, using the thermal wind equation, indicate a strong winter circumpolar vortex, with maximum winds of 190 m/s at mid northern latitudes near 300 km. Above this level, the vortex decays. Curiously, the stratospheric zonal winds and temperatures in both hemispheres are symmetric about a pole that is offset from the surface pole by about 4 deg. The cause of this is not well understood, but it may reflect the response of a cyclostrophic circulation to the onset between the equator, where the distance to the rotation axis is greatest, and the seasonally varying subsolar latitude. The mean meridional circulation can be inferred from the temperature field and the meridional distribution of organic molecules and condensates and hazes. Both the warm temperatures near 400 km and the enhanced concentration of several organic molecules suggest subsidence in the north polar region during winter and early spring. Stratospheric condensates are localized at high northern latitudes, with a sharp cut-off near 50 deg N. Titan's winter polar vortex appears to share many of the same characteristics of isolating high and low-latitude air masses as do the winter polar vortices on Earth that envelop the ozone holes. Global mapping of temperatures, winds, and composition in the troposphere, by contrast, is incomplete. The few suitable discrete clouds that have been found for tracking indicate smaller velocities than aloft, consistent: with the Huygens measurements, Along the descent trajectory, the Huygens measurements indicate eastward zonal winds down to 7 km, where they shift westward, and then eastward again below 1 km dawn to the surface. The low-latitude dune fields seen in Cassini RADAR images have been interpreted as longitudinal dunes occurring in a mean eastward zonal wind. This is not like Earth, where the low-latitude winds are westward above the surface. Because the net zonal-mean time-averaged torque exerted by the surface on the atmosphere should vanish, there must be westward flow over part of the surface; the question is where and when. The meridional contrast in tropospheric temperatures deduced from radio occultations at low, mid, and high latitudes. is small, approximately 5 K at the tropopause and approximately 3 K at the surface. This implies efficient heat transport, probably by axisymmetric meridional circulations. The effect of the methane "hydrological" cycle on the atmospheric circulation is not well constrained by existing measurements, Understanding the mature of the surface-atmosphere coupling will be critical to elucidating the atmospheric transports of momentum, heat, and volatiles.

Flasar, F. M.↗
NASA NTRS2012

The Structure and Emission Model of the Relativistic Jet in the Quasar 3C 279 Inferred From Radio To High-Energy Gamma-Ray Observations in 2008-2010

We present time-resolved broad-band observations of the quasar 3C 279 obtained from multiwavelength campaigns conducted during the first two years of the Fermi Gamma-ray Space Telescope mission. While investigating the previously reported gamma-ray/optical flare accompanied by a change in optical polarization, we found that the optical emission appears delayed with respect to the gamma-ray emission by about 10 days. X-ray observations reveal a pair of 'isolated' flares separated. by approx. 90 days, with only weak gamma-ray/optical counterparts. The spectral structure measured by Spitzer reveals a synchrotron component peaking in the mid-infrared band with a sharp break at the far-infrared band during the gamma-ray flare, while the peak appears in the mm/sub-mm band in the low state. Selected spectral energy distributions are fitted with leptonic models including Comptonization of external radiation produced in a dusty torus or the broad-line region. Adopting the interpretation of the polarization swing involving propagation of the emitting region along a curved trajectory, we can explain the evolution of the broad-band spectra during the gamma-ray flaring event by a shift of its location from approx. 1 pc to approx. 4 pc from the central black hole. On the other hand, if the gamma-ray flare is generated instead at sub-pc distance from the central black hole, the far-infrared break can be explained by synchrotron self-absorption. We also model the low spectral state, dominated by the mm/sub-mm peaking synchrotron component, and suggest that the corresponding inverse-Compton component explains the steady X-ray emission.

Source record↗
DOE OSTI2026

Artificial Intelligence for Event Reconstruction and Higgs Physics at CMS and Future Colliders

This dissertation charts a trajectory in which advances in artificial intelligence (AI) play a central role in pushing the high-energy physics frontier, complementing progress driven by higher collision energies and larger colliders. The discovery potential of the LHC and future colliders relies on accurate reconstruction of increasingly complex particle collision events. In the CMS experiment, this task is performed by the particle-flow (PF) algorithm. This dissertation presents the first implementation of a machine-learning-based particle-flow (MLPF) reconstruction in the CMS detector based on transformer architectures. In simulated top quark--antiquark pair (ttbar) events under LHC Run~3 (2023--2024) conditions, MLPF improves jet energy resolution by 10--20\% compared to standard PF for jets with transverse momentum between 30--100\GeV. Runtime performance is evaluated using simulated multijet events, with a median inference time of 20\unit{ms} per event on an NVIDIA L4 GPU, compa red to approximately 110\unit{ms} for standard PF. The MLPF algorithm is also validated on Run~3 collision data, representing the first data-validated ML-based reconstruction pipeline at any LHC experiment. We then extend MLPF toward future electron--positron colliders and introduce the first full-simulation cross-detector transfer learning workflow for PF reconstruction. The model is pre-trained on simulated events from the Compact Linear Collider detector (CLICdet) and fine-tuned on the CLIC-like detector (CLD) proposed for the Future Circular Collider (FCC). This approach achieves up to a 40\% improvement in jet energy resolution over rule-based reconstruction while reducing the required training dataset size by an order of magnitude, demonstrating the potential of AI to accelerate detector development and optimization. This dissertation also demonstrates how modern AI techniques enhance the sensitivity of LHC physics analyses. A CMS search for highly Lorentz-boosted Higgs bosons decaying to \textrm{W} boson pairs is presented, focusing on the single-lepton final state. A dedicated fine-tuning strategy for \ParT yields an approximately 70\% increase in expected sensitivity relative to the baseline model. The analysis uses proton--proton collision data at a center-of-mass energy of \ensuremath{\sqrt{s}=13\TeV} collected by CMS between 2016 and 2018, corresponding to an integrated luminosity of 138\ensuremath{\ \mathrm{fb}^{-1}}. The expected significance of the search is $1.86\sigma$, with an observed signal strength of $-0.19^{+0.48}_{-0.46}$. Finally, explainable AI techniques are applied to the MLPF and \ParticleNet algorithms using layerwise relevance propagation, showing that both models base their predictions on physically meaningful features consistent with our physics intuition. Together, these results demonstrate how advanced AI methods can enhance reconstruction, analysis sensitivity, and interpretability, shaping the next era of experimental parti cle physics.

Mokhtar, Farouk [UC, San Diego]↗
NASA NTRS1999

Processes Affecting Tropospheric Ozone Inferred from Ozonesonde and Other Tracer Data from the R/V R H Brown Atlantic Cruise (37N-34S) in January-February 1999

During the Aerosols-99 trans-Atlantic cruise from Norfolk, VA, to Cape Town, South Africa, 22 ozonesondes were launched from the NOAA R/V R H Brown between 17 Jan and 6 Feb 1999, with all sondes but one reaching 30 km. A composite of ozone profiles along the transect shows high free tropospheric ozone (up to 100 ppbv at 9 km) between 5N and 20S, a coherent feature straddling either side of the ITCZ. Latitudinal variations of tropospheric ozone are interpreted using correlative measurements of surface ozone, CO, water vapor, and aerosol optical thickness (column absorbance) measured from the ship. Elevated ozone in the lower troposphere results from photochemical reactions of precursors emitted by biomass burning north of the ITCZ. However, the greatest ozone mixing ratios are in the mid-troposphere south of the ITCZ, which gives evidence of interhemispheric transport. Column-integrated tropospheric ozone, 35 DU from 0-16 km, agrees with that derived from the TOMS satellite by the modified-residual method [Thompson and Hudson, 1999]. NCEP wind fields, ship-launched radiosondes and back trajectories are consistent with a picture of recirculating air parcels centered in the tropical Atlantic region which is identified with the maximum wave-one amplitude in total ozone seen in sondes and by satellite.

Thompson, Anne M.↗
NASA NTRS2015

HSRL-2 Observations of Aerosol Variability During an Aerosol Build-up Event in Houston and Comparisons With WRF-Chem

The NASA Langley airborne multi-wavelength High Spectral Resolution Lidar (HSRL-2) provides vertical distribution of aerosol optical properties as curtains of aerosol extinction, backscatter and depolarization along the flight track, plus intensive properties that are used to infer aerosol type and external mixing of types. Deployed aboard the NASA Langley King Air on the DISCOVER-AQ field mission in Houston in September 2013, HSRL-2 flew a pattern that included 18 ground sites, repeated four times a day, coordinated with a suite of airborne in situ measurements. The horizontally and vertically resolved curtains of HSRL-2 measurements give an unparalleled view of the spatial and temporal variability of aerosol, which provide broad context for interpreting other measurements and models. Detailed comparisons of aerosol extinction are made with the WRF-Chem chemical transport model along the HSRL-2 flight path. The period from Sept. 11-14 is notable for a large aerosol build-up and persistent smoke layers. We investigate the aerosol properties using the vertically resolved HSRL-2 measurements and aerosol typing analysis plus WRFChem model tracers and back trajectories, and modeling of humidification effects.

Burton, Sharon P.↗
DOE OSTI2022

Locating fast-varying line disturbances with the frequency mismatch

In an attempt to provide an efficient method for line disturbance identification in complex networks of diffusively coupled agents, we recently proposed to leverage the frequency mismatch. The frequency mismatch filters out the intricate combination of interactions induced by the network structure and quantifies to what extent the trajectory of each agent is affected by the disturbance. In this previous work, we provided some analytical evidence of its efficiency when the perturbation is assumed to be slow. Here we claim that the frequency mismatch performs actually well for most disturbance regimes. This is shown through a series of simulations and is backed up by an analytical argument. Therefore, we argue that the frequency mismatch is an efficient and elegant tool for line disturbance location in complex networks of diffusively coupled agents.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2010

Plasma IMS Composition Measurements for Europa and the Other Galilean Moons

NASA and ESA are planning the joint Europa Jupiter System Mission (EJSM) to the Jupiter system with specific emphasis to Europa and Ganymede, respectively. The Japanese Space Agency is also planning an orbiter mission to explore Jupiter's magnetosphere and the Galilean satellites. For NASA's Jupiter Europa Orbiter (JEO) we are developing the 3D Ion Mass Spectrometer (IMS) with two main goals which can also be applied to the other Galilean moons, 1) measure the plasma interaction between Europa and Jupiter's magnetosphere and 2) infer the 4 pi surface composition to trace elemental and significant isotopic levels. The first goal supports the magnetometer (MAG) measurements, primarily directed at detection of Europa's sub-surface ocean, while the second gives information about transfer of material between the Galilean moons, and between the moon surfaces and subsurface layers putatively including oceans. The measurement of the interactions for all the Galilean moons can be used to trace the in situ ion measurements of pickup ions back to either Europa's or Ganymede's surface from the respectively orbiting spacecraft. The IMS instrument, being developed under NASA's Astrobiology Instrument Development Program, would maximally achieve plasma measurement requirements for JEO and EJSM while moving forward our knowledge of Jupiter system composition and source processes to far higher levels than previously envisaged. The composition of the global surfaces of Europa and Ganymede can be inferred from the measurement of ejected neutrals and pick-up ions using at minimum an in situ payload including MAG and IMS also fully capable of meeting Level 1 mission requirements for ocean detection and survey. Elemental and isotopic analysis of potentially extruded oceanic materials at the moon surfaces would further support the ocean objectives. These measurements should be made from a polar orbiting spacecraft about Europa or Ganymede at height 100 km. The ejecta produced by sputtering of the surfaces of Europa and Ganymede has been shown to be representative of the surface composition. Level 2 science on surface geology and composition can then be further enhanced by addition of the following: 3D Ion Neutral Mass Spectrometer (INNS), 3D plasma electron spectrometer (ELS), and hot plasma energetic particle instrument. The measurement approach is to alternate between times measuring pickup ions and times measuring plasma and magnetic field parameters along the spacecraft trajectory. By measuring the pickup ion energy, arrival direction and mass-per-charge, the ion can be traced back along the ejection trajectory to the approximate area of origin if the 3-D electric field and magnetic field are known. In situ observations of plasma flow velocities and vector magnetic fields can be used to determine the local convective electric field (E = -VXB) along the spacecraft trajectory. By combining this information with models of the magnetospheric interaction with Europa, one can generate 3D maps of the electric and magnetic field and compute the trajectories of the pickup ions back to the surface or exospheric points of origin. In the case of Ganymede there is the additional complexity of its own internal dipole magnetic field, while Io's volcanic activity introduces the complexity of a highly structured denser atmosphere. Callisto with its less globally extended exosphere will have a simpler interaction than for Europa (i.e., more like our moon). We will discuss these differences in light of the above proposed technique. Finally, the INNS observations and neutral exosphere models are needed to estimate production rates of pickup ions. The hot plasma measurements are needed to correct for sputtering rates which can be time dependent and electron plasma observations for electron impact ionization rates. Instrument characteristics, field-of-view requirements, modes of operation and effects of radiation on instrument functionality will be discussed.

Sittler, Edward↗
DOE OSTI2026

Data-driven particle dynamics: Structure-preserving coarse-graining for emergent behavior in non-equilibrium systems

Multiscale systems are ubiquitous in science and technology, but are notoriously challenging to simulate as short spatiotemporal scales must be appropriately linked to emergent bulk physics. When expensive high-dimensional dynamical systems are coarse-grained into low-dimensional models, the entropic loss of information leads to emergent physics which are dissipative, history-dependent, and stochastic. To machine learn coarse-grained dynamics from time-series observations of particle trajectories, we propose a framework using the metriplectic bracket formalism that preserves these properties by construction; most notably, the framework guarantees discrete notions of the first and second laws of thermodynamics, conservation of momentum, and a discrete fluctuation-dissipation balance crucial for capturing non-equilibrium statistics. We introduce the mathematical framework abstractly before specializing to a particle discretization. As labels are generally unavailable for entropic state variables, we introduce a novel self-supervised learning strategy to identify emergent structural variables. We validate the method on benchmark systems and demonstrate its utility on two challenging examples: (1) coarse-graining star polymers at challenging levels of coarse-graining while preserving non-equilibrium statistics, and (2) learning models from high-speed video of colloidal suspensions that capture coupling between local rearrangement events and emergent stochastic dynamics. We provide open-source implementations in both PyTorch and LAMMPS, enabling large-scale inference and extensibility to diverse particle-based systems.

Computational Engineering, Finance, and Science (c↗
NASA NTRS1983

Nuclear cascades in Saturn's rings - Cosmic ray albedo neutron decay and origins of trapped protons in the inner magnetosphere

The nearly equatorial trajectory of the Pioneer 11 spacecraft through Saturn's high energy proton radiation belts and under the main A-B-C rings provided a unique opportunity to study the radial dependence of the greater than 30 MeV proton intensities in the belts in terms of models for secondary nucleon production by cosmic ray interactions in the rings, in situ proton injection in the radiation belts by neutron beta decay, magnetospheric diffusion, and absorption by planetary rings and satellites. Maximum trapped proton intensities measured by Pioneer 11 in the radiation belts are compared with calculated intensities and found consistent with trapping times of roughly 40 years and a radial diffusion coefficient of about 10 to the -15th L to the 9th R sub s squared/s. Differential energy spectra proportional to E to the -2 estimated from integral measurements of trapped photons with E greater than 100 MeV are consistent with the beta decay model, but an inferred turndown of the spectra toward lower energies and reported integral proton anisotropies of a specified form both indicate the need for more realistic calculations of the neutron source from the rings and the radiation belt loss processes.

Cooper, J. F.↗
NASA NTRS2011

Plasma IMS Composition Measurements for Europa, Ganymede, and the Jovian System

NASA and ESA are now planning a reduced version of the joint Europa Jupiter System Mission (EJSM), potentially including a radically descoped Jupiter Europa Orbiter (JEO) but still with magnetometer and plasma instruments. Similar field and plasma instrumentation would also reside on ESA's Jupiter Ganymede Orbiter (JGO), which conceivably could carry out multiple flybys of Europa before entering orbit at Ganymede. We are developing the 3D Ion Mass Spectrometer (IMS) designed to measure both major and minor ion species within the high radiation environment of Jupiter's magnetosphere and the icy Galilean moons. The IMS covers the energy range from 10 eV to 30 keY, wide field-of-view (FOV) capability and 10-60 sec time resolution for major ions. This instrument has two main goals: 1) measure the plasma interaction between Europa and Jupiter's magnetosphere and 2) infer the global surface composition to trace elemental and significant isotopic levels; these goals are also applicable for in-situ measurements at Ganymede and Callisto, and remotely everywhere via the iogenic plasma for 10. The first goal supports the magnetometer (MAG) measurements, primarily directed at detection of Europa's sub-surface ocean, while the second goal gives information about transfer of material between the Galilean moons, e.g. mainly from 10 to the other moons, and further allows detection of oceanic materials emergent to the moon surfaces from subsurface layers putatively including salt water oceans. Outgassed exospheric materials are probed by the IMS by measuring pickup ions accelerated up to spacecraft altitudes of approximately 100-200 km in electric fields extending through the local magnetospheric environment and moon exosphere to the surface. Our 3D hybrid kinetic model of the moon-magnetosphere interaction is used to construct a global model of electric and magnetic fields for tracing of pickup ion trajectories back to the sources at approximate surface resolution of 100 km. We show that Europa's exospheric ionosphere is dominated by pickup ions with energies of 100-1000 eV. We also expect field aligned polar ion outflows driven by ionospheric electrons via the polarization electric field at Europa; the IMS will observe such outflows and thus sample the ionosphere below spacecraft orbit altitude approximately 100 km. Based on previous Ganymede studies, we also comment on IMS applications to a Ganymede orbiter. The IMS and the Europa interaction model are respectively being developed with support from NASA's Astrobiology Instrument Development (ASTID) and Outer Planets Research (OPR) programs.

Sittler, E. C., Jr.↗
NASA NTRS2011

Plasma IMS Composition Measurements for Europa, Ganymede, and the Jovian Systems

NASA and ESA are now planning a reduced version of the joint Europa Jupiter System Mission (EJSM), potentially including a radically descoped Jupiter Europa Orbiter (JEO) but still with magnetometer and plasma instruments. Similar field and plasma instrumentation would also reside on ESA's Jupiter Ganymede Orbiter (JGO), which conceivably could carry out multiple flybys of Europa before entering orbit at Ganymede. We are developing the 3D Ion Mass Spectrometer (IMS) designed to measure both major and minor ion species within the high radiation environment of Jupiter s magnetosphere and the icy Galilean moons. The IMS covers the energy range from 10 eV to 30 keV, wide field-ofview (FOV) capability and 10-60 sec time resolution for major ions. This instrument has two main goals: 1) measure the plasma interaction between Europa and Jupiter s magnetosphere and 2) infer the global surface composition to trace elemental and significant isotopic levels; these goals are also applicable for in-situ measurements at Ganymede and Callisto, and remotely everywhere via the iogenic plasma for Io. The first goal supports the magnetometer (MAG) measurements, primarily directed at detection of Europa's sub-surface ocean, while the second goal gives information about transfer of material between the Galilean moons, e.g. mainly from Io to the other moons, and further allows detection of oceanic materials emergent to the moon surfaces from subsurface layers putatively including salt water oceans. Outgassed exospheric materials are probed by the IMS by measuring pickup ions accelerated up to spacecraft altitudes of approximately 100-200 km in electric fields extending through the local magnetospheric environment and moon exosphere to the surface. Our 3D hybrid kinetic model of the moon-magnetosphere interaction is used to construct a global model of electric and magnetic fields for tracing of pickup ion trajectories back to the sources at approximate surface resolution of 100 km. We show that Europa's exospheric ionosphere is dominated by pickup ions with energies of 100-1000 eV. We also expect field aligned polar ion outflows driven by ionospheric electrons via the polarization electric field at Europa; the IMS will observe such outflows and thus sample the ionosphere below spacecraft orbit altitude approximately 100 km. Based on previous Ganymede studies, we also comment on IMS applications to a Ganymede orbiter. The IMS and the Europa interaction model are respectively being developed with support from NASA's Astrobiology Instrument Development (ASTID) and Outer Planets Research (OPR) programs.

Sittler, E.↗
NASA NTRS2016

Using Natural Language to Enhance Mission Effectiveness

The availability of highly capable, yet relatively cheap, unmanned aerial vehicles (UAVs) is opening up new areas of use for hobbyists and for professional-related activities. The driving function of this research is allowing a non-UAV pilot, an operator, to define and manage a mission. This paper describes the preliminary usability measures of an interface that allows an operator to define the mission using speech to make inputs. An experiment was conducted to begin to enumerate the efficacy and user acceptance of using voice commands to define a multi-UAV mission and to provide high-level vehicle control commands such as "takeoff." The primary independent variable was input type - voice or mouse. The primary dependent variables consisted of the correctness of the mission parameter inputs and the time needed to make all inputs. Other dependent variables included NASA-TLX workload ratings and subjective ratings on a final questionnaire. The experiment required each subject to fill in an online form that contained comparable required information that would be needed for a package dispatcher to deliver packages. For each run, subjects typed in a simple numeric code for the package code. They then defined the initial starting position, the delivery location, and the return location using either pull-down menus or voice input. Voice input was accomplished using CMU Sphinx4-5prealpha for speech recognition. They then inputted the length of the package. These were the option fields. The subject had the system "Calculate Trajectory" and then "Takeoff" once the trajectory was calculated. Later, the subject used "Land" to finish the run. After the voice and mouse input blocked runs, subjects completed a NASA-TLX. At the conclusion of all runs, subjects completed a questionnaire asking them about their experience in inputting the mission parameters, and starting and stopping the mission using mouse and voice input. In general, the usability of voice commands is acceptable. With a relatively well-defined and simple vocabulary, the operator can input the vast majority of the mission parameters using simple, intuitive voice commands. However, voice input may be more applicable to initial mission specification rather than for critical commands such as the need to land immediately due to time and feedback constraints. It would also be convenient to retrieve relevant mission information using voice input. Therefore, further on-going research is looking at using intent from operator utterances to provide the relevant mission information to the operator. The information displayed will be inferred from the operator's utterances just before key phrases are spoken. Linguistic analysis of the context of verbal communication provides insight into the intended meaning of commonly heard phrases such as "What's it doing now?" Analyzing the semantic sphere surrounding these common phrases enables us to predict the operator's intent and supply the operator's desired information to the interface. This paper also describes preliminary investigations into the generation of the semantic space of UAV operation and the success at providing information to the interface based on the operator's utterances.

Trujillo, Anna C.↗
NASA NTRS2012

Adaptive Sampling of Time Series During Remote Exploration

This work deals with the challenge of online adaptive data collection in a time series. A remote sensor or explorer agent adapts its rate of data collection in order to track anomalous events while obeying constraints on time and power. This problem is challenging because the agent has limited visibility (all its datapoints lie in the past) and limited control (it can only decide when to collect its next datapoint). This problem is treated from an information-theoretic perspective, fitting a probabilistic model to collected data and optimizing the future sampling strategy to maximize information gain. The performance characteristics of stationary and nonstationary Gaussian process models are compared. Self-throttling sensors could benefit environmental sensor networks and monitoring as well as robotic exploration. Explorer agents can improve performance by adjusting their data collection rate, preserving scarce power or bandwidth resources during uninteresting times while fully covering anomalous events of interest. For example, a remote earthquake sensor could conserve power by limiting its measurements during normal conditions and increasing its cadence during rare earthquake events. A similar capability could improve sensor platforms traversing a fixed trajectory, such as an exploration rover transect or a deep space flyby. These agents can adapt observation times to improve sample coverage during moments of rapid change. An adaptive sampling approach couples sensor autonomy, instrument interpretation, and sampling. The challenge is addressed as an active learning problem, which already has extensive theoretical treatment in the statistics and machine learning literature. A statistical Gaussian process (GP) model is employed to guide sample decisions that maximize information gain. Nonsta tion - ary (e.g., time-varying) covariance relationships permit the system to represent and track local anomalies, in contrast with current GP approaches. Most common GP models are stationary, e.g., the covariance relationships are time-invariant. In such cases, information gain is independent of previously collected data, and the optimal solution can always be computed in advance. Information-optimal sampling of a stationary GP time series thus reduces to even spacing, and such models are not appropriate for tracking localized anomalies. Additionally, GP model inference can be computationally expensive.

Thompson, David R.↗
DOE OSTI2021

Topographical Controls on Hillslope-Scale Hydrology Drive Shrub Distributions on the Seward Peninsula, Alaska

Observations indicate shrubs are expanding across the Arctic tundra, mainly on hillslopes and primarily in response to climate warming. However, the impact topography exerts on hydrology, nutrient dynamics, and plant growth can make untangling the mechanisms behind shrub expansion difficult. We examined the role topography plays in determining shrub expansion by applying a coupled transect version of a mechanistic ecosystem model (ecosys) in a tundra hillslope site in the Seward Peninsula, Alaska. Modeled biomass of the dominant plant functional types agreed well with field measurements (R 2 = 0.89) and accurately represented shrub expansion over the past 30 years inferred from satellite observations. In the well-drained crest position, canopy water potential and plant nitrogen (N) uptake was modeled to be low from plant and microbial water stress. Intermediate soil water content in the mid-slope position enhanced mineralization and plant N uptake, increasing shrub biomass. The deciduous shrub growth in the mid-slope position was further enhanced by symbiotic N 2 fixation primed by increased root carbon allocation. The gentle slope in the poorly drained lower-slope position resulted in saturated soil conditions that reduced soil O 2 concentrations, leading to lower root O 2 uptake and lower nutrient uptake and plant biomass. Additionally, a simulation that removed topographical interconnectivity between grid cells resulted in (1) a 28% underestimate of mean shrub biomass and (2) over or underestimated shrub productivity at the various hillslope positions. Our results indicate that land models need to account for hillslope-scale coupled surface and subsurface hydrology to accurately predict current plant distributions and future trajectories in Arctic ecosystems.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2003

The Livingstone Model of a Main Propulsion System

Livingstone is a discrete, propositional logic-based inference engine that has been used for diagnosis of physical systems. We present a component-based model of a Main Propulsion System (MPS) and say how it is used with Livingstone (L2) in order to implement a diagnostic system for integrated vehicle health management (IVHM) for the Propulsion IVHM Technology Experiment (PITEX). We start by discussing the process of conceptualizing such a model. We describe graphical tools that facilitated the generation of the model. The model is composed of components (which map onto physical components), connections between components and constraints. A component is specified by variables, with a set of discrete, qualitative values for each variable in its local nominal and failure modes. For each mode, the model specifies the component's behavior and transitions. We describe the MPS components' nominal and fault modes and associated Livingstone variables and data structures. Given this model, and observed external commands and observations from the system, Livingstone tracks the state of the MPS over discrete time-steps by choosing trajectories that are consistent with observations. We briefly discuss how the compiled model fits into the overall PITEX architecture. Finally we summarize our modeling experience, discuss advantages and disadvantages of our approach, and suggest enhancements to the modeling process.

Bajwa, Anupa↗
NASA NTRS2003

Analytic Solution to the Problem of Aircraft Electric Field Mill Calibration

It is by no means a simple task to retrieve storm electric fields from an aircraft instrumented with electric field mill sensors. The presence of the aircraft distorts the ambient field in a complicated way. Before retrievals of the storm field can be made, the field mill measurement system must be "calibrated". In other words, a relationship between impressed (i.e., ambient) electric field and mill output must be established. If this relationship can be determined, it is mathematically inverted so that ambient field can be inferred from the mill outputs. Previous studies have primarily focused on linear theories where the relationship between ambient field and mill output is described by a "calibration matrix" M. Each element of the matrix describes how a particular component of the ambient field is enhanced by the aircraft. For example the product M(sub ix), E(sub x), is the contribution of the E(sub x) field to the i(th) mill output. Similarly, net aircraft charge (described by a "charge field component" E(sub q)) contributes an amount M(sub iq)E(sub q) to the output of the i(th) sensor. The central difficulty in obtaining M stems from the fact that the impressed field (E(sub x), E(sub y), E(sub z), E(sub q) is not known but is instead estimated. Typically, the aircraft is flown through a series of roll and pitch maneuvers in fair weather, and the values of the fair weather field and aircraft charge are estimated at each point along the aircraft trajectory. These initial estimates are often highly inadequate, but several investigators have improved the estimates by implementing various (ad hoc) iterative methods. Unfortunately, none of the iterative methods guarantee absolute convergence to correct values (i.e., absolute convergence to correct values has not been rigorously proven). In this work, the mathematical problem is solved directly by analytic means. For m mills installed on an arbitrary aircraft, it is shown that it is possible to solve for a single 2m-vector that provides all other needed variables (i.e., the unknown fair weather field, the unknown aircraft charge, and the unknown matrix M). Numerical tests of the solution, effects of measurement errors, and studies of solution non-uniqueness are ongoing as of this writing.

Koshak, William↗
DOE OSTI2021

Measurement of Volatile Compounds for Real-Time Analysis of Soil Microbial Metabolic Response to Simulated Snowmelt

Snowmelt dynamics are a significant determinant of microbial metabolism in soil and regulate global biogeochemical cycles of carbon and nutrients by creating seasonal variations in soil redox and nutrient pools. With an increasing concern that climate change accelerates both snowmelt timing and rate, obtaining an accurate characterization of microbial response to snowmelt is important for understanding biogeochemical cycles intertwined with soil. However, observing microbial metabolism and its dynamics non-destructively remains a major challenge for systems such as soil. Microbial volatile compounds (mVCs) emitted from soil represent information-dense signatures and when assayed non-destructively using state-of-the-art instrumentation such as Proton Transfer Reaction-Time of Flight-Mass Spectrometry (PTR-TOF-MS) provide time resolved insights into the metabolism of active microbiomes. In this study, we used PTR-TOF-MS to investigate the metabolic trajectory of microbiomes from a subalpine forest soil, and their response to a simulated wet-up event akin to snowmelt. Using an information theory approach based on the partitioning of mutual information, we identified mVC metabolite pairs with robust interactions, including those that were non-linear and with time lags. The biological context for these mVC interactions was evaluated by projecting the connections onto the Kyoto Encyclopedia of Genes and Genomes (KEGG) network of known metabolic pathways. Simulated snowmelt resulted in a rapid increase in the production of trimethylamine (TMA) suggesting that anaerobic degradation of quaternary amine osmo/cryoprotectants, such as glycine betaine, may be important contributors to this resource pulse. Unique and synergistic connections between intermediates of methylotrophic pathways such as dimethylamine, formaldehyde and methanol were observed upon wet-up and indicate that the initial pulse of TMA was likely transformed into these intermediates by methylotrophs. Increases in ammonia oxidation signatures (transformation of hydroxylamine to nitrite) were observed in parallel, and while the relative role of nitrifiers or methylotrophs cannot be confirmed, the inferred connection to TMA oxidation suggests either a direct or indirect coupling between these processes. Overall, it appears that such mVC time-series from PTR-TOF-MS combined with causal inference represents an attractive approach to non-destructively observe soil microbial metabolism and its response to environmental perturbation.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2006

On the Rates of Coronal Mass Ejections: Remote Solar and In Situ Observations

We compare the rates of coronal mass ejections (CMEs) as inferred from remote solar observations and interplanetary CMEs (ICMEs) as inferred from in situ observations at both 1 AU and Ulyssses from 1996 through 2004. We also distinguish between those ICMEs that contain a magnetic cloud (MC) and those that do not. While the rates of CMEs and ICMEs track each other well at solar minimum, they diverge significantly in early 1998, during the ascending phase of the solar cycle, with the remote solar observations yielding approximately 20 times more events than are seen at 1 AU. This divergence persists through 2004. A similar divergence occurs between MCs and non-MC ICMEs. We argue that these divergences are due to the birth of midlatitude active regions, which are the sites of a distinct population of CMEs, only partially intercepted by Earth, and we present a simple geometric argument showing that the CME and ICME rates are consistent with one another. We also acknowledge contributions from (1) an increased rate of high-latitude CMEs and (2) focusing effects from the global solar field. While our analysis, coupled with numerical modeling results, generally supports the interpretation that whether one observes a MC within an ICME is sensitive to the trajectory of the spacecraft through the ICME (i.e., an observational selection effect), one result directly contradicts it. Specifically, we find no systematic offset between the latitudinal origin of ICMEs that contain MCs at 1 AU in the ecliptic plane and that of those that do not.

Riley, Pete↗