Crustal accretion and dynamic feedback on mantle melting of a ridge centred plume: The Iceland case

A major consequence of the interaction of a plume with an oceanic ridge is the enhanced melt production and associated crust generation. In the case of Iceland crustal thickness as large as 20 to 40 km has been reported. Crustal seismic velocities are high, and have to be explained by thermal or chemical effects. In the first part of the paper we address the question whether extraction of melt out of the plume beneath a slowly spreading ridge and deposition of extracted basalt volumes at the surface produces a dynamic feedback mechanism on mantle melting. To study this question we solve the convection equations for a ridge centred plume with non-Newtonian rheology including melting, melt extraction associated with deposition of cold crust at the surface of the model, and using a simplified approach for compaction. The assumption of cold crust is justified if the thickness of each deposited basaltic layer is less than roughly 1 km. Depending on the buoyancy flux of the plume, crustal thicknesses between 10 and 40 km are modelled, showing characteristic dipping structures resembling the rift-ward dipping basaltic layers of East- and Western Iceland. Comparing the resulting crustal thickness and magma generation rate with models in which the dynamic effect of crust deposition has been suppressed indicates, that melt generation beneath a slowly spreading ridge is considerably damped by the dynamic feedback mechanism if the plume buoyancy flux exceeds 400 to 600 kg/s. Based on the observed crustal thickness of Iceland our models predict a plume buoyancy flux of 1140 kg/s.In the second part we study the accretion of the Icelandic crust by a thermo-mechanical model in more detail based on the Navier–Stokes-, the heat transport and the mass conservation equations including volumetric sources. Hot (1200 °C) molten crustal material is injected into the newly forming crust with a constant rate at different crustal source regions: a) deep, widespread emplacement of dykes and sills including crustal underplating, b) magma chambers at shallow to mid-crustal level, and c) surface extrusions and intrusions in fissure swarms at shallow depth connected to volcanic centres. We identify the material from the different source regions by a marker approach. Varying the relative dominance of these source regions, characteristic crustal structures evolve, showing shallow dipping upper crustal layers with dip angles between 10 and 15°. The thermal structure of the crust varies between cold crust (shallow-source region dominating) and hot crust (deep-source region dominating). We use observations of maximum depth of seismicity to constrain the depth of the 650 °C isotherm and seismological inferences on the lower crust to constrain temperatures in that region. The best agreement with our models is achieved for crust formation dominated by deep dykes and underplating with a considerable influence of magma chamber accretion.