Partial melting in an upwelling mantle column

Decompression melting of hot upwelling rock in the mantle creates a region of partial melt comprising a porous solid matrix through which magma rises buoyantly. Magma transport and the compensating matrix deformation are commonly described by two-phase compaction models, but melt production is less often incorporated. Melting is driven by the necessity to maintain thermodynamic equilibrium between mineral grains in the partial melt; the position and amount of partial melting that occur are thus thermodynamically determined. We present a consistent model for the ascent of a one-dimensional column of rock and provide solutions that reveal where and how much partial melting occurs, the positions of the boundaries of the partial melt being determined by conserving energy across them. Thermodynamic equilibrium of the boundary between partial melt and the solid lithosphere requires a boundary condition on the effective pressure (solid pressure minus melt pressure), which suggests that large effective stresses, and hence fracture, are likely to occur near the base of the lithosphere. Matrix compaction, melt separation and temperature in the partially molten region are all dependent on the effective pressure, a fact that can lead to interesting oscillatory boundary-layer structures.