[1] The scaling of plate tectonic convection is investigated by simulating thermal convection with pseudoplastic rheology and strongly temperature‐dependent viscosity. The effect of mantle melting is also explored with additional depth‐dependent viscosity. Heat flow scaling can be constructed with only two parameters, the internal Rayleigh number and the lithospheric viscosity contrast, the lat...