Quaternary Science Reviews xxx (2012) 1e10 Contents lists available at SciVerse ScienceDirect Quaternary Science Reviews journal homepage: www.elsevier.com/locate/quascirev Ice sheet retreat dynamics inferred from glacial morphology of the central Pine Island Bay Trough, West Antarctica Martin Jakobsson a, *, John B. Anderson b, Frank O. Nitsche c, Richard Gyllencreutz a, Alexandra E. Kirshner b, Nina Kirchner d, Matthew O’Regan e, Rezwan Mohammad a, Björn Eriksson a a Department of Geological Sciences, Stockholm University, 106 91 Stockholm, Sweden Department of Earth Sciences, Rice University, Houston, TX 77005, USA Lamont-Doherty Earth Observatory of Columbia University, Route 9 W, Palisades, NY 10964, USA d Department of Physical Geography and Quaternary Geology, Stockholm University, 106 91 Stockholm, Sweden e School of Earth and Ocean Sciences, Cardiff University, Wales, United Kingdom b c a r t i c l e i n f o a b s t r a c t Article history: Received 26 September 2011 Received in revised form 23 December 2011 Accepted 27 December 2011 Available online xxx Pine Island Glacier drains portions of the West Antarctic Ice Sheet into the Amundsen Sea. During the Last Glacial Maximum the glacier extended nearly 500 km from its present location onto the outer continental shelf. Unusually restricted sea-ice cover during the austral summer of 2010 allowed for a systematic multibeam swath-bathymetric and chirp sonar survey of the mid-shelf section of Pine Island Trough. The mapped glacial landforms reveal new information about the paleo-Pine Island Ice Stream’s dynamic retreat from the mid-shelf area and confirm previous suggestion of a retreat in distinct steps. The periods of grounding line stability during the overall retreat phase are marked by sediment accumulations, i.e. grounding zone wedges. These wedges are here mapped in sufficient detail to characterize spatial dimensions and estimate the volume of deposited sediment. Considering a range of sediment flux rates from the paleo-Pine Island Ice Stream we estimate that the largest and most clearly defined grounding zone wedge, located at about 73 S in the surveyed area, took between 600 and 2000 years to form. The ice stream retreated landward of this wedge before 12.3 cal ka BP. The swath-bathymetric imagery of landforms in Pine Island Trough includes glacial features that suggest that retreat between periods of grounding line stability may be associated with episodes of ice shelf break-up. The depths of grounding line wedges decrease in a landward direction, from 740 to 670 m, and record elevation of the grounding line as it stepped landward. In all, the grounding line elevation varied by only w80 m over a distance of just over 100 km, implying a low ice sheet profile during retreat. Finally, we revisited seismic reflection profile NB9902, acquired along Pine Island Trough in 1999, in combination with the newly acquired swathbathymetric imagery from 2010. Together these data show that the ice stream paused during its retreat to form grounding zone wedges at an area in central Pine Island Trough where a high in dipping bedrock strata exists and the glacial trough is narrow, forming a bathymetric “bottle neck”. Ó 2011 Elsevier Ltd. All rights reserved. Keywords: West Antarctica Ice streams Grounding zone wedge Multibeam 1. Introduction The West Antarctic Ice Sheet (WAIS) comprises approximately 10% of the entire Antarctic Ice Sheet volume (Lythe et al., 2001). It is an ice sheet with its base mainly grounded below sea level, making it particularly sensitive to environmental change and potentially unstable. Sea level would rise more than 3 m if the marine-based part of WAIS disintegrated (Bamber et al., 2009). Slightly less than two thirds of WAIS drains into the Ross and RonneeFilchner Ice shelves * Corresponding author. Tel.: þ46 8 164719. E-mail address: [email protected] (M. Jakobsson). while the remaining part drains into the Amundsen Sea through Pine Island and Thwaites glaciers (Fig. 1) (Rignot et al., 2008). These latter two glaciers were considered by Hughes (1981) as the WAIS “weak underbelly”, referring to their direct exposure to the ocean and lack of protective physiographic barriers or large buttressing ice shelves. In addition, the bedrock of these glaciers deepens significantly inland, adding to their vulnerability (Holt et al., 2006; Vaughan et al., 2006). Recent observations indeed show rapid thinning and near constant acceleration of the Pine Island Glacier (PIG) (e.g. Joughin et al., 2003). The thinning is suggested to be caused by intrusion of warmer ocean water beneath the relatively small ice shelf in front of the grounded glacier (Jenkins et al., 2010; Jacobs et al., 2011). Inflow of warmer water has recently also been documented west of Pine Island Bay 0277-3791/$ e see front matter Ó 2011 Elsevier Ltd. All rights reserved. doi:10.1016/j.quascirev.2011.12.017 Please cite this article in press as: Jakobsson, M., et al., Ice sheet retreat dynamics inferred from glacial morphology of the central Pine Island Bay Trough, West Antarctica, Quaternary Science Reviews (2012), doi:10.1016/j.quascirev.2011.12.017 2 M. Jakobsson et al. / Quaternary Science Reviews xxx (2012) 1e10 Fig. 1. Overview map showing tracklines of IB Oden during expedition OSO0910 (black lines), the general bathymetry of Pine Island Bay continental shelf (Nitsche et al., 2007) and locations of grounding zone wedges GZW1-5. Red dots show locations of sound velocity profiles (CTD, XCTD, XBT, and XSV, see text) used to calibrate the multibeam system. The yellow line show chirp profile PI-125A and the orange line seismic profile NBP9902, both shown in Fig. 6. PIT ¼ Pine Island Trough; PITE ¼ Pine Island Trough East; PITW ¼ Pine Island Trough West. The black box outlines the detailed map of acquired multibeam data shown in Fig. 2a. (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.) (PIB) in the glacial trough offshore of western PIB and Getz and Dotson ice shelves (Walker et al., 2007; Wåhlin et al., 2010). The current behavior of the PIG highlights the need to better understand its grounding line retreat history; how fast has it retreated since the Last Glacial Maximum (LGM) and what are the main drivers of rapid retreat? There is now compelling geomorphic and chronostratigraphic evidence for the WAIS having extended across the outer continental shelf during the LGM (Anderson et al., 2002; Lowe and Anderson, 2002; Ó Cofaigh et al., 2002; Evans et al., 2006; Hillenbrand et al., 2009; Mosola and Anderson, 2006; Smith et al., 2011). In the Amundsen Sea, the present shelf bathymetry is dominated by crossshelf glacial troughs (Evans et al., 2006; Nitsche et al., 2007) marking the paths of former ice streams that advanced across the shelf during past glacial maxima (Lowe and Anderson, 2002; Graham et al., 2010b). One of these troughs is the Pine IslandeThwaites Trough (PIT) (Fig. 1). Here, the Pine Island ice stream extended to the continental shelf edge and delivered water-saturated sediments generating turbidity currents that formed gully-channel systems on the uppermost continental slope (Dowdeswell et al., 2006a). This implies that the Pine Island Glacier extended more than 500 km from its present position during LGM. The central part of PIT has been sparsely mapped due to persistent sea-ice cover (Nitsche et al., 2007), although sets of grounding zone wedges were identified previously, marking still-stands in the ice retreat from the LGM position (Lowe and Anderson, 2002; Graham et al., 2010a) (Fig. 1). Prior to now, few constraints on the dynamics of the retreat history existed. Based on 14C-dates of calcareous foraminifera retrieved from sediment cores, Lowe and Anderson (2002) suggested that the ice retreated from the outer shelf before 16 ka (14C yr), but halted at about the latitude of Burke Island where it remained long enough to produce a massive grounding zone wedge (Fig. 1). The Oden-Southern-Ocean 2009/2010 (OSO0910) expedition with Swedish icebreaker IB Oden carried out swath-bathymetric mapping, chirp sonar profiling, oceanographic station work and coring in Pine Island Trough. Since the trough area was virtually ice free in the austral summer of 2010, the swath-bathymetric mapping and chirp sonar profiling were conducted as a systematic survey covering 4140 km2 of the mid-shelf section of the PIT (Fig. 1). Here we present the detailed bathymetric data and chirp sonar profiles along with a previously acquired seismic reflection profile that provide new insights into the glacial dynamics during the deglaciation of the continental shelf. The landforms imaged in PIT indicate stepwise retreat of Pine Island Ice Stream from the midshelf that was punctuated by periods of grounding line stability and minor advance followed by episodes of rapid retreat. This is similar to episodic retreat of ice streams in a number of Antarctic and Arctic cross-shelf troughs (Dowdeswell et al., 2008). Complementary to the geomorphic results presented here, Kirshner et al. (submitted for publication) provide new age constraints and sedimentary evidence for the changing glacial setting during grounding line retreat. These independent data sets support episodic post-LGM grounding line retreat within central PIT between 16.4 and 12.3 k cal yr BP. Please cite this article in press as: Jakobsson, M., et al., Ice sheet retreat dynamics inferred from glacial morphology of the central Pine Island Bay Trough, West Antarctica, Quaternary Science Reviews (2012), doi:10.1016/j.quascirev.2011.12.017 M. Jakobsson et al. / Quaternary Science Reviews xxx (2012) 1e10 2. Materials and methods 2.1. Multibeam swath-bathymetric mapping The swath-bathymetric data presented in this work were collected using a hull-mounted Kongsberg 12 kHz EM122 1 1 multibeam echo sounder. This multibeam system includes a Seatex Seapath 200 unit for integration of GPS navigation, heading and attitude information. The motion sensor consists of Seatex’s MRU5. EM122 is capable of producing 432 beams per ping and apply multipinging in the intermediate water depths. This implies that the transmit fan is duplicated after a small offset in along track tilt resulting in a doubling of ping density. Due to the ice protection of the EM122 transceivers, the useable across-track angular coverage is reduced to less than 65 65 resulting in a swath-width of typically threeefour times the water depth. The favorable ice conditions during our survey allowed swath mapping of the central part of PIT along a regular survey track (Fig. 1). This has not previously been possible due to severe sea-ice conditions in this part of the bay (e.g. Evans et al., 2006; Graham et al., 2010a). Line density was maintained with a swath overlap between about 30 and 100%. Sound velocity control was achieved through regular CTD (Conductivity, Temperature, Depth) stations supplemented with XCTD (Expendable CTD), XBT (Expendable Bathy Thermograph) and XSV (Expendable Sound Velocimeter) casts. During the 12 day-long survey in PIT, 46 stations with sound velocity control were acquired. All data were processed using the software Fledermaus and gridded to a horizontal resolution of 20 20 m. Seafloor morphology was interpreted in the 3D environment of Fledermaus and maps were subsequently produced in the GIS software Geomedia Professional. 2.2. Chirp sonar sub-bottom profiling The IB Oden sub-bottom profiler is a SBP 120 3 3 chirp sonar integrated with the multibeam through the use of the same receiving array and the capability of using the multibeam center beam depth to automatically adjust the acquisition window. During the OSO0910 survey the chirp sonar was operated continuously using a 2.5e7 kHz and 35 ms long pulse. The ping interval was set to a fixed rate of 500 ms in order to facilitate post processing of acoustic water column data acquired with the EM122 for physical oceanographic studies. Further information about the chirp sonar settings is found in the OSO0010 cruise report (Anderson and Jakobsson, 2010). 2.3. Seismic reflection profiling The seismic reflection profile used in this study was collected during a 1999 austral summer cruise (NBP9902) of the RV Nathaniel B. Palmer. Parts of the profile were published by Lowe and Anderson (2002, 2003). A 210 in3 GI air gun was used as a seismic source and a single-channel streamer was used as a receiver. The data were recorded with an Elics seismic acquisition system and band-pass filtered. The vertical resolution of the seismic profile is on the order of 10 m. Henceforth, the full seismic profile is referred to as Line NBP9902 (Fig. 1). 3 mode of formation and significance to ice dynamics during recession. Fig. 3 shows all identified landforms and their relationship to each other. These relationships will be further discussed after the individual landforms are described. 3.1. Grounding zone wedges Sedimentary wedges, or grounding zone wedges (GZW), that occur within troughs on the Antarctic continental shelf are interpreted to represent grounding line deposits formed during pauses in ice sheet’s retreat from the continental shelf (Bart and Anderson, 1996; Anderson, 1999; Mosola and Anderson, 2006; Dowdeswell et al., 2008). Their occurrence within troughs that typically occur seaward of modern ice streams and association with mega-scale glacial lineations and “soft” water saturated till indicates that they are the products of deposition by fast-flowing ice streams (Alley et al., 1989). Larger wedges are commonly 50e100 m thick and considered to be formed by a broad range of processes, including sub-glacial sedimentation and deposition via sediment gravity flows debouching from the grounding line (Anderson, 1999; Christoffersen et al., 2010). The occurrence of two GZWs within the central PIT was originally demonstrated by using swath bathymetry data and seismic reflection profiles (Lowe and Anderson, 2003). Subsequently, Graham et al. (2010a) constrained the location of the front of the larger GZW using data from two nearby multibeam swaths and referred to it as GZW5 (Fig. 1). In total, they identified five GZWs (GZW5-1; Fig. 1). GZW1 and 2 are located on the outer shelf and GZW3-5 within the central trough (Fig. 1). We will use the same nomenclature when further describing the GZWs that fall within our survey area. The OSO0910 survey includes the full extent of GZW3, 4 and 5 (Figs. 2e4). The data show that the seafloor morphology at GZW3 is actually comprised of a series of back-stepping transverse ridges that deepen in a landward direction (Fig. 5). These ridges were in part mapped by Graham et al. (2010a,b). However, it is difficult to exclude the possibility that these transverse ridges overly a small and morphologically subdued GZW. Together, the transverse ridges and grounding zone wedges record the changing elevation of the grounding line in a landward direction. The transverse ridges occur in water depths ranging from 720 to 740 m. GZW4 is located in 705e726 m and GZW5 in 726 to 670 m water depths. We further note that GZW5 appears to partly overly GZW4, the two wedges together spanning about 38 km along the axis of the trough (Fig. 4). GZW4 has a more subdued bathymetric expression relative to GZW5, with a maximum thickness of approximately 20 m (Fig. 4). GZW5 extends approximately 22 km along the trough axis, 12 km across the trough and has a maximum thickness of about 56 m when assuming its base at a water depth of 726 m (Fig. 4). However, the wedge width across the PIT is difficult to pinpoint as its outer edges, along the shallow sides of the trough, have been eroded by icebergs. The inferred width uses the outer trough walls as limits since it is evident that the wedge extends over the entire trough. The sediment volume of GZW5 within these horizontal constraints, calculated using the multibeam data, is approximately 12 km3. 3.2. Streamlined sub-glacial landforms 3. Results Several types of sedimentary landforms were produced when the ice sheet was grounded in PIT during the LGM and the subsequent retreat. Because sedimentation rates have been very low since icesheet retreat, these landforms are expressed in detail in the seafloor morphology (Fig. 2). The relatively high resolution of our multibeam swath bathymetry and pristine condition of these features allows us to describe them in detail and to interpret their The surveyed part of PIT contains streamlined ridge-groove features that conform to the morphological description of megascale glacial lineations (MSGL) (Clark, 1993; Clark et al., 2003; Wellner et al., 2006). Distinct sets of MSGLs occur north and south of GZWs 4 and 5 (Fig. 2cee). The landward set trends about 350 , parallel to the trough axis, while the seaward set trends 340e345 and oblique to the trough axis (Fig. 3). Cross-cutting relations show that the seaward set is older than the landward set. Individual MSGLs Please cite this article in press as: Jakobsson, M., et al., Ice sheet retreat dynamics inferred from glacial morphology of the central Pine Island Bay Trough, West Antarctica, Quaternary Science Reviews (2012), doi:10.1016/j.quascirev.2011.12.017 4 M. Jakobsson et al. / Quaternary Science Reviews xxx (2012) 1e10 Fig. 2. Main geomorphic features of the PIB central trough. a) Iceberg plow ridges; b) Corrugation ridges that overprint large linear to curvilinear furrows; c) Isolated corrugation ridges and MSGLs superimposed on GZW4; d) GZW5 overprinted by iceberg plowmarks; e) MSGLs and corrugation ridges in the southern part of central PIT landward of GZW5. can be followed for more than 15 km in both areas and range from >500 to 150 m in ridge-to-ridge width and have ridge heights of about 2 to >6 m. North of GZW5, the MSGLs are overprinted by iceberg furrows at a water depth above approximately 670 m, while south of GZW5 iceberg furrows generally overprint MSGLs to a depth of about 750 m. 3.3. Large curvilinear e linear furrows North of GZW4, linear to curvilinear lineations occur in water depths deeper (690 and 720 m) than the MSGLs and are aligned parallel to the axis of the trough (Fig. 2b). Some furrows have more pronounced relief with >20 m ridge height. These lineations were Please cite this article in press as: Jakobsson, M., et al., Ice sheet retreat dynamics inferred from glacial morphology of the central Pine Island Bay Trough, West Antarctica, Quaternary Science Reviews (2012), doi:10.1016/j.quascirev.2011.12.017 M. Jakobsson et al. / Quaternary Science Reviews xxx (2012) 1e10 5 GZW5 and shallower than w750 m south of this wedge (Fig. 2). These features are generally not traceable for more than 2 km before they change directions and their relief is generally <10 m. The seabed morphology of GZW5 is also dominated by iceberg furrows and plowmarks (Fig. 4). Deep-keeled Antarctic icebergs are typically calved from fast-flowing ice streams and outlet glaciers that have only limited floating margins, whereas more extensive ice shelves often produce bergs with keels no deeper than about 300 m (Dowdeswell and Bamber, 2007). 3.5. Ridges One of the most conspicuous geomorphic features in the area consists of extremely regular and small sedimentary ridges previously described by Jakobsson et al. (2011) and named corrugation ridges. These are 1e2 m from trough to crest and separated by about 60e200 m (Fig. 2b). Within the surveyed area, corrugation ridges are mostly constrained to the northern part of the trough (Figs. 2 and 3). However, some isolated sets of corrugation ridges occur also within the MSGLs (Fig. 3). The glacial dynamic connection between the large furrows, the corrugation ridges, and the iceberg plow ridges was interpreted by Jakobsson et al. (2011) and is further discussed below. 3.6. Transverse ridges (moraines) The surface of GZW3 is dominated by a series of ridges with their spatial extensions transverse to the axis of PIT (Fig. 5). They range from a few meters to 6 m high and occur between a few hundred meters to >1000 m apart. Some of these ridges bend and continue along the PIT at the eastern most part of the trough, near the outer wall (Fig. 5). 3.7. Seismic stratigraphy previously described by Jakobsson et al. (2011) and interpreted to be formed from an armada of grounded large icebergs discharged during an ice shelf break-up and likely associated grounding line retreat of the Pine Island Ice Stream. These features predate the MSGLs north of the GZW4 as indicated from cross-cutting relationships. Seismic profile NBP9902 extends along the deepest central part of PIT in the southern section of the surveyed area (Figs. 1 and 6). From the southern part of GZW5 and northward, the seismic line continues slightly westward of the deepest part of the trough. The seismic stratigraphy of GZW5 is thus imaged along its extension in a slightly oblique angle. Fig. 6 displays the section between shot point 8100 and 1600 along with chirp sonar profile PI-125A, which is the profile most closely following NBP9902. The rest of profile NBP9902 is shown in Supplementary data. The north-dipping strata underlying the trough have been eroded by the ice sheet, resulting in a prominent unconformity that is an amalgamation of several erosion surfaces visible only on the outer shelf (Lowe and Anderson, 2002). The older strata below the unconformity dip seaward at progressively lower angles north of GZW5 and near the shelf break the strata are nearly horizontal (Supplementary data). Unfortunately, a lack of crossing lines makes it difficult to fully access the inclination of the strata and the degree to which the trough has shifted in it position on the shelf with time. GZW5 is not characterized by a specific visible acoustic layering, so it is not possible to determine if there are several wedges stacked upon one another. Landwards of GZW5 there is a very thin conformable unit visible above the unconformity. The chirp profile provides limited additional information due to poor acoustic penetration. Nevertheless, a thin (<3 ms equivalent to w2 m) conformable sediment drape can be distinguished in many places at closer inspection (Fig. 6d) and sediment cores confirm the existence of these surface sediments. 3.4. Iceberg plowmarks 4. Discussion and interpretation Randomly oriented iceberg plowmarks dominate the seafloor morphology in water depths shallower than w670 m north of The landforms imaged by our survey provide a more detailed view on the WAIS retreat from the central PIT than previously revealed by Fig. 3. Interpretation of the glaciogenic geomorphic features of the central PIT and their relationship to each another. Note that the different sets of MSGLs have different orientations. Please cite this article in press as: Jakobsson, M., et al., Ice sheet retreat dynamics inferred from glacial morphology of the central Pine Island Bay Trough, West Antarctica, Quaternary Science Reviews (2012), doi:10.1016/j.quascirev.2011.12.017 6 M. Jakobsson et al. / Quaternary Science Reviews xxx (2012) 1e10 Fig. 4. 3D-view of GZWs 4 and 5. The bathymetric profile between A and B shows the wedge heights from an estimated base depth of 726 m. the more sparse data from the region (Lowe and Anderson, 2002; Graham et al., 2010a). Previous mapping of the outer shelf has suggested that the major PIT bifurcates into two paths, both reaching the shelf edge (Evans et al., 2006; Dowdeswell et al., 2006a). Following previous authors the two outer troughs of these ice stream pathways are henceforth referred to as PITW and PITE (Fig.1). Considering these troughs together with previously mapped landforms, three conceptual models for the Pine IslandeThwaites paleo-ice stream flow during the LGM were presented by Graham et al. (2010a,b). Their preferred model suggests that the different flow directions recorded in the two troughs were formed at different times and that the ice stream of PITW predates that of PITE (Fig. 7). Our new data do not allow us to constrain the age of the outer shelf troughs. It is also worth noting that the sole evidence for ice sheet grounding at the shelf break is the occurrence of gullies on the upper slope, which are argued to be the product of sediment-laden water flowing from beneath the margin of the ice sheet (Anderson, 1999; Dowdeswell et al., 2006a). However, our data do reveal a complex history of grounding line retreat from the central part of the trough (Fig. 7). The mapped landforms indicate that this retreat was marked by one episode of ice shelf break-up (Jakobsson et al., 2011) while the sediment stratigraphy suggests an additional episode of ice shelf breakup (Kirshner et al., submitted for publication). Furthermore, prolonged pauses in grounding line retreat occurred that were associated with build-up of grounding line wedges, possibly ice shelf regrowth, and an episode of ice stream re-advance. The relative timing of these events is based on the depth of geomorphic features that record grounding line depths and superposition of these features. The deepest part of the main trough north of GZW5-3 is oriented towards PITE and contains corrugation ridges (Figs. 2 and 3). These features are interpreted to have been formed during an ice shelf Fig. 5. 3D-view showing that GZW3 overprinted by back-stepping transverse ridges (moraines) that decrease in height in a landward direction as shown in the bathymetric profile from A to B. Please cite this article in press as: Jakobsson, M., et al., Ice sheet retreat dynamics inferred from glacial morphology of the central Pine Island Bay Trough, West Antarctica, Quaternary Science Reviews (2012), doi:10.1016/j.quascirev.2011.12.017 M. Jakobsson et al. / Quaternary Science Reviews xxx (2012) 1e10 7 Fig. 6. Seismic profile NBP9902 acquired in 1999 from RV Nathaniel B. Palmer (Lowe and Anderson, 2002) and chirp profile PI-125A acquired from IB Oden in 2010. The locations of the profiles are shown in Fig. 1. Panel A) displays the section between shot points 1600 and 8100 and panel B) includes an interpretation of this section with the most prominent reflectors marked. Reflectors indicating dipping bedrock strata are marked with red lines and the erosional unconformity is marked by a stippled purple line. The location of GZW5 is shaded in gray. Note that the seismic line extends west of the bathymetric expression of GZW4. Panel C shows chirp profile PI-125A that occurs between shot points 2000 and 6000. A zoom in of this profile is displayed in panel D. (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.) break-up and likely associated grounding line retreat (Jakobsson et al., 2011). This hypothesis implies that an ice shelf fringed the Pine Island Ice Stream while it extended across the outer-mid shelf. We cannot exclude that this ice shelf broke up in similar fashion as Fig. 7. Proposed LGM reconstruction for the outer Pine Island continental shelf based on results of Graham et al. (2010a,b). Initially, an ice stream reached the shelf edge at PITW. This ice stream was abandoned and a new one established to the east (PITE). the Larsen A and B ice shelves in 1995 and 2002 (MacAyeal et al., 2003), although in the case of the PIB break-up, the ice shelf must have disintegrated all the way back to the grounding line and even a bit beyond as the oriented features occur in water depths down to 650 m. The mechanism for the corrugation ridge formation involves sediment being squeezed into ridges by the trailing edge of a portion of a broken-up ice shelf rising and settling to the seafloor under tidal influence as it drifts seaward (Jakobsson et al., 2011). Eventually, at a pace of 60e200 m/day, the armada of large icebergs reached the outer more shallow part of the central PIT where the icebergs grounded to form iceberg plow ridges (Fig. 2a). It should be noted that the grounding line retreat due to the ice shelf break-up was asymmetric and appears to have occurred only in the western part of the trough. The series of landward stepping transverse ridges that occur between 720 and 740 m water depth may have been formed as the ice-stream retreated in this eastern part of the trough following the break-up of the ice shelf (Fig. 5). These small transverse ridges may be similar to De Geer moraines (Lindén and Möller, 2005). Furthermore, similar transverse ridges have elsewhere been demonstrated to be formed annually in ice contact settings (Ottesen and Dowdeswell, 2006). If the transverse ridges mapped in this work are annual formations, the average retreat rate ranged from hundreds of meters to more than 1000 m per year. Eventually this retreat stopped and a prolonged pause in grounding line retreat resulted in the formation of GZW4 when the grounding line was located at water depths Please cite this article in press as: Jakobsson, M., et al., Ice sheet retreat dynamics inferred from glacial morphology of the central Pine Island Bay Trough, West Antarctica, Quaternary Science Reviews (2012), doi:10.1016/j.quascirev.2011.12.017 8 M. Jakobsson et al. / Quaternary Science Reviews xxx (2012) 1e10 between 726 and 706 m (Fig. 5). Evident from the MSGL north set is that this still-stand was followed by an ice stream re-advance (Figs. 3 and 8b). As the grounding line was elevated to the level of GZW4 the flow was no longer controlled by the outer trough and the ice stream was free to follow a more westward path. Notably, the Graham et al. (2010a,b) model would suggest that the ice stream would actually have followed a more easterly path as it retreated, which is not consistent with our data (Fig. 7). However, since there is no coherent data coverage and due to increasing overprinting of iceberg scours west of our survey area it is not possible to determine the extent of the re-advance. The ice stream may only have advanced over GZW4 and a little distance beyond, at least as far as we mapped the MSGL north set. From its GZW 4 location, the grounding line stepped landward to the location of GZW5, where the depth of the grounding line was elevated by deposition of the wedge to between 726 and 670 m (Fig. 8c). The precise timing of the formation of a GZW is difficult to determine as it requires that sediment units constraining both the beginning and end of wedge formation are identified and dated. Our radiocarbon data are not sufficient to constrain the exact timing of GZW formation (Kirshner et al., submitted for publication). However, the GZW50 s morphology is mapped in detail allowing us to calculate that the wedge’s sediment volume per meter of ice-stream width is approximately 1 106 m3, i.e. the same volume as previously concluded per meter width of GZW5 by Graham et al. (2010a,b). The volume per meter ice stream-width, together with estimations of sediment flux per year of the Pine Island and Thwaites Ice Streams, allows us to evaluate the time required to construct GZW5 and, thus, the time of one of Pine Island Ice Stream’s still-stands during its general retreat from the central trough. If we assume that the transverse moraines located on GZW3 (Fig. 5) are annual retreat features, it is possible to estimate the sediment flux of the Pine Island ice stream when it was active at this location of PIT. The three larger ridges R1-R3 range in volume from w1260e1620 m3 m1 while the three smaller ones R4eR6 range from w200e300 m3 m1, calculated for a 1 m wide section along profile AeB shown in Fig. 5. Thus, these ridges suggest yearly sediment fluxes of between 200 and 1620 m3 a1 m1 (meter ice stream width). The recent study of the Kamb Ice Stream in Ross embayment suggests that sediment wedges near grounding lines in that area were formed mostly by melt out of basal debris (Christoffersen et al., 2010). Assuming different scenarios for Kamb Ice Stream’s active and stagnant phases over time and ice flow speeds between 500 m a1 (active) and 0 m a1 (stagnant), Christoffersen et al. (2010) calculated average sediment fluxes between 525 and 875 m3 a1 m1. This is within range of the sediment fluxes of about w100e800 m3 a1 m1 estimated for the Marguerite Bay paleo-ice stream, which flowed across the western continental shelf of the Antarctic Peninsula (Dowdeswell et al., 2004). Considering that the mapped landforms in PIT all seem to suggest a highly dynamic ice stream behavior and that the geology of the outer shelf is characterized by easily eroded sedimentary strata (Anderson, 1999), low sediment fluxes near 100 m3 a1 m1 seems unlikely. If we consider flux rates between 500 and 1650 m3 a1 m1, it would have taken between 600 and 2000 years to form GZW5 (Fig. 9). The Pine Island Ice Stream retreated from the location of GZW5 before 12.3 cal ka BP, a timing inferred from dated foraminifera in a sediment core retrieved just landward of the wedge (Fig. 8d; Kirshner et al., submitted for publication). A major question is why the Pine Island ice stream changed path during the re-advance that followed the ice shelf break-up, grounding line retreat and subsequent still-stand and formation of Fig. 8. Ice retreat from the Pine Island outer and central continental shelf based on previous work and the interpreted glaciogenic features mapped in this study. Map a) shows the oldest reconstruction while d is the youngest when the ice finally withdraw from the area at about 12.3 cal ka BP. Please cite this article in press as: Jakobsson, M., et al., Ice sheet retreat dynamics inferred from glacial morphology of the central Pine Island Bay Trough, West Antarctica, Quaternary Science Reviews (2012), doi:10.1016/j.quascirev.2011.12.017 M. Jakobsson et al. / Quaternary Science Reviews xxx (2012) 1e10 9 Acknowledgment The expedition OSO0910 was carried out as collaboration between Swedish Polar Research Secretariat, the Swedish Research Council and the US National Science Foundation (NSF). We thank the Captain and Crew of the IB Oden. Financial support was received from the Swedish Research Council (VR), the Swedish Royal Academy of Sciences through a grant financed by the Knut and Alice Wallenberg Foundation and the Bert Bolin Centre for Climate Research at Stockholm University. F. Nitsche was supported by NSF grant ANT-0838735 and Anderson by NSF/ARRA grant ANT-0837925. We thank one anonymous reviewer and Alastair Graham for insightful comments that improved the manuscript. Appendix. Supplementary data Supplementary data associated with this article can be found, in the online version, at doi:10.1016/j.quascirev.2011.12.017. References Fig. 9. Calculated formation time of GZW5 based on its volume and assuming sediment flux rates from the Pine Island paleo-ice stream between 500 and 1650 m3 a1 m1. The range of flux rates estimated by Dowdeswell et al. (2004) and Christoffersen et al., 2010 are inferred as a reference (denoted D. and Ch. respectively). GZW4? Even if the re-advance was minor, it still changed the flow path as evident from the MSGL north set (Fig. 3). It has been suggested that trough morphology and underlying geology influence ice stream flow (Dowdeswell et al., 2006b). The one available seismic profile (NBP9902) does not provide a clear answer to whether the geology in PIB could have influenced ice stream behavior (Fig. 6). However, the locations of GZWs 5e3, and thus the location for grounding line stability, do coincide with a high in the dipping strata that may indicate more resistant layers, and thus suggests that this may constitute an area of greater bed resistance to the ice stream flow. In addition, both GZW4 and GZW5 occur in locations of ‘bottle-necking’ in the trough. This, plus the fact that lineation direction occurs at a location where the trough orientation changes, highlights the importance of trough bathymetry on ice stream flow direction. 5. Conclusions In conclusion, the glacial retreat from the mid-shelf section of PIT was highly dynamic, involving both ice shelf break-up and associated grounding line retreat. The retreat appears to have been episodic in accordance with typical retreat behavior previously suggested for the last deglaciation of polar continental shelves (Dowdeswell et al., 2008). The ice retreated from the mid-shelf section of PIT first at about 12.3 cal ka BP to a location southward of our surveyed area. This landward grounding line location was likely coincident with an underlying bedrock change from seaward dipping sedimentary strata to crystalline bedrock, which has been shown to be the location of a notable change in glacial landforms (Lowe and Anderson, 2002, 2003). The depths of grounding line wedges decrease in a landward direction, from 740 to 670 m, and record elevation of the grounding line as it stepped landward. In all, the grounding line elevation varied by only w80 m over a distance of just over 100 km, likely implying a low ice sheet profile. The calculated volume of GZW5 suggests that the grounding line in the mid-shelf section of PIT was stabile, before it finally retreated from the area, during a time period between about 600 and 2000 years if sediment flux rates between 500 and 1650 m3 a1 m1 are considered. Alley, R.B., Blankenship, D.D., Rooney, S.T., Bentley, C.R., 1989. Sedimentation beneath ice shelves e the view from ice stream B. Marine Geology 85, 101e120. Anderson, J., Jakobsson, M., OSO0910 Shipboard Scientific Party, 2010. Oden Southern Ocean 0910: Cruise Report, Meddelanden Från Stockholms Universitets Institution För Gelogiska Vetenskaper. Department of Geological Sciences, Stockholm, 134 pp. Anderson, J.B., Shipp, S.S., Lowe, A.L., Wellner, J.S., Mosola, A.B., 2002. The Antarctic ice sheet during the Last Glacial Maximum and its subsequent retreat history: a review. Quaternary Science Reviews 21, 49e70. Anderson, J.B., 1999. Antarctic Marine Geology. Cambridge University Press, Cambridge. Bamber, J.L., Riva, R.E.M., Vermeersen, B.L.A., LeBrocq, A.M., 2009. Reassessment of the Potential sea-level rise from a collapse of the West Antarctic ice sheet. Science 324, 901e903. Bart, P.J., Anderson, J.B., 1996. Seismic expression of depositional sequences associated with expansion and contraction of ice sheets on the northwestern Antarctic Peninsula continental shelf. In: De Batist, M., Jacobs, P. (Eds.), Geology of Siliclastic Shelf Seas. Geological Society, New York, pp. 171e186. Christoffersen, P., Tulaczyk, S., Behar, A., 2010. Basal ice sequences in Antarctic ice stream: exposure of past hydrologic conditions and a principal mode of sediment transfer. Journal of Geophysical Research 115, F03034. Clark, C.D., Tulaczyk, M., Stokes, C.R., Canals, M., 2003. A groove-ploughing theory for the production of mega-scale glacial lineations, and implications for icestream mechanics. Journal of Glaciology 49, 240e256. Clark, C.D., 1993. Mega-scale glacial lineations and cross-cutting ice-flow landforms. Earth Surface Processes and Landforms 18, 1e29. Dowdeswell, J.A., Bamber, J.L., 2007. Keel depths of modern Antarctic icebergs and implications for sea-floor scouring in the geological record. Marine Geology 243, 120e131. Dowdeswell, J.A., Ó Cofaigh, C., Pudsey, C.J., 2004. Thickness and extent of the subglacial till layer beneath an Antarctic paleo-ice stream. Geology 32, 13e16. Dowdeswell, J.A., Evans, J., Ó Cofaigh, C., Anderson, J.B., 2006a. Morphology and sedimentary processes on the continental slope off Pine Island Bay, Amundsen Sea, West Antarctica. Geological Society of America Bulletin 118, 606e619. Dowdeswell, J.A., Ottesen, D., Rise, L., 2006b. Flow switching and large-scale deposition by ice streams draining former ice sheets. Geology 34, 313e316. Dowdeswell, J.A., Ottesen, D., Evans, J., Ó Cofaigh, C., Anderson, J.B., 2008. Submarine glacial landforms and rates of ice-stream collapse. Geology 36, 819e822. Evans, J., Dowdeswell, J.A., Ó Cofaigh, C., Benham, T.J., Anderson, J., 2006. Extent and dynamics of the West Antarctic ice sheet on the outer continental shelf of Pine Island bay during the last glaciation. Marine Geology 230, 53e72. Graham, A.G.C., Larter, R.D., Gohl, K., Dowdeswell, J.A., Hillenbrand, C.-D., Smith, J.A., Evans, J., Kuhn, G., Deen, T., 2010a. Flow and retreat of the Late Quaternary Pine IslandeThwaites palaeo-ice stream, West Antarctica. Journal of Geophysical Research 115, F03025. Graham, A.G.C., Nitsche, F.O., Larter, R.D., 2010b. An improved bathymetry compilation for the Bellingshausen Sea, Antarctica, to inform ice-sheet and ocean models. The Cryosphere Discuss 4, 2079e2101. Hillenbrand, C.D., Ehrmann, W., Larter, R.D., Benetti, S., Dowdeswell, J.A., Ó Cofaigh, C., Graham, A.G.C., Grobe, H., 2009. Clay mineral provenance of sediments in the southern Bellingshausen Sea reveals drainage changes of the West Antarctic ice sheet during the Late Quaternary. Marine Geology 265, 1e18. Holt, J.W., Blankenship, D.D., Morse, D.L., Young, D.A., Peters, M.E., Kempf, S.D., Richter, T.G., Vaughan, D.G., Corr, H.F.J., 2006. New boundary conditions for the Please cite this article in press as: Jakobsson, M., et al., Ice sheet retreat dynamics inferred from glacial morphology of the central Pine Island Bay Trough, West Antarctica, Quaternary Science Reviews (2012), doi:10.1016/j.quascirev.2011.12.017 10 M. Jakobsson et al. / Quaternary Science Reviews xxx (2012) 1e10 West Antarctic ice sheet: subglacial topography of the Thwaites and Smith glacier catchments. Geophysical Research Letters 33, L09502. Hughes, T.J., 1981. The weak underbelly of the West Antarctic ice-Sheet. Journal of Glaciology 27, 518e525. Jacobs, S.S., Jenkins, A., Giulivi, C.F., Dutrieux, P., 2011. Stronger ocean circulation and increased melting under Pine Island Glacier ice shelf. Nature Geoscience 4, 519e523. Jakobsson, M., Anderson, J.B., Nitsche, F.O., Dowdeswell, J.A., Gyllencreutz, R., Kirchner, N., O’Regan, M.A., Alley, R.B., Anandakrishnan, S., Mohammad, R., Eriksson, B., Fernandez, R., Kirshner, A., Minzoni, R., Stolldorf, T., Majewski, W., 2011. Geological record of ice shelf breakup and grounding line retreat, Pine Island bay, West Antarctica. Geology 39, 691e694. Jenkins, A., Dutrieux, P., Jacobs, S.S., McPhail, S.D., Perrett, J.R., Webb, A.T., White, D., 2010. Observations beneath Pine Island glacier in West Antarctica and implications for its retreat. Nature Geoscience 3, 468e472. Joughin, I., Rignot, E., Rosanova, C.E., Lucchitta, B.K., Bohlander, J., 2003. Timing of recent accelerations of Pine Island glacier, Antarctica. Geophysical Research Letters 30, 1706. Kirshner, A.E., Anderson, J.B., Jakobsson, M., O’Regan, M., Majewski, W., Nitsche, F.O., submitted for publication. Post-LGM deglaciation in Pine Island Bay, West Antarctica. Quaternary Science Reviews. Lindén, M., Möller, P., 2005. Marginal formation of De Geer moraines and their implications to the dynamics of grounding-line recession. Journal of Quaternary Science 20, 113e133. Lowe, A.L., Anderson, J.B., 2002. Reconstruction of the West Antarctic ice sheet in Pine Island Bay during the Last Glacial Maximum and its subsequent retreat history. Quaternary Science Reviews 21, 1879e1897. Lowe, A.L., Anderson, J.B., 2003. Evidence for abundant subglacial meltwater beneath the paleo-ice sheet in Pine Island Bay, Antarctica. Journal of Glaciology 49, 125e138. Lythe, M.B., Vaughan, D.G., Consortium, B., 2001. BEDMAP: a new ice thickness and subglacial topographic model of Antarctica. Journal of Geophysical Research 106, 11335e11352. MacAyeal, D.R., Scambos, T.A., Hulbe, C.L., Fahnestock, M.A., 2003. Catastrophic iceshelf break-up by an ice-shelf-fragment-capsize mechanism. Journal of Glaciology 49, 22e36. Mosola, A.B., Anderson, J.B., 2006. Expansion and rapid retreat of the West Antarctic Ice Sheet in eastern Ross Sea: possible consequence of over-extended ice streams? Quaternary Science Reviews 25, 2177e2196. Nitsche, F.O., Jacobs, S.S., Larter, R.D., Gohl, K., 2007. Bathymetry of the Amundsen Sea continental shelf: implications for geology, oceanography, and glaciology. Gechemistry Geophysics Geosystems 8, Q10009. Ó Cofaigh, C., Pudsey, C.J., Dowdeswell, J.A., Morris, P., 2002. Evolution of subglacial bedforms along a paleo-ice stream, Antarctic Peninsula continental shelf. Geophysical Research Letters 29, 1e4. Ottesen, D., Dowdeswell, J.A., 2006. Assemblages of submarine landforms produced by tidewater glaciers in Svalbard. Journal of Geophysical Research 111. Rignot, E., Bamber, J.L., van den Broeke, M.R., Davis, C., Li, Y., van de Berg, W.J., van Meijgaard, E., 2008. Recent Antarctic ice mass loss from radar interferometry and regional climate modelling. Nature Geoscience 1, 106e110. Smith, J.A., Hillenbrand, C.-D., Kuhn, G., Larter, R.D., Graham, A.G.C., Ehrmann, W., Moreton, S.G., Forwick, M., 2011. Deglacial history of the West Antarctic ice sheet in the western Amundsen Sea embayment. Quaternary Science Reviews 30, 488e505. Vaughan, D.G., Corr, H.F.J., Ferraccioli, F., Frearson, N., O’Hare, A., Mach, D., Holt, J.W., Blankenship, D.D., Morse, D.L., Young, D.A., 2006. New boundary conditions for the West Antarctic ice sheet: subglacial topography beneath Pine Island glacier. Geophysical Research Letters 33, L09501. Wåhlin, A.K., Yuan, X., Björk, G., Nohr, C., 2010. Inflow of Warm Circumpolar Deep water in the Central Amundsen Shelf. Journal of Physical Oceanography 40, 1427e1433. Walker, D.P., Brandon, M.A., Jenkins, A., Allen, J.T., Dowdeswell, J.A., Evans, J., 2007. Oceanic heat transport onto the Amundsen Sea shelf through a submarine glacial trough. Geophysical Research Letters 34. Wellner, J.S., Heroy, D.C., Anderson, J.B., 2006. The death mask of the Antarctic ice sheet: comparison of glacial geomorphic features across the continental shelf. Geomorphology 75, 157e171. Please cite this article in press as: Jakobsson, M., et al., Ice sheet retreat dynamics inferred from glacial morphology of the central Pine Island Bay Trough, West Antarctica, Quaternary Science Reviews (2012), doi:10.1016/j.quascirev.2011.12.017
© Copyright 2026 Paperzz