A model for evolving crust on 4 Vesta through combined compositional and thermal modelling

Combined phase equilibrium and thermal modelling has been used to investigate the evolution of asteroid 4 Vesta. Orthopyroxene compositions of 200 natural diogenite meteorites are used as a basis for constructing a staged mantle melting model for Vesta, which is then used to develop a staged thermal evolution model. Our pMELTS models find that removal of 15–20% of a mean eucrite component from an initial Vestan mantle composition allows a second stage of melting that crystallises low-calcium orthopyroxenes that match the observed compositions of those in natural diogenites, whereas single stage melting produces orthopyroxenes that are too calcic. Using the compositions generated by the pMELTS modelling, THERMOCALC models were created for an initial Vestan mantle composition and an evolved composition generated by a melt extraction stage. These models suggest that melt production for second-stage diogenite generation required considerably hotter temperatures (>1340 °C) than for eucrites (<1240 °C). Staged and layered thermal evolution models developed using these composition and temperature constraints, based on the decay of ²⁶Al and ⁶⁰Fe, suggest that Vesta accreted 1.50 to 1.75 Myr after calcium-aluminium inclusion (CAI) formation. Earlier accretion results in conditions that are inconsistent with the petrology of the HED meteorites, whereas later accretion predicts temperatures that are insufficient to produce diogenites. We suggest that upward migration of ²⁶Al-rich melt initially created a convecting shallow magma ocean of <20 km depth that rapidly crystallised to form a ²⁶Al-rich eucritic crust that acted as a hot insulating lid. The second stage of crust formation began once the depleted mantle residue reached high enough temperatures to produce diogenite-forming magmas. These results further support the view that diogenites likely formed as crustal intrusions rather than as magma ocean cumulates.