Flow regimes as a function of the unmolten and partially molten domains Rayleigh numbers (Ra_UM and Ra_PM), after 20 Myr of basal heating of a 45 km thick continental crust containing white buoyant and low-viscosity material and black heavy and high-viscosity material (inclusions). The ambient domain remains motionless as long as RaPM < 200 (with at most local segregation of some inclusions one with respect to a neighbour), diapirism occurs when 200 < RaPM <3000, and the suspension regime initiates when RaUM > 10 (convection). The layering regime occurs at low RaUM and RaPM > 3000 : convection is not vigorous enough so that the buoyant inclusions can stack at ca. 20 km depth (after Louis-Napoléon et al., 2022).