Melt Migration in a Heterogeneous Mantle
The generation and segregation of melt from the earth's interior is a fundamental, but only partially understood, process. Melt segregation involves interacting physical processes on a wide range of length scales. Melt, generated at the millimeter scale of individual mineral grains, migrates through porous solid mantle, perhaps in an interconnected network of channels that forms at scales of tens of meters. This melt, driven by buoyancy forces and stresses that develop in solid mantle, flows on scales of hundreds of kilometers (Figure 1). It has been suggested that high-porosity pyroxene-free dunite channels may serve as pathways for efficient extraction of melt formed in the deep mantle. With postdocs Alan Schemenz and Marc Hesse, and Profs. Marc Parmentier and Jan Hesthaven (Applied Math, Brown University), we examined the dynamics of reactive dissolution along a solubility gradient in an upwelling and viscously deformable mantle column. Three dynamics regimes are identified through linear stability analysis (Hesse et al., 2011): channel regime where melt localizes into high porosity channels, wave regime where melt migrates in the form of compaction-dissolution waves, and a stable regime. The compaction-dissolution waves are different from the solitary waves discussed in the magma dynamics literature in that the wave speed of the former is smaller than that of the buoyancy-driven melt migration. A parameter study shows that compaction-dissolution waves may play an important role for melt migration beneath mid-ocean ridges.
Brown faculty collaborators:Jan Hesthaven
Other project collaborators:Marc Hesse, Alan Schiemenz, Seshu Tirupathi, Chenguang Sun, Conroy Baltzell
Figure 1. Multiscales of melt generation and segregation in the mantle. The largest scale defines the setting in which melting occurs. At the smallest scale partial melting occurs by chemical reactions between different mineral grains. The intermediate scale that controls the rate of melt transport can be broadly described as a porous block with multiple scales of porosity. High porosity melt channels consisting primarily of the mineral olivine are thought to be a consequence of rising melt dissolving minerals clinopyroxene (cpx) and orthopyroxene (opx) while precipitating olivine. Here D refers to the rock dunite that comprises almost entirely of olivine with small amount of cpx and spinel. H represents harzburgite composed of olivine + opx ± cpx ± spinel.
Figure 2. Steady-state distributions of (a) porosity and (d) the soluble mineral opx in a short upwelling column. (b) Close-up view of the porosity field in the highlighted area near the top of the wider channel. (c) Close-up view of the opx field in the highlighted region near the bottom of the wider dunite channel. Dark blue regions in Figure 1d mark the opx-free dunite channels. White curves in Figures 1b and 1c are streamlines of the melt. The two dunite channels were initiated by two small but sustained Gaussian perturbations in porosity at the bottom of the upwelling column. Steady state is established after one solid overturn (t ≥ 2). Porosity and opx modal abundance are normalized to their respective background values (2% and 15%, respectively) at the bottom of the upwelling column.
Figure 3. Distributions of (ad) porosity and (eh) the soluble mineral opx at four selected times in an upwelling column. Dark blue regions in Figures 2f2h mark the opx-free dunite channels. Curves in Figures 2e2h are melt streamlines. Dimensions of the waves are outlined by the dashed rectangles in Figures 2b and 2d. Porosity and opx abundance are normalized to their background values (1% and 20%, respectively) at the bottom of the upwelling column. Time is measured in units of solid upwelling time which is the time needed to advance the solid by one compaction length. One solid overturn takes 6.75 time units in this example.
Figure 4. Distributions of (ac) 86Sr/87Sr and (df) 143Nd/144Nd isotopic ratios in the melt and solid at three selected times in an upwelling column. The isotopic ratios were calculated using the porosity, opx, and velocity fields from the case shown in Figure 3 while neglecting diffusion and dispersion in the melt. The white lines in each panel are melt streamlines.Back to Yan Liang's Brown Research Profile