Radiogenically heated lava-lake crust model

For a well-mixed lava lake of depth $H$ and a thin conductive crust of thickness $a$, put $\Delta T=T_m-T_s$, and use the convective flux $K(\overline{T}-T_m)$ with constant heat transfer coefficient $K>0$. The Stefan condition and liquid heat balance give, to leading order in $a/H$,Here latent heat $L_f$ is per unit mass and $k=\rho c_p\kappa$. Crust heating and changes of liquid depth are neglected. The linear conductive profile uses a quasi-steady approximation; thinness alone does not justify neglecting crust thermal storage. During initial growth, a small Stefan number supplies a sufficient rapid-adjustment regime. Define $h=a/H$, $\theta =(\overline{T}-T_m)/\Delta T$, and $\tau =\kappa t/H^2$. Scaling by the thermal diffusion time $H^2/\kappa$ and $\Delta T$ gives $h^{\prime }=S(1/h-B\theta )$ and ${\theta }^{\prime }=R-B\theta$, where $B=KH/k$, $R=Q{H}^2/(c_p\kappa \Delta T)$, and $S=c_p\Delta T/L_f$ is the Stefan number. Thus $\theta =R/B+({\theta }_0-R/B)e^{-B\tau }$ exactly. The reduction applies while $a/H\ll 1$, the conductive profile is approximately linear, and the thermal convection closure remains valid; it does not model complete freezing of a deep lake.

In the radiogenically heated lava-lake crust model, the thickness derivative vanishes at $h_e=1/(B\theta )$ when the liquid superheat $\theta$ is positive. This moving curve balances thermal conduction and convective heat flux; it is not itself the evolving crust. Cooling raises the curve and keeps an initially absent crust below it. Warming lowers the curve, allowing lava-crust overshoot under warming. At zero superheat the stationary thickness is infinite, so its finite formula applies only after positive superheat develops.

In the radiogenically heated lava-lake crust model, initially cool liquid warms exponentially and its instantaneous stationary lava-crust thickness $1/(B\theta )$ decreases. The actual crust can cross this moving value, reach a maximum, and melt back to $h_{\ast }=1/R$. Overshoot is not automatic. Linearizing about equilibrium givesIf $B<S{R}^2$ and ${\theta }_0<R/B$, the slowly decaying forced response has positive sign, so the crust eventually approaches equilibrium from above; an initially absent crust must then have overshot. In contrast, initially hot liquid cools, its instantaneous stationary thickness increases, and an initially absent crust cannot cross it from below: at a proposed first crossing its own derivative is zero but the stationary thickness has positive derivative. The hot-start crust therefore grows monotonically. These comparisons assume nonnegative liquid superheat.

In the radiogenically heated lava-lake crust model, temperature stays constant precisely when ${\theta }_0=R/B$. Then $h^{\prime }=S(1/h-R)$, and an initially absent crust satisfies the displayed law for $0\le h<1/R$. It follows by integrating $h/(1-Rh)$ from zero to $h$. The crust starts as $\sqrt{2S\tau }$ and tends to $1/R$ only at infinite time. The equilibrium remains thin only if $R\gg 1$.