Buoyant melting instabilities beneath extending lithosphere: 1. Numerical models

Buoyant decompression melting instabilities in regions of partially molten upper mantle have been proposed to be an important process that might account for some characteristics of intraplate volcanism on Earth and other terrestrial planets. The instability is driven by variations in the melting rate within a partially molten layer whenever a relative decrease in density accompanies decompression melting of ascending mantle. Here, the development of buoyant decompression melting instabilities in a plane layer of passively upwelling and partially melting mantle beneath diffusely extending lithosphere is studied using numerical convection models covering a wide range of physical parameters. We find that the occurrence and nature of these instabilities in such a scenario is strongly affected by the rate of extension and melt percolation, as well as depth distribution of solid density variations arising from melt depletion. In some cases, instabilities do not occur during extension, but only develop after extension has slowed or stopped completely. This behavior creates two pulses of magma generation due to passive upwelling accompanying extension followed by the subsequent instability and is favored by a faster rate of extension, higher mantle viscosity, higher rate of melt percolation, and smaller amount of solid residuum depletion‐derived buoyancy. Larger degrees of solid density changes accompanying melt depletion can enhance the instability of partially molten mantle during extension but decrease the cumulative volume of generated melt. This kind of behavior modifies the conventional expectation of spatially and temporally correlated volcanism and extension and may lend insight into the observed increase in localized volcanic activity following Miocene Basin and Range extension in the western United States

Buoyant melting instabilities beneath extending lithosphere: 2. Linear analysis

In a companion paper, numerical models reveal that buoyant melting instabilities can occur beneath extending lithosphere for a sufficiently small mantle viscosity, extension rate, and rate of melt percolation. However, in some cases, instabilities do not develop during extension but only occur after extension slows or stops. These results are suggestive of a critical behavior in the onset of these kinds of instabilities and motivate a linear analysis to study the onset of instability in a partially melting, passively upwelling plane layer of mantle beneath extending lithosphere. The model we employ includes the effects of buoyancy arising from thermal expansion, the presence of a retained fraction of partial melt, and depletion of the solid by melt extraction. We find a critical behavior in the onset of instability controlled by melt retention buoyancy that is characterized by a “Rayleigh” number M, such that M must exceed some critical value M_(crit) which depends on the efficiency of Stokes rise of a partially molten body relative to the rate of background percolation. Comparison of this theory to the numerical results in the companion paper yields a close quantitative agreement. We also find that solid depletion buoyancy can either stabilize or destabilize a partially melting layer, depending upon both the distribution of preexisting depletion and the magnitude of density changes with depth. This theory is compared with previous studies of buoyant melting instabilities beneath mid‐ocean ridges where similar behavior was reported, and it suggests that the stability of passively upwelling, partially melting mantle underlying both narrow and wide rift settings is controlled by similar processes

Buoyant melting instabilities beneath extending lithosphere: 2. Linear analysis

In a companion paper, numerical models reveal that buoyant melting instabilities can occur beneath extending lithosphere for a sufficiently small mantle viscosity, extension rate, and rate of melt percolation. However, in some cases, instabilities do not develop during extension but only occur after extension slows or stops. These results are suggestive of a critical behavior in the onset of these kinds of instabilities and motivate a linear analysis to study the onset of instability in a partially melting, passively upwelling plane layer of mantle beneath extending lithosphere. The model we employ includes the effects of buoyancy arising from thermal expansion, the presence of a retained fraction of partial melt, and depletion of the solid by melt extraction. We find a critical behavior in the onset of instability controlled by melt retention buoyancy that is characterized by a “Rayleigh” number M, such that M must exceed some critical value M_(crit) which depends on the efficiency of Stokes rise of a partially molten body relative to the rate of background percolation. Comparison of this theory to the numerical results in the companion paper yields a close quantitative agreement. We also find that solid depletion buoyancy can either stabilize or destabilize a partially melting layer, depending upon both the distribution of preexisting depletion and the magnitude of density changes with depth. This theory is compared with previous studies of buoyant melting instabilities beneath mid‐ocean ridges where similar behavior was reported, and it suggests that the stability of passively upwelling, partially melting mantle underlying both narrow and wide rift settings is controlled by similar processes

Buoyant melting instabilities beneath extending lithosphere: 1. Numerical models

Buoyant decompression melting instabilities in regions of partially molten upper mantle have been proposed to be an important process that might account for some characteristics of intraplate volcanism on Earth and other terrestrial planets. The instability is driven by variations in the melting rate within a partially molten layer whenever a relative decrease in density accompanies decompression melting of ascending mantle. Here, the development of buoyant decompression melting instabilities in a plane layer of passively upwelling and partially melting mantle beneath diffusely extending lithosphere is studied using numerical convection models covering a wide range of physical parameters. We find that the occurrence and nature of these instabilities in such a scenario is strongly affected by the rate of extension and melt percolation, as well as depth distribution of solid density variations arising from melt depletion. In some cases, instabilities do not occur during extension, but only develop after extension has slowed or stopped completely. This behavior creates two pulses of magma generation due to passive upwelling accompanying extension followed by the subsequent instability and is favored by a faster rate of extension, higher mantle viscosity, higher rate of melt percolation, and smaller amount of solid residuum depletion‐derived buoyancy. Larger degrees of solid density changes accompanying melt depletion can enhance the instability of partially molten mantle during extension but decrease the cumulative volume of generated melt. This kind of behavior modifies the conventional expectation of spatially and temporally correlated volcanism and extension and may lend insight into the observed increase in localized volcanic activity following Miocene Basin and Range extension in the western United States

Persistence of strong silica-enriched domains in the Earth's lower mantle

The composition of the lower mantle—comprising 56% of Earth’s volume—remains poorly constrained. Among the major elements, Mg/Si ratios ranging from ∼0.9–1.1, such as in rocky Solar-System building blocks (or chondrites), to ∼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 homogenous 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 associated with a 20-fold change in viscosity, such as due to the dominance of intrinsically strong (Mg, Fe)SiO3–bridgmanite in low-Mg/Si domains, is sufficient to prevent efficient mantle mixing, even on large scales. Models predict that intrinsically strong domains stabilize mantle convection patterns, and coherently persist at depths of about 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 centres, and geophysical changes in seismic-tomography patterns, radial viscosity, rising plumes and sinking slabs near 1,000 km depth. Moreover, these ancient structures may provide a reservoir to host primordial geochemical signatures

Some dynamical consequences of partial melting in Earth's deep mantle

International audienceWe use regional scale numerical models of mantle convection to investigate the simple hypothesis that seismically anomalous thin patches of Earth's lowermost mantle, termed “ultralow-velocity zones” or ULVZ, are derived from partial melting of ordinary mantle. The models span the lower 500 km of Earth's mantle, employ both temperature and melting-related contributions to buoyancy, and include a cold randomly moving downwelling introduced from above to maintain a thermal boundary layer. Partial melting of ordinary mantle introduces a ubiquitous partially molten layer above an isothermal core–mantle boundary as a consequence of its isothermal and isobaric conditions, although it naturally develops variations in thickness greater than two orders of magnitude, with the thickest portions occurring at the base of upwelling plumes and a thin layer elsewhere. We find that only a dense partially molten mixture produces partial melt distributions that are compatible with seismic observations of ULVZ; however, if such a melt percolates downward a dense basal liquid layer accumulates above the core–mantle boundary. The apparent requirement of a volumetrically dense and non-percolating melt phase in the lowermost mantle presents serious problems for the hypothesis that ULVZ arise from melting of ordinary mantle, and suggests that such features likely form as a consequence of more complex processes. Furthermore, these considerations suggest that the solidus of ordinary mantle is a reasonable upper bound on the present day temperature of the CMB