Supplementary material to "Calendar effects on surface air temperature and precipitation based on model-ensemble equilibrium and transient simulations from PMIP4 and PACMEDY"

Abstract. Numerical modelling 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 is from the same position along the Earth's orbit. Most models use the classical "fixed-length" 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 conclude that: (a) The largest cooling bias occurs in autumn when the classical calendar is applied for the mid-Holocene and last interglacial. (b) The sign of the temperature anomalies between the Last Interglacial and pre-industrial in boreal autumn can be reversed after the switch from 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 "fixed-length" seasonal cycle is applied. (d) Finally, correcting the calendar based on the monthly model results can reduce the biases to a large extent, but not completely eliminate them. In addition, we examine the calendar effects in 3 transient simulations for 6–0 ka by AWI-ESM, MPI-ESM, and IPSL. We find significant discrepancies between adjusted and unadjusted temperature values over ice-free 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 "fixed-length" 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. One important fact to be noted here is that the discrepancy in seasonality under different calendars is an analysis bias and is highly depends on the choice of the reference position/date (usually the vernal equinox, which is set to 31th March) on the Earth's ellipse around the sun. Different model groups may apply different reference dates, so ensuring a consistent reference date and seasonal definition is key when we compare results across multiple models.

[1]  C. Brierley,et al.  Analysing the PMIP4-CMIP6 collection: a workflow and tool (pmip_p2fvar_analyzer v1) , 2022, Geoscientific Model Development.

[2]  P. Valdes,et al.  The PMIP4 Last Glacial Maximum experiments: preliminary results and comparison with the PMIP3 simulations , 2021 .

[3]  B. Otto‐Bliesner,et al.  Large-scale features of Last Interglacial climate: results from evaluating the lig127k simulations for the Coupled Model Intercomparison Project (CMIP6)–Paleoclimate Modeling Intercomparison Project (PMIP4) , 2021 .

[4]  C. Brierley,et al.  Analyzing the PMIP4-CMIP6 ensemble: a workflow and tool (pmip_p2fvar_analyzer v1) , 2021 .

[5]  H. Curtis Orbital position as a function of time , 2021, Orbital Mechanics for Engineering Students.

[6]  M. Toohey,et al.  Global temperature modes shed light on the Holocene temperature conundrum , 2020, Nature Communications.

[7]  K. Wyser,et al.  Supplementary material to "Simulating the mid-Holocene, Last Interglacial and mid-Pliocene climate with EC-Earth3-LR" , 2020 .

[8]  Zhenghui Xie,et al.  The Flexible Global Ocean‐Atmosphere‐Land System Model Grid‐Point Version 3 (FGOALS‐g3): Description and Evaluation , 2020, Journal of Advances in Modeling Earth Systems.

[9]  J. Dufresne,et al.  Implementation of the CMIP6 Forcing Data in the IPSL‐CM6A‐LR Model , 2020, Journal of Advances in Modeling Earth Systems.

[10]  D. Sidorenko,et al.  Early-Holocene simulations using different forcings and resolutions in AWI-ESM , 2020, The Holocene.

[11]  Bo Sun,et al.  A multi-model CMIP6-PMIP4 study of Arctic sea ice at 127 ka: sea ice data compilation and model differences , 2021, Climate of the Past.

[12]  A. LeGrande,et al.  Large-scale features and evaluation of the PMIP4-CMIP6 midHolocene simulations , 2020, Climate of the Past.

[13]  M. Losch,et al.  Evaluation of FESOM2.0 Coupled to ECHAM6.3: Preindustrial and HighResMIP Simulations , 2019, Journal of Advances in Modeling Earth Systems.

[14]  J. Aerts,et al.  Agreement between reconstructed and modeled boreal precipitation of the Last Interglacial , 2019, Science Advances.

[15]  M. Mills,et al.  High Climate Sensitivity in the Community Earth System Model Version 2 (CESM2) , 2019, Geophysical Research Letters.

[16]  Guoxiong Wu,et al.  CAS FGOALS-f3-L Model Datasets for CMIP6 Historical Atmospheric Model Intercomparison Project Simulation , 2019, Advances in Atmospheric Sciences.

[17]  P. Braconnot,et al.  Strengths and challenges for transient Mid- to Late Holocene simulations with dynamical vegetation , 2019, Climate of the Past.

[18]  M. Toohey,et al.  What was the source of the atmospheric CO2 increase during the Holocene? , 2019, Biogeosciences.

[19]  S. V. Emelina,et al.  Simulation of the modern climate using the INM-CM48 climate model , 2018, Russian Journal of Numerical Analysis and Mathematical Modelling.

[20]  P. Bartlein,et al.  Paleo calendar-effect adjustments in time-slice and transient climate-model simulations (PaleoCalAdjust v1.0): impact and strategies for data analysis , 2018, Geoscientific Model Development.

[21]  Michel Crucifix,et al.  The PMIP4 contribution to CMIP6 – Part 1: Overview and over-arching analysis plan , 2018 .

[22]  Bin Wang,et al.  The NUIST Earth System Model (NESM) version 3: description and preliminary evaluation , 2017, Geoscientific Model Development.

[23]  O. Boucher,et al.  IPSL IPSL-CM6A-LR model output prepared for CMIP6 CMIP , 2018 .

[24]  Michel Crucifix,et al.  The PMIP 4 contribution to CMIP 6 – Part 1 : Overview and overarching analysis plan , 2018 .

[25]  J. Schmitt,et al.  A 156 kyr smoothed history of the atmospheric greenhouse gases CO 2 , CH 4 , and N 2 O and their radiative forcing , 2017 .

[26]  Xiaoping Zhou,et al.  The NUIST Earth System Model (NESM) version 3: description and preliminary evaluation , 2017 .

[27]  W. Lipscomb,et al.  The PMIP4 contribution to CMIP6 – Part 2: Two interglacials, scientific objective and experimental design for Holocene and Last Interglacial simulations , 2016 .

[28]  G. Lohmann,et al.  Simulated response of the mid-Holocene Atlantic meridional overturning circulation in ECHAM6-FESOM/MPIOM , 2016 .

[29]  T. Semmler,et al.  Towards multi-resolution global climate modeling with ECHAM6-FESOM. Part II: climate variability , 2013, Climate Dynamics.

[30]  G. Lohmann,et al.  Greenland Ice Sheet influence on Last Interglacial climate: global sensitivity studies performed with an atmosphere–ocean general circulation model , 2015 .

[31]  Thomas F. Stocker,et al.  Revision of the EPICA Dome C CO2 record from 800 to 600 kyr before present , 2015 .

[32]  D. Jiang,et al.  Mid-Holocene global monsoon area and precipitation from PMIP simulations , 2015, Climate Dynamics.

[33]  J. Schröter,et al.  Towards multi-resolution global climate modeling with ECHAM6–FESOM. Part I: model formulation and mean climate , 2015, Climate Dynamics.

[34]  J. N. Tokis A Solution of Kepler’s Equation , 2014 .

[35]  K. Cook,et al.  West African monsoon demise: Climatology, interannual variations, and relationship to seasonal rainfall , 2014 .

[36]  A. Timmermann,et al.  The Holocene temperature conundrum , 2014, Proceedings of the National Academy of Sciences.

[37]  F. Joos,et al.  A reconstruction of atmospheric carbon dioxide and its stable carbon isotopic composition from the penultimate glacial maximum to the last glacial inception , 2013 .

[38]  Shaun A Marcott,et al.  A Reconstruction of Regional and Global Temperature for the Past 11,300 Years , 2013, Science.

[39]  A. Berger,et al.  Modelling the climatic diversity of the warm interglacials , 2012 .

[40]  A. Berger,et al.  The last interglacial (Eemian) climate simulated by LOVECLIM and CCSM3 , 2012 .

[41]  B. Otto‐Bliesner,et al.  A multi-model assessment of last interglacial temperatures , 2012 .

[42]  J. Kutzbach,et al.  Calendar effect on phase study in paleoclimate transient simulation with orbital forcing , 2011 .

[43]  R. S. Thompson,et al.  Pollen-based continental climate reconstructions at 6 and 21 ka: a global synthesis , 2011 .

[44]  T. Stocker,et al.  Atmospheric nitrous oxide during the last 140,000 years , 2010 .

[45]  C. Turney,et al.  Does the Agulhas Current amplify global temperatures during super‐interglacials? , 2010 .

[46]  V. Masson‐Delmotte,et al.  TALDICE-1 age scale of the Talos Dome deep ice core, East Antarctica , 2010 .

[47]  T. Stocker,et al.  Glacial–interglacial and millennial-scale variations in the atmospheric nitrous oxide concentration during the last 800,000 years , 2010 .

[48]  J. Jungclaus,et al.  Effects of orbital forcing on atmosphere and ocean heat transports in Holocene and Eemian climate simulations with a comprehensive Earth system model , 2009 .

[49]  A. Timmermann,et al.  On the definition of seasons in paleoclimate simulations with orbital forcing , 2008 .

[50]  T. Stocker,et al.  Orbital and millennial-scale features of atmospheric CH4 over the past 800,000 years , 2008, Nature.

[51]  F. Joos,et al.  Rates of change in natural and anthropogenic radiative forcing over the past 20,000 years , 2008, Proceedings of the National Academy of Sciences.

[52]  T. Fichefet,et al.  Results of PMIP2 coupled simulations of the Mid-Holocene and Last Glacial Maximum – Part 2: feedbacks with emphasis on the location of the ITCZ and mid- and high latitudes heat budget , 2007 .

[53]  J. Tison,et al.  One-to-one coupling of glacial climate variability in Greenland and Antarctica. , 2006 .

[54]  Epica Community Members One-to-one coupling of glacial climate variability in Greenland and Antarctica , 2006, Nature.

[55]  T. Stocker,et al.  Atmospheric Methane and Nitrous Oxide of the Late Pleistocene from Antarctic Ice Cores , 2005, Science.

[56]  Kenji Kawamura,et al.  Evidence for substantial accumulation rate variability in Antarctica during the Holocene, through synchronization of CO 2 in the Taylor Dome, Dome C and DML ice cores , 2004 .

[57]  N. Noblet-Ducoudré,et al.  Some mechanisms of mid-Holocene climate change in Europe, inferred from comparing PMIP models to data , 2004 .

[58]  S. Janicot,et al.  The West African Monsoon Dynamics. Part II: The “Preonset” and “Onset” of the Summer Monsoon , 2003 .

[59]  J. Kutzbach,et al.  Mid-Holocene climates of the Americas: a dynamical response to changed seasonality , 2003 .

[60]  D. Pollard,et al.  A calendar conversion method for monthly mean paleoclimate model output with orbital forcing , 2002 .

[61]  W. Broecker,et al.  Last Interglacial Climates , 2002, Quaternary Research.

[62]  T. Stocker,et al.  High‐resolution Holocene N2O ice core record and its relationship with CH4 and CO2 , 2002 .

[63]  T. Stocker,et al.  Atmospheric CO2 concentrations over the last glacial termination. , 2001, Science.

[64]  Daryl R. Myers,et al.  Mean preserving algorithm for smoothly interpolating averaged data , 2001 .

[65]  Pascale Braconnot,et al.  Sensitivity of paleoclimate simulation results to season definitions , 1997 .

[66]  Schwander,et al.  Atmospheric CO2 . . . , 1992 .

[67]  E. S. Epstein On Obtaining Daily Climatological Values from Monthly Means , 1991 .

[68]  André Berger,et al.  Insolation values for the climate of the last 10 , 1991 .

[69]  J. M. A. Danby,et al.  The solution of Kepler's equation, I , 1983 .

[70]  André Berger,et al.  Long-term variations of daily insolation and Quaternary climatic changes , 1978 .

[71]  André Berger,et al.  Long-term variations of the Earth's orbital elements , 1977 .