Glacial isostatic adjustment (GIA) has a stabilizing effect on the evolution of the Antarctic ice sheet by reducing the grounding line migration following ice melt. The timescale and strength of this feedback depends on the spatially varying viscosity of the Earth's mantle. Most studies assume a relatively long and laterally homogenous response time of the bedrock. However, the mantle viscosity is spatially variable, with a high mantle viscosity beneath East Antarctica and a low mantle viscosity beneath West Antarctica. For this study, we have developed a new method to couple a 3D GIA model and an ice sheet model to study the interaction between the solid Earth and the Antarctic ice sheet during the last glacial cycle. With this method, the ice sheet model and GIA model exchange ice thickness and bedrock elevation during a fully coupled transient experiment. The feedback effect is taken into account with a high temporal resolution, where the coupling time steps between the ice sheet and GIA model are 5000 years over the glaciation phase and vary between 500 and 1000 years over the deglaciation phase of the last glacial cycle. During each coupling time step, the bedrock elevation is adjusted at every ice sheet model time step, and the deformation is computed for a linearly changing ice load. We applied the method using the ice sheet model ANICE and a 3D GIA finite element model. We used results from a regional seismic model for Antarctica embedded in the global seismic model SMEAN2 to determine the patterns in the mantle viscosity. The results of simulations over the last glacial cycle show that differences in mantle viscosity of an order of magnitude can lead to differences in the grounding line position up to 700gkm and to differences in ice thickness of the order of 2gkm for the present day near the Ross Embayment. These results underline and quantify the importance of including local GIA feedback effects in ice sheet models when simulating the Antarctic ice sheet evolution over the last glacial cycle.

@en

The Earth’s surface and interior deform due to a changing load of the Antarctic Ice Sheet (AIS) during the last glacial cycle, called Glacial Isostatic Adjustment (GIA). This deformation changes the surface height of the ice sheet and indirectly the groundling line position. These changes in surface height and grounding line position influence the evolution of the AIS and consequently, again the load on the Earth’s surface. As a result, GIA operates as a negative feedback loop and could stabilize the evolution of the AIS. This feedback maybe particularly relevant for relatively low viscosities of the mantle in West Antarctica which lead to a relatively fast response time of the bedrock due to changes in the West Antarctic Ice Sheet loading. Most studies capture this process by ignoring lateral variations in the viscosity of the mantle and the stabilizing GIA feedback loop. Here we present a new method to couple an ice sheet model to a GIA model at a variable timestep in the order of a thousand years. Several experiments have been done using different radial and lateral varying rheologies for simulations of the last glacial cycle. It is shown that the effect of including lateral variations and accounting for the stabilizing GIA feedback is up to 80 kilometers for the grounding line position and 400 meters for the ice thickness. The largest differences are observed close to the grounding line of the Ronne ice shelf and at several locations in East Antarctica. The total ice volume of the AIS increases by 0.5 percent over 5000 years when including the 3D GIA feedback loops in the coupled model. These results quantify the local importance of including GIA feedback effects in ice dynamic models when simulating the Antarctic Ice Sheet evolution over the full glacial cycle. @en

Many processes affect sea level near the coast. In this paper, we discuss the major uncertainties in coastal sea-level projections from a process-based perspective, at different spatial and temporal scales, and provide an outlook on how these uncertainties may be reduced. Uncertainty in centennial global sea-level rise is dominated by the ice sheet contributions. Geographical variations in projected sea-level change arise mainly from dynamical patterns in the ocean response and other geophysical processes. Finally, the uncertainties in the short-duration extreme sea-level events are controlled by near coastal processes, storms and tides.@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

Sea-level indicators dated to the Last Interglacial, or Marine Isotope Stage (MIS) 5e, have a twofold value. First, they can be used to constrain the melting of Greenland and Antarctic Ice Sheets in response to global warming scenarios. Second, they can be used to calculate the vertical crustal rates at active margins. For both applications, the contribution of glacio- and hydro-isostatic adjustment (GIA) to vertical displacement of sea-level indicators must be calculated. In this paper, we re-assess MIS 5e sea-level indicators at 11 Mediterranean sites that have been generally considered tectonically stable or affected by mild tectonics. These are found within a range of elevations of 2–10 m above modern mean sea level. Four sites are characterized by two separate sea-level stands, which suggest a two-step sea-level highstand during MIS 5e. Comparing field data with numerical modeling we show that (i) GIA is an important contributor to the spatial and temporal variability of the sea-level highstand during MIS 5e, (ii) the isostatic imbalance from the melting of the MIS 6 ice sheet can produce a >2.0 m sea-level highstand, and (iii) a two-step melting phase for the Greenland and Antarctic Ice Sheets reduces the differences between observations and predictions. Our results show that assumptions of tectonic stability on the basis of the MIS 5e records carry intrinsically large uncertainties, stemming either from uncertainties in field data and GIA models. The latter are propagated to either Holocene or Pleistocene sea-level reconstructions if tectonic rates are considered linear through time.

@en

We evaluate modelled Greenland ice sheet (GrIS) near-surface climate, surface energy balance (SEB) and surface mass balance (SMB) from the updated regional climate model RACMO2 (1958-2016). The new model version, referred to as RACMO2.3p2, incorporates updated glacier outlines, topography and ice albedo fields. Parameters in the cloud scheme governing the conversion of cloud condensate into precipitation have been tuned to correct inland snowfall underestimation: snow properties are modified to reduce drifting snow and melt production in the ice sheet percolation zone. The ice albedo prescribed in the updated model is lower at the ice sheet margins, increasing ice melt locally. RACMO2.3p2 shows good agreement compared to in situ meteorological data and point SEB/SMB measurements, and better resolves the spatial patterns and temporal variability of SMB compared with the previous model version, notably in the north-east, south-east and along the K-transect in south-western Greenland. This new model version provides updated, high-resolution gridded fields of the GrIS present-day climate and SMB, and will be used for projections of the GrIS climate and SMB in response to a future climate scenario in a forthcoming study.

@en

Ice flow forced by gravity is governed by the full Stokes (FS) equations, which are computationally expensive to solve due to the nonlinearity introduced by the rheology. Therefore, approximations to the FS equations are commonly used, especially when modeling a marine ice sheet (ice sheet, ice shelf, and/or ice stream) for 103 years or longer. The shallow ice approximation (SIA) and shallow shelf approximation (SSA) are commonly used but are accurate only for certain parts of an ice sheet. Here, we report a novel way of iteratively coupling FS and SSA that has been implemented in Elmer/Ice and applied to conceptual marine ice sheets. The FS-SSA coupling appears to be very accurate; the relative error in velocity compared to FS is below 0.5 % for diagnostic runs and below 5 % for prognostic runs. Results for grounding line dynamics obtained with the FS-SSA coupling are similar to those obtained from an FS model in an experiment with a periodical temperature forcing over 3000 years that induces grounding line advance and retreat. The rapid convergence of the FS-SSA coupling shows a large potential for reducing computation time, such that modeling a marine ice sheet for thousands of years should become feasible in the near future. Despite inefficient matrix assembly in the current implementation, computation time is reduced by 32 %, when the coupling is applied to a 3-D ice shelf..

@en

Isostatic response of the bedrock, or glacial isostatic adjustment (GIA) in included in most ice sheet models. This is important because the surface elevation determines the mass balance and thereby implicitly also the strength of the mass balance feedback where higher surface elevation yields lower temperatures implying less melt and vice versa. Usually a single relaxation time or a set of relaxation times is used to model the response everywhere on Earth or at least for an entire ice sheet. In reality the viscosity in the Earth’s mantle, and hence the relaxation time experienced by the ice, varies with location. Seismic studies indicate that several regions that were covered by ice during the last glacial cycle are underlain by mantle in which viscosity varies with orders of magnitude, such as Antarctica and North America. The question is whether such a variation of viscosity influences ice evolution. Several GIA models exist that can deal with 3D viscosity, but their large computation times make it nearly impossible to couple them to ice sheet models. Here we use the ANICE ice-sheet model (de Boer et al. 2013) with a simple bedrock-relaxation model in which a different relaxation time is used for separate regions. A temperature anomaly is applied to grow a schematic ice sheet on a flat earth, with other forcing mechanisms neglected. It is shown that in locations with a fast relaxation time of 300 years the equilibrium ice sheet is significantly thinner and narrower but also ice thickness in neighbouring regions (with the more standard relaxation time of 3000 years) is affected. @en

The glacial history of Antarctica during the most recent Milankovitch cycles is poorly constrained relative to the Northern Hemisphere. As a consequence, the contribution of mass changes in the Antarctic ice sheet to global sea-level change and the prediction of its future evolution remain uncertain. The process of Glacial Isostatic Adjustment (GIA) represents the ongoing response of the solid Earth to the Late-Pleistocene deglaciation and, therefore, provides information about Antarctic glacial history. Moreover, insufficient knowledge of GIA hampers the determination of present-day changes in the Antarctic mass balance through satellite gravity measurements. Previous studies have laid the theoretical foundation for distinguishing between signals of ongoing GIA and contemporary ice mass change through the combination of satellite gravimetry and satellite altimetry. This distinction is made possible by the fact the GIA-induced changes (involving relatively dense rock) will produce a different combination of topography and gravity change than those produced by variations in ice or firn thickness (due to the lower density of these materials); however, no conclusive results have been produced to date. Here we show that, by combining laser altimetry and gravity data from the ICESat and GRACE satellite missions over the period March 2003-March 2008, the GIA contribution can indeed be isolated. The inferred GIA signal over the Antarctic continent, which represents the first result derived from direct observations by satellite techniques, strongly supports Late-Pleistocene ice models derived from glacio-geologic studies. The GIA impact on GRACE-derived estimates of mass balance is found to be 100 ± 67 Gt/yr.

@en