Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamical core”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Prediction of cccDNA dynamics in hepatitis B patients by a combination of serum surrogate markers

Quantification of intrahepatic covalently closed circular DNA (cccDNA) is a key for evaluating an elimination of hepatitis B virus (HBV) in infected patients. However, quantifying cccDNA requires invasive methods such as a liver biopsy, which makes it impractical to access the dynamics of cccDNA in patients. Although HBV RNA and HBV core-related antigens (HBcrAg) have been proposed as surrogate markers for evaluating cccDNA activity, they do not necessarily estimate the amount of cccDNA. Here, we employed a recently developed multiscale mathematical model describing intra- and intercellular viral propagation and applied it in HBV-infected patients under treatment. We developed a model that can predict intracellular HBV dynamics by use of extracellular viral markers, including HBsAg, HBV DNA, and HBcrAg in peripheral blood. Importantly, the model prediction of the amount of cccDNA in patients over time was confirmed to be well correlated with the data for quantified cccDNA by paired liver biopsy. Thus, our method combining classic and emerging surrogate markers enables us to predict the decay dynamics of cccDNA in patients undergoing treatment.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Heterogeneity and Hysteresis in the Polymer Collapse of Single Core–Shell Stimuli-Responsive Plasmonic Nanohybrids

Broad application of polymeric stimuli-responsive smart nanohybrids requires understanding the mechanisms governing active control. Ensemble techniques have identified inhomogeneous polymer collapse in microgels that potentially arise from heterogeneous interchain interactions and differences in core size. A single-particle examination would establish the influence of core size and internal polymer network heterogeneity on local interactions that contribute to the observed inhomogeneous polymer collapse dynamics of nanohybrids. Using single-particle dark-field spectroscopy, we investigated the complex polymer collapse profiles of core–shell plasmonic nanohybrids comprising thermoresponsive poly(N-isopropylacrylamide) (pNIPAM)-encapsulated gold nanorods (AuNRs). We report that the polymer collapse behavior was independent of the core size. For thinner polymer shells, we observed hysteresis in the collapse of AuNR@pNIPAMs, likely related to local pNIPAM aggregation due to interchain hydrogen bonding. For thicker polymer shells, we observed a broad polymer collapse distribution that we attributed to a two-step phase transition that arises from a polymer network density gradient. Furthermore, our single-particle approach relates the internal heterogeneity of the polymer network of nanohybrids to the mechanisms underlying heterogeneous phase transitions that traditional, ensemble-averaged approaches are unable to discern.

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

Effects of iron on the lattice thermal conductivity of Earth’s deep mantle and implications for mantle dynamics

We discuss how iron may critically influence the physical properties and thermochemical structures of Earth’s lower mantle. Its effects on thermal conductivity, with possible consequences on heat transfer and mantle dynamics, however, remain largely unknown. We measured the lattice thermal conductivity of lower-mantle ferropericlase to 120 GPa using the ultrafast optical pump-probe technique in a diamond anvil cell. The thermal conductivity of ferropericlase with 56% iron significantly drops by a factor of 1.8 across the spin transition around 53 GPa, while that with 8–10% iron increases monotonically with pressure, causing an enhanced iron substitution effect in the low-spin state. Combined with bridgmanite data, modeling of our results provides a self-consistent radial profile of lower-mantle thermal conductivity, which is dominated by pressure, temperature, and iron effects, and shows a twofold increase from top to bottom of the lower mantle. Such increase in thermal conductivity may delay the cooling of the core, while its decrease with iron content may enhance the dynamics of large low shear-wave velocity provinces. Our findings further show that, if hot and strongly enriched in iron, the seismic ultralow velocity zones have exceptionally low conductivity, thus delaying their cooling.

58 GEOSCIENCES↗
DOE OSTI2016

Reducing I/O variability using dynamic I/O path characterization in petascale storage systems

In petascale systems with a million CPU cores, scalable and consistent I/O performance is becoming increasingly difficult to sustain mainly because of I/O variability. Furthermore, the I/O variability is caused by concurrently running processes/jobs competing for I/O or a RAID rebuild when a disk drive fails. We present a mechanism that stripes across a selected subset of I/O nodes with the lightest workload at runtime to achieve the highest I/O bandwidth available in the system. In this paper, we propose a probing mechanism to enable application-level dynamic file striping to mitigate I/O variability. We also implement the proposed mechanism in the high-level I/O library that enables memory-to-file data layout transformation and allows transparent file partitioning using subfiling. Subfiling is a technique that partitions data into a set of files of smaller size and manages file access to them, making data to be treated as a single, normal file to users. Here, we demonstrate that our bandwidth probing mechanism can successfully identify temporally slower I/O nodes without noticeable runtime overhead. Experimental results on NERSC’s systems also show that our approach isolates I/O variability effectively on shared systems and improves overall collective I/O performance with less variation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Insights on dramatic radial fluctuations in track formation by energetic ions

We discuss the insights on the unexpected dramatic radial variations in the ion tracks formed by energetic ion (2.3 GeV 208 Pb) irradiation at a constant electronic energy-loss (~42 keV/nm) in pyrochlore structured Gd 2 TiZrO 7 . Though previous studies have shown track formation and average track diameter measurements, this work brings further clarity on why quantitative analysis of ion track formation in Gd 2 Ti x Zr (1-x) O 7 systems can be more complicated than the currently accepted behavior for ion tracks. The ion track profile is usually considered to be diametrically uniform at constant values of the electronic energy-loss. This study shows the diameter variations to be as large as ~40% within an extremely short incremental track length of ~20 nm. Our molecular dynamics simulations show that these fluctuations in diameter of amorphous core and overall track diameter are attributed to (i) the stochastic nature of inelastic energy loss along the track and (ii) the random substitution of Ti atoms by Zr atoms on the B-site in the pyrochlore lattice. Furthermore, the partial substitution of Ti by Zr increases the favorability of the defect-fluorite structure formation over amorphous phase stochastically, by introducing localized inhomogeneity in atomic structure, density and strain.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Time-resolved atomic-resolution Brownian tomography of single nanocrystals reveals size-dependent dynamics

Atomic-resolution structure identification of nanocrystals by graphene liquid cell electron microscopy (GLC-EM) has revealed that small, solubilized platinum nanocrystals consist of an ordered crystalline core surrounded by mobile surface atoms, which dissociate during oxidative etching, resulting in distinct temporal structural states. Requirements imposed by the 3D reconstruction algorithm limit the number of structural states that can be resolved. We introduce a regularized 3D reconstruction algorithm that exploits the redundancy inherent in the experimental data, allowing us to improve the time resolution. Our developments provide a comprehensive molecular picture at unprecedented spatial and temporal resolution of the nonlinear, linear, and fluctuating dynamic phenomena that single nanocrystals undergo during the GLC-EM experiment. We determined atomic structures of 66 temporal structural states, extracted from 15 time trajectories of individual nanocrystals. Large (478 to 698 atoms) and small (<300 atoms) nanocrystals show etching that preserves a stable core, whereas mid-sized (351 to 571 atoms) nanocrystals present dynamics that change the coordination of the core.

Meana-Pañeda, Rubén↗
DOE OSTI2024

Dynamical friction in dark matter superfluids: The evolution of black hole binaries

Abstract The theory of superfluid dark matter is characterized by self-interacting sub-eV particles that thermalize and condense to form a superfluid core in galaxies. Massive black holes at the center of galaxies, however, modify the dark matter distribution and result in a density enhancement in their vicinity known as dark matter spikes. The presence of these spikes affects the evolution of binary systems by modifying their gravitational wave emission and inducing dynamical friction effects on the orbiting bodies. In this work, we assess the role of dynamical friction for bodies moving through a superfluid core enhanced by a central massive black hole. As a first step, we compute the dynamical friction force experienced by bodies moving in a circular orbit. Then, we estimate the gravitational wave dephasing of the binary, showing that the effect of the superfluid drag force is beyond the reach of space-based experiments like LISA, contrarily to collisionless dark matter, therefore providing an opportunity to distinguish these dark matter models.

Astronomy & Astrophysics↗
DOE OSTI2016

Column-by-column observation of dislocation motion in CdTe: Dynamic scanning transmission electron microscopy

The dynamics of partial dislocations in CdTe have been observed at the atomic scale using aberration-corrected scanning transmission electron microscopy (STEM), allowing the mobility of different dislocations to be directly compared: Cd-core Shockley partial dislocations are more mobile than Te-core partials, and dislocation cores with unpaired columns have higher mobility than those without unpaired columns. The dynamic imaging also provides insight into the process by which the dislocations glide. Dislocations with dangling bonds on unpaired columns are found to be more mobile because the dangling bonds mediate the bond exchanges required for the dislocations to move. Furthermore, a screw dislocation has been resolved to dissociate into a Shockley partial-dislocation pair along two different directions, revealing a way for the screw dislocation to glide in the material. The results show that dynamic STEM imaging has the potential to uncover the details of dislocation motion not easily accessible by other means. This research was sponsored by the U.S. DOE, Office of Energy Efficiency and Renewable Energy, Foundational Program to Advance Cell Efficiency (F-PACE, DE-FOA-0000492), (CL, YLW, NP, YFY, SJP), the Office of DOE-BES, Materials Science and Engineering Division (ARL), DOE Grant No. DE-FG02-09ER46554 (YYZ, STP), and a user project supported by ORNL's Center for Nanophase Materials Sciences (CNMS), which is also sponsored by DOE-BES. Supercomputer time was provided by the National Center for Supercomputing Applications which is supported by the DOE Office of Science under Contract No. DE-AC02-05CH11231, and the Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by National Science Foundation Grant No. ACI-1053575. Chen Li [No. 656378] and Timothy Pennycook [No.655760] are currently funded by the European Union's Marie Sklodowska-Curie Grants.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Are superluminous supernovae and long GRBs the products of dynamical processes in young dense star clusters?

Superluminous supernovae (SLSNe) occur almost exclusively in small galaxies (Small/Large Magellanic Cloud (SMC/LMC)-like or smaller), and the few SLSNe observed in larger star-forming galaxies always occur close to the nuclei of their hosts. Another type of peculiar and highly energetic supernovae are the broad-line Type Ic SNe (SN Ic-BL) that are associated with long-duration gamma-ray bursts (LGRBs). Also these have a strong preference for occurring in small (SMC/LMC-like or smaller) star-forming galaxies, and in these galaxies LGRBs always occur in the brightest spots. Studies of nearby star-forming galaxies that are similar to the hosts of LGRBs show that these brightest spots are giant H II regions produced by massive dense young star clusters with many hundreds of O- and Wolf-Rayet-type stars. Such dense young clusters are also found in abundance within a few hundred parsecs from the nucleus of larger galaxies like our own. We argue that the SLSNe and the SNe Ic-BL/LGRBs are exclusive products of two types of dynamical interactions in dense young star clusters. In our model the high angular momentum of the collapsing stellar cores required for the engines of an SN Ic-BL results from the post-main-sequence mergers of dynamically produced cluster binaries with almost equal-mass components. The merger produces a critically rotating single helium star with sufficient angular momentum to produce an LGRB; the observed 'metal aversion' of LGRBs is a natural consequence of the model. We argue that, on the other hand, SLSNe could be the products of runaway multiple collisions in dense clusters, and we present (and quantize) plausible scenarios of how the different types of SLSNe can be produced.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Protostar l1455 IRS1: a rotating disk connecting to a filamentary network

We conducted IRAM-30 m C{sup 18}O (2–1) and SMA 1.3 mm continuum {sup 12}CO (2–1) and C{sup 18}O (2–1) observations toward the Class 0/I protostar L1455 IRS1 in Perseus. The IRAM-30 m C{sup 18}O results show IRS1 in a dense 0.05 pc core with a mass of 0.54 M {sub ⊙}, connecting to a filamentary structure. Inside the dense core, compact components of 350 au and 1500 au are detected in the SMA 1.3 mm continuum and C{sup 18}O, with a velocity gradient in the latter one perpendicular to a bipolar outflow in {sup 12}CO, likely tracing a rotational motion. We measure a rotational velocity profile ∝r{sup −0.75} that becomes shallower at a turning radius of ∼200 au, which is approximately the radius of the 1.3 mm continuum component. These results hint at the presence of a Keplerian disk with a radius <200 au around L1455 IRS1 with a protostellar mass of about 0.28 M {sub ⊙}. We derive a core rotation that is about one order of magnitude faster than expected. A significant velocity gradient along a filament toward IRS1 indicates that this filament is dynamically important, providing a gas reservoir and possibly responsible for the faster-than-average core rotation. Previous polarimetric observations show a magnetic field aligned with the outflow axis and perpendicular to the associated filament on a 0.1 pc scale, while on the inner 1000 au scale, the field becomes perpendicular to the outflow axis. This change in magnetic field orientations is consistent with our estimated increase in rotational energy from large to small scales that overcomes the magnetic field energy, wrapping the field lines and aligning them with the disk velocity gradient. These results are discussed in the context of the interplay between filament, magnetic field, and gas kinematics from large to small scales. Possible emerging trends are explored with a sample of 8 Class 0/I protostars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Multiscale simulations of electron and ion dynamics in self-irradiated silicon

The interaction of energetic ions with the electronic and ionic system of target materials is an interesting but challenging multiscale problem, and understanding of the early stages after impact of heavy, initially charged ions is particularly poor. At the same time, energy deposition during these early stages determines later formation of damage cascades. In this work, we address the multiscale character by combining real-time time-dependent density functional theory for electron dynamics with molecular dynamics simulations of damage cascades. Our first-principles simulations prove that core electrons affect electronic stopping and have an unexpected influence on the charge state of the projectile. We show that this effect is absent for light projectiles, but dominates the stopping physics for heavy projectiles. By parametrizing an inelastic energy loss friction term in the molecular dynamics simulations using our first-principles results, we also show a qualitative influence of electronic stopping physics on radiation-damage cascades.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Dust Formation and Growth in Supernovae Explosions

Dust is a multi-messenger probe able to explore the dynamics of the Supernova Engine. The yield and composition of grains formed in the supernova explosion are dependent on the initial explosion characteristics of the progenitor and on the explosion mechanism. Using dust as a prove, the dynamics of dust formation and growth is modeled in a diverse set of Core-Collapse Supernovae Explosions (CCSNe). The progenitor mass, explosion energy engine type, and composition are varied across 72 models. The explosion is evolved using a 1-D Lagrangian Hydro-dynamical code (hydrocodes) to roughly 3 years. A multigrain nucleation and dust growth model is used to post process the hydrocode results and model grain formation and growth. The results show an earlier onset of grain formation, a larger silicate grain abundance, and smaller grain sizes in high explosion energy models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Spatial access and resource limitations control carbon mineralization in soils

Core-scale soil carbon fluxes are ultimately regulated by pore-scale dynamics of substrate availability and microbial access. These are constrained by physicochemical and biochemical phenomena (e.g. spatial access and hydrologic connectivity, physical occlusion, adsorption-desorption with mineral surfaces, nutrient and resource limitations). We conducted an experiment to determine how spatial access and resource limitations influence core-scale SOM mineralization, and how these are regulated by antecedent moisture conditions. Intact soil cores were incubated at field-moist vs. drought conditions, after which they were saturated from above (to simulate precipitation) or below (to simulate groundwater recharge). Soluble C (acetate) and N (nitrate) forms were added to some cores during the rewetting process to alleviate potential nutrient limitations. Respiration showed evidence of C-limitation, as C amendments increased the cumulative CO 2 evolved. Drought and rewetting increased soil respiration, and these soils also exhibited increased complex aromatic molecules in porewater. This newly available substrate appeared to alleviate nutrient limitations on respiration, because there were no further respiration increases with subsequent C and N amendments. We expected that respiration would be influenced by wetting direction, as simulated precipitation would mobilize C from the surface. However, as a main effect, this response was seen only in the C-amended soils, indicating that surface-C may not have been bioavailable. At the pore scale (porewater samples), compared to the baseline soils, drought and the C, N amendments caused a net loss of identified molecules when the soils were rewet from below, whereas wetting from above caused a net increase in identified molecules. This indicates that as soils were wet from below, the fresh inputs simulated the C-and N-limited microbial populations present deeper in the soil profile. Furthermore, our experiment highlights the complex and interactive role of antecedent moisture conditions, wetting direction, and resource limitations in driving core-scale C fluxes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Data for: Spatial access and resource limitations control carbon mineralization in soils

This dataset contains data and code used for the paper "Spatial access and resource limitations control carbon mineralization in soils", https://doi.org/10.1016/j.soilbio.2021.108427. Core-scale soil carbon fluxes are ultimately regulated by pore-scale dynamics of substrate availability and microbial access. These are constrained by physicochemical and biochemical phenomena (e.g. spatial access and hydrologic connectivity, physical occlusion, adsorption-desorption with mineral surfaces, nutrient and resource limitations). We conducted an experiment to determine how spatial access and resource limitations influence core-scale water-soluble soil organic matter (SOM) mineralization, and how these are regulated by antecedent moisture conditions. Intact soil cores were incubated at field-moist vs. drought conditions, after which they were saturated from above (to simulate precipitation) or below (to simulate groundwater recharge). Soluble carbon (acetate) and nitrogen (nitrate) forms were added to some cores during the rewetting process to alleviate potential nutrient limitations. Soil respiration was measured during the incubation, after which pore water was extracted from the saturated soils and analyzed for water soluble organic carbon concentrations and characterization. Our results showed that carbon (C) amendments increased the cumulative carbon dioxide (CO2) evolved from the soil cores, suggesting that the soils were C-limited. Drought and rewetting increased soil respiration, and there was a greater abundance of complex aromatic molecules in pore waters sampled from these soils. This newly available substrate appeared to alleviate nutrient limitations on respiration, because there were no further respiration increases with subsequent C and N amendments. We had hypothesized that respiration would be influenced by wetting direction, as simulated precipitation would mobilize C from the surface. However, as a main effect, this response was seen only in the C-amended soils, indicating that surface-C may not have been bioavailable. At the pore scale (pore water samples), drought and the C, N amendments caused a net loss of identified molecules when the soils were rewet from below, whereas wetting from above caused a net increase in identified molecules, suggesting that fresh inputs stimulated the C-and N-limited microbial populations present deeper in the soil profile. Our experiment highlights the complex and interactive role of antecedent moisture conditions, wetting direction, and resource limitations in driving core-scale C fluxes.This dataset contains a compressed (.zip) archive of the data and R scripts used for this manuscript. The dataset includes files in .csv format, which can be accessed and processed using MS Excel or R. This archive can also be accessed on GitHub at https://github.com/kaizadp/TES_spatial_access_2021 (DOI: 10.5281/zenodo.5522938).

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2001

Longtime Durability of PMR-15 Matrix Polymer at 204, 260, 288, and 316 C

Isothermal weight loss studies at the Glenn (Lewis) Research Center were conducted at four temperatures (204, 260, 288, and 316 C) with specimens of varied geometric shapes to investigate the mechanisms involved in the thermal degradation of PMR-15. Both neat resin behavior and composite behaviors were studied. Two points of interest in these studies are the role(s) of oxygen in the mechanisms involved in the thermo-oxidative degradation of these composite materials and the dimensional changes that occur during their useable lifetime. Specimen dimensional changes and surface layer growth were measured and recorded. It was shown that physical and chemical changes take place as a function of time and location in PMR-15 neat resin and composites as aging takes place in air at elevated temperatures. These changes initiate at the outer surfaces of both materials and progress inward following the oxygen as it proceeds by diffusion into the central core of each material. Microstructural changes cause changes in density, material shrinkage (strains), glass transition temperature, dimension, dynamic shear modulus, and compression properties. These changes also occur slowly dividing the polymer material into two distinct parts: a visibly undamaged core section between two visibly damaged surface layers. The surface layer has a significant effect on compression properties of thinner specimens, but the visibly undamaged core material controls these properties for specimens having eight or more plies. It was demonstrated that there are three different mechanisms involved in the degradation of PMR-15 during aging at elevated temperatures. These are a weight gain, a small weight fraction bulk material weight loss, and a large mass fraction weight loss concentrated at the surface of the polymer or composite. At the higher temperatures (260 C and above), the surface loss predominates. Below 260 C, the surface loss and the bulk core loss become more equivalent. Between 175 and 260 C, the initial weight change is due to a weight gain mechanism with a visible lifetime that diminishes as the aging temperature increases.

Bowles, Kenneth J.↗
DOE OSTI2016

Optical transients powered by magnetars: Dynamics, LIGHT curves, and transition to the nebular phase

Millisecond magnetars can be formed via several channels: core collapse of massive stars, accretion-induced collapse of white dwarfs (WDs), double WD mergers, double neutron star (NS) mergers, and WD-NS mergers. Because the mass of ejecta from these channels could be quite different, their light curves are also expected to be diverse. We evaluate the dynamic evolution of optical transients powered by millisecond magnetars. We find that the magnetar with a short spin-down timescale converts its rotational energy mostly into the kinetic energy of the transient, while the energy of a magnetar with a long spin-down timescale goes into radiation of the transient. This leads us to speculate that hypernovae could be powered by magnetars with short spin-down timescales. At late times the optical transients will gradually evolve into a nebular phase because of the photospheric recession. We treat the photosphere and nebula separately because their radiation mechanisms are different. In some cases the ejecta could be light enough that the magnetar can accelerate it to a relativistic speed. It is well known that the peak luminosity of a supernova (SN) occurs when the luminosity is equal to the instantaneous energy input rate, as shown by Arnett. We show that photospheric recession and relativistic motion can modify this law. The photospheric recession always leads to a delay of the peak time t{sub pk} relative to the time t{sub ×} at which the SN luminosity equals the instantaneous energy input rate. Relativistic motion, however, may change this result significantly.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Observing crystal nucleation in four dimensions using atomic electron tomography

Nucleation plays a critical role in many physical and biological phenomena that range from crystallization, melting and evaporation to the formation of clouds and the initiation of neurodegenerative diseases. However, nucleation is a challenging process to study in experiments, especially in the early stage when several atoms or molecules start to form a new phase from a parent phase. A number of experimental and computational methods have been used to investigate nucleation processes, but experimental determination of the three-dimensional atomic structure and the dynamics of early stage nuclei has been unachievable. Here we use atomic electron tomography to study early stage nucleation in four dimensions (4D: that is, including time) at atomic resolution. Using FePt nanoparticles as a model system, we find that early stage nuclei are irregularly shaped, each has a core of one to a few atoms with the maximum order parameter, and the order parameter gradient points from the core to the boundary of the nucleus. We capture the structure and dynamics of the same nuclei undergoing growth, fluctuation, dissolution, merging and/or division, which are regulated by the order parameter distribution and its gradient. These experimental observations are corroborated by molecular dynamics simulations of heterogeneous and homogeneous nucleation in liquid–solid phase transitions of Pt. Our experimental and molecular dynamics results differ from classical nucleation theory, indicating that a theory beyond this is needed to describe early stage nucleation at the atomic scale. Looking forward, we anticipate that the reported approach will open the door to the study of many fundamental problems in materials science, nanoscience, condensed matter physics and chemistry, such as phase transition, atomic diffusion, grain boundary dynamics, interface motion, defect dynamics and surface reconstruction with 4D atomic resolution.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Water foraging with dynamic roots in E3SM; The role of roots in terrestrial ecosystem memory on intermediate timescales (Final Technical Report)

Terrestrial ecosystems can show sustained responses to stress events that may last for years after the initial perturbation. This phenomenon is referred to as legacy or memory and it emerges from a set of ecosystem processes that are poorly captured by Earth System and Land Surface Models. Consequently, these models struggle to predict the impact that extreme climate events have on surface energy, carbon and hydrological exchange once a stressor such as drought has relaxed or in response to repeated exposure to stress. In this project, we set out to test how the addition of dynamic root profiles in land surface models impact the capacity for Earth System Models, such as the Department of Energy’s E3SM, to capture realistic legacy effects. Observations have shown that vegetation shifts their root profiles during stress events to forage for water and these altered root profiles may take years to relax back to the initial state. We hypothesized that the alteration of the belowground root structures may be a source of terrestrial ecosystem legacy that is missing from models. To test this idea we undertook three core activities. (1) We developed a global-scale analysis of the impact of dynamic roots on ecosystem legacy building from a recently developed dynamic root module. (2) We developed new root dynamics in the DOE’s Energy Land Model (ELM) that allowed not only root profiles to shift their profile but also simultaneously alter carbon allocation to fine root pools. We then tested the impact of these new dynamics with intensive sensitivity analysis at four long-term AmeriFlux sites. (3) We undertook detailed isotopic analysis of tree rings from these AmeriFlux sites to assess ecophysiological and ecohydrological legacy to stress events that can be used to benchmark the sensitivity experiments with ELM. From these core activities, we report the following key findings. Firstly, on a global scale, the addition of dynamic roots led to chronically water-stressed ecosystems recovering faster to climate stress while wetter ecosystem showed enhanced legacy. This is because across ecosystems, stress events were almost universally associated with water shortages that led to the development of deeper root profiles. These deeper root profiles proved beneficial for recovery from drought stress. While the root dynamics did not universally improve the modeled representation of legacy it showed complex transient dynamics that emerge from the addition of dynamic roots. Secondly, the sensitivity analysis illustrated long term shifts in rooting depths away from the prescribed default profile suggesting that initializing of root profiles could benefit from spin-up simulations that converge on locally optimized root profiles. In addition, by enabling dynamic allocation some of the sites predicted unrealistically low allocation to roots (and high allocation to leaves). While these changes did not dramatically alter modeled gross primary productivity, they illustrated that without more sophisticated root processes in models, there is little penalty to dramatically disinvest in roots. Thirdly, the isotopic analysis of tree rings showed highly distinct legacy responses across species and sites. For example, T. canadensis showed reduced transpiration the year after stress events illustrated by sustained elevated $\delta^{18}$O. In contrast, A saccharum displayed elevated $\delta^{13}$C associated with reduced stomatal conductance in response to the previous years’ stress event. These geochemical signatures of legacy were present despite tree growth returning to normal the year after the stress. The results show the importance of species-level dynamics in legacy that are absent in modeling that assumes common traits within plant functional types. In summary, the work here established new avenues to explore root dynamics in models while illustrating how additional root processes are needed before dynamic carbon allocation can be implemented. Lastly, this project provided mentorship to a postdoctoral fellow, training for an early career scientist and multiple undergraduate students recruited from a minority serving institution.

54 ENVIRONMENTAL SCIENCES↗