Results for “parameter inference” — Page 6

DOE OSTI2025

MINE: a new way to design genetics experiments for discovery

Abstract The Maximally Informative Next Experiment or MINE is a new experimental design approach for experiments, such as those in omics, in which the number of effects or parameters p greatly exceeds the number of samples n (p > n). Classical experimental design presumes n > p for inference about parameters and its application to p > n can lead to over-fitting. To overcome p > n, MINE is an ensemble method, which makes predictions about future experiments from an existing ensemble of models consistent with available data in order to select the most informative next experiment. Its advantages are in exploration of the data for new relationships with n < p and being able to integrate smaller and more tractable experiments to replace adaptively one large classic experiment as discoveries are made. Thus, using MINE is model-guided and adaptive over time in a large omics study. Here, MINE is illustrated in two distinct multiyear experiments, one involving genetic networks in Neurospora crassa and a second one involving a genome-wide association study in Sorghum bicolor as a comparison to classic experimental design in an agricultural setting.

Biochemistry & Molecular Biology↗

NASA NTRS2019

A NICER view of PSR J0030+0541: Millisecond Pulsar Parameter Estimation

We report on Bayesian parameter estimation of the mass and equatorial radius of the millisecond pulsar PSRJ0030+0451, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray spectral timing event data. We perform relativistic ray-tracing of thermal emission from hot regions of the pulsar’s surface. We assume two distinct hot regions based on two clear pulsed components in the phase-folded pulse-profile data; we explore a number of forms (morphologies and topologies) for each hot region, inferring their parameters in addition to the stellar mass and radius. For the family of models considered, the evidence (prior predictive probability of the data) strongly favors a model that permits both hot regions to be located in the same rotational hemisphere. Models wherein both hot regions are assumed to be simply connected circular single-temperature spots, in particular those where the spots are assumed to be reflection-symmetric with respect to the stellar origin, are strongly disfavored. For the inferred configuration, one hot region subtends an angular extent of only a few degrees (in spherical coordinates with origin at the stellar center) and we are insensitive to other structural details; the second hot region is far more azimuthally extended in the form of a narrow arc, thus requiring a larger number of parameters to describe. The inferred mass M and equatorial radius Req are, respectively, -1.34+ M 0.16 0.15 and-12.71+ km 1.19 1.14 , while the compactness = -GM R c 0.156+ eq 2 0.010 0.008 is more tightly constrained; the credible interval bounds reported here are approximately the 16% and 84% quantiles in marginal posterior mass.

T E Riley↗

DOE OSTI2025

Generalized fiducial inference on differentiable manifolds

We introduce a novel approach to inference on parameters that take values in a Riemannian manifold embedded in a Euclidean space. Parameter spaces of this form are ubiquitous across many fields, including chemistry, physics, computer graphics, and geology. Here, this new approach uses generalized fiducial inference (GFI) to obtain a posterior-like distribution on the manifold, without needing to know local parameterizations that map to the constrained space from an unconstrained Euclidean space. Using mathematical tools from Riemannian geometry, we construct a constrained generalized fiducial distribution (CGFD). A Bernstein-von Mises-type result for the CGFD, which provides intuition for how the desirable asymptotic qualities of the unconstrained generalized fiducial distribution are inherited by the CGFD, is provided. To illustrate the practical use of the CGFD, we provide a proof-of-concept example in the context of a linear logspline density estimation problem, and demonstrate that CGFD-based confidence sets exhibit desirable coverage properties via simulation. As an application, we fit a CGFD to COVID-19 case count data from North Carolina, USA.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2022

Optimal Estimation Framework for Ocean Color Atmospheric Correction and Pixel-level Uncertainty Quantification

Ocean color remote sensing requires compensation for atmospheric scattering and absorption (aerosol, Rayleigh, and trace gases), referred to as atmospheric correction (AC). AC allows inference of parameters such as spectrally resolved remote sensing reflectance ( R rs )(λ) ; sr 1 ) at the ocean surface from the top-of-atmosphere reflectance. Often, the uncertainty of this process is not fully explored. Bayesian inference techniques provide a simultaneous AC and uncertainty assessment via a full posterior distribution of the relevant variables, given the prior distribution of those variables and the radiative transfer (RT) likelihood function. Given uncertainties in the algorithm inputs, the Bayesian framework enables better constraints on the AC process by using the complete spectral information compared to traditional approaches that use only a subset of bands for AC. This paper investigates a Bayesian inference research method (Optimal Estimation, OE) for ocean color AC by simultaneously retrieving atmospheric and ocean properties using all visible and near-infrared spectral bands. The OE algorithm analytically approximates the posterior distribution of parameters based on normality assumptions and provides a potentially viable operational algorithm with a reduced computational expense. We developed a Neural Network (NN) RT forward model look-up-table-based emulator to increase algorithm efficiency further and thus speed up the likelihood computations. We then applied the OE algorithm to synthetic data and observations from the MODerate resolution Imaging Spectroradiometer (MODIS) on NASA’s Aqua spacecraft. We compared the R rs )(λ) retrieval and its uncertainty estimates from the OE method with in-situ validation data from the SeaWiFS Bio-optical Archive and Storage System (SeaBASS) and Aerosol Robotic Network Ocean Color (AERONET-OC) datasets. The OE algorithm improved R rs )(λ) estimates relative to the NASA standard operational algorithm by improving all statistical metrics at 443, 555, and 667 nm. Unphysical negative R rs )(λ) , which often appear in complex water conditions, was reduced by a factor of 3. The OE-derived pixel-level R rs )(λ) uncertainty estimates were also assessed relative to in-situ data and were shown to have skill.

Atmospheric correction↗

NASA NTRS1998

Mapping the Solar Wind from its Source Region into the Outer Corona

Knowledge of the radial variation of the plasma conditions in the coronal source region of the solar wind is essential to exploring coronal heating and solar wind acceleration mechanisms. The goal of the present proposal is to determine as many plasma parameters in that region as possible by coordinating different observational techniques, such as Interplanetary Scintillation Observations, spectral line intensity observations, polarization brightness measurements and X-ray observations. The inferred plasma parameters are then used to constrain solar wind models.

Esser, Ruth↗

NASA NTRS1981

Inference of sea surface temperature, near surface wind, and atmospheric water by Fourier analysis of Scanning Multichannel Microwave Radiometer data

The Scanning Multichannel Microwave Radiometer measures thermal microwave emission from the earth in both polarizations at wavelengths of 0.8, 1.4, 1.7, 2.8 and 4.6 cm. Similar instruments were launched on Nimbus 7 and Seasat. Both spatial resolution on the earth and relative sensitivity to different geophysical parameters change with wavelength. Therefore, spatial Fourier components of geophysical parameters are inferred from the corresponding Fourier components of the radiometer measurements, taking into account the different dependence of signal-to-noise ratio on spatial frequency for each radiometer wavelength. The geophysical parameters are sea surface temperature, near-surface wind speed, integrated water vapor mass, integrated liquid water mass, and the product of rainfall rate with height of the rain layer. The capabilities and limitations of the inversion method are illustrated by means of data from the North Atlantic and from tropical storms.

Rosenkranz, P. W.↗

NASA NTRS1997

Mapping the Solar Wind from its Source Region into the Outer Corona

Knowledge of the radial variation of the plasma conditions in the coronal source region of the solar wind is essential to exploring coronal heating and solar wind acceleration mechanisms. The goal of the proposal was to determine as many plasma parameters in the solar wind acceleration region and beyond as possible by coordinating different observational techniques, such as Interplanetary Scintillation Observations, spectral line intensity observations, polarization brightness measurements and X-ray observations. The inferred plasma parameters were then used to constrain solar wind models.

Esser, Ruth↗

NASA NTRS1995

Relating Radar Backscatter to Boreal Forest Canopy Parameters

During the first intensive field campaign for the Boreal Ecosystem Atmospheric Study (BOREAS) in August of 1993, AIRSAR data were acquired over the entire BOREAS study area. On of the objectives of the AIRSAR deployment was to examine the sensitivity of the radar signal to vegetation type and biomass distribution in the boreal forest and to develop algorithms for inferring vegetation parameters. During the experiment a set of ground measurements were also made to support the AIRSAR data analysis. The dominant stands in the study area consist of black spruce, young jack pine, old jack pine, and aspen. These stands represent a wide range of biomass and canopy architectural variations which can be distinguished in SAR images.

radar backscatter↗

NASA NTRS1997

Gaseous abundances and methane supersaturation in Titan's troposphere

Various properties of Titan's troposphere are inferred from an analysis of Voyager I infrared spectrometer (IRIS) data between 200 and 600/ cm. Two homogeneous spectral averages acquired at widely separated emission angles are chosen for the analysis. Both data sets are associated with northern low latitudes very close to that of the radio science ingress occultation point. Solutions require simultaneous nonlinear least-squares fits to the two IRIS data sets, coupled with iteration of the radio occultation refractivity data. Values and associated 1-sigma uncertainties of several parameters are inferred from our analysis. These include mole fractions for molecular hydrogen (approx. 0.0011), argon (small), and methane near the surface ( approx. 0.057). Solutions are also obtained for the hydrogen parafraction (close to equilibrium, with considerable uncertainty), air temperature near the surface ( approx. 93 K), surface surface temperature discontinuity (approx. 1 K), and maximum degree of methane supersaturation in the upper troposphere (approx. 1.5). Actual values for the above-mentioned parameters depend on the amount of ethane cloud near the tropopause. There is no evidence for methane clouds in the upper troposphere, nor is their presence compatible with large degrees of supersaturation. A wave number dependence for the stratospheric haze opacity is inferred similar to that found for a polymeric residue created in laboratory discharge experiments. This haze appears to be uniformly distributed with latitude between altitudes of 40 and 160 km, provided those nighttime data at southern high latitudes that are discounted.

Samuelson, Robert E.↗

NASA NTRS2004

Implicit and Explicit Spacecraft Gyro Calibration

This paper presents a comparison between two approaches to sensor calibration. According to one approach, called explicit, an estimator compares the sensor readings to reference readings, and uses the difference between the two to estimate the calibration parameters. According to the other approach, called implicit, the sensor error is integrated to form a different entity, which is then compared with a reference quantity of this entity, and the calibration parameters are inferred from the difference. In particular this paper presents the comparison between these approaches when applied to in-flight spacecraft gyro calibration. Reference spacecraft rate is needed for gyro calibration when using the explicit approach; however, such reference rates are not readily available for in-flight calibration. Therefore the calibration parameter-estimator is expanded to include the estimation of that reference rate, which is based on attitude measurements in the form of attitude-quaternion. A comparison between the two approaches is made using simulated data. It is concluded that the performances of the two approaches are basically comparable. Sensitivity tests indicate that the explicit filter results are essentially insensitive to variations in given spacecraft dynamics model parameters.

Bar-Itzhack, Itzhack Y.↗

NASA NTRSDate not supplied

Highlights of NASA’s Orbital Debris Program Office In Situ and Laboratory Measurements

NASA’s Orbital Debris Program Office (ODPO) maintains various returned spacecraft materials, capabilities, and facilities used for in situ and laboratory measurements that directly support orbital debris environmental models. In situ measurements include the analysis of exposed and returned hardware surfaces. These surfaces serve as passive sensors for the small-sized micrometeoroid and orbital debris (MMOD) flux below the sensitivity of ground-based radar and optical sensors. Various instruments and techniques are used to determine the size and depth of selected impact features, and – if feasible – the composition of the projectile material. Analysis of the impactor residues enables the differentiation of MM and OD for debris below 1 mm to support modeling the OD environment. In addition, projectiles identified as OD can be further differentiated in low-, medium-, and high- density impactors based on chemical analyses. In addition to in situ measurements, the ODPO has also worked in collaboration with the U.S. Space Force Space Systems Command (formerly the U.S. Air Force Space and Missile Systems Center), the Aerospace Corporation, and the University of Florida on a laboratory-based hypervelocity impact test, DebriSat, conducted at the Air Force Arnold Engineering Development Complex in 2014. The resulting data from this impact test series are being analyzed to assess the fragments’ sizes/masses, materials/densities, shapes, and other parameters of interest. The DebriSat project provides the data needed to update NASA’s breakup models and size estimation models using the simulated orbital breakup of a modern, low Earth orbit spacecraft. Ultimately, over 200,000 fragments from this impact test will be stored at NASA Johnson Space Center (JSC) and further analyzed by the ODPO. This project will also use machine learning techniques to infer physical parameters of fragments embedded in the soft-catch foam used in the impact experiment. Applied to X-ray imagery of the foam panels, these techniques promise to minimize human-in-the-loop processes for fragment extraction and physical characterization. A brief overview of this project and data collected will be presented. Lastly, the ODPO provides various capabilities hosted at NASA JSC for optical inspections and measurements using a variety of techniques and scientific instrumentation to support both in situ and laboratory measurements. The ODPO’s Optical Measurement Center (OMC) is an advanced facility for photometric and spectroscopic laboratory measurements of targets, including fragments from the DebriSat project. The OMC simulates telescopic observations by using space-like illumination conditions and source-target-sensor orientation techniques. Additionally, the OMC is uniquely equipped to acquire pseudo-bidirectional reflectance distribution data for broadband photometric measurements, thus removing aspect angle dependencies that can affect target size estimates using the optical size estimation model. Narrow-band surface material characterization using spectroscopic instrumentation gives insight into how the albedo parameter – also important in the optical size estimation model – may vary depending on the state of the material. The OMC also performs simulations of photometric measurements using optical ray-tracing software to model the OMC optical throughput. In addition to the OMC, the ODPO houses a start-of-the-art Fragment Analysis Facility that uses multiple microscopic inspection instruments to support in situ measurements and material characterization. An overview of both facilities will be highlighted in this paper.

Orbital Debris↗

NASA NTRSDate not supplied

Highlights of NASA’s Orbital Debris Program Office In Situ and Laboratory Measurements

NASA’s Orbital Debris Program Office (ODPO) maintains various returned spacecraft materials, capabilities, and facilities used for in situ and laboratory measurements that directly support orbital debris environmental models. In situ measurements include the analysis of exposed and returned hardware surfaces. These surfaces serve as passive sensors for the small-sized micrometeoroid and orbital debris (MMOD) flux below the sensitivity of ground-based radar and optical sensors. Various instruments and techniques are used to determine the size and depth of selected impact features, and – if feasible – the composition of the projectile material. Analysis of the impactor residues enables the differentiation of MM and OD for debris below 1 mm to support modeling the OD environment. In addition, projectiles identified as OD can be further differentiated in low-, medium-, and high- density impactors based on chemical analyses. In addition to in situ measurements, the ODPO has also worked in collaboration with the U.S. Space Force Space Systems Command (formerly the U.S. Air Force Space and Missile Systems Center), the Aerospace Corporation, and the University of Florida on a laboratory-based hypervelocity impact test, DebriSat, conducted at the Air Force Arnold Engineering Development Complex in 2014. The resulting data from this impact test series are being analyzed to assess the fragments’ sizes/masses, materials/densities, shapes, and other parameters of interest. The DebriSat project provides the data needed to update NASA’s breakup models and size estimation models using the simulated orbital breakup of a modern, low Earth orbit spacecraft. Ultimately, over 200,000 fragments from this impact test will be stored at NASA Johnson Space Center (JSC) and further analyzed by the ODPO. This project will also use machine learning techniques to infer physical parameters of fragments embedded in the soft-catch foam used in the impact experiment. Applied to X-ray imagery of the foam panels, these techniques promise to minimize human-in-the-loop processes for fragment extraction and physical characterization. A brief overview of this project and data collected will be presented. Lastly, the ODPO provides various capabilities hosted at NASA JSC for optical inspections and measurements using a variety of techniques and scientific instrumentation to support both in situ and laboratory measurements. The ODPO’s Optical Measurement Center (OMC) is an advanced facility for photometric and spectroscopic laboratory measurements of targets, including fragments from the DebriSat project. The OMC simulates telescopic observations by using space-like illumination conditions and source-target-sensor orientation techniques. Additionally, the OMC is uniquely equipped to acquire pseudo-bidirectional reflectance distribution data for broadband photometric measurements, thus removing aspect angle dependencies that can affect target size estimates using the optical size estimation model. Narrow-band surface material characterization using spectroscopic instrumentation gives insight into how the albedo parameter – also important in the optical size estimation model – may vary depending on the state of the material. The OMC also performs simulations of photometric measurements using optical ray-tracing software to model the OMC optical throughput. In addition to the OMC, the ODPO houses a start-of-the-art Fragment Analysis Facility that uses multiple microscopic inspection instruments to support in situ measurements and material characterization. An overview of both facilities will be highlighted in this paper.

Orbital Debris↗

NASA NTRS1995

Application of Space Shuttle photography to studies of upper ocean dynamics

Three studies have been conducted using space shuttle imagery to explain the dynamics behavior of internal waves in the Atlantic and Indian Oceans and to derive tide-related parameters for Delaware Bay. By interpreting space shuttle photographs taken during mission STS-40, a total of 34 internal wave packets on the continental shelf of the Middle Atlantic Bight have been recognized. Using the finite-depth theory we derived that the maximum amplitude of solitons is 5.6 m, the phase speed 0.42 m/s, and the period 23.8 min. Deep-ocean internal waves in the western equatorial Indian Ocean on photographs taken during mission STS-44 were also interpreted and analyzed. The internal waves occurred in the form of a multisoliton packet in which there are about a dozen solitons. The average wavelength of the solitons is 1.8 +/- 0.5 km. The crest lines are mostly straight and reach as long as 100 km. The distance between two adjacent packets is about 66 km. Using the deepwater soliton theory, we derived that the mean amplitude of the solitons is 25 m, the nonlinear phase speed 1.7 m/s, and the average period 18 min. For both cases, the semidiural tides are the principal generating mechanism. The tide-related parameters of Delaware Bay were derived from space shuttle time-series photographs taken during mission STS-40. The water area in the bay were measured from interpretation maps of the photographs. The corresponding tidal levels were calculated using the exposure time. From these data, an approximate function relating the water area to the tidal level at a reference point was determined. Then, the water areas of the Delaware Bay at mean high water (MHW) and mean low water (MLW), below 0 m, for the tidal zone, and the tidal flux were inferred. All parameters derived were reasonable and compared well with results of previous investigations.

Zheng, Quanan↗

DOE OSTI2025

Multifrequency Models of Black Hole Photon Rings from Low-luminosity Accretion Disks

Images of black holes encode both astrophysical and gravitational properties. Detecting highly lensed features in images can differentiate between these two effects. We present an accretion disk emission model coupled to the Adaptive Analytical Ray Tracing (AART) code that allows a fast parameter space exploration of black hole photon ring images produced from synchrotron emission from 10 to 670 GHz. As an application, we systematically study several disk models and compute their total flux density, average radii, and optical depth. The model parameters are chosen around fiducial values calibrated to general relativistic magnetohydrodynamic (GRMHD) simulations and observations of M87*. For the parameter space studied, we characterize the transition between optically thin and thick regimes and the frequency at which the first photon ring is observable. Our results highlight the need for careful definitions of photon ring radius in the image domain, as in certain models the highly lensed photon ring is dimmer than the direct emission at certain angles. We find that at low frequencies the ring radii are set by the electron temperature, while at higher frequencies the magnetic field strength plays a more significant role, demonstrating how multifrequency analysis can also be used to infer plasma parameters. Lastly, we show how our implementation can qualitatively reproduce multifrequency black hole images from GRMHD simulations when adding time variability to our disk model through Gaussian random fields. This approach provides a new method for simulating observations from the Event Horizon Telescope and the proposed Black Hole Explorer space mission.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2007

Ring-Down Spectroscopy for Characterizing a CW Raman Laser

.A relatively simple technique for characterizing an all-resonant intracavity continuous-wave (CW) solid-state Raman laser involves the use of ring-down spectroscopy. As used here, characterizing signifies determining such parameters as threshold pump power, Raman gain, conversion efficiency, and quality factors (Q values) of the pump and Stokes cavity modes. Heretofore, in order to characterize resonant-cavity-based Raman lasers, it has usually been necessary to manipulate the frequencies and power levels of pump lasers and, in each case, to take several sets of measurements. In cases involving ultra-high-Q resonators, it also has been desirable to lock pump lasers to resonator modes to ensure the quality of measurement data. Simpler techniques could be useful. In the present ring-down spectroscopic technique, one infers the parameters of interest from the decay of the laser out of its steady state. This technique does not require changing the power or frequency of the pump laser or locking the pump laser to the resonator mode. The technique is based on a theoretical analysis of what happens when the pump laser is abruptly switched off after the Raman generation reaches the steady state. The analysis starts with differential equations for the evolution of the amplitudes of the pump and Stokes electric fields, leading to solutions for the power levels of the pump and Stokes fields as functions of time and of the aforementioned parameters. Among other things, these solutions show how the ring-down time depends, to some extent, on the electromagnetic energy accumulated in the cavity. The solutions are readily converted to relatively simple equations for the parameters as functions of quantities that can be determined from measurements of the time-dependent power levels. For example, the steady-state intracavity conversion efficiency is given by G1/G2 1 and the threshold power is given by Pin(G2/G1)2, where Pin is the steady-state input pump power immediately prior to abrupt switch-off, G1 is the initial rate of decay of the pump field, and G2 is the final rate of decay of the pump field. Hence, it is possible to determine all the parameters from a single ring-down scan, provided that the measurements taken in that scan are sufficiently accurate and complete.

Matsko, Andrey↗

NASA NTRS1994

On the rms anisotropy at 7 deg and 10 deg observed in the COBE-DMR two year sky maps

The frequency-independent rms temperature fluctuations determined from the Cosmic Background Explorer-Differential Microwave Radiometer (COBE-DMR) two-year sky maps are used to infer the parameter Q(sub rms-PS), which characterizes the normalization of power-law models of primordial cosmological temperature anisotropy, for a forced fit to a scale-invariant Harrison-Zel'dovich (n = 1) spectral model. Using a joint analysis of the 7 deg and 10 deg 'cross'-rms derived from both the 53 and 90 GHz sky maps, we find Q(sub rms-PS) = 17.0(sub -2.1 sup +2.5) micro Kelvin when the low quadrupole is included, and Q(sub rms-PS) = 19.4(sub -2.1 sup +2.3) micro Kelvin excluding the quadrupole. These results are consistent with the n = 1 fits from more sensitive methods. The effect of the low quadrupole derived from the COBE-DMR data on the inferred Q(sub rms-PS) normalization is investigated. A bias to lower Q(sub rms-PS) is found when the quadrupole is included. The higher normalization for a forced n = 1 fit is then favored by the cross-rms technique.

Banday, A. J.↗

DOE OSTI2026

Improving the stability and performance of MagLIF implosions by applying dielectric coatings and increasing applied B z , fuel preheat, and load current

We report two magnetized liner inertial fusion (MagLIF) experiments that produced record thermonuclear D–D neutron yields of 2.11×10 13 and 2.33×10 13 . These yields are about a factor of two higher than previous MagLIF results. The experiments achieved ion temperatures of 3.0 and 3.3 keV and stagnation pressures of 1.6 and 1.3 Gbar. The inferred Lawson parameters were χ=0.2 and 0.1, which are the largest reported for MagLIF. The performance increase used a high-aspect-ratio beryllium liner with a dielectric coating and modest increases in preheat energy (∼2.2 kJ), peak current (18.5 MA), and axial magnetic field (15 T). Three-dimensional HYDRA simulations are consistent with the measured liner dynamics and fusion outputs. These results indicate a pathway to higher-yield MagLIF designs using coated, high-aspect-ratio liners and improved input parameters. Simulations further suggest that adding an ice fuel layer could increase yield by up to a factor of 2.5 by reducing liner convergence, instability feedthrough, and mix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Bayesian calibration of bubble size dynamics applied to CO 2 gas fermenters

To accelerate the scale-up of gaseous CO 2 fermentation reactors, computational models need to predict gas-to-liquid mass transfer which requires capturing the bubble size dynamics, i.e. bubble breakup and coalescence. However, the applicability of existing models beyond air–water mixtures remains to be established. Here, an inverse modeling approach, accelerated with a neural network surrogate, calibrates the breakup and coalescence closure models, that are used in class methods for population balance modeling (PBM). The calibration is performed based on experimental results obtained in a CO 2 -air–water-coflowing bubble column reactor. Bayesian inference is used to account for noise in the experimental dataset and bias in the simulation results. To accurately capture gas holdup and interphase mass transfer, the results show that the breakage rate needs to be increased by one order of magnitude. In conclusion, the inferred model parameters are then used on a separate configuration and shown to also improve bubble size distribution predictions.

09 BIOMASS FUELS↗