ABSTRACT The monitoring of the regional (provincial) economic situation is of particular importance due to the high level of heterogeneity and interdependences among different territories. Although econometric models allow for spatial and serial correlation of various kinds, the limited availability of territorial data restricts the set of relevant predictors at a more disaggregated level, especially for gross domestic product (GDP). Combining data from different sources at NUTS-3 level, this paper evaluates the predictive performance of a spatial dynamic panel data model with individual fixed effects and some relevant exogenous regressors, by using data on total gross value added (GVA) for 103 Italian provinces over the period 2000–2016. A comparison with nested panel sub-specifications as well as pure temporal autoregressive specifications has also been included. The main finding is that the spatial dynamic specification increases forecast accuracy more than its competitors throughout the out-of-sample, recognising an important role played by both space and time. However, when temporal cointegration is detected, the random-walk specification is still to be preferred in some cases even in the presence of short panels.