LIG Climate Research Articles

Southern Ocean sea surface temperature synthesis: Part 2. Penultimate glacial and last interglacial

The last interglacial (LIG: ∼130 to 115 thousand years before present) is often used as an analogue for near-future climate warming. Antarctic Ice Sheet response to LIG warming is of particular interest, because of its implications for sea level rise. Comparison between LIG climate simulations and proxy-based reconstructions of Southern Ocean sea surface temperature (SST) remains challenging, due to high uncertainties in both reconstructions and simulations. In this two-part study, the accompanying paper (Part 1) addressed uncertainties in the SST reconstructions by evaluating proxies relevant to Southern Ocean SST, and made recommendations for which proxies and respective calibrations are most reliable on glacial-interglacial time scales in this region. In the second part (this paper), we now apply these recommendations to a synthesis of Southern Ocean SST over the penultimate glacial and LIG. Similar to previous studies, we find that LIG warming at 40°S to 60°S reached 1.6 ± 1.1 °C (annual mean) or 1.9 ± 1.3 °C (austral summer: JFM) relative to present. Annual/summer cooling in the penultimate glacial maximum reached −3.6 ± 1.0 °C/−4.0 ± 1.2 °C, similar to the last glacial maximum. Compared with the previous LIG SST syntheses, our reported uncertainties more strongly reflect geographic variability and dating errors, as we have reduced errors in the individual temperature reconstructions and do not date records by aligning peaks in their SST. However, the reconstruction errors are still important, and we do not recommend detailed interpretation of temperature records from small numbers of sites. Instead, comparisons of our new synthesis with model simulations should focus only on the regional average.

Contrasting northern and southern European winter climate trends during the Last Interglacial

Abstract The Last Interglacial (LIG; 130–115 ka) is an important test bed for climate science as an instance of significantly warmer than preindustrial global temperatures. However, LIG climate patterns remain poorly resolved, especially for winter, affected by a suite of strong feedbacks such as changes in sea-ice cover in the high latitudes. We present a synthesis of winter temperature and precipitation proxy data from the Atlantic seaboard of Europe, spanning from southern Iberia to the Arctic. Our data reveal distinct, opposite latitudinal climate trends, including warming winters seen in the European Arctic while cooling and drying occurred in southwest Europe over the LIG. Climate model simulations for 130 and 120 ka suggest these contrasting climate patterns were affected by a shift toward an atmospheric circulation regime with an enhanced meridional pressure gradient and strengthened midlatitude westerlies, leading to a strong reduction in precipitation across southern Europe.

Potential Himalayan community turnover through the Late Pleistocene

Prevailing models for the high biodiversity in tropical mountains assume that organisms can survive through past climate change by performing short elevational movements to track suitable habitats. However, dramatic Pleistocene climatic oscillations could also lead to species turnover but have widely been ignored. Here, we used ecological niche modelling (ENM) of 288 passerine species in the Himalayas to test the effect of climate change during the Last Interglacial Period (LIG), the Last Glacial Maximum (LGM) and the present day. The ENM analyses hindcasted species persistence through climate change from the LGM to the present day but likely showed a high degree of species turnover (e.g. 32.6–46.2%) from the LIG to the LGM. Further elevational dynamic reconstructions demonstrated that species might survive these two periods of climate change by upward and downward shifts, respectively. Statistical analyses of climatic variables showed increased climatic variability in the Himalayas during the LIG, which might have caused community turnover. The severe evolutionary consequence of the LIG climate in the Himalayas contrasts with the paradigm of the climatic optimum in Europe and North America and suggests potential geography-dependent effects of past climate change. More importantly, our results demonstrate that dramatic historical climate change might overwhelm the buffering effect of elevational heterogeneity, which should be considered when investigating the origin of tropical montane biodiversity.

Holocene and Last Interglacial climate of the Faroe Islands from sedimentary plant wax hydrogen and carbon isotopes

The Last Interglacial period (LIG) is Earth's most recent globally warm period and is analogous in some ways to projected future global warming. However, questions remain regarding the state of the climate during the LIG in the North Atlantic, a region that is extremely sensitive to changes in oceanic and atmospheric circulation. Here, we present hydrogen and carbon isotope (δD and δ13C) records from a suite of plant wax biomarkers preserved in Holocene and LIG lacustrine sediments from the North Atlantic Faroe Islands and interpret them as qualitative proxies for temperature and hydroclimate variability. These data are used to directly compare LIG and Holocene climate using the same proxy approaches from the same terrestrial location. Measuring multiple isotopes on multiple types of waxes elucidates the sources of homologous plant waxes. We deduce that the δD values of long-chain n-alkanes (C27–C33) and mid-chain n-alkanes (C23–C25) in these sedimentary archives reflect leaf water and lake water δD values, respectively, while the δD values for both long-chain (C28–C30) and mid-chain n-alkanoic acids (C24–C26) primarily represent lake water δD values. Plant wax-inferred δD values of precipitation during the early Holocene (10,100 to 8,200 cal yr BP) are ∼35‰ more positive than late Holocene values, and decline over the Holocene. δD-inferred hydrologic change and δ13C -inferred plant water use efficiency both indicate that the Faroe Islands became drier throughout the Holocene. Comparison with measurements from LIG plant waxes indicates that late LIG in the Faroe Islands was hydrologically similar to the early-to mid-Holocene (8,200 to 4,000 cal yr BP), with enriched precipitation isotopes and reduced evapotranspiration indicating a warmer, wetter environment.

The penultimate deglaciation: protocol for Paleoclimate Modelling Intercomparison Project (PMIP) phase 4 transient numerical simulations between 140 and 127 ka, version 1.0

Abstract. The penultimate deglaciation (PDG, ∼138–128 thousand years before present, hereafter ka) is the transition from the penultimate glacial maximum (PGM) to the Last Interglacial (LIG, ∼129–116 ka). The LIG stands out as one of the warmest interglacials of the last 800 000 years (hereafter kyr), with high-latitude temperature warmer than today and global sea level likely higher by at least 6 m. Considering the transient nature of the Earth system, the LIG climate and ice-sheet evolution were certainly influenced by the changes occurring during the penultimate deglaciation. It is thus important to investigate, with coupled atmosphere–ocean general circulation models (AOGCMs), the climate and environmental response to the large changes in boundary conditions (i.e. orbital configuration, atmospheric greenhouse gas concentrations, ice-sheet geometry and associated meltwater fluxes) occurring during the penultimate deglaciation. A deglaciation working group has recently been set up as part of the Paleoclimate Modelling Intercomparison Project (PMIP) phase 4, with a protocol to perform transient simulations of the last deglaciation (19–11 ka; although the protocol covers 26–0 ka). Similar to the last deglaciation, the disintegration of continental ice sheets during the penultimate deglaciation led to significant changes in the oceanic circulation during Heinrich Stadial 11 (∼136–129 ka). However, the two deglaciations bear significant differences in magnitude and temporal evolution of climate and environmental changes. Here, as part of the Past Global Changes (PAGES)-PMIP working group on Quaternary interglacials (QUIGS), we propose a protocol to perform transient simulations of the penultimate deglaciation under the auspices of PMIP4. This design includes time-varying changes in orbital forcing, greenhouse gas concentrations, continental ice sheets as well as freshwater input from the disintegration of continental ice sheets. This experiment is designed for AOGCMs to assess the coupled response of the climate system to all forcings. Additional sensitivity experiments are proposed to evaluate the response to each forcing. Finally, a selection of paleo-records representing different parts of the climate system is presented, providing an appropriate benchmark for upcoming model–data comparisons across the penultimate deglaciation.

The penultimate deglaciation: protocol for PMIP4 transient numerical simulations between 140 and 127 ka, version 1.0

Abstract. The penultimate deglaciation (PDG, ~ 138–128 thousand years before present, hereafter ka) is the transition from the penultimate glacial maximum to the Last Interglacial (LIG, ~ 129–116 ka). The LIG stands out as one of the warmest interglacials of the last 800 ka, with high-latitude temperature warmer than today and global sea level likely higher by at least 6 meters. Considering the transient nature of the Earth system, the LIG climate and ice-sheets evolution were certainly influenced by the changes occurring during the penultimate deglaciation. It is thus important to investigate, with coupled Atmosphere-Ocean General Circulation Models (AOGCMs), the climate and environmental response to the large changes in boundary conditions (i.e. orbital configuration, atmospheric greenhouse gas concentrations, ice-sheet geometry, and associated meltwater fluxes) occurring during the penultimate deglaciation. A deglaciation working group has recently been set up as part of the Paleoclimate Modelling Intercomparison Project (PMIP) phase 4, with a protocol to perform transient simulations of the last deglaciation (19–11 ka; although the protocol covers 26–0 ka). Similar to the last deglaciation, the disintegration of continental ice-sheets during the penultimate deglaciation led to significant changes in the oceanic circulation during Heinrich Stadial 11 (~ 136–129 ka). However, the two deglaciations bear significant differences in magnitude and temporal evolution of climate and environmental changes. Here, as part of the PAGES-PMIP working group on Quaternary Interglacials, we propose a protocol to perform transient simulations of the penultimate deglaciation under the auspices of PMIP4. This design includes time-varying changes in orbital forcing, greenhouse gas concentrations, continental ice-sheets as well as freshwater input from the disintegration of continental ice-sheets. This experiment is designed for AOGCMs to assess the coupled response of the climate system to all forcings. Additional sensitivity experiments are proposed to evaluate the response to each forcing. Finally, a selection of paleo records representing different parts of the climate system is presented, providing an appropriate benchmark for upcoming model-data comparisons across the penultimate deglaciation.

Recent advances on the dynamical representation and our understanding of the warmer‑than‑present last interglacial climate

The Last Interglacial (LIG, ~129-116 thousand years Before Present, hereafter ka) represents an ideal case study to understand the climate mechanisms at play under a warmer-than-present climate. However a spatio-temporal representation of the LIG climatic changes remains difficult to obtain, mainly because aligning paleoclimatic records from various archives (i.e. polar ice cores, marine sediments, speleothems) from around the globe is challenging. Here we summarize recent studies that highlight how the coupling of a synthesis of surface air and sea temperature records (above polar ice sheets and from the North Atlantic and Southern Ocean respectively) associated with harmonized chronologies and of appropriate climate model experiments improved our spatio-temporal representation of the LIG high-latitude climate evolution, and our understanding of the mechanisms at play, especially at the beginning of the LIG. In particular, we describe commonly-used record alignment strategies for marine sediment cores and we show that age discrepancies larger than 4 ka can exist between the timescales inferred from the different approaches. Providing harmonized chronologies when comparing multiple records is thus essential and we propose a new high latitude LIG data synthesis based on coherent record time scales together with associated time slices at 130, 125, 120 and 115 ka of surface temperature anomalies relative to present-day. The results provide the first robust evidence for asynchronous surface temperature evolutions at the LIG onset across the world and also enable one to identify important missing processes in state-of-the-art model climate simulations to reproduce the early LIG climate evolution. Our integrated model-data approach shows that a freshwater input into the North Atlantic (due to the Northern Hemisphere ice sheet early melting) needs to be accounted for in addition to the orbital and greenhouse gas concentration forcing in climate simulations, to explain the evolution of the early LIG climate.

Critical evaluation of climate syntheses to benchmark CMIP6/PMIP4 127 ka Last Interglacial simulations in the high-latitude regions

The Last Interglacial (LIG, ∼129-116 thousand years ago, ka) represents an excellent case study to investigate the response of sensitive components of the Earth System and mechanisms of high-latitude amplification to a climate warmer than present-day. The Paleoclimate Model Intercomparison Project (Phase 4, hereafter referred as PMIP4) and the Coupled Model Intercomparison Project (Phase 6, hereafter referred as CMIP6) are coordinating the design of (1) a LIG Tier 1 equilibrium simulation to simulate the climate response at 127 ka, a time interval associated with a strong orbital forcing and greenhouse gas concentrations close to preindustrial levels and (2) associated Tier 2 sensitivity experiments to examine the role of the ocean, vegetation and dust feedbacks in modulating the response to this orbital forcing.Evaluating the capability of the CMIP6/PMIP4 models to reproduce the 127 ka polar and sub-polar climate will require appropriate data-based benchmarks which are currently missing. Based on a recent data synthesis that offers the first spatio-temporal representation of high-latitude (i.e. poleward of 40°N and 40°S) surface temperature evolution during the LIG, we produce a new 126–128 ka time slab, hereafter named 127 ka time slice. This 127 ka time slice represents surface temperature anomalies relative to preindustrial and is associated with quantitative estimates of the uncertainties related to relative dating and surface temperature reconstruction methods. It illustrates warmer-than-preindustrial conditions in the high-latitude regions of both hemispheres. In particular, summer sea surface temperatures (SST) in the North Atlantic region were on average 1.1 °C (with a standard error of the mean of 0.7 °C) warmer relative to preindustrial and 1.8 °C (with a standard error of the mean of 0.8 °C) in the Southern Ocean. In Antarctica, average 127 ka annual surface air temperature was 2.2 °C (with a standard error of the mean of 1.4 °C) warmer compared to preindustrial.We provide a critical evaluation of the latest LIG surface climate compilations that are available for evaluating LIG climate model experiments. We discuss in particular our new 127 ka time-slice in the context of existing LIG surface temperature time-slices. We also compare the 127 ka time slice with the ones published for the 125 and 130 ka time intervals and we discuss the potential and limits of a data-based time slice at 127 ka in the context of the upcoming coordinated modeling exercise. Finally we provide guidance on the use of the available LIG climate compilations for future model-data comparison exercises in the framework of the upcoming CMIP6/PMIP4 127 ka experiments. We do not recommend the use of LIG peak warmth-centered syntheses. Instead we promote the use of the most recent syntheses that are based on coherent chronologies between paleoclimatic records and provide spatio-temporal reconstruction of the LIG climate. In particular, we recommend using our new 127 ka data-based time slice in model-data comparison studies with a focus on the high-latitude climate.

Pleistocene Niche Stability and Lineage Diversification in the Subtropical Spider Araneus omnicolor (Araneidae)

The influence of Quaternary climate oscillations on the diversification of the South American fauna is being increasingly explored. However, most of these studies have focused on taxa that are endemic to tropical environments, and relatively few have treated organisms restricted to subtropical biomes. Here we used an integrative phylogeographical framework to investigate the effects of these climate events on the ecological niche and genetic patterns of the subtropical orb-weaver spider Araneus omnicolor (Araneidae). We analyzed the mitochondrial (Cytochrome Oxidase I, COI) and nuclear (Internal Transcribed Subunit II, ITS2) DNA of 130 individuals throughout the species’ range, and generated distribution models in three different climate scenarios [present, Last Glacial Maximum (LGM), and Last Interglacial Maximum (LIG)]. Additionally, we used an Approximate Bayesian Computation (ABC) approach to compare possible demographic scenarios and select the hypothesis that better explains the genetic patterns of A. omnicolor. We obtained high haplotype diversity but low nucleotide variation among sequences. The population structure and demographic analyses showed discrepancies between markers, suggesting male-biased dispersal in the species. The time-calibrated COI phylogenetic inference showed a recent diversification of lineages (Middle/Late Pleistocene), while the paleoclimate modeling indicated niche stability since ~120 Kya. The ABC results agreed with the niche models, supporting a panmictic population as the most likely historical scenario for the species. These results indicate that A. omnicolor experienced no niche or population reductions during the Late Pleistocene, despite the intense landscape modifications that occurred in the subtropical region, and that other factors beside LGM and LIG climate oscillations might have contributed to the demographic history of this species. This pattern may be related to the high dispersal ability and wide environmental tolerance of A. omnicolor, highlighting the need for more phylogeographical studies with invertebrates and other generalist taxa, in order to understand the effects of Quaternary climate changes on Neotropical biodiversity.

The Greenland Ice Sheet at the peak of warming during the previous Interglacial

The Last Interglacial (LIG or the Eemian) between ca. 130 and 115 kyr BP is probably the best analogue for future climate warming for which increasingly better proxy data are becoming available. The volume of the Greenland Ice Sheet (GrIS) during this period is of particular interest to better assess how much and how fast sea-level can rise in a future Earth undergoing gradual climatic warming. Sea-level during the LIG is inferred to have been up to 9 meter higher than today, but contribution of the GrIS into this rise remains unclear. Various ice-sheet modeling studies have come up with a very broad range of the LIG volume loss by the GrIS to between 60 cm and 6 m of equivalent sea-level rise. This wide range is explained by the sensitivity of GrIS models to the imposed climatic conditions and to poor knowledge of the LIG climate itself in terms of the magnitude and precise timing of the maximum warming, as well as in terms of spatial and annual patterns. To partially circumvent these uncertainties we made use of the newest temperature record over the Central Greenland reconstructed from the isotopic composition of the recently obtained NEEM ice core containing undisturbed LIG segment to build the climatic forcing of the model. The NEEM record unequivocally indicates times of the start and of the end of the LIG warming in Greenland as well as the duration of the warmest time period within the Eemian. Using a three-dimensional thermomechanical ice-sheet model, we produced an ensemble of possible LIG configurations by varying only four key parameters for temperature, precipitation rate, surface melting magnitude and melting pattern within realistic bounds. The outcome of a series of the numerical experiments is a variety of glaciologically consistent GrIS geometries corresponding to a wide range of possible «climates». To constrain the ensemble of GrIS geometries, we used data inferred from 5 Greenland ice cores such as the presence or absence of LIG ice, borehole temperature and isotopic composition. Lagrangian backtracing of particles was used to calculate non-climatic biases in isotopic records introduced by horizontal advection, systematic latitudinal contrast and local elevation changes. Comparison of model-generated ice-core characteristics with the observed data enabled to narrow down the ensemble to a bound on the GrIS contribution to the LIG sea-level rise of between 1.3 and 2.9 m with the best choice of 1.8–2.2 m. This conclusion in general supports the point of view about the modest GrIS contribution to global sea level rise during the LIG.

Temporal and spatial structure of multi-millennial temperature changes at high latitudes during the Last Interglacial

The Last Interglacial (LIG, 129–116 thousand of years BP, ka) represents a test bed for climate model feedbacks in warmer-than-present high latitude regions. However, mainly because aligning different palaeoclimatic archives and from different parts of the world is not trivial, a spatio-temporal picture of LIG temperature changes is difficult to obtain.Here, we have selected 47 polar ice core and sub-polar marine sediment records and developed a strategy to align them onto the recent AICC2012 ice core chronology. We provide the first compilation of high-latitude temperature changes across the LIG associated with a coherent temporal framework built between ice core and marine sediment records. Our new data synthesis highlights non-synchronous maximum temperature changes between the two hemispheres with the Southern Ocean and Antarctica records showing an early warming compared to North Atlantic records. We also observe warmer than present-day conditions that occur for a longer time period in southern high latitudes than in northern high latitudes. Finally, the amplitude of temperature changes at high northern latitudes is larger compared to high southern latitude temperature changes recorded at the onset and the demise of the LIG.We have also compiled four data-based time slices with temperature anomalies (compared to present-day conditions) at 115 ka, 120 ka, 125 ka and 130 ka and quantitatively estimated temperature uncertainties that include relative dating errors. This provides an improved benchmark for performing more robust model-data comparison. The surface temperature simulated by two General Circulation Models (CCSM3 and HadCM3) for 130 ka and 125 ka is compared to the corresponding time slice data synthesis. This comparison shows that the models predict warmer than present conditions earlier than documented in the North Atlantic, while neither model is able to produce the reconstructed early Southern Ocean and Antarctic warming. Our results highlight the importance of producing a sequence of time slices rather than one single time slice averaging the LIG climate conditions.

Simulating last interglacial climate with NorESM: role of insolation and greenhouse gases in the timing of peak warmth

Abstract. The last interglacial (LIG, ~130–116 ka, ka = 1000 yr ago) is characterized by high-latitude warming and is therefore often considered as a possible analogue for future warming. However, in contrast to predicted future greenhouse warming, the LIG climate is largely governed by variations in insolation. Greenhouse gas (GHG) concentrations were relatively stable and similar to pre-industrial values, with the exception of the early LIG when, on average, GHGs were slightly lower. We performed six time-slice simulations with the low-resolution version of the Norwegian Earth System Model covering the LIG. In four simulations only the orbital forcing was changed. In two other simulations, representing the early LIG, additionally the GHG forcing was reduced. With these simulations we investigate (1) the different effects of GHG versus insolation forcing on the temperatures during the LIG; (2) whether reduced GHGs can explain the low temperatures reconstructed for the North Atlantic; and (3) the timing of the observed LIG peak warmth. Our simulations show that the insolation forcing results in seasonal and hemispheric differences in temperature. In contrast, a reduction in the GHG forcing causes a global and seasonal-independent cooling. Furthermore, we compare modelled temperatures with proxy-based LIG sea-surface temperatures along a transect in the North Atlantic. The modelled North Atlantic summer sea-surface temperatures capture the general trend of the reconstructed summer temperatures, with low values in the early LIG, a peak around 125 ka, and a steady decrease towards the end of the LIG. Simulations with reduced GHG forcing improve the model–data fit as they show lower temperatures in the early LIG. Furthermore we show that the timing of maximum summer and winter surface temperatures is in line with the local summer and winter insolation maximum at most latitudes. Two regions where the maximum local insolation and temperature do not occur at the same time are Antarctica and the Southern Ocean. The austral summer insolation has a late maximum at ~115 ka. In contrast the austral summer temperatures in Antarctica show maxima at both ~130 ka and ~115 ka, and the Southern Ocean temperatures peak only at ~130 ka. This is probably due to the integrating effect of the ocean, storing heat from other seasons and resulting in relatively warm austral summer temperatures. Reducing the GHG concentrations in the early LIG (125 and 130 ka) results in a similar timing of peak warmth, except over Antarctica. There, the lower austral summer temperatures at 130 ka shift the maximum warmth to a single peak at 115 ka.

Early Last Interglacial Greenland Ice Sheet melting and a sustained period of meridional overturning weakening: a model analysis of the uncertainties

Proxy-data suggest that the Last Interglacial (LIG; ~130–116 ka BP) climate was characterized by higher temperatures, a partially melted Greenland Ice Sheet (GIS) and a changed Atlantic meridional overturning circulation (AMOC). Notwithstanding the uncertainties in LIG palaeoclimatic reconstructions, this setting potentially provides an opportunity to evaluate the relation between GIS melt and the AMOC as simulated by climate models. However, first we need to assess the extent to which a causal relation between early LIG GIS melt and the weakened AMOC is plausible. With a series of transient LIG climate simulations with the LOVECLIM earth system model, we quantify the importance of the major known uncertainties involved in early LIG GIS melt scenarios. Based on this we construct a specific scenario that is within the parameter space of uncertainties and show that it is physically consistent that early LIG GIS melting kept the AMOC weakened. Notwithstanding, this scenario is at the extreme end of the parameter space. Assuming that proxy-based reconstructions of early LIG AMOC weakening offer a realistic representation of its past state, this indicates that either (1) the AMOC weakening was caused by other forcings than early LIG GIS melt or (2) the early LIG AMOC was less stable than indicated by our simulations and a small amount of GIS melt was sufficient to keep the AMOC in the weak state of a bi-stable regime. We argue that more intensive research is required because of the high potential of the early LIG to evaluate model performance in relation to the AMOC response to GIS melt.

Sensitivity of the North Atlantic climate to Greenland Ice Sheet melting during the Last Interglacial

Abstract. During the Last Interglacial (LIG; ~130 000 yr BP), part of the Greenland Ice Sheet (GIS) melted due to a warmer than present-day climate. However, the impact of this melting on the LIG climate in the North Atlantic region is relatively unknown. Using the LOVECLIM Earth system model of intermediate complexity, we have systematically tested the sensitivity of the LIG climate to increased freshwater runoff from the GIS. In addition, experiments have been performed to investigate the impact of an idealized reduction of both surface elevation and extent of the GIS on the LIG climate. Based on changes in the maximum sea-ice cover and the strength of the overturning circulation, three regimes have been identified, which are characterized by a specific pattern of surface temperature change in the North Atlantic region. By comparing the simulated deep ocean circulation with proxy-based reconstructions, the most realistic simulated climate could be discerned. The resulting climate is characterized by a shutdown of deep convection and a subsequent ~4 °C cooling in the Labrador Sea. Furthermore, a cooling of ~1 °C over the North Atlantic Ocean between 40° N and 70° N is seen. The prescribed reduction in surface elevation and extent of the GIS results in a local warming of up to 4 °C and amplifies the freshwater-forced reduction in deep convection and the resultant cooling in the Nordic Seas. A further comparison of simulated summer temperatures with both continental and oceanic proxy records reveals that the partial melting of the GIS during the LIG could have delayed maximum summer temperatures in the western part of the North Atlantic region relative to the insolation maximum.