---
title: Stratified Lunar Magma Ocean
url: https://www.emergentmind.com/topics/stratified-lunar-magma-ocean
type: topic
---

# Stratified Lunar Magma Ocean

A stratified lunar magma ocean is the chemically and density layered state produced as the Moon’s early global magma ocean crystallized, generating compositionally distinct cumulates and residual products rather than a homogeneous mantle. In the lunar literature, the term refers both to the primary sequence of olivine- and pyroxene-dominated early cumulates, plagioclase flotation products, late Fe–Ti-rich cumulates, KREEP-type residual materials, and an anorthositic crust, and to the subsequent partial overturn, retention, and remelting of those layers that shaped the present mantle and late mare volcanism [1809.01629][2104.06847].

## 1. Fractional crystallization and primary cumulate layering

The canonical stratified-LMO framework treats crystallization as a nearly fractional, batch-fractional, or mixed equilibrium–fractional process following liquidus–solidus phase relations for a chondritic bulk lunar composition. The residual melt evolves along a magma-ocean adiabat, and the identity and proportions of crystallizing phases are read from experimental or thermodynamic phase diagrams; densities and related physical properties are then computed from phase equations of state and combined by mass or volume fractions [1809.01629].

The crystallization sequence inferred for the LMO is strongly ordered. Early Mg-rich olivine and then orthopyroxene crystallize and sink, producing a lower cumulate pile with only small density gradients. Plagioclase appears later and, being less dense than the melt, floats to form the anorthositic crust. Near the end of solidification, progressive Fe enrichment of the residual melt leads to ilmenite-bearing cumulates that are denser than the underlying mantle and therefore gravitationally unstable. In related mantle-stratigraphy models, the final residual melts are KREEP-type cumulates and the last flotation products are the anorthositic crust [1809.01629][2104.06847].

| Stage of solidification | Dominant phase(s) | Dynamical outcome |
|---|---|---|
| 0–20% | Mg-rich olivine (Fo\(_{90}\)) | Sinks; builds early lower cumulates |
| 20–70% | Orthopyroxene | Sinks; enlarges lower cumulate pile |
| 70–80% | Plagioclase (An\(_{98-100}\)) | Floats; forms anorthositic lid/crust |
| 80–95% | Fe–Ti oxides, including ilmenite | Produces dense late cumulates |
| Final residual melt | KREEP-type cumulates | Concentrates incompatible elements |

This sequence yields a first-order vertical stratigraphy in which buoyant plagioclase-derived crust overlies mafic cumulates, while the densest late products tend to migrate downward. The resulting mantle is therefore stratified both chemically and rheologically, not merely compositionally.

## 2. Mantle reservoirs, FeO dependence, and present-day structure

A detailed cumulate-reservoir formulation divides the mantle into early-stage lower-mantle cumulates dominated by olivine \(\pm\) orthopyroxene, mid-stage upper-mantle clinopyroxene-rich cumulates, late-stage ilmenite-bearing cumulates (IBC), KREEP-type residual cumulates, and an anorthositic crust fixed at 40 km in the model configuration. The lower-mantle cumulates are Fe-poor and Ti-poor, the upper mantle is modestly more Fe rich, and the IBC are the densest and most Fe–Ti enriched of the mantle reservoirs [2104.06847].

The properties of these reservoirs depend strongly on bulk silicate Moon FeO. For the lower mantle,
\[
\rho_{\rm LM} \simeq 3380~{\rm kg/m^3} + (11~{\rm kg/m^3~per~wt\%FeO_{BSM}})\cdot({\rm FeO}_{\rm BSM}-8),
\]
with \({\rm TiO_2}_{\rm LM}\lesssim 0.1\). For the upper mantle,
\[
\rho_{\rm UM} \simeq 3340~{\rm kg/m^3} + (11~{\rm kg/m^3~per~wt\%FeO_{BSM}})\cdot({\rm FeO}_{\rm BSM}-8),
\]
with \({\rm TiO_2}_{\rm UM}\lesssim 0.2\). The IBC obey the strongest FeO dependence:
\[
t_{\rm IBC}({\rm km}) \simeq 3.9\,({\rm FeO}_{\rm BSM}) - 19.0,
\]
so that \(t_{\rm IBC}\approx 11.6\) km at \(8\) wt% FeO and \(t_{\rm IBC}\approx 31.1\) km at \(13\) wt% FeO. Their density at formation depth is
\[
\rho_{\rm IBC} \simeq 3750~{\rm kg/m^3} + (14~{\rm kg/m^3~per~wt\%FeO_{BSM}})\cdot({\rm FeO}_{\rm BSM}-8).
\]
KREEP cumulates are volumetrically minor, with thickness \(\lesssim 5\) km and density \(\simeq 3200\)–\(3300~{\rm kg/m^3}\) [2104.06847].

These relations are used together with selenodetic and seismic constraints. The same modeling indicates that the bulk silicate Moon contains about \(9.4\)–\(10.9\) weight percent FeO, with a lowermost limit of \(8.3\) weight percent and an uppermost limit of \(11.9\) weight percent. It also implies incomplete sinking of the IBC: only \(20\)–\(60\) percent of the original IBC mass might have reached the core-mantle boundary, while the rest either remained near its formation depth or was mixed into the middle mantle. In that interpretation, the present-day mantle retains both primary fractionation signatures and overturn-generated heterogeneity, including a high-density, low-velocity zone at the base of the mantle of about \(170\) km thickness [2104.06847].

## 3. Fe–Mg isotopes and the Chang’E 5 evidence for hybrid cumulate remelting

Direct evidence that late lunar volcanism sampled stratified LMO reservoirs comes from the Chang’E 5 mare basalts, whose Fe and Mg isotopic compositions, petrology, and trace-element systematics were interpreted as recording a hybrid mantle source. The isotope notation is
\[
\delta^{56}{\rm Fe} \equiv \left[\frac{(^{56}{\rm Fe}/^{54}{\rm Fe})_{\rm sample}}{(^{56}{\rm Fe}/^{54}{\rm Fe})_{\rm IRMM\text{-}014}} - 1\right]\times 10^3,
\]
and
\[
\delta^{26}{\rm Mg} \equiv \left[\frac{(^{26}{\rm Mg}/^{24}{\rm Mg})_{\rm sample}}{(^{26}{\rm Mg}/^{24}{\rm Mg})_{\rm DSM\text{-}3}} - 1\right]\times 10^3.
\]
Measured CE-5 values are \(\delta^{56}{\rm Fe}=+0.161\pm0.010\)‰ and \(\delta^{26}{\rm Mg}=-0.264\pm0.015\)‰. Modeled or previously measured LMO-cumulate endmembers are \(\delta^{56}{\rm Fe}\approx 0\)‰ and \(\delta^{26}{\rm Mg}\approx -0.30\)‰ for early-stage olivine-dominated cumulates, versus \(\delta^{56}{\rm Fe}\approx +0.30\)‰ and \(\delta^{26}{\rm Mg}\approx -0.45\)‰ for late-stage clinopyroxene + ilmenite-rich cumulates. The CE-5 Fe isotope value lies between low-Ti and high-Ti mare basalt fields, whereas the Mg isotope value remains within the range of early-cumulate melts, pointing to a mixed source rather than a single cumulate lithology [2302.10507].

Petrologically, the CE-5 clasts have Mg# \(\simeq 29\), the lowest of all returned mare basalts, and display a sub-ophitic groundmass of pyroxene + olivine + plagioclase + minor ilmenite. They also contain symplectites of olivine + silica or olivine + Si–K glass derived from pyroxferroite breakdown, together with “swiss-cheese” intergrowths of fayalite + Si–K glass. Olivine compositions of Fo\(_{52-58}\), in equilibrium with the whole-rock Mg#, indicate no significant crystal accumulation, so the bulk clast chemistry approximates a primary liquid. Their trace-element pattern is KREEP-like, with \({\rm La}_N=158\), \({\rm La/Yb}_N\approx 2.8\), and a deep Eu anomaly of \({\rm Eu/Eu^*}\approx 0.5\), which was argued to require a source enriched in late-stage plagioclase-depleted melt rather than simple olivine fractionation alone [2302.10507].

The isotopic mixing relation is written
\[
\delta X_{\rm mix}=f\,\delta X_{\rm early}+(1-f)\,\delta X_{\rm late},
\]
where \(f\) is the mass fraction of the early-cumulate melt. Substituting the CE-5 \(\delta^{56}{\rm Fe}\) and \(\delta^{26}{\rm Mg}\) values and the endmember compositions yields \(f\approx 0.70\)–\(0.80\), corresponding to a \(20\)–\(30\) percent contribution of late-stage cumulate melt. In the specific geometric interpretation given for an \(\sim 600\) km deep LMO, this is reproduced by about \(70\%\) melt of a \(50\%\)-solid cumulate, mainly olivine + orthopyroxene, and \(30\%\) melt of a \(98\%\)-solid cumulate, clinopyroxene + ilmenite rich, originating at about \(300\) km and \(\ll 100\) km depth, respectively. Fractional crystallization or impact-related evaporation was argued not to generate the observed combined isotopic and geochemical signature [2302.10507].

## 4. Overturn physics, mixing, and the survival of stratification

The persistence or destruction of LMO layering depends on the competition between thermal convection, compositional buoyancy, chemical diffusion, and rheology. A standard density expansion used in magma-ocean and cumulate models is
\[
\rho(P,T,X_i)=\rho_0+\left(\frac{\partial \rho}{\partial P}\right)_T(P-P_0)-\alpha\rho_0(T-T_0)+\sum_i\left(\frac{\partial \rho}{\partial X_i}\right)(X_i-X_{i0}),
\]
with representative values \(\alpha\sim 3\times 10^{-5}~{\rm K^{-1}}\) and \(\left(\partial \rho/\partial P\right)_T\sim 10^{-6}~{\rm bar^{-1}}\). The vigor of overturn is commonly parameterized by
\[
Ra=\frac{\alpha g \Delta T d^3}{\kappa \nu},
\]
where \(g\simeq 1.6~{\rm m\,s^{-2}}\), \(\kappa\simeq 10^{-6}~{\rm m^2\,s^{-1}}\), and convection requires \(Ra>Ra_c\simeq 10^3\)–\(10^4\). For a compositional density contrast \(\Delta\rho\simeq 400~{\rm kg\,m^{-3}}\), the buoyancy flux per unit area is
\[
B=\frac{g\,\Delta\rho}{\rho_0}\,\kappa.
\]
Chemical diffusion follows
\[
\frac{\partial C}{\partial t}=D\nabla^2C,\qquad D=D_0\exp\!\left(-\frac{E_a}{RT}\right),
\]
with \(D_0\sim 10^{-9}\)–\(10^{-10}~{\rm m^2\,s^{-1}}\) and \(E_a\sim 200~{\rm kJ\,mol^{-1}}\), while rheology may be written in Arrhenius form as
\[
\dot\varepsilon=A\,\sigma^n\exp\!\left(-\frac{E_v}{RT}\right),
\]
with \(E_v\sim 200\)–\(400~{\rm kJ\,mol^{-1}}\) [1809.01629].

These relations support a time-dependent rather than static view of stratification. Rapid initial cooling can crystallize the magma ocean to \(\sim 80\%\) by volume in \(t_0\simeq 10^3\) yr, whereas formation of the anorthositic lid delays the final \(\sim 20\%\) by an extra \(\sim 10\) Myr. In models including strong tidal heating, full solidification can extend to \(\sim 200\) Myr. A \(10\) km thick dense ilmenite layer has a Stokes-type settling time of \(\simeq 10^5\)–\(10^6\) yr for cumulus viscosities of \(\eta_{\rm cumulus}\simeq 10^{17}\)–\(10^{20}\) Pa s, and if solid-state convection starts early, mixing can reduce heterogeneity over \(\sim 10^7\)–\(10^8\) yr, although the ilmenite-rich layer may persist as a dense root [1809.01629].

The CE-5 model places partial remelting within this dynamical setting. Dense late cumulates sink or overturn into a still-warm mantle, lower the local solidus, and permit partial melting of mixtures of early olivine-rich and late clinopyroxene–ilmenite cumulates. Melt compositions are then governed by batch melting,
\[
C_{\rm melt}= \frac{C_{\rm source}}{D+F(1-D)},
\]
where \(F\) is melt fraction and \(D=\sum X_iD_i\) is the bulk partition coefficient. Isotopic evolution during equilibrium crystallization is described using Rayleigh-type relations such as
\[
\delta_{\rm melt}=\left[\frac{\delta_0+1000}{f_{\rm melt}}\right]-1000,
\]
with equilibrium fractionation factors for Fe scaling approximately as \(1/T^2\). In this framework, stratification is both preserved and reworked: preserved in the sense that late cumulates remain compositionally distinct, and reworked in the sense that overturn and remelting generate hybrid magmas [2302.10507].

## 5. A two-layer proto-Earth–Theia stratification model

A separate model extends the idea of a stratified LMO back to the immediate post-impact state by proposing that the magma ocean inherited two dominantly different silicate melts, one from proto-Earth’s mantle and one from Theia’s, and that these unmixed by density before solidification. In that scenario, the newborn Moon was effectively molten to a depth up to \(\sim 1000\)–\(1200\) km, corresponding to \(0\)–\(3.8\) GPa, and disk or lunar temperatures were \(\simeq 2000\)–\(5000\) K. The proto-Earth silicate melt is represented by \(44.9\) wt% \({\rm SiO_2}\), \(37.7\) wt% \({\rm MgO}\), \(7.93\) wt% \({\rm FeO}\), and \(4.44\) wt% \({\rm Al_2O_3}\), whereas the Theia melt is represented by \(43.90\) wt% \({\rm SiO_2}\), \(31.66\) wt% \({\rm MgO}\), \(18.8\) wt% \({\rm FeO}\), and \(3.15\) wt% \({\rm Al_2O_3}\) [2509.06519].

Because the assumed Theia mantle is more Fe rich, its melt is denser: \(\rho_{\rm Theia}\approx 2.78~{\rm g\,cm^{-3}}\) at \(0\) GPa and \(2000^\circ\)C, compared with \(\rho_{\rm proto\mbox{-}Earth}\approx 2.67~{\rm g\,cm^{-3}}\). For a molten Theia-derived blob of radius \(r\sim 1\) km, equating drag and buoyancy gives
\[
F_d=\frac{1}{2}C_d\rho_E\pi r^2 v^2,\qquad 
F_g=\frac{4}{3}\pi r^3(\rho_T-\rho_E)g,
\]
and hence
\[
v\simeq \sqrt{\frac{8r(\rho_T-\rho_E)g}{3C_d\rho_E}}\simeq 20~{\rm m\,s^{-1}}.
\]
At lunar radius \(R\simeq 1700\) km, the inferred sinking time is \(t_{\rm sink}\simeq R/v\simeq 23\) h. By contrast, diffusion in silicate melt with \(D\sim 10^{-11}~{\rm m^2\,s^{-1}}\) over \(L\sim 1\) m gives
\[
t_{\rm mix}\simeq \frac{L^2}{D}\sim 3\times 10^3~{\rm yr},
\]
much longer than \(t_{\rm sink}\). The model therefore predicts prompt segregation into a dense lower Theia-derived melt layer and an upper proto-Earth melt layer [2509.06519].

Fractional crystallization of the two layers was modeled with alphaMELTS (pMELTS) at \(1.5\), \(2.0\), \(2.5\), and \(3.0\) GPa. At \(2\) GPa, orthopyroxene and clinopyroxene from the upper layer are less dense than the lower melt and float upward across the interface; by the time these phases are exhausted, \(\sim 27\%\) of the upper-layer mass sits atop the interface as solids. Between \(1549^\circ\)C and \(1240^\circ\)C, olivine from the upper melt also floats, adding \(\sim 16\%\) of the upper-layer mass. In the lower layer, olivine is neutrally buoyant at \(\sim 3\) GPa, but orthopyroxene above \(1531^\circ\)C and clinopyroxene down to \(\sim 1258^\circ\)C are less dense than the lower melt and float into the upper cumulate pile, contributing \(\sim 13\%\) of the lower-layer mass. In this proposed solution to the lunar isotopic crisis, the final solid Moon has an upper cumulate mantle derived almost exclusively from proto-Earth material, while isotopically distinct Theia material is sequestered in the lower mantle and core [2509.06519].

## 6. Synthesis, misconceptions, and active points of interpretation

The available stratified-LMO models agree that lunar magma-ocean solidification generated sharp chemical and density gradients, but they differ on the extent to which those gradients survived intact. Reviews and interior-structure models emphasize a bottom-up fractional-crystallization sequence followed by partial overturn, incomplete IBC sinking, and retention of major lithological boundaries; the CE-5 basalt study emphasizes later remelting of mixed early and late cumulates; and the two-layer isotopic-crisis model adds a pre-solidification density separation between proto-Earth and Theia-derived melts [1809.01629][2104.06847][2302.10507][2509.06519].

One recurrent misconception is that stratification implies a permanently static stack of pristine layers. The dynamical calculations instead indicate rapid crystallization, conductive delay by an anorthositic lid, gravitational instability of late dense cumulates, and mixing or entrainment over \(10^7\)–\(10^8\) yr, even though dense ilmenite-bearing material may remain concentrated near the core-mantle boundary. Another misconception is that KREEP-like signatures in young mare basalts require large admixture of KREEP. The CE-5 interpretation argues that a hybrid mantle source containing both early- and late-stage LMO cumulates can reproduce the combination of low Mg#, intermediate \(\delta^{56}{\rm Fe}\), early-cumulate-like \(\delta^{26}{\rm Mg}\), and KREEP-like REE systematics without requiring large admixture of KREEP [2302.10507].

Taken together, these results suggest that “stratified lunar magma ocean” is best understood as a coupled petrological and dynamical framework. In its narrow sense, it denotes the phase-equilibria-driven cumulate architecture generated during bottom-up solidification. In a broader sense, it includes the flotation of plagioclase, the sinking and partial survival of IBC, the formation and preservation of KREEP-rich reservoirs, the possibility of compositional two-layering inherited from the Moon-forming impact, and the later selective remelting of mixed cumulate packages. Within that framework, the anorthositic crust, mantle seismic layering, basal dense material, the Procellarum KREEP Terrane, and the geochemistry of the youngest mare basalts are all treated as linked consequences of early lunar differentiation rather than isolated phenomena [1809.01629][2104.06847].

Source: https://www.emergentmind.com/topics/stratified-lunar-magma-ocean