Dynamics of mantle flow and melting at a ridge-centered hotspot: Iceland and the Mid-Atlantic Ridge

The dynamics of mantle flow and melting of a ridge-centered plume were investigated with three-dimensional variable-viscosity numerical models, focusing on three buoyancy sources: temperature, melt depletion, and melt retention. The width, W, to which a plume spreads along a ridge axis, depends on plume volume flux,Q, full spreading rate, U, buoyancy number, B, and ambient/plume viscosity contrast γ. When all melting effects are considered, our numerical results are best parameterized by W = 2.37(Q/U) 1/2>(By) ^0.04. Thermal buoyancy is first-order in controlling along-axis plume spreading while latent heat loss due to melting, and depletion and retention buoyancy forces contribute second-order effects. We propose two end-member models for the Iceland plume beneath the Mid-Atlantic Ridge (MAR). The first has a broad plume source with temperature anomaly ΔT _p of 75°C, radius, a, of 300 km, andQof 1.2 × 10 ⁷ km ³/my. The second is of a narrower and hotter plume source with ΔT _p of 170°C, a radius of 60 km, andQof 2.1 × 10 ⁶ km ³/my. The broad plume source predicts successfully the observed seismic crustal thickness, topographic, and gravity anomalies along the MAR, but predicts an along-axis geochemical plume width substantially broader than that suggested by the observed ⁸⁷Sr/ ⁸⁶Sr anomaly. The narrow plume source model predicts successfully the total excess crustal production rate along the MAR (2.5 × 10 ⁵km ³/my) and a geochemical width consistent with that of the ⁸⁷Sr/ ⁸⁶Sr anomaly, but it requires substantial along-axis melt transport to explain the observed along-axis variations in crustal thickness, bathymetry, and gravity. Calculations suggest that lateral plume dispersion may be radially symmetric rather than channelled along the ridge axis and that the topographic swell, which is elongated along the Reykjanes Ridge, may be due to rapid off-axis subsidence associated with lithospheric cooling superimposed on a broader hotspot swell. The two plume source models predict seismic P-wave velocity reductions of 0.5 - 2% in the center of the plume, producing travel time delays of 0.2 - 1.2 s. Predicted P-wave delay times for the narrow plume source model are more consistent with recent seismic observations beneath Iceland, suggesting that this model may be more representative of the Iceland plume.