The nanodevices often work in a thermal environment and lead to the viscoelastic and thermoelastic coupling which becomes more complicated when further coupled with the size effects. A spatial-temporal fractional order model is proposed to study the dynamic behavior of thermoelastic nanoplate in the present work. The integer order differential is extended to the fractional order differential with the integer order differential as a special case. This extension makes the modeling of small-scale effects of mechanical behavior and generalized heat conduction more flexible. Non-local strain gradient elasticity and non-local heat conduction are used to describe the spatial non-local effects while the fractional order differential to describe the complicated history-dependent effects or the temporal nonlocal effects. The standard solid viscoelastic model and generalized heat conduction model with the fractional order differentials lead to the fractional-order differential governing equations. The fractional order differential equation is solved by the Laplace transform method, and the analytical solution of the response is expressed in terms of the Mittag-Leffler function. In order to verify the reliability of the analytical solution, the numerical solution is also provided for comparison. Based on the numerical results of dynamic response under the step load, the effects of thermoelastic fractional order parameters, viscoelastic fractional order parameters, and small-scale parameters are discussed. The new spatial-temporal fractional order model is a natural extension of the traditional integer order model which can be recovered from the present model.