Spatial maps of extreme precipitation are crucial in flood prevention. Withthe aim of producing maps of precipitation return levels, we propose a novelapproach to model a collection of spatially distributed time series where theasymptotic assumption, typical of the traditional extreme value theory, isrelaxed. We introduce a Bayesian hierarchical model that accounts for thepossible underlying variability in the distribution of event magnitudes andoccurrences, which are described through latent temporal and spatial processes.Spatial dependence is characterized by geographical covariates and effects notfully described by the covariates are captured by spatial structure in the hierar-chies. The performance of the approach is illustrated through simulation studiesand an application to daily rainfall extremes across North Carolina (USA). Theresults show that we significantly reduce the estimation uncertainty with respectto state of the art techniques.