Search

Abstract The widespread High Arctic Large Igneous Province (HALIP) exhibits prolonged melting over more than 50 Myr, an observation that is difficult to reconcile with the classic view that large igneous provinces (LIPs) originate from melting in plume heads. Hence, the suggested plume‐related origin and classification of HALIP as a LIP have been questioned. Here, we use numerical models that include melting and melt migration to investigate a rising plume interacting with lithosphere of variable thickness, that is, a basin‐to‐craton setting applicable to the Arctic. Models reveal that melt migration introduces significant spatial and temporal variations in melt volumes and pulses of melt production, including protracted melting for at least about 30–40 Myr, because of the dynamic feedback between migrating melt and local lithosphere thinning. For HALIP, plume material deflected from underneath the Greenland craton can re‐activate melting zones below the previously plume‐influenced Sverdrup Basin after a melt‐free period of about 10–15 Myr, even though the plume is already ∼500 km away. Hence, actively melting zones do not necessarily represent the location of the deeper plume stem at a given time, especially for secondary pulses. Additional processes such as (minor) plume flux variations or local lithospheric extension may alter the timing and volume of HALIP pulses, but are to first order not required to reproduce the long‐lived and multi‐pulse magmatism of HALIP. Since melting zones are always plume‐fed, we would expect HALIP magmatism to exhibit plume‐related trace element signatures throughout time, potentially shifting from mostly tholeiitic toward more alkalic compositions.

Abstract Melt transport across the ductile mantle is essential for oceanic crust formation or intraplate volcanism. However, mechanisms of melt migration and associated chemical interaction between melt and solid mantle remain unclear. Here, we present a thermo‐hydro‐mechanical‐chemical (THMC) model for melt migration coupled to chemical differentiation. We consider melt migration by porosity waves and a chemical system of forsterite‐fayalite‐silica. We solve the one‐dimensional (1D) THMC model numerically using the finite difference method. Variables, such as solid and melt densities or MgO and SiO2 mass concentrations, are functions of pressure (P), temperature (T), and total silica mass fraction (CTSiO2). These variables are pre‐computed with Gibbs energy minimization and their variations with evolving P, T, and CTSiO2 are implemented in the THMC model. We consider P and T conditions relevant around the lithosphere‐asthenosphere boundary. Systematic 1D simulations quantify the impact of initial distributions of porosity and CTSiO2 on the melt velocity. Larger perturbations of CTSiO2 cause larger melt velocities. An adiabatic or conductive geotherm cause fundamentally different vertical variations of densities and concentrations, and an adiabatic geotherm generates higher melt velocities. We quantify differences between melt transport (considering incompatible tracers), major element transport and porosity evolution. Melt transport is significant in the models. We also quantify the relative importance of four porosity variation mechanisms: (a) mechanical compaction and decompaction, (b) density variation, (c) compositional variation, and (d) solid‐melt mass exchange. In the models, (de)compaction dominates the porosity variation. We further discuss preliminary results of 2D THMC simulations showing blob‐like and channel‐like porosity waves.

Abstract To investigate channelization during migration of a reactive melt, we performed a series of Darcy‐type experiments in which an alkali basalt infiltrated partially molten harzburgites and lherzolites at a confining pressure of 300 MPa, temperatures of 1200°C and 1250°C, and pore pressure gradients of ~2 to 60 MPa/mm. We compare our results to those from previously published experiments performed on wehrlites. In all experiments, irrespective of the exact mineralogy, a planar reaction layer composed of olivine + melt developed in which all of the pyroxene was consumed. Under specific conditions controlled primarily by the melt flow velocity, finger‐like channels composed of olivine + melt also developed. In wehrlites, these reaction infiltration instabilities formed at 1200°C and 1250°C at pressure gradients >25 and >5 MPa/mm, respectively. In harzburgites, channelization occurred only at 1250°C at a pressure gradient of 35 MPa/mm. In lherzolites, a planar melt‐filled vein developed at 1250°C; no finger‐like channels formed under a pressure gradient of ~25 MPa/mm. Both the finger‐like channels and the planar vein led to very efficient extraction of melt from the reservoir. Channelization established large compositional variations over short distances in the crystalized phases as well as in the local melt and greatly enhanced the abundance of the reaction product, olivine, similar to dunite channels in the Earth. The range of chemical‐mechanical responses displayed by this array of compositions provides a set of targets for reactive transport and mechanical modeling studies.