Spatial and temporal Bayesian hierarchical model over large domains with application to Holocene Sea surface temperature reconstruction in the equatorial Pacific