Asiaq Greenland Survey

ArcticDEM is an NGA-NSF public-private initiative to automatically produce a high-resolution, high quality, digital surface model (DSM) of the Arctic using optical stereo imagery, high-performance computing, and open source photogrammetry software. This record represents both the version 3 mosaics and the version 3 strip DEMs.

Abstract This paper provides a snapshot of the permafrost thermal state in the Nordic area obtained during the International Polar Year (IPY) 2007–2009. Several intensive research campaigns were undertaken within a variety of projects in the Nordic countries to obtain this snapshot. We demonstrate for Scandinavia that both lowland permafrost in palsas and peat plateaus, and large areas of permafrost in the mountains are at temperatures close to 0°C, which makes them sensitive to climatic changes. In Svalbard and northeast Greenland, and also in the highest parts of the mountains in the rest of the Nordic area, the permafrost is somewhat colder, but still only a few degrees below the freezing point. The observations presented from the network of boreholes, more than half of which were established during the IPY, provide an important baseline to assess how future predicted climatic changes may affect the permafrost thermal state in the Nordic area. Time series of active‐layer thickness and permafrost temperature conditions in the Nordic area, which are generally only 10 years in length, show generally increasing active‐layer depths and rising permafrost temperatures. Copyright © 2010 John Wiley & Sons, Ltd.

Abstract The grey wolf ( Canis lupus ) was the first species to give rise to a domestic population, and they remained widespread throughout the last Ice Age when many other large mammal species went extinct. Little is known, however, about the history and possible extinction of past wolf populations or when and where the wolf progenitors of the present-day dog lineage ( Canis familiaris ) lived 1–8 . Here we analysed 72 ancient wolf genomes spanning the last 100,000 years from Europe, Siberia and North America. We found that wolf populations were highly connected throughout the Late Pleistocene, with levels of differentiation an order of magnitude lower than they are today. This population connectivity allowed us to detect natural selection across the time series, including rapid fixation of mutations in the gene IFT88 40,000–30,000 years ago. We show that dogs are overall more closely related to ancient wolves from eastern Eurasia than to those from western Eurasia, suggesting a domestication process in the east. However, we also found that dogs in the Near East and Africa derive up to half of their ancestry from a distinct population related to modern southwest Eurasian wolves, reflecting either an independent domestication process or admixture from local wolves. None of the analysed ancient wolf genomes is a direct match for either of these dog ancestries, meaning that the exact progenitor populations remain to be located.

Land surface temperature (LST) is one of the sources of input data for modeling land surface processes. The Landsat satellite series is the only operational mission with more than 30 years of archived thermal infrared imagery from which we can retrieve LST. Unfortunately, stray light artifacts were observed in Landsat-8 TIRS data, mostly affecting Band 11, currently making the split-window technique impractical for retrieving surface temperature without requiring atmospheric data. In this study, a single-channel methodology to retrieve surface temperature from Landsat TM and ETM+ was improved to retrieve LST from Landsat-8 TIRS Band 10 using near-surface air temperature (Ta) and integrated atmospheric column water vapor (w) as input data. This improved methodology was parameterized and successfully evaluated with simulated data from a global and robust radiosonde database and validated with in situ data from four flux tower sites under different types of vegetation and snow cover in 44 Landsat-8 scenes. Evaluation results using simulated data showed that the inclusion of Ta together with w within a single-channel scheme improves LST retrieval, yielding lower errors and less bias than models based only on w. The new proposed LST retrieval model, developed with both w and Ta, yielded overall errors on the order of 1 K and a bias of −0.5 K validated against in situ data, providing a better performance than other models parameterized using w and Ta or only w models that yielded higher error and bias.

Abstract. Glacier inventories provide the basis for further studies on mass balance and volume change, relevant for local hydrological issues as well as for global calculation of sea level rise. In this study, a new Austrian glacier inventory has been compiled, updating data from 1969 (GI 1) and 1998 (GI 2) based on high-resolution lidar digital elevation models (DEMs) and orthophotos dating from 2004 to 2012 (GI 3). To expand the time series of digital glacier inventories in the past, the glacier outlines of the Little Ice Age maximum state (LIA) have been digitalized based on the lidar DEM and orthophotos. The resulting glacier area for GI 3 of 415.11 ± 11.18 km2 is 44% of the LIA area. The annual relative area losses are 0.3% yr−1 for the ~119-year period GI LIA to GI 1 with one period with major glacier advances in the 1920s. From GI 1 to GI 2 (29 years, one advance period of variable length in the 1980s) glacier area decreased by 0.6% yr−1 and from GI 2 to GI 3 (10 years, no advance period) by 1.2% yr−1. Regional variability of the annual relative area loss is highest in the latest period, ranging from 0.3 to 6.19% yr−1. The mean glacier size decreased from 0.69 km2 (GI 1) to 0.46 km2 (GI 3), with 47% of the glaciers being smaller than 0.1 km2 in GI 3 (22%).

Abstract. Greenland runoff, from ice mass loss and increasing rainfall, is increasing. That runoff, as discharge, impacts the physical, chemical, and biological properties of the adjacent fjords. However, where and when the discharge occurs is not readily available in an open database. Here we provide data sets of high-resolution Greenland hydrologic outlets, basins, and streams, as well as a daily 1958 through 2019 time series of Greenland liquid water discharge for each outlet. The data include 24 507 ice marginal outlets and upstream basins and 29 635 land coast outlets and upstream basins, derived from the 100 m ArcticDEM and 150 m BedMachine. At each outlet there are daily discharge data for 22 645 d – ice sheet runoff routed subglacially to ice margin outlets and land runoff routed to coast outlets – from two regional climate models (RCMs; MAR and RACMO). Our sensitivity study of how outlet location changes for every inland cell based on subglacial routing assumptions shows that most inland cells where runoff occurs are not highly sensitive to those routing assumptions, and outflow location does not move far. We compare RCM results with 10 gauges from streams with discharge rates spanning 4 orders of magnitude. Results show that for daily discharge at the individual basin scale the 5 % to 95 % prediction interval between modeled discharge and observations generally falls within plus or minus a factor of 5 (half an order of magnitude, or +500 %/-80 %). Results from this study are available at https://doi.org/10.22008/promice/freshwater (Mankoff, 2020a) and code is available at http://github.com/mankoff/freshwater (last access: 6 November 2020) (Mankoff, 2020b).

Growing season conditions are widely recognized as the main driver for tundra shrub radial growth, but the effects of winter warming and snow remain an open question. Here, we present a more than 100 years long Betula nana ring-width chronology from Disko Island in western Greenland that demonstrates a highly significant and positive growth response to both summer and winter air temperatures during the past century. The importance of winter temperatures for Betula nana growth is especially pronounced during the periods from 1910-1930 to 1990-2011 that were dominated by significant winter warming. To explain the strong winter importance on growth, we assessed the importance of different environmental factors using site-specific measurements from 1991 to 2011 of soil temperatures, sea ice coverage, precipitation and snow depths. The results show a strong positive growth response to the amount of thawing and growing degree-days as well as to winter and spring soil temperatures. In addition to these direct effects, a strong negative growth response to sea ice extent was identified, indicating a possible link between local sea ice conditions, local climate variations and Betula nana growth rates. Data also reveal a clear shift within the last 20 years from a period with thick snow depths (1991-1996) and a positive effect on Betula nana radial growth, to a period (1997-2011) with generally very shallow snow depths and no significant growth response towards snow. During this period, winter and spring soil temperatures have increased significantly suggesting that the most recent increase in Betula nana radial growth is primarily triggered by warmer winter and spring air temperatures causing earlier snowmelt that allows the soils to drain and warm quicker. The presented results may help to explain the recently observed 'greening of the Arctic' which may further accelerate in future years due to both direct and indirect effects of winter warming.

Active rock glaciers represent the best visual expression of mountain permafrost that can be mapped and monitored directly using remotely sensed data. Active rock glaciers are bodies that consist of a perennially frozen ice/rock mixture and express a distinct flow-like morphology indicating downslope permafrost creep movement. Annual rates of motion have ranged from a few millimeters to several meters per year, varying within the annual cycle, from year to year, as well as at the decennial time scale. During the last decade, in situ observations in the European Alps have shown that active rock glaciers are responding almost synchronously to inter-annual and decennial changes in ground temperature, suggesting that the relative changes of their kinematics are a general indicator of the evolution of mountain permafrost conditions. Here, we used satellite radar interferometry (InSAR) to monitor the rate of motion of various active rock glaciers in the Swiss Alps, Qeqertarsuaq (Western Greenland), and the semiarid Andes of South America. Velocity time series computed with Sentinel-1 SAR images, regularly acquired since 2014, every six days over Europe and Greenland and every 12 days over the Andes, show annual fluctuations, with higher velocities at the end of the summer. A JERS-1 image pair of 1996 and stacks of very high-resolution SAR images from TerraSAR-X and Cosmo-SkyMed from 2008 to 2017 were analyzed using InSAR and offset tracking over the Western Swiss Alps in order to extend the main observation period of our study. A quantitative assessment of the accuracy of InSAR and offset tracking was performed by comparison with in situ methods. Our results for the three different study regions demonstrate that Sentinel-1 InSAR can complement worldwide in situ measurements of active rock glacier kinematics.

Abstract A degree-day model developed for parameterizing melt rates on the Greenland ice sheet is adapted to the climatic conditions on glaciers in southern Norway. The model is calibrated by means of observed average mass-balance-elevation relationships (1963–90) for three glaciers in a west-east transect in southern Norway and 30 year normals (1961–90) of temperature and precipitation observed at nearby climate stations. The calibration gives a surprisingly small variation of the model parameters (degree-day factors for snow-and ice-melt, and precipitation-elevation gradient) from one glacier to another. The derived values of the parameters are used to estimate the change of the mass-balance-elevation relationship for different climatic scenarios. The study indicates that a low-lying glacier in the maritime, high-precipitation environment near the Atlantic coast is more sensitive to both temperature and precipitation changes than the high elevated glaciers in the dry, more continental climate farther away from the coast. However, all of the glaciers studied will lose mass in a warmer climate, unless the warming is accompanied by a dramatic increase in the precipitation of 25–40% deg −1 warming.

In the Arctic domain, there is evidence for very long distance transport of pollen grains from boreal forests to tundra environments. However, the sampling protocol used in earlier studies did not allow the determination of precise timing of the deposition of exotic grains. The ongoing “Epilobe” project monitors the deposition of exotic grains on the western and eastern coasts of Greenland utilizing wind filters, which are changed weekly throughout the entire year. A combination of the identification of tree pollen grains, the dates of deposition of the pollen on the filters, and the distribution map of the trees identified, allows the selection of modeled backward trajectories of air parcels responsible for the capture of the grains in the growing area of northeastern North America and their long‐distance transport to Greenland. A survey of data obtained from four stations, analyzed during 2004 and 2005, indicates the occurrence of a general pattern every spring, which follows the main cyclone tracks reaching this Arctic region.

Grey wolves (Canis lupus) are one of the few large terrestrial carnivores that have maintained a wide geographical distribution across the Northern Hemisphere throughout the Pleistocene and Holocene. Recent genetic studies have suggested that, despite this continuous presence, major demographic changes occurred in wolf populations between the Late Pleistocene and early Holocene, and that extant wolves trace their ancestry to a single Late Pleistocene population. Both the geographical origin of this ancestral population and how it became widespread remain unknown. Here, we used a spatially and temporally explicit modelling framework to analyse a data set of 90 modern and 45 ancient mitochondrial wolf genomes from across the Northern Hemisphere, spanning the last 50,000 years. Our results suggest that contemporary wolf populations trace their ancestry to an expansion from Beringia at the end of the Last Glacial Maximum, and that this process was most likely driven by Late Pleistocene ecological fluctuations that occurred across the Northern Hemisphere. This study provides direct ancient genetic evidence that long-range migration has played an important role in the population history of a large carnivore, and provides insight into how wolves survived the wave of megafaunal extinctions at the end of the last glaciation. Moreover, because Late Pleistocene grey wolves were the likely source from which all modern dogs trace their origins, the demographic history described in this study has fundamental implications for understanding the geographical origin of the dog.

Recirculation of wood ash from energy production to forest soil improves the sustainability of this energy production form as recycled wood ash contains nutrients that otherwise would be lost at harvest. In addition, wood-ash is beneficial to many soils due to its inherent acid-neutralizing capabilities. However, wood ash has several ecosystem-perturbing effects like increased soil pH and pore water electrical conductivity both known to strongly impact soil bacterial numbers and community composition. Studies investigating soil bacterial community responses to wood ash application remain sparse and the available results are ambiguous and remain at a general taxonomic level. Here we investigate the response of bacterial communities in a spruce forest soil to wood ash addition corresponding to 0, 5, 22 and 167 t wood ash ha-1. We used culture-based enumerations of general bacteria, Pseudomonas and sporeforming bacteria combined with 16S rRNA gene amplicon sequencing to valuate soil bacterial responses to wood ash application. Results showed that wood ash addition strongly increased soil pH and electrical conductivity. Soil pH increased from acidic through neutral at 22 t ha-1 to alkaline at 167 t ha-1. Bacterial numbers significantly increased up to a wood ash dose of 22 t ha-1 followed by significant decrease at 167 t ha-1 wood ash. The soil bacterial community composition changed after wood ash application with copiotrophic bacteria responding positively up to a wood ash dose of 22 t ha-1 while the adverse effect was seen for oligotrophic bacteria. Marked changes in bacterial community composition occurred at a wood ash dose of 167 t ha-1 with a single alkaliphilic genus dominating. Additionally, spore-formers became abundant at an ash dose of 167 t ha-1 whereas this was not the case at lower ash doses. Lastly, bacterial richness and diversity strongly decreased with increasing amount of wood ash applied. All of the observed bacterial responses can be directly explained by the wood ash induced changes in pH, electrical conductivity and the addition of wood ash inherent nutrients.

This study examines the relationships linking methanesulfonic acid (MSA, arising from the atmospheric oxidation of the biogenic dimethylsulfide, DMS) in atmospheric aerosol, satellite-derived chlorophyll a (Chl-a), and oceanic primary production (PP), also as a function of sea ice melting (SIM) and extension of the ice free area in the marginal ice zone (IF-MIZ) in the Arctic. MSA was determined in PM10 samples collected over the period 2010–2012 at two Arctic sites, Ny Ålesund (78.9°N, 11.9°E), Svalbard islands, and Thule Air Base (76.5°N, 68.8°W), Greenland. PP is calculated by means of a bio-optical, physiologically based, semi-analytical model in the potential source areas located in the surrounding oceanic regions (Barents and Greenland Seas for Ny Ålesund, and Baffin Bay for Thule). Chl-a peaks in May in the Barents sea and in the Baffin Bay, and has maxima in June in the Greenland sea; PP follows the same seasonal pattern of Chl-a, although the differences in absolute values of PP in the three seas during the blooms are less marked than for Chl-a. MSA shows a better correlation with PP than with Chl-a, besides, the source intensity (expressed by PP) is able to explain more than 30% of the MSA variability at the two sites; the other factors explaining the MSA variability are taxonomic differences in the phytoplanktonic assemblages, and transport processes from the DMS source areas to the sampling sites. The taxonomic differences are also evident from the slopes of the correlation plots between MSA and PP: similar slopes (in the range 34.2–36.2 ng m−3of MSA/(gC m−2 d−1)) are found for the correlation between MSA at Ny Ålesund and PP in Barents Sea, and between MSA at Thule and PP in the Baffin Bay; conversely, the slope of the correlation between MSA at Ny Ålesund and PP in the Greenland Sea in summer is smaller (16.7 ng m−3of MSA/(gC m−2 d−1)). This is due to the fact that DMS emission from the Barents Sea and Baffin Bay is mainly related to the MIZ diatoms, which are prolific DMS producers, whereas in the Greenland Sea the DMS peak is related to an offshore pelagic bloom where low-DMS producer species are present. The sea ice dynamic plays a key role in determining MSA concentration in the Arctic, and a good correlation between MSA and SIM (slope = 39 ng m−3 of MSA/106 km2 SIM) and between MSA and IF-MIZ (slope = 56 ng m−3 of MSA/106 km2 IF-MIZ) is found for the cases attributable to bloomings of diatoms in the MIZ. Such relationships are calculated by combining the data sets from the two sites and suggest that PP is related to sea ice melting and to the extension of marginal sea ice areas, and that these factors are the main drivers for MSA concentrations at the considered Arctic sites.

Low-land permafrost areas are subject to intense freeze-thaw cycles and characterized by remarkable surface displacement. We used Sentinel-1 SAR interferometry (InSAR) in order to analyse the summer surface displacement over four spots in the Arctic and Antarctica since 2015. Choosing floodplain or outcrop areas as the reference for the InSAR relative deformation measurements, we found maximum subsidence of about 3 to 10 cm during the thawing season with generally high spatial variability. Sentinel-1 time-series of interferograms with 6–12 day time intervals highlight that subsidence is often occurring rather quickly within roughly one month in early summer. Intercomparison of summer subsidence from Sentinel-1 in 2017 with TerraSAR-X in 2013 over part of the Lena River Delta (Russia) shows a high spatial agreement between both SAR systems. A comparison with in-situ measurements for the summer of 2014 over the Lena River Delta indicates a pronounced downward movement of several centimetres in both cases but does not reveal a spatial correspondence between InSAR and local in-situ measurements. For the reconstruction of longer time-series of deformation, yearly Sentinel-1 interferograms from the end of the summer were considered. However, in order to infer an effective subsidence of the surface through melting of excess ice layers over multi-annual scales with Sentinel-1, a longer observation time period is necessary.

Abstract Freshwater runoff to fjords with marine-terminating glaciers along the Greenland Ice Sheet margin has an impact on fjord circulation and potentially ice sheet mass balance through increasing heat transport to the glacier front. Here, the authors use the high-resolution (5.5 km) HIRHAM5 regional climate model, allowing high detail in topography and surface types, to estimate freshwater input to Godthåbsfjord in southwest Greenland. Model output is compared to hydrometeorological observations and, while simulated daily variability in temperature and downwelling radiation shows high correlation with observations (typically >0.9), there are biases that impact the results. In particular, overestimated albedo leads to underestimation of melt and runoff at low elevations. In the model simulation (1991–2012), the ice sheet experiences increasing energy input from the surface turbulent heat flux (up to elevations of 2000 m) and shortwave radiation (at all elevations). Southerly wind anomalies and declining cloudiness due to an increase in atmospheric pressure over north Greenland contribute to increased summer melt. This results in declining surface mass balance (SMB), increasing surface runoff, and upward shift of the equilibrium line altitude. SMB is reconstructed back to 1890 though regression between simulated SMB and observed temperature and precipitation, with added uncertainty in the period 1890–1952 because of possible inhomogeneity in the precipitation record. SMB as low as in recent years appears to have occurred before, most notably around 1930, 1950, and 1960. While previous low SMBs were mainly caused by low accumulation, those around 1930 and in the 2000s are mainly due to warming.

The supply of freshwater to fjord systems in Greenland is increasing as a result of climate change-induced acceleration in ice sheet melt. However, insight into the marine implications of the melt water is impaired by lack of observations demonstrating the fate of freshwater along the Greenland coast and providing evaluation basis for ocean models. Here we present 13 years of summer measurements along a 120 km transect in Young Sound, Northeast Greenland and show that sub-surface coastal waters are decreasing in salinity with an average rate of 0.12 ± 0.05 per year. This is the first observational evidence of a significant freshening on decadal scale of the waters surrounding the ice sheet and comes from a region where ice sheet melt has been less significant. It implies that ice sheet dynamics in Northeast Greenland could be of key importance as freshwater is retained in southward flowing coastal currents thus reducing density of water masses influencing major deep water formation areas in the Subarctic Atlantic Ocean. Ultimately, the observed freshening could have implications for the Atlantic meridional overturning circulation.

Abstract Soils are warming as air temperatures rise across the Arctic and Boreal region concurrent with the expansion of tall-statured shrubs and trees in the tundra. Changes in vegetation structure and function are expected to alter soil thermal regimes, thereby modifying climate feedbacks related to permafrost thaw and carbon cycling. However, current understanding of vegetation impacts on soil temperature is limited to local or regional scales and lacks the generality necessary to predict soil warming and permafrost stability on a pan-Arctic scale. Here we synthesize shallow soil and air temperature observations with broad spatial and temporal coverage collected across 106 sites representing nine different vegetation types in the permafrost region. We showed ecosystems with tall-statured shrubs and trees (>40 cm) have warmer shallow soils than those with short-statured tundra vegetation when normalized to a constant air temperature. In tree and tall shrub vegetation types, cooler temperatures in the warm season do not lead to cooler mean annual soil temperature indicating that ground thermal regimes in the cold-season rather than the warm-season are most critical for predicting soil warming in ecosystems underlain by permafrost. Our results suggest that the expansion of tall shrubs and trees into tundra regions can amplify shallow soil warming, and could increase the potential for increased seasonal thaw depth and increase soil carbon cycling rates and lead to increased carbon dioxide loss and further permafrost thaw.

Abstract Firn-density variations have been studied in the lower accumulation area of the Greenland ice sheet (1440-1620 m a.s.l.) near Pâkitsoq, West Greenland. The main control on density in the near-surface firn layer (of 5-10 m thickness) is the formation of ice layers by the refreezing of meltwater that reaches depths of 2-4 m below the surface. The density variations are described by the ratio of annual surface melt M to the annual accumulation C . The ratio M/C is about 0.6 at the run-off limit (at about 1400 m a.s.l. in the study area) where refreezing of meltwater transforms snow into impermeable ice. The mean density of near-surface firn decreases with elevation, reflecting a decrease in melt with elevation. There is a surprising decrease in firn density at depths of more than about 4 m below the 1991 summer surface, which reflects lower melt rates and/or higher accumulation in the early 1980s and late 1970s when this firn was passing through the surface layer. The formation of such low-density firn may have partially contributed to the 1978-85 thickening of the ice sheet observed by satellite-radar altimetry. Near-surface firn density is therefore very sensitive to climate change and might be an attractive target for climate monitoring.

ABSTRACT Glacier surface mass-balance measurements on Greenland started more than a century ago, but no compilation exists of the observations from the ablation area of the ice sheet and local glaciers. Such data could be used in the evaluation of modelled surface mass balance, or to document changes in glacier melt independently from model output. Here, we present a comprehensive database of Greenland glacier surface mass-balance observations from the ablation area of the ice sheet and local glaciers. The database spans the 123 a from 1892 to 2015, contains a total of ~3000 measurements from 46 sites, and is openly accessible through the PROMICE web portal ( http://www.promice.dk ). For each measurement we provide X, Y and Z coordinates, starting and ending dates as well as quality flags. We give sources for each entry and for all metadata. Two thirds of the data were collected from grey literature and unpublished archive documents. Roughly 60% of the measurements were performed by the Geological Survey of Denmark and Greenland (GEUS, previously GGU). The data cover all regions of Greenland except for the southernmost part of the east coast, but also emphasize the importance of long-term time series of which there are only two exceeding 20 a. We use the data to analyse uncertainties in point measurements of surface mass balance, as well as to estimate surface mass-balance profiles for most regions of Greenland.

In the age of the Anthropocene, the ocean has typically been viewed as a sink for pollution. Pollution is varied, ranging from human-made plastics and pharmaceutical compounds, to human-altered abiotic factors, such as sediment and nutrient runoff. As global population, wealth and resource consumption continue to grow, so too does the amount of potential pollution produced. This presents us with a grand challenge which requires interdisciplinary knowledge to solve. There is sufficient data on the human health, social, economic, and environmental risks of marine pollution, resulting in increased awareness and motivation to address this global challenge, however a significant lag exists when implementing strategies to address this issue. This review draws upon the expertise of 17 experts from the fields of social sciences, marine science, visual arts, and Traditional and First Nations Knowledge Holders to present two futures; the Business-As-Usual, based on current trends and observations of growing marine pollution, and a More Sustainable Future, which imagines what our ocean could look like if we implemented current knowledge and technologies. We identify priority actions that governments, industry and consumers can implement at pollution sources, vectors and sinks, over the next decade to reduce marine pollution and steer us towards the More Sustainable Future. Supplementary Information: The online version contains supplementary material available at 10.1007/s11160-021-09674-8.