Search NASA⌕ Search

SEARCH · Search NASA

Results for “Planetary differentiation”

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 OSTI2025

The variety and origin of materials accreted by Bennu’s parent asteroid

The first bodies to form in the Solar System acquired their materials from stars, the presolar molecular cloud and the protoplanetary disk. Asteroids that have not undergone planetary differentiation retain evidence of these primary accreted materials. However, geologic processes such as hydrothermal alteration can dramatically change their bulk mineralogy, isotopic compositions and chemistry. Here we analyse the elemental and isotopic compositions of samples from asteroid Bennu to uncover the sources and types of material accreted by its parent body. We show that some primary accreted materials escaped the extensive aqueous alteration that occurred on the parent asteroid, including presolar grains from ancient stars, organic matter from the outer Solar System or molecular cloud, refractory solids that formed close to the Sun, and dust enriched in neutron-rich Ti isotopes. We find Bennu to be richer in isotopically anomalous organic matter, anhydrous silicates, and light isotopes of K and Zn than its closest compositional counterparts, asteroid Ryugu and Ivuna-type (CI) carbonaceous chondrite meteorites. We propose that the parent bodies of Bennu, Ryugu and CI chondrites formed from a common but spatially and/or temporally heterogeneous reservoir of materials in the outer protoplanetary disk.

Kuiper belt↗
DOE OSTI2024

Regularized Differentiation for Bioburden Density Estimation in Planetary Protection

In this paper, we propose and investigate the performance of two novel shrinkage estimators for bioburden density estimation in planetary protection. The estimators are based on the regularized differentiation of a cumulative count of colony forming units collected throughout the data collecting session or the life cycle of the entire mission. The regularized differentiation recasts the problem of bioburden density estimation as a linear least squares problem. The least squares problem is then solved through regularization techniques, such as truncated singular value decomposition and penalized least squares. The regularization is necessary to avoid noise amplification during the differentiation of noisy data. The two regularization estimators are compared with four other commonly used estimators to simultaneously evaluate the means of multivariable independent Poisson distributions: the maximum likelihood, noninformative Bayes estimator with Jeffreys prior, Empirical Bayes using conjugate gamma-Poisson model with gamma parameters selected by method of moments, and the Clevenson-Zidek estimator. It is shown through computer-simulated data that the regularized differentiation based on ridge regression has the smallest mean-squared error among all estimators. The analysis of shrinkage mechanism implemented by regularized differentiation is performed, and it is shown that the regularized differentiation amounts to performing a weighted averaging of all the samples. The weights are determined by the regularization parameter automatically selected by the L-curve technique. Since the method of least squares makes no distributional assumptions about the data, it presents an attractive technique for bioburden density estimation when there are concerns about the misspecification of the distributional model. The paper concludes with the analysis of the bioburden data collected during InSight mission and directions for future work.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Tungsten and molybdenum isotopic evidence for an impact origin of pallasites

The origin of pallasites—stony-iron meteorites mainly composed of olivine and Fe-Ni metal—is debated and proposed formation scenarios broadly range from models that explain pallasite formation by internal processes in the mantle of a differentiated planetesimal to those that involve impact–induced mixing of core and mantle materials. Here, the origin of pallasites is examined by studying the nebular source regions of their precursor material using Mo isotopes and their history of metal-silicate segregation using Hf-W chronometry. We report new Mo and W isotopic data for a large suite of pallasite metal samples, alongside Pt isotope data to quantify superimposed cosmic ray exposure effects. Most main-group pallasites exhibit uniform pre-exposure 182 W and Mo isotopic compositions that bear an excellent similarity to those of IIIAB iron meteorites. Four main-group pallasites and the IIIAB iron Thunda have more radiogenic pre-exposure 182 W compositions, but display the same Mo isotopic composition as other main-group pallasites and IIIAB irons. This strong chronological and genetic link strongly suggests that main-group pallasite metal originated in the IIIAB parent body core. This, combined with prior Pd-Ag chronometric evidence for an early collisional disruption of the IIIAB parent body, implies that main-group pallasites formed by impact–induced mixing of metal and silicates rather than by an internal process on the IIIAB parent body. This mixing led to elevated 182 W compositions in some pallasites, which are best accounted for by partial re-equilibration of IIIAB metal with radiogenic 182 W from the colliding body. Altogether, our results support models that explain main-group pallasite formation by injection of pallasite metal into the mantle of another differentiated body, implying that pallasite silicates did not primarily derive from the IIIAB mantle, but instead from that of the colliding body.

58 GEOSCIENCES↗
DOE OSTI2022

Ca-isotopes as a robust tracer of magmatic differentiation

The large mass difference (~10%) between the two most abundant isotopes of calcium, 40 Ca and 44 Ca, gives Ca great potential in tracking mass-dependent fractionation during magmatic processes. Resolvable Ca-isotope fractionation during fractional crystallization of magma, particularly by feldspar in evolved melts, has been theoretically inferred but not robustly tested in nature. To further explore the effects of magmatic differentiation on Ca-isotope systematics, we studied the late-Permian alkaline igneous suite of the Øyangen Caldera, Oslo Rift, Norway, consisting of volcanic and intrusive units ranging from basanitic to rhyolitic compositions. Major and trace element variations and modeling demonstrate that the main series of samples (N = 21), including basanites, ring-dyke syenites, and central-dome syenites, likely documents a co-genetic and closed-system fractional crystallization sequence. Our data show minimal δ 44/40 Ca variation (< 0.05 ‰) in the intermediate magma and a marked increase in δ 44/40 Ca in the felsic magma of the Øyangen Caldera (from 0.62 ± 0.02 ‰ to 1.15 ± 0.03 ‰ relative to Ca standard, SRM915a). The systematic increase is best explained by equilibrium isotopic fractionation dominated by alkali feldspar in the fractionating mineral assemblage. This is further supported by strong correlations between δ 44/40 Ca, CaO, and Eu/Eu* in the main-series samples. Implementing a Monte Carlo approach, isotopic modeling of the liquid line of descent using Rayleigh fractionation is highly consistent with the observed Ca-isotope evolution. For the first time, we confirm prominent Ca stable isotope fractionation in felsic-stage differentiation of alkaline magma and constrain the isotope fractionation factors of plagioclase and K-feldspar. Integrated with extant estimations on mineral fractionation factors from the literature, our results suggest increasing fractionation effects of rock-forming minerals with decreasing Ca content. Finally, the affirmation of significant Ca-isotope fractionation in alkaline magma by feldspar empowers the application of Ca as a versatile tracer of crustal evolution, allowing further tests in other magmatic conditions across various planetary objects.

58 GEOSCIENCES↗
DOE OSTI2020

Onset of magma ocean solidification on Mars inferred from Mn-Cr chronometry

The mantle of Mars probably differentiated through the crystallization of a magma ocean during the first tens of million years (Ma) of Solar System evolution. However, the exact timescale of large-scale silicate differentiation of the martian mantle is debated, and in particular, it remains unclear when differentiation commenced. In this paper, we applied the short-lived 53 Mn- 53 Cr system to martian meteorites in order to date the onset of large-scale mantle differentiation on Mars. The new Cr isotope data demonstrate that martian meteorites exhibit no resolvable radiogenic 53 Cr variations, and instead have a uniform +20.2±1.2 (95% conf.) parts-per-million excess in 53 Cr/ 52 Cr relative to the terrestrial mantle. The investigated groups of martian meteorites are lithologically varied and derive from diverse mantle sources that probably had variable Mn/Cr. Hence, the lack of 53 Cr variability among martian meteorites demonstrates that silicate differentiation on Mars occurred after the extinction of 53 Mn. Provided that the sources of the martian meteorites have Mn/Cr variations that are typical of the terrestrial planets, this result implies that the onset of large-scale silicate differentiation must have occurred later than 20±5 Ma after Solar System formation. The onset of silicate differentiation on Mars inferred here is significantly later than time estimates for segregation of the martian core which conservatively occurred within <10 Ma after Solar System formation. Thus, the new Mn-Cr data imply that there was a small, but resolvable, time gap of at least 5 Ma between core formation and magma ocean solidification on Mars. If the age of core segregation is taken at face value, our results imply that the martian magma ocean remained mostly molten over several Ma. This inferred longevity of the magma ocean is inconsistent with thermal models predicting rapid (<1 Ma) solidification of the martian magma ocean. Although there is currently no unique solution to this conundrum, our results can potentially be explained by a protracted history of impact bombardment that delayed differentiation in a shallow magma ocean on Mars, or perhaps more readily, by the presence of an early and dense atmosphere that acted as an insulator and prevented the magma ocean from cooling quickly.

58 GEOSCIENCES↗
DOE OSTI2025

A tale of two planets: Disparate evolutionary models for Mars inferred from radiogenic isotope compositions of Martian meteorites

The radiogenic isotopic compositions of basaltic Martian meteorites (shergottites) and clinopyroxene/olivine cumulate meteorites (nakhlite/chassignites) are used to define the global evolution of Mars. However, the two main groups of meteorites demonstrate that their sources underwent divergent styles of magmatic evolution. The shergottites portray a planet that differentiated ~4.52 billion years ago via solidification of a magma ocean, producing incompatible element-depleted and -enriched reservoirs that remained isolated until melt production. In contrast, the reservoir from which the nakhlite/chassignites derive may have formed earlier, produced melts that fractionated Sm/Nd and Hf/W differently, was compositionally less variable, and experienced a significantly more complex history following primordial differentiation than the shergottite sources. The disparate histories recorded by these two groups of meteorites elucidate important questions that could be addressed by acquiring additional samples. Obtaining samples that shared the isotopic systematics of the shergottites would provide confidence that extrapolating the primordial differentiation history of Mars from shergottite radiogenic isotope systematics is reasonable. Returned samples from Mars will also constrain the physical locations of the meteorite source regions, providing insights into the general structure of the Martian mantle. In addition, they will help constrain the phases present in the martian mantle during melting and the conditions under which they are stable. Finally, identifying an evolved lithology that satisfies the geochemical and isotopic constraints placed on the incompatible element-enriched endmember observed in the shergottites would define the nature of magmatic evolution on Mars and whether it is more akin to processes on the Earth or the Moon.

Mars↗
DOE OSTI2022

Formation, cooling history and age of impact events on the IIE iron parent body: Evidence from the Miles meteorite

Most iron meteorites formed in planetary cores during differentiation, but the IIE iron meteorites have chemical and physical features that are inconsistent with this origin. By combining mineral chemistry, mineral modes and three-dimensional petrography, we reconstruct the bulk chemistry of the felsic silicate-bearing Miles IIE iron meteorite and demonstrate that the silicate inclusion compositions are similar to partial melts produced experimentally from an H chondrite composition. We use the reconstructed bulk composition, mineralogy and thermodynamic modelling to show that melting above ~ 1200°C under reducing conditions formed metal (Fe-Ni alloy) and felsic silicate partial melts. Upon cooling, the melts crystallized Mg-rich pyroxenes, Na- and K-rich feldspars, and tridymite. Importantly, this mechanism enriches cosmochemically volatile elements (i.e., those with a 50% condensation temperature of ~430–830°C, like Na and K) to the level found in the felsic silicate inclusions. The presence of crystallographically disordered srilankite (only stable above 1160°C) and an absence of Widmanstätten texture require both high peak temperatures and rapid cooling, which cannot be explained by core formation. Instead, they point to small melt volumes, a transient heat pulse, and small thermal mass, and imply efficient physical segregation of silicate and metallic melts through buoyancy separation followed by rapid cooling that arrested the separation of metal and silicate liquid phases. In situ 207 Pb/ 206 Pb age of 4542.3 ± 4.0 Ma in Zr-oxide and phosphate minerals dates the melting event that formed the silicate inclusions. This age aligns with the earliest ages found in other IIE iron meteorite silicates and requires a heating event ~25 million years after the solar system formed. We found 39 Ar/ 40 Ar ages of 3495 ± 52 Ma (low-T) and 4303 ± 7 Ma (high-T) in a K-feldspar grain, with the 3495 Ma age aligning with later thermal events recorded in other IIE iron meteorites. Dating reveals the complex petrogenetic and thermal history of Miles and the IIE iron meteorites. This is the first IIE iron meteorite found to record evidence of heating at 4.5 and 3.5 Ga likely from impact events. We propose that high-velocity impact(s) into an iron-rich, porous chondritic parent body at ~4.54 Ga produced immiscible metal and silicate melts that cooled rapidly and trapped low density silicate inclusions within high density metal. Other IIE irons that formed at lower peak temperatures (900–1000°C) contain chondritic silicate inclusions and relict chondrules, supporting this conceptual model. In conclusion, this petrogenesis is consistent with thermodynamic modelling, experimental data and the wide range of peak temperatures and cooling rates observed in the IIE iron meteorites.

58 GEOSCIENCES↗
DOE OSTI2025

Sound Velocity and Grüneisen Parameter in Shock‐Melted Silica at Deep Earth Pressures

Silica is a primary component of rocky planet interiors and its melt properties are important for understanding planetary formation and differentiation, magma oceans, and the deep mantle. Although well understood in the solid state, the high-pressure behavior of liquid silica is poorly constrained at lower mantle pressures. Using laser interferometry to measure shock wave profiles, we report measured stress-density states and longitudinal sound speeds in shock-synthesized stishovite, from fused silica staring material, and across the solid-liquid phase boundary up to 154 GPa. Our results constrain completion of melt at 80 GPa and show that at pressures relevant to the deep mantles of Earth-sized, rocky planets, the Grüneisen parameter for liquid silica increases with compression. This finding is consistent with a continuous increase in Si-O coordination above six for liquid silica at core-mantle boundary relevant pressures and temperatures.

Ocampo, I. K. [Princeton University, NJ (United St↗
DOE OSTI2021

Detecting Atmospheric Molecules of Nontransiting Temperate Terrestrial Exoplanets Using High-resolution Spectroscopy in the Mid-infrared Domain

Motivated by the development of high-dispersion spectrographs in the mid-infrared (MIR) range, we study their application to the atmospheric characterization of nearby nontransiting temperate terrestrial planets around M-type stars. We examine the detectability of CO{sub 2}, H{sub 2}O, N{sub 2}O, and O{sub 3 }features in high-resolution planetary thermal emission spectra at 12–18 μm assuming an Earth-like profile and a simplified thermal structure. The molecular line width of such planets can be comparable to or broader than the Doppler shift due to the planetary orbital motion. Given the likely difficulty in knowing the high-resolution MIR spectrum of the host star with sufficient accuracy, we propose observing the target system at two quadrature phases and extracting the differential spectra as the planetary signal. In this case, the signals can be substantially suppressed compared with the case where the host star spectrum is perfectly known, as some parts of the spectral features do not remain in the differential spectra. Despite this self-subtraction, the CO{sub 2} and H{sub 2}O features of nearby (≲5 pc) systems with mid-/late-M host stars would be feasible with a 6.5 m class cryogenic space telescope, and orbital inclination could also be constrained for some of them. For CO{sub 2} and N{sub 2}O in a 1 bar Earth-like atmosphere, this method would be sensitive when the mixing ratio is 1–10{sup 3} ppm. The detectability of molecules except O{sub 3} is not significantly improved when the spectral resolution is higher than R≳10,000, although the constraint on the orbital inclination is improved. This study provides some benchmark cases useful for assessing the value of MIR high-resolution spectroscopy in terms of characterization of potentially habitable planets.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Temperature and Density on the Forsterite Liquid‐Vapor Phase Boundary

Abstract The physical processes during planet formation span a large range of pressures and temperatures. Giant impacts, such as the one that formed the Moon, achieve peak pressures of 100s of GPa. The peak shock states generate sufficient entropy such that subsequent decompression to low pressures intersects the liquid‐vapor phase boundary. The entire shock‐and‐release thermodynamic path must be calculated accurately in order to predict the post‐impact structures of planetary bodies. Forsterite (Mg 2 SiO 4 ) is a commonly used mineral to represent the mantles of differentiated bodies in hydrocode models of planetary collisions. Here, we performed shock experiments on the Sandia Z Machine to obtain the density and temperature of the liquid branch of the liquid‐vapor phase boundary of forsterite. This work is combined with previous work constraining pressure, density, temperature, and entropy of the forsterite principal Hugoniot. We find that the vapor curves in previous forsterite equation of state models used in giant impacts vary substantially from our experimental results, and we compare our results to a recently updated equation of state. We have also found that due to under‐predicted entropy production on the principal Hugoniot and elevated temperatures of the liquid vapor phase boundary of these past models, past impact studies may have underestimated vapor production. Furthermore, our results provide experimental support to the idea that giant impacts can transform much of the mantles of rocky planets into supercritical fluids.

58 GEOSCIENCES↗
DOE OSTI2021

Bayesian constraints on the origin and geology of exoplanetary material using a population of externally polluted white dwarfs

ABSTRACT White dwarfs that have accreted planetary bodies are a powerful probe of the bulk composition of exoplanetary material. In this paper, we present a Bayesian model to explain the abundances observed in the atmospheres of 202 DZ white dwarfs by considering the heating, geochemical differentiation, and collisional processes experienced by the planetary bodies accreted, as well as gravitational sinking. The majority (>60 per cent) of systems are consistent with the accretion of primitive material. We attribute the small spread in refractory abundances observed to a similar spread in the initial planet-forming material, as seen in the compositions of nearby stars. A range in Na abundances in the pollutant material is attributed to a range in formation temperatures from below 1000 K to higher than 1400 K, suggesting that pollutant material arrives in white dwarf atmospheres from a variety of radial locations. We also find that Solar System-like differentiation is common place in exoplanetary systems. Extreme siderophile (Fe, Ni, or Cr) abundances in eight systems require the accretion of a core-rich fragment of a larger differentiated body to at least a 3σ significance, whilst one system shows evidence that it accreted a crust-rich fragment. In systems where the abundances suggest that accretion has finished (13/202), the total mass accreted can be calculated. The 13 systems are estimated to have accreted masses ranging from the mass of the Moon to half that of Vesta. Our analysis suggests that accretion continues for 11 Myrs on average.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Contrasting biological production trends over land and ocean

Terrestrial and marine ecosystems constitute the primary components of the Earth’s biosphere, yet their photosynthetic productions are typically studied separately, which limits understanding of planetary carbon uptake and biosphere health. Here, in this study, using multiple satellite-derived products, we identify contrasting net primary production (NPP) trends between land and ocean, probably reflecting their differential sensitivity to climate warming, especially in tropical regions. Planetary NPP shows an overall increase of 0.11 ± 0.13 PgC yr −1 (P = 0.05) from 2003 to 2021, driven by a significant terrestrial enhancement of 0.20 ± 0.07 PgC yr −1 (P < 0.001) and partially offset by an oceanic decline of −0.12 ± 0.12 PgC yr −1 (P = 0.07). While land contributes to the strong upwards NPP trend, the interannual variability in global NPP is predominantly driven by the ocean, especially during strong El Niño–Southern Oscillation events. Our findings highlight the resilience and potential vulnerability of biosphere primary productivity in a warming climate, calling for integrated land–ocean monitoring and assessment to support climate mitigation initiatives.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Structural evolution of liquid silicates under conditions in Super-Earth interiors

Molten silicates at depth are crucial for planetary evolution, yet their local structure and physical properties under extreme conditions remain elusive due to experimental challenges. In this study, we utilize in situ X-ray diffraction (XRD) at the Matter in Extreme Conditions (MEC) end-station of the Linear Coherent Linac Source (LCLS) at SLAC National Accelerator Laboratory to investigate liquid silicates. Using an ultrabright X-ray source and a high-power optical laser, we probed the local atomic arrangement of shock-compressed liquid (Mg,Fe)SiO 3 with varying Fe content, at pressures from 81(9) to 385(40) GPa. We compared these findings to ab initio molecular dynamics simulations under similar conditions. Results indicate continuous densification of theO-O and Mg-Si networks beyond Earth’s interior pressure range, potentially altering melt properties at extreme conditions. This could have significant implications for early planetary evolution, leading to notable differences in differentiation processes between smaller rocky planets, such as Earth and Venus, and super-Earths, which are exoplanets withmasses nearly three times that of Earth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Water‐Induced Diamond Formation at Earth's Core‐Mantle Boundary

Abstract The carbon and water cycles in the Earth's interior are linked to key planetary processes, such as mantle melting, degassing, chemical differentiation, and advection. However, the role of water in the carbon exchange between the mantle and core is not well known. Here, we show experimental results of a reaction between Fe 3 C and H 2 O at pressures and temperatures of the deep mantle and core‐mantle boundary (CMB). The reaction produces diamond, FeO, and FeH x , suggesting that water can liberate carbon from the core in the form of diamond (“core carbon extraction”) while the core gains hydrogen, if subducted water reaches to the CMB. Therefore, Earth's deep water and carbon cycles can be linked. The extracted core carbon can explain a significant amount of the present‐day mantle carbon. Also, if diamond can be collected by mantle flow in the region, it can result in unusually high seismic‐velocity structures.

Ko, Byeongkwan↗
DOE OSTI2025

A small core in Vesta inferred from Dawn’s observations

Vesta’s large-scale interior structure had previously been constrained primarily using the gravity and shape data from the Dawn mission. However, these data alone still allow a wide range of possibilities for the differentiation state of the body. The moment of inertia is arguably the most diagnostic parameter related to the radial density distribution of a planetary body, making it crucial for assessing the body’s state of internal differentiation. Determining the moment of inertia requires additional measurements of the amplitudes of small rotational motions, such as precession and nutation. Here we report an updated estimate of the moment of inertia of Vesta inferred from Dawn’s Doppler tracking via the Deep Space Network and onboard imaging data. The recovered value for Vesta’s normalized polar moment of inertia is $\overline{C}$/MR 2 = 0.4208 ± 0.0047 (where M is the mass of Vesta and R is the reference radius), which is only 6.6% lower than the homogeneous value of 0.4505. This value, combined with the gravity field and global shape, suggests that Vesta’s interior has limited density stratification beneath its howardite–eucrite–diogenite-dominated crust. We propose two possible origin scenarios that are consistent with the observed constraints. In the first scenario, Vesta’s interior did not undergo full differentiation due to late accretion. In the second scenario, Vesta originated as an impact remnant of a larger differentiated body re-accreted with non-chondritic bulk composition produced from a catastrophic impact. Vesta did not experience complete differentiation in either scenario, suggesting that its current state reflects a complex interplay between its accretion timing, thermal evolution, redistribution of 26 Al bearing melt and/or impact processes.

CNEOS 2014-01-08 bolide↗
DOE OSTI2024

MITgcm-AD v2: Open source tangent linear and adjoint modeling framework for the oceans and atmosphere enabled by the Automatic Differentiation tool Tapenade

The Massachusetts Institute of Technology General Circulation Model (MITgcm) is widely used by the climate science community to simulate planetary atmosphere and ocean circulations. A defining feature of the MITgcm is that it has been developed to be compatible with an algorithmic differentiation (AD) tool, TAF, enabling the generation of tangent-linear and adjoint models. These provide gradient information which enables dynamics-based sensitivity and attribution studies, state and parameter estimation, and rigorous uncertainty quantification. Importantly, gradient information is essential for computing comprehensive sensitivities and performing efficient large-scale data assimilation, ensuring that observations collected from satellites and in-situ measuring instruments can be effectively used to optimize a large uncertain control space. As a result, the MITgcm forms the dynamical core of a key data assimilation product employed by the physical oceanography research community: Estimating the Circulation and Climate of the Ocean (ECCO) state estimate. Although MITgcm and ECCO are used extensively within the research community, the AD tool TAF is proprietary and hence inaccessible to a large proportion of these users. The new version 2 (MITgcm-AD v2) framework introduced here is based on the source-to-source AD tool Tapenade, which has recently been open-sourced. Another feature of Tapenade is that it stores required variables by default (instead of recomputing them) which simplifies the implementation of efficient, AD-compatible code. The framework has been integrated with the MITgcm model’s main branch and is now freely available.

Adjoints↗
DOE OSTI2024

Viscosity Measurements at High Pressures: A Critical Appraisal of Corrections to Stokes' Law

Abstract Fluids and melts in planetary interiors significantly influence geodynamic processes from volcanism to global‐scale differentiation. The roles of these geofluids depend on their viscosities ( η ). Constraining geofluid η at relevant pressures and temperatures relies on laboratory‐based measurements and is most widely done using Stokes' Law viscometry with falling spheres. Yet small sample chambers required by high‐pressure experiments introduce significant drag on the spheres. Several correction schemes are available for Stokes' Law but there is no consensus on the best scheme(s) for high‐pressure experiments. We completed high‐pressure experiments to test the effects of (a) the relative size of the sphere diameter to the chamber diameter and (b) the top and bottom of the chamber, that is, the ends, on the sphere velocities. We examined the influence of current correction schemes on the estimated viscosity using Monte Carlo simulations. We also compared previous viscometry work on various geofluids in different experimental setups/geometries. We find the common schemes for Stokes' Law produce statistically distinct values of η . When inertia of the sphere is negligible, the most appropriate scheme may be the Faxén correction for the chamber walls. Correction for drag due to the chamber ends depends on the precision in the sinking distance and may be ineffective with decreasing sphere size. Combining the wall and end corrections may overcorrect η . We also suggest the uncertainty in η is best captured by the correction rather than propagated errors from experimental parameters. We develop an overlying view of Stokes' Law viscometry at high pressures.

Geochemistry & Geophysics↗
DOE OSTI2026

As above, not so below: Ion fractionation in planetary analog ices

The geophysical evolution and astrobiological potential of ocean worlds are indelibly linked to the chemical compositions of their oceans and ice shells. In the absence of direct measurements, empirical estimates of subsurface ocean compositions have relied on the assumption that the ionic compositions of ice shell surfaces are representative of their underlying ocean compositions. Here, we present experimental results demonstrating that ion fractionation—the differential entrainment of ion species in forming ices—is likely a prevalent process on ocean worlds, suggesting that planetary ice shell compositions do not directly reflect their underlying ocean compositions. We measure in-ice depletions and amplifications of relative ion concentrations ranging between −40 and +77%, compared to the parent fluid composition. Although this may complicate the interpretation of spacecraft data, ion fractionation provides a mechanism for generating compositionally diverse ices that could help explain the geological complexity of planetary ice shells.

58 GEOSCIENCES↗