Modern and fossil pollen data are widely used in paleoenvironmental research to characterize past environmental changes in a given location. However, their discrete and discontinuous nature can limit the inferences that can be made from them. Deriving continuous spatial maps of the pollen presence from point-based datasets would enable more robust regional characterization of such past changes. To address this problem, we propose a comprehensive collection of European pollen presence maps including 194 pollen taxa derived from the interpolation of pollen data from the Eurasian Modern Pollen Database (EMPD v2) restricted to the Euro-Mediterranean Basin. To do so, we developed an automatic Kriging-based interpolation workflow to select an optimal geostatistical model describing the spatial variability for each taxon. The output of the interpolation model consists of a series of multivariate predictive maps of Europe at 25 km scale, showing the occurrence probability of pollen taxa, the predicted presence based on diverse probability thresholds, and the interpolation uncertainty for each taxon. Combined visual inspections of the maps and systematic cross-validation tests demonstrated that the ensemble of predictions is reliable even in data-scarce regions, with a relatively low uncertainty, and robust to complex and non-stationary pollen distributions. The maps, freely distributed as GeoTIFF files (https://doi.org/10.5281/zenodo.10015695, Oriani et al., 2023), are proposed as a ready-to-use tool for spatial paleoenvironmental characterization. Since the interpolation model only uses the coordinates of the observation to spatialize the data, the model can also be employed with fossil pollen records (or other presence/absence indicators), thus enabling the spatial characterization of past changes, and possibly their subsequent use for quantitative paleoclimate reconstructions.
