Simulated temperature changes as a result of Antarctic ice sheet growth and CO2 reduction

During the Middle Miocene climate transition about 14 million years ago, the Antarctic ice sheet expanded to near-modern volume. Surprisingly, this ice sheet growth was accompanied by a warming in the surface waters of the Southern Ocean, whereas a slight deep-water temperature increase was delayed by more than 200 thousand years. Here we use a coupled atmosphere-ocean model to assess the relative effects of changes in atmospheric CO2 concentration and ice sheet growth on regional and global temperatures. In the simulations, changes in the wind field associated with the growth of the ice sheet induce changes in ocean circulation, deep-water formation and sea-ice cover that result in sea surface warming and deep-water cooling in large swaths of the Atlantic and Indian ocean sectors of the Southern Ocean. We interpret these changes as the dominant ocean surface response to a 100-thousand-year phase of massive ice growth in Antarctica. A rise in global annual mean temperatures is also seen in response to increased Antarctic ice surface elevation. In contrast, the longer-term surface and deep-water temperature trends are dominated by changes in atmospheric CO2 concentration. We therefore conclude that the climatic and oceanographic impacts of the Miocene expansion of the Antarctic ice sheet are governed by a complex interplay between wind field, ocean circulation and the sea-ice system.

Lagrangian back-tracing results for Kohnen station; output of time dependent numerical experiment with a 3D nested Antarctic ice sheet model

A nested ice flow model was developed for eastern Dronning Maud Land to assist with the dating and interpretation of the EDML deep ice core. The model consists of a high-resolution higher-order ice dynamic flow model that was nested into a comprehensive 3-D thermomechanical model of the whole Antarctic ice sheet. As the drill site is on a flank position the calculations specifically take into account the effects of horizontal advection as deeper ice in the core originated from higher inland. First the regional velocity field and ice sheet geometry is obtained from a forward experiment over the last 8 glacial cycles. The result is subsequently employed in a Lagrangian backtracing algorithm to provide particle paths back to their time and place of deposition. The procedure directly yields the depth-age distribution, surface conditions at particle origin, and a suite of relevant parameters such as initial annual layer thickness. This paper discusses the method and the main results of the experiment, including the ice core chronology, the non-climatic corrections needed to extract the climatic part of the signal, and the thinning function. The focus is on the upper 89% of the ice core (appr. 170 kyears) as the dating below that is increasingly less robust owing to the unknown value of the geothermal heat flux. It is found that the temperature biases resulting from variations of surface elevation are up to half of the magnitude of the climatic changes themselves.

Planimetric ice-flow convergence and flow-orthonormal strain rate across the Antarctic Ice Sheet, with links to model result files

Fast-flowing ice streams discharge most of the ice from the interior of the Antarctic Ice Sheet coastward. Understanding how their tributary organisation is governed and evolves is essential for developing reliable models of the ice sheet's response to climate change. Despite much research on ice-stream mechanics, this problem is unsolved, because the complexity of flow within and across the tributary networks has hardly been interrogated. Here I present the first map of planimetric flow convergence across the ice sheet, calculated from satellite measurements of ice surface velocity, and use it to explore this complexity. The convergence map of Antarctica elucidates how ice-stream tributaries draw ice from the interior. It also reveals curvilinear zones of convergence along lateral shear margins of streaming, and abundant convergence ripples associated with nonlinear ice rheology and changes in bed topography and friction. Flow convergence on ice-stream tributaries and their feeding zones is markedly uneven, and interspersed with divergence at distances of the order of kilometres. For individual drainage basins as well as the ice sheet as a whole, the range of convergence and divergence decreases systematically with flow speed, implying that fast flow cannot converge or diverge as much as slow flow. I therefore deduce that flow in ice-stream networks is subject to mechanical regulation that limits flow-orthonormal strain rates. These properties and the gridded data of convergence and flow-orthonormal strain rate in this archive provide targets for ice- sheet simulations and motivate more research into the origin and dynamics of tributarization.

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.

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.

Corrected ICESat altimetry data, surface mass balance, and firn elevation change on Antarctic ice shelves

Accurate prediction of global sea-level rise requires that we understand the cause of recent, widespread and intensifying glacier acceleration along Antarctic ice-sheet coastal margins. Floating ice shelves buttress the flow of grounded tributary glaciers and their thickness and extent are particularly susceptible to changes in both climate and ocean forcing. Recent ice-shelf collapse led to retreat and acceleration of several glaciers on the Antarctic Peninsula. However, the extent and magnitude of ice-shelf thickness change, its causes and its link to glacier flow rate are so poorly understood that its influence on the future of the ice sheets cannot yet be predicted. Here we use satellite laser altimetry and modelling of the surface firn layer to reveal for the first time the circum-Antarctic pattern of ice-shelf thinning through increased basal melt. We deduce that this increased melt is the primary driver of Antarctic ice-sheet loss, through a reduction in buttressing of the adjacent ice sheet that has led to accelerated glacier flow. The highest thinning rates (~7 m/a) occur where warm water at depth can access thick ice shelves via submarine troughs crossing the continental shelf. Wind forcing could explain the dominant patterns of both basal melting and the surface melting and collapse of Antarctic ice shelves, through ocean upwelling in the Amundsen and Bellingshausen Seas and atmospheric warming on the Antarctic Peninsula. This implies that climate forcing through changing winds influences Antarctic Ice Sheet mass balance, and hence global sea-level, on annual to decadal timescales.

Water stable oxygen isotope records from six Antarctic ice core sites for the present and last interglacial periods

We compare the present and last interglacial periods as recorded in Antarctic water stable isotope records now available at various temporal resolutions from six East Antarctic ice cores: Vostok, Taylor Dome, EPICA Dome C (EDC), EPICA Dronning Maud Land (EDML), Dome Fuji and the recent TALDICE ice core from Talos Dome. We first review the different modern site characteristics in terms of ice flow, meteorological conditions, precipitation intermittency and moisture origin, as depicted by meteorological data, atmospheric reanalyses and Lagrangian moisture source diagnostics. These different factors can indeed alter the relationships between temperature and water stable isotopes. Using five records with sufficient resolution on the EDC3 age scale, common features are quantified through principal component analyses. Consistent with instrumental records and atmospheric model results, the ice core data depict rather coherent and homogenous patterns in East Antarctica during the last two interglacials. Across the East Antarctic plateau, regional differences, with respect to the common East Antarctic signal, appear to have similar patterns during the current and last interglacials. We identify two abrupt shifts in isotopic records during the glacial inception at TALDICE and EDML, likely caused by regional sea ice expansion. These regional differences are discussed in terms of moisture origin and in terms of past changes in local elevation histories, which are compared to ice sheet model results. Our results suggest that elevation changes may contribute significantly to inter-site differences. These elevation changes may be underestimated by current ice sheet models

(Appendix A) Geochemical composition of snow and eolian mineral dust from Berkner Island ice sheet, Antarctica

The Sr and Nd isotopic composition of dust extracted from recent snow layers at the top of Berkner Island ice sheet (located within the Filchner-Ronne Ice Shelf at the southern end of the Weddell Sea) enables us, for the first time, to document dust provenance in Antarctica outside the East Antarctic Plateau (EAP) where all previous studies based on isotopic fingerprinting were carried out. Berkner dust displays an overall crust-like isotopic signature, characterized by more radiogenic 87Sr/86Sr and much less radiogenic 143Nd/144Nd compared to dust deposited on the EAP during glacial periods. Differences with EAP interglacial dust are not as marked but still significant, indicating that present-day Berkner dust provenance is distinct, at least to some extent, from that of the dust reaching the EAP. The fourteen snow-pit sub-seasonal samples that were obtained span a two-year period (2002-2003) and their dust Sr and Nd isotopic composition reveals that multiple sources are at play over a yearly time period. Southern South America, Patagonia in particular, likely accounts for part of the observed spring/summer dust deposition maxima, when isotopic composition is shifted towards 'younger' isotopic signatures. In the spring, possible additional inputs from Australian sources would also be supported by the data. Most of the year, however, the measured isotopic signatures would be best explained by a sustained background supply from putative local sources in East Antarctica, which carry old-crust-like isotopic fingerprints. Whether the restricted East Antarctic ice-free areas produce sufficient eolian material has yet to be substantiated however. The fact that large (> 5 µm) particles represent a significant fraction of the samples throughout the entire time-series supports scenarios that involve contributions from proximal sources, either in Patagonia and/or Antarctica (possibly including snow-free areas in the Antarctic Peninsula and other areas as well). This also indicates that additional dust transport, which does not reach the EAP, must occur at low-tropospheric levels to this coastal sector of Antarctica.

(Table 1) GRACE-derived mass balance of the Greenland ice sheet 2002-2010

We re-evaluate the Greenland mass balance for the recent period using low-pass Independent Component Analysis (ICA) post-processing of the Level-2 GRACE data (2002-2010) from different official providers (UTCSR, JPL, GFZ) and confirm the present important ice mass loss in the range of -70 and -90 Gt/y of this ice sheet, due to negative contributions of the glaciers on the east coast. We highlight the high interannual variability of mass variations of the Greenland Ice Sheet (GrIS), especially the recent deceleration of ice loss in 2009-2010, once seasonal cycles are robustly removed by Seasonal Trend Loess (STL) decomposition. Interannual variability leads to varying trend estimates depending on the considered time span. Correction of post-glacial rebound effects on ice mass trend estimates represents no more than 8 Gt/y over the whole ice sheet. We also investigate possible climatic causes that can explain these ice mass interannual variations, as strong correlations between GRACE-based mass balance and atmosphere/ocean parallels are established: (1) changes in snow accumulation, and (2) the influence of inputs of warm ocean water that periodically accelerate the calving of glaciers in coastal regions and, feed-back effects of coastal water cooling by fresh currents from glaciers melting. These results suggest that the Greenland mass balance is driven by coastal sea surface temperature at time scales shorter than accumulation.

Elevation lines of modelled ice sheet extension during last glacial maximum

A numerical ice-sheet model was used to reconstruct the Late Weichselian glaciation of the Eurasian High Arctic, between Franz Josef Land and Severnaya Zemlya. An ice sheet was developed over the entire Eurasian High Arctic so that ice flow from the central Barents and Kara seas toward the northern Russian Arctic could be accounted for. An inverse approach to modeling was utilized, where ice-sheet results were forced to be compatible with geological information indicating ice-free conditions over the Taymyr Peninsula during the Late Weichselian. The model indicates complete glaciation of the Barents and Kara seas and predicts a "maximum-sized" ice sheet for the Late Weichselian Russian High Arctic. In this scenario, full-glacial conditions are characterized by a 1500-m-thick ice mass over the Barents Sea, from which ice flowed to the north and west within several bathymetric troughs as large ice streams. In contrast to this reconstruction, a "minimum" model of glaciation involves restricted glaciation in the Kara Sea, where the ice thickness is only 300 m in the south and which is free of ice in the north across Severnaya Zemlya. Our maximum reconstruction is compatible with geological information that indicates complete glaciation of the Barents Sea. However, geological data from Severnaya Zemlya suggest our minimum model is more relevant further east. This, in turn, implies a strong paleoclimatic gradient to colder and drier conditions eastward across the Eurasian Arctic during the Late Weichselian.

Glaciological characteristics and carbon dioxide record of ice cores from Antarctica

We report a decadally resolved record of atmospheric CO2 concentration for the last 1000 years, obtained from the West Antarctic Ice Sheet (WAIS) Divide shallow ice core. The most prominent feature of the pre-industrial period is a rapid ~7 ppm decrease of CO2 in a span of ~20-50 years at ~1600 A.D. This observation confirms the timing of an abrupt atmospheric CO2 decrease of ~10 ppm observed for that time period in the Law Dome ice core CO2 records, but the true magnitude of the decrease remains unclear. Atmospheric CO2 variations over the time period 1000-1800 A.D. are statistically correlated with northern hemispheric climate and tropical Indo-Pacific sea surface temperature. However, the exact relationship between CO2 and climate remains elusive due to regional climate variations and/or uneven geographical data density of paleoclimate records. We observe small differences of 0 ~2% (0 ~ 6 ppm) among the high-precision CO2 records from the Law Dome, EPICA Dronning Maud Land and WAIS Divide Antarctic ice cores. However, those records share common trends of CO2 change on centennial to multicentennial time scales, and clearly show that atmospheric CO2 has been increasing above preindustrial levels since ~1850 A.D.