Modelling the extension of heterogeneous hot lithosphere