AbstractThe spatiotemporal variations of the climate of the last millennium provide the context for understanding the recent climate. However, there are few 1000‐year‐long spatiotemporally continuous climate reconstructions at hemispheric to global scales. Here, a millennial, spatiotemporally continuous and annually resolved air temperature dataset for the Northern Hemisphere (NH) was reconstructed using paleoclimate data assimilation (PDA) approach based on an expanded tree ring proxy dataset and artificial neural network‐based proxy system models. Verifications show the PDA‐based reconstruction agrees well with instrumental temperature, previous proxy‐based reconstructions and the Last Millennium Climate Reanalysis version 2 (LMRv2) dataset in space and time domains. Since the PDA‐based reconstruction is spatiotemporally continuous fields without missing value over the whole NH, which is a distinct advantage of the PDA‐based reconstruction over other proxy‐based NH‐averaged temperature reconstructions, we focussed on analysing the spatial variations of the NH temperature during the last millennium, for example, the primary modes of temperature variabilities and the spatial variabilities of the NH temperature in three typical climatic epochs of the millennium. Results show that the first‐leading mode (explained variance of 72.50%) of the NH temperature variability is a spatially consistent mode, while the second‐leading mode (explained variance of 10.60%) exhibits a different pattern that is different from the first‐leading empirical orthogonal function mode, particularly in the Arctic, which shows a dipole structure anomaly. The 20th century was the warmest, with 94.17% of the NH displaying positive anomalies. The 11th century was the second warmest century in the record, with 70.26% of the NH showing a positive anomaly. The 19th century was the coldest, with negative anomalies covering 86.97% of the NH. Northwest North America and the Northern Europe‐Barents Sea are two distinctly anomalous centres during these three centuries. Overall, the reconstructed air temperatures provide a new and reliable dataset for studying the millennial NH temperature variations.