Linear Covariance Analysis Framework for Aerospace Vehicle Trajectory Modeling and Parametric Design.
Abstract not provided.
Heflin, Lucas↗
SEARCH · Search NASA
Results for “linear covariance”
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.
Linear Covariance Analysis Framework for Aerospace Vehicle Trajectory Modeling and Parametric Design.
Abstract not provided.
Whitten, William↗
Theoretical Understanding of the Linear Relationship between Convective Updrafts and Cloud-Base Height for Shallow Cumulus Clouds. Part II: Continental Conditions
This is the Part II of a two-part study that seeks a theoretical understanding of an empirical relationship for shallow cumulus clouds: subcloud updraft velocity covaries linearly with the cloud-base height. This work focuses on continental cumulus clouds that are more strongly forced by surface fluxes and more deviated from equilibrium than those over oceans (Part I). We use a simple analytical model for shallow cumulus that is well tested against a high-resolution (25 m in the horizontal) large-eddy simulation model. Consistent with a conventional idea, we find that surface Bowen ratio is the key variable that regulates the covariability of both parameters: under the same solar insolation, a drier surface allows for stronger buoyancy flux, triggering stronger convection that deepens the subcloud layer. We find that the slope of the Bowen-ratio-regulated relationship between the two parameters (defined as λ) is dependent on both the local time and the stability of the lower free atmosphere. The value of λ decreases with time exponentially from sunrise to early afternoon and linearly from early afternoon to sunset. The value of λ is larger in a more stable atmosphere. In addition, continental λ in the early afternoon more than doubles the oceanic λ. Validation of the theoretical results against ground observations over the Southern Great Plains shows a reasonable agreement. Physical mechanisms underlying the findings are explained from the perspective of different time scales at which updrafts and cloud-base height respond to a surface flux forcing.
54 ENVIRONMENTAL SCIENCES↗
Compositional zoning of the Otowi Member of the Bandelier Tuff, Valles caldera, New Mexico, USA
The Otowi Member of the Bandelier Tuff erupted at ca. 1.60 Ma from the Valles caldera (New Mexico, USA). It consists of as much as 400 km 3 (dense rock equivalent) of strongly differentiated high-silica rhyolite and shows systematic upward variations in crystallinity, mineral chemistry, and trace element concentrations through its thickness, but the major element composition is almost constant and is near the low-pressure granite minimum. Incompatible trace elements in whole pumice fragments and glasses show well-correlated linear covariations. Upward zoning to lower abundances of incompatible trace elements is accompanied by development of overgrowths on quartz and alkali feldspar, although earlier-formed interiors of quartz and feldspar have near-constant compositions throughout the tuff, modified by cation diffusion in the case of feldspar. Melt inclusions in remnant quartz cores show diverse Pb isotope ratios, pointing to a wide range of distinct protoliths that contributed rhyolitic melt to the Otowi magma. Mineral thermometers suggest a modest temperature gradient through the melt body, perhaps of 40 °C, at the time of eruption. Chemical, textural, and mineralogical variations and volume-composition relations through the tuff are consistent with an origin for zoning by melting of a high-crystallinity cumulate layer beneath cognate supernatant liquid to produce denser, remobilized liquid of accumulative composition (i.e., the “modified mush model”). Melting may have occurred in several episodes. The latest of these episodes, probably thousands of years prior to eruption, introduced new rhyolitic liquid into the system and was associated with a thermal excursion, recorded in core compositions of pyroxene, during which much of the earlier crystal mass was dissolved. This left inherited cores and interiors of accumulated quartz and feldspar mantled with new growth having less-evolved compositions (higher Ti, Sr, and Ba). Changing solubility of zircon during cumulate melting produced a reversal of Zr concentrations. There is no clear petrologic evidence of a recharge eruption trigger; nonetheless, compositional zoning resulted mainly from repeated recharge-induced remobilization of quartz-feldspar cumulate. The Otowi system was built, evolved, and modified by several events over the course of a few hundred thousand years.
Wolff, J. A.↗
A data science approach for analysis and reconstruction of spinodal-like composition fields in irradiated FeCrAl alloys
A statistical method for the analysis of continuously distributed data representative of composition fluctuations in irradiated FeCrAl alloys acquired using Energy Dispersive X-ray Spectroscopy (EDS) method is presented. Using probability distribution functions, direct and cross-covariances between the elemental compositions, the effects of alloy composition and irradiation dose were investigated on the spatial distribution and length scale of composition fluctuations at the nanoscale. We have observed that, for neutron-irradiated FeCrAl alloys, the distribution of Fe and Cr followed a left-skewed and right-skewed distribution, respectively for all (average) alloy compositions and irradiation doses. The analysis also revealed enhanced spatial gradients in the elemental compositions at higher irradiation dose. Direct and cross-covariance estimates of the experimental data were also utilized for reconstruction of composition data through fitting it to a parametric form of the covariance functions. Linear Model of Coregionalization was used to determine the parameters of the covariance functions. Subsequently, a spectral method was utilized for simulating a realization of the alloy compositions. Close correspondence was observed between the experimental and the reconstructed data which was analyzed using probability distribution functions and covariance functions. Composition space of the experimental and reconstructed data and dislocation velocities as a function of applied stress and line directions over the entire composition maps were also examined.
36 MATERIALS SCIENCE↗
Size-Dependent Characteristics of Surface-Rooted Three-Dimensional Convective Objects in Continental Shallow Cumulus Simulations
A segmentation algorithm is applied to high resolution simulations of shallow continental convection to identify individual convective 3D objects within the convective boundary layer, in order to investigate which properties of the objects vary with the object width. The study analyses the geometry of the objects, along with their profiles of vertical velocity and total water, to assess various assumptions often used in spectral mass-flux convection schemes. The methodology of this paper is unique in that we use (a) a novel application of the watershed algorithm to detect individual objects in the constantly evolving continental boundary layer efficiently, and (b) an unprecedentedly large number of simulations being analyzed. In total, 26 days of LASSO simulations at the Atmospheric Radiation Measurement Southern Great Plains site are analyzed, yielding roughly one million objects. Plume-like surface-rooted objects are found to be omnipresent, the vertical extent of which is strongly dependent on the object width. The vertical velocity and moisture anomaly profiles of the widest objects are roughly consistent with the classic buoyancy-driven rising plume model. The kinematic and thermodynamic properties of the objects vary with object width. This width dependence is strongest above cloud base, but much weaker below. Finally the impact of neglecting the contribution of covariances to the vertical moisture flux, which is commonly used in mass-flux parameterizations, is investigated. The average effect of neglecting covariances increases linearly with object width, leading to a 20% flux underestimation for 2 km wide objects. Implications of the results for spectral convection scheme development are briefly discussed.
54 ENVIRONMENTAL SCIENCES↗
Eigen-decomposition of covariance matrices: An application to the BAO linear point
Not Available
Lee, Jaemyoung (Jason) (ORCID:0000000166339793)↗
Allosteric prediction via convolutional neural networks and protein structural and dynamical features
Allostery is the phenomenon whereby a binding event or covalent modification at one site in a protein modulates function at a distal site, thus changing a protein’s functional state. As such, it is a ubiquitous aspect of protein functional regulation. Computationally predicting allosteric states is important as part of the broader challenge of functional annotation, but it also has practical implications for drug development, as targeting an allosteric site often affords greater specificity compared with targeting an orthosteric site. This study introduces a machine learning approach to predict the allosteric functional state using the small G-protein KRas as the model system, due to its implication in many types of cancer and being well studied as a result with many x-ray crystallographic structures of KRas available with different mutations and ligands bound. Using structural and dynamical features that can be cast as images, namely interatomic distances, contact maps, covariance, and mutual information, supervised learning was performed using convolutional neural networks. Two pretrained convolutional neural network architectures, GoogLeNet and ResNet18, were fine-tuned to classify KRas into active or inactive states based on these features. Across training regimes, atomic contact maps emerged as the most effective structural feature, whereas linearized mutual information outperformed covariance in capturing dynamical correlations relevant to allostery. Models achieved significant validation accuracy, with atomic contact maps yielding up to 90% accuracy. In conclusion, the findings suggest that integrating global structural rearrangements and correlated motion patterns with deep learning can reliably predict protein allosteric states, offering a promising framework for understanding allosteric regulation and developing targeted therapeutics.
Rajeshwar T., Rajitha [Oak Ridge National Laborato↗
Large-Scale Inference of Multivariate Regression for Heavy-Tailed and Asymmetric Data
Large-scale multivariate regression is a fundamental statistical tool with a wide range of applications. Here, this study considers the problem of simultaneously testing a large number of general linear hypotheses, encompassing covariate-effect analysis, analysis of variance, and model comparisons. The challenge that accompanies a large number of tests is the ubiquitous presence of heavy-tailed and/or highly skewed measurement noise, which is the main reason for the failure of conventional least squares-based methods. For large-scale multivariate regression, we develop a set of robust inference methods to explore data features such as heavy tailedness and skewness, which are not visible to least squares methods. The new testing procedure is based on the data-adaptive Huber regression and a new covariance estimator of regression estimates. Under mild conditions, we show that our methods produce consistent estimates of the false discovery proportion. Extensive numerical experiments and an empirical study on quantitative linguistics demonstrate the advantage of the proposed method over many state-of-the-art methods when the data are generated from heavy-tailed and/or skewed distributions.
97 MATHEMATICS AND COMPUTING↗
Nonstationary multivariate Gaussian processes for electronic health records
Not Available
97 MATHEMATICS AND COMPUTING↗
Regression of exchangeable relational arrays
Relational arrays represent measures of association between pairs of actors, often in varied contexts or over time. Trade flows between countries, financial transactions between individuals, contact frequencies between school children in classrooms and dynamic protein-protein interactions are all examples of relational arrays. Elements of a relational array are often modelled as a linear function of observable covariates. Uncertainty estimates for regression coefficient estimators, and ideally the coefficient estimators themselves, must account for dependence between elements of the array, e.g., relations involving the same actor. Existing estimators of standard errors that recognize such relational dependence rely on estimating extremely complex, heterogeneous structure across actors. This paper develops a new class of parsimonious coefficient and standard error estimators for regressions of relational arrays. Here we leverage an exchangeability assumption to derive standard error estimators that pool information across actors, and are substantially more accurate than existing estimators in a variety of settings. This exchangeability assumption is pervasive in network and array models in the statistics literature, but not previously considered when adjusting for dependence in a regression setting with relational data. We demonstrate improvements in inference theoretically, via a simulation study, and by analysis of a dataset involving international trade.
97 MATHEMATICS AND COMPUTING↗
Modeling Systematic Discrepancies in Nuclear-Data Measurements with Machine Learning
Nuclear data and their associated co-variances are constantly being reevaluated as techniques improve and as new experimental data, as well as nuclearmodel developments, emerge. A standard technique used to evaluate mean values in nuclear data, and their associated covariances, is the generalized linear least squares (GLLS) method. Aligning with recent efforts to incorporate measurement features into nuclear data evaluation, we augment GLLS by including a linear term which attempts to predict potential systematic discrepancies in experimental data as related to the measurement features. Due to the general nature of this augmentation, we are able to apply this evaluation to three key observables: neutroninduced fission cross sections, the average prompt neutron multiplicity, and the prompt-fission neutron spectrum of ²³⁹Pu.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
Bayesian Monte Carlo Evaluation Framework for Cross Sections Nuclear Data and Integral Benchmark Experiments
The new Bayesian Monte Carlo (MC) evaluation framework described in this abstract has been conceived as an attempt to improve nuclear data evaluations of differential crosssection data by removing the following two approximations conventionally employed for nuclear data evaluations: all probability density functions (PDFs) of all data and model parameters, both prior and posterior, are assumed to be normal (i.e., Gaussian) PDFs, and all uncertainties and covariances are propagated using a linear approximation. With these approximations removed, the Bayesian MC (BMC) framework could be used to account for nonlinear effects and would enable improved evaluations of differential cross sections and IBE data that are presently performed based on the assumptions itemized above. The BMC would also improve upon the uniform sampling of IBE parameters from within ranges defined by their evaluated uncertainties.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Chiral properties of the nucleon interpolating current and θ-dependent observables
We revisit the chiral properties of nucleon interpolating currents, and show that of the two leading order currents j 1 and j 2 , only two linear combinations j 1 ± j 2 transform covariantly under the anomalous U(1) A symmetry. As a result, calculations of quantities which vanish by symmetry in the chiral limit may produce unphysical results if carried out with different linear combinations of the currents. This includes observables such as electric dipole moments, induced by the quantum chromodynamics (QCD) parameter θ, and the θ-dependence of the nucleon mass. For completeness, we also exhibit the leading order results for nucleon electric dipole moments (d n,p ) induced by θ, and the nucleon magnetic moments (μ n,p ), when calculated using QCD sum rules for both the covariant choices of the nucleon interpolating current. The results in each channel, conveniently expressed as the ratios, d n,p /μ n,p , are numerically consistent, and reflect the required physical dependence on θ.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Assessing variation in faecal glucocorticoid concentrations in gray whales exposed to anthropogenic stressors
Understanding how individual animals respond to stressors behaviourally and physiologically is a critical step towards quantifying long-term population consequences and informing management efforts. Glucocorticoid (GC) metabolite accumulation in various matrices provides an integrated measure of adrenal activation in baleen whales and could thus be used to investigate physiological changes following exposure to stressors. In this study, we measured GC concentrations in faecal samples of Pacific Coast Feeding Group (PCFG) gray whales (Eschrichtius robustus) collected over seven consecutive years to assess the association between GC content and metrics of exposure to sound levels and vessel traffic at different temporal scales, while controlling for contextual variables such as sex, reproductive status, age, body condition, year, time of year and location. We develop a Bayesian Generalized Additive Modelling approach that accommodates the many complexities of these data, including non-linear variation in hormone concentrations, missing covariate values, repeated samples, sampling variability and some hormone concentrations below the limit of detection. Estimated relationships showed large variability, but emerging patterns indicate a strong context-dependency of physiological variation, depending on sex, body condition and proximity to a port. Our results highlight the need to control for baseline hormone variation related to context, which otherwise can obscure the functional relationship between faecal GCs and stressor exposure. Therefore, extensive data collection to determine sources of baseline variation in well-studied populations, such as PCFG gray whales, could shed light on cetacean stress physiology and be used to extend applicability to less-well-studied taxa. GC analyses may offer greatest utility when employed as part of a suite of markers that, in aggregate, provide a multivariate measure of physiological status, better informing estimates of individuals’ health and ultimately the consequences of anthropogenic stressors on populations.
59 BASIC BIOLOGICAL SCIENCES↗
Taylor approximation variance reduction for approximation errors in PDE-constrained Bayesian inverse problems
In numerous applications, surrogate models are used as a replacement for accurate parameter-to-observable mappings when solving large-scale inverse problems governed by partial differential equations (PDEs). The surrogate model may be a computationally cheaper alternative to the accurate parameter-to-observable mappings and/or may ignore additional unknowns or sources of uncertainty. The Bayesian approximation error (BAE) approach provides a means to account for the induced uncertainties and approximation errors, i.e. the errors between the accurate parameter-to-observable mapping and the surrogate. The statistics of these errors are, however, in general unknown a priori, and are thus calculated using Monte Carlo sampling. Although the sampling is typically carried out offline, i.e. before considering the data, the process can still represent a computational bottleneck. In this work, we develop a scalable computational approach for reducing the costs associated with the sampling stage of the BAE approach. Specifically, we consider the Taylor expansion of the accurate and surrogate forward models with respect to the uncertain parameter fields either as a control variate for variance reduction or as a means to directly and efficiently approximate the mean and covariance of the approximation errors. We propose efficient methods for evaluating the expressions for the mean and covariance of the Taylor approximations based on linear(-ized) PDE solves. Furthermore, the proposed approach is independent of the dimension of the uncertain parameter, depending instead on the intrinsic dimension of the data, ensuring scalability to high-dimensional problems. The potential benefits of the proposed approach are demonstrated for two high-dimensional inverse problems governed by PDE examples, namely for the estimation of a distributed Robin boundary coefficient in a linear diffusion problem, and for a coefficient estimation problem governed by a nonlinear diffusion problem.
Bayesian approximation error↗
Bayesian Monte Carlo Evaluation Framework for Imperfect Data and Models [Abstract]
Nuclear data evaluation methods conventionally make the following assumptions: prior and posterior probability distribution functions (PDFs) of all model parameters and data are normal (Gaussian); the linear approximation is sufficiently accurate for minimization of a cost function (even for non-linear models); and that both the model (of, e.g., neutron cross section) and experimental data (including their covariance data) are perfect. These assumptions are inherent to the well-known generalized linear least squares (GLLS) minimization method commonly used for evaluations of resolved resonance region (RRR) neutron cross sections. However, these assumptions are often not justified due to the presence of non-normal PDFs, non-linear models (e.g. R -matrix formalism), and inherent imperfections in data and models (e.g. discrepant data sets, discrepancies between the previous evaluation and newly measured data, or imperfect covariance data). We remove the said assumptions in a mathematical framework of Bayes’ theorem, and implement it using the Metropolis-Hastings Monte Carlo method. Parameters of a new kind are introduced to parameterize inherent imperfections, e.g. , any discrepancies between the theoretical model and measured data. These new parameters enable evaluators to quantify their expert judgement about any discrepancies or imperfections in a reproducible manner. We demonstrate the framework with an ongoing evaluation of 233 U in the eV region using the ENDF-B/VIII library and transmission data measured by Guber, et al. , and compare the posterior parameters to those obtained by conventional evaluation methods.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Leveraging design of experiments to build chemometric models for the quantification of uranium (VI) and HNO3 by Raman spectroscopy
Partial least squares regression (PLSR) and support vector regression (SVR) models were optimized for the quantification of U(VI) (10–320 g L −1 ) and HNO 3 (0.6–6 M) by Raman spectroscopy with optimized calibration sets chosen by optimal design of experiments. The designed approach effectively minimized the number of samples in the calibration set for PLSR and SVR by selecting sample concentrations with a quadratic process model, despite complex confounding and covarying spectral features in the spectra. The top PLS2 model resulted in percent root mean square errors of prediction for U(VI), HNO 3 , and NO 3 − of 3.7%, 3.6%, and 2.9%, respectively. PLS1 models performed similarly despite modeling an analyte with a majority linear response (i.e., uranyl symmetric stretch) and another with more covarying vibrational modes (i.e., HNO 3 ). Partial least squares (PLS) model loadings and regression coefficients were evaluated to better understand the relationship between weaker Raman bands and covarying spectral features. Support vector machine models outperformed PLS1 models, resulting in percent root mean square error of prediction values for U(VI) and HNO 3 of 1.5% and 3.1%, respectively. The optimal nonlinear SVR model was trained using a similar number of samples (11) compared with the PLSR model, even though PLS is a linear modeling approach. The generic D-optimal design presented in this work provides a robust statistical framework for selecting training set samples in disparate two-factor systems. This approach reinforces Raman spectroscopy for the quantification of species relevant to the nuclear fuel cycle and provides a robust chemometric modeling approach to bolster online monitoring in challenging process environments.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗