A Bayesian hierarchical approach for spatial analysis of climate model bias in multi-model ensembles