Search NASA⌕ Search

SEARCH · Search NASA

Results for “Mechanistic Modeling”

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 271 records · Page 15

Recent Trends in Global Ocean Chlorophyll

Recent analyses of SeaWiFS data have shown that global ocean chlorophyll has increased more than 5% since 1998. The North Pacific ocean basin has increased nearly 19%. To understand the causes of these trends we have applied the newly developed NASA Ocean Biogeochemical Assimilation Model (OBAM), which is driven in mechanistic fashion by surface winds, sea surface temperature, atmospheric iron deposition, sea ice, and surface irradiance. The mode1 utilizes chlorophyll from SeaWiFS in a daily assimilation. The model has in place many of the climatic variables that can be expected to produce the changes observed in SeaWiFS data. Th~s enables us to diagnose the model performance, the assimilation performance, and possible causes for the increase in chlorophyll.

Gregg, Watson↗

Microbial inoculants and invasions: a call to action

Microbial inoculants are increasingly used for beneficial purposes in agriculture, bioremediation, and medicine, but they can carry risks of generating invasive microbes. Here, we present a roadmap for guarding against these invasions, proposing developing (i) coherent mechanistic understandings of how microbial inoculants can effect invasions, (ii) predictive models forecasting microbial invasion risks, and (iii) effective management strategies. To guide mechanistic understandings, we distill 17 guiding hypotheses. For predictive modeling, we highlight data collection needs and qualitative approaches. For management strategies, we stress the importance of accurately weighing the risks against benefits. The unified approach presented here provides a route toward an effective research and management infrastructure for microbial inoculants in order to avoid potentially catastrophic microbial invasions.

invasive species↗

Applicability of 1D site response analysis to shallow sedimentary basins: A critical evaluation through physics-based 3D ground motion simulations

One-dimensional site response analysis (1D SRA) remains the standard practice in considering the effect of local soil deposits and predicting site-specific ground motions, although its range of applicability to realistic seismic wavefields is still in question. In this 1D approach, horizontal and vertical ground shaking are assumed to be induced by vertically propagating shear and compressional waves, respectively. A recent study based on analytical two-dimensional (2D) plane waves and simple point source earthquake simulations has shown two mechanistic limitations in this 1D modelling technique for general inclined seismic waves, that is, systematic over-prediction of the vertical motion and wave trapping in the 1D soil column. In this article, we evaluate in detail the applicability of this 1D modelling approach to realistic three-dimensional (3D) simulated seismic wavefields in shallow sedimentary basins. Linear-viscoelastic 1D SRA predictions using two types of input motions that are commonly used in practice—rock outcrop and in-column motions, are compared with the reference true site response results from 3D earthquake simulations in terms of various measures in the frequency and time domain. It is shown that the horizontal motion in the 3D seismic wavefield exhibits dominant shear wave propagation phenomenon, while the vertical motion is a combined effect of compressional and shear waves and can be over-predicted by the 1D approach when the incident seismic waves are inclined. Direct evidence of the wave refraction process that leads to the vertical motion over-prediction is provided. 1D SRA with in-column inputs can yield motions that have significantly longer duration compared to the true 3D site response solution due to trapped waves, casting in doubt the frequent need for increased soil damping in existing site studies to compensate for wave attenuation due to scattering alone. Sensitivity investigation on the increase of soil profile damping by a multiplier D mul shows D mul values compatible with those found in the literature for both horizontal and vertical motions. It is shown that the level of D mul optimized for a best match of the spectral acceleration is dependent on the characteristic of the input motion and a larger D mul is typically required for the vertical component. In contrast, 1D SRA with outcrop motions predicts motions with shorter significant duration due to its inability to capture the basin-edge generated surface waves. A suite of ground motion simulations was performed to assess the sensitivity of the observations to the basin geologic structure including the velocity gradient, rock-basin impedance contrast and basin depth. The analysis results show that the accuracy of the simplified 1D procedure is dependent on the wavefield composition of both the input motions and the true 3D site response solution. While the horizontal motions in shallow sedimentary basins can, to the first order, be reasonably captured by the simplified 1D approach, 1D SRA for the vertical component is in general not reliable and contributions from inclined shear waves should be accounted for in site-specific evaluation of the vertical design ground motion.

1D assumption↗

Selective phosphate removal with manganese oxide composite anion exchange membranes in membrane capacitive deionization

The discharge of excessive phosphorous into water bodies can lead to serious eutrophication threatening aquatic ecosystem. Membrane capacitive deionization (MCDI) is an effective platform for deionizing aqueous streams; however, conventional MCDI is unable to selectively remove targeted ions from a liquid mixture. Here, in this work, we fabricated manganese oxide composite anion exchange membranes (AEMs) for MCDI to enhance phosphate removal selectivity from sodium chloride-sodium dihydrogen phosphate (10:1 M ratio) aqueous mixtures. We systematically investigated several critical factors, such as constant current or voltage operation, applied voltage amount, process stream pH, and manganese oxide (Mn 2 O 3 ) content in the AEM, on phosphate removal efficiency and phosphate selectivity. A trade-off was observed between phosphate removal and selectivity when increasing the cell voltage. Under the best conditions, a MCDI unit with a 20 wt% Mn 2 O 3 composite AEM and a bipolar membrane facilitated high phosphate removal efficiency of ≥ 31.8 % and a phosphate over chloride selectivity of 1.1 while showing stability for at least 30 cycles. To help understand how Mn 2 O 3 composite AEM boosts phosphate selectivity, static electronic structure calculations were performed, and they revelated that hydrogen phosphate absorption on Mn 2 O 3 composite AEM was 314 kcal/mol more exothermic than that on pristine AEM while chloride adsorption on Mn 2 O 3 composite AEM was 2.2 kcal/mol less exothermic than that on a pristine AEM. Overall, this work presents an effective strategy for selectively removing phosphate from model wastewater solutions and the mechanistic understanding that governs ion selectivity in composite ion-exchange membranes used in MCDI.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanisms of Three-Dimensional Solid-Phase Epitaxial Crystallization of Strontium Titanate

Strontium titanate (SrTiO 3 , STO) is a complex metal oxide with a cubic perovskite crystal structure. Due to its easily described and understood crystal structure in the cubic phase, STO is an ideal model system for exploring the mechanistic details of solid-phase epitaxy (SPE) in complex oxides. SPE is a crystallization approach that aims to guide crystal growth at low homologous temperatures to achieve targeted microstructures. Beyond planar thin films, SPE can also exploit the addition of a chemically inert, noncrystallizing, amorphous obstacle in the path of crystallization to generate complex three-dimensional structures. The introduction of this mask fundamentally alters the SPE process, inducing a transition from two- to three-dimensional geometries and from vertical to lateral crystal growth under the influence of the crystal/mask/amorphous boundary. Using a combination of molecular dynamics simulations and experiments, we identify several unique phenomena in the nanoscale growth behaviors in both conventional (unmasked) and masked SPE. Examining conventional SPE of STO, we find that crystallization at the interface is strongly correlated to, and potentially driven by, density fluctuations in the region of the amorphous STO near the crystalline/amorphous interface with a strong facet dependence. In the masked case, we find that the crystalline growth front becomes nonplanar near contact with the mask. We also observe a minimum vertical growth requirement prior to lateral crystallization. Both phenomena depend on the relative bulk and interfacial free energies of the three-phase (crystal/mask/amorphous) system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Real-Space Constrained Density Functional Theory Investigation of Site-Specific, Interfacial Charge Recombination Dynamics Across the Au Nanoparticle/TiO 2 Heterojunction

Au nanoparticle (NP)/TiO 2 heterojunction is a representative system to study interfacial charge transfer in photocatalysis and photovoltaics, where suppressing recombination from TiO 2 to Au can enhance hot carrier extraction. We apply real-space constrained density functional theory (CDFT) with Marcus theory to quantify charge recombination time scales across Au/TiO 2 . This approach enables direct control and visualization of charge-separated states, aligning with site-specific probes like time-resolved X-ray photoelectron spectroscopy (trXPS). We find that the charge-separated state features a bipolaron, with recombination dominated by TiO 2 LUMO to Au HOMO transitions, primarily at interfacial Au sites. Marcus rate predictions are benchmarked with surface hopping methods, quantifying differences in time scales and computational efficiency. Lastly, we examine how the Au cluster size affects the free energy change (ΔG) and reorganization energy (λ), explaining trends in closed-shell systems and highlighting challenges for open-shell extrapolations. Overall, CDFT + Marcus theory provides efficient, mechanistically transparent interfacial charge transfer modeling, and we clearly defined its applicability and limitation.

Glenna, Drew M. [Univ. of Idaho, Idaho Falls, ID (↗

Thermally driven surface phase separation in intermetallic alloys

Intermetallic compounds are widely recognized for their high-temperature phase stability and resistance to composition and structural changes. However, we reveal athermally activated bulk-to-surface mass exchange mechanism that drives surface phase separation, resulting in the formation of surface precipitates with distinct composition andstructure from the bulk matrix. Using the archetypal β-NiAl system, we show that asymmetries in vacancy formation energies between Ni and Al atoms induce preferential Ni segregation to the surface, forming Ni-rich γ'-Ni 3 Al precipitates. By integrating in-situ electron microscopy, synchrotron X-ray absorption spectroscopy and first-principles computational modeling, we establish a direct mechanistic connection between bulk thermal defect dynamics, surface compositional evolution, and phase segregation behavior. This bulk-surface coupling mechanism can be a driver of surface phase separation in multicomponent alloys under thermal stress. In conclusion, these results refine the thermodynamic boundaries of intermetallic stability and provide insights into managingthe performance and durability of intermetallic alloys for demanding high-temperature applications.

36 MATERIALS SCIENCE↗

Mesoscale density variability in the mesosphere and thermosphere: Effects of vertical flow accelerations

A mechanistic one dimensional numerical (iteration) model was developed which can be used to simulate specific types of mesoscale atmospheric density (and pressure) variability in the mesosphere and the thermosphere, namely those due to waves and those due to vertical flow accelerations. The model was developed with the idea that it could be used as a supplement to the TGCMs (thermospheric general circulation models) since such models have a very limited ability to model phenomena on small spatial scales. The simplest case to consider was the integration upward through a time averaged, height independent, horizontally divergent flow field. Vertical winds were initialized at the lower boundary using the Ekman pumping theory over flat terrain. The results of the computations are summarized.

Revelle, D. O.↗

Kinetics and Mechanism of the Reaction of Hydoxyl Radicals with Acetonitrile under Atmospheric Conditions

The pulsed laser photolysis-pulsed laser induced fluorescence technique has been employed to determine absolute rate coefficients for the reaction OH + CH3CN (1) and its isotopic variants, OH + CD3CN (2), OD + CH3CN (3), and OD + CD3CN (4). Reactions 1 and 2 were studied as a function of pressure and temperature in N2, N2/O2, and He buffer gases. In the absence of O2 all four reactions displayed well-behaved kinetics with exponential OH decays and pseudo-first rate constants which were proportional to substrate concentration. Data obtained in N2 over the range 50-700 Torr at 298 K are consistent with k(sub 1), showing a small pressure dependence. The Arrhenius expression obtained by averaging data at all pressures in k(sub 1)(T) = (1.1(sup +0.5)/(sub -0.3)) x 10(exp -12) exp[(-1130 +/- 90)/T] cu cm /(molecule s). The kinetics of reaction 2 are found to be pressure dependent with k(sub 2) (298 K) increasing from (1.21 +/- 0.12) x 10(exp -14) to (2.16 +/- 0.11) x 10(exp -14) cm(exp 3)/ (molecule s) over the pressure range 50-700 Torr of N2 at 298 K. Data at pressures greater than 600 Torr give k(sub 2)(T) = (9.4((sup +13.4)(sub -5.0))) x 10(exp -13) exp[(-1180 +/- 250)/T] cu cm/(molecule s). The rates of reactions 3 and 4 are found to be independent of pressure over the range 50-700 Torr of N2 with 298 K rate coefficient given by k(sub 3) =(3.18 +/- 0.40) x 10(exp -14) cu cm/(molecule s) and k(sub 4) = (2.25 +/-0.28) x 10(exp -14) cu cm/(molecule s). In the presence of O2 each reaction shows complex (non-pseudo-first-order) kinetic behavior and/or an apparent decrease in the observed rate constant with increasing [O2], indicating the presence of significant OH or OD regeneration. Observation of regeneration of OH in (2) and OD in (3) is indicative of a reaction channel which proceeds via addition followed by reaction of the adduct, or one of its decomposition products, with O2. The observed OH and OD decay profiles have been modeled by using a simple mechanistic scheme to extract kinetic information about the adduct reations with O2 and branching ratios for OH regeneration. A plausible mechanism for OH regeneration in (2) involves OH addition to the nitrogen atom followed by O2 addition to the cyano carbon atom, isomeriazation and decomposition to D2CO + DOCN + OH. Our results suggest that the OH + CH3CN reaction occurs via a complex mechanism involving both bimolecular and termolecular pathways, analogous to the mechanisms for the the important atmospheric reactions of OH with CO and HNO3.

Hynes, A. J.↗

Terrestrial Models and Global Change: Challenges for the Future

A wide variety of models have illustrated the potential importance of terrestrial biological feedbacks on climate and climate change, yet our ability to make precise predictions is severely limited, due to a high degree of uncertainty. In this paper, after briefly reviewing current models, we present challenges for new terrestrial models and introduce a simple mechanistic approach that may complement existing approaches.

Hurtt, George C.↗

Airborne Dust Plumes Lofted by Dislodged Ice Blocks at Russell Crater, Mars

Linear dune gullies on poleward-facing Martian slopes are enigmatic. Formation by CO2-ice block or snow cornice falls has been proposed based on optical imagery of bright, high-albedo features inside gully channels. Because these features often resemble patchy frost residue rather than three-dimensional blocks, more evidence is needed to support the ice-block formation mechanism. Satellite imagery captured two simultaneous airborne plumes with in-channel sources at the Russell crater megadune, thrust up, and dispersed outward along the path of linear dune gullies. We use spectral data analyses, climatic analyses of bolometric temperatures, and thermal modeling to further develop the mechanistic framework for linear dune gully development. Basal sublimation and CO2 gas venting likely cause CO2-ice-block detachment and falls from gully alcoves in southern early spring, accompanied by ice-block off-gassing and saltation of sands and coarse silts that are redeposited around gully channels, and lofting of sublimation lag (coarse dust/silt) into airborne plumes.

Cynthia Dinwiddie↗

Dynamical Sensitivity in Response to a Wide Range of Abrupt CO 2 Forcings

An improved understanding of dynamical variability is important for projecting future changes in extratropical weather and the interaction between the extratropical troposphere and the tropics, stratosphere, cryosphere and ocean. Despite their simplicity, the abrupt 2x- and 4xCO 2 forcing simulations from the Coupled Model Intercomparison Project (CMIP) Phase 6 DECK experiments enable a mechanistic look into the forcing and feedback response characteristics of models that can be unambiguously attributed to an increase in carbon dioxide concentrations. Thus, while typically used to evaluate the climate sensitivity in models, with a primary focus on global surface temperature change, here we focus on several measures of extratropical variability, including projected changes in the storm tracks and in stratospheric polar vortex variability. Results are primarily based on simulations produced using low- and high-top versions of the NASA Goddard Institute for Space Studies Model (ModelE) but results are also presented from the larger CMIP6 multi-model ensemble. In addition to the 2x- and 4xCO 2 simulations, we also explore the linearity of the response of extratropical dynamical variability in ModelE to varying levels of CO 2 spanning the range 1/8-8xCO 2 . In particular, we show that the expansion of the jet streams varies nonlinearly with increasing CO 2 , especially in the Northern Hemisphere, and that this can be interpreted in terms of variations in ocean heat transport. The impact of composition feedbacks on changes in variability is also discussed.

dynamical variability↗

Tree Growth Enhancement Drives a Persistent Biomass Gain in Unmanaged Temperate Forests

While enhanced tree growth over the last decades has been reported in forests across the globe, it remains unclear whether it drives persistent biomass increases of forest stands, particularly in mature forests. Enhanced tree growth and stand-level biomass are often linked with a simultaneous increase in density-driven mortality and a reduction in tree longevity. Identifying empirical evidence regarding the balance between these processes is challenging due to the confounding effects of stand history, management, and environmental changes. Here, we investigate the link between growth and biomass via the negative relationship between average tree size and stand density (tree number per area). We find increasing stand density for a given mean tree size in unmanaged closed-canopy forests in Switzerland over the past six decades and a positive relationship between tree growth and stand density across forest plots—qualitatively consistent with our simulations using a mechanistic, cohort-resolving ecosystem model (BiomeE). Model simulations show that, in the absence of other disturbances, enhanced tree growth persistently increases biomass stocks despite simultaneous decreases in carbon residence time and tree longevity. However, the magnitude of simulated biomass changes for a given growth enhancement critically depends on the shape of the mortality functions. Our analyses reconcile reports of growth-induced reductions of tree longevity with model predictions of persistent biomass increases, and with our finding of trends toward denser forests in response to growth—also in mature stands.

biomass stocks↗

Towards a liana plant functional type for vegetation models

Lianas (woody climbers) are crucial components of tropical forests and they have been increasingly recognized to have profound effects on tropical forest carbon dynamics. Despite their importance, lianas' representation in vegetation models remains limited, partly due to the complexity of liana-tree dynamics and the diversity in liana life history strategies. This paper provides a comprehensive review of advances and challenges for mechanistically representing lianas in forest ecosystem models and a proposed path towards effectively representing lianas in these models. Defining a liana plant functional type is a significant challenge because of the high morphological and physiological diversity amongst liana species, and because of their structural association with trees. Here, we identify critical liana traits that likely should contribute to establishing a liana plant functional type, along with key processes to properly represent lianas in ecosystem models. Subsequently, we discuss a variety of possible liana implementation strategies with their associated strengths, limitations, computational costs and data requirements. A fundamental redesign of the tree-centric demographic vegetation models seems appropriate to accommodate the unique growth and competition strategies of lianas. We illustrate the potential of such models with a single-site case study where we disentangle putative mechanisms of liana increasing abundance. Furthermore, we underscore the critical need for comprehensive liana demographic and functional data (including long-term, physiological, and pantropical observations) for the qualitative implementation and evaluation in the proposed modeling efforts. Currently, there is a scarcity of liana data and the data that do exist have a neotropical bias. We finally introduce a new liana functional trait database that can centralize existing liana trait data, incentivize improved data gathering and thus facilitate model development and scientific analyses.

54 ENVIRONMENTAL SCIENCES↗

Comparison of kinetic models for atom recombination on high-temperature reusable surface insulation

Five kinetic models are compared for their ability to predict recombination coefficients for oxygen and nitrogen atoms over high-temperature reusable surface insulation (HRSI). Four of the models are derived using Rideal-Eley or Langmuir-Hinshelwood catalytic mechanisms to describe the reaction sequence. The fifth model is an empirical expression that offers certain features unattainable through mechanistic description. The results showed that a four-parameter model, with temperature as the only variable, works best with data currently available. The model describes recombination coefficients for oxygen and nitrogen atoms for temperatures from 300 to 1800 K. Kinetic models, with atom concentrations, demonstrate the influence of atom concentration on recombination coefficients. These models can be used for the prediction of heating rates due to catalytic recombination during re-entry or aerobraking maneuvers. The work further demonstrates a requirement for more recombination experiments in the temperature ranges of 300-1000 K, and 1500-1850 K, with deliberate concentration variation to verify model requirements.

Willey, Ronald J.↗

A Unified Framework to Reconcile Different Approaches of Modeling Transpiration Response to Water Stress: Plant Hydraulics, Supply Demand Balance, and Empirical Soil Water Stress Function

Plant responses to water stress is a major uncertainty to predicting terrestrial ecosystem sensitivity to drought. Different approaches have been developed to represent plant water stress. Empirical approaches (the empirical soil water stress (or Beta) function and the supply-demand balance scheme) have been widely used for many decades; more mechanistic based approaches, that is, plant hydraulic models (PHMs), were increasingly adopted in the past decade. However, the relationships between them—and their underlying connections to physical processes—are not sufficiently understood. This limited understanding hinders informed decisions on the necessary complexities needed for different applications, with empirical approaches being mechanistically insufficient, and PHMs often being too complex to constrain. Here we introduce a unified framework for modeling transpiration responses to water stress, within which we demonstrate that empirical approaches are special cases of the full PHM, when the plant hydraulic parameters satisfy certain conditions. We further evaluate their response differences and identify the associated physical processes. Finally, we propose a methodology for assessing the necessity of added complexities of the PHM under various climatic conditions and ecosystem types, with case studies in three typical ecosystems: a humid Midwestern cropland, a semi-arid evergreen needleleaf forest, and an arid grassland. Notably, Beta function overestimates transpiration when VPD is high due to its lack of constraints from hydraulic transport and is therefore insufficient in high VPD environments. With the unified framework, we envision researchers can better understand the mechanistic bases of and the relationships between different approaches and make more informed choices.

54 ENVIRONMENTAL SCIENCES↗

Elucidating Phosphate and Cadmium Cosorption Mechanisms on Mineral Surfaces with Direct Spectroscopic and Modeling Evidence

The simultaneous sorption of cations and anions at the mineral–water interface can substantially alter their individual sorption characteristics; however, this phenomenon lacks a mechanistic understanding. Our study provides direct spectroscopic and modeling evidence of the molecular cosorption mechanisms of the cadmium ion (Cd 2+ ) and phosphate (P) on goethite and layered manganese (Mn) oxide of birnessite, through in situ attenuated total reflection Fourier-transform infrared (ATR-FTIR), P K-edge X-ray absorption near-edge structure (XANES) spectroscopy, and surface complexation modeling. Phosphate synergistically cosorbed with Cd on goethite predominantly through P-bridged ternary complexes (≡Fe–P–Cd) and electrostatic interactions at wide pH conditions. Likewise, P and Cd exhibited synergistic cosorption on birnessite by forming P-bridged ternary complexes (≡Mn–P–Cd) and weak competitive sorption at the layer edge sites. As pH and Cd loading increased, the surface P species transitioned from a binary complex to a ternary complex and/or Cd 3 (PO 4 ) 2 precipitate for both goethite and birnessite. Compared to that in solution at pH 8, the formation of Cd 3 (PO 4 ) 2 was inhibited by the presence of goethite and birnessite, ascribed to the specific adsorption of P and Cd, more pronounced in birnessite due to the stronger sorption of Cd at its vacant sites. Finally, the discovered cosorption mechanisms of P and Cd have important implications for understanding and predicting their mobility and availability in Cd-contaminated settings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanistic Insights into the Demethylation of Lignin‐Derived Structures Using Protic Ionic Liquids: A Density Functional Theory Study

Lignin valorization is restricted by the stability of its methoxy groups, creating a critical need for efficient demethylation strategies. Here, in this study, density functional theory (DFT) is employed to dissect the mechanistic pathways of demethylation in lignin model compounds, guaiacol and syringol, using protic ionic liquids (PILs) that act as both solvent and catalyst. Conductor‐like Screening Model for Real Solvents (COSMO‐RS) analysis identifies monoethanolammonium acetate ([MEOA][Ace]) as the most promising medium, attributed to its strong hydrogen bonding network and solvation ability. By integrating implicit and explicit solvation models, it is revealed that an acid‐catalyzed hydrolytic mechanism governs demethylation, with PILs stabilizing crucial transition states and intermediates. Complementary electronic structure evaluations, including highest occupied molecular orbital–lowest unoccupied molecular orbital (HOMO‐LUMO) gap analysis, charge distribution, and electrostatic potential mapping, demonstrate how PILs lower energetic barriers and enhance reactivity. To simulate realistic environments, this study is extended to lignin dimer complexes with varying water content, uncovering how water‐bridged solvation reverses demethylation preference from guaiacyl to syringyl units. This mechanistic shift aligns with experimental observations showing faster S‐unit reactivity in hydrated systems. Together, these findings provide atomic‐level insight into lignin demethylation dynamics and highlight how tuning acid concentration in PILs can accelerate kinetics, enabling a rational pathway toward next‐generation biomass conversion technologies.

COSMO-RS↗