We present probabilistic projections for spatial patterns of future temperature change using a multivariate Bayesian analysis. The methodology is applied to the output from 21 global coupled climate models used for the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. The statistical technique is based on the assumption that spatial patterns of climate change can be separated into a large scale signal related to the true forced climate change and a small scale signal due to model bias and variability. The different scales are represented via dimension reduction techniques in a hierarchical Bayesian model. Posterior probabilities are obtained with a Markov chain Monte Carlo simulation. We show that with 66% (90%) probability 79% (48%) of the land areas warm by more than 2°C by the end of the century for the SRES A1B scenario.