Only 5 southern Greenland shelf edge glaciations since the early Pliocene.

Much uncertainty exists about the history of the Greenland Ice Sheet (GIS), particularly as to the frequency of extreme shelf edge glaciations. Because the last glaciation removed most of the record of earlier GIS extent on land and shelf exploration of the older GIS glacial record requires the use of deeper marine archives. Here we present seismic evidence for the frequency of extreme shelf edge glaciations offshore southwest Greenland. Our findings reveal that since the GIS formation only 5 glaciations were characterized by an ice sheet covering the entire shelf of southern Greenland. We estimate an age of around 4.5 million years (my) for the oldest episode and found that such extreme GIS expansions may have occur here only 3 times within the past c. 1.5 my. We thus conclude that the first large shelf edge glaciation of southern Greenland did occur prior to the Pliocene warmth epoch.

Documentation of ice-rafted debris (IRD) demonstrates the existence of Northern Hemisphere ice sheets back to at least 10 my ago12, though glaciers on Greenland may have existence already 20 my earlier3. As the Greenland ice core record is limited to the last glacial cycle4, the history of the GIS before this time is only known from scattered and fragmentary geological information. Only in very few places older glacial deposits are found on land or on the inner shelf56. In contrast, the outer Greenland shelf and slope are characterized by large prograding wedges of glacigenic sediments deposited by multiple ice streams7, documenting the existence of former large-scale ice sheet expansions. The greater part of the wedges is formed by down-slope mass-flow deposits triggered every time the ice reached the (paleo-) shelf edge, the so-called glacigenic debris flows (GDF)8. The wedges are the cumulative and amalgamated result of multiple shelf glaciations. Identification and discrimination of GDF units from individual glaciations is difficult, not only on the shelf, but also in upper slope areas where extreme ice sheet conditions led to calving of giant, deep-draft icebergs scouring the seabed down to about 1000 m of water depth9. In contrast, in the adjacent deep water areas away from direct ice contact, records from the most extreme shelf edge glaciations may be preserved owing to a long run-out potential of the GDF's8. A first prerequisite for tracing here extreme shelf-edge glaciation through time is a marked difference between glacial and inter-glacial depositional environments. We have found that the basinal area off the glaciated southwest Greenland margin bordering the Davis Strait Drift Complex (DSDC10) toward west fulfils these conditions (Fig. 1a). In extreme glacial periods this basin is dominated by the GDF deposition, in contrast to DSDC current dominated deposition during interglacial periods.

Figure 1

(a) Bathymetry of the study area together with the outline of the main body of the Davis Strait Drift Complex (DSDC, orange), and location of the seismic section shown in Fig. 2. The basin in question forms the seabed depression between the shelf break and the drift crest. The insert map shows the regional setting of the study area and present-day extent of the Greenland Ice Sheet. (b) Study area with the grid of seismic lines used for this study. Also shown is the location of exploration well Qulleq-1 which provided the base for the seismo-stratigraphic framework as previously published10. The green lines mark the location of the seismic examples shown in Figs. 2 and 3. The inset map shows the pattern of the present-day main current systems, both shallow (small arrows) and deep water (thick arrow). WGC = West Greenland Current; DWBC = Deep Western Boundary Current. For further details, see Ref. 11. Bathymetry compiled by GEUS. Maps produced by GEUS.

The thermohaline circulation regime is known to fluctuate with the glacial-interglacial cycles11. The present-day oceanographic setting of the study area is characterized by significant boundary current activity both at shallow depth (north flowing West Greenland Current) and greater depth (Fig. 1b, inset map)11. This current system has been responsible for the formation of the DSDC since mid Miocene times10. An oceanographic setting similar to the present-day may be expected in interglacial times, while a current slow-down and change of the circulation pattern occurred in glacial times512.

Results

Based on a dense grid of conventional 2D seismic data (Fig. 1b; see Methods), our study is the first to report the frequency of extreme southern Greenland shelf edge glaciations. Stratigraphic information has been derived from a published seismo-stratigraphy of the DSDC10. The architectural development of the basin deposits is studied by identifying intervals representative for a prevalence of either current- or glacial-related processes. The seismic identification of the two resulting depositional types is based on reflection pattern variations (see Methods). We present here examples of multi-channel seismic profiles (Figs. 2, 3) showing the alternation of units mainly representative of a contourite facies and intercalated several major GDF units. The number of major GDF units in southern Davis Strait is constant (five), whereas at least six or seven major GDF units can be observed in sections more to the north. Common for all GDF units are the timing of the appearance of the oldest one (GDF 1), which based on the regional seismo-stratigraphic model10 is assigned to have an early Pliocene age.

Figure 2

Overview (A) and zoom (B) of seismic section crossing the basin offshore Southwest Greenland, extending from the glacigenic prograding wedge deposits on the Greenland shelf in the east to the current-dominated Davis Strait Drift Complex (DSDC) deposits in central Davis Strait in the west. The seismic section illustrates the chaotic and semi-transparent appearance of the five major glacigenic debris flow (GDF) units (blue) identified within the DSDC deposits (orange). The five GDF units are labelled 1 to 5 from bottom to top. The seismic section documents that GDF depositional processes can extend over distances of at least 60 km from the base of the slope. For location – see Fig. 1.

Figure 3

Overview (A) and zoom (B) of seismic section crossing the basin offshore Southwest Greenland, extending from the glacigenic prograding wedge deposits on the Greenland shelf in the east to the current-dominated Davis Strait Drift Complex (DSDC) deposits in central Davis Strait in the west. Location of the seismic section is shown in Fig. 1. For seismic facies interpretation, see Fig. 2. This section illustrates that a similar seismic scenario with five GDF units is still present further to the north. In some cases the GDF units suggest a 2-stage formation as indicated by an intercalated continues reflection horizon (e.g. Unit 3 in the lower panel B). This we attribute to ice margin fluctuations across the shelf edge associated with a stadial-interstadial climate shift.

Discussion

An early Pliocene age of GDF 1 is supported by records of ice-rafted materials with Greenlandic provenance documented by ODP drilling Leg 10513 and comparable to the timing of glacier formation on Iceland14, but notably postdates the first occurrence of drop stones off south-eastern Greenland around 7 my ago15. The Neogene drop stone sequence from this area has been interpreted to reflect the results of three major glaciations of southeast Greenland between 2.0 and 5.3 my ago. This is supported by the results from ODP Site 987 on the east Greenland Scoresby Sund Fan, revealing glacial debris flows back to about 5 my ago16. Based on our seismo-stratigraphy10, and supported by the age of the first major coarse-sand peak recorded at ODP Site 918 off Southeast Greenland17, we propose an age of near 4.6 my for GDF 1. This age is close to the timing of a negative excursion observed in the time-adjusted global sea level curve1819 suggesting a prominent glaciation event at Northern and Southern high latitudes. These various sources of evidence make us conclude that the first extreme shelf edge glaciation of southern Greenland can be dated between 4.4 and 4.6 my. A full and extreme southern Greenland glaciation at that time does, however, not imply that a permanent ice sheet necessarily has existed here since. Major GIS reductions can be assumed during subsequent Pliocene warming between 3.6 my and 2.6 my20 and sometime later in the Pleistocene21.

As for the youngest GDF 5 unit, we tentatively assign an age corresponding to the Saale glaciation (c. 0.16 my; Marine Isotope Stage 6).This glaciation was characterized by an extreme northern hemisphere ice sheet expansion22, and must have been responsible for the production of giant, deep-drafted icebergs originating from southern Greenland ice streams23. For information on younger, post-Saale deposits we may refer to the Methods section below. In between these oldest and youngest GDF episodes, glacial advances had resulted in a significant sediment input to the shelf and upper slope, but during this long period only three glaciations produced an extensive ice sheet responsible for the deposition of basinal GDF units. Using our seismo-stratigraphic model we can not assign with high confidence a more specific age to each of these three GDF units. Still, based on recent GIS modelling24 we suggest that GDF 2 could originate from the earliest Quaternary glaciation era (c. 2.5 my). A shelf edge glaciation during this period could explain the increased accumulation of coarse sand off southeast Greenland dated at ODP site 918 between ca. 2.5 and 2.0 my ago15. For reason of climate warmth both the Mid-Pliocene and early Pleistocene2025 may be ruled out for age consideration of one of the remaining GDF units (GDF 3 and 4). The same may apply to the ‘early Brunhes’ era, ca. 0.40–0.75 my, when inter-hemispheric climate asymmetry appears to have been larger, and climate belt and ocean fronts were displaced to the north262728. During this period ice-free interglacial conditions with more abundant vegetation may have been widespread in southern Greenland21. Based on our seismo-stratigraphic model we may infer that GDF 3 was deposited after about 1.5 my ago, but further data lack for assigning a more precise age.

The seismic records (Figs. 2, 3) suggest a close timing for GDF 4 and GDF 5, which could imply that these may represent an Elster (GDF 4) – Saale (GDF 5) complex. Such a close timing may also be inferred from the coarse sand mass accumulation record of ODP site 91815 showing continuously high accumulation rates during the latter period.

Regarding the GDF frequency, it is noteworthy that seismic records of the East Greenland margin further north display six units that likewise were identified as representing Plio-Pleistocene phases of extensive ice sheets29. This larger number (six) fits well with our records from shelf areas further north in Davis Strait. An important factor which presumably affects the large-scale GIS morphology and extent of shelf glaciation is (subsurface) advection of warm, Atlantic-derived Irminger Sea Water. This mainly influences southern Greenland shelf waters south of Denmark Strait and Davis Strait, respectively11. It has been a matter of debate how far the initiation of the North Atlantic Current system has had an impact on large-scale GIS expansion. Also mechanisms such as a decline in atmospheric CO230 and changed atmospheric circulation patterns due to tectonic uplift of the Rocky Mountains and Himalaya31 have been proposed. The latter changes may have affected the cyclone tracks and precipitation patterns over the North Atlantic. It would be beyond the scope of this study to discuss these possible trigger mechanisms of the early Pliocene GIS expansion in more detail. Within this context it is, however, worth noting that a recent study shows an early Pliocene maximum ocean temperature similar as today32.

Our findings show that although the late Neogene and Quaternary oxygen isotope record indicate altogether more than 30 glacial stages, only relatively few (five) were characterized by a GIS geometry covering the entire shelf offshore southern Greenland. Greenland shelf areas more to the north are likely to have experienced a larger number of shelf edge glaciations. We estimate an age of around 4.5 my for the oldest southern Greenland shelf edge glaciation and conclude that such extreme GIS expansions did occur here only three times within the past 1.5 my, the last event probably corresponding to the Saale glaciation (Marine Isotope Stage 6). Our findings demonstrate the large natural variability of GIS dynamic behaviour with extreme boundary conditions that apparently have changed through time. Grasping the entire range of natural variability is a precondition for improving modelling of future GIS behaviour.

Methods

A large digital database of conventional 2D seismic reflection data, with a line spacing varying between 5 and 10 km, was used for the study. The database comprises data from a series of individual surveys carried out since the 1980's. Stratigraphic information from the exploration well Qulleq-1 was used for seismic-to-well calibration.

We distinguish between a parallel high-amplitude seismic reflection pattern and a chaotic and semi-transparent reflection pattern. The parallel reflection pattern originates primarily from bottom current deposition, but can also indicate turbidite deposition. To differentiate between these two types of deposits is not always easy. However, in both cases there is seismically a marked difference with the GDF units that are characterized by a chaotic seismic reflection patterns or a semi-transparent appearance without any distinct continuous internal reflections. In contrast to melt-water generated turbidite deposits, the GDF units can be primarily related to sedimentary processes associated with shelf edge glaciation, providing thus a sedimentary record of such extreme glaciations. Further, the GDF can be distinguished from non-glacigenic slide and slump deposits by the lack of up-slope slide scars. Evidence for the latter is entirely lacking in our study area.

The seismic vertical resolution is ca. 10 m, which implies that an upper, post-Saale sediment unit less than about 10 m thick cannot be recognized on the seismic records. Sediment core data confirm the presence of such a post-Saale unit in the basin33, which for the LGM yields a well-stratified sequence clearly lacking any GDF signature. This is supported by studies on the Southwest Greenland shelf that demonstrate that during the LGM the position of the GIS was on the mid shelf6.

Author Contributions

T.N. is responsible for the seismic interpretation and A.K. for chronological aspects and climate-related issues. Both authors contributed equally to the conclusion and writing of the paper.