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.

Origins of ultralow velocity zones through slab-derived metallic melt

Understanding the ultralow velocity zones (ULVZs) places constraints on the chemical composition and thermal structure of deep Earth and provides critical information on the dynamics of large-scale mantle convection, but their origin has remained enigmatic for decades. Recent studies suggest that metallic iron and carbon are produced in subducted slabs when they sink beyond a depth of 250 km. Here we show that the eutectic melting curve of the iron-carbon system crosses the current geotherm near Earth's core-mantle boundary, suggesting that dense metallic melt may form in the lowermost mantle. If concentrated into isolated patches, such melt could produce the seismically observed density and velocity features of ULVZs. Depending on the wetting behavior of the metallic melt, the resultant ULVZs may be short-lived domains that are replenished or regenerated through subduction, or long-lasting regions containing both metallic and silicate melts. Slab-derived metallic melt may produce another type of ULVZ that escapes core sequestration by reacting with the mantle to form iron-rich postbridgmanite or ferropericlase. The hypotheses connect peculiar features near Earth's core-mantle boundary to subduction of the oceanic lithosphere through the deep carbon cycle.

Structure and density of basaltic melts at mantle conditions from first-principles simulations

The origin and stability of deep-mantle melts, and the magmatic processes at different times of Earth's history are controlled by the physical properties of constituent silicate liquids. Here we report density functional theory-based simulations of model basalt, hydrous model basalt and near-MORB to assess the effects of iron and water on the melt structure and density, respectively. Our results suggest that as pressure increases, all types of coordination between major cations and anions strongly increase, and the water speciation changes from isolated species to extended forms. These structural changes are responsible for rapid initial melt densification on compression thereby making these basaltic melts possibly buoyantly stable at one or more depths. Our finding that the melt-water system is ideal (nearly zero volume of mixing) and miscible (negative enthalpy of mixing) over most of the mantle conditions strengthens the idea of potential water enrichment of deep-mantle melts and early magma ocean.

Probing failure susceptibilities of earthquake faults using small-quake tidal correlations

Mitigating the devastating economic and humanitarian impact of large earthquakes requires signals for forecasting seismic events. Daily tide stresses were previously thought to be insufficient for use as such a signal. Recently, however, they have been found to correlate significantly with small earthquakes, just before large earthquakes occur. Here we present a simple earthquake model to investigate whether correlations between daily tidal stresses and small earthquakes provide information about the likelihood of impending large earthquakes. The model predicts that intervals of significant correlations between small earthquakes and ongoing low-amplitude periodic stresses indicate increased fault susceptibility to large earthquake generation. The results agree with the recent observations of large earthquakes preceded by time periods of significant correlations between smaller events and daily tide stresses. We anticipate that incorporating experimentally determined parameters and fault-specific details into the model may provide new tools for extracting improved probabilities of impending large earthquakes.