Tectonic Plates & Mantle Convection

Structure of the Earth

Earth is not a uniform solid. Its interior is layered like an onion, each layer with distinct composition, temperature and mechanical behaviour:

The lithosphere — the crust plus the rigid uppermost mantle — sits atop the asthenosphere, a partially molten, mechanically weak zone in the upper mantle roughly 80–200 km deep. It is on this soft layer that the tectonic plates "float" and slide.

Mantle Convection

Although solid, the mantle behaves like an extremely viscous fluid on geological timescales (millions of years). Its viscosity is approximately 10²¹ Pa·s — about 10²⁴ times more viscous than water, yet slow thermal convection is still possible.

The mechanism is simple thermodynamics: hot material deep near the core is less dense and rises; cooler, denser material near the base of the lithosphere sinks. This sets up convection cells that drag the overlying plates with them.

The governing equations are the Stokes equations coupled with the heat transport equation. For incompressible flow in the mantle:

Where u is velocity, P pressure, η viscosity, ρ density, T temperature and κ thermal diffusivity. The first two equations say the mantle flows to balance pressure and viscous forces against gravity — inertia is negligible at these scales. The third equation says temperature is advected by the flow and diffused by conduction.

The dimensionless Rayleigh number determines whether convection occurs:

Where α is thermal expansion coefficient, ΔT the temperature difference across the layer and d the layer thickness. Convection begins when Ra > ~1000. For the mantle, Ra ≈ 10⁷, so convection is very vigorous by that measure.

Types of Plate Boundaries

The 17 major tectonic plates interact at three kinds of boundary, each producing distinct geological features:

Divergent Boundaries

Plates pull apart. The gap is filled by upwelling magma that solidifies to form new oceanic crust — a process called seafloor spreading. The Mid-Atlantic Ridge is the world's longest mountain chain: 16 000 km of volcanic ridge running down the centre of the Atlantic. Shallow earthquakes and basaltic volcanism are characteristic.

Convergent Boundaries

Plates collide. The denser oceanic plate is forced under the lighter continental plate in a process called subduction. The subducting slab sinks into the mantle, dragging the surface plate with it (slab pull, see below). Subduction zones produce deep-focus earthquakes (down to 700 km), explosive arc volcanoes and accretionary wedges of scraped sediment. The Pacific "Ring of Fire" marks a chain of subduction zones.

Where two continental plates collide, neither subducts easily (both are buoyant). Instead, the crust crumples upward into mountain ranges. The Himalayas formed this way when India collided with Asia ~40–50 million years ago.

Transform Boundaries

Plates slide horizontally past each other. Stress builds as irregularities lock the fault, eventually releasing in large earthquakes when the lock breaks. The San Andreas Fault in California is a classic example — the Pacific Plate moves northwest relative to the North American Plate at ~5 cm/year.

Driving Forces

What actually moves the plates? The answer involves several competing forces:

• Slab pull (dominant, ~90%): Old, cold oceanic lithosphere is denser than the underlying mantle. At subduction zones it sinks under its own weight, pulling the rest of the plate behind it like a tablecloth. This is the strongest single driver.
• Ridge push (~10%): At mid-ocean ridges, hot new crust is elevated and gravitationally slides downhill away from the ridge. A smaller but real contribution.
• Mantle drag: Convection in the asthenosphere exerts viscous drag on the base of plates. Current research suggests it can both drive and resist plate motion depending on the geometry of the convection cell.
• Basal suction: Upwelling mantle plumes (hot spots) under continents can create a suction effect that tears plates apart, initiating rifting.

How the Simulation Works

The browser-based tectonic plate simulation uses a simplified but physically motivated approach. Rather than solving the full Navier–Stokes equations (computationally prohibitive at real dimensions), it uses three coupled models:

1. Voronoi Plate Segmentation

The globe surface is divided into a Voronoi diagram of N plates. Each plate is a convex polygon of mesh triangles. The Voronoi seeds move on the sphere according to velocity vectors sampled from the mantle convection field below.

2. Mantle Convection via 2-D Grid

Below the crust, a 256 × 256 grid stores temperature and velocity. Each frame:

1. Advect temperature: T_new[i] = T_old[i − v*dt] (bilinear interpolation)
2. Diffuse temperature: 5-point Laplacian stencil with κ = 0.01
3. Compute buoyancy velocity from temperature gradient using simplified Stokes solver

The result is a velocity field that drives the plate seeds.

3. Boundary Deformation

At each boundary edge between plates, the relative velocity vector determines the boundary type:

• Divergent: gap opens → new triangles inserted, coloured red (hot crust)
• Convergent: plates overlap → triangles removed, surface buckles upward by offsetting vertex Y positions along the boundary normal
• Transform: vertices slide tangentially, boundary irregularity accumulates as "earthquake energy" released as a particle burst when threshold exceeded

Mountain building is approximated by raising vertex heights proportional to convergence velocity over time, then applying Gaussian smoothing to blend the terrain. The result is visually convincing without being geologically accurate.

Mesh Update Algorithm

Every 60 simulation frames the mesh is re-triangulated via constrained Delaunay triangulation, which prevents degenerate triangles from accumulating at plate boundaries. The key insight is that most vertices stay fixed — only boundary vertices within one plate width of an active edge need retriangulation:

function updateBoundaryMesh(plateA, plateB, velocity) {
  const normal   = boundary.normal;
  const relSpeed = velocity.dot(normal);          // + = convergent, - = divergent
  const dt       = FRAME_TIME * SIM_SPEED;

  if (relSpeed > 0) {
    // Convergent: raise terrain along boundary
    boundary.vertices.forEach(v => {
      v.y += relSpeed * dt * MOUNTAIN_SCALE;
      v.y  = Math.min(v.y, MAX_HEIGHT);
    });
    smoothTerrain(plateA, 3);
    smoothTerrain(plateB, 3);
  } else {
    // Divergent: insert new vertices, assign young-crust colour
    const newVerts = interpolateBoundary(boundary, -relSpeed * dt);
    insertVertices(plateA, newVerts, HOT_CRUST_COLOR);
  }
}

"We can think of the mantle as a giant lava lamp — heated at the bottom by the core, cooled at the top by the surface. The wax blobs drifting up and down are the mantle plumes and downwellings, and the lamp's outer skin is our slowly drifting continents."

Try the Simulation

The interactive simulation lets you watch convection cells form, plates drift and collide, mountains rise and rifts open — all in real time. You can adjust mantle viscosity, core heat flux and plate count, and watch the supercontinent cycle play out in fast-forward.