Articles | Volume 18, issue 3
Clim. Past, 18, 465–483, 2022

Special issue: Reconstructing Southern Ocean sea-ice dynamics on glacial-to-historical...

Clim. Past, 18, 465–483, 2022

Research article 14 Mar 2022

Research article | 14 Mar 2022

Sea ice changes in the southwest Pacific sector of the Southern Ocean during the last 140 000 years

Sea ice changes in the southwest Pacific sector of the Southern Ocean during the last 140 000 years
Jacob Jones1, Karen E. Kohfeld1,2, Helen Bostock3,4, Xavier Crosta5, Melanie Liston6, Gavin Dunbar6, Zanna Chase7, Amy Leventer8, Harris Anderson7, and Geraldine Jacobsen9 Jacob Jones et al.
  • 1School of Resource and Environmental Management, Simon Fraser University, Burnaby, Canada
  • 2School of Environmental Science, Simon Fraser University, Burnaby, Canada
  • 3School of Earth and Environmental Sciences, The University of Queensland, Brisbane, Australia
  • 4National Institute of Water and Atmospheric Research (NIWA), Wellington, New Zealand
  • 5Université de Bordeaux, CNRS, EPHE, UMR 5805 EPOC, Pessac, France
  • 6Antarctic Research Centre, Victoria University of Wellington, Wellington, New Zealand
  • 7Institute of Marine and Antarctic Studies, University of Tasmania, Hobart, Australia
  • 8Geology Department, Colgate University, Hamilton, NY, USA
  • 9Centre for Accelerator Science, Australian Nuclear Science and Technology Organisation, Lucas Heights, NSW, Australia

Correspondence: Jacob Jones (


Sea ice expansion in the Southern Ocean is believed to have contributed to glacial–interglacial atmospheric CO2 variability by inhibiting air–sea gas exchange and influencing the ocean's meridional overturning circulation. However, limited data on past sea ice coverage over the last 140 ka (a complete glacial cycle) have hindered our ability to link sea ice expansion to oceanic processes that affect atmospheric CO2 concentration. Assessments of past sea ice coverage using diatom assemblages have primarily focused on the Last Glacial Maximum (∼21 ka) to Holocene, with few quantitative reconstructions extending to the onset of glacial Termination II (∼135 ka). Here we provide new estimates of winter sea ice concentrations (WSIC) and summer sea surface temperatures (SSST) for a full glacial–interglacial cycle from the southwestern Pacific sector of the Southern Ocean using the modern analog technique (MAT) on fossil diatom assemblages from deep-sea core TAN1302-96. We examine how the timing of changes in sea ice coverage relates to ocean circulation changes and previously proposed mechanisms of early glacial CO2 drawdown. We then place SSST estimates within the context of regional SSST records to better understand how these surface temperature changes may be influencing oceanic CO2 uptake. We find that winter sea ice was absent over the core site during the early glacial period until MIS 4 (∼65 ka), suggesting that sea ice may not have been a major contributor to early glacial CO2 drawdown. Sea ice expansion throughout the glacial–interglacial cycle, however, appears to coincide with observed regional reductions in Antarctic Intermediate Water production and subduction, suggesting that sea ice may have influenced intermediate ocean circulation changes. We observe an early glacial (MIS 5d) weakening of meridional SST gradients between 42 and 59 S throughout the region, which may have contributed to early reductions in atmospheric CO2 concentrations through its impact on air–sea gas exchange.

1 Introduction

Antarctic sea ice has been suggested to have played a key role in glacial–interglacial atmospheric CO2 variability (e.g., Stephens and Keeling, 2000; Ferrari et al., 2014; Kohfeld and Chase, 2017; Stein et al., 2020). Sea ice has been dynamically linked to several processes that promote deep ocean carbon sequestration, namely by (1) reducing deep ocean outgassing by ice-induced “capping” and surface water stratification (Stephens and Keeling, 2000; Rutgers van der Loeff et al., 2014) and (2) influencing ocean circulation through water mass formation and deep-sea stratification, leading to reduced diapycnal mixing and reduced CO2 exchange between the surface and deep ocean (Toggweiler, 1999; Bouttes et al., 2010; Ferrari et al., 2014). Numerical modelling studies have shown that sea-ice-induced capping, stratification, and reduced vertical mixing may be able to account for a significant portion of the total CO2 variability on glacial–interglacial timescales (between 40–80 ppm) (Stephens and Keeling, 2000; Galbraith and de Lavergne, 2018; Marzocchi and Jansen, 2019; Stein et al., 2020). However, debate continues surrounding the timing and magnitude of sea ice impacts on glacial-scale carbon sequestration (e.g., Morales Maquede and Rahmstorf, 2002; Archer et al., 2003; Sun and Matsumoto, 2010; Kohfeld and Chase, 2017).

Past Antarctic sea ice coverage has been estimated primarily through diatom-based reconstructions, with most work focusing on the Last Glacial Maximum (LGM), specifically the EPILOG time slice as outlined in Mix et al. (2001), corresponding to 23 to 19 thousand years before present (ka, calibrated backwards from 1950). During the LGM, these reconstructions suggest that winter sea ice expanded by 7–10 latitude (depending on the sector of the Southern Ocean), which corresponds to substantial expansion of total winter sea ice coverage compared to modern observations (Gersonde et al., 2005; Benz et al., 2016; Lhardy et al., 2021). Currently, only a handful of studies provide quantitative sea ice coverage estimates back to the penultimate glaciation, Marine Isotope Stage (MIS) 6 (∼194 to 135 ka) (Gersonde and Zielinksi, 2000; Crosta et al., 2004; Schneider-Mor et al., 2012; Esper and Gersonde 2014a; Ghadi et al., 2020). These studies primarily cover the Atlantic sector, with only one published sea ice record from each of the Indian (SK200-33 from Ghadi et al., 2020), eastern Pacific (PS58/271-1 from Esper and Gersonde, 2014a), and southwestern Pacific sectors (SO136-111 from Crosta et al., 2004). These glacial–interglacial sea ice records show heterogeneity between sectors in both timing and coverage. While the Antarctic Zone (AZ) in the Atlantic sector experienced early sea ice advance corresponding to MIS 5d cooling (i.e., 115 to 105 ka) (Gersonde and Zielinksi, 2000; Bianchi and Gersonde, 2002; Esper and Gersonde, 2014a), the Indian and Pacific sector cores in the AZ show only minor sea ice advances during this time (Crosta et al., 2004; Ghadi et al., 2020). The lack of spatial and temporal resolution has resulted in significant uncertainty in our ability to evaluate the timing and magnitude of sea ice change during a full glacial cycle across the Southern Ocean, and to link sea ice to glacial–interglacial CO2 variability.

This paper provides new winter sea ice concentration (WSIC) and summer sea surface temperature (SSST) estimates for the southwestern Pacific sector of the Southern Ocean over the last 140 ka. WSIC, which is a grid-scale observation of the mean state fraction of ocean area that is covered by sea ice over the sample period, and SSST estimates are produced by applying the modern analog technique (MAT) to fossil diatom assemblages from sediment core TAN1302-96 (59.09 S, 157.05 E, water depth 3099 m). We place this record within the context of sea ice and SSST changes from the region using previously published records from SO136-111 (56.66 S, 160.23 E, water depth 3912 m), which has recalculated WSIC and SSST estimates presented in this study, and nearby marine core E27-23 (59.61 S, 155.23 E, water depth 3182 m) (Ferry et al., 2015). Using these records, we compare the timing of sea ice expansion to early glacial–interglacial CO2 variability to test the hypothesis that the initial CO2 drawdown (∼115 to 100 ka) resulted from reduced air–sea gas exchange in response to sea ice capping and surface water stratification (Kohfeld and Chase, 2017). We then consider alternative oceanic drivers of early atmospheric CO2 variability and place our SSST estimates within the context of other studies to examine how regional cooling and a weakening in meridional SST gradients might affect air–sea disequilibrium and early CO2 drawdown (Khatiwala et al., 2019). Finally, we compare our WSIC estimates with regional reconstructions of Antarctic Intermediate Water (AAIW) production and subduction variability using previously published carbon isotope analyses on benthic foraminifera from intermediate to deep-water depths in the southwest Pacific sector of the Southern Ocean, to test the hypothesis that sea ice expansion is dynamically linked to AAIW production and variability (Ronge et al., 2015).

2 Methods

2.1 Study site and age determination

We reconstruct diatom-based WSIC and SSST using marine sediment core TAN1302-96 (59.09 S, 157.05 E, water depth 3099 m) (Fig. 1). The 364 cm core was collected in March 2013 using a gravity corer during the return of the RV Tangaroa from the Mertz Polynya in eastern Antarctica (Williams, 2013). The core is situated in the western Pacific sector of the Southern Ocean, on the southwestern side of the Macquarie Ridge, approximately 3–4 south of the average position of the Polar Front (PF) at 157 E (Sokolov and Rintoul, 2009).

Figure 1Map of the southwestern Pacific sector of the Southern Ocean including the study site, TAN1302-96 (blue circle), and additional published cores providing sea ice extent data, SO136-111 and E27-23 (green circles); SST reconstructions (red circles); and δ13C of benthic foraminifera (yellow circles). Note that some cores may not appear present in the figure because of their proximity to other cores. Data for all cores are provided in Table 2. Dashed lines show the average location of the subtropical and polar fronts (Smith et al., 2013; Bostock et al., 2015), and red and blue lines show mean positions of modern summer sea ice (SSI) and winter sea ice (WSI) extents, respectively (Reynolds et al., 2002, 2007).

The age model for TAN1302-96 (Figs. 2 and 3) was based on a combination of radiocarbon dating of mixed foraminiferal assemblages and stable oxygen isotope stratigraphy on Neogloboquadrina pachyderma (180–250 µm). Seven accelerator mass spectrometry (AMS) 14C samples were collected (Table A1 in Appendix A) and consisted of mixed assemblages of planktonic foraminifera (N. pachyderma and Globigerina bulloides, >250µm). Three of the seven radiocarbon samples (NZA 57105, 57109, and 61429) were previously published in Prebble et al. (2017), and four additional samples (OZX 517-520) were added to improve the dating reliability (Table A1 in Appendix A). OZX 519 and OZX 520 produced dates that were not distinguishable from the background (>57.5 ka) and were subsequently excluded from the age model. The TAN1302-96 oxygen isotopes were run at the National Institute of Water and Atmospheric Research (NIWA) using the Kiel IV individual acid-on-sample device and analysed using Finnigan MAT 252 mass spectrometer. The precision is ±0.07 % for δ18O and ±0.05 % for δ13C.

Figure 2Age model of TAN1302-96. Red circles indicate the depth of AMS 14C samples, and yellow circles indicate tie points between the TAN1302-96 oxygen isotope stratigraphy and the LR04 benthic stack (Lisiecki and Raymo, 2005). Two radiocarbon dates, OZX 519 and 520 (at 130 and 170 cm, respectively), were not included in the age model as they produced dates that were NDFB (not distinguishable from the background).


Figure 3Age model of TAN1302-96. Tie points are depicted as yellow dots, and grey shading represents associated uncertainty between tie points. The age model used a marine reservoir calibration of 1000±100 years.


The age model was constructed using the “Undatable” MATLAB software (Lougheed and Obrochta, 2018) by bootstrapping at 10 % and using an xfactor of 0.1 (Lougheed and Obrochta, 2019), which scales Gaussian distributions of sediment accumulation uncertainty (Table A2 in Appendix A). Below 100 cm, nine tie points were selected at positions of maximum change in δ18O and were correlated to the LR04 benthic stack (Lisiecki and Raymo, 2005) (Fig. 2; Table A2 in Appendix A). Uncertainty associated with stratigraphic correlation to the LR04 stack has been estimated to be ±4 ka (Lisiecki and Raymo, 2005). We used a conservative marine reservoir age (MRA) for radiocarbon calibration of 1000±100 years, in line with regional estimates in Paterne et al. (2019) and modelled estimates by Butzin et al. (2017, 2020). The age model shows that TAN1302-96 extends to at least 140 ka, capturing a full glacial–interglacial cycle. Linear sedimentation rates in TAN1302-96 were observed to be higher during interglacial periods, averaging ∼3.5 cm ka−1, compared to glacial periods, averaging ∼2.5 cm ka−1. It is worth noting that there can be significant MRA variability over time due to changes in ocean ventilation, sea ice coverage, and wind strength, specifically in the polar high latitudes (Heaton et al., 2020), and as a result, caution should be taken when interpreting the precision of radiocarbon dates. For more information on age model construction and selection, refer to the Supplement.

2.2 Diatom analysis

TAN1302-96 was sampled every 3–4 cm throughout the core except between 130–180 cm, where samples were collected every 10 cm due to limited availability of sample materials (Table A3 in Appendix A). Diatom slide preparation followed two procedures. The first approach approximated the methods outlined in Renberg (1990), while the second followed the protocol outlined in Warnock and Scherer (2015). To ensure there were no biases between preparation techniques, results from each technique were first visually compared followed by a comparison of sample means (see Fig. B1 in Appendix B). No biases in the data were observed between methods.

The first procedure was conducted at Victoria University of Wellington and Simon Fraser University on samples every 10 cm throughout the core. Sediment samples contained high concentrations of diatoms with little carbonaceous or terrigenous materials, so no dissolving aids were used. Instead, approximately 50 mg of sediment was weighed, placed into a 50 mL centrifuge tube, and topped up with 40 mL of deionized water. Samples were then manually shaken to disaggregate sediment, followed by a 10 s mechanical stir using a vortex machine. Samples were then left to settle for 25 s. A total of 0.25 mL of the solution was then pipetted onto a microscope slide from a consistent depth, where it was left to dry overnight. Once the sample had dried, coverslips were permanently mounted to the slide using Permount, a high refractive index mountant. Slides were redone if they contained too many diatoms and identification was not possible, or if they contained too few diatoms (generally <40 specimens per transect). Sediment sample weight was adjusted to achieve the desired dilution.

The second procedure was conducted at Colgate University on samples every 3–4 cm throughout the core. Oven-dried samples were placed into a 20 mL vial with 1–2 mL of 10 % H2O2 and left to react for up to several days, followed by a brief (2–3 s) ultrasonic bath to disaggregate samples. The diatom solution was then added into a settling chamber, where microscope coverslips were placed on stages to collect settling diatoms. The chamber was gradually emptied through an attached spigot, and samples were evaporated overnight. Cover slips were permanently mounted onto the slides with Norland Optical Adhesive 61, a mounting medium with a high refractive index.

Diatom identification was conducted at Simon Fraser University using a Leica Leitz DMBRE light microscope using standard microscopy techniques. Following transverses, a minimum of 300 individual diatoms were identified at 1000× magnification from each sample throughout the core. Individuals were counted towards the total only if they represented at least one-half of the specimen so that fragmented diatoms were not counted twice. Identification was conducted to the highest taxonomic level possible, either to the species or species-group level. Taxonomic identification was conducted using numerous identification materials, including (but not limited to) Fenner et al. (1976), Fryxell and Hasle (1976, 1980), Johansen and Fryxell (1985), Hasle and Syversten (1997), Cefarelli et al. (2010), and Wilks and Armand (2017). The relative abundances were calculated by dividing the number of identified specimens of a particular species by the total number of identified diatoms from the sample. Based on previously established taxonomic groups (Crosta et al., 2004), diatoms were grouped into one of three categories based on temperature preference and sea ice tolerance. The following main taxonomic groups were used (Table 1):

  1. Sea ice group. This represents diatoms that thrive in or near the sea ice margin in SSTs generally ranging from −1 to 1 C.

  2. Permanently Open Ocean Zone (POOZ). This represents diatoms that thrive in open ocean conditions, with SSTs generally ranging from ∼2 to 10 C.

  3. Sub-Antarctic Zone (SAZ). This represents diatoms that thrive in warmer sub-Antarctic waters, with SSTs generally ranging from 11 to 14 C.

Table 1Species comprising each of the diatom taxonomic groups (updated from Crosta et al., 2004).

Download Print Version | Download XLSX

2.3 Modern analog technique

Past WSIC and SSST (January to March) were estimated for TAN1302-96 and recalculated for SO136-111 by applying the modern analog technique (MAT) to the fossil diatom assemblages, as outlined in Crosta et al. (1998, 2020). Summer (January to March) SST was estimated because it is considered to be a better explanatory variable than spring or annual SST (Esper et al., 2010; Esper and Gersonde, 2014b). The MAT reference database used for this analysis is comprised of 249 modern core top samples (analogs) located primarily in the Atlantic and Indian sectors from ∼40 S to the Antarctic coast. The age of the core tops included in the reference database have been assessed through radiocarbon and/or isotope stratigraphy when possible. Core tops were visually evaluated for selective diatom dissolution, so it is believed that sub-modern assemblages contain well-preserved and unbiased specimens. Modern SSST and WSIC were interpolated from the reference core locations using a 1×1 grid from the World Ocean Atlas (Locarnini et al., 2013) through the Ocean Data View (Schlitzer, 2005). The MAT was applied using the “bioindic” package (Guiot and de Vernal, 2011) through the R platform. Fossil diatom assemblages were compared to the modern analogs using 33 species or species groups to identify the five most similar modern analogs using both the LOG and CHORD distance. The dissimilarity threshold, above which the fossil assemblages are considered to be too dissimilar to the modern dataset, is fixed at the first quartile of random distances (Crosta et al., 2020). The reconstructed SSST and WSIC are the distance-weighted mean of the climate values associated with the selected modern analog (Guiot et al., 1993; Ghadi et al., 2020). Both MAT approaches produce an R2 value of 0.96 and a root mean square error of prediction (RMSEP) of ∼1C for SSST and an R2 of 0.93 and a RMSEP of 10 % for WSIC (Ghadi et al., 2020). As outlined in Ferry et al. (2015), we consider <15 % WSIC to represent an absence of winter sea ice, 15 %–40 % WSIC as present but unconsolidated, and >40 % to represent consolidated winter sea ice.

2.4 Additional core data

We use additional published marine cores from the southwestern Pacific throughout this analysis (Table 2), for WSIC comparisons (E27-23), %AAIW calculations (MD06-2990/SO136-003, MD06-2986, and MD97-2120), and regional SST gradient comparisons (SO136-003, FR1/94-GC3, ODP 1119-181, DSDP 594, and Q200).

Table 2Additional data on published marine cores used throughout this analysis.

Download Print Version | Download XLSX

3 Results

3.1 TAN1302-96 diatom assemblage results

In this core, 51 different species or species groups were identified, of which 33 were used in the transfer function. These 33 species represent >82 % of the total diatom assemblages (mean of 92 %). Permanently Open Ocean Zone (POOZ) diatoms made up the largest proportion of diatoms identified, representing between 72 %–91 % of the assemblage (Fig. 4), with higher values observed during warmer interstadial periods of MIS 1, 3, and 5. Sea ice diatoms made up the second most abundant group, representing between 0.5 %–7.5 % of the assemblage, with higher values observed during cooler stadial periods (MIS 2, 4, and 6). The Sub-Antarctic Zone group had relatively low abundances, with higher values occurring during warmer interstadial periods (MIS 5 and the Holocene) and briefly during MIS 4 at ∼65 ka.

Figure 4Diatom assemblages results from TAN1302-96 separated into percentage contribution from each taxonomic group (sea ice group, POOZ, and SAZ; see Table 1) over a full glacial–interglacial cycle. Using the modern analog technique (MAT), winter sea ice concentration (WSIC) and summer sea surface temperature (SSST) were estimated and compared against the δ18O signature of TAN1302-96.


3.2 TAN1302-96 SSST and WSIC estimates

There were no non-analog conditions observed in TAN1302-96 samples, and all estimates were calculated on five analogs. Estimates of SSST and WSIC from both LOG and CHORD MAT outputs produced similar results (Fig. 4). During Termination II, SSST began to rise from ∼1C at 140 ka (MIS 6) to ∼4.5C at 132 ka (MIS 5e–6 boundary). This warming corresponded with a decrease in WSIC from 48 % to approximately 0 % over the same time periods (Fig. 4). Reconstructed SSSTs were variable throughout MIS 5e, reaching a maximum value of ∼4.5C at 118 ka, after which they declined throughout MIS 5. During this period of SSST decline, winter sea ice was largely absent, punctuated by brief periods during which sea ice was present but unconsolidated (WSIC =15 % and 17 % at 105 and 85 ka, respectively). During MIS 4 (71 to 57 ka), SSST cooled to between roughly 1 and 3 C, and sea ice expanded to 36 %, such that it was present but unconsolidated for intervals of a few thousand years. SSST increased slightly from 1.5 C at 61 ka (during MIS 4) to ∼2.5C at 50 ka (during MIS 3), followed by a general cooling trend into MIS 2. Sea ice appears to have been largely absent during MIS 3 (57 to 29 ka), although sampling resolution is low, but increased rapidly to 48 % cover during MIS 2 where winter sea ice was consolidated over the core site. During MIS 2, SSST cooled to a minimum of <1C at 24.5 ka. After 18 ka, the site rapidly transitioned from cool, ice-covered conditions to warmer, ice-free winter conditions during the early deglaciation. This warming was interrupted by a brief cooling around 13.5 ka, following which SSSTs quickly reached their maximum values of ∼5C at 11.5 ka and remained relatively high throughout the rest of the Holocene. Winter sea ice was not present during the Holocene.

3.3 SO136-111 SSST and WSIC recalculation

In core SO136-111, the 33 species included in the transfer function represent values >79 % of the total diatom assemblages (mean of 91 %). There were no non-analog conditions observed in SO136-111 samples, and all estimates were calculated on five analogs. Recalculated estimates of SSST and WSIC from both LOG and CHORD MAT outputs produced similar results for SO136-111 (Fig. 5a, d). During Termination II, SSST rose from ∼2C at 137 ka (MIS 6) to a maximum value of 6 C at 125 ka (MIS5e), corresponding to a rapid decline in WSIC from 37 % to ∼0 % during the same period. SSST remained relatively high (between 4 and 5 C) from 125 ka until 115 ka where they declined to ∼2C. SSST remained variable from 110 ka until ∼40 ka, fluctuating between ∼2 and 4 C. Winter sea ice was largely absent during MIS 5, with a brief period where sea ice was present but unconsolidated (WSIC =17 % at 84 ka). Beginning at ∼76 ka, WSIC began to increase and continued throughout early MIS 4 to a maximum 36 % at 69 ka. WSIC remained present but unconsolidated throughout most of MIS 3 and 2 with brief periods of absence (WSIC ≤15 %) lasting a few thousand years. SSST and WSIC reached their coolest values and highest concentration at 24.5 ka before SSST increased to ∼5C and stabilized throughout the Holocene, while WSIC declined to virtually 0 % throughout the same period.

4 Discussion

4.1 Regional SSST and WSIC estimates

The new WSIC and SSST estimates from TAN1302-96 and recalculated WSIC estimates from SO136-111 show a coherent regional pattern (Fig. 5). TAN1302-96 shows slightly higher concentrations during MIS 2 (maximum WSIC =48 % at 24.5 ka) and 4 (maximum WSIC =37 % at 65 ka) compared with SO136-111 (maximum WSIC =35 % at 24.5 ka and 36 % at 68 ka, respectively), which can be explained by a more poleward position of TAN1302-96 relative to SO136-111. The estimates between cores differ during MIS 3, with seemingly lower WSIC in TAN1302-96 than in SO136-111, which might result from the low sampling resolution in TAN1302-96 during this period. Overall, these cores show a highly similar and coherent history of sea ice over the last 140 ka.

Figure 5(a) WSIC estimates using MAT from SO136-111 (recalculated in this study; see Appendix D); (b) WSIC estimates using GAM from E27-23 (Ferry et al., 2015); (c) WSIC estimates using MAT from TAN1302-96 (this study); (d) SSST estimates using MAT from TAN1302-96 (solid black line) and recalculated SSST for SO136-111 (black dotted line); (e) Antarctic atmospheric CO2 concentrations over 140 ka (Bereiter et al., 2015); (f) δ13C data from nearby cores MD06-2990/SO136-003, MD97-2120, and MD06-2986 (Ronge et al., 2015); (g) %Antarctic Intermediate Water (%AAIW) as calculated in Ronge et al. (2015), which tracks when core MD97-2120 was bathed primarily by AAIW (green) or Upper Circumpolar Deep Water (UCDW) (blue).


When compared with E27-23 (Fig. 5b), which is located only ∼120 km to the southwest of TAN1302-96 (Fig. 1), the TAN1302-96 core shows lower estimates of WSIC, especially during MIS 3. During early and mid-MIS 2, both cores show similar WSIC estimates, while later in MIS 2 (∼17 ka), E27-23 reports a maximum WSIC of 72 % compared to only 22 % at TAN1302-96. A discrepancy between estimates is also observed during the Holocene, with E27-23 reporting sea ice estimates of up to nearly 50 % during the mid-Holocene (∼6 ka), while TAN1302-96 experienced values well below the RMSEP of 10 %.

Possible explanations for the observed differences in WSIC estimates include (1) differences in statistical applications, (2) lateral sediment redistribution, (3) differences in laboratory protocols, (4) differences in diatom identification/counting methodology, and (5) selective diatom dissolution. Of these explanations, we believe that (1) and (2) are the most likely candidates and are discussed below (for further discussion on 3, 4, and 5, see Appendix C).

The first possible explanation is the use of different statistical applications. Ferry et al. (2015) used a generalized additive model (GAM) to estimate WSIC for both E27-23 and SO136-111, while we have used the MAT for TAN1302-96 and SO136-111. A simple comparison of WSIC estimates between the results in Ferry et al. (2015) and our recalculated WSIC estimates for SO136-111 can provide insights into the magnitude of estimation differences. Generally speaking, the GAM estimation produced higher WSIC estimates than the MAT (e.g., ∼50 % WSIC at 23 ka while the MAT produced ∼37 % for the same time period); however, we believe it is unlikely that statistical approaches alone could explain a larger difference (i.e., 50 %) between E27-23 and TAN1302-96.

The second possible explanation involves lateral sediment redistribution and focusing by the ACC. We estimated sediment focusing for E27-23 using 230Th data from Bradtmiller et al. (2009) together with dry bulk density estimated using calcium carbonate content (Froelich, 1991). Both sedimentation rates and focusing factors (FF) for the E27-23 are relatively high (maximum ≅35 cm ka−1 and 26, respectively) during the LGM and Holocene, which could influence the reliability of WSIC and SSST estimation (see Fig. C1 in Appendix C). Several peaks in focusing occurring around 16, 12, and 3 ka appear to closely correspond to periods of peak WSIC (∼67 %, ∼54 %, and ∼35 %, respectively), suggesting a possible link. Lateral redistribution could artificially increase or decrease relative abundances of some diatom groups, which could lead to over- or under-estimations of sea ice coverage. Thorium analysis for TAN1302-96 is beyond the scope of this study; however, future work could help address this uncertainty.

Although we are unable to identify the specific cause of the differences, we suggest considering the results from all cores when drawing conclusions of regional sea ice history.

4.2 The role of sea ice on early CO2 drawdown

Kohfeld and Chase (2017) hypothesized that the initial drawdown of atmospheric CO2 (∼35 ppm) during the glacial inception of MIS 5d (∼115 to 100 ka) was primarily driven by sea ice capping and a corresponding stratification of surface waters, which reduced the CO2 outgassing of upwelled carbon-rich waters. This hypothesis is supported by several lines of evidence, including (1) sea salt sodium (ssNA) archived in Antarctic ice cores, suggesting sea ice expansion near the Antarctic continent (Wolff et al., 2010); (2) δ15N proxy data from the central Pacific sector of the Southern Ocean, suggesting increased stratification south of the modern-day Antarctic Polar Front (Studer et al., 2015); and (3) diatom assemblages in the Permanently Open Ocean Zone (POOZ) of the Atlantic sector, suggesting a slight cooling and northward expansion of sea ice during MIS 5d (Bianchi and Gersonde, 2002). Our data address this hypothesis by providing insights into early sea ice expansion into the polar frontal zone of the western Pacific sector.

Our data show that, in contrast to the Atlantic sector (Bianchi and Gersonde, 2002), there does not appear to be any evidence of sea ice expansion in the southwestern Pacific during MIS 5d at either the TAN1302-96 or SO136-111 core sites (Fig. 5). Unfortunately, the lack of spatially extensive quantitative records extending back to Termination II limits our ability to estimate the timing and magnitude of sea ice changes for regions poleward of 59 S in the southwestern Pacific. We anticipate, however, that an advance in the sea ice edge, consistent with those outlined in Bianchi and Gersonde (2002), likely would have reduced local SST as the sea ice edge advanced closer to the core site. Indeed, the TAN1302-96 SSST record does show a decrease to ∼2C (observed at 108 ka), which quickly rebounded to ∼4C by ∼102 ka (Fig. 5). However, this SSST drop occurred roughly 7 ka after the initial CO2 reduction, suggesting that the CO2 drawdown event and local SSST reduction may not be linked. Thus, while we cannot rule out the possibility of modest sea ice advances or consolidation of pre-existing sea ice (particularly to the south of the core sites), the quantitative WSI and SSST reconstructions suggest that sea ice cover over our core site was limited during glacial inception.

Given that sea ice was not at its maximum extent during the early glacial, it stands to reason that any reductions to air–sea gas exchange in response to the hypothetically expanded sea ice would not have been at its maximum impact either. Previous modelling work has suggested that the maximum impact of sea ice expansion on glacial–interglacial atmospheric CO2 reductions ranged from 5 to 14 ppm (Kohfeld and Ridgwell, 2009). More recent modelling studies are consistent with this range, suggesting a 10 ppm reduction (Stein et al., 2020), while some studies even suggest a possible increase in atmospheric CO2 concentrations due to sea ice expansion (Khatiwala et al., 2019). Furthermore, Stein et al. (2020) suggest that the effects of sea ice capping would have taken place after changes in deep ocean stratification had occurred and would have contributed to CO2 drawdown later during the mid-glacial period. These model results, when combined with our data, suggest that even if modest sea ice advances did take place during the early glacial (i.e., MIS 5d), their impacts on CO2 variability likely would have been modest, ultimately casting doubt on the hypothesis that early glacial CO2 reductions of 35 ppm can be linked solely to the capping and stratification effects of sea ice expansion.

4.3 Other potential contributors to early glacial CO2 variability

The changes observed in WSIC and SSST from TAN1302-96 suggest that sea ice expansion was likely not extensive enough early in the glacial cycle for a sea ice capping effect to be solely responsible for early atmospheric CO2 drawdown. This leaves open the question of what may have contributed to early drawdown of atmospheric CO2. In terms of the ocean's role, we highlight three contenders: (1) a potentially non-linear response between sea ice coverage and CO2 sequestration potential, (2) links between sea ice expansion and early changes in global ocean overturning, and (3) the impact of cooling on air–sea disequilibrium in the Southern Ocean.

The first possible explanation considers that not all sea ice has the same capacity to facilitate or inhibit air–sea gas exchange. We previously suggested that because sea ice was not at its maximum extent during MIS 5d, the contribution of sea ice on CO2 sequestration would likely not be at its maximum extent either. However, this assumes a linear relationship between sea ice coverage and CO2 sequestration potential. We know that different sea ice properties, such as thickness and temperature, determine overall porosity, with thicker and colder sea ice being less porous and more effective at reducing air–sea gas exchange compared to thinner and warmer sea ice (Delille et al., 2014). It is therefore possible that if modest sea ice advances took place closer to the Antarctic continent (and were therefore not captured by TAN1302-96), they may have been more effective at reducing CO2 outgassing either by experiencing some type of reorganization or consolidation, or through a change in properties such as temperature or thickness. It is also possible that sea ice coverage over some regions leads to more effective capping, while in other regions sea ice growth contributes only to marginal reductions in air–sea gas exchange. This, theoretically, could point to a non-linear response between sea ice expansion and CO2 sequestration potential, and thus modest sea ice growth around the Antarctic continent could have contributed in part to the ∼35 ppm initial CO2 drawdown event. While this is theoretical and cannot be adequately addressed in this analysis, it is worthy of deeper consideration.

The second possible explanation involves changes in the global overturning circulation. Kohfeld and Chase (2017) previously examined the timing of changes in δ13C of benthic foraminifera solely from the Atlantic basin and observed that the largest changes in the Atlantic Meridional Overturning Circulation (AMOC) coincided with the mid-glacial reductions in atmospheric CO2 changes mentioned above. Subsequent work of O'Neill et al. (2021) examined whole-ocean changes in δ13C of benthic foraminifera and noted that the separation between δ13C values of abyssal and deep ocean waters – and therefore the isolation of the abyssal ocean – was actually initiated between MIS 5d and MIS 5a (114 to 71 ka). Evidence for early changes in abyssal circulation and reductions in deep-ocean overturning have also been detected in Indian Ocean δ13C records (Govin et al., 2009). More recently, Indian Ocean εNd records (Williams et al., 2021) have suggested that the abyssal ocean may have responded to sea ice changes around the Antarctic continent early in the glacial cycle, with colder and more saline AABW forming as sea ice expanded near the continent. If indications of an early-glacial response in the global ocean circulation in the Indo-Pacific are correct, these data may also point to an elevated importance of sea ice near the Antarctic continent in triggering early, deep-ocean overturning changes.

The third possible explanation involves changes in surface ocean temperature gradients in the Southern Ocean, and how they could influence air–sea gas exchange. Several recent studies have pointed to the importance of changes to air–sea disequilibrium as a key contributor to CO2 uptake in the Southern Ocean (Eggleston and Galbraith, 2018; Marzocchi and Jansen, 2019; Khatiwala et al., 2019). Khatiwala et al. (2019) suggested that modelling studies have traditionally underrepresented (or neglected) the role of air–sea disequilibrium in amplifying the impact of cooling on potential CO2 sequestration in the middle to high southern latitudes during glacial periods. They argue that when the full effects of air–sea disequilibrium are considered, ocean cooling can result in a 44 ppm decrease due to temperature-based solubility effects alone. They attributed this increased impact of SST to a reduction in sea-surface temperature gradients explicitly in polar mid-latitude regions (roughly between 40 and 60 north and south). If we compare the SST gradients in the southwest Pacific sector over the last glacial–interglacial cycle (Fig. 6), we see an early cooling response between MIS 5e–d corresponding to roughly half of the full glacial cooling, specifically in the cores located south of the modern STF. While they did not quantify them, Bianchi and Gersonde (2002) also described a weakening of meridional SST gradients between the Subantarctic and Antarctic Zones during MIS 5d in the Atlantic sector. Although this analysis is based on sparse data, our SSST reconstructions are consistent with the notion that surface ocean cooling, a weakening of meridional SST gradients, and changes to the overall air–sea disequilibrium could be responsible for at least some portion of the early CO2 drawdown. Further SST estimates from the region, and from the global ocean, are needed to substantiate this hypothesis.

Figure 6SST estimates from seven cores located in the southwestern Pacific. SST used were five-point averages (depending on sampling resolution) taken at MIS peaks and median dates in accordance with boundaries outlined in Lisiecki and Raymo (2005). Due to the complex circulation and frontal structures in the region, cores were plotted in ± distance from the average position of the modern STF. Cores used include SO136-003 (SSTs calculated from alkenones, Pelejero et al., 2006), FR1/94-GC3 (alkenones, Pelejero et al., 2006), ODP 181-1119 (PF-MAT, Hayward et al., 2008), DSDP594 (PF-MAT, Schaefer et al., 2005), Q200 (PF-MAT, Weaver et al., 1998), SO136-111 (D-MAT, Crosta et al., 2004), and TAN1302-96 (D-MAT; this study). The blue band represents the modern STF zone, while the red dotted line represents the southern shift in the STF during MIS 5e (Cortese et al., 2013).


4.4 Sea ice expansion and ocean circulation

Although the TAN1302-96 WSIC record suggests that sea ice was largely absent at the core site until the mid-glacial (∼65 ka), the observed changes in sea ice could have modulated regional fluctuations in Antarctic Intermediate Water (AAIW) subduction throughout the glacial–interglacial cycle. The annual growth and decay of Antarctic sea ice plays a critical role in regional water mass formation. Brine rejection results in net buoyancy loss in regions of sea ice formation, while subsequent melt results in freshwater inputs and net buoyancy gains near the ice margin (Shin et al., 2003; Pellichero et al., 2018). This increased freshwater input and buoyancy gain near the ice margin can hinder AAIW subduction, with direct and indirect impacts on both the upper and lower branches of the meridional overturning circulation (Pellichero et al., 2018).

Previous research has used δ13C in benthic foraminifera to track changes in the depth of the interface between AAIW and Upper Circumpolar Deep Water (UCDW) (Pahnke and Zahn, 2005; Ronge et al., 2015). Low δ13C values are linked to high nutrient concentrations found at depths below ∼1500 m in the UCDW, and higher δ13C values are associated with the shallower AAIW waters (Fig. 5). Marine sediment core MD97-2120 (45.535 S, 174.9403 E, core depth 1210 m) was retrieved from a water depth near the interface between the AAIW and UCDW water masses (Pahnke and Zahn, 2005). Over the last glacial–interglacial cycle, fluctuations in the benthic δ13C values from MD97-2120 suggest that the core site was intermittently bathed in AAIW and UCDW, and that the vertical extent of AAIW fluctuated throughout the last glacial–interglacial cycle. Ronge et al. (2015) used the δ13C values from MD97-2120 and other core sites to quantify the contributions of AAIW to the waters overlying MD97-2120 (%AAIW, Appendix D). These results suggest that during warm periods, MD97- 2120 exhibited more positive δ13C values, corresponding to higher %AAIW, while cooler periods exhibited more negative values, corresponding to lower %AAIW (Fig. 5). This suggests that during cooler periods, the AAIW-UCDW interface shoaled, reducing the total volume of AAIW and indirectly causing an expansion of UCDW (Ronge et al., 2015).

Our comparison between %AAIW and regional WSIC estimates suggest a strong link between the two (Fig. 5). Specifically, we observe that AAIW shoaled and UCDW expanded (i.e., %AAIW is low) during periods when sea ice expansion occurred. In contrast, during periods of low WSIC, a reduced seasonal sea ice cycle, and warmer summer sea surface temperatures (e.g., MIS 5e), %AAIW is observed to be high. This correlation supports the idea that increased concentrations of regional sea ice resulted in a substantial summer freshwater flux into the AAIW source region. This regional freshening likely promoted a shallower subduction of AAIW and a corresponding volumetric expansion of UCDW, which can be seen by the isotopic offset of the δ13C values between the reference cores, and also by the increased carbonate dissolution in MD97-2120 during glacial periods (Fig. 7) (Pahnke et al., 2003; Ronge et al., 2015). These findings directly link sea ice proxy records to observed changes in ocean circulation and water mass geometry.

Figure 7Schematic of changes in southwestern Pacific sector sea ice coverage and water mass geometry between interglacial and glacial stages. Panel (a) depicts interglacial conditions where sea ice coverage is minimal and freshwater input from summer sea ice melt is low. This lack of freshwater input allows AAIW to subduct to deeper depths and bath core MD97-2120, capturing the higher δ13C signature of the overlying AAIW waters. The AAIW-UCDW interface (red dashed line) is located beneath MD97-2120. CO2 outgassing is occurring as a carbon-rich Circumpolar Deep Water upwell near Antarctica. Panel (b) depicts glacial conditions where sea ice expansion has occurred beyond TAN1302-96, increasing brine rejection, and stabilizing the water column. As a result of the increased sea ice growth, subsequent summer melt increases the freshwater flux into the AAIW source region and increases AAIW buoyancy. This buoyancy gain shoals the AAIW-UCDW interface above core MD97-2120, causing the core site to be bathed in low δ13C UCDW. The shoaling of AAIW causes an indirect expansion of CDW, increasing the glacial carbon stocks of the deep ocean, while sea ice reduces CO2 outgassing via the capping mechanism.


In addition to its influence on regional freshwater forcing and AAIW reductions, these sea ice changes may also coincide with larger-scale deep ocean circulation changes. The most dramatic increases in winter sea ice observed in TAN1302-96 and SO136-111, along with changes in %AAIW, are initiated during MIS 4. These shifts also correspond to basin-wide changes in benthic δ13C values in the Atlantic Ocean that suggest a shoaling in the AMOC during MIS 4 (Oliver et al., 2010; Kohfeld and Chase, 2017). Changes in deep ocean circulation are also recorded in εNd isotope data in the Indian sector of the Southern Ocean (Wilson et al., 2015), suggesting extensive reductions in the AMOC during this period. Recent modelling literature (Marzocchi and Jansen, 2019; Stein et al., 2020) suggests that sea ice formation directly impacts marine carbon storage by increasing density stratification and reducing diapycnal mixing, especially in simulations where brine rejection is enhanced near the Antarctic continental slope and open ocean vertical mixing (and subsequent CO2 outgassing) is reduced (Bouttes et al., 2010, 2012; Menviel et al., 2012). These simulations suggest a resulting CO2 sequestration of 20–40 ppm into the deep ocean.

Taken collectively, the available data show that sea ice expansion, AAIW-UCDW shoaling, changes in the AMOC, and a decrease in atmospheric CO2 all occur concomitantly during MIS 4 (Fig. 5). It appears likely, therefore, that sea ice expansion during this time influenced intermediate water density gradients through increased freshening and consequent shoaling of AAIW, which may also have increased the efficiency of the carbon pump and increased CO2 uptake by phytoplankton (Sigman et al., 2021). This appears to have occurred while simultaneously influencing deep-ocean density, and therefore stratification, through brine rejection and enhanced deep water formation, which ultimately lead to decreased ventilation (Abernathey et al., 2016). These changes in ocean stratification, combined with the sea ice “capping” mechanism, appear to agree with both the recent modelling efforts (Stein et al., 2020) and observed proxy data and fit well within the hypothesis that mid-glacial CO2 variability was primarily the result of a more sluggish overturning circulation (Kohfeld and Chase, 2017).

5 Summary and conclusion

This study presents new WSIC and SSST estimates from marine core TAN1302-96, located in the southwestern Pacific sector of the Southern Ocean. We find that the WSIC remained low during the early glacial cycle (130 to 70 ka), expanded during the middle glacial cycle (∼65 ka), and reached its maximum just prior to the LGM (∼24.5 ka). These results largely agree with nearby core SO136-111 but display some differences in WSIC magnitude with E27-23. This discrepancy may be explained by differences in statistical applications and/or lateral sediment redistribution, although more analysis is required to determine the exact cause(s).

The lack of changes in SSST and the absence of winter sea ice over the core site during the early glacial suggests that the sea ice capping mechanism and corresponding surface stratification in this region is an unlikely cause for early CO2 drawdown, and that alternative hypotheses should be considered when evaluating the mechanism(s) responsible for the initial drawdown. More specifically, we consider the impact of changes in SSST gradients between ∼40 to 60 S and support the idea that changes in air–sea disequilibrium associated with reduced sea-surface temperature gradients could be a potential mechanism that contributed to early glacial reductions in atmospheric CO2 concentrations (Khatiwala et al., 2019). Another key consideration is the potentially non-linear response between sea ice expansion and CO2 sequestration potential (i.e., that not all sea ice is equal in its capacity to sequester carbon). More analyses are required to adequately address this.

We also observe a strong link between regional sea ice concentrations and vertical fluctuations in the AAIW-UCDW interface. Regional sea ice expansion appears to coincide with the shoaling of AAIW, likely due to the freshwater flux from summer sea ice melt increasing buoyancy in the AAIW formation region. Furthermore, major sea ice expansion and AAIW shoaling occurs during the middle of the glacial cycle and is coincident with previously recognized shoaling in AMOC and mid-glacial atmospheric CO2 reductions, suggesting a mechanistic link between sea ice and ocean circulation.

In conclusion, this paper has focused exclusively on sea ice as a driver of physical changes, but we recognize that these changes in sea ice will be accompanied by multiple processes that interact and compete with each other. Marzocchi and Jansen (2019) note that teasing apart the individual components of CO2 fluctuations is complicated because of interactions between sea ice capping, air–sea disequilibrium, AABW formation rates, and the biological pump. We recognize that these processes may not act independently and, as such, have contributed new data to help advance our collective understanding of the role of sea ice on influencing atmospheric CO2 variability on a glacial–interglacial timescale.

Appendix A: Age model and sampling depths

Table A1Radiocarbon dates taken from TAN1302-96. NDFB = not distinguishable from the background.

Download Print Version | Download XLSX

Table A2Tie points used in construction of the TAN1302-96 age model.

Download Print Version | Download XLSX

Table A3Sample depth and corresponding age. Diatom slides using Method 1 used even-numbered sediment samples (e.g., 10, 20, 30), while diatom slides using Method 2 used odd-numbered sediment samples (e.g., 53, 87).

* Indicates the sample was calculated based on linear sedimentation rates.

Download Print Version | Download XLSX

Appendix B: Diatom slide preparation comparison

Figure B1Results from diatom slide preparation methods 1 and 2. No notable differences or biases were observed between the two different methods.


Appendix C: TAN1302-96 and E27-23 comparison

Potential causes for WSIC estimate differences

The third potential cause for the observed differences between TAN1302-96 and E27-23 WSIC estimates is through the cumulative effects of different laboratory protocols. While it is difficult to determine precisely how much different laboratory protocols could influence the results, we cannot exclude this explanation as a possible contributor to differences in WSIC.

The fourth potential cause for differences in WSIC estimates between E27-23 and TAN1302-96 are differences in counting and identification methods. We believe this is an unlikely cause for the differences observed between E27-23 and TAN1302-96 primarily because of the magnitude of counting discrepancies required to cause a difference of 50 % WSIC estimates between the two cores. The close coupling of WSIC estimates between TAN1302-96 and SO136-111 over the entire glacial–interglacial cycle supports that a fundamental issue relating to taxonomic identification and/or methodology is an unlikely explanation for the observed WSIC differences.

Finally, the fifth potential cause of differing WSIC estimates is selective diatom preservation (e.g., Pichon et al., 1992; Ragueneau et al., 2000). The similarities between TAN1302-96 and SO136-111 WSIC estimates, along with independent indicators in cores E27-23 and TAN1302-96, suggest that this is unlikely. For E27-23, Bradtmiller et al. (2009) used the consistent relationship between 231Pa/230Th ratios and opal fluxes to suggest that dissolution remained relatively constant between the LGM and Holocene periods. In TAN1302-96, we assigned a semi-quantitative diatom preservation value between 1 (extreme dissolution) and 4 (virtually perfect preservation) for each counted specimen. The average preservation of diatoms for the entire core was 3.38±0.13, with no observed bias based on sedimentation rate or MIS. This assessment, although semi-qualitative, suggests that preservation remained relatively constant (and good) throughout TAN1302-96 and is therefore unlikely to cause large differences in WSIC between the two cores.

Figure C1Preliminary focusing factor (FF) values for E27-23. These results suggest notable lateral sediment redistribution over the last 26 ka, requiring further analysis (Bradtmiller et al., 2009).


Appendix D: %AAIW calculation

The calculation of %AAIW in this study is the same as was used in Ronge et al. (2015):


All core information for MD97-2120, MD06-2986, and MD06-2990 can be found through the original publication.

Data availability

All data has been published on Pangaea and can be found at (Jones et al., 2021).

Author contributions

Study conception and design was completed by KK and HB. Data collection was completed by JJ, KK, HB, XC, ML, GD, ZC, and AL. Data analysis and the interpretation of results was completed by JJ, KK, HB, XC, ZC, AL, HA, and GJ. Draft manuscript preparation and editing was completed by JJ, KK, HB, XC, GD, ZC, AL, HA, and GJ. All authors reviewed the results and approved the final version of the paper.

Competing interests

The contact author has declared that neither they nor their co-authors have any competing interests.


Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Special issue statement

This article is part of the special issue “Reconstructing Southern Ocean sea-ice dynamics on glacial-to-historical timescales”. It is not associated with a conference.


A special thanks to Rachel Meyne (Colgate University), who assisted with slide preparation, Maureen Soon (University of British Columbia), who assisted with opal concentration measurements, and Marlow Pellatt (Parks Canada), who assisted with project conceptualization and guidance. The TAN1302-96 core was collected during the TAN1302 RV Tangaroa voyage to the Mertz Polynya. We would like to thank the voyage leader Mike Williams and Captain Evan Solly and the crew, technicians, and scientists involved in the TAN1302 voyage.

Financial support

This research has been supported by the Canadian Natural Sciences and Engineering Research Council Grant’s Discovery Grant (grant no. RGPIN342251), the Australian Nuclear Science and Technology (ANSTO) grant (grant no. AP11676), and the Australian Research Council’s Discovery Projects funding scheme (grant no. DP180102357). Additional funding for travel and workshop collaboration was provided to Jacob Jones by a Past Global Changes (PAGES) grant (grant no. C-SIDE WS_163) to the Cycles of Sea Ice Dynamics in the Earth System (C-SIDE) working group.

Review statement

This paper was edited by Bjørg Risebrobakken and reviewed by Diana Krawczyk and two anonymous referees.


Abernathey, R. P., Cerovecki, I., Holland, P. R., Newsom, E., Mazloff, M., and Talley, L. D.: Water-mass transformation by sea ice in the upper branch of the Southern Ocean overturning, Nat. Geosci., 9, 596–601,, 2016.

Archer, D. E., Martin, P. A., Milovich, J., Brovkin, V., Plattner, G. K., and Ashendel, C.: Model sensitivity in the effect of Antarctic sea ice and stratification on atmospheric pCO2, Paleoceanography, 18, 1012,, 2003.

Barrows, T., Juggins, S., De Deckker, P., Calvo, E., and Pelejero, C.: Long-term sea surface temperature and climate change in the Australian–New Zealand region, Paleoceanography, 22, PA2215,, 2007.

Benz, V., Esper, O., Gersonde, R., Lamy, F., and Tiedemann, R.: Last Glacial Maximum sea surface temperature and sea-ice extent in the Pacific sector of the Southern Ocean, Quaternary Sci. Rev., 146, 216–237,, 2016.

Bianchi, C. and Gersonde, R.: The Southern Ocean surface between Marine Isotope Stages 6 and 5d: shape and timing of climate changes, Palaeogeogr. Palaeocl., 187, 151–177,, 2002.

Bereiter, B., Eggleston, S., Schmitt, J., Nehrbass-Ahles, C., Stocker, T.F., Fischer, H., Kipfstuhl, S., and Chappellaz, J.: Revision of the EPICA Dome C CO2 record from 800 to 600 kyr before present, Geophys. Res. Lett., 42, 542–549,, 2015.

Bostock, H., Hayward, B., Neil, H., Sabaa, A., and Scott, G.: Changes in the position of the Subtropical Front south of New Zealand since the last glacial period, Paleoceanography, 30, 824–844,, 2015.

Bouttes, N., Paillard, D., and Roche, D. M.: Impact of brine-induced stratification on the glacial carbon cycle, Clim. Past, 6, 575–589,, 2010.

Bouttes, N., Paillard, D., Roche, D. M., Waelbroeck, C., Kageyama, M., Lourantou, A., Michel, E., and Bopp, L.: Impact of oceanic processes on the carbon cycle during the last termination, Clim. Past, 8, 149–170,, 2012

Bradtmiller, L. I., Anderson, R. F., Fleisher, M. Q., and Burckle, L. H.: Comparing glacial and Holocene opal fluxes in the Pacific sector of the Southern Ocean, Paleoceanography, 24, PA2214,, 2009.

Butzin, M., Köhler, P., and Lohmann, G.: Marine radiocarbon reservoir age simulations for the past 50,000 years, Geophys. Res. Lett., 44, 8473–8480,, 2017.

Butzin, M., Heaton, T.J., Köhler, P., and Lohmann, G.: A short note on marine reservoir age simulations used in INTCAL20, Radiocarbon, 62, 865–871,, 2020.

Cefarelli, A. O., Ferrario, M. E., Almandoz, G. O., Atencio, A. G., Akselman, R., and Vernet, M.: Diversity of the diatom genus Fragilariopsis in the Argentine Sea and Antarctic waters: morphology, distribution, and abundance, Polar Biol., 33, 1463–1484,, 2010.

Cortese, G., Dunbar, G. B., Carter, L., Scott, G., Bostock, H., Bowen, M., Crundwell, M., Hayward, B. W., Howard, W., Martínez, J. I., Moy, A., Neil, H., Sabaa, A., and Sturm, A.: Southwest Pacific Ocean response to a warmer world: Insights from Marine Isotope Stage 5e, Paleoceanography, 28, 585–598,, 2013.

Crosta, X., Pichon, J.-J., and Burckle, L. H.: Application of modern analog technique to marine Antarctic diatoms: reconstruction of maximum sea-ice extent at the Last Glacial Maximum, Paleoceanography, 13, 284–297,, 1998.

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., 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, 101861,, 2020.

Delille, B., Vancoppenolle, M., Geilfus, N. X., Tilbrook, B., Lannuzel, D., Schoemann, V., Becquevort, S., Carnat, G., Delille, D., Lancelot, C., Chou, L., Dieckmann, G. S., and Tison, J. L.: Southern Ocean CO2 sink: The contribution of sea ice, J. Geophys. Res., 119, 6340–3655,, 2014.

Eggleston, S. and Galbraith, E. D.: The devil's in the disequilibrium: multi-component analysis of dissolved carbon and oxygen changes under a broad range of forcings in a general circulation model, Biogeosciences, 15, 3761–3777,, 2018.

Esper, O. and Gersonde, R.: New tools for the reconstruction of Pleistocene Antarctic Sea ice, Palaeogeogr. Palaeocl., 399, 260–283,, 2014a.

Esper, O. and Gersonde, R.: Quaternary surface water temperature estimations: New diatom transfer functions for the Southern Ocean, Palaeogeogr. Palaeocl., 414, 1–19,, 2014b.

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.

Fenner, J., Schrader, H., and Wienigk, H. (Eds.): Diatom Phytoplankton Studies in the Southern Pacific Ocean, Composition and Correlation to the Antarctic Convergence and Its Paleoecological Significance, Geologisch-Paläontologisches Institut and Museum der Universität Kiel, Germany, 757–813,, 1976.

Ferrari, R., Jansen, M. F., Adkins, J.F., Burke, A., Stewart, A. L., and Thompson, A. F.: Antarctic sea ice control on ocean circulation in present and glacial times, P. Natl. Acad. Sci. USA, 111, 8753–8758,, 2014.

Ferry, A. J., Crosta, X., Quilty, P. G., Fink, D., Howard, W., and Armand, L. K.: First records of winter sea ice concentration in the southwest Pacific sector of the Southern Ocean, Paleoceanography, 30, 1525–1539,, 2015.

Froelich, P. N.: Biogenic opal and carbonate accumulation rates in the Subantarctic South Atlantic: The late Neogene of Meteor Rise site 704, Proceedings of the Ocean Drilling Program, Scientific Results, 120, 515–549, 1991.

Fryxell, G. A. and Hasle, G. R.: The genus Thalassiosira: some species with a modified ring of central strutted processes, Nova Hedwigia Beihefte, 54, 67–98, 1976.

Fryxell, G. A. and Hasle, G. R.: The marine diatom Thalassiosira oestrupii: structure, taxonomy and distribution, Am. J. Bot., 67, 804–814, 1980.

Galbraith, E. and de Lavergne, C.: Response of a comprehensive climate model to a broad range of external forcings: relevance for deep ocean ventilation and the development of late Cenozoic ice ages, Clim. Dynam., 52, 653–679,, 2019.

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.

Gersonde, R., Crosta, X., Abelmann, A., and Armand, L.: Sea-surface temperature and sea ice distribution of the Southern Ocean at the EPILOG last Glacial Maximum–a circum-Antarctic view based on siliceous microfossil records, Quaternary Sci. Rev., 24, 869–896,, 2005.

Ghadi, P., Nair, A., Crosta, X., Mohan, R., Manoj, M. C., and Meloth, T.: Antarctic sea-ice and palaeoproductivity variation over the last 156,000 years in the Indian sector of Southern Ocean, Mar. Micropaleontol., 160, 101894,, 2020.

Govin, A., Michel, E., Labeyrie, L., Waelbroeck, C., Dewilde, F., and Jansen, E.: Evidence for northward expansion of Antarctic Bottom Water mass in the Southern Ocean during the last glacial inception, Paleoceanography, 24, PA1202,, 2009.

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., Chceddadi, R., David, F., Ponel, P., and Reille, M.: The climate of western Europe during the last Glacial/Interglacial cycle derived from pollen and insect remains, Palaeogeogr. Palaeocl., 103, 73–93,, 1993.

Hasle, G. R. and Syvertsen, E. E.: Marine diatoms, in: Identifying marine phytoplankton, edited by: Tomas, C. R., Academic Press, 5–385,, 1997.

Hayward, B. W., Scott, G. H., Crundwell, M. P., Kennett, J. P., Carter, L., Neil, H. L., Sabaa, A. T., Wilson, K., Rodger, J. S., Schaefer, G., Grenfell, H. R., and Li, Q.: The effect of submerged plateaux on Pleistocene gyral circulation and sea-surface temperatures in the Southwest Pacific, Global Planet. Change, 63, 309–316,, 2008.

Heaton, T., Köhler, P., Butzin, M., Bard, E., Reimer, R., Austin, W., Bronk Ramsey, C., Grootes, P., Hughen, K., Kromer, B., Reimer, P., Adkins, J., Burke, A., Cook, M., Olsen, J., and Skinner, L.: Marine20–The Marine Radiocarbon Age Calibration Curve (0–55,000 cal BP), Radiocarbon, 62, 779–820,, 2020.

Johansen, J. R. and Fryxell, G. A.: The genus Thalassiosira (Bacillariophyceae): studies on species occurring south of the Antarctic Convergence Zone, Deep Sea Research Part B. Oceanographic Literature Review, 32, 1050,, 1985.

Jones, J., Kohfeld, K. E., Bostock, H., Crosta, X., Liston, M., Dunbar, G., Chase, Z., Levente, A., Anderson, H., and Jacobsen, G.: Radiocabon, d18O, winter sea ice, and summer sea surface temperature estimates from sediment cores TAN1302-96 and SO136-111, PANGAEA [data set],, 2021.

Khatiwala, S., Schmittner, A., and Muglia, J.: Air-sea disequilibrium enhances ocean carbon storage during glacial periods, Science Advances, 5, eaaw4981,, 2019.

Kohfeld, K. E. and Chase, Z.: Temporal evolution of mechanisms controlling ocean carbon uptake during the last glacial cycle, Earth Planet. Sc. Lett., 472, 206–215,, 2017.

Kohfeld, K. E. and Ridgwll, A.: Glacial-Interglacial Variability in Atmospheric CO2, in: Surface Ocean-Lower Atmospheric Processes, edited by: Quéré, C. L. and Saltzman, E. S., Washington DC, USA,, 2009.

Kowalski, E. and Meyers, P.: Glacial–interglacial variations in Quaternary production of marine organic matter at DSDP Site 594, Chatham Rise, southeastern New Zealand margin, Mar. Geol., 140, 249–263,, 1997.

Lhardy, F., Bouttes, N., Roche, D. M., Crosta, X., Waelbroeck, C., and Paillard, D.: Impact of Southern Ocean surface conditions on deep ocean circulation during the LGM: a model analysis, Clim. Past, 17, 1139–1159,, 2021.

Lisiecki, L. E. and Raymo, M. E.: A Pliocene–Pleistocene stack of 57 globally distributed benthic δ18O records, Paleoceanography, 20, PA1003,, 2005.

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

Lougheed, B. C. and Obrochta, S. P.: Age-depth modelling in Matlab, Zenodo [code],, 2018.

Lougheed, B. C. and Obrochta, S. P.: A Rapid, Deterministic Age-Depth Modeling Routine for Geological Sequences With Inherent Depth Uncertainty, Paleoceanography and Paleoclimatology, 34, 122–133,, 2019.

Marzocchi, A. and Jansen, M. F.: Global cooling linked to increased glacial carbon storage via changes in Antarctic sea ice, Nat. Geosci., 12, 1001–1005,, 2019.

Menviel, L., Joos, F., and Ritz, S.: Simulating atmospheric CO2, 13C and the marine carbon cycle during the Last Glacial–Interglacial cycle: possible role for a deepening of the mean remineralization depth and an increase in the oceanic nutrient inventory, Quaternary Sci. Rev., 56, 46–68,, 2012.

Mix, A. C., Bard, E., and Schneider, R.: Environmental processes of the ice age: land, oceans, glaciers (EPILOG), Quaternary Sci. Rev., 20, 627–657,, 2001.

Morales Maqueda, M. A. and Rahmstorf, S.: Did Antarctic sea-ice expansioncause glacial CO2 decline?, Geophys. Res. Lett., 29, 11-1–11-3,, 2002.

Nelson, C. S., Hendy, C. H., Jarrett, G. R., and Cuthbertson, A. M.: Near-synchroneity of New Zealand alpine glaciations and Northern Hemisphere continental glaciations during the past 750 kyr, Nature, 318, 361–363,, 1985.

Oliver, K. I. C., Hoogakker, B. A. A., Crowhurst, S., Henderson, G. M., Rickaby, R. E. M., Edwards, N. R., and Elderfield, H.: A synthesis of marine sediment core δ13C data over the last 150 000 years, Clim. Past, 6, 645–673,, 2010.

O'Neill, C. M., Hogg, A. McC., Ellwood, M. J., Opdyke, B. N., and Eggins, S. M.: Sequential changes in ocean circulation and biological export productivity during the last glacial–interglacial cycle: a model–data study, Clim. Past, 17, 171–201,, 2021.

Pahnke, K. and Zahn, R.: Southern Hemisphere Water Mass Conversion Linked with North Atlantic Climate Variability, Science, 307, 1741–1746,, 2005.

Pahnke, K., Zahn, R., Elderfield, H., and Schulz, M.: 340,000-year centennial-scale marine record of Southern Hemisphere climatic oscillation, Science, 301, 948–952,, 2003.

Paterne, M., Michel, E., and Héros, V. : Variability of marine 14C reservoir ages in the Southern Ocean highlighting circulation changes between 1910 and 1950, Earth Planet. Sc. Lett., 511, 99–104,, 2019.

Pellichero, V., Sallée, J. B., Chapman, C., and Downes, S.: The Southern Ocean meridional overturning in the sea-ice sector is driven by freshwater fluxes, Nat. Commun., 9, 1789,, 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.

Prebble, J. G., Bostock, H. C., Cortese, G., Lorrey, A. M., Hayward, B. W., Calvo, E., Northcote, L. C., Scott, G. H., and Neil, H. L.: Evidence for a Holocene Climatic Optimum in the southwest Pacific: A multiproxy study, Paleoceanography, 32, 763–779,, 2017.

Ragueneau, O., Tréguer, P., Leynaert, A., Anderson, R. F., Brzezinski, M. A., DeMaster, D.J., Dugdale, R. C., Dymond, J., Fischer, G., François, R., Heinze, C., Maier-Reimer, E., Martin-Jézéquel, V., Nelson, D. M., and Quéguiner, B.: A review of the Si cycle in the modern ocean: recent progress and missing gaps in the application of biogenic opal as a paleoproductivity proxy, Global Planet. Change, 26, 317–365,, 2000.

Renberg, I.: A procedure for preparing large sets of diatom slides from sediment cores, J. Paleolimnol., 4, 87–90,, 1990.

Reynolds, R., Rayner, N., Smith, T., Stokes, D., and Wang, W.: An Improved In Situ and Satellite SST Analysis for Climate, J. Climate, 15, 1609–1625,<1609:AIISAS>2.0.CO;2, 2002.

Reynolds, R., Smith, T., Chunying, L., Chelton, D., Casey, K., and Schlax, M.: Daily High-Resolution-Blended Analyses for Sea Surface Temperature, J. Climate, 20, 5473–5496,, 2007.

Ronge, T. A., Steph, S., Tiedemann, R., Prange, M., Merkel, U., Nürnberg, D., and Kuhn, G.: Pushing the boundaries: Glacial/interglacial variability of intermediate and deep waters in the southwest Pacific over the last 350,000 years, Paleoceanography, 30, 23–38,, 2015.

Rutgers van der Loeff, M. M., Cassar, N., Nicolaus, M., Rabe, B., and Stimac, I.: The influence of sea ice cover on air-sea gas exchange estimated with radon-222 profiles, J. Geophys. Res.-Oceans, 119, 2735–2751,, 2014.

Schaefer, G., Rodger, J. S., Hayward, B. W., Kennett, J. P., Sabaa, A. T., and Scott, G. H.: Planktic foraminiferal and sea surface temperature record during the last 1 Myr across the Subtropical Front, Southwest Pacific, Mar. Micropaleontol., 54, 191–212,, 2005.

Schlitzer, R.: Interactive analysis and vizualization of geoscience data with Ocean Data View, Comput. Geosci., 28, 1211–1218,, 2005.

Schneider Mor, A., Yam, R., Bianchi, C., Kunz-Pirrung, M., Gersonde, R., and Shemesh, A.: Variable sequence of events during the past seven terminations in two deep-sea cores from the Southern Ocean, Quaternary Res., 77, 317–325,, 2012.

Shin, S. I., Liu, Z., Otto-Bliesner, B., Kutzbach, J., and Vavrus, S. J.: Southern Ocean sea-ice control of the glacial North Atlantic thermohaline circulation, Geophys. Res, Lett., 30, 1096,, 2003.

Sigman, D., Fripiat, F., Studer, A. S., Kemeny, P. C., Martínez-García, A., Hain, M. P., Ai, X., Wang, X., Ren, H., and Haug, G. H.: The Southern Ocean during the ice ages: A review of the Antarctic surface isolation hypothesis, with comparison to the North Pacific, Quaternary Sci. Rev., 254, 106732,, 2021.

Smith, R. O., Vennell, R., Bostock, H. C., and Williams, M. J.: Interaction of the subtropical front with topography around southern New Zealand, Deep-Sea Res. Pt. I, 76, 13–26,, 2013.

Sokolov, S. and Rintoul, S.: Circumpolar structure and distribution of the Antarctic Circumpolar Current fronts: 2. Variability and relationship to sea surface height, J. Geophys. Res., 114, C11019,, 2009.

Stein, K., Timmermann, A., Kwon, E. Y., and Friedrich, T.: Timing and magnitude of Southern Ocean sea ice/carbon cycle feedbacks, P. Natl. Acad Sci. USA, 117, 4498–4504,, 2020.

Stephens, B. B. and Keeling, R. F.: The influence of Antarctic sea ice on glacial–interglacial CO2 variations, Nature, 404, 171–174,, 2000.

Studer, A. S., Sigman, D. M., Martínez-García, A., Benz, V., Winckler, G., Kuhn, G., Esper, O., Lamy, F., Jaccard, S. L., Wacker, L., Oleynik, S., Gersonde, R., and Haug, G. H.: Antarctic Zone nutrient conditions during the last two glacial cycles, Paleoceanography, 30, 845–862,, 2015.

Sun, X. and Matsumoto, K.: Effects of sea ice on atmospheric pCO2: A revised view and implications for glacial and future climates, J. Geophys. Res., 115, G02015,, 2010.

Toggweiler, J. R.: Variation of atmospheric CO2 by ventilation of the ocean's deepest water, Paleoceanography, 14, 571–588,, 1999.

Warnock, J. P. and Scherer, R. P.: A revised method for determining the absolute abundance of diatoms, J. Paleolimnol., 53, 157–163,, 2015.

Weaver, A., Eby, M., Fanning, A. F., and Wiebe, E. C.: Simulated influence of carbon dioxide, orbital forcing and ice sheets on the climate of the Last Glacial Maximum, Nature, 394, 847–853,, 1998.

Wilks, J. V. and Armand, L. K.: Diversity and taxonomic identification of Shionodiscus spp. in the Australian sector of the Subantarctic Zone, Diatom Res., 32, 295–307,, 2017.

Williams, M. J.: Voyage Report TAN1302, Mertz Polynya, Tech. Rep., National Institute of Water and Atmospheric Research (NIWA), Wellington, 2013.

Williams, T. J., Martin, E. E., Sikes, E., Starr, A., Umling, N. E., and Glaubke, R.: Neodymium isotope evidence for coupled Southern Ocean circulation and Antarctic climate throughout the last 118,000 years, Quaternary Sci. Rev., 260, 106915,, 2021.

Wilson, D. J., Piotrowski, A. M., Galy, A., and Banakar, V. K.: Interhemispheric controls on deep ocean circulation and carbon chemistry during the last two glacial cycles, Paleoceanography, 30, 621–641,, 2015.

Wilson, K., Hayward, B. W., Sabaa, A. T., Scott, G. H., and Kennett, J. P.: A one-million-year history of a north-south segment of the Subtropical Front, east of New Zealand, Paleoceanography, 20, PA2004,, 2005.

Wolff, E. W., Barbante, C., Becagli, S., Bigler, M., Boutron, C. F., Castellano, E., de Angelis, M., Federer, U., Fischer, H., Fundel, F., Hansson, M., Hutterli, M., Jonsell, U., Karlin, T., Kaufmann, P., Lambert, F., Littot, G. C., Mulvaney, R., Röthlisberger, R., and Wegner, A.: Changes in environment over the last 800,000 years from chemical analysis of the EPICA Dome C ice core, Quaternary Sci. Rev., 29, 285–295,, 2010.

Short summary
We provide new winter sea ice and summer sea surface temperature estimates for marine core TAN1302-96 (59° S, 157° E) in the Southern Ocean. We find that sea ice was not consolidated over the core site until ~65 ka and therefore believe that sea ice may not have been a major contributor to early glacial CO2 drawdown. Sea ice does appear to have coincided with Antarctic Intermediate Water production and subduction, suggesting it may have influenced intermediate ocean circulation changes.