The current surface boundary condition solves for surface temperature and atmospheric flux, given the magma temperature. This assumes a fixed lid thickness, and a constant dT/dz in the conductive boundary layer.
In reality, the thickness of the boundary layer would evolve with time. We should think about ways to incorporate a solution of both surface temperature and lid thickness, given the mantle temperature and the atmospheric flux calculation.
The current surface boundary condition solves for surface temperature and atmospheric flux, given the magma temperature. This assumes a fixed lid thickness, and a constant dT/dz in the conductive boundary layer.
In reality, the thickness of the boundary layer would evolve with time. We should think about ways to incorporate a solution of both surface temperature and lid thickness, given the mantle temperature and the atmospheric flux calculation.
See this recent paper: https://iopscience.iop.org/article/10.3847/PSJ/ad4351