---
title: 'Lava: Planetary Molten Silicate Dynamics'
url: https://www.emergentmind.com/topics/lava
type: topic
---

# Lava: Planetary Molten Silicate Dynamics

Lava, in planetary science, designates the manifestation of molten silicates at or near the surfaces of planets and satellites. The scope of "lava" spans thin near-surface flows and lakes observed on contemporary volcanic terrains, through to hemispherical or global "magma oceans"—molten silicate layers hundreds to thousands of kilometers deep—which can dominate the early evolution or persistent states of terrestrial worlds and exoplanets. Lava plays a central role in planetary differentiation, crust formation, atmosphere–interior exchange, and surface–atmosphere observables. Investigations leverage laboratory rheology, thermochemical modeling, direct geological observation, computational hazard mapping, and exoplanetary emission spectroscopy to constrain the underlying physics and evolutionary pathways.

## 1. Definition, Physical Regimes, and Key Parameters

"Lava" refers to silicate melt erupted onto a planetary surface, while "magma ocean" denotes an interior layer where temperature exceeds the silicate liquidus, yielding a fully molten shell of depth $D \sim 10^2$–$10^3$ km. Terrestrial lava flows and lakes typically have thickness $<1$ km and surface temperatures $T_s \sim 1,400$–$1,700$ K. Key physical parameters are:

- Melting/liquidus temperature: $T_m \approx 1,300$ K at 1 bar (peridotite), rising to $\sim2,000$ K at $P \sim 10^2$ GPa.
- Latent heat of crystallization: $L \approx 2 \times 10^6$ J kg$^{-1}$.
- Melt fraction $\phi$ denotes the liquid volume/mass fraction in partially crystallized rock.
- Viscosity $\eta$ for silicate melts can span $10^{-1}$–$10^3$ Pa$\cdot$s, depending on composition and temperature, with superliquidus Mercury analog lavas measured at $4$–$16$ Pa$\cdot$s in the $1,600$–$1,736$ K range [1708.04810], and lower limits near $10^{-1}$ Pa$\cdot$s in silicate vaporization regimes [2012.07337].
- For magma oceans, the Rayleigh number $Ra \sim 10^{18}$ ensures "hard-turbulent" convection.

## 2. Energetics, Dynamics, Crystallization, and Solidification

The thermal evolution of lava and magma-ocean systems is governed by the interplay of radiative cooling, convective transport, latent heat effects during crystallization, and, for planetary-scale melt layers, internal and external heat sources.

- The energy balance for a magma ocean per unit area is:
  \[
  \frac{dE}{dt} = -q_\mathrm{rad} - q_\mathrm{conv} + \rho L \frac{dD_\mathrm{cryst}}{dt}
  \]
  where $q_\mathrm{rad}$ is radiative loss ($\epsilon \sigma T_s^4$), $q_\mathrm{conv}$ is convective flux, and the final term accounts for latent heat released as depth $dD_\mathrm{cryst}$ crystallizes [2012.07337, 2409.11459].

Crystallization proceeds via either batch (equilibrium) or fractional (with segregated phases) processes. The solidus and liquidus are pressure- and composition-dependent:
  \[
  T_\mathrm{sol}(P) \simeq T_0 + m_\mathrm{sol} P, \quad T_\mathrm{liq}(P) \simeq T_1 + m_\mathrm{liq} P
  \]
  where slopes $m_\mathrm{sol} \sim 10$–$15$ K GPa$^{-1}$, $m_\mathrm{liq} \sim 15$–$20$ K GPa$^{-1}$. Convection ceases ("lock-up") near melt fractions $\phi \approx 0.4$–$0.6$.

For planetary lavas, high melt mobility is facilitated by low viscosity, rapid turbulent flow (Re $\gg 2,000$), and high effusion rates. For example, Mercury’s northern volcanic plains exhibit channel velocities $u = 4.8$–$7.2$ m s$^{-1}$ (depths $h = 1$–$10$ m, $\eta_m \approx 13$ Pa$\cdot$s, slope $5^\circ$), and require effusion rates $Q > 10^4$ m$^3$ s$^{-1}$ to traverse $>100$ km before quenching [1708.04810].

## 3. Atmospheric Coupling, Evaporation, and Surface–Atmosphere Exchange

Lava/magma-ocean surfaces exchange mass and energy with their atmospheres via volatile outgassing and rock vaporization, especially on "lava planets"—rocky exoplanets with substellar temperatures $T_s \gtrsim 2,000$ K.

- Vapor pressures $p_i(T_s)$ for surface constituents are well-described by:
  \[
  p_i(T_s) = p_{0,i} \exp\left(-\frac{L_{v,i}}{R T_s}\right)
  \]
  where $L_{v,i}$ is latent heat of evaporation and $R$ the gas constant. Representative $p_i$ at $T_s = 3,000$ K: $p_\mathrm{Na} \approx 1.4 \times 10^4$ Pa, $p_\mathrm{SiO} \approx 5.3 \times 10^3$ Pa, $p_\mathrm{SiO_2} \approx 2.4 \times 10^2$ Pa [2010.14101].

- The resultant silicate vapor atmospheres display significant compositional and scale-height variations: for K2-141b, scale heights $h_\mathrm{Na} \approx 35$ km, $h_\mathrm{SiO} \approx 18$ km, $h_\mathrm{SiO_2} \approx 13$ km at the limb [2010.14101].
- Atmospheric mass-loading, driven by net evaporation, must be balanced by return melt flow in the magma ocean. Steady-state SiO/SiO$_2$ atmospheres require day–night ocean currents $v_o \sim 10^{-4}$ m s$^{-1}$, while Na-rich atmospheres would require over $10^{-2}$ m s$^{-1}$—not sustainable, implying long-term surface fractionation [2010.14101].

Thermal and compositional coupling between the melt layer and atmosphere sets observable properties and drives planetary evolution [2012.07337, 2407.21111].

## 4. Observational Manifestations and Modeling

Lava and associated atmospheres express themselves via distinct observables on Solar System bodies and exoplanets.

- Infrared thermal emission is the primary signature; for lava worlds, secondary-eclipse depths and phase curves yield $T_s$, constrain day–night heat transport, and can reveal atmospheric bands (e.g., Si–O at 9–11 $\mu$m, CO$_2$, Na features) [2012.07337, 2407.21111].
- Clouds in silicate atmospheres (e.g., SiO condensation at lava-planet terminators) reduce surface temperatures by 100–200 K, though their effect on emergent infrared spectra is mainly to flatten the brightness temperature continuum rather than introduce strong lines [2407.21111].
- On terrestrial planets and satellites, lava lakes such as Erta Ale (Ethiopia) and Kīlauea Iki (Hawaii) exhibit convection, crust cycling, and measurable surface heat fluxes (e.g., $\sim20$ W m$^{-2}$ for Erta Ale), serving as scaling analogs for deeper magma-ocean convection [2012.07337].
- On Io, persistent melting is sustained by tidal heating; Loki Patera's convecting magma sea is $\sim200$ km wide with $\sim1$–$3$ W m$^{-2}$ heat flux [2012.07337].

Modern computational tools model lava emplacement as stochastic or probabilistic processes. The Flowy code employs the MrLavaLoba method, simulating lava emplacement as stochastic sequences of elliptical "lobes," achieving computational speedups of $100$–$400\times$ over previous codes with high physical fidelity for hazard and thickness mapping [2405.20144].

## 5. Lava Emplacement: Experimental Constraints and Modeling

Rheological and emplacement experiments empirically constrain flow laws and landscape modification by lava.

- Mercury NVP analog lavas display moderate superliquidus viscosities ($4$–$16$ Pa$\cdot$s for $1,600$–$1,736$ K), and crystal content–dependent, shear-thinning behavior subliquidus (e.g., a decrease in $\eta$ by 1 log unit as shear rate increases from $0.1$ to $5$ s$^{-1}$) [1708.04810].
- Turbulent flow velocities and critical effusion rates required to avoid premature cooling and quenching are derived from numerical models coupling flow kinematics with thermal budget constraints (e.g., maintaining $T$ above liquidus over 100 km requires $Q > 10^4$ m$^3$ s$^{-1}$ for channel widths $\sim100$ m).
- Probabilistic hazard mapping as in Flowy involves Monte Carlo–like ensembles of flows, tracking inundation and thickness fields over digital elevation models (DEMs), with volume-conserving rasterization and convergence metrics (e.g., coverage index $\mu > 0.99$ after 60 runs for $N_\mathrm{thresh} = 4$) [2405.20144].

## 6. Lava Worlds and Exoplanets: Evolution, Feedbacks, and Detection

In ultra-short-period exoplanets, persistent magma oceans and their interactions with irradiation and tidal forces drive both observational signatures and thermal evolution.

- Tidally locked “lava planets” develop hemispherical dichotomy: dayside maintains a thick, hot magma ocean with fixed $T_s$ (e.g., $3,000$ K), while the nightside cools, crystallizes and solidifies over $10^8$–$10^9$ yr without net energy transfer or tidal heating [2409.11459].
- Sustaining a molten nightside requires either: (i) horizontal convection carrying $>$20% of intercepted stellar power (necessitating $\eta_{ns} < 10^{-3}$ Pa$\cdot$s), or (ii) tidal heating rates $q_T \gtrsim 8 \times 10^{-4}$ W/kg in the mush, corresponding to orbital eccentricities $e \gtrsim 7 \times 10^{-3}$ [2409.11459].
- If these criteria are not met, positive feedbacks—viscosity rise and mush mass loss—accelerate nightside freezing.
- Night–day phase curve observations of exoplanets (e.g., K2-141b) provide direct diagnostics of the efficiency of day–night heat transport and thus the physical state of lava and planetary interiors [2409.11459, 2010.14101, 2012.07337].
- Cloud radiative effects in partial atmospheres have been predicted to chill the dayside surface by up to 200 K without substantially altering banded spectral signatures, particularly in systems like HD 213885b and HD 20329b; detection requires absolute spectral precision of 10–20 ppm [2407.21111].

## 7. Knowledge Frontiers and Future Opportunities

Outstanding research needs and frontiers include:

- Rheological properties and phase diagrams of ultra-high-$T$, volatile-poor, and silicate-vaporizing melts at $>2,000$ K and $>100$ GPa [2012.07337].
- Detailed evaporation kinetics, multi-phase convection modeling, and coupled atmosphere–interior chemistry under intense irradiation.
- Direct measurement and mapping of lava world

Source: https://www.emergentmind.com/topics/lava