Small-scale mineralogical heterogeneity from variations in phase assemblages in the transition zone and D″ layer predicted by convection modelling

Abstract

Small-scale heterogeneity in the deep mantle is concentrated in the upper-mantle transition zone (TZ), in the depth range 410–660 km and also at the bottom 250 km D″ region. This encourages a more detailed investigation of the potential for seismic reflectivity imaging by modelling heterogeneous structures in mantle convection models including phase transitions of the TZ and D″ regions. We applied finite elements with variable spacing near the boundary layers in 2-D cylindrical geometry that allow for sufficient spatial resolution. We investigated several models including an extended Boussinesq (EBA) model, focused on the D″ region, and a compressible (ALA) model for the TZ region. The results for the D″ region show typical lens-shaped structures of post-perovskite (PPV) embedded in the perovskite (PV) background mantle, where the thickness of the lenses, at 200–400 km, strongly depends on the Clapeyron slope of the PV-PPV transition. A second phase transition (double crossing) occurs in case the core temperature is higher than the intercept temperature T _i. Our phase-dependent rheology results in contrasting effective viscosity between PV and PPV. Our model results reveal a distinctly clear mechanical weakening of the PPV lenses with about an order of magnitude lower viscosity. The shear wave-speed distributions computed from our convection results are strongly correlated with the heterogeneous distribution of the mineral phase. Gradients in the seismic wave-speed that are the target of seismological reflectivity imaging are clearly revealed. The wave-speed results show a clear resolution of the top and bottom interfaces of the PPV lenses. Our ALA model for the TZ is based on a thermodynamical model for the magnesium end-member of an olivine-pyroxene mantle. The model predicts a much more complex distribution of mineral phases, compared to our D″ results, in agreement with the greater number of mineral phases involved in the olivine-pyroxene phase diagram for the P, T conditions of the transition zone. Near cold downwelling flows representing subducting lithospheric slabs, where the local geotherm can differ by up to 1 000 K from the horizontal average, and small-scale lateral variations in the mineral phases can occur.

Keywords

mantle convection / phase transition / post-perovskite / mineralogical heterogeneity

References

[1] Ammann M. W., Brodholt J. P., Wookey J., . First-Principles Constraints on Diffusion in Lower-Mantle Minerals and a Weak D″ Layer. Nature, 2010, 465(7297): 462-465.

[2] Boehler R.. High-Pressure Experiments and the Phase Diagram of Lower Mantle and Core Materials. Reviews of Geophysics, 2000, 38(2): 221-245.

[3] Cadek O., Fleitout L.. Effect of Lateral Viscosity Variations in the Core-Mantle Boundary Region on Predictions of the Long-Wavelength Geoid. Stud. Geophys. Geod., 2006, 50(2): 217-232.

[4] Cao Q., Wang P., van der Hilst R. D., . Imaging the Upper Mantle Transition Zone with a Generalized Radon Transform of SS Precursors. Phys. Earth Planet. Inter., 2010, 180(1–2): 80-91.

[5] Catalli K., Shim S. H., Prakapenka V.. Thickness and Claeyron Slope of the Post-Perovskite Boundary. Nature, 2009, 462(7274): 782-785.

[6] Christensen A. U. R., Yuen D. A.. Layered Convection Induced by Phase Transitions. J. Geophys. Res., 1985, 90: 10291-10300.

[7] Connolly J. A. D.. Computation of Phase Equilibria by Linear Programming: A Tool for Geodynamic Modeling and Its Application to Subduction Zone Decarbonation. Earth Planet. Sci. Lett., 2005, 236(1–2): 524-541.

[8] Daessler R., Yuen D. A.. The Effects of Phase Transition Kinetics on Subducting Slabs. Geophys. Res. Lett., 1993, 20(23): 2603-2606.

[9] de Hoop M. V., Smith H., Uhlmann G., . Seismic Imaging with the Generalized Radon Transform: A Curvelet Transform Perspective. Inverse Problems, 2009, 25 2 025005

[10] Hernlund J. W., Thomas C., Tackley P. J.. A Doubling of the Post-Perovskite Phase Boundary and Structure of the Earth’s Lowermost Mantle. Nature, 2005, 434(7035): 882-886.

[11] Hunt S. A., Weidner D. J., Li L., . Weakening of Calcium Iridate during Its Transformation from Perovskite to Post-Perovskite. Nature Geoscience, 2009, 2(11): 794-797.

[12] Jacobs M. H. G., de Jong B. H. W. S.. Placing Constraints on Phase Equilibria and Thermophysical Properties in the System MgO-SiO₂ by a Thermodynamically Consistent Vibrational Method. Geochimica et Cosmochimica Acta, 2007, 71(14): 3630-3655.

[13] Jacobs, M. H. G., van den Berg, A. P., 2011. Complex Phase Distribution and Seismic Velocity Structure of the Transition Zone: Convection Model Predictions for a Magnesium-Endmember Olivine-Pyroxene Mantle. Phys. Earth Planet. Inter. (Accepted)

[14] Jarvis G. T., McKenzie D. P.. Convection in a Compressible Fluid with Infinite Prandtl Number. J. Fluid Mech., 1980, 96: 515-583.

[15] Karato S. I.. Deformation of Earth Materials: An Introduction to the Rheology of Solid Earth, 2008, Cambridge: Cambridge University Press

[16] Karato S. I.. The Influence of Anisotropic Diffusion on the High-Temperature Creep of a Polycrystalline Aggregate. Phys. Earth Planet. Inter., 2010, 183(3–4): 468-472.

[17] King S. D., Lee C., van Keken P. E., . A Community Benchmark for 2-D Cartesian Compressible Convection in the Earth’s Mantle. Geophys. J. Int., 2010, 180(1): 73-87.

[18] Lay T., Hernlund J., Garnero E. J., . A Post-Perovskite Lens and D″ Heat Flux beneath the Central Pacific. Science, 2006, 314(5803): 1272-1276.

[19] Murakami M., Hirose K., Kawamura K., . Post-Perovskite Phase Transition in MgSiO₃. Science, 2004, 304(5672): 855-858.

[20] Nakagawa T., Tackley P. J.. Influence of Initial CMB Temperature and Other Parameters on the Thermal Evolution of Earth’s Core Resulting from Thermochemical Spherical Mantle Convection. Geochemistry, Geophysics, Geosystems, 2010, 11 Q06001

[21] Schubert G., Turcotte D. L., Olson P.. Mantle Convection in the Earth and Planets, 2001, Cambridge: Cambridge University Press

[22] Steinbach V., Hansen U., Ebel A.. Compressible Convection in the Earth’s Mantle: A Comparison of Different Approaches. Geophys. Res. Lett., 1989, 16: 633-636.

[23] Tosi N., Cadek O., Martinec Z., . Is the Long-Wavelength Geoid Sensitive to the Presence of Post-Perovskite above the Core-Mantle Boundary?. Geophys. Res. Lett., 2009, 36 L05303

[24] Vacher P., Spakman W., Wortel M. J. R.. Numerical Tests on the Seismic Visibility of Metastable Minerals in Subduction Zones. Earth Planet. Sci. Lett., 1999, 170(3): 335-349.

[25] van den Berg A. P., van Keken P. E., Yuen D. A.. The Effects of a Composite Non-Newtonian and Newtonian Rheology on Mantle Convection. Geophysical Journal International, 1993, 115(1): 62-78.

[26] van den Berg A. P., de Hoop M. V., Yuen D. A., . Geodynamical Modeling and Multiscale Seismic Expression of Thermo-chemical Heterogeneity and Phase Transitions in the Lowermost Mantle. Phys. Earth Planet. Int., 2010, 180(3–4): 244-257.

[27] van der Hilst R. D., de Hoop M. V., Wang P., . Seismostratigraphy and Thermal Structure of Earth’s Core-Mantle Boundary Region. Science, 2007, 315(5820): 1813-1817.

[28] van Hunen J., van den Berg A. P., Vlaar N. J.. On the Role of Subducting Oceanic Plateaus in the Development of Shallow Flat Subduction. Tectonophysics, 2002, 352(3–4): 317-33.