Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “invariance”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Growth and Performance of High-Quality SWCNT Forests on Inconel Foils as Lithium-Ion Battery Anodes

Large-scale production of vertically aligned single-walled carbon nanotubes (VA-SWCNTs) on metal foils promises to enable technological advancements in many fields, from functional composites to energy storage to thermal interfaces. In this work, we demonstrate growth of high-quality (G/D > 6, average diameters ~ 2–3 nm, densities > 10 12 cm –2 ) VA-SWCNTs on Inconel metal for use as a lithium-ion battery (LIB) anode. Scale-up of SWCNT growth on Inconel 625 to 100 cm 2 exhibits nearly invariant CNT structural properties, even when synthesis is performed near atmospheric pressure, and this robustness is attributed to a growth kinetic regime dominated by the carbon precursor diffusion in the bulk gas mixture. SWCNT forests produced on large-area metal substrates at close to atmospheric pressure possess a combination of structural features that are among the best demonstrated so far in the literature for growth on metal foils. Leveraging these achievements for energy applications, we demonstrate a VA-SWCNT LIB anode with capacity >1200 mAh/g at 1.0C and stable cycling beyond 300 cycles. Furthermore, this robust synthesis of high-quality VA-SWCNTs on metal foils presents a promising route toward mass production of high-performance CNT devices for a broad range of applications.

Energy↗
DOE OSTI2025

Growth of Hexagonal Boron Nitride from Molten Nickel Solutions: A Reactive Molecular Dynamics Study

Metal flux methods are excellent for synthesizing high-quality hexagonal boron nitride (hBN) crystals, but the atomic mechanisms of hBN nucleation and growth in these systems are poorly understood and difficult to probe experimentally. Here, we harness classical reactive molecular dynamics (ReaxFF) to unravel the mechanisms of hBN synthesis from liquid nickel solvent over time scales up to 30 ns. These simulations mimic experimental conditions by including relatively large liquid nickel slabs containing dissolved boron and a molecular nitrogen gas phase. Overall, the reaction takes place almost exclusively on the surface of the liquid nickel, owing to the low solubility of nitrogen in bulk nickel and the intermediate species’ preference for the metal–gas interface. The formation of hBN invariably begins by reaction of dinitrogen with nickel-solvated boron atoms at the surface, forming intermediate N–N–B species, which typically evolve into B–N–B units through a short-lived intermediate where a single nitrogen atom is coordinated by one nitrogen and two boron atoms. The resulting B–N–B units, in turn, coalesce with growing hBN nuclei and carry nitrogen between hBN nanocrystals in an Ostwald ripening process. The amount of hBN produced on the tens of nanosecond time scale depends critically on the boron concentration, while having a much weaker dependence on the N 2 pressure for the regime considered (N 2 pressures of 2.5–10 MPa, Ni–B solutions with 6–12% boron by atom fraction). The highest rate of hBN formation occurs at the lowest temperature considered (1750 K, just above the melting point of nickel), while no hBN sheets are formed at 2000 K or above. An analysis of the transition pathways for nitrogen atoms shows that the final step, incorporation of small B–N motifs into larger hBN sheets, is the rate-limiting step in the regimes considered. While raising the temperature from 1750 to 2000 K has little effect on the formation of intermediates (N–N–B, B–N–B, etc.), the lack of large hBN sheets at temperatures >1900 K is explained by decreased probability of the final step and increased probability of breakup of hBN into B–N motifs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Mechanism of Chlorine-Mediated Electrochemical Ethylene Oxidation in Saline Water

Chlorine as a redox mediator allows for the selective oxidation of ethylene to 2-chloroethanol, which converts to ethylene oxide in alkaline aqueous electrolyte. This strategy utilizes abundant saline water as an electrolyte and source of oxygen atoms for functionalization. We present a mechanistic study of ethylene oxidation in saline water using cobalt oxide nanoparticle catalysts. Electrochemical kinetic analysis and in situ X-ray absorption spectroscopy suggest that the resting state of the catalyst and the rate-determining step differ for the chlorine evolution reaction in the presence and absence of ethylene. In 0.6 M NaCl pH 8 electrolyte, which resembles seawater, the average current density was ~60 mA/cm 2 with a Faradaic efficiency of ~41% toward ethylene functionalization. The use of synthetic and natural seawater achieved Faradaic efficiencies above 70%, while the partial current toward the product remained invariant. Further conversion of the initial product 2-chloroethanol into ethylene glycol was also demonstrated. Here, we present a broader vision of harnessing saline water in electrochemical functionalization of organic molecules and coproduction of hydrogen.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

How a Formate Dehydrogenase Responds to Oxygen: Unexpected O 2 Insensitivity of an Enzyme Harboring Tungstopterin, Selenocysteine, and [4Fe–4S] Clusters

The reversible two-electron interconversion of formate and CO 2 is catalyzed by both nonmetallo- and metallo-formate dehydrogenases (FDHs). The latter group comprises molybdenum- or tungsten-containing enzymes with the metal coordinated by two equivalents of a pyranopterin cofactor, a cysteinyl or selenocysteinyl (Sec) ligand supplied by the polypeptide, and a catalytically essential terminal sulfido ligand. In addition, these biocatalysts incorporate one or more [4Fe–4S] clusters for facilitating long-distance electron transfer. However, an interesting dichotomy arises when attempting to understand how the metallo-FDHs react with O 2 . Whereas existing scholarship portrays these enzymes as being unable to perform in air due to extreme O 2 lability of their metal centers, studies dating as far back as the 1930s emphasize that some of these systems exhibit formate oxidase (FOX) activity, coupling formate oxidation to O 2 reduction. Therefore, to reconcile these conflicting views, we explored context-dependent functional linkages between metallo-FDHs and their cognate electron acceptors within the same organism vis-à-vis catalysis under atmospheric O 2 . Here, we report the discovery and characterization of an O 2 -insensitive FDH2 from the sulfate-reducing bacterium Desulfovibrio vulgaris Hildenborough (DvH) that ligates tungsten, Sec, and four [4Fe–4S] clusters. By advancing a robust expression platform for its recombinant production, we eliminate both the requirement of nitrate or azide during purification and reductive activation with thiols and/or formate prior to catalysis. Because the distinctive spectral signatures of formate-reduced DvH-FDH2 remain invariant under anaerobic and aerobic conditions, we benchmarked the enzyme activity in air, identifying CO 2 as the catalytic product. Full reaction progress curve analysis discloses a high catalytic efficiency when probed with a high-potential artificial electron acceptor. Furthermore, we show that DvH-FDH2 enables near-stoichiometric hydrogen peroxide production without superoxide release to achieve O 2 insensitivity. Notably, simultaneous electron transfer to cytochrome c and O 2 reveals that metal-based electron bifurcation is operational in this system. Taken together, our work proves the co-occurrence of redox bifurcated FDH and FOX activities within a metalloenzyme scaffold. These findings set the stage for uncovering previously unknown O 2 -insensitive flavin-based electron bifurcation mechanisms, as well as for developing authentic formate/air biofuel cells, engineering O 2 -stable FDHs and biohybrid metallocatalysts, and discerning formate bioenergetics of gut microbiota.

13C NMR↗
DOE OSTI2026

Quantifying the Impact of Metal Population Distribution in MFI-Supported Mo Catalysts for Methane Dehydroaromatization

Precise evaluation of intrinsic kinetic behavior in Mo/MFI catalysts for methane dehydroaromatization (MDA) is confounded by variations in Mo dispersion and speciation, which are influenced by metal loading and zeolite acidity. This work advances the utility of H 2 -temperature programmed reduction (H 2 -TPR) for characterizing Mo/MFI catalysts by enabling quantitative comparison of MoO x populations distinguished by reduction behavior and linked to initial catalytic performance. A comprehensive H 2 reduction pathway is established through systematic H 2 - TPR studies varying catalyst composition (1–10 wt % Mo, Si/Al = 15, 40, ∞), supplemented by UV-Raman spectroscopy, X-ray powder diffraction, N 2 physisorption, NH 3 -TPD, and advanced spectroscopic analysis (in situ XAS with principal component analysis/multivariate curve resolution—alternating least squares). Two low-temperature H 2 -TPR regions capture distinct Mo populations undergoing initial Mo(VI)→Mo(IV) reduction: a lower-temperature population (Mo-RI) associated primarily with highly dispersed, anchored MoO x species expected to predominantly reside within MFI channels, and a higher-temperature population (Mo-RII) corresponding to a broader set of MoO x species that becomes increasingly bulk-like/extrazeolitic at higher Mo loading. Quantification of these populations provides practical, kinetically relevant descriptors for comparing initial MDA rates across catalysts with varying Mo loading, Si/Al ratio, and MoO x heterogeneity. Application of lower-temperature Mo-RI estimates to kinetic measurements reveals a minimum threshold of ∼0.12 × 10 –3 mol Mo-RI/g cat , above which initial forward benzene rates normalized to this population converge despite differences in metal loading and zeolite Brønsted acidity. This threshold coincides with a transition toward a common C 2 -mediated benzene-forming regime, as indicated by approach-to-equilibrium analysis of methane-to-ethane, ethane-to-ethylene, and ethylene-to-benzene reaction steps. Above this threshold, normalized initial benzene rates are nearly invariant with increasing Mo-RII/Mo-RI population ratio, indicating that excess Mo-RII populations, including bulk-like/extrazeolitic MoO x domains present at higher loading, do not measurably suppress benzene formation associated with anchored and mostly channel-confined Mo population under the initial-rate conditions examined.

Mo/MFI↗
DOE OSTI2025

Experimental and Computational Evaluation of Nicotinamide Cofactor Biomimetics

Oxidoreductase enzymes are widely used biocatalysts due to their high enantioselectivity and broad substrate compatibility in useful transformations. Many oxidoreductases require nicotinamide cofactors (i.e., NAD(P)H). To replace this costly natural cofactor, synthetic nicotinamide cofactor biomimetics (NCBs) offer different shapes, binding affinities, and reducing potentials that exceed the capabilities of wild-type NAD(P)H. However, the ill-defined structure–activity relationships (SARs) of various NCBs slow rationally guided innovation, such as customized reducing potentials. Here, we dissect two essential elements of NCB design, holding the nicotinamide invariant. First, the linker length between the nicotinamide and an unconjugated aromatic ring uncovered unexpected benefits to redox activity for two or three carbon linkers. Second, substitution on this unconjugated aryl group (Ring 2) might not be expected to affect activity. However, SAR trends demonstrate substantial benefits to reductive potential conferred by electron-donating functionalities on Ring 2. Furthermore, catalysis by two enzymes demonstrates enzyme-dependent tolerance or sensitivity to the NCB structures. Density functional theory (DFT) and computational modeling provide a theoretical framework to understand and build upon these observations. Ring 2 reaches up to the nicotinamide to stabilize its positive charge after oxidation through π–π stacking and charge transfer. Thus, the systematic examination of NCB’s stability, electrochemical redox potentials, and kinetics uncovers trends for the improved design of NCBs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Distinguishing Models for Mixed Halide Lead Perovskite Photosegregation via Terminal Halide Stoichiometry

Photoinduced halide segregation in mixed halide hybrid perovskites [i.e., APb(I 1–x Br x ) 3 ] represents an intrinsic instability that impedes their commercialization. Resolving this issue requires developing a microscopic understanding of the phenomenon. Key to this is distinguishing existing models of halide photosegregation by comparing their corresponding predictions to experiment. In this work, we test the temperature dependency of predicted perovskite terminal stoichiometries (x terminal ), following photosegregation, in single and double cation mixed halide perovskites. Our results show a general temperature invariance of x terminal . This largely supports the idea that band gap difference between parent and halide-segregated perovskite phases drives photosegregation. Beyond this, careful examination of temperature- and halide composition-dependent emission energies suggests an alternative explanation for the enhanced photostability of mixed cation/mixed halide perovskites, linked to band gap energy differences between parent and phase-segregated perovskite phases. Together, our results make inroads in clarifying the microscopic origin of mixed halide perovskite photosegregation and pave the way toward approaches that stabilize their use in applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Extrinsic Nonlinear Kerr Rotation in Topological Materials under a Magnetic Field

Topological properties in quantum materials are often governed by symmetry and tuned by crystal structure and external fields, and hence, symmetry-sensitive nonlinear optical measurements in a magnetic field are a valuable probe. Here, we report nonlinear magneto-optical second harmonic generation (SHG) studies of nonmagnetic topological materials including bilayer WTe 2 , monolayer WSe 2 , and bulk TaAs. The polarization-resolved patterns of optical SHG under a magnetic field show nonlinear Kerr rotation in these time-reversal symmetric materials. For materials with 3-fold rotational symmetric lattice structure, the SHG polarization pattern rotates just slightly in a magnetic field, whereas in those with mirror or 2-fold rotational symmetry, the SHG polarization pattern rotates greatly and distorts. Further, these different magneto-SHG characters can be understood by considering the superposition of the magnetic field-induced time-noninvariant nonlinear optical tensor and the crystal-structure-based time-invariant counterpart. The situation is further clarified by scrutinizing the Faraday rotation, whose subtle interplay with crystal symmetry accounts for the diverse behavior of the extrinsic nonlinear Kerr rotation in different materials. Our work illustrates the application of magneto-SHG techniques to directly probe nontrivial topological properties, and underlines the importance of minimizing extrinsic nonlinear Kerr rotation in polarization-resolved magneto-optical studies.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Symmetry-Protected Moiré Band Engineering and Enhanced Electron–Phonon Coupling in Xe/Bi 2 Se 3 Superlattices: Path to Topological Superconductivity

Observation of superconductivity, magnetism, and correlated insulating phases driven by the moiré potential in twisted graphene bilayer has opened the exciting new field of “twistronics”. Even richer physics is expected if moiré superlattice could be generated on topological insulators; however, until now, experimental studies have been scarce. Here, we demonstrate topological moirés generated by adsorbing a monolayer of noble gas on a topological insulator. By angle-resolved photoemission spectroscopy, we show that the moiré potential replicates the topological surface state and affects it in a way fundamentally different from the trivial states. Replicated Dirac cones generally avoid crossings, except at the time-reversal invariant momenta that remain gapless. This creates van Hove singularities at the moiré Brillouin zone corners, providing the mechanism of enhancing correlations. Indeed, we observe a strong enhancement of the electron–phonon coupling strength that, if properly tuned, might lead to topological superconductivity and Majorana Fermions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Magnetism in Mixed Valence, Defect, Cubic Perovskites: BaIn 1– x Fe x O 2.5+δ , x = 0.25, 0.50, and 0.75. Local and Average Structures

The series BaIn 1–x Fe x O 2.5+δ , x = 0.25, 0.50, and 0.75, has been prepared under air-fired and argon-fired conditions and studied using X-ray diffraction, d.c. and a.c. susceptibility, Mössbauer spectroscopy, neutron diffraction, X-ray near edge absorption spectroscopy (XANES), and X-ray pair distribution (PDF) methods. While Ba 2 In 2 O 5 (BaInO 2.5 ) crystallizes in an ordered brownmillerite structure, Ibm2, and Ba 2 Fe 2 O 5 (BaFeO 2.5 ) crystallizes in a complex monoclinic structure, P2 1 /c, showing seven Fe 3+ sites with tetrahedral, square planar, and octahedral environments, all phases studied here crystallize in the cubic perovskite structure, Pm$\bar{}3}$m, with long-range disorder on the small cation and oxygen sites. 57 Fe Mössbauer studies indicate a mixed valency, Fe 4+ /Fe 3+ , for both the air-fired and argon-fired samples. The increased Fe 3+ content for the argon-fired samples is reflected in increased cubic cell constants and in the increased Mössbauer fraction. It appears that the Pm$\bar{}3}$m phases are only metastable when fired in argon. From a slightly modified percolation theory for a primitive cubic lattice (taking into account the presence of random O atom vacancies), long-range spin order is permitted for the x = 0.50 and 0.75 phases. Instead, the d.c. susceptibility shows only zero-field-cooled (ZFC) and field-cooled (FC) divergences at ~6 K [5 K] for x = 0.50 and at ~22 K [21 K] for x = 0.75, with values for the argon-fired samples in [ ]. Neutron diffraction data for the air-fired samples confirm the absence of long-range magnetic order at any studied temperature. For the air-fired x = 0.50, a.c. susceptibility data show a frequency-dependent χ'(max) and spin glass behavior, while for x = 0.75, χ'(max) is invariant with frequency, ruling out either a spin glass or a superparamagnetic ground state. These behaviors are discussed in terms of competing Fe 3+ –Fe 3+ antiferromagnetic exchange and ferromagnetic Fe 3+ –Fe 4+ exchange. The PDF and 57 Fe Mössbauer data indicate a local structure at short interatomic distances, which deviates strongly from the average Pm$\bar{}3}$m model. Fe Mössbauer, PDF, and XANES data show a systematic dependence on x and indicate that the Fe 3+ sites are largely fourfold-coordinated and Fe 4+ sites are fivefold- or sixfold-coordinated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Cyclometallated Co(III) Complexes with Lowest-Energy Charge Transfer Excited States Accessible with Visible Light

The Co(III) complexes, cis-[Co(ppy) 2 (L)]PF 6 , where ppy = 2-phenylpyridine and L = bpy (2,2′-bipyridine; 1), phen (1,10-phenanthroline; 2), and DAP (1,12-diazaperylene; 3), are reported and their photophysical properties were investigated to evaluate their potential as sensitizers for applications that include solar energy conversion schemes and photoredox catalysis. Calculations show that cyclometallation in the cis-[Co(ppy) 2 (L)]PF 6 series affords strong Co(dπ)/ppy(π) orbital interactions that result in a Co/ppy(π*) highest occupied molecular orbital (HOMO) and a lowest unoccupied molecular orbital (LUMO) localized on the diimine ligand, L(π*). Complexes 1–3 exhibit relatively invariant oxidation potentials, whereas the reduction event is dependent on the identity of the diimine ligand, L, consistent with the theoretical predictions. For 3 a broad Co/ppy(π*) → L(π*) metal/ligand-to-ligand charge transfer (ML-LCT) absorption band is observed in CH 3 CN with a maxima at 507 nm, extending beyond 600 nm. Upon excitation of the 1ML-LCT transition, transient absorption features consistent with the population of a 3 ML-LCT excited state with lifetimes, τ, of 3.0 ps, 4.6 and 42 ps for 1, 2 and 3 in CH 3 CN respectively are observed. Upon irradiation with 505 nm, 3 is able to reduce methyl viologen (MV 2+ ), an electron acceptor commonly in photocatalytic schemes. Here, to our knowledge, 3 represents the first heteroleptic molecular Co(III) complex that combines cyclometallation with a diimine ligand with lowest-lying metal-to-ligand charge transfer excited states able to undergo photoinduced charge transfer with low-energy green light. As such, the structural design of 3 represents an important step toward d6 photosensitizers based on earth abundant metals.

Energy↗
DOE OSTI2025

Reversible Interfacial Hydride Transfer as a Complementary Tool To Measure Molecular Hydricity

Hydride transfer is an essential elementary reaction across the chemical value chain, but there are limited methods available for quantifying thermodynamic hydricity (ΔG H- ), particularly among main group reagents. Herein, we exploit facile H 2 activation and reversible hydride transfer from a metal surface to a molecular reagent, the net hydrogen reduction reaction (HRR), to develop a potentiometric method for quantifying ΔG H- of main group reagents recalcitrant to conventional methods. HRR potentiometry is first validated with a benzimidazole-based hydride donor and then applied to uncover the impact of the reaction environment on hydricity. For a benzimidazole-based hydride donor, HRR equilibrium potentials are roughly invariant across solvents, indicating that the solvent dependence of its hydricity largely reflects the differential solvation of H - across media. For formate, HRR potentials and corresponding hydricities depend strongly on water content. Additionally, for borohydrides, HRR potentiometry reveals that effective hydricity values are strongly influenced by Lewis acid-base adduct formation with the hydride acceptor but are minimally influenced by the countercation. Together with these studies, the advantages, limitations, and practical considerations of the HRR potentiometry method are discussed, highlighting the power of this methodology as a complementary tool to measure molecular hydricity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Driving Force Dependent Photoinduced Charge Transfer Dynamics in Polymer-Wrapped Semiconducting Single-Walled Carbon Nanotubes

Here, we investigate the thermodynamic driving-force dependences of photoinduced charge separation (CS) and subsequent charge transfer dynamics in single-walled carbon nanotube (SWNT)–perylenediimide (PDI) donor–acceptor (D–A) superstructures. Pump–probe spectroscopy reveals that [SWNT (•+)n ]-(PDI –• ) n CS states form on an ∼100 fs time scale following photoexcitation; these dynamics are invariant across an ∼400 mV driving force range, indicating that SWNT hole polaron formation time scales are determined by nanotube lattice and solvent relaxation. These CS states feature SWNT hole polarons adjacent to (geminate) and nearby (nongeminate) PDI radical anions. Analysis of the free energy dependence for charge recombination (CR) of [SWNT •+ ] geminate -(PDI –• ) CS states highlights an ∼2 meV value for D–A electronic coupling (HAB) and ∼0.93 eV for the total reorganization energy (λ T ). A corresponding driving force dependence of the CR dynamics for [SWNT •+ ] nongeminate -(PDI –• ) CS states indicates a diminished H AB value (∼0.6 meV) and a larger λ T (∼1.1 eV), consistent with larger transfer distances. SWNT excitons that persist following photoinduced CS drive photooxidation of PDI –• components of [SWNT (•+)n ]-(PDI –• ) n CS states ( 1 SWNT* + PDI •– → PDI + SWNT •– ); this reaction manifests a significantly reduced λ T value (∼0.67 eV) as the initially prepared SWNT reduced state bears the character of a conduction band injected electron ([SWNT •– ] CB ). This reaction thus gives rise to relaxed, nongeminate SWNT electron and hole polarons on the same nanotube; these polarons react on a 10 2 ps time scale independent of the electronic structure of these SWNT-PDI superstructures.

anions↗
DOE OSTI2020

Development of Future Heatwaves for Different Hazard Thresholds

In 2018 and 2019, heatwaves set all-time temperature records around the world and caused adverse effects on human health, agriculture, natural ecosystems, and infrastructure. Often, severe impacts relate to the joint spatial and temporal extent of the heatwaves, but most research so far focuses either on spatial or temporal attributes of heatwaves. Furthermore, sensitivity of heatwaves characteristics to the choice of the heatwave thresholds in a warming climate are rarely discussed. Here, we analyze the largest spatiotemporal moderate heatwaves—that is, three-dimensional (space-time) clusters of hot days—in simulations of global climate models. We use three different hazard thresholds to define a hot day: fixed thresholds (time-invariant climatological thresholds), seasonally moving thresholds based on changes in the summer means, and fully moving thresholds (hot days defined relative to the future climatology). We find a substantial increase of spatiotemporally contiguous moderate heatwaves with global warming using fixed thresholds, whereas changes for the other two hazard thresholds are much less pronounced. In particular, no or very little changes in the overall magnitude, spatial extent, and duration are detected when heatwaves are defined relative to the future climatology using a temporally fully moving threshold. This suggests a dominant contribution of thermodynamic compared to dynamic effects in global climate model simulations. The similarity between seasonally moving and fully moving thresholds indicates that seasonal mean warming alone can explain large parts of the warming of extremes. The strong sensitivity of simulated future heatwaves to hazard thresholds should be considered in the projections of potential future heat-related impacts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

River Dynamics Control Transit Time Distributions and Biogeochemical Reactions in a Dam-Regulated River Corridor

Residence Time Distributions (RTDs) exert an important control on biogeochemical translation in the hydrological systems. Previous tracer studies have revealed that RTDs often followed time-invariant exponential, lognormal, or heavy-tailed RTDs that have power-law behavior for long tails in headwater or low-order streams. However, there is a recognition that RTDs can be more complicated and time-variable in response to dynamic hydrological forcing. Here, we use particle tracking to estimate RTDs along the Hanford Reach of the Columbia River and to quantify the influences of river stage fluctuations. Particle tracking is conducted using the velocity field from high-resolution 3D groundwater flow simulations. The effects of hydrological forcing on the residence time distribution were evaluated by varying river flow boundary conditions and releasing particles in different time windows. Our results revealed that dynamic stage fluctuations created rapidly changing losing-gaining conditions in the river and led to highly transient RTDs, which contributes to multiple modes of RTDs. Dam-induced high-frequency (sub-daily) flow variation contributes to the short-time (sub-daily) component of the RTDs. Deviation of the reactant consumption under the single-mode assumption compared to the multimodal RTDs is relatively small (~5%) and appears when the Damköhler number was close to one.More specifically, high-frequency (daily) dam-induced stage variations have higher impacts on biogeochemical reactions with faster reaction rate, and can potentially increase the reactant consumption by 26.68%. Based on these findings, we suggest that current river basin models could be improved by including bank storage and more complex RTDs influenced by both short and long term river stage fluctuations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Modification of the Loss Cone for Energetic Particles in the Earth's Inner Magnetosphere

The behavior of energetic particles in the Earth's dipole and stretched magnetic-field models is explored with a focus on the shift of the atmospheric loss cone away from the magnetic-field direction and on non-adiabatic behavior occurring when the particle gyroradius becomes comparable with the gradient scale length of the local magnetic field. It is shown that the equatorial loss cone is aberrated away from the magnetic-field direction (pitch angle of 0°) by the perpendicular drift of a charged particle (which is referred to as the “Mozer transform” described in Mozer (1966, https://doi.org/10.1029/JZ071i011p02701). It is found that the Mozer coordinate transformation in pitch-angle/gyrophase-angle space better organizes the behavior of bounce times, mirror altitudes, drift speeds, and first adiabatic invariant. It also describes the loss-cone shift accurately for a certain range of values of the adiabaticity parameter ϵ, defined by the ratio of the particle gyroradius to the local radius of curvature of the magnetic field. For particles with larger ϵ, the Mozer transform (evaluated with the standard gradient-curvature drift) breaks down and the “central trajectory” theory can be used to calculate the angular shift of the loss cone, which now includes an Earthward shift. When ϵ becomes relatively large, stochastic field line curvature (FLC) scattering occurs: we show the intimate connection between the strength of FLC scattering and the loss-cone shift. Finally, by comparing ion orbits in the dipole and Tsyganenko (Ts89 and Ts04) magnetic fields, we conclude that the loss cone of ring-current ions is significantly modified during geomagnetically active times.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The Impact of Crop Rotation and Spatially Varying Crop Parameters in the E3SM Land Model (ELMv2)

Abstract Earth System Models (ESMs) are increasingly representing agriculture due to its impact on biogeochemical cycles, local and regional climate, and fundamental importance for human society. Realistic large scale simulations may require spatially varying crop parameters that capture crop growth at various scales and among different cultivars, as well as common crop management practices, but their importance is uncertain, and they are often not represented in ESMs. In this study, we examine the impact of using constant versus spatially varying crop parameters using a novel, realistic crop rotation scenario in the Energy Exascale Earth System Model (E3SM) Land Model version 2 (ELMv2). We implemented crop rotation by using ELMv2's dynamic land unit capability, and then calibrated and validated the model against observations collected at three AmeriFlux sites in the US Midwest with corn soybean rotation. The calibrated model closely captured the magnitude and observed seasonality of carbon and energy fluxes across crops and sites. We performed regional simulations for the US Midwest using the calibrated model and found that spatially varying only a few crop parameters across the region, as opposed to using constant parameters, had a large impact, with the carbon fluxes and energy fluxes both varying by up to 40%. These results imply that large scale ESM simulations using spatially invariant crop parameters may result in biased energy and carbon fluxes estimation from agricultural land, and underline the importance of improving human‐earth systems interactions in ESMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Evaluating the Effects of Precipitation and Evapotranspiration on Soil Moisture Variability Within CMIP5 Using SMAP and ERA5 Data

=The effects of precipitation (Pr) and evapotranspiration (ET) on surface soil moisture (SSM) play an essential role in the land-atmosphere system. Here we evaluate multimodel differences of these effects within the Coupled Model Intercomparison Project Phase 5 (CMIP5) compared to Soil Moisture Active Passive (SMAP) products and ECMWF Reanalysis v5 (ERA5) as references in a frequency domain. The variability of SSM, Pr, and ET within three frequency bands (1/7 ~ 1/30 days -1 , 1/30 ~ 1/90 days -1 , and 1/90 ~ 1/365 days -1 ) after normalization is quantified using Fourier transform. We analyze the impact of ET and Pr on SSM variability based on a transfer function assuming that these variables form a linear time-invariant (LTI) system. For the total effects of ET and Pr on SSM variability, the CMIP5 estimations are smaller than the reference data in the two higher frequency bands and are larger than the reference data in the lowest frequency band. Besides, the effects on SSM by Pr and ET are found to be different across the three frequency bands. In each frequency band, the variability of the factor that dominates SSM (i.e., Pr or ET) from CMIP5 is smaller than that from the references. This study identifies the spatiotemporal distribution of differences between CMIP5 models and references (SMAP and ERA5) in simulating ET and Pr effects on SSM within three frequency bands. This study provides insightful information on how soil moisture variability is affected by varying precipitation and evapotranspiration at different time scales within Earth System Models.

54 ENVIRONMENTAL SCIENCES↗