Patterns of changing surface climate variability from the Last Glacial Maximum to present in transient model simulations

As of 2023, global mean temperature has risen by about 1.45±0.12 °C with respect to the 1850–1900 pre-industrial (PI) baseline according to the World Meteorological Organization. This rise constitutes the first period of substantial global warming since the Last Deglaciation, when global temperatures rose over several millennia by about 4.0–7.0 °C according to proxy reconstructions. Similar levels of warming could be reached in the coming centuries considering current and possible future emissions. Such warming causes widespread changes in the climate system, of which the mean state provides only an incomplete picture. Instead, fluctuations around the mean and in higher-order statistics need to be considered. Indeed, climate's variability and the distributions of climate variables change with warming, impacting, for example, ecosystems and the frequency and intensity of extremes. However, previous investigations of climate variability focus mostly on measures such as variance, or standard deviation, and on quasi-equilibrium states such as the Holocene or Last Glacial Maximum (LGM). Changes in the tails of distributions of climate variables and transition periods such as the Last Deglaciation remain largely unexplored. Therefore, we investigate changes of climate variability on annual to millennial timescales in 15 transient climate model simulations of the Last Deglaciation. This ensemble consists of models of varying complexity, from an energy balance model to Earth system models (ESMs), and includes sensitivity experiments, which differ only in terms of their underlying ice sheet reconstruction, meltwater protocol, or consideration of volcanic forcing. The ensemble simulates an increase in global mean temperature of 3.0–6.6 °C between the LGM and Holocene. Against this backdrop, we examine whether common patterns of variability emerge in the ensemble. To this end, we compare the variability in surface climate during the LGM, Deglaciation, and Holocene by estimating and analyzing the distributions and power ...