Search

Marine-terminating outlet glacier terminus traces, mapped from satellite and aerial imagery, have been used extensively in understanding how outlet glaciers adjust to climate change variability over a range of timescales. Numerous studies have digitized termini manually, but this process is labor intensive, and no consistent approach exists. A lack of coordination leads to duplication of efforts, particularly for Greenland, which is a major scientific research focus. At the same time, machine learning techniques are rapidly making progress in their ability to automate accurate extraction of glacier termini, with promising developments across a number of optical and synthetic aperture radar (SAR) satellite sensors. These techniques rely on high-quality, manually digitized terminus traces to be used as training data for robust automatic traces. Here we present a database of manually digitized terminus traces for machine learning and scientific applications. These data have been collected, cleaned, assigned with appropriate metadata including image scenes, and compiled so they can be easily accessed by scientists. The TermPicks data set includes 39 060 individual terminus traces for 278 glaciers with a mean of 136 ± 190 and median of 93 of traces per glacier. Across all glaciers, 32 567 dates have been digitized, of which 4467 have traces from more than one author, and there is a duplication rate of 17 %. We find a median error of ∼ 100 m among manually traced termini. Most traces are obtained after 1999, when Landsat 7 was launched. We also provide an overview of an updated version of the Google Earth Engine Digitization Tool (GEEDiT), which has been developed specifically for future manual picking of the Greenland Ice Sheet.

Climate warming is rapidly reshaping the Arctic cryosphere and ocean conditions, with consequences for sea ice and pelagic productivity patterns affecting the entire marine food web. To predict how ongoing changes will impact Arctic marine ecosystems, concerted effort from various disciplines is required. Here, we contribute multi‐decadal reconstructions of changes in diatom production and sea‐ice conditions in relation to Holocene climate and ocean conditions off northwest Greenland. Our multiproxy study includes diatoms, sea‐ice biomarkers (IP25 and HBI III) and geochemical tracers (TOC [total organic carbon], TOC:TN [total nitrogen], δ13C, δ15N) from a sediment core record spanning the last c. 9,000 years. Our results suggest that the balance between the outflow of polar water from the Arctic, and input of Atlantic water from the Irminger Current into the West Greenland Current is a key factor in controlling sea‐ice conditions, and both diatom phenology and production in northeastern Baffin Bay. Our proxy record notably shows that changes in sea‐surface conditions initially forced by Neoglacial cooling were dynamically amplified by the shift in the dominant phase of the Arctic Oscillation (AO) mode that occurred at c. 3,000 yr BP, and caused drastic changes in community composition and a decline in diatom production at the study site. In the future, with projected dominant‐positive AO conditions favored by Arctic warming, increased water column stratification may counteract the positive effect of a longer open‐water growth season and negatively impact diatom production., We present the first continuous and high‐resolution diatom record from northeastern Baffin Bay spanning the last c. 9,000 years. Our data show that the variability in sea ice and primary production on the northwest Greenland shelf is linked to modulations in millennial‐scale Arctic Oscillation. The phenology and diatom community structure is largely controlled by water‐column stratification and sea‐ice conditions.

Baffin Bay is a semi-enclosed basin connecting the Arctic Ocean and the western North Atlantic, thus making out a significant pathway for heat exchange. Here we reconstruct the alternating advection of relatively warmer and saline Atlantic waters versus the incursion of colder Arctic water masses entering Baffin Bay through the multiple gateways in the Canadian Arctic Archipelago and the Nares Strait during the Holocene. We carried out benthic foraminiferal assemblage analyses, X-ray fluorescence scanning, and radiocarbon dating of a 738 cm long marine sediment core retrieved from eastern Baffin Bay near Upernavik, Greenland (Core AMD14-204C; 987 m water depth). Results reveal that eastern Baffin Bay was subjected to several oceanographic changes during the last 9.2 kyr. Waning deglacial conditions with enhanced meltwater influxes and an extensive sea-ice cover prevailed in eastern Baffin Bay from 9.2 to 7.9 ka. A transition towards bottom water amelioration is recorded at 7.9 ka by increased advection of Atlantic water masses, encompassing the Holocene Thermal Maximum. A cold period with growing sea-ice cover at 6.7 ka interrupts the overall warm subsurface water conditions, promoted by a weaker northward flow of Atlantic waters. The onset of the neoglaciation at ca. 2.9 ka is marked by an abrupt transition towards a benthic fauna dominated by agglutinated species, likely in part explained by a reduction of the influx of Atlantic Water, allowing an increased influx of the cold, corrosive Baffin Bay Deep Water originating from the Arctic Ocean to enter Baffin Bay through the Nares Strait. These cold subsurface water conditions persisted throughout the Late Holocene, only interrupted by short-lived warmings superimposed on this cooling trend.

Subglacial roughness can be determined at a variety of length scales from radio-echo sounding (RES) data either via statistical analysis of topography or inferred from basal radar scattering. Past studies have demonstrated that subglacial terrain exhibits self-affine (power law) roughness scaling behaviour, but existing radar scattering models do not take this into account. Here, using RES data from northern Greenland, we introduce a self-affine statistical framework that enables a consistent integration of topographic-scale roughness with the electromagnetic theory of radar scattering. We demonstrate that the degree of radar scattering, quantified using the waveform abruptness (pulse peakiness), is topographically controlled by the Hurst (roughness power law) exponent. Notably, specular bed reflections are associated with a lower Hurst exponent, with diffuse scattering associated with a higher Hurst exponent. Abrupt waveforms (specular reflections) have previously been used as a RES diagnostic for basal water, and to test this assumption we compare our radar scattering map with a recent prediction for the basal thermal state. We demonstrate that the majority of thawed regions (above pressure melting point) exhibit a diffuse scattering signature, which is in contradiction to the prior approach. Self-affine statistics provide a generalised model for subglacial terrain and can improve our understanding of the relationship between basal properties and ice-sheet dynamics.

Satellite imagery represents a unique opportunity to quantify the spatial and temporal changes of glaciers world-wide. Glacier velocity has been measured from repeat satellite scenes for decades now, yet a range of satellite missions, feature tracking programs, and user approaches have made it a laborious task. To date, there has been no tool developed that would allow a user to obtain displacement maps of any specified glacier simply by establishing the key temporal, spatial and feature tracking parameters. This work presents the application and development of a unique, semi-automatic, open-source, flexible processing toolbox for the retrieval of displacement maps with a focus on obtaining glacier surface velocities. SenDiT combines the download, pre-processing, feature tracking, and postprocessing of the highest resolution Sentinel-2A and Sentinel-2B satellite images into a semi-automatic toolbox, leaving a user with a set of rasterized and georeferenced glacier flow magnitude and direction maps for their further analyses. The solution is freely available and is tailored so that non-glaciologists and people with limited geographic information system (GIS) knowledge can also benefit from it. The system can be used to provide both regional and global sets of ice velocities. The system was tested and applied on a range of glaciers in mainland Norway, Iceland, Greenland and New Zealand. It was also tested on areas of stable terrain in Libya and Australia, where sources of error involved in the feature tracking using Sentinel-2 imagery are thoroughly described and quantified., This work is a contribution to the project Copernicus bretjeneste (Contract NIT.06.15.5), which is partly funded by the Norwegian Space Centre through Copernicus.

The Greenland Ice Sheet has experienced significant mass loss in recent years. A substantial component of this is attributable to the retreat of marine-terminating outlet glaciers, which lose mass through increases in calving, submarine melting and terrestrial meltwater discharge. In terms of iceberg production, Jakobshavn Isbræ is the largest marine-terminating glacier in Greenland, yet relatively little is known about its history before the first glacier margin observations in 1851. Two marine sediment cores obtained 15 and 19 km northwest from the mouth of Jakobshavn Isfjord were analysed to reconstruct the past behaviour of Jakobshavn Isbræ and to investigate the response of the glacier system to ocean forcing. These records provide long-term (~2000) context for assessing the significance of the rapid changes in glacier stability over the last century. The X-ray imagery and high-resolution grain size analysis from both cores reveal distinct multi-centennial-scale changes in the flux of iceberg-rafted debris (IRD) from Jakobshavn Isbræ. Foraminiferal analysis shows that variability in the relatively warm West Greenland Current (WGC) may have been an important driver of calving activity at Jakobshavn Isbræ. We find that iceberg rafting and WGC inflow were relatively high from onset of the record, at 60 BC, until AD 1100. Subsequently, the inflow of the WGC into Disko Bugt decreased. This was accompanied by a dramatic reduction in IRD from AD 1500 to 1850, which is attributed to the establishment of a floating ice tongue. We also show that ocean warming in the 20th century is part of a longer-term warming trend in the WGC which started at around AD 1700. Finally, these new records underline the complexity of glaciomarine sediments; IRD variability was driven by the inflow of the WGC but was also modulated by a complex interplay of air temperature, sea-ice coverage and ice margin proximity.

There is widespread, but often indirect, evidence that a significant fraction of the bed beneath the Greenland Ice Sheet is thawed (at or above the pressure melting point for ice). This includes the beds of major outlet glaciers and their tributaries and a large area around the NorthGRIP borehole in the ice-sheet interior. The ice-sheet-scale distribution of basal water is, however, poorly constrained by existing observations. In principle, airborne radio-echo sounding (RES) enables the detection of basal water from bed-echo reflectivity, but unambiguous mapping is limited by uncertainty in signal attenuation within the ice. Here we introduce a new, RES diagnostic for basal water that is associated with wet–dry transitions in bed material: bed-echo reflectivity variability. This technique acts as a form of edge detector and is a sufficient, but not necessary, criteria for basal water. However, the technique has the advantage of being attenuation insensitive and suited to combined analysis of over a decade of Operation IceBridge survey data.The basal water predictions are compared with existing analyses of the basal thermal state (frozen and thawed beds) and geothermal heat flux. In addition to the outlet glaciers, we demonstrate widespread water storage in the northern and eastern interior. Notably, we observe a quasilinear corridor of basal water extending from NorthGRIP to Petermann Glacier that spatially correlates with elevated heat flux predicted by a recent magnetic model. Finally, with a general aim to stimulate regional- and process-specific investigations, the basal water predictions are compared with bed topography, subglacial flow paths and ice-sheet motion. The basal water distribution, and its relationship with the thermal state, provides a new constraint for numerical models.

Understanding feedbacks between the Greenland Ice Sheet (GrIS) and the Atlantic Meridional Overturning Circulation (AMOC) is crucial for reducing uncertainties over future sea level and ocean circulation change. Reconstructing past GrIS dynamics can extend the observational record and elucidate mechanisms that operate on multi-decadal timescales. We report a highly-constrained last glacial vertical profile of cosmogenic isotope exposure ages from Sermilik Fjord, a marine-terminating ice stream in the southeast sector of the GrIS. Our reconstruction reveals substantial ice-mass loss throughout the Younger Dryas (12.9-11.7 ka), a period of marked atmospheric and sea-surface cooling. Earth-system modelling reveals that southern GrIS marginal melt was likely driven by strengthening of the Irminger Current at depth due to a weakening of the AMOC during the Younger Dryas. This change in North Atlantic circulation appears to have drawn warm subsurface waters to southeast Greenland despite markedly cooler sea surface temperatures, enhancing thermal erosion at the grounding lines of palaeo ice-streams, supporting interpretation of regional marine-sediment cores. Given current rates of GrIS meltwater input into the North Atlantic and the vulnerability of major ice streams to water temperature changes at the grounding line, this mechanism has important implications for future AMOC changes and northern hemisphere heat transport.

The transition from the last ice age to the present-day interglacial was interrupted by the Younger Dryas (YD) cold period. While many studies exist on this climate event, only few include high-resolution marine records that span the YD. In order to better understand the interactions between ocean, atmosphere and ice sheet stability during the YD, more high-resolution proxy records from the Arctic, located proximal to ice sheet outlet glaciers, are required. Here we present the first diatom-based high-resolution quantitative reconstruction of sea surface conditions from central-eastern Baffin Bay, covering the period 14.0–10.2 kyr BP. Our record reveals warmer sea surface conditions and strong interactions between the ocean and the West Greenland ice margin during the YD. These warmer conditions were caused by increased Atlantic-sourced water inflow combined with amplified seasonality. Our results emphasize the importance of the ocean for ice sheet stability under the current changing climate., We would like thank the captain and all members of the cruise MSM09/2. We also thank O. Hyttinen, University of Helsinki, for her help with grain size distribution analysis. Funding from the Finnish Graduate School in Geology is kindly acknowledged as is the funding from the Academy of Finland (A.E.K. Ojala, QUAL-project 259343). We are grateful to F.M. Nick, H. Machguth, T. Luoto, and K. Pauli for insightful discussions.

Jakobshavn Isbrae is a major marine terminating outlet glacier of the western Greenland Ice Sheet, which has been undergoing widespread acceleration and strong mass loss since the disintegration of its floating ice tongue in the late 1990s. The underlying mechanisms are poorly understood despite a wealth in observational and modelling studies. This doctoral thesis analyses the dynamic changes of Jakobshavn Isbrae using the Ice Sheet System Model (ISSM), a state-of-the-art finite-element ice flow model. Two missing model features for 1) the modelling the polythermal regime of glaciers and ice sheets, and 2) the dynamic evolution of its horizontal calving front position are designed and implemented into ISSM. A three-dimensional, thermodynamically coupled model of Jakobshavn Isbrae is set up and calibrated using modern observational data products. Low basal drag in the trough under the ice stream requires that its high driving stress is balanced by lateral drag in the shear margins, which allows for high flow velocities, as the ice viscosity is strain-rate-dependent. The developed modules are applied to the glacier model, which captures 90% of the observed changes from 1985 to 2015. Analysis of the model results reveals that calving front retreat is able to trigger widespread inland acceleration due to a rheological ice viscosity drop in the shear margins. Thermal feedbacks contribute 5 to 10% to the total acceleration. The study shows that Jakobshavn Isbrae will continue to contribute to eustatic sea level rise for at least the next century due to ongoing geometry adjustment to the new calving front position. Future fields of research include deriving a suitable calving rate parametrisation for large-scale ice flow models, a material law for temperate ice with a microscopic water content larger than 1%, and technical refinements of the modules developed for this thesis.

Calving glaciers contribute substantially to sea level rise, but they are challenging to represent in models. Fine resolution is required for continental-scale models to accurately resolve calving dynamics, and in many cases glacier geometry is too complicated to be adequately reflected by more simplified models. Flowline models are able to resolve flow along the main branch of a glacier, but many of those in current use either ignore tributaries entirely or parameterize their effect using a measure of “equivalent width.” Here we present a simple method to simulate terminus advance and retreat for an interacting network of glacier branches, based on a model extending Nye's (1953) perfect plastic flow approximation to calving glaciers. We apply the method to case studies of four marine-terminating glaciers: Jakobshavn Isbræ and Helheim Glacier of Greenland, and Columbia and Hubbard Glaciers of Alaska. Given bed topography and upstream elevation history, our method reproduces observed patterns of terminus advance and retreat in all cases, as well as centerline profiles for all branches.

One of the most prominent manifestations of the ongoing climate warming is the retreat of glaciers and ice sheets around the world. Retreating glaciers result in the formation of new ponds and lakes, which are available for colonization. The gradual appearance of these new habitat patches allows us to determine to what extent the composition of asexual Daphnia (water flea) populations is affected by environmental drivers vs. dispersal limitation. Here, we used a landscape genetics approach to assess the processes structuring the clonal composition of species in the D. pulex species complex that have colonized periglacial habitats created by ice-sheet retreat in western Greenland. We analysed 61 populations from a young ( 150 years) of lakes and ponds. We identified 42 asexual clones that varied widely in spatial distribution. Beta-diversity was higher among older than among younger systems. Lineage sorting by the environment explained 14% of the variation in clonal composition whereas the pure effect of geographical distance was very small and statistically insignificant (R-adj(2) = 0.010, P = 0.085). Dispersal limitation did not seem important, even among young habitat patches. The observation of several tens of clones colonizing the area combined with environmentally driven clonal composition of populations illustrates that population assembly of asexual species in the Arctic is structured by environmental gradients reflecting differences in the ecology of clones.

International audience; Monitoring the evolution of polar glaciers, ice caps and ice streams is of utmost importance because they constitute a good indicator of global climate change and contribute significantly to ongoing sea level rise. Accurate topographic surveys are particularly relevant as they reflect the geometric evolution of ice masses. Unfortunately, the precision and/or spatial coverage of current satellite missions (radar altimetry, ICESat) or field surveys are generally insufficient. Improving our knowledge of the topography of Polar Regions is the goal of the SPIRIT (SPOT 5 stereoscopic survey of Polar Ice: Reference Images and Topographies) international polar year (IPY) project. SPIRIT will allow (1) the acquisition of a large archive of SPOT 5 stereoscopic images covering most polar ice masses and, (2) the delivery of digital terrain models (DTM) to the scientific community. Here, we present the architecture of this project and the coverage achieved over northern and southern polar areas during the first year of IPY (July 2007 to April 2008). We also provide the first accuracy assessments of the SPIRIT DTMs. Over Jakobshavn Isbrae (West Greenland), SPIRIT elevations are within ± 6 m of ICESat elevations for 90% of the data. Some comparisons with ICESat profiles over Devon ice cap (Canada), St Elias Mountains (Alaska) and west Svalbard confirm the good overall quality of the SPIRIT DTMs although large errors are observed in the flat accumulation area of Devon ice cap. We then demonstrate the potential of SPIRIT DTMs for mapping glacier elevation changes. The comparison of summer-2007 SPIRIT DTMs with October-2003 ICESat profiles shows that the thinning of Jakobshavn Isbrae (by 30 to 40 m in 4 years) is restricted to the fast glacier trunk. The thinning of the coastal part of the ice stream (by over 100 m) and the retreat of its calving front (by up to 10 km) are clearly depicted by comparing the SPIRIT DTM to an ASTER April-2003 DTM.

We present a new record of ice thickness change, reconstructed at nearly 100,000 sites on the Greenland Ice Sheet (GrIS) from laser altimetry measurements spanning the period 1993-2012, partitioned into changes due to surface mass balance (SMB) and ice dynamics. We estimate a mean annual GrIS mass loss of 243 ± 18 Gt ⋅ y(-1), equivalent to 0.68 mm ⋅ y(-1) sea level rise (SLR) for 2003-2009. Dynamic thinning contributed 48%, with the largest rates occurring in 2004-2006, followed by a gradual decrease balanced by accelerating SMB loss. The spatial pattern of dynamic mass loss changed over this time as dynamic thinning rapidly decreased in southeast Greenland but slowly increased in the southwest, north, and northeast regions. Most outlet glaciers have been thinning during the last two decades, interrupted by episodes of decreasing thinning or even thickening. Dynamics of the major outlet glaciers dominated the mass loss from larger drainage basins, and simultaneous changes over distances up to 500 km are detected, indicating climate control. However, the intricate spatiotemporal pattern of dynamic thickness change suggests that, regardless of the forcing responsible for initial glacier acceleration and thinning, the response of individual glaciers is modulated by local conditions. Recent projections of dynamic contributions from the entire GrIS to SLR have been based on the extrapolation of four major outlet glaciers. Considering the observed complexity, we question how well these four glaciers represent all of Greenland's outlet glaciers.

Observations over the past decade show significant ice loss associated with the speed-up of glaciers in southeast Greenland from 2003, followed by a deceleration from 2006. These short-term, episodic, dynamic perturbations have a major impact on the mass balance on the decadal scale. To improve the projection of future sea level rise, a long-term data record that reveals the mass balance beyond such episodic events is required. Here, we extend the observational record of marginal thinning of Helheim and Kangerdlugssuaq glaciers from 10 to more than 80 years. We show that, although the frontal portion of Helheim Glacier thinned by more than 100 m between 2003 and 2006, it thickened by more than 50 m during the previous two decades. In contrast, Kangerdlugssuaq Glacier underwent minor thinning of 40–50 m from 1981 to 1998 and major thinning of more than 100 m after 2003. Extending the record back to the end of the Little Ice Age (prior to 1930) shows no thinning of Helheim Glacier from its maximum extent during the Little Ice Age to 1981, while Kangerdlugssuaq Glacier underwent substantial thinning of 230 to 265 m. Comparison of sub-surface water temperature anomalies and variations in air temperature to records of thickness and velocity change suggest that both glaciers are highly sensitive to short-term atmospheric and ocean forcing, and respond very quickly to small fluctuations. On century timescales, however, multiple external parameters (e.g. outlet glacier shape) may dominate the mass change. These findings suggest that special care must be taken in the projection of future dynamic ice loss.

Here we present a 100 year long proxy record for the renewal rate of the subsurface ocean waters in Sermilik Fjord at the edge of Helheim Glacier, based on investigations of two sediment cores (ER11-24 and ER11-25) obtained from the head of the fjord. By calculating the mean sortable silt ((SS) over bar) in current-sorted melt water plume sediments we find that episodes of increased water renewal rates lasting 3-5 years coincide with a positive North Atlantic Oscillation (NAO) index. This is not surprising as low pressure systems and northeasterly storms are observed more frequently along the east coast of Greenland during positive NAO years as a result of the northward shift in the North Atlantic storm track. Previous studies of sediment cores obtained from the mid-region of the fjord showed that Helheim Glacier destabilization coincides with a negative NAO index. Therefore we conclude that inter-annual variability in storm-induced flushing of Sermilik Fjord and thus the water renewal rate towards the glacier margin is not the controlling factor for inter-annual variability in Helheim Glacier destabilization. Such knowledge may have implications on predictive model studies of ice-ocean interactions and glacier behavior. (C) 2014 Elsevier Ltd. All rights reserved.

A multiproxy study of paleoceanographic and climatic changes in Holsteinsborg Dyb on the West Greenland shelf shows that the area is sensitive to variations in the mixing between the two components of the West Greenland Current, i.e., the Atlantic Water of the Irminger Current and the Polar Water from the East Greenland Current. A seismic profile illustrates the general development in the area since the deglaciation, and analyses of benthic foraminifera, stable isotopes (delta O-18 and delta C-13), and X-ray fluorescence elemental counts in two cores, covering ca. 6650 and 1100 years, respectively, provide details on the oceanographic changes during the mid- and late Holocene. The chronology is based on Pb-210 and AMS C-14 age determinations. During the time interval of 6650-1600 cal. yr BP, the bottom waters are characterized by the presence of modified Atlantic Water, with some influence of meltwater, either local or regional. Depleted delta O-18 values and relatively high Si/Fe ratio suggest warmer conditions from 6650 to 3800 cal. yr BP, corresponding in time to part of the Holocene Thermal Maximum, followed by a general slight cooling with some fluctuation after 3800 cal. yr BP. At 1600 cal. yr BP, there is an abrupt change in almost all investigated proxies. The foraminiferal assemblage and the delta O-18 values suggest cold, low-salinity bottom-water conditions, and there is indication of seasonal sea-ice cover between 1400 and 1200 cal. yr BP, corresponding in time to the Dark Ages Cold Period and to a minimum in solar irradiance. From ca. 1150 to 580 cal. yr BP, a time interval equivalent to the historical Medieval Warm Period (MWP), the foraminiferal fauna remains indicative of cold, low-salinity bottom waters, but without any indication of sea-ice cover. The influence of Polar Water/meltwater at the sea floor during the MWP is thought to be a result of enhanced mixing of the two components of the West Greenland Current caused by an increase in wind activity during generally positive North Atlantic Oscillation conditions, maybe combined with a displacement of the Irminger Current pathway to the outer shelf. During the Little Ice Age (LIA), there is a strong influence of modified Atlantic Water at the sea floor, related to sea-ice formation and stratification of the water column. After AD 1900, the benthic foraminiferal assemblage suggests a slight increase in Polar Water influence and thereby mixing, even though the environment is not much different from that during the LIA. (C) 2013 Elsevier B.V. All rights reserved.

A marine sediment core from Vaigat in Disko Bugt, West Greenland, has been analysed in terms of lithology, dinoflagellate cysts and foraminifera in order to evaluate the influence of oceanographic variability on West Greenland glacier stability. The data show that during the past 5200 years the Atlantic foraminiferal abundance in the subsurface waters of the West Greenland Current (WGC) episodically increased, indicating periods of increases in the inflow of subsurface warm Atlantic water at 2000—1500 cal. yr BP and 1300 cal. yr BP as well as periods of less pronounced increased bottom-water temperatures around 4700—4000 cal. yr BP, 3100—2800, 2600, 1000—800, 500—400, and at 200 cal. yr. The sedimentological and dinoflagellate cyst data indicate that these episodes with enhanced advection of Irminger Sea-derived waters are accompanied by increased iceberg rafting, which we link to increased iceberg calving in relation to destabilization of the Jakobshavn Isbrae. The long-term trend in the data documents the end of a late-Holocene Thermal Maximum between 5200 and 4300 cal. yr BP and a final onset of the Neoglaciation at 3500 cal. yr BP. Increased responses of the iceberg rafting after 3500 cal. yr BP, reflects a westward/seaward advance of the glacier margin in relation to onset of Neoglaciation and a development of the glacier into a floating tongue after 2000 cal. yr BP. A comparison of our record with a record from the eastern North Atlantic indicates that a NAO-like anomaly pattern between subsurface waters in West Greenland and atmospheric temperature in the Eastern North Atlantic may have been operating during most of the late Holocene. However, during the past 1000 years the NAO signal may have weakened as some other mode of climate variability overprints the anti-phase climate signal in this region.

During recent summers (2007–2012), several surface melt records were broken over the Greenland Ice Sheet (GrIS). The extreme summer melt resulted in part from a persistent negative phase of the North Atlantic Oscillation (NAO), favoring warmer atmospheric conditions than normal over the GrIS. Simultaneously, large anomalies in sea ice cover (SIC) and sea surface temperature (SST) were observed in the North Atlantic, suggesting a possible connection. To assess the direct impact of 2007–2012 SIC and SST anomalies on GrIS surface mass balance (SMB), a set of sensitivity experiments was carried out with the regional climate model MAR forced by ERA-Interim. These simulations suggest that perturbations in SST and SIC in the seas surrounding Greenland do not considerably impact GrIS SMB, as a result of the katabatic wind blocking effect. These offshore-directed winds prevent oceanic near-surface air, influenced by SIC and SST anomalies, from penetrating far inland. Therefore, the ice sheet SMB response is restricted to coastal regions, where katabatic winds cease. A topic for further investigation is how anomalies in SIC and SST might have indirectly affected the surface melt by changing the general circulation in the North Atlantic region, hence favoring more frequent warm air advection towards the GrIS.