The giant impact hypothesis for Moon formation successfully explains the dynamic properties of the Earth-Moon system but remains challenged by the similarity of isotopic fingerprints of the terrestrial and lunar mantles. Moreover, recent geochemical evidence suggests that the Earth's mantle preserves ancient (or "primordial") heterogeneity that predates the Moon-forming giant impact. Using a new hydrodynamical method, we here show that Moon-forming giant impacts lead to a stratified starting condition for the evolution of the terrestrial mantle. The upper layer of the Earth is compositionally similar to the disk, out of which the Moon evolves, whereas the lower layer preserves proto-Earth characteristics. As long as this predicted compositional stratification can at least partially be preserved over the subsequent billions of years of Earth mantle convection, the compositional similarity between the Moon and the accessible Earth's mantle is a natural outcome of realistic and high-probability Moon-forming impact scenarios. The preservation of primordial heterogeneity in the modern Earth not only reconciles geochemical constraints but is also consistent with recent geophysical observations. Furthermore, for significant preservation of a proto-Earth reservoir, the bulk composition of the Earth-Moon system may be systematically shifted towards chondritic values.

The composition of the lower mantle $-$ comprising 56% of Earth's volume $-$ remains poorly constrained. Among the major elements, Mg/Si ratios ranging from $\sim$0.9$-$1.1, such as in rocky solar-system building blocks (or chondrites), to $\sim$1.2$-$1.3, such as in upper-mantle rocks (or pyrolite), have been proposed. Geophysical evidence for subducted lithosphere deep in the mantle has been interpreted in terms of efficient mixing and thus homogeneous Mg/Si across most of the mantle. However, previous models did not consider the effects of variable Mg/Si on the viscosity and mixing efficiency of lower-mantle rocks. Here, we use geodynamic models to show that large-scale heterogeneity with viscosity variations of $\sim$20$\times$, such as due to the dominance of intrinsically strong (Mg,Fe)SiO$_3-$bridgmanite in low-Mg/Si domains, are sufficient to prevent efficient mantle mixing, even on large scales. Models predict that intrinsically strong domains stabilize degree-two mantle-convection patterns, and coherently persist at depths of $\sim$1,000$-$2,200 km up to the present-day, separated by relatively narrow up-/downwelling conduits of pyrolitic material. The stable manifestation of such "bridgmanite-enriched ancient mantle structures" (BEAMS) may reconcile the geographical fixity of deep-rooted mantle-upwelling centers, and fundamental geophysical changes near 1,000 km depth (e.g. in terms of seismic-tomography patterns, radial viscosity increase, lateral deflections of rising plumes and sinking slabs). Moreover, these ancient structures may provide a reservoir to host primordial geochemical signatures.

SiO$_2$ may have been expelled from the core directly following core formation in the early stages of Earth's accretion and onwards through the present day. On account of SiO$_2$'s low density with respect to both the core and the lowermost mantle, we examine the process of SiO$_2$ accumulation at the core-mantle boundary (CMB) and its incorporation into the mantle by buoyant rise. Today, if SiO$_2$ is 100-10000 times more viscous than lower mantle material, the dimensions of SiO$_2$ diapirs formed by the viscous Rayleigh-Taylor instability at the CMB would cause them to be swept into the mantle as inclusions of 100 m - 10 km diameter. Under early Earth conditions of rapid heat loss after core formation, SiO$_2$ diapirs of ~1 km diameter could have risen independently of mantle flow to their level of neutral buoyancy in the mantle, trapping them there due to a combination of intrinsically high viscosity and neutral buoyancy. We examine the SiO$_2$ yield by assuming Si+O saturation at the conditions found at the base of a magma ocean and find that for a range of conditions, dispersed bodies could reach as high as 8.5 vol.% in parts of the lower mantle. At such low concentration, their effect on aggregate seismic wavespeeds is within observational seismology uncertainty. However, their presence can account for small-scale scattering in the lower mantle due to the bodies' large velocity contrast. We conclude that the shallow lower mantle (700-1500 km depth) could harbor SiO$_2$ released in early Earth times.

Terrestrial planets are thought to experience episode(s) of large-scale melting early in their history. Fractionation during magma-ocean freezing leads to unstable stratification within the related cumulate layers due to progressive iron enrichment upward, but the effects of incremental cumulate overturns during MO crystallization remain to be explored. Here, we use geodynamic models with a moving-boundary approach to study convection and mixing within the growing cumulate layer, and thereafter within the fully crystallized mantle. For fractional crystallization, cumulates are efficiently stirred due to subsequent incremental overturns, except for strongly iron-enriched late-stage cumulates, which persist as a stably stratified layer at the base of the mantle for billions of years. Less extreme crystallization scenarios can lead to somewhat more subtle stratification. In any case, the long-term preservation of at least a thin layer of extremely enriched cumulates with Fe#>0.4, as predicted by all our models, is inconsistent with seismic constraints. Based on scaling relationships, however, we infer that final-stage Fe-rich magma-ocean cumulates originally formed near the surface should have overturned as small diapirs, and hence undergone melting and reaction with the host rock during sinking. The resulting moderately iron-enriched metasomatized/hybrid rock assemblages should have accumulated at the base of the mantle, potentially fed an intermittent basal magma ocean, and be preserved through the present-day. Such moderately iron-enriched rock assemblages can reconcile the physical properties of the large low shear-wave velocity provinces in the presentday lower mantle. Thus, we reveal Hadean melting and rock-reaction processes by integrating magmaocean crystallization models with the seismic-tomography snapshot.