Search

The projected decay of the Antarctic Ice Sheet (AIS) over the coming centuries will lead to uplift of the Earth's surface due to Glacial Isostatic Adjustment (GIA). GIA slows down grounding line migration and therefore has a stabilizing effect on the ice sheet evolution. GIA acts on timescales of decades to centennial depending on the magnitude of the mantle viscosity. The mantle viscosity is several orders of magnitude higher in East Antarctica than in West Antarctica and varies with one order of magnitude within West Antarctica. Studies of the AIS evolution over the last glacial cycle have shown that including lateral variations of the Earth's mantle viscosity can lead to 1.5-kilometer thicker ice in West Antarctica at present day. However, current projections do not include GIA, or they use a laterally homogeneous GIA model. One study applied a uniform high mantle viscosity under East Antarctica and a uniform low mantle viscosity under West Antarctica and showed that, on longer timescales of hundreds of years, mass loss projections of Antarctica may be underestimated because spatially uniform GIA models overestimate the stabilizing effect of GIA across East Antarctica. We developed a coupled GIA - ice-sheet model using the ice-sheet model IMAU-ICE, and a 3D GIA finite element model that includes lateral mantle viscosity variations, and a seismic model to determine the patterns of the viscosity. The results of projections for two IPCC scenarios show that including lateral variations in the Earth's mantle viscosity leads to local ice thickness differences of up to 600 meters in West Antarctica 2300. The results underline and quantify the importance of including this local feedback effect in ice-sheet models when projecting the long-term sea level contribution from Antarctica.

Both the Greenland and Antarctic ice sheets are expected to experience substantial mass loss in the case of unmitigated anthropogenic climate change. The exact rate of future mass loss under high warming scenarios remains uncertain, depending strongly on physical quantities that are difficult to constrain from observations, such as basal sliding and sub-shelf melt. We apply a novel model initialisation protocol, that combines elements from existing approaches such as the equilibrium spin-up, basal inversion, and palaeo spin-up, to models of both the Greenland and Antarctic ice sheets. We show the results in term of sea-level projections including the uncertainties, under different warming scenarios, following the ISMIP6 protocol.This abstract is a companion to “On the initialisation of ice sheet models: equilibrium assumptions, thermal memory, and present-day states” by Bernales et al. We hope that, if both abstracts are lucky enough to be accepted, the conveners can program the two talks in sequence.

A significant portion of the spread in future projections of ice sheet volume changes is attributed to uncertainties in their present-day state, and the way this state is represented in ice-sheet models. The scientific literature already contains a variety of classic initialisation approaches used by modelling groups around the globe, each with its own advantages and limitations. We propose a generalised protocol that allows for the quantification of the impact of individual initialisation choices, such as steady-state assumptions, the inclusion of internal paleoclimatic thermal signals, sea level and glacial isostatic effects, and calibration methods. We then apply this protocol to an ensemble of multi-millennia model spin-ups of the present-day Greenland and Antarctic ice sheets and show the importance of the choices made during initialisation.[This abstract is a companion to “Sensitivity of future projections of ice sheet retreat to initial conditions” by Berends et al. We hope that, if both abstracts are lucky enough to be accepted, the conveners can program the two talks in sequence.]

In geodynamic studies, most Finite-Element (FE) models in the commercial FE software Abaqus use elastic foundations at internal boundaries. This method works well for incompressible and so-called material-compressible material parameters but it is unclear if it works sufficiently well for implementing compressibility, especially in a 3D spherical model. The latter is of importance in investigations of glacial isostatic adjustment (GIA). A possible alternative method is based on a combination of explicit gravity loading with non-linear geometry (NLGEOM parameter in Abaqus) (Hampel et al., 2019). This method would remove the need to make a stress transformation to get the correct GIA stresses, and automatically accounts for the change in internal buoyancy forces that arises when allowing for compression, according to the Abaqus Documentation. We compared the method for (in)compressible flat (~half-space) FE models with existing numerical half-space and spherical (in)compressible codes and tested the applicability of this method in a spherical FE model. We confirm that this method works for multi-layer incompressible flat FE models. We furthermore notice that horizontal displacement rates of incompressible flat FE models match those of spherical incompressible GIA models below the current GNSS (Global Navigation Satellite System) measurement accuracy of 0.2-0.3 mm/a, but only for ice sheets that are smaller than 450 km in extent. For compressible models, disagreements in the vertical displacement rates are found between the flat NLGEOM model and the compressible Normal Mode code ICEAGE (Kaufmann, 2004). An extension of the NLGEOM-gravity method to a spherical FE model, where gravity must be implemented in the form of body forces combined with initial stress, leads to a divergence of the solution when viscous behaviour is turned on. We thus conclude that the applicability of the NLGEOM method is so far limited to flat FE models, and in GIA investigations for flat models the applicability further depends on the size of the load (ice sheet, glacier). References: Hampel, A., Lüke, J., Krause, T., & Hetzel, R., 2019. Finite-element modelling of glacial isostatic ad-justment (GIA): Use of elastic foundations at material boundaries versus the geometrically non-linearformulation, Computers & geosciences, 122, 1–14. Kaufmann, G. (2004). Program Package ICEAGE, Version 2004. Manuscript. Institut für Geophysik der Universität Göttingen.

Changes in ice load over time deform the Earth’s crust and mantle. This effect, Glacial Isostatic Adjustment (GIA), induces vertical deformation of the bedrock of the Antarctic continent and impacts the grounding line position which is critical for the dynamical state of the Antarctic Ice Sheet (AIS). GIA introduces a negative feedback and stabilizes the ice sheet evolution, hence GIA modelling is important for transient studies. Most ice dynamic models use a two-layer flat Earth approach with a laterally homogenous relaxation time or a layered Earth approach with a laterally homogenous viscosity (1D) to compute the bedrock deformation. However, viscosity of the Earth’s interior varies laterally (3D) and radially with several orders of magnitude across the Antarctic continent. Here we present a new coupled 3D GIA – ice dynamic model which can run over hundred thousands of years with a resolution of 500 years. The method is applied using various 1D and 3D rheologies. Results show that the present-day ice volume is 3 % lower when using a 1D viscosity of 1021 Pa·s than using a 3D viscosity. However, local differences in grounding line position maybe up to a hundred kilometres around the Ronne and the Ross Ice Shelfs, and ice thickness differences are up to a kilometre for present day conditions when comparing 1D rheologies and 3D rheologies. The difference between the use of various 3D rheologies is significantly smaller. These results underline and quantify the importance of including local GIA feedback effects in ice dynamic models when simulating the Antarctic Ice Sheet evolution over the Last Glacial Cycle.

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.

Plastic pollution in the marine environment is an urgent global environmental challenge. Land-based plastics, emitted into the ocean through rivers, are believed to be the main source of marine plastic litter. According to the latest model-based estimates, most riverine plastics are emitted in Asia. However, the exact amount of global riverine plastic emission remains uncertain due to a severe lack of observation. Field-based studies are rare in numbers, focused on rivers in Europe and North America and used strongly varying data collection methods. We present a harmonized assessment of floating macroplastic transport from observations at 24 locations in rivers in seven countries in Europe and Asia. Visual counting and debris sampling were used to assess (1) magnitude of plastic transport, (2) the spatial distribution across the river width, and (3) the plastic polymer composition. Several waterways in Indonesia and Vietnam contain up to four orders of magnitude more plastic than waterways in Italy, France, and The Netherlands in terms of plastic items per hour. We present a first transcontinental overview of plastic transport, providing observational evidence that, for the sampled rivers, Asian rivers transport considerably more plastics towards the ocean. New insights are presented in the magnitude, composition, and spatiotemporal variation of riverine plastic debris. We emphasize the urgent need for more long-term monitoring efforts. Accurate data on riverine plastic debris are extremely important to improve global and local modeling approaches and to optimize prevention and collection strategies.

Over glacial-interglacial cycles, the evolution of an ice sheet is influenced by Glacial isostatic adjustment (GIA) via two negative feedback loops. Firstly, vertical bedrock deformation due to a changing ice load alters ice-sheet surface elevation. For example, an increasing ice load leads to a lower bedrock elevation that lowers ice-sheet surface elevation. This will increase surface melting of the ice sheet, following an increase of atmospheric temperature at lower elevations. Secondly, bedrock deformation will change the height of the grounding line of the ice sheet. For example, a lowering bedrock height following ice-sheet advance increases the melt due to ocean water that in turn leads to a retreat of the grounding line and a slow-down of ice-sheet advance. GIA is mainly determined by the viscosity of the interior of the solid Earth which is radially and laterally varying. Underneath the Antarctic ice sheet, there are relatively low viscosities in West Antarctica and higher viscosities in East Antarctica, in turn affecting the response time of the above mentioned feedbacks. However, most ice-dynamical models do not consider the lateral variations of the viscosity in the GIA feedback loops when simulating the evolution of the Antarctic ice sheet. The method developed by Gomez et al. (2018) includes the feedback between GIA and ice-sheet evolution and alternates between simulations of the two models where each simulation covers the full time period. We presents a different method to couple ANICE, a 3-D ice-sheet model, to a 3-D GIA finite element model. In this method the model computations alternates between the ice-sheet and GIA model until convergence of the result occurs at each timestep. We simulate the evolution of the Antarctic ice sheet from 120 000 years ago to the present. The results of the coupled simulation will be discussed and compared to results of the uncoupled ice-sheet model (using an ELRA GIA model) and the method developed by Gomez et al. (2018).