SEARCH · Search NASA
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.
An analytic model approach is applied to several simple atmospheric inversion problems. This method gives a sharp determination of aerosol size distribution parameters. It is shown that this analytic approach, together with ground level point sampling data measurements, can be used to infer information on the tropospheric ozone profile.
Observations from magnetic field and plasma experiments on ISEE 3 are used to identify and to analyze 16 crossings of the distant magnetotail neutral sheet. Four methods of analyses have been employed to measure the extent of the magnetotail twisting, which is found to be much greater than what was previously observed or, the theoretically inferred twist. A twist of more than 90 deg has frequently been observed on both open and closed magnetic field lines. In addition, the location of the magnetic neutral line of the magnetotail has been found to move over the distance interval of 115 to 220 Earth radii. Some features of these results suggest that a modification of the present solar wind-geomagnetic field interaction model may be needed.
Various killing mechanisms have been suggested to contribute to the mass extinctions at the KT boundary, including severe, global deterioration of the atmosphere and hydrosphere due to SO(x) released from heavily shocked, sulfate-bearing target rocks. The devolatilization of anhydrite is predominantly inferred from thermodynamic considerations and lacks experimental confirmation. To date, the experimentally determined shock behavior of anhydrite is limited to solid-state effects employing X-ray diffraction methods. The present report employs additional methods to characterize experimentally shocked anhydrite.
NASA has been collecting massive amounts of remote sensing data about Earth's systems for more than a decade. Missions are selected to be complementary in quantities measured, retrieval techniques, and sampling characteristics, so these datasets are highly synergistic. To fully exploit this, a rigorous methodology for combining data with heterogeneous sampling characteristics is required. For scientific purposes, the methodology must also provide quantitative measures of uncertainty that propagate input-data uncertainty appropriately. We view this as a statistical inference problem. The true but notdirectly- observed quantities form a vector-valued field continuous in space and time. Our goal is to infer those true values or some function of them, and provide to uncertainty quantification for those inferences. We use a spatiotemporal statistical model that relates the unobserved quantities of interest at point-level to the spatially aggregated, observed data. We describe and illustrate our method using CO2 data from two NASA data sets.
We present a detection of the unnormalized skewness (T(sup ~)(sup 2)(n(circumflex)) induced by the thermal Sunyaev-Zel'dovich (tSZ) effect in filtered Atacama Cosmology Telescope (ACT) 148 GHz cosmic microwave background temperature maps. Contamination due to infrared and radio sources is minimized by template subtraction of resolved sources and by constructing a mask using outlying values in the 218 GHz (tSZ-null) ACT maps. We measure (T(sup ~)(sup 3) (n(circumflex)) = -31 plus or minus 6 micro-K(sup 3) (measurement error only) or plus or minus 14 micro-K(sup 3) (including cosmic variance error) in the filtered ACT data, a 5sigma detection. We show that the skewness is a sensitive probe of sigma(sub 8), and use analytic calculations and tSZ simulations to obtain cosmological constraints from this measurement. From this signal alone we infer a value of sigma(sub 8) = 0.78 sup +0.03 sub -0.04 (68% C.L.) sup +0.05 sub -0.16. Our results demonstrate that measurements of nonGaussianity can be a useful method for characterizing the tSZ effect and extracting the underlying cosmological information.
Stratospheric aerosols play key roles in the chemistry and radiation balance of the atmosphere and are a key input parameter for global chemistry and climate models. The degree to which aerosols impact chemistry and radiation balance depends primarily on the relative abundance of different sized particles within the sample volume, often referred to as the particle size distribution (PSD). If the PSD is accurately known then other key modeling parameters (e.g., surface area density and effective radius) can be derived. Historically, occultation observations from orbital instruments such as SAGE III/ISS have been used to infer these PSD parameters by inverting the extinction coefficient spectra. However, past efforts routinely failed to account for measurement uncertainty and lacked a rigorous estimate of the inferred PSD uncertainty. We carried out a series of simulations to evaluate the accuracy of these inferences and, for every valid SAGE III/ISS extinction spectrum, determined the range of PSD parameters that fell withing the bounds of the extinction error bars. Special application of this method was applied to estimate the impact of the 2022 Hunga Tonga eruption had on particle size distributions.
Noise bias is a significant source of systematic error in weak gravitational lensing measurements that must be corrected to satisfy the stringent standards of modern imaging surveys in the era of precision cosmology. This paper reviews the analytical noise bias correction method and provides analytical derivations demonstrating that we can recover shear to its second order using the ‘renoising’ noise bias correction approach introduced by METACALIBRATION. We implement this analytical noise bias correction within the AnaCal shear estimation framework and propose several enhancements to the noise bias correction algorithm. We evaluate the improved AnaCal using simulations designed to replicate Rubin Legacy Survey of Space and Time (LSST) imaging data. These simulations feature semi-realistic galaxies and stars, complete with representative distributions of magnitudes and Galactic spatial density. We conduct tests under various observational challenges, including cosmic rays, defective CCD columns, bright star saturation, bleed trails, and spatially variable point spread functions. Our results indicate a multiplicative bias in weak lensing shear recovery of less than a few tenths of a per cent, meeting LSST Dark Energy Science Collaboration requirements without requiring calibration from external image simulations. Additionally, our algorithm achieves rapid processing, handling one galaxy in less than a millisecond.
The majority of research on efficient and scalable algorithms in computational science and engineering has focused on the forward problem : given parameter inputs, solve the governing equations to determine output quantities of interest. In contrast, in this paper, we consider the broader question: given a (large-scale) model containing uncertain parameters, (possibly) noisy observational data, and a prediction quantity of interest, how do we construct efficient and scalable algorithms to (1) infer the model parameters from the data (the deterministic inverse problem ), (2) quantify the uncertainty in the inferred parameters (the Bayesian inference problem), and (3) propagate the resulting uncertain parameters through the model to issue predictions with quantified uncertainties (the forward uncertainty propagation problem )? We present efficient and scalable algorithms for this end-to-end, data-to-prediction process under the Gaussian approximation and in the context of modeling the flow of the Antarctic ice sheet and its effect on loss of grounded ice to the ocean. The ice is modeled as a viscous, incompressible, creeping, shear-thinning fluid. The observational data come from satellite measurements of surface ice flow velocity, and the uncertain parameter field to be inferred is the basal sliding parameter, represented by a heterogeneous coefficient in a Robin boundary condition at the base of the ice sheet. The prediction quantity of interest is the present-day ice mass flux from the Antarctic continent to the ocean. We show that the work required for executing this data-to-prediction process—measured in number of forward (and adjoint) ice sheet model solves—is independent of the state dimension, parameter dimension, data dimension, and the number of processor cores. The key to achieving this dimension independence is to exploit the fact that, despite their large size, the observational data typically provide only sparse information on model parameters. This property can be exploited to construct a low rank approximation of the linearized parameter-to-observable map via randomized SVD methods and adjoint-based actions of Hessians of the data misfit functional.
Background: One-nucleon knockout reactions provide insightful information on the single-particle structure of nuclei. When applied to one-neutron halo nuclei, they are purely peripheral, suggesting that they could be properly modeled by describing the projectile within a halo effective field theory (halo-EFT). Purpose: We reanalyze the one-neutron knockout measurements of 11 Be and 15 C —both one-neutron halo nuclei—on beryllium at about 60 MeV/nucleon. We consider halo-EFT descriptions of these nuclei which already provide excellent agreement with breakup and transfer data. Method: Here we include a halo-EFT description of the projectile within an eikonal-based model of the reaction and compare its outcome to existing data. Results: Excellent agreement with experiment is found for both nuclei. The asymptotic normalization coefficients inferred from this comparison confirm predictions from ab initio nuclear-structure calculations and values deduced from transfer data. Conclusions: Halo-EFT can be reliably used to analyze one-neutron knockout reactions measured for halo nuclei and test predictions from state-of-the-art nuclear structure models on these experimental data.
Abstract With machine learning applications now spanning a variety of computational tasks, multi-user shared computing facilities are devoting a rapidly increasing proportion of their resources to such algorithms. Graph neural networks (GNNs), for example, have provided astounding improvements in extracting complex signatures from data and are now widely used in a variety of applications, such as particle jet classification in high energy physics (HEP). However, GNNs also come with an enormous computational penalty that requires the use of GPUs to maintain reasonable throughput. At shared computing facilities, such as those used by physicists at Fermi National Accelerator Laboratory (Fermilab), methodical resource allocation and high throughput at the many-user scale are key to ensuring that resources are being used as efficiently as possible. These facilities, however, primarily provide CPU-only nodes, which proves detrimental to time-to-insight and computational throughput for workflows that include machine learning inference. In this work, we describe how a shared computing facility can use the NVIDIA Triton Inference Server to optimize its resource allocation and computing structure, recovering high throughput while scaling out to multiple users by massively parallelizing their machine learning inference. To demonstrate the effectiveness of this system in a realistic multi-user environment, we use the Fermilab Elastic Analysis Facility augmented with the Triton Inference Server to provide scalable and high-throughput access to a HEP-specific GNN and report on the outcome.
Unveiling the relationship between phylogeny and function of the microbiome is crucial to determine its contribution to ecosystem functioning. However, while there is a considerable amount of information on microbial phylogenetic diversity, our understanding of its relationship to functional diversity is still scarce. Here we predicted the total microbiome functions of bacteria and fungi on Earth using the total known functions from level 3 of KEGG Orthology by modelling the increase of functions with increasing diversity of bacteria or fungi. For bacteria and fungi, the unsaturated model described the data significantly better (for both P?<2.2e-16), suggesting the presence of two types of functions. Widespread functions ubiquitous in every living organism that make up two thirds of our current knowledge of microbiome functions are separated from rare functions from specialised enzymes present in only a few species. Given previous estimates on species richness, we predicted a global total of 35.5 million functions in bacteria and 3.2 million in fungi; of which only 0.02% and 0.14% are known today. Our approach highlights the necessity of novel and more sophisticated methods to unveil the entirety of rare functions to fully understand the involvement of the microbiome in ecosystem functioning. SIGNIFICANCE: The functionality of and within a microbial community is generally inferred based on the taxonomic annotation of the organism. However, our understanding of functional diversity and how it relates to taxonomy is still limited. Here we predict the total microbiome functionality in bacteria and fungi on Earth using known and annotated protein-coding sequences in species accumulation curves. Our estimates reveal that the majority of functionality (>99%) could be assigned to yet unknown and rare functions, highlighting that our current knowledge is incomplete and functional inference is thus lackluster.
This article investigates the applicability of a recently proposed, nonlinear sparse Bayesian learning (NSBL) algorithm to identify and estimate the complex aerodynamics of limit cycle oscillations. NSBL provides a semi-analytical framework for determining the data-optimal sparse model nested within a (potentially) over-parameterized model. This is particularly relevant to nonlinear dynamical systems where modelling approaches involve the use of physics-based and data-driven components. In such cases, the data-driven components, where analytical descriptions of the physical processes are not readily available, are often prone to overfitting, meaning that the empirical aspects of these models will often involve the calibration of an unnecessarily large number of parameters. While an overparameterized model may fit the observed data well, such models may be inadequate for making predictions in regimes that are different from those wherein the data were recorded. In view of this, it is desirable to not only calibrate the model parameters, but also identify the optimal compromise between data fit and model complexity. In this article, we exhibit the optimal model discovery for an aeroelastic system wherein the structural dynamics are well-known and described by a differential equation model, coupled with a semi-empirical aerodynamic model for laminar separation flutter, resulting in low-amplitude limit cycle oscillations (LCO). To illustrate the performance of the algorithm, in this article, we use synthetic data and demonstrate the ability of the algorithm to correctly rediscover the optimal model and model parameters, given a known data-generating model. The synthetic data are generated from a forward simulation of a known differential equation model with parameters selected so as to mimic the dynamics observed in wind-tunnel experiments. Subsequently, we demonstrate the performance of the algorithm for model selection using noisy LCO data from wind tunnel experiments. As there is no ground truth available for the experimental data case, we provide a comparison between NSBL and Bayesian model selection to validate the results, and demonstrate the use of NSBL as an efficient alternative to traditional methods.
Digital elevation maps obtained using TanDEM-X and Pleiades data combined with newly obtained surface age estimates using cosmogenic radionuclide (CRN) and optically simulated luminescence (OSL) methods are used to quantify the slip rate along the western section of the Altyn Tagh Fault in southern Xinjiang. The reconstruction of the conical shape of massive alluvial fans inferred to be from the Eemian (115±7 ka) from CRN dating shows consistent left-lateral offsets of 300±20 m, yielding a slip rate of 2.6±0.3 mm/year. Successive episodes of incision have left cut terraces inset in wide canyons, 10–25 m below the fans' surface. The incision was followed by the deposition of a broad terrace of early Holocene age, which is reincised by modern stream channels. Near the village of Shanxili, a 200 m wide valley is partially dammed by a shutter ridge displaced by the fault. A fill terrace deposited upstream from the ridge has an OSL age of 8.8±0.6 ka. The 23±2 m offset of the riser incising the terrace indicates a minimum postdepositional movement on the fault, yielding a Holocene rate of 2.6±0.5 mm/year, consistent with the 115-ka average slip rate. Furthermore, scarp degradation analysis using mass diffusion reveals a nonlinear relationship between fault displacement and degradation coefficient along the progressively exposed fault scarp, a pattern suggesting either seismic clustering or variable diffusion rate since the Eemian. Together with the Gozha Co-Longmu Co fault to the south, the Karakax section of the Altyn Tagh Fault contributes to the eastward movement of the western corner of Tibet.
A unified, Bayesian inference of midplane electron temperature and density profiles using both Thomson scattering (TS) and interferometric data is presented. Beyond the Bayesian nature of the analysis, novel features of the inference are the use of a Gaussian process prior to infer a mollification length-scale of inferred profiles and the use of Gauss-Laguerre quadratures to directly calculate the depolarisation term associated with the TS forward model. Results are presented from an application of the method to data from the high resolution TS system on the Mega-Ampere Spherical Tokamak, along with a comparison to profiles coming from the standard analysis carried out on that system.
We discuss methods to analyze a quantum Gibbs family in the ultra-cold regime where the norm closure of the Gibbs family fails due to discontinuities of the maximum-entropy inference. The current discussion of maximum-entropy inference and irreducible correlation in the area of quantum phase transitions is a major motivation for this research. We extend a representation of the irreducible correlation from finite temperatures to absolute zero.
Scanning near-field optical microscopy is one of the most effective techniques for spectroscopy of nanoscale systems. However, inferring optical constants from the measured near-field signal can be challenging because of a complicated and highly nonlinear interaction between the scanned probe and the sample. Conventional fitting methods applied to this problem often suffer from the lack of convergence or require human intervention. Here, we develop an alternative approach where the optical parameter extraction is automated by a deep learning network. The network provides an initial estimate that is subsequently refined by a traditional fitting algorithm. In conclusion, we show that this method demonstrates superior accuracy, stability against noise, and computational speed when applied to simulated near-field spectra.
Differential cross section measurements are the currency of scientific exchange in particle and nuclear physics. A key challenge for these analyses is the correction for detector distortions, known as deconvolution or unfolding. Binned unfolding of cross section measurements traditionally rely on the regularized inversion of the response matrix that represents the detector response, mapping pre-detector (`particle level') observables to post-detector (`detector level') observables. In this paper we introduce Neural Posterior Unfolding, a modern, Bayesian approach that leverages normalizing flows for unfolding. By using normalizing flows for neural posterior estimation, NPU offers several key advantages including implicit regularization through the neural network architecture, fast amortized inference that eliminates the need for repeated retraining, and direct access to the full uncertainty in the unfolded result. In addition to introducing NPU, we implement a classical Bayesian unfolding method called Fully Bayesian Unfolding (FBU) in modern Python so it can also be studied. These tools are validated on simple Gaussian examples and then tested on simulated jet substructure examples from the Large Hadron Collider (LHC). We find that the Bayesian methods are effective and worth additional development to be analysis ready for cross section measurements at the LHC and beyond.
A solar active region (AR) is a three-dimensional (3D) magnetic structure formed in the convection zone, whose property is fundamentally important for determining the coronal structure and solar activity when emerged. However, our knowledge of the detailed 3D structure prior to its emergence is rather poor, largely limited by the low cadence and sensitivity of previous instruments. Here, using the 45 s high-cadence observations from the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory, we are able for the first time to reconstruct a 3D data cube and infer the detailed subsurface magnetic structure of NOAA AR 11158, and to characterize its magnetic connectivity and topology. This task is accomplished with the aid of the image-stacking method and advanced 3D visualization. We find that the AR consists of two major bipoles or four major polarities. Each polarity in 3D shows interesting tree-like structure, i.e., while the root of the polarity appears as a single tree-trunk-like tube, the top of the polarity has multiple branches consisting of smaller and thinner flux tubes which connect to the branches of the opposite polarity that is similarly fragmented. The roots of the four polarities align well along a straight line, while the top branches are slightly non-coplanar. Our observations suggest that an active region, even appearing highly complicated on the surface, may originate from a simple straight flux tube that undergoes both horizontal and vertical bifurcation processes during its rise through the convection zone.