Carbon dioxide in igneous petrogenesis: II. Fluid dynamics of mantle metasomatism

Petrographic, fluid inclusion, geochemical and isotopic evidence from xenoliths in alkali basalts suggests that low-viscosity fluids rich in O-H-C, dissolved silicates and especially the incompatible elements may ascend, decompress and precipitate crystalline phases and/or induce partial fusion in the upper mantle. Such mantle metasomatic fluids (MMF) may be important in generating isotopic heterogeneity and in transporting and focusing mantle heat. In order to model the movement of MMF, the ordinary differential equations governing the variation ofP, T, ascent velocity and fluid density of a compressible, viscous, single-phase (H₂O or CO₂) non-reacting fluid ascending through a vertical crack of constant width have been solved. A large number of numerical simulations were carried out in which the significant factors affecting flow behavior (thermodynamic and transport fluid properties, roughness and width of cracks, geothermal gradient, initial conditions, etc.) were systematically varied. The calculations show that: (1) MMF tends to move at uniform rates following a short period of rapid initial acceleration, (2) MMF ascends nearly isothermally, (3) MMF acts as an efficient heat transfer agent; numerical experiments show that transport of heat into regions undergoing metasomatism can lead to partial fusion. The heat transported by movement of MMF averaged over the age of the Earth is sufficient to generate about 0.1 km³ of basaltic magma per year, which is approximately equal to the production rate of alkaline magma. If an intense period of mantle degassing occured early in the history of the Earth, the transport of heat and mass (K, U, Rb, LREE) by migrating fluids might have been important.