Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “data 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 937 records · Page 52

Analysis of hydrogen Lyman-alpha observations of the coma of Comet P/Halley near the perihelion

The pioneer Venus Orbiter Ultraviolet Spectrometer measurements of the Lyman-alpha intensity of atomic hydrogen excited by solar resonance scattering in the coma of Comet P/Halley acquired from December 28, 1985, to January 6, 1986, and from January 31, 1986, to March 6, 1986, are simulated with the Monte Carlo Particle Trajectory Model corrected for optical depth effects. Spatially detailed comparisons between data and model show excellent agreement and are used to infer that the highest cometary activity may not be at perihelion, but about 2 1/2 weeks before. An improved set of H2O production rates is presented for the period of time that the spectrometer was observing and is found to be consistent with the rates from other types of measurements. The apparent discrepancy between Stewart (1987) in early March and International Ultraviolet Explorer OH derived rates is resolved. The problem with the conversion of 18-cm OH radio brightness to H2O production rates is also discussed.

Smyth, William H.↗

Optical fiber sensors for damage analysis in aerospace materials

Under this grant, fiber optic sensors were investigated for use in the nondestructive evaluation of aging aircraft. Specifically, optical fiber sensors for detection and location of impacts on a surface, and for detection of corrosion in metals were developed. The use of neural networks was investigated for determining impact location by processing the output of a network of fiberoptic strain sensors distributed on a surface. This approach employs triangulation to determine location by comparing the arrival times at several sensors, of the acoustic signal generated by the impact. For this study, a neural network simulator running on a personal computer was used to train a network using a back-propagation algorithm. Fiber optic extrinsic Fabry-Perot interferometer (EFPI) strain sensors are attached to or embedded in the surface, so that stress waves emanating from an impact can be detected. The ability of the network to determine impact location by time-or-arrival of acoustic signals was assessed by comparing network outputs with actual experimental results using impacts on a panel instrumented with optical fiber sensors. Using the neural network to process the sensor outputs, the impact location can be inferred to centimeter range accuracy directly from the arrival time data. In addition, the network can be trained to determine impact location, regardless of material anisotropy. Results demonstrate that a back-propagation network identifies impact location for an anisotropic graphite/bismaleimide plate with the same accuracy as that for an isotropic aluminum plate. Two different approaches were investigated for the development of fiber optic sensors for corrosion detection in metals, both utilizing optical fiber sensors with metal coatings. In the first approach, an extrinsic Fabry-Perot interferometric fiber optic strain sensor was placed under tensile stress, and while in the resulting strained position, a thick coating of metal was applied. Due to an increase in the quantity of material, the sensor does not return to its original position upon removal of the applied stress, and some residual strain is maintained within the sensor element. As the metal thickness decreases due to corrosion, this strain is released, providing the sensing mechanism for corrosion detection. In the second approach, photosensitive optical fibers with long period Bragg gratings in the core were coated with metal. The Bragg gratings serve to couple core modes at discrete wavelengths to cladding modes. Since cladding modes interact with the metal coating surrounding the fiber cladding, the specific wavelengths coupled from core to cladding depend on the refractive index of the metal coating. Therefore, as the metal corrodes, the resulting change in index of the coating may be measured by measuring the change in wavelength of the coupled mode. Results demonstrate that both approaches can be successfully used to track the loss in metal coating on the optical fiber sensors due to corrosion.

Schindler, Paul↗

Climatology, Natural Cycles, and Modes of Interannual Variability of the Great Plains Low-Level Jet as Assimilated by the GEOS-1 Data Analysis System

Despite the fact that the low-level jet of the southern Great Plains (the GPLLJ) of the U.S. is primarily a nocturnal phenomenon that virtually vanishes during the daylight hours, it is one of the most persistent and stable features of the low-level continental flow during the warm-season months, May through August. We have first used significant-level data to validate the skill of the GEOS-1 Data Assimilation System (DAS) in realistically detecting this jet and inferring its structure and evolution. We have then carried out a 15-year reanalysis with the GEOS-1 DAS to determine and validate its climatology and mean diurnal cycle and to study its interannual variability. Interannual variability of the GPLLJ is much smaller than mean diurnal and random intraseasonal variability and comparable in magnitude, but not location, to mean seasonal variability. There are three maxima of interannual low-level meridional flow variability of the GPLLJ over the upper Great Plains, southeastern Texas, and the western Gulf of Mexico. Cross-sectional profiles of mean southerly wind through the Texas maximum remain relatively stable and recognizable from year to year with only its eastward flank showing significant variability. This variability, however, exhibits a distinct, biennial oscillation during the first six years of the reanalysis period and only then. Each of the three variability maxima corresponds to a spatially coherent, jet-like pattern of low-level flow interannual variability. There are three prominent modes of interannual. variability. These include the intermittent biennial oscillation (IBO), local to the Texas maximum. Its signal is evident in surface pressure, surface temperature, ground wetness and upper air flow, as well. A larger-scale continental convergence pattern (CCP) of covariance, exhibiting strong anti-correlation between the flow near the Texas and the upper Great Plains variability maxima, is revealed only when the IBO is removed from the interannual time series. A third, subtropical mode of covariance is associated with the Gulf of Mexico variability maximum. Significant interannual anti-correlations of the southeasterly flow over the Arizona/New Mexico region with the CCP and the subtropical mode are enhanced when restricted to the month of July. These anti-correlations may relate to an observed out-of-phase precipitation relationship between the Great Plains and the southwestern U.S.. The typical duration of interannual low-level meridional wind anomalies within a given season increases over the continent with decreasing latitude from two to three weeks over the upper Great Plains to six to seven weeks over eastern Texas.

Helfand, H. M.↗

Mars Exploration Rovers as Virtual Instruments for Determination of Terrain Roughness and Physical Properties

The two 2003 Mars Exploration Rovers (MERs), in combination with the Athena Payload, will be used as virtual instrument systems to infer terrain properties during traverses, in addition to using the rover wheels to excavate trenches, exposing subsurface materials for remote and in-situ observations. The MERs are being modeled using finite element-based rover system transfer functions that utilize the distribution of masses associated with the vehicle, together with suspension and wheel dynamics, to infer surface roughness and mechanical properties from traverse time series data containing vehicle yaw, pitch, roll, encoder counts, and motor currents. These analyses will be supplemented with imaging and other Athena Payload measurements. The approach is being validated using Sojourner data, the FIDO rover, and experiments with MER testbed vehicles. In addition to conducting traverse science and associated analyses, trenches will be excavated by the MERs to depths of approximately 10-20 cm by locking all but one of the front wheels and rotating that wheel backwards so that the excavated material is piled up on the side of the trench away from the vehicle. Soil cohesion and angle of internal friction will be determined from the trench telemetry data. Emission spectroscopy and in-situ observations will be made using the Athena payload before and after imaging. Trenching and observational protocols have been developed using Sojourner results; data from the FIDO rover, including trenches dug into sand, mud cracks, and weakly indurated bedrock; and experiments with MER testbed rovers. Particular attention will be focused on Mini-TES measurements designed to determine the abundance and state of subsurface water (e.g. hydrated, in zeolites, residual pore ice?) predicted to be present from Odyssey GRS/NS/HEND data.

R E Arvidson↗

On Iron Enrichment, Star Formation, and Type Ia Supernovae in Galaxy Clusters

The nature of star formation and Type Ia supernovae (SNIa) in galaxies in the field and in rich galaxy clusters are contrasted by juxtaposing the buildup of heavy metals in the universe inferred from observed star formation and supernovae rate histories with data on the evolution of Fe abundances in the intracluster medium (ICM). Models for the chemical evolution of Fe in these environments are constructed, subject to observational constraints, for this purpose. While models with a mean delay for SNIa of 3 Gyr and standard initial mass function (IMF) are fully consistent with observations in the field, cluster Fe enrichment immediately tracked a rapid, top-heavy phase of star formation - although transport of Fe into the ICM may have been more prolonged and star formation likely continued beyond redshift 1. The means of this prompt enrichment consisted of SNII yielding greater than or equal to 0.1 solar mass per explosion (if the SNIa rate normalization is scaled down from its value in the field according to the relative number of candidate progenitor stars in the 3 - 8 solar mass range) and/or SNIa with short delay times originating during the rapid star formation epoch. Star formation is greater than 3 times more efficient in rich clusters than in the field, mitigating the overcooling problem in numerical cluster simulations. Both the fraction of baryons cycled through stars, and the fraction of the total present-day stellar mass in the form of stellar remnants, are substantially greater in clusters than in the field.

Loewenstein, Michael↗

Theory of Radiation Transfer in Neutron Star Atmospheres

The possibility for direct investigation of thermal emission from isolated neutron stars opened about a quarter of century ago with the launch of the first X-ray observatories Einstein and EXOSAT stimulated developing models of the neutron star surface radiation which began at the end of 80's. Confronting observational data with theoretical models of thermal emission allows one to infer the surface temperatures, magnetic fields, chemical composition, and neutron star masses and radii. This information, supplemented with the model equations of state and neutron star cooling models, provides an opportunity to understand the fundamental properties of the superdense matter in the stars' interiors. Almost all available models are based on the assumption that thermal radiation emitted by a neutron star is formed in the superficial star's layers--atmosphere. The neutron star atmospheres are very different from those of usual stars due to the immense gravity and huge magnetic fields. In this presentation we review the current status of the neutron star atmosphere modeling, present most important results, discuss problems and possible future developments.

Zavlin, Vyacheslav↗

Analysis of Terrestrial Water Storage Changes from GRACE and GLDAS

Since March 2002, the Gravity Recovery and Climate Experiment (GRACE) has provided first estimates of land water storage variations by monitoring the time-variable component of Earth's gravity field. Here we characterize spatial-temporal variations in terrestrial water storage changes (TWSC) from GRACE and compare them to those simulated with the Global Land Data Assimilation System (GLDAS). Additionally, we use GLDAS simulations to infer how TWSC is partitioned into snow, canopy water and soil water components, and to understand how variations in the hydrologic fluxes act to enhance or dissipate the stores. Results quantify the range of GRACE-derived storage changes during the studied period and place them in the context of seasonal variations in global climate and hydrologic extremes including drought and flood, by impacting land memory processes. The role of the largest continental river basins as major locations for freshwater redistribution is highlighted. GRACE-based storage changes are in good agreement with those obtained from GLDAS simulations. Analysis of GLDAS-simulated TWSC illustrates several key characteristics of spatial and temporal land water storage variations. Global averages of TWSC were partitioned nearly equally between soil moisture and snow water equivalent, while zonal averages of TWSC revealed the importance of soil moisture storage at low latitudes and snow storage at high latitudes. Evapotranspiration plays a key role in dissipating globally averaged terrestrial water storage. Latitudinal averages showed how precipitation dominates TWSC variations in the tropics, evapotranspiration is most effective in the midlatitudes, and snowmelt runoff is a key dissipating flux at high latitudes. Results have implications for monitoring water storage response to climate variability and change, and for constraining land model hydrology simulations.

Syed, Tajdarul H.↗

An Analysis of AERONET Aerosol Absorption Properties and Classifications Representative of Aerosol Source Regions

Partitioning of mineral dust, pollution, smoke, and mixtures using remote sensing techniques can help improve accuracy of satellite retrievals and assessments of the aerosol radiative impact on climate. Spectral aerosol optical depth (tau) and single scattering albedo (omega (sub 0) ) from Aerosol Robotic Network (AERONET) measurements are used to form absorption [i.e., omega (sub 0) and absorption Angstrom exponent (alpha(sub abs))] and size [i.e., extinction Angstrom exponent (alpha(sub ext)) and fine mode fraction of tau] relationships to infer dominant aerosol types. Using the long-term AERONET data set (1999-2010), 19 sites are grouped by aerosol type based on known source regions to: (1) determine the average omega (sub 0) and alpha(sub abs) at each site (expanding upon previous work); (2) perform a sensitivity study on alpha(sub abs) by varying the spectral omega (sub 0); and (3) test the ability of each absorption and size relationship to distinguish aerosol types. The spectral omega (sub 0) averages indicate slightly more aerosol absorption (i.e., a 0.0 < delta omega (sub 0) <= 0.02 decrease) than in previous work and optical mixtures of pollution and smoke with dust show stronger absorption than dust alone. Frequency distributions of alpha(sub abs) show significant overlap among aerosol type categories and at least 10% of the alpha(sub abs) retrievals in each category are below 1.0. Perturbing the spectral omega (sub 0) by +/- 0.03 induces significant alpha(sub abs) changes from the unperturbed value by at least approx. +/- 0.6 for Dust, approx. +/-0.2 for Mixed, and approx. +/-0.1 for Urban/Industrial and Biomass Burning. The omega (sub 0)440nm and alpha(sub ext) 440-870nm relationship shows the best separation among aerosol type clusters, providing a simple technique for determining aerosol type from surface- and future space-based instrumentation.

Giles, David M.↗

Geophysical Information from Advanced Sounder Infrared Spectral Radiance

Advanced satellite sensors are tasked with improving global observations of the Earth's atmosphere, clouds, and surface to enable enhancements in weather prediction, climate monitoring capability, and environmental change detection. Satisfying this type of improvement for inferred geophysical information from these observations requires optimal usage of data from current systems as well as enhancements to future sensors. This presentation addresses the information content present in infrared spectral radiance from advanced atmospheric sounders with an emphasis on knowledge of thermodynamic state and trace species. Results of trade-off studies conducted to evaluate the impact of spectral resolution, spectral coverage, instrument noise, and a priori knowledge on remote sensing system information content will be discussed. A focus is placed on information achievable from the Atmospheric InfraRed Sounder (AIRS) on the NASA EOS Aqua satellite in orbit since 2002, the Infrared Atmospheric Sounding Interferometer (IASI) aboard MetOp-A since 2006, and the Cross-track Infrared Sounder (CrIS) instrument aboard the NPP and JPSS series of satellites which began 28 October 2011.

Larar, Allen M.↗

A Corrected Formulation of the Multilayer Model (MLM) for Inferring Gaseous Dry Deposition to Vegetated Surfaces

The Multilayer Model (MLM) has been used for many years to infer dry deposition fluxes from measured trace species concentrations and standard meteorological measurements for national networks in the U.S., including the U.S. Environmental Protection Agency's Clean Air Status and Trends Network (CASTNet). MLM utilizes a resistance analogy to calculate deposition velocities appropriate for whole vegetative canopies, while employing a multilayer integration to account for vertically varying meteorology, canopy morphology and radiative transfer within the canopy. However, the MLM formulation, as it was originally presented and as it has been subsequently employed, contains a non-physical representation related to the leaf-level quasi-laminar boundary layer resistance that affects the calculation of the total canopy resistance. In this note, the non-physical representation of the canopy resistance as originally formulated in MLM is discussed and a revised, physically consistent, formulation is suggested as a replacement. The revised canopy resistance formulation reduces estimates of HNO3 deposition velocities by as much as 38% during mid-day as compared to values generated by the original formulation. Inferred deposition velocities for SO2 and O3 are not significantly altered by the change in formulation (less than 3%). Inferred deposition loadings of oxidized and total nitrogen from CASTNet data may be reduced by 10-20% and 5-10%, respectively, for the Eastern U. S. when employing the revised formulation of MLM as compared to the original formulation.

vegetated surfaces↗

AeroCom Phase III Multi-Model Evaluation of the Aerosol Life Cycle and Optical Properties Using Ground and Space-Based Remote Sensing as Well as Surface In Situ Observations

Within the framework of the AeroCom (Aerosol Comparisons between Observations and Models) initiative, the state-of-the-art modelling of aerosol optical properties is assessed from 14 global models participating in the phase III control experiment (AP3). The models are similar to CMIP6/AerChemMIP Earth System Models (ESMs) and provide a robust multi-model ensemble. Inter-model spread of aerosol species lifetimes and emissions appears to be similar to that of mass extinction coefficients (MECs), suggesting that aerosol optical depth (AOD) uncertainties are associated with a broad spectrum of parameterised aerosol processes. Total AOD is approximately the same as in AeroCom phase I (AP1) simulations. However, we find a 50 % decrease in the optical depth (OD) of black carbon (BC), attributable to a combination of decreased emissions and lifetimes. Relative contributions from sea salt (SS) and dust (DU) have shifted from being approximately equal in AP1 to SS contributing about 2∕3 of the natural AOD in AP3. This shift is linked with a decrease in DU mass burden, a lower DU MEC, and a slight decrease in DU lifetime, suggesting coarser DU particle sizes in AP3 compared to AP1. Relative to observations, the AP3 ensemble median and most of the participating models underestimate all aerosol optical properties investigated, that is, total AOD as well as fine and coarse AOD (AODf, AODc), Ångström exponent (AE), dry surface scattering (SCdry), and absorption (ACdry) coefficients. Compared to AERONET, the models underestimate total AOD by ca. 21 % ± 20 % (as inferred from the ensemble median and interquartile range). Against satellite data, the ensemble AOD biases range from −37 % (MODIS-Terra) to −16 % (MERGED-FMI, a multi-satellite AOD product), which we explain by differences between individual satellites and AERONET measurements themselves. Correlation coefficients (R) between model and observation AOD records are generally high (R>0.75), suggesting that the models are capable of capturing spatio-temporal variations in AOD. We find a much larger underestimate in coarse AODc (∼ −45 % ± 25 %) than in fine AODf (∼ −15 % ± 25 %) with slightly increased inter-model spread compared to total AOD. These results indicate problems in the modelling of DU and SS. The AODc bias is likely due to missing DU over continental land masses (particularly over the United States, SE Asia, and S. America), while marine AERONET sites and the AATSR SU satellite data suggest more moderate oceanic biases in AODc. Column AEs are underestimated by about 10 % ± 16 %. For situations in which measurements show AE > 2, models underestimate AERONET AE by ca. 35 %. In contrast, all models (but one) exhibit large overestimates in AE when coarse aerosol dominates (bias ca. +140 % if observed AE < 0.5). Simulated AE does not span the observed AE variability. These results indicate that models overestimate particle size (or underestimate the fine-mode fraction) for fine-dominated aerosol and underestimate size (or overestimate the fine-mode fraction) for coarse-dominated aerosol. This must have implications for lifetime, water uptake, scattering enhancement, and the aerosol radiative effect, which we can not quantify at this moment. Comparison against Global Atmosphere Watch (GAW) in situ data results in mean bias and inter-model variations of −35 % ± 25 % and −20 % ± 18 % for SCdry and ACdry, respectively. The larger underestimate of SCdry than ACdry suggests the models will simulate an aerosol single scattering albedo that is too low. The larger underestimate of SCdry than ambient air AOD is consistent with recent findings that models overestimate scattering enhancement due to hygroscopic growth. The broadly consistent negative bias in AOD and surface scattering suggests an underestimate of aerosol radiative effects in current global aerosol models. Considerable inter-model diversity in the simulated optical properties is often found in regions that are, unfortunately, not or only sparsely covered by ground-based observations. This includes, for instance, the Sahara, Amazonia, central Australia, and the South Pacific. This highlights the need for a better site coverage in the observations, which would enable us to better assess the models, but also the performance of satellite products in these regions. Using fine-mode AOD as a proxy for present-day aerosol forcing estimates, our results suggest that models underestimate aerosol forcing by ca. −15 %, however, with a considerably large interquartile range, suggesting a spread between −35 % and +10 %.

space-based remote sensing↗

The Atmospheric Response to High Nonthermal Electron Beam Fluxes in Solar Flares. I. Modeling the Brightest NUV Footpoints in the X1 Solar Flare of 2014 March 29

The 2014 March 29 X1 solar flare (SOL20140329T17:48) produced bright continuum emission in the far- and near-ultraviolet (NUV) and highly asymmetric chromospheric emission lines, providing long-sought constraints on the heating mechanisms of the lower atmosphere in solar flares. We analyze the continuum and emission line data from the Interface Region Imaging Spectrograph (IRIS) of the brightest flaring magnetic footpoints in this flare. We compare the NUV spectra of the brightest pixels to new radiative-hydrodynamic predictions calculated with the RADYN code using constraints on a nonthermal electron beam inferred from the collisional thick-target modeling of hard X-ray data from Reuven Ramaty High Energy Solar Spectroscopic Imager. We show that the atmospheric response to a high beam flux density satisfactorily achieves the observed continuum brightness in the NUV. The NUV continuum emission in this flare is consistent with hydrogen (Balmer) recombination radiation that originates from low optical depth in a dense chromospheric condensation and from the stationary beam-heated layers just below the condensation. A model producing two flaring regions (a condensation and stationary layers) in the lower atmosphere is also consistent with the asymmetric Fe II chromospheric emission line profiles observed in the impulsive phase.

Kowalski, Adam F.↗

A test of alternative Caribbean plate relative motion models

The new NUVEL-1 data set for global relative plate motions is used here to discriminate between the two prevailing models for Caribbean plate motion. One model, by Jordan (1975), assumes that North America-Caribbean motion is reflected by the spreading rate inferred from magnetic anomalies at the Cayman Spreading Center and the azimuths of nearby transforms. The other model, by Sykes et al. (1982), uses rates and azimuths inferred from the geometry of the Lesser Antilles Wadati-Benioff zone. Overall, it is found that the data fit the Jordan geometry better, that the data used in global plate motion models are more suitable than rates and azimuths inferred from the geometry of the Wadati-Benioff zone for determining relative motions, and that incorporation of all relevant plate boundaries is essential.

Stein, Seth↗

NuSTAR Detection of X-Rays Concurrent with Gamma-Rays in the Nova V5855 Sgr

We report the first detection of hard (>10 keV) X-ray emission simultaneous with gamma-rays in a nova eruption. Observations of the nova V5855 Sgr carried out with the NuSTAR satellite on Day 12 of the eruption revealed faint, highly absorbed thermal X-rays. The extreme equivalent hydrogen column density toward the X-ray emitting region (∼3 × 10(exp 24 per sq.cm) indicates that the shock producing the X-rays was deeply embedded within the nova ejecta. The slope of the X-ray spectrum favors a thermal origin for the bulk of the emission, and the constraints of the temperature in the shocked region suggest a shock velocity compatible with the ejecta velocities inferred from optical spectroscopy. While we do not claim the detection of nonthermal X-rays, the data do not allow us to rule out an additional, fainter component dominating at energies above 20 keV, for which we obtained upper limits. The inferred luminosity of the thermal X-rays is too low to be consistent with the gamma-ray luminosities if both are powered by the same shock under standard assumptions regarding the efficiency of nonthermal particle acceleration and the temperature distribution of the shocked gas.

Nelson, Thomas↗

The state of the atmosphere as inferred from the FGGE satellite observing systems during SOP-1

The statistical properties, and coverage, of satellite temperature sounding data are described. Tropical regions are observed every two days, extratropics from one to four times a day. Oceans are covered two to three times a day. Asynoptic coverage is comparable to the U.S. rawinsonde network twice daily coverage. Lack of ground truth for data sparse areas makes accuracy difficult to assess. The rms differences of layer mean temperatures obtained from collocating rawinsonde observations with satellite temperature profiles in space and time differ from rms differences of layer mean satellite temperature soundings. The FGGE satellite systems can infer the three dimensional motion field and improve the representation of the large scale state of the atmosphere.

Halem, M.↗

Empirical relationships among atmospheric variables from rawinsonde and field data as surrogates for AVIRIS measurements: Estimation of regional land surface evapotranspiration

Empirical relationships between variables are ways of securing estimates of quantities difficult to measure by remote sensing methods. The use of empirical functions was explored between: (1) atmospheric column moisture abundance W (gm H2O/cm(sup 2) and surface absolute water vapor density rho(q-bar) (gm H2O/cm(sup 3), with rho density of moist air (gm/cm(sup 3), q-bar specific humidity (gm H2O/gm moist air), and (2) column abundance and surface moisture flux E (gm H2O/(cm(sup 2)sec)) to infer regional evapotranspiration from Airborne Visible/Infrared Imaging Spectrometers (AVIRIS) water vapor mapping data. AVIRIS provides, via analysis of atmospheric water absorption features, estimates of column moisture abundance at very high mapping rate (at approximately 100 km(sup 2)/40 sec) over large areas at 20 m ground resolution.

Conel, James E.↗

A statistical model for interpreting computerized dynamic posturography data

Computerized dynamic posturography (CDP) is widely used for assessment of altered balance control. CDP trials are quantified using the equilibrium score (ES), which ranges from zero to 100, as a decreasing function of peak sway angle. The problem of how best to model and analyze ESs from a controlled study is considered. The ES often exhibits a skewed distribution in repeated trials, which can lead to incorrect inference when applying standard regression or analysis of variance models. Furthermore, CDP trials are terminated when a patient loses balance. In these situations, the ES is not observable, but is assigned the lowest possible score--zero. As a result, the response variable has a mixed discrete-continuous distribution, further compromising inference obtained by standard statistical methods. Here, we develop alternative methodology for analyzing ESs under a stochastic model extending the ES to a continuous latent random variable that always exists, but is unobserved in the event of a fall. Loss of balance occurs conditionally, with probability depending on the realized latent ES. After fitting the model by a form of quasi-maximum-likelihood, one may perform statistical inference to assess the effects of explanatory variables. An example is provided, using data from the NIH/NIA Baltimore Longitudinal Study on Aging.

NASA Discipline Neuroscience↗