DOI: 10.5281/zenodo.17406367
title: Density and Viscosity Variations Due to Plume Melting of a Bilithologic Mantle: Implications for Asthenosphere and Hotspot Swell Root Dynamics

description: This dataset comprises the output data from Perple_X thermodynamic computations and numerical simulations, accompanied by the MATLAB scripts necessary to reproduce all figures presented in this paper:

Shao and Morgan Density and viscosity variations due to plume melting of a bilithologic mantle: Implications for asthenosphere and hotspot swell root dynamics, under review in Earth and Planetary Science Letters.

% ----------------------------------------------------------------

1. DataCodes_Thermodynamic (include the data from perple_X and m-files)

% ----------------------------------------------------------------

% m-file

Fig4_PEffectDensity: produce figure 4

 

% file folds contain original Perple_X output files

Pe_DP10% : the original Perple_X output files (density and mineralogy data) for a 10%-depleted peridotite, calculated under increasing temperature at constant pressures of 1-3 GPa.

Pe_KLB1 : the original Perple_X output files (density and mineralogy data) for a KLB-1 peridotite, calculated under increasing temperature at constant pressures of 1-3 GPa.

Py_G2     : the original Perple_X output files (density and mineralogy data) for a G2 pyroxenite, calculated under increasing temperature at constant pressures of 3-7 GPa.

Py_MIX1G : the original Perple_X output files (density and mineralogy data) for a MIX1G pyroxenite, calculated under increasing temperature at constant pressures of 3-7 GPa.

 

% file folds contain processed data

Results_F_rhoF_percent : This fold contains the processed data (interim and final results) alongside the MATLAB scripts used for analysis.

                   F_DeltaDensFPercent_per.m : produce Figure S5

                   F_DeltaDensFPercent_pyx.m : produce Figure S6,S7

 

ThermalExpansivity_P : Perple_X output thermal expansivity codes and files, calculated under increasing pressures and at constant temperature (1573K)

                   Figure_S1_ThermalExpan_P.m : produce Figure S1

 



% ----------------------------------------------------------------

2. Codes_NumericalSimulation (a m-file and a database)

% ----------------------------------------------------------------

WMCM_V3_REACT_Mix : The primary numerical code developed and employed in this study.

database2                        : The fitted relationship between lithologic density variations and melt extraction (a necessary input for the numerical codes).

 


% ----------------------------------------------------------------

3. Data_NumericalSimulation (MATLAB files, *.mat)

% ----------------------------------------------------------------

DT_X0_MIX1G        : Pure 10%-depleted peridotite with an initial temperatures of 1300-1450 oC (Figure 6,S8).

DT_X8218_MIX1G : 82% 10%-depleted peridotite + 18% MIX1G pyroxenite with an initial temperatures of 1300-1450 oC (Figure 6,S8).

G2                           : Bilitholigic mantle containing 1-30% G2 pyroxenite under the conditions of initial temperatures of 1300-1700 oC (Figure 2,3,5b,5c,S4).

MIX1G                     : Bilitholigic mantle containing 1-30% MIX1G pyroxenite under the conditions of initial temperatures of 1300-1700 oC (Figure 5b,c).

T1500_DVis           : 80% 10%-depleted peridotite + 20% pyroxenite with an initial temperatures of 1500 oC, test the viscosity prefactor (Figure 7).

T1500_P                 : 80% 10%-depleted peridotite + 20% pyroxenite with an initial temperatures of 1500 oC, test the effect of pressure from therodynamic results (Figure 5a)

T1500_Xsdry400. : 80% 10%-depleted peridotite + 20% pyroxenite with an initial temperatures of 1500 oC, peridotite is dry and pyroxenite contains 400 ppm water (Figure S3)

 

% ----------------------------------------------------------------

4. PostCodesNumericalSimulation (MATLAB scipts, *.m)

% ----------------------------------------------------------------

Fig2_DensVis1DPlot                    : produce figure 2

Fig3_Comp_DensVis1DPlot_v2 : produce figure 3

Fig5_PlotMain2                           : produce figure 5

Fig6a_S8                                     : produce figure 6a and S8

Fig6b                                           : produce figure 6b

Fig7_Comp_Vis                           : produce figure 7

FigS2_PerDepl                           : produce figure S2

FigS3_dryper                              : produce figure S3

Figure_S9_Homogenization_KLB1 : produce figure S9