Density and Viscosity Variations Due to Plume Melting of a Bilithologic Mantle: Implications for Asthenosphere and Hotspot Swell Root Dynamics

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