Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational”

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,783 records · Page 99
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Controlling the Morphology of PEDOT:PSS Blend Films with Pre-Deposition Solution Composition and Deposition Technique

Understanding the relationships between morphology, fabrication processes, and thermoelectric performance in conducting polymers is essential to the development of high- efficiency organic thermoelectrics as an alternative to commonly used rare metals. Altering the film fabrication process of poly(3,4-ethylenedioxythiophene):poly(styrene sulfonate) (PEDOT:PSS) with the addition of high boiling solvents to the precast solution improves the electrical conductivity and significantly increases its Seebeck value. Neutron scattering monitors the changes in the atomic, nanoscale, and mesoscale morphologies of PEDOT:PSS thin films with the addition of dimethyl sulfoxide (DMSO) to the aqueous solution prior to film formation and with varying fabrication procedures. The neutron scattering results show a decrease in the deuterated PSS domain size along with systematic variations in PEDOT fibril assemblies in the final blend film with the addition of DMSO to the pre-deposition solution. These structural modifications indicate that the increase in conductivity of PEDOT:PSS blends with addition of DMSO reported in the literature can be ascribed to the disruption of solvated PEDOT assemblies by the DMSO, forming smaller PSS domains in the pre-deposition solution and allowing smoother film formation. These improvements are observed significantly with the addition of just 1% DMSO but continue to modestly improve with the addition of up to 5% DMSO to the PEDOT:PSS blend pre-deposition solution. The fact that the variations in the measured morphology are independent of whether the films were deposited by spin or ultrasonic spray casting methods emphasizes the crucial importance of the structure of the blend in the pre-deposition solution in determining the final thin film blend morphology.

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

Influence of Grafting Density and the Ionic Environment on the Structure of Zwitterionic Brushes

Zwitterionic polymer brushes possess a high potential for applications as surface coatings, e.g., in antifouling applications. Their complex association behavior due to the coexistence of oppositely charged groups on the same monomeric unit allows for a broad variation of polarity, hydrophilicity, and the antipolyelectrolyte effect with the variation of the surrounding environment. In this study, planar polysulfobetaine brushes were investigated with neutron reflectometry (NR) and neutron spin–echo spectroscopy under grazing incidence (GINSES) to explore the brush structure and its inner dynamics perpendicular to the substrate. In particular, the effects of the substrate–initiator system and ionic strength on the structure were investigated with neutron reflectometry and a swelling by a factor of 3–5 was found in the presence of aqueous NaCl and MgCl 2 solutions compared to the dry state. During the data analysis, the applicability of a model-free evaluation was also demonstrated for the investigated polysulfobetaine brushes. Furthermore, it was found with GINSES that in contrast to neutral brushes, these polysulfobetaine brushes do not show the typical concentration fluctuations in the nanosecond time range when partially swollen with salt-free water.

36 MATERIALS SCIENCE↗
DOE OSTI2020

High-Throughput Characterization of (Fe x Co 1–x ) 3 O 4 Thin-Film Composition Spreads

In this study, thin-film continuous composition spreads of Fe–Co–O were fabricated by reactive cosputtering from elemental Fe and Co targets in reactive Ar/O 2 atmosphere using deposition temperatures ranging from 300 to 700 °C. Fused silica and platinized Si/SiO 2 strips were used as substrates. Ti and Ta were investigated as adhesion layer for Pt and the fabrication of the Fe–Co–O films. The thin-film composition spreads were characterized by high-throughput electron-dispersive X-ray spectroscopy, X-ray diffraction, X-ray photoelectron spectroscopy, atomic force microscopy, scanning electron microscopy, and optical transmission spectroscopy. The Fe-content ranged from 28 to 72 at. %. The spinel phases Fe 2 CoO 4 and FeCo 2 O 4 could be synthesized and stabilized at all deposition temperatures with a continuous variation in spinel composition in between. The dependence of the film surface microstructure on the deposition temperature and the composition was mapped. Moreover, the band gap values, ranging from 2.41 eV for FeCo 2 O 4 to 2.74 eV for Fe 2 CoO 4 , show a continuous variation with the composition.

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

Investigating Particle Phase State Dependencies during a Historic Pacific Northwest Wildfire Event

The function and lifetime of an atmospheric particle are greatly dependent upon its phase state. Studying particle phase states from specific ambient source types is often challenging due to uncontrolled environmental conditions and multiple source contributions. However, a historic Pacific Northwest wildfire event occurred in September 2020, providing a rare opportunity to study a pseudo-single particle type ambient sample (i.e., wildfire) over time. Sampling was performed in Richland, Washington, capturing aged particles from shifting wildfire sources. In addition to source variability, temporal variation in wildfire particle phase state was observed from single particle analyses. Despite this, single particle elemental compositions were homogenous over time (=87% carbonaceous) with no significant temporal variation in organic functional group abundances within individual particles either (e.g., 4% relative standard deviation in alkene contributions). While particle composition is indeed known to affect phase state, water uptake appeared to be the most significant contributor to phase state variability amongst single particles here, as particle aspect ratios (related to phase state) linearly correlated with ambient relative humidity during sampling (R2 = 0.65). Therefore, while limited to a case study, meteorological considerations are found to convey greater comparative importance than particle composition for aged particle phase state predictions.

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

Confinement and Processing Can Alter the Morphology and Periodicity of Bottlebrush Block Copolymers in Thin Films

Bottlebrush block copolymers (BBCPs) are intriguing architectural variations on linear BCPs with highly tunable structure. Confinement can have a significant impact on polymer assembly, giving rise to changes in morphology, assembly kinetics, and properties like the glass transition. Given that confinement leads to significant changes in the persistence length of bottlebrush homopolymers, it is reasonable to expect that BBCPs will see significant changes in their structure and periodicity relative to the bulk morphology. Understanding how confinement influences assembly will be important for designing BBCPs for thin film applications including membranes, integrated photonic structures, and potentially BCP lithography. In this work, in order to study the effects of confinement on BBCP conformation and morphology, a blade coating was used to prepare films with continuous variation in film thickness. Unlike thin films of linear BCPs, islands/holes were not observed, and instead mixtures of parallel and perpendicular morphologies emerge after annealing. The lamellar periodicity ($L_0$) of the morphologies is found to be thickness dependent, increasing $L_0$ with decreasing film thickness for blade coated films. Films coated out of tetrahydrofuran (THF) resulted in a single well-defined lamellar periodicity, verified through atomic force microscopy (AFM) and grazing incidence small-angle X-ray scattering (GISAXS), which increases dramatically from the bulk value (30.6 nm) and continues to increase as the film thickness decreases. The largest observed $L_0$ was 65.5 nm, and this closely approaches the estimated upper limit of 67 nm corresponding to a fully extended backbone in a bilayer arrangement. Films coated out of propylene glycol methyl ether acetate (PGMEA) resulted in a mixture of perpendicular lamellae and a smaller, likely cylindrical morphology. The lamellar portion of the film shows the same thickness dependence as the lamellae observed in the THF coated films. The scaling of the lamellar $L_0$ with respect to film thickness follows predictions for confined semiflexible polymers with weak excluded volume interactions and can be related to models for confinement of DNA. Spin coated films shows the same reduction in periodicity, although at very different film thicknesses. This result suggests that the material has shallow free-energy barriers to transitioning between different $L_0$ and morphologies, a property that could be taken advantage of for patterning diverse structures with a single material.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Temperature-Dependent Properties of Molten Li 2 BeF 4 Salt Using Ab Initio Molecular Dynamics

Molten lithium tetrafluoroberyllate (Li 2 BeF 4 ) salt, also known as FLiBe, with a 2:1 mixture of LiF and BeF 2 is being proposed as a coolant and solvent in advanced nuclear reactor designs, such as the molten salt reactor or the fluoride salt cooled high-temperature reactor. We present the results on the structure and properties of FLiBe over a wide range of temperatures, 0–2000 K, from high-throughput ab initio molecular dynamics simulation using a supercell model of 504 atoms. The variations in the local structures of solid and liquid FLiBe with temperature are discussed in terms of a pair distribution function, coordination number, and bond angle distribution. The temperature-dependent electronic structure and optical and mechanical properties of FLiBe are calculated. The optical and mechanical property results are reported for the first time. The results above and below the melting temperature (~732 K) are compared with the experimental data and with data for crystalline FLiBe. The electronic structure and interatomic bonding results are discussed in correlation with the mechanical strength. A novel concept of total bond order density (TBOD), an important quantum mechanical parameter, is used to characterize the internal cohesion and strength in the simulated models. The results show a variation in the rate of change in properties in solid and liquid phases with anomalous behavior across the melting region. The observed trend is the decrease in mechanical strength, band gap, and TBOD in a nonlinear fashion as a function of temperature. The refractive index shows a surprising minimum at 850 K, among the tested temperatures, which lies above the melting point. These findings provide a new platform to understand the interplay between the temperature-dependent structures and properties of FLiBe salt.

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

Comparative Pore Structure and Dynamics for Bacterial Microcompartment Shell Protein Assemblies in Sheets or Shells

Bacterial microcompartments (BMCs) are protein-bound organelles found in some bacteria that encapsulate enzymes for enhanced catalytic activity. These compartments spatially sequester enzymes within semipermeable shell proteins, analogous to many membrane-bound organelles. The shell proteins assemble into multimeric tiles; hexamers, trimers, and pentamers, and these tiles self-assemble into larger assemblies with icosahedral symmetry. While icosahedral shells are the predominant form in vivo , the tiles can also form nanoscale cylinders or sheets. The individual multimeric tiles feature central pores that are key to regulating transport across the protein shell. Our primary interest is to quantify pore shape changes in response to alternative component morphologies at the nanoscale. We used molecular modeling tools to develop atomically detailed models for both planar sheets of tiles and curved structures representative of the complete shells found in vivo . Subsequently, these models were animated using classical molecular dynamics simulations. From the resulting trajectories, we analyzed the overall structural stability, water accessibility to individual residues, water residence time, and pore geometry for the hexameric and trimeric protein tiles from the Haliangium ochraceu m model BMC shell. These exhaustive analyses suggest no substantial variation in pore structure or solvent accessibility between the flat and curved shell geometries. We additionally compare our analysis to hydroxyl radical footprinting data to serve as a check against our simulation results, highlighting specific residues where water molecules are bound for a long time. Although with little variation in morphology or water interaction, we propose that the planar and capsular morphology can be used interchangeably when studying permeability through BMC pores.

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

Experimental Investigation of Uranium and Iron Condensation from High-Temperature Plasma Conditions

We used a plasma flow reactor (PFR) to generate synthetic fallout nanoparticles from vapor-phase condensation of two different input concentrations of uranium and iron analytes (U/Fe = 1:1 and 1:2). Synthetic fallout from complex chemical matrices (e.g., mixtures of U and Fe) has not been generated in this setup before and allows for the observation of relative condensation and fractionation of nuclear debris. Experiments were conducted under two different temperature histories with variations in particle flow patterns along the PFR. Transmission electron microscopy (TEM) observation and analysis of the nanoparticles revealed variations in speciation of uranium oxides (UO 2 and α-UO 3 ) depending on the competition between flow mixing and oxygen sequestration by iron. A ternary metal oxide, UFeO 4 , was observed in addition to iron oxides (e.g., FeO and Fe 3 O 4 ), which suggests that fallout models should account for chemical speciation of ternary metal oxides (i.e., UFeO 4 ) and their relative condensation behaviors in addition to those of singular metal oxides (e.g., FeO and UO 2 ). X-ray Energy Dispersive Spectroscopy (EDS) elemental maps showed that some particles had Fe-rich cores surrounded by U-rich regions. This suggests that either condensed U oxides coagulate onto molten Fe oxides or that the increase in iron analyte concentration might drive the system toward a higher degree of supersaturation, leading to earlier formation of iron oxide particles and providing an energetically favored pathway for nucleation of uranium oxides around iron oxide particles.

and nuclear chemistry↗
DOE OSTI2023

Cyclic versus Linear Alkylammonium Cations: Preventing Phase Transitions at Operational Temperatures in 2D Perovskites

Two-dimensional lead halide perovskites offer numerous attractive features for optoelectronics owing to their soft, deformable lattices and high degree of chemical tunability. While alteration of the metal and halide ions gives rise to significant modification of the bandgap energy, the organic spacer cations offer in-roads to tuning phase behavior and more subtle functionalities in ways that remain to be understood. Here, we study six variations of 2D perovskites changing only the organic spacer cations and demonstrate that these components intrinsically impact material response in important ways such as altering crystallographic structure, temperature-induced phase transitions, and photoluminescence emission. Two-dimensional perovskites containing commonly utilized aliphatic linear spacers, such as butylammonium, undergo phase transitions near room temperature. These transitions and temperature changes induce spacer-dependent variations in the emission spectra. Conversely, 2D perovskites comprising cyclic aliphatic spacers, such as cyclobutylammonium, are found to lack first-order phase transitions. These cyclic molecules are more sterically hindered within the crystal lattice, leading to temperature-induced contraction or expansion along certain crystallographic planes but no other significant thermal effects; additionally, they undergo changes in their emission spectra that cannot be explained by simple thermal expansion. Finally, given the similarities in the dielectric and chemical makeup of this set of six alkylammonium molecules, these results are unexpected and suggest a large structural and thermal phase space via spacer manipulation that could lead to improved 2D perovskite functionalization.

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

Precise Modulation of CO 2 Sorption in Ti 8 Ce 2 –Oxo Clusters: Elucidating Lewis Acidity of the Ce Metal Sites and Structural Flexibility

Investigating the structure-property correlation in porous materials is a fundamental and consistent focus in various scientific domains, especially within sorption research. Metal oxide clusters with capping ligands, characterized by intrinsic cavities formed through specific solid-state packing, demonstrate significant potential as versatile platforms for sorption investigations due to their precisely tunable atomic structures and inherent long-range order. This study presents a series of Ti 8 Ce 2 -oxo clusters with subtle variations in coordinated linkers and explores their sorption behavior. Notably, Ti 8 Ce 2 -BA (BA denotes benzoic acid) manifests a distinctive two-step profile during CO 2 adsorption, accompanied by a hysteresis loop. This observation marks a new instance within the metal oxide cluster field. Of intrigue, the presence of unsaturated Ce(IV) sites was found to be correlated with the stepped sorption property. Moreover, the introduction of an electrophilic fluorine atom, positioned ortho or para to the benzoic acid, facilitated precise control over gate pressure and stepped sorption quantities. Advanced in-situ techniques systematically unraveled the underlying mechanism behind this unique sorption behavior. The findings elucidate that robust Lewis base-acid interactions are established between CO 2 molecules and Ce ions, consequently altering the conformation of coordinated linkers. Conversely, the F atoms primarily contribute to gate pressure variation by influencing the Lewis acidity of the Ce sites. This research advances the understanding in fabricating metal-oxo clusters with structural flexibility and provides profound insights into their host-guest interaction motifs. Furthermore, these insights hold substantial promise across diverse fields and offer valuable guidance for future adsorbent designs grounded in fundamental theories of structure-property relationships.

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

High-Throughput Discovery Illuminates Design Principles and Limits for Long-Lived Charged Species in Organic Electrolytes

The chemical stability of charged molecules in all-organic redox flow batteries (RFBs) is required for the prolonged operation of these devices. Molecular engineering and electrolyte optimization are used to mitigate parasitic reactions and extend the lifetimes of the charge carriers. However, how much can structural variation extend the lifetime? To probe this query, we designed a high-throughput kinetic study of the radical cation of N-methylphenothiazinium, guided by statistical sampling and learning algorithms. Using Argonne’s autonomous discovery facility, we conducted over 6,000 kinetic experiments with robotic sample preparation, parallel kinetic measurements, and machine learning inputs, testing 188 solvent molecules selected from a space of over 540 candidates from 11 chemical classes. Algorithmic selections guided us to stable solvent candidates, which were further tested in high concentration with and without supporting electrolyte. Our findings reveal the inherent difficulty of exceeding the current state of the art through solvent variation. The desired stability is statistically rare and poorly predictable. Among the many tested, only three solvents significantly outperformed our baseline, acetonitrile─and none by more than a factor of 3─suggesting a general challenge in achieving the necessary techno-economic targets. Furthermore, we suggest that self-discharge through solvent homolysis is the cause of the observed limitations. Several structural motifs contribute to >1,000 h half-life stability including molecular simplicity, symmetry, oxidation complement, and strategic fluorination. Importantly, this workflow establishes effective assays for diagnosing and predicting oxidative stress for highly stable liquid electrolytes in all batteries.

Batteries↗
DOE OSTI2025

Chromium versus Aluminum: Impact of Nickel Alloy Composition and Interfacial Kinetics on High-Temperature Passivating Oxide Formation

High-temperature corrosion resistance depends critically on the formation of a passivating surface oxide, which is highly sensitive to alloy composition and structure. Such details often elude experimental investigation, and simplified analytical models fail to provide a truly chemical view of passivating oxide evolution. Here, we explicitly compare the fundamental chemistry of Cr and Al as prototypical passivating elements in Ni alloys by directly simulating competing reaction and diffusion processes within the oxide film using kinetic Monte Carlo and density functional theory. We find that the origin and expression of passivating behavior during early-stage thermal oxidation are qualitatively different between the two alloy systems. Ni–Cr alloys feature a sudden onset of passivation associated with a sharp phase transition upon Cr enrichment that directly couples oxidation kinetics to phase transformation behavior. In contrast, Ni–Al alloys display more continuous oxide phase variation with Al enrichment, ultimately resulting in a lower composition threshold for passivation and a thinner passivating layer. In addition, we elucidate the nonobvious role of metal exchange within the alloy near the oxide boundary, which fundamentally alters film composition and passivating behavior. Furthermore, our results have key implications for engineering improved corrosion-resistant alloys, both in terms of compositional variation and processing.

Alloys↗
DOE OSTI2026

Photoluminescent Cu(I) HETPHENs Featuring Bulky Alkyl-Substituted Phenanthrolines

Successful excited-state-enhancing substituent effect strategies applied in homoleptic Cu(I) metal-to-ligand charge transfer (MLCT) chromophores have been adapted to the heteroleptic phenanthroline (HETPHEN) platform, leveraging 2,9-mesityl-1,10-phenanthroline (mesPhen) in conjunction with a series of seven distinct 2,9- and 2,3,4,7,8,9-substituted phenanthrolines. The newly conceived CuHETPHEN complexes feature MLCT lifetimes ranging from 62 to 443 ns, all of which exhibit unprecedented room-temperature photoluminescence and demonstrate quantitative adherence to the energy gap law. TD-DFT calculations successfully modeled the systematic variation in electronic transition intensities, accounting for the experimentally observed UV–vis absorption bands across the entire series of molecules while providing detailed structural explanations for the variations in spectral profiles. These heteroleptic Cu(I) diimine chromophores exhibit thermally activated delayed fluorescence (TADF), with delayed fluorescence occurring between closely spaced 1 MLCT and 3 MLCT excited manifolds, having energetic separations of ΔE = 713–1009 cm –1 , as demonstrated here for the first time. Nanosecond and ultrafast transient absorption spectroscopy verified the MLCT nature of these excited states and extracted the time constants for the initial pseudo-Jahn–Teller distortion and intersystem crossing throughout the entire series of molecules. Lastly, the triplet photosensitization properties of these Cu(I) HETPHENs are demonstrated in a number of model photochemical transformations.

Chromophores↗
DOE OSTI2020

Detection Limits and Near‐Field Ground Motions of Fast and Slow Earthquakes

Abstract We investigate theoretical limits to detection of fast and slow seismic events, and spatial variations of ground motion expected from M 6 earthquakes at short epicentral distances. The analyses are based on synthetic velocity seismograms calculated with the discrete wavenumber method assuming seismic velocities and attenuation properties of the crust in Southern California. The examined source properties include different magnitudes ( M −1.0 to M 6.0), static stress drops (0.1–10 MPa), and slow and fast ruptures (0.1–0.9 of shear wave velocity). For the M 6 events we also consider variations in rise times producing crack‐ and pulse‐type events and different rupture directivities. Slow events produce ground motion with considerably lower amplitude than corresponding regular fast earthquakes with the same magnitude, and hence are significantly more difficult to detect. The static stress drop and slip rise time also affect the maximum radiated seismic motion, and hence event detectability. Apart from geometrical factors, the saturation and depletion of seismic ground motion at short epicentral distances stem from radiation pattern, earthquake size (magnitude, stress drop), and rupture directivity. The rupture velocity, rise time, and directivity affect significantly the spatial pattern of the ground motions. The results can help optimizing detection of slow and fast small earthquakes and understand the spatial distribution of ground motion generated by large events.

58 GEOSCIENCES↗
DOE OSTI2020

Stable Methane Isotopologues From Northern Lakes Suggest That Ebullition Is Dominated by Sub‐Lake Scale Processes

Abstract Stable isotopes have emerged as popular study targets when investigating emission of methane (CH 4 ) from lakes. Yet little is known on how isotopic patterns conform to variations in emission magnitudes—a highly relevant question. Here, we present a large multiyear data set on stable isotopes of CH 4 ebullition (bubbling) from three small adjacent subarctic lakes. The δ 13 C‐CH 4 and δD‐CH 4 range from −78.4‰ to −53.1‰ and from −369.8‰ to −218.8‰, respectively, and vary greatly among the lakes. The signatures suggest dominant hydrogenotrophic methanogenesis, particularly in the deep zones, but there are also signals of seemingly acetoclastic production in some high fluxing shallow areas, possibly fueled by in situ vegetation, but in‐sediment anaerobic CH 4 oxidation cannot be ruled out as an alternative cause. The observed patterns, however, are not consistent across the lakes. Neither do they correspond to the spatiotemporal variations in the measured bubble CH 4 fluxes. Patterns of acetoclastic and hydrogenotrophic production plus oxidation demonstrate that gains and losses of sediment CH 4 are dominated by sub‐lake scale processes. The δD‐CH 4 in the bubbles was significantly different depending on measurement month, likely due to evaporation effects. On a larger scale, our isotopic data, combined with those from other lakes, show a significant difference in bubble δD‐CH 4 between postglacial and thermokarst lakes, an important result for emission inventories. Although this characteristic theoretically assists in source partitioning studies, most hypothetical future shifts in δD‐CH 4 due to high‐latitude lake area or production pathway are too small to lead to atmospheric changes detectable with current technology.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Kilometer-Scale Hydrologic Exchange Flows in a Gravel Bed River Corridor and Their Implications to Solute Migration

In this work, a well-characterized field site along a major river corridor was used to investigate the dynamic pathways and impacts of subsurface hydrogeologic structures on kilometer-scale hydrologic exchange flows between river water and groundwater. A aqueous uranium (U aq ) plume exists at the site that discharges to the Columbia River. We performed temporally intensive monitoring of specific conductance (SpC) and U aq concentrations within the plume at increasing distances from the river shoreline and at locations both within and outside a presumed subsurface pathway of lateral hydrologic exchange. Specific conductance and U aq were utilized as in-situ tracers of hydrologic exchange and associated groundwater-surface water mixing. Monitoring data was collected for an approximate two year period during which seasonal river elevations varied by more than 2 m causing distinct events of river water intrusion and retreat from the near-shore aquifer, and large temporal variations in well-water analyte concentrations. Statistical clustering analysis was applied to identify well groupings exhibiting common time series behaviors. PFLOTRAN was implemented within a modeling domain of 3.3 x 1.9 km with dynamic flow boundary conditions and accurate subsurface topography for the basal aquitard underlying the aquifer. In-silico calculations of river water tracer behavior were performed to assess predominant flow directions and velocities influenced by river stage variations and variable aquitard topography, and their collective influences on SpC and U aq concentrations. We demonstrate complicated multidirectional flow behaviors at the kilometer scale resulting from hydrologic exchange and subsurface structure that influence plume dynamics and determine the overall system behavior.

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 OSTI2020

Comparison of equilibrium climate sensitivity estimates from slab ocean, 150-year, and longer simulations

We compare equilibrium climate sensitivity (ECS) estimates from pairs of long (= 800-year) control and abruptly quadrupled CO2 simulations with shorter (150, 300 year) coupled atmosphere-ocean simulations and Slab Ocean Models (SOM). Consistent with previous work, ECS estimates from shorter coupled simulations based on annual averages for years 1-150 underestimate those from SOM (-8% ± 13%) and long (-14% ± 8%) simulations. Analysis of only years 21-150 improved agreement with SOM (-2% ± 14%) and long (-8% ± 10%) estimates. Use of pentadal averages for years 51-150 results in improved agreement with long simulations (-4% ± 11%). While ECS estimates from current generation US models based on SOM and coupled annual averages of years 1-150 range from 2.6°C to 5.3°C, estimates based longer simulations of the same models range from 3.2°C to 7.0°C. Such variations between methods argues for caution in comparison and interpretation of ECS estimates across models.

Dunne, John P.↗