We constructed a three-dimensional numerical model of thermal convection in the upper mantle of the Earth in spherical variables using an artificial compressibility method. Results of three-dimensional modeling of convection beneath the Eurasia cratons are presented. The calculation results illustrate the structure of convective flows
