Significance New information from the ANDRILL-2A drill core and a complementary ice sheet modeling study show that polar climate and Antarctic ice sheet (AIS) margins were highly dynamic during the early to mid-Miocene. Changes in extent of the AIS inferred by these studies suggest that high southern latitudes were sensitive to relatively small changes in atmospheric CO 2 (between 280 and 500 ppm). Importantly, reconstructions through intervals of peak warmth indicate that the AIS retreated beyond its terrestrial margin under atmospheric CO 2 conditions that were similar to those projected for the coming centuries.

Abstract Geological records from the Antarctic margin offer direct evidence of environmental variability at high southern latitudes and provide insight regarding ice sheet sensitivity to past climate change. The early to mid-Miocene (23–14 Mya) is a compelling interval to study as global temperatures and atmospheric CO 2 concentrations were similar to those projected for coming centuries. Importantly, this time interval includes the Miocene Climatic Optimum, a period of global warmth during which average surface temperatures were 3–4 °C higher than today. Miocene sediments in the ANDRILL-2A drill core from the Western Ross Sea, Antarctica, indicate that the Antarctic ice sheet (AIS) was highly variable through this key time interval. A multiproxy dataset derived from the core identifies four distinct environmental motifs based on changes in sedimentary facies, fossil assemblages, geochemistry, and paleotemperature. Four major disconformities in the drill core coincide with regional seismic discontinuities and reflect transient expansion of grounded ice across the Ross Sea. They correlate with major positive shifts in benthic oxygen isotope records and generally coincide with intervals when atmospheric CO 2 concentrations were at or below preindustrial levels (∼280 ppm). Five intervals reflect ice sheet minima and air temperatures warm enough for substantial ice mass loss during episodes of high (∼500 ppm) atmospheric CO 2 . These new drill core data and associated ice sheet modeling experiments indicate that polar climate and the AIS were highly sensitive to relatively small changes in atmospheric CO 2 during the early to mid-Miocene.

Knowledge regarding Antarctic ice sheet (AIS) response to warming climate is of fundamental importance due to the role ice sheets play in global sea level change. Paleoenvironmental records from Earth’s past offer a means to examine AIS variability under climatic conditions that were similar to today and those projected for the next several decades (1, 2). In this respect, the early to mid-Miocene is a compelling interval to study as proxy reconstructions of atmospheric CO 2 , albeit uncertain, suggest that concentrations generally varied between preindustrial levels (PAL = 280 ppm) and values at or above 500 ppm (3⇓⇓⇓⇓⇓–9). Additionally, global mean surface temperature during peak Miocene warmth was up to 3–4° higher than today (10), similar to best-estimate temperatures expected by 2100 under the highest projected greenhouse gas concentration pathway (Representative Concentration Pathway 8.5) (1, 2). Finally, Miocene geography was similar to today (11), and major circum-Antarctic oceanic and atmospheric circulation patterns that dominate the modern Southern Ocean were well established (12, 13).

Much of our understanding regarding AIS history through the early to mid-Miocene comes from far-field records from deep ocean basins. Benthic oxygen and carbon isotope proxies for global paleoclimate suggest that early to mid-Miocene climate and glacial environments were highly variable (14⇓⇓⇓⇓–19). These records include evidence for major transient glacial episodes and sea level fall, intervals of relative ice sheet stability, and periods of climatic warmth with major ice sheet retreat and sea level rise. Furthermore, reconstructions from sedimentary sequences on the Marion Plateau, offshore northeast Australia (20), and New Jersey margin (21) suggest sea level varied by up to 100 m. Episodes of sea level maxima (+40 m) suggest loss of Antarctica’s marine-based ice sheets that, at present, occupy much of West Antarctica and large portions of East Antarctica (22), as well as substantial loss of mass from Antarctica’s terrestrial ice sheets. Episodes of maximum sea level fall (up to −60 m) suggest that the AIS occasionally grew and advanced across continental shelves.

Geological records proximal to Antarctica’s coastal margin provide direct evidence of past ice sheet variability in response to changing global climate. The ANDRILL (AND)-2A drill core, a 1,138-m-long stratigraphic archive of climate and ice sheet variability from the McMurdo Sound sector of the western Ross Sea (77°45.488′S, 165°16.613′E), was recovered by drilling from an ∼8.5-m-thick floating sea ice platform in 380 m of water, located ∼30 km off the coast of Southern Victoria Land (SVL) (Fig. 1A) (23). The drill core comprises lower Miocene to Quaternary glacial-marine strata deposited in the steadily subsiding Victoria Land Basin (VLB) (24). Paleogeography was broadly similar to today, although continental shelves in the Ross and Weddell seas were likely shallower (SI Text and Fig. S1). Recovered core sediments, and the proxies they contain, allow us to assess past ice sheet dynamics along a coastal margin influenced by ice flowing from East Antarctica and across the West Antarctic continental shelf. Through analysis of an integrated proxy environmental dataset, we derive a new environmental reconstruction and combine this with a suite of global environmental data to establish a history of AIS response to global climate events and episodes during the early to mid-Miocene. This integrated dataset allows us to evaluate key drivers of high latitude climate and ice sheet variability between 21 and 13 Ma. Outcomes from research reported here and in a companion ice sheet modeling study (25) suggest the AIS advanced across continental shelves during cold orbital configurations and retreated well inland of the coast under warm orbits. This large range of AIS variability occurred under a relatively low range in atmospheric CO 2 concentration (∼280–500 ppm) and indicates the Antarctic environment was highly sensitive during the early to mid-Miocene.

Fig. 1. (A) Map showing Ross Sea area (Inset) and AND-2A drill site (white box indicates approximate area for schematic reconstructions shown below). (B) Stratigraphic summary of lower 925 m of AND-2A (214.13–1,138.54 mbsf) showing 61 sedimentary cycles. Glacial proximity curve tracks relative position of the grounding line through ice-contact (I), ice-marginal (IM), ice-proximal (P), ice-distal (D), and open marine (O) environments. Continuously acquired datasets include magnetic susceptibility and Niobium (Nb) XRF-CS counts. CIA (curve and bar) indicates arid (<50, blue) and less arid (>60, green) conditions. Intervals of peak palynomorph concentration shown by orange boxes. Foraminifera assemblages include cold water/ice marginal benthic species (red circles) and cool water planktonic species (blue circles). Blue bars, sea ice diatoms; M, intervals with well-preserved molluscs. Sea water temperature estimates based on TEX 86 L (black circles) and Δ 47 (blue circles, open, less well-preserved specimens). Environmental Motif curve based on the proxy environmental dataset. (C) Schematic reconstructions of region around AND-2A showing likely conditions for each environmental motif (I–IV).

Fig. S1. (Top Left) Ross Sea regional topography (ice sheet removed). Red square indicates area of sketch at right, which shows major geographic and topographic features in the region of the AND-2A drill hole (dashed lines indicate Late Neogene volcanic features). (Bottom Left) Miocene topography is based on a linear interpolation between an Eocene/Oligocene reconstruction (49) and an isostatically rebounded version of BEDMAP-2 (22). Sketch at right shows reconstruction of early Miocene topography in the region of the drill site.

Summary Our analysis of the AND-2A drill core and synthesis with regional and global data show that the early to mid-Miocene Antarctic coastal climate was highly variable. During relatively short-lived intervals of peak warmth, summer land surface air temperature was at least 10 °C, tundra vegetation extended to locations 80 km inland (34, 40), surface water temperatures in the Ross Sea were between 6 °C and 10 °C, and the AIS retreated inland. During intermittent intervals of peak cold climate, vegetation vanished and the AIS grew and advanced into the marine environment, expanding across the continental shelf. Whereas glacial cycles were paced by orbital variability through the early to mid-Miocene, maximum ice sheet retreat occurred when atmospheric CO 2 was ≥500 ppm and maximum advance when CO 2 was ≤280 ppm (Fig. 2 and Fig. S5) (4). New numerical ice sheet simulations (25) also show that the Miocene AIS expanded across the continental shelf when atmospheric CO 2 was low (280 ppm) and retreated well inland of the coast when CO 2 was high (500 ppm). These ice sheet proximal data and model simulations support inferences from benthic deep sea records that suggest the global climate system and AIS were highly sensitive during the mid-Miocene (14⇓–16). These results are consistent with observations and numerical climate and ice sheet simulations based on the warm Pliocene (41⇓⇓–44), which indicate that sustained levels of atmospheric CO 2 >400 ppm may represent a stability threshold for marine-based portions of the West and East Antarctic ice sheets. Furthermore, outcomes from our complementary drill core analysis and ice sheet modeling indicate that Antarctica’s terrestrial ice sheets were vulnerable when atmospheric CO 2 concentrations last exceeded 500 ppm. Given current atmospheric CO 2 levels have risen above 400 ppm (45) and are projected to go higher (2), paleoclimate reconstructions such as this one for the early to mid-Miocene imply an element of inevitability to future polar warming, Antarctic ice sheet retreat, and sea level rise.

Methods and Materials Methods are presented in detail in SI Text and ref. 23. AND-2A was described using standard sedimentological techniques to produce detailed stratigraphic logs (27). An age model for the core (Fig. S2) uses magnetostratigraphy, biostratigraphy, 87Sr/86Sr dating of macrofossils, and 40Ar-39Ar ages on lava clasts and tephra layers to correlate rock units to the Global Polarity Timescale (46). Assemblages of fossil pollen, dinoflagellates, diatoms, foraminifera, and molluscs were used to constrain paleoenvironmental conditions. A standard suite of continuous physical properties was collected on whole and split core and in the borehole. Whole rock inorganic geochemical data were collected at high sampling resolution using an X-ray fluorescence core scanner (XRF-CS). Additional chemical data were collected from discrete bulk sediment samples at lower resolution to provide calibration points for near-continuous noninvasive sampling obtained via XRF-CS. Concentrations of Al 2 O 3 , Na 2 O, CaO, K 2 O, P 2 O 5 , and total organic carbon (TOC) from bulk sediment samples were used to calculate the CIA (47). Samples for TEXL 86 were prepared at Utrecht University and LC-MS analyses performed at the Royal Netherlands Institute for Sea Research. Samples for Δ 47 (clumped isotopes) were prepared and analyzed at the California Institute of Technology.

SI Text Miocene Paleogeography. Miocene paleogeography of the broader Ross Sea region is uncertain due to a relative paucity of geological data. However, reconstructions based on geodynamic modeling, backstripping, and erosion models indicate that much of West Antarctica was subaerial in the Eocene and early Oligocene (49), and, although subsidence history is unclear, seismic data suggest that structural ridges including the Coulman and Central highs formed relatively shallow marine ridges that were sites for ice cap growth during the late Oligocene-early Miocene. Geomorphologic data indicate the main outlet valleys dissecting the TAM were formed by rivers during the Eocene and were modified subsequently and/or incised via glacial processes under both temperate and polar climatic conditions (50). Exposed landforms and sediments preserved at high elevation are relict Miocene features that underwent little modification since middle Miocene time (50). Uplift has also been minimal since ∼14 Ma, with rates of 100–200 m/My during the Late Neogene (51). A 15-Ma paleogeographic reconstruction for the Ross Sea region (Fig. S1) is based on a linear interpolation between reconstructed topography at the E/O boundary (49) and an isostatically rebounded version of modern topography (22). Present-day geography in the vicinity of the AND-2A drillsite is dominated by the TAM to the west (Fig. S1), where mountain peaks in the Royal Society Range reach heights of >4,000 m. Uplift of this coastal mountain range began in the early Eocene (∼55 Ma) at an average rate of 100 m/My and with subsequent removal (denudation) of 4.5–5 km of rock and sediment along the coastal wedge (52). Major Neogene alkalic volcanic features including Mount Discovery, Minna Bluff, Mount Morning, and Ross Island occur to the south and east of the drill site. This volcanic landscape is relatively young with most of these volcanic edifices postdating the AND-2A record. However, rock outcrop on Mount Morning have been dated at 18.7 and 14.6 Ma and indicate that volcanic features associated with the McMurdo Volcanic Group occurred in the region of the Koetlitz Glacier during the early to Middle Miocene (53). Age Model. An age model for the core (Fig. S2) uses magnetostratigraphy (23), biostratigraphy (23), 87Sr/86Sr dating of macrofossils (54), and 40Ar-39Ar ages on volcanic clasts and tephra layers (55) to correlate rock units to the Global Polarity Timescale (46). Key events and tie points are summarized in Table S1. Biostratigraphic datums include first and last appearance of diatoms. The maximum and minimum age range reported for these datums are derived from the total, average, and hybrid range models of ref. 56, but we use the hybrid range model ages for tie points (marked by base of arrows in Fig. S2 and included in the text below). The section from 1,138.54 (the base of AND-2A) to 986.94 mbsf is characterized by normal polarity (magnetozone N8). A radiometric date of 20.01 ± 0.12 Ma on pumice at 1,093 mbsf ties this normal interval to C6An.1n and the magnetic polarity reversal (MPR) at 986.94 mbsf to C6An.1n/C6r. Unconformity U1 occurs at the base of sedimentary cycle 14 and in the middle of magnetostratigraphic unit R7. We correlate the MPR at 958.59 mbsf (6.84 m above U1) to the C6r/C6n boundary based in part on radiometric dates of 19.61 ± 0.17 and 9.56 ± 0.12 Ma on pumice at 953.94 and 953.28 mbsf, respectively, that tie magnetozone N7 to C6n. This correlation suggests that ∼200 ky is missing in U1. A radiometric age of 18.82 ± 0.15 Ma on pumice at 831.66 mbsf is consistent with our correlation of N7 with C6n and indicates that the MPR at 783.69 mbsf is the C6n/C5Er boundary (18.748 Ma). U2 is located 8.75 m above the C6n/C5Er boundary at the base of cycle 22, which also contains the first occurrence (FO) of diatom Thalassiosira praefraga [first appearance datum (FAD): 18.58] and is characterized by a reversed polarity zone (R6). The last occurrence (LO) of T. nansenii (LAD: 17.5 Ma) at 763.84 mbsf indicates that sediments above this datum must be younger than ∼17.5 Ma and suggests that the MPRs at 759.25 and 749.4 mbsf are correlated to C5Dr.1n/C5Dr.2r (17.74 Ma) and C5Dr.1r/C5Dr.1n (17.717 Ma) respectively. We note that this correlation requires that specimens of T. praefraga at 736.2 mbsf and pumice clasts at 709 mbsf have been reworked upcore. The diatom specimens occur in the middle of diamictite unit and pumice clasts occur in a sandstone unit. Reworking is certainly plausible. Furthermore, alternate solutions to the age model require downhole reworking of younger specimens, which is highly unlikely given the nature of this compact drill core. Therefore, we conclude the correlation presented herein is the most parsimonious. We correlate the MPR at 722.45 mbsf with C5Dn/C5Dr.1r (17.533 Ma) based on the radiometric age of 17.5 ± 0.14 Ma on pumice from a tuff layer at 640.13 mbsf and infer that the dominantly normal magnetic stratigraphy between 722.45 and 581.34 mbsf correlates with C5Dn. Therefore, the MPR at 581.34 mbsf correlates with C5Cr/C5Dn (17.235 Ma). This solution is supported by the radiometric age of 17.2 ± 0.17 Ma on a lava clast at 564.92 mbsf and the occurrence of reversed magnetozones (R3.1 and R3.2) between 581.34 and 482.69 mbsf that likely correlate to the long reversed chron C5Cr (17.235–16.721 Ma). We infer that the C5Cn.3n/C5Cr boundary occurs at 482.69 mbsf but note that the boundary may lie anywhere between 482.69 and 460.62 mbsf as magnetostratigraphic measurements through this interval are ambiguous. The interval from 460.62 and 448.59 mbsf is characterized by normal polarity (magnetozone N3) and is correlated to C5Cn.3n. We tie the MPR at 448.59 to C5Cn.2r/C5Cn.3n (16.543 Ma). Reversed polarity interval R2.2 is tied to C5Cn.2r (16.472–16.543 Ma), and we correlate the ambiguous magnetozone from 439.23 to 429.56 with normal Chron C5Cn.2n (16.303–16.472 Ma) and reversed magnetozone R2.3 with C5Cn.1r. This solution is supported by the radiometric date of 16.64 ± 0.37 Ma on a lava clast at 440.83 mbsf and the FO of Denticulopsis maccollumii (FAD: 16.77 Ma) at 434.35 mbsf. This solution requires that the MPR at 412.71 mbsf correlates with C5Cn.1n/C5Cn.1r (16.268 Ma). We therefore infer that magnetozones N2.1A and B and N2.2 that occur between 412.71 and 328.52 mbsf, and the ambiguous magnetozones between the normal units, correlate with C5Cn.1n (15.974–16.268 Ma). This correlation is supported by strontium dates of 16.05 ± 0.35 Ma on an Adamusium shell fragment at 366.8 mbsf and radiometric date of 16.01 ± 0.17 Ma on a lava clast at 358.11 mbsf. Reversed polarity interval R2.0A and B must therefore correlate with C5Br but the exact location of the C5Cn.1n /C5Br boundary cannot be identified as the magnetic data are ambiguous. The absence of Actinocyclus ingens from the section below U2 suggests that the sediments immediately below the unconformity are at least older than 15.52 Ma (hybrid range model age for the FAD of A. ingens) and possibly older than 15.96 Ma (total range model age for the FAD of A. ingens). If we construct a line of correlation that accommodates all of the constraints outlined above, the C5Cn.1n/C5Br boundary would occur at the base of R2.0B. Two major unconformities (U4 and U3) bracket the top 48.44 m of the section described in this paper. Rocks below the U4 disconformity contain the LO of diatom Nitzschia sp. 17 (230 mbsf, LAD: 14.94–13.7 Ma) and FO of Actinocyclus ingens var. nodus (237.03, FAD: 14.42–14.57 Ma), which indicate the sequence is at least older than 13.7 Ma. Magnetic polarity is uncertain through much of this section of AND-2A, although an interval of normal polarity occurs between 246.31 and 253.92 mbsf. We correlate this with C5ADn based on the FO of A. ingens var. nodus, which is located 9 m above the normal interval. We also infer that the slope of the LOC (sedimentation rate) for the sequence of rocks between U4 and U3 is similar to that defined immediately below U3. Therefore, the rocks most likely range in age from 14.6 to 14.4 Ma. U4 occurs between cycles 61 and 62 (214.13 mbsf). The absence of diatom taxa including A. ingens var. nodus (LAD: 12.57 Ma) and Denticulopsis lauta (LAD: 11.08 Ma) at 205 mbsf suggests that the overlying rocks are no older than 12.57 Ma and are possibly younger than 11.08 Ma. Radiometric ages of 11.363 ± 0.072 and 11.43 ± 0.46 Ma on crystals from reworked lava clasts at 127.50 and 129.96 mbsf, respectively, support this interpretation. Based on data described above, we infer the event(s) causing U4 must have occurred at or after 14.4 Ma (possible after 13.7 Ma based on LAD of Nitzschia sp. 17) and before 11 Ma. Environmental Proxies and Interpretation. A summary of the suite of data and environmental proxies obtained from the AND-2A core is presented in Fig. S3. Details regarding these data and the environmental interpretations generated from the integrated proxy dataset are summarized below. Sedimentology. Detailed sedimentological description and facies analysis provide the fundamental stratigraphic framework on which our paleoenvironmental interpretations are based. Nine nonvolcanic sedimentary facies are identified and range from diatomaceous siltstone to massive diamictite and matrix supported diamictite (26, 27), which reflect ice-distal to ice proximal (and occasional subglacial) settings, respectively. The Miocene lithologic succession is divided into sixty-one sequences on the basis of the repetitive vertical occurrence of characteristic facies that reflect cycles of advance and retreat of glaciers through the Transantarctic Mountains (27). Distinctive facies associations (26) and stratigraphic motifs (27) are distinguished based on diversity and abundance of different facies that reflect varying glacial environmental conditions. Whereas details differ between climatic interpretations based on the facies association scheme (26) and stratigraphic motif designations (27), respectively, the long-term patterns of warming and cooling are compatible. Both interpretations suggest that diamictite-dominated sequences in the lowermost 250 m of the core (1,138.54–901.54 mbsf) reflect cold, subpolar, ice-proximal conditions. A 50-m-thick sequence of thin stratified diamictite and abundant mudstone (1,042.55–996.69 mbsf) in the middle of the diamictite-dominated interval indicates a period of ice-distal conditions and inferred warmer climate (27). A 100-m-thick mudstone dominated interval (901.54–774.94 mbsf) has been inferred to reflect deposition within a hemipelagic shelf environment with minimal glacial influence and a warmer temperate climate (26, 27). However, in this paper we argue that additional proxy climate data indicate that this interval reflects a cold climate and that the sediments were deposited beneath semipermanent floating ice shelf or land-fast sea ice. An abrupt return to ice-proximal environments at the drill site is reflected in a 137-m-thick sequence (774.94–637.96 mbsf) that is dominated by massive and stratified diamictite. Passchier et al. (26) suggested that the facies associations between ∼259 and ∼648 mbsf are characteristic of open-shelf sedimentation under melt-water influenced settings and an inferred warmer climate. Stratified clast-rich diamictite units were deposited during glacial peaks and laminated and bioturbated mudstone intervals were deposited from sediment-rich plumes during interglacial periods. Massive and stratified diamictite dominated units above 259 mbsf indicate a return to ice-proximal environments and colder climate. Fielding et al. (27) differ somewhat in their facies classification in this section of the core. They identified a lithologically diverse composite sequence and mudrock-dominated units from 637.96 to 389.03 mbsf and suggested that this sequence was deposited under a high-latitude temperate glacial regime. Diamictite-dominated units reflect a return to a subpolar glacial regime at 389.03 mbsf and a transition to cold subpolar/polar conditions at 296.34 mbsf. Regardless of the minor differences in composite unit boundaries, both lithologic interpretations imply that the sequence between 640 and 390 mbsf records an extended interval of diminished glacial influence at the drill site and that a transition to colder more ice-proximal environments occurs between 390 and 259 mbsf. Most of the lithologic units in AND-2A contain gravel-sized clasts that were eroded from the TAM and transported to the drill site during advance of grounded ice and via melt-out from floating ice shelves and icebergs. The quantity and composition of clasts change in parallel with sedimentary facies and vary in abundance from their near absence (in cycles 9–12) to an interval with more than 1,000 clasts per meter. Changes in clast composition reflect the varying influence of different glacial drainages and have been used to infer periods of local vs. distal glacial sources (57). For example, the presence of foliated granodiorites and tonalities and medium-grade metamorphic rocks indicate provenance from basement rocks exposed in the Royal Society Range between the Blue and Koetlitz glaciers (57). In contrast, clast assemblages containing a distinctive monzogranite together with hornblende-biotite paragneiss and medium to high-grade metamorphic rocks are likely derived from rock exposure in the Britannia Range between the Carlyon and Byrd glaciers. Therefore, the distribution of distinctive basement clast assemblages in the AND-2A core reflect provenance from three regions including the Royal Society Range, the area between the Mulock and Skelton glaciers, and an area between the Carlyon and Darwin glaciers. Changes in relative heavy mineral abundance delineate three intervals within the early to middle Miocene section of AND-2A (58). Interval I (Fig. S3) comprises three highly variable samples that are inferred to reflect provenance from rocks exposed at locations spanning the central TAM south of Byrd Glacier to the flanks of Mount Morning immediately south of the drill site. Of particular relevance is the high abundance of carbonate minerals in a sample at 599.02 mbsf, which suggests that glaciers were eroding limestone exposed at the coast of the central TAM. Interval II (Fig. S3) is dominated by heavy minerals derived from Ferrar Dolerite and high grade metamorphic rock exposures that suggest a local source in the vicinity of the Koettlitz Glacier. The presence of kyanite at several discrete horizons also suggests that ice periodically delivered sediment derived from rocks exposed in the Nimrod Glacier catchment. Increased abundance of heavy minerals with a Ferrar Group “fingerprint” indicates that glaciers had retreated from the coastal margin and were actively eroding dolerite sills exposed inland (58). Interval III (Fig. S3) is characterized by a decrease in the relative abundance of minerals derived from the Ferrar Group and high-grade metamorphic rocks, an increase in Granite Harbor intrusive rocks, and the reappearance of carbonate minerals. These data suggest that the EAIS margin readvanced, and outlet glaciers were actively eroding rocks exposed along the TAM front closer to the present-day coastline. Paleontology. Whereas AND-2A is generally fossil poor, terrestrial flora and marine flora and fauna occur throughout the succession and are abundant in several discrete intervals (23). Assemblages of fossil foraminifera are dominated by calcareous benthic taxa. Intervals containing abundant and well-preserved forams extend to ∼800 mbsf, and biogenic carbonate is present at several deeper horizons (30). Assemblages contain taxa that indicate paleo-water depths between 50 and 150 m (although depths up to 400 m are also possible). Notably, the interval above ∼350 mbsf contains the first occurrence of the cold water planktonic species Neogloboquadrina pachyderma. Fossil diatoms generally occur in trace levels to the bottom of the core and are only abundant in four discrete intervals above 450 mbsf. The presence of Paralia sulcata between 430 and 434 mbsf indicates relatively shallow water depths (<100 m) at the drill site, which is consistent with foram data. Thalassionema nitzschioides is also recovered in this interval and indicates a cool water seasonally ice-free marine environment. A unique diatomite unit occurs between ∼310 and 312 mbsf. This interval is dominated by Chaetoceras resting spores and contains the first occurrence of sea-ice indicators, including Fragilariopsis truncata and Synedropsis cheethamii, and taxa associated with ice-edge environments. T. nitzschioides also occurs in the diatomite unit, which indicates that seasonally ice free conditions existed at or nearby the drill site. Upcore diatom-rich intervals at 231–237 and 283–294 mbsf consist of massive, clast poor diamictite that was deposited in ice-proximal environments. Diatom assemblages recovered from these intervals are similar to those in the diatomite, although abundances are generally lower. Marine palynomorph abundance and assemblage composition also fluctuates throughout the core (23). Many reworked Paleogene dinocyst specimens are present below ∼857 mbsf. In situ dinocysts occur in discrete intervals throughout the core but specimens are most abundant above 430 mbsf. This in-situ dinocyst assemblage is dominated by Operculodinim centrocarpum and Pyxidinopsis braboi (32). High concentrations of these taxa occur within a diatom-bearing muddy diamictite at ∼431 mbsf and peak concentrations occur in the diatomite at 312–310 mbsf. Operculodinium centrocarpum is an extant dinocyst that is distributed across a broad range of environments, but today it occurs in low abundance south of the Subantarctic Front. High numbers of both O. centrocarpum and P. braboi in these two distinct intervals in AND-2A indicate seasonally ice-free conditions (32). Abundance and diversity of pollen and spores also vary throughout the core, although assemblage composition is relatively constant and comprises bryophyte spores, pteridophyte spores, and both conifer and angiosperm pollen. In general, the spore-pollen assemblages represent mossy tundra vegetation dominated by shrub podocarps and Nothofagus with a variety of subordinate plants dependent upon site microclimate and physiography (23, 32, 33). These data suggest that AND-2A records a general climate that was cooler than the modern austral polar-alpine tree limit in which the January mean temperature rarely exceeds 10 °C. However, an increase in Podocarpites pollen and peak counts of Nothofagidites pollen occur within the biosiliceous-bearing muddy diamictite between 432 and 434 mbsf and the diatomaceous mudstone between 312 and 310 mbsf and is concomitant with the increase in diatom and dinocyst abundance described above. This significant increase in terrestrial pollen may reflect an interval of elevated atmospheric temperature (>10 °C January mean), which allowed proliferation of woody plants (23) and possibly growth of tree-like forms (31). Increases in abundance of these pollen taxa also occur in discrete intervals within the lower sections of the core, although a comprehensive study through this interval has yet to be completed. Finally, more than 600 marine macrofossil-bearing horizons occur throughout the core. Most of these fossil-bearing intervals are dominated by polychaete tubes but fragments of bryozoa, echinoderms, cirripeds, brachiopods, and a relatively diverse molluscan fauna are also present (23, 29). Well-preserved bivalves have been recovered from several discrete intervals and include an unusually diverse assemblage of six pectinid species that represent free-lying and attached forms. Of particular note is a “pectinid shellbed” at 430 mbsf, which contains the highest concentration of valves in the core (29). Most of the shells in this interval are crushed together at high angles, suggesting the deposit records a death assemblage that was likely transported from a shallow marine/subtidal environment and deposited in deeper water offshore. The unusually diverse pectinid assemblage reflects warmer temperatures in the Miocene than at present. Today, diverse large scallop assemblages are typically encountered in subpolar to temperate regions and do not occur further south than on the Magellanic shelf (57°S) and subantarctic waters around Campbell and Auckland Islands (55°S). Sea surface temperatures at these locations have an annual range between ∼5 °C and 7 °C warmer than in the Ross Sea today (29). Physical properties and inorganic geochemistry. A standard suite of continuous physical properties measurements were collected on whole and split core (23). Measurements of magnetic susceptibility reach their highest values (>200) within several relatively long intervals in the lower section of the core: 1,105–1,040, 990–950, and 850–775 mbsf. Shorter intervals occur from 990 to 980, 755 to 745, and 570 to 550, 480–470, and 445–435 mbsf. Most of the high susceptibility intervals generally correspond to units with high volcanic clast content. They also typically coincide with other proxy data that suggest cooler and drier climate including low CIA values and low TEX 86 L derived SWTs. These data suggest that intervals characterized by high magnetic susceptibility were deposited during episodes of concentrated glacial erosion of local volcanic rocks under cool/cold climatic conditions. Whole rock inorganic geochemical data were collected from the AND-2A core at high sampling resolution using an XRF-CS (23). Additional data were collected from discrete bulk sediment samples at lower resolution (1–3 m) to provide calibration points for the near-continuous noninvasive samples obtained via XRF-CS. Percentage composition of major elements including TiO 2 , Al 2 O 3 , CaO, K2O, Fe 2 O 3 , Ce, Cr, Ni, Nb, and Zr were determined for the entire core. These data provide an additional tool to identify likely source rocks for glacial-marine sediments in the core, which can help constrain glacial reconstructions. For example, Niobium content is significantly higher in McMurdo Volcanic Group rocks (up to 200 ppm) (59) than in basement rocks exposed in the TAM. Therefore, intervals of high Nb content in the core reflect increased proportions of sediment derived from McMurdo Volcanic Group rocks exposed in the Mount Morning/Koetlitz Glacier region. Inorganic geochemical data from bulk sediment samples has been used to calculate the CIA, which essentially provides a proxy for feldspar weathering but is generally used as a proxy for weathering of silicate rocks (47). The index is based on measured concentrations of Al 2 O 3 , Na 2 O, CaO, and K 2 O (on salt-free samples) and has been corrected for CaO associated with carbonate and phosphate minerals. Values <50 indicate that no, or only minimal, chemical weathering occurred, whereas CIA values >70 reflect significant chemical alteration (60). CIA values in AND-2A range between 46 and 68. CIA values determined from unweathered exposures of Ferrar Dolerite, McMurdo Volcanic Group, and Granite Harbor Intrusive Complex are typically around 50 or less. Therefore, intervals with low CIA values (<52) reflect a sediment source comprised of largely unweathered source rocks and likely indicate a cold and relatively dry climate with minimal time for source rocks to undergo chemical alteration. For instance, several intervals with low CIA values coincide with high magnetic susceptibility and high Nb XRF-CS counts. Together these data suggest that sediment was sourced from relatively unweathered McMurdo Volcanic Group rocks exposed in the region of Mount Morning. In contrast, the interval between ∼400 and ∼315 mbsf is characterized by low CIA values, magnetic susceptibility, and Nb concentration. These data indicate that sediment was derived from unweathered outcrop of Ferrar Dolerite and/or Granite Harbor Instrusives under cool and dry climatic conditions. Intervals with CIA values >60 are also typically associated with low magnetic susceptibility and Nb concentration, which likely reflect sediment derived from more weathered outcrops of non-MVG rocks under warmer climate with more moisture available for silicate weathering. Lipid biomarkers. Analyses of lipid biomarkers recovered from sedimentary sequences offer an ever-increasing range of proxy data for environmental reconstructions, which include sea surface water temperature (upper 200 m) (SWT), mean air temperature (MAT), and regional precipitation. Notably, proxies based on organic compounds have become a particularly important tool for high latitude reconstructions of SWTs through the Neogene as the sediments at these latitudes contain low amounts of carbonate material used in more traditional approaches. A lipid-based paleothermometer, TEX 86 (61, 62), was applied to 49 samples in the AND-2A record to reconstruct SSTs (Fig. S3 and Table S3). We used the low temperature (0–200 m depth) calibration TEX 86 L (63), which is best suited for the polar oceans and Antarctic waters as it takes into consideration the depth in the water column where Thaumarchaeota (the organisms on which TEX 86 paleothermometry is based) live year round (64). Sea surface (upper 200 m) water temperature (SWT) is derived using the following equation: SWT = 50.8 × TE X 86 L + 36.1 ( r 2 = 0.87 ; n = 396 ) . The calibration error is ±2.8 °C (62). Results using an alternative calibration to that described in ref. 33 are also shown in Table S3. This calibration produces results with a similar trend but slightly different absolute temperatures. Furthermore, SWTs calculated using the Bayesian calibration of ref. 65 produce values that are several degrees higher, but the downcore trend remains similar (data not shown). Samples with a branched vs. isoprenoid tetraether (BIT) index >0.3 (66) indicate high soil input and were not used to interpret SWTs. Samples with insufficient signal or poor peak resolution were discarded. We also consider the value of the sample at ∼968 msbf unreliable as the sample occurs at an unconformity. Calculated SWTs vary throughout the early to middle Miocene section of the core and range from −1.4 to 7.0 ± 2.8 °C. Several intervals reveal persistent cold temperatures (SSTs <0 to ∼2 °C), similar to those recorded in the Ross Sea today; these include an extended section from ∼939 to ∼700 mbsf. The warmest SST values between 6 °C and 7 °C occur in discrete intervals at 433, 311, and 285 mbsf. SSTs higher than 5 °C are rare below 500 mbsf and more common above this level. A summary of TEX 86 L SWT data for EMs II, III, and IV is presented in Fig. S4 and highlights distinct differences between each motif. Carbonate stable isotope ratios. The stable isotopic analysis of carbonates from sediments is widely used to constrain paleoenvironmental conditions. We measured carbonate δ18O, δ13C, and mass-47 (Δ 47 ) values in 37 samples of well-preserved to reasonably well-preserved carbonate macrofossils to determine SWT. Δ 47 values form the basis for a new kind of stable isotope paleothermometer for carbonate minerals that is rooted in the same statistical thermodynamics as Urey’s carbonate-water thermometer, but that circumvents the need to know the δ18O of water. Δ 47 or “clumped” isotope thermometry is based on the principle that the pairing or clumping of heavy isotopes into bonds with each other in molecules is a temperature-dependent phenomenon. In carbonates, rare heavy isotopes (13C, 18O) occur in a pool of light abundant isotopes. Carbonate clumped isotope thermometry examines the proportion of 13C and 18O that are bound to each other within the carbonate mineral lattice. The extent of ordering is predicted to be temperature dependent based on theoretical calculations, with disorder increasing as temperature increases. The basis for the clumping of these heavy isotopes into bonds with each other is thought to be the thermodynamically controlled exchange of stable isotopes among isotopologues of carbonate ion. The equilibrium constant for the reaction C 18 12 O 16 O 2 2 − + C 16 13 O 3 2 − = C 18 13 O 16 O 2 2 − + C 16 12 O 3 2 − forms the theoretical basis for carbonate clumped isotope thermometry, with the doubly substituted species (or heavy isotope clump) slightly more stable than the other isotopologues, and a progressively more random distribution of heavy isotopes among all possible isotopologues preferentially favored with increasing temperatures. So in an equilibrium precipitate, if the equilibrium constant for the above reaction is known and the abundance of all four isotopic species is measured, the temperature can be calculated independent of the isotopic composition of the water in which the carbonate grew. In addition, clumped isotope-derived temperatures can be combined with carbonate δ18O to estimate water δ18O, by applying an appropriate carbonate-water oxygen-isotope fractionation relationship. In practice, carbonate clumped isotope thermometry is based on analyses of 13C18O16O in CO 2 (which has a mass of 47 amu) that is produced by acid digestion of CaCO 3 . Fractionation during acid digestion is ∼0.2 per mil (‰) and is a predictable function of acid-digestion temperature. The abundance of the doubly substituted isotopologue 13C18O16O is reported using the variable ∆ 47 . ∆ 47 is defined as the ‰ enrichment of 13C18O16O above the amount expected for a random distribution of isotopes. The methods and instrumentation used are identical to those used in a previous study (67). Briefly, samples were digested at Caltech at either 25 °C in a two-legged reaction vessel or on a 90 °C common acid bath system with automated cryogenic and chromatographic purification. An acid digestion fractionation factor of 0.092‰ was used to correct 90 °C-reacted samples to 25 °C. Measurements were made on a Finnegan MAT 253 gas source mass spectrometer at Caltech with 1 SE precision as good as ±0.009‰. The total analysis time (including peak centering, background measurement, and pressure balancing) ranged from 4 and 8 h. Measurements were made to yield a 16-V signal for mass 44, with peak centering, background measurement, and pressure balancing before each acquisition. These conditions were selected to ensure stable ion currents and minimize time-dependent fractionation of sample and reference gas reservoirs. Each day we analyzed a heated gas composed of CO 2 with a stochastic distribution of isotopes between isotopologues. Gases with different bulk δ18O and δ13C ratios in quartz breakseals were heated to 1,000 °C for 2 h and then quenched at room temperature. These standard gases (or heated gases) were then purified and analyzed using the same protocol as sample gases. Our analysis were conducted before the proposition of a new reference frame for ∆ 47 measurements based on normalization to a reference gas equilibrated with water, and therefore samples were projected onto the absolute reference frame using carbonate standard data. Temperatures were calculated using a new mollusk calibration equation (67), which was generated using the same methods and instrumentation described herein. For calculations of equilibrium 18O/16O ratios in calcite (and to calculate δ18O-calcification temperatures), we used a published relationship from Kim and O’Neil (68) to describe calcite-water fractionation. This calibration line is used in this study as it spans a temperature range (−1 to 30 °C) appropriate for work on the Antarctic and includes measurements from modern specimens collected in shallow Southern Ocean sites including the Ross Sea and Ushuaia. This calibration was also generated using the same 90 °C common acid bath digestion and automated sample cleanup apparatus and mass spectrometer systems that were used to make measurements on the ANDRILL samples. Well-preserved specimens were identified based on visual appearance and previously published Sr iosotope data and elemental ratios (54, 69). In addition, we measured other samples that were visibly well preserved and yielded plausible seawater δ18O values (i.e., nonbrine water δ18O) (70). These samples were classified as possibly well preserved. Reconstructed water δ18O values for all samples range from 0‰ to 2‰ (Table S4). Miocene ∆ 47 -SWT estimates (Table S4) are similar to those based on paleoecology and TEX 86 paleothermometry and match published model results (31).

Acknowledgments We thank three anonymous reviewers. We acknowledge the Antarctic Geological Drilling (ANDRILL) Program, which acquired the unique geological section examined in this study. Scientific research was supported by New Zealand Ministry of Business Innovation and Employment Contracts C05X0410 and C05X1001, the US National Science Foundation (Cooperative Agreement 0342484 to the University of Nebraska–Lincoln), the Italian Antarctic Research Programme, the German Research Foundation, the Alfred Wegener Institute for Polar and Marine Research (Helmholtz Association of German Research Centres), and New Zealand Antarctic Research Institute Grant NZARI 2013-1.

Footnotes Author contributions: R.L., D.H., F.F., E.G., R.D., and D.P. designed research; R.L., D.H., F.F., F.S., R.T., H.v.E., E.G., G.K., A.T., R.D., C.F., B.F., M.O., D.P., S.S., F.T., S.W., V.W., G.A., K.P., T.P., M.T., and S.S.T. performed research; D.H., F.S., R.T., H.v.E., E.G., G.K., A.T., C.F., B.F., N.G., R.M., T.N., M.O., D.P., S.S., F.T., S.W., V.W., G.A., K.P., T.P., M.T., and S.S.T. analyzed data; R.L. and D.H. wrote the paper; N.G. contributed to interpretation of data and glacial dynamics; and R.M. and T.N. contributed to interpretation of data.

The authors declare no conflict of interest.

This article is a PNAS Direct Submission.

See Commentary on page 3419.

This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1516030113/-/DCSupplemental.