Cold Extremes during Dansgaard-Oeschger Oscillations
Abstract. This paper studies the statistics of extreme cold air surface temperatures in a climate that experiences abrupt changes. We employ a CCSM4 simulation of Last Glacial Maximum conditions that exhibits rapid switching between stadial and interstadial states as data. Non-stationary linear Generalized Extreme Value (GEV) distributions are fitted to find the regions that show the most prominent changes and associate them with physical processes. The results are then compared with a non-linear model, which gives a more detailed picture of how the parameters change as a function of the AMOC strength. While different regions of the world show different extremal behaviours, many regions show an approximately linear relationship between the parameters of the GEV distributions and the strength of the AMOC within the stadial and interstadial states, with some non-linear oscillation or jump where the transition between them takes place. Physical processes such as the expansion and retreat of the sea ice and the relative changes in the strength of the currents are shown to impact the magnitude and variability of the extremes, with significant changes observed in the three parameters of the GEV. They also create teleconnections that are compared whenever possible with various proxies for the temperature of the air and the water. We show how mapping the parameters of the GEV distributions into the AMOC strength gives a way to compare results between different climate models and different climate states. Comparisons, however, need to pay heed to the dynamical characteristics of the state and the location to be meaningful.
The authors study the extreme value statistics of a low resolution CCSM4 simulation with Last Glacial Maximum conditions, containing Dansgaard-Oeschger (DO) events. To model the minima the authors employ a non-stationary generalized extreme value (GEV) distribution, where the location, scale and shape parameters are expressed using linear predictors of the AMOC strength as a covariate. Different configurations of this model are fitted independently to all local time series, and the dependence of the parameters on AMOC strength is further explored by fitting a stationary GEV model on a rolling window. Notable regions provide insight into the relationship between the AMOC strength and the distribution of annual minimum temperatures, and are highlighted and discussed.
The manuscript presents a thorough analysis of how cold extremes respond to abrupt climate transitions, an interesting topic that is well within the scope of ESD. However, I have one main concern regarding the applicability of extreme value theory to the dataset, and I therefore recommend a major revision before the manuscript can be considered for publication. My main concern, together with some smaller comments and technical corrections, is detailed below.
While the dataset used is almost 8000 years long, it is sampled at monthly resolution. This means that the number of observations in each block is 12, which is very small. In addition, seasonal effects would likely make annual minima cluster near the coldest winter months, further reducing the effective number of observations in each block. My primary concern is therefore that extreme value theory might not be appropriate for this dataset, since the GEV distribution is an asymptotic result.
If it is feasible, I would recommend repeating the analysis using a CCSM4 simulation output with higher temporal resolution (e.g. daily). Otherwise, I believe it is necessary to validate the GEV assumptions, for example by showing which months tend to provide the annual minimum across the domain (to check the effective block size) and by providing goodness-of-fit diagnostics (e.g. Q-Q plots) for a selection of grid points, and to compare the GEV with alternative distributions such as the normal or skew-normal (e.g. by AIC or another appropriate criterion). If the GEV model does not fit well, it might be necessary to consider alternative models, and to reinterpret the results.
Other comments:
Technical corrections: