SEARCH · Engineering Papers

Results for “spatial heterogeneity”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2026

A Multi-Region SEIR Model Incorporating Inter-County Mobility and Time-Dependent Transmission Dynamics: Application to COVID-19 Disease Outbreak Data in North Carolina.

Classical infectious disease compartmental models typically do not incorporate spatial heterogeneity or mobility. We develop a multi-region susceptible-exposed-infected-recovered (SEIR) model in which disease dynamics are coupled via inter-region mobility and the transmission rate is both region and time dependent. We calibrate the model using rolling averages of daily COVID-19 data in all 100 North Carolina counties. Mobility parameters are prescribed using daily inter-county commuter data. The number of transmission rate parameters is substantially reduced by hypothesizing that the dynamics correlate with county-level population density. Parameter estimation is carried out using several objective functions with error terms at different scales. An additive combination of least squares error at the county-level and the state-level, along with a quadratic transmission rate polynomial, yields the lowest overall error at both spatial scales. The calibrated model is used to simulate regional effects of perturbing disease transmission rates in adjacent counties and to illustrate effects of the state’s mobility infrastructure on disease dynamics and spread for a new disease outbreak.

COVID-19 modeling↗

DOE OSTI2026

Permafrost Thaw, Uneven Subsidence and Projected Drying of Ice-wedge Polygon Tundra: Modeling Archive

This dataset is a model archive of the paper Permafrost Thaw, Uneven Subsidence and Projected Drying of Ice-wedge Polygon Tundra (in prep) to support a modeling study investigating how projected increases in Arctic temperature and precipitation will jointly influence hydrologic conditions in ice-rich tundra landscapes. With this dataset, this study is to address the research question: Will Arctic tundra landscapes become wetter or drier with increasing precipitation and temperature in the future when thaw-induced ground subsidence and associated microtopographic evolution are represented? The simulations focus on ice-wedge polygon tundra, a widespread form of ice-rich permafrost terrain that is highly sensitive to thaw-driven landscape change. This dataset contains model input and output data for four study watersheds in Alaska: Anaktuvuk, Utqiagvik (formerly Barrow), Brooks Foothills, and Prudhoe Bay. Simulations were performed using the Advanced Terrestrial Simulator (ATS, v1.5), a physics-rich integrated surface–subsurface hydrologic model. For each watershed, ten modeling cases were performed representing two landscape evolution conditions (with subsidence and without subsidence) combined with five climate forcing scenarios derived from Shared Socioeconomic Pathways (SSP5, SSP5 with precipitation trend, SSP2, SSP2 with precipitation trend, and SSP2 with double precipitation trend). Particularly, for each watershed under the forcing SSP2 with precipitation trend, there are two additional simulations considering spatially heterogeneous subsidence distributions: one assumes randomly distributed scaling and the other includes elevation dependent distribution scaling. These simulations span 1980–2099 and include spin-up runs (1980–2009) followed by transient projections (2010–2099). To facilitate reproducibility of simulations, all datasets are organized by watershed. For each study watershed, the dataset contains: (1) Pre-partitioned mesh files for 32-core modeling (.par.32.XX), located in EACH_WATERSHED/mesh/basin; and also a non-partitioned mesh file (.exo) located in EACH_WATERSHED/mesh; (2) Climate forcings corresponding to the five SSP scenarios (.h5), located in EACH_WATERSHED/data; (3) Final states (.h5) from column spin-up modeling used to initialize historical watershed-scale spin-up runs from 1980 to 2009, located in EACH_WATERSHED/PreSpinupHistorical; (4) Final states (.h5) of historical watershed-scale spin-up runs from 1980 to 2009 used to initialize projection runs, located in EACH_WATERSHED/Spinup_daymetERA5; (5) ATS modeling input files (.xml), located in EACH_WATERSHED/EACH_SIMULATION_SCENARIO/inputfiles; (6) ATS modeling output files (.dat), located in in EACH_WATERSHED/EACH_SIMULATION_SCENARIO/combined_obs; (7) For the Brooks Foothills watershed, additional spatial model outputs are provided (.h5) for selected years (2033 and 2093) used to generate spatial figures in this study, located in Brooksfoothills/EACH_SIMULATION_SCENARIO/results-WITH/WITHOUT_SUBSIDENCE-year2033/2093. All data files with suffix .h5 can be accessible through Python h5py, and all data files with suffix of .dat can be imported by Python pandas. Mesh file with .exo can be visualized through Paraview or read by Python netCDF. The Next-Generation Ecosystem Experiments in the Arctic (NGEE Arctic) project is a research effort to reduce uncertainty in the Department of Energy’s Energy Exascale Earth System Model (E3SM) by developing a predictive understanding of Arctic tundra ecosystems underlain by permafrost and to quantify feedbacks from the Arctic tundra to the Earth system. NGEE Arctic is supported by the Department of Energy's Office of Biological and Environmental Research. Over Phases 1–3, observations made by the NGEE Arctic team across a gradient of permafrost landscapes in Arctic Alaska improved the representation of tundra processes in the land surface component of E3SM (the E3SM Land Model, ELM). Model improvements emphasized unique aspects of permafrost environments and explored reductions in model complexity while retaining predictive power. The Arctic-informed ELM developed by NGEE Arctic has been used to make novel predictions on processes ranging from permafrost thaw to soil biogeochemical cycling to Earth system feedbacks associated with the unique characteristics of tundra plants. In Phase 4, the NGEE Arctic team is evaluating our new predictive understanding under novel conditions across the Arctic domain. In collaboration with partners at long-term pan-Arctic research sites we are examining whether an Arctic-informed ELM can faithfully simulate interactions among surface and subsurface processes at site, regional, and pan-Arctic scales. In turn, we are using variety of tools to dynamically extend and evaluate ELM inference, with an emphasis on data synthesis and pan-Arctic model evaluation, reintegration of code with an evolving E3SM, scaling across heterogeneous Arctic landscapes, and the appropriate representation of the impacts of increasingly frequent Arctic disturbances.

EARTH SCIENCE > ATMOSPHERE > PRECIPITATION↗

DOE OSTI2020

A Bezier Curve Informed Melt Pool Geometry to Model Additive Manufacturing Microstructures Using SPPARKS

Additive manufacturing is a transformative technology with the potential to manufacture designs which traditional subtractive machining methods cannot. Additive manufacturing offers fast builds at near final desired geometry; however, material properties and variability from part to part remain a challenge for certification and qualification of metallic components. AM induced metallic microstructures are spatially heterogeneous and highly process dependent. Engineering properties such as strength and toughness are significantly affected by microstructure morphologies resulting from the manufacturing process Linking process parameters to microstructures and ultimately to the dynamic response of AM materials is critical to certifying and qualifying AM built parts and components and improving the performance of AM materials. The AM fabrication process is characterized by building parts layer by layer using a selective laser melt process guided by a computer. A laser selectively scans and melts metal according to a designated geometry. As the laser scans, metal melts, fuses, and solidifies forming the final geometry in a layerwise fashion. As the laser heat source moves away, the metal cools and solidifies forming metallic microstructures. This work describes a microstructure modeling application implemented in the SPPARKS kinetic Monte Carlo computational framework for simulating the resulting microstructures. The application uses Bzier curves and surfaces to model the melt pool surface and spatial temperature profile induced by moving the laser heat source; it simulates the melting and fusing of metal at the laser hot spot and microstructure formation and evolution when the laser moves away. The geometry of the melt pool is quite flexible and we explore effects of variances in model parameters on simulated microstructures.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Hidden Features: How Subsurface and Landscape Heterogeneity Govern Hydrologic Connectivity and Stream Chemistry in a Montane Watershed

ABSTRACT Hydrologic connectivity is defined as the connection among stores of water within a watershed and controls the flux of water and solutes from the subsurface to the stream. Hydrologic connectivity is difficult to quantify because it is goverened by heterogeniety in subsurface storage and permeability and responds to seasonal changes in precipitation inputs and subsurface moisture conditions. How interannual climate variability impacts hydrologic connectivity, and thus stream flow generation and chemistry, remains unclear. Using a rare, four‐year synoptic stream chemistry dataset, we evaluated shifts in stream chemistry and stream flow source of Coal Creek, a montane, headwater tributary of the Upper Colorado River. We leveraged compositional principal component analysis and end‐member mixing to evaluate how seasonal and interannual variation in subsurface moisture conditions impacts stream chemistry. Overall, three main findings emerged from this work. First, three geochemically distinct end members were identified that constrained stream flow chemistry: reach inflows, and quick and slow flow groundwater contributions. Reach inflows were impacted by historic base and precious metal mine inputs. Bedrock fractures facilitated much of the transport of quick flow groundwater and higher‐storage subsurface features (e.g., alluvial fans) facilitated the transport of slow flow groundwater. Second, the contributions of different end members to the stream changed over the summer. In early summer, stream flow was composed of all three end members, while in late summer, it was composed predominantly of reach inflows and slow flow groundwater. Finally, we observed minimal differences in proportional composition in stream chemistry across all four years, indicating seasonal variability in subsurface moisture and spatial heterogeneity in landscape and geologic features had a greater influence than interannual climate fluctuation on hydrologic connectivity and stream water chemistry. These findings indicate that mechanisms controlling solute transport (e.g., hydrologic connectivity and flow path activation) may be resilient (i.e., able to rebound after perturbations) to predicted increases in climate variability. By establishing a framework for assessing compositional stream chemistry across variable hydrologic and subsurface moisture conditions, our study offers a method to evaluate watershed biogeochemical resilience to variations in hydrometeorological conditions.

Johnson, Keira [College of Earth, Ocean, and Atmos↗

DOE OSTI2022

The global decline in the sensitivity of vegetation productivity to precipitation from 2001 to 2018

The sensitivity of vegetation productivity to precipitation (S ppt ) is a key metric for understanding the variations in vegetation productivity under changing precipitation and predicting future changes in ecosystem functions. However, a comprehensive assessment of S ppt over all the global land is lacking. Here, we investigated spatial patterns and temporal changes of S ppt across the global land from 2001 to 2018 with multiple streams of satellite observations. We found consistent spatial patterns of S ppt with different satellite products: S ppt was highest in dry regions while low in humid regions. Grassland and shrubland showed the highest S ppt , and evergreen needle-leaf forest and wetland showed the lowest. Temporally, S ppt showed a generally declining trend over the past two decades (p < .05), yet with clear spatial heterogeneities. The decline in S ppt was especially noticeable in North America and Europe, likely due to the increase in precipitation. In central Russia and Australia, however, S ppt showed an increasing trend. Biome-wise, most ecosystem types exhibited significant decrease in S ppt , while grassland, evergreen broadleaf forest, and mixed forest showed slight increases or non-significant changes in S ppt . Our finding of the overall decline in S ppt implies a potential stabilization mechanism for ecosystem productivity under climate change. However, the revealed S ppt increase for some regions and ecosystem types, in particular global grasslands, suggests that grasslands might be increasingly vulnerable to climatic variability with continuing global climate change.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Discharge Rate‐Driven Li 2 O 2 Growth Exhibits Unconventional Morphology Trends in Solid‐State Li‐O 2 Batteries

Solid-state lithium oxygen batteries (LOBs) are known for their enhanced safety, higher electrochemical stability, and improved energy density compared to liquid-state LOBs. However, the investigation of solid-state LOBs is limited with little understanding of their discharge and charge processes. In this work, a polymer-based solid-state LOB is used to investigate the effect of discharge rate on lithium peroxide (Li 2 O 2 ) formation, the oxygen evolution reaction (OER), and cycle performance. Notably, we observe a counterintuitive trend: Li 2 O 2 particle size increases with increasing discharge current density, in contrast to liquid systems. This behavior arises from inherent space charge layers that restrict Li⁺ transport under high current, and spatially heterogeneous active sites at the solid electrolyte–cathode interface, directly evidenced by small angle X-ray scattering (SAXS), which govern nucleation accessibility and promote site-selective Li 2 O 2 growth. Furthermore, higher current densities improve ORR and OER efficiency but accelerate anode degradation, while lower currents promote side reactions. These opposing effects result in a trade-off that defines an optimal discharge rate (0.1 mA cm -2 ) for maximizing cycle life. This study provides a new mechanistic perspective on discharge-driven processes in solid-state LOBs and offers practical guidelines for performance optimization in future high-energy battery systems.

Discharge current density↗

DOE OSTI2026

Spatiotemporal Dynamics of the Relative Abundance of Soil Nutrient‐Degrading Enzyme‐Encoding Genes Across Continental US Ecoregions

Understanding the spatiotemporal patterns in the relative abundance of soil extracellular enzyme‐encoding genes is critical for predicting microbial responses to environmental change and their potential role in nutrient cycling. Yet, integrating novel metagenomic observations with spatiotemporal environmental gradients to infer regional patterns and future trajectories has remained unclear. To address this gap, we applied a machine learning (ML) approach, integrating soil metagenomic data with environmental variables—soil properties, topography, vegetation, and climate—to predict the relative abundance of enzyme‐encoding genes for soil carbon (C), nitrogen (N), and phosphorus (P) across surface soils of the continental United States. We assessed potential responses under future emission scenarios (SSP2‐4.5 and SSP5‐8.5) by comparing a baseline (1985–2014) to a future period (2071–2100). The ML model explained 57%–63% of baseline variation. Precipitation was identified as the most influential factor for the relative abundance of C‐ and N‐degrading enzyme‐encoding genes, while slope length, representing horizontal distance that water can travel downslope, was the primary driver for P‐degrading enzyme‐encoding genes abundance. Projections revealed spatially heterogeneous shifts across continental US ecoregions: the relative abundance of C‐ and N‐degrading enzyme‐encoding genes decreased in wetter ecoregions and increased in drier ecoregions under future climate, while P‐degrading enzyme‐encoding genes abundance decreased significantly in semiarid and Mediterranean ecoregions. This study demonstrates the utility of metagenomic data for mapping soil genetic potential and predicting its regional response to environmental change, to inform ecosystem management strategies.

extracellular enzyme-encoding genes↗

DOE OSTI2024

Localization of proteins involved in the biogenesis and repair of the photosynthetic apparatus to thylakoid subdomains in Arabidopsis

Abstract Thylakoid membranes in chloroplasts and cyanobacteria harbor the multisubunit protein complexes that catalyze the light reactions of photosynthesis. In plant chloroplasts, the thylakoid membrane system comprises a highly organized network with several subcompartments that differ in composition and morphology: grana stacks, unstacked stromal lamellae, and grana margins at the interface between stacked and unstacked regions. The localization of components of the photosynthetic apparatus among these subcompartments has been well characterized. However, less is known about the localization of proteins involved in the biogenesis and repair of the photosynthetic apparatus, the partitioning of proteins between two recently resolved components of the traditional margin fraction (refined margins and curvature), and the effects of light on these features. In this study, we analyzed the partitioning of numerous thylakoid biogenesis and repair factors among grana, curvature, refined margin, and stromal lamellae fractions of Arabidopsis thylakoid membranes, comparing the results from illuminated and dark‐adapted plants. Several proteins previously shown to localize to a margin fraction partitioned in varying ways among the resolved curvature and refined margin fractions. For example, the ALB3 insertase and FtsH protease involved in photosystem II (PSII) repair were concentrated in the refined margin fraction, whereas TAT translocon subunits and proteins involved in early steps in photosystem assembly were concentrated in the curvature fraction. By contrast, two photosystem assembly factors that facilitate late assembly steps were depleted from the curvature fraction. The enrichment of the PSII subunit OE23/PsbP in the curvature fraction set it apart from other PSII subunits, supporting the previous conjecture that OE23/PsbP assists in PSII biogenesis and/or repair. The PSII assembly factor PAM68 partitioned differently among thylakoid fractions from dark‐adapted plants and illuminated plants and was the only analyzed protein to convincingly do so. These results demonstrate an unanticipated spatial heterogeneity of photosystem biogenesis and repair functions in thylakoid membranes and reveal the curvature fraction to be a focal point of early photosystem biogenesis.

Chotewutmontri, Prakitchai↗

DOE OSTI2023

Improving evapotranspiration computation with electrical resistivity tomography in a maize field

Abstract Hydrogeophysical methods have been increasingly used to study subsurface soil–water dynamics, yet their application beyond the soil compartment or the quantitative link to soil hydraulic properties remains limited. To examine how these methods can inform model‐based evapotranspiration (ET) calculation under varying soil water conditions, we conducted a pilot‐scale field study at an experimental maize plot with manipulated irrigation treatments. Our goal was to develop a workflow for (1) acquiring and inverting field electrical resistivity tomography (ERT) data, (2) correlating ERT to soil hydraulic properties, (3) spatially characterizing soil water stress that feeds into ET modeling (the FAO‐56 model), and (4) evaluating the performance of ERT‐based ET computation. Our results showed that ERT was able to capture decimeter‐scale soil water content (SWC) dynamics from root water uptake and irrigation manipulation and the contrast of soil water stress between deficiently and fully irrigated maize. We also demonstrated the flexibility of using ERT to spatially integrate soil water stress in the soil volume of interest, which could be adjusted based on different crops and plot layouts. The integration of the ERT datasets into ET modeling provided insights into the spatial heterogeneity of the subsurface that has been challenging for point‐based sensing, which can further our understanding of the hydraulic dynamics in the soil‐plant‐atmosphere continuum.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

The Observed and Projected Changes of Global Monsoons: Current Status and Future Perspectives

The global monsoon system, encompassing the Asian-Australian, African, and American monsoons, sustains two-thirds of the world’s population by regulating water resources and agriculture. Monsoon anomalies pose severe risks, including floods and droughts. Recent research associated with the implementation of the Global Monsoons Model Intercomparison Project under the umbrella of CMIP6 has advanced our understanding of its historical variability and driving mechanisms. Observational data reveal a 20th-century shift: increased rainfall pre-1950s, followed by aridification and partial recovery post-1980s, driven by both internal variability (e.g., Atlantic Multidecadal Oscillation) and external forcings (greenhouse gases, aerosols), while ENSO drives interannual variability through ocean-atmosphere interactions. Future projections under greenhouse forcing suggest long-term monsoon intensification, though regional disparities and model uncertainties persist. Models indicate robust trends but struggle to quantify extremes, where thermodynamic effects (warming-induced moisture rise) uniformly boost heavy rainfall, while dynamical shifts (circulation changes) create spatial heterogeneity. Volcanic eruptions and proposed solar radiation modification (SRM) further complicate predictions: tropical eruptions suppress monsoons, whereas high-latitude events alter cross-equatorial flows, highlighting unresolved feedbacks. The emergent constraint approach is booming in terms of correcting future projections and reducing uncertainty with respect to the global monsoons. Critical challenges remain. Model biases and sparse 20th-century observational data hinder accurate attribution. The interplay between natural variability and anthropogenic forcings, along with nonlinear extreme precipitation risks under warming, demands deeper mechanistic insights. Additionally, SRM’s regional impacts and hemispheric monsoon interactions require systematic evaluation. Addressing these gaps necessitates enhanced observational networks, refined climate models, and interdisciplinary efforts to disentangle multiscale drivers, ultimately improving resilience strategies for monsoon-dependent regions.

climate extreme events↗

DOE OSTI2025

Radiochronometric discordance in cast uranium metal: a multi-laboratory intercomparison exercise

Here, the model age of a nuclear material is crucial in nuclear forensic analysis. Uranium metals with complex production histories often exhibit discordant model ages from the 230 Th– 234 U and 231 Pa– 235 U chronometers. Recent studies involving targeted uranium metal castings have enhanced our understanding of decay product behavior during casting, aiding nuclear forensic interpretation. Building on this prior work, forensics laboratories at Atomic Weapons Establishment (AWE), Lawrence Livermore National Laboratory (LLNL), and Los Alamos National Laboratory (LANL) conducted an interlaboratory comparison to investigate spatial heterogeneity in uranium metal cast under controlled conditions. Each laboratory measured samples of a mixed feedstock and its corresponding cast product. This work furthers our understanding of discordant model ages and the use of discordance as a signature to enhance confidence in interpretations of radiochronometric data for nuclear forensics.

230Th/234U↗

DOE OSTI2023

Ultrasound and X-ray Cross-Characterization of a Graded Impedance Impactor used for Shock-Ramp Compression Experiments

Abstract In this work we perform ultrasound measurements on an impedance graded impactor made by tape casting magnesium, copper, and tungsten. We also destructively extract small representative samples from the part for complementary characterization with x-ray computed tomography. Combining the two data sets enables direct assignment of some of the measured ultrasound features to specific material characteristics identified by x-ray tomography. Our results demonstrate how ultrasound inspection, informed by x-ray computed tomography, can be used to identify sub-millimeter material amalgamations and spatial heterogeneities in this graded material.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Direct observations and characterization of crack closure during microstructurally small fatigue crack growth via in-situ high-energy X-ray characterization

This study provides direct observation of the crack closure mechanism of a naturally occurring, tortuous, 3D microstructurally small fatigue crack (SFC) in additively manufactured Inconel 718. In-situ non-destructive characterization is performed using high-energy X-ray diffraction techniques to capture the evolution of the 3D microstructure and micromechanical response in the vicinity of the crack front. Based on the stress state of twelve grains analyzed at the crack tip, the crack closure events of the SFC front was found to be spatially heterogeneous with respect to loading progression governed by the local stress state of the grains (specifically the stress reversal from compression to tension). From this analysis, three grains that displayed different degrees of crack closure were further investigated, based on the orientation of the crack relative to the grains and the associated modality of crack growth. The stress normal to the crack plane and the associated degree of Mode I crack behavior were correlated with events of the crack opening earlier during the loading cycle. This was further corroborated by diffraction spot spreading analysis that quantified the crystallographic lattice distortion caused by the opening crack. Additionally, the detailed characterization of the opening behavior of the grains located at the crack tip and their associated states of stress elucidates the mechanism governing crack closure and will inform future modeling efforts of this phenomenon.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Site-specific field management adaptation is key to feeding the world in the 21st century

Rapid climate change and growing population threaten global food security across the globe. Several studies have proposed early planting, increased irrigation, and increased fertilizer applications as climate adaptation strategies, yet none have considered combined and site-specific field management strategies as a comprehensive solution. Here, we analyzed non-irrigated wheat yield responses to climate change and field management adaptation using a mechanistic crop model evaluated against observed global non-irrigated wheat-yield over 3-year intervals spanning 13 years at 3749 sites (RMSE = 36 gC m -2 ). Early planting with later-maturing varieties provided the most benefit to future yields among the proposed field management adaptation strategies. Improved water use efficiency from increased CO 2 led to relatively low benefits of additional irrigation. Here we estimated that spatially heterogeneous adaption strategies had the potential to improve global wheat yields by 91% by 2100 compared to the present day. The yield improvements from combined field adaptation strategies were larger than the sum of improvements from the individual strategies. These synergistic benefits were shown to result from complementary processes regulating nutrient and water uptake, physiological tolerance to heat stress, and internal carbon and nutrient cycling.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Dust-planetary boundary layer interactions amplified by entrainment and advections

Mineral dust contributes to more than half of the global aerosol loading. However, the radiative impacts of dust aerosols on planetary boundary layer (PBL) structure have not been explored sufficiently. During a typical dust storm event over Tarim Basin, dust aerosols exhibit a well-mixed distribution during the daytime in spite of a shallow layer of dust particles accumulated at higher altitudes. By contrast, nocturnal dust plumes are located near the surface due to stable stratification. We demonstrate that these differentiated vertical distributions determine the spatial heterogeneity of dust loading, radiative fluxes and PBL height variations. Dust aerosols cause daytime PBL suppression and nighttime PBL promotion through modulating surface and atmospheric radiative budgets. Specifically, dust-induced cooling effect within PBL directly inhibits the daytime PBL development. PBL suppression effect is then amplified by entrainment processes resulting in excessively low PBL height, especially for dust particles below but near the PBL top. Dust plumes weaken updrafts from PBL and downdrafts of the free atmosphere, which further reduce the entrainment mixing through attenuating horizontal and vertical advection, and eventually amplify PBL suppression. At night, near-surface dust aerosols stimulate a warm and unstable lower atmosphere, generate warm advection heating and promote the PBL development. Our study highlights the importance of specifying entrainment parameters and nighttime advection activities in quantifying the dust-PBL interactions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Spatiotemporal 4D Whole-cell Modeling of a Minimal Autotroph Reveals Central Carbon Metabolism Regulated Locally by Protein Megacomplexes via Post-translational Modifications under Light Disturbance

Photosynthetic microorganisms rely on multiple pathways in central carbon metabolism to adapt to fluctuating light and energy availability across diel cycles. Mechanistic insight into the regulatory dynamics of this adaptation requires integrating processes spanning disparate timescales, from rapid redox-dependent post-translational modifications (PTMs) to slower changes in protein expression and metabolic pathway usage. To address this complexity beyond genome-based inference and traditional modeling, we develop a whole-cell four-dimensional (3D + time) model of the marine cyanobacterium Prochlorococcus marinus MED4 that explicitly represents the spatial organization of enzymatic and molecular processes in central carbon metabolism under light perturbation. We employ a perturbation-based research design to experimentally generate time-series, multi-omics measurements that provide molecular descriptors and cryo-ET derived 3D segmented volumes as constraints for this dynamic 4D framework. The integration of experiments and modeling across defined light regimes enables quantitative validation of system-level responses and forecasting under distinct light disturbances. We test the hypothesis that light-dependent redox PTMs regulating the structural assembly of a protein megacomplex, the “dark complex,” modulate metabolic flux at a conserved regulatory node of the Calvin–Benson cycle (CBC) in cyanobacteria. Our model shows that subcellular spatial organization buffers rapid light-induced changes in thylakoid reaction rates, which are followed by redox-PTM-mediated sequestration or release of CBC enzymes in the dark complex, ultimately impacting carbon fixation dynamics within carboxysomes. Comparison with an equivalently parameterized well-mixed stochastic model demonstrates that post-translational regulation not only buffers transcriptional noise and diffusion-driven fluctuations but also stabilizes phenotypic outcomes, underscoring the importance of spatial heterogeneity in phenotypic robustness. This ability to probe adaptive, spatiotemporally resolved mechanisms in photosynthetic machinery and central carbon metabolism addresses a critical gap in genotype-to-phenotype inference and expands modeling and design capabilities for understudied or genetically intractable autotrophs such as P. marinus MED4.

Johnson, Connah G.↗

DOE OSTI2022

Plutonium mobilization from contaminated estuarine sediments, Esk Estuary (UK)

Since 1952, liquid radioactive effluent containing 238-242 Pu, 241 Am, 237 Np, 137 Cs, and 99 Tc has been released with authorization from the Sellafield nuclear complex (UK) into the Irish Sea. This represents the largest source of plutonium (Pu) discharged in all western Europe, with 276 kg having been released. In the Eastern Irish Sea, the majority of the transuranic activity has settled into an area of sediments (Mudpatch) located off the Cumbrian coast. Radionuclides from the Mudpatch have been re-dispersed via particulate transport in fine-grained estuarine and intertidal sediments to the North-East Irish Sea, including the intertidal saltmarsh located at the mouth of the Esk Estuary. Saltmarshes are highly dynamic systems which are vulnerable to external agents (sea level change, erosion, sediment supply, and freshwater inputs), and their stability remains uncertain under current sea level rise projections and possible increases in storm activity. In this work, we examined factors affecting Pu mobility in contaminated sediments collected from the Esk Estuary by conducting leaching experiments under both anoxic and oxic conditions. Here, leaching experiments were conducted over a 9-month period and were periodically sampled to determine solution phase Pu via multicollector-inductively coupled plasma-mass spectrometry (MC-ICP-MS), and to measure redox indicators (Eh, pH and extractable Fe(II)). Microbial community composition was also characterized in the sediments, and at the beginning and end of the anoxic/oxic experiments. Results show that: 1) Pu leaching is about three times greater in solutions leached under anoxic conditions compared to oxic conditions, 2) the sediment slurry microbial communities shift as conditions change from anoxic to oxic, 3) Pu leaching is enhanced in the shallow sediments (0–10 cm depth), and 4) the magnitude of Pu leached from sediments is not correlated with total Pu, indicating that the biogeochemistry of sediment-associated Pu is spatially heterogeneous. These findings provide constraints on the stability of redox sensitive Pu in biogeochemically dynamic/transient environments on a timescale of months and suggests that anoxic conditions can enhance Pu mobility in estuarine systems.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2026

Machining-inducted residual strain and ductile-brittle transition in sapphire revealed by X-ray nanoprobe diffraction microscopy

Hard X-ray nanoprobe diffraction microscopy is employed to investigate machining-induced residual lattice strain in single-crystal sapphire following ultra-precision orthogonal plunge cutting. In this study, tilt-series and single-angle acquisition strategies are systematically compared to evaluate their capability for resolving lattice variations associated with the machining process. An analytical framework is developed to enable the efficient and reliable extraction of strains, which serve as strain-based indicators of changes in the dominant deformation mechanism. The measured high-dimensional strain exhibits spatial heterogeneity and systematic evolution along the cutting path, consistent with deformation-regime changes from ductile response to crack-prone behavior.

Liang, Rui↗