Mantle Convection Simulation

simulator intermediate ~10 min
Loading simulation...
3 convection cells — Surface velocity: 4.2 cm/yr. Heat flow: 72 mW/m². Nu = 12.4.

At Rayleigh number 10^7 with 2500 K temperature contrast and 60% internal heating, convection produces approximately 3 cell pairs with surface velocity 4.2 cm/yr and heat flow 72 mW/m2.

Formula

Ra = α g ΔT D³ / (κ ν)
Nu ~ 0.294 Ra^(1/3) (Nusselt-Rayleigh scaling)
v_rms ~ (κ/D) Ra^(2/3) (velocity scaling)

Thermal Convection in Earth's Mantle

Earth's mantle, despite being solid rock, flows like an extremely viscous fluid over geological time scales. Heat from the core and radioactive decay creates buoyancy forces that drive convection: hot, less-dense material rises while cold, dense material sinks. This process transports over 70% of Earth's internal heat to the surface and drives plate tectonics.

The Rayleigh Number

The Rayleigh number is the key dimensionless parameter governing convection. It quantifies the competition between thermal buoyancy (promoting flow) and viscous resistance plus thermal diffusion (inhibiting flow). When Ra exceeds a critical value, convection initiates as steady rolls. At higher Ra values, the flow becomes increasingly complex, time-dependent, and eventually chaotic.

Boundary Layers and Plumes

Convection in the mantle develops thin thermal boundary layers at the top (lithosphere) and bottom (D" layer above the core-mantle boundary). These boundary layers can become gravitationally unstable, spawning cold downwellings (subducting slabs) from the top and hot upwellings (mantle plumes) from the bottom. The asymmetry between these features reflects the role of internal heating from radioactive decay.

Heat Flow and Nusselt Number

The Nusselt number compares actual convective heat transport to what pure conduction would achieve. For mantle-like Rayleigh numbers, Nu ~ Ra^(1/3), meaning convection enhances heat loss by an order of magnitude. Earth's mean surface heat flow of ~87 mW/m² reflects this convective enhancement, with higher values at mid-ocean ridges and lower values in continental shield regions.

FAQ

What drives mantle convection?

Mantle convection is driven by thermal buoyancy. Heat from the core and radioactive decay in the mantle creates density differences: hot material rises and cold material sinks. The Rayleigh number, proportional to temperature contrast and inversely proportional to viscosity, determines whether convection occurs and how vigorous it is.

What is the Rayleigh number?

The Rayleigh number Ra = alpha * g * dT * D^3 / (kappa * nu) compares buoyancy forces to viscous and thermal diffusion. For Ra above the critical value (~657-1100 depending on boundary conditions), convection begins. Earth's mantle Ra is estimated at 10^6-10^8, indicating vigorous convection.

How fast does the mantle convect?

Mantle convection velocities are typically 1-10 cm/year, comparable to plate tectonic speeds. The mantle takes hundreds of millions of years to overturn once. Despite these slow speeds, the enormous scale (2900 km depth) and high viscosity make convection the dominant heat transport mechanism.

What is the relationship between mantle convection and plate tectonics?

Plate tectonics is the surface expression of mantle convection. Subducting slabs are cold downwelling limbs, mid-ocean ridges are passive upwelling responses to plate separation, and hotspot volcanoes may be fed by deep mantle plumes. The coupling between convection and plate motion remains an active research area.

Sources

Embed

<iframe src="https://homo-deus.com/lab/geophysics/mantle-convection/embed" width="100%" height="400" frameborder="0"></iframe>
View source on GitHub