ALBERT

Description: We use the BISICLES adaptive mesh ice sheet model to carry out one, two, and three century simulations of the fast-flowing ice streams of the West Antarctic Ice Sheet. Each of the simulations begins with a geometry and velocity close to present day observations, and evolves according to variation in meteoric ice accumulation, ice shelf melting, and mesh resolution. Future changes in accumulation and melt rates range from no change, through anomalies computed by atmosphere and ocean models driven by the E1 and A1B emissions scenarios, to spatially uniform melt rates anomalies that remove most of the ice shelves over a few centuries. We find that variation in the resulting ice dynamics is dominated by the choice of initial conditions, ice shelf melt rate and mesh resolution, although ice accumulation affects the net change in volume above flotation to a similar degree. Given sufficient melt rates, we compute grounding line retreat over hundreds of kilometers in every major ice stream, but the ocean models do not predict such melt rates outside of the Amundsen Sea Embayment until after 2100. Sensitivity to mesh resolution is spurious, and we find that sub-kilometer resolution is needed along most regions of the grounding line to avoid systematic under-estimates of the retreat rate, although resolution requirements are more stringent in some regions – for example the Amundsen Sea Embayment – than others – such as the Möller and Institute ice streams.

Description: The largest uncertainty in projections of future sea-level change still results from the potentially changing dynamical ice discharge from Antarctica. While ice discharge can alter through a number of processes, basal ice-shelf melting induced by a warming ocean has been identified as a major if not the major cause for possible additional ice flow across the grounding line. Here we derive dynamic ice-sheet response functions for basal ice-shelf melting using experiments carried out within the Sea-level Response to Ice Sheet Evolution (SeaRISE) intercomparison project with five different Antarctic ice-sheet models. As used here these response functions provide separate contributions for four different Antarctic drainage regions. Under the assumptions of linear-response theory we project future ice-discharge for each model, each region and each of the four Representative Concentration Pathways (RCP) using oceanic temperatures from 19 comprehensive climate models of the Coupled Model Intercomparison Project, CMIP-5, and two ocean models from the EU-project Ice2Sea. Uncertainty in the climatic forcing, the oceanic response and the ice-model differences is combined into an uncertainty range of future Antarctic ice-discharge induced from basal ice-shelf melt. The additional ice-loss (Table 6) is clearly scenario-dependent and results in a median of 0.07 m (66%-range: 0.04–0.10 m; 90%-range: −0.01–0.26 m) of global sea-level equivalent for the low-emission RCP-2.6 scenario and yields 0.1 m (66%-range: 0.06–0.14 m; 90%-range: −0.01–0.45 m) for the strongest RCP-8.5. If only models with an explicit representation of ice-shelves are taken into account the scenario dependence remains and the values change to: 0.05 m (66%-range: 0.03–0.08 m) for RCP-2.6 and 0.07 m (66%-range: 0.04–0.11 m) for RCP-8.5. These results were obtained using a time delay between the surface warming signal and the subsurface oceanic warming as observed in the CMIP-5 models. Without this time delay the ranges for all ice-models changes to 0.10 m (66%-range: 0.07–0.12 m; 90%-range: 0.01–0.28 m) for RCP-2.6 and 0.15 m (66%-range: 0.10–0.21 m; 90%-range: 0.02–0.53 m) for RCP-8.5. All probability distributions as provided in Fig. 10 are highly skewed towards high values.

Description: We use the BISICLES adaptive mesh ice sheet model to carry out one, two, and three century simulations of the fast-flowing ice streams of the West Antarctic Ice Sheet, deploying sub-kilometer resolution around the grounding line since coarser resolution results in substantial underestimation of the response. Each of the simulations begins with a geometry and velocity close to present-day observations, and evolves according to variation in meteoric ice accumulation rates and oceanic ice shelf melt rates. Future changes in accumulation and melt rates range from no change, through anomalies computed by atmosphere and ocean models driven by the E1 and A1B emissions scenarios, to spatially uniform melt rate anomalies that remove most of the ice shelves over a few centuries. We find that variation in the resulting ice dynamics is dominated by the choice of initial conditions and ice shelf melt rate and mesh resolution, although ice accumulation affects the net change in volume above flotation to a similar degree. Given sufficient melt rates, we compute grounding line retreat over hundreds of kilometers in every major ice stream, but the ocean models do not predict such melt rates outside of the Amundsen Sea Embayment until after 2100. Within the Amundsen Sea Embayment the largest single source of variability is the onset of sustained retreat in Thwaites Glacier, which can triple the rate of eustatic sea level rise.

Description: We analyzed hydrographic data from the northwestern Weddell Sea continental shelf of the three austral winters 1989, 1997, and 2006 and two summers following the last winter cruise. During summer a thermal front exists at ~64° S separating cold southern waters from warm northern waters that have similar characteristics as the deep waters of the central basin of the Bransfield Strait. In winter, the whole continental shelf exhibits southern characteristics with high Neon (Ne) concentrations, indicating a significant input of glacial melt water. The comparison of the winter data from the shallow shelf off the tip of the Antarctic Peninsula, spanning a period of 17 yr, shows a salinity decrease of 0.09 for the whole water column, which has a residence time of

Description: We analysed hydrographic data from the northwestern Weddell Sea continental shelf of three austral winters (1989, 1997 and 2006) and two summers following the last winter cruise. During summer a thermal front exists at ~64° S separating cold southern waters from warm northern waters that have similar characteristics as the deep waters of the central basin of the Bransfield Strait. In winter, the whole continental shelf exhibits southern characteristics with high Neon (Ne) concentrations, indicating a significant input of glacial melt water. The comparison of the winter data at the tip of the Antarctic Peninsula, spanning a period of 17 years, shows a salinity decrease of 0.09 for the whole water column. We interpret this freshening as a reduction in salt input to the water masses being advected northward on the western Weddell Sea continental shelf. Possible causes for the reduced winter salinification are a southward retreat of the summer sea ice edge together with more precipitation in this sector. However, the latter might have happened in conjunction with an increase in ice shelf mass loss, counteracting an enhanced salt input due to sea ice formation in coastal areas formerly occupied by Larsen A and B ice shelves.

Description: Sub-ice shelf circulation and freezing/melting rates in ocean general circulation models depend critically on an accurate and consistent representation of cavity geometry. Existing global or pan-Antarctic topography data sets have turned out to contain various inconsistencies and inaccuracies. The goal of this work is to compile independent regional surveys and maps into a global data set. We use the S-2004 global 1-min bathymetry as the backbone and add an improved version of the BEDMAP topography (ALBMAP bedrock topography) for an area that roughly coincides with the Antarctic continental shelf. The position of the merging line is individually chosen in different sectors in order to capture the best of both data sets. High-resolution gridded data for ice shelf topography and cavity geometry of the Amery, Fimbul, Filchner-Ronne, Larsen C and George VI Ice Shelves, and for Pine Island Glacier are carefully merged into the ambient ice and ocean topographies. Multibeam survey data for bathymetry in the former Larsen B cavity and the southeastern Bellingshausen Sea have been obtained from the data centers of Alfred Wegener Institute (AWI), British Antarctic Survey (BAS) and Lamont-Doherty Earth Observatory (LDEO), gridded, and blended into the existing bathymetry map. The resulting global 1-min Refined Topography data set (RTopo-1) contains self-consistent maps for upper and lower ice surface heights, bedrock topography, and surface type (open ocean, grounded ice, floating ice, bare land surface). The data set is available in NetCDF format from the PANGAEA database at doi:10.1594/pangaea.741917.

Description: The largest uncertainty in projections of future sea-level change results from the potentially changing dynamical ice discharge from Antarctica. Basal ice-shelf melting induced by a warming ocean has been identified as a major cause for additional ice flow across the grounding line. Here we derive dynamic ice-sheet response functions for basal ice-shelf melting for four different Antarctic drainage regions using experiments from the Sea-level Response to Ice Sheet Evolution (SeaRISE) intercomparison project with five different Antarctic ice-sheet models. Under the assumptions of linear-response theory we project future ice-discharge for each model, each region and each of the four Representative Concentration Pathways (RCP) using oceanic temperatures from 19 comprehensive climate models of the Coupled Model Intercomparison Project, CMIP-5, and two ocean models from the EU-project Ice2Sea. The uncertainty in the climatic forcing, the oceanic response and the ice-model response is combined into an uncertainty range of future Antarctic ice-discharge induced from basal ice-shelf melt. The uncertainty range we derived for the Antarctic contribution to global sea-level rise from 1992 to 2011 is in full agreement with the observed contribution for this period if we use the three ice-sheet models with an explicit representation of ice-shelf dynamics and account for the time delayed warming of the oceanic subsurface compared with the surface air temperature. The median of the additional ice-loss for the 21st century (Table 6) is 0.07 m (66%-range: 0.02–0.14 m; 90%-range: 0.0–0.23 m) of global sea-level equivalent for the low-emission RCP-2.6 scenario and 0.09 m (66%-range: 0.04–0.21 m; 90%-range: 0.01–0.37 m) for the strongest RCP-8.5 if models with explicit ice-shelf representation are applied. These results were obtained using a time delay between the surface warming signal and the subsurface oceanic warming as observed in the CMIP-5 models. Without this time delay the values increase to 0.09 m (66%-range: 0.04–0.17 m; 90%-range: 0.02–0.25 m) for RCP-2.6 and 0.15 m (66%-range: 0.07–0.28 m; 90%-range: 0.04–0.43 m) for RCP-8.5. Our results are scenario dependent which is most visible in the upper percentiles of the distribution, i.e. highest contributions to sea level rise. All probability distributions, as provided in Fig. 12, are highly skewed towards high values. The applied ice-sheet models are coarse-resolution with limitations in the representation of grounding-line motion. However, we find the main uncertainty to be introduced by the external forcing to the ice-sheets, i.e. the climatic and oceanic uncertainty dominate. The scaling coefficients for the four different drainage basins provide valuable information for further assessments of future Antarctic ice discharge.

Description: The largest uncertainty in projections of future sea-level change results from the potentially changing dynamical ice discharge from Antarctica. Basal ice-shelf melting induced by a warming ocean has been identified as a major cause for additional ice flow across the grounding line. Here we attempt to estimate the uncertainty range of future ice discharge from Antarctica by combining uncertainty in the climatic forcing, the oceanic response and the ice-sheet model response. The uncertainty in the global mean temperature increase is obtained from historically constrained emulations with the MAGICC-6.0 (Model for the Assessment of Greenhouse gas Induced Climate Change) model. The oceanic forcing is derived from scaling of the subsurface with the atmospheric warming from 19 comprehensive climate models of the Coupled Model Intercomparison Project (CMIP-5) and two ocean models from the EU-project Ice2Sea. The dynamic ice-sheet response is derived from linear response functions for basal ice-shelf melting for four different Antarctic drainage regions using experiments from the Sea-level Response to Ice Sheet Evolution (SeaRISE) intercomparison project with five different Antarctic ice-sheet models. The resulting uncertainty range for the historic Antarctic contribution to global sea-level rise from 1992 to 2011 agrees with the observed contribution for this period if we use the three ice-sheet models with an explicit representation of ice-shelf dynamics and account for the time-delayed warming of the oceanic subsurface compared to the surface air temperature. The median of the additional ice loss for the 21st century is computed to 0.07 m (66% range: 0.02–0.14 m; 90% range: 0.0–0.23 m) of global sea-level equivalent for the low-emission RCP-2.6 (Representative Concentration Pathway) scenario and 0.09 m (66% range: 0.04–0.21 m; 90% range: 0.01–0.37 m) for the strongest RCP-8.5. Assuming no time delay between the atmospheric warming and the oceanic subsurface, these values increase to 0.09 m (66% range: 0.04–0.17 m; 90% range: 0.02–0.25 m) for RCP-2.6 and 0.15 m (66% range: 0.07–0.28 m; 90% range: 0.04–0.43 m) for RCP-8.5. All probability distributions are highly skewed towards high values. The applied ice-sheet models are coarse resolution with limitations in the representation of grounding-line motion. Within the constraints of the applied methods, the uncertainty induced from different ice-sheet models is smaller than that induced by the external forcing to the ice sheets.

Description: Sub-ice shelf circulation and freezing/melting rates in ocean general circulation models depend critically on an accurate and consistent representation of cavity geometry. Existing global or pan-Antarctic data sets have turned out to contain various inconsistencies and inaccuracies. The goal of this work is to compile independent regional fields into a global data set. We use the S-2004 global 1-min bathymetry as the backbone and add an improved version of the BEDMAP topography (ALBMAP bedrock topography) for an area that roughly coincides with the Antarctic continental shelf. The position of the merging line is individually chosen in different sectors in order to get the best out of each data set. High-resolution gridded data for ice shelf topography and cavity geometry of the Amery, Fimbul, Filchner-Ronne, Larsen C and George VI Ice Shelves, and for Pine Island Glacier are carefully merged into the ambient ice and ocean topographies. Multibeam survey data for bathymetry in the former Larsen B cavity and the southeastern Bellingshausen Sea have been obtained from the data centers of Alfred Wegener Institute (AWI), British Antarctic Survey (BAS) and Lamont-Doherty Earth Observatory (LDEO), gridded, and blended into the existing bathymetry map. The resulting global 1-min topography data set (RTopo-1) contains maps for upper and lower ice surface heights, bedrock topography, and consistent masks for open ocean, grounded ice, floating ice, and bare land surface. The data set is available in NetCDF format from the PANGAEA database at doi:10.1594/pangaea.741917.

Description: The Finite-Element Sea-Ice Model, used as a component of the Finite-Element Sea ice Ocean Model, is presented. Version 2 includes the elastic-viscous-plastic (EVP) and viscous-plastic (VP) solvers and employs a flux corrected transport algorithm to advect the ice and snow mean thicknesses and concentration. The EVP part also includes a modified approach proposed recently by Bouillon et al., which is characterized by an improved stability compared to the standard EVP approach. The model is formulated on unstructured triangular meshes. It assumes a collocated placement of ice velocities, mean thicknesses and concentration at mesh vertices, and relies on a piecewise-linear (P1) continuous elements. Simple tests for the modified EVP and VP solvers are presented to show that they may produce very close results provided the number of iterations is sufficiently high.