Periodic Reporting for period 1 - MoVEMENT (Mobility of Volatiles in the Earth’s Mantle by Experimental and Numerical Technics)
Período documentado: 2020-02-01 hasta 2022-01-31
Volatiles cycle has a leading place in the evolution and fate of a planet, such as typically by defining the habitability conditions at the surface. On Earth, CO2 and H2O are the two most abundant volatiles in the upper mantle, playing a critical role on melting of mantle rocks, and triggering the production of magmatic liquids whose composition ranges from carbonatites (i.e. low silica content) deep in the mantle to basanites-basalts (i.e. high silica content) closer to the surface. These melts are the major hosts of carbon (and hydrogen) and their mobility through the mantle controls the distribution of carbon in the mantle, and at a larger scale the global deep carbon cycle with fluxes between surficial and deep reservoirs.
The aim of the project MoVEMENT is to combine fundamental constraints on the physical properties, namely density and viscosity, of CO2-H2O-bearing melts with complex modelling to gain a quantitative understanding of the mobility and geophysical signature of these important phases during mantle magmatic processes. Better constraints on these processes will lead in turn to a new understanding of volatiles circulation and recycling in the deep Earth, and their impact in surficial processes.
Especially, we identified three main research questions to be answered by the MoVEMENT project:
1. What is the density and viscosity of carbonatite-kimberlite-basalt melts in the upper mantle (WP1)? What are the effects of pressure, temperature, composition (e.g. change in SiO2 and H2O content) on these properties?
2. What is the most appropriate parameterization to describe the physical properties (density, viscosity) as a function of compositional variations in mantle melts with P-T (WP2)?
3. What is the mobility of volatile-bearing melts in the upper mantle (WP3)? What is the geophysical signature of melts in the upper mantle and their role in mantle geodynamics?
To answer these questions, we have provided new data on the viscosity of carbonate-rich melts by classical Molecular Dynamic (MD) simulations and further designed, implemented and applied thermodynamics-based models predicting density and viscosity of mantle melts covering a broad range of relevant pressures, temperatures and compositions. Further, we applied these models to predict geophysical signatures of such melts within the Earth’s mantle.
The aim of the project MoVEMENT is to combine fundamental constraints on the physical properties, namely density and viscosity, of CO2-H2O-bearing melts with complex modelling to gain a quantitative understanding of the mobility and geophysical signature of these important phases during mantle magmatic processes. Better constraints on these processes will lead in turn to a new understanding of volatiles circulation and recycling in the deep Earth, and their impact in surficial processes.
Especially, we identified three main research questions to be answered by the MoVEMENT project:
1. What is the density and viscosity of carbonatite-kimberlite-basalt melts in the upper mantle (WP1)? What are the effects of pressure, temperature, composition (e.g. change in SiO2 and H2O content) on these properties?
2. What is the most appropriate parameterization to describe the physical properties (density, viscosity) as a function of compositional variations in mantle melts with P-T (WP2)?
3. What is the mobility of volatile-bearing melts in the upper mantle (WP3)? What is the geophysical signature of melts in the upper mantle and their role in mantle geodynamics?
To answer these questions, we have provided new data on the viscosity of carbonate-rich melts by classical Molecular Dynamic (MD) simulations and further designed, implemented and applied thermodynamics-based models predicting density and viscosity of mantle melts covering a broad range of relevant pressures, temperatures and compositions. Further, we applied these models to predict geophysical signatures of such melts within the Earth’s mantle.
The first step of the project was to develop a continuous model for the density of carbonate melts in the system MgO-CaO-Na2O-K2O-Li2O-H2O-CO2. Such a model was missing, in spite of its unique potential to investigate a diversity of carbon-mediated magmatic processes, from carbonatitic eruptions to carbonate melts migration and crystal fractionation, and to constrain the seismic signature of carbonate melts in the upper mantle/transition zone. The present density model is calibrated to span most of these conditions of interest. After having tested different equations of state (EoS), the pressure-volume relations of each oxide component are ultimately described by employing a simple Murnaghan EOS based on 5 adjustable parameters, and describing the partial molar volume and compressibility of the different oxide components. The calibrated density model reproduces well the calibration dataset and predicts densities with an average relative error lower than ±5%. Combined with reasonable assumptions for the volumetric properties of SiO2 and FeO, the applications of the model to density and seismic P-wave velocities Vp predictions in complex carbonatite systems issued from partial melting carbonated lithologies provide new insights into their migration and geophysical signature. Some of the mean observations issued from our calculations are:
1- Both density and VP contrasts between carbonatitic melts and the host solid mantle diminish with depth, which implies larger effect of carbonate-rich melts on the seismic velocity signature of the Earth’s upper mantle at shallower depths;
2- We pinpoint that carbonate melts are seismically faster than silicate melts, which together with their higher electrical conductivity, are promising diagnostic features to identify molten phases at depth and to discriminate between carbonate and silicate melts in the Earth’s upper mantle using geophysical observations.
In a second step, we developed a continuous model for the viscosity of mantle melts, spanning the whole compositional spectrum of melts prevailing within the Earth’s upper mantle. The model is calibrated to be used in most conditions of interest for the crust and upper mantle. We have developed a new non-Arrhenian three-parameters model to predict the contrasted viscosity properties of mantle melts, from carbonate-rich to silicate-rich compositions, inspired of a reformulation of the MYEGA model. The model is parameterized in a way to include effects of pressure-temperature-composition, and is still under construction while most of the calibration and related results have been obtained so far, and a publication is now under preparation.
Viscosity of mantle melts with various compositions, from carbonatites rich in CO2 to basalts enriched in SiO2, will be calculated and presented in the publication, along with estimations for the mobility of melts within the host rocks.
The database used for the development of the continuous viscosity model has benefited from the viscosity points of mixed alkali and alkaline-earth melts at upper mantle conditions, i.e. up to 13 GPa and 1000 – 2000 K, obtained by MD simulations. The latter were performed by a former PhD student in the host group, Dr. Xenia Ritter, in collaboration with Profs. B. Guillot and N. Sator from Sorbonne University (Paris, France). I have participated to this study which provides first evidence for non-Arrhenian temperature dependent viscosity of carbonate melts at upper mantle pressures, which may be a behaviour extensible to other carbonate melt compositions. To facilitate the integration of the new data into carbonate mobility and geodynamic modelling, I employed the data to parameterize a model describing the pressure-temperature dependence of the viscosity of alkali(ne) carbonate melts at upper mantle conditions.
The density model and applications are reported in a publication currently in revision (moderate revisions) in ‘Chemical Geology’, where I am the first author, as well as at different conferences: GeoMünster – Annual Conference of the German Geological Society and the German Mineralogical Society, Münster, Germany; EMPG 2021, Potsdam, Germany (online); Goldschmidt 2021, Lyon, France (online); and the Japan Geoscience Union 2022. The viscosity model and applications will be submitted to ‘Earth and Planetary Science Letters’.
1- Both density and VP contrasts between carbonatitic melts and the host solid mantle diminish with depth, which implies larger effect of carbonate-rich melts on the seismic velocity signature of the Earth’s upper mantle at shallower depths;
2- We pinpoint that carbonate melts are seismically faster than silicate melts, which together with their higher electrical conductivity, are promising diagnostic features to identify molten phases at depth and to discriminate between carbonate and silicate melts in the Earth’s upper mantle using geophysical observations.
In a second step, we developed a continuous model for the viscosity of mantle melts, spanning the whole compositional spectrum of melts prevailing within the Earth’s upper mantle. The model is calibrated to be used in most conditions of interest for the crust and upper mantle. We have developed a new non-Arrhenian three-parameters model to predict the contrasted viscosity properties of mantle melts, from carbonate-rich to silicate-rich compositions, inspired of a reformulation of the MYEGA model. The model is parameterized in a way to include effects of pressure-temperature-composition, and is still under construction while most of the calibration and related results have been obtained so far, and a publication is now under preparation.
Viscosity of mantle melts with various compositions, from carbonatites rich in CO2 to basalts enriched in SiO2, will be calculated and presented in the publication, along with estimations for the mobility of melts within the host rocks.
The database used for the development of the continuous viscosity model has benefited from the viscosity points of mixed alkali and alkaline-earth melts at upper mantle conditions, i.e. up to 13 GPa and 1000 – 2000 K, obtained by MD simulations. The latter were performed by a former PhD student in the host group, Dr. Xenia Ritter, in collaboration with Profs. B. Guillot and N. Sator from Sorbonne University (Paris, France). I have participated to this study which provides first evidence for non-Arrhenian temperature dependent viscosity of carbonate melts at upper mantle pressures, which may be a behaviour extensible to other carbonate melt compositions. To facilitate the integration of the new data into carbonate mobility and geodynamic modelling, I employed the data to parameterize a model describing the pressure-temperature dependence of the viscosity of alkali(ne) carbonate melts at upper mantle conditions.
The density model and applications are reported in a publication currently in revision (moderate revisions) in ‘Chemical Geology’, where I am the first author, as well as at different conferences: GeoMünster – Annual Conference of the German Geological Society and the German Mineralogical Society, Münster, Germany; EMPG 2021, Potsdam, Germany (online); Goldschmidt 2021, Lyon, France (online); and the Japan Geoscience Union 2022. The viscosity model and applications will be submitted to ‘Earth and Planetary Science Letters’.
Models and associated results/interpretations/messages are, or will be, made available to the scientific community through different publication in peer-reviewed journals listed above. A main objective is to further supplement these models into a more global framework where a comprehensive description of the physicochemical properties of melts and their host rocks will be calculated. Such outputs will bring quantitative constraints on different fundamental mechanisms ruling the Earth interior, e.g. the mobility of the melts and their impact on geophysical (seismicity, electrical conductivity), geodynamic (Earth’s plate tectonic engine) and chemical/biological (volatiles cycle) properties.