Description: Highlights: • Multi-centennial oscillation with 100–200 years periods is evident in proxy data and model simulations during the Holocene. • Multi-centennial oscillation is a global signal and is more significant in the Northern Hemisphere high latitudes. • None of the external forcings is found to be the sole driver of the multi-centennial variability. • It indicates the multi-centennial oscillation may be due to potential internal drivers and essential feedbacks. Abstract: Variability on centennial to multi-centennial timescales is mentioned as a feature in reconstructions of the Holocene climate. As more long transient model simulations with complex climate models become available and efforts have been made to compile large proxy databases, there is now a unique opportunity to study multi-centennial variability with greater detail and a large amount of data than earlier. This paper presents a spectral analysis of transient Holocene simulations from 9 models and 120 proxy records to find the common signals related to oscillation periods and geographic dependencies and discuss the implications for the potential driving mechanisms. Multi-centennial variability is significant in most proxy records, with the dominant oscillation periods around 120–130 years and an average of 240 years. Spectra of model-based global mean temperature (GMT) agree well with proxy evidence with significant multi-centennial variability in all simulations with the dominant oscillation periods around 120–150 years. It indicates a comparatively good agreement between model and proxy data. A lack of latitudinal dependencies in terms of oscillation period is found in both the model and proxy data. However, all model simulations have the highest spectral density distributed over the Northern hemisphere high latitudes, which could indicate a particular variability sensitivity or potential driving mechanisms in this region. Five models also have differentiated forcings simulations with various combinations of forcing agents. Significant multi-centennial variability with oscillation periods between 100 and 200 years is found in all forcing scenarios, including those with only orbital forcing. The different forcings induce some variability in the system. Yet, none appear to be the predominant driver based on the spectral analysis. Solar irradiance has long been hypothesized to be a primary driver of multi-centennial variability. However, all the simulations without this forcing have shown significant multi-centennial variability. The results then indicate that internal mechanisms operate on multi-centennial timescales, and the North Atlantic-Arctic is a region of interest for this aspect.

Description: Numerical modeling enables a comprehensive understanding not only of the Earth's system today, but also of the past. To date, a significant amount of time and effort has been devoted to paleoclimate modeling and analysis, which involves the latest and most advanced Paleoclimate Modelling Intercomparison Project phase 4 (PMIP4). The definition of seasonality, which is influenced by slow variations in the Earth's orbital parameters, plays a key role in determining the calculated seasonal cycle of the climate. In contrast to the classical calendar used today, where the lengths of the months and seasons are fixed, the angular calendar calculates the lengths of the months and seasons according to a fixed number of degrees along the Earth's orbit. When comparing simulation results for different time intervals, it is essential to account for the angular calendar to ensure that the data for comparison are from the same position along the Earth's orbit. Most models use the classical calendar, which can lead to strong distortions of the monthly and seasonal values, especially for the climate of the past. Here, by analyzing daily outputs from multiple PMIP4 model simulations, we examine calendar effects on surface air temperature and precipitation under mid-Holocene, Last Interglacial, and pre-industrial climate conditions. We came to the following conclusions. (a) The largest cooling bias occurs in boreal autumn when the classical calendar is applied for the mid-Holocene and Last Interglacial, due to the fact that the vernal equinox is fixed on 21 March. (b) The sign of the temperature anomalies between the Last Interglacial and pre-industrial in boreal autumn can be reversed after the switch from the classical to angular calendar, particularly over the Northern Hemisphere continents. (c) Precipitation over West Africa is overestimated in boreal summer and underestimated in boreal autumn when the classical seasonal cycle is applied. (d) Finally, month-length adjusted values for surface air temperature and precipitation are very similar to the day-length adjusted values, and therefore correcting the calendar based on the monthly model results can largely reduce the artificial bias. In addition, we examine the calendar effects in three transient simulations for 6–0 ka by AWI-ESM, MPI-ESM, and IPSL-CM. We find significant discrepancies between adjusted and unadjusted temperature values over continents for both hemispheres in boreal autumn, while for other seasons the deviations are relatively small. A drying bias can be found in the summer monsoon precipitation in Africa (in the classical calendar), whereby the magnitude of bias becomes smaller over time. Overall, our study underlines the importance of the application of calendar transformation in the analysis of climate simulations. Neglecting the calendar effects could lead to a profound artificial distortion of the calculated seasonal cycle of surface air temperature and precipitation.