A geodynamic model of mantle density heterogeneity