Present day shallow and deep density anomaly distribution predicted by models A (panel a₂) and B (panel b₂), with crust density of 2670 kg m⁻³. Dashed black lines indicates the 300 and 1500 K isotherms. Panels (a₁) and (b₁) shows the corresponding surface gravity anomalies predicted by, respectively, model A and model B (continuum black lines), compared to the gravity anomalies due to solely thermal contraction and assuming no crust–mantle density stratification.