Although there has been significant progress in the seismic imaging of mantle heterogeneity,the outstanding issue that remains to be resolved is the unknown distribution of mantle temperature anomalies in the distant geological past that give rise to the present-day anomalies inferred by global tomography models.To address this question,we present 3-D convection models in compressible and self-gravitating mantle initialised by different hypothetical temperature patterns.A notable feature of our forward convection modelling is the use of self-consistent coupling of the motion of surface tectonic plates to the underlying mantle flow,without imposing prescribed surface velocities(i.e.,plate-like boundary condition).As an approximation for the surface mechanical conditions before plate tectonics began to operate we employ the no-slip(rigid) boundary condition.A rigid boundary condition demonstrates that the initial thermally-dominated structure is preserved,and its geographical location is fixed during the evolution of mantle flow.Considering the impact of different assumed surface boundary conditions(rigid and plate-like) on the evolution of thermal heterogeneity in the mantle we suggest that the intrinsic buoyancy of seven superplumes is most-likely resolved in the tomographic images of present-day mantle thermal structure.Our convection simulations with a plate-like boundary condition reveal that the evolution of an initial cold anomaly beneath the Java-Indonesian trench system yields a long-term,stable pattern of thermal heterogeneity in the lowermost mantle that resembles the presentday Large Low Shear Velocity Provinces(LLSVPs),especially below the Pacific.The evolution of subduction zones may be,however,influenced by the mantle-wide flow driven by deeply-rooted and longlived superplumes since Archean times.These convection models also detect the intrinsic buoyancy of the Perm Anomaly that has been identified as a unique slow feature distinct from the two principal LLSVPs.We find there is no need for dense chemical ’piles’ in the lower mantle to generate a stable distribution of temperature anomalies that are correlated to the LLSVPs and the Perm Anomaly.Our tomography-based convection simulations also demonstrate that intraplate volcanism in the south-east Pacific may be interpreted in terms of shallow small-scale convection triggered by a superplume beneath the East Pacific Rise. Although there has been significant progress in the seismic imaging of mantle heterogeneity, the outstanding issue that remains to be resolved is the unknown distribution of mantle temperature anomalies in the distant geological past that give rise to to present-day anomalies inferred by global tomography models. To address this question, we present 3-D convection models in compressible and self-gravitating mantle initialised by different hypothetical temperature patterns. A notable feature of our forward convection modeling is the use of self-consistent coupling of the motion of surface tectonic plates to the underlying mantle flow, without imposing prescribed surface velocities (ie, plate-like boundary conditions). As an approximation for the surface mechanical conditions before plate tectonics began to operate we employ the no-slip (rigid) boundary condition. A rigid boundary condition demonstrates that the initial thermally-dominated structure is preserved, and its geographical location is fixed during the evolution of mantle flow. Copleidering the impact of different assumed surface boundary conditions (rigid and plate-like) on the evolution of thermal heterogeneity in the mantle we suggest that the intrinsic buoyancy of seven superplumes is most-likely resolved in the tomographic images of present-day mantle thermal structure. Our convection simulations with a plate-like boundary condition reveal that the evolution of an initial cold anomaly beneath the Java-Indonesian trench system yields a long-term, stable pattern of thermal heterogeneity in the lowermost mantle that resembles the presentday Large Low Shear Velocity Provinces (LLSVPs), especially below the Pacific.The evolution of subduction zones may be, however, influenced by the mantle-wide flow driven by deeply-rooted and longlived superplumes since Archean times.These convection models also detect the intrinsic buoyancy of the Perm Anomaly that has been identified as a unique slow feature distinct from the two principal LLSVPs.We find there is no need for dense chemical ’piles’ in the lower mantle to generate stable distribution of temperature anomalies that are correlated to the LLSVPs and the Perm Anomaly. Our tomography-based convection simulations also demonstrate that intraplate volcanism in the south-east Pacific may be interpreted in terms of shallow small-scale convection triggered by a superplume beneath the East Pacific Rise.