Search NASA⌕ Search

SEARCH · Search NASA

Results for “statistical sampling techniques”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Designing an Optimal Sensor Network via Minimizing Information Loss

Optimal experimental design is a classic topic in statistics, with many well-studied problems, applications, and solutions. The design problem we study is the placement of sensors to monitor spatiotemporal processes, explicitly accounting for the temporal dimension in our modeling and optimization. We observe that recent advancements in computational sciences often yield large datasets based on physics-based simulations, which are rarely leveraged in experimental design. We introduce a novel model-based sensor placement criterion, along with a highly-efficient optimization algorithm, which integrates physics-based simulations and Bayesian experimental design principles to identify sensor networks that “minimize information loss” from simulated data. Our technique relies on sparse variational inference and (separable) Gauss-Markov priors, and thus may adapt many techniques from Bayesian experimental design. We validate our method through a case study monitoring air temperature in Phoenix, Arizona, using state-of-the-art physics-based simulations. Our results show our framework to be superior to random or quasi-random sampling, particularly with a limited number of sensors. We conclude by discussing practical considerations and implications of our framework, including more complex modeling tools and real-world deployments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

The development and application of the stirred‐reactor coupon analysis (SRCA) test method

A new technique, termed the stirred‐reactor coupon analysis (SRCA) method, has been developed to measure the rate of glass dissolution in forward‐rate conditions. Monolithic glass coupons are partially masked with an inert material before placement in a large volume of well‐mixed solution with known chemistry and temperature for a predetermined duration. After the test, the mask is removed, and the difference in step height between the protected area and the exposed corroded portions of the sample coupon is measured to determine the extent of glass dissolution. The step height is converted to a rate measurement using the test duration and glass density. Test parameters such as sample surface preparation and test duration were evaluated to determine their effects on the measured rates. Additionally, results from an interlaboratory study (ILS) consisting of 12 laboratories from 11 different institutions are presented, where each laboratory performed 12 independent tests. When removing experimental outlier data, the 95% reproducibility limits for the SRCA method has no statistical difference with previously published standardized test methods used to determine the forward rate of glass dissolution. Overall, this paper describes steps necessary to perform the test method and provides the statistical calculations to evaluate test accuracy.

chemical durability↗
NASA NTRS2000

Effect of scanline orientation on ventricular flow propagation: assessment using high frame-rate color Doppler echocardiography

Color M-mode echocardiography has recently been utilized to describe diastolic flow propagation velocity (Vp) in the left ventricle. While increasing temporal resolution from 15 to 200 Hz, this M-mode technique requires the user to select a single scanline, potentially limiting quantification of Vp due to the complex three-dimensional inflow pattern. We previously performed computational fluid dynamics simulations to demonstrate the insignificance of the scanline orientation, however geometric complexity was limited. The purpose of this study was to utilize high frame-rate 2D color Doppler images to investigate the importance of scanline selection in patients for the quantification of Vp. 2D color Doppler images were digitally acquired at 50 frames/s in 6 subjects from the apical 4-chamber window (System 5, GE/Vingmed, Milwaukee, WI). Vp was determined for a set of scanlines positioned through 5 locations across the mitral annulus (from the anterior to posterior mitral annulus). An analysis of variance was performed to examine the differences in Vp as a function of scanline position. Vp was not effected by scanline position in sampled locations from the center of the mitral valve towards the posterior annulus. Although not statistically significant, there was a trend to slower propagation velocities on the anterior side of the valve (60.8 +/- 16.7 vs. 54.4 +/- 13.6 cm/s). This study clinically validates our previous numerical experiment showing that Vp is insensitive to small perturbations of the scanline through the mitral valve. However, further investigation is necessary to examine the impact of ventricular geometry in pathologies including dilated cardiomyopathy.

Non-NASA Center↗
DOE OSTI2025

The DESI One-Percent Survey: Modelling the clustering and halo occupation of all four DESI tracers with U CHUU

We present results from a set of mock lightcones for the DESI One-Percent Survey, created from the UCHUU simulation. This 8 h −3 Gpc 3 N-body simulation comprises 2.1 trillion particles and provides high-resolution dark matter (sub)haloes in the framework of the Planck-based ΛCDM cosmology. Employing the subhalo abundance matching (SHAM) technique, we populated the UCHUU (sub)haloes with all four DESI tracers – Bright Galaxy Survey (BGS), luminous red galaxies (LRGs), emission line galaxies (ELGs), and quasars (QSOs) – to z = 2.1. Our method accounts for redshift evolution as well as the clustering dependence on luminosity and stellar mass. The two-point clustering statistics of the DESI One-Percent Survey generally agree with predictions from UCHUU across scales ranging from 0.3 h −1 Mpc to 100 h −1 Mpc for the BGS and across scales ranging from 5 h −1 Mpc to 100 h −1 Mpc for the other tracers. We observed some differences in clustering statistics that can be attributed to incompleteness of the massive end of the stellar mass function of LRGs, our use of a simplified galaxy-halo connection model for ELGs and QSOs, and cosmic variance. We find that at the high precision of UCHUU, the shape of the halo occupation distribution (HOD) of the BGS and LRG samples is smaller bias values, likely due to cosmic variance. The bias dependence on absolute magnitude, stellar mass, and redshift aligns with that of previous surveys. These results provide DESI with tools to generate high-fidelity lightcones for the remainder of the survey and enhance our understanding of the galaxy-halo connection.

cosmology↗
NASA NTRS2018

Radar Measurements of the LEO Orbital Debris Environment

Access to space and the preservation of the near-Earth space environment is of critical significance. Increased interest in issues surrounding space traffic management and the continued assessment and discussion of orbital debris at the United Nations Committee on the Peaceful Uses of Outer Space (COPUOS) illustrates the significance of the topic of orbital debris. There are currently over 20,000 tracked objects in the publicly available satellite catalog on Space-Track.org. The catalog is maintained by the US Air Force Space Command using a network of optical and radar ground-based sensors and is believed to be complete for a characteristic size of 10 cm or larger in low Earth orbit (LEO). Based on the work of the NASA Orbital Debris Program Office (ODPO) over approximately the past 40 years it is understood that the small debris population (< 10 cm) increases by orders of magnitude as characteristic size decreases. With population estimates ranging from 500,000 to 1,000,000 small debris objects (> 5 mm) in orbit, it is currently not practical to track and maintain precision orbits on every object. Instead the NASA ODPO uses powerful ground-based radars to sample the low Earth Orbit (LEO) environment and assign approximate orbits to each detection. This poses an interesting signal processing challenge as we are trying to detect the smallest objects possible on the edge of the radar's sensitivity. For approximately the last 30 years, NASA ODPO has partnered with the Massachusetts Institute of Technology Lincoln Laboratory (MIT/LL) to utilize the Haystack Ultra-wideband Satellite Imaging Radar (HUSIR - formerly the Long-Range Imaging Radar or simply Haystack) and the Haystack Auxiliary (HAX) radar to collect orbital debris radar data. Additionally, the ODPO collaborates with the NASA Jet Propulsion Laboratory (JPL) to use the Goldstone Solar System Radar. The orbital debris detections from these radars serve as inputs for statistical risk models used by the human spaceflight and satellite communities to assess risk to spacecraft posed by orbital debris. In this paper, we will describe the history of orbital debris radar measurements conducted by NASA, provide an overview of current radar measurements techniques and facilities, discuss the signal processing software used for orbital debris measurements and the inference of debris size and orbital parameters from these measurements, and discuss how orbital debris radar measurements are validated for use in models that are used throughout the aerospace industry.

Blackwell, Chris↗
DOE OSTI2024

Probabilistic inference of the structure and orbit of Milky Way satellites with semi-analytic modelling

Semi-analytic modelling furnishes an efficient avenue for characterizing dark matter haloes associated with satellites of Milky Way-like systems, as it easily accounts for uncertainties arising from halo-to-halo variance, the orbital disruption of satellites, baryonic feedback, and the stellar-to-halo mass (SMHM) relation. We use the SatGen semi-analytic satellite generator, which incorporates both empirical models of the galaxy–halo connection as well as analytic prescriptions for the orbital evolution of these satellites after accretion onto a host to create large samples of Milky Way-like systems and their satellites. By selecting satellites in the sample that match observed properties of a particular dwarf galaxy, we can infer arbitrary properties of the satellite galaxy within the cold dark matter paradigm. For the Milky Way’s classical dwarfs, we provide inferred values (with associated uncertainties) for the maximum circular velocity v max and the radius r max at which it occurs, varying over two choices of baryonic feedback model and two prescriptions for the SMHM relation. While simple empirical scaling relations can recover the median inferred value for v max and r max , this approach provides realistic correlated uncertainties and aids interpretability. We also demonstrate how the internal properties of a satellite’s dark matter profile correlate with its orbit, and we show that it is difficult to reproduce observations of the Fornax dwarf without strong baryonic feedback. Furthermore, the technique developed in this work is flexible in its application of observational data and can leverage arbitrary information about the satellite galaxies to make inferences about their dark matter haloes and population statistics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Robust Design Under Uncertainty in Quantum Error Mitigation

Error mitigation techniques are crucial to achieving near-term quantum advantage. Classical postprocessing of quantum computation outcomes is a popular approach for error mitigation, which includes methods, such as zero noise extrapolation, virtual distillation, and learning-based error mitigation. However, these techniques have limitations due to the propagation of uncertainty resulting from the finite shot number of a quantum measurement. In this work, we introduce general and unbiased methods for quantifying the uncertainty and error of error-mitigated observables based on the strategic sampling of error mitigation outcomes. We then extend our approach to demonstrate the optimization of performance and robustness of error mitigation under uncertainty. To illustrate our methods, we apply them to zero noise extrapolation and Clifford date regression in the ground state of the XY model simulated using depolarizing and International Business Machines Corporation (IBM) Toronto noise models, respectively. In particular, we optimize the choice of noise levels and the allocation of shots for zero noise extrapolation and the distribution of the training circuits for Clifford data regression. While our methods are readily applicable to any postprocessing-based error mitigation approach, in practice they must not be prohibitively expensive—even though they perform optimizations of the error mitigation hyperparameters requiring sampling of a statistical distribution of error mitigation outcomes. By leveraging surrogate-based optimization, we show that our methods can efficiently perform optimal design for a zero noise extrapolation implementation. We then further demonstrate the transferability of learned zero noise extrapolation hyperparameters to other similar circuits.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Sampling Size Optimization for Bioburden Density Estimation in Planetary Protection

Planetary protection (PP) is a discipline that focuses on minimizing the biological contamination of spacecraft to ensure compliance with international policy. Precise estimation of bioburden - the total number of microbes in or on spacecraft hardware – and the bioburden density are of utmost importance for PP. Such estimation is the way concordance with requirements is demonstrated, and it is critical for quantifying the potential risk of inadvertently contaminating other planetary bodies. Although a suite of molecular techniques have been used to thoroughly characterize and profile the microbiome of various cleanroom environments and spacecraft, the gold standard remains the physical enumeration of microbes via culturing of samples directly taken from spacecraft and associated surfaces. However, due to technical, budgetary, and programmatic constraints, only a manageable portion (around 10%) of the entire spacecraft surface is directly sampled with cotton swabs or wipes. To generate the bioburden current best estimate (CBE) for components not directly verifiable, the accepted approach is to apply a NASA-defined bioburden estimate based on the components’ manufacturing or assembly environment. This approach utilizes a prespecified bioburden density estimation that applies a maximum value across the total surface area of the specified component. For hardware components that underwent similar assembly processes, an implied bioburden is adopted for all components, based on a direct verification of a representative component within the same lot. Once all components have a CBE, the bioburden estimates are generated. In previous publication [ 1], we have shown that statistical risks quantifying the accuracy of the estimates for sampled, prespecified, and implied components can be derived and ranked. For mean squared error (MSE) function, the risks are available analytically and hence a cost function can be obtained to optimize the risks with respect to the sampling area and sampling cost. Since the sampling area and sampling cost are two complimentary variables, their sum will have a well-defined minimum. This paper presents the multivariate optimization of the integrated risk of an empirical Bayes estimator to determine the optimal sampling schedule for a given number of components. It is assumed that given a number of components, N, the bioburden density for each component can either be sampled, implied, or prespecified. The multivariate optimization searches through different options to sample, imply or prespecify the bioburden density for a component, and account for the component’s surface area and cost of sampling. The idea of the optimization is based on the observation that the statistical risk of using an estimator is a monotonically decreasing function of the sampled area. The larger the sampled area, the lower the risk of using the estimator as the estimator becomes more and more accurate as the sampling area increases. On the other hand, the cost of sampling is monotonically increasing as the sampled surface grows. This makes the risk and total cost of sampling complimentary variables which can be counterbalanced to achieve an optimal overall value with respect to the sampled surface. In this paper, the integrated risk has been used to quantify the accuracy of the estimator. This risk has been selected because it depends on neither the true value of the parameter nor on the collected data. The cost of each sample was also available to obtain the total cost of sampling of N components. The paper will present the results based on computer-simulated data as well as the data collected during the InSight mission. The computer-simulated data have N components with randomly generated total areas and each component assigned to one of the three categories according to the method of estimating of bioburden density: sampled, implied, or prespecified. The cost of sampling is also available. The cost of sampling is estimated based on a cost model provided by the planetary protection group at JPL. For this paper, the overall cost was assumed to be a linear function of exposure. The optimization process finds the allocation of the components to the three categories that minimizes the tradeoff between integrated risk and total cost. For the InSight data, a set of components is selected representing all three categories, and optimization is performed to determine if the performed allocation was optimal or if a better allocation could have been obtained. To the best of our knowledge, this work is the first attempt not only perform an accurate estimation of bioburden density but also do it in an optimal way.

97 - MATHEMATICS AND COMPUTING↗
NASA NTRS2016

Characterization of Carbon Dioxide Washout Measurement Techniques in the Mark-III Space Suit

A space suit must provide adequate carbon dioxide (CO2) washout inside the helmet to prevent symptoms of hypercapnia. In the past, an oronasal mask has been used to measure the inspired air of suited subjects to determine a space suit's CO2 washout capability. While sufficient for super-ambient pressure testing of space suits, the oronasal mask fails to meet several human factors and operational criterion needed for future sub-ambient pressure testing (e.g. compatibility with a Valsalva device). This paper describes the evaluation of a nasal cannula as a device for measuring inspired air within a space suit. Eight test subjects were tasked with walking on a treadmill or operating an arm ergometer to achieve target metabolic rates of 1000, 2000, and 3000 British thermal units per hour (BTU/hr), at flow rates of 2, 4, and 6 actual cubic feet per minute (ACFM). Each test configuration was conducted twice, with subjects instructed to breathe either through their nose only, or however they felt comfortable. Test data shows that the nasal cannula provides more statistically consistent data across test subjects than the oronasal mask used in previous tests. The data also shows that inhaling/exhaling through only the nose provides a lower sample variance than a normal breathing style. Nose-only breathing reports better CO2 washout due to several possible reasons, including a decreased respiratory rate, an increased tidal volume, and because nose-only breathing directs all of the exhaled CO2 down and away from the oronasal region. The test subjects in this study provided feedback that the nasal cannula is comfortable and can be used with the Valsalva device.

Meginnis, Ian M.↗
DOE OSTI2025

Evaluating User Errors and Temporal Trends in Marine Fish Communities Using 360-Degree Underwater Photography

The use of environmental DNA (eDNA) sampling has been proposed as a complementary method to monitor fish species in marine environments, offering a non-invasive and potentially more efficient approach to marine species observations. eDNA monitoring could be especially useful in and around sites targeted for marine energy generation as these regions need regular monitoring that would be impractical with traditional techniques. Before we can fully rely upon eDNA, we must first verify its accuracy against other proven methods, such as the use of underwater photography. In this study, I deployed a 360-degree camera in the tidal channel of Sequim Bay once a month during several hours overlapping slack tide. I investigated how having multiple people identify and count fish on underwater images could affect the overall results. Using chi square tests in R, I compared my fish identifications and counts to those made by another intern on the same images recorded in August. I found significant differences in the number of species identified and the total individual counts between the two different datasets. I also tested the statistical differences in both Shannon diversity and Pielou evenness indices between the August, September, and November camera deployments using a Hutcheson t-test. Only one significant difference was found in the Shannon index comparisons, and none were found between the Pielou evenness comparisons. These findings show that if multiple identifiers are used to process underwater images, quality control checks must be made to reduce the potential for error. This also points toward the possibility to leverage more advanced image analysis processes, such as automated image analysis software. The findings from this study also show that the dynamics of marine fish communities can vary over a few months; however, further analysis is needed to determine the extent of the seasonal changes in Sequim Bay.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2011

Standardization of a Volumetric Displacement Measurement for Two-Body Abrasion Scratch Test Data Analysis

A limitation has been identified in the existing test standards used for making controlled, two-body abrasion scratch measurements based solely on the width of the resultant score on the surface of the material. A new, more robust method is proposed for analyzing a surface scratch that takes into account the full three-dimensional profile of the displaced material. To accomplish this, a set of four volume- displacement metrics was systematically defined by normalizing the overall surface profile to denote statistically the area of relevance, termed the Zone of Interaction. From this baseline, depth of the trough and height of the plowed material are factored into the overall deformation assessment. Proof-of-concept data were collected and analyzed to demonstrate the performance of this proposed methodology. This technique takes advantage of advanced imaging capabilities that allow resolution of the scratched surface to be quantified in greater detail than was previously achievable. When reviewing existing data analysis techniques for conducting two-body abrasive scratch tests, it was found that the ASTM International Standard G 171 specified a generic metric based only on visually determined scratch width as a way to compare abraded materials. A limitation to this method was identified in that the scratch width is based on optical surface measurements, manually defined by approximating the boundaries, but does not consider the three-dimensional volume of material that was displaced. With large, potentially irregular deformations occurring on softer materials, it becomes unclear where to systematically determine the scratch width. Specifically, surface scratches on different samples may look the same from a top view, resulting in an identical scratch width measurement, but may vary in actual penetration depth and/or plowing deformation. Therefore, two different scratch profiles would be measured as having identical abrasion properties, although they differ significantly.

Street, K. W. Jr.↗
DOE OSTI2025

An In Situ , Automated High-Explosives Aging Method Utilizing Two-Dimensional Gas Chromatography–Mass Spectrometry

Understanding chemical changes that occur in high explosives as they age is of great importance to the safe employment and storage of these compounds. Traditional methods of aging high explosives even under accelerated aging conditions are time intensive with durations on the order of months to years. The nature of traditional aging analyses reduces each sample to a snapshot data point often separated widely in time, requiring many assumptions as to how the degradation products develop. Further complicating matters, several analytical techniques are typically employed for each sample analysis in order to ascertain an entire picture of the decomposition pathways. To address these shortcomings with existing methods, a new method of accelerated aging of high explosives utilizing comprehensive two-dimensional gas chromatography coupled to high-resolution mass spectrometry (GC × GC-HRMS) was developed using 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (CL-20) as a model compound for method development. This in situ automated method reduces the time scale of aging to a matter of hours using the inlet of the GC × GC as the aging vessel. GC × GC in combination with HRMS allowed for the collection of both evolved gases and other decomposition products produced during the entire aging process in real time with HRMS providing far greater certainty in identification of explosives aging products. Additionally, this method allowed for a higher throughput of samples with greatly simplified sample preparation. Chemometric analysis of the GC × GC-HRMS data set via the alteration analysis (ALA) enabled discovery of statistically significant chemical changes providing insight into the variation of decomposition pathways with varying aging temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2002

Transit Detection with a Distributed Network of Telescopes

The discovery since 1995 of more than 80 planets around nearby solar-like stars and the photometric detection of a transit of the planet orbiting HD 209458 (producing a more than 1% drop in brightness that lasts 3 hours) has heralded a new era in astronomy. It has now been demonstrated that small telescopes equipped with sensitive and stable electronic detectors can produce fundamental scientific discoveries regarding the frequency and nature of planets outside the solar system. The modest equipment requirements for the measurement of extrasolar planetary transits are achieved by commercial small aperture telescopes and CCD imagers common among amateur astronomers. With equipment already in hand and armed with target lists, observing techniques and software procedures developed b NASA's Ames Research Center and the University of California at Santa Cruz, non-professional astronomers can contribute significantly to the study of planets around others stars. Statistical analyses of the population of parent stars of the known extrasolar planets indicate that approximately one in ten metal-rich stars should harbor a short-period planet. Given the ten percent chance that a given short-period planet displays transits, we therefore expect that approximately 1% of the most metal rich stars will have a planetary companion detectable by this project. A catalog of 206 highly metal rich nearby F, G and K stars has been compiled, and this catalog will provide a rich source of targets. In addition, main sequence F, G, K and M stars identified to have "transit-like" features in the Hipparcos satellite photometry archive will also be monitored. A commercially available "amateur grade" telescope/CCD/software system acquired late during the 2001 "transit season" for HID 209458 has achieved 0.47% RMS precision for 13 minute time sampling from a suburban backyard under less than ideal observing conditions and a realistic range of airmass values.

Castellano, T.↗
NASA NTRS2003

Vector-Ordering Filter Procedure for Data Reduction

The vector-ordering filter (VOF) technique involves a procedure for sampling a large population of data vectors to select a subset of data vectors that fully characterize the state space of the large population. The VOF technique enables a large reduction of the volume of data that must be handled in the automated monitoring system and method discussed in the two immediately preceding articles. In so doing, the VOF technique enables the development of data-driven mathematical models of a monitored asset from sets of data that would otherwise exceed the memory capacities of conventional engineering computers. Data-driven mathematical models have been shown to offer high fidelity for purposes of control and monitoring of assets. In practice, a collection of asset-operating observations is acquired with the intention that the collection contain observations characteristic of the full dynamic range of operation of the asset. Often, such a collection contains an extremely large number of observations, many of which are redundant. The VOF technique fills the need for a means to extract, from the original collection of observational data, a reduced data matrix that excludes redundant data while maintaining the full statistical character and dynamic range of the original data. The reduced data matrix can then be used as the input data for development of a mathematical model of the monitored asset, or as training data for a neural-network substitute for an explicit mathematical model of the asset. Alternatively, the reduced data matrix can, itself, be used directly as a mathematical model of the monitored asset, as is commonly done in multivariate state-estimation techniques. The original data are collected from the asset over a range of operating states and are put in matrix form. Each column vector in the original data matrix represents the signal values acquired at a particular operational state of the asset. Thus, the number of columns of the original data matrix equals the number of observed states and the number of rows in this matrix equals the number of signals acquired at each observation. In the VOF technique, one extracts the reduced data matrix from the original data matrix through the selection of a representative subset of the column (state) vectors.

Bickford, Randall L.↗
NASA NTRS1999

The NASA/MSFC Global Reference Atmospheric Model: 1999 Version (GRAM-99)

The latest version of Global Reference Atmospheric Model (GRAM-99) is presented and discussed. GRAM-99 uses either (binary) Global Upper Air Climatic Atlas (GUACA) or (ASCII) Global Gridded Upper Air Statistics (GGUAS) CD-ROM data sets, for 0-27 km altitudes. As with earlier versions, GRAM-99 provides complete geographical and altitude coverage for each month of the year. GRAM-99 uses a specially-developed data set, based on Middle Atmosphere Program (MAP) data, for 20-120 km altitudes, and NASA's 1999 version Marshall Engineering Thermosphere (MET-99) model for heights above 90 km. Fairing techniques assure smooth transition in overlap height ranges (20-27 km and 90-120 km). GRAM-99 includes water vapor and 11 other atmospheric constituents (O3, N2O, CO, CH4, CO2, N2, O2, O, A, He and H). A variable-scale perturbation model provides both large-scale (wave) and small-scale (stochastic) deviations from mean values for thermodynamic variables and horizontal and vertical wind components. The small-scale perturbation model includes improvements in representing intermittency ("patchiness"). A major new feature is an option to substitute Range Reference Atmosphere (RRA) data for conventional GRAM climatology when a trajectory passes sufficiently near any RRA site. A complete user's guide for running the program, plus sample input and output, is provided. An example is provided for how to incorporate GRAM-99 as subroutines in other programs (e.g., trajectory codes).

Justus, C. G.↗
NASA NTRS1997

Ultrasonic wave velocity measurement in small polymeric and cortical bone specimens

A system was refined for the determination of the bulk ultrasonic wave propagation velocity in small cortical bone specimens. Longitudinal and shear wave propagations were measured using ceramic, piezoelectric 20 and 5 MHz transducers, respectively. Results of the pulse transmission technique were refined via the measurement of the system delay time. The precision and accuracy of the system were quantified using small specimens of polyoxymethylene, polystyrene-butadiene, and high-density polyethylene. These polymeric materials had known acoustic properties, similarity of propagation velocities to cortical bone, and minimal sample inhomogeneity. Dependence of longitudinal and transverse specimen dimensions upon propagation times was quantified. To confirm the consistency of longitudinal wave propagation in small cortical bone specimens (< 1.0 mm), cut-down specimens were prepared from a normal rat femur. Finally, cortical samples were prepared from each of ten normal rat femora, and Young's moduli (Eii), shear moduli (Gij), and Poisson ratios (Vij) were measured. For all specimens (bone, polyoxymethylene, polystyrene-butadiene, and high-density polyethylene), strong linear correlations (R2 > 0.997) were maintained between propagation time and distance throughout the size ranges down to less than 0.4 mm. Results for polyoxymethylene, polystyrene-butadiene, and high-density polyethylene were accurate to within 5 percent of reported literature values. Measurement repeatability (precision) improved with an increase in the wave transmission distance (propagating dimension). No statistically significant effect due to the transverse dimension was detected.

Non-NASA Center↗
NASA NTRS2022

GNSS-RO Deep Refraction Signals from Moist Marine Atmospheric Boundary Layer (MABL)

The marine atmospheric boundary layer (MABL) has a profound impact on sensible heat and moisture exchanges between the surface and the free troposphere. The goal of this study is to develop an alternative technique for retrieving MABL-specific humidity (q) using GNSS-RO data in deep-refracted signals. The GNSS-RO signal amplitude (i.e., signal-to-noise ratio or SNR) at the deep straight-line height (H(SL)) was been found to be strongly impacted by water vapor within the MABL. This study presents a statistical analysis to empirically relate the normalized SNR (S(RO)) at deep HSL to the MABL q at 950 hPa (~400 m). When compared to the ERA5 reanalysis data, a good linear q–S(RO) relationship is found with the deep H(SL) S(RO) data, but careful treatments of receiver noise, SNR normalization, and receiver orbital altitude are required. We attribute the good q–S(RO) correlation to the strong refraction from a uniform, horizontally stratiform and dynamically quiet MABL water vapor layer. Ducting and diffraction/interference by this layer help to enhance the S(RO) amplitude at deep H(SL). Potential MABL water vapor retrieval can be further developed to take advantage of a higher number of S(RO) measurements in the MABL compared to the Level-2 products. A better sampled diurnal variation of the MABL q is demonstrated with the S(RO) data over the Southeast Pacific (SEP) and the Northeast Pacific (NEP) regions, which appear to be consistent with the low cloud amount variations reported in previous studies.

diurnal variation↗
NASA NTRS1991

On the estimation of climatological Z-R relationships

A statistical framework for climatological Z-R parameter estimation is developed and simulation experiments are conducted to examine sampling properties of the estimators. Both parametric and nonparametric models are considered. For parametric models, it is shown that Z-R parameters can be estimated by maximum likelihood, a procedure with optimal large sample properties. A general nonparametric framework for climatological Z-R estimation is also developed. Nonparametric procedures are attractive because of their flexibility in dealing with certain types of measurement errors common to radar data. Simulation experiments show that even under favorable assumptions on error characteristics of radar and raingages, large datasets are required to obtain accurate Z-R parameter estimates. Another important conclusion is that estimation results are generally quite sensitive to radar and raingage measurement thresholds. For fixed sample size, the simulation results can be used to provide quantitative assessments of the accuracy of Z-R model parameter estimates. These results are particularly useful for error analysis of precipitation products that are derived using climatological Z-R relations. One example is the large-area rainfall estimates derived using the height-area rainfall threshold (HART) technique.

Krajewski, Witold F.↗