Schematic representation of the formation of Naxos’s nested domes inspired by the numerical experiments. Partial melting of the crust allows upward migration of buoyant material (melt, magma, felsic partially molten layers represented in white) and accumulation of dense material (residual solids, cumulates, mafic partially molten layers in dark brown). The upward migration of the melting front is increasing the thickness of the convective domain with time, increasing the Rayleigh number and permitting larger convection cells to form and causing an additional decrease in viscosity and density. The first order dome is controlled by the size of the convection cells (drawn in pink) and the second order domes form within the low-density material aggregated at the top of the convection cells (drawn in white).