Search NASA⌕ Search

SEARCH · Search NASA

Results for “Earth’s accretion”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

The Non-carbonaceous–Carbonaceous Meteorite Dichotomy

The isotopic dichotomy between non-carbonaceous (NC) and carbonaceous (CC) meteorites indicates that meteorite parent bodies derive from two genetically distinct reservoirs, which presumably were located inside (NC) and outside (CC) the orbit of Jupiter and remained isolated from each other for the first few million years of the solar system. Here we review the discovery of the NC–CC dichotomy and its implications for understanding the early history of the solar system, including the formation of Jupiter, the dynamics of terrestrial planet formation, and the origin and nature of Earth’s building blocks. The isotopic difference between the NC and CC reservoirs is probably inherited from the solar system’s parental molecular cloud and has been maintained through the rapid formation of Jupiter that prevented significant exchange of material from inside (NC) and outside (CC) its orbit. The growth and/or migration of Jupiter resulted in inward scattering of CC bodies, which accounts for the co-occurrence of NC and CC bodies in the present-day asteroid belt and the delivery of presumably volatile-rich CC bodies to the growing terrestrial planets. Earth’s primitive mantle, at least for siderophile elements like Mo, has a mixed NC–CC composition, indicating that Earth accreted CC bodies during the final stages of its growth, perhaps through the Moon-forming giant impactor. The late-stage accretion of CC bodies to Earth is sufficient to account for the entire budget of Earth’s water and highly volatile species.

58 GEOSCIENCES↗
DOE OSTI2022

The Lithophile Element Budget of Earth's Core

Abstract The relative composition of Earth's core and mantle were set during core formation. By determining how elements partition between metal and silicate at high pressures and temperatures, measurements of the mantle composition and geophysical observations of the core can be used to understand the mechanisms by which Earth formed. Here we present the results of metal‐silicate partitioning experiments for a range of nominally lithophile elements (Al, Ca, K, Mg, O, Si, Th, and U) and S to 85 GPa and up to 5400 K. With our results and a compilation of literature data, we developed a parameterization for partitioning that accounts for compositional dependencies in both the metal and silicate phases. Using this parameterization in a range of planetary growth models, we find that, in general, lithophile element partitioning into the metallic phase is enhanced at high temperatures. The relative abundances of FeO, SiO 2 , and MgO in the mantle vary significantly between planetary growth models, and the mantle abundances of these elements can be used to provide important constraints on Earth's accretion. To match Earth's core mass and mantle composition, Earth's building blocks must have been enriched in Fe and depleted in Si compared with CI chondrites. Finally, too little Mg, Si, and O are partitioned into the core for precipitation of oxides to be a major source of energy for the geodynamo. In contrast, several ppb of U can be partitioned into the core at high temperatures, and this energy source must be accounted for in thermal evolution models.

Chidester, B. A.↗
DOE OSTI2022

Did Earth Eat Its Leftovers? Impact Ejecta as a Component of the Late Veneer

The presence of highly siderophile elements in Earth's mantle indicates that a small percentage of Earth's mass was delivered after the last giant impact in a stage of "late accretion." There is ongoing debate about the nature of late-accreted material and the sizes of late-accreted bodies. Earth appears isotopically most similar to enstatite chondrites and achondrites. It has been suggested that late accretion must have been dominated by enstatite-like bodies that originated in the inner disk, rather than ordinary or carbonaceous chondrites. Here we examine the provenances of "leftover" planetesimals present in the inner disk in the late stages of accretion simulations. Dynamically excited planet formation produces planets and embryos with similar provenances, suggesting that the Moon-forming impactor may have had a stable isotope composition very similar to the proto-Earth. Commonly, some planetesimal-sized bodies with similar provenances to the Earth-like planets are left at the end of the main stage of growth. The most chemically similar planetesimals are typically fragments of protoplanets ejected millions of years earlier. If these similar-provenance bodies are later accreted by the planet, they will represent late-accreted mass that naturally matches Earth's composition. The planetesimal-sized bodies that exist during the giant impact phase can have large core mass fractions, with core provenances similar to the proto-Earth. These bodies are an important potential source for highly siderophile elements. The range of core fractions in leftover planetesimals complicates simple inferences as to the mass and origin of late accretion based on the highly siderophile elements in the mantle.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

No 182W evidence for early Moon formation

The Moon-forming giant impact was probably the last major event in Earth’s accretion, so dating this event is critical to determine the timeline of terrestrial planet formation. Recently, Thiemens et al. used the short-lived 182 Hf– 182 W system to argue that the Moon formed within the first 60 Myr of Solar System history. Here we demonstrate, however, that mixing processes during and after the giant impact modified the 182W compositions of the Earth and Moon, which hampers the use of the Hf–W system to date the Moon. Our results show that the lunar 182 W record is fully consistent with a recently determined, younger age of the Moon of 142 ± 25 Myr after Solar System formation.

58 GEOSCIENCES↗
DOE OSTI2022

A dry ancient plume mantle from noble gas isotopes

Primordial volatiles were delivered to terrestrial reservoirs during Earth’s accretion, and the mantle plume source is thought to have retained a greater proportion of primordial volatiles compared with the upper mantle. This study shows that mantle He, Ne, and Xe isotopes require that the plume mantle had low concentrations of volatiles like Xe and H 2 O at the end of accretion compared with the upper mantle. A lower extent of mantle processing alone is not sufficient to explain plume noble gas signatures. Ratios of primordial isotopes are used to determine proportions of solar, chondritic, and regassed atmospheric volatiles in the plume mantle and upper mantle. The regassed Ne flux exceeds the regassed Xe flux but has a small impact on the mantle Ne budget. Pairing primordial isotopes with radiogenic systems gives an absolute concentration of 130 Xe in the plume source of ∼1.5 × 10 7 atoms 130 Xe/g at the end of accretion, ∼4 times less than that determined for the ancient upper mantle. A record of limited accretion of volatile-rich solids thus survives in the He-Ne-Xe signatures of mantle rocks today. A primordial viscosity contrast originating from a factor of ∼4 to ∼250 times lower H 2 O concentration in the plume mantle compared with the upper mantle may explain (a) why giant impacts that triggered whole mantle magma oceans did not homogenize the growing planet, (b) why the plume mantle has experienced less processing by partial melting over Earth’s history, and (c) how early-formed isotopic heterogeneities may have survived ∼4.5 Gy of solid-state mantle convection.

58 GEOSCIENCES↗
DOE OSTI2021

Wide Dust Gaps in Protoplanetary Disks Induced by Eccentric Planets: A Mass-eccentricity Degeneracy

The tidal perturbation of embedded protoplanets on their natal disks has been widely attributed to be the cause of gap-ring structures in submillimeter images of protoplanetary disks around T Tauri stars. Numerical simulations of this process have been used to propose scaling of characteristic dust-gap width/gap-ring distance with respect to planet mass. Applying such scaling to analyze observed gap samples yields a continuous mass distribution for a rich population of hypothetical planets in the range of several Earth to Jupiter masses. In contrast, the conventional core-accretion scenario of planet formation predicts a bimodal mass function due to (1) the onset of runaway gas accretion above ~20 Earth masses and (2) suppression of accretion induced by gap opening. Here, we examine the dust disk response to the tidal perturbation of eccentric planets as a possible resolution of this paradox. Based on simulated gas and dust distributions, we show the gap-ring separation of Neptune-mass planets with small eccentricities might become comparable to that induced by Saturn-mass planets on circular orbits. Furthermore, this degeneracy may obliterate the discrepancy between the theoretical bimodal mass distribution and the observed continuous gap width distribution. Despite damping due to planet–disk interaction, modest eccentricity may be sustained either in the outer regions of relatively thick disks or through resonant excitation among multiple super Earths. Moreover, the ring-like dust distribution induced by planets with small eccentricities is axisymmetric even in low viscosity environments, consistent with the paucity of vortices in Atacama Large Millimeter/submillimeter Array images.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The gallium isotopic composition of the Moon

Here, in this study, we present new Ga isotope data from a suite of 28 mare basalts and lunar highland rocks. The δ 71 Ga values of these samples range from -0.10 to +0.66‰ (where δ 71 Ga is the relative difference between the 71 Ga/ 69 Ga ratio of a sample and the Ga-IPGP standard), which is an order of magnitude more heterogeneous than δ 71 Ga values in terrestrial magmatic rocks. The cause of this isotopic heterogeneity must be established to estimate the bulk δ 71 Ga value of the Moon. In general, low-Ti basalts and ferroan anorthosite suite (FAS) rocks have δ 71 Ga values that are lower than high-Ti basalts and KREEP-rich rocks. The observation that rocks derived from later forming LMO cumulates have higher δ 71 Ga values suggests that Ga isotopes are fractionated by processes that operate within the chemically evolving LMO, rather than localized degassing or volatile redistribution. Correlations between indices of plagioclase removal from the LMO (e.g. Eu/Eu*) with Ga isotope ratios suggest that a Δ 71 Ga plagioclase-melt of -0.3‰, (where Δ 71 Ga plagioclase-melt is the isotopic fractionation associated with crystallization of plagioclase from a melt), could drive the observed isotopic fractionation in high-Ti mare basalts and KREEP-rich rocks. This would be consistent with the observation that FAS rocks have δ 71 Ga values that are lower than mare basalts. However, the addition of KREEP-like material into the mare basalt source regions would not contribute enough Ga to perturb the isotopic composition outside of analytical uncertainty. Thus, basalts derived from early formed LMO cumulates such as those from Apollo 15, would preserve light Ga isotopic compositions despite containing modest amounts of urKREEP. We estimate that the δ 71 Ga value of the LMO was ~0.14‰ prior to the onset of plagioclase crystallization and extraction. Whether this δ 71 Ga value is representative of the initial BSM cannot be ascertained from the current dataset. It remains plausible that the Moon accreted with a heavier Ga isotopic composition than the Earth. Alternatively, the Moon and Earth could have accreted with similar isotopic compositions (BSE = 0.00 ± 0.06‰, Kato et al., 2017) and volatile loss drove the LMO to higher δ 71 Ga values prior to formation of the lunar crust.

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

Noble gas insights into early impact delivery and volcanic outgassing to Earth's atmosphere: A limited role for the continental crust

Earth's atmosphere, crust and mantle have evolved together through continuous geochemical exchange throughout Earth's history. Constraints on the transport of volatile elements and compounds between these reservoirs are crucial for understanding how Earth could have stayed habitable for extended periods of time. In this work, we present a new forward model of He, Ne and Ar in the mantle, crust and atmosphere. We explore concentrations of noble gases at the end of accretion, bulk silicate Earth K/U ratios, crustal growth scenarios, and upper mantle processing rates throughout Earth's history. We search for parameter combinations that simultaneously satisfy observational constraints on present-day mantle 4 He/ 3 He and 40 Ar/ 36 Ar (sensitive to mantle outgassing and continental crust growth, which depletes the mantle of U, Th and K), atmospheric 20 Ne/ 22 Ne (which tracks the mix of outgassed vs. delivered Ne in the atmosphere), and atmospheric 40 Ar/ 36 Ar in the past and today (sensitive to volatile delivery, mantle outgassing and crustal growth). Leveraging this intertwined set of noble gas abundances and isotopic compositions yields new constraints on initial noble gas abundances in the mantle and on the proportions of atmospheric volatiles originating from delivery by impact degassing, mantle outgassing, and degassing of the continental crust. We find that atmospheric Ar isotopic evolution is primarily sensitive to the mantle processing rate history; the atmospheric Ar isotopic record should therefore not be used to reconstruct continental crust growth, but instead provides valuable insights into mantle processing rates. Our model predicts a measurably low 20 Ne/ 22 Ne ratio of ~9.7 in Archean atmospheric samples. Most of atmospheric primordial 36 Ar was directly delivered by chondritic bodies and not transferred to the atmosphere during an intense early episode of mantle outgassing. Nitrogen delivered by impact degassing could account for the present-day atmospheric nitrogen inventory.

58 GEOSCIENCES↗
DOE OSTI2020

The Preservation of Super-Earths and the Emergence of Gas Giants after Their Progenitor Cores Have Entered the Pebble-isolation Phase

The omnipresence of super-Earths suggests that they are able to be retained in natal disks around low-mass stars, whereas exoplanets' mass distributions indicate that some cores have transformed into gas giants through runaway gas accretion at ≳1 au from solar-type stars. In this paper, we show that transition to runaway gas accretion by cores may be self-impeded by an increase of the grain opacity in their envelope after they have acquired sufficient mass (typically ~ 10M ⊕ ) to enter a pebble-isolation phase. The accumulation of approximately millimeter- to meter-sized pebbles in their migration barriers enhances their local fragmentation rates. The freshly produced submillimeter grains pass through the barrier, elevate the effective dust opacity, and reduce the radiative flux in the core's envelope. These effects alone are adequate to suppress the transition to runaway accretion and preserve super-Earths in the stellar proximity (~0.1 au), albeit entropy advection between the envelope and the disk can further reduce the accretion rate. At an intermediate distance (~1 au) from their host stars, the escalation in the dust opacity dominates over entropy advection in stalling the transition to runaway accretion for marginally pebble-isolated cores. Beyond a few astronomical unit, the transformation of more-massive cores to gas giants is reachable before severe depletion of disk gas. This requirement can be satisfied either in extended disks with large scale height via orderly accretion of migrating pebbles or through the mergers of oligarchic protoplanetary embryos, and can account for the correlated occurrence of long-period gas giants and close-in super-Earths.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Coastal Wetland Carbon Sequestration in a Warmer Climate (Final Report)

Coastal wetlands are global hotspots of carbon storage and locations where carbon and nitrogen cycles have a disproportionately large impact on land, water, and air in comparison to the area they occupy. The extremely high rates of carbon sequestration in these systems are the result of complex feedbacks between vegetation and the physical environment. The plant and microbial interactions that drive these feedbacks are absent in ecosystem- and global-scale Earth System models. The goal of this proposal was to address the uncertain response of the coastal carbon sink to climate change, a critical knowledge gap in coastal carbon research and a significant barrier to incorporating coastal wetlands into Earth System Models. The Salt Marsh Accretion Response to Temperature eXperiment (SMARTX) was designed to understand the ecosystem-scale consequences of warming and elevated CO₂ (eCO₂) in tidal wetlands at the coastal terrestrial-aquatic interface. We successfully designed and built a novel whole-ecosystem warming and eCO₂ experiment in a coastal wetland. The gradient design of the warming treatments (+0, +1.7, +3.4, +5.1 ⁰C) allowed us to discover unexpected non-linear and non-additive responses arising from plant-microbe interactions. We concluded that these non-linear responses to warming were the result of plants and microbes responding to temperature at different thresholds, with important consequences for forecasting terrestrial ecosystem responses to climate change and understanding global trends.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Microphysically modified magnetosonic modes in collisionless, high-β plasmas

With the support of hybrid-kinetic simulations and analytic theory, we describe the nonlinear behaviour of long-wavelength non-propagating (NP) modes and fast magnetosonic waves in high-β collisionless plasmas, with particular attention to their excitation of and reaction to kinetic micro-instabilities. The perpendicularly pressure balanced polarization of NP modes produces an excess of perpendicular pressure over parallel pressure in regions where the plasma β is increased. For mode amplitudes |δB/B 0 |≳0.3, this excess excites the mirror instability. Particle scattering off these micro-scale mirrors frustrates the nonlinear saturation of transit-time damping, ensuring that large-amplitude NP modes continue their decay to small amplitudes. At asymptotically large wavelengths, we predict that the mirror-induced scattering will be large enough to interrupt transit-time damping entirely, isotropizing the pressure perturbations and morphing the collisionless NP mode into the magnetohydrodynamic (MHD) entropy mode. In fast waves, a fluctuating pressure anisotropy drives both mirror and firehose instabilities when the wave amplitude satisfies |δB/B 0 |≳2β -1 . The induced particle scattering leads to delayed shock formation and MHD-like wave dynamics. Taken alongside prior work on self-interrupting Alfvén waves and self-sustaining ion-acoustic waves, our results establish a foundation for new theories of electromagnetic turbulence in low-collisionality, high-β plasmas such as the intracluster medium, radiatively inefficient accretion flows and the near-Earth solar wind.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Understanding the eco-geomorphologic feedback of coastal marsh under sea level rise: vegetation dynamic representations, processes interaction, and parametric sensitivity

A growing number of coastal eco-geomorphologic modeling studies have been conducted to understand coastal marsh evolution under sea-level rise (SLR). Although these models quantify marsh topographic change as a function of sedimentation and erosion, their representations of vegetation dynamics that control organic sedimentation differ. How vegetation dynamic schemes contribute to simulation outcomes is not well quantified. Additionally, the sensitivity of modeling outcomes to parameter selection in the available formulations has not been rigorously tested to date, especially under the influence of an accelerating SLR. In this paper, we used a coastal eco-geomorphologic model with different vegetation dynamic schemes to investigate the eco-geomorphologic feedbacks of coastal marshes and parametric sensitivity under SLR scenarios. We found that marsh platform relief increased with sea level rise rate. The simulations with different vegetation schemes exhibited different spatial-temporal variations in elevation and biomass. The nonlinear Spartina scheme presented the most resilient prediction with generally the highest marsh accretion and vegetation biomass, and the least elevation relief under SLR. But the linear Spartina scheme predicts the lowest unvegetated-vegetated ratio. We also found that vegetation-related parameters and sediment diffusivity, which were not well measured or discussed in previous studies, were identified as some of the most critical parameters. Additionally, the model sensitivity to vegetation-related parameters increased with SLR rates. The identified most sensitive parameters may inform how to appropriately choose modeling representations of key processes and parameters for different coastal marsh landscapes under SLR, and demonstrate the importance of future field measurements of these key parameters.

58 GEOSCIENCES↗
DOE OSTI2020

River bank erosion and accretion rates, and planform metrics of the East River, downstream of Gothic Colorado over intervals between 1955 and 2015

This dataset provides the tabular summary of analysis of an alluvial floodplain reach of the East River, downstream of Gothic, CO near Crested Butte. The measurements include bank erosion and accretion rates, channel width, bank curvature, and the aspect/orientation of the river banks. The measurements were derived from binary masks of the location of the river channel from aerial and satellite imagery collected on seven dates over a sixty-year timespan: 1955, 1973, 1983, 1990, 2001, 2011, and 2015. The masks were analyzed using the Spatially Continuous Riverbank Erosion and Accretion Measurements (SCREAM) software detailed in Rowland et al. 2016. The masks used in this analysis can be found an accompanying dataset (DOI: ). Rates of change along the East River were measured over a total of 8-time intervals: 1955-1973; 1955-2015; 1973-1983; 1983-1990; 1990-2001; 2001-2011; 2001-2015; and 2011-2015.In files with “summary” in the name, the data is provided at a pixel level, where each mapped bank pixel has an associated erosion or accretion value, a channel width, a curvature value, and an aspect each river and time period will have an individual file. Files with “Segments” in name provide data that is averaged along segments of the rivers. These segments are approximately 10 channel widths in length. In addition to erosion and accretion rates, the segment-based results include area measurements of erosion and accretion, islands, and channels. The number of islands is also included.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Ion and Electron Acoustic Bursts during Anti-Parallel Reconnection Driven by Lasers

Magnetic reconnection converts magnetic energy into thermal and kinetic energy in plasma. Among the numerous candidate mechanisms, ion acoustic instabilities driven by the relative drift between ions and electrons (or equivalently, electric current) have been suggested to play a critical role in dissipating magnetic energy in collisionless plasmas. However, their existence and effectiveness during reconnection have not been well understood due to ion Landau damping and difficulties in resolving the Debye length scale in the laboratory. We report a sudden onset of ion acoustic bursts measured by collective Thomson scattering in the exhaust of anti-parallel magnetically driven reconnection using high-power lasers. The ion acoustic bursts are followed by electron acoustic bursts with electron heating and bulk acceleration. We reproduce these observations with one- and two-dimensional particle-in-cell simulations in which an electron outflow jet drives ion acoustic instabilities, forming double layers. These layers induce electron two-stream instabilities that generate electron acoustic bursts and energize electrons. Our results demonstrate the importance of ion and electron acoustic dynamics during reconnection when ion Landau damping is ineffective, a condition applicable to a range of astrophysical plasmas including near-Earth space, stellar flares and black hole accretion engines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Origin of 182 W Anomalies in Ocean Island Basalts

Ocean island basalts (OIB) show variable 182 W deficits that have been attributed to either early differentiation of the mantle or core-mantle interaction. However, 182 W variations may also reflect nucleosynthetic isotope heterogeneity inherited from Earth's building material, which would be evident from correlated 182 W and 183 W anomalies. Some datasets for OIB indeed show hints for such correlated variations, meaning that a nucleosynthetic origin of W isotope anomalies in OIB cannot be excluded. We report high-precision W isotope data for OIB from Samoa and Hawaii, which confirm previously reported 182 W deficits for these samples, but also demonstrate that none of these samples have resolvable 183 W anomalies. These data therefore rule out a nucleosynthetic origin of the 182 W deficits in OIB, which most likely reflect the entrainment of either core material or an overabundance of late-accreted materials within OIB mantle sources. If these processes occurred over Earth's history, they may have also been responsible for shifting the 182 W composition of the bulk mantle to its modern-day value. We also report Mo isotope data for some Hawaiian OIB, which reveal no resolved nucleosynthetic Mo isotopic anomalies. This is consistent with inheritance of 182 W deficits in OIB from the addition of either core or late-accreted material, but only if these materials have a non-carbonaceous (NC) meteorite-like heritage. As such, these data rule out significant contributions of carbonaceous chondrite (CC)-like materials to either Earth's core or late accretion.

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

Biogeochemistry simulations for the Salt Marsh Accretion Response to Temperature eXperiment (SMARTX)

Coastal ecosystems have been largely ignored in Earth system models but are important zones for carbon and nutrient processing. Interactions between water, microbes, soil, sediments, and vegetation are important for mechanistic representation of coastal processes and ecosystem function. To investigate the role of these feedbacks, we used a reactive transport model (PFLOTRAN) that has the capability to be connected to the Energy Exascale Earth System Model (E3SM). PFLOTRAN was used to incorporate redox reactions and track chemical species important for coastal ecosystems as well as define simple representations of vegetation dynamics. Our goal was to incorporate oxygen flux, salinity, pH, sulfur cycling, and methane production along with plant-mediated transport of gases and tidal flux. Using porewater profile and incubation data for model calibration and evaluation, we were able to create depth-resolved biogeochemical soil profiles for saltmarsh habitat and use this updated representation to simulate direct and indirect effects of elevated CO2 and temperature on subsurface biogeochemical cycling. We found that simply changing the partial pressure of CO2 or increasing temperature in the model did not fully reproduce observed changes in the porewater profile, but the inclusion of plant or microbial responses to CO2 and temperature manipulations was more accurate in representing porewater concentrations. This indicates the importance of characterizing tightly coupled vegetation-subsurface processes for developing predictive understanding and the need for measurement of plant-soil interactions on the same time scale to understand how hotspots or moments are generated.Included in this data package are PFLOTRAN input (PFLOTRAN input files and chemical database) files for simulating single column biogeochemistry, root, and tide interactions at the Global Change Research Wetland (Kirkpatrick Marsh; Edgewater, MD). The biogeochemical network includes soil organic matter decomposition, nitrogen, iron, and sulfur cycling, and methanogenesis. Reduced species can be oxidized and plant processes include oxygen and nutrient priming, methane release, and nutrient uptake.Inputs:TAI_database.dat - geochemical database for reactions, more information on database structure and variables can be found here https://www.pflotran.org/documentation/user_guide/cards/pages/geochemical_database.htmlswamp.in - input file for biogeochemical network in PFLOTRANswamp_eCO2.in - input file for biogeochemical network in PFLOTRAN with input gas partial pressures/concentrations adjusted for elevated CO2 treatmentsOutputs:swamp_obs_0.tec - hourly porewater concentrations from from multiple depths in the soil columnswamp_eCO2_obs_0.tec - hourly porewater concentrations from multiple depths in the soil column for elevated CO2 treatmentsPFLOTRAN code access: https://github.com/fmyuan/pflotran-elm-interface.git

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Activity model for 36 elements in Fe-Ni-Si-S-C liquids with application to terrestrial planet accretion and mantle geochemistry: New data for Ru, Re, Pt, Os, Ti, Nb, and Ta

Understanding siderophile element partitioning between metal and silicate melts under diverse conditions can be used to place important constraints on the materials and conditions of planetary accretion and core formation, as well as post core formation processes. However, the effects of Si on the partitioning and activity coefficients for these elements are not well known, despite Si likely being one of the dominant light elements in Earth’s core. To address this gap in understanding, we have undertaken a systematic study of the highly siderophile elements Re, Pt, Os, and Ru, and the refractory lithophile elements Nb, Ta and Ti at 1600 °C and 1 GPa, to derive epsilon interaction parameters for these elements in FeSi metallic liquids. Positive epsilon interaction parameters were measured for Nb, Ta, Ti, Ru, Re, Pt, and Os, indicating that dissolved Si in Fe liquids causes a decrease in their metal/silicate partition coefficients (or ‘silicophobic’ behavior). Furthermore, ε$_{Re, Os, or Ru}^{Si}$ > ε$_{Re, Os, or Ru}^{S}$ which means Si causes a larger decrease in D(metal/silicate) than S, and the chalcophile behavior expected from some elements will be completely masked by the presence of Si in a metallic liquid. The new parameters are used to update an activity model that now includes 36 siderophile elements in Fe-Ni-Si-S-C liquids (27 trace elements considered here). Systematic assessment of these 27 elements shows which have the strongest affinity for Si, C, and S, and also how activity coefficients for these elements would vary during accretion and core formation in Earth, Mars, and Mercury of widely differing fO 2 and core compositional conditions. The activity model is combined with new partitioning expressions for Mo, W, Cr, Re Ru, Pt, and Os and applied to aspects of post core formation mantle geochemistry of Earth, Mars, and Mercury. Our updated expressions show that the BSE Mo/W ratio can easily be achieved with metal/silicate partitioning during growth of the Earth, whereas Re, Os and Ru become lower than and highly fractionated compared with BSE values during core formation and accretion, and thus nearly 99% of their BSE abundances are likely contributed by late accretion. Ru isotopes should be a very good indicator of the source material for the late accretion. The high Pt/Os and Re/Os developed in a deepening magma ocean during the growth of the Earth, indicates 186 Os and 187 Os isotopes could be coupled if this ancient material remained isolated and subsequently became entrained in mantle plumes and measured in surficial lavas. The extent to which this occurred will be limited by the low Os content of this ancient material, thus requiring mixing as a major component in plume sources. Martian mantle Hf/W ratio stays low during accretion and core formation modelling, suggesting that W isotope anomalies are more likely due to solid/liquid silicate fractionation than to core formation. Finally, Ti contents measured by MESSENGER at Mercury’s surface can be explained by segregation of either a metallic core (IW-6 to -8) or metallic core + sulfide (IW-4 to -7.5) followed by mantle melting.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Comparison of Polarized Radiative Transfer Codes Used by the EHT Collaboration

Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT Collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I , Q , U , and V , respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗