Interglacials of the last 800 , 000 years

Interglacials, including the present (Holocene) period, are warm, low land ice extent (high sea level), end-members of glacial cycles. Based on a sea level definition, we identify eleven interglacials in the last 800,000 years, a result that is robust to alternative definitions. Data compilations suggest that despite spatial heterogeneity, Marine Isotope Stages (MIS) 5e (last interglacial) and 11c (~400 ka ago) were globally strong (warm), while MIS 13a (~500 ka ago) was cool at many locations. A step change in strength of interglacials at 450 ka is apparent only in atmospheric CO2 and in Antarctic and deep ocean temperature. The onset of an interglacial (glacial termination) seems to require a reducing precession parameter (increasing Northern Hemisphere summer insolation), but this condition alone is insufficient. Terminations involve rapid, nonlinear, reactions of ice volume, CO2, and temperature to external astronomical forcing. The precise timing of events may be modulated by millennial-scale climate change that can lead to a contrasting timing of maximum interglacial intensity in each hemisphere. A variety of temporal trends is observed, such that maxima in the main records are observed either early or late in different interglacials. The end of an interglacial (glacial inception) is a slower process involving a global sequence of changes. Interglacials have been typically 10–30 ka long. The combination of minimal reduction in northern summer insolation over the next few orbital cycles, owing to low eccentricity, and high atmospheric greenhouse gas concentrations implies that the next glacial inception is many tens of millennia in the future.


CITE THIS VERSION
http:// hdl.handle.net/2078.1/175429 --DOI : 10.1002/2015RG000482 Le dépôt institutionnel DIAL est destiné au dépôt et à la diffusion de documents scientifiques émanant des membres de l'UCLouvain. Toute utilisation de ce document à des fins lucratives ou commerciales est strictement interdite. L'utilisateur s'engage à respecter les droits d'auteur liés à ce document, principalement le droit à l'intégrité de l'oeuvre et le droit à la paternité. La politique complète de copyright est disponible sur la page Copyright policy DIAL is an institutional repository for the deposit and dissemination of scientific documents from UCLouvain members. Usage of this document for profit or commercial purposes is stricly prohibited. User agrees to respect copyright about this document, mainly text integrity and source mention. Full content of copyright policy is available at Copyright policy become warmer in the next century and remain warm for many millennia to come [Intergovernmental Panel on Climate Change (IPCC), 2013]. This makes particularly pertinent the study of periods in which at least sectors of the Earth system may have been "warmer" than they are currently.
Our understanding of the climate of our current interglacial is good in a static sense: given the boundary conditions (solar irradiance, astronomical characteristics, geography, ice cover, and concentrations of long-lived greenhouse gases such as CO 2 ) and the resulting radiation balance of the planet. However, we cannot currently give a satisfactory explanation of how the climate system evolved to this state: why we live in an interglacial at this time, or indeed why we live in an interglacial embedded in 100 ka rather than 40 ka cycles.
Finally the entries into interglacials (and exits from them) are fascinating examples of nonlinear behavior of the Earth system that need to be understood. For these reasons-importance of warm periods, need to understand how our current climate arose, fascinating climate dynamics, and recent availability of significant amounts of data and model simulations-it is timely to review and synthesize existing knowledge about the warm periods (interglacials) of the last 800 ka.
After identifying what an interglacial is, this review aims to summarize the common features and differences between them. This allows us to highlight particular interglacials that may enlighten us about the climate system and to describe current knowledge of the factors that lead into and out of an interglacial. The review ends with a discussion about the likely extent of the present interglacial, both in the absence and presence of human influence on the climate system.

Nomenclature and History
Terrestrial, marine, and ice core sequences deposited since the intensification of NH glaciation,~2.6 Ma ago have traditionally been divided on the basis of climate changes inferred from their lithological, biological, and geochemical signature. The fundamental units of this geologic-climate classification [American Commission on Stratigraphic Nomenclature, 1961] are (i) glaciations (or glacials), representing episodes during which extensive glaciers developed, and (ii) interglaciations (or interglacials), representing episodes during which the climate was incompatible with the wide extent of glaciers. Further subdivisions include stades (stadials) and interstades (interstadials), reflecting secondary advances and recessions or standstills of glaciers, respectively.
The term "interglacial" derives from the position of temperate-climate sediments interbedded between glacial sediments [Heer, 1865;Geikie, 1874]. In the higher midlatitudes of Europe and North America the concept of an interglacial as an interval of time during which conditions were at least as warm as the present interglacial (Holocene) was traditionally based on floral evidence, indicating the development of temperate deciduous forest [e.g., Jessen and Milthers, 1928;Fairbridge, 1972;Woillard, 1978;West, 1984;Gibbard and West, 2000]. By comparison, the term interstadial has been used to describe a period that was either too short or too cold to allow the development of temperate deciduous forest. Similar distinctions were also adopted for faunal evidence [e.g., Coope, 1977].  [Lisiecki and Raymo, 2005]. Oxygen isotopes in benthic foraminifera represent a combination of changes in the oxygen isotope content of the water (dominated by ice volume) and deepwater temperature. This plot shows the onset of cyclicity, when Northern Hemisphere ice started to wax and wane, at~2.6 Ma ago (arrow marks the official start of the Quaternary), and the strengthening and lengthening of cycles across the mid-Pleistocene transition (horizontal bar).

10.1002/2015RG000482
In deep-sea sediments, the terms glacial and interglacial were first applied to even-and odd-numbered oxygen isotope stages, respectively, defined by variations in the oxygen isotope ratio (δ 18 O) in planktonic foraminifera and first thought to reflect changes in sea surface temperature [Emiliani, 1955]. However, the concept of a marine isotopic stratigraphy is based on the principle that changes in the oxygen isotopic composition of the oceans primarily arises from changes in the volume of ice stored on the continents [Shackleton, 1967] and that given the slow mixing of the deep ocean, such a signal would be global on timescales of thousands of years. Shackleton and Opdyke [1973] proposed the formalization of Marine Isotope Stages (MIS) 1-22, using the δ 18 O in foraminifera in core V28-238 from the Eastern Equatorial Pacific, a practice which has been gradually extended over the entire Pleistocene, Pliocene, and beyond. Marine isotope stages may be further divided into substages (e.g., MIS 5 into 5a-5e [Shackleton, 1969]), characterized by local minima and maxima in δ 18 O values known as marine isotopic events [Prell et al., 1986], such that MIS 5a contains event 5.1, MIS 5b contains event 5.2, and so on. The decimal event notation has occasionally been used to denote substages, but the two systems are not interchangeable. An event refers to a point in time rather than an interval; in a detailed record there will be space between successive events, whereas substages have boundaries and together fill the stratigraphic record and the time that it represents [Shackleton, 2006]. To avoid confusion, please note that we use a recently recommended substage nomenclature [Railsback et al., 2015], which contrasts with some literature usage, particularly for MIS 9 and MIS 15.
A major development in the understanding of the correspondence between marine and terrestrial stages was the realization that the Eemian (Last) Interglacial of northwest Europe is not equivalent to the whole of MIS 5, but rather approximates to a substage (MIS 5e) within it [Shackleton, 1969]. In line with the terrestrial evidence, other isotopically depleted marine substages corresponded to interstadials associated with the same interglacial complex (so that MIS 5c and 5a would broadly correspond to the Amersfoort/Brörup and Odderade interstadials, respectively, of the northwest European stratigraphy). Joint foraminiferal isotopic and pollen analyses from marine sequences have not only confirmed such correlations but also showed that the marine and terrestrial stage boundaries are not necessarily time parallel, reflecting geographical asynchroneity and differences in response times of, e.g., vegetation and ice [Sánchez Goñi et al., 1999;Shackleton et al., 2002Shackleton et al., , 2003Tzedakis et al., 2004].

Definition of Interglacial
The population of interglacials must be defined before any attempt is made to discuss their individual characteristics, including the statistics of their timing and spacing. Given that interglacial complexes (odd-numbered Marine Isotope Stages, such as MIS 5, but excluding MIS3 which is not now considered to include an interglacial period) contain a number of temperate intervals within them, a crucial question is how to categorize interglacials and distinguish them from interstadials. In some conceptual models [Parrenin and Paillard, 2003], the climate system occupies (and jumps between) a glacial and deglacial (or interglacial) state, with different governing equations in each. While a bifurcation structure may indeed be relevant for describing the dynamics of glacial-interglacial cycles, identifying it from the available paleo-observations is complex [Imbrie et al., 2011]. Distinguishing interglacial from interstadial climates is nonetheless important, not only on formal stratigraphic grounds but also in climatic context. Though temperate in climatic signature, interstadials are still part of a glacial world with excess ice in the Northern Hemisphere outside Greenland. Thus, for example, the Amersfoort/Brörup and Odderade interstadials in northwestern European nomenclature (roughly MIS 5a and 5c) are part of the Early Glacial (MIS 5d-5a) in this region, which reflects not only their stratigraphic position but also their overall climatic character compared to the Eemian.
Central to any definition of an interglacial is the choice of an appropriate criterion that can be applied consistently through time. A number of approaches are possible, each with inherent advantages and disadvantages (Table 1). A recurring issue is that distinctions between interglacials and interstadials depend on the metric used. It is desirable that whatever metric is chosen has global significance, such that interglacials are phenomena that can be considered widespread (probably of global extent), even if their regional expression is neither globally uniform or synchronous. Some definitions that have been used are based on temperature. Temperature records exist at many sites but are essentially local, and therefore, the definition of an interglacial based on temperature would vary according to the location. This will remain the case until there are sufficient data with accurate chronologies, such that some kind of global average temperature can be constructed across multiple glacial/interglacial cycles. Only recently has a first such reconstruction [Shakun et al., 2012;Marcott et al., 2013] become available for the last 22 ka and none is available for earlier interglacials.
We might also be tempted to define interglacials on the basis of the external forcing. At the multimillennial timescale, the dominant external forcing of Earth's climate is the astronomical forcing. This drives the climate system dynamics which, lead, for example, to variations in atmospheric CO 2 concentration and in the mass of the ice sheets. These in turn may be viewed as feedbacks, or as "internal" forcing factors of the atmospheric and oceanic conditions all over the globe. Increased concentrations of CO 2 , in particular, are a characteristic of interglacials. However, the relationships between astronomical forcing, CO 2 and the mass of continental ice are nonlinear and complex. In particular, eccentricity, often cited as a possible origin of the 100 ka periodicity characteristic of glacial-interglacial cycles does not appear in the spectrum of regional insolation changes; it merely modulates the amplitude of insolation changes. Combinations of obliquity and precession predict the occurrence of some interglacials [Kukla et al., 1981], but as longer paleoclimate records have appeared, it has become clear that such simple formulae do not correctly predict the full roster observed in the data. It is therefore not currently practical to use astronomical forcing to define interglacials.
As discussed earlier, interglacials have been defined as episodes during which global climate was incompatible with the wide extent of glaciers [American Commission on Stratigraphic Nomenclature, 1961]. While remembering that there is still a large volume of ice in Antarctica during every interglacial, this nonetheless Lack of objective criteria to recognize such a bifurcation from the existing records "As warm or warmer than the Holocene" Unambiguous in that a demarcation value is set (e.g., Holocene peak) In absence of sufficient data to calculate global mean temperature, definition will vary with location.
With this definition, no interglacials between 800 and 450 ka ago Fixed threshold (similar to above, but cutoff value is not determined by Holocene values) Unambiguous for given metric Depends on metric used (e.g., temperature, sea level, and CO 2 ).
Depending on where threshold is set, it can exclude "cool" interglacials or could include "interstadials" (e.g., MIS 5a and 5c) With specific reference to Antarctic records, depends on whether "overshoots" are included in the interglacial (see section 7.1) Varying threshold, for example, different pre-MBE versus post-MBE values, or one for each g-ig cycle If amplitude of signal is genuinely nonstationary when subjected to similar controls, then this method addresses that.
Subjective choice of time period over which threshold can change-makes unproven assumption about nonstationarity.
(alternatively, segments of time series are normalized and a fixed threshold, e.g., one standard deviation is applied across) With specific reference to Antarctic records, depends on whether "overshoots" are included in the interglacial (see section 7.1) Interglacial defined as the interval following a glacial termination Good for calculating repeat times of g-ig cycles Shifts problem to defining what a termination is.
Interglacial defined as the most prominent peak(s) within each odd-numbered marine isotopic complex Allows for more than one interglacial within one complex Depends on metric used (e.g., temperature, sea level, and CO 2 ).
Makes an unwarranted assumption that original numbering system satisfactorily defines interglacial recurrence. Interglacials are characterized by absence of NH ice outside Greenland; different interglacials must be separated by lowering of sea level below a set threshold.
Uses a metric (global sea level/ice volume) that integrates global climate effects Uncertainties with sea level reconstructions Consistent with stratigraphic concepts Subjective choice of remaining NH ice that can be tolerated and of ice volume increase to separate interglacials.
Allows for more than one interglacial within one complex a It should be noted that some of these definitions are also appropriate for estimating interglacial duration. The difference, however, is that the identification of interglacials using a threshold approach usually is made with reference to its acme (peak), whereas estimating duration does not necessarily consider only peak conditions.

10.1002/2015RG000482
emphasizes that the fundamental concept underlying the terminology of an interglacial is that of low continental ice volume, a consequence of integrated global climate effects, which would also be observed as high sea level or less positive seawater δ 18 O. This is in line with early definitions of interglacials based on the marine transgression [e.g., Harting, 1874]. The slow response of ice sheets would mean that the sea level highstand or isotopic minimum may be attained relatively late compared to increases in warmth or climate forcing. Here we clearly distinguish the use of such values to define whether an interglacial is occurring, from any measure used to define when the interglacial starts and ends. As an example, for the current interglacial, sea level close to present, a value we might use to define MIS 1 as an interglacial, occurs as late as~7-8 ka ago. However, having determined that an interglacial is occurring, its start can be defined in other ways [Walker et al., 2009].
A eustatic, or global mean, sea level curve would therefore provide a useful measure to address this issue, but in practice direct sea level determinations are often complicated by local isostatic and gravitational effects and dating uncertainties, are discontinuous, and only cover the most recent part of the record [e.g., Thompson and Goldstein, 2006]. The δ 18 O from benthic foraminifera is sometimes used as a sea level proxy, but it is overprinted by local deepwater temperature and hydrographic effects, which can obscure the glacioeustatic signal [Chappell and Shackleton, 1986;Skinner and Shackleton, 2005]; averaging many globally distributed δ 18 O benthic records as in the stacked record of Lisiecki and Raymo [2005] (hereafter referred to as LR04) does however remove some of the regional variability. Various approaches to isolate the sea level component of δ 18 O benthic [Shackleton, 1987;McManus et al., 1999;Shackleton, 2000;Waelbroeck et al., 2002] have been attempted, but until recently none of them have been extended to cover the last 800 ka. A modeled sea level curve by Bintanja and van de Wal [2008] extends over the last 3 Ma but depends on uncertain assumptions about deepwater temperatures and their coupling with atmospheric temperatures. A record using shallow-infaunal benthic foraminifera to deconvolve deepwater temperatures (using Mg/Ca ratios) and the δ 18 O of seawater (δ 18 O sw ) over the last 1.5 Ma has recently become available from Ocean Drilling Program (ODP) site 1123 on the Chatham Rise, east of New Zealand [Elderfield et al., 2012]. With hydrographic effects considered minimal at this site, the δ 18 O sw provides an important insight into ice volume changes over this interval, although the propagated error of ±0.2‰ in calculated δ 18 O sw estimated by Elderfield et al. [2012] means that the uncertainty is ±20 m sea level equivalent. A continuous sea level reconstruction based on planktonic δ 18 O records and a hydrological model for the Red Sea is available for the last 520 ka [Rohling et al., 2009]. A related method covering a longer time period, using a hydrological model for the Mediterranean [Rohling et al., 2014], gives a further useful overview, but with a large propagated uncertainty in sea level, similar to that of Elderfield et al. [2012]. Finally, recent work [Shakun et al., 2015] has used a global compilation of paired planktonic δ 18 O and sea surface temperatures in an attempt to deduce δ 18 O of surface seawater, interpreted to reflect mainly ice volume and sea level. However, it is hard to use this record to compare the strength of interglacials because the authors had to make a correction for a long-term trend, believed to be unrelated to sea level.
Given that the LR04 stack is affected by deepwater temperature and hydrographic effects of variable influence through time, we first examine those records where the sea level component has been (mostly) isolated. The ODP Site 1123 [Elderfield et al., 2012], the Red Sea/Mediterranean [Rohling et al., 2009[Rohling et al., , 2014, and the detrended surface seawater stack [Shakun et al., 2015] records all suggest that within the uncertainties of the reconstructions (±20 m), highstands during the most prominent isotopic substages of the last 800 ka reached a level similar to present ( Figure 2). Assuming that the East Antarctic ice sheet was not significantly smaller than today, this observation points to a more specific definition of an interglacial as an interval within which the distribution of Northern Hemisphere ice resembled the present (0 ± 20 m), i.e., there was little Northern Hemisphere ice outside Greenland, with periods of significantly greater ice volume (sea level passing below approximately À50 m compared to present) before and after the interglacial period. The second half of the definition is required in order to separate periods of interglacial character that are interrupted by clearly glacial conditions. However, it has consequences for the number of interglacials in MIS 15 and 7. The distinction between an interglacial and an interstadial under this definition lies in the absence/presence of significant Northern Hemisphere ice outside Greenland. An example of this distinction is MIS 5. Evidence from the Norwegian Sea and western Scandinavia indicates ice-free conditions during the Eemian (interglacial), but expansion of mountain glaciers during the Amersfoort/Brørup (MIS 5c) and Odderade (MIS 5a) interstadials leading to increased background values in the relative abundance of ice-rafted detritus in the Norwegian Sea during these intervals [Baumann et al., 1995]. Tzedakis et al. [2009] assigned an interglacial status to the most prominent temperate (low δ 18 O) interval within each odd-numbered marine isotope stage of the LR04 benthic δ 18 O stack [Lisiecki and Raymo, 2005] (with the exception of the historical artifact of MIS 3). On this basis, MIS 5e, 9e, 11c, 13a, 17c, and 19c (and by default MIS 1) are unambiguously the most prominent intervals within each MIS (Figure 2), but choosing the most prominent interval in MIS 7 and MIS 15 is difficult. Yin andBerger [2010, 2012] chose a similar set of interglacials to simulate, using the expectation of a 100 ka cycle to justify including MIS 7e but not 7a-7c, and including MIS 15a but not 15e. Using our sea level definition, MIS 7a-7c and 7e, and also MIS 15a and 15e are all equally prominent within their respective stages, and as we do not assume a particular periodicity, they are all included in our interglacial roster, which therefore comprises MIS 1, 5e, 7a-7c (as a single interglacial), 7e, 9e, 11c, 13a, 15a, 15e, 17c, and 19c.
Although a definition based on sea level, as we have suggested, is attractive in its simplicity, its application is hampered by uncertainties that remain in the values of sea level provided by indirect but continuous methods. In particular, direct sea level evidence suggests [Dutton and Lambeck, 2012;Kopp et al., 2013] that sea level during MIS 5e was higher than present. However, the (ODP 1123) derived sea level record [Elderfield et al., 2012] places the MIS 5e highstand at the edge of the lowest range of our proposed threshold (À20 m), and at a similar level in MIS 5a, with MIS 5c higher than both ( Figure 2). Had we only the ODP1123 data [Elderfield et al., 2012], then MIS 5a and MIS 5c would have been identified as interglacials on a par with MIS 5e, but we use knowledge from other sea level data sets [Dutton and Lambeck, 2012;Grant et al., 2012;Kopp et al., 2013] to distinguish MIS 5e as an interglacial and MIS 5a and 5c as interstadials. This, however, casts doubt on the reliability of the reconstruction in earlier periods (especially before 450 ka ago) when alternative sea level evidence is weaker. Given this uncertainty, we also test the sensitivity of our interglacial assignments by using simple thresholds in the uncorrected benthic oxygen isotope record (in this case the LR04 stack [Lisiecki and Raymo, 2005]), and in the ice core record of CO 2 (chosen as a global-scale indicator which, together with ice volume/extent, affects the global radiative forcing).
Our sensitivity studies with the LR04 benthic isotope stack [Lisiecki and Raymo, 2005] would produce the same set of interglacials for any threshold between 3.5 and 3.73‰ (to qualify as an interglacial, values lower than the threshold must be achieved). Any threshold set lower than 3.5‰ would exclude first MIS 17c and very quickly all other stages except MIS 1, 5e, 9e, and 11c. At thresholds higher than 3.73‰, several other stages, starting with MIS 5c, would become interglacials. The very wide range of thresholds (3.5-3.73‰) that lead to the same answer suggest that this solution is robust, i.e., that new data are unlikely to make it obsolete.
We need also to assess the criterion for how deep the surrounding glacials must be to decide that a peak qualifies as a separate interglacial. Using LR04, it would in theory be possible to set a second threshold O of benthic foraminifera in the LR04 stack [Lisiecki and Raymo, 2005]; (bottom) probability maximum Red Sea relative sea level [Rohling et al., 2009] (orange line); and deconvolved δ 18 O of seawater [Elderfield et al., 2012] (grey line). Marine Isotopic substages of interglacial status are indicated (as well as interstadials MIS 5a and 5c). Note that there is some confusion in the literature about the lettering of substages in MIS 9 and 15. Here we use the substage numbering recommended in a recent paper [Railsback et al., 2015], such that the oldest substages in MIS 9 and 15 are 9e and 15e. This contrasts with the usage in some recent papers [Tzedakis et al., 2012b]. Dashed lines mark values discussed in the text that might be used in the definition of interglacials.
for the depth of the valley between interglacials, such that MIS 13 and 15 are separate interglacials, while MIS 7a-7e and MIS 15a-15e each form single very long interglacials. However, the range of values for such a second threshold is very narrow (~0.1‰), suggesting that it has no robust validity. In addition, it would imply inclusion in an interglacial of isotopic values (during MIS 7d and 15b) that are generally characteristic of periods that are clearly part of the glacial world. A final argument against this solution is that two of the derived sea level records [Elderfield et al., 2012;Shakun et al., 2015] place sea level lower in MIS 7d than in MIS 14. Thus, unless we are willing to describe MIS 13-15 as a very long continuous interglacial, we are forced to set the threshold such that MIS 7ac and 7e (similarly 15a and 15e) are separate interglacials.
As a final sensitivity study, we investigate the ice core CO 2 record [Lüthi et al., 2008;Bereiter et al., 2015]; because of its importance as a radiative forcer at global scale, CO 2 may plausibly be one of the controlling factors on interglacial strength in a range of other records (see section 6.7). If the threshold is 260 ppm or higher, only one interglacial (MIS 19) is recognized before MIS 11. However, as the threshold is lowered other interglacials qualify: by 250 ppm, all the interglacials in our roster (including 7c and 15a) have qualified, with the exception of MIS 17c. Periods of clearly interstadial nature in other records get drawn in for lower thresholds, starting with MIS 5a. The low values of CO 2 in MIS 17c were already noted [Lüthi et al., 2008] and persist even in the revised data set recently published [Bereiter et al., 2015]. Perhaps most importantly for this sensitivity study, there is again no reasonable criterion using CO 2 that qualifies MIS 13a as an interglacial without also qualifying MIS 7ac and 15a as separate interglacials.
The scheme we have used allows for the occurrence of more than one interglacial within the same complex and does not require an interglacial to follow a traditionally numbered glacial termination. Only a definition that insists on a minimum temporal spacing (of order 100 ka) between interglacials could force a more traditional list of interglacials. While we recognize that the definition applied here will lead to some inconsistency in relation to traditional nomenclature, we feel it is justified by its objectivity and robustness against a range of definitions and thresholds; the traditional assumptions on the other hand will lead to misleading conclusions about the spacing and timing of interglacials. It is important to reiterate that this is a working definition based on a global metric and as such it has global stratigraphic significance. Interglacial stages may continue to be defined in local contexts using local signals (e.g., local temperature or forest development), but they do not carry general applicability. One cautionary note to add is that our definition is clearly dependent on the context of the last 800 ka and would not necessarily be suitable for the early Quaternary, even though we recognize from benthic isotopes that periods with interglacial character existed then, nor for earlier periods containing alternations of icy and ice-free intervals.

Methods
Before attempting to assess the nature, intensity, and timing of interglacials, it is necessary to establish and summarize the main methods used to measure interglacial climate characteristics, to place them on an age scale, and to model the processes involved.

Stratigraphy and Chronology
Regardless of what definition is adopted for an interglacial, the study of past interglacials must be underpinned by a robust stratigraphy and chronology. Comparing climate signals in different archives is challenging and requires precise stratigraphic correlation of marine, ice, and terrestrial records, including those providing information about sea level and atmospheric composition. Study of the initiation, termination, and duration of interglacials as well as rates of processes is entirely dependent on the accuracy of the timescale. The age model controls the rate, shape, and duration of interglacial periods and absolute ages are needed to discuss the phasing of proxies relative to astronomical forcing. Sampling frequency as well as chronological precision can affect the degree of variability and peak values (amplitude) of the interglacial signal.
The long ice core records have relied mostly on the Antarctic Ice Core Chronology 2012 (AICC2012) based on information from several ice cores [Bazin et al., 2013], and its predecessor, the European Project for Ice Coring in Antarctica (EPICA) Dome C, third generation (EDC3) chronology that was derived by inversion of a snow accumulation and mechanical ice flow model using a set of independently dated horizons (age markers) along the core [Parrenin et al., 2007]. Note that all ages shown in figures in this paper are nominally relative to 1950 as present. The uncertainty of the EDC3 age model is estimated to be ±3 ka at 100 ka ago and Reviews of Geophysics 10.1002/2015RG000482 increasing to ±6 ka from 130 ka to the base of the core. AICC2012 [Bazin et al., 2013] shows only small differences from EDC3 and is well within the quoted age uncertainty. The duration of the last four interglacial periods is not affected by more than 5% using AICC2012 instead of EDC3 [Bazin et al., 2013]. One powerful alternative method for developing ice core chronologies is orbital tuning of the O 2 /N 2 ratio to local insolation forcing [Bender, 2002;Kawamura et al., 2007]. This method has the advantage that it dates the ice (which holds the climate signal), rather than the enclosed air, which can be considerably younger because the air bubbles close only at depths of typically 60-100 m). With an assumption that O 2 /N 2 is in phase with local summer insolation, this method provides a link to absolute ages; however, the precision of this assumption has been questioned . Vostok O 2 /N 2 data are included in AICC2012 as one component of the age model, but the Dome Fuji (DFO2006) age model relies mainly on this tuning, and a smaller estimated uncertainty (up to 2.5 ka) is quoted. Comparison of EDC3 and Dome Fuji chronologies reveals reasonable agreement for the last 350 ka (within the estimated uncertainty of EDC3, but not always within the uncertainty quoted for DFO2006).
Most deep-sea records are dated by correlating benthic δ 18 O signals to the LR04 stack [Lisiecki and Raymo, 2005]. The LR04 stack is tuned to the output of an ice volume model [Imbrie and Imbrie, 1980] that, in turn, is driven by insolation, assuming fixed response times. The estimated error for LR04 is ±4 ka for the past million years. One of the caveats associated with the LR04 timescale is its reliance on assumptions associated with astronomical forcing of climate. The depth-derived ages of Huybers and Wunsch [2004] are almost free of such assumptions but are estimated to be accurate only to within ±9 ka. The differences between LR04 and other age models fall within the uncertainties of these error estimates (Figure 3). (see discussion below) and agrees closely with the timing of Antarctic temperature change ( Figure 3). This likely reflects the late response of ice volume compared to temperature in the benthic δ 18 O signal. The lead of warming Antarctic and deep-sea temperatures over (NH) ice volume change underscores the problem of how to define the start of interglacial periods because the early warming may not be considered to be part of the interglacial period if defined by ice volume reduction.
Speleothems and corals provide an opportunity for absolute dating of the timing of terminations [e.g., Cheng et al., 2009] and the duration of interglacials by radiometric U series dating, and to determine the phasing of climate change relative to insolation forcing. If age markers from speleothems can be identified as unequivocally synchronous to features in ice cores and/or marine sediment cores, then the speleothem chronology can be exported to other archives [e.g., Barker et al., 2011;Grant et al., 2012] and used to evaluate the accuracy of independent chronologies. Full use of such methods awaits firm proof that rapid jumps in local speleothem δ 18 O are synchronous with rapid jumps in ice core temperature or atmospheric methane and marine sea surface temperature (SST) records. Assumptions of this nature have been disputed [e.g., Blaauw et al., 2010]. Barker et al. [2011] placed the EPICA Dome C (EDC) record on a "Speleo-Age" timescale for the last 400 ka by first differentiating the Antarctic δD record, according to a conceptual model of interhemispheric millennial temperature change, to produce a synthetic Greenland record, and then correlating cold events in Greenland with "weak monsoon events" in the detrended speleothem δ 18 O record. The ages agree to within ±2-3 ka with EDC3 from present to MIS 11 where the difference increases to 5 ka. A similar approach has been applied to Iberian Margin sediments where planktonic δ 18 O and SST have been used to synchronize millennial-scale events among marine sediment core, Antarctic ice core, and speleothem records [Barker et al., 2011;Hodell et al., 2013]. Differences between the speleo-synchronized and LR04 age models from the Iberian Margin are less than ±5 ka for the last 400 ka.
While speleothem records can be radiometrically dated, the age models of other long terrestrial records have been derived by tuning lithostratigraphical or biostratigraphical variations to an astronomical target (which may be a local insolation target) or to marine records [Prokopenko et al., 2006;Tzedakis et al., 2006;Melles et al., 2012;Torres et al., 2013]. While this can be done rather precisely, it involves assumptions about synchroneity that could be in error by several thousand years. Unless these assumptions can be tested (for example, through the use of pollen in marine cores, of tephra layers that can be shown to have the same origin, or through the use of dated paleomagnetic boundaries), it is impossible to assess the phasing of change between different processes during relatively rapid changes that occur at terminations in particular.
The accuracy of timescales directly influences our interpretations of the causes of interglacial periods. For example, the errors associated with the ages for the onset and end of an interglacial period (regardless of defining criteria) are large relative to the accuracy of astronomical solutions, which we take to have no error for the last 800 ka. The absolute chronological error (of order ±5 ka for ice cores and marine records) is also large relative to the duration of interglacial periods, which ranges from~10 to 30 ka. The age uncertainty (about one quarter of a precession cycle) provides considerable leeway in interpretation of how the onset and end of interglacials align with astronomical forcing, underscoring the importance of reducing chronological error in the future.

Proxies
Temperature, greenhouse gas (GHG) concentrations, and ice volume (sea level) have been the three reconstructed climate parameters used most often for studying past interglacial periods, although changes in the hydrologic cycle are also important especially at low latitudes and changes in ocean circulation are also highly relevant. Clearly, it is not appropriate here to review all the methods used in paleoclimate reconstruction; however, we briefly summarize the main methods used to derive the records that appear most often in this paper.
In the polar regions, temperature is derived from the ratios of water isotopes (δ 18 O, δD) in ice cores. This has a strong physical basis as a temperature proxy [Jouzel et al., 1997], but there remains uncertainty about its calibration, particularly for higher temperatures [Sime et al., 2009]. Additional methods to constrain temperature in ice cores have been used [e.g., Buizert et al., 2014] in Greenland ice cores, but so far the water isotopes are the only method used to reconstruct temperatures in the long Antarctic ice cores that extend through several glacial cycles. SSTs are estimated from marine sediment cores using a range of geochemical and biological proxies: fossil assemblages (using foraminifera, coccolithophores, radiolarians, diatoms, or dinocysts), alkenone, and other biomarkers (U k 37 , U k′ 37 , TEX 86 ), Mg/Ca ratios and "clumped" isotopologs of carbon and oxygen.

10.1002/2015RG000482
Each of these methods comes with uncertainty about calibration and about the precise interpretation of the proxy (for example, whether it indicates an annual or seasonal signal). Where methods have been compared, discrepancies are found in some regions [Nurnberg et al., 2000;de Vernal et al., 2006;de Garidel-Thoron et al., 2007]. This means that it is difficult to construct consistent spatial time slices of absolute temperature. However, the relative signal within a single proxy at a single site should be much more reliable, so that anomaly data can be used with some confidence to assess the relative strength of different interglacials, or the timing of maxima and transitions at a given location. For terrestrial sites, pollen assemblages may sometimes be used to estimate temperature; many other commonly described terrestrial proxies are influenced by factors other than temperature.
GHG concentrations (particularly, CO 2 and CH 4 ) for the last 800 ka can be determined directly from measurements in the air within ice cores [e.g., Loulergue et al., 2008;Lüthi et al., 2008]. The uncertainty in the concentrations measured in this way is very small (a few parts per million for CO 2 and about 10 parts per billion for CH 4 ), although methodology is important for obtaining correct values, especially in deep ice [Bereiter et al., 2015]. However, in the cores for which old ice is available, there is an inherent limit to the resolution (of several centuries), and an uncertainty in the age relative to the temperature proxy in the same core, because of changes in delta age, the age difference between the ice and the younger air that diffuses from the surface and is trapped within the ice matrix only at depth [e.g., Parrenin et al., 2013]. CH 4 changes result from changes in wetlands in the tropics and high northern latitudes and thus are often used as an indicator in the gas phase of the timing of Northern Hemisphere and tropical climate change relative to other signals.
Estimates of relative sea level change at individual sites and dates can be made using a range of methods (already partly discussed in section 2.2). While dated corals are of particular importance, to obtain a continuous (over multiple glacial cycles) and global, eustatic, sea level record, marine sediments remain important. The δ 18 O of benthic foraminifera is affected not only by changes in continental ice volume but also by deepwater temperature and local hydrographic (water mass) changes. On timescales significantly greater than the mixing time of the oceans (10 3 years), the δ 18 O calcite signal represents changes in the mean state of the "global ocean" as local/regional hydrographic changes that exist in the ocean are diminished [Shackleton, 1967;Imbrie et al., 1984;Skinner and Shackleton, 2006]. If an independent estimate of temperature is derived, the seawater δ 18 O component (i.e., δ 18 O seawater ) of the benthic δ 18 O signal can be calculated. As discussed in section 2.2, one approach has been to stack benthic δ 18 O records assuming that local and regional hydrographic signals are averaged out, thereby leaving the common global δ 18 O signal [Imbrie et al., 1984;Prell et al., 1986;Lisiecki and Raymo, 2005]. In reality, however, the relative contributions of temperature and ice volume still remain unconstrained in a stacked record. Another approach has been to measure benthic δ 18 O in cores from regions where deepwater temperature is close to freezing [e.g., Labeyrie et al., 1987]. Others have attempted to model δ 18 O seawater and sea level from the benthic δ 18 O stack [Bintanja et al., 2005;Bintanja and van de Wal, 2008] or use a sea level transfer function to convert benthic δ 18 O to relative sea level equivalent [Waelbroeck et al., 2002]. Mg/Ca ratios of foraminifera have the potential to reconstruct temperature and, when combined with foraminiferal δ 18 O measurements, provide an estimate of changes in the δ 18 O of seawater [Elderfield et al., , 2012. This approach was taken to deconvolve δ 18 O seawater and deepwater temperature at Site ODP1123 in the deep SW Pacific. A completely separate approach for obtaining a continuous (albeit local, and therefore relative) sea level record has been used in a series of papers using marine cores from the Red Sea. This relies on an assumption that δ 18 O in foraminifera in these cores is controlled mainly by changes in salinity induced by evaporation, whose influence varies depending on the depth of seawater above a shallow sill that controls exchange with the open ocean. The Red Sea records so far extend only to just over 500 ka ago [Rohling et al., 2009], although a longer record using a similar method for the Mediterranean has been published recently [Rohling et al., 2014]. Different sea level estimates are compared in Figure 2.

Models
A range of different types of model has been used to study past interglacials. Experiments can be carried out with a range of underlying purposes, for example: testing models that are used to simulate the present and future climate, understanding the dynamics of termination and inception of glacials (i.e., the start and end of interglacials), or understanding the whole dynamics of glacial-interglacial cycles.
By far the most common approach is to carry out snapshot simulations of particular time slices, using models ranging from full ESMs and atmosphere-ocean general circulation models (GCMs), in some cases identical to those used to study present and future climate, through to EMICs. The majority of such simulations have been  , which focused on time slices at 130, 128, 125, and 115 ka. Such simulations are driven by boundary conditions distinct from those of the Holocene. For most interglacial simulations, ice sheet extent is assumed similar to that of today, as are other components such as vegetation and aerosol concentration. Therefore, the main changes are in astronomical parameters (which can be calculated precisely) and GHG concentrations (which are derived from ice core data); these can be modified as appropriate for the time period under study. A very few studies have carried out snapshots for earlier interglacials Herold et al., 2012;Yin and Berger, 2012;Milker et al., 2013;Muri et al., 2013;Yin and Berger, 2015]. The evaluation of snapshot simulations is severely hampered by the scarcity of data compilations for previous interglacials. As an example, the PMIP3 LIG intercomparison was evaluated against a compilation that consisted of a single value of temperature for each site [Turney and Jones, 2010] for the whole interglacial, which clearly makes the separate evaluation of snapshots at 125 and 130 ka impossible. Work is underway to improve this situation for a number of interglacials, including MIS 5e [Capron et al., 2014] and MIS 11 [Milker et al., 2013;Kleinen et al., 2014].

Reviews of Geophysics
To evaluate the mechanisms at the start and the end of interglacials, and eventually to understand the interactions during interglacials between ice sheet history, climate, vegetation, and greenhouse gases, transient simulations are essential. However, due to limitations of computer resources, it remains challenging to run even simplified GCMs for many thousands of years (remember that a full simulation of the last interglacial would need to run at least from 135 to 110 ka). A few EMICs, and GCMs running at low resolution (and in some cases with accelerated astronomical forcing), have been deployed across the last interglacial (130-115 ka) [Bakker et al., 2013], and transient simulations of earlier interglacials (in this case 9e, 11c, and 19c) are starting to appear [Yin and Berger, 2015]. Another development that may allow a detailed analysis of the sensitivity of climate (including that of interglacials) to different boundary conditions (such as orbital parameters, CO 2 , and ice sheet size) is the use of an emulator that is based on, and attempts to mimic, a GCM or EMIC [Araya-Melo et al., 2015;Bounceur et al., 2015].
At least two EMICs have been run across the full 800 ka period we are considering in this paper, one with prescribed [Holden et al., 2010[Holden et al., , 2011, the other with interactive [Ganopolski and Calov, 2011], ice sheets, and both with prescribed CO 2 concentration. While these models necessarily have extremely simplified dynamics and other simplifications, they do allow the exploration of the influence of different factors (such as GHG concentrations, insolation, and dust) on the timing and occurrence of interglacials. An alternative way to look at what influences glacial cycles is to drive a state of the art ice sheet model with parameterized climates, which are themselves developed from GCM simulations [Abe-Ouchi et al., 2013]. This latter study covered a period of the last 400 ka. A range of studies have investigated the loss and growth of ice sheets across one or more glacial cycles using a variety of simplified forcings [e.g., Berger et al., 1998b].
A final class of models that should be mentioned seeks to understand the dynamics of glacial-interglacial cycles through simplified mathematical concepts [Paillard, 1998;Crucifix, 2013, and references therein]. Although these conceptual models are addressing an issue much wider than just interglacial climate, we include them here because they form part of the discussion in sections 7 and 8.

Forcings
Because of the effect of components of the climate system with a long response time, the state of the climate at any given time depends not only on the instantaneous forcing but also the forcing history. This section focuses on the instantaneous forcing during interglacials, while the impact of forcing history during the preceding glacial is considered in later sections. This will provide a context for discussions about the relative strength of interglacials, their timing, and the trends that occur during them.

Insolation Characteristics due to Astronomical Changes
The largest difference in forcing between and during different interglacials lies in the latitudinal and seasonal pattern of incoming solar radiation [Yin andBerger, 2012, 2015], which is controlled by the three astronomical parameters precession, obliquity, and eccentricity ( Figure 4). Here we consider first the instantaneous forcing at a characteristic time in each interglacial, and then the trends and patterns of astronomical forcing through the interglacial.

Reviews of Geophysics
10.1002/2015RG000482 Yin and Berger [2010] discussed the insolation patterns at the precession minimum (NH summer insolation maximum) near to what they considered to be the "peak" of each interglacial. In the absence of robust data sets for global temperature or ice volume, they identified the period of interest using the LR04 benthic δ 18 O stack [Lisiecki and Raymo, 2005]. They chose this data set because it has some global character, and peak values (isotopic minima) are easily identified. Within the chronological uncertainty, the LR04 dates for the interglacial peaks (shown on Figure 4) are generally similar to those for the Antarctic temperature proxy from the EDC3 or AICC2012 Antarctic ice core age models; however, we note that all these age models are partly dependent on astronomical tuning, so that the exact phasing between astronomical parameters and peak interglacial should be treated cautiously.
The structure of the latitudinal and seasonal distribution of insolation depends upon the precessional parameter and the obliquity [Berger, 1978;Berger and Loutre, 1991]. When precession minima (i.e., NH summer occurs at perihelion (when Earth is closest to the Sun in its elliptical orbit)) and obliquity maxima are in phase, they strengthen the NH insolation during boreal summer, especially at high latitudes. According to Milanković theory (Milankovitch [1941] and Milanković [1998] (translation)), boreal summer insolation controls the growth and shrinkage of ice sheets, so one might expect the interglacial minimum in ice volume to occur sometime after the insolation maximum. This timing issue will be discussed further in the sections about glacial terminations and inceptions. Yin and Berger [2010] noted that the δ 18 O minima of many interglacials occur while precession is rising (NH summer moving from perihelion toward aphelion), and in six cases (of 11) the isotopic minimum occurs about one fourth of a precessional cycle after the precessional minimum (NH summer at perihelion) ( Figure 4). This timing would be consistent with the expected response time of the melting ice sheets to the peak forcing.
However, such a relationship should be treated cautiously as it is partly related to the astronomical tuning of the LR04 stack, and the phase relationship of the second control on LR04 (i.e., deepwater temperature) is not directly addressed by Milanković theory. Furthermore, there are also several interglacials where this timing is not seen. The minimum in δ 18 O also always occurs within or just later than a quarter of an obliquity cycle after the obliquity maximum, with the notable exception of MIS 7c and 7e.
It is not obvious that the astronomical factors controlling the timing of interglacials should also control their strength. Of the two factors governing the benthic δ 18 O signal, ice volume (and hence boreal insolation) may have a limited effect on interglacial strength, if most NH ice disappears in every interglacial [Elderfield et al., 2012]. The insolation control on deepwater temperature is not clear. It is therefore unsurprising that (middle) Precession parameter (red, negative plotted high) with eccentricity as an envelope (blue); obliquity (in degrees).
(bottom) calculated insolation at 65°N [Berger and Loutre, 1992], for the summer solstice. Dashed vertical lines show the LR04 isotopic minima referred to in the text.

Reviews of Geophysics
10.1002/2015RG000482 ( Figure 4) there is no strict general relationship valid for all the interglacials of the last 800 ka between the strength of the interglacials and either the strength of the nearest precessional peak or the phase of obliquity maximum and precession minimum .
The latitudinal-seasonal distribution of daily insolation for a single time slice is usually used in model snapshot simulations. Compared to the average of recent interglacials, MIS 1, MIS 11c, and MIS 19c are characterized by less (more) insolation over the Earth during boreal (austral) summer because of weak precession caused by low eccentricity. In contrast, MIS 5e, MIS 7c, MIS 15a, MIS 15e, and MIS 17c are characterized by more (less) insolation over the Earth during boreal (austral) summer due to their high eccentricity. Due to its low obliquity, MIS 7e shows low insolation over the summer hemispheres with a deep minimum over the poles, and higher than average insolation over the winter hemispheres with a maximum in the midlatitudes. MIS 9e has high obliquity and a unique insolation pattern with high insolation centered over the poles in both summer hemispheres. For MIS 13a, the precession peak (506 ka) occurs so long before the δ 18 O peak (491 ka) that it does not make sense to consider it in the same way .
A full analysis must consider the insolation distribution as a function of time, as illustrated in Figure 5 for MIS 5e for an example. This interglacial shows the common features of many others [Yin and Berger, 2015] (the exceptions are those where the interglacial peak does not occur during rising precession (MIS 13a, 15e, and 17c) and falling obliquity (MIS 7c)). For the eight interglacials with rising precession index, a positive anomaly in insolation at the NH summer solstice ( Figure 5a) occurs before the interglacial peaks (based on LR04), covers the whole Earth and lasts more than 10 ka. It is followed by a negative anomaly a few thousands of years later. The amplitude of the negative and positive anomalies varies with interglacial and is controlled by the amplitude of precession variation (which in turn is determined by eccentricity), and at  [Berger, 1978]. The y axis is latitude in Figures 5a, 5c, and 5d and month in Figure 5b; the x axis is age in ka.

10.1002/2015RG000482
high latitude also by obliquity. The time evolution of the mean summer irradiance is very similar to that of the insolation at the summer solstice except that its maximum occurs 1 or 2 ka later and is shifted toward lower latitudes (Figure 5c). This underlines the problem raised by the choice of the insolation parameter that is used to "tune" or to compare with proxy records, the timing being different from one parameter to the other by a non-negligible amount.
A seasonal maximum in insolation at 65°N is seen before the interglacial peaks in LR04 (Figure 5b), shifting from the month of April to July over a period of more than 10 ka. The same pattern occurs, with some differences in timing, for all the interglacials with rising precession. After the peaks of the interglacials (in LR04), the maximum insolation occurs in boreal fall and winter, and the minimum in boreal spring and summer.
The values in Figures 5a-5c are controlled mainly by precession, while the total annual irradiation (Figure 5d) is only a function of obliquity . As a result of the falling obliquity at the interglacial peak, the time evolution of total annual irradiation exhibits also a more or less similar pattern for all interglacials, except MIS 7c. The interglacial peak (in LR04) ( Figure 5d) is preceded by a positive anomaly in the high latitudes of both hemispheres centered over the poles, lasting more than 10 ka, and followed by a negative anomaly. The differences between the interglacials and today and between the interglacials themselves are mainly in the high latitudes and are very small between 60°N and 60°S. Therefore, the high latitudes are considered critical regions for explaining the differences between the interglacials. An interesting feature is that the 10 ka long positive anomalies of the high latitudes occurring before the interglacial peaks tend on average to be larger since the mid-Brunhes (MB) than before, related to the larger amplitude of obliquity after the MB associated with the 1.3 Ma periodicity in the amplitude modulation of obliquity [Berger et al., 1998a]. To what degree this difference in obliquity and therefore in the total energy received at high latitudes explains any systematic differences between the pre-MB and post-MB interglacials needs further investigation.
In addition to the insolation distributions discussed above, other combinations of insolation factors might also deserve attention for the determination of the differences between the interglacials. The mean annual insolation, and its gradient, might be important for moisture transport and ocean physics [Loutre et al., 2004], while the latitudinal insolation gradient may induce a systematic difference between the pre-MB and the post-MB interglacials in the Southern Ocean ventilation and deep-sea temperature [Yin, 2013].

Ice Sheets and Greenhouse Gases
In general, and as discussed earlier, our definition assumes that the Northern Hemisphere ice sheets were similar to today for at least part of each interglacial, i.e., with significant ice only in Greenland. Most studies have also assumed that the Antarctic Ice Sheet was at a similar size to today in each previous interglacial. Thus, changes in ice sheet forcing are generally ignored when comparing interglacials. However, given the high sea levels [Dutton et al., 2015] that have been inferred for MIS 5e [Kopp et al., 2009] and MIS 11 [Raymo and Mitrovica, 2012], a significant reduction in these interglacials of either or both of the Greenland Ice Sheet Overpeck et al., 2006;de Vernal and Hillaire-Marcel, 2008] or West Antarctic Ice Sheet [Scherer et al., 1998] cannot be discounted and might impose a significant additional forcing to regional climate [Holden et al., 2011].
The GHG concentrations reached during different interglacials show considerable variation ( Figure 6). In particular, the interglacials between 450 and 800 ka ago have significantly lower concentrations of both CO 2 and CH 4 than do the later interglacials. The radiative forcing from greenhouse gases between an interglacial with relatively high GHG concentrations (280 ppm, 700 ppb) and the ones with lower values is around 1 W m À2 . This radiative forcing is sufficient, all other things being equal, to drive a global temperature drop of just under 1°C (for typical values of climate sensitivity [IPCC, 2013]).

The Diversity and Structure of Interglacials
The following sections will consider key aspects of interglacials: their strength, timing, shape, variability, and length. However, it is first important to have an overview of the collection of interglacials occurring in the last 800 ka, and of their anatomy and characteristics, in order to introduce the terminology we use later.
To illustrate the terminology, we use the last interglacial (MIS 5e, identified with the Eemian in NW Europe and the Sangamonian in North America) as an example (Figure 7). The figure shows benthic δ 18 O, relative Reviews of Geophysics 10.1002/2015RG000482 sea level as estimated by the "Red Sea method," and Antarctic deuterium (temperature proxy, representing an example of a Southern Hemisphere (SH) signal) and methane concentration (representing a signal based in the NH).
Moving forward through time, we first observe the end of the preceding glacial (MIS 6). In this case, as in many but not all other glacials, the glacial maximum (lowest sea level/coldest temperatures) was reached just before the termination. The glacial termination is characterized by several thousand years of rising sea level and temperature. A brief early maximum is seen in some parameters, notably Antarctic deuterium (representing temperature), and the transient peak is sometimes referred to as an "overshoot." Caution is needed in the use of this latter term, as it has been used to refer to brief centennial-scale rises of greenhouse gas concentrations above the interglacial value that followed, and to the much longer millennial-scale peak in temperatures seen for Antarctic δD at the start of the interglacial in Figure 7. We will discuss the likely origin of this early maximum in the section on terminations. However since it clearly provides a particularly warm climate at least in Antarctica, we consider it a part of the local interglacial when assessing maximum amplitudes. Methane clearly displays a different pattern of change, and this reflects the fact that maximum values for a given interglacial are achieved at different times in each record. After the initial maximum, a period of relatively stable interglacial climate is seen in MIS 5e (but not in all interglacials) and then values start to descend (generally more slowly than they rise during the termination): the first period of descent is the glacial inception. We also note that full glacial conditions are achieved only slowly, with cooling continuing well beyond the period shown in Figure 7. In this case, the interstadials MIS 5c and 5a intervene before strong glacial conditions arrive only in MIS 4,~74 ka ago. Figure 6. Concentrations of atmospheric greenhouse gases over the last 800 ka, bassed on measurements in Antarctic ice cores. Data for CO 2 [Lüthi et al., 2008] as recently revised and updated [Bereiter et al., 2015] and CH 4 [Loulergue et al., 2008], both plotted on the AICC2012 age model. Horizontal red lines are drawn at values considered typical of the preindustrial atmosphere.  [Jouzel et al., 2007] (representing Antarctic temperature, blue) and CH 4 [Loulergue et al., 2008] (red) on the AICC2012 age scale. Note that although there are not large age offsets, this plot should not be used to determine phasing between parameters as no attempt at synchronizing the three age scales used has been made.

Reviews of Geophysics
10.1002/2015RG000482 Figure 8 shows, for each interglacial, a range of key climate data sets. Each data set is plotted on its own age scale (i.e., there is no attempt at synchronizing), so spurious phasings may result. Each is presented with the same scales on each x and y axis, so that it is immediately possible to recognize long interglacials as opposed to short ones, and ones that are strong in a particular property. In each figure, the four panels relate to ( Figure  Antarctic air temperature (EPICA Dome C)).

The Intensity of Interglacials
Numerous records from the marine and terrestrial (including ice core) realm show a pattern in which the basic glacialinterglacial sequence is seen and in which peaks and troughs can easily be aligned (but not necessarily synchronized) with marine isotope stages (MIS) identified in marine benthic isotope records, such as the LR04 stack [Lisiecki and Raymo, 2005]. This indicates that, however, we define them, interglacial states are global phenomena, expressed in almost every type Figure 8. Presentation of key data sets for each interglacial discussed in this paper. In each case after this first expanded exemplar, the x axis (time) and y axis for each interglacial have approximately the same scaling to allow comparison between interglacials. Each plot runs (right to left) from the preceding glacial period through the termination, across the interglacial peak, and into the inception of the following glacial period (except MIS 1 which stops at the present).
(a) Obliquity and precession [Laskar et al., 2004]. (b) The LR04 benthic oxygen isotope stack [Lisiecki and Raymo, 2005], and the sea level estimate based on separating deepwater temperature and water isotopic content from a record at site ODP1123 (South Pacific) [Elderfield et al., 2012].
(c) The ice core record of CO 2 [Lüthi et al., 2008;Bereiter et al., 2015] and CH 4 [Loulergue et al., 2008]. (d) The temperature estimated from δD at EPICA Dome C, Antarctica [Jouzel et al., 2007], and SSTs at the North Atlantic site ODP982 [Lawrence et al., 2009]. This first set of panels is for MIS 1.  [Cheng et al., 2012] expresses strong precessional variability, but the interglacials (as we have defined them) do not stand out as unusual.) However, the intensity and precise timing of that expression varies between records and the different spatial and environmental patterns seen in each of them are important evidence to test against our understanding of the mechanisms behind them. In this section, we focus on the intensity or strength of interglacials.

Reviews of Geophysics
In some cases we can relate the proxy record we use with a single climate parameter (for example, the ice core measurements of CO 2 or Mg/Ca in marine records representing near surface temperature); in others (for example, the biosilica content of sediments in Lake Baikal) the association is more complex but still the glacial-interglacial pattern is seen. In each record one can define the intensity of the interglacial based on the extent to which the measured parameter changes in the interglacial direction (higher CO 2 , more arboreal pollen, less dust, etc.): this allows us to use a range of proxies without necessarily having a clear or unique interpretation of its meaning. In assessing the intensity of each interglacial, we therefore have a range of measures representing different aspects of the Earth system and different geographical regions.
A further issue to consider is whether one defines the intensity based on the maximum value achieved, or on the average over the entire interglacial or over a specified period. There are advantages and disadvantages to each approach: 1. Comparing the identical time period in each record would allow us to build the kind of snapshot climate reconstructions that are useful for comparison with most climate model outputs. However, the accuracy of the chronology in many records is variable, and our ability to synchronize records in any of the earlier interglacials is too low to allow such an approach in practice. 2. An average over the entire interglacial would probably be the most intuitive approach, allowing us to assess the integrated strength of an interglacial without assuming that the strength occurred simultaneously at every location. However, the period of averaging, and hence the result, is very dependent on the definitions used to determine the timing of the start and end of an interglacial, as discussed in section 7.5. In addition, this method carries an underlying assumption that there is a stable interglacial period over which it is possible to define stationary statistics, whereas in some parameters and interglacials no such period exists: rather, values approach a peak and immediately start to descend toward glacial values again. 3. Using the maximum value obtained during a particular interglacial puts focus on what may be a short anomalous period caused by local factors. This approach also suffers from the fact that we have records of varying resolution available, implying that we could find a maximum that is resolved in one record and miss it because it is unresolved in another. We must also be very aware that the maximum value may occur at different dates in different records, so that comparisons must not be used incorrectly. Nonetheless, it is very easy to define the maximum value objectively, it represents the value that the human eye is drawn to and is of environmental importance because it is the most intense condition that other aspects of the system (including ecological ones) have had to endure. 4. Finally, we could assess the amplitude of change over the termination. This has the advantage that it can be defined statistically and pinned to a particular time period in which the forcing may be determined. It shares with the maximum value criterion the disadvantage that termination may not be simultaneous in all records. Additionally, it is sensitive to the magnitude of the preceding glaciation and places emphasis on only one part of the interglacial, reducing its relevance for aspects such as sea level and ecological impact.

Methodology
In this paper we will follow and build on the approach and compilation used in earlier work [Lang and Wolff, 2011]. To ensure that we are comparing the same marine isotopic substage (but with no expectation that we have synchronized records to within a few kiloannum), marine records are all aligned to the LR04 age model [Lisiecki and Raymo, 2005]; ice core records are aligned to the AICC2012 age scale [Bazin et al., 2013]; terrestrial records are on their own age scale, which often implicitly uses a similar tuning target as LR04. We then take the maximum (or minimum, as appropriate) value that occurs within the interglacial (i.e., for practical reasons we use method 3 above): the only exception is that where resolution is better than 1000 years, we take the maximum value after applying a 1000 year smoothing, to avoid the influence of noisy outliers. More complex methods were attempted in previous work [Lang and Wolff, 2011], but as they gave similar results, we have decided to take the simplest approach here. It is important to note that the resulting indices or maps of interglacial intensity do not represent snapshots of a single time slice, since the maximum value does not always occur at the same date in different records. In some cases we will note where maxima occur particularly and noticeably asynchronously in different records, but will leave a more general discussion of the trends during interglacials to section 7.2.
Data have been included only if they meet the following criteria: 1. Records must cover the entire 800 ka period to allow a consistent comparison to be made between interglacials. We have relaxed this criterion for two 700 ka marine records that supply information in a region (western South Pacific) that is otherwise poorly represented. 2. Records must be continuous, with no hiatuses or data gaps, at least during interglacials. 3. The resolution must be better than 3 ka. However, in order to include a reasonable geographical range of SSTs, we relaxed this rule to just over 4 ka in some cases for such records. Reviews of Geophysics 10.1002/2015RG000482

Data Compilation
There are a few differences (Figures 9 and 10) from a recent compilation [Lang and Wolff, 2011], on which we have built.
First, a few additional marine records of SST have been included, allowing better, geographical coverage. There are now 17 SST records (of which two fall short of the full 800 ka period). Antarctic temperature is still represented by only one record (EPICA Dome C): although a second long ice core record, covering 720 ka from Dome F, distant from Dome C on the East Antarctic plateau, exists, the data are as yet unpublished. However, we note that the shorter published record [Kawamura et al., 2007;Sime et al., 2009] suggests that the relative intensity of different interglacials at Dome F is similar (though not identical) to that at Dome C and that this likely represents an Antarctic-wide pattern.
A second difference is that in the present discussion we have more strongly emphasized temperature (and especially SST) records; the earlier paper [Lang and Wolff, 2011] devoted a lot of attention to benthic and planktonic marine isotope records. While that study did show some minor differences in pattern between sites, it was clear that the common overprint (due to ice volume and for benthic records, deepwater temperature) on such records is sufficiently strong that they basically tell a single story. Here we therefore only display the LR04 stack but now supplemented with a new record that splits the isotopic signal at one site into a deepwater temperature (based on Mg/Ca) and seawater δ 18 O (ice volume) component [Elderfield et al., 2012].
In the terrestrial realm, the previous compilation [Lang and Wolff, 2011] included stacked loess records from China, biogenic silica from Lake Baikal and arboreal tree pollen from Tenaghi Philippon. Here we have added data for the lacustrine sediments of Funza in South America and for Lake El'gygytgyn in the Russian Arctic. For Lake El'gygytgyn [Melles et al., 2012], temperature data have been deduced from pollen for some periods in the record but are not yet available through 800 ka, so we are forced to use the Si/Ti data as an indicator of interglacial strength without attaching any particular climatic interpretation to them. For Funza, a new age model [Torres et al., 2013] has facilitated use of the data, but it remains difficult to know which data provide a climate signal. While a relationship is observed between arboreal pollen fraction, tree line altitude and temperature [Hooghiemstra and Ran, 1994], total arboreal pollen percentages are certainly biased by the fact that Quercus arrived in the region only at the mid-Brunhes. It is not obvious how to compare data before and after this arrival: we have calculated arboreal pollen percentage after removing Quercus but urge caution in interpreting the resulting record, particularly across the MB.
The pollen record from Lake Van (Turkey) offers a chance to discuss conditions in south western Asia but only extends to MIS 15b. We therefore discuss it but do not include it in Figure 9. For the Chinese loess record, the earlier compilation [Lang and Wolff, 2011] used magnetic susceptibility and mass accumulation rate: both parameters showed a pattern of glacial-interglacial variability that became indistinct in parts of the last 800 ka, and it is not clear that they gave a meaningful indication of interglacial intensity. It has been argued that indices of chemical weathering may be more useful, being controlled largely by the summer monsoon related moisture and temperature [Guo et al., 2009]. Here we have therefore used the ratio of free to total Fe 2 O 3 (Fed/Fet, representing weathering) and the redness index (believed to reflect soil temperature) at Xifeng in northern China [Guo et al., 2009]. However, we have additionally used a new frequencydependent record of magnetic susceptibility [Hao et al., 2012] that does show a glacialinterglacial pattern and that is interpreted as a record of East Asian summer monsoon strength. The same study also presents a grain size parameter that is interpreted as a record of winter monsoon strength. We refer to this record but do not include it in our tabular compilation because the minimum values of the parameter (indicating weakest winter monsoon) generally occur after the end of the interglacial (as judged in other records) and are therefore not directly indicative of interglacial amplitude.
A further long, continuous record is available in data from lake sediments at Heqing Basin, China, where an Indian summer monsoon index was compiled . This record shows glacial-interglacial variability in the part of the record from 900 to 350 ka, but little pattern after that, and does not include MIS 5e or MIS 1. We therefore discuss it but have not included it in our table.
A final continuous and well-resolved 800 ka record is available from the sapropel record obtained in the eastern Mediterranean, at sites ODP 967 and 968 [Ziegler et al., 2010;Konijnendijk et al., 2014]. Following the authors, we interpret the Ti/Al ratio as an indicator for monsoon strength in central to north Africa, and we therefore treat it as a terrestrial record. Like many speleothem records, it is dominated by precessional rather than glacial-interglacial variability.
Finally, there is a small change in methodology. In the earlier work, the maximum immediately after the glacial termination was used as the intensity value. Here we simply take the highest (or lowest) value that occurs during the interglacial, making no judgment about which part of the interglacial it occurs in.
6.3. Findings-Large-Scale Measures (Ice Volume, Deep Ocean Temperature, and CO 2 ) We start by considering some measures (Figure 9) that have a wide significance as an indicator, forcing and/or feedback at global scale: benthic oxygen isotopes, deepwater temperature (using site ODP1123, east of New Zealand), seawater oxygen isotope content (which is equivalent to global ice volume or sea level), and atmospheric CO 2 .
In the LR04 benthic oxygen isotope stack, MIS 5e, 11c, 1, and 9e (in that order) are considerably stronger than the remaining interglacials. MIS 7e and 7c are of comparable strength to all the interglacials that occurred before 450 ka ago, with MIS 17 as the weakest. Although there are small variations in rank order, this general pattern is confirmed in each of the numerous benthic and planktonic isotope records that were examined earlier [Lang and Wolff, 2011]: MIS 5e and 11c are consistently the strongest interglacials, with MIS 1 and 9e just behind, and MIS 17 and 13 are most commonly the weakest interglacial in the last 800 ka. Note that in the earlier compilation the values for MIS 13 were taken from MIS 13c, making it stand out as a particularly cold stage; in the new compilation, the inclusion of stronger values in MIS 13a, which is considered here as part of the continuing MIS 13 interglacial, means that it is not always the weakest.
Although the benthic isotope record from ODP Site 1123 has very similar intensities to that of the LR04 stack, the pattern becomes more complex when the benthic isotopes are split into their two components, deepwater temperature and seawater δ 18 O [Elderfield et al., 2012]. For seawater δ 18 O [Elderfield et al., 2012], MIS 11 stands out as having the lowest value (implying lowest ice volume), with MIS 9e and, surprisingly, MIS 17 next in rank. Most other interglacials (whether before or after 450 ka) show similar minimum values, comparable to those of the Holocene. One interpretation is that this is a threshold effect: all interglacials were warm enough to lose most NH ice volume outside Greenland, and higher values (for example, in MIS 11) imply further loss from Antarctica and/or Greenland. However, the imprecision of the ice volume calculation (using Mg/Ca values to correct the benthic δ 18 O) means that the latter inference is not robust; this is illustrated by the relatively low inferred sea level for MIS 5e, for which abundant other evidence suggests a level higher than today [Thompson and Goldstein, 2006;Dutton and Lambeck, 2012;Kopp et al., 2013].
Other attempts to partition benthic δ 18 O, such that sea level/ice volume can be derived from it [Bintanja et al., 2005;Siddall et al., 2010], used only a single measurement series to infer two outputs and therefore had to make assumptions that were not necessary in the work of Elderfield et al. [2012]. However, the latter suffers from the fact that the joint precision of the derived seawater δ 18 O is poor (±20 m for single data points). To test how robust the derived interglacials strengths are, we can look at independent estimates of sea level, albeit not covering 800 ka continuously. Continuous measurements over the last 500 ka using the "Red Sea" method [Rohling et al., 2009] suggest that MIS 5e has the highest sea level in that period; MIS 9e does not show high values, and even MIS 11c is not unusual compared to the other interglacials. Similar results emerge from globally distributed studies of foraminifera δ 18 O [McManus et al., 1999;Hodell et al., 2000;King and Howard, 2000;McManus et al., 2003]. Estimates using coral terraces and other indicators certainly suggest that sea level was higher than today in MIS 5e [Kopp et al., 2009;Dutton and Lambeck, 2012]. Some less direct evidence suggests the same for MIS 11 [Raymo and Mitrovica, 2012], although there continues to be more debate about sea level in this earlier interglacial [e.g., Hearty et al., 1999;Olson and Hearty, 2009;Bowen, 2010;Raymo and Mitrovica, 2012]. One way to explain particularly low ice volume/high sea level during interglacials would be through melting of all or part of the Greenland ice sheet. There is only indirect evidence about this: measurements of pollen in a marine core off south Greenland [de Vernal and Hillaire-Marcel, 2008] suggest a much reduced ice sheet particularly in MIS 11 and, more surprisingly, in MIS 13a. The result for MIS 11 is supported by geochemical provenancing of silt deposited south of the Greenland margin, suggesting near complete deglaciation of southern Greenland [Reyes et al., 2014]. Relatively, pollen concentrations in MIS 5e were only slightly raised, and other interglacials showed low values similar to those of today. The MIS 5e result is also supported by the geochemistry of silt, which in this case suggests that only limited retreat  [Colville et al., 2011]. In summary, ice volume was not very different in each interglacial of the last 800 ka, and there is no particular evidence for a change at the MB. However, there appears to be relatively strong evidence for somewhat higher sea level and lower ice volume during MIS 5e and some evidence suggesting the same for MIS 11. The low ice volume deduced for MIS 9 [Elderfield et al., 2012] is not yet supported by other studies, while the reduced Greenland ice sheet inferred in MIS 13 [de Vernal and Hillaire-Marcel, 2008] does not appear to have led to high sea level or low overall ice volume.
For deepwater temperature, where there are values deduced directly from Mg/Ca ratios at a South Pacific site [Elderfield et al., 2012] MIS 11c, 7e, 5e, 1, and 9e have significantly warmer values than MIS 7c or any of the interglacials before 450 ka ago. We note that a different ranking would be derived for the only other site where bottom water temperature has been measured, DSDP 607 [Sosdian and Rosenthal, 2009] in the North Atlantic. As well as representing a different water mass, the latter study used epifaunal rather than infaunal benthic foraminifera, and the reliability of deepwater temperature estimates in this case has been questioned [Yu and Broecker, 2010] and defended [Sosdian and Rosenthal, 2010].
The same interglacials emerge as the highest in atmospheric CO 2 , although the exact rank order is different, with MIS 9e showing the highest value at nearly 300 ppm. All of MIS 1, 5e, 7e, 9e, and 11c exceed 275 ppm, while none of MIS 7c, 13, 15a, 15e, or 17c exceed 260 ppm. MIS 19c stands out as an interglacial in the earlier part of the record that has a moderately high concentration in the newly corrected compilation [Bereiter et al., 2015], while MIS 17c stands out in the case of CO 2 for its particularly low values, reaching a maximum of 244 ppm.

Findings-SST and Surface Air Temperatures
We now consider temperature estimates from around the world (Figures 10 and 11), considering SSTs derived from Mg/Ca, alkenones or foraminifer assemblages from marine sediments, and surface air temperatures (SAT) derived from water isotopes for Antarctic ice. A pattern that seems consistent across most of the world is of particular warmth during MIS 5e and to a lesser extent MIS 11c, and of particularly cold conditions during MIS 13. There are apparent exceptions to this first-order rule; for example, at ODP1123 in the SW Pacific, MIS 5e is relatively cool and MIS 13a not unusually so. Nonetheless high temperatures are observed, e.g., for MIS 5e in Antarctica (EDC), the South Atlantic (ODP1090), SE Pacific (GeoB3388-1 and PS75/034-2), Eastern equatorial Pacific (ODP846), western equatorial Pacific (ODP806B, ODP1146, and MD97-2140), NE Pacific (ODP1020), and the North Atlantic (ODP982 and ODP607). Only at site ODP552 (which is close to ODP982 in the Atlantic), and at site ODP1123 mentioned above, does MIS 5e rank lower. MIS 5e is also an exceptionally warm interglacial in records from regions that are not represented in our compilation because they do not meet our criteria: for example, MIS 5e has the warmest SSTs of any interglacial of the last 600 ka [Martrat et al., 2007;Rodrigues et al., 2011] in records from the Portuguese Margin on the east of the North Atlantic, while MIS 11c, 5e, and 9e stand out as considerably warmer than other interglacials [Caley et al., 2011] in a record from the east coast of South Africa (MD96-2048, Indian Ocean) that lacks resolution for a quantitative analysis.

10.1002/2015RG000482
A similar global intensity seems to be observed for MIS 11c, but there are some notable exceptions in the table as well as in SST records that do not reach back 800 ka [e.g., McManus et al., 1999;Martrat et al., 2007] and in some isotopic records [Helmke and Bauch, 2003]. The only exceptions to the cold conditions of MIS 13 in the continuous records we used occur at a couple of tropical sites where the differences between interglacials are anyway small. However we note that some shorter SST records show a relatively warm MIS 13 [e.g., McManus et al., 1999], and some areas of the world have little or no data.
Other interglacials show warmth at only some sites. MIS 9e is the warmest interglacial of all in Antarctica, mirroring the high CO 2 concentration. Similar warmth is also seen at some other sites (such as site ODP1123 in the SW Pacific, site ODP1146 in the eastern equatorial Pacific, and site ODP982 in the north Atlantic). However, there is no consistent geographical pattern to warmth in this interglacial and because this interglacial has an exceptionally short peak, we must exercise caution in comparing records with low resolution. MIS 7e, though warm in Antarctica and at site ODP1123 in the SW Pacific, is otherwise cool, comparable to MIS 15a and 15e, 17c, and 19c. At least in the Northern Hemisphere, MIS 7c is generally warmer than MIS 7e.

Findings-Terrestrial Records
There are few terrestrial records we can use in this synthesis. Ice core CH 4 may be viewed as a signal integrating parts of the NH terrestrial biosphere. It shows a pattern in which MIS 9e is strongest, while MIS 11c, 19c, 1, and 5e stand significantly (more than 25 ppb higher concentration) stronger than MIS 7e, 15e, and 17c. MIS 13a again emerges as the weakest interglacial (Figure 6a).
For arboreal pollen at Tenaghi Philippon in Greece [Tzedakis et al., 2006], MIS 11c and 1 are strongest and 17c and 19c are weakest. Going eastward, arboreal pollen at Lake Van (Turkey, not included in Figure 9) is not especially prominent in MIS 11, and instead MIS 5e (followed by MIS 9 and MIS 7c) appears particularly strong. For biogenic silica in Lake Baikal [Prokopenko et al., 2006], MIS 11c and 15a are the strongest interglacials and MIS 13a and 7e the weakest. For Si/Ti at Lake El'gygytgyn in the Russian Arctic [Melles et al., 2012], MIS 11c is hugely stronger than any other interglacial. We counsel caution in interpreting the Funza [Torres et al., 2013] (arboreal pollen minus Quercus) data because of the changing structure of the tree population but have presented the data for completeness.
With the Xifeng loess indices used here, MIS 13a, 11c, and 5e are particularly intense [Guo et al., 2009]. The unusual strength of MIS 13 in the East Asian monsoon region is supported by the magnetic susceptibility stack (summer monsoon proxy, strong in interglacials, strongest in MIS 13) [Hao et al., 2012]; the winter monsoon proxy (not shown in the table; weakest winter monsoon in interglacials and inceptions) also reaches its most extreme (weak winter monsoon) values during MIS 13. The intensity of MIS 13a in these records contrasts with its weakness in the methane record, also assumed to have a strong East Asian control.
The Indian monsoon index reported from lake sediments at Heqing  is hard to interpret in terms of interglacial strength. The Indian summer monsoon index was, at face value, strongest in MIS 15e, and very weak in all interglacials since MIS 10. However, the authors of this study interpreted the variability in terms of a combination of Northern and Southern Hemisphere controls that lead to a complex pattern that is not expected to relate to other climate indices in any simple way. Interestingly, the authors also present  a winter temperature proxy (percent Tsuga) for the site; in this proxy, MIS 11c is the warmest interglacial, and MIS 13a, MIS 7e, and MIS 7c are the coolest, more in line with the pattern seen in marine temperatures.
The Mediterranean sapropel record (Ti/Al at ODP sites 967/968) suggests that rainfall in the African monsoonal belt was particularly intense in early MIS 13a, and particularly weak in MIS 19c. However, we emphasize that the record is dominated by precession and also shows intense monsoonal peaks in glacial periods, so this record cannot really be considered to be providing a measure of interglacial intensity.

Synthesis of Interglacial Intensity
The main conclusions do not differ from those of an earlier synthesis. Across the whole range of records, there is a tendency for more intense interglacials in the period after 450 ka compared to that from 800 to 450 ka. However, as others have emphasized [Candy et al., 2010;Candy and McClymont, 2013] this is by no means a general rule: MIS 7c and in most cases 7e is often of similar strength to earlier interglacials and each of MIS 13a, 15a and 15e, 17c, and 19c show strength comparable to the later interglacials in one or more important records (MIS 13a shows a strong loess weathering signal, MIS 15a and 15e are particularly strong at Baikal, MIS 17c is particularly strong in seawater δ 18 O and in two North Atlantic SST records, and MIS 19c shows high CH 4 values). For the midlatitude Atlantic, it has been argued, using a range of shorter SST records along with British terrestrial evidence that there is no significant shift even in mean interglacial strength at 450 ka [Candy and McClymont, 2013], though more data are needed to confirm whether this is a true regional finding.
In any case, it can be agreed that there is not an exclusive step change for most properties at the MB (mid-Brunhes), as MIS 7e (and certainly MIS 7c) show a widespread weakness that is more characteristic of the pre-MB events. If we exclude MIS 7c, then a step change to persistently stronger values after 450 ka can be inferred only for CO 2 , Antarctic temperature and deepwater temperature. As an example, for CO 2 , the new interglacial range after the MB (279-299 ppm), excluding MIS 7 ac, is about 30 ppm higher than that seen before it (244-269 ppm). The general strength of MIS 5e and the widespread weakness of MIS 13a are of almost global nature, extending to temperatures nearly everywhere. A very strong MIS 5e and very weak MIS 13a are also observed in the first attempt at a global SST stack [Shakun et al., 2015]. The one characteristic in which MIS 13 shows unusual strength is the summer monsoon recorded in Chinese loess. The patchy nature of warmth in MIS 9e may arise in part from the short length of this interglacial compared to the resolution of many records. Its overwhelming strength in CO 2 , CH 4 , Antarctic temperature, and in the low δ 18 O seawater values suggest that this interglacial may warrant further detailed study.
A final point is that MIS 15a and 7c, which we newly introduce as interglacials, stand comparison in all respects with their peers. MIS 15a and 15e are of comparable strength in most records. MIS 7c is weaker than 7e in many but not all records, but still as strong as MIS 13a.

Understanding Interglacial Intensity
The mechanisms for the observed pattern of interglacial intensities can be explored in more detail by model simulations. Conceptual models are able to some extent to simulate the sea level amplitude of interglacials [Parrenin and Paillard, 2003], but more physically based models must be used to understand the mechanisms. However, few transient simulations designed to explore this issue have not yet been carried out, except for MIS 5e (discussed later), so that studies to date have generally simply used peak forcing in snapshot experiments. The individual contributions of insolation and CO 2 to the interglacial climates of the past 800 ka were quantified by Yin and Berger [2012] through snapshot simulations with the model LOVECLIM (and did not include MIS 7a-7c and 15e). They carried out simulations using astronomical parameters at the times when NH summer occurs at perihelion just preceding the interglacial peaks of the LR04 stack, combined with the maximum GHG concentrations occurring around the isotopic peak. Their simulations show that, globally, MIS 9e, 5e, and 11c are indeed expected to be the warmest interglacials, MIS 1 and 19c are very similar and slightly above the average of the last nine interglacials, MIS 7e and 15a should be cool interglacials, and MIS 13a and 17c the coolest. This is in broad agreement with the qualitative conclusions of the data. However, the differences in their intensity result from different contributions of insolation and CO 2 . According to the model , the insolation and CO 2 of MIS 5e and 9e both contribute to a warming, but MIS 11 is a warm interglacial only because of its high CO 2 concentration. The insolation and CO 2 of MIS 13a and 17c both contribute to a cooling. MIS 7e is a cool interglacial because its insolation-induced cooling beats its CO 2 -induced warming, and MIS 15a is due to the reverse situation.
The relative contributions of insolation and CO 2 to the modeled climate are different from one region to another and from one climatic variable to another . The modeled anomalies of the annual mean temperatures averaged over the globe and over the southern high latitudes are mainly controlled by CO 2 and perhaps obliquity, whereas over the northern high latitudes, regional insolation plays a dominant role. In low latitudes, insolation (and precession) definitely dominates changes in Reviews of Geophysics 10.1002/2015RG000482 monsoon precipitation and therefore in vegetation. For example, due to high eccentricity, MIS 15a is the simulated wettest and most vegetated interglacial at least at low latitudes, followed by MIS 5e and MIS 17. The relative roles of eccentricity, precession, and obliquity on temperature, precipitation, sea ice, and vegetation at different latitudes have been assessed and quantified using snapshot  and transient simulations [Yin and Berger, 2015] and is shown to vary between climatic variables and latitudes. The varying importance of insolation, CO 2 , and the three astronomical parameters among variables and latitudes may explain at least partly the regional and temporal diversity of the interglacials.
In addition to EMIC studies, a similar intercomparison [Herold et al., 2012] of the climate response to the combined effect of insolation and CO 2 for five warm interglacials (MIS 1, 5e, 9e, 11c, and 19c) used the ESM CCSM3. Their results show that the greatest variation between these interglacials occurs in sea ice margins and in the regions where the insolation difference is the largest. At northern middle latitudes, large insolation variations and extensive land masses lead to large temperature differences between the interglacials during boreal summer. However, the mean annual temperature exhibits relatively little change over most of the continents but shows large variation in far North Pacific, Arctic, and around Wilkes Land. These model results reveal the complex interplay between insolation and CO 2 in different regions and for different variables. However, when comparing them with proxy reconstructions, one should keep in mind that 1. Ice sheets were fixed in these simulations, and therefore, the possible impact of ice sheet changes were not taken into account. Although the difference in minimum ice sheet size between the interglacials might be small according to the sea water δ 18 O [Elderfield et al., 2012], early peaks in some locations (e.g., the early maximum in Antarctic temperature seen in Figure 7) might still be strongly influenced by both ice sheets and freshwater inputs. 2. All the experiments described are snapshot simulations; therefore, the possible impact of interglacial duration and deglacial history on intensity is not taken into account, and they do not represent the peak insolation forcing ( Figure 5) at particular latitudes. 3. The climate response to insolation change has strong seasonal, as well as regional, dependence. It is therefore essential to consider the seasonality of the signals recorded in proxies before comparing them with modeling results.
As far as the MB is concerned, one modeling study  shows that an MB clearly appears in the modeled variables dominated by CO 2 (such as the global annual mean temperature and the southern high-latitude temperature), acting particularly in boreal winter and the SH . There is no clear MB in modeled precipitation, vegetation, and the northern high-latitude temperature which are all dominated by insolation . A related study [Yin, 2013] suggests that insolation alone can induce stronger southern westerly winds and weaker stratification during several of the pre-MB Reviews of Geophysics 10.1002/2015RG000482 interglacials (MIS 13a, 15a, and 17c but also 7e). These lead to stronger Southern Ocean ventilation and a cooler deep ocean. A cooler pre-MB deep ocean is observed in the data [Elderfield et al., 2012]; the model simulations [Yin, 2013] suggest that this cooling results from a combination of insolation and lower CO 2 concentration. This work implies that the apparent shift at 450 ka resulted from a series of individual interglacial responses to various combinations of insolation conditions. To the extent that there is a change in behavior at 450 ka, it is essential mainly to understand the change in peak CO 2 concentrations.
While these model experiments are a promising start, the one interglacial (MIS 5e) for which multiple experiments are available suggests that we should interpret the results with caution. Within the Paleoclimate Model Intercomparison Project (PMIP3), time slice climate simulations have been run by 13 modeling groups with a hierarchy of climate models (AOGCMs and EMICs, different resolutions) forced with the astronomical and GHG changes of 130-125 ka .
These simulations can be compared to global quantitative data syntheses of maximum annual surface temperatures for MIS 5e [Turney and Jones, 2010;McKay et al., 2011]. A caveat is that these data syntheses implicitly treat the warmest phases as globally synchronous [Cortese et al., 2013], and this emphasizes the need to reconstruct individual time slices within the interglacial [Capron et al., 2014]. Much higher mean annual temperatures than preindustrial are reconstructed for Europe, the Arctic coastal regions of Alaska and Siberia, Greenland, and the Atlantic Ocean north of 40°N. The multimodel mean simulates the reconstructed pattern of NH annual warming [Bakker et al., 2013]. The magnitude of the warming over the NH high-latitude continental regions, though, is underestimated in the multimodel mean. This could be partly due to the lack of interactive vegetation in many of the models: vegetation has been shown to have a positive feedback with surface temperature, associated with changes in land surface albedo and plant transpiration [Foley et al., 1994;Swann et al., 2010]. Strong warming at proxy sites near the Arctic Ocean also suggests a role for sea ice [Schurgers et al., 2007;Fischer and Jungclaus, 2010;Otto-Bliesner et al., 2013]. The East Antarctic ice cores record much higher annual temperatures than present early in MIS 5e, which are not present in simulations. Sensitivity simulations that account for changes in the East Antarctic ice sheet elevation [Bradley et al., 2013] [Capron et al., 2014] also emphasizes that the Northern Hemisphere temperature signal is very different at 130 ka and 125 ka in the data (in contrast to the models), which again might be partly explained by a bipolar seesaw response. These results emphasize that we cannot expect good agreement in the magnitude of warming at particular times and sites in earlier interglacials without accounting for the transient evolution of the forcings and range of internal feedbacks. The MIS 5e work also points to the need for reconstructions of other parameters such as precipitation and sea ice to further constrain model performance. As an example, recent work [Kleinen et al., 2014] comparing climate and vegetation against the few available terrestrial records for time slices in MIS 11, shows general agreement between models and reconstructions, although reconstructed precipitation changes are often larger than those simulated by models.

Glacial Intensity
In previous work [Lang and Wolff, 2011], comparative glacial intensities (looking at extreme values occurring in the period just before termination) were also surveyed and we will not repeat this analysis here. It has been suggested [Berger and Wefer, 2003;Lang and Wolff, 2011] that, particularly in marine isotope records (including LR04), the strength of interglacial response partly mirrors the strength of the preceding glacial maximum, already implicated in determining the timing of terminations [Raymo, 1997]. This is an interesting concept that might perhaps be related to ideas that will be developed later, that the ice volume during the glacial maximum acts as a threshold for termination to occur at all. However, while it is an intriguing idea, it has one prominent failure (in that the very strong glacial MIS 16 is followed by one of the weakest interglacials (15e)). While it is valuable to consider that the size of the glacial maximum ice sheet likely determines the volume and rate of freshwater forcing, and therefore has an influence on the shape and regional pattern of the following interglacial, there does not appear to be a straightforward predictive rule.

The Structure and Timing of Interglacials
Having gained an overview of the interglacials of the last 800 ka, we now dissect aspects of them in more detail, following a chronological sequence. First, we discuss the way that the onset of interglacials (which equates Figure 10. Maximum values for each interglacial for temperature records.

10.1002/2015RG000482
effectively to glacial termination) occurs, up to the point at which full interglacial conditions exist. We then examine the trends and variability of climate during full interglacial conditions. After that we review the end of interglacials, which roughly equates to glacial inception. In each section, considerations of timing will be important-what determines when an interglacial starts and when it ends? We will bring all this together in a separate section about the duration of interglacials.

The Onset of Interglacials (Glacial Terminations) 7.1.1. Interglacial Onsets in the Paleoclimate Record
The onset of each interglacial (Figure 8) is characterized by a change, relatively rapid in the context of the length of a glacial cycle, in all the properties associated with an interglacial: rises in sea level, associated falls in marine oxygen isotope ratios, increased SSTs and air temperatures (including those inferred from Antarctic ice cores), increases in CO 2 and methane concentrations, and so on. Glacial cycles of the 100 ka world are often referred to as having a sawtooth pattern, reflecting the rather sharp onset of interglacials compared to the slow descent back to full glacial conditions.
Obviously interglacial onset coincides approximately with glacial termination: a glacial "Termination" is strictly defined as the midpoint of rapid transition between late Pleistocene glacial and interglacial climate conditions, as recorded in marine stable oxygen isotopes [Broecker and Van Donk, 1970]. Terminations therefore represent key stratigraphic boundaries in the Pleistocene MIS stratigraphy [Shackleton and Opdyke, 1973]. However, the glacial terminations assigned simply Roman numerals are only those associated with the onset of odd-numbered MIS, and so there was traditionally no numbered termination associated with interglacials at MIS 7c or MIS 15a. Nonetheless, there is a sharp and rapid jump in the marine oxygen isotope record at the MIS 7d/7c and 15b/15a transitions, of a magnitude larger than that of Termination VI (MIS 14/13), and accompanied by changes in all the usual properties ( Figure 8). This was partly recognized by previous authors [Huybers and Wunsch, 2005;Cheng et al., 2009] who numbered the onset of MIS 7c as T3a or T-IIIA, distinguishing it from the onset of MIS 7e (T3b or just T-III). The period of rapid isotopic change (greater than 0.1‰/ka) in the LR04 stack [Lisiecki and Raymo, 2005] typically spans 4-9 ka, sometimes including an interval with a lower rate of change. We note that 4-9 ka is also the typical transition duration for Antarctic ice core isotope (representing temperature) records [Röthlisberger et al., 2008]. Figure 11. Maps showing peak interglacial SSTs and EDC temperature for MIS 5e, 7e, 11c, and 13a. For each location the colors follow those in Figure 10, i.e., the site has a dark red dot in the warmest interglacial, dark blue in the coldest.

10.1002/2015RG000482
On timescales that are longer than the mixing time of the global ocean (e.g., several thousand years), the MIS stratigraphy can be interpreted to reflect glacioeustatic sea level (i.e., total land ice volume) changes [Shackleton, 1967], and terminations therefore reflect periods of relatively rapid deglacial glacioeustatic sea level rise, marking the onset of an interglacial sea level high stand. However, on shorter (centennial to millennial) timescales, terminations are in fact diachronous climatostratigraphic boundaries: when exactly the marine oxygen isotope Termination happens will depend on where it is observed in the ocean, due to the finite (and potentially variable) mixing time of the ocean, resulting in local temporal offsets between MIS boundaries of up to~4 millennia [Skinner and Shackleton, 2005;Lisiecki and Raymo, 2009;Primeau and Deleersnijder, 2009].
It is important to note that "Terminations" (defined from marine oxygen isotopes) are not identical to, nor synchronous with, other markers of the transition between glacial and interglacial climates (in the remainder of this section we will, however, refer to terminations meaning the transition in all the properties that change at the glacial-interglacial transition). A classic illustration of this is provided by the distinction between the onset of MIS 5e (defined in marine oxygen isotopes) and of the Eemian interglacial (defined on the basis of terrestrial vegetation change) [Shackleton, 1969;Shackleton et al., 2003;Sánchez Goñi et al., 2005]. A particularly unambiguous example is provided by the trace gas records measured in Antarctic ice cores (Figure 8c), and therefore, on the same chronology; CO 2 typically jumps by 50-100 ppm over a period of several thousand years, more or less synchronously with rising Antarctic temperature [Pedro et al., 2012;Parrenin et al., 2013], while CH 4 makes its transition to interglacial conditions within a few decades toward the end of the CO 2 increase (and at least in the case of Termination I, synchronously with Greenland temperature [Rosen et al., 2014]). Another example comes from comparing changes in SST and benthic δ 18 O in the same (South Atlantic) core, where the benthic signal lags at TI and TV by several ka [Vazquez Riveiros et al., 2010]. In general, however, it is challenging to synchronize records from different locations and different archives beyond the radiocarbon era, so it is only for Termination I that a clear sequence of events can be defined [Shakun et al., 2012].
Milanković theory would suggest that terminations (and associated interglacial onset) should occur on the rising limb of NH summer insolation. As not every precessional cycle leads to an interglacial, there must be another factor that leads to interglacials occurring in some precessional cycles and not others [Imbrie and Imbrie, 1980]. This will be discussed further below. However, we can assess the timing of the transition to an interglacial, relative to an absolute astronomical cycle, in some parameters. It is important to avoid circularity, and therefore, records that have been tuned using an assumption about the temporal control of the astronomical signal on the record being considered. Independent constraints on the timing of interglacial onset are provided (a) in ice core records, using ice core O 2 /N 2 ratios which, for reasons unrelated to those controlling the climate signals, are assumed to reflect local insolation [Kawamura et al., 2007] and (b) in corals and speleothems by uranium series dating [Wang et al., 2001;Thompson and Goldstein, 2005] of rapid changes in oxygen isotope content which are assumed to mirror rapid changes in other proxies such as ice core methane. These confirm the broad association between terminations and insolation changes that has been inferred from the similar frequency structure of climate and insolation time series [Shackleton, 2000]. In particular, within dating uncertainties of the independent methods, the last four interglacial onsets in properties such as Antarctic temperature, CO 2 , and marine oxygen isotopes coincide with the period during which NH daily summer insolation is rising [Kawamura et al., 2007;Cheng et al., 2009], and with changes in any other correlative aspect of the global insolation field [Huybers and Denton, 2008]. Because of its faster pacing, precession tends to dominate in discussions about timing, but it is important to note that high obliquity should also be conducive to glacial termination [Drysdale et al., 2009], and in most cases termination indeed occurs during a time of high obliquity [Huybers and Wunsch, 2005]. It remains an open question, whether the precise timing of a termination, relative to its insolation forcing, may influence the character of the ensuing interglacial [Tzedakis et al., 2012b]. Furthermore, it is unlikely that the climate system responds only to insolation at a specific latitude and season.
In seeking to determine how interglacials arise (i.e., how deglaciations are triggered) and more specifically whether or not differences between terminations lead to differences between interglacials, it is probably most useful to consider separately the deglacial evolution (precise timing and character) of various key global climate parameters, such as: ice volume, regional surface temperatures, low-latitude precipitation/monsoon

10.1002/2015RG000482
systems, greenhouse gases, and the marine carbonate system. All of these parameters show commonalities between successive terminations, but key differences also emerge, most notably in terms of their amplitude and/or duration.
One important commonality between terminations (of the last 800 ka) is that they appear consistently to include millennial-scale events with a particular regional fingerprint. This fingerprint is that of the so-called bipolar seesaw [Broecker, 1998;Schmittner et al., 2003;Stocker and Johnsen, 2003]. High-latitude SH (Antarctic) temperatures begin gradually to rise prior to a more abrupt rise in high-latitude NH (Greenland and North Atlantic) temperatures, after which Antarctic temperatures gradually fall again ( Figure 12). The canonical theory for this "asymmetry" between Southern and Northern Hemisphere high-latitude temperatures is that it reflects the impact of a "collapse" or at least significant reduction of the Atlantic overturning circulation on cross-equatorial meridional heat transport [Broecker, 1998;McManus et al., 1999;Schmittner et al., 2003;Stocker and Johnsen, 2003]. Thus, a significant freshwater discharge into the North Atlantic at the onset of deglaciation would maintain cold conditions in the northern high latitudes while contributing to warming in the southern high latitudes. Additional global impacts would arise from these "terminal seesaw events," including perturbations of the low-latitude atmospheric convection [Cheng et al., 2009] and the marine carbon cycle [Anderson et al., 2009;Skinner et al., 2013]. A number of papers have noted evidence of an increase in Agulhas leakage (involving transfer of warm, salty water from the Indian to Atlantic Ocean) during each glacial termination [Peeters et al., 2004;Caley et al., 2012]. While this has been proposed as a possible cause of the resumption of strong Atlantic meridional overturning circulation (AMOC) during interglacials, other studies would imply that a reduction in meltwater forcing is alone sufficient to explain the resumption [Liu et al., 2009]. Further work is needed to assess the interactions and cause-effect relationship, between AMOC strength, frontal movements, and Agulhas leakage.
The occurrence of terminal seesaw events raises important issues for the definition and study of interglacial onset. It has been suggested that terminations may simply start as a millennial southern warming (bipolar seesaw tilted south, probably induced by northern meltwater) that through the operation of feedbacks runs out of control until interglacial onset becomes inevitable  and that terminations are intrinsically associated with millennial-scale events [Sima et al., 2004;Broecker et al., 2010]. In such a scenario, while the conditions for interglacial onset are under astronomical control, the precise (millennial-scale) timing of termination should be variable, as the actual termination is overprinted by the millennial-scale events. Such a model also implies that there may be a continuum from small Antarctic Isotopic Maxima [EPICA Community Members, 2006] through to full terminations, with some events having an appearance of "aborted terminations." At a practical level, the impacts of the bipolar seesaw on high-latitude temperatures and on global GHG concentrations mean that peak values of these parameters might be overprinted and dominated by an "independent" millennial climate dynamic or anomaly or may even fall outside of the interglacial interval defined by ice volume [Tzedakis et al., 2012a]. Thus, many important peaks in the Antarctic ice core records of temperature and of GHG concentrations (especially the early interglacial maxima of MIS 5e, 7e, and 9e, for example) must be treated with care when attempting to link peak climatic conditions during different interglacials to peak radiative forcings, for example. This underlines the more general point that Antarctic temperature trends, especially on millennial but perhaps also on longer timescales, are driven by local to regional processes and can differ significantly from global average surface temperature trends [Shakun et al., 2012].

Simulating and Explaining Interglacial Onset
From a climate dynamics perspective, glacial terminations are a manifestation of the strongly nonlinear nature of the Earth system [MacAyeal, 1979;Paillard, 2001], and of the ways in which a seasonal/regional forcing can lead, though feedbacks, to a long-term/global response. While the rate of sea level rise during the middle-late Pleistocene scales broadly with the concurrent intensity of NH insolation forcing, terminations are characterized by anomalously high rates of sea level rise ( Figure 13). The exact presentation of this finding in Figure 13 is determined by the scaling and smoothing we have used. However, it is interesting that the high rate of sea level rise (at least based on the benthic isotope record) at the onset of MIS 7c and 15a can be seen mainly as a response to a very strong insolation forcing (these periods have the strongest forcings in the record). In contrast, at most traditionally numbered terminations, high rates of change are observed only because they stand well above the scaled forcing (implying high levels of nonlinearity). This suggests Reviews of Geophysics 10.1002/2015RG000482 that these two interglacials, each occurring after a very short cold period, may be the result of a response to a very strong insolation forcing, rather than a nonlinear response to a weaker forcing.
From an energy balance perspective, two key amplifiers of the deglaciation process, both internal to the climate system as opposed to the external forcing, are the atmospheric greenhouse gas concentration (that of CO 2 in particular) and land albedo [e.g., Köhler et al., 2010]. The latter will have contributed to glacial terminations via an ice-albedo feedback on temperature. In principle, the trigger for this feedback could stem partly from a marine ice sheet instability mechanism [Weertman, 1974;Schoof, 2007], which would lead to the rapid collapse of a part of an ice sheet and enhanced ice discharge into the ocean. Recent explanations have emphasized the observation that termination generally only occurs when a large ice volume (or low sea level) has been reached, and this was also the basis for successful conceptual models [Parrenin and Paillard, 2003]. This underlies, for example, the proposal that an ice-lithosphere feedback, leading to a nonlinear tendency toward ice sheet instability when the ice sheet is large, could be at play [Abe-Ouchi et al., 2013]. It has also been proposed that intensive deposition of aeolian dust over the NH ice sheets can lead to lowering of snow and ice albedo and thus significantly enhance surface melt beyond a threshold in dust deposition [Peltier and Marshall, 1995;Ganopolski et al., 2010], with strong dust deposition apparently occurring only when ice sheets are large. Climate-ice sheet models of different complexity are able to simulate the timing of glacial terminations successfully with prescribed astronomical variations and time-dependent CO 2 concentration Bonelli et al., 2009] or even with a constant, sufficiently low, CO 2 concentration [Ganopolski and Calov, 2011;Abe-Ouchi et al., 2013]. The idea that large ice volume is required for termination to occur seems to be breached for the onset of MIS 15a, 13, and 7c. However, for 15a and 7c, as discussed earlier, very strong insolation forcing seems to be overcoming such a requirement; one could think of that forcing as aborting the glacial that would otherwise be expected, returning the system to an interglacial state. The onset of MIS 13 remains enigmatic, but still the paradigm that large ice volume induces nonlinear feedbacks seems strong. However, these studies leave open the question of exactly which ice sheet instability mechanism has been the most important for terminations during the last 800 ka.  [Jouzel et al., 2007], Antarctica (grey), and methane concentration [Loulergue et al., 2008] (blue), whose fast changes track those of Greenland temperature. The high-pass-filtered records (filtered to retain frequencies from 1 to 10 ka) clearly show the importance of a hemispherically different millennial-scale variability in each termination.

10.1002/2015RG000482
Explanations of the greenhouse gas amplifier also remain equivocal, though it would appear that the primary source of CO 2 to the atmosphere during deglaciation must have been the deep ocean [Sigman et al., 2010] and that the primary locations for ocean-atmosphere CO 2 exchange were likely in the Southern Ocean [Anderson et al., 2009;Skinner et al., 2013] and/or the North Pacific [Galbraith et al., 2007]. Changes in biological export productivity, ocean carbonate chemistry and the overturning circulation (and physical properties) of the ocean will have all played important roles in mediating deglacial CO 2 rise and interglacial CO 2 levels [Kohfeld and Ridgwell, 2009]. Some papers also suggest that the CO 2 plays an active role in triggering deglaciation, with low sea level inducing changes in southern circulation that releases carbon that has accumulated in the deep ocean [Paillard and Parrenin, 2004;Bouttes et al., 2010;Paillard, 2015]. However, no quantitative physical explanation has yet been advanced for the different CO 2 levels and trends that were experienced during each interglacial of the last 800 ka [Wolff et al., 2005]. The exact mechanism linking insolation forcing and deglacial CO 2 rise remains unclear.
The millennial-scale regional climate perturbations that occur across terminations (i.e., terminal seesaw events) could be part of the link between CO 2 rise and deglaciation. This would rely on their combined impacts on the marine carbon cycle (i.e., the CO 2 amplifier) and polar temperatures (i.e., the ice-albedo amplifier). Thus, a large freshwater flux into the North Atlantic caused by rapid (astronomically induced) melting of the NH ice sheets during glacial terminations would have strongly affected the Atlantic meridional overturning circulation (AMOC) [McManus et al., 2004], which in turn would have caused anomalous cooling in the Northern Hemisphere and anomalous warming in the Southern Hemisphere via the thermal "bipolar seesaw" mechanism [Broecker, 1998;Schmittner et al., 2003;Stocker and Johnsen, 2003;Shakun et al., 2012] (Figure 12). Although the first of these effects may have acted as a negative feedback on European deglaciation (via reduced air temperatures), it may also have helped to trigger the destabilization of parts of the North American ice sheets via subsurface warming in the North Atlantic, ice-shelf destabilization and enhanced ice discharge into the ocean [Alvarez-Solas et al., 2011]. This latter process (a positive feedback on freshwater forcing in the North Atlantic) could help to explain why particularly large ice-rafting events (Heinrich events) would be a persistent feature of glacial terminations [McManus et al., 1999;Venz et al., 1999;Hodell et al., 2008;Vazquez Riveiros et al., 2013].
In the SH, rapid warming during terminal seesaw events would produce an "overshoot" in Antarctic temperatures relative to the unperturbed equilibrium climate state [Ganopolski and Roche, 2009] (Figure 12). This transient "excess" Antarctic warmth may have contributed to a complete disintegration of the West Antarctic ice sheet during some interglacials [Marino et al., 2015] and could help to explain much higher than Holocene Antarctic temperatures during early phases of MIS 11c, 9e, and 5e, for example, [Holden et al., 2010]. Another important SH impact of terminal seesaw events would be to trigger a rapid pulse of CO 2 release from the Southern Ocean, for example, via their impact on winds [Anderson et al., 2009] or on sea ice/upper ocean stratification [Stephens and Keeling, 2000;Skinner et al., 2014]. This effect is consistent with existing atmospheric and marine radiocarbon records [Hughen et al., 2006;Skinner et al., 2010;Burke and Robinson, 2012;Skinner et al., 2013] and atmospheric stable carbon isotope records [Lourantou et al., 2010], and would be a key contributor to deglacial CO 2 rise. While numerical models Figure 13. Interglacial onsets over the last 800 ka (termination numbers TI to TIX mark traditionally numbered glacial terminations). (bottom) Benthic LR04 isotope record (blue) [Lisiecki and Raymo, 2005] and ice core CO 2 (red) [Lüthi et al., 2008;Bereiter et al., 2015]. (top) The rate of change of the benthic isotope ratio (light blue), overlying a plot of Northern Hemisphere insolation forcing (black). The rate (related partly to the rate of sea level change) broadly tracks insolation forcing throughout the last 800 ka but stands out above the scaled insolation curve at glacial terminations.

10.1002/2015RG000482
disagree on the exact mechanisms at play, many are able to produce a significant increase in atmospheric CO 2 following a shutdown of the AMOC [Marchal et al., 1998;Schmittner and Galbraith, 2008;Menviel et al., 2014], though they do not produce the very high rates and magnitudes of CO 2 rise observed at intervals during the last deglaciation [Marcott et al., 2014].
In summary, terminations clearly represent a strongly nonlinear response of the global climate system to regional changes in the seasonality of solar radiation. Although incompletely understood, this nonlinearity appears to emerge primarily from instabilities in the ice sheets (triggering an ice-albedo feedback) and instabilities in the carbon cycle (triggering an enhanced greenhouse effect). Furthermore, it would appear that transient millennial-scale regional climate anomalies could represent a crucial bridging mechanism between the regional influence of astronomical insolation forcing, and global climate and carbon cycle impacts. This would underline the possible role of short-term regional climate perturbations in longer-term global climate transitions, and would raise the question of whether or not the intensity, the duration, or the stability of an interglacial period might depend on the pattern of change during the termination that preceded it. This remains an important and unresolved question.

Trends and Variability During Interglacials 7.2.1. Intra-Interglacial Climate Variability
During glacial periods, the dominant variability is at millennial scale, documented over the last glacial period through the identification of Heinrich events [Heinrich, 1988] in marine sediment cores, and Dansgaard-Oeschger [Dansgaard et al., 1993] and Antarctic Isotopic Maximum events [EPICA Community Members, 2006] in ice cores. Long marine, terrestrial, and Antarctic ice core records have subsequently revealed the recurrence of such millennial-scale features in all glacial periods [McManus et al., 1999;Jouzel et al., 2007;Wang et al., 2008].
No comparable millennial events have been identified under warm climate conditions once the major phase of ice melt is complete, and interglacial periods have commonly been considered as less variable than glacial ones. For example, one of the most conspicuous features of Greenland ice cores is the presence of large, abrupt millennial-scale variability during the glacial period compared to the relatively subdued climate variability under full interglacial conditions. This pattern of enhanced suborbital climate variability during glacials and its suppression during full interglacial periods is a consistent pattern in North Atlantic sediment cores throughout the Pleistocene [McManus et al., 1999;Hodell et al., 2013;Barker et al., 2015;Hodell et al., 2015].
However, from the well-documented Holocene, it is clear that there is significant variability at centennialmillennial timescales [e.g., Wanner et al., 2008;Marcott et al., 2013], although the evidence that this variability is spatially coherent is less clear. In addition to such variability, each interglacial exhibits differences in multimillennial climatic trends controlling the observed "shape" of each past interglacial.
Exploring this intra-interglacial variability is important for placing Holocene climate into context, for detection and attribution of climate changes and for projecting a future where anthropogenic forcings will be added to natural ones and internal variability. In addition, the trends during interglacials are a crucial factor in understanding the relevance of the interglacial intensities discussed in section 6 and may play a role in the timing of the subsequent inception. Documenting intra-interglacial climate variability, however, remains challenging because of the limited temporal resolution of available records and the difficulty of comparing signals between proxies and between interglacials where relative and absolute chronologies may have large uncertainties.
Moreover, the response time specific to each climatic component may induce different delays in the imprint of climatic variations (for example, the time of ocean overturning being several centuries against a few months for tropospheric mixing). For the same reason combined with the spatial heterogeneity of climate forcings and feedbacks, local, or regional proxy records (e.g., temperature, SST, and pollen) have to be differentiated from globally integrated signals such as sea level or atmospheric greenhouse gas composition. The characteristic averaging time of proxy signals can also result in a loss of information by the smoothing of highest-frequency signals due to inherent processes (i.e., the bioturbation in sediments may smooth the records by several thousand years at some sites; the timescale of bubble enclosure in ice core sites with low snow accumulation rates may smooth GHG records by several centuries [e.g., Spahni et al., 2003]).

10.1002/2015RG000482
The description of interglacial variability is therefore inherently limited by the timescales resolved by proxy records and their spatial representativeness.

Millennial to Submillennial-Scale Variability During Past Interglacials
The climate variability of the last~12 ka has been well preserved in many paleoclimatic records and described at high temporal resolution, up to decadal scale (see Wanner et al. [2008] for a synthesis). Persistent significant millennial-scale to multidecadal-scale oscillations have emerged from those Holocene records and millennial to multicentennial events ( Figure 14)-first detected in the North Atlantic sector [O'Brien et al., 1995;Bond et al., 2001;Debret et al., 2007]-have been documented at a global scale in many records [Mayewski et al., 2004;Wanner et al., 2011]. The lack of paleoclimatic data at sufficient time resolution for past interglacial periods has long prevented such climatic features being identified within earlier interglacials. The emergence of climatic records resolved at decadal or better scale for MIS 5e Bigler et al., 2010;Pol et al., 2014] and MIS 11c [Koutsodendris et al., 2011;Pol et al., 2011] starts to allow the identification of such events in a few locations (Figure 14), thus confirming evidence from Europe [Martrat et al., 2004;Allen and Huntley, 2009], the Mediterranean [Martrat et al., 2004;Sprovieri et al., 2006;Mangili et al., 2007;Milner et al., 2013], and the North Atlantic [Bauch and Kandiano, 2007], without, however, providing a global framework. High-resolution records for earlier interglacials, such as MIS 19c, are starting to emerge [Ferretti et al., 2015], and the retrieval of high-resolution records from, for example, annually laminated and partially annually laminated lakes [Stockhecke et al., 2014] raises the prospect of future studies down to subcentennial scales.
Understanding the causes of intra-interglacial variability requires the knowledge of past natural forcings. Using available records over the last 12 ka, solar activity [Steinhilber et al., 2012] and volcanic eruptions [Castellano et al., 2005;Gao et al., 2008] have been investigated as potential drivers of the Holocene millennialscale to submillennial-scale climate oscillations. This has led to the conclusion that observed frequency patterns of the current interglacial period cannot be simply explained by external forcings [Wanner et al., 2008]. In parallel, model studies have demonstrated the possibility for the climate system to exhibit internal modes of variability at those key timescales Park and Latif, 2008], which not only add to the usual forced variability but can also act as amplifying factors [Jongma et al., 2007]. Holocene multicentennial to millennial variability has also been reported to result from fluctuations in ocean overturning in response to changes in freshwater fluxes [Renssen et al., , 2012Mathiot et al., 2013].
The lack of solar or volcanic forcing records for earlier interglacials obviously limits the exploration of drivers of intra-interglacial variability. Similarly, no systematic modeling study is yet available that allows investigation of the characteristics of simulated internal climate variability under different astronomical contexts. The combination of such modeling studies with additional high-resolution data from a range of sites and interglacials is needed to assess whether climate variability at millennial and submillennial timescales is itself dependent on the background climate. This was for instance recently suggested in the case of an East Antarctic ice core record, where multicentennial to millennial climate variability appeared stronger during the warm MIS 5e phase than during the Holocene [Pol et al., 2014].

Multimillennial Trends and Variability
In contrast with submillennial-scale variability, multimillennial climate trends can be observed in numerous records covering the last 800 ka. We use the paleoclimatic records of Figure 8 to draw general patterns for each interglacial, within the available spatial proxy coverage. Here we focus on just a few data sets. The marine benthic isotope record [Lisiecki and Raymo, 2005] serves as a proxy of global sea level changes along with deepwater temperature. The CH 4 [Loulergue et al., 2008] and CO 2 [Lüthi et al., 2008] records define changes in the GHG radiative forcing of the last 800 ka. Note that changes in CH 4 concentration also reflect the impacts of NH climate on methane emissions. Antarctic temperatures inferred from the EDC δD signal [Jouzel et al., 2007] have been related to global-scale temperature changes with a polar amplification factor of~2 at glacial-interglacial scale [Masson-Delmotte et al., 2010a], consistent with recent findings [Elderfield et al., 2012;Milker et al., 2013].
Based on those data and multimillennial timescales, interglacial periods can be classified into two different structures [e.g., Tzedakis et al., 2012aTzedakis et al., , 2012b]. The first group (MIS 5e, 7e, 9e, and 19) is characterized by highest values reached at the beginning of the interglacial (i.e., soon after glacial termination), particularly in the EDC δD and GHG records. The relatively short peak in these properties (Antarctic temperature, CO 2 , and CH 4 ) Reviews of Geophysics 10.1002/2015RG000482 is sometimes described as an overshoot (as discussed earlier) and could be a result of millennial-scale changes discussed in section 7.1. The marine isotopic/sea level maximum typically occurs soon after the peak in these properties and then most parameters slowly decline toward inception.
A second group (MIS 11c, 13a, and 17c) shows gradually increasing trends in many properties, followed by maximum values late in the interglacial intervals. In MIS 11c and 17c, only about half the inferred change in sea level (based on site ODP1123) occurs in the period of~5 ka of sharp sea level rise that would generally be thought of as the termination, with the remaining rise occurring over the following 15-20 ka [Rohling et al., 2010]. Other properties show a positive but weaker trend. MIS 13 requires a special comment as Termination VI, consisting of the sharpest change in marine isotope values (LR04) in this period of the record and the main jump in Antarctic temperature and CO 2 , occurs at about 530 ka. However, sea level is inferred to have reached its maximum only after 500 ka, separated from the termination by a distinct dip in sea level and other properties, and by almost two precessional cycles. This description is partly a result of our criteria that defined MIS 13 as a single interglacial, but the absence of a further sharp rise in sea level that could be interpreted as a termination does suggest that the pattern is indeed unusual. We note that there is no obvious insolation trend to explain these two classifications, and indeed, sea level apparently stabilizes or even continues to rise through an entire half-cycle of reducing NH summer insolation at MIS 17, 13, and 11. However, there may be systematic differences in the phasing of obliquity and precession maxima that should be further explored.
MIS 15a, 15e, and 7a-7c show either rather flat patterns without significant trends or different trends in different key parameters, making them difficult to classify. In some parameters (such as Antarctic temperature and CO 2 ), there is a distinct dip between MIS 7c and 7a, but this is by no means obvious in, for example, LR04.
As examples of local or regional nuances, we, however, note that the transition between interglacial and glacial states is not always obvious in some records, as observed between MIS 19 and 18 and between MIS 9 and 8 in local SST (in North Atlantic [Ruddiman et al., 1989] or Southern Ocean [Martinez-Garcia et al., 2009] sectors for instance) or pollen records [Tzedakis et al., 2006;Fletcher et al., 2013]. In the same SST records, the return to  Wanner et al. [2008]. Detected events with respect to methodology displayed in Figures 14c and 14f are labeled according to Pol et al. [2014]. (c, f) Detrended signals obtained by the subtraction of the red signals of Figures 14a and 14d from the black signals of Figures 14b and 14e. Black dashed lines represent the respective ± 1σ levels of the MIS 1 and 5e data distributions. Extreme isotopic values which do not range within the ± 1σ interval are highlighted by red (warm excursions) and blue points (cold excursions). Significant interglacial climatic subevents are taken as sequences of two consecutive opposite excursions.

MIS 1 Versus Past Interglacials
As shown in the EDC δD and CO 2 signal, MIS 1 shows a small optimum and might therefore belong to the first group of interglacials, together with MIS 5e, 7e, 9e, and 19c (we note that MIS 11c also shows a small early peak in these properties, but after a small decreasing trend, it eventually trended upward over a long period [Rohling et al., 2010]). Early Antarctic optima of MIS 5e and 1 have both been related to the bipolar seesaw in respect to Greenland or North Atlantic records [Masson-Delmotte et al., 2010b;Marino et al., 2015]. This supports the hypothesis that Antarctic early maxima during interglacial periods are responses to the NH deglaciation history.
However, analogies between MIS 1 and earlier interglacial periods have also been investigated on the basis of astronomical contexts. MIS 11c has long been considered as a partial analog due to its low eccentricity values (cf. section 4.1), implying weak precessional forcing [Droxler et al., 2003;Loutre and Berger, 2003;McManus et al., 2003]. MIS 19c also offers an astronomical configuration comparable to that of MIS 1, including (unlike MIS 11c) similar phasing between obliquity and precession [Pol et al., 2010;Tzedakis, 2010]. Although not often mentioned in this context, a number of other interglacials (MIS 5e, 9e, 15a, and 15e) also have this phasing , but with varying amounts of precessional power and obliquity amplitude.
These astronomical comparisons have assumed particular relevance because of the observation that CO 2 and CH 4 concentrations have been increasing for the last 8 ka and 5 ka, respectively. It has been observed that this is in contrast to decreases in similar astronomical contexts in other interglacials (especially MIS 19) and suggested that it is driven by land use changes due to deforestation and agriculture [Ruddiman, 2003[Ruddiman, , 2007. However, this "early anthropogenic" hypothesis has been challenged [e.g., Joos et al., 2004;Claussen et al., 2005;Broecker and Stocker, 2006;Elsig et al., 2009] (see also the special journal issue introduced by Ruddiman et al. [2011]). During deglaciations, atmospheric CO 2 usually rises by 40-80 ppm, as a consequence of carbon released by oceans into the atmosphere. The CO 2 shows a decreasing multimillennial trend during most interglacials (including MIS 19, often considered the best astronomical analog) after the initial posttermination peak. The observation of the increasing trend in MIS 1, coupled with new estimates of possible anthropogenic carbon release [Ruddiman, 2013], provides motivation for the early anthropogenic hypothesis. There are two clear exceptions to the decreasing trend apart from MIS 1 (20 ppm increase): MIS 11c (10 ppm increase) and MIS 15e (20 ppm), both discussed above as partial astronomical analogs in different senses. This emphasizes the point that the relationship between astronomical parameters and CO 2 trends is expected to be indirect and complex. The magnitude and isotopic composition of the MIS 1 CO 2 trend can be explained by purely natural mechanisms related to ocean carbonate chemistry and terrestrial carbon, including peats [Elsig et al., 2009;Kleinen et al., 2010;Schmitt et al., 2012]. Small differences between land and ocean processes can drive either positive or negative multimillennial atmospheric CO 2 trends, the overall result depending on astronomical context-which affects continental climate and associated land storage-and deglaciation history (including the pattern of bipolar seesaw changes).
For methane, a natural explanation for the Holocene trend has also been proposed on the basis of modeling experiments , while new measurements of the interpolar difference in CH 4 through the Holocene [Mitchell et al., 2013] suggest that both a natural southern source and a (probably anthropogenic) northern source have been increasing together. In summary, while natural mechanisms may be capable of explaining the MIS 1 observations, they are not yet diagnosing why MIS 1 behaves differently to other interglacials; as a result, the extent to which anthropogenic activities have influenced greenhouse gas concentrations over the last several millennia remains open.
Recent (instrumental, last century) changes have been compared to the broader context of MIS 1 temperature variations. Continental temperature reconstructions spanning the past 1500 years have highlighted the exceptional worldwide structure of the late twentieth century warming, compared to all previous multidecadal intervals in the period studied [IPCC, 2013;Pages 2 k Consortium, 2013]. Further work [Marcott et al., 2013] suggests that global temperatures in recent decades are approaching those of the early Holocene optimum; however, care must be taken when comparing decadally resolved temperatures with less well-resolved proxy data from earlier periods. The same study suggests clearly that projected year 2100 global temperatures [IPCC, 2013] will be warmer than those reconstructed during all of MIS 1 under all emissions scenarios considered.

Simulating Interglacial Dynamics
While the Mid-Holocene (6 ka) has long been a target for snapshot paleoclimate modeling intercomparisons [Braconnot et al., 2012] major efforts have been recently made (e.g., http://pmip3.lsce.ipsl.fr) to perform multimillennial transient simulations for MIS 1 (8-0 ka ago) and 5e (130-115 ka ago), which complement the MIS 5e time slice simulations discussed earlier . The transient experiments used nine climate models of different complexities. Those simulations were used to test model abilities in reproducing the astronomically driven interglacial temperature trends [Bakker et al., 2013[Bakker et al., , 2014. Consistent with the expected seasonal impact of climatic precession changes, simulated trends during the warmest months of the year appear systematically larger than annual mean temperature trends. Overall, simulated temperatures also show smaller variations for MIS 1 than for MIS 5e, as expected from the different rates of seasonal insolation changes in the studied intervals. While decreasing temperature trends appear robust within periods and models for middle to high latitudes of the Northern Hemisphere, no clear signals emerge in the Southern Hemisphere. The intermodel spread is moreover particularly large in monsoon areas. As regards MIS 5e, model results appear difficult to reconcile with the magnitude of temperature changes reconstructed from Greenland ice cores [NEEM Community Members, 2013] and do not capture the warming depicted in Antarctic ice core data [Masson-Delmotte et al., 2011].
Transient experiments covering at least one precessional cycle have also been made to simulate the climate of five warm interglacials (MIS 1,5e,9e,11c,and 19c) with the model LOVECLIM [Yin and Berger, 2015]. The authors show that the phase relationship between precession and obliquity is important in shaping the trends of regional temperature within each individual interglacial. Although other factors, such as the aftereffects of meltwater forcing may also be involved, this result might partly explain the differences in the internal structure between the interglacials as observed in proxy records discussed in section 7.2.3.
Specific sensitivity modeling studies have de facto emphasized the key role of Laurentide and possibly West Antarctic fresh water fluxes to correctly simulate the polar multimillennial early Holocene temperature variability patterns [Renssen et al., 2009;Mathiot et al., 2013]. Similarly, past changes in Antarctic topography have been suggested to explain the MIS 5e Antarctic warming [Holden et al., 2010;Otto-Bliesner et al., 2013]. Those indirect lines of evidence suggest that ice sheet-climate interactions play a key role on interglacial climate and its variability for climatic conditions 1-2°C warmer than during the preindustrial period. More generally, understanding the part of the intra-interglacial variability that is predictable and a response to radiative perturbations and the part that is linked to the internal behavior of ocean and atmosphere dynamics is a major challenge and requires additional paleoclimate information, simulations, and new modeling methods such as data assimilation [Goosse et al., 2006]. Assessing the magnitude of this internal intra-interglacial variability component over different timescales is particularly important in the context of ongoing global warming.

Glacial Inception
Each peak interglacial interval of the Pleistocene gave way to a subsequent glaciation. Due to the sawtooth character of the most recent ice age cycles [Broecker and Van Donk, 1970], the end of the peak interglacial conditions often occurred long before the glacial maximum of the same cycle. Yet the long-term trend following the peak interglacial within each large climate cycle of the last 0.8 Ma can be seen as leading toward increased glaciation as each glacial proceeds, despite significant precessional and millennial oscillations. For this reason, the discussion here will consider the end of the peak interglacial conditions as the initial stage of the subsequent glaciation and thus also as the glacial inception.
Following the Milanković theory (Milankovitch [1941] and Milanković [1998] (translation)) the succession of glacial and interglacial periods is immediately related to the amount of Northern Hemisphere summer insolation at latitudes where ice sheets would nucleate, because it controls the annual mass balance of snow. The direct effects of changes in insolation are amplified by the feedback caused by the higher albedo of snow and ice compared to land and sea. It is now understood that glacial inception is a dynamical process during which environmental changes (including changes in GHG concentrations) gradually affect climate in a way that further promotes glaciation [Berger, 1992;Imbrie et al., 1992].
The decline in northern summer insolation is the combined result of changes in the eccentricity-modulated precession and obliquity. Precession is particularly important for processes sensitive to maximum (solstice)
Along with insolation, existing data indicate that nearly every aspect of the Earth system undergoes changes during glacial inception. Ice begins to accumulate on land, causing global sea level to fall [Thompson and Goldstein, 2006]. The atmospheric concentration of CO 2 and other GHG declines [Barnola et al., 1987;Chappellaz et al., 1990;Petit et al., 1999;Lüthi et al., 2008;Bereiter et al., 2012]. Dust deposition increases [EPICA Community Members, 2004;Winckler et al., 2008], reflecting stronger winds, increased aridity, or both. Rain belts migrate and monsoons weaken [Yuan et al., 2004;Wang et al., 2008]. Millennial-scale variability starts to appear, most visibly around the North Atlantic Landais et al., 2006;Oppo et al., 2006;Bauch et al., 2011;Capron et al., 2012;Mokeddem et al., 2014]. Declining NH summer insolation reduces the amount of accumulated growing-season warmth, critical for the survival of boreal trees, and leads to expansion of tundra [Tzedakis, 2003]. Vegetation in temperate regions also changes, with, for example, latitudinal shifts in vegetation in Europe [Sánchez Goñi et al., 2005]. Loess deposition expands at the expense of soils [Kukla et al., 1990;Hao et al., 2012]. The expansion of ice on land and sea, the increase in dust loading in the atmosphere, loess over soils, and the shift in vegetation zones all contribute to an increase in the planetary albedo. Atmospheric, continental, and sea surface temperatures decrease. Although this cooling is neither globally synchronous nor homogeneous in magnitude [Ruddiman and McIntyre, 1984;Imbrie et al., 1992;Landais et al., 2006], it eventually represents a net reduction of the Earth's surface temperature [e.g., CLIMAP Project Members, 1984], which also extends into the deep ocean [Labeyrie et al., 1987;Shackleton, 1987;Martin et al., 2002;Elderfield et al., 2012].
Nearly all of these changes appear to accompany each glacial inception [Imbrie et al., 1992], although they do not necessarily occur in a fixed sequence or with consistent timing relative to insolation. This may be due to the nonstationary nature of glacial onset [Tzedakis et al., 2012b] or because each climatic component and location is influenced by a unique combination of local and indirect forcing effects of insolation, temperature gradients, seasonality, heat transport, etc. Some internal systems such as the monsoon [Yuan et al., 2004;Wang et al., 2008;Cheng et al., 2009;Govin et al., 2014] appear to respond rapidly to declining insolation. Other systems, such as changes in the strength and position of the westerly winds and the tropical rain belts, may be more likely to respond to ice sheet growth or hemispheric thermal gradients, as they do late in glaciation [e.g., Manabe and Broccoli, 1985;Peterson et al., 2000;Wang et al., 2004;Toggweiler et al., 2006;McGee et al., 2014], although their behavior during glacial inception is currently poorly documented. Important changes emanate from the high northern latitudes, where initial ice growth, vegetation shifts and albedo feedback are accompanied by dynamical changes in the ocean and atmosphere [Imbrie et al., 1992;Landais et al., 2006;Bauch et al., 2011;Capron et al., 2012;Mokeddem et al., 2014;Pol et al., 2014]. During the most recent and best studied glacial inception (from MIS 5e to MIS 5d), the early increase in sea ice export from the Arctic weakens and constricts the subpolar gyre circulation in the northwest Atlantic [Born et al., 2010], enhancing northward penetration of warm, saline subtropical waters toward the northeast Atlantic and into the Nordic Seas [Risebrobakken et al., 2007;Born et al., 2010Born et al., , 2011Mokeddem et al., 2014]. This in turn enhances convection and deepwater production in those seas and provides a negative feedback on sea ice production that contributes to a progressive millennial-scale oceanographic oscillation as insolation declines [Mokeddem et al., 2014]. The intermittent, yet persistent, production of deepwater maintains an overturning circulation throughout the glacial inception [Hall et al., 1998;McManus et al., 2002;Guihou et al., 2011] that contributes to extended warmth in the North Atlantic [Ruddiman and McIntyre, 1979;Ruddiman et al., 1980;McManus et al., 1994;Kukla et al., 1997;Shackleton et al., 2003]. This warmth persists despite, and may even enhance, ice sheet growth on surrounding continents [Ruddiman and McIntyre, 1979;McManus et al., 2002;Wang and Mysak, 2002]. Variability in the extent of southern deep waters occurs at the same time [Govin Reviews of Geophysics 10.1002/2015RG000482 et al., 2009, along with surface cooling [Ruddiman and McIntyre, 1984] and changes in the winds and sea ice around Antarctica [Toggweiler et al., 2006;Govin et al., 2009].
Northern ice sheet growth in this configuration appears to proceed rapidly, reaching the sea in the form of tidewater glaciers relatively early [McManus et al., 1994]. The ocean-cryosphere interaction of icebergs and meltwater then influences deep convection and the overturning circulation [McManus et al., 1999;Oppo et al., 2006;Mokeddem et al., 2014], initiating the bipolar seesaw of millennial-scale temperature changes that characterizes glaciations [Broecker, 1998;Landais et al., 2006;Barker et al., 2011;Capron et al., 2012]. A threshold for this response to initial ice growth is estimated to be 30-50 m of sea level equivalent [Chapman and Shackleton, 1999;McManus et al., 1999;Schulz et al., 1999] and may be used to define the end of the peak interglacial [Tzedakis et al., 2012b]. By this time, ice growth is sufficient and accompanied by reduced GHG concentrations and increased planetary albedo so that, despite subsequent increases in northern summer insolation, a return to full interglacial conditions is prevented, and glaciation generally proceeds progressively with only interspersed interstadial intervals until the conditions are eventually met for the next termination, as described earlier.
Reductions in CO 2 concentration clearly form an important part of the feedbacks that drive the system from an interglacial toward glacial maximum and prevent return to interglacial conditions. However, it may be less important for the onset of inception. In contrast to the situation in terminations, much of the CO 2 reduction in some inceptions occurs rather late compared to temperature decreases, at least in Antarctica. As an example [Schneider et al., 2013], in the inception following MIS 5e, CO 2 reduces from around 280 to 240 ppm only after 115 ka, when the inception reduction in Antarctic deuterium is already half completed. A similar phasing of change is seen in several other inceptions (see, e.g., end of MIS 11 in Figure 8).
The modern understanding of the glacial inception process also relies on advances in climate system theory, which have been made possible through the combination of conceptual modeling and experiments with numerical models of various levels of complexity. Different experiment designs are needed, both to quantify the effects of the astronomical forcing, possibly in combination with greenhouse gas concentration changes, on the snow mass balance [Royer et al., 1983;Harvey, 1988] and to identify the dynamics of the glacial inception over several thousands of years [Crowley et al., 1992;Gallee et al., 1992].
Consistent with Milanković theory, the effect of variations in summer insolation on snow and ice melt continues to be identified as at least one of the dominant controls on the annual mass balance of ice [Vettoretti and Peltier, 2011]. The localization of early nucleation sites is determined by a combination of orography, summer temperature and snow accumulation that point to the Hudson Bay and the Kewatin regions as the most likely candidates [Calov et al., 2005b;Otieno and Bromwich, 2009;Gregory et al., 2012]. At the hemispheric scale, summer temperature appears to be the most important driver of glacial inception and varies as a result of the combination of direct radiative effects amplified locally by the snow albedo feedback, along with significant remote effects associated with ocean cooling, sea ice expansion, and vegetation feedbacks [de Noblet et al., 1996;Crucifix and Loutre, 2002;Meissner et al., 2003;Kageyama et al., 2004;Wang et al., 2005]. The detailed evolution of individual ice sheets, however, may be largely determined by regional effects associated with changes in atmosphere and ocean circulation; hence, the importance of controlling as much as possible warm or cold biases in climate models for glacial inception studies [Vettoretti and Peltier, 2011]. Calov et al. [2005b] reported significant sensitivity of the Fennoscandian ice sheet on Norwegian Sea circulation patterns, involving significant variations in both snow accumulation and ablation. Models generally confirm the above-quoted interpretations of observations, according to which ocean circulation and sea ice dynamics significantly influence the regional expression and global propagation of glacial inception. For example, based on experiments with the Institut Pierre-Simon Laplace (IPSL) model, it was proposed [Born et al., 2010] that sea ice export toward the Canadian Archipelago, during the early stage of hemispheric cooling at the end of the Eemian, may disrupt Labrador Sea convection and favor Canadian ice inception, while Norwegian Sea convection is being maintained and delays the growth of the Fennoscandian ice sheet. As a consequence, the European glacial inception would occur later than the American inception, which, in contrast, may actually be accelerated by the evaporative moisture source of subtropical waters moving northward to balance the deepwater production in the Nordic Seas. Several models display continued or even enhanced North Atlantic deep circulation as glacial inception proceeds [Meissner and Gerdes, 2002;Wang and Mysak, 2002;Govin et al., 2012;Jochum et al., 2012].

10.1002/2015RG000482
While it is now clear that astronomical forcing affects the timing of glacial inception, a full understanding of the process must also explain why the climate system does not usually return to an interglacial state at the following insolation maximum. This is not covered by the original Milanković theory. A successful theory should also explain why the timing of glacial inception is not as tightly coupled to the northern summer insolation minimum as deglaciation is tied to the maximum [Tzedakis et al., 2012b].
Two conceptual frameworks may provide appropriate explanations. The first one is based on the notion of a "tipping point." It assumes that interglacial and glacial states are each dynamically quasi-stable. At some stage, the positive feedbacks induced by the early glacial inception process dominate the dynamics of the system, so that the latter runs toward ever stronger glaciation and reaches the glacial quasi-stable state. Mathematically, this hypothesis can be expressed by positing the existence of a bifurcation point in the climate system that is being crossed under the influence of the astronomical forcing. There are several elements supporting this scenario. First, ice sheet models consistently present a bifurcation structure that is compatible with this explanation. Oerlemans [1980], for example, showed the coexistence of at least two stable states in a 1-D ice sheet model, and this result was confirmed at least with a range of ice sheet atmosphere models [Calov et al., 2005a;Crucifix, 2011;Abe-Ouchi et al., 2013]. The tipping point model also gets some indirect support from the data. In particular, the bifurcation hypothesis is expected to leave an imprint on the dynamics of glaciation, which should be marked by an acceleration phase followed by a relaxation phase [Calov et al., 2009]. This is compatible with sea level observations over the last glacial inception. The approach of a bifurcation should also translate into specific dynamical phenomena, such as an increase in autocovariance (sometimes referred to as critical slow down) and an increase in variability [Dakos et al., 2008]. Increased variability, at least, is documented in Antarctic ice core records at the end of interglacial periods [Pol et al., 2011[Pol et al., , 2014 and may also be evident in some surface and deep ocean records Galaasen et al., 2014;Mokeddem et al., 2014].
In the second framework, the dynamics of the deglaciation prime the system for the following glacial inception, possibly with a time lag that depends on the amplitude, millennial patterns and rate of the preceding termination, creating the possibility of self-sustained glacial cycles for which the astronomical forcing is only a time control or a pacemaker. The idea was extensively developed by Saltzman and coworkers [Saltzman and Maasch, 1988;Maasch and Saltzman, 1990] and has recently been revisited [Crucifix, 2012]. Carbon cycle dynamics were stressed [Saltzman and Maasch, 1988] as an essential element for the generation of such a limit cycle. It has been suggested more recently that Southern Ocean dynamics controlling ocean carbon storage and release [Bouttes et al., 2010] or processes associated with the carbonate balance biogeochemistry [Rickaby et al., 2010] may indeed be involved, although the literature is not definitive about the importance of such mechanisms. We note in particular that in the CLIMBER model, effects associated with changes in sea surface temperatures and in the volume of bottom water of southern origin dominate biogeochemical effects during the glacial inception process [Brovkin et al., 2012].
This distinction between tipping point and limit cycle has consequences for understanding the factors that control the duration and the end of an interglacial. The tipping point model implies the existence of forcing thresholds beyond which inception occurs. Hence, the problem of identifying the controls on the duration of past interglacials or predicting the timing of the next glacial inception relates conceptually to the identification of this threshold [Petit et al., 1999;Berger and Loutre, 2002;McManus et al., 2003], even if, as already noted, the distinct effects of precession and obliquity on the distribution of insolation are such that defining a threshold on a specific insolation quantity will always be an oversimplification [Crucifix, 2011;Tzedakis et al., 2012b]. The self-sustained limit cycle model gives a clearer role to long-memory effects, so that the timing of each glacial inception is also conditioned in part by aspects of the preceding deglaciation, and may thus be intrinsically less predictable [Crucifix, 2013].

Spacing of Interglacials
While spectral analysis of Pleistocene climate records produces a significant spectral peak at a period centered on 100 ka [Hays et al., 1976;Lisiecki and Raymo, 2007], the averaging nature of this method is not appropriate for discussion of the spacing between interglacials (or terminations). In the LR04 record, the spacing of the interglacial peak value (as shown by dashed lines in Figure 4) is 84 ka , 86 ka (17c-15e), 35 ka (15e-15a), 84 ka (15a-13a), 86 ka (13a-11c), 76 ka (11c-9e), 90 ka (9e-7e), 22 ka (7e-7c), 94 ka (7c-5e), and 121 ka (5e-1). It is noteworthy that none of the intervals are within any reasonable uncertainty of 100 ka nor would they be so if we skipped MIS 15a or 7c (this would produce intervals of 119 ka Reviews of Geophysics 10.1002/2015RG000482 (15e-13a) and 116 ka (7e-5e)). It therefore appears that intervals cluster around 80-90 ka and perhaps 110-120 ka, as well as the shorter intervals within MIS 7 and 15. This could be related to the fact that the eccentricity period normally discussed as 100 ka is actually a group of periods of which 95 ka and 123 ka have the largest amplitude [Berger, 1977]. It also matches the possibility that the repeat times for glacial cycles may be related to multiples of obliquity or precession cycles [Raymo, 1997;Ridgwell et al., 1999;Huybers and Wunsch, 2005;Raymo and Huybers, 2008].

Duration of Interglacials
Knowledge of the duration of past interglacial periods provides valuable insights into the climate system dynamics, and is important from the perspective of understanding how natural processes, modified by human influence on atmospheric composition, will control the future evolution of the Holocene. It is also fundamental for comparing the interglacials and better understanding their diversity. As discussed previously, no single, continuous, integrative, and absolutely dated indicator exists for past interglacial intervals. Instead, a wide range of evidence exists for many climatic components and processes at local, regional, and global scales. The duration of interglacial conditions in each of these may be estimated, first, in order to consider how the timescales may differ between components of the climate system and second,, to evaluate the different duration of each past interglacial in the same record or records.
We emphasize that the question of whether an interglacial exists (discussed in section 2) may be a separate one from that of defining its length. In particular we proposed that an interglacial requires that sea level approaches modern levels (within ±20 m). However, it would then be extremely restrictive to say that the interglacial starts only when sea level reaches zero: taking MIS 5e as an example, Antarctic temperature and CO 2 reached values generally considered characteristic of interglacials well before such a level is reached.
The first question to be discussed is whether there is such a thing as an objective interglacial duration at all. In some conceptual models of glacial cycles [Paillard, 1998;Parrenin and Paillard, 2003], the system switches from a glacial to interglacial (or deglacial) state based on insolation and ice volume criteria, thus clearly defining the start and end of the interglacial. As already noted, defining such an interglacial state from the data is not straightforward.
We have already noted that the most recent glacial terminations seem to be synchronized (within uncertainty of a few kiloannum) to a decreasing precession index (i.e., rising northern summer insolation). On theoretical grounds, we would expect inception to occur during rising precession, i.e., approximately 0.5 (or 1.5, 2.5, …) precession cycles later. Obliquity may also play a significant role in determining how many precession cycles an interglacial includes [Tzedakis et al., 2012b]. Assessments of interglacial duration can be made in the context of this expectation.
Recognizing that there is no simple definition for the length of an interglacial, we calculate in Table 2 the length of each interglacial using a range of methods, and three different measurement sets, each represented by a precise data set and considered to have wide significance: the marine sediment benthic oxygen isotope stack (LR04) [Lisiecki and Raymo, 2005], atmospheric CO 2 concentration [Bereiter et al., 2015], and the Antarctic temperature proxy, deuterium [Jouzel et al., 2007]. The methods we have used in Table 2 are as follows: 1. Fixed predetermined thresholds. This method has the advantage of simplicity. However, thresholds based on Holocene values, as might be implied by the idea that interglacials are periods of comparable warmth to the present, inevitably fail for the weaker interglacials. As an example, the typically quoted post-MB interglacial CO 2 value of 280 ppm, if used as a threshold, would provide a duration for only three interglacials., We have therefore chosen less stringent subjective threshold values. Specifically, we used (a) LR04 δ 18 O value of 3.6, chosen because it is in the middle of the range of values used in our sensitivity study for the identification of interglacials; (b) CO 2 from Antarctic ice with a threshold at 260 ppm (the minimum observed in the Holocene); and (c) δD from the EPICA Dome C ice core with a threshold at À403‰, which was used previously because it is the minimum 300 year average for the Holocene. 2. Half-transition boundaries. In this method, the minimum value in the glacial maximum immediately before glacial termination, and the maximum value immediately following termination are identified. The average of these values is then the transition value determined independently for each interglacial. An interglacial is defined as starting when the measure last passes through the transition value,

10.1002/2015RG000482
and ending at the earliest point that the measure passes back down through the value. The difference between the two times is the duration. (Note that we chose this method over the variant that defines the inception half-transition as a different value bisecting the maximum and minimum value at the start and end of the inception itself.) This method has the advantage that it ensures that a value can be determined for each period we have defined as an interglacial and that it uses the rapid termination (which might be considered the point at which the system changes mode). However, where inception is slow, it produces long interglacials that are controlled by the speed of glaciation rather than the period of low ice volume or high temperature. The method was applied to the same three data sets as for the threshold method. 3. A third series of methods uses data sets that have been normalized so that they have a threshold that is adapted to different time periods with different amplitudes of variability. For one set of estimates, the threshold (variable on the raw data but fixed on the normalized records) was 1 standard deviation (SD) above the mean for each individual cycle (in this case the normalization was carried out across the periods of order 100 ka encompassing MIS 2-5, 6-7, etc., not the shorter time period representing, e.g., MIS 7d/7e). With this method the threshold can vary strongly from one interglacial to the next, and the interglacial length is adapted to the climate in which it is embedded. In a second test we calculated the threshold using a 400 ka interval around the relevant interglacial, meaning that the threshold varies slowly with time, and again set the threshold to 1 SD above the mean. These methods identify the warmest 17% of the record within a given window and ensure that there are still identified interglacials when the amplitude of cycles is smaller. In both these methods, the averaging interval and SD threshold are subjective choices. 4. A final method identifies the cessation and onset of millennial-scale variability believed to arise from freshwater input that can only occur when there are substantial ice sheets. As applied here, the start of the interglacial is at the terminal oscillation of the bipolar seesaw and the end of the interglacial is set to 3 ka before the first identified reactivation of the millennial switch during inception [Tzedakis et al., 2012b]. Advantages of this method are that it is process based, and it places the start of the Holocene at the end of the Younger Dryas, exactly where the official definition of the start of the Holocene has been placed [Walker et al., 2009]. However, it is not always simple to identify millennial activity, especially on inception. 7.5.1. Effects of Methodology and Measured Parameter A first examination of the resulting durations (Table 2 and Figure 15) highlights some of the issues with using different methods. As we have discussed, the ice age cycles of the first half of the period (~400-800 ka) were generally smaller in magnitude and did not reach the peak interglacial values of the last four cycles prior to the Holocene [Shackleton, 1987;EPICA Community Members, 2004]. Using fixed thresholds, unless they are placed very low, can imply that there were therefore no, few, or very short, interglacials within the interval prior to MIS 11. Half-threshold values avoid this issue, ensuring a finite length for all defined interglacials but tend to give greater durations than other methods. For some of the weaker interglacials (especially MIS 13), they imply a particularly long duration and cause interglacial length to accumulate in periods with values that would be strongly glacial in other cycles. The results of the normalizing methods depend on the criteria used, and particularly on the number of SD used for the threshold: with lower values than we have used, very long durations are again found for the weaker interglacials. The process-based approach of using the occurrence of the bipolar seesaw is unaffected by baseline shifts [Tzedakis et al., 2012b]. There is more scatter between methods for the earlier interglacials, reflecting different philosophies as to whether interglacial length is defined by changes relative to a baseline or by absolute conditions experienced.
For the same methodology, there appears to be little systematic difference between calculations using LR04, CO 2 , or Antarctic temperature proxy. For example, using the half-transition method, each of the three proxies gives a longer duration than the others for at least two interglacials; the same result is found using the normalized (400 ka interval) method. 7.5.2. Evaluation of Interglacial Duration Both in their individual estimates and particularly when compared to the Holocene duration using the same method, a number of consistent distinctions emerge among the previous interglacials (Table 2 and Figure 15). In

10.1002/2015RG000482
comparing durations with that of the Holocene (MIS 1), it must be remembered that the Holocene has not yet ended: under none of the criteria has it passed a threshold for inception.
Interglacials generally have a duration of~10-30 ka, consistent with the timescale of the astronomical forcing, but significantly less than the duration of the more glacial portion of climate cycles of the past 800 ka. The elapsed duration of the Holocene is generally estimated to be 10-15 ka, comparable to one half of a precession cycle. Few previous interglacials were shorter than this, and some were similar in duration. The most recent complete interglacial, MIS 5e, was longer, with a length approximately 1.5 times the length of the Holocene to date. MIS 9e is generally comparable in length, while MIS 7e is the shortest interglacial peak under most criteria used.
MIS 11c appears to be an exceptionally long interglacial. On all criteria, it is much longer than the Holocene to date, with a range of 20 ka to 40 ka depending on the individual method, and generally 2-3 times the Holocene duration using the same method. The MIS 11c duration is also longer in comparison to other interglacials, none of which (except MIS 13, see below) have mean durations twice as long as the Holocene. The exceptional length of MIS 11 has been noted previously [Berger and Loutre, 1996;Droxler et al., 2003;Loutre and Berger, 2003;McManus et al., 2003;EPICA Community Members, 2004;Rohling et al., 2010], with potentially important implications for the natural continuation of the Holocene [Berger and Loutre, 2002;McManus et al., 2003], which (until the last century) shares similar GHG concentrations and astronomical configuration, particularly the nearly circular orbit and damped precession of MIS 11. Results of transient model simulations [Yin and Berger, 2015] have been used to suggest that the length of MIS 11 may be influenced by the near antiphase between obliquity and precession as well as the weak precession itself. The other interglacial occurring at a time of weak precession (MIS 19c) has a length much shorter than MIS 11c, and in the mean and median, comparable to the current length of the Holocene; however, while MIS 11c is longer than either MIS 19c or elapsed MIS 1 in all methods, the relative length of MIS 19c compared to MIS 1 differs between methods.   [Lisiecki and Raymo, 2005], atmospheric CO 2 concentration from Antarctic ice cores [Lüthi et al., 2008;Bereiter et al., 2015], and δ-deuterium for the EPICA Dome C ice core (EDC deut) [Jouzel et al., 2007]. Methods are as described in the text, where T is threshold method, HT is half-threshold method, and "norm x ka, 1 SD" represents a normalization over an x ka window, with the threshold set at 1 standard deviation above the mean.
"na" means that the method was not applicable; "bt" means that values never reached the threshold for this interglacial. For purposes of averaging, bt was counted as zero, while na values were disregarded.

10.1002/2015RG000482
MIS 13 is exceptional in the fact that it is a particularly weak interglacial, so that in some definitions the length of the interglacial is zero, and by several others it is short. On the other hand, the time between the apparent termination event and the subsequent maximum in warmth or sea level is particularly long, so that in some criteria it stands out as a particularly long interglacial, encompassing an additional precession cycle and isotopic substage. Possibly this combination resulting in an unusually long period of intermediate warmth explains the surprising finding that Greenland may have been partially deglaciated in MIS 13 [de Vernal and Hillaire-Marcel, 2008].

Present and Future Interglacials
This review has identified and described the interglacials that have occurred over the last 800 ka and discussed some of the factors that may have controlled the diversity of interglacials we have observed. Important remaining questions relate to the climate of the present and the future. The first question is that of determinism and predictability. Given the (known) astronomical forcing and the feedbacks that are present, is the occurrence and character of interglacials (in terms of timing, strength, and shape) predictable? In other words, starting from 800 ka ago with only astronomical forcing prescribed, is it inevitable that we would find ourselves in today's interglacial climate following the same sequence of glacial and interglacials that has occurred? Several additional critical questions follow from this first one and relate to the timing of the next glacial inception. What is the natural progression of the current interglacial? Most previous interglacials have been observed to be as long as or longer than this one. Is it possible to predict when the next glacial inception would have occurred if greenhouse gases had followed a "natural" trajectory? Has human activity already played a role in prolonging this transition? If CO 2 concentrations continue to rise due to human activity as anticipated, by how long will the inception be delayed? We address these questions briefly in the concluding sections.

Predictability
It is widely accepted that astronomical variations pace the sequence of glacial-interglacial cycles. The appearance of the main astronomical periods in palaeoclimate records [Hays et al., 1976] and the evidence that recent terminations line up with periods of decreasing precession support this [Cheng et al., 2009]. This leads to a view that the ice age cycles are to a large degree deterministic and predictable. A number of different conceptual models are capable of producing the qualitative structure of the ice age cycles [Maasch and Saltzman, 1990;Paillard, 1998;Tziperman et al., 2006] and, with an appropriate choice of parameters, the timing of the observed sequence over the last 800 ka. In most models, the initial conditions do not affect the timing and amplitude of subsequent glacial cycles beyond a horizon of 200 ka [Parrenin and Paillard, 2012]. However, in some models, small changes in the parameters or stochastic perturbations can perturb the sequence, leading to later inception or to termination in a different precessional cycle [Paillard, 2001;Crucifix, 2013]. EMICs can produce a good Figure 15. Estimates for interglacial duration in paleoclimate records based on various methods for identifying each interglacial interval. (see Table 2). Interglacial intervals are labeled here using the marine isotope stage (MIS) nomenclature. Data sets are the benthic isotope stack (LR04) [Lisiecki and Raymo, 2005], atmospheric CO 2 concentration from Antarctic ice cores [Lüthi et al., 2008;Bereiter et al., 2015], and δ deuterium for the EPICA Dome C ice core (EDC Deut) [Jouzel et al., 2007]. Methods are as described in the text, where T is threshold method, HT is half-threshold method, and Norm x ka, 1SD represents a normalization over an x ka window, with the threshold set at 1 standard deviation above the mean. The mean and median are shown as bold symbols.

10.1002/2015RG000482
representation of the climate of the last 800 ka, when GHG and ice sheets are prescribed [Holden et al., 2010;Ganopolski and Calov, 2011]. Although it has been proposed that glacial cycles will still occur with constant CO 2 [Ganopolski and Calov, 2011;Abe-Ouchi et al., 2013], the exact shape and timing of interglacial occurrence is clearly dependent on both the CO 2 and other aspects of the model parameters. At this stage it therefore remains unclear whether the timing and shape of interglacials is predictable or not, and over what length of time [Parrenin and Paillard, 2012;Crucifix, 2013].

Next Glacial Inception
Most model experiments of the end of an interglacial were performed for the last glacial inception (circa 115 ka ago) when summer insolation was significantly lower than at present. However, over the last million years several glacial inceptions occurred under astronomical configuration similar to the present one. Based on model experiments, the glaciation threshold depends not only on insolation but also on the atmospheric CO 2 content [Archer and Ganopolski, 2005]. Models of different complexity and time series analyses of proxy data have been used to investigate the response to orbital forcing in the future for a range of atmospheric CO 2 levels. These results show consistently, that a glacial inception is unlikely to happen within the next approximate 50 ka (when the next strong drop in Northern Hemisphere summer insolation occurs) if either atmospheric CO 2 concentration remains above 300 ppm or cumulative carbon emissions exceed 1000 Pg C [Loutre and Berger, 2000;Archer and Ganopolski, 2005;Cochelin et al., 2006]. Only for an atmospheric CO 2 content below the preindustrial level may a glaciation occur within the next 10 ka [Loutre and Berger, 2000;Cochelin et al., 2006;Kutzbach et al., 2011;Vettoretti and Peltier, 2011;Tzedakis et al., 2012a]. This can be used to explain why glacial inception occurred during MIS 19c under astronomical configuration rather similar to the present one but at lower CO 2 concentration (~240 ppm). Given the continued anthropogenic CO 2 emissions, glacial inception is very unlikely to occur in the next 50 ka, because the timescale for CO 2 and temperature reduction toward unperturbed values in the absence of active removal is very long [IPCC, 2013], and only weak precessional forcing occurs in the next two precessional cycles.

Conclusions
Interglacials of the last 800 ka represent the warm, low ice extent (high sea level) end-member of the climate system during this period. There has been no single definition of interglacials that satisfies all the attributes we associate with them, but a working definition based on a sea level close to present is becoming practical with the emergence of new sea level data sets and is supported by sensitivity studies on other data sets. Using such a definition (in which there is little or no ice in the Northern Hemisphere outside Greenland), we have identified 11 interglacials in the last 800 ka. Using the median values for duration in Table 2, interglacials represent about 20% of the past 800 ka (but note that for some of the definitions included in the calculation, such a percentage was prescribed by the method).
Although there is a lack of available terrestrial data sets, there are now sufficient well-resolved long records, particularly of SSTs, that it is possible to assess the pattern of each interglacial. There is, however, no apparent template for the shape, amplitude, duration, or spatial pattern of an interglacial. They all start with a sea level rise that is rapid compared to the length of a glacial cycle, and highly nonlinear in relation to the external forcing. The interglacials that are seen in the ocean as MIS 5e and 11c show a global strength (based mainly on temperature records), while MIS 13a is generally weaker. There is, however, no simple astronomical cause for the pattern of amplitude between interglacials, which seems to arise at least partly from the pattern of CO 2 . This emphasizes the need to better understand and model the carbon cycle across glacial cycles.
Snapshot, and a few transient, model runs, using boundary conditions appropriate for previous interglacials, are now becoming available. The majority of them are for the last interglacial (MIS 5e), and while there are some areas of agreement with data, it is already apparent that the interhemispheric pattern of climate across that interglacial cannot be explained unless significant millennial bipolar seesaw activity (presumably due to changes in AMOC) is taken into account. This poses a challenge for models (as the correct forcing will be hard to deduce) and offers hints for understanding how astronomical forcing and millennial-scale events interact to produce glacial terminations.

10.1002/2015RG000482
Glacial inception involves a sequence of events that still needs to be better documented for most interglacials prior to the last one. There is not yet a clear understanding of what controls the length of each interglacial, and there is no perfect analog for the present interglacial. However, inception at preindustrial greenhouse gas concentrations seems unlikely without a strong precession cycle, suggesting a long interglacial ahead of us. Anthropogenic additions to these gases further increase the likelihood of a long interglacial.
Within this paper we have identified a number of key concepts, which also give rise to some key questions that remain to be solved. We summarize these here: 1. A definition of interglacials that respects the origin of the term but also the implications of modern data sets would be that they are those climate states with little or no Northern Hemisphere ice outside Greenland. However, available data sets do not yet provide a robust and consistent identification of these intervals. a. Can a sea level record be constructed with sufficient precision and confidence that such a definition becomes a practical reality? 2. While the intensity of different interglacials varies geographically and among the different climate indices, the concentration of CO 2 appears to play a key role in influencing the relative strength of each interglacials. a. What determines the CO 2 concentration and trend during each interglacial? 3. High ice volume seems to be a strong indicator for the subsequent occurrence of an interglacial within a particular precession cycle, although strong insolation forcing may also be capable of returning the climate system to its interglacial state.
a. What are the mechanistic links (ice-lithosphere interaction, dust, ocean circulation and carbon cycle, other) that drive this observation? b. How does this apply to interglacials before 800 ka? 4. Millennial variability, involving a bipolar seesaw signature, seems crucial for the relative timing of changes in different climatic parameters and at different locations during glacial terminations and interglacial onset.
a. To what extent does stochastic variability (such as changes in volcanic and solar forcing, and internal variability) affect the precise timing of interglacial onset and demise? b. To what extent is the pattern of millennial variability at interglacial onset also important for trends during each subsequent interglacial? c. And alternatively, to what extent is the relative phasing of different orbital parameters responsible? 5. It remains a challenge, only partly addressed in the last three concepts, to quantitatively link causal factors (either primary astronomical forcing or secondary factors such as ice volume and CO 2 concentration) to observed interglacial properties (intensity, timing, and duration). a. This nonetheless remains an aim, which can be pursued as further data make new data analyses possible; its feasibility is, however, closely related to the next point. 6. It is possible to create models of varying complexity that do a reasonable job of reproducing the pattern of interglacial occurrence in the last 800 ka.
a. Does this mean that interglacials are predictable? b. Are such models sufficiently constrained by the data and by understanding of processes that we can confidently predict the timing of the end of the current interglacial for a range of CO 2 concentration scenarios?
A number of more technical research issues need to be solved if we are to make further breakthroughs in understanding recent interglacials. The data allow a broad overview of the pattern of each interglacial. However, chronological advances, both in assessing absolute ages relative to astronomical forcing and in aligning different proxies and locations, are essential if we are to assess the dynamics of interglacials and their termination and inception. The paucity of terrestrial records precludes the assessment of many important aspects of climate. Although sea level is quasi-similar at each interglacial apex, better knowledge of the state of the Greenland and Antarctic ice sheets in stronger interglacials could be critical. Identifying the controls on intra-interglacial variability remains a challenge. It appears that inclusion of realistic freshwater forcing is essential for simulating terminations and the interhemispheric pattern in climate models. Furthermore, while the pattern of glacial-interglacial cycles is well-understood if CO 2 is prescribed, we will not be able to say we understand it until realistic modeling of the carbon cycle over these cycles is achieved. Finally, it remains unrealistic to run the latest state of the art ESMs over the full 800 ka. While this must be a goal for future model Reviews of Geophysics 10.1002/2015RG000482 development, conceptual understanding of glacial cycles may also involve better integration of some of the ideas addressed above into simpler models.