A new multivariable benchmark for Last Glacial Maximum climate simulations

Abstract. We present a new global reconstruction of seasonal climates at the Last Glacial Maximum (LGM, 21 000 years BP) made using 3-D variational data assimilation with pollen-based site reconstructions of six climate variables and the ensemble average of the PMIP3—CMIP5 simulations as a prior (initial estimate of LGM climate). We assume that the correlation matrix of the uncertainties in the prior is both spatially and temporally Gaussian, in order to produce a climate reconstruction that is smoothed both from month to month and from grid cell to grid cell. The pollen-based reconstructions include mean annual temperature (MAT), mean temperature of the coldest month (MTCO), mean temperature of the warmest month (MTWA), growing season warmth as measured by growing degree days above a baseline of 5  ∘ C (GDD 5 ), mean annual precipitation (MAP), and a moisture index (MI), which is the ratio of MAP to mean annual potential evapotranspiration. Different variables are reconstructed at different sites, but our approach both preserves seasonal relationships and allows a more complete set of seasonal climate variables to be derived at each location. We further account for the ecophysiological effects of low atmospheric carbon dioxide concentration on vegetation in making reconstructions of MAP and MI. This adjustment results in the reconstruction of wetter climates than might otherwise be inferred from the vegetation composition. Finally, by comparing the uncertainty contribution to the final reconstruction, we provide confidence intervals on these reconstructions and delimit geographical regions for which the palaeodata provide no information to constrain the climate reconstructions. The new reconstructions will provide a benchmark created using clear and defined mathematical procedures that can be used for evaluation of the PMIP4–CMIP6 entry-card LGM simulations and are available at https://doi.org/10.17864/1947.244 (Cleator et al., 2020b).