Methods for resolving the origin of large igneous provinces from crustal seismology

Abstract

We present a new quantitative framework to understand the process of mantle melting based on the velocity structure of igneous crust. Our approach focuses on the lower crustal section, which is expected to be least affected by porosity and seawater alteration, especially for thick igneous crust. Our methodology is thus best for constraining the origin of large igneous provinces. First, a quantitative relation between bulk crustal velocity and mantle melting parameters is established on the basis of data from mantle melting experiments. Second, we show how lower crustal velocity can be used to place bounds on the expected range of bulk crustal velocity, despite ambiguity in crustal emplacement processes. By modeling fractional crystallization processes at a range of crustal pressures, these bounds are derived as a function of the proportion of lower versus upper crust. Finally, a simple mantle melting model is constructed to illustrate the effects of potential temperature, active upwelling, and a preexisting lithospheric lid on predicted crustal thickness and velocity. As an example, this new interpretation is applied to a seismic transect across the southeast Greenland margin to constrain mantle dynamics during the opening of the North Atlantic. Some complicating factors in seismological inference on mantle melting process, such as the possibilities of subcrustal igneous fractionation and mantle source heterogeneity, are also discussed with this example.