Palaeo-ice sheets are important analogues for understanding contemporary ice sheets, offering a record of ice sheet behaviour that spans millennia. There are two main approaches to reconstructing palaeo-ice sheets. Empirical reconstructions use the available glacial geological and chronological evidence to estimate ice sheet extent and dynamics but lack direct consideration of ice physics. In contrast, numerically modelled simulations implement ice physics, but often lack direct quantitative comparison with empirical evidence. Despite being long identified as a fruitful scientific endeavour, few ice sheet reconstructions attempt to reconcile the empirical and model-based approaches. To achieve this goal, model-data comparison procedures are required. Here, we compare three numerically modelled simulations of the former British–Irish Ice Sheet with the following lines of evidence: (a) position and shape of former margin positions, recorded by moraines; (b) former ice-flow direction and flow-switching, recorded by flowsets of subglacial bedforms; and (c) the timing of ice-free conditions, recorded by geochronological data. These model–data comparisons provide a useful framework for quantifying the degree of fit between numerical model simulations and empirical constraints. Such tools are vital for reconciling numerical modelling and empirical evidence, the combination of which will lead to more robust palaeo-ice sheet reconstructions with greater explicative and ultimately predictive power.

@en

Spinning up a highly complex, coupled Earth system model (ESM) is a time consuming and computationally demanding exercise. For models with interactive ice sheet components, this becomes a major challenge, as ice sheets are sensitive to bidirectional feedback processes and equilibrate over glacial timescales of up to many millennia. This work describes and demonstrates a computationally tractable, iterative procedure for spinning up a contemporary, highly complex ESM that includes an interactive ice sheet component. The procedure alternates between a computationally expensive coupled configuration and a computationally cheaper configuration where the atmospheric component is replaced by a data model. By periodically regenerating atmospheric forcing consistent with the coupled system, the data atmosphere remains adequately constrained to ensure that the broader model state evolves realistically. The applicability of the method is demonstrated by spinning up the preindustrial climate in the Community Earth System Model Version 2 (CESM2), coupled to the Community Ice Sheet Model Version 2 (CISM2) over Greenland. The equilibrium climate state is similar to the control climate from a coupled simulation with a prescribed Greenland ice sheet, indicating that the iterative procedure is consistent with a traditional spin-up approach without interactive ice sheets. These results suggest that the iterative method presented here provides a faster and computationally cheaper method for spinning up a highly complex ESM, with or without interactive ice sheet components. The method described here has been used to develop the climate/ice sheet initial conditions for transient, ice sheet-enabled simulations with CESM2-CISM2 in the Coupled Model Intercomparison Project Phase 6 (CMIP6).

@en

The Greenland ice sheet (GrIS) is now losing mass at a rate of 0.7 mm of sea level rise (SLR) per year. Here we explore future GrIS evolution and interactions with global and regional climate under high greenhouse gas forcing with the Community Earth System Model version 2.1 (CESM2.1), which includes an interactive ice sheet component (the Community Ice Sheet Model v2.1 [CISM2.1]) and an advanced energy balance-based calculation of surface melt. We run an idealized 350-year scenario in which atmospheric CO₂ concentration increases by 1% annually until reaching four times pre-industrial values at year 140, after which it is held fixed. The global mean temperature increases by 5.2 and 8.5 K by years 131–150 and 331–350, respectively. The projected GrIS contribution to global mean SLR is 107 mm by year 150 and 1,140 mm by year 350. The rate of SLR increases from 2 mm yr⁻¹ at year 150 to almost 7 mm yr⁻¹ by year 350. The accelerated mass loss is caused by rapidly increasing surface melt as the ablation area expands, with associated albedo feedback and increased sensible and latent heat fluxes. This acceleration occurs for a global warming of approximately 4.2 K with respect to pre-industrial and is in part explained by the quasi-parabolic shape of the ice sheet, which favors rapid expansion of the ablation area as it approaches the interior “plateau.”.

@en

We describe and evaluate version 2.1 of the Community Ice Sheet Model (CISM). CISM is a parallel, 3-D thermomechanical model, written mainly in Fortran, that solves equations for the momentum balance and the thickness and temperature evolution of ice sheets. CISM's velocity solver incorporates a hierarchy of Stokes flow approximations, including shallow-shelf, depth-integrated higher order, and 3-D higher order. CISM also includes a suite of test cases, links to third-party solver libraries, and parameterizations of physical processes such as basal sliding, iceberg calving, and sub-ice-shelf melting. The model has been verified for standard test problems, including the Ice Sheet Model Intercomparison Project for Higher-Order Models (ISMIP-HOM) experiments, and has participated in the initMIP-Greenland initialization experiment. In multimillennial simulations with modern climate forcing on a 4 km grid, CISM reaches a steady state that is broadly consistent with observed flow patterns of the Greenland ice sheet. CISM has been integrated into version 2.0 of the Community Earth System Model, where it is being used for Greenland simulations under past, present, and future climates. The code is open-source with extensive documentation and remains under active development.@en

This paper summarises developments in understanding sea level change during the Quaternary in Scotland since the publication of the Quaternary of Scotland Geological Conservation Review volume in 1993. We present a review of progress in methodology, particularly in the study of sediments in isolation basins and estuaries as well as in techniques in the field and laboratory, which have together disclosed greater detail in the record of relative sea level (RSL) change than was available in 1993. However, progress in determining the record of RSL change varies in different areas. Studies of sediments and stratigraphy offshore on the continental shelf have increased greatly, but the record of RSL change there remains patchy. Studies onshore have resulted in improvements in the knowledge of rock shorelines, including the processes by which they are formed, but much remains to be understood. Studies of Late Devensian and Holocene RSLs around present coasts have improved knowledge of both the extent and age range of the evidence. The record of RSL change on the W and NW coasts has disclosed a much longer dated RSL record than was available before 1993, possibly with evidence of Meltwater Pulse 1A, while studies in estuaries on the E and SW coasts have disclosed widespread and consistent fluctuations in Holocene RSLs. Evidence for the meltwater pulse associated with the Early Holocene discharge of Lakes Agassiz–Ojibway in N America has been found on both E and W coasts. The effects of the impact of storminess, in particular in cliff-top storm deposits, have been widely identified. Further information on the Holocene Storegga Slide tsunami has enabled a better understanding of the event, but evidence for other tsunami events on Scottish coasts remains uncertain. Methodological developments have led to new reconstructions of RSL change for the last 2000 years, utilising state-of-the-art GIA models and alongside coastal biostratigraphy to determine trends to compare with modern tide gauge and documentary evidence. Developments in GIA modelling have provided valuable information on patterns of land uplift during and following deglaciation. The studies undertaken raise a number of research questions which will require addressing in future work.

@en

Observational evidence, including offshore moraines and sediment cores, confirm that at the Last Glacial Maximum (LGM) the Greenland ice sheet (GrIS) expanded to a significantly larger spatial extent than seen at present, grounding into Baffin Bay and out onto the continental shelf break. Given this larger spatial extent and its close proximity to the neighbouring Laurentide Ice Sheet (LIS) and Innuitian Ice Sheet (IIS), it is likely these ice sheets will have had a strong non-local influence on the spatial and temporal behaviour of the GrIS. Most previous paleo ice-sheet modelling simulations recreated an ice sheet that either did not extend out onto the continental shelf or utilized a simplified marine ice parameterization which did not fully include the effect of ice shelves or neglected the sensitivity of the GrIS to this non-local bedrock signal from the surrounding ice sheets. In this paper, we investigated the evolution of the GrIS over the two most recent glacial-interglacial cycles (240 ka BP to the present day) using the ice-sheet-ice-shelf model IMAU-ICE. We investigated the solid earth influence of the LIS and IIS via an offline relative sea level (RSL) forcing generated by a glacial isostatic adjustment (GIA) model. The RSL forcing governed the spatial and temporal pattern of sub-ice-shelf melting via changes in the water depth below the ice shelves. In the ensemble of simulations, at the glacial maximums, the GrIS coalesced with the IIS to the north and expanded to the continental shelf break to the southwest but remained too restricted to the northeast. In terms of the global mean sea level contribution, at the Last Interglacial (LIG) and LGM the ice sheet added 1.46 and -2.59 m, respectively. This LGM contribution by the GrIS is considerably higher (∼1.26 m) than most previous studies whereas the contribution to the LIG highstand is lower (∼0.7 m). The spatial and temporal behaviour of the northern margin was highly variable in all simulations, controlled by the sub-ice-shelf melting which was dictated by the RSL forcing and the glacial history of the IIS and LIS. In contrast, the southwestern part of the ice sheet was insensitive to these forcings, with a uniform response in all simulations controlled by the surface air temperature, derived from ice cores.

@en

The new sea-level database for Britain and Ireland contains >2100 data points from 86 regions and records relative sea-level (RSL) changes over the last 20 ka and across elevations ranging from ∼+40 to −55 m. It reveals radically different patterns of RSL as we move from regions near the centre of the Celtic ice sheet at the last glacial maximum to regions near and beyond the ice limits. Validated sea-level index points and limiting data show good agreement with the broad patterns of RSL change predicted by current glacial isostatic adjustment (GIA) models. The index points show no consistent pattern of synchronous coastal advance and retreat across different regions, ∼100–500 km scale, indicating that within-estuary processes, rather than decimetre- and centennial-scale oscillations in sea level, produce major controls on the temporal pattern of horizontal shifts in coastal sedimentary environments. Comparisons between the database and GIA model predictions for multiple regions provide potentially powerful constraints on various characteristics of global GIA models, including the magnitude of MWP1A, the final deglaciation of the Laurentide ice sheet and the continued melting of Antarctica after 7 ka BP.

@en

The spatial variability of Holocene relative sea level (RSL) in the South China Sea is unknown, with data restricted to Thailand, the Malay Peninsula, and a few other isolated sites. In this study, we present new continuous RSL records for Borneo using surveyed and U–Th dated coral microatolls from four sites in western Sarawak. The record spans 450 years of RSL from 7450 to 7000 yr BP. Our data suggest that RSL was higher than present and rapid RSL rise had ceased by 7450 yr BP. We compare these RSL reconstructions with a regional model of glacial-isostatic adjustment (GIA). The RSL reconstructions from three sites off the coast of Sarawak show a spatial gradient opposite to that predicted by the GIA model. This disagreement can best be explained by tectonic deformation since 7000 yr BP, which was previously unrecognized. We propose vertical land motion of 0.7–1.45 m due to slip on the Serabang fault, which runs between our four sites. This slip may have occurred in response to the loading of the Sunda Shelf by rising sea level.

@en

Tidal marshes rank among Earth's vulnerable ecosystems, which will retreat if future rates of relative sea-level rise (RSLR) exceed marshes' ability to accrete vertically. Here, we assess the limits to marsh vulnerability by analyzing >780 Holocene reconstructions of tidal marsh evolution in Great Britain. These reconstructions include both transgressive (tidal marsh retreat) and regressive (tidal marsh expansion) contacts. The probability of a marsh retreat was conditional upon Holocene rates of RSLR, which varied between -7.7 and 15.2 mm/yr. Holocene records indicate that marshes are nine times more likely to retreat than expand when RSLR rates are ≥7.1 mm/yr. Coupling estimated probabilities of marsh retreat with projections of future RSLR suggests a major risk of tidal marsh loss in the twenty-first century. All of Great Britain has a >80% probability of a marsh retreat under Representative Concentration Pathway (RCP) 8.5 by 2100, with areas of southern and eastern England achieving this probability by 2040.

@en

Ice sheets are a major component of the Earth System, however they are not yet interactively coupled to most global climate models. Here we present past achievements in this front with the CESM1.0 version as well as first results and challenges with the upcoming CESM2.0, where the Community Ice Sheet Model 2.1 is bi-directionally coupled to the atmospheric and ocean components, as opposed to only one-way coupling in CESM1.0. In both CESM versions, the surface mass balance (SMB) is calculated in the land component (CLM) with explicit albedo and refreezing calculations, and downscaled to the ice sheet model resolution via elevation classes and bi-linear (horizontal) and linear (vertical) interpolations. A major highlight of CESM is that the most important coupling process between ice sheet and atmosphere, the albedo feedback, is explicitly modeled, as opposed to state-of-the-art parameterizations of albedo and/or surface melt. Regarding future Greenland ice sheet projections with CESM1, we summarize our results on: 1) doubled end-of-the-century melt rates under RCP8.5 from increased incoming thermal radiation and turbulent fluxes despite decreased incoming solar radiation over Greenland from more clouds, 2) increase in SMB variability due to reduced accumulation to ablation area ration, 3) bimodal emergence of an anthropogenic signal on the SMB due to both increasing melt and snow accumulation, 4) reduction in ice discharge from marginal thinning. We also outline work in progress in preparation for our contribution to the Ice Sheet Model Intercomparison Project 6 (ISMIP6) and paleo-research on the last deglaciation, e.g., on model initialization, improved SMB calculation, evaluation of CESM2.0 climate over Greenland, and parameter optimization of the higher-order CISM2.1. @en

The ice streams and outlet glaciers of the Greenland Ice sheet (GrIS) are dynamic active parts of the ice sheets; with the recent accelerated contribution to global sea level rise attributed to retreat and thinning from many marine-terminating outlet glaciers. For example, across the Northern margin the outlet glaciers have been seen to be retreating, and increasing in speed. (Hill, et al, 2017). Over longer time scales geomorphological evidence from previously glaciated regions, such as below the Laurentide Ice Sheet (Stokes, Nature, 2016) highlighted the very spatial dynamic nature of the paleo ice streams. In this study, we aim to investigate the variability, both spatially and temporally of the GrIS ice streams and outlet glaciers to try and understand the mechanisms that control the ice stream formation and stability. We perform an ensemble of sensitivity experiments using the new updated version 2.1 of Community Ice sheet model (CISM) for the Greenland ice sheet, at a 4kmx4km resolution. In all simulations, we force the ice sheet with the surface mass Balance (SMB) taken from RACMO2.3 with the aim of recreating a steady-state present day ice sheet. There have been a number of improvements since version2 which we will adopt: a higher-order velocity solver ‘vertically integrated approximation to the Stokes flow law’ taken from Goldberg, 2011; the introduction of a new pseudo-plastic power law sliding scheme, that relates changes in the basal shear stress to the till yield stress and a pseudo-plastic flow law exponent and finally a simple model of basal hydrology which is used to calculate the changes in till yield stress. This basal hydrology model relates changes in rheology of subglacial sediments (defined by a till friction angle) to changes in basal water pressure to produce a till yield stress that varies both in space and time. This new sliding law and basal hydrology model have a number of@en

Sea-level rise is a global problem, yet to forecast future changes, we must understand how and why relative sea level (RSL) varied in the past, on local to global scales. In East and Southeast Asia, details of Holocene RSL are poorly understood. Here we present two independent high-resolution RSL proxy records from Belitung Island on the Sunda Shelf. These records capture spatial variations in glacial isostatic adjustment and paleotidal range, yet both reveal a RSL history between 6850 and 6500 cal years BP that includes two 0.6 m fluctuations, with rates of RSL change reaching 13±4 mm per year (2σ). Observations along the south coast of China, although of a lower resolution, reveal fluctuations similar in amplitude and timing to those on the Sunda Shelf. The consistency of the Southeast Asian records, from sites 2,600 km apart, suggests that the records reflect regional changes in RSL that are unprecedented in modern times.

@en