Southwestern Greenland Research Articles

Modelling lateral meltwater flow and superimposed ice formation atop Greenland's near-surface ice slabs

Abstract At high elevations on the Greenland ice sheet meltwater percolates and refreezes in place, and hence does not contribute to mass loss. However, meltwater generation and associated surface runoff is occurring from increasingly higher altitudes, causing changes in firn stratigraphy that have led to the presence of near-surface ice slabs. These ice slabs force meltwater to flow laterally instead of percolating downwards. Here we present a simple, physics-based quasi-2-D model to simulate lateral meltwater runoff and superimposed ice (SI) formation on top of ice slabs. Using an Eulerian Darcy flow scheme, the model calculates how far meltwater can travel within a melt season and when it appears at the snow surface. Results show that lateral flow is a highly efficient runoff mechanism, as lateral outflow exceeds locally generated meltwater in all model gridcells, with total meltwater discharge sometimes reaching more than 30 times the average amount of in situ generated melt. SI formation, an important process in the formation and thickening of the ice slabs, can retain up to 40% of the available meltwater, and generally delays the appearance of visible runoff. Validating the model against field- or remote-sensing data remains challenging, but the results presented here are a first step towards a more comprehensive understanding and description of the hydrological system in the accumulation zone of the southwestern Greenland ice sheet.

Amide groups in 3.7 billion years old liquid inclusions

Carbon with depleted d13C (down to − 25.6‰ VPDB) found in > 3.7 billion year old metamorphic sediments from the Isua Supracrustal Belt, Southwestern Greenland, has been proposed to represent the oldest remains of life on Earth. Graphitic inclusions within garnet porphyroblasts from this locality have been shown to associate with elements consistent with biogenic remains. In this report, we focus on certain liquid inclusions found in the Isua garnets, characterizing their chemical composition using atomic force microscopy, AFM-based infrared spectroscopy, optical photothermal infrared spectroscopy, Raman spectroscopy, and time-of-flight secondary ion mass spectrometry. Our results show that the liquid inclusions contain functional groups consisting of carbon, nitrogen, and oxygen in a configuration similar to amide functional groups. We suspect that the amide groups formed from N, O and C-containing volatile components that were released from the original kerogenous material enclosed in the garnets, as this was graphitized during thermal maturation. This is consistent with the observed inclusion assemblage of solid graphitic and viscous fluid inclusions alike. Our observations are compatible with the inclusions forming from biogenic precursor material, and when considered alongside previous reports on the carbonaceous material in the Isua metamorphic sediments, these and our study collectively indicate that the carbonaceous material in the Isua metasediments represents the oldest traces of life on Earth.

Ice Darkening on Turner Glacier in Auyuittuq National Park, Nunavut

The glaciers of the Canadian Arctic are currently the third largest contributor to global sea level rise due to enhanced Arctic warming driving increased glacier melt and runoff (Zemp et al., 2019). Recent modelling suggests that the glaciers of the southern Qikiqtaaluk region, specifically Baffin and Bylot Islands, are 70% more sensitive to 1°C of warming (Brice et al., 2018). However, these studies lack field-based data to contextualize the drivers of glacier change in this region. Light absorbing particles (LAPs) such as dust, ash, and algae lower the albedo of snow and ice, which is the fraction of light a surface reflects, leading to an increase in the absorption of incoming radiation (Aubry-Wake et al., 2022). On glaciers in western Canada (Engstrom et al., 2022), the European Alps (Di Mauro et al., 2020) and south-western Greenland (Cook et al., 2020), positive feedback mechanisms have been discovered where summer snow and ice melt, augmented by LAPs, leads to the localized release of nutrients in ice and the subsequent enhanced growth of algae. Preliminary remote sensing analysis revealing a dark purple hue on glaciers in the Akshayuk Pass region of Auyuittuq National Park in Baffin Island, NU indicate these processes are likely occurring. The question of how this phenomenon is impacting glacier albedo, and therefore glacier melt which has implications for mass balance and regional sea level rise in the Canadian High Arctic has yet to be answered. Field-based sampling of LAPs on Turner Glacier, characterized through light microscopy and scanning electron microscopy, combined with in-situ hyperspectral spectroradiometer measurements (320nm-1100nm) of the sampled LAPs, and fine resolution (1cm) UAV surveys of sampling sites, are integral to develop a field-calibrated remote sensing approach to monitor LAP accumulation on the glacier surface and to assess the associated melt potential.

Rock avalanches in northeastern Baffin Island, Canada: understanding low occurrence amid high hazard potential

Rock avalanches in fjord environments can cause direct catastrophic damage and trigger secondary submarine landslides and tsunamis. These are well-documented in Greenland, Norway, and Alaska but have gone largely unreported in the extensive fjord terrain of the eastern Canadian Arctic. We provide the first inventory of rock avalanche deposits in northeastern Baffin Island—a region characterized by moderate to high seismic hazard, steep and high-walled fjords and glacial valleys, active deglaciation, and observed climate warming. Over a broad study area of ~60,000 km2, one sixth of the terrain had sufficient slope height and gradient to potentially generate rock avalanches. Within that hazard zone, we identified eight rock avalanche deposits at six locations. Only three rock avalanche deposits at two locations are dated, using aerial imagery (1958-present), to the last century while five deposits at four locations are inferred as syn- to post-glacial, likely occurring shortly after local debuttressing. These total numbers fall well below documented inventories from Greenland, Norway, and Alaska. We hypothesize that (1) continuous permafrost persists throughout this region and continues to act as a stabilizing factor and (2) rock mass quality is high in areas of most extreme relief contrast within the study region relative to analogous high-latitude fjord systems such as those in southwestern Greenland. We suggest that Baffin Island is currently in a period of quasi-stability that follows the intense instability during initial deglaciation, yet precedes the higher anticipated slope instability that may occur during permafrost degradation.

Air temperature changes in SW Greenland in the second half of the 18th century

Abstract. The thermal conditions of south-western Greenland in the second half of the 18th century were estimated using two unique series of meteorological observations. The first series (Neu-Herrnhut, 1 September 1767 to 22 July 1768, hereinafter 1767–1768) is the oldest long-term series of instrumental measurements of air temperature available for Greenland. The second (Godthaab, September 1784 to June 1792) contains the most significant and reliable data for Greenland for the study period. The quality-controlled and corrected data were used to calculate daily, monthly, seasonal and yearly means. The daily means were further used to calculate day-to-day temperature variability (DDTV), thermal seasons, growing degree days (GDDs), the air thawing index (ATI), positive degree days (PDDs) and air freezing index (AFI) degree days. Air temperature in Godthaab (now Nuuk) was, on average, warmer than the present day (1991–2020) in 1767–1768 and colder in 1784–1792. Compared to the early 20th-century Arctic warming (ETCAW) period, the data for the two sub-periods show that the late 18th century was as warm or even warmer. Except winter 1767/1768, winters and springs in the study period were longer, while summers and autumns were shorter than at present. The analysed climate indices usually did not exceed the maximum and minimum values from 1991–2020. Mean daily air temperature in the studied historical periods rarely exceed ±2 SD (standard deviation) of the long-term mean calculated for the contemporary period. Air temperature distribution was usually close to normal, in both historical and contemporary periods.

Evaluating middle to late Holocene climate variability from δ18O of aquatic invertebrate remains in southwestern Greenland

Holocene paleotemperature records from ice-free margins of Greenland suggest spatiotemporal heterogeneity in the timing of key climate transitions following retreat of the Greenland Ice Sheet. However, the paucity of high-resolution climate records from this region limits our understanding of the timing, spatial variation, and drivers of Holocene climate variability. Here, we present a middle to late Holocene climate reconstruction based on δ18O of aquatic invertebrate (chironomid and cladoceran) remains and sediment geochemistry from lake sediments to elucidate millennial to submillennial-scale drivers of climate in southwestern Greenland. We establish that our site, Arrowhead Lake, is a through-flowing basin with an integrated signal of annual meteoric water with no evidence of significant evaporative enrichment. Core sediment geochemistry captures the transition from minerogenic-dominated to organic-rich sediments ∼8 cal ka BP. Chitin-based δ18O measurements of both chironomid head capsules and ephippia exhibit similar temporal variability through the record, broadly decreasing from 8.2 cal ka BP to present, concurrent with decreasing summer insolation and regional Neoglacial cooling. A simplified assumption that temperature is the sole driver of changes in isotopes of precipitation yields an estimated average temperature change from the middle to late Holocene at Arrowhead Lake of 2.1–3.0 °C, comparable to existing estimates of middle to late Holocene temperature decline in western Greenland. Depleted δ18O values after 3.3 cal ka BP correspond to the onset of Neoglacial cooling and independent evidence for a slowdown of the West Greenland Current, and suggests that both temperature and sea-ice variability were important controls on isotopic variability at this site. Highly depleted isotopic values at 8.2 cal ka BP followed by high isotopic variability in the early part of the middle Holocene suggests that our record may capture regional 8.2 ka cooling. This new isotope record provides insight into climate influences during the middle and late Holocene along the western Greenland margin that are important for understanding how ongoing melt of the Greenland Ice Sheet and amplified warming in the Arctic may impact the climate system.

Lipid-based paleoecological and biogeochemical reconstruction of Store Saltsø, an extreme lacustrine system in SW Greenland

Polar lakes harbour a unique biogeochemistry that reflects the implications of climatic fluctuations against a susceptible yet extreme environment. In addition to polar, Store Saltsø (Kangerlussuaq, southwestern Greenland) is an endorheic lake with alkaline and oligotrophic waters that host a distinctive ecology adapted to live in such particular physico-chemical and environmental conditions. By exploring the sedimentary record of Store Saltsø at a molecular and compound-specific isotopic level, we were able to understand its ecology and biogeochemical evolution upon climate change. We employed lipid biomarkers to identify biological sources and metabolic traits in different environmental samples (shore terrace, sediment core, and white precipitates at the shore), and their succession over time to reconstruct the lake paleobiology. Different molecular ratios and geochemical proxies provided further insights toward the evolution of environmental conditions in the frame of the deglaciation history of Kangerlussuaq. The relative abundance of terrestrial (i.e., plant derived) biomarkers (odd long-chain n-alkanes, even long-chain n-alkanols, and phytosterols) in the upper half of the shore terrace versus the relatively more present aquatic biomarkers (botryococcenes and long-chain alkenones) in its lower half revealed higher lake water levels in the past. Moreover, the virtual absence of organics in the deepest section of the sediment core (32–29 cm depth) suggested that the lake did not yet exist at the northwestern shore of Store Saltsø ∼5100 years ago. According to the relative abundance of lipid biomarkers detected in the adjacent section above (29–25 cm depth), we hypothesize that the northwestern shore of Store Saltsø formed ∼4900 years ago. By combining the molecular and compound-specific isotopic analysis of lipids in a ∼360 cm sedimentary sequence, we recreated the paleobiology and evolution of an extreme lacustrine environment suitable for the study of the limits of life and the effects of climate warming.

Holocene southwest Greenland ice sheet behavior constrained by sea-level modeling

The melting of the Greenland ice sheet (GrIS) is a major contributor to past and future global sea-level rise. Understanding the response of the GrIS to times in the past when temperatures were as warm or warmer than today offers insights into its current and future response to climate change. In the southwest sector, the GrIS retreated inland beyond its current margin during the (at least regionally) warmer-than-present mid-Holocene, before it readvanced to the historical maximum position during the Little Ice Age. This was then followed by a slight retreat to its current position. To investigate the timing and magnitude of southwest GrIS retreat and readvance in response to Holocene warmth, we model the response of the solid Earth and local relative sea level (RSL) to past ice sheet change. We test a suite of eleven ice sheet scenarios that are based on ICE-6G_C but are modified in the timing and magnitude of ice retreat and readvance and pair them with four different viscoelastic Earth structures. We compare model predictions to observations of paleo sea level, present-day sea-level change, and present-day vertical land motion (VLM) around Nuuk, Greenland. We find that the modeled timing and magnitude of the Holocene retreat and readvance have a significant impact on modern sea-level change and VLM in Nuuk. Models that assume a readvance approaching the southwest GrIS’ historical maximum between 2 and 1 ka are most consistent with observations. The RSL response, however, is less sensitive to the timing of the minimum GrIS extent. Nonetheless, better data-model fits are generally obtained when the minimum ice sheet extent is reached between 5 and 3 ka, within the tested range of 6-3 ka. Comparing this timing to local and regional records of temperature and ice-sheet change suggest that the evolution of the southwestern GrIS presented here was in-phase with the likely evolution of southwestern GrIS mass balance through the Holocene. Our results have implications for future ice sheet modeling studies targeting southwestern Greenland by providing additional constraints and strengthening existing ones. Moreover, this work provides a deeper understanding of the interactions between the climate and the cryosphere and thus of future ice sheet change.

Improved monitoring of subglacial lake activity in Greenland

Abstract. Subglacial lakes form beneath ice sheets and ice caps if water is available and if bedrock and surface topography are able to retain the water. On a regional scale, the lakes modulate the timing and rate of freshwater flow through the subglacial system to the ocean by acting as reservoirs. More than 100 hydrologically active subglacial lakes that drain and recharge periodically have been documented under the Antarctic Ice Sheet, while only approximately 20 active lakes have been identified in Greenland. Active lakes may be identified by local changes in ice topography caused by the drainage or recharge of the lake beneath the ice. The small size of the Greenlandic subglacial lakes puts additional demands on mapping capabilities to resolve the evolving surface topography in sufficient detail to record their temporal behaviour. Here, we explore the potential for using CryoSat-2 swath-processed data, together with TanDEM-X digital elevation models, to improve the monitoring capabilities of active subglacial lakes in Greenland. We focus on four subglacial lakes previously described in the literature and combine the data with ArcticDEMs to obtain improved measurements of the evolution of these four lakes. We find that with careful tuning of the swath processor and filtering of the output data, the inclusion of these data, together with the TanDEM-X data, provides important information on lake activity, documenting, for example, that the ice surface collapse basin on Flade Isblink Ice Cap was 50 % (30 m) deeper than previously recorded. We also present evidence of a new, active subglacial lake in southwestern Greenland, which is located close to an already known lake. Both lakes probably drained within 1 month in the summer of 2012, which suggests either that they are hydrologically connected or that the drainages were independently triggered by extensive surface melt. If the hydrological connection is confirmed, this would to our knowledge be the first indication of hydrologically connected subglacial lakes in Greenland.

An automated algorithm to retrieve the location and depth of supraglacial lakes from ICESat-2 ATL03 data

Supraglacial lakes are widespread on the surface of the Antarctic Ice Sheet and Greenland Ice Sheet due to global warming. Estimating the location and depth of supraglacial lakes is critical for evaluating the impact of surface meltwater on ice dynamics. The accuracy of traditional image-based methods is limited, especially when supraglacial lakes are deep. Recently, Ice, Clouds, and land Elevation Satellite-2 data have been used to develop algorithms to measure the depth of supraglacial lakes. The ICESat-2-based methods are currently either semiautomated or sensitive to noise, which makes them less efficient and less accurate in detecting supraglacial lakes. In this study, we develop an automated algorithm for the retrieval of the location and depth of supraglacial lakes using ICESat-2 ATL03 data. We investigated 10 supraglacial lakes as case studies over Antarctica and Greenland (4 lakes on the Amery Ice Shelf, Antarctica and 6 lakes in southwestern Greenland). The results show that our algorithm can detect supraglacial lakes at depths up to 8.25 m. The RMSE for the lake depth retrievals in Antarctica is approximately 0.30 m and that of the lakes in Greenland is approximately 0.32 m. Our algorithm reveals high robustness and accuracy, especially in shallow areas. Additionally, the fully automated algorithm does not require prior knowledge and guarantees efficiency when there is an extensive study area. This result implies that our algorithm has the potential to address large amounts of data and may even be useful to estimate meltwater volumes across the surface of the entire ice sheet.

Herbivore Diversity Helps Maintain Arctic Tundra Diversity

A long-term experiment in southwestern Greenland reveals that the presence of musk oxen and caribou helps stave off declines in Arctic tundra diversity brought on by climate change.

Spatial and temporal differences in surface and subsurface meltwater distribution over Greenland ice sheet using multi-frequency passive microwave observations

Increasingly larger portions of the Greenland ice sheet are undergoing seasonal melting-refreeze cycles due to global climate warming. The cycle begins with the arrival of high temperatures and increased solar radiation in the spring and summer seasons generating meltwater on the ice sheet's surface. Meltwater percolates to deeper ice layers, either refreezing within the firn, creating longer-term meltwater pockets (firn aquifers), or generating peripheral runoff. Depending on the location and climate, the refreeze duration, the depth of infiltration, and meltwater persistence are temporally and spatially complex. Our recent study showed that multi-frequency passive microwave measurements in the 1.4 GHz to 36.5 GHz range effectively distinguished seasonal meltwater between the immediate surface and deeper firn layers at an experiment site in the accumulation zone of the southwestern Greenland ice sheet. Here, we further explored the vertically and horizontally polarized multi-frequency melt response at the pan-Greenland scale. We employed 1.4 GHz brightness temperature (TB) measurements from the NASA Soil Moisture Active Passive (SMAP) satellite and 6.9, 10.7, 18.9, and 36.5 GHz TB measurements from the JAXA Global Change Observation Mission-Water Shizuku (GCOM-W) satellite. The results show that the frequency-dependent response was consistent across the ice sheet. The multi-frequency melt indications match with lasting seasonal subsurface meltwater with delayed refreezing compared to the surface. These results suggest persistent seasonal subsurface meltwater occurrences that are spatially and temporally significant but concealed from the high-frequency observations. Retrieving the meltwater evolution in snow and firn presents a complex problem; this work represents an initial step toward developing an ice-sheet-wide algorithm for more comprehensive retrieval of the meltwater profile.

The impact of surface melt rate and catchment characteristics on Greenland Ice Sheet moulin inputs

Abstract. The supraglacial drainage system of the Greenland Ice Sheet, in combination with surface melt rate, controls the rate of water flow into moulins, a major driver of subglacial water pressure. We apply the Subaerial Drainage System (SaDS) model, a physically based surface meltwater flow model, to a ∼20×27km2 catchment on the southwestern Greenland Ice Sheet for 4 years of melt forcing (2011, 2012, 2015, and 2016) to (1) examine the relationship between surface melt rate and the rate, diurnal amplitude, and timing of surface inputs to moulins; (2) compare SaDS to contemporary models; and (3) present a framework for selecting appropriate supraglacial drainage models for different modelling objectives. We find that variations in the rate and timing of modelled moulin inputs related to the development of supraglacial channels are relatively more important in years with low melt volumes than years with high melt volumes. We suggest that a process-resolving supraglacial hydrology model (e.g., SaDS) should be considered when modelling outcomes are sensitive to subdiurnal and long-term seasonal changes in the rate of discharge into moulins.

Sea-level rise in Southwest Greenland as a contributor to Viking abandonment

The first records of Greenland Vikings date to 985 CE. Archaeological evidence yields insight into how Vikings lived, yet drivers of their disappearance in the 15th century remain enigmatic. Research suggests a combination of environmental and socioeconomic factors, and the climatic shift from the Medieval Warm Period (~900 to 1250 CE) to the Little Ice Age (~1250 to 1900 CE) may have forced them to abandon Greenland. Glacial geomorphology and paleoclimate research suggest that the Southern Greenland Ice Sheet readvanced during Viking occupation, peaking in the Little Ice Age. Counterintuitively, the readvance caused sea-level rise near the ice margin due to increased gravitational attraction toward the ice sheet and crustal subsidence. We estimate ice growth in Southwestern Greenland using geomorphological indicators and lake core data from previous literature. We calculate the effect of ice growth on regional sea level by applying our ice history to a geophysical model of sea level with a resolution of ~1 km across Southwestern Greenland and compare the results to archaeological evidence. The results indicate that sea level rose up to ~3.3 m outside the glaciation zone during Viking settlement, producing shoreline retreat of hundreds of meters. Sea-level rise was progressive and encompassed the entire Eastern Settlement. Moreover, pervasive flooding would have forced abandonment of many coastal sites. These processes likely contributed to the suite of vulnerabilities that led to Viking abandonment of Greenland. Sea-level change thus represents an integral, missing element of the Viking story.

A regime shift in North Pacific annual mean sea surface temperature in 2013/14

This study revealed a new regime shift in the North Pacific annual mean sea surface temperature (SST) in 2013/14, robust in multiple SST data sets. This regime shift shows a horseshoe pattern with warming SST along the North American western coast and toward the Central Pacific. It suggests a phase reversal of the Victoria Mode (VM) and it seems unrelated to the shift in the Pacific Decadal Oscillation (PDO) and Interdecadal Pacific Oscillation (IPO). Associated with this regime shift, the annual surface air temperature significantly enhanced over the high-latitude landmass after 2014, whereas it declined over northeastern Canada and southwestern Greenland. The annual precipitation generally increased (decreased) in the Northern (Southern) Hemisphere landmass monsoon regions. There are no regime shifts in Equatorial Central-East Pacific SST and the intensity and locations of the Aleutian Low around 2013/14. These facts suggest the regime shift of North Pacific SST in 2013/14 may not originate from the decadal variations of the midlatitude atmosphere and tropical ocean. The extensive influence of this regime shift in 2013/14 on the atmospheric, oceanic and ecological systems warrants further investigation.

Whole genome population structure of North Atlantic kelp confirms high-latitude glacial refugia.

Coastal refugia during the Last Glacial Maximum (~21,000 years ago) have been hypothesized at high latitudes in the North Atlantic, suggesting marine populations persisted through cycles of glaciation and are potentially adapted to local environments. Here, whole-genome sequencing was used to test whether North Atlantic marine coastal populations of the kelp Alaria esculenta survived in the area of southwestern Greenland during the Last Glacial Maximum. We present the first annotated genome for A. esculenta and call variant positions in 54 individuals from populations in Atlantic Canada, Greenland, Faroe Islands, Norway and Ireland. Differentiation across populations was reflected in ~1.9 million single nucleotide polymorphisms, which further revealed mixed ancestry in the Faroe Islands individuals between putative Greenlandic and European lineages. Time-calibrated organellar phylogenies suggested Greenlandic populations were established during the last interglacial period more than 100,000 years ago, and that the Faroe Islands population was probably established following the Last Glacial Maximum. Patterns in population statistics, including nucleotide diversity, minor allele frequencies, heterozygosity and linkage disequilibrium decay, nonetheless suggested glaciation reduced Canadian Atlantic and Greenlandic populations to small effective sizes during the most recent glaciation. Functional differentiation was further reflected in exon read coverage, which revealed expansions unique to Greenland in 337 exons representing 162 genes, and a modest degree of exon loss (103 exons from 56 genes). Altogether, our genomic results provide strong evidence that A. esculenta populations were resilient to past climatic fluctuations related to glaciations and that high-latitude populations are potentially already adapted to local conditions as a result.

In situ measurements of meltwater flow through snow and firn in the accumulation zone of the SW Greenland Ice Sheet

Abstract. The Greenland Ice Sheet is losing mass, part of which is caused by increasing runoff. The location of the runoff limit, the highest elevation from which meltwater finds its way off the ice sheet, plays an important role in the surface mass balance of the ice sheet. The recently observed rise in runoff area might be related to an increasing amount of refreezing: ice layer development in the firn reduces vertical percolation and promotes lateral runoff. To investigate meltwater flow near the runoff limit in the accumulation zone on the southwestern Greenland Ice Sheet, we carried out in situ measurements of hydrological processes and properties of firn and snow. The hydraulic conductivity of icy firn in pre-melt conditions measured using a portable lysimeter ranges from 0.17 to 12.8 m h−1, with flow predominantly occurring through preferential flow fingers. Lateral flow velocities of meltwater on top of the near-surface ice slab, measured at the peak of the melt season by salt dilution and tracer experiments, range from 1.3 to 15.1 m h−1. With these lateral flow velocities, the distance between the slush limit, the highest elevation where liquid water is visible on the ice sheet surface, and the runoff limit could be roughly 4 km in regions where near-surface ice slabs are present. These measurements are a first step towards an integrated set of hydrological properties of firn on the SW Greenland Ice Sheet and show evidence that meltwater runoff may occur from elevations above the visible runoff area.

Lu–Hf, Sm–Nd, and U–Pb isotopic coupling and decoupling in apatite

Apatite is a useful geochronological tool due to its common occurrence and incorporation of a variety of radioactive parent isotopes during crystallisation. However, an understanding of the geological conditions recorded by each of these isotopic systems in apatite is necessary for their effective application to geological problems. In this study, U–Pb, Sm–Nd, and Lu–Hf dates for apatite grains from samples of ∼3.0 Ga TTG gneisses in the Akia terrane of the North Atlantic Craton (south-western Greenland) were obtained in order to better understand the conditions under which isotopic re-equilibration in apatite is achieved. This isotopic data records at least two thermotectonic events subsequent to the initial magmatic crystallisation of these rocks. Whereas the Lu–Hf and Sm–Nd isotopic systems in these apatites record a regional high temperature metamorphic event at ∼2.7 Ga, the U–Pb system records a milder thermal event at ∼1.7 Ga. This decoupling of the U–Pb from the Sm–Nd and Lu–Hf systems in apatite reflects a difference in closure temperature. The regional thermal effects of the ∼1.7 Ga Nagssugtoqidian orogeny in adjacent parts of the North Atlantic Craton resulted in partial to complete reequilibration of the U–Pb systems in apatite across the Akia terrane but were insufficient to cause open system behaviour of Sm–Nd or Lu–Hf. The indistinguishable apatite Sm–Nd and Lu–Hf ages produced in each analysed sample suggest that these systems are similarly susceptible to isotopic re-equilibration. The older ages produced by the Sm–Nd and Lu–Hf systems relative to U–Pb in apatite demonstrate that the former are more robust and thus more likely to faithfully record high temperature events than the U–Pb system.

Morphological variability and diagnostics of two sculpins Icelus bicornis and I. spatula (Cottiformes: Cottidae) from the Arctic

The morphological variability of Icelus spatula Gilbert et Burke, 1912 and I. bicornis Reinhardt, 1840 was studied in order to clarify the diagnostic characters of these poorly studied fishes. Allometry, geographic variability and sexual dimorphism of species were evaluated by statistical methods based on 24 counts and 30 morphometric indices. Samples from 5 Arctic regions were studied: southwestern Greenland, Spitsbergen, Barents Sea, Laptev and East Siberian seas. An analysis of the variability of diagnostic characters showed that the following can be used to identify species with high certainty: 1) the shape of the urogenital papilla (males); 2) the presence/absence of spinules on the lower-posterior margin of scales of the lateral line below the pore of the canal and the degree of development of spinules on the dorsal margin of these scales; 3) the presence of scales at the base of the caudal fin; 4) presence/absence of bony granules on the skin between the lateral line and dorsal row of scales. Other characters (the number of scales in the lateral line and in rows on the body behind the base of the pectoral fin, the length of the nuchal spines) are more variable. No clinal changes in meristic characters were found over the studied areas of ranges of the two species. Greenland population of I. bicornis differs markedly from the samples from the Barents and Laptev seas, which suggests the possibility of the existence of groups of subspecies rank in these waters. An updated key for identification of juveniles, females and males of I. bicornis and I. spatula is provided.

Simulating the Holocene deglaciation across a marine-terminating portion of southwestern Greenland in response to marine and atmospheric forcings

Abstract. Numerical simulations of the Greenland Ice Sheet (GrIS) over geologic timescales can greatly improve our knowledge of the critical factors driving GrIS demise during climatically warm periods, which has clear relevance for better predicting GrIS behavior over the upcoming centuries. To assess the fidelity of these modeling efforts, however, observational constraints of past ice sheet change are needed. Across southwestern Greenland, geologic records detail Holocene ice retreat across both terrestrial-based and marine-terminating environments, providing an ideal opportunity to rigorously benchmark model simulations against geologic reconstructions of ice sheet change. Here, we present regional ice sheet modeling results using the Ice-sheet and Sea-level System Model (ISSM) of Holocene ice sheet history across an extensive fjord region in southwestern Greenland covering the landscape around the Kangiata Nunaata Sermia (KNS) glacier and extending outward along the 200 km Nuup Kangerula (Godthåbsfjord). Our simulations, forced by reconstructions of Holocene climate and recently implemented calving laws, assess the sensitivity of ice retreat across the KNS region to atmospheric and oceanic forcing. Our simulations reveal that the geologically reconstructed ice retreat across the terrestrial landscape in the study area was likely driven by fluctuations in surface mass balance in response to Early Holocene warming – and was likely not influenced significantly by the response of adjacent outlet glaciers to calving and ocean-induced melting. The impact of ice calving within fjords, however, plays a significant role by enhancing ice discharge at the terminus, leading to interior thinning up to the ice divide that is consistent with reconstructed magnitudes of Early Holocene ice thinning. Our results, benchmarked against geologic constraints of past ice-margin change, suggest that while calving did not strongly influence Holocene ice-margin migration across terrestrial portions of the KNS forefield, it strongly impacted regional mass loss. While these results imply that the implementation and resolution of ice calving in paleo-ice-flow models is important towards making more robust estimations of past ice mass change, they also illustrate the importance these processes have on contemporary and future long-term ice mass change across similar fjord-dominated regions of the GrIS.