Hostname: page-component-74656f9884-j7cmm Total loading time: 0 Render date: 2026-11-07T12:03:33.852Z Has data issue: false hasContentIssue false

Asymmetric ice-sheet retreat pattern around northern Scotland revealed by marine geophysical surveys

Published online by Cambridge University Press:  05 August 2015

Tom Bradwell
Affiliation:
Biological and Environmental Sciences, University of Stirling, Stirling FK9 4LA, UK British Geological Survey, Murchison House, West Mains Road, Edinburgh EH9 3LA, UK
Martyn Stoker
Affiliation:
British Geological Survey, Murchison House, West Mains Road, Edinburgh EH9 3LA, UK
Rights & Permissions [Opens in a new window]

Abstract

This study uses marine geophysical data, principally single-beam and high-resolution multibeam echo sounder bathymetry, combined with seismic sub-bottom profiles, and existing Quaternary geological information, to map the glacial geomorphology of a large area of seafloor (∼50,000 km2) on the continental shelf around northern Scotland, from west of Lewis to north of the Orkney Islands. Our new mapping reveals the detailed pattern of submarine glacial landforms, predominantly moraines, relating to ice sheets that covered Scotland and much of the continental shelf during the Late Weichselian glaciation and earlier in the Mid to Late Pleistocene. The reconstructed retreat pattern based on geomorphological evidence highlights the large number of different retreat stages and the asymmetric, non-uniform evolution of this ice sheet sector during Late Weichselian deglaciation. Time-equivalent ice-front reconstructions show that marine sectors of the ice sheet, such as the Minch, changed their geometry significantly, perhaps rapidly; whilst other sectors remained relatively unchanged and stable. We suggest that this behaviour, governed principally by bed topography/bathymetry and ice dynamics, led to reorganisation of the Late Weichselian ice sheet as it retreated back to two main ice centres: one in Western Scotland and the other over Orkney and Shetland. This retreat pattern suggests relatively early deglaciation of NW Lewis (ca. 25 ka BP) and the mountains of far NW Scotland – the latter possibly forming a substantial ice-free land corridor. Our reconstructions differ from most previous syntheses, but are strongly supported by the independently-mapped offshore Quaternary succession and key onshore dating constraints.

Figure 0

Figure 1 Extent of the study area in the northern UK, showing general topography and bathymetry around northern Scotland. Bathymetry coloured to emphasise depth variations on the continental shelf; water depths off the continental shelf (>200 m) are not shown. Bathymetric features in italic font; islands and headlands in roman font.

Figure 1

Figure 2 Echo sounder bathymetry data coverage on the UK continental shelf within the study area (as of December 2013). Both single-beam and multibeam data sets were used in this work. Single-beam echo sounder data are part of a global dataset managed and compiled by Olex AS (Norway) and licensed for scientific research. Multibeam echo sounder data were primarily collected by the Maritime and Coastguard Agency between 2005 and 2011 and form part of ongoing UK Hydrographic Office surveys in UK waters.

Figure 2

Figure 3 (A) Distribution of Quaternary formations on the NW UK Continental shelf, determined from seismostratigraphy and seabed cores (modified from Stoker et al.1993). Where the Quaternary succession is thin, patchy or cannot be resolved in seismic profiles, it remains undifferentiated; this results in some artificial boundaries. Inset maps show (B) available marine geophysical data and (C) seabed cores, on which these original interpretations were made. Black lines are single-channel seismic profiles, mostly collected by the British Geological Survey between 1970 and 1990. Green dots are existing BGS boreholes (variable depth penetration, typically 20–300 m); orange dots are BGS vibrocores (max. 6 m penetration).

Figure 3

Figure 4 Quaternary stratigraphy of (A) the West Shetland Shelf and (B) the Hebrides Shelf, derived from seismic profiles (from Stoker et al.1993). (C) Schematic showing stratigraphic relationship between Pleistocene units on the mid-outer shelf, highlighting their angular discordance. See Figure 1 for lines of section.

Figure 4

Figure 5 Bathymetry of study area. Image (top left) shows shelf-wide digital surface model compiled from various bathymetric datasets (single-beam and multibeam); surface model is hill-shaded with illumination from the NW. These data have been used to map the seabed glacial geomorphology around northern Scotland in detail. Bathymetric transects (1–10) show cross-profile (2-D) morphology of selected ridges on the continental shelf (named and described in this study). Long dashed lines show boundaries between three sub-areas, defined for convenience: (i) the Minch and Hebrides Shelf; (ii) Central Sector; (iii) West Shetland Shelf and Orkney–Shetland Platform.

Figure 5

Figure 6 Distribution of all ice-marginal landforms mapped in this study on the NW UK Continental Shelf. Areas of bedrock at seabed are also shown. For more detailed geomorphology of sub-areas, see Figures 6–15.

Figure 6

Figure 7 (A) Hill-shaded multibeam bathymetry data, and (B) bathymetric cross-profiles of North Minch Ridge on the mid shelf, ∼40 km NW of Cape Wrath, in present-day water depths ca. 80 m. Note the irregular pitted surface, suggestive of iceberg keel marks and scours. Contains Maritime and Coastguard Agency MBES data.

Figure 7

Figure 8 BGS sparker profile (upper panel) and seismostratigraphic interpretation (lower panel) of the Eddrachillis Ridge, eastern part of the Minch. This large 30–40 m-high ridge and the smaller ridges inboard to the SE have not been formally defined within the existing Quaternary stratigraphic framework. The collective evidence suggests they are relatively late-stage (MIS2) subaqueous moraines recording grounded ice-sheet margin oscillations offshore mainland NW Scotland. For line of section, see Figure 1.

Figure 8

Figure 9 Seabed moraines in the Summer Isles region, NW Scotland (mapped and studied previously by Stoker et al.2006, 2009; Bradwell et al.2008b): (A) hill-shaded multibeam bathymetric image showing suite of transverse ridges (moraines) between Tanera Mor and Carn nan Sgeir. [Multibeam echo sounder data acquisition by BGS]; (B) hill-shaded multibeam bathymetric image of seabed moraines ∼10 km NW of Tanera Mor. Note the small delicate de Geer moraines (m) overprinting the larger (older) set of regional moraines (RM). Location of bathymetric profile (lower panel) shown by line; (C) BGS seismic reflection (boomer) profile of subaqueous moraines between Tanera Mor and Carn nan Sgeir. Note the asymmetric cross-profile. Rockhead/diamicton contact mapped where acoustically resolvable. Line of seismic profile shown in (A). Abbreviations: SBM=seabed multiple; BT=bottom tracking pulse. Modified from Bradwell et al.2008b.

Figure 9

Figure 10 Seabed moraines in the central sector, SW of Nun Rock High: (A) hill-shaded multibeam bathymetry data showing area of dense seabed ridges interpreted as ice-sheet moraines; (B) outline geomorphological map of seabed ridges. Note how the morphology of zigzag moraine (ZZM) overprints adjacent (older) ridges, suggesting a re-advance of the ice sheet margin; (C) hill-shaded greyscale bathymetric surface model, illuminated from the NW, used to map ridge morphology and crest lines in detail; (D) bathymetric cross-profiles of closely spaced, well developed seabed moraines SW of Nun Rock; profile lines shown in (A). For bathymetric colour ramp, see Figure 7. Contains Maritime and Coastguard Agency MBES data.

Figure 10

Figure 11 Seabed moraines in the central sector, S of Nun Rock High: (A) hill-shaded multibeam bathymetric image showing suites of transverse ridges between Nun Rock and Cape Wrath; (B) outline geomorphological map of ridges (large and small) interpreted as ice sheet moraines; stipple is bedrock at seabed; (C) bathymetric transects perpendicular to ridge crests showing cross-profile (2-D) morphology of selected ridges. Depths and distances in metres. Lines of transects shown in (A); (D) BGS seismic reflection (airgun) profile across seabed moraines S of Nun Rock. Line of profile shown in (A) and (B). Note the Pleistocene sediment package thickening in the vicinity of the Nun Rock moraine; sediment cover is thin or absent in places to the east. For bathymetric colour ramp, see Figure 7. Contains Maritime and Coastguard Agency MBES data. Abbreviations: NRM=Nun Rock Moraine; RM1, RM2=regional moraines; m1, m2=smaller recessional moraines; SBM=seabed multiple.

Figure 11

Figure 12 Seabed moraines off the Sutherland coast due N of Faraidh head: (A) hill-shaded multibeam bathymetric image showing the area of complex multi-phase glacial geomorphology; (B) outline geomorphological map of transverse seabed ridges interpreted as ice-marginal landforms (moraines). Larger, older ridges (grey) are regional moraines relating to ice sheet retreat generally from west to east; smaller, younger ridges (black, FHM) relate to a later phase of advance and retreat of an ice mass sourced to the S, in NW Sutherland; (C) BGS seismic reflection (airgun) profile across seabed moraines off Faraidh Head. Line of survey shown in (A) & (B). Pleistocene sediment package is locally up to 20 m thick in ridges, although sequence is undivided. Exact stratigraphic relationship between regional moraines and Faraidh Head Moraines is uncertain, owing to low resolution of airgun profile; however, geomorphology shows the latter, more pronounced, ridges superimposed on the former, broader ridges, in places. Abbreviations: FHM=Faraidh Head Moraines; RM=Regional Moraines; SBM=Seabed multiple. For bathymetric colour ramp, see Figure 6. Contains Maritime and Coastguard Agency MBES data.

Figure 12

Figure 13 Seabed moraines off the N Caithness coast: (A, B) hill-shaded single-beam echo sounder bathymetric images (Olex dataset); arrows highlight prominent arcuate sediment ridge (RR) abutting or overprinting older more subdued forms (MR): (A) illuminated from the NE; (B) illuminated from the NW. (C) BGS seismic reflection (airgun) profile across seabed offshore Caithness. Line of survey shown in (A), perpendicular to crest of ridges. Note the Pleistocene sediment package thickening to form a distinctive double crested ridge – the Reay Moraine. Steep-sided transverse sediment ridges to the E are similar to large Holocene sand waves seen elsewhere on the Continental Shelf; although seismic reflection data suggest some may be draped over Pleistocene (glacial) sediment cores. Abbreviations: MR=Melvich Ridge; RR=Reay Ridge; SBM=Seabed multiple.

Figure 13

Figure 14 Seabed moraines on the mid-shelf NW of Orkney: (A) hill-shaded echo sounder bathymetric surface model illuminated from the NW; dashed line shows join between single-beam and multibeam data; (B) hill-shaded multibeam bathymetry data showing the prominent Westray Loop Moraine (WLM) complex in detail (area of image shown in (A)). Lower panel: bathymetric profiles perpendicular to ridge complex showing different cross-profile (2-D) morphology of moraine at various points. Depths and distances in metres; (C) regional bathymetric transect showing depth profile of mid-shelf. Note the reverse (inshore-deepening) slope and location of main moraine complexes; (D) BGS seismic reflection (sparker) profile across Westray Loop Moraine and adjacent basin. Line of seismic profile shown in (A). Note the disrupted reflectors within the moraine ridge – possible evidence of glaciotectonic deformation. Contains Maritime and Coastguard Agency MBES data.

Figure 14

Figure 15 Seabed moraines on the western flank of the Orkney–Shetland Platform: (A) hill-shaded multibeam bathymetry data showing suite of long transverse ridges W of Westray; (B) outline map of transverse seabed ridges interpreted as subaqueous recessional moraines. Red arrow indicates direction of ice-sheet retreat inferred from moraine morphology. Bedrock at seabed is stippled; (C) oblique view of seabed moraines, looking NW, showing line of bathymetric profile; generated in Fledermaus software. Lower panel: bathymetric profile (in metres) perpendicular to crest line of moraines along line X–X'; (D) morphometric analysis of moraines (from W to E; 40 cross-profiles). Grey profiles show end members: smallest (and most asymmetric) and largest (almost symmetrical) ridge. Red profile is more typical. [Note: other profiles removed to aid clarity]; (E) extract of seabed slope model (derived from bathymetric xyz data) for area shown in (A). Warm colours are steeper slopes. Note the generally higher slope angles on east-facing slopes, highlighting moraine asymmetry. Steep west-facing slopes occur locally. [Diagonal stripes are data and processing artefacts.] Contains Maritime and Coastguard Agency MBES data.

Figure 15

Figure 16 Seabed moraines on the Orkney–Shetland Platform: (A) hill-shaded multibeam bathymetry data showing suites of transverse ridges in the North Sound, offshore Papa Westray, Orkney. Onshore elevation model: hill-shaded NEXTMap Britain digital surface model; (B) summary map of transverse seabed ridges interpreted as subaqueous recessional (de Geer) moraines (black, grey lines). Red arrows indicate direction of ice-sheet retreat inferred from moraine morphology; large sandbanks are shown in blue; bedrock at seabed is stippled. Note how the moraines chart retreat in two different directions, suggesting separation of a palaeo-ice front around submerged bedrock highs (marked). Dashed line shows extent of multibeam data; (C) bathymetric profiles (in metres) perpendicular to crest line of moraines, showing typical slope asymmetry. Profile lines are shown in (A). Contains Maritime and Coastguard Agency MBES data.

Figure 16

Figure 17 Seabed moraines SE of Orkney: (A) hill-shaded multibeam bathymetry data showing suite of transverse ridges ∼40 km E of Duncansby Head, Outer Moray Firth; (B) summary map of transverse seabed ridges interpreted as subaqueous recessional (de Geer) moraines (black lines). Red arrow indicates direction of ice-sheet retreat inferred from moraine morphology. Note the absence of moraines in deeper water areas (>80 m); (C) bathymetric profile (in metres) perpendicular to crest line of moraines showing typical slope asymmetry; profile line shown in (A); (D) oblique view of seabed ∼25 km ENE of moraine suite in (A), looking N; shows moraines and marine bedforms at strongly different, almost perpendicular, orientations. Lines of bathymetric profiles are also shown. Lower panel: bathymetric profiles (in metres) perpendicular to transverse features, showing the markedly different cross-profile of de Geer moraines (steeper slopes; greater relief) vs. marine bedforms. Note different scales on x and y axes. Contains Maritime and Coastguard Agency MBES data.

Figure 17

Figure 18 Reconstructed ice-sheet margins around northern Scotland, based on new geomorphological seabed mapping. Coloured lines show key ice-margin positions (Stages 1–10) during retreat. Solid lines are based on landform evidence (moraines); dashed lines are inferred or extrapolated. Key published terrestrial cosmogenic nuclide (TCN) exposure ages (from Phillips et al.2008; Bradwell et al.2008a; Everest et al.2013) are also shown. TCN ages are means of several samples. Note the strongly asymmetric pattern of deglaciation relative to the shelf break – exemplified in Stage 7 – with ice-sheet retreat back to the coastline in NW mainland Scotland, at the same time as ice-sheet oscillations on the mid shelf W of Orkney and on the outer shelf W of Shetland. This reconstructed pattern of British–Irish Ice Sheet retreat differs from current models (cf. Bradwell et al.2008a; Clark et al.2012), and points to a large ice mass centred over Orkney–Shetland at a relatively late stage during Weichselian deglaciation. This model forms a new framework which will be tested by future work.

Figure 18

Figure 19 Distribution of seabed moraines mapped from echo sounder bathymetry data (this study) overlaid on the distribution of Quaternary Formations on the NW UK Continental Shelf, mapped from seismostratigraphy and borehole data (modified from Stoker et al.1993). Ice-marginal landforms are from Figure 5; grey polygons (large); thin black lines (small) moraines. Note the clear agreement between the two independently derived maps: the generally eastward-younging Pleistocene sequence, north of 59°N, matches well with the pattern of seabed moraines, indicating predominantly eastward or southeastward ice sheet retreat.