Glacial Isostatic Adjustment Models Research Articles

Constraining Upper Mantle Viscosity Using Temperature and Water Content Inferred From Seismic and Magnetotelluric Data

AbstractMantle viscosity controls a variety of geodynamic processes such as glacial isostatic adjustment (GIA), but it is poorly constrained because it cannot be measured directly from geophysical measurements. Here we develop a method that calculates viscosity using empirical viscosity flow laws coupled with mantle parameters (temperature and water content) inferred from seismic and magnetotelluric (MT) observations. We find that combining geophysical constraints allows us to place significantly tighter bounds on viscosity estimates compared to using seismic or MT observations alone. In particular, electrical conductivity inferred from MT data can determine whether upper mantle minerals are hydrated, which is important for viscosity reduction. Additionally, we show that rock composition should be considered when estimating viscosity from geophysical data because composition directly affects seismic velocity and electrical conductivity. Therefore, unknown composition increases uncertainty in temperature and water content, and makes viscosity more uncertain. Furthermore, calculations that assume pure thermal control of seismic velocity may misinterpret compositional variations as temperature, producing erroneous interpretations of mantle temperature and viscosity. Stress and grain size also affect the viscosity and its associated uncertainty, particularly via their controls on deformation regime. Dislocation creep is associated with larger viscosity uncertainties than diffusion creep. Overall, mantle viscosity can be estimated best when both seismic and MT data are available and the mantle composition, grain size and stress can be estimated. Collecting additional MT data probably offers the greatest opportunity to improve geodynamic or GIA models that rely on viscosity estimates.

Exceptionally stable preindustrial sea level inferred from the western Mediterranean Sea.

An accurate record of preindustrial (pre-1900 CE) sea level is necessary to contextualize modern global mean sea level (GMSL) rise with respect to natural variability. Precisely dated phreatic overgrowths on speleothems (POS) provide detailed rates of Late Holocene sea-level rise in Mallorca. Statistical analysis indicates that sea level rose locally by 0.12 to 0.31 m (95% confidence) from 3.26 to 2.84 thousand years (ka) ago (2σ) and remained within 0.08 m (95% confidence) of preindustrial levels from 2.84 ka to 1900 CE. This sea-level history is consistent with glacial isostatic adjustment models adopting relatively weak upper mantle viscosities of ~1020 Pa s. There is virtual certainty (>0.999 probability) that the average GMSL rise since 1900 CE has exceeded even the high average rate of sea-level rise between 3.26 and 2.84 ka inferred from the POS record. We conclude that modern GMSL rise is anomalous relative to any natural variability in ice volumes over the past 4000 years.

Last Interglacial sea-level data points from Northwest Europe

Abstract. Abundant numbers of sites and studies exist in NW Europe that document the geographically and geomorphologically diverse coastal record from the Last Interglacial (Eemian, Ipswichian, Marine Isotope Stage 5e). This paper summarises a database of 146 known Last Interglacial sea-level data points from in and around the North Sea (35 entries in the Netherlands, 10 Belgium, 23 in Germany, 17 in Denmark, 9 in Britain) and the English Channel (24 entries for the British and 25 for the French side, 3 on the Channel Isles) believed to be a representative and fairly complete inventory and assessment from ∼80 published sites. The geographic distribution (∼1500 km SW–NE) across the near field of the Scandinavian and British ice sheets and the attention paid to relative and numeric age control are assets of the NW European database. The research history of Last Interglacial coastal environments and sea level for this area is long, methodically diverse and spread through regional literature in several languages. Our review and database compilation effort drew from the original regional literature and paid particular attention to distinguishing between sea-level index points (SLIPs) and marine and terrestrial limiting points. We also incorporated an updated quantification of background rates of basin subsidence for the central and eastern North Sea region, utilising revised mapping of the base Quaternary, to correct for significant basin subsidence in this depocentre. As a result of subsidence, lagoonal and estuarine Last Interglacial shorelines of the Netherlands and the German Bight are preserved below the surface. In contrast, Last Interglacial shorelines along the English Channel are encountered above modern sea level. This paper describes the dominant sea-level indicators from the region compliant with the WALIS database structure and referenced to original data sources (https://doi.org/10.5281/zenodo.6478094, Cohen et al., 2021). The sea-level proxies are mostly obtained from locations with good lithostratigraphic, morphostratigraphic and biostratigraphical constraints. Most continental European sites have chronostratigraphic age control, notably through regional pollen association zones with duration estimates. In all regions, many SLIPs and limiting points have further independent age control from luminescence, uranium series, amino acid racemisation and electron spin resonance dating techniques. Main foreseen usage of this database for the near-field region of the European ice sheets is in glacial isostatic adjustment modelling and fingerprinting Last Interglacial ice sheet melt.

Revisiting sea-level budget by considering all potential impact factors for global mean sea-level change estimation

Accurate estimates of global sea-level change from the observations of Altimetry, Argo and Gravity Recovery and Climate Experiment (GRACE) and GRACE Follow-on (GRACE-FO) are of great value for investigating the global sea-level budget. In this study, we analyzed the global sea-level change over the period from January 2005 to December 2019 by considering all potential impact factors, i.e. three factors for Altimetry observations (two Altimetry products, ocean bottom deformation (OBD) and glacial isostatic adjustment (GIA)), three factors for Argo observations (four Argo products, salinity product error and deep-ocean steric sea-level change), and seven factors for GRACE/GRACE-FO observations including three official RL06 solutions, five spatial filtering methods, three GIA models, two C20 (degree 2 order 0) products, Geocenter motion, GAD field and global mass conservation. The seven impact factors of GRACE/GRACE-FO observations lead to ninety combinations for the post-procession of global mean barystatic sea-level change estimation, whose rates range from 2.00 to 2.45 mm/year. The total uncertainty of global barystatic sea-level change rate is ± 0.27 mm/year at the 95% confidence level, estimated as the standard deviation of the differences between the different datasets constituting the ensembles. The statistical results show that the preferred GIA model developed by Caron et al. in 2018 can improve the closure of the global sea-level budget by 0.20–0.30 mm/year, which is comparable with that of neglecting the halosteric component. About 30.8% of total combinations (GRACE/GRACE-FO plus Argo) can close the global sea-level budget within 1-sigma (0.23 mm/year) of Altimetry observations, 88.9% within 2-sigma. Once the adopted factors including GRACE/GRACE-FO solutions from Center for Space Research (CSR), Caron18 GIA model, SWENSON filtering and Argo product from China Second Institute of Oceanography, the linear trend of global sterodynamic sea-level change derived from GRACE/GRACE-FO plus Argo observations is 3.85 ± 0.14 mm/year, nearly closed to 3.90 ± 0.23 mm/year of Altimetry observations.

Resolving glacial isostatic adjustment (GIA) in response to modern and future ice loss at marine grounding lines in West Antarctica

Abstract. Accurate glacial isostatic adjustment (GIA) modelling in the cryosphere is required for interpreting satellite, geophysical and geological records and for assessing the feedbacks of Earth deformation and sea-level change on marine ice-sheet grounding lines. GIA modelling in areas of active ice loss in West Antarctica is particularly challenging because the ice is underlain by laterally varying mantle viscosities that are up to several orders of magnitude lower than the global average, leading to a faster and more localised response of the solid Earth to ongoing and future ice-sheet retreat and necessitating GIA models that incorporate 3-D viscoelastic Earth structure. Improvements to GIA models allow for computation of the viscoelastic response of the Earth to surface ice loading at sub-kilometre resolution, and ice-sheet models and observational products now provide the inputs to GIA models at comparably unprecedented detail. However, the resolution required to accurately capture GIA in models remains poorly understood, and high-resolution calculations come at heavy computational expense. We adopt a 3-D GIA model with a range of Earth structure models based on recent seismic tomography and geodetic data to perform a comprehensive analysis of the influence of grid resolution on predictions of GIA in the Amundsen Sea Embayment (ASE) in West Antarctica. Through idealised sensitivity testing down to sub-kilometre resolution with spatially isolated ice loading changes, we find that a grid resolution of ∼ 13 of the radius of the load or higher is required to accurately capture the elastic response of the Earth. However, when we consider more realistic, spatially coherent ice loss scenarios based on modern observational records and future ice-sheet model projections and adopt a viscoelastic Earth, we find that predicted deformation and sea-level change along the grounding line converge to within 5 % with grid resolutions of 7.5 km or higher, and to within 2 % for grid resolutions of 3.75 km and higher, even when the input ice model is on a 1 km grid. Furthermore, we show that low mantle viscosities beneath the ASE lead to viscous deformation that contributes to the instrumental record on decadal timescales and equals or dominates over elastic effects by the end of the 21st century. Our findings suggest that for the range of resolutions of 1.9–15 km that we considered, the error due to adopting a coarser grid in this region is negligible compared to the effect of neglecting viscous effects and the uncertainty in the adopted mantle viscosity structure.

Antarctic Basal Water Storage Variation Inferred from Multi-Source Satellite Observation and Relevant Models

Antarctic basal water storage variation (BWSV) refers to mass changes of basal water beneath the Antarctic ice sheet (AIS). Identifying these variations is critical for understanding Antarctic basal hydrology variations and basal heat conduction, yet they are rarely accessible due to a lack of direct observation. This paper proposes a layered gravity density forward/inversion iteration method to investigate Antarctic BWSV based on multi-source satellite observations and relevant models. During 2003–2009, BWSV increased at an average rate of 43 ± 23 Gt/yr, which accounts for 29% of the previously documented total mass loss rate (−76 ± 20 Gt/yr) of AIS. Major uncertainty arises from satellite gravimetry, satellite altimetry, the glacial isostatic adjustment (GIA) model, and the modelled basal melting rate. We find that increases in basal water mainly occurred in regions with widespread active subglacial lakes, such as the Rockefeller Plateau, Siple Coast, Institute Ice Stream regions, and marginal regions of East Antarctic Ice Sheet (EAIS), which indicates the increased water storage in these active subglacial lakes, despite the frequent water drainage events. The Amundsen Sea coast experienced a significant loss during the same period, which is attributed to the basal meltwater discharging into the Amundsen Sea through basal channels.

First a wudd, and syne a sea: postglacial coastal change of Scotland recalled in ancient stories

ABSTRACT The question of whether culturally-grounded stories (myths and legends) have any value in understanding past environmental changes is uncertain. Focused on stories that plausibly recall movements of the land-sea interface in Scotland in postglacial times, this paper summarises details of 11 stories that refer to submergence and 7 stories that refer to emergence. Most submergence stories are confined to the Outer Hebrides and include those recalling when it was possible to walk between places that are now islands. Emergence stories come from around the Scottish coast and include some from the Inner Hebrides. The agreement in direction of movement (submergence or emergence) with models of postglacial landscape change is almost perfect. By comparing submergence/emergence magnitudes to histories of relative sea-level change from glacial-isostatic adjustment models, it is possible to estimate ages for all the submergence stories to at least 2107­-8695 years BP and all the emergence stories to 674-7120 years BP. Land-uplift rates calculated from emergence stories agree with those from palaeo-shoreline analysis. As is becoming increasingly clear for other places to which ancient culturally-grounded stories about environmental change refer, these Scottish stories likely represent residues of millennia-old observations of coastal change. This study should encourage further investigations of ancient Scottish narratives.

What are the drivers of Caspian Sea level variation during the late Quaternary?

Quaternary Caspian Sea level variations depended on geophysical processes (affecting the opening and closing of gateways and basin size/shape) and hydro-climatological processes (affecting water balance). Disentangling the drivers of past Caspian Sea level variation, as well as the mechanisms by which they impacted the Caspian Sea level variation, is much debated. In this study we examine the relative impacts of hydroclimatic change, ice-sheet accumulation and melt, and isostatic adjustment on Caspian Sea level change. We performed model analysis of ice-sheet and hydroclimate impacts on Caspian Sea level and compared these with newly collated published palaeo-Caspian sea level data for the last glacial cycle. We used palaeoclimate model simulations from a global coupled ocean-atmosphere-vegetation climate model, HadCM3, and ice-sheet data from the ICE-6G_C glacial isostatic adjustment model. Our results show that ice-sheet meltwater during the last glacial cycle played a vital role in Caspian Sea level variations, which is in agreement with hypotheses based on palaeo-Caspian Sea level information. The effect was directly linked to the reorganization and expansion of the Caspian Sea palaeo-drainage system resulting from topographic change. The combined contributions from meltwater and runoff from the expanded basin area were primary factors in the Caspian Sea transgression during the deglaciation period between 20 and 15 kyr BP. Their impact on the evolution of Caspian Sea level lasted until around 13 kyr BP. Millennial scale events (Heinrich events and the Younger Dryas) negatively impacted the surface water budget of the Caspian Sea but their influence on Caspian Sea level variation was short-lived and was outweighed by the massive combined meltwater and runoff contribution over the expanded basin.

Postglacial relative sea level change in Norway

We present the first comprehensive postglacial relative sea-level (RSL) database for the Norwegian coast from Oslo to the Kola Peninsula. The database spans the last 20 kyrs and is composed of 413 index points and 610 limiting data points derived from raised beaches, archeological data, glaciomarine terraces, and sedimentary indicators from isolation basins, salt marshes, and peat bogs. The data are quality controlled, assigned standardized indicative meanings, and recalibrated to current standards. We use an ensemble of Bayesian statistical models, trained on sea-level index points and weighted by their fit to both index and sea-level limiting points, to assess the spatiotemporal patterns of Norwegian RSL change. Continuous RSL fall driven by isostatic rebound in response to Eurasian ice sheet collapse dominates the RSL signal at every inland location in Norway. A first transgression (episode of RSL rise), which occurred in southwest Norway during the Younger Dryas (14–11.7 ka), increased RSL by as much as 15 m in some locations. A second transgression, named the Tapes transgression, occurred during the early-mid Holocene between 10 and 5.5 ka. The spatiotemporal model ensemble constrains the timing, amplitude, and spatial distribution of the Younger Dryas and Tapes transgressions. Based on our modeling results, we speculate that the Tapes Transgression was the result of global mean sea level rise temporarily outpacing isostasy-driven RSL fall, while the Younger Dryas Transgression was likely the result of local ice sheet readvance combined with low viscosity asthenosphere and weak lithosphere in the region. We also describe the effects of peripheral bulge migration on Norwegian RSL, which caused increased RSL in the early Holocene and a delayed Tapes transgression. We show that postglacial RSL data in Norway contain complex spatiotemporal patterns of nearfield RSL change that can best be estimated by combining a high-quality data compilation with glacial isostatic adjustment modeling via a robust statistical model.

SeeLevelViz: A simple data science tool for dynamic visualization of shoreline displacement caused by sea-level change

We present SeeLevelViz, a free, open-source program written in Python for making interactive visualizations of relative sea-level change in landscapes and shorelines. The accurate reconstruction of shoreline positions is a crucial factor in coastal palaeolandscape studies, particularly in areas where the coast is fronted by islands, since the separation of islands from the mainland drives important ecological and sociocultural outcomes. This program creates accurate time-slice reconstructions of shoreline positions and palaeolandscapes when the user provides two components: 1) a digital elevation model of the target region (including currently submerged areas), and 2) a simple spreadsheet of relative sea-level elevations at different dates derived either from a glacio isostatic adjustment model of relative sea-level change, or from observed past sea-level data points. The tool is presented using the eastern coast of the Adriatic Sea in the Mediterranean as a test case, since this region has a complex coastline articulation due to combined geological and geomorphological factors. In this area, like in many other Mediterranean coastal areas, the separation of islands from the mainland following the last glacial maximum and throughout the Holocene has occurred in connection with important phases of the development, particularly of Mesolithic and Neolithic cultures, influencing human migrations and the spread of seafaring techniques. Reliable palaeolandscape reconstructions at different time slices are thus crucial for supporting archaeological interpretation. Flexibile and user-friendly, SeeLevelViz can compliment reconstructions of coastal landscape changes either based on glacial isostatic adjustment models or on relative palaeo-sea- level evidence, since simple, interactive visualizations are a powerful technique for understanding spatial time-series data, both for the interpretation phase of research, and for presentation to colleagues and the public. The program can be modified or used freely for papers, presentations, etc. By crediting and citing this article.

Palaeo-Shoreline Configuration of the Adventure Plateau (Sicilian Channel) at the Last Glacial Maximum

The Adventure Plateau, located in the NW sector of the Sicilian Channel, experienced several episodes of exposure/erosion and subsequent drowning, with the most recent occurring after the Last Glacial Maximum (LGM). Unlike other parts of the Sicilian Channel, the Adventure Plateau is relatively tectonically stable and is therefore best suitable for reconstructing its coastal configuration before the post-LGM marine transgression. Here, we use high-resolution seismic data to identify and map the palaeo-coastline at the LGM on the basis of the internal architecture of the prograding wedges (i.e., the location of the subaqueous clinoform rollover point) and the erosional markers such as the subaerial unconformities and the wave ravinement surfaces. These data, which show an extreme variability in the palaeo-morphology of the coastal margins of the Adventure Plateau, have been complemented with vintage seismic profiles in order to entirely cover its perimeter. The mapped LGM coastline has then been compared to predictions from glacial isostatic adjustment (GIA) modeling, which considers the horizontal migration of the shorelines in response to sea level rise and to Earth’s rotational and deformational effects associated with deglaciation. The two shorelines (i.e., the coastline derived from the marine data interpretation and the one derived from the GIA model) are in good agreement at 21 kyears BP, although some discrepancies occur in the southern part of the plateau, where the seabed slope is extremely gentle, which makes the clinoform rollover points and the buried erosional unconformities difficult to detect. After 20 kyears BP, an acceleration in the rate of the sea level rise occurred. The results of this study indicate the importance of comparing experimental data with model predictions in order to refine and calibrate boundary parameters and to gain a better picture of the evolution of sea level rise over various time scales.

Influence of 3D Earth Structure on Glacial Isostatic Adjustment in the Russian Arctic

AbstractAnalyses of glacial isostatic adjustment (GIA) and deglacial relative sea‐level (RSL) change in the Russian Arctic deliver important insights into the Earth's viscosity structure and the deglaciation history of the Eurasian ice sheet complex. Here, we validate the 1D GIA models ICE‐6G_C (VM5a) and ICE‐7G_NA (VM7) and select new 3D GIA models in the Russian Arctic against a quality‐controlled deglacial RSL database of >500 sea‐level data points from 24 regions. Both 1D models correspond to the RSL data along the southern coast of the Barents Sea and Franz Josef Land from ∼11 ka BP to present but show notable misfits (>50 m at 10 ka BP) with the White Sea data. We find 3D model predictions of deglacial RSL resolve most of the misfits with the observed data for the White Sea while retaining comparable fits in other regions of the Russian Arctic. Our results further reveal: (a) RSL in the western Russian Arctic is sensitive to elastic lithosphere with lateral thickness variation and 3D viscosity structure in the upper mantle; and (b) RSL in the whole Russian Arctic is less sensitive to 3D viscosity structure in the lower mantle compared to the upper mantle. The 3D models reveal a compromise in the upper mantle between the background viscosity and scaling factor to best fit the RSL data, which needs to be considered in future 3D GIA studies.

GPS Rates of Vertical Bedrock Motion Suggest Late Holocene Ice‐Sheet Readvance in a Critical Sector of East Antarctica

AbstractWe investigate present‐day bedrock vertical motion using new Global Positioning System (GPS) timeseries from the Totten‐Denman glacier region of East Antarctica (∼77°–120°E) where models of glacial isostatic adjustment (GIA) disagree, glaciers are likely losing mass, and few data constraints on GIA exist. We show that varying surface mass balance loading (SMBL) is a dominant signal, contributing random‐walk‐like noise to GPS timeseries across Antarctica. In the study region, it induces site velocity biases of up to ∼+1 mm/yr over 2010–2020. After correcting for SMBL displacement and GPS common mode error, subsidence is evident at all sites aside from the Totten Glacier region where uplift is ∼1.5 mm/yr. Uplift near the Totten Glacier is consistent with late Holocene ice retreat while the widespread subsidence further west suggests possible late Holocene readvance of the region’s ice sheet, in broad agreement with limited glacial geological data and highlighting the need for sampling beneath the current ice sheet.

Assessment of Contemporary Antarctic GIA Models Using High-Precision GPS Time Series

Past redistributions of the Earth’s mass resulting from the Earth’s viscoelastic response to the cycle of deglaciation and glaciation reflect the process known as glacial isostatic adjustment (GIA). GPS data are effective at constraining GIA velocities, provided that these data are accurate, have adequate spatial coverage, and account for competing geophysical processes, including the elastic loading of ice/snow ablation/accumulation. GPS solutions are significantly affected by common mode errors (CMEs) and the choice of optimal noise model, and they are contaminated by other geophysical signals due primarily to the Earth’s elastic response. Here, independent component analysis is used to remove the CMEs, and the Akaike information criterion is used to determine the optimal noise model for 79 GPS stations in Antarctica, primarily distributed across West Antarctica and the Antarctic Peninsula. Next, a high-resolution surface mass variation model is used to correct for elastic deformation. Finally, we use the improved GPS solution to assess the accuracy of seven contemporary GIA forward models in Antarctica. The results show that the maximal GPS crustal displacement velocity deviations reach 4.0 mm yr−1, and the mean variation is 0.4 mm yr−1 after removing CMEs and implementing the noise analysis. All GIA model-predicted velocities are found to systematically underestimate the GPS-observed velocities in the Amundsen Sea Embayment. Additionally, the GPS vertical velocities on the North Antarctic Peninsula are larger than those on the South Antarctic Peninsula, and most of the forward models underestimate the GIA impact on the Antarctic Peninsula.

An approach for constraining mantle viscosities through assimilation of palaeo sea level data into a glacial isostatic adjustment model

Abstract. Glacial isostatic adjustment is largely governed by the rheological properties of the Earth's mantle. Large mass redistributions in the ocean–cryosphere system and the subsequent response of the viscoelastic Earth have led to dramatic sea level changes in the past. This process is ongoing, and in order to understand and predict current and future sea level changes, the knowledge of mantle properties such as viscosity is essential. In this study, we present a method to obtain estimates of mantle viscosities by the assimilation of relative sea level rates of change into a viscoelastic model of the lithosphere and mantle. We set up a particle filter with probabilistic resampling. In an identical twin experiment, we show that mantle viscosities can be recovered in a glacial isostatic adjustment model of a simple three-layer Earth structure consisting of an elastic lithosphere and two mantle layers of different viscosity. We investigate the ensemble behaviour on different parameters in the following three set-ups: (1) global observations data set since last glacial maximum with different ensemble initialisations and observation uncertainties, (2) regional observations from Fennoscandia or Laurentide/Greenland only, and (3) limiting the observation period to 10 ka until the present. We show that the recovery is successful in all cases if the target parameter values are properly sampled by the initial ensemble probability distribution. This even includes cases in which the target viscosity values are located far in the tail of the initial ensemble probability distribution. Experiments show that the method is successful if enough near-field observations are available. This makes it work best for a period after substantial deglaciation until the present when the number of sea level indicators is relatively high.

A GNSS velocity field for crustal deformation studies: The influence of glacial isostatic adjustment on plate motion models

SUMMARYThe two main causes of global-scale secular deformation of the Earth are tectonic plate motion and glacial isostatic adjustment (GIA). We create a bespoke global 3D GNSS surface velocity field ‘NCL20’ to investigate tectonic plate motion and the effect of GIA on plate motion models (PMMs), drawing on a set of 1D and 3D GIA model predictions. The main motivation for creating NCL20 is to include a larger number of GNSS sites in the most GIA-affected areas of investigation, namely North America, Europe, and Antarctica. We do this using the IGS repro2 data and other similarly processed GNSS data sets. Our final GNSS velocity field has horizontal uncertainties mostly within ±0.5 mm yr–1 and vertical uncertainties mostly within ±1 mm yr–1 (at 95 per cent confidence), which make it suitable for testing GIA models. We generate a suite of 117 global GIA model predictions by combining three different ice history models (ICE-5G, ICE-6G and W12) with a range of 1D and 3D Earth models. By subtracting this ensemble from the GNSS velocity field, we identify and compare a range of PMMs which are expected to be unaffected by GIA. Our method allows us to include GNSS sites that are typically excluded from PMM estimations due to their location in GIA-affected regions. We demonstrate that significant GIA-related horizontal motion outside of the rapidly uplifting regions may bias PMMs if left uncorrected. Based on their ability to explain the observed surface velocity field, a group of best-performing GIA models is selected for three regions of interest: North America, Europe and Antarctica. The range of GIA predictions in each best-performing group is assumed to represent the uncertainty in regional GIA modelling insofar as it can be constrained by present-day geodetic velocities. In the horizontal component, we note that 3D GIA models show more variation in the direction of predicted velocities than 1D GIA models, confirming that horizontal velocities are strongly sensitive to lateral variations in Earth structure. Furthermore, for Antarctica the variation in predicted GIA vertical velocities suggests that the total GIA contribution to annual gravimetric mass change ranges from –3 to 22 Gt yr–1 depending on which of the best-performing GIA models is used.

Investigating the roles of relative sea-level change and glacio-isostatic adjustment on the retreat of a marine based ice stream in NW Scotland

The record of ice-sheet demise since the last glacial maximum (LGM) provides an opportunity to test the relative importance of instability mechanisms, including relative sea-level (RSL) change, controlling ice-sheet retreat. Here we examine the record of RSL changes accompanying the retreat of the Minch Ice Stream (MnIS) of northwest Scotland during the deglaciation following the LGM as well as use the record to provide additional age constraints on a local late-glacial readvance known as the Wester Ross Readvance. We use new and existing records of RSL change obtained from isolation basins in Wester Ross along the flanks of the former MnIS to test available glacial-isostatic adjustment (GIA) predictions of the deglacial RSL history for the region. Using these GIA model predictions we examine the nature of RSL change across the retreating front of the MnIS through the early deglaciation. Our new radiocarbon ages from these basins confirm the timing of deglaciation within the inner trough of the former MnIS as well as refines the age of the Wester Ross Readvance, both established by earlier cosmogenic-based studies. We find that the Wester Ross Readvance culminated around 15.8 ± 0.1 ka, slightly earlier than recent suggestions. Near Gairloch, Wester Ross, RSL fell from a marine limit ∼20 m above present at ∼16.1–16.5 ka. Three isolation basins record RSL fall over the following ∼0.8 ka allowing a comparison between GIA predictions and RSL observations. Our new analyses suggest that the rate of RSL rise increased at the ice front, in concert with the MnIS encountering a landward sloping bed potentially aiding the rapid retreat of the MnIS from 17.6 to 16.4 ka BP. This observation suggests that GIA during deglaciation does not necessarily induce a stabilizing RSL change to marine-based ice streams as some models have suggested. Along indented ice margins, the RSL field at the front of individual ice streams may be governed by the regional GIA signal driven by the ice sheet as a whole, rather than the local ice front. In addition, the stabilizing impact of post-glacial rebound is dependent on an Earth rheology weak enough to respond quickly to the ice-sheet retreat. In the case of the MnIS, the RSL experienced at the front of the ice stream was likely governed by the earlier ice mass extent, the larger ice masses lying to the east and south of the highly indented ice front, and the relatively strong Earth rheology beneath the British Isles. Thus, the geometry of the ice sheet margins, such as those in Greenland and Antarctica today, and the Earth rheology beneath them need to be taken into account when considering the stabilizing impact of post-glacial rebound on marine ice sheet retreat.

Rapid retreat of the southwestern Laurentide Ice Sheet during the Bølling-Allerød interval

Abstract The timing of Laurentide Ice Sheet deglaciation along its southwestern margin controlled the evolution of large glacial lakes and has implications for human migration into the Americas. Accurate reconstruction of the ice sheet's retreat also constrains glacial isostatic adjustment models and is important for understanding ice-sheet sensitivity to climate forcing. Despite its significance, retreat of the southwestern Laurentide Ice Sheet (SWLIS) is poorly constrained by minimum-limiting 14C data. We present 26 new cosmogenic 10Be exposure ages spanning the western Interior Plains, Canada. Using a Bayesian framework, we combine these data with geomorphic mapping, 10Be, and high-quality minimum-limiting 14C ages to provide an updated chronology. This dataset presents an internally consistent retreat record and indicates that the initial detachment of the SWLIS from its convergence with the Cordilleran Ice Sheet began by ca. 15.0 ka, concurrent with or slightly prior to the onset of the Bølling-Allerød interval (14.7–12.9 ka) and retreated >1200 km to its Younger Dryas (YD) position in ~2500 yr. Ice-sheet stabilization at the Cree Lake Moraine facilitated a meltwater drainage route to the Arctic from glacial Lake Agassiz within the YD, but not necessarily at the beginning. Our record of deglaciation and new YD constraints demonstrate deglaciation of the Interior Plains was ~60% faster than suggested by minimum 14C constraints alone. Numerical modeling of this rapid retreat estimates a loss of ~3.7 m of sea-level equivalent from the SWLIS during the Bølling-Allerød interval.

The Impact of a 3‐D Earth Structure on Glacial Isostatic Adjustment in Southeast Alaska Following the Little Ice Age

AbstractIn Southeast Alaska, extreme uplift rates are primarily caused by glacial isostatic adjustment (GIA), as a result of ice thickness changes from the Little Ice Age to the present combined with a low‐viscosity asthenosphere. Previous GIA models adopted a 1‐D Earth structure. However, the actual Earth structure is likely more complex due to the long history of subduction and tectonism and the transition from a continental to an oceanic plate. Seismic evidence indeed shows a laterally heterogenous Earth structure. In this study, a numeral model is constructed for Southeast Alaska, which allows for the inclusion of lateral viscosity variations. The viscosity follows from scaling relationships between seismic velocity anomalies and viscosity variations. We use this scaling relationship to constrain the thermal effect on seismic variations and investigate the importance of lateral viscosity variations. We find that a thermal contribution to seismic anomalies of 10% is required to explain the GIA observations. This implies that non‐thermal effects control seismic anomaly variations in the shallow upper mantle. Due to the regional geologic history, it is likely that hydration of the mantle impact both viscosity and seismic velocity. The best‐fit model has a background viscosity of 5.0 × 1019 Pa‐s, and viscosities at ∼80 km depth range from 1.8 × 1019 to 4.5 × 1019 Pa‐s. A 1‐D averaged version of the 3‐D model performed slightly better, however, the two models were statistically equivalent within a 2σ measurement uncertainty. Thus, lateral viscosity variations do not contribute significantly to the uplift rates measured with the current accuracy and distribution of sites.

A review of the composition and chemistry of peridotite mantle xenoliths in volcanic rocks from Antarctica and their relevance to petrological and geophysical models for the lithospheric mantle

Abstract This chapter reviews the geochemistry and petrology of mantle peridotite xenoliths from across Antarctica, including parameters that are of most relevance to geophysical studies. This Memoir is the first time such a complete overview of the chemistry of Antarctic mantle xenoliths has been available and Antarctica should no longer be the ignored continent in studies of mantle xenoliths in volcanic rocks. Xenoliths indicate that the chemistry, heat flow and water content of the Antarctic lithospheric mantle varies regionally at scales of one to thousands of kilometres. The prevalence of variability in xenoliths suggests that the Antarctic mantle is ubiquitously heterogeneous. This has important, yet unquantified, implications for interpreting geophysical data and for reference Earth models used in Antarctic geophysical studies. Information about and interpretations of Antarctic mantle xenoliths can be linked to studies from once adjacent continental blocks in Africa, India, Australia, New Zealand and South America. Together, this can improve understanding of the mantle contribution to glacial isostatic adjustment and geodynamic models to show how the Antarctic mantle fits with adjacent continents in the puzzle of lithospheric blocks. Numerous, fundamental and important research questions remain unanswered making further study of the Antarctic mantle an exciting prospect for future research.