190 likes | 325 Views
Chemical consequences of perovskite fractionation from an ultramafic liquid with application to the dynamics of a basal magma ocean . Leah Ziegler, Hongluo Zhang, Colin Jackson, Matt Jackson, Dave Stegman. Team BMO:. What is a basal magma ocean?.
E N D
Chemical consequences of perovskite fractionation from an ultramafic liquid with application to the dynamics of a basal magma ocean Leah Ziegler, Hongluo Zhang, Colin Jackson, Matt Jackson, Dave Stegman Team BMO:
What is a basal magma ocean? - a dense melt layer that forms at the bottom of the mantle – likely linked to a whole mantle magma ocean Labrosse et al., 2007
Why a basal magma ocean? -poorly understood lower mantle structures: Garnero and McNamara, 2008 -missing geochemical reservoirs: Labrosse et al., 2007
How would a basal magma ocean form? w/ Fe w/o Fe -Curvature of mantle liquidus can lead to crystallization initiating at mid-mantle depths -Density crossover between liquid and solid at mid-mantle depths is possible – depends on Kd-Fe Figures: Stixrude et al., 2009
Key parameters for evaluating BMO: Density of liquid and solid with crystallization • f(Xliq, Xsol, XU+Th+K, time) Concentration of REE and heat producers (U+Th) associated with BMO products • Both Require: - major element partition coefficients (Kd) - trace element partition coefficients
Step 1: Compile Kd database for Mg-Pv, Ca-Pv, & Fe-MgO -Database compiled at CIDER 2012 -Details of database: Total number of studies: 19 Experiments including -- Mg-Pv: 31 -- Ca-Pv: 10 – Fe-MgO:12 Pressure Range: 23 – 86 GPa Temperature Range: 2300 – 3500°C Database will be published with initial reporting of findings
Effect of Al content on the solubility of Fe in MgPv & Kd-TE right: Frost and Langenhorst, 2002 left: Liebske et al., 2005 Importance of charge coupled substitutions
Step 2: Parameterize Kds for Mg-Pv Models for Mg-Pv major elements: Predicted Observed Equal weight, multiple linear regression
Step 2: Parameterize Kds for Mg-Pv Models for Mg-Pv trace elements: Predicted Observed Similar TE substitutions as pyroxene
Step 3: Apply Kds for Mg-PvBMO crystallization Team BMO met in October at BU: - Refined Kd parameterizations • Developed stoichiometric Mg-Pv crystallization model • Quantified uncertainties in crystallization model predictions using montecarlo approach • Accounts for covariation in model parameters • Thank you Don Forsyth
Model Results Liquid evolutionwith crystallization (batch): Bulk: McDonough & Sun BSE, 40 GPa Majors: cation mol fraction, Trace: PPM*
Model Results KdFewith crystallization (BSE bulk, 40 GPa): Al-bearing Al-free Results from 78.5 GPa: Andrault et al., 2012 Al + Si causing offset between Nomura and Andrault ?
Model Results Removal of heat producing elements and REE from liquid (batch): High Al in MgPv makes it significant reservoir for TEs Still working to incorporate full uncertainties into TE Kds
Preliminary Findings -Al partitions equally between MgPV and liquid -REE+U+Th have similar substitutions into pyx and MgPv - Al AND Ca are important -Fe is a moderately incompatible element • Fe is strongly controlled by Al, but other effects appear significant (Siliq?)
Current Directions Post-AGU Team BMO Meet up: - Develop new crystallization model • Weaver and Langmuir (1990) • Working with StephaneEscrig to implement his code - stoichiometry - self-consistent saturation of phases - Predict past Fe-MgO or Ca-Pv saturation - Couple chemical model to dynamical model
Crystallization Model System of Equations: -Maintains stoichiometry on Si and Mg site
Effect of Al content on Kd-Fe at high P & T Al-bearing Al-free Andraultet al., 2012