Holocene climate variations in the western Antarctic Peninsula

Introduction Conclusions References

The late Holocene (the last ∼ 3000 yr) was characterized by more variable temperatures and increased sea ice presence, accompanied by reduced local primary productivity likely in response to a shorter growing season compared to the early or mid-Holocene.The stepwise increase in annual sea ice duration over the last 7000 yr might have been influenced by decreasing mean annual and spring insolation despite an increasing summer insolation.We postulate that in addition to precessional changes in insolation, seasonal variability, via changes in the strength of the circumpolar Westerlies and upwelling activity, was further amplified by the increasing frequency/amplitude of El Ni ño-Southern Oscillation (ENSO).However, between 4000 and 2100 yr BP, the lack of correlation between ENSO and climate variability in the WAP suggests that other climatic factors might have been more important in controlling WAP climate at this time.

Introduction
The Antarctic seasonal sea ice cycle impacts ocean-atmosphere heat and gas fluxes (Anderson et al., 2009), the formation of deep and intermediate waters that participate in driving the thermohaline circulation (e.g.Orsi et al., 2002), biogeochemical cycles and local biological production (Marinov et al., 2006;Sarmiento et al., 2004;Toggweiler et al., 2006).The Antarctic seasonal sea ice cycle is strongly affected by recent warming, especially in the western Antarctic sector where the seasonal duration of sea ice has decreased by over 85 ± 20 days between 1979 and 2004 (Stammerjohn et al., 2008).In the WAP, global warming has had a stronger and faster effect than in any other region of Antarctica (Vaughan et al., 2003;Steig et al., 2009).For instance, the WAP has experienced a winter air warming trend of almost 6 • C since the 1950's (King et al., 1994), which have been associated with increasingly reduced duration of winter sea ice cover (Stammerjohn et al., 2008), and the widespread retreat of WAP ice shelves (Doake and Vaughan, 1991;Rott et al., 1996;Vaughan and Doake, 1996).Subsurface ocean warming of +0.3 • C (the 0-300 m upper layer) over the last century (Levitus et al., 2000) strongly affects ice shelf stability in some sectors of West Antarctica and has the potential to significantly increase global sea-level (King, 1994;Joughin and Alley, 2011;Rignot, 2006).Bentley et al. (2009) provide a comprehensive review of Holocene paleoclimate records from the Antarctic Peninsula, including data from ice cores, and lacustrine and marine sediment cores.Similarly, Taylor and Sjunneskog (2002) and Allen et al. (2010) summarized existing proxy data sets.These records, which show some regional variability in the timing and expression of the climate signal, document two warm periods, the first during the early Holocene and the second, during the "MidHolocene Hypsithermal" (Bentley et al., 2009).Recent marine geochemical data from the Palmer Deep (WAP) based on the TEX 86 proxy, similarly document a relatively warm early Holocene (Shevenell et al., 2011), with average sea surface temperature (SST) of 3.7 ± 2.2 • C between 11 800 and 9000 yr BP and some peaks reaching +10  by millennial-scale variability of SST superimposed on an overall cooling of 3-4 • C throughout the Holocene.More recently, a newly published ice core record from James Ross Island, Eastern Antarctic Peninsula (Mulvaney et al., 2012) also documents a warm early Holocene, with atmospheric temperatures approximately +1.3 ± 0.3 • C warmer than today.The ice core record revealed relatively cooling temperatures of ∼ 1 • C from ∼ 9200 to 5000 yr BP, with a slowdown or slightly warmer temperatures from ∼ 5000 to 3000 yr BP corresponding to the "MidHolocene Hypsithermal", before a second cooling step of ∼ 1 • C over the late Holocene.
Here, we provide a new multi-proxy record from a core located in the Palmer Deep Basin (WAP) on the continental open shelf, south of Anvers Island.We combine more traditional diatom assemblage data with more recent developments in multi-proxy analytical work, including a temperature record based on the well-calibrated TEX L 86 paleothermometer, and a high-resolution record of specific diatom biomarker lipids, the highly branched isoprenoid (HBI) alkenes, which are directly related to sea ice variability.Combining these records has enabled us to assess Holocene climatic changes in the central part of the WAP, to compare our new data with other records from the region, including those from the nearby ODP Site 1098, also in the Palmer Deep (Domack et al., 2001;Leventer et al., 1996Leventer et al., , 2002;;Shevenell and Kennett, 2002;Sjunneskog and Taylor, 2002;Taylor and Sjunneskog, 2002;Ishman and Sperling, 2002;Shevenell et al., 2001), and to test hypotheses presented previously to explain WAP climate fluctuations over the last 9000 yr BP (insolation and ENSO) (e.g.Shevenell et al., 2011).

Oceanographic setting
The northern WAP is influenced by the Southern Ocean Westerlies and the Antarctic Circumpolar Current (ACC) (Fig. 1).The ACC mainly transports Circumpolar Deep Water (CDW) (Sievers and Nowlin, 1984), which is a relatively warm water mass and com-Figures

Back Close
Full Deep Water with a salinity maximum (Rintoul et al., 2001).CDW is derived from a mixture of deep waters arising from the North Atlantic, Pacific and Indian Oceans as well as waters formed in the Antarctic region.Due to Ekman pumping, CDW upwells onto the shelf when strong Westerlies are displaced south and penetrates the WAP shelf via deep glacial troughs (Klinck et al., 2004).In addition, strong offshore winds during the sea ice retreat season cause the relatively warm and salty UCDW (1.5 , 34.6-34.73 psu (Moffat et al., 2008)) to mix with Antarctic Surface Water (AASW).AASW forms in the upper 100-150 m of the water column and is mostly controlled by ice formation and melting, as well as wind forcing (Dierssen et al., 2002).AASW differs from UCDW by its cold temperatures and relatively low salinity (between −1.8 and 1 • C, 33.0-33.7 psu) (Smith et al., 1999).During the ice-free season, the modified UCDW upwells heat and nutrients, which are important for biological productivity (e.g.Moffat et al., 2008).The vertical transport of relatively warm UCDW is believed to prevent high rates of sea ice production along the WAP shelf (Martinson et al., 2008).Satellite monitoring demonstrates large spatial and temporal variability of sea ice along the WAP (Stammerjohn and Smith, 1996) and significant changes in sea ice duration over the past several decades (Stammerjohn et al., 2008a;2008b).While mean winter sea ice limits are located between 60-63 • S, the average summer sea ice boundary lies between 67-69 • S, 3-4 • south of Anvers Island (Fig. 1).There is little ice coverage during the austral summer, between January and April, with a minimum in March (Stammerjohn and Smith, 1996).Sea ice grows during the austral autumn and reaches its maximum cover in August (Stammerjohn and Smith, 1996), before melting between September and November.During the spring sea ice melt season, when temperatures warm, injections of fresh water, from the melting sea ice, stratify the water column, trap nutrients and phytoplankton in the light-rich photic zone and promote widespread primary production dominated by spring bloom diatom assemblages (Annett et al., 2010).
Later in the summer season, brine-rich waters are mixed, by the Westerlies and associated extratropical cyclonic flow, with deeper and warmer upwelled waters, which renew the nutrient pool and promote a second phytoplankton bloom (Moffat et al., 2008).Figures

Back Close
Full  1).We also analyzed 19 sediment samples of the core drilled at the latter site for comparison with the record of Shevenell et al. (2011).Sedimentary structures and physical properties of the JPC-10 core are discussed in Domack et al. (2003) where they showed that sediments are not disturbed by turbidites.

Age model
The age model for JPC-10 was primarily based on the correlation of its magnetic susceptibility (MS) with the MS recorded at ODP Site 1098, which in turn has a robust agemodel derived from 54 14 C dates (Domack et al., 2001).The two MS records show remarkably similar trends throughout the Holocene (Fig. 2).Correlation of the MS records of the two cores was further controlled by two carbonate-based radiocarbon dates from JPC-10.These two 14 C dates at 1250 cm and 1323 cm were corrected using a 1230 yr regional age reservoir (Domack et al., 2001) and calibrated using the Calib 6.1.1 radiocarbon calibration program, thus corresponding to 8404 and 9040 yr BP, respectively.
The good correlation between the NBP9903 JPC-10 and ODP Site 1098 allows for the reconstruction of Holocene variability at a centennial scale.According to the stratigraphy used, the top of the JPC-10 core corresponds to ∼ 270 yr BP within the age model uncertainties.The sedimentation rate (SR) is 0.1-0.15cm yr −1 between 9000 and 8000 yr BP and again for the last 5000 yr BP, while during the intervening period, between 8000-5000 yr BP, the SR is higher, ranging between 0.2-0.3cm yr −1 .Figures

Back Close
Full

Biomarkers extraction
For isoprenoid glycerol dibiphytanyl glycerol tetraethers (GDGTs) and highly branched isoprenoids (HBI) alkenes, lipids were first extracted at the Laboratoire d'Oc éanographie et du Climat: Experimentations et Approches Num ériques (LOCEAN) using a mixture of 9 mL CH 2 Cl 2 /MeOH (2 : 1, v : v) to which internal standards (7 hexyl nonadecane, 9 octyl heptadecene and androstanol) were added.Several sonication and centrifugation steps were applied in order to properly extract the desired compounds.After drying with N 2 at 35 • C, the total lipid extract was fractionated over a silica column into an apolar and a polar fraction using 3 mL hexane and 6 mL CH 2 Cl 2 /MeOH (1 : 1, v : v), respectively.

TEX L 86
The TEX 86 (TetraEther Index of tetraethers with 86 carbon atoms) proxy is based on the relative distribution of Thaumarchaeotal lipids or GDGTs (Schouten et al., 2002).
For temperature reconstruction, we measured TEX L 86 at a centennial scale resolution.
TEX L 86 is the modified form of TEX 86 and is recommended for application in polar oceans (Kim et al., 2010).In contrast to a previous study in the WAP (Shevenell et al., 2011) which adapted the original TEX 86 proxy and performed a local calibration, we used TEX L 86 as a proxy for reconstructing paleotemperature at the JPC-10 core site.Correlation of the TEX 86 index with SST < 10 • C, i.e. in the polar oceans, shows that changes in TEX 86 are relatively minor with temperature and therefore this index is likely unsuitable for polar regions (Kim et al., 2008).Kim et al. (2010)  GDGT-1, GDGT-2, and GDGT-3 indicate GDGTs containing 1, 2, and 3 cyclopentane moieties, respectively.Since Thaumarchaeota are in low abundance in the Antarctic summer surface water (the ∼ 0-45 m layer of low salinity water mass) but more abundant in winter in the ∼ 45-105 m depth interval of cold, salty water (i.e. the summer remnant of the previous winter, surface-mixed layer) (Kalanetra et al., 2009), we used the TEX L 86 calibration established against the integrated 0-200 m water depth temperatures (see Kim et al., 2012): It is also worthwhile to note that Thaumarchaeota are found in Antarctic sea ice (Cowie et al., 2011), albeit in low abundance and below the sea ice (Alonso-S áez et al., 2012;Grzymski et al., 2012).This shows that their production may occur even below the freezing point of water (i.e.< 0 • C).The application of TEX L 86 on the JPC-10 core results in a 1.2 • C temperature estimate for the most recent sediment (see results), which is within the range of the present annual mean temperature between 0-200 m at this location (Martinson et al., 2008).Thus, this proxy seems to be appropriate for the reconstruction of temperature in sediment records from polar oceans with an associated residual standard error of ± 2.8 • C (Kim et al., 2012).In order to depict a regional temperature pattern in the WAP, we also analyzed 19 samples from core ODP 1098 (Hole B) for which Shevenell et al. (2011) already reported TEX 86 temperature estimates.

Back Close
Full eluted isocratically with 99 % A and 10 % B for 5 min, followed by a linear gradient to 16 % B in 45 min, where A = hexane and B = hexane/propanol (9 : 1, v : v).Flow rate was 0.2 mL min −1 .After each analysis the column was cleaned by back-flushing hexane/propanol (9 : 1, v : v) at 0.2 mL min −1 for 10 min.Detection was achieved using atmospheric pressure positive ion chemical ionization mass spectrometry (APCI-MS).
Fractional abundances of each isoprenoid GDGT were obtained by normalizing each peak area to the summed area of all six isoprenoid GDGTs.

Highly Branched Isoprenoids (HBIs)
The concentrations of di and tri-unsaturated HBI lipids (respectively, C 25:2 and C 25:3 HBI alkenes), which are diatom-specific biomarkers (Volkman et al., 1994), were determined at a decadal scale resolution, within age uncertainties, to reconstruct variations in past sea ice extent.It has been demonstrated that a di-unsaturated isomer (diene) is mainly synthesized by sea ice diatom species (Johns et al., 1999;Mass é et al., 2011) whilst a tri-unsaturated HBI (triene) is produced by some open water diatom species (Belt et al., 2000;Mass é et al., 2011).Recent studies have highlighted the potential for the application of HBIs as sea ice indicators in the Southern Ocean (Barbara et al., 2010;Denis et al., 2010;Mass é et al., 2011) and their coherent evolution with diatom assemblages (Barbara et al., 2010;Denis et al., 2010).Furthermore, HBI structures appear to be relatively resistant to degradation (Robson and Rowland, 1988) although they might undergo rapid sulfurization under highly anoxic conditions (e.g. in the Antarctic Ellis fjord; Sinninghe Damst é et al., 2007).While some sedimentary sulfur species were identified at certain depths in JPC-10 sediments, GC-MS analyses did not identify any organic sulfur compounds with an HBI carbon skeleton.We are therefore relatively confident that degradation is minimal throughout the core.HBIs were obtained from the apolar fraction by fractionation over a silica column using hexane as eluent following the procedures reported by Belt et al. (2007) andMass é et al. (2011).After removing the solvent with N 2 at 35 • C, elemental sulfur was removed using the TBA (Tetrabutylammonium) sulfite method (Jensen et al., 1977;Riis and Babel, 1999).The obtained hydrocarbon fraction was analyzed within an Agilent 7890A gas chomatograph (GC) fitted with 30 m fused silica Agilent J&C GC column (0.25 mm i.d., 0.25 µm film thickness), coupled to a Agilent 5975C Series mass selective detector (MSD).The oven temperature in the GC was initially programmed at 40 • C, then increased by 10 • C min −1 up to 320 • C where it was held for 12 min.Spectra were then collected using the Agilent MS-Chemstation software.Individual HBIs were identified on the basis of comparison between their GC retention times and mass spectra with those of previously authenticated HBIs (e.g.Johns et al., 1999) and peaks of the two compounds were integrated using the Mass Hunter software.Values are expressed as concentration relative to the internal standards.

Diatoms
Diatom distribution and abundance in core JPC-10 was defined at a centennial scale resolution using the Crosta and Koc ¸(2007) counting convention.Details regarding slide preparation and diatom identification are described in Crosta et al. (2004).Although the abundance of each diatom species was estimated as a percentage (%) of the Chaetoceros-free diatom assemblage, the total diatom abundance, Chaetoceros Introduction

Conclusions References
Tables Figures

Back Close
Full resting spores (CRS) and non-CRS assemblage absolute abundance were expressed as Mv g −1 (millions of valves per gram of sediment).Species recording less than 1 % in every given sample were removed to decrease the effect of noise in the dataset, with the exception of several benthic species, which were grouped in a special benthic group.Eucampia specific counts, evaluating the relative abundances of the Eucampia antarctica var.recta ("polar" variety) versus Eucampia antarctica var.antarctica ("subpolar" variety) and terminal versus intercalary valves were completed by counting 100 specimens per slide if abundance permitted, or by counting all specimens on a slide.
Here, we selected a number of diatom indicator species and groups of indicator species that reflect specific environments.In the Antarctic Peninsula, diatom species distribution follows a well-marked latitudinal gradient.CRS mostly develop in the northern part of the Antarctic Peninsula (Buffen et al., 2007;Pike et al., 2008), where the area has a longer ice-free season, and usually bloom during periods of most intense water column stratification (Crosta et al., 2007(Crosta et al., , 2008;;Denis et al., 2006;Leventer et al., 1996).CRS are usually associated with high productivity (Crosta et al., 2007(Crosta et al., , 2008;;Denis et al., 2006;Leventer et al., 1996Leventer et al., , 2002) ) and constitute most of the preserved diatom species found in the JPC-10 sediments accounting for more than 90 % of the total diatom species.The dominant CRS exerts a strong control on diatom absolute abundance.We therefore distinguished the total diatom abundance with and without CRS.We performed CRS-free counts in order to insure that potentially important variations in the relative abundance of more minor diatom indicator species were not dampened by the overwhelming presence of CRS (e.g.Allen et al., 2005).
Other environmentally diagnostic and numerically abundant diatom species include cold variety of Thalassiosira antarctica which grows at the sea ice edge when the duration of seasonal sea ice is greatest (Buffen et al., 2007;Pike et al., 2009).Fragilariopsis curta, combined here with the very similar species Fragilariopsis cylindrus, is found close to sea ice (Buffen et al., 2007;Pike et al., 2008) and is used to reflect sea ice distribution.This ecological association has been widely applied by a number of authors (Crosta et al., 1998;Gersonde and Zielinski, 2000;Leventer, 1998;Gersonde Introduction

Conclusions References
Tables Figures

Back Close
Full  , 2005;Allen et al., 2011;Collins et al., 2012).Currently, both F. curta and F. cylindrus are found to be more abundant in the southern WAP (Pike et al., 2008).Eucampia Antarctica, also common in JPC-10, has been documented to occur as two varieties, termed "sub-polar" and "polar" that are distinguished by valve symmetry (Fryxell, 1989;Fryxell and Prasad, 1990), and generally occur either north or south, respectively, of the Antarctic Convergence Zone.In addition, a paleoenvironmental proxy based on the ratio of terminal to intercalary valves of Eucampia has been used to record oscillations of winter sea ice (Kaczmarska et al., 1993;Whitehead et al., 2005).The abundance of Fragilariopsis kerguelensis and Thalassiosira lentiginosa, which presently reside off shore, and share a close relationship with the ACC, are used as indicators of open ocean conditions (Crosta et al., 2004), and might highlight summer sea-surface temperatures (Crosta et al., 2007).Finally, the sea ice related diatom group is compared to the so-called benthic group (composed of Amphora, Cocconeis sp., Gomphonema, Licmophora and Navicula sp.) which include cryophilic diatoms that live in near-coastal sea ice, and are commonly found during spring and summer blooms (Krebs, 1983).

Results
In JPC-10, the Holocene TEX L 86 -derived temperature record (0-200 water depth) reaches its maximum (> ± 3 • C) at the bottom of the core, ∼ 9000 yr BP (Fig. 3b).Following this peak, temperatures decreased to ∼ +1.5 • C, until ∼ 7000 yr BP, when a stepwise decrease occurred, followed by a declining trend.From ∼ 7000 to 4200 yr BP, temperatures declined slightly from +1.6 to +0.7 • C, before experiencing a slowdown between ∼ 4200 and 3000 yr BP and a second cooling phase, reaching a low (∼ +0. the general trends of both records (JPC-10 and ODP Site 1098) are in good agreement (Fig. 4b), both in terms of overall temperature range and in highlighting the abrupt decrease in temperature from ∼ 9000 to 8200 yr BP.BIT values varied between 0.01 and 0.02 for both marine records (Fig. 4c).The low BIT values indicated that the input of soil organic matter (Hopmans et al., 2004) to our core sites was very low and, thus, did not influence our TEX L 86 record (cf.Weijers et al., 2006).The present-day vertical temperature profile in the study area shows a surface mixed layer (∼ 30 m) with an average annual mean temperature of −0.The HBI record (the ratio of D/T) is characterized by a two-step increase.Prior to ∼ 7000 yr BP the ratio is uniformly low, close to 0, before increasing to values of 1-2 throughout the period of mid-Holocene cooling (Fig. 3).The ratio showed lower and more stable values (< 1) prior to and during the late Holocene (3700 to 1500 yr BP), before increasing again to higher values and a peak > 5 around 1000 yr BP, its highest values on average over the last 1200 yr BP.The mid-Holocene cooling and late Holocene patterns resemble each other (increasing from low to high values), albeit with greater D/T ratios during the late Holocene.
The total diatom and CRS absolute abundance evolved similarly over the last 9000 yr (Fig. 3).While abundance with values close to those found prior to 7000 yr BP.Absolute abundance of diatoms other than Chaetoceros (Non-CRS assemblage data) increased from 20-40 Mv g −1 from 7000 yr BP, which is consistent with the total diatom abundance data, but then remained relatively high until the most recent part of the record (Figs. 3 and  5).CRS represent 80-90 % of the total abundance during the early Holocene, and increased slightly to 90-95 % during the mid-Holocene, prior to declining to ∼ 70 % of the total abundance since 3000 yr (Fig. 5).The trend of F. curta and F. cylindrus showed a similar profile with low values, less than 10 % of the non-CRS total abundance, before 7000 yr BP, a net increase towards higher values during the mid-Holocene, reaching nearly 50 % at 5500 yr, followed by a decrease to values around 10-20 % over the last 1000 yr BP.In comparison, the cold variety of T. antarctica exhibited an overall increasing trend throughout the Holocene, rising from 2 to 12 %.T. antarctica relative abundance decreased slightly at around 3000 yr BP before experiencing a strong increase over the last 3000 yr BP.In contrast, E. antarctica was a significant component of the non-CRS total abundance between 9000 and 6000 yr BP (20-40 %) with a dominance of the sub-polar variety of the species, but then declined sharply at 6000 yr BP, representing a minor component throughout the remainder of the Holocene (mean 3.6 %).This decline was accompanied by the switch to a polar form of the species.In addition, the ratio of terminal/intercalary Eucampia valves was uniformly low until 6000 yr BP and then became higher and more variable.F. kerguelensis and T. lentiginosa occurred in relatively high abundances (20-30 %) before 7000 yr BP and progressively declined towards lower values, ranging between 10-20 %, over the course of the Holocene.This trend is opposite to the benthic group record, which revealed the lowest abundance during the early Holocene warmth, less than 3 %, before increasing to > 3 % from ∼ 7000 to 4000 yr BP and declining thereafter.

CPD Figures Back Close
Full  2011) (Fig. 4a).These differences are firstly due to the different indices used (TEX L 86 in our study and TEX 86 in the study of Shevenell et al., 2011) as well as the different calibrations.We used the more extensive core-top data set of Kim et al. (2010) for our study, which contains far more data points from polar oceans than the calibration established by Kim et al. (2008) and used in Shevenell et al. (2011).Furthermore, we did not calibrate the TEX 86 index against sea-surface temperature but the mean temperature integrated over the 0-200 m water depth range as outlined in Kim et al. (2012).Importantly, our ODP Site 1098 TEX L 86 record is fully consistent with the one from JPC-10 (Fig. 4b), thus suggesting that they represent a regional record of temperature variability.

The early Holocene
Between ∼ 10 000 and 7000 yr BP, TEX L 86 -derived temperatures were warmer than those during the mid and late Holocene (Figs. 3 and 4).With the exception of a cooling event centered ∼ 8200 yr BP, which might be related to a global cooling phenomenon (Barber et al., 1999;Ellison et al., 2006;Morril and Jacobsen, 2005;Rohling and P älike, 2005), temperatures varied from 1 • C to 3 • C, which is similar to modern observations under a warming climate (Martinson et al., 2008).Warm upper ocean temperatures are supported by the uniformly low relative abundance of the colder variety of T. antarctica (T1 variety) between 9000 and 7000 yr BP (Fig. 5).These relatively warm temperatures

CPD Figures Back Close
Full were probably accompanied by reduced seasonal sea ice presence as illustrated by the very low D/T ratio (< 03) (Fig. 3), the low relative abundances of diatoms groups associated with sea ice, including F. curta and F. cylindrus, the benthic diatom group, and the low ratio of terminal/intercalary Eucampia valves.The resulting prolonged open ocean conditions likely allowed for glacial melt and thus the production of stratification-related diatom species, such as E. antarctica and CRS, and open ocean species such as F. kerguelensis and T. lentiginosa (Fig. 5).
The early Holocene warm period (∼ 12 000-7000 yr BP) has been recorded throughout the AP from both marine (Allen et al., 2010;Heroy et al., 2008;Sjunneskog and Taylor, 2002;Taylor et al., 2001;Taylor and Sjunneskog, 2002) and ice core records (Mulvaney et al., 2012), which report warmer conditions and/or weaker sea ice presence than today.High temperatures in our reconstructed record coincided with maximum annual and spring insolation since the onset of the Holocene (Fig. 3).Warmerthan-modern winter and spring conditions would have limited sea ice persistence and prolonged the local sea ice-free season.Given the link between insolation and sea ice, we infer that the latter likely retreated early in the spring and advanced late in the autumn during the early Holocene warmth, as supported by diatom assemblages.This would be consistent with warmer summer temperatures, estimated from model simulations to be slightly higher (0.5 • C) than today (Renssen et al., 2005), and the reconstructed atmospheric temperatures about +1 • C warmer than present days (Mulvaney et al., 2012).
A longer and warmer summer season, alongside a weaker sea ice presence, would have also likely favored the southward migration of the Westerlies as previously hypothesized (Shevenell et al., 2011).The southward migration of the Westerlies might have in turn increased the upwelling of warm UCDW and its transport onto the shelf, thus leading to higher ocean temperatures of between 1 and 2 • C. In addition, increased winddriven coastal upwelling along with the reduced amount of near-shore sea ice cover might have enhanced the transport of nutrient-rich waters to the surface and hence allowed the renewal of the nutrient pool in the euphotic zone.Long seasonal periods of Introduction

Conclusions References
Tables Figures

Back Close
Full stratified conditions, followed by strong upwelling activity probably contributed in concert to stimulate diatom blooms during the early Holocene warmth along the WAP as recorded for instance here, in the Palmer Deep Basin.

The mid-Holocene
From ∼ 7000 to 3000 yr BP, the TEX L 86 -derived temperature record indicated a slow cooling from +2 to +0.5 • C.This trend parallels decreasing annual and spring insolation, and increasing summer insolation (Renssen et al., 2005).The transition to a colder upper water column was also reflected in the overall increase in the relative abundance of T. antarctica T1 species at this time.This period of transition was concomitant with a pronounced increase in the D/T ratio (up to 2) (Fig. 3), suggesting an increasing seasonal presence of spring sea ice.The increasing abundance of diatoms species diagnostic of sea ice presence (F.curta and F. cyclindrus, Eucampia and the benthic diatom group) further support this scenario (Fig. 5) and suggest that regional sea ice persisted for longer throughout the year between ∼ 7000 and 4200 yr BP.Reduced open ocean conditions are supported further by the decreasing abundance of open ocean diatom indicators, F. kerguelensis and T. lentiginosa.Proximity to a seasonally retreating sea ice edge likely promoted conditions conducive to CRS growth, also observed by Sjunneskog and Taylor (2002) in the Palmer Deep, through increased periods of sea ice melt-induced upper water column stratification.The strongly stratified conditions associated with an extended sea ice season generated higher productivity than permanently open ocean conditions as reflected in the total diatom abundance record.The mid-Holocene cooling ended with a shift in all proxy (for instance, increasing TEX L 86derived temperatures and D/T ratio, and decreasing total diatom abundance) records between ∼ 4200 and 3000 yr BP.

Conclusions References
Tables Figures

Back Close
Full our TEX L 86 -derived temperature record points to a gradual cooling trend during the mid-Holocene, probably as the result of colder annual atmospheric temperatures and shorter summer duration associated with colder and longer winter and spring seasons, in turn driven by decreasing winter and spring insolation (Fig. 3) It is therefore likely that decreasing insolation in the winter and spring resulted in the development of a more extensive sea ice season which restricted the duration of the summer season, despite increasing summer insolation.This interpretation is supported by modeled temperatures that indicate a relatively continuous cooling trend of the annual and spring temperatures of 1 and 2 • C, respectively, over West Antarctica between ∼ 9000 yr BP and today (Renssen et al., 2005).This is also consistent with recently reconstructed atmospheric temperatures that show a similar profile during this time interval (Mulvaney et al., 2012).
During the mid-Holocene, sea ice probably retreated later in the spring and advanced earlier in the autumn than previous Holocene periods.The melting of extensive sea ice would have led to stronger stratification of the water column, likely favoring the growth of Chaetoceros species and enhancing overall primary productivity in the Palmer Deep basin.A similar increase in marine primary productivity was also reported in Lallemand Fjord (Taylor et al., 2001;Shevenell et al., 1996) and Marguerite Bay (Allen et al., 2010), although a decrease in CRS was reported for the former region.This could be related to the more southerly and coastal location of Lallemand Fjord, which may have contributed to stronger sea ice development that prevented the occurrence of large blooms of Chaetoceros.The extension of sea ice might have displaced the Westerlies northward, and the resulting reduced upwelling of UCDW might have led to a diminished transport of warm waters masses over the WAP shelf.Due to the longer seasonal persistence of sea ice along the shore and the reduced effects of the katabatic and synoptic winds, decreased coastal upwelling would have acted to limit the intrusion of warm waters to the surface, thus further decreasing temperatures.Introduction

Conclusions References
Tables Figures

Back Close
Full

The late Holocene
The late Holocene was characterized by more variable TEX L 86 -derived temperatures, ranging between ± 0.3 • C and 2.5 • C with the maximum reached at ∼ 1000 yr BP (Fig. 3).The T. antarctica T1 data also exhibit large swings, superimposed onto a general cooling trend as suggested by the overall increased average relative abundance of this colder water variety of T. antarctica (Fig. 5).The most abrupt shift in the ratio of D/T (up to 5) occurred at ∼ 2000 yr BP and implies that the late Holocene corresponded to a period of increased sea ice presence (Fig. 3).This is further supported by the continued high relative abundances of F. curta and F. cylindrus, as well as the benthic diatom group (Fig. 5).Conversely, the typically open ocean diatoms F. kerguelensis and T. lentiginosa reached their lowest abundance during this period indicating the shortest seasonal open ocean period throughout the Holocene.Furthermore, the CRS and the total diatom abundance declined during this period (Fig. 5) a trend unexpected considering the extent of upper water column stratification anticipated from the melting of such extensive sea ice cover.It has been shown that high surface water stratification generally leads to high productivity and, especially, high CRS production in the western Weddell Sea (Buffen et al., 2007).It is therefore possible that decreased water column stratification due to wind action reduced diatom abundance during this period as already evidenced in Bransfield Strait (Heroy et al., 2008).
The late Holocene (∼ 3000-600 yr BP) has been extensively described as a period of drastic environmental conditions changes in the AP that might be related to ocean cooling (for example, Allen et al., 2010;Bj örk et al., 1996a;Domack and McClennen, 1996;Heroy et al., 2008;Taylor et al., 2001;Taylor and Sjunneskog, 2002;Sjunneskog and Taylor, 2002) and atmospheric cooling (Mulvaney et al., 2012).During this period, the HBI record from JPC-10 suggests the greatest sea ice duration since the onset Introduction

Conclusions References
Tables Figures

Back Close
Full We suggest that despite maximum summer insolation, summer duration was probably its shortest since the early Holocene and colder and longer winter and spring seasons would have likely promoted longer sea ice seasons, thus supporting previous work in East Antarctica (Pike et al., 2009).
In contrast to the icy environment implied by diatom assemblage and HBI data, the JPC-10 TEX-derived temperature record suggests relatively warmer conditions during the ice-free months, and thus a rather rapid transition from a cold spring to a warm summer.A regime shift of this nature would have likely hindered a long, stable period of stratification and curtailed the associated growing season, as evidenced by the decreasing trends in CRS and total diatom abundance.A shorter growing season in the WAP is further supported by similar results from Marguerite Bay (Allen et al., 2010).However, we believe it is unlikely that variations in seasonal insolation alone were responsible for the climatic changes observed over the WAP during the mid to late Holocene.For example, despite low and relatively stable annual and spring insolation since 2000 yr BP, HBI and sea ice-related diatom records continued to show large amplitude changes throughout this period.Furthermore, the simulated atmospheric temperatures indicate summer cooling and relatively stable winter and spring conditions (Renssen et al., 2005) which is consistent with the atmospheric cooling trend reconstructed in the eastern Antarctic Peninsula (Mulvaney et al., 2012), while our TEX L 86 -derived temperatures illustrate a slight annual warming (+1 • C) and the D/T ratio and diatom assemblages report winter and spring cooling.An additional explanation is therefore required.

The influence of ENSO in the WAP
An intriguing aspect of the D/T ratio from JPC-10 is the increased range of variability from ∼ 4500 to 4000 yr BP and again from ∼ 1200 yr BP to the top of the core.The was less variable and uniformly low, respectively.The question is not only what caused the changes in the average D/T ratio, but also the cause of higher frequency variability.Similar high climate variability has been found at low latitudes since ∼ 5000 yr BP and related to the intensification of centennial scale ENSO activity (Chazen et al., 2009;Conroy et al., 2008;Moy et al., 2002).Recent work additionally demonstrates the impact of ENSO on Antarctic climate (Kwok and Comiso, 2002;Turner et al., 2004;Yuan and Martinson, 2000;Shevenell et al., 2011).Strength and frequency of ENSO vary in line with changes in east-west temperature and pressure gradients along the equatorial Pacific.During El Ni ño events, greater heat flux in the eastern Pacific propagates to the Antarctic coastal regions in the Pacific Ocean sector through Rossby wave trains, resulting in warmer conditions and reduced sea ice (Yuan et al., 2004).Conversely, during La Ni ña events, when temperatures are much lower in the eastern than western equatorial Pacific, heat flux is diminished and colder conditions are correlated with more sea ice in the Pacific Ocean sector of Antarctica (Yuan et al., 2004).Temperatures and sea ice around Antarctica seem to be more sensitive to the cold La Ni ña events than the warm El Ni ño events and respond with greater amplitude to the former (Yuan and Martinson, 2001).However, these two climatic states also induce a shift in the Westerlies, with the circumpolar wind belt transposed further south during La Ni ña events, which would in turn cause an increase in the upwelling of warm UCDW, thereafter propagating the El Ni ño and La Ni ña signal into the warming WAP surface waters (Martinson et al., 2008;Willmott et al., 2010).Steig et al. (2012)  peak to peak correlation between the ENSO and HBI data is not as clear for the remainder of the record but generally higher D/T ratios are in phase with greater ENSO activity (Fig. 6).From ∼ 7000 to 4500 yr BP sea ice extent and variability demonstrate increasing trends, with mid-Holocene peaks in D/T ratios and variability observed between ∼ 4500 to 4000 yr BP.The ENSO records from Conroy et al. (2006) and Moy et al. (2002) similarly record increasing variability and frequency of ENSO at this time.
Between ∼ 4000 and 2100 yr BP, no strong sea ice variations are recorded in phase with low ENSO activity recorded at the Equator (e.g.Moy et al., 2002) and Peru margin (Makou et al., 2011).However, a phase of increasing ENSO frequency was conversely recorded in the Galapagos area (Conroy et al., 2008).Finally, over the last 2100 yr BP, increasing ENSO frequency concomitant with the appearance of stronger La Ni ña events (Makou et al., 2011) is coincident with the greatest variability in the D/T ratio.Given these data, we suggest, in line with Shevenell et al. (2011), that ENSO variability may have been one of the main drivers of the sea ice seasonal cycle at centennial to millennial timescales in the WAP, with increasing influence in the later Holocene.
In particular, we note the data of Das and Alley (2008) that documents increased melt frequency at Siple Dome through the late Holocene; they suggest the role of increased storm activity and the inland movement of maritime warm and moist air, perhaps related to ENSO variability.However, considering the lack of correlation between our sea ice data and ENSO variability between ∼ 4000 and 2100 yr BP, we further propose that other factors also linked to the Antarctic Dipole (the Southern Annular Mode) and other low-latitude climate oscillations (e.g.Indian or Atlantic Ocean Dipole), may also play a role in controlling sea ice in the Antarctic Peninsula.Clearly, the role of ENSO and tropical forcing of Southern Ocean climate needs greater investigation.
Our results indicate that sea ice persisted longer during phases of stronger La Ni ña events (the last 2100 yr) while upper ocean temperatures became higher as shown by the increasing D/T ratio TEX L 86 -derived temperature from JPC-10, respectively.During La Ni ña events, it is likely that a southernmost position of the Westerlies would be accompanied by a reinforcement of the upwelling of warm waters, which

CPD Introduction
Full would in turn affect the WAP surface waters and preclude thick and long seasonal sea ice cover (Martinson et al., 2008;Willmott et al., 2010).The HBI record is not consistent with such a scenario, suggesting instead that the pronounced presence of sea ice coincided with our reconstructed high TEX L 86 -derived temperatures.In contrast, increasing abundance of cold water diatom species (e.g.T. antarctica) indicated that colder conditions characterized the winter, potentially promoting pronounced sea ice extent.Such a pattern can be validated if we reasonably suppose that the increasing ENSO frequency associated with more La Ni ña events (Makou et al., 2011) might have led to an increase in seasonal variability with greater amplitude between the warm summer and cold winter seasons.Given the low annual insolation, higher ENSO frequency and stronger La Ni ña events during the past 2000 yr BP, the WAP climatic variability might have been governed by an alternation of "extremely" warm and cold temperatures, which would explain the contrasting relatively high temperatures and enhanced sea ice cover recorded in JPC-10 and reported elsewhere (Denis et al., 2010).Despite increased sea ice during winter and late spring, warmer temperatures and stronger wind activity during the summer could have led to more rapid sea ice retreat and a stronger upwelling of warm waters, which together may prevent long periods of the stratified conditions that are conducive to diatom and particularly CRS blooms.

Summary and conclusions
In this study, we examined the WAP climate variability over the last 9000 yr BP using a suite of proxies, including diatom assemblage data and biomarkers, to reconstruct the variations in sea ice temperature and diatom productivity, from a core located in the Palmer Deep Basin.Overall, the data document a cooling trend and gradually increasing sea ice extent throughout the Holocene, potentially in response to decreasing annual and spring insolation.This lends support to recent suggestions that decreasing regional annual and seasonal mean insolation were the main drivers of temperature Introduction

Conclusions References
Tables Figures

Back Close
Full and sea ice changes occurring in the Southern Hemisphere over the last climatic cycle (Huybers and Denton, 2009).
Early Holocene warmth corresponded to a longer sea ice-free season with warm upper ocean temperatures probably linked to the upwelling of warm UCDW.Coeval with declining insolation, the midHolocene was in contrast characterized by a cooling trend, longer seasonal sea ice duration and extensive marginal ice zone conditions, which in turn resulted in maximum regional productivity.During the late Holocene, particularly over the last 2000 yr BP, the WAP experienced a major climatic step, and was characterized by an exceptionally long sea ice season which persisted until late in spring season.
During the Holocene, we suggest the role of low-latitude climate forcing on sea ice and local oceanic conditions.When comparing lower latitude proxy records of ENSO frequency with our multi-proxy data, we suggest that increasing ENSO frequency and decreasing insolation might have together, amplified the seasonal contrast between a warm summer and a cold winter.An exception is made for the period between 4000 and 2100 yr BP when neither HBIs nor TEX L 86 records showed a strong correspondence with ENSO variability.We propose the operation of additional influencing factors related to the West and East Antarctic dipoles (SAM), northern high latitudes (e.g.AMOC) and/or other tropical/subtropical climate influences (e.g.Indian Ocean Dipole) expressed through oceanic and atmospheric teleconnections during this period.Thus, further high-resolution investigations around Antarctica, in particularly in the eastern sector, are required for a better understanding of the processes controlling WAP climate during this period of time.Introduction

Conclusions References
Tables Figures

Back Close
Full Age (years BP) Fig. 5. Chaetoceros-free diatom relative abundance (%) of selected environmental indicator species from core JPC-10.On the right panel, the grey curve indicates the ratio between polar and subpolar E. antarctica species, while in black is represented the terminal and intercalary valves ratio.
Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | records by analyzing sediment slices of the Jumbo Piston Core 10 (JPC-10), collected during the United States Antarctic Program Cruise 03 in 1999 by the R/V Nathaniel B. Palmer.NBP9903 JPC-10 core, 13.33 m long and mainly composed of diatomaceous mud, was obtained from the Palmer Deep basin (64 • 53 S, 64 • 12 W, 905 m water depth), south of the Anvers Island, close to the ODP (Ocean Drilling Program) Site 1098 (64 • 51 S, 64 • 12 W, 1010 m water depth) (Fig.
Discussion Paper | Discussion Paper | Discussion Paper | readdressed the relationship of GDGTs with SST leading to the definition of a new GDGT index, TEX L 86 TEX L 86 = log [GDGT-2] [GDGT-1] + [GDGT-2] + [GDGT-3Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | . Determination of the ratio of the diene and triene (D/T) HBI alkenes, combined with diatom assemblage data, provides an indication of the relative input of sea ice algae and open water phytoplankton to the sediment.This new proxy has recently been applied to the reconstruction of Antarctic paleoenvironments and used to monitor the onset of deglaciation in Prydz Bay (Barbara et al., 2010), and the sea ice evolution through the Holocene in the Ad élie Basin (Denis et al., 2010).Bacterial degradation is considered to have had minor effects given the good preservation of HBI over the last 60 000 yr (Collins et al., unpublished data) Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | et al.
3 • C) at ∼ 3000 yr BP.Although the last 3000 yr BP were characterized by variable temperatures with a mean of +1.0 • C, the temperature record exhibited a peak at ∼ 1000 yr BP, with values close to +2.5 • C. In comparison, our TEX L 86 record from ODP Site 1098 Hole B covering the last 12 000 yr BP revealed a temperature range from +3 • C to +0.3 • C (Fig. 4a).Although the temporal resolution of our ODP Site 1098 record is relatively low, Discussion Paper | Discussion Paper | Discussion Paper | 200 m water depth layer with an annual mean temperature of around +0.3 • C(Levitus and Boyer, 1994).Hence, the estimated temperatures of the upper most sediments analyzed (+1.2 • C for JPC-10 and +1.0 • C for ODP1098) are well above the regional freezing point of seawater (−1.7•C,Martinson et al., 2008) and are only +1.0 • C warmer than the annual mean temperature of the 0-200 m water layer.Nevertheless, care has to be taken in interpreting the absolute values of TEX L 86 -derived reconstructions due to the relatively large calibration error (± 2.8 • C) Relative changes in TEX L 86derived temperature can be viewed with more confidence.
the early Holocene warmth corresponded to low diatom and CRS abundance, around 200 Mv g −1 , the mid-Holocene cooling was represented by the highest diatom abundance with two peaks of 1200 and 900 Mv g −1 at around 4800 and 4000 yr BP, respectively.The Late Holocene was marked by a decline in total diatom Figures Back Close Full Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | JPC-10 and ODP Site 1098 Hole B (Fig.4b) show a consistent picture of temperature variation in Palmer Deep Basin.However, our results give different temperature estimates in terms of amplitude and variations to the TEX 86 -SST values published byShevenell et al. ( Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | applicability of this proxy as an indicator of sea ice suggests greater variability of sea ice extent during those time intervals, and in contrast a lesser variability during the periods ∼ 7000-4500 yr BP, ∼ 4000-1200 yr BP and ∼ 9000 to 7000 yr BP, when the D/T ratio Figures Back Close Full Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | also highlight the potential importance of forcing from the Tropical Pacific in regulating the flow of UCDW onto the West Antarctic continental shelf.Here we compare the sea ice-related HBI data from the Palmer Deep, with a suite of ENSO proxy records, including those of Conroy et al. (2008) from the Galapagos, Moy et al. (2002) from the Ecuadorian Andes and Makou et al. (2011) from ODP Site 1228D on the Peruvian continental margin (Fig. 6).All the records document minimum frequency of ENSO events in the early Holocene, coincident with the time of the most uniform, and lowest D/T ratios, indicative of uniformly reduced sea ice cover.The Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |

Fig. 1 .
Fig. 1.Location of the JPC-10 core and the ODP Site 1098 in Western Antarctic Peninsula (WAP).The dashed yellow line represents the limit of sea ice presence during the summer while the solid green line delineates the maximum sea ice extent during the winter.Major oceanic currents: ACC, Antarctic Circumpolar Current and APCC, Antarctic Peninsula Circumpolar Countercurrent which is derived from the Weddell Sea surface waters.