Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model Coupling”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Patterns of Precipitation Change and Climatological Uncertainty among CMIP5 Models, with a Focus on the Midlatitude Pacific Storm Track

Projections of modeled precipitation ( P ) change in global warming scenarios demonstrate marked intermodel disagreement at regional scales. Empirical orthogonal functions (EOFs) and maximum covariance analysis (MCA) are used to diagnose spatial patterns of disagreement in the simulated climatology and end-of-century P changes in phase 5 of the Coupled Model Intercomparison Project (CMIP5) archive. The term principal uncertainty pattern (PUP) is used for any robust mode calculated when applying these techniques to a multimodel ensemble. For selected domains in the tropics, leading PUPs highlight features at the margins of convection zones and in the Pacific cold tongue. The midlatitude Pacific storm track is emphasized given its relevance to wintertime P projections over western North America. The first storm-track PUP identifies a sensitive region of disagreement in P increases over the eastern midlatitude Pacific where the storm track terminates, related to uncertainty in an eastward extension of the climatological jet. The second PUP portrays uncertainty in a zonally asymmetric meridional shift of storm-track P , related to uncertainty in the extent of a poleward jet shift in the western Pacific. Both of these modes appear to arise primarily from intermodel differences in the response to radiative forcing, distinct from sampling of internal variability. The leading storm-track PUPs for P and zonal wind change exhibit similarities to the leading uncertainty patterns for the historical climatology, indicating important and parallel sensitivities in the eastern Pacific storm-track terminus region. However, expansion coefficients for climatological uncertainties tend to be weakly correlated with those for end-of-century change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2013

Effective field theory of weakly coupled inflationary models

The application of Effective Field Theory (EFT) methods to inflation has taken a central role in our current understanding of the very early universe. The EFT perspective has been particularly useful in analyzing the self-interactions determining the evolution of co-moving curvature perturbations (Goldstone boson modes) and their influence on low-energy observables. However, the standard EFT formalism, to lowest order in spacetime differential operators, does not provide the most general parametrization of a theory that remains weakly coupled throughout the entire low-energy regime. Here we study the EFT formulation by including spacetime differential operators implying a scale dependence of the Goldstone boson self-interactions and its dispersion relation. These operators are shown to arise naturally from the low-energy interaction of the Goldstone boson with heavy fields that have been integrated out. We find that the EFT then stays weakly coupled all the way up to the cutoff scale at which ultraviolet degrees of freedom become operative. This opens up a regime of new physics where the dispersion relation is dominated by a quadratic dependence on the momentum ω ∼ p{sup 2}. In addition, provided that modes crossed the Hubble scale within this energy range, the predictions of inflationary observables — including non-Gaussian signatures — are significantly affected by the new scales characterizing it.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Modeling optical coupling of plasmons and inhomogeneously broadened emitters

Optically coupling quantum emitters to nanoparticles provides the foundation for many plasmonic applications. Including quantum mechanical effects within the calculations can be crucial for designing new devices, but classical approximations are sometimes sufficient. Comprehending how the classical and quantum mechanical descriptions of quantum emitters alter their calculated optical response will lead to a better understanding of how to design devices. In this work, we describe how the semiclassical Maxwell-Liouville method can be used to calculate the optical response from inhomogeneously broadened states. After describing the Maxwell-Liouville algorithm, we use the method to study the photon echoes from quantum dots and compare the results against analytical models. We then modify the quantum dot’s state distribution to match a PbS 850 nm quantum dot’s absorption spectra to see how the complete quasi-band structure affects their coupling to gold nanoislands. Finally, we compare the results with previously published work to demonstrate where the complete quantum dot description is necessary.

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

Dynamics of Johari-Goldstein β relaxation and its universal relation to α relaxation in bulk metallic glasses by mechanical spectroscopy

The dynamics of the Johari–Goldstein (JG) β relaxation and the α relaxation in bulk metallic glasses (MGs) has been investigated by using mechanical spectroscopy combined with the Coupling Model. The β relaxations of MGs exhibit different behaviors such as peaks, humps, and excess wings due to the different fluctuations of the chemical interactions among the constituting atoms. A universal correlation between the β relaxation and the α relaxation is generally found by their activation energies and relaxation times as well as the non-exponentiality parameter of the α relaxation, which can be predicted quantitatively from the Coupling Model. Based on the quasi-point defects theory, a correlation factor χ shows a broad peak along with the β relaxation, suggesting that the concentration and the correlation degree of the string-like configurations involved in the β relaxation vary with increasing temperature, which challenges the previous view that the system is in an iso-configuration state below T{sub g} and may shed new light on the nature of the JG β relaxation in metallic glasses.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Origin of Large Effective Phonon Magnetic Moments in Monolayer MoS 2

Recent helicity-resolved magneto-Raman spectroscopy measurement demonstrates large effective phonon magnetic moments of ~2.5 μ B in monolayer MoS 2 , highlighting resonant excitation of bright excitons as a feasible route to activate Γ-point circularly polarized phonons in transition metal dichalcogenides. However, a microscopic picture of this intriguing phenomenon remains lacking. In this work, we show that an orbital transition between the split conduction bands (Δ 0 = 4 meV) of MoS 2 couples to the doubly degenerate E" phonon mode (Ω 0 = 33 meV), forming two hybridized states. Our phononic and electronic Raman scattering measurements capture these two states: (i) one with predominantly phonon contribution in the helicity-switched channels and (ii) one with primarily orbital contribution in the helicity-conserved channels. An orbital-phonon coupling model successfully reproduces the large effective magnetic moments of the circularly polarized phonons and explains their thermodynamic properties. Strikingly, the Raman mode from the orbital transition is superimposed on a strong quasi-elastic scattering background, indicating the presence of spin fluctuations. As a result, the electrons excited to the conduction bands through the exciton exhibit paramagnetic behavior although MoS 2 is generally considered as a nonmagnetic material. By depositing nanometer-thickness nickel thin films on monolayer MoS 2 , we tune the electronic structure so that the A exciton perfectly overlaps with the 633 nm laser. The optimization of resonance excitation leads to pronounced tunability of the orbital-phonon hybridized states. Our results generalize the orbital-phonon coupling model of effective phonon magnetic moments to material systems beyond the paramagnets and magnets.

2D transition metal dichalcogenide↗
DOE OSTI2014

Right-handed neutrino production rate at T>160 GeV

The production rate of right-handed neutrinos from a Standard Model plasma at a temperature above a hundred GeV has previously been evaluated up to NLO in Standard Model couplings (g∼2/3) in relativistic (M∼πT) and non-relativistic regimes (M≫πT), and up to LO in an ultrarelativistic regime (M≲gT). The last result necessitates an all-orders resummation of the loop expansion, accounting for multiple soft scatterings of the nearly light-like particles participating in 1↔2 reactions. In this paper we suggest how the regimes can be interpolated into a result applicable for any right-handed neutrino mass and at all temperatures above 160 GeV. The results can also be used for determining the lepton number washout rate in models containing right-handed neutrinos. Numerical results are given in a tabulated form permitting for their incorporation into leptogenesis codes. We note that due to effects from soft Higgs bosons there is a narrow intermediate regime around M∼g{sup 1/2}T in which our interpolation is phenomenological and a more precise study would be welcome.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Energetics and the Role of Defects in Fe(II)-Catalyzed Goethite Recrystallization from Molecular Simulations

Goethite is one of the most stable and common iron (III) minerals at the Earth’s near surface. However, recent isotope- tracer studies suggest that goethite continuously recrystallizes in the presence of aqueous Fe(II) ions. Such studies often indicate the presence of two regimes of atom exchange kinetics, a rapid stage assigned to reactive defect sites initially available at particle surfaces, followed by slower continuous exchange. An autocatalytic solid-state electron conduction model coupling Fe(II) oxidative adsorption to its reductive release at spatially distinct sites has been proposed, but the thermodynamic driving force has yet to be pinpointed. Here, using a novel hybrid/reactive molecular simulation method, for goethite (110) surfaces at circumneutral pH, we rigorously tested whether surface free energy minimization, including examining the role of structural defects, is sufficient to overcome the activation energy for interfacial electron transfer and conduction. The simulations quantitatively show that: (i). on smooth stable surfaces the available thermal energy at dynamic equilibrium is sufficient to sustain the slow continuous regime of atom exchange kinetics via short intra-surface electron conduction pathways of 1-2 nm (3-5 Fe site hops); (ii). in this slower regime, the model converges to atom exchange kinetics of 10 -5 Fe s -1 cm -2 , a rate recently deduced from stochastic modeling of experimental data and linked to the reductive dissolution rate of goethite; (iii). the driving force for initially rough defective goethite surfaces to smoothen can accelerate atom exchange to an extent quantitatively consistent with that observed in the initial rapid stage, in this case accessing conduction pathways of up to 8 nm. The findings suggest that the interaction of Fe(II) with initially defective goethite surfaces can drive, by the conduction model, atom exchange that is capable of recrystallizing the interiors of nanoscale particles, and that, closer to equilibrium on smooth surfaces, slower atom exchange continues in perpetuity but likely involving only the outermost atomic layers.

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

Quantification of marine aerosol subgrid variability and its correlation with clouds based on high‐resolution regional modeling

Abstract One limitation of most global climate models (GCMs) is that with the horizontal resolutions they typically employ, they cannot resolve the subgrid variability (SGV) of clouds and aerosols, adding extra uncertainties to the aerosol radiative forcing estimation. To inform the development of an aerosol subgrid variability parameterization, here we analyze the aerosol SGV over the southern Pacific Ocean simulated by the high‐resolution Weather Research and Forecasting model coupled to Chemistry. We find that within a typical GCM grid, the aerosol mass subgrid standard deviation is 15% of the grid‐box mean mass near the surface on a 1 month mean basis. The fraction can increase to 50% in the free troposphere. The relationships between the sea‐salt mass concentration, meteorological variables, and sea‐salt emission rate are investigated in both the clear and cloudy portion. Under clear‐sky conditions, marine aerosol subgrid standard deviation is highly correlated with the standard deviations of vertical velocity, cloud water mixing ratio, and sea‐salt emission rates near the surface. It is also strongly connected to the grid box mean aerosol in the free troposphere (between 2 km and 4 km). In the cloudy area, interstitial sea‐salt aerosol mass concentrations are smaller, but higher correlation is found between the subgrid standard deviations of aerosol mass and vertical velocity. Additionally, we find that decreasing the model grid resolution can reduce the marine aerosol SGV but strengthen the correlations between the aerosol SGV and the total water mixing ratio (sum of water vapor, cloud liquid, and cloud ice mixing ratios).

Lin, Guangxing↗
DOE OSTI2024

The DECOVALEX international collaboration on modeling of coupled subsurface processes and its contribution to confidence building in radioactive waste disposal

Abstract The long-lived radiotoxicity of the high-level radioactive waste generated by nuclear power plants requires safe isolation from the biosphere for many hundreds of thousands of years. An international consensus has emerged that such isolation can best be provided by disposal in mined geologic repositories, a strategy that today is pursued by most countries dealing with radioactive waste. However, the need to predict the performance of such repositories over very long time periods generates large uncertainties that have to be accounted for in safety assessments. The findings from such safety assessments need to be conveyed to all stakeholders in a clear way, such that public confidence in geologic disposal solutions can be achieved. It is suggested here that close international collaboration on the technical aspects of geologic waste disposal has helped, and will continue to help, building trust and increasing confidence. This paper discusses a particular international collaboration initiative referred to as DECOVALEX, which brings together multiple teams and disciplines to collectively tackle complex experimental and modeling challenges related to geologic disposal. By describing how DECOVALEX works and by providing joint research examples, a case is made that such international collaboration contributes to knowledge transfer and confidence building in radioactive waste disposal science.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Reducing uncertainty of polar to midlatitude linkages using DOE’s E3SM in a coordinated model-experiment setting

This project brought DOE’s climate modeling effort with the Energy Exascale Earth System Model into the Polar Amplification Model Intercomparison Project (PAMIP), which is part of the sixth and latest Coupled Model Intercomparison Project, CMIP6. PAMIP examines the causes and consequences of polar amplification, when external forcing results in a larger temperature increase in high latitudes than the global average, in a coordinated set of model experiments with a common modeling protocol. Our teams from UC Irvine and the University of Toronto have designed, carried out, analyzed, and disseminated PAMIP output from the Energy Exascale Earth System Model (E3SM) and the Community Earth System Model (CESM). PAMIP’s ongoing stream of significant new results have advanced progress in the community’s understanding and led to new outstanding research questions that have motivated further work. PAMIP has led to improved consensus on the atmospheric response to sea ice loss. The important finding is that for a similar sea ice anomaly forcing, the simulated atmospheric response in the troposphere is remarkably consistent among the 16 models’ runs analyzed. The zonal-mean tropospheric response consists of a very robust equatorward shift of the westerly flow in mid-latitudes. However, while the multi-model mean response is robust, it has a weak amplitude relative to internal variability. We identified a weakness in the models (including E3SM) in terms of their response to sea-ice forcing that is related to eddy forcing (or nonlinear dynamical effects) at mid- to high latitudes. In fact, E3SM is an outlier in terms of the models participating, and in that sense, it turned out to be a vital participant model. We found that reductions in energy transport due sea-ice loss and involving dry air only are compensated by increases in moist energy transport from warmer sea surface temperature in midlatitudes. This leads to a large spread in energy transport into the Arctic and is a potential source of spread in Arctic amplification. We identified an important role that climate modes, including tropical modes of variability (El Nino and the Southern Oscillation (ENSO); the Quasi-biennial Oscillation (QBO)) play in the response to sea-ice anomalies, including in ocean coupled experiments. Similarly, we identified and quantified the contribution of sea-ice thickness to the atmospheric response compared to the response to sea-ice extent only. We found that it is important to run large ensembles and even with an ensemble size of 100 simulations the response is largely influenced by internal variability. We demonstrated convincingly that Ural blocking, not sea-ice loss, provides the weakening of the stratospheric polar vortex in fall/early winter and a negative phase of the North Atlantic Oscillation that can last for up to two months. However, sea-ice anomalies can influence the background flow so that the response to Ural blocking is more persistent under low sea-ice conditions in the Barents/Kara Sea than high sea ice Atmospheric model hierarchies that progressively add individual processes have a long history in providing dynamical insight for modeling the atmosphere. Similarly, coupled model hierarchies that progressively add individual ocean processes can provide insights into the workings of the coupled climate system, however such hierarchies have not been available except for a non-dynamical slab ocean model. Because of the missing processes, surface flux corrections must be added to produce a target climate. In this project, we managed to overcome this problem and develop a globally coupled ocean model hierarchy in CESM that can turn on and off the processes of mixed-layer entrainment and Ekman flow. We used the hierarchy to study the impact of Arctic sea-ice loss on the climate system. We find that the effect of mixed-layer entrainment on ocean heat uptake influences the atmospheric circulation by shifting the latitudinal positions of the mid-latitude westerly jet and the Intertropical Convergence Zone (ITCZ). In quadrupled CO 2 experiments, we studied how air-sea coupling affects the response of tropical rainfall under global warming. In order to identify the importance of individual ocean processes, we used the hierarchy of ocean models to separate the effects of seasonal mixed-layer entrainment, wind-driven Ekman flows, and frictional flows. We showed that including Ekman and frictional flows allows our simulation to produce the Pacific Ocean's enhanced equatorial warming pattern and equatorward ITCZ contraction noted in previous climate simulations. We also showed that the frictional flow, which has yet to receive much attention, is as important as the Ekman flow in generating equatorial heat convergence.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Modeling the Coupled Chemo-Thermo-Mechanical Behavior of Amorphous Polymer Networks.

Amorphous polymers exhibit a rich landscape of time-dependent behavior including viscoelasticity, structural relaxation, and viscoplasticity. These time-dependent mechanisms can be exploited to achieve shape-memory behavior, which allows the material to store a programmed deformed shape indefinitely and to recover entirely the undeformed shape in response to specific environmental stimulus. The shape-memory performance of amorphous polymers depends on the coordination of multiple physical mechanisms, and considerable opportunities exist to tailor the polymer structure and shape-memory programming procedure to achieve the desired performance. The goal of this project was to use a combination of theoretical, numerical and experimental methods to investigate the effect of shape memory programming, thermo-mechanical properties, and physical and environmental aging on the shape memory performance. Physical and environmental aging occurs during storage and through exposure to solvents, such as water, and can significantly alter the viscoelastic behavior and shape memory behavior of amorphous polymers. This project – executed primarily by Professor Thao Nguyen and Graduate Student Rui Xiao at Johns Hopkins University in support of a DOE/NNSA Presidential Early Career Award in Science and Engineering (PECASE) – developed a theoretical framework for chemothermo- mechanical behavior of amorphous polymers to model the effects of physical aging and solvent-induced environmental factors on their thermoviscoelastic behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The fully coupled regionally refined model of E3SM version 2: overview of the atmosphere, land, and river results

Abstract. This paper provides an overview of the United States (US) Department of Energy's (DOE's) Energy Exascale Earth System Model version 2 (E3SMv2) fully coupled regionally refined model (RRM) and documents the overall atmosphere, land, and river results from the Coupled Model Intercomparison Project 6 (CMIP6) DECK (Diagnosis, Evaluation, and Characterization of Klima) and historical simulations – a first-of-its-kind set of climate production simulations using RRM. The North American (NA) RRM (NARRM) is developed as the high-resolution configuration of E3SMv2 with the primary goal of more explicitly addressing DOE's mission needs regarding impacts to the US energy sector facing Earth system changes. The NARRM features finer horizontal resolution grids centered over NA, consisting of 25→100 km atmosphere and land, a 0.125∘ river-routing model, and 14→60 km ocean and sea ice. By design, the computational cost of NARRM is ∼3× of the uniform low-resolution (LR) model at 100 km but only ∼ 10 %–20 % of a globally uniform high-resolution model at 25 km. A novel hybrid time step strategy for the atmosphere is key for NARRM to achieve improved climate simulation fidelity within the high-resolution patch without sacrificing the overall global performance. The global climate, including climatology, time series, sensitivity, and feedback, is confirmed to be largely identical between NARRM and LR as quantified with typical climate metrics. Over the refined NA area, NARRM is generally superior to LR, including for precipitation and clouds over the contiguous US (CONUS), summertime marine stratocumulus clouds off the coast of California, liquid and ice phase clouds near the North Pole region, extratropical cyclones, and spatial variability in land hydrological processes. The improvements over land are related to the better-resolved topography in NARRM, whereas those over ocean are attributable to the improved air–sea interactions with finer grids for both atmosphere and ocean and sea ice. Some features appear insensitive to the resolution change analyzed here, for instance the diurnal propagation of organized mesoscale convective systems over CONUS and the warm-season land–atmosphere coupling at the southern Great Plains. In summary, our study presents a realistically efficient approach to leverage the fully coupled RRM framework for a standard Earth system model release and high-resolution climate production simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Understanding the changes in ductility and Poisson's ratio of metallic glasses during annealing from microscopic dynamics

In the paper K. L. Ngai et al., [J. Chem. 140, 044511 (2014)], the empirical correlation of ductility with the Poisson's ratio, ν{sub Poisson}, found in metallic glasses was theoretically explained by microscopic dynamic processes which link on the one hand ductility, and on the other hand the Poisson's ratio. Specifically, the dynamic processes are the primitive relaxation in the Coupling Model which is the precursor of the Johari–Goldstein β-relaxation, and the caged atoms dynamics characterized by the effective Debye–Waller factor f{sub 0} or equivalently the nearly constant loss (NCL) in susceptibility. All these processes and the parameters characterizing them are accessible experimentally except f{sub 0} or the NCL of caged atoms; thus, so far, the experimental verification of the explanation of the correlation between ductility and Poisson's ratio is incomplete. In the experimental part of this paper, we report dynamic mechanical measurement of the NCL of the metallic glass La{sub 60}Ni{sub 15}Al{sub 25} as-cast, and the changes by annealing at temperature below T{sub g}. The observed monotonic decrease of the NCL with aging time, reflecting the corresponding increase of f{sub 0}, correlates with the decrease of ν{sub Poisson}. This is important observation because such measurements, not made before, provide the missing link in confirming by experiment the explanation of the correlation of ductility with ν{sub Poisson}. On aging the metallic glass, also observed in the isochronal loss spectra is the shift of the β-relaxation to higher temperatures and reduction of the relaxation strength. These concomitant changes of the β-relaxation and NCL are the root cause of embrittlement by aging the metallic glass. The NCL of caged atoms is terminated by the onset of the primitive relaxation in the Coupling Model, which is generally supported by experiments. From this relation, the monotonic decrease of the NCL with aging time is caused by the slowing down of the primitive relaxation and β-relaxation on annealing, and vice versa.

36 MATERIALS SCIENCE↗
DOE OSTI2013

DeepCwind Consortium Research Program: January 15, 2010 - March 31, 2013 (Final Technical Report)

This is the final technical report for the U.S. Department of Energy-funded program, DE-0002981: DeepCwind Consortium Research Program. The project objective was the partial validation of coupled models and optimization of materials for offshore wind structures. The United States has a great opportunity to harness an indigenous abundant renewable energy resource: offshore wind. In 2010, the National Renewable Energy Laboratory (NREL) estimated there to be over 4,000 GW of potential offshore wind energy found within 50 nautical miles of the US coastlines (Musial and Ram, 2010). The US Energy Information Administration reported the total annual US electric energy generation in 2010 was 4,120 billion kilowatt-hours (equivalent to 470 GW) (US EIA, 2011), slightly more than 10% of the potential offshore wind resource. In addition, deep water offshore wind is the dominant US ocean energy resource available comprising 75% of the total assessed ocean energy resource as compared to wave and tidal resources (Musial, 2008). Through these assessments it is clear offshore wind can be a major contributor to US energy supplies. The caveat to capturing offshore wind along many parts of the US coast is deep water. Nearly 60%, or 2,450 GW, of the estimated US offshore wind resource is located in water depths of 60 m or more (Musial and Ram, 2010). At water depths over 60 m building fixed offshore wind turbine foundations, such as those found in Europe, is likely economically infeasible (Musial et al., 2006). Therefore floating wind turbine technology is seen as the best option for extracting a majority of the US offshore wind energy resource. Volume 1 - Test Site; Volume 2 - Coupled Models; and Volume 3 - Composite Materials.

17 WIND ENERGY↗
DOE OSTI2024

Computational modeling of coupled mechanical damage and electrochemistry in ternary oxide composite electrodes

Performance degradation of ternary layered oxide cathodes largely originates from their loss of structural integrity in cyclic usage. Mechanical damage, such as intergranular fracture of the active particles, is not only a mechanical cleavage process but also interferes with electrochemical kinetics such as infiltration of liquid electrolyte, surface corrosion of the constituent primary particles, and may eventually isolate the primary grains from the electron conducting network. Here, in this work, we develop a computational framework that integrates electrochemistry of a LiNi x Mn y Co 1−x−y O 2 (NMC) composite cathode with mechanical damage of the active particles. To fully examine the intricate chemomechanical behavior of the electrode, we evaluate the effects of the anisotropic material properties, the influence of mechanical potential on Li transport, and the concurrent intergranular fracture and electrolyte penetration along the grain boundaries upon multiple cycles. Electrolyte infiltration benefits capacity retention but aggravates further mechanical damage by corrosion. Structural failure mostly occurs in the first charging due to the anisotropic mechanical strain between the primary grains, while the resulting damage remains stable in the later few cycles. The results are consistent with experimental observations and the integration of electrochemistry and mechanical failure enables a step further understanding of the complex mechanism of battery degradation.

Battery degradation↗
DOE OSTI2013

Subsurface Uranium Fate and Transport: Integrated Experiments and Modeling of Coupled Biogeochemical Mechanisms of Nanocrystalline Uraninite Oxidation by Fe(III)-(hydr)oxides - Project Final Report

Subsurface bacteria including sulfate reducing bacteria (SRB) reduce soluble U(VI) to insoluble U(IV) with subsequent precipitation of UO 2 . We have shown that SRB reduce U(VI) to nanometer-sized UO 2 particles (1-5 nm) which are both intra- and extracellular, with UO 2 inside the cell likely physically shielded from subsequent oxidation processes. We evaluated the UO 2 nanoparticles produced by Desulfovibrio desulfuricans G20 under growth and non-growth conditions in the presence of lactate or pyruvate and sulfate, thiosulfate, or fumarate, using ultrafiltration and HR-TEM. Results showed that a significant mass fraction of bioreduced U (35-60%) existed as a mobile phase when the initial concentration of U(VI) was 160 µM. Further experiments with different initial U(VI) concentrations (25 - 900 M) in MTM with PIPES or bicarbonate buffers indicated that aggregation of uraninite depended on the initial concentrations of U(VI) and type of buffer. It is known that under some conditions SRB-mediated UO 2 nanocrystals can be reoxidized (and thus remobilized) by Fe(III)-(hydr)oxides, common constituents of soils and sediments. To elucidate the mechanism of UO 2 reoxidation by Fe(III) (hydr)oxides, we studied the impact of Fe and U chelating compounds (citrate, NTA, and EDTA) on reoxidation rates. Experiments were conducted in anaerobic batch systems in PIPES buffer. Results showed EDTA significantly accelerated UO 2 reoxidation with an initial rate of 9.5 M day-1 for ferrihydrite. In all cases, bicarbonate increased the rate and extent of UO 2 reoxidation with ferrihydrite. The highest rate of UO 2 reoxidation occurred when the chelator promoted UO 2 and Fe(III) (hydr)oxide dissolution as demonstrated with EDTA. When UO 2 dissolution did not occur, UO 2 reoxidation likely proceeded through an aqueous Fe(III) intermediate as observed for both NTA and citrate. To complement to these laboratory studies, we collected U-bearing samples from a surface seep at the Rifle field site and have measured elevated U concentrations in oxic iron-rich sediments. To translate experimental results into numerical analysis of U fate and transport, a reaction network was developed based on Sani et al. (2004) to simulate U(VI) bioreduction with concomitant UO 2 reoxidation in the presence of hematite or ferrihydrite. The reduction phase considers SRB reduction (using lactate) with the reductive dissolution of Fe(III) solids, which is set to be microbially mediated as well as abiotically driven by sulfide. Model results show the oxidation of HS– by Fe(III) directly competes with UO 2 reoxidation as Fe(III) oxidizes HS– preferentially over UO 2 . The majority of Fe reduction is predicted to be abiotic, with ferrihydrite becoming fully consumed by reaction with sulfide. Predicted total dissolved carbonate concentrations from the degradation of lactate are elevated (log(pCO 2 ) ~ –1) and, in the hematite system, yield close to two orders-of-magnitude higher U(VI) concentrations than under initial carbonate concentrations of 3 mM. Modeling of U(VI) bioreduction with concomitant reoxidation of UO 2 in the presence of ferrihydrite was also extended to a two-dimensional field-scale groundwater flow and biogeochemically reactive transport model for the South Oyster site in eastern Virginia. This model was developed to simulate the field-scale immobilization and subsequent reoxidation of U by a biologically mediated reaction network.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Trans-Pacific transport and evolution of aerosols: Evaluation of quasi-global WRF-Chem simulation with multiple observations

A fully coupled meteorology-chemistry model (WRF-Chem, the Weather Research and Forecasting model coupled with chemistry) has been configured to conduct quasi-global simulation for 5 years (2010–2014) and evaluated with multiple observation data sets for the first time. The evaluation focuses on the simulation over the trans-Pacific transport region using various reanalysis and observational data sets for meteorological fields and aerosol properties. The simulation generally captures the overall spatial and seasonal variability of satellite retrieved aerosol optical depth (AOD) and absorbing AOD (AAOD) over the Pacific that is determined by the outflow of pollutants and dust and the emissions of marine aerosols. The assessment of simulated extinction Ångström exponent (EAE) indicates that the model generally reproduces the variability of aerosol size distributions as seen by satellites. In addition, the vertical profile of aerosol extinction and its seasonality over the Pacific are also well simulated. The difference between the simulation and satellite retrievals can be mainly attributed to model biases in estimating marine aerosol emissions as well as the satellite sampling and retrieval uncertainties. Compared with the surface measurements over the western USA, the model reasonably simulates the observed magnitude and seasonality of dust, sulfate, and nitrate surface concentrations, but significantly underestimates the peak surface concentrations of carbonaceous aerosol likely due to model biases in the spatial and temporal variability of biomass burning emissions and secondary organic aerosol (SOA) production. A sensitivity simulation shows that the trans-Pacific transported dust, sulfate, and nitrate can make significant contribution to surface concentrations over the rural areas of the western USA, while the peaks of carbonaceous aerosol surface concentrations are dominated by the North American emissions. Both the retrievals and simulation show small interannual variability of aerosol characteristics for 2010–2014 averaged over three Pacific sub-regions. Furthermore, the evaluation in this study demonstrates that the WRF-Chem quasi-global simulation can be used for investigating trans-Pacific transport of aerosols and providing reasonable inflow chemical boundaries for the western USA, allowing one to further understand the impact of transported pollutants on the regional air quality and climate with high-resolution nested regional modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Interannual to decadal climate variability of sea salt aerosols in the coupled climate model CESM1.0

This study reported herein examines multi-year climate variability associated with sea salt aerosols and their contribution to the variability of shortwave cloud forcing (SWCF) using a 150-year simulation for pre-industrial conditions of the Community Earth System Model version 1.0 (CESM1). The results suggest that changes in sea salt and related cloud and radiative properties on interannual timescales are dominated by the ENSO cycle. Sea salt variability on longer (interdecadal) timescales is associated with low-frequency Pacific ocean variability similar to the interdecadal Pacific Oscillation (IPO), but does not show a statistically significant spectral peak. A multivariate regression suggests that sea salt aerosol variability may contribute to SWCF variability in the tropical Pacific, explaining up to 25-35% of the variance in that region. Elsewhere, there is only a small aerosol influence on SWCF through modifying cloud droplet number and liquid water path that contributes to the change of cloud effective radius and cloud optical depth (and hence cloud albedo), producing a multi-year aerosol-cloud-wind interaction.

54 ENVIRONMENTAL SCIENCES↗