Search NASA⌕ Search

SEARCH · Search NASA

Results for “Sensitivity index function”

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 19 records

Stochastic Framework for Optimal Control of Planetary Reentry Trajectories Under Multilevel Uncertainties

We present a novel stochastic optimal control framework that accounts for various types of uncertainties, with application to reentry trajectory planning. The formulation of the optimal trajectory control problem is presented in the context of an indirect method where a functional objective associated with the terminal vehicle speed is to be minimized. Uncertain input parameters in the optimal trajectory control model, including aerodynamic parameters and initial and terminal conditions, are modeled as aleatory random variables, while the statistical parameters of these aleatory distributions are themselves random variables. The parametric and model uncertainties are simultaneously propagated through an extended polynomial chaos expansion (EPCE) formalism. Several metrics are described to evaluate response statistics and presented as insightful tools for robust decision making. Specifically, the response probability density function (PDF) reflecting influence of both epistemic and aleatory uncertainties is obtained. By sampling over the random variables representing model error, an ensemble of response PDFs is generated and the associated failure probability is estimated as a random variable with its own polynomial chaos expansion. Besides, the sensitivity index functions of response PDF with respect to the statistical parameters are evaluated. Coupling parametric and model uncertainties within the EPCE framework leads to a robust and efficient paradigm for multilevel uncertainty propagation and PDF characterization in general optimal control problems.

Engineering↗

A functional global sensitivity measure and efficient reliability sensitivity analysis with respect to statistical parameters

Sensitivity analysis and reliability assessment are two important aspects of structural and system safety. Epistemic uncertainty with respect to probabilistic model of input parameters due to lack of knowledge is present in many scarce-data applications and complicates the characterization of uncertainty in model response. In this article, we present two importance measures to evaluate the impact of distribution parameters on the probability distribution function (PDF) of the output and the failure probability. The epistemic uncertainty associated with the distribution parameters is modeled as random variables. Additionally, a modified extended polynomial chaos expansion (MEPCE) approach is introduced in which aleatory and epistemic random variables are modeled and propagated simultaneously while allowing the separate assessment for any single epistemic variable. A MEPCE-based kernel density estimation (KDE) construction provides a composite map from each epistemic variable to the response PDF. The functional global sensitivity index of the PDF with respect to the distribution parameters is thus derived, as a function of output, which is both more informative and more efficient than standard scalar sensitivity measures. Reliability sensitivity indices can be readily evaluated by integrating the global sensitivity index function over the failure zone. Three illustrative examples are used to demonstrate the proposed methodology.

42 ENGINEERING↗

Process Interactions Can Change Process Ranking in a Coupled Complex System Under Process Model and Parametric Uncertainty

For a complex hydrologic system with multiple processes and process interactions, global sensitivity analysis is often used to identify important or influential parameters for model development and improvement. The identification is complicated by process model uncertainty, when a system process can be represented by multiple process models. This study develops a new total-effect process sensitivity index to identify influential processes under model uncertainty. This is done by extending Sobol's total-effect parameter sensitivity index for one system model to total-effect process sensitivity index for multiple system models to account for uncertainty in process models and model parameters. The total-effect process sensitivity index includes not only the first-order process sensitivity index for measuring the importance of individual processes but also higher-order indices that account for process interactions. The total-effect process sensitivity index can identify an influential process that itself and its interactions with other processes influence a model output. Here, the total-effect process sensitivity index is applied to two numerical examples: (a) Sobol's G*-functions with analytical solutions of first-order and total-effect process sensitivity indices, and (b) groundwater flow models with interactions between recharge, geology, and snowmelt processes. The second evaluation shows that, due to second-order and higher-order process interactions, the first-order and total-effect process sensitivity indices give different process ranking. It is thus necessary to estimate both first-order and total-effect process sensitivity indices to appreciate the difference between the first-order impact of a process alone and the overall total-effect impact of the process itself and its interactions with other processes on a model output.

54 ENVIRONMENTAL SCIENCES↗

Probing plasma physics with spectral index maps of accreting black holes on event horizon scales

The Event Horizon Telescope (EHT) collaboration has produced the first resolved images of the supermassive black holes at the centre of our galaxy and at the centre of the elliptical galaxy M87. As both technology and analysis pipelines improve, it will soon become possible to produce spectral index maps of black hole accretion flows on event horizon scales. Here, in this work, we predict spectral index maps of both M87* and Sgr A* by applying the general relativistic radiative transfer (GRRT) code IPOLE to a suite of general relativistic magnetohydrodynamic (GRMHD) simulations. We analytically show that the spectral index increases with increasing magnetic field strength, electron temperature, and optical depth. Consequently, spectral index maps grow more negative with increasing radius in almost all models, since all of these quantities tend to be maximized near the event horizon. Additionally, photon ring geodesics exhibit more positive spectral indices, since they sample the innermost regions of the accretion flow with the most extreme plasma conditions. Spectral index maps are sensitive to highly uncertain plasma heating prescriptions (the electron temperature and distribution function). However, if our understanding of these aspects of plasma physics can be tightened, even the spatially unresolved spectral index around 230 GHz can be used to discriminate between models. In particular, Standard and Normal Evolution (SANE) flows tend to exhibit more negative spectral indices than Magnetically Arrested Disc (MAD) flows due to differences in the characteristic magnetic field strength and temperature of emitting plasma.

(magnetohydrodynamics) MHD↗

Files and scripts to support manuscript Needham et al. Canopy Gradients of Respiration

This dataset includes the parameter files, relevant output files, and scripts to perform analysis with Jupyter notebooks that support the manuscript Needham et al 2025 “Canopy Gradients of Respiration Drive Plant Carbon Budgets and Leaf Area Index.” We add functionality to the Functionally Assembled Terrestrial Ecosystem Simulator (FATES) to allow flexible vertical gradients of leaf maintenance respiration (Rdark) and maximum carboxylation rate (Vcmax) through the canopy. We test the sensitivity of FATES to canopy gradients in Rdark, both in global simulations to assess broad scale impacts on leaf area index (LAI) and vegetation carbon, and in single site simulations where we assess impacts on plant functional type (PFT) competitive dynamics. Parameter files are netcdf files that can be converted to human readable .cdl files using NCO tools. Analysis scripts are Jupyter notebook files. These can be opened and run using the open source Jupyter notebook software. Model outputs are netcdf files.

54 ENVIRONMENTAL SCIENCES↗

Vertical canopy gradients of respiration drive plant carbon budgets and leaf area index

Despite its importance for determining global carbon fluxes, leaf respiration remains poorly constrained in land surface models (LSMs). We tested the sensitivity of the Energy Exascale Earth System Model Land Model – Functionally Assembled Terrestrial Ecosystem Simulator (ELM-FATES) to variation in the canopy gradients of leaf maintenance respiration (R dark ). We ran global and point simulations varying the canopy gradient of R dark to explore the impacts on forest structure, composition, and carbon cycling. In global simulations, steeper canopy gradients of R dark lead to increased understory survival and leaf biomass. Leaf area index (LAI) increased up to 77% in tropical regions compared with the default parameterization, improving alignment with remotely sensed benchmarks. Global vegetation carbon varied from 308 Pg C to 449 Pg C across the ensemble. In tropical forest simulations, steeper gradients of R dark had a large impact on successional dynamics. Results show the importance of canopy gradients in leaf traits and fluxes for determining plant carbon budgets and emergent ecosystem properties such as competitive dynamics, LAI, and vegetation carbon. The high-model sensitivity to canopy gradients in R dark highlights the need for more observations of how leaf traits and fluxes vary along light micro-environments to inform critical dynamics in LSMs.

59 BASIC BIOLOGICAL SCIENCES↗

The role of coverage effects on the structure–sensitivity of formic acid electrooxidation on Pd surfaces

Direct formic acid (FA) fuel cells (DFAFCs) with Pd anode catalysts are promising for small portable applications. During the electrooxidation of formic acid (FAO), CO accumulates on the Pd electrocatalyst, thereby decreasing its performance. In this work, we use density functional theory (DFT; RPBE) to elucidate the significance of coverage effects on the structure–sensitivity of FAO over four low–index Pd surfaces: Pd(1 1 1), Pd(1 0 0), Pd(1 1 0), and Pd(2 1 1). We construct coverage–dependent Gibbs free energy diagrams and in agreement with experiment, predict that Pd(1 0 0) and Pd(2 1 1) are the most active surfaces. Additionally, we find that the presence of CO spectator species can alter the preferred electrooxidation pathway under different applied potentials. Furthermore, such pathway transitions appear to be critical for rationalizing experimental trends regarding electrocatalytic activity and are only captured when coverage effects are accounted for. Insights derived from this work may facilitate the design of improved FAO electrocatalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reviving the carbon sink: The influence of moderate wind disturbance in a secondary temperate mixed forest

Global forests are increasingly exposed to climate-driven perturbations, which may in turn alter their climate mitigation potential. As tropical cyclones expand poleward due to climate warming, wind disturbances in temperate forests have become increasingly frequent. The consequences of moderate wind disturbances remain poorly understood, hindering efforts to quantify their role in the global carbon cycle. Here, we used 16 years of continuous eddy covariance and biometric measurements to investigate the impacts of moderate wind disturbances on the structure and carbon sink dynamics of a temperate forest in Northeast China. Following Typhoon Maysak in 2020, the mortality of large trees (particularly the aging pioneer species) increased ninefold, whereas that of small trees decreased by nearly half. Both stand basal area and leaf area index were reduced between 2019 and 2023, with aging pioneer tree species being more vulnerable than mid-to-late species to wind disturbances. Shifts in species composition altered the environmental sensitivity of forest carbon sink function. Unexpectedly, wind disturbances reversed the declining trends in net ecosystem production and ecosystem carbon use efficiency of this secondary forest. A novel composite structural indicator—the standardized leaf area index (the maximum leaf area supported by per basal area of the stand)—provided robust predictions (R 2 > 0.4) of carbon sink dynamics throughout the study period. The selective removal of less efficient pioneer trees accelerated succession and reversed the aging-related decline in forest carbon sink strength and carbon use efficiency. In conclusion, these findings highlight the potential role of moderate wind disturbances in enhancing forest carbon sink function and offer a framework for understanding, assessing, and predicting forest carbon dynamics under increasing disturbance frequencies driven by climate change.

Carbon sink↗

Simultaneous measurement of surface velocity and plasma density with interferometric velocimetry

The apparent velocity measured by an interferometric surface velocimeter is a function of both the surface velocity and the time derivative of the refractive index along the measurement path. We employed this dual sensitivity to simultaneously measure km/s surface velocities and 10 18 cm -3 average plasma densities with combined VISAR (velocity interferometer system for any reflector) and PDV (photonic Doppler velocimetry) measurements in experiments performed on the Z Pulsed Power Facility. We detail the governing equations, associated assumptions, and analysis specifics and show that the surface velocity can be extracted without knowledge of the specific plasma density profile.

47 OTHER INSTRUMENTATION↗

Accurate prediction of global-density-dependent range-separation parameters based on machine learning

In this work, we develop an accurate and efficient XGBoost machine learning model for predicting the global-density-dependent range-separation parameter, ωGDD, for long-range corrected functional (LRC)-ωPBE. This ωGDDML model has been built using a wide range of systems (11 466 complexes, ten different elements, and up to 139 heavy atoms) with fingerprints for the local atomic environment and histograms of distances for the long-range atomic correlation for mapping the quantum mechanical range-separation values. The promising performance on the testing set with 7046 complexes shows a mean absolute error of 0.001 117 a0−1 and only five systems (0.07%) with an absolute error larger than 0.01 a0−1, which indicates the good transferability of our ωGDDML model. In addition, the only required input to obtain ωGDDML is the Cartesian coordinates without electronic structure calculations, thereby enabling rapid predictions. LRC-ωPBE(ωGDDML) is used to predict polarizabilities for a series of oligomers, where polarizabilities are sensitive to the asymptotic density decay and are crucial in a variety of applications, including the calculations of dispersion corrections and refractive index, and surpasses the performance of all other popular density functionals except for the non-tuned LRC-ωPBE. Finally, LRC-ωPBE (ωGDDML) combined with (extended) symmetry-adapted perturbation theory is used in calculating noncovalent interactions to further show that the traditional ab initio system-specific tuning procedure can be bypassed. The present study not only provides an accurate and efficient way to determine the range-separation parameter for LRC-ωPBE but also shows the synergistic benefits of fusing the power of physically inspired density functional LRC-ωPBE and the data-driven ωGDDML model.

Chemistry↗

Dynamic subcanopy leaf traits drive resistance of net primary production across a disturbance severity gradient

Across the globe, the forest carbon sink is increasingly vulnerable to an expanding array of low- to moderate-severity disturbances. However, some forest ecosystems exhibit functional resistance (i.e., the capacity of ecosystems to continue functioning as usual) following disturbances such as extreme weather events and insect or fungal pathogen outbreaks. Unlike severe disturbances (e.g., stand-replacing wildfires), moderate severity disturbances do not always result in near-term declines in forest production because of the potential for compensatory growth, including enhanced subcanopy production. Community-wide shifts in subcanopy plant functional traits, prompted by disturbance-driven environmental change, may play a key mechanistic role in resisting declines in net primary production (NPP) up to thresholds of canopy loss. However, the temporal dynamics of these shifts, as well as the upper limits of disturbance for which subcanopy production can compensate, remain poorly characterized. In this study, we leverage a 4-year dataset from an experimental forest disturbance in northern Michigan to assess subcanopy community trait shifts as well as their utility in predicting ecosystem NPP resistance across a wide range of implemented disturbance severities. Through mechanical girdling of stems, we achieved a gradient of severity from 0% (i.e., control) to 45, 65, and 85% targeted gross canopy defoliation, replicated across four landscape ecosystems broadly representative of the Upper Great Lakes ecoregion. We found that three of four examined subcanopy community weighted mean (CWM) traits including leaf photosynthetic rate (p = 0.04), stomatal conductance (p = 0.07), and the red edge normalized difference vegetation index (p < 0.0001) shifted rapidly following disturbance but before widespread changes in subcanopy light environment triggered by canopy tree mortality. Surprisingly, stimulated subcanopy production fully compensated for upper canopy losses across our gradient of experimental severities, achieving complete resistance (i.e., no significant interannual differences from control) of whole ecosystem NPP even in the 85% disturbance treatment. Additionally, we identified a probable mechanistic switch from nutrient-driven to light-driven trait shifts as disturbance progressed. Our findings suggest that remotely sensed traits such as the red edge normalized difference vegetation index (reNDVI) could be particularly sensitive and robust predictors of production response to disturbance, even across compositionally diverse forests. The potential of leaf spectral indices to predict post-disturbance functional resistance is promising given the capabilities of airborne to satellite remote sensing. We conclude that dynamic functional trait shifts following disturbance can be used to predict production response across a wide range of disturbance severities.

resistance↗

BeyondPlanck: XV. Limits on large-scale polarized anomalous microwave emission from Planck LFI and WMAP

We constrained the level of polarized anomalous microwave emission (AME) on large angular scales using Planck Low-Frequency Instrument (LFI) and WMAP polarization data within a Bayesian cosmic microwave background (CMB) analysis framework. We modeled synchrotron emission with a power-law spectral energy distribution, as well as the sum of AME and thermal dust emission through linear regression with the Planck High-Frequency Instrument (HFI) 353 GHz data. This template-based dust emission model allowed us to constrain the level of polarized AME while making minimal assumptions on its frequency dependence. We neglected CMB fluctuations, but show through simulations that these fluctuations have a minor impact on the results. We find that the resulting AME polarization fraction confidence limit is sensitive to the polarized synchrotron spectral index prior. In addition, for prior means β s < -3.1 we find an upper limit of p AME max ≲ 0.6% (95% confidence). In contrast, for means β s = -3.0, we find a nominal detection of p AME = 2.5 ± 1.0% (95% confidence). These data are thus not strong enough to simultaneously and robustly constrain both polarized synchrotron emission and AME, and our main result is therefore a constraint on the AME polarization fraction explicitly as a function of β s . Combining the current Planck and WMAP observations with measurements from high-sensitivity low-frequency experiments such as C-BASS and QUIJOTE will be critical to improve these limits further.

79 ASTRONOMY AND ASTROPHYSICS↗

Machine-Learning Functional Zonation Approach for Characterizing Terrestrial–Aquatic Interfaces: Application to Lake Erie

Ecosystems at coastal terrestrial–aquatic interfaces play a significant role in global biogeochemical cycles. In this study, we aimed to characterize coastal wetlands with particular focus on the co-variability between plant dynamics, topography, soil, and other environmental factors. We proposed a functional zonation approach based on machine learning clustering to identify the spatial regions, i.e., zones that capture these co-varied properties. This approach was applied to publicly available datasets along Lake Erie, in the Great Lakes Region. We investigated the heterogeneity of coastal ecosystem structures as a function of along-shore distance and transverse distance, based on the spatial data layers, including topography, wetland vegetation cover, and the time series of Landsat’s enhanced vegetation index (EVI) between 1990 and 2020. Results showed that the topographic metrics (elevation and slope), soil texture, and plant productivity influence the spatial distribution of wetland land-covers (emergent and phragmites). These results highlight a natural organization along the transverse axis, where the elevation and the EVI increase further away from the coastline. In addition, the clustering analysis allowed us to identify regions with distinct environmental characteristics, as well as the ones that are more sensitive to interannual lake-level variations.

47 OTHER INSTRUMENTATION↗

Detecting Multipartite Entanglement Patterns Using Single-Particle Green’s Functions

Here, we present a protocol for detecting multipartite entanglement in itinerant many-body electronic systems using single-particle Green’s functions. To achieve this, we first establish a connection between the quantum Fisher information and single-particle Green’s functions by constructing a set of witness operators built out of single electron creation and destruction operators in a doubled system. This set of witness operators is indexed by a momentum k. We compute the quantum Fisher information for these witness operators and show that for thermal ensembles it can be expressed as an autoconvolution of the single-particle spectral function. We then apply our framework to a one-dimensional fermionic system to showcase its effectiveness in detecting entanglement in itinerant electron models. We observe that the detected entanglement level is sensitive to the wave vector associated with witness operator. Our protocol will permit detecting entanglement in many-body systems using scanning tunneling microscopy and angle-resolved photoemission spectroscopy, two spectroscopies that measure the single-particle Green’s function. It offers the prospect of the experimental detection of entanglement through spectroscopies beyond the established route of measuring the dynamical spin response.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Mineral Deposition on the Rough Walls of a Fracture

Modeling carbonate growth in fractures and pores is important for understanding carbon sequestration in the environment or when supersaturated solutions are injected into rocks. Here, we study the simple but nontrivial problem of calcite growth on fractures with rough walls of the same mineral using kinetic Monte Carlo simulations of attachment and detachment of molecules and scaling approaches. First, we consider wedge-shaped fracture walls whose upper terraces are in the same low-energy planes and show that the valleys are slowly filled by the propagation of parallel monolayer steps in the wedge sides. The growth ceases when the walls reach these low-energy configurations so that a gap between the walls may not be filled. Second, we consider fracture walls with equally separated monolayer steps (vicinal surfaces with roughness below 1 nm) and show that growth by step propagation will eventually clog the fracture gap. In both cases, scaling approaches predict the times to attain the final configurations as a function of the initial geometry and the step-propagation velocity, which is set by the saturation index. The same reasoning applied to a random wall geometry shows that step propagation leads to lateral filling of surface valleys until the wall reaches the low-energy crystalline plane that has the smallest initial density of molecules. Thus, the final configurations of the fracture walls are much more sensitive to the crystallography than to the roughness or the local curvature. The framework developed here may be used to determine those configurations, the times to reach them, and the mass of deposited mineral. Effects of transport limitations are discussed when the fracture gap is significantly narrowed.

calcite↗

High-throughput identification of novel heat tolerance genes via genome-wide pooled mutant screens in the model green alga Chlamydomonas reinhardtii

Different high temperatures adversely affect crop and algal yields with various responses in photosynthetic cells. The list of genes required for thermotolerance remains elusive. Additionally, it is unclear how carbon source availability affects heat responses in plants and algae. Here, we utilized the insertional, indexed, genome-saturating mutant library of the unicellular, eukaryotic green alga Chlamydomonas reinhardtii to perform genome-wide, quantitative, pooled screens under moderate (35°C) or acute (40°C) high temperatures with or without organic carbon sources. We identified heat-sensitive mutants based on quantitative growth rates and identified putative heat tolerance genes (HTGs). By triangulating HTGs with heat-induced transcripts or proteins in wildtype cultures and MapMan functional annotations, we presented a high/medium-confidence list of 933 Chlamydomonas genes with putative roles in heat tolerance. Triangulated HTGs include those with known thermotolerance roles and novel genes with little or no functional annotation. About 50% of these high-confidence HTGs in Chlamydomonas have orthologs in green lineage organisms, including crop species. Arabidopsis thaliana mutants deficient in the ortholog of a high-confidence Chlamydomonas HTG were also heat sensitive. This work expands our knowledge of heat responses in photosynthetic cells and provides engineering targets to improve thermotolerance in algae and crops.

59 BASIC BIOLOGICAL SCIENCES↗

Chemical Kinetics Underlying the Sooting Tendency and Auto-Ignition Characteristics of Linear, Branched, and Cyclic Ether Compounds

Biofuels present opportunities for improving the performance and reducing emissions from internal combustion engines by incorporating oxygenated functional groups to the fuels. Among various oxygenates, ethers have been recognized as promising candidates for an alternative to conventional diesel fuel owing to their higher reactivity and lower sooting tendency. The detailed guidelines for designing ethers, however, have not been fully discussed, even though their combustion characteristics are sensitive to the molecular structure. This study was devoted to exploring the structure-property relationships, particularly focusing on the cetane number and yield sooting index, using five linear, branched, and cyclic ethers: di-amyl-ether, 4-butoxy-heptane, 3,3-dimethyl-oxetane, 2-ethyl-4-methyl-1,3-dioxolane, and 2-isopropyl-4-methyl-1,3-dioxolane. First, we examined the chemical kinetics underlying the sooting tendency of the test fuels. The combustion product distribution was measured from flow reactor experiments at 750-1100 K, F=3, at atmospheric pressure. As a result, it was revealed that the sooting tendency is closely related to the size of hydrocarbon intermediates in the high-temperature regime (>1000 K); that is, larger hydrocarbons lead to more soot precursor formation. The underlying chemistry determining the size of the hydrocarbons from the tested fuels was analyzed using reaction pathway analysis and quantum mechanics calculations, which showed that the branched and cyclic ether structures form abundant C3-C4 compounds. Moreover, the auto-ignition characteristics of the test fuels were studied using the flow reactor at low-temperature (400-700 K) and F=1. We found a clear difference in the combustion-product distribution from high and low reactivity fuels, which was then correlated to the systematic analysis of the key reaction energy barriers with the varying molecular structure.

BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗

Imaging Photonic Resonances within an All‐Dielectric Metasurface via Photoelectron Emission Microscopy

Dielectric metasurfaces, through volume‐type photonic resonances, enable precise control of light‐matter interactions for applications including imaging, holography, and sensing. The application space of dielectric metasurfaces has extended from infrared to visible wavelengths by incorporating high refractive index materials, such as titanium dioxide (TiO 2 ). Understanding the fundamental and fabrication limits for these applications requires metrology with nanoscale resolution, sensitivity to electromagnetic fields within the meta‐atom volume, and far‐field excitation. In this work, photoelectron emission microscopy (PEEM) is used to image field distributions of photonic resonances in a TiO 2 metasurface excited with far‐field, visible‐wavelength illumination. The local volumetric field variations within the meta‐atoms are analyzed as a function of illumination angle and polarization by comparing photoelectron images to finite‐difference time‐domain simulations. This study determines the inelastic mean free path of very low‐energy (<1 eV) photoelectrons to be 35 ± 10 nm, which is comparable to the meta‐atom height thereby highlighting PEEM sensitivity to resonances within the volume. Additionally, the simulations reveal high sensitivity of PEEM images to an in‐plane component of the illumination k ‐vector. These results demonstrate that photoelectron imaging with subwavelength resolution offers unique advantages for examining light‐matter interactions in volume‐type (as opposed to surface) photonic modes within dielectric nanophotonic structures.

electron inelastic mean free path↗