Relative sea‐level observations in western Scotland since the Last Glacial Maximum for testing models of glacial isostatic land movements and ice‐sheet reconstructions

Observations of relative sea‐level change and local deglaciation in western Scotland provide critical constraints for modelling glacio‐isostatic rebound in northern Britain over the last 18 000 years. The longest records come from Skye, Arisaig and Knapdale with a shorter, Holocene, record from Kintail. Biostratigraphic (diatom, pollen, dinoflagellate, foraminifera and thecamoebian), lithological and radiocarbon analyses provide age and elevation parameters for each sea‐level index point. All four sites reveal relative sea‐level change that is highly non‐monotonic in time as the local vertical component of glacio‐isostatic rebound and eustasy (or global meltwater influx) dominate at different periods. Copyright © 2006 John Wiley & Sons, Ltd.     

Influence of glacial isostatic adjustment upon current sea level variations in the Mediterranean

Present-day sea level variations in the Mediterranean depend on various factors, including recent climatic forcing, tectonic activity, anthropogenic effects, and glacio-isostatic adjustment. The latter is governed by mantle rheology and the spatio–temporal distribution of the late-Pleistocene ice sheets and it is expected to produce a long-wavelength pattern of sea level variations across the Mediterranean, mostly determined by the response of the solid earth and of the geoid to loading effects of melt water since the end of deglaciation. Modeling glacio-isostatic effects in this region is necessary for a correct interpretation of tide gauge and GPS time-series, and thereby to constrain both the present-day climate-related sea level rise and regional or local geological, tectonic and human-driven displacements. By an exhaustive exploration of the parameter space of mantle rheology and ice sheet chronologies, in this work we outline upper and lower bounds on the current rate of sea level variation associated with glacial isostatic adjustment in the Mediterranean. This may contribute to a full assessment of coastal vulnerability by sea level rise on a regional and local scale.     

Relative sea‐level changes,glacial isostatic modelling and ice‐sheet reconstructions from the British Isles since the Last Glacial Maximum

While contributing <1 m equivalent eustatic sea‐level rise the British Isles ice sheet produced glacio‐isostatic rebound in northern Britain of similar magnitude to eustatic sea‐level change, or global meltwater influx, over the last 18 000 years. The resulting spatially variable relative sea‐level changes combine with observations from far‐field locations to produce a rigorous test for quantitative models of glacial isostatic adjustment, local ice‐sheet history and global meltwater influx. After a review of the attributes of relative sea‐level observations significant for constraining large‐scale models of the isostatic adjustment process we summarise long records of relative sea‐level change from the British Isles and far‐field locations. We give an overview of different global theoretical models of the isostatic adjustment process before presenting intercomparisons of observed and predicted relative sea levels at sites in the British Isles and far‐field for a range of Earth and ice model parameters in order to demonstrate model sensitivity and the resolving power available from using evidence from the British Isles. For the first time we show a good degree of fit between relative sea‐level observations and predictions that are based upon global Earth and ice model parameters, independently derived from analysis of far‐field data, with a terrain‐corrected model of the British Isles ice sheet that includes extensive glaciation of the North Sea and western continental shelf, that does not assume isostatic equilibrium at the Last Glacial Maximum and keeps to trimline constraints of ice surface elevation. We do not attempt to identify a unique solution for the model lithosphere thickness parameter or the local‐scale detail of the ice model in order to provide a fit for all sites, but argue that the next stage should be to incorporate an ice‐sheet model that is based on quantitative, glaciological model simulations. We hope that this paper will stimulate this debate and help to integrate research in glacial geomorphology, glaciology, sea‐level change, Earth rheology and quantitative modelling. Copyright © 2006 John Wiley & Sons, Ltd.     

Lateglacial and Holocene relative sea‐level changes and first evidence for the Storegga tsunami in Sutherland,Scotland

We reconstruct one of the longest relative sea‐level (RSL) records in north‐west Europe from the north coast of mainland Scotland, using data collected from three sites in Loch Eriboll (Sutherland) that we combine with other studies from the region. Following deglaciation, RSL fell from a Lateglacial highstand of +6?8 m OD (Ordnance Datum = ca. mean sea level) at ca. 15 k cal a BP to below present, then rose to an early Holocene highstand and remained at ca. +1 m OD between ca. 7 and 3 k cal a BP, before falling to present. We find no evidence for significant differential Holocene glacio‐isostatic adjustment between sites on the north‐west (Lochinver, Loch Laxford), north (Loch Eriboll) and north‐east (Wick) coast of mainland Scotland. This suggests that the region was rapidly deglaciated and there was little difference in ice loads across the region. From one site at the head of Loch Eriboll we report the most westerly sedimentary evidence for the early Holocene Storegga tsunami on the Scottish mainland. The presence of the Storegga tsunami in Loch Eriboll is predicted by a tsunami wave model, which suggests that the tsunami impacted the entire north coast of Scotland and probably also the Atlantic coastline of north‐west Scotland.

     

An improved glacial isostatic adjustment model for the British Isles

The glacial isostatic adjustment (GIA) of the British Isles is of interest due to the constraints that can be provided on key model parameters such as the global meltwater signal, local ice sheet history and viscoelastic earth structure. A number of recent studies have modelled relative sea‐level (RSL) data from this region to constrain model parameters. As indicated in these studies, the sensitivity of these data to both local and global parameters results in a highly non‐unique problem. This study aims to address this inherent non‐uniqueness by combining a previously published British–Irish ice model that is based on the most recent geomorphological data with a new global ice sheet model that provides an accurate prediction of eustatic sea‐level change. In addition, constraints from Global Positioning System (GPS) measurements of present‐day vertical land motion are considered alongside the entirety of RSL data from both Great Britain and Ireland. A model solution is found that provides a high‐quality fit to both the RSL data and the GPS data. Within the range of earth viscosity values considered, the optimal data model fits were achieved with a relatively thin lithosphere (71 km), upper mantle viscosities in the range 4–6 × 10²⁰ Pa s and lower mantle viscosities ≥ 3 × 10²² Pa s. Copyright © 2011 John Wiley & Sons, Ltd.     

Late Quaternary sea‐level change and evolution of Belfast Lough,Northern Ireland: new offshore evidence and implications for sea‐level reconstruction

The interplay of eustatic and isostatic factors causes complex relative sea‐level (RSL) histories, particularly in paraglacial settings. In this context the past record of RSL is important in understanding ice‐sheet history, earth rheology and resulting glacio‐isostatic adjustment. Field data to develop sea‐level reconstructions are often limited to shallow depths and uncertainty exists as to the veracity of modelled sea‐level curves. We use seismic stratigraphy, 39 vibrocores and 26 radiocarbon dates to investigate the deglacial history of Belfast Lough, Northern Ireland, and reconstruct past RSL. A typical sequence of till, glacimarine and Holocene sediments is preserved. Two sea‐level lowstands (both max. ?40 m) are recorded at c. 13.5 and 11.5k cal a bp . Each is followed by a rapid transgression and subsequent periods of RSL stability. The first transgression coincides temporally with a late stage of Meltwater Pulse 1a and the RSL stability occurred between c. 13.0 and c. 12.2k cal a bp (Younger Dryas). The second still/slowstand occurred between c. 10.3 and c. 11.5k cal a bp . Our data provide constraints on the direction and timing of RSL change during deglaciation. Application of the Depth of Closure concept adds an error term to sea‐level reconstructions based on seismic stratigraphic reconstructions.     

Mid‐ to late Holocene relative sea‐level change in Poole Harbour,southern England

A foraminiferal transfer function for mean tide level (MTL) is used in combination with AMS radiocarbon dated material to construct a record of relative sea‐level (RSL) change from Poole Harbour, southern Britain. These new data, based on multiple cores from duplicate sites, indicate four phases of change during the last 5000 cal. (calendar) yr: (i) rising RSL between ca. 4700 cal. yr BP and ca. 2400 cal. yr BP; (ii) stable to falling RSL from ca. 2400 cal. yr BP until ca. 1200 cal. yr BP; (iii) a brief rise in RSL from ca. 1200 cal. yr BP to ca. 900 cal. yr BP, followed by a period of stability; (iv) a recent increase in the rate of RSL rise from ca. 400–200 cal. yr BP until the present day. In addition, they suggest that the region has experienced long‐term crustal subsidence at a rate of 0.5 mm C¹⁴ yr^?1. Although this can account for the overall rise in MTL observed during the past 2500 yr, it fails to explain the changes in the rate of rise during this period. This implies that the phases of RSL change recorded in the marshes of Poole Harbour reflect tidal range variations or ‘eustatic’ fluctuations in sea‐level. Copyright © 2001 John Wiley & Sons, Ltd.     

Postglacial relative sea‐level observations from Ireland and their role in glacial rebound modelling

The British Isles have been the focus of a number of recent modelling studies owing to the existence of a high‐quality sea‐level dataset for this region and the suitability of these data for constraining shallow earth viscosity structure, local to regional ice sheet histories and the magnitude/timing of global meltwater signals. Until recently, the paucity of both glaciological and relative sea‐level (RSL) data from Ireland has meant that the majority of these glacial isostatic adjustment (GIA) modelling studies of the British Isles region have tended to concentrate on reconstructing ice cover over Britain. However, the recent development of a sea‐level database for Ireland along with emergence of new glaciological data on the spatial extent, thickness and deglacial chronology of the Irish Ice Sheet means it is now possible to revisit this region of the British Isles. Here, we employ these new data to constrain the evolution of the Irish Ice Sheet. We find that in order to reconcile differences between model predictions and RSL evidence, a thick, spatially extensive ice sheet of ～600–700 m over much of north and central Ireland is required at the LGM with very rapid deglaciation after 21 k cal. yr BP. Copyright © 2007 John Wiley & Sons, Ltd.     

Late Devensian and Holocene relative sea‐level change in North Wales,UK

Deglacial sea‐level index points defining relative sea‐level (RSL) change are critical for testing glacial isostatic adjustment (GIA) model output. Only a few observations are available from North Wales and until recently these provided a poor fit to GIA model output for the British‐Irish Ice Sheet. We present results of an integrated offshore geophysical (seismic reflection), coring (drilling rig), sedimentological, micropalaeontological (foraminifera), biostratigraphical (palynology) and geochronological (AMS ¹⁴C) investigation into a sequence of multiple peat/organic sediment horizons interbedded within a thick estuarine–marine sequence of minerogenic clay‐silts to silty sands from the NE Menai Strait, North Wales. Ten new sea‐level index points and nine new limiting dates from the Devensian Late‐glacial and early Holocene are integrated with twelve pre‐existing Holocene sea‐level index points and one limiting point from North Wales to generate a regional RSL record. This record is similar to the most recent GIA predictions for North Wales RSL change, supporting either greater ice load and later deglaciation than in the GIA predictions generated before 2004, or a modified eustatic function. There is no evidence for a mid‐Holocene highstand. Tidally corrected RSL data indicate initial breaching of the Menai Strait between 8.8 and 8.4 ka BP to form a tidal causeway, with final submergence between 5.8 and 4.6 ka BP. Final breaching converted the NE Menai Strait from a flood‐dominated estuary into a high energy ebb tidal delta with extensive tidal scouring of pre‐existing Late‐glacial and Holocene sequences. The study confirms the value of utilising offshore drilling/coring technology to recover sea‐level records which relate to intervals when rates of both eustatic and isostatic change were at their greatest, and therefore of most value for constraining GIA models. Copyright © 2011 John Wiley & Sons, Ltd.     

Late Holocene great earthquakes and relative sea‐level change at Kenai,southern Alaska

Kenai, located on the west coast of the Kenai Peninsula, Alaska, subsided during the great earthquake of AD 1964. Regional land subsidence is recorded within the estuarine stratigraphy as peat overlain by tidal silt and clay. Reconstructions using quantitative diatom transfer functions estimate co‐seismic subsidence (relative sea‐level rise) between 0.28±0.28 m and 0.70±0.28 m followed by rapid post‐seismic recovery. Stratigraphy records an earlier co‐seismic event as a second peat‐silt couplet, dated to ～1500–1400 cal. yr BP with 1.14±0.28 m subsidence. Two decimetre‐scale relative sea‐level rises are more likely the result of glacio‐isostatic responses to late Holocene and Little Ice Age glacier expansions rather than to co‐seismic subsidence during great earthquakes. Comparison with other sites around Cook Inlet, at Girdwood and Ocean View, helps in constructing regional patterns of land‐level change associated with three great earthquakes, AD 1964, ～950–850 cal. yr BP and ～1500–1400 cal. yr BP. Each earthquake has a different spatial pattern of co‐seismic subsidence which indicates that assessment of seismic hazard in southern Alaska requires an understanding of multiple great earthquakes, not only the most recent. All three earthquakes show a pre‐seismic phase of gradual land subsidence that marked the end of relative land uplift caused by inter‐seismic strain accumulation. Copyright © 2005 John Wiley & Sons, Ltd.     

Relative sea‐level change and postglacial isostatic adjustment along the coast of south Devon,United Kingdom

Previous sea‐level studies suggest that southwest Britain has the fastest subsiding coastline in the United Kingdom, but tide‐gauge data, GPS and gravity measurements and geophysical models show little evidence of anomalous subsidence in this region. In this paper we present 15 new sea‐level index points from four coastal barrier systems in south Devon. Eight are from compaction‐free basal sediments and others were corrected for autocompaction. Our data suggest that relative sea level along the south Devon coastline has risen by 21 ± 4 m during the past 9000 years. Sea‐level rise slowed during the middle and late Holocene and a rise of 8 ± 1 m has occurred since ca. 7000 cal. yr BP. Anomalous ages for many rejected points are attributed to sediment reworking during barrier transgression. The relative sea‐level history during the early and middle Holocene shows a good fit with geophysical model predictions, but the geological and modelled data diverge in the later Holocene. Unlike the geophysical models, sea‐level index points cannot differentiate between late Holocene relative sea‐level histories of south Devon and southwest Cornwall. It is suggested that this discrepancy can be resolved by obtaining additional high‐quality sea‐level index points covering the past 4000 years. Copyright © 2008 John Wiley & Sons, Ltd.     

Lateglacial to Holocene relative sea‐level changes in the Stykkishólmur area,northern Snæfellsnes,Iceland

Comparatively little research has been undertaken on relative sea‐level (RSL) change in western Iceland. This paper presents the results of diatom, tephrochronological and radiocarbon analyses on six isolation basins and two coastal lowland sediment cores from the Stykkishólmur area, northern Snæfellsnes, western Iceland. The analyses provide a reconstruction of Lateglacial to mid‐Holocene RSL changes in the region. The marine limit is measured to 65–69 m above sea level (asl), with formation being estimated at 13.5 cal ka BP. RSL fall initially occurred rapidly following marine limit formation, until ca. 12.6 cal ka BP, when the rate of RSL fall decreased. RSL fell below present in the Stykkishólmur area during the early Holocene (by ca. 10 cal ka BP). The rates of RSL change noted in the Stykkishólmur area demonstrate lesser ice thicknesses in Snæfellsnes than Vestfirðir during the Younger Dryas, when viewed in the regional context. Consequently, the data provide an insight into patterns of glacio‐isostatic adjustment surrounding Breiðafjörður, a hypothesized major ice stream at the Last Glacial Maximum.     

A tephrochronology for the Lateglacial palynological record of the Endinger Bruch (Vorpommern,north‐east Germany)

The tephrostratigraphy of lake sediments in the Endinger Bruch provides the first robust age model for the Lateglacial palynological records of Vorpommern (north‐east Germany). Cryptotephra investigations revealed six tephra layers within sediments spanning from Open vegetation phase I (～Bølling, ～15 ka) to the Early Holocene Betula/Pinus forest phase (～Pre‐boreal, ～10.5 ka). Four of these layers have been correlated with previously described tephra layers found in sites across Europe. The Laacher See Tephra (Eifel Volcanic Field) is present in very high concentrations within sediments of the Lateglacial Betula (/Pinus) forest phase (～Allerød). The Vedde Ash (Iceland) lies midway through Open vegetation phase III (～Younger Dryas). The Hässeldalen and the Askja tephras (Iceland) lie in the Early Holocene Betula/Pinus forest phase (～Preboreal). These tephra layers have independently derived age estimates, which have been imported into the Endinger Bruch record. Furthermore, the layers facilitate direct correlation of the regional vegetation record with other palaeoenvironmental archives, which contain one or more of the same tephra layers, from Greenland to Southern Europe. In doing this, localized variations are confirmed in some aspects of the pollen stratigraphy; however, transitions between the main vegetation phases appear to occur synchronously (within centennial errors) with the equivalent environmental transitions observed in sites across the European continent. Copyright © 2011 John Wiley & Sons, Ltd.     

Evaluation of a numerical model of the British–Irish ice sheet using relative sea‐level data: implications for the interpretation of trimline observations

The glacial isostatic adjustment (GIA) of the British Isles is complex due to the interplay between local and non‐local signals. A number of recent studies have modelled the GIA response of the British Isles using relative sea‐level data. This study extends these previous analyses by using output from a numerical glaciological model as input to a GIA model. This is a necessary step towards more realistic GIA models, and although there have been similar studies for the major late Pleistocene ice sheets, this is the first study to do so for the British Isles. We test three reconstructions, classed as ‘minimal’, ‘median’ and ‘maximal’ in terms of their volume at maximum extent, and find it is possible to obtain good data–model fits. The minimal reconstruction is clearly preferred by the sea‐level data. The ice reconstructions tested were not constrained by geomorphological information of past ice extent (lateral and vertical). As a consequence, the reconstructions extend further than much of this information suggests, particularly in terms of ice thickness. It is notable, however, that the reconstructions produce good fits to many sea‐level data from central, mountainous regions (e.g. Scottish highlands), which lends support to the suggestion that trimlines, often used as an constraint on the palaeo ice surface location, are in fact features formed at the transition from warm‐ to cold‐based ice and so mark a minimum constraint on the ice surface altitude. Based on data–model misfits, suggestions are made for improving the ice model reconstructions. However, in many locations, the cause of the misfit could be due to non‐local, predominantly Fennoscandian ice and so interpretation is not straightforward. As a result, we suggest that future analyses of this type consider models and observations for both Fennoscandia and the British Isles. Copyright © 2012 John Wiley & Sons, Ltd.     

A new Holocene relative sea‐level curve for Terra Nova Bay,Victoria Land,Antarctica

More than 100 radiocarbon dates of penguin guano and remains, shells and seal skin afford ages for raised beaches adjacent to Terra Nova Bay, Antarctica. These dates permit construction of a new relative sea‐level curve that bears on the timing of deglaciation. Recession of the Ross Sea ice‐sheet grounding line from Terra Nova Bay occurred no earlier than 7200 ¹⁴C yr (8000 cal. yr) BP. Retreat along the Victoria Land coast may have been rapid, possibly contributing to eustatic sea‐level rise centred at ca. 7600 cal. yr BP. The presence of a significant amount of ice remaining in the Ross Sea Embayment in Holocene time lessens the chance that Antarctica contributed significantly to meltwater pulse 1A several thousand years earlier. Copyright © 2004 John Wiley & Sons, Ltd.     

Holocene relative sea‐level changes in Harris,Outer Hebrides,Scotland, UK

Evidence for relative sea‐level changes during the middle and late Holocene is examined from two locations on the Atlantic coast of Harris, Outer Hebrides, Scotland, using morphological mapping and survey, stratigraphical, grain size and diatom analysis, and radiocarbon dating. The earliest event identified is a marine flood, which occurred after 7982–8348 cal. a (7370 ± 80 ¹⁴C a) BP, when the sea crossed a threshold lying at ?0.08 m Ordnance Datum Newlyn (OD) (?2.17 m mean high water springs (MHWS)) before withdrawing. This could have been due to a storm or to the Holocene Storegga Slide tsunami. By 6407–6122 cal. a (5500 ± 60 ¹⁴C a) BP, relative sea levels had begun to fall from a sandflat surface with an indicated MHWS level of between 0.08 and ?1.96 m (?2.01 to ?4.05 m). This fall reached between ?0.30 and ?2.35 m (?2.39 to ?4.44 m) after 5841–5050 cal. a (4760 ± 130 ¹⁴C a) BP, but was succeeded by a relative sea‐level rise which reached between 0.54 and ?1.57 m (?1.55 to ?3.66 m) by 5450–4861 cal. a (4500 ± 100 ¹⁴C a) BP. This rise continued, possibly with an interruption, until a second sandflat surface was reached between 2.34 and ?0.26 m (0.25 to ?2.35 m) between 2952–3375 cal. a (3000 ± 80 ¹⁴C a) and 1948–2325 cal. a (2130 ± 70 ¹⁴C a) BP, before present levels were reached. The regressive episode from the earliest sandflat is correlated with the abandonment of the Main Postglacial Shoreline. It is maintained that the fluctuations in relative sea level recorded can be correlated with similar events elsewhere on the periphery of the glacio‐isostatic centre and may therefore reflect secular changes in nearshore sea surface levels. Despite published evidence from trim lines of differential ice sheet loading across the area, no evidence of variations in uplift between the locations concerned could be found. Copyright © 2009 John Wiley & Sons, Ltd.     

High‐resolution reconstruction of a coastal barrier system: impact of Holocene sea‐level change

This study presents a detailed reconstruction of the sedimentary effects of Holocene sea‐level rise on a modern coastal barrier system. Increasing concern over the evolution of coastal barrier systems due to future accelerated rates of sea‐level rise calls for a better understanding of coastal barrier response to sea‐level changes. The complex evolution and sequence stratigraphic framework of the investigated coastal barrier system is reconstructed using facies analysis, high‐resolution optically stimulated luminescence and radiocarbon dating. During the formation of the coastal barrier system starting 8 to 7 ka rapid relative sea‐level rise outpaced sediment accumulation. Not before rates of relative sea‐level rise had decreased to ca 2 mm yr^?1 did sediment accumulation outpace sea‐level rise. From ca 5·5 ka, rates of regionally averaged sediment accumulation increased to 4·3 mm yr^?1 and the back‐barrier basin was filled in. This increase in sediment accumulation resulted from retreat of the barrier island and probably also due to formation of a tidal inlet close to the study area. Continued transgression and shoreface retreat created a distinct hiatus and wave ravinement surface in the seaward part of the coastal barrier system before the barrier shoreline stabilized between 5·0 ka and 4·5 ka. Back‐barrier shoreline erosion due to sediment starvation in the back‐barrier basin was pronounced from 4·5 to 2·5 ka but, in the last 2·5 kyr, barrier sedimentation has kept up with and outpaced sea‐level. In the last 0·4 kyr the coastal barrier system has been prograding episodically. Sediment accumulation shows considerable variation, with periods of rapid sediment deposition and periods of non‐deposition or erosion resulting in a highly punctuated sediment record. The study demonstrates how core‐based facies interpretations supported by a high‐resolution chronology and a well‐documented sea‐level history allow identification of depositional environments, erosion surfaces and hiatuses within a very homogeneous stratigraphy, and allow a detailed temporal reconstruction of a coastal barrier system in relation to sea‐level rise and sediment supply.     

Scale considerations in using diatoms as indicators of sea‐level change: lessons from Alaska

This paper assesses variations in quantitative reconstructions of late Holocene relative sea‐level (RSL) change arising from using modern diatom datasets from different spatial scales, applied to case studies from Alaska. We investigate the implications of model choice in transfer functions using local‐, sub‐regional‐ and regional‐scale modern training sets, and produce recommendations on the creation and selection of modern datasets for reconstructing RSL change over Holocene timescales in tidal marsh environments comparable with those in Alaska. We show that regional modern training sets perform best in terms of providing fossil samples with good modern analogues, and in producing reconstructions that most closely match observations, where these are available. Local training sets are frequently insufficient to provide fossil samples with good modern analogues and may over‐estimate the precision of RSL reconstructions. This is particularly apparent when reconstructing RSL change for periods beyond the last century. For reconstructing RSL change we recommend using regional modern training sets enhanced by local samples. Copyright © 2012 John Wiley & Sons, Ltd.