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.

Prior to intrusion, magma migrates through the crustal plumbing system that likely contains layers or columns of crystal mush. To better understand the behavior of the crustal magmatic system during magmatic unrest, it is important to examine the process of melt migration within the crystal mush and the associated evolution in pressure and temperature. In this study I use an analytical model to explore the characteristics of transport of melt, pressure, and heat through an idealized crystal mush layer/column under uniaxial strain condition. The model invokes a thermo-poro-viscoelastic rheology and uses a frequency-domain method to explore two scenarios of magmatic unrest: harmonic perturbation of fluid pressure, and step-rise in fluid pressure at a source location. Several factors influence the transport of melt, pressure and heat, including the thermal-mechanical coupling arising from the mush rheology, the advection of heat by melt flows, the competition between thermal diffusivity and poroelastic diffusivity, and the viscoelastic relaxation of the crystalline framework. One key finding is the development of transport asymmetry: when a background temperature gradient exists, the transport properties become different for propagation along the background thermal gradient and propagation against the background thermal gradient. Analysis on an endmember case shows that the transport asymmetry is associated to the competition between the diffusion and advection of pore pressure, which determines a Peclet number that depends on the temperature difference across the mush and the thermal expansion coefficients. Because the temperature in magma mushes in the crust likely increase with depth, the observed propagation asymmetry suggests some intrinsic difference between a bottom-up vs. a top-down triggering mechanism for magmatic unrest. The results from this study highlight the importance for further exploration for a more complete description of the transport properties in the crystal mush.

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.

Core formation has left a lasting geochemical signature on the Earth. In order to constrain the composition of the Earth we must fully understand the processes by which newly formed Earth, and the bodies which accreted to it, differentiated. Percolation of iron-rich melt through solid silicate has been invoked as a mechanism for differentiation and core formation in terrestrial bodies in the early solar system. However, to date the contribution of percolation to core formation cannot be assessed due to the absence of data on Fe-rich melt migration velocities. Here we use a novel experimental design to investigate textural changes in an analog system, Au melt in polycrystalline h-BN, at 3 GPa, relevant to core formation in the early solar system. Using a combination of high resolution, in-situ X-ray tomography and fast 2-D radiographic imaging, we obtain the first direct data on melt migration velocities at high PT. Melt migration is highly variable and episodic, driven by variations in differential pressure during melt migration and matrix compaction. Smaller scale melt processes, representing migration of melt along pre-existing melt networks, give comparatively fast velocities of 0.6–60 μms−1. Ex-situ experiments are used to compare melt networks in analog systems to Fe-rich melt in silicates. Two competing processes for melt migration are percolation of melt along grain boundaries, and hydraulic fracturing induced by melt injection. Typically, both processes are noted in experimental and natural systems, although the relative importance of each mechanism is variable. Using a simple model for melt flow through a porous media, migration velocities determined here account for full differentiation of Earth-sized bodies within 101–103 Myr, for submicron diameter melt bands, or within a few Myr or micron-sized melt bands. This is consistent with rapid timescales inferred from geochemistry for core formation in planetesimals, implying that percolation may have had an important contribution to core differentiation in the Earth.