Journal cover Journal topic
Climate of the Past An interactive open-access journal of the European Geosciences Union
Journal topic
CP | Articles | Volume 16, issue 4
Clim. Past, 16, 1451–1467, 2020
© Author(s) 2020. This work is distributed under
the Creative Commons Attribution 4.0 License.
Clim. Past, 16, 1451–1467, 2020
© Author(s) 2020. This work is distributed under
the Creative Commons Attribution 4.0 License.

Research article 07 Aug 2020

Research article | 07 Aug 2020

Sea surface temperature in the Indian sector of the Southern Ocean over the Late Glacial and Holocene

Sea surface temperature in the Indian sector of the Southern Ocean over the Late Glacial and Holocene
Lisa Claire Orme1,2, Xavier Crosta3, Arto Miettinen1,4, Dmitry V. Divine1, Katrine Husum1, Elisabeth Isaksson1, Lukas Wacker5, Rahul Mohan6, Olivier Ther3, and Minoru Ikehara7 Lisa Claire Orme et al.
  • 1Norwegian Polar Institute, Tromsø, 9296, Norway
  • 2ICARUS, Department of Geography, Maynooth University, Maynooth, Ireland
  • 3UMR 5805 EPOC, Université de Bordeaux, Bordeaux, France
  • 4Ecosystems and Environment Research Programme, University of Helsinki, Helsinki, Finland
  • 5Department of Physics, ETH Zürich, 8093 Zurich, Switzerland
  • 6National Centre for Polar and Ocean Research (NCPOR), Ministry of Earth Sciences, Vasco-da-Gama, Goa 403804, India
  • 7Center for Advanced Marine Core Research, Kochi University, Nankoku, 783-8502, Japan

Correspondence: Lisa Claire Orme (

Back to toptop

Centennial- and millennial-scale variability of Southern Ocean temperature over the Holocene is poorly known, due to both short instrumental records and sparsely distributed high-resolution temperature reconstructions, with evidence for past temperature variations in the region coming mainly from ice core records. Here we present a high-resolution (∼60 year), diatom-based sea surface temperature (SST) reconstruction from the western Indian sector of the Southern Ocean that spans the interval 14.2 to 1.0 ka (calibrated kiloyears before present). During the late deglaciation, the new SST record shows cool temperatures at 14.2–12.9 ka and gradual warming between 12.9 and 11.6 ka in phase with atmospheric temperature evolution. This supports the evolution of the Southern Ocean SST during the deglaciation being linked with a complex combination of processes and drivers associated with reorganisations of atmospheric and oceanic circulation patterns. Specifically, we suggest that Southern Ocean surface warming coincided, within the dating uncertainties, with the reconstructed slowdown of the Atlantic Meridional Overturning Circulation (AMOC), rising atmospheric CO2 levels, changes in the southern westerly winds and enhanced upwelling. During the Holocene the record shows warm and stable temperatures from 11.6 to 8.7 ka followed by a slight cooling and greater variability from 8.7 to 1 ka, with a quasi-periodic variability of 200–260 years identified by spectral analysis. We suggest that the increased variability during the mid- to late Holocene reflects the establishment of centennial variability in SST connected with changes in the high-latitude atmospheric circulation and Southern Ocean convection.

1 Introduction
Back to toptop

Research into the sequence of events associated with the deglacial period has highlighted that there were contrasting patterns of millennial-scale warming and cooling between the southern and northern high latitudes, considered to have resulted from global-scale processes, namely fluctuations in the strength of the Atlantic Meridional Overturning Circulation (AMOC) and latitudinal shifts in the major atmospheric circulation cells. For example, during the Younger Dryas (12.9–11.7 ka; Rasmussen et al., 2014) a prolonged weakening of the AMOC reduced the northward heat transport in the Atlantic allowing heat to accumulate in the Southern Ocean causing a “bipolar seesaw” characterised by opposite ocean temperature anomalies between the two hemispheres (Denton et al., 2010, and references therein). These temperature changes caused a rapid alteration of the atmospheric circulation, with a southward shift of the Intertropical Convergence Zone (ITCZ) and southern westerly winds (SWWs). A combination of processes including enhanced eddy transport of heat across the Antarctic Circumpolar Current (ACC), upwelling of warm water and dissipation of sea ice then resulted in the warming of the Southern Ocean and Antarctica (Levermann et al., 2007; Screen et al., 2009; Denton et al., 2010). Finally, the increased upwelling in the Southern Ocean enhanced CO2 outgassing and thus warming globally. Despite many of these processes occurring within the Southern Ocean the strongest evidence for the patterns and magnitude of southern high-latitude warming during the Younger Dryas comes from Antarctic ice cores (e.g. Stenni et al., 2011). Although a number of records from the Southern Ocean show warming during the Younger Dryas, the timing, magnitude and duration of warming are far from consistent between records (Bianchi and Gersonde, 2004; Divine et al., 2010; Siani et al., 2013; Xiao et al., 2016). In the same vein, the cooler Antarctic Cold Reversal (ACR; 14.7–12.7 ka; Stenni et al., 2011), preceding the Younger Dryas, similarly displays a strong regional heterogeneity that is not fully understood (Pedro et al., 2015). Acquiring high-resolution proxy records of past sea surface temperature (SST) is therefore important for establishing the spatio-temporal patterns of warming in the Southern Ocean and their link with Antarctic temperature variability.

During the Holocene the majority of records from the open ocean of the Atlantic and Indian sectors of the Southern Ocean show that the early Holocene, defined as the period 11.7–8.2 ka (Walker et al., 2018), was warmer than the subsequent mid- to late Holocene (e.g. Bianchi and Gersonde, 2004; Nielsen et al., 2004; Anderson et al., 2009; Divine et al., 2010; Xiao et al., 2016). The suggested causes of the temperature evolution during the Holocene include variations in Southern Hemisphere insolation (e.g. Shevenell et al., 2011; Xiao et al., 2016), in the interplay of Southern and Northern Hemisphere insolation (e.g. Nielsen et al., 2004), and in the weakening of the AMOC in response to Northern Hemisphere deglaciation (Renssen et al., 2010).

There is also evidence of centennial variability of southern high-latitude climate through the Holocene as shown by some palaeoclimate records, including those from the Southern Ocean (e.g. Leventer et al., 1996; Crosta et al., 2007; Katsuki et al., 2012) and those reflecting atmospheric circulation changes such as the Southern Annular Mode (SAM) (e.g. Moreno et al., 2014; Turney et al., 2016). Observational and modelling evidence has suggested that centennial fluctuations in deep convection in the Southern Ocean (Weddell Sea) lead to variability in the surface water temperature of the Southern Ocean and the strength of the atmospheric circulation (Latif et al., 2013). Others have linked the oscillations to AMOC variability (Debret et al., 2009) or solar activity (e.g. Leventer et al., 1996; Crosta et al., 2007; Turney et al., 2016). To understand further how the Southern Ocean temperature has varied during the Holocene, there is a need for high-resolution records capturing centennial- to multi-centennial-scale variations.

We here present a new diatom-based SST reconstruction from the western Indian sector of the Southern Ocean spanning the last 14.2 kyr with an average temporal resolution of 60 years, to show the patterns, timing and magnitude of temperature variability through the late deglaciation and the Holocene.

2 Study area
Back to toptop

The Conrad Rise is located in the west Indian sector of the Southern Ocean (around 54  S, 40 E, Fig. 1) close to the Atlantic Sector. The rise is at <3000 m depth, with two sea mounts (Ob Bank and Lena Seamount) at depths of 300 m (Ansorge et al., 2008). The Conrad Rise is situated within the diatom ooze belt, which circles Antarctica, extending from the polar front in the north to the late winter or spring sea ice extent in the south (Diekmann, 2007).

Figure 1Location of core COR1GC in the western Indian sector of the Southern Ocean. (a) Map of the Southern Ocean and Antarctica. Key terrestrial ice and marine sediment cores mentioned in the text and in Fig. 7 are labelled: A: TN057-13-PC4 (Anderson et al., 2009; Divine et al., 2010); B: TN057-17-PC1 (Nielsen et al., 2004); C: MD07-3088 (Siani et al., 2013); D: COR1GC (this study); E: EDML (EPICA Community Members, 2006);F: Epica Dome C (EDC) (Jouzel et al., 2007). The mean location of the southern Antarctic Circumpolar Current Front (sACCF; dotted line), the polar front (PF; black line) and the Subantarctic Front (SAF; grey line) according to Orsi et al. (1995) are shown. (b) Bathymetric map of the Conrad Rise adapted from Ansorge et al. (2008) including areas of high absolute geostrophic velocity in red shading and direction of flow shown by red arrows.

Today the site is located within the Permanent Open Ocean Zone (POOZ), with the mean winter sea ice limit lying to the south at 59 S (Park et al., 1998; Parkinson and Cavalieri, 2012). The mean summer SST at 10 m depth is 3.5 C (Locarnini et al., 2013), while summer measurements show that at the Conrad Rise the surface mixed layer extends to 70–80 m depth with temperatures of 2–2.5 C and salinity of 33.90–33.95 psu (Park et al., 1998; Anilkumar et al., 2006). The surface oceanographic characteristics are likely influenced by the atmospheric climate and circulation. The core of the SWW was situated north of the Conrad Rise at approximately 52 S during the period 1979–2010 CE (Swart and Fyfe, 2012). Across the Southern Ocean positive SAM anomalies cause stronger, southward-shifted westerlies, and this also alters SST, causing cooling in the Antarctic Zone due to stronger Ekman transport of cold Antarctic surface water to the north (Lovenduski and Gruber, 2005).

The Conrad Rise is situated within the ACC and close to the polar front, which is defined by the maximum northern extent of the 2 C isotherm between 100 and 300 m. Some modelling and observations support the idea that two branches of the polar front flow to the north and south of the Conrad Rise, transporting 35 and 25 Sv of water, respectively, with low flow velocities over the rise itself (Fig. 1b; Pollard and Read, 2001; Ansorge et al., 2008; Sokolov and Rintoul, 2009). However, assessments of the polar front position based on the identification of a subsurface temperature minimum indicate that it is in fact situated to the north, passing between the Conrad Rise and Crozet Islands at approximately 50 S (Pauthenet et al., 2018). Topography has a strong influence on the position and form of the ACC in this region. At the Conrad Rise, branches of the polar front become topographically trapped and therefore have low spatial variability through time (Pollard and Read, 2001; Pauthenet et al., 2018).

3 Materials and methods
Back to toptop

3.1 Core sampling and sediment description

The 2.48 m long core KH-10-7 COR1GC was retrieved from the Conrad Rise (54.2673 S, 39.7663 E; Fig. 1; water depth of 2834 m) using a gravity corer during the expedition KH-10-7 on board the R/V Hakuho-Maru in 2010.

The sediment consisted of homogenous diatom ooze; from 2.48 to 2.2 m the sediment colour was greenish grey, while above this the sediment had a pale yellow colour. At approximately 5 cm depth there was a dark grey horizon.

3.2 Age model

Fifteen AMS radiocarbon dates were obtained on mono-specific samples of the planktic foraminifera Neogloboquadrina pachyderma sinistral. The core chronology was constructed in Bacon version 2.3.3 (Blaauw and Christen, 2011), an age–depth modelling technique that uses Bayesian statistics to reconstruct accumulation histories for deposits, by combining radiocarbon and other dates with prior information. This analysis used the Marine13 calibration dataset (Reimer et al., 2013). The shape of sedimentation rate and memory strength parameters used in the Bacon modelling were set to 1.5 and 4, respectively.

To correct the raw radiocarbon ages we used a constant reservoir age of 890±100 years recommended for this region (Bard, 1988; Butzin et al., 2005). However, we acknowledge the possibility that the reservoir age may have varied through time, particularly during the deglaciation when the ventilation of deep waters occurred. In the Atlantic sector of the Southern Ocean, the comparison of radiocarbon ages and 226Ra decay in barite suggests a change in regional reservoir ages between the early and mid- to late Holocene of ∼800 years (van Beek et al., 2002). In the Pacific sector, Holocene reservoir ages have been estimated in the Ross Sea through a comparison of 14C and U−Th dating of corals, which showed limited changes in reservoir age during the last 6 kyr (Hall et al., 2010). In the east Pacific sector (core MD07-3088, Fig. 1), a comparison between radiocarbon and tephra dating of marine sediments suggests that at 14 and 11.5 ka the reservoir age was 975–920 years, while in the mid- to late Holocene it was 800 years (Siani et al., 2013), indicating a change of just 100–175 years. While this study was located north of the Subantarctic Front in a different oceanographic setting to the Conrad Rise, the MD07-3088 study site is bathed in waters originating from the Southern Ocean that are deflected north along the coast of South America, supporting the idea that changes here reflect Southern Ocean reservoir age variability. Given the information delivered by these remote studies and the lack of evidence from the Indian sector of the Southern Ocean, we believe it is more sensible to apply a constant reservoir age correction but acknowledge that there is greater age uncertainty in the late deglacial section of the COR1GC reconstruction.

3.3 Diatom analysis

Diatom analysis was conducted at ∼1 cm intervals. Slide preparation and counting of over 300 diatom valves per sample followed standard methods (Crosta and Koç, 2007). Summer SST (at 10 m depth) and winter sea ice concentration (WSIC) were estimated using the modern analogue technique (MAT) applied to diatom assemblages (Crosta et al., 2020). This method uses a modern database composed of 249 surface sediment samples. Modern summer (January–February–March) SSTs from 10 m depth were interpolated on a 1 grid from the World Ocean Atlas 2013 (Locarnini et al., 2013) using Ocean Data View (Schlitzer, 2014). Winter sea ice concentration (WSIC), represented by the month of September when sea ice is at its maximum extent around Antarctica (Parkinson, 2019), was interpolated on a 1 grid from the numerical atlas of Schweitzer (1995). The MAT was implemented from the “bioindic” package (Guiot and de Vernal, 2011) built on the R platform (, last access: 1 August 2019) and used the relative abundances of 32 diatom species and the chord distance, as a similarity metric, to select the five most similar modern analogues. The threshold above which modern analogues are supposed to be too dissimilar to the fossil assemblage is fixed at the first quartile of random distances on the validation or modern dataset. Quantitative estimates of summer SST and WSIC are a similarity-weighted mean of the SST or WSIC associated with the selected modern analogues (Guiot et al., 1993). This method yields an R2 of 0.96 and a root mean square error of prediction (RMSEP) of ∼1C for summer SST and an R2 of 0.92 and an RMSEP of ∼10 % for WSIC.

The relative abundance of different diatom species and predictions of SST and WSIC may be altered by differential dissolution and preservation between species. While temperature is the dominant control on diatom species assemblages (Pichon et al., 1992; Esper et al., 2010), the dissolution of weakly silicified species, which are often those from cold waters, could potentially result in warmer reconstructed temperatures or an underestimation of sea ice extent (Xiao et al., 2016). However, as the dissolution of species is strongest in the Subantarctic Zone (SAZ) and sea ice zone of the Southern Ocean (Esper et al., 2010) and a dissolution minimum (of 10 %) is centred on 55 S (Pichon et al., 1992), the Conrad Rise is a region likely to be less influenced by diatom dissolution than locations at both lower and higher latitudes in the Southern Ocean. Furthermore, several studies have shown that the surface water gradients in most abundant diatom species are preserved in the surface sediments (Gersonde and Zielinski, 2000; Armand et al., 2008a, b; Cárdenas et al., 2019). The strong association between ocean surface conditions and diatom assemblages within the sediment is supported by several statistical techniques that have been used to quantitatively infer past surface conditions using diatom analysis (e.g. Gersonde and Zielinski, 2000; Crosta et al., 2004; Esper and Gersonde, 2014).

The total diatom flux (DF), or diatom accumulation rate, was calculated to provide an approximation of productivity through time, following the method described in Romero et al. (2015). To calculate the DF, the number of valves per gram of dry sediment (diatom absolute abundance) was multiplied by the bulk density of wet sediment and the sedimentation rate. The resulting value was then divided by 2 to represent the number of diatoms rather than valves. The sedimentation rates and diatom abundances were calculated for core COR1GC; however, because density was not measured on core COR1GC the wet bulk density measurements for core COR1bPC were instead used (Oiwane et al., 2014), as this core was recovered from the same location. The density measurements were transferred from COR1bPC to COR1GC using the chronologies of the two cores.

3.4 SiZer analysis

SiZer (Significance of Zero Crossings of the Derivative) (Chaudhuri and Marron, 1999) is a scale–space technique that was applied to explore statistically significant features in the reconstructed SST. A key idea in SiZer is that significant features are found at different timescales, that is, at different levels of data smoothing. This makes it particularly useful in palaeoclimate studies since the salient features in a time series may depend heavily on the time horizon on which it is analysed. This method has been used previously in a number of palaeoclimate studies (e.g. Divine et al., 2010).

In SiZer the observed data are viewed at varying levels of resolution, while the notion of scale is controlled through the bandwidth h in the local linear kernel estimator. For each scale h and time t of the signal, a test is performed to see whether the smoothed data series has a local derivative significantly different from 0, in other words, to see if the slope at a specific time point for a given scale is significantly different from 0. The Gaussian kernel estimator embedded in SiZer does not require an analysed time series to be evenly spaced; the method is therefore applied to the data directly without any prior resampling. SiZer visualises the output of the analysis in a feature map where the results are displayed as a function of time and scale.

3.5 Spectral analysis

Many geophysical time series have distinctive red noise characteristics that can be modelled by a first-order autoregressive (AR1) process. REDFIT spectral analysis (Schulz and Mudelsee, 2002) was carried out on the Holocene section of the SST reconstruction to establish whether quasi-periodic variability existed. This method of spectral analysis was selected as it has been designed for unevenly spaced palaeoclimate data, therefore avoiding the need to resample and interpolate the data, which can introduce bias in the inferred spectrum (Schulz and Mudelsee, 2002). The analysis was conducted using the REDFIT tool integrated in the PAST 3.25 software (Hammer et al., 2001) using the Welch window with two overlapping segments (Schulz and Mudelsee, 2002). The appropriateness of the AR1 model to describe the analysed data was tested using a nonparametric run test (Bendat and Piersol, 1986) embedded in the package.

4 Results
Back to toptop

According to the constructed chronology the core covers the period from 14.2 to 1 ka and has a mean sedimentation rate of ∼20 cm kyr−1 (Table 1; Fig. 2). The sedimentation rate decreased from 14 to 12 ka when it stabilised at its lowest values until 10 ka. It subsequently increased over the course of the Holocene to reach values >20 cm kyr−1 during the late Holocene.

Table 1AMS radiocarbon dates and calibrated age estimates.

Download Print Version | Download XLSX

Figure 2(a) Age–depth model for core COR1GC. The age model is constrained by 15 calibrated radiocarbon dates, the continuous line is the median age and shaded area represents the 95 % confidence interval of the modelled core chronology. (b) Sedimentation rate estimates for COR1GC based on median modelled ages for AMS 14C dates.


The diatom assemblage of COR1GC is dominated by pelagic open-ocean taxa, particularly Fragilariopsis kerguelensis and Thalassiosira lentiginosa (Zielinski and Gersonde, 1997; Crosta et al., 2005), with accompanying species typical of the Polar Front Zone (PFZ) and POOZ (Fig. 3). Species from the PFZ (F. kerguelensis, Thalassionema nitzschioides var. lanceolata) and SAZ (Shionodiscus oestrupii, previously known as Thalassiosira oestrupii; Crosta et al., 2005; Romero et al., 2005) were more abundant in the interval 12–9.5 ka, with lowest relative abundances at 14.2–12 ka and 4–1 ka. Conversely, POOZ species abundances (especially Thalassiothrix antarctica, Thalassiosira gracilis, Fragilariopsis separanda and Fragilariopsis rhombica; Zielinski and Gersonde, 1997; Crosta et al., 2005) were lower at 14.2–9.5 ka, after which point they gradually increased in proportion over the Holocene. Diatom species associated with sea ice, including Fragilariopsis ritscheri, Fragilariopsis curta, Fragilariopsis cylindrus, Porosira glacialis, Fragilariopsis obliquecostata and Thalassiosira tumida (Zielinski and Gersonde, 1997; Armand et al., 2005), had individually low average abundances of ≤0.3 %. These species have been grouped together (Fig. 3). The sea ice species were present in low abundances through the record but were slightly more abundant between 14.2 and 12.6, 7.6 and 6.2, and 3.9 and 1 ka and lower from 12.6 to 7.6 and 6.2 to 3.9 ka.

Figure 3COR1GC diatom abundances as percentages of total sample and the total diatom flux. Note that the species are plotted on different x-axis scales.


The diatom flux record of COR1GC (Fig. 3) shows a range from ∼2000 to 14 000×106 diatoms cm2 kyr−1. From 14.2 to 12 ka the diatom flux decreased (from ∼12 000 to ×106 diatoms cm2 kyr−1). Concentrations remained low between 12 and 10.3 ka followed by higher fluxes with an average of 7200×106 diatoms cm2 kyr−1 between 10.3 and 5.2 ka. After 5.2 ka fluxes decreased, with an average of 4500×106 diatoms cm2 kyr−1 between 5.2 and 1 ka.

The SST reconstruction for COR1GC shows that during the 14.2–12.9 ka period the average temperature was 2.9 C with a range of 1.8–4.1 C and standard deviation of 0.6 C (Fig. 4), which was followed by a phase of warming between 12.9 and 11.6 ka when temperatures increased from 2.1 to 5 C. Between 11.6 and 8.7 ka SST was on average 4.3 C, with a range of 3.3–5 C and a standard deviation of 0.4 C. After 8.7 ka the SST record shows a slightly lower average of 3.9 C, a range of 2.2–5.2 C, a standard deviation of 0.7 C and a slight cooling trend. The record shows that during the slightly warmer interval at 11.6–8.7 ka the temperature varied less than after 8.7 ka, indicating a potentially more stable climate during the early Holocene and a more variable climate during the mid- to late Holocene. SiZer analysis flags three events as statistically significant: the two short-term cooling events at ∼8.2 ka and ∼2.2 ka and the warming trend during the deglaciation, the latter of which is marked as significant when identifying variability over a broad range of timescales (Fig. 5). Conversely, SiZer did not find the cooling through the Holocene to be significant. Spectral analysis conducted using the REDFIT method shows that during the Holocene there were quasi-periodic oscillations of 220 and 260 years, significant above the 95 % false alarm level (Fig. 6).

Figure 4Reconstructed mean summer (January–March) sea surface temperature and mean winter (September) sea ice concentrations for core COR1GC. Key periods discussed are highlighted, along with significant centennial cool events at ca. 8.2 and 2.2 ka, as identified by the SiZer analysis (Fig. 5). The modern average summer SST at 10 m depth is 3.5 C (Locarnini et al., 2013) marked on the graph, while the modern winter sea ice concentration (WSIC) is marked as 0 as the sea ice limit is at 59 S (Park et al., 1998; Parkinson and Cavalieri, 2012).


Figure 5SiZer analysis of the COR1GC SST reconstruction. (a) Family plot of the reconstructed SST. The green dots represent the raw data. Blue lines show the family of smoothings obtained by the local linear kernel estimator for a range of bandwidths (scales h). (b) A SiZer map, given as a function of location (time t) and scale h. A significantly positive (increased SST) derivative is flagged as red, while a significantly negative (decreased SST) derivative is flagged as blue. The colour purple is used at locations where the derivative is not significantly different from 0. The colour dark grey is used to indicate that too few data points are available to make a correct inference. The distance between the two dotted lines in the cone-shaped curve for a horizontal line in the SiZer plot can be interpreted as the scale for that level of resolution.


Figure 6REDFIT spectral analysis of the COR1GC SST reconstruction during the period 12–1 ka. The 95 % significance level (Chi2) is shown by the dashed black line. The time series is fitted to an AR(1) red noise model (dashed red line).


Summer sea ice did not reach the Conrad Rise through the last 14.2 kyr (reconstruction not shown) and winter sea ice concentrations were very low and below the RMSEP throughout the period (Fig. 4). These data however suggest slightly more extensive sea ice cover between 14 and 12 ka and after 8 ka, and absent winter sea ice between 12.2 and 8 ka when SSTs were ∼1C higher than today. Relatively low values of the inferred winter sea ice concentration suggest only sporadic and relatively short-term (perhaps decadal- to multi-decadal-scale) expansions of winter sea ice to the coring location throughout the mid- to late Holocene.

5 Discussion
Back to toptop

5.1 Southern high-latitude temperature

The millennial trends in the COR1GC temperature record are similar to those of SST records from the Atlantic sector of the Southern Ocean (Nielsen et al., 2004; Anderson et al., 2009; Divine et al., 2010; Xiao et al., 2016; Fig. 7a, b) within the age uncertainty of each record. Many of these show cooler temperatures at ∼14.2–12.9 ka during the ACR (14.7–12.7 ka; Stenni et al., 2011), gradual warming during the Younger Dryas, and high temperatures between 11.6 and 8.7 ka during the early Holocene, followed by a cooling trend thereafter. Another SST reconstruction from the Conrad Rise (PS2606-6; Fig, 7c; Xiao et al., 2016), with lower temporal resolution than the COR1GC record, shows warm temperatures during the ACR and the early Holocene along with a slight cooling during the Younger Dryas, at odds with the other records from the POOZ. We suggest that the opposite trend in the PS2606-6 record during the late deglaciation could be related to either chronological issues, as it is the only SST record that is out of phase during this period, or the result of different reconstruction techniques. Furthermore, while most records show cooling over the Holocene, either gradually or as an abrupt cooling at the end of the early Holocene (Xiao et al., 2016), there are some differences between records. For example, the cooling in the COR1GC is small compared to other records as the change in mean temperature from the early to mid- to late Holocene was ∼0.4C, although superimposed on this long-term trend were centennial events with 1–2 C cooler temperatures during the mid- to late Holocene (Fig. 4). We note, on the other hand, that core TN057-17 close to the polar front in the Atlantic sector shows a warming trend since ∼4 ka (Fig. 7; Nielsen et al., 2004), in contrast to a nearby low-resolution record (cores PS1654/OPD1093), which indicates that SSTs were relatively stable over the late Holocene (Xiao et al., 2016). Despite many records having similar long-term patterns, these results therefore suggest some regional heterogeneity of temperature evolution at the centennial to multi-millennial scales through the late deglacial and Holocene.

Figure 7Reconstructions reflecting temperature since 14.2 ka from the Southern Ocean and Antarctica at locations shown in Fig. 1. (a) TN057-17-PC1 diatom-based SST reconstruction (Nielsen et al., 2004); (b) TN057-13-PC4 diatom-based SST reconstruction (Anderson et al., 2009; Divine et al., 2010); (c) PS2606-6 diatom-based SST reconstruction from the Conrad Rise (Xiao et al., 2016); (d) COR1GC SST reconstruction (this study) from the Conrad Rise; (e) EDML ice core δ18O reconstruction (EPICA Community Members, 2006); (f) EDC ice core temperature reconstruction (Jouzel et al., 2007).


The records from the Bouvet Island region, in the eastern Atlantic sector of the Southern Ocean, and the Conrad Rise (Fig. 7a–d; cores TN057-17-PC1, TN057-13-PC4 and PS2606-6; Xiao et al., 2016) show that the magnitude of warming between the ACR and early Holocene was 1–2 C. Nevertheless, there are variations in the exact magnitudes of the temperature fluctuations between the different records. Some previous records suggested that the SST during the ACR and the mid- to late Holocene were similar (Nielsen et al., 2004; Anderson et al., 2009; Divine et al., 2010; Xiao et al., 2016; Fig. 7a–c). Our new record from core COR1GC however shows SSTs were 1 C lower during the ACR compared to the mid- to late Holocene, in agreement with reconstructions of air temperature from ice cores (Fig. 7e, f). The contrasting findings between SST reconstructions may be explained by the reconstructed temperatures being close to the prediction error of SST transfer functions, which is ∼1C in this study and 0.86 C in Xiao et al. (2016).

The millennial temperature evolution of the COR1GC record during the deglacial and early Holocene closely resembles the patterns of atmospheric temperature and sea ice extent, reconstructed from Antarctic ice core water isotope records and sodium flux records, respectively (Fig. 7e, f; EPICA Community Members, 2006; Fischer et al., 2007; Jouzel et al., 2007; Stenni et al. 2011). This is particularly evident when comparing the smoothed and normalised reconstructions of temperature (δ18O) and regional sea ice extent (ssNa) from the closest ice core, EDML, with the COR1GC SST record (Supplement). The start of the COR1GC record at 14.2 ka is after the beginning of the ACR (at 14.7 ka; Stenni et al., 2011) meaning that the timing of the onset of the ACR at the Conrad Rise cannot be established. However despite some degree of dating uncertainty the warming from 12.9 to 11.6 ka appears closely aligned with the timing of warming identified in the ice core records at 12.7–11.9 ka (Stenni et al., 2011). A good synchronicity between Southern Ocean SST and Antarctic temperature evolution over the deglacial was also shown by a planktic foraminifera-based SST reconstruction from core MD07-3088 from the eastern Pacific sector, which like the COR1GC record showed cold temperatures during the ACR, followed by gradual warming between ca. 12.7 and 11.5 ka, leading to slightly warmer conditions through the Holocene (Siani et al., 2013; Haddam et al., 2018). The MD07-3088 record was constrained by both radiocarbon dates and tephra horizons minimising dating uncertainties (Siani et al., 2013), which supports the chronological framework of core COR1GC and the new SST record. The marginally warmer early Holocene temperatures that followed at 12–8.7 ka are observed widely in ice core, terrestrial and coastal records from across Antarctica (Stenni et al., 2011; Verleyen et al., 2011; Fig. 7e, f). Together this evidence suggests that the sea surface and atmospheric temperatures across the southern high latitudes varied synchronously during the deglacial to Holocene transition.

5.2 Causes of temperature variability

The SST variations may have been due to a number of interrelated factors, including in response to atmospheric temperature changes, variations in the AMOC, changes in the strength and position of the SWW, and shifts in the location of the ACC fronts relative to the core site. Here these causes will be discussed in relation to the Younger Dryas and Holocene climate patterns, with the suggested causes summarised in Table 2.

Table 2Summary table of the suggested forcings of the core COR1GC SST reconstruction, as proposed in Sect. 5.2.

Download Print Version | Download XLSX

5.2.1 Younger Dryas

Southern high-latitude warming during the Younger Dryas is widely considered to be the result of a slowdown in the AMOC (McManus et al., 2004; Fig. 8a), which caused reduced northward heat transport and the accumulation of heat in the Southern Hemisphere (Broecker, 1998; Stocker and Johnsen, 2003). It is considered that this hemispheric temperature change caused a southward shift of the ITCZ (Fig. 8b; Hughen et al., 1996), which in turn caused poleward-shifted and strengthened SWW (Denton et al., 2010; Pedro et al., 2018). Evidence from the Southern Ocean supports the idea that the SWWs were stronger over a longer interval from 14 to 11 ka but only shifted southward to a more poleward position at 12.5 ka (Fig. 8c; Fletcher and Moreno, 2011; Saunders et al., 2018).

Figure 8Summary of key reconstructions showing the late deglacial to Holocene variations in the AMOC, atmospheric circulation and Southern Ocean. (a) AMOC reconstruction based on Pa231∕Th230 from a core from the western subtropical North Atlantic (McManus et al., 2004); (b) greyscale measurements from cores from the Cariaco Basin related to productivity, interpreted as reflecting ITCZ position (Hughen et al., 1996); (c) summary of reconstructed SWW strength based on a collation of palaeoclimate data from the southern mid-latitudes (41–52 S; Fletcher and Moreno, 2011); (d) upwelling in the Atlantic sector of Southern Ocean (core TN057-13PC4) based on opal flux (Anderson et al., 2009; Divine et al., 2010); (e) productivity at the Conrad Rise inferred from the diatom flux in COR1GC; * measured in units: 109 diatoms cm−2 kyr−1; (f) EDC CO2 reconstruction (Monnin et al., 2001); (g) COR1GC reconstruction of SST (this study).


The changes to the SWW during the Younger Dryas are thought to have warmed the Southern Ocean in a number of ways, including through the southward migration of the polar front, causing local warming in regions currently to the south of the polar front (e.g. Barker et al., 2009) and through strengthening of the eddy transport of heat across the ACC (e.g. Pedro et al., 2018) causing a more widespread warming of the POOZ. Although there is a lack of palaeo-proxy evidence for stronger eddy transport of heat across the ACC, this mechanism has been identified in both model simulations and modern observations (Meredith and Hogg, 2006; Screen et al., 2009). It is considered that this initial warming of the Southern Ocean resulted in reduced sea ice in the Indian and Atlantic sectors at 12 ka (Bianchi and Gersonde, 2004; Iizuka et al., 2008), promoting further ocean warming in the newly ice-free region (Denton et al., 2010; Xiao et al., 2016; Pedro et al., 2018) and warming of the lower troposphere and Antarctica by ice-albedo feedbacks (Pedro et al., 2018). Greater upwelling, caused by stronger winds and a reduction in sea ice extent, has been shown by higher opal deposition to the south of the polar front in the Southern Ocean (Atlantic, Indian and Pacific sectors) through the period 12.7–11.5 ka (Fig. 8d; Dezileau et al., 2003; Anderson et al., 2009; Siani et al., 2013). The combination of greater upwelling and reduced sea ice allowed more outgassing of CO2 to the atmosphere (Fig. 8f), causing warming globally (Monnin et al., 2001; Anderson et al., 2009; Clark et al., 2012; Saunders et al., 2018).

The warming and reduced sea ice shown by the COR1GC reconstructions during the Younger Dryas (Figs. 3 and 4) coincide with many of these changes, including the AMOC weakening, strengthened and southward-shifted SWW, increased upwelling, reduced sea ice, rising CO2, and warming over Antarctica (Fig. 8; Supplement). However, there is no evidence of increased productivity and therefore upwelling at the Conrad Rise, as diatom fluxes instead decreased through this period (Fig. 8e). Although the marine records have sources of age uncertainty, the similarity in the timing of changes in marine and terrestrial records supports warming being a result of processes connected with the interhemispheric oceanic and atmospheric reorganisation (Fig. 8).

One challenge is to establish the relative importance of each of the outlined causes of warming because although changes occurred synchronously they may not have all contributed equally to the warming of the Southern Ocean. For example, the Conrad Rise is located close to the modern polar front and therefore a southward shift of the SWW during the Younger Dryas may have pushed the polar front to the south, leading to warming. Indeed, a southward shift is indicated by an increase in the polar front species Thalassionema nitzschioides var. lanceolata at ca. 12 ka (Fig. 3). Although observations and models show that the fronts of the polar front are topographically constrained in this region (Pollard and Read, 2001; Graham et al., 2012; Pauthenet et al., 2018), it is possible that larger magnitude changes in the SWW were able to divert the polar front to the south of the Conrad Rise. The relative importance of CO2 in contributing to the warming of Antarctica and the Southern Ocean, which resulted from the above-described changes in circulation and Southern Ocean conditions, is also not clear. Principle component analysis on a synthesis of global temperature reconstructions spanning the last deglaciation indicate that temperatures globally were primarily controlled by greenhouse gas concentrations (Clark et al., 2012) and temperature over Antarctica rose synchronously with CO2 increases during the Younger Dryas (Parrenin et al., 2013). The synchronicity between the CO2 record and some of the temperature records (COR1GC, MD07-3088 and ice core records) over the deglaciation, suggests that the CO2 increase, following the initial warming and sea ice loss, contributed to the amplitude of the deglacial warming and subsequent early Holocene temperatures through a strong positive feedback on the climate system (Past Interglacials Working Group, 2016).

5.2.2 Holocene

The warmth during the early Holocene (12–8.7 ka) in the southern high latitudes following the Younger Dryas warming can be attributed to high spring insolation causing a longer summer season (Shevenell et al., 2011; Etourneau et al., 2013), which would have enhanced the duration of summer warming and reduced the duration of winter cooling of the open ocean. Insolation-driven changes to the SWW may have also influenced the SST. Reconstructions support the idea that the SWW was weaker generally between 11 and 7 ka (Fletcher and Moreno, 2011; Saunders et al., 2018; Fig. 8c), shifted southward (Lamy et al., 2010) and that atmospheric circulation was more meridional (Mayewski et al., 2013). The more meridional circulation may have led to a greater poleward penetration of warm air masses and therefore warming over the Southern Ocean. A general weakening of the SWW may have reduced the amount of northward Ekman transport of cold water (e.g. Hall and Visbeck, 2002; Lovenduski and Gruber, 2005), again leading to warming further north. However at the latitude of the Conrad Rise (54 S), it is likely the SWW strengthened, as the southward SWW caused stronger winds in South America at ca. 53 S (Lamy et al., 2010), and this may have potentially increased heat loss from the ocean surface, damping the warming.

At ca. 8.2 ka the SiZer analysis indicated that there was a significant cool interval in the COR1GC record (Fig. 5), while the diatom record by Katsuki et al. (2012) from the Conrad Rise also showed a peak in cold species at 8 ka. The duration and timing of the SST cooling observed in the COR1GC record (from 8.5 to 7.9 ka) coincides with an AMOC reduction and North Atlantic cooling associated with the 8.2 event (Ellison et al., 2006). The AMOC reduction should have caused a warming of the Southern Ocean via the bipolar seesaw. However a modelling study suggests that cooler North Atlantic Deep Water (NADW) was formed in the North Atlantic, due to the retreat of the Laurentide Ice Sheet, and this upwelled in the Southern Ocean causing a cooling in the POOZ between 9 and 7 ka (Renssen et al., 2010). The 8.2 event has not generally been observed in records from the Southern Hemisphere (Alley and Ágústsdóttir, 2005). While there are some indications of a longer period with lower temperatures in Antarctic ice core records at this time (Fig. 7e, f; EPICA Community Members, 2006; Jouzel et al., 2007) and evidence of colder SST at the Antarctic Peninsula (Etourneau et al., 2019), the few other high-resolution SST records from the open ocean do not show a cool event (Fig. 7; TN057-13PC4; TN057-17TC; Nielsen et al., 2004; Anderson et al., 2009; Divine et al., 2010). Therefore although there appears to have been a cooling on the Conrad Rise at this time, it cannot be reliably associated with the 8.2 ka event in the Northern Hemisphere without additional evidence from other high-resolution marine sediment records.

From the mid- to late Holocene the spring southern high-latitude insolation decreased (Laskar et al., 2004), causing a shorter summer duration and cooling, which is coherent with the northward shift of the SWW (Lamy et al., 2001). Together this can explain the slight cooling in the COR1GC record, the increased frequency of cold events and increase in sea ice (Fig. 4). The winter and summer latitudinal insolation gradient also increased (Laskar et al., 2004; Divine et al., 2010), causing the SWW to strengthen (e.g. Saunders et al., 2018), which may have caused cooling by promoting enhanced northward Ekman transport of cooler water from the south towards the Conrad Rise (e.g. Hall and Visbeck, 2002; Lovenduski and Gruber, 2005). The strengthened winds may also have increased upwelling and therefore productivity; however, at the Conrad Rise this does not seem to have occurred as the estimated diatom flux decreased, particularly after 5 ka (Fig. 3). Potentially, a strengthening of the SWW to the north, rather than over the Conrad Rise, may explain why upwelling and productivity did not increase, which agrees with evidence from the Atlantic sector that during the mid- to late Holocene productivity decreased south of the polar front (Anderson et al., 2009).

The COR1GC SST record also shows increasingly variable SSTs through the mid- to late Holocene after 8 ka (Fig. 4). Another record from the Conrad Rise found an increased frequency of cold events from 5.5 ka onwards (Katsuki et al., 2012), and the closest ice core record, EDML, shows increasing variability in the δ18O record from ca. 6 ka onwards (EPICA Community Members, 2006; Fischer et al., 2007; Fig. 7e). Spectral analysis conducted on the COR1GC SST record identifies significant quasi-periodicities of ∼220 to 260 years (Fig. 6) during the last 12 000 years. Similar quasi-periodicities of 200–300 years have previously been observed in Southern Ocean palaeoceanographic records, including those reflecting SST, sea ice and productivity changes (Leventer et al., 1996; Bárcena et al., 1998; Nielsen et al., 2004; Crosta et al., 2007). Furthermore, quasi-periodicities of 250 years have been identified in a record of westerly wind intensity from the Falkland Islands (Turney et al., 2016), and multi-centennial variations were similarly identified in a 3000-year reconstruction of the SAM (Moreno et al., 2014). These similar patterns of cyclicity in records reflecting both ocean and atmospheric conditions suggest that they were caused by (1) internal variability and ocean–atmosphere coupling or (2) ocean and atmosphere responses to a shared forcing. In support of the first hypothesis, model simulations have shown similar 200–300-year oscillations in Southern Ocean temperature and atmospheric circulation (Latif et al., 2013). In a simulation by the Kiel Climate Model, centennial episodes of strong (weak) deep convection in the Weddell Sea caused warm (cold) SSTs and reduced (expanded) sea ice, which altered atmospheric pressure patterns and led to weakened (strengthened) atmospheric circulation (Latif et al., 2013). In this case, internal oscillations in the climate system have a hemispheric impact, which may explain the similar variability observed in the COR1GC record and other records. In support of the second hypothesis, external forcing, such as cycles in solar activity, may have caused these changes (Leventer et al., 1996; Bárcena et al., 1998; Crosta et al., 2007; Turney et al., 2016), with changes in solar irradiance directly altering ocean SST, and potentially regulating the SAM and associated changes in atmospheric and oceanic circulation, as suggested by some modern evidence that shows a connection between the SAM and the 11-year solar cycle (e.g. Kuroda and Kodera, 2005; Kuroda, 2018). While it is yet not possible to establish the cause of this cycle, the enhanced variability in the mid- to late Holocene appears to reflect the establishment of modern coupled ocean–atmosphere relationships over the Southern Ocean.

6 Conclusions
Back to toptop

The COR1GC reconstruction of SST variability at the Conrad Rise in the west Indian Sector of the Southern Ocean spans the period from 14.2 to 1 ka with an average resolution of ∼60 years. SSTs were cool (∼2.9C) from 14.2 to 12.9 ka during the ACR before temperatures increased to 5 C between 12.9 and 11.6 ka during the Younger Dryas. During the Holocene SSTs were more stable and marginally warmer between 12 and 9 ka before a significant, centennial-scale cooling event occurred at ca. 8.2 ka, followed by greater SST variability through the period 8–1 ka.

The timing and duration of warming in the COR1GC record at 12.9–11.6 ka reflects warming also shown by a precisely dated marine record of SST from the Pacific (Siani et al., 2013) and Antarctic ice core temperature records, supporting the idea that atmospheric and ocean temperatures varied in phase. The COR1GC warming also coincides with palaeo-evidence for a weakened AMOC, rising CO2 levels, greater upwelling in the Atlantic sector of the Southern Ocean, southward displacement and strengthening of the SWW, and reduced sea ice, each of which may have contributed to warming in the Southern Ocean. The results support the idea that the warming was initiated by an interhemispheric oceanic and atmospheric reorganisation during the Younger Dryas but do not allow an attribution of the relative importance of each process to warming.

For the Holocene period, the most striking change in the COR1GC record is the switch from stable and slightly warmer conditions to variable and slightly cooler conditions at 8.7 ka, potentially linked with a global reorganisation of the climate system. It is suggested that the early Holocene warmth may have resulted from higher spring insolation increasing the duration of heat accumulation in the Southern Ocean during the spring–summer season or as a result of changes to latitudinal heat transport as the SWW weakened. During the mid- to late Holocene the enhanced variability may represent the establishment of modern coupled atmosphere–ocean relationships, whereby the SAM and SSTs in the Southern Ocean vary in phase over centennial timescales. This is supported by the identification of a quasi-periodic oscillation of ∼200–260 years in the Holocene section of the COR1GC SST reconstruction, which mirrors similar centennial variability identified in reconstructions of the SAM and westerly winds.

Data availability
Back to toptop
Data availability. 

The complete dataset can be obtained at (Orme et al., 2020a) and (Orme et al., 2020b).

Back to toptop

The supplement related to this article is available online at:

Author contributions
Back to toptop
Author contributions. 

LCO performed the diatom analysis and wrote the paper. XC produced the SST record and did training and writing and editing of the paper. DVD did SiZer analysis. AM, DVD, KH, EI, and RM edited the paper. AM, KH, EI, XC, RM, and DVD conceptualised and supervised the project. LW and OT conducted laboratory analysis. MI provided core material.

Competing interests
Back to toptop
Competing interests. 

The authors declare that they have no conflict of interest.

Back to toptop

This work was funded by the Research Council of Norway and the Ministry of Earth Science, Earth System Science Organization, through the OCTEL project. We would like to thank Svetlana Divina and Linda Rossignol-Malaizé for preparation of the foraminifera samples for radiocarbon dating. Gravity core sampling and subsample preparation were partly supported by JSPS KAKENHI. Rahul Mohan would like to thank the Director, NCPOR, India, for the support to the project, and this is NCPOR Contribution No. J-18/2020-21.

Financial support
Back to toptop
Financial support. 

This research has been supported by the Research Council of Norway (grant no. 248776/E10), the Ministry of Earth Science, Earth System Science Organization (grant no. MoES/Indo-Nor/PS-2/2015), and the JSPS KAKENHI (grant nos. 23244102, 17H06318).

Review statement
Back to toptop
Review statement. 

This paper was edited by Erin McClymont and reviewed by two anonymous referees.

Back to toptop

Alley, R. B. and Ágústsdóttir, A. M.: The 8k event: Cause and consequences of a major Holocene abrupt climate change, Quaternary Sci. Rev., 24, 1123–1149,, 2005. 

Anderson, R. F., Ali, S., Bradtmiller, L. I., Nielsen, S. H. H., Fleisher, M. Q., Anderson, B. E., and Burckle, L. H.: Wind-driven upwelling in the southern ocean and the deglacial rise in atmospheric CO2, Science, 323, 1443–1448,, 2009. 

Anilkumar, N., Luis, A. J., Somayajulu, Y. K., Babu, V. R., Dash, M. K., Pednekar, S. M., Babu, K. N., Sudhakar, M., and Pandey, P. C.: Fronts, water masses and heat content variability in the Western Indian sector of the Southern Ocean during austral summer 2004, J. Marine Syst., 63, 20–34,, 2006. 

Ansorge, I. J., Roman, R., Durgadoo, J. V., Ryan, P. G., Dlamini, L., Gebhardt, Z., Rainier, S., Smith, M., Mtonsti, T., and Lutjeharms, J. R. E.: The first oceanographic survey of the Conrad Rise, S. Afr. J. Sci., 104, 333–336, 2008. 

Armand, L. K., Crosta, X., Romero, O., and Pichon, J.-J.: The Biogeography of Major Diatom Taxa in Southern Ocean Sediments: 1. Sea Ice Related Species, Palaeogeogr. Palaeocl., 223, 93-126,, 2005. 

Armand, L. K., Cornet-Barthaux, V., Mosseri, J., and Quéguiner, B.: Late summer diatom biomass and community structure on and around the naturally iron-fertilised Kerguelen Plateau in the Southern Ocean, Deep Sea Res. II, 55, 653–676,, 2008a. 

Armand, L. K., Crosta, X., Quéguiner, B., Mosseri, J., and Garcia, N.: Diatoms preserved in surface sediments of the northeastern Kerguelen Plateau, Deep Sea Res. II, 55, 677–692,, 2008b. 

Bárcena, M. A., Gersonde, R., Ledesma, S., Fabrés, J., Calafat, A. M., Canals, M., Sierro, F. J., and Flores, J. A.: Record of Holocene glacial oscillations in Bransfield Basin as revealed by siliceous microfossil assemblages, Antarct. Sci., 10, 269–285,, 1998. 

Bard, E.: Correction of accelerator mass spectrometry 14C ages measured in planktonic foraminifera: paleoceanographic implications, Paleoceanogr. Paleocl., 3, 635–645, 1988. 

Barker, S., Vautravers, M. J., Pike, J., Hall, I. R., Broecker, W. S., Diz, P., and Knorr, G.: Interhemispheric Atlantic seesaw response during the last deglaciation, Nature, 457, 1097–1102,, 2009. 

Bendat, J. S. and Piersol, A. G.: Decomposition of wave forces into linear and non-linear components, J. Sound Vib., 106, 391–408,, 1986. 

Bianchi, C. and Gersonde, R.: Climate evolution at the last deglaciation: The role of the Southern Ocean, Earth Planet. Sci. Lett., 228, 407–424,, 2004. 

Blaauw, M. and Christen, J. A.: Flexible paleoclimate age-depth models using an autoregressive gamma process, Bayesian Anal., 6, 457–474,, 2011. 

Broecker, W. S.: Paleocean circulation during the last deglaciation: a bipolar seesaw?, Paleoceanogr., 13, 119–121,, 1998. 

Butzin, M., Prange, M., and Lohmann, G.: Radiocarbon simulations for the glacial ocean: the effects of wind stress, Southern Ocean sea ice and Heinrich events, Earth Planet. Sci. Lett., 235, 45–61,, 2005. 

Chaudhuri, P. and Marron, J. S.: SiZer for exploration of structures in curves, J. Am. Stat. Assoc., 94, 807–823,, 1999. 

Cárdenas, P., Lange, C. B., Vernet, M., Esper, O., Srain, B., Vorrath, M.-E., Ehrhardt, S., Müller, J., Kuhn, G., Arz, H. W., Lembke-Jene, L., and Lamy, F.: Biogeochemical proxies and diatoms in surface sediments across the Drake Passage reflect oceanic domains and frontal systems in the region, Prog. Oceano., 174, 72–88,, 2019. 

Clark, P. U., Shakun, J. D., Baker, P. A., Bartlein, P. J., Brewer, S., Brook, E., Carlson, A. E., Cheng, H., Kaufman, D. S., Liu, Z., Marchitto, T. M., Mix, A. C., Morrill, C., Otto-Bliesner, B. L., Pahnke, K., Russell, J. M., Whitlock, C., Adkins, J. F., Blois, J. L., Clark, J., Colman, S. M., Curry, W. B., Flower, B. P., He, F., Johnson, T. C., Lynch-Stieglitz, J., Markgraf, V., McManus, J., Mitrovica, J. X., Moreno, P. I., and Williams, J. W.: Global climate evolution during the last deglaciation, P. Natl. Acad. Sci. USA, 109, E1134–E1142,, 2012. 

Crosta, X. and Koç, N.: Chapter eight diatoms: From micropaleontology to isotope geochemistry, Develop. Mar. Geol., 1, 327–369,, 2007. 

Crosta, X., Sturm, A., Armand, L., and Pichon, J. J.: Late Quaternary sea ice history in the Indian sector of the Southern Ocean as recorded by diatom assemblages, Mar. Micropaleontol., 50, 209–223,, 2004. 

Crosta, X., Romero, O., Armand, L. K., and Pichon, J.-J.: The Biogeography of Major Diatom Taxa in Southern Ocean Sediments: 2. Open Ocean Related Species, Palaeogeogr. Palaeocl., 223, 66–92,, 2005. 

Crosta, X., Debret, M., Denis, D., Courty, M. A., and Ther, O.: Holocene long- and short-term climate changes off Adélie Land, East Antarctica, Geochem., Geophys., Geosyst., 8, Q11009,, 2007. 

Crosta, X., Shukla, S. K., Ther, O., Ikehara, M., Yamane, M., and Yokoyama, Y.: Last Abundant Appearance Datum of Hemidiscus Karstenii Driven by Climate Change, Mar. Micropaleontol, 157, 10186,, 2020. 

Debret, M., Sebag, D., Crosta, X., Massei, N., Petit, J. R., Chapron, E., and Bout-Roumazeilles, V.: Evidence from wavelet analysis for a mid-Holocene transition in global climate forcing, Quaternary Sci. Rev., 28, 2675–2688,, 2009. 

Denton, G. H., Anderson, R. F., Toggweiler, J. R., Edwards, R. L., Schaefer, J. M., and Putnam, A. E.: The Last Glacial Termination, Science, 328, 1652 LP-1656,, 2010. 

Dezileau, L., Bareille, G., and Reyss, J. L.: The 231Pa/230Th ratio as a proxy for past changes in opal fluxes in the Indian sector of the Southern Ocean, Mar. Chem., 81, 105–117,, 2003. 

Diekmann, B.: Sedimentary Patterns in the Late Quaternary Southern Ocean, Deep Sea Res Part II, 54, 2350–2366,, 2007. 

Divine, D. V., Koç, N., Isaksson, E., Nielsen, S., Crosta, X., and Godtliebsen, F.: Holocene Antarctic climate variability from ice and marine sediment cores: Insights on ocean-atmosphere interaction, Quaternary Sci. Rev., 29, 303–312,, 2010. 

Ellison, C. R. W., Chapman, M. R., and Hall, I. R.: Surface and deep ocean interactions during the cold climate event 8200 years ago, Science, 312, 1929–1932,, 2006. 

EPICA Community Members: One-to-one coupling of glacial climate variability in Greenland and Antarctica, Nature, 444, 195–198,, 2006. 

Esper, O. and Gersonde, R.: Quaternary Surface Water Temperature Estimations: New Diatom Transfer Functions for the Southern Ocean, Palaeogeogr. Palaeocl., 414, 1–19,, 2014. 

Esper, O., Gersonde, R., and Kadagies, N.: Diatom Distribution in Southeastern Pacific Surface Sediments and Their Relationship to Modern Environmental Variables, Palaeogeogr. Palaeocl., 287, 1–27,, 2010. 

Etourneau, J., Collins, L. G., Willmott, V., Kim, J.-H., Barbara, L., Leventer, A., Schouten, S., Sinninghe Damsté, J. S., Bianchini, A., Klein, V., Crosta, X., and Massé, G.: Holocene climate variations in the western Antarctic Peninsula: evidence for sea ice extent predominantly controlled by changes in insolation and ENSO variability, Clim. Past, 9, 1431–1446,, 2013. 

Etourneau, J., Sgubin, G., Crosta, X., Swingedouw, D., Willmott, V., Barbara, L., Houssais, M-N, Schouten, S., Damsté, J. S. S., Goosse, H., Escutia, C., Crespin, J., Massé, G., and Kim, J.-H.: Ocean Temperature Impact on Ice Shelf Extent in the Eastern Antarctic Peninsula, Nat. Comm., 10, 1431–1446,, 2019. 

Fischer, H., Fundel, F., Ruth, U., Twarloh, B., Wegner, A., Udisti, R., Becagli, S., Castellano, E., Morganti, A., Severi, M., Wolff, E., Littot, G., Röthlisberger, R., Mulvaney, R., Hutterli, M. A., Kaufmann, P., Federer, U., Lambert, F., Bigler, M., Hansson, M., Jonsell, U., de Angelis, M., Boutron, C., Siggaard-Andersen, M.-L., Steffensen, J. P., Barbante, C., Gaspari, V., Gabrielli, P., and Wagenbach,, D.: Reconstruction of millennial changes in dust emission, transport and regional sea ice coverage using the deep EPICA ice cores from the Atlantic and Indian Ocean sector of Antarctica, Earth Planet. Sc. Lett., 260, 340–354,, 2007. 

Fletcher, M. S. and Moreno, P. I.: Zonally symmetric changes in the strength and position of the Southern Westerlies drove atmospheric CO2 variations over the past 14 k.y, Geology, 39, 419–422,, 2011. 

Gersonde, R. and Zielinski, U.: The reconstruction of Late Quaternary Antarctic sea-ice distribution – the use of diatoms as a proxy for sea ice, Palaeogeogr. Palaeocl., 162, 263–286, 2000. 

Graham, R. M., De Boer, A. M., Heywood, K. J., Chapman, M. R., and Stevens, D. P.: Southern Ocean fronts: Controlled by wind or topography? J. Geophys. Res.-Oceans, 117, C08018,, 2012. 

Guiot, J. and de Vernal, A.: Is spatial autocorrelation introducing biases in the apparent accuracy of paleoclimatic reconstructions?, Quaternary Sci. Rev., 30, 1965–1972,, 2011. 

Guiot, J., de Beaulieu, J. L., Cheddadi, R., David, F., Ponel, P., and Reille, M.: The climate in Western Europe during the last Glacial/Interglacial cycle derived from pollen and insect remains, Palaeogeogr. Palaeocl., 103, 73–93,, 1993. 

Haddam, N. A., Siani, G., Michel, E., Kaiser, J., Lamy, F., Duchamp-Alphonse, S., Hefter, J., Braconnot, P., Dewilde, F., Isgüder, G., Tisnerat-Laborde, N., Thil, F., Durand, N., and Kissel, C.: Changes in latitudinal sea surface temperature gradients along the Southern Chilean margin since the last glacial, Quaternary Sci. Rev., 194, 62–76,, 2018. 

Hall, A. and Visbeck, M.: Synchronous Variability in the Southern Hemisphere Atmosphere, Sea Ice, and Ocean Resulting from the Annular Mode, J. Climate, 15, 3043–3057,<3043:SVITSH>2.0.CO;2, 2002. 

Hall, B. L., Henderson, G. M., Baroni, C., and Kellogg, T. B.: Constant Holocene Southern-Ocean 14C Reservoir Ages and Ice-Shelf Flow Rates, Earth Planet. Sci. Lett., 296, 115–123,, 2010. 

Hammer, Ø., Harper, D. A. T., and Ryan, P. D.: PAST: paleontological statistics software package for education and data analysis, Palaeontologia Electronica, 4, 4, available at: (last access: 6 August 2020), 2001. 

Hughen, K. A., Overpeck, J. T., Peterson, L. C., and Trumbore, S.: Rapid climate changes in the tropical Atlantic region during the last deglaciation, Nature, 380, 51–54,, 1996. 

Iizuka, Y., Hondoh, T., and Fujii, Y.: Antarctic sea ice extent during the Holocene reconstructed from inland ice core evidence, J. Geophys. Res.-Atmos., 113, 1–8,, 2008. 

Jouzel, J., Masson-Delmotte, V., Cattani, O., Dreyfus, G., Falourd, S., Hoffmann, G., Minster, B., Nouet, J., Barnola, J. M., Chappellaz, J., Fischer, H., Gallet, J. C., Johnsen, S., Leuenberger, M., Loulergue, L., Luethi, D., Oerter, H., Parrenin, F., Raisbeck, G., Raynaud, D., Schilt, A., Schwander, J., Selmo, E., Souchez, R., Spahni, R., Stauffer, B., Steffensen, J. P., Stenni, B., Stocker, T. F., Tison, J. L., Werner, M., and Wolff, E. W.: Orbital and millennial Antarctic climate variability over the past 800,000 years, Science, 317, 793–796,, 2007. 

Katsuki, K., Ikehara, M., Yokoyama, Y., Yamane, M., and Khim, B. K.: Holocene migration of oceanic front systems over the Conrad Rise in the Indian Sector of the Southern Ocean, J. Quaternary Sci., 27, 203–210,, 2012. 

Kuroda, Y.: On the Origin of the Solar Cycle Modulation of the Southern Annular Mode, J. Geophys. Res.-Atmos., 123, 1959–1969,, 2018. 

Kuroda, Y. and Kodera, K.: Solar Cycle Modulation of the Southern Annular Mode, Geophys. Res. Lett., 32, L13802,, 2005. 

Lamy, F., Hebbeln, D., Röhl, U., and Wefer, G.: Holocene rainfall variability in southern Chile: a marine record of latitudinal shifts of the Southern Westerlies, Earth Planet. Sci. Lett., 185, 369–382,, 2001. 

Lamy, F., Kilian, R., Arz, H. W., Francois, J.-P., Kaiser, J., Prange, M., and Steinke, T.: Holocene changes in the position and intensity of the southern westerly wind belt, Nat. Geosci., 3, 695–699, 2010. 

Laskar, J., Robutel, P., Joutel, F., Gastineau, M., Correia, A., C., M., and Levrard, B.: A Long-Term Numerical Solution for the Insolation Quantities of the Earth, Astron. Astrophys., 428, 261–85,, 2004. 

Latif, M., Martin, T., and Park, W.: Southern Ocean sector centennial climate variability and recent decadal trends, J. Climate, 26, 7767–7782,, 2013. 

Leventer, A., Domack, E. W., Ishman, S. E., Brachfeld, S., McClennen, C. E., and Manley, P.: Productivity cycles of 200–300 years in the Antarctic Peninsula region: understanding linkages among the sun, atmosphere, oceans, sea ice, and biota, Geol. Soc. Am. Bulletin, 108, 1626–1644,<1626:PCOYIT>2.3.CO;2, 1996. 

Levermann, A., Schewe, J., and Montoya, M.: Lack of bipolar see-saw in response to Southern Ocean wind reduction, Geophys. Res. Lett., 34, 1–5,, 2007. 

Locarnini, R. A., Mishonov, A. V, Antonov, J. I., Boyer, T. P., Garcia, H. E., Baranova, O. K., Zweng, C. R., Paver, J. R., Reagan, D. R., Johnson, D. R., Hamilton, M., and Seidov, D.: World ocean atlas 2013. Volume 1, Temperature, edited by: Levitus, S. and Mishonov, A., NOAA Atlas NESDIS 73, Washington, DC, USA, 40 pp., 2013. 

Lovenduski, N. S. and Gruber, N.: Impact of the Southern Annular Mode on Southern Ocean circulation and biology, Geophys. Res. Lett., 32, 1–4,, 2005. 

Mayewski, P.A., Maasch, K.A., Dixon, D., Sneed, S.B., Oglesby, R., Korotkikh, E., Potocki, M., Grigholm, B., Kreutz, K., Kurbatov, A.V., and Spaulding, N.: West Antarctica's sensitivity to natural and human-forced climate change over the Holocene, J. Quaternary Sci., 28, 40–48,, 2013. 

McManus, J. F., Francois, R., Gherardi, J.-M., Keigwin, L. D., and Brown-Leger, S.: Collapse and rapid resumption of Atlantic meridional circulation linked to deglacial climate changes, Nature, 428, 834–837,, 2004. 

Meredith, M. P. and Hogg, A. M.: Circumpolar response of Southern Ocean eddy activity to a change in the Southern Annular Mode, Geophys. Res. Lett., 33, 2–5,, 2006. 

Monnin, E., Indermühle, A., Dällenbach, A., Flückiger, J., Stauffer, B., Stocker, T. F., Raynaud, D., and Barnola, J.-M.: Atmospheric CO2 Concentrations over the Last Glacial Termination, Science, 291, 112 LP-114,, 2001. 

Moreno, P. I., Vilanova, I., Villa-Martínez, R., Garreaud, R. D., Rojas, M., and De Pol-Holz, R.: Southern annular mode-like changes in southwestern Patagonia at centennial timescales over the last three millennia, Nat. Commun., 5, 4375,, 2014. 

Nielsen, S. H. H., Koç, N., and Crosta, X.: Holocene climate in the Atlantic sector of the Southern Ocean: Controlled by insolation or oceanic circulation?, Geology, 32, 317–320,, 2004. 

Oiwane, H., Ikehara, M., Suganuma, Y., Miura, H., Nakamura, Y., Sato, T., Nogi, Y., Yamane, M., and Yokoyama, Y.: Sediment Waves on the Conrad Rise , Southern Indian Ocean?: Implications for the Migration History of the Antarctic Circumpolar Current, Mar. Geol., 348, 27–36,, 2014. 

Orme, L. C., Crosta, X., Miettinen, A., Divine, D., Husum, K., Isaksson, E., Wacker, L., Mohan, R., Ther, O., and Ikehara, M.: Diatom percentage abundances for sediment core KH-10-7 COR1GC from the Conrad Rise, Southern Ocean, PANGAEA,, 2020a. 

Orme, L. C., Crosta, X., Miettinen, A., Divine, D., Husum, K., Isaksson, E., Wacker, L., Mohan, R., Ther, O., and Ikehara, M.: Sea surface temperature reconstruction for sediment core KH-10-7 COR1GC from the Conrad Rise, Southern Ocean, PANGAEA,, 2020b. 

Orsi, A. H., Whitworth III, T., and Nowlin Jr, W. D.: On the meridional extent and fronts of the Antarctic Circumpolar Current, Deep Sea Res. Part I, 42, 641–673,, 1995. 

Park, Y.-H., Charriaud, E., and Fieux, M.: Thermohaline structure of the Antarctic Surface Water/Winter Water in the Indian sector of the Southern Ocean, J. Mar. Syst., 17, 5–23,, 1998. 

Parkinson, C. L.: A 40-y record reveals gradual Antarctic sea ice increases followed by decreases at rates far exceeding the rates seen in the Arctic, Proc. Nat. Am. Soc., 116, 14414–14423, 2019. 

Parkinson, C. L. and Cavalieri, D. J.: Antarctic sea ice variability and trends, 1979–2010, The Cryosphere, 6, 871–880,, 2012. 

Parrenin, F., Masson-Delmotte, V., Köhler, P., Raynaud, D., Paillard, D., Schwander, J., Barbante, C., Landais, A., Wegner, A., and Jouzel, J.: Synchronous change of atmospheric CO2 and Antarctic temperature during the last deglacial warming, Science, 339, 1060–1063,, 2013. 

Past Interglacials Working Group: Interglacials of the Last 800,000 Years, Rev. Geophys., 54, 162–219,, 2016. 

Pauthenet, E., Roquet, F., Madec, G., Guinet, C., Hindell, M., McMahon, C. R., Harcourt, R., and Nerini, D.: Seasonal meandering of the Polar Front upstream of the Kerguelen Plateau, Geophys. Res. Lett., 45, 9774–9781,, 2018. 

Pedro, J. B., Bostock, H., C., Bitz, C., M., He, F., Vandergoes, M., J., Steig, E., J., Chase, B., M., Krause, C., E., Rasmussen, S., O., Markle, B., R., and Cortese, G.: The Spatial Extent and Dynamics of the Antarctic Cold Reversal, Nat. Geosci., 9, 51–55,, 2015. 

Pedro, J. B., Jochum, M., Buizert, C., He, F., Barker, S., and Rasmussen, S. O.: Beyond the bipolar seesaw?: Toward a process understanding of interhemispheric coupling, Quaternary Sci. Rev., 192, 27–46,, 2018. 

Pichon, J.-J., Bareille, G., Labracherie, M., Labeyrie, L., D., Baudrimont, A., and Turon, J.-L.: Quantification of the Biogenic Silica Dissolution in Southern Ocean Sediments, Quaternary Res., 37, 361–378,, 1992. 

Pollard, R. T. and Read, J. F.: Circulation pathways and transports of the Southern Ocean in the vicinity of the Southwest Indian Ridge, J. Geophys. Res.-Oceans, 106, 2881–2898,, 2001. 

Rasmussen, S., O., Bigler, M., Blockley, S., P., Blunier, T., Buchardt, S., L., Clausen, H., B., Cvijanovic, I., Dahl-Jensen, D., Johnsen, S., J., and Fischer, H.: A Stratigraphic Framework for Abrupt Climatic Changes during the Last Glacial Period Based on Three Synchronized Greenland Ice-Core Records: Refining and Extending the INTIMATE Event Stratigraphy, Quaternary Sci. Rev., 106, 14–28,, 2014. 

Reimer, P. J., Bard, E., Bayliss, A., Beck, J. W., Blackwell, P. G., Ramsey, C. B., Buck, C. E., Cheng, H., Edwards, R. L, Friedrich, M., Grootes, P. M., Guilderson, T. P., Haflidason, H., Hajdas, I., Hatté, C., Heaton, T. J., Hoffmann, D. L., Hogg, A. G., Hughen, K. A., Kaiser, K. F., Kromer, B., Manning, S. W., Niu, M., Reimer, R. W., Richards, D. A., Scott, E. M., Southon, J. R.., Staff, R. A., Turney, C. S. M., and van der Plicht, J.: IntCal13 and Marine13 radiocarbon age calibration curves 0–50,000 years cal BP, Radiocarbon, 55, 1869–1887,, 2013. 

Renssen, H., Goosse, H., Crosta, X., and Roche, D. M.: Early Holocene aurentide ice sheet deglaciation causes cooling in the high-latitude southern hemisphere through oceanic teleconnection, Paleoceanography, 25, 1–15,, 2010. 

Romero, O., E, Armand, L., K., Crosta, X., and Pichon, J-J.: The Biogeography of Major Diatom Taxa in Southern Ocean Surface Sediments: 3. Tropical/Subtropical Species, Palaeogeogr. Palaeocl., 223, 49–65,, 2005. 

Romero, O., E., Crosta, X., Kim, J.-H., Pichevin, L., and Crespin, J.: Rapid Longitudinal Migrations of the Filament Front off Namibia (SE Atlantic) during the Past 70 kyr, Global Planet. Change, 125, 1–12,, 2015. 

Saunders, K. M., Roberts, S. J., Perren, B., Butz, C., Sime, L., Davies, S., Van Nieuwenhuyze, W., Grosjean, M., and Hodgson, D. A.: Holocene dynamics of the Southern Hemisphere westerly winds and possible links to CO2 outgassing, Nat. Geosci., 11, 650–655,, 2018. 

Schlitzer, R.: Ocean Data View, available at: (last access: January 2019), 2014. 

Schulz, M. and Mudelsee, M.: REDFIT?: estimating red-noise spectra directly from unevenly spaced paleoclimatic time series, Comput. Geosci., 28, 421–426,, 2002. 

Schweitzer, P. N.: Monthly averaged polar sea-ice concentration. U.S. Geological Survey Digital Data Series, Virginia,, 1995. 

Screen, J. A., Gillett, N. P., Stevens, D. P., Marshall, G. J., and Roscoe, H. K.: The role of eddies in the Southern Ocean temperature response to the Southern Annular Mode, J. Climate, 22, 806–818,, 2009. 

Shevenell, A. E., Ingalls, A. E., Domack, E. W., and Kelly, C.: Holocene Southern Ocean surface temperature variability west of the Antarctic Peninsula, Nature, 470, 250–254,, 2011. 

Siani, G., Michel, E., Pol-holz, R. De, Devries, T., Lamy, F., Isguder, G., Dewilde, F., and Lourantou, A.: Carbon isotope records reveal precise timing of enhanced Southern Ocean upwelling during the last deglaciation, Nat. Commun., 4, 2758,, 2013. 

Sokolov, S. and Rintoul, S. R.: Circumpolar structure and distribution of the Antarctic Circumpolar Current fronts: 1. Mean circumpolar paths, J. Geophys. Res.-Oceans, 114, C11018,, 2009. 

Stenni, B., Buiron, D., Frezzotti, M., Albani, S., Barbante, C., Bard, E., Barnola, J. M., Baroni, M., Baumgartner, M., Bonazza, M., Capron, E., Castellano, E., Chappellaz, J., Delmonte, B., Falourd, S., Genoni, L., Iacumin, P., Jouzel, J., Kipfstuhl, S., Landais, A., Lemieux-Dudon, B., Maggi, V., Masson-Delmotte, V., Mazzola, C., Minster, B., Montagnat, M., Mulvaney, R., Narcisi, B., Oerter, H., Parrenin, F., Petit, J. R., Ritz, C., Scarchilli, C., Schilt, A., Schüpbach, S., Schwander, J., Selmo, E., Severi, M., Stocker T. F., and Udisti, R.: Expression of the bipolar see-saw in Antarctic climate records during the last deglaciation, Nat. Geosci., 4, 46–49,, 2011. 

Stocker, T. F. and Johnsen, S. J.: A minimum thermodynamic model for the bipolar seesaw, Paleoceanography, 18, 1087,, 2003.  

Swart, N. C. and Fyfe, J. C.: Observed and simulated changes in the Southern Hemisphere surface westerly wind-stress, Geophys. Res. Lett., 39, 6–11,, 2012. 

Turney, C. S. M., Jones, R. T., Fogwill, C., Hatton, J., Williams, A. N., Hogg, A., Thomas, Z. A., Palmer, J., Mooney, S., and Reimer, R. W.: A 250-year periodicity in Southern Hemisphere westerly winds over the last 2600 years, Clim. Past, 12, 189–200,, 2016. 

van Beek, P., Reyss, J.-L., Paterne, M., Gersonde, R., van der Loeff, M. R., and Kuhn, G.: 226Ra in barite: Absolute dating of Holocene Southern Ocean sediments and reconstruction of sea-surface reservoir ages, Geology, 30, 731–734,<0731:RIBADO>2.0.CO;2, 2002. 

Verleyen, E., Hodgson, D. A., Sabbe, K., Cremer, H., Emslie, S. D., Gibson, J., Hall, B., Imura, S., Kudoh, S., Marshall, G. J., McMinn, A., Melles, M., Newman, L., Roberts, D., Roberts, S. J., Singh, S. M., Sterken, M., Tavernier, I., Verkulich, S., Van de Vyver, E., Van Nieuwenhuyze, W., Wagner, B., and Vyverman, W.: Post-glacial regional climate variability along the East Antarctic coastal margin—evidence from shallow marine and coastal terrestrial records, Earth-Sci. Rev., 104, 199–212,, 2011. 

Walker, M., Head, M., J., Berkelhammer, M., Björck, S., Cheng, H., Cwynar, L., Fisher, D., Gkinis, V., Long, A., and Lowe, J.: Formal Ratification of the Subdivision of the Holocene Series/Epoch (Quaternary System/Period): Two New Global Boundary Stratotype Sections and Points (GSSPs) and Three New Stages/Subseries, Episodes, 41, 213–223,, 2018. 

Xiao, W., Esper, O., and Gersonde, R.: Last Glacial - Holocene climate variability in the Atlantic sector of the Southern Ocean, Quaternary Sci. Rev., 135, 115–137,, 2016. 

Zielinski, U. and Gersonde. R.: Diatom Distribution in Southern Ocean Surface Sediments (Atlantic Sector): Implications for Paleoenvironmental Reconstructions, Palaeogeogr. Palaeocl., 129, 213–250,, 1997. 

Publications Copernicus
Short summary
A record of past sea temperature in the Indian sector of the Southern Ocean, spanning the last 14 200 years, has been developed by analysis of fossil diatoms in marine sediment. During the late deglaciation the reconstructed temperature changes were highly similar to those over Antarctica, most likely due to a reorganisation of global ocean and atmospheric circulation. During the last 11 600 years temperatures gradually cooled and became increasingly variable.
A record of past sea temperature in the Indian sector of the Southern Ocean, spanning the last...