Melting systematics in mid-ocean ridge basalts: Application of a plagioclase-spinel melting model to global variations in major element chemistry and crustal thickness

Abstract

We present a new model for anhydrous melting in the spinel and plagioclase stability fields that provides enhanced predictive capabilities for the major element compositional variability found in mid-ocean ridge basalts (MORBs). The model is built on the formulation of Kinzler and Grove (1992) and Kinzler (1997) but incorporates new experimental data collected since these calibrations. The melting model is coupled to geodynamic simulations of mantle flow and mid-ocean ridge temperature structure to investigate global variations in MORB chemistry and crustal thickness as a function of mantle potential temperature, spreading rate, mantle composition, and the pattern(s) of melt migration. While the initiation of melting is controlled by mantle temperature, the cessation of melting is primarily determined by spreading rate, which controls the thickness of the lithospheric lid, and not by the exhaustion of clinopyroxene. Spreading rate has the greatest influence on MORB compositions at slow to ultraslow spreading rates (<2 cm/yr half rate), where the thermal boundary layer becomes thicker than the oceanic crust. A key aspect of our approach is that we incorporate evidence from both MORB major element compositions and seismically determined crustal thicknesses to constrain global variations in mantle melting parameters. Specifically, we show that to explain the global data set of crustal thickness, Na[subscript 8], Fe[subscript 8], Si[subscript 8], Ca[subscript 8]/Al[subscript 8], and K[subscript 8]/Ti[subscript 8] (oxides normalized to 8 wt % MgO) require a relatively narrow zone over which melts are pooled to the ridge axis. In all cases, our preferred model involves melt transport to the ridge axis over relatively short horizontal length scales (~25 km). This implies that although melting occurs over a wide region beneath the ridge axis, up to 20–40% of the total melt volume is not extracted and will eventually refreeze and refertilize the lithosphere. We find that the temperature range required to explain the global geochemical and geophysical data sets is 1300°C to 1450°C. Finally, a small subset of the global data is best modeled as melts of a depleted mantle source composition (e.g., depleted MORB mantle—2% melt).