Influence of magmatism on mantle cooling, surface heat flow and Urey ratio

Two-dimensional thermo-chemical mantle convection simulations are used to investigate the influence of melting-inducted differentiation on the thermal evolution of Earth's mantle, focussing in particular on matching the present-day surface heat flow and the ‘Urey ratio’. The influence of internal heating rate, initial mantle temperature and partitioning of heat-producing elements into basaltic crust are studied. High initial mantle temperatures, which are expected following Earth's accretion, cause major differences in early mantle thermo-chemical structures, but by the present-day surface heat flux and internal structures are indistinguishable from cases with a low initial temperature. Assuming three different values of mantle heat production that vary by more than a factor of two results in small differences in present-day heat flow, as does assuming different partitioning ratios of heat-producing elements into crust. Indeed, all of the cases presented here, regardless of exact parameters, have approximately Earth's present-day heat flow, with substantial fractions coming from the core and from mantle cooling. As a consequence of the model present-day surface heat flow varying only slightly with parameters, the Urey ratio (the ratio of total heat production to the total surface heat flow) is highly dependent on the amount of internal heat production, and due to the large uncertainty in this, the Urey ratio is considered to be a much poorer constraint on thermal evolution than the heat flow. The range of present-day Urey ratio observed in simulations here is about 0.3 to 0.5, which is consistent with observational and geochemical constraints (Jaupart et al., 2007). Magmatic heat transport contributes an upper bound of 9% to Earth's present-day heat loss but a much higher fraction at earlier times—often more than convective heat loss—so neglecting this causes an overestimation of the Urey ratio. Magmatic heat transport also plays an important role in mantle cooling. Considering these points, it is important to include magmatic effects when attempting to understand the thermal evolution of the Earth.

► We investigate the plausibility of constraint of thermal evolution of Earth's mantle. ► The heat transfer by magmatism is important for understanding mantle cooling. ► The Urey ratio is not good constraint for understanding thermal evolution. ► The surface heat flow would be better than the Urey ratio. ► Current mantle structure cannot be constrained for the early Earth's state.