Search NASA⌕ Search

SEARCH · Search NASA

Results for “Operator 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 235 records · Page 13

An Investigation of the Kinematic and Microphysical Control of Lightning Rate, Extent and NOX Production using DC3 Observations and the NASA Lightning Nitrogen Oxides Model (LNOM)

The Deep Convective Clouds and Chemistry (DC3) experiment seeks to quantify the relationship between storm physics, lightning characteristics and the production of nitrogen oxides via lightning (LNOx). The focus of this study is to investigate the kinematic and microphysical control of lightning properties, particularly those that may govern LNOx production, such as flash rate, type and extent across Alabama during DC3. Prior studies have demonstrated that lightning flash rate and type is correlated to kinematic and microphysical properties in the mixed‐phase region of thunderstorms such as updraft volume and graupel mass. More study is required to generalize these relationships in a wide variety of storm modes and meteorological conditions. Less is known about the co‐evolving relationship between storm physics, morphology and three‐dimensional flash extent, despite its importance for LNOx production. To address this conceptual gap, the NASA Lightning Nitrogen Oxides Model (LNOM) is applied to North Alabama Lightning Mapping Array (NALMA) and Vaisala National Lightning Detection NetworkTM (NLDN) observations following ordinary convective cells through their lifecycle. LNOM provides estimates of flash rate, flash type, channel length distributions, lightning segment altitude distributions (SADs) and lightning NOx production profiles. For this study, LNOM is applied in a Lagrangian sense to multicell thunderstorms over Northern Alabama on two days during DC3 (21 May and 11 June 2012) in which aircraft observations of NOx are available for comparison. The LNOM lightning characteristics and LNOX production estimates are compared to the evolution of updraft and precipitation properties inferred from dual‐Doppler and polarimetric radar analyses applied to observations from a nearby radar network, including the UAH Advanced Radar for Meteorological and Operational Research (ARMOR). Given complex multicell evolution, particular attention is paid to storm morphology, cell mergers and possible dynamical, microphysical and electrical interaction of individual cells when testing various hypotheses.

Carey, Lawrence↗

Comparing Two Approaches for Assessing Observation Impact

Langland and Baker introduced an approach to assess the impact of observations on the forecasts. In that approach, a state-space aspect of the forecast is defined and a procedure is derived ultimately relating changes in the aspect with changes in the observing system. Some features of the state-space approach are to be noted: the typical choice of forecast aspect is rather subjective and leads to incomplete assessment of the observing system, it requires availability of a verification state that is in practice correlated with the forecast, and it involves the adjoint operator of the entire data assimilation system and is thus constrained by the validity of this operator. This article revisits the topic of observation impacts from the perspective of estimation theory. An observation-space metric is used to allow inferring observation impact on the forecasts without the limitations just mentioned. Using differences of observation-minus-forecast residuals obtained from consecutive forecasts leads to the following advantages: (i) it suggests a rather natural choice of forecast aspect that directly links to the data assimilation procedure, (ii) it avoids introducing undesirable correlations in the forecast aspect since verification is done against the observations, and (iii) it does not involve linearization and use of adjoints. The observation-space approach has the additional advantage of being nearly cost free and very simple to implement. In its simplest form it reduces to evaluating the statistics of observationminus- background and observation-minus-analysis residuals with traditional methods. Illustrations comparing the approaches are given using the NASA Goddard Earth Observing System.

Variational Analysis↗

Radial diffusion in Saturn's radiation belts - A modeling analysis assuming satellite and ring E absorption

A modeling analysis is carried out of six experimental phase space density profiles for nearly equatorially mirroring protons using methods based on the approach of Thomsen et al. (1977). The form of the time-averaged radial diffusion coefficient D(L) that gives an optimal fit to the experimental profiles is determined under the assumption that simple satellite plus Ring E absorption of inwardly diffusing particles and steady-state radial diffusion are the dominant physical processes affecting the proton data in the L range that is modeled. An extension of the single-satellite model employed by Thomsen et al. to a model that includes multisatellite and ring absorption is described, and the procedures adopted for estimating characteristic satellite and ring absorption times are defined. The results obtained in applying three representative solid-body absorption models to evaluate D(L) in the range where L is between 4 and 16 are reported, and a study is made of the sensitivity of the preferred amplitude and L dependence for D(L) to the assumed model parameters. The inferred form of D(L) is then compared with that which would be predicted if various proposed physical mechanisms for driving magnetospheric radial diffusion are operative at Saturn.

Hood, L. L.↗

Portfolios in Stochastic Local Search: Efficiently Computing Most Probable Explanations in Bayesian Networks

Portfolio methods support the combination of different algorithms and heuristics, including stochastic local search (SLS) heuristics, and have been identified as a promising approach to solve computationally hard problems. While successful in experiments, theoretical foundations and analytical results for portfolio-based SLS heuristics are less developed. This article aims to improve the understanding of the role of portfolios of heuristics in SLS. We emphasize the problem of computing most probable explanations (MPEs) in Bayesian networks (BNs). Algorithmically, we discuss a portfolio-based SLS algorithm for MPE computation, Stochastic Greedy Search (SGS). SGS supports the integration of different initialization operators (or initialization heuristics) and different search operators (greedy and noisy heuristics), thereby enabling new analytical and experimental results. Analytically, we introduce a novel Markov chain model tailored to portfolio-based SLS algorithms including SGS, thereby enabling us to analytically form expected hitting time results that explain empirical run time results. For a specific BN, we show the benefit of using a homogenous initialization portfolio. To further illustrate the portfolio approach, we consider novel additive search heuristics for handling determinism in the form of zero entries in conditional probability tables in BNs. Our additive approach adds rather than multiplies probabilities when computing the utility of an explanation. We motivate the additive measure by studying the dramatic impact of zero entries in conditional probability tables on the number of zero-probability explanations, which again complicates the search process. We consider the relationship between MAXSAT and MPE, and show that additive utility (or gain) is a generalization, to the probabilistic setting, of MAXSAT utility (or gain) used in the celebrated GSAT and WalkSAT algorithms and their descendants. Utilizing our Markov chain framework, we show that expected hitting time is a rational function - i.e. a ratio of two polynomials - of the probability of applying an additive search operator. Experimentally, we report on synthetically generated BNs as well as BNs from applications, and compare SGSs performance to that of Hugin, which performs BN inference by compilation to and propagation in clique trees. On synthetic networks, SGS speeds up computation by approximately two orders of magnitude compared to Hugin. In application networks, our approach is highly competitive in Bayesian networks with a high degree of determinism. In addition to showing that stochastic local search can be competitive with clique tree clustering, our empirical results provide an improved understanding of the circumstances under which portfolio-based SLS outperforms clique tree clustering and vice versa.

Mengshoel, Ole J.↗

Pulsed Electric Propulsion Thrust Stand Calibration Method

The evaluation of the performance of any propulsion device requires the accurate measurement of thrust. While chemical rocket thrust is typically measured using a load cell, the low thrust levels associated with electric propulsion (EP) systems necessitate the use of much more sensitive measurement techniques. The design and development of electric propulsion thrust stands that employ a conventional hanging pendulum arm connected to a balance mechanism consisting of a secondary arm and variable linkage have been reported in recent publications by Polzin et al. These works focused on performing steady-state thrust measurements and employed a static analysis of the thrust stand response. In the present work, we present a calibration method and data that will permit pulsed thrust measurements using the Variable Amplitude Hanging Pendulum with Extended Range (VAHPER) thrust stand. Pulsed thrust measurements are challenging in general because the pulsed thrust (impulse bit) occurs over a short timescale (typically 1 micros to 1 millisecond) and cannot be resolved directly. Consequently, the imparted impulse bit must be inferred through observation of the change in thrust stand motion effected by the pulse. Pulsed thrust measurements have typically only consisted of single-shot operation. In the present work, we discuss repetition-rate pulsed thruster operation and describe a method to perform these measurements. The thrust stand response can be modeled as a spring-mass-damper system with a repetitive delta forcing function to represent the impulsive action of the thruster.

Wong, Andrea R.↗

Automatic Facial Expression Recognition and Operator Functional State

The prevalence of human error in safety-critical occupations remains a major challenge to mission success despite increasing automation in control processes. Although various methods have been proposed to prevent incidences of human error, none of these have been developed to employ the detection and regulation of Operator Functional State (OFS), or the optimal condition of the operator while performing a task, in work environments due to drawbacks such as obtrusiveness and impracticality. A video-based system with the ability to infer an individual's emotional state from facial feature patterning mitigates some of the problems associated with other methods of detecting OFS, like obtrusiveness and impracticality in integration with the mission environment. This paper explores the utility of facial expression recognition as a technology for inferring OFS by first expounding on the intricacies of OFS and the scientific background behind emotion and its relationship with an individual's state. Then, descriptions of the feedback loop and the emotion protocols proposed for the facial recognition program are explained. A basic version of the facial expression recognition program uses Haar classifiers and OpenCV libraries to automatically locate key facial landmarks during a live video stream. Various methods of creating facial expression recognition software are reviewed to guide future extensions of the program. The paper concludes with an examination of the steps necessary in the research of emotion and recommendations for the creation of an automatic facial expression recognition program for use in real-time, safety-critical missions.

Blanson, Nina↗

Automatic Facial Expression Recognition and Operator Functional State

The prevalence of human error in safety-critical occupations remains a major challenge to mission success despite increasing automation in control processes. Although various methods have been proposed to prevent incidences of human error, none of these have been developed to employ the detection and regulation of Operator Functional State (OFS), or the optimal condition of the operator while performing a task, in work environments due to drawbacks such as obtrusiveness and impracticality. A video-based system with the ability to infer an individual's emotional state from facial feature patterning mitigates some of the problems associated with other methods of detecting OFS, like obtrusiveness and impracticality in integration with the mission environment. This paper explores the utility of facial expression recognition as a technology for inferring OFS by first expounding on the intricacies of OFS and the scientific background behind emotion and its relationship with an individual's state. Then, descriptions of the feedback loop and the emotion protocols proposed for the facial recognition program are explained. A basic version of the facial expression recognition program uses Haar classifiers and OpenCV libraries to automatically locate key facial landmarks during a live video stream. Various methods of creating facial expression recognition software are reviewed to guide future extensions of the program. The paper concludes with an examination of the steps necessary in the research of emotion and recommendations for the creation of an automatic facial expression recognition program for use in real-time, safety-critical missions

Blanson, Nina↗

The Kinematic and Microphysical Control of Lightning Rate, Extent and NOX Production

The Deep Convective Clouds and Chemistry (DC3) experiment seeks to quantify the relationship between storm physics, lightning characteristics and the production of nitrogen oxides via lightning (LNOx). The focus of this study is to investigate the kinematic and microphysical control of lightning properties, particularly those that may govern LNOx production, such as flash rate, type and extent across Alabama during DC3. Prior studies have demonstrated that lightning flash rate and type is correlated to kinematic and microphysical properties in the mixed-phase region of thunderstorms such as updraft volume and graupel mass. More study is required to generalize these relationships in a wide variety of storm modes and meteorological conditions. Less is known about the co-evolving relationship between storm physics, morphology and three-dimensional flash extent, despite its importance for LNOx production. To address this conceptual gap, the NASA Lightning Nitrogen Oxides Model (LNOM) is applied to North Alabama Lightning Mapping Array (NALMA) and Vaisala National Lightning Detection Network(TM) (NLDN) observations following ordinary convective cells through their lifecycle. LNOM provides estimates of flash rate, flash type, channel length distributions, lightning segment altitude distributions (SADs) and lightning NOx production profiles. For this study, LNOM is applied in a Lagrangian sense to multicell thunderstorms over Northern Alabama on two days during DC3 (21 May and 11 June 2012) in which aircraft observations of NOx are available for comparison. The LNOM lightning characteristics and LNOX production estimates are compared to the evolution of updraft and precipitation properties inferred from dual-Doppler and polarimetric radar analyses applied to observations from a nearby radar network, including the UAH Advanced Radar for Meteorological and Operational Research (ARMOR). Given complex multicell evolution, particular attention is paid to storm morphology, cell mergers and possible dynamical, microphysical and electrical interaction of individual cells when testing various hypotheses.

Carey, Lawrence↗

Quantum optimal control of superconducting qubits based on machine-learning characterization

Implementing fast and high-fidelity quantum operations using open-loop quantum optimal control relies on having an accurate model of the quantum dynamics. Any deviations between this model and the complete dynamics of the device, such as the presence of spurious modes or pulse distortions, can degrade the performance of optimal controls in practice. Here, we propose an experimentally simple approach to realize optimal quantum controls tailored to the device parameters and environment while specifically characterizing this quantum system. Concretely, we use physics-inspired machine learning to infer an accurate model of the dynamics from experimentally available data and then optimize our experimental controls on this trained model. We show the power and feasibility of this approach by optimizing arbitrary single-qubit operations in detailed numerical simulations of a superconducting transmon qubit. Furthermore, we demonstrate that this framework produces an accurate description of the device dynamics under arbitrary controls, together with the precise pulses achieving arbitrary single-qubit gates with a high fidelity of ∼99.99%.

Artificial neural networks↗

Developing Dual Polarization Applications For 45th Weather Squadron's (45 WS) New Weather Radar: A Cooperative Project With The National Space Science and Technology Center (NSSTC)

A new weather radar is being acquired for use in support of America s space program at Cape Canaveral Air Force Station, NASA Kennedy Space Center, and Patrick AFB on the east coast of central Florida. This new radar includes dual polarization capability, which has not been available to 45 WS previously. The 45 WS has teamed with NSSTC with funding from NASA Marshall Spaceflight Flight Center to improve their use of this new dual polarization capability when it is implemented operationally. The project goals include developing a temperature profile adaptive scan strategy, developing training materials, and developing forecast techniques and tools using dual polarization products. The temperature profile adaptive scan strategy will provide the scan angles that provide the optimal compromise between volume scan rate, vertical resolution, phenomena detection, data quality, and reduced cone-of-silence for the 45 WS mission. The mission requirements include outstanding detection of low level boundaries for thunderstorm prediction, excellent vertical resolution in the atmosphere electrification layer between 0 C and -20 C for lightning forecasting and Lightning Launch Commit Criteria evaluation, good detection of anvil clouds for Lightning Launch Commit Criteria evaluation, reduced cone-of-silence, fast volume scans, and many samples per pulse for good data quality. The training materials will emphasize the appropriate applications most important to the 45 WS mission. These include forecasting the onset and cessation of lightning, forecasting convective winds, and hopefully the inference of electrical fields in clouds. The training materials will focus on annotated radar imagery based on products available to the 45 WS. Other examples will include time sequenced radar products without annotation to simulate radar operations. This will reinforce the forecast concepts and also allow testing of the forecasters. The new dual polarization techniques and tools will focus on the appropriate applications for the 45 WS mission. These include forecasting the onset of lightning, the cessation of lightning, convective winds, and hopefully the inference of electrical fields in clouds. This presentation will report on the results achieved so far in the project.

Roeder, W.P.↗

Near Real-Time Probabilistic Damage Diagnosis Using Surrogate Modeling and High Performance Computing

This work investigates novel approaches to probabilistic damage diagnosis that utilize surrogate modeling and high performance computing (HPC) to achieve substantial computational speedup. Motivated by Digital Twin, a structural health management (SHM) paradigm that integrates vehicle-specific characteristics with continual in-situ damage diagnosis and prognosis, the methods studied herein yield near real-time damage assessments that could enable monitoring of a vehicle's health while it is operating (i.e. online SHM). High-fidelity modeling and uncertainty quantification (UQ), both critical to Digital Twin, are incorporated using finite element method simulations and Bayesian inference, respectively. The crux of the proposed Bayesian diagnosis methods, however, is the reformulation of the numerical sampling algorithms (e.g. Markov chain Monte Carlo) used to generate the resulting probabilistic damage estimates. To this end, three distinct methods are demonstrated for rapid sampling that utilize surrogate modeling and exploit various degrees of parallelism for leveraging HPC. The accuracy and computational efficiency of the methods are compared on the problem of strain-based crack identification in thin plates. While each approach has inherent problem-specific strengths and weaknesses, all approaches are shown to provide accurate probabilistic damage diagnoses and several orders of magnitude computational speedup relative to a baseline Bayesian diagnosis implementation.

Warner, James E.↗

Graph neural network for neutrino physics event reconstruction

Liquid argon time projection chamber (LArTPC) detector technology offers a wealth of high-resolution information on particle interactions, and leveraging that information to its full potential requires sophisticated automated reconstruction techniques. Here, this article describes NUGRAPH 2, a graph neural network for low-level reconstruction of simulated neutrino interactions in a LArTPC detector. Simulated neutrino interactions in the MicroBooNE detector geometry are described as heterogeneous graphs, with energy depositions on each detector plane forming nodes on planar subgraphs. The network utilizes a multihead attention message-passing mechanism to perform background filtering and semantic labeling on these graph nodes, identifying those associated with the primary physics interaction with 98.0% efficiency and labeling them according to particle type with 94.9% efficiency. The network operates directly on detector observables across multiple two-dimensional representations but utilizes a three-dimensional-context-aware mechanism to encourage consistency between these representations. Model inference takes 0.12 s / event on a CPU and 0.005 s / event batched on a GPU. This architecture is designed to be a general-purpose solution for particle reconstruction in neutrino physics, with the potential for deployment across a broad range of detector technologies, and offers a core convolution engine that can be leveraged for a variety of tasks beyond the two described in this paper.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Basic data requirements for microwave radiometer systems

Microwave radiometry has emerged over the last two decades to become an integral part of the field of environmental remote sensing. Numerous investigations were conducted to evaluate the use of microwave radiometry for atmospheric, oceanographic, hydrological, and geological applications. Remote sensing of the earth using microwave radiometry began in 1968 by the Soviet satellite Cosmos 243, which included four microwave radiometers (Ulably, 1981). Since then, microwave radiometers were included onboard many spacecraft, and were used to infer many physical parameters. Some of the basic concepts of radiometric emission and measurement will be discussed. Several radiometer systems are presented and an overview of their operation is discussed. From the description of the radiometer operation the data stream required from the radiometer and the general type of algorithm required for the measurement is discussed.

Lawrence, R. W.↗

Measurement and interpretation of crustal deformation rates associated with postglacial rebound

This project involves obtaining GPS measurements in Scandinavia and using the measurements to estimate the viscosity profile of the earth's mantle and to correct tide-gauge measurements for the rebound effect. We report on several aspects of this project. The DSGS was not scheduled to be reoccupied with DOSE receivers during the report period. The permanent network set up by Onsala Space Observatory continues to operate, and the data are being evaluated. An important technical advance we intend for this project is to use the full three dimensional site velocity information for inferring geophysical parameters. During the report period, two papers have been been accepted for publication in the Journal of Geophysical Research and will be published in April. Reprints of these papers are contained in the Appendix.

Davis, James L.↗

Implementation and Comparison of Acoustic Travel-Time Measurement Procedures for the Helioseismic and Magnetic Imager Time-Distance Helioseismology Pipeline

The Helioseismic and Magnetic Imager (HMI) instrument on board the Solar Dynamics Observatory (SDO) satellite is designed to produce high-resolution Doppler velocity maps of oscillations at the solar surface with high temporal cadence. To take advantage of these high-quality oscillation data, a time-distance helioseismology pipeline has been implemented at the Joint Science Operations Center (JSOC) at Stanford University. The aim of this pipeline is to generate maps of acoustic travel times from oscillations on the solar surface, and to infer subsurface 3D flow velocities and sound-speed perturbations. The wave travel times are measured from cross covariances of the observed solar oscillation signals. For implementation into the pipeline we have investigated three different travel-time definitions developed in time-distance helioseismology: a Gabor wavelet fitting (Kosovichev and Duvall, 1997), a minimization relative to a reference cross-covariance function (Gizon and Birch, 2002), and a linearized version of the minimization method (Gizon and Birch, 2004). Using Doppler velocity data from the Michelson Doppler Imager (MDI) instrument on board SOHO, we tested and compared these definitions for the mean and difference travel-time perturbations measured from reciprocal signals. Although all three procedures return similar travel times in a quiet Sun region, the method of Gizon and Birch (2004) gives travel times that are significantly different from the others in a magnetic (active) region. Thus, for the pipeline implementation we chose the procedures of Kosovichev and Duvall (1997) and Gizon and Birch (2002). We investigated the relationships among these three travel-time definitions, their sensitivities to fitting parameters, and estimated the random errors they produce

Couvidat, S.↗

The Kinematic and Microphysical Control of Storm Integrated Lightning Flash Extent

The objective of this preliminary study is to investigate the kinematic and microphysical control of lightning properties, particularly those that may govern the production of nitrogen oxides (NOx) in thunderstorms, such as flash rate, type and extent. The mixed-phase region is where the noninductive charging (NIC) process is thought to generate most storm electrification during rebounding collisions between ice particles in the presence of supercooled water. As a result, prior radar-based studies have demonstrated that lightning flash rate is well correlated to kinematic and microphysical properties in the mixed-phase region of thunderstorms such as updraft volume, graupel mass, or ice mass flux. There is also some evidence that lightning type is associated with the convective state. Intracloud (IC) lightning tends to dominate during the updraft accumulation of precipitation ice mass while cloud-to-ground (CG) lightning is more numerous during the downdraft-driven descent of radar echo associated with graupel and hail. More study is required to generalize these relationships, especially regarding lightning type, in a wide variety of storm modes and meteorological conditions. Less is known about the co-evolving relationship between storm kinematics, microphysics, morphology and three-dimensional flash extent, despite its importance for lightning NOx production. To address this conceptual gap, the NASA MSFC Lightning Nitrogen Oxides Model (LNOM) is applied to North Alabama Lightning Mapping Array (NALMA) and Vaisala National Lightning Detection NetworkTM (NLDN) observations following ordinary convective cells through their lifecycle. LNOM provides estimates of flash type, channel length distributions, lightning segment altitude distributions (SADs) and lightning NOx production profiles. For this study, LNOM is applied in a Lagrangian sense to well isolated convective cells on 3 April 2007 (single cell and multi-cell hailstorm, non-severe multicell) and 6 July 2007 (non-severe multi-cell) over Northern Alabama. The LNOM lightning characteristics are compared to the evolution of updraft and precipitation properties inferred from dual-Doppler and polarimetric radar analyses applied to observations from a nearby Doppler radar network, including the UA Huntsville Advanced Radar for Meteorological and Operational Research (ARMOR, C-band, polarimetric). The LNOM estimated SAD and lightning NOx production profiles are placed in the context of radar derived profiles of vertical motion, precipitation types and amounts. Finally, these analyses are used to determine if storm integrated flash channel extent is as well correlated to volumetric updraft and precipitation ice characteristics in the mixed phase region as flash rate for these individual convective cells.

Carey, Lawrence↗

Implementation and Comparison of Acoustic Travel-Time Measurement Procedures for the Solar Dynamics Observatory-Helioseismic and Magnetic Imager Time-Distance Helioseismology Pipeline

The Helioseismic and Magnetic Imager (HMI) instrument onboard the Solar Dynamics Observatory (SDO) satellite is designed to produce high-resolution Doppler-velocity maps of oscillations at the solar surface with high temporal cadence. To take advantage of these high-quality oscillation data, a time - distance helioseismology pipeline (Zhao et al., Solar Phys. submitted, 2010) has been implemented at the Joint Science Operations Center (JSOC) at Stanford University. The aim of this pipeline is to generate maps of acoustic travel times from oscillations on the solar surface, and to infer subsurface 3D flow velocities and sound-speed perturbations. The wave travel times are measured from cross-covariances of the observed solar oscillation signals. For implementation into the pipeline we have investigated three different travel-time definitions developed in time - distance helioseismology: a Gabor-wavelet fitting (Kosovichev and Duvall, SCORE'96: Solar Convection and Oscillations and Their Relationship, ASSL, Dordrecht, 241, 1997), a minimization relative to a reference cross-covariance function (Gizon and Birch, Astrophys. J. 571, 966, 2002), and a linearized version of the minimization method (Gizon and Birch, Astrophys. J. 614, 472, 2004). Using Doppler-velocity data from the Michelson Doppler Imager (MDI) instrument onboard SOHO, we tested and compared these definitions for the mean and difference traveltime perturbations measured from reciprocal signals. Although all three procedures return similar travel times in a quiet-Sun region, the method of Gizon and Birch (Astrophys. J. 614, 472, 2004) gives travel times that are significantly different from the others in a magnetic (active) region. Thus, for the pipeline implementation we chose the procedures of Kosovichev and Duvall (SCORE'96: Solar Convection and Oscillations and Their Relationship, ASSL, Dordrecht, 241, 1997) and Gizon and Birch (Astrophys. J. 571, 966, 2002). We investigated the relationships among these three travel-time definitions, their sensitivities to fitting parameters, and estimated the random errors that they produce.

HMI↗