Mass balance and acceleration of the Greenland ice sheet

Within the last decade, the Greenland ice sheet (GrIS) and its surroundings have experienced record high surface temperatures (Mote, 2007, doi:10.1029/2007GL031976; Box et al., 2010), ice sheet melt extent (Fettweis et al., 2011, doi:10.5194/tc-5-359-2011) and record-low summer sea-ice extent (Nghiem et al., 2007, doi:10.1029/2007GL031138). Using three independent data sets, we derive, for the first time, consistent ice-mass trends and temporal variations within seven major drainage basins from gravity fields from the Gravity Recovery and Climate Experiment (GRACE; Tapley et al., 2004, doi:10.1029/2004GL019920), surface-ice velocities from Inteferometric Synthetic Aperture Radar (InSAR; Rignot and Kanagaratnam, 2006, doi:10.1126/science.1121381) together with output of the regional atmospheric climate modelling (RACMO2/ GR; Ettema et al., 2009, doi:10.1029/2009GL038110), and surface-elevation changes from the Ice, cloud and land elevation satellite (ICESat; Sorensen et al., 2011, doi:10.5194/tc-5-173-2011). We show that changing ice discharge (D), surface melting and subsequent run-off (M/R) and precipitation (P) all contribute, in a complex and regionally variable interplay, to the increasingly negative mass balance of the GrIS observed within the last decade. Interannual variability in P along the northwest and west coasts of the GrIS largely explains the apparent regional mass loss increase during 2002-2010, and obscures increasing M/R and D since the 1990s. In winter 2002/2003 and 2008/2009, accumulation anomalies in the east and southeast temporarily outweighed the losses by M/R and D that prevailed during 2003-2008, and after summer 2010. Overall, for all basins of the GrIS, the decadal variability of anomalies in P, M/R and D between 1958 and 2010 (w.r.t. 1961-1990) was significantly exceeded by the regional trends observed during the GRACE period (2002-2011).

Reconstructed ice-sheet limit in the Kara and Barents Sea area at about 50,000 years ago (Markhida Line)

Glacial landforms in northern Russia, from the Timan Ridge in the west to the east of the Urals, have been mapped by aerial photographs and satellite images supported by field observations. An east-west trending belt of fresh hummock-and-lake glaciokarst landscapes has been traced to the north of 67°N. The southern boundary of these landscapes is called the Markhida Line, which is interpreted as a nearly synchronous limit of the last ice sheet that affected this region. The hummocky landscapes are subdivided into three types according to the stage of postglacial modification: Markhida, Harbei and Halmer. The Halmer landscape on the Uralian piedmont in the east is the freshest, whereas the westernmost Markhida landscape is more eroded. The west- east gradient in morphology is considered to be a result of the time-transgressive melting of stagnant glacier ice and of the underlying permafrost. The pattern of ice-pushed ridges and other directional features reflects a dominant ice flow direction from the Kara Sea shelf. Traces of ice movement from the central Barents Sea are only discernible in the Pechora River left bank area west of 50°E. In the Polar Urals the horseshoe-shaped end moraines at altitudes of up to 560 m a.s.l. reflect ice movement up-valley from the Kara Ice Sheet, indicating the absence of a contemporaneous ice dome in the mountains. The Markhida moraines, superimposed onto the Eemian strata, represent the maximum ice sheet extent in the western part of the Pechora Basin during the Weichselian. The Markhida Line truncates the huge arcs of the Laya-Adzva and Rogovaya ice-pushed ridges protruding to the south. The latter moraines therefore reflect an older ice advance, probably also of Weichselian age. Still farther south, fluvially dissected morainic plateaus without lakes are of pre-Eemian age, because they plunge northwards under marine Eemian sediments. Shorelines of the large ice-dammed Lake Komi, identified between 90 and 110 m a.s.l. in the areas south of the Markhida Line, are radiocarbon dated to be older than 45 ka. The shorelines, incised into the Laya-Adzva moraines, morphologically interfinger with the Markhida moraines, indicating that the last ice advance onto the Russian mainland reached the Markhida Line during the Middle or Early Weichselian, before 45 ka ago.

Greenland and Antarctic ice sheet topography, cavity geometry, and global bathymetry (RTopo-2), links to NetCDF files

The ocean plays an important role in modulating the mass balance of the polar ice sheets by interacting with the ice shelves in Antarctica and with the marine-terminating outlet glaciers in Greenland. Given that the flux of warm water onto the continental shelf and into the sub-ice cavities is steered by complex bathymetry, a detailed topography data set is an essential ingredient for models that address ice-ocean interaction. We followed the spirit of the global RTopo-1 data set and compiled consistent maps of global ocean bathymetry, upper and lower ice surface topographies and global surface height on a spherical grid with now 30-arc seconds resolution. We used the General Bathymetric Chart of the Oceans (GEBCO, 2014) as the backbone and added the International Bathymetric Chart of the Arctic Ocean version 3 (IBCAOv3) and the Interna- tional Bathymetric Chart of the Southern Ocean (IBCSO) version 1. While RTopo-1 primarily aimed at a good and consistent representation of the Antarctic ice sheet, ice shelves and sub-ice cavities, RTopo-2 now also contains ice topographies of the Greenland ice sheet and outlet glaciers. In particular, we aimed at a good representation of the fjord and shelf bathymetry sur- rounding the Greenland continent. We corrected data from earlier gridded products in the areas of Petermann Glacier, Hagen Bræ and Sermilik Fjord assuming that sub-ice and fjord bathymetries roughly follow plausible Last Glacial Maximum ice flow patterns. For the continental shelf off northeast Greenland and the floating ice tongue of Nioghalvfjerdsfjorden Glacier at about 79°N, we incorporated a high-resolution digital bathymetry model considering original multibeam survey data for the region. Radar data for surface topographies of the floating ice tongues of Nioghalvfjerdsfjorden Glacier and Zachariæ Isstrøm have been obtained from the data centers of Technical University of Denmark (DTU), Operation Icebridge (NASA/NSF) and Alfred Wegener Institute (AWI). For the Antarctic ice sheet/ice shelves, RTopo-2 largely relies on the Bedmap-2 product but applies corrections for the geometry of Getz, Abbot and Fimbul ice shelf cavities.

Sulfate and nitrate firn concentrations on the Greenland ice sheet

Intercomparison of three new chemical ice core records from northern Greenland (covering the time span from approximately 1500 A.D. to present) with previously published records for southern and central Greenland reveals a uniform timing of anthropogenic changes in sulfate and nitrate firn concentrations over the entire ice sheet. The anthropogenic sulfate increase started around 1890, was interrupted by a transient decrease in the 1930s, and has resumed a major increase since 1950. Since the late 1970s though, a significant 30% decline in Greenland sulfate firn levels can be documented. The maximum anthropogenic increase in northern Greenland sulfate firn concentrations (up to 200-230 ppb) is 2-3 times larger than in southern and central Greenland. Nitrate records show an essentially steady increase since 1950 and, documented for the first time, a slight reduction during most recent years. Maximum nitrate firn levels of 100-130 ppb exceed the preindustrial background by 100% all over the Greenland ice sheet. Comparison with anthropogenic SO2 and NO x emission records indicates that the major increase in sulfate firn concentrations since 1950 can be attributed to Eurasian sources, while firn levels during the first half of this century appear to be dominated by North American emissions. A stronger North American source contribution is indicated over the entire 20th century in the case of nitrate. Applying a macroscopic deposition model separate time series for wet and dry deposition were derived which revealed a close correspondence of wet deposited sulfate with the timing of U.S. emissions, while the temporal evolution of Eurasian emissions is mainly reflected in the dry sulfate deposition record. During this century wet sulfate deposition increased by a factor of two while the total dry sulfate deposition flux increased by more than 500%. Wet and dry nitrate deposition both increased by 100% during the same period.

Antarctic ice sheet topography, cavity geometry, and global bathymetry (RTopo 1.0.5-beta)

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-minute bathymetry as the backbone and add an improved version of the BEDMAP topography for an area that roughly coincides with the Antarctic continental shelf. Locations of the merging line have been carefully adjusted in order to get the best out of each data set. High-resolution gridded data for upper and lower ice surface topography and cavity geometry of the Amery, Fimbul, Filchner-Ronne, Larsen C and George VI Ice Shelves, and for Pine Island Glacier have been 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 again carefully merged into the existing bathymetry map. The global 1-minute dataset (RTopo-1 Version 1.0.5) has been split into two netCDF files. The first contains digital maps for global bedrock topography, ice bottom topography, and surface elevation. The second contains the auxiliary maps for data sources and the surface type mask. A regional subset that covers all variables for the region south of 50 deg S is also available in netCDF format. Datasets for the locations of grounding and coast lines are provided in ASCII format.

(Table 1) Mass loss of the Greenland Ice Sheet

Using satellite radar interferometry observations of Greenland, we detected widespread glacier acceleration below 66° north between 1996 and 2000, which rapidly expanded to 70° north in 2005. Accelerated ice discharge in the west and particularly in the east doubled the ice sheet mass deficit in the last decade from 90 to 220 cubic kilometers per year. As more glaciers accelerate farther north, the contribution of Greenland to sea-level rise will continue to increase.