Investigating Magma Ocean Solidification on Earth Through Laser-Heated Diamond Anvil Cell Experiments

We carried out a series of silicate fractional crystallization experiments at lower mantle pressures using the laser-heated diamond anvil cell. Phase relations and the compositional evolution of the cotectic melt and equilibrium solids along the liquid line of descent were determined and used to assemble the melting phase diagram. In a pyrolitic magma ocean, the first mineral to crystallize in the deep mantle is iron-depleted calcium-bearing bridgmanite. From the phase diagram, we estimate that the initial 33%-36% of the magma ocean will crystallize to form such a buoyant bridgmanite. Substantial calcium solubility in bridgmanite is observed up to 129 GPa, and significantly delays the crystallization of the calcium silicate perovskite phase during magma ocean solidification. Residual melts are strongly iron-enriched as crystallization proceeds, making them denser than any of the coexisting solids at deep mantle conditions, thus supporting the terrestrial basal magma ocean hypothesis (Labrosse et al., 2007).

Counter-Gradient Momentum Transport Through Subtropical Shallow Convection in ICON-LEM Simulations

It is well known that subtropical shallow convection transports heat and water vapor upwards from the surface. It is less clear if it also transports horizontal momentum upwards to significantly affect the trade winds in which it is embedded. We utilize unique multiday large-eddy simulations run over the tropical Atlantic with ICON-LEM to investigate the character of shallow convective momentum transport (CMT). For a typical trade-wind profile during boreal winter, CMT acts as an apparent friction to decelerate the north-easterly flow. This effect maximizes below the cloud base while in the cloud layer, friction is very small, although present over a relatively deep layer. In the cloud layer, the zonal component of the momentum flux is counter-gradient and penetrates deeper than reported in traditional shallow cumulus LES cases. The transport through conditionally sampled convective updrafts and downdrafts explains a weak friction effect, but not the counter-gradient flux near the cloud tops. The analysis of the momentum flux budget reveals that, in the cloud layer, the counter-gradient flux is driven by convectively triggered nonhydrostatic pressure-gradients and horizontal circulations surrounding the clouds. A model set-up with large domain size and realistic boundary conditions is necessary to resolve these effects.

How Equivalent Are Equivalent Porous Media?

Geoenergy and geoengineering applications usually involve fluid injection into and production from fractured media. Accounting for fractures is important because of the strong poromechanical coupling that ties pore pressure changes and deformation. A possible approach to the problem uses equivalent porous media to reduce the computational cost and model complexity instead of explicitly including fractures in the models. We investigate the validity of this simplification by comparing these two approaches. Simulation results show that pore pressure distribution significantly differs between the two approaches even when both are calibrated to predict identical values at the injection and production wells. Additionally, changes in fracture stability are not well captured with the equivalent porous medium. We conclude that explicitly accounting for fractures in numerical models may be necessary under some circumstances to perform reliable coupled thermohydromechanical simulations, which could be used in conjunction with other tools for induced seismicity forecasting.

InterCarb: A Community Effort to Improve Interlaboratory Standardization of the Carbonate Clumped Isotope Thermometer Using Carbonate Standards

Increased use and improved methodology of carbonate clumped isotope thermometry has greatly enhanced our ability to interrogate a suite of Earth-system processes. However, interlaboratory discrepancies in quantifying carbonate clumped isotope (Δ47) measurements persist, and their specific sources remain unclear. To address interlaboratory differences, we first provide consensus values from the clumped isotope community for four carbonate standards relative to heated and equilibrated gases with 1,819 individual analyses from 10 laboratories. Then we analyzed the four carbonate standards along with three additional standards, spanning a broad range of δ47 and Δ47 values, for a total of 5,329 analyses on 25 individual mass spectrometers from 22 different laboratories. Treating three of the materials as known standards and the other four as unknowns, we find that the use of carbonate reference materials is a robust method for standardization that yields interlaboratory discrepancies entirely consistent with intralaboratory analytical uncertainties. Carbonate reference materials, along with measurement and data processing practices described herein, provide the carbonate clumped isotope community with a robust approach to achieve interlaboratory agreement as we continue to use and improve this powerful geochemical tool. We propose that carbonate clumped isotope data normalized to the carbonate reference materials described in this publication should be reported as Δ47 (I-CDES) values for Intercarb-Carbon Dioxide Equilibrium Scale.

IASI-Derived Sea Surface Temperature Data Set for Climate Studies

Sea surface temperature (SST) is an essential climate variable, that is directly used in climate monitoring. Although satellite measurements can offer continuous global coverage, obtaining a long-term homogeneous satellite-derived SST data set suitable for climate studies based on a single instrument is still a challenge. In this work, we assess a homogeneous SST data set derived from reprocessed Infrared Atmospheric Sounding Interferometer (IASI) level-1 (L1C) radiance data. The SST is computed using Planck's Law and simple atmospheric corrections. We assess the data set using the ERA5 reanalysis and the EUMETSAT-released IASI level-2 SST product. Over the entire period, the reprocessed IASI SST shows a mean global difference with ERA5 close to zero, a mean absolute bias under 0.5°C, with a SD of difference around 0.3°C and a correlation coefficient over 0.99. In addition, the reprocessed data set shows a stable bias and SD, which is an advantage for climate studies. The interannual variability and trends were compared with other SST data sets: ERA5, Hadley Centre's SST (HadISST), and NOAA's Optimal Interpolation SST Analysis (OISSTv2). We found that the reprocessed SST data set is able to capture the patterns of interannual variability well, showing the same areas of high interannual variability (>1.5°C), including over the tropical Pacific in January corresponding to the El Niño Southern Oscillation. Although the period studied is relatively short, we demonstrate that the IASI data set reproduces the same trend patterns found in the other data sets (i.e., cooling trend in the North Atlantic, warming trend over the Mediterranean).

Impact of Remineralization Profile Shape on the Air-Sea Carbon Balance

The ocean's "biological pump" significantly modulates atmospheric carbon dioxide levels. However, the complexity and variability of processes involved introduces uncertainty in interpretation of transient observations and future climate projections. Much research has focused on "parametric uncertainty," particularly determining the exponent(s) of a power-law relationship of sinking particle flux with depth. Varying this relationship's functional form introduces additional "structural uncertainty." We use an ocean biogeochemistry model substituting six alternative remineralization profiles fit to a reference power-law curve, to systematically characterize structural uncertainty, which, in atmospheric pCO2 terms, is roughly 50% of parametric uncertainty associated with varying the power-law exponent within its plausible global range, and similar to uncertainty associated with regional variation in power-law exponents. The substantial contribution of structural uncertainty to total uncertainty highlights the need to improve characterization of biological pump processes, and compare the performance of different profiles within Earth System Models to obtain better constrained climate projections.

The Viscous and Ohmic Damping of the Earth's Free Core Nutation

The cause for the damping of the Earth's free core nutation (FCN) and the free inner core nutation eigenmodes has been a matter of debate since the earliest reliable estimations from nutation observations were made available. Numerical studies are difficult given the extreme values of some of the parameters associated with the Earth's fluid outer core, where important energy dissipation mechanisms can take place. We present a fully 3D numerical model for the FCN capable of describing accurately viscous and Ohmic dissipation processes taking place in the bulk of the fluid core as well as in the boundary layers. We find an asymptotic regime, appropriate for Earth's parameters, where viscous and Ohmic processes contribute to the total damping, with the dissipation taking place almost exclusively in the boundary layers. By matching the observed nutational damping, we infer an enhanced effective viscosity matching and validating methods from previous studies. We suggest that turbulence caused by the Earth's precession can be a source for the enhanced viscosity.

Global Modeled Sinking Characteristics of Biofouled Microplastic

Microplastic debris ending up at the sea surface has become a known major environmental issue. However, how microplastic particles move and when they sink in the ocean remains largely unknown. Here, we model microplastic subject to biofouling (algal growth on a substrate) to estimate sinking timescales and the time to reach the depth where particles stop sinking. We combine NEMO-MEDUSA 2.0 output, that represents hydrodynamic and biological properties of seawater, with a particle-tracking framework. Different sizes and densities of particles (for different types of plastic) are simulated, showing that the global distribution of sinking timescales is largely size-dependent as opposed to density-dependent. The smallest particles we simulate (0.1 μm) start sinking almost immediately around the globe and their trajectories take the longest time to reach their first sinking depth (relative to larger particles). In oligotrophic subtropical gyres with low algal concentrations, particles between 1 and 0.01 mm do not sink within the simulation time of 90 days. This suggests that in addition to the comparatively well-known physical processes, biological processes might also contribute to the accumulation of floating plastic (of 1-0.01 mm) in subtropical gyres. Particles of 1 μm in the gyres start sinking largely due to vertical advection, whereas in the equatorial Pacific they are more dependent on biofouling. The qualitative impacts of seasonality on sinking timescales are small, however, localized sooner sinking due to spring algal blooms is seen. This study maps processes that affect the sinking of virtual microplastic globally, which could ultimately impact the ocean plastic budget.

Scenarios of Human Responses to Unprecedented Social-Environmental Extreme Events

In a rapidly changing world, what is today an unprecedented extreme may soon become the norm. As a result, extreme-related disasters are expected to become more frequent and intense. This will have widespread socio-economic consequences and affect the ability of different societal groups to recover from and adapt to rapidly changing environmental conditions. Therefore, there is the need to decipher the relation between genesis of unprecedented events, accumulation and distribution of risk, and recovery trajectories across different societal groups. Here, we develop an analytical approach to unravel the complexity of future extremes and multiscalar societal responses-from households to national governments and from immediate impacts to longer term recovery. This requires creating new forms of knowledge that integrate analyses of the past-that is, structural causes and political processes of risk accumulation and differentiated recovery trajectories-with plausible scenarios of future environmental extremes grounded in the event-specific literature. We specifically seek to combine the physical characteristics of the extremes with examinations of how culture, politics, power, and policy visions shape societal responses to unprecedented events, and interpret the events as social-environmental extremes. This new approach, at the nexus between social and natural sciences, has the concrete advantage of providing an impact-focused vision of future social-environmental risks, beyond what is achievable within conventional disciplinary boundaries. In this paper, we focus on extreme flooding events and the societal responses they elicit. However, our approach is flexible and applicable to a wide range of extreme events. We see it as the first building block of a new field of research, allowing for novel and integrated theoretical explanations and forecasting of social-environmental extremes.

Single Particle Multipole Expansions From Micromagnetic Tomography

Micromagnetic tomography aims at reconstructing large numbers of individual magnetizations of magnetic particles from combining high-resolution magnetic scanning techniques with micro X-ray computed tomography (microCT). Previous work demonstrated that dipole moments can be robustly inferred, and mathematical analysis showed that the potential field of each particle is uniquely determined. Here, we describe a mathematical procedure to recover higher orders of the magnetic potential of the individual magnetic particles in terms of their spherical harmonic expansions (SHE). We test this approach on data from scanning superconducting quantum interference device microscopy and microCT of a reference sample. For particles with high signal-to-noise ratio of the magnetic scan we demonstrate that SHE up to order n = 3 can be robustly recovered. This additional level of detail restricts the possible internal magnetization structures of the particles and provides valuable rock magnetic information with respect to their stability and reliability as paleomagnetic remanence carriers. Micromagnetic tomography therefore enables a new approach for detailed rock magnetic studies on large ensembles of individual particles.

Buoyancy and Structure of Volatile-Rich Silicate Melts

The early Earth was marked by at least one global magma ocean. Melt buoyancy played a major role for its evolution. Here we model the composition of the magma ocean using a six-component pyrolite melt, to which we add volatiles in the form of carbon as molecular CO or CO2 and hydrogen as molecular H2O or through substitution for magnesium. We compute the density relations from first-principles molecular dynamics simulations. We find that the addition of volatiles renders all the melts more buoyant compared to the reference volatile-free pyrolite. The effect is pressure dependent, largest at the surface, decreasing to about 20 GPa, and remaining roughly constant to 135 GPa. The increased buoyancy would have enhanced convection and turbulence, and thus promoted the chemical exchanges of the magma ocean with the early atmosphere. We determine the partial molar volume of both H2O and CO2 throughout the magma ocean conditions. We find a pronounced dependence with temperature at low pressures, whereas at megabar pressures the partial molar volumes are independent of temperature. At all pressures, the polymerization of the silicate melt is strongly affected by the amount of oxygen added to the system while being only weakly affected by the specific type of volatile added.

Vertically Resolved Magma Ocean-Protoatmosphere Evolution: H2, H2O, CO2, CH4, CO, O2, and N2 as Primary Absorbers

The earliest atmospheres of rocky planets originate from extensive volatile release during magma ocean epochs that occur during assembly of the planet. These establish the initial distribution of the major volatile elements between different chemical reservoirs that subsequently evolve via geological cycles. Current theoretical techniques are limited in exploring the anticipated range of compositional and thermal scenarios of early planetary evolution, even though these are of prime importance to aid astronomical inferences on the environmental context and geological history of extrasolar planets. Here, we present a coupled numerical framework that links an evolutionary, vertically resolved model of the planetary silicate mantle with a radiative-convective model of the atmosphere. Using this method, we investigate the early evolution of idealized Earth-sized rocky planets with end-member, clear-sky atmospheres dominated by either H2, H2O, CO2, CH4, CO, O2, or N2. We find central metrics of early planetary evolution, such as energy gradient, sequence of mantle solidification, surface pressure, or vertical stratification of the atmosphere, to be intimately controlled by the dominant volatile and outgassing history of the planet. Thermal sequences fall into three general classes with increasing cooling timescale: CO, N2, and O2 with minimal effect, H2O, CO2, and CH4 with intermediate influence, and H2 with several orders of magnitude increase in solidification time and atmosphere vertical stratification. Our numerical experiments exemplify the capabilities of the presented modeling framework and link the interior and atmospheric evolution of rocky exoplanets with multiwavelength astronomical observations.

Resource Colimitation Drives Competition Between Phytoplankton and Bacteria in the Southern Ocean

Across the Southern Ocean, phytoplankton growth is governed by iron and light, while bacterial growth is regulated by iron and labile dissolved organic carbon (LDOC). We use a mechanistic model to examine how competition for iron between phytoplankton and bacteria responds to changes in iron, light, and LDOC. Consistent with experimental evidence, increasing iron and light encourages phytoplankton dominance, while increasing LDOC and decreasing light favors bacterial dominance. Under elevated LDOC, bacteria can outcompete phytoplankton for iron, most easily under lower iron. Simulations reveal that bacteria are major iron consumers and suggest that luxury storage plays a key role in competitive iron uptake. Under seasonal conditions typical of the Southern Ocean, sources of LDOC besides phytoplankton exudation modulate the strength of competitive interactions. Continued investigations on the competitive fitness of bacteria in driving changes in primary production in iron-limited systems will be invaluable in refining these results.

The Nature and Origins of Sub-Neptune Size Planets

Planets intermediate in size between the Earth and Neptune, and orbiting closer to their host stars than Mercury does the Sun, are the most common type of planet revealed by exoplanet surveys over the last quarter century. Results from NASA's Kepler mission have revealed a bimodality in the radius distribution of these objects, with a relative underabundance of planets between 1.5 and 2.0R ⊕ . This bimodality suggests that sub-Neptunes are mostly rocky planets that were born with primary atmospheres a few percent by mass accreted from the protoplanetary nebula. Planets above the radius gap were able to retain their atmospheres ("gas-rich super-Earths"), while planets below the radius gap lost their atmospheres and are stripped cores ("true super-Earths"). The mechanism that drives atmospheric loss for these planets remains an outstanding question, with photoevaporation and core-powered mass loss being the prime candidates. As with the mass-loss mechanism, there are two contenders for the origins of the solids in sub-Neptune planets: the migration model involves the growth and migration of embryos from beyond the ice line, while the drift model involves inward-drifting pebbles that coagulate to form planets close-in. Atmospheric studies have the potential to break degeneracies in interior structure models and place additional constraints on the origins of these planets. However, most atmospheric characterization efforts have been confounded by aerosols. Observations with upcoming facilities are expected to finally reveal the atmospheric compositions of these worlds, which are arguably the first fundamentally new type of planetary object identified from the study of exoplanets.