A 2-million-year-old ecosystem in Greenland uncovered by environmental DNA.

Late Pliocene and Early Pleistocene epochs 3.6 to 0.8 million years ago1 had climates resembling those forecasted under future warming2. Palaeoclimatic records show strong polar amplification with mean annual temperatures of 11-19 °C above contemporary values3,4. The biological communities inhabiting the Arctic during this time remain poorly known because fossils are rare5. Here we report an ancient environmental DNA6 (eDNA) record describing the rich plant and animal assemblages of the Kap København Formation in North Greenland, dated to around two million years ago. The record shows an open boreal forest ecosystem with mixed vegetation of poplar, birch and thuja trees, as well as a variety of Arctic and boreal shrubs and herbs, many of which had not previously been detected at the site from macrofossil and pollen records. The DNA record confirms the presence of hare and mitochondrial DNA from animals including mastodons, reindeer, rodents and geese, all ancestral to their present-day and late Pleistocene relatives. The presence of marine species including horseshoe crab and green algae support a warmer climate than today. The reconstructed ecosystem has no modern analogue. The survival of such ancient eDNA probably relates to its binding to mineral surfaces. Our findings open new areas of genetic research, demonstrating that it is possible to track the ecology and evolution of biological communities from two million years ago using ancient eDNA.

Extensive inland thinning and speed-up of Northeast Greenland Ice Stream.

Over the past two decades, ice loss from the Greenland ice sheet (GrIS) has increased owing to enhanced surface melting and ice discharge to the ocean1-5. Whether continuing increased ice loss will accelerate further, and by how much, remains contentious6-9. A main contributor to future ice loss is the Northeast Greenland Ice Stream (NEGIS), Greenland's largest basin and a prominent feature of fast-flowing ice that reaches the interior of the GrIS10-12. Owing to its topographic setting, this sector is vulnerable to rapid retreat, leading to unstable conditions similar to those in the marine-based setting of ice streams in Antarctica13-20. Here we show that extensive speed-up and thinning triggered by frontal changes in 2012 have already propagated more than 200 km inland. We use unique global navigation satellite system (GNSS) observations, combined with surface elevation changes and surface speeds obtained from satellite data, to select the correct basal conditions to be used in ice flow numerical models, which we then use for future simulations. Our model results indicate that this marine-based sector alone will contribute 13.5-15.5 mm sea-level rise by 2100 (equivalent to the contribution of the entire ice sheet over the past 50 years) and will cause precipitous changes in the coming century. This study shows that measurements of subtle changes in the ice speed and elevation inland help to constrain numerical models of the future mass balance and higher-end projections show better agreement with observations.

Hypothesis: Metformin is a potential reproductive toxicant.

Metformin is the first-line oral treatment for type 2 diabetes mellitus and is prescribed to more than 150 million people worldwide. Metformin's effect as a glucose-lowering drug is well documented but the precise mechanism of action is unknown. A recent finding of an association between paternal metformin treatment and increased numbers of genital birth defects in sons and a tendency towards a skewed secondary sex ratio with less male offspring prompted us to focus on other evidence of reproductive side effects of this drug. Metformin in humans is documented to reduce the circulating level of testosterone in both men and women. In experimental animal models, metformin exposure in utero induced sex-specific reproductive changes in adult rat male offspring with reduced fertility manifested as a 30% decrease in litter size and metformin exposure to fish, induced intersex documented in testicular tissue. Metformin is excreted unchanged into urine and feces and is present in wastewater and even in the effluent of wastewater treatment plants from where it spreads to rivers, lakes, and drinking water. It is documented to be present in numerous freshwater samples throughout the world - and even in drinking water. We here present the hypothesis that metformin needs to be considered a potential reproductive toxicant for humans, and probably also for wildlife. There is an urgent need for studies exploring the association between metformin exposure and reproductive outcomes in humans, experimental animals, and aquatic wildlife.

A Late Paleocene age for Greenland's Hiawatha impact structure.

The ~31-km-wide Hiawatha structure, located beneath Hiawatha Glacier in northwestern Greenland, has been proposed as an impact structure that may have formed after the Pleistocene inception of the Greenland Ice Sheet. To date the structure, we conducted 40Ar/39Ar analyses on glaciofluvial sand and U-Pb analyses on zircon separated from glaciofluvial pebbles of impact melt rock, all sampled immediately downstream of Hiawatha Glacier. Unshocked zircon in the impact melt rocks dates to ~1915 million years (Ma), consistent with felsic intrusions found in local bedrock. The 40Ar/39Ar data indicate Late Paleocene resetting and shocked zircon dates to 57.99 ± 0.54 Ma, which we interpret as the impact age. Consequently, the Hiawatha impact structure far predates Pleistocene glaciation and is unrelated to either the Paleocene-Eocene Thermal Maximum or flood basalt volcanism in east Greenland. However, it was contemporaneous with the Paleocene Carbon Isotope Maximum, although the impact's exact paleoenvironmental and climatic significance awaits further investigation.

Late Quaternary dynamics of Arctic biota from ancient environmental genomics.

During the last glacial-interglacial cycle, Arctic biotas experienced substantial climatic changes, yet the nature, extent and rate of their responses are not fully understood1-8. Here we report a large-scale environmental DNA metagenomic study of ancient plant and mammal communities, analysing 535 permafrost and lake sediment samples from across the Arctic spanning the past 50,000 years. Furthermore, we present 1,541 contemporary plant genome assemblies that were generated as reference sequences. Our study provides several insights into the long-term dynamics of the Arctic biota at the circumpolar and regional scales. Our key findings include: (1) a relatively homogeneous steppe-tundra flora dominated the Arctic during the Last Glacial Maximum, followed by regional divergence of vegetation during the Holocene epoch; (2) certain grazing animals consistently co-occurred in space and time; (3) humans appear to have been a minor factor in driving animal distributions; (4) higher effective precipitation, as well as an increase in the proportion of wetland plants, show negative effects on animal diversity; (5) the persistence of the steppe-tundra vegetation in northern Siberia enabled the late survival of several now-extinct megafauna species, including the woolly mammoth until 3.9 ± 0.2 thousand years ago (ka) and the woolly rhinoceros until 9.8 ± 0.2 ka; and (6) phylogenetic analysis of mammoth environmental DNA reveals a previously unsampled mitochondrial lineage. Our findings highlight the power of ancient environmental metagenomics analyses to advance understanding of population histories and long-term ecological dynamics.

Centennial response of Greenland's three largest outlet glaciers.

The Greenland Ice Sheet is the largest land ice contributor to sea level rise. This will continue in the future but at an uncertain rate and observational estimates are limited to the last few decades. Understanding the long-term glacier response to external forcing is key to improving projections. Here we use historical photographs to calculate ice loss from 1880-2012 for Jakobshavn, Helheim, and Kangerlussuaq glacier. We estimate ice loss corresponding to a sea level rise of 8.1 ± 1.1 millimetres from these three glaciers. Projections of mass loss for these glaciers, using the worst-case scenario, Representative Concentration Pathways 8.5, suggest a sea level contribution of 9.1-14.9 mm by 2100. RCP8.5 implies an additional global temperature increase of 3.7 °C by 2100, approximately four times larger than that which has taken place since 1880. We infer that projections forced by RCP8.5 underestimate glacier mass loss which could exceed this worst-case scenario.

A large impact crater beneath Hiawatha Glacier in northwest Greenland.

We report the discovery of a large impact crater beneath Hiawatha Glacier in northwest Greenland. From airborne radar surveys, we identify a 31-kilometer-wide, circular bedrock depression beneath up to a kilometer of ice. This depression has an elevated rim that cross-cuts tributary subglacial channels and a subdued central uplift that appears to be actively eroding. From ground investigations of the deglaciated foreland, we identify overprinted structures within Precambrian bedrock along the ice margin that strike tangent to the subglacial rim. Glaciofluvial sediment from the largest river draining the crater contains shocked quartz and other impact-related grains. Geochemical analysis of this sediment indicates that the impactor was a fractionated iron asteroid, which must have been more than a kilometer wide to produce the identified crater. Radiostratigraphy of the ice in the crater shows that the Holocene ice is continuous and conformable, but all deeper and older ice appears to be debris rich or heavily disturbed. The age of this impact crater is presently unknown, but from our geological and geophysical evidence, we conclude that it is unlikely to predate the Pleistocene inception of the Greenland Ice Sheet.

Parasitic infections and resource economy of Danish Iron Age settlement through ancient DNA sequencing.

In this study, we screen archaeological soil samples by microscopy and analyse the samples by next generation sequencing to obtain results with parasites at species level and untargeted findings of plant and animal DNA. Three separate sediment layers of an ancient man-made pond in Hoby, Denmark, ranging from 100 BC to 200 AD, were analysed by microscopy for presence of intestinal worm eggs and DNA analysis were performed to identify intestinal worms and dietary components. Ancient DNA of parasites, domestic animals and edible plants revealed a change in use of the pond over time reflecting the household practice in the adjacent Iron Age settlement. The most abundant parasite found belonged to the Ascaris genus, which was not possible to type at species level. For all sediment layers the presence of eggs of the human whipworm Trichuris trichiura and the beef tapeworm Taenia saginata suggests continuous disposal of human faeces in the pond. Moreover, the continuous findings of T. saginata further imply beef consumption and may suggest that cattle were living in the immediate surrounding of the site throughout the period. Findings of additional host-specific parasites suggest fluctuating presence of other domestic animals over time: Trichuris suis (pig), Parascaris univalens (horse), Taenia hydatigena (dog and sheep). Likewise, alternating occurrence of aDNA of edible plants may suggest changes in agricultural practices. Moreover, the composition of aDNA of parasites, plants and vertebrates suggests a significant change in the use of the ancient pond over a period of three centuries.

DNA evidence of bowhead whale exploitation by Greenlandic Paleo-Inuit 4,000 years ago.

The demographic history of Greenland is characterized by recurrent migrations and extinctions since the first humans arrived 4,500 years ago. Our current understanding of these extinct cultures relies primarily on preserved fossils found in their archaeological deposits, which hold valuable information on past subsistence practices. However, some exploited taxa, though economically important, comprise only a small fraction of these sub-fossil assemblages. Here we reconstruct a comprehensive record of past subsistence economies in Greenland by sequencing ancient DNA from four well-described midden deposits. Our results confirm that the species found in the fossil record, like harp seal and ringed seal, were a vital part of Inuit subsistence, but also add a new dimension with evidence that caribou, walrus and whale species played a more prominent role for the survival of Paleo-Inuit cultures than previously reported. Most notably, we report evidence of bowhead whale exploitation by the Saqqaq culture 4,000 years ago.

Geodetic measurements reveal similarities between post-Last Glacial Maximum and present-day mass loss from the Greenland ice sheet.

Accurate quantification of the millennial-scale mass balance of the Greenland ice sheet (GrIS) and its contribution to global sea-level rise remain challenging because of sparse in situ observations in key regions. Glacial isostatic adjustment (GIA) is the ongoing response of the solid Earth to ice and ocean load changes occurring since the Last Glacial Maximum (LGM; ~21 thousand years ago) and may be used to constrain the GrIS deglaciation history. We use data from the Greenland Global Positioning System network to directly measure GIA and estimate basin-wide mass changes since the LGM. Unpredicted, large GIA uplift rates of +12 mm/year are found in southeast Greenland. These rates are due to low upper mantle viscosity in the region, from when Greenland passed over the Iceland hot spot about 40 million years ago. This region of concentrated soft rheology has a profound influence on reconstructing the deglaciation history of Greenland. We reevaluate the evolution of the GrIS since LGM and obtain a loss of 1.5-m sea-level equivalent from the northwest and southeast. These same sectors are dominating modern mass loss. We suggest that the present destabilization of these marine-based sectors may increase sea level for centuries to come. Our new deglaciation history and GIA uplift estimates suggest that studies that use the Gravity Recovery and Climate Experiment satellite mission to infer present-day changes in the GrIS may have erroneously corrected for GIA and underestimated the mass loss by about 20 gigatons/year.

Postglacial viability and colonization in North America's ice-free corridor.

During the Last Glacial Maximum, continental ice sheets isolated Beringia (northeast Siberia and northwest North America) from unglaciated North America. By around 15 to 14 thousand calibrated radiocarbon years before present (cal. kyr bp), glacial retreat opened an approximately 1,500-km-long corridor between the ice sheets. It remains unclear when plants and animals colonized this corridor and it became biologically viable for human migration. We obtained radiocarbon dates, pollen, macrofossils and metagenomic DNA from lake sediment cores in a bottleneck portion of the corridor. We find evidence of steppe vegetation, bison and mammoth by approximately 12.6 cal. kyr bp, followed by open forest, with evidence of moose and elk at about 11.5 cal. kyr bp, and boreal forest approximately 10 cal. kyr bp. Our findings reveal that the first Americans, whether Clovis or earlier groups in unglaciated North America before 12.6 cal. kyr bp, are unlikely to have travelled by this route into the Americas. However, later groups may have used this north-south passageway.

Digital elevation model and orthophotographs of Greenland based on aerial photographs from 1978-1987.

Digital Elevation Models (DEMs) play a prominent role in glaciological studies for the mass balance of glaciers and ice sheets. By providing a time snapshot of glacier geometry, DEMs are crucial for most glacier evolution modelling studies, but are also important for cryospheric modelling in general. We present a historical medium-resolution DEM and orthophotographs that consistently cover the entire surroundings and margins of the Greenland Ice Sheet 1978-1987. About 3,500 aerial photographs of Greenland are combined with field surveyed geodetic ground control to produce a 25 m gridded DEM and a 2 m black-and-white digital orthophotograph. Supporting data consist of a reliability mask and a photo footprint coverage with recording dates. Through one internal and two external validation tests, this DEM shows an accuracy better than 10 m horizontally and 6 m vertically while the precision is better than 4 m. This dataset proved successful for topographical mapping and geodetic mass balance. Other uses include control and calibration of remotely sensed data such as imagery or InSAR velocity maps.

Holocene ice marginal fluctuations of the Qassimiut lobe in South Greenland.

Knowledge about the Holocene evolution of the Greenland ice sheet (GrIS) is important to put the recent observations of ice loss into a longer-term perspective. In this study, we use six new threshold lake records supplemented with two existing lake records to reconstruct the Holocene ice marginal fluctuations of the Qassimiut lobe (QL) - one of the most dynamic parts of the GrIS in South Greenland. Times when the ice margin was close to present extent are characterized by clastic input from the glacier meltwater, whereas periods when the ice margin was behind its present day extent comprise organic-rich sediments. We find that the overall pattern suggests that the central part of the ice lobe in low-lying areas experienced the most prolonged ice retreat from ~9-0.4 cal. ka BP, whereas the more distal parts of the ice lobe at higher elevation re-advanced and remained close to the present extent during the Neoglacial between ~4.4 and 1.8 cal. ka BP. These results demonstrate that the QL was primarily driven by Holocene climate changes, but also emphasises the role of local topography on the ice marginal fluctuations.

Spatial and temporal distribution of mass loss from the Greenland Ice Sheet since AD 1900.

The response of the Greenland Ice Sheet (GIS) to changes in temperature during the twentieth century remains contentious, largely owing to difficulties in estimating the spatial and temporal distribution of ice mass changes before 1992, when Greenland-wide observations first became available. The only previous estimates of change during the twentieth century are based on empirical modelling and energy balance modelling. Consequently, no observation-based estimates of the contribution from the GIS to the global-mean sea level budget before 1990 are included in the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Here we calculate spatial ice mass loss around the entire GIS from 1900 to the present using aerial imagery from the 1980s. This allows accurate high-resolution mapping of geomorphic features related to the maximum extent of the GIS during the Little Ice Age at the end of the nineteenth century. We estimate the total ice mass loss and its spatial distribution for three periods: 1900-1983 (75.1 ± 29.4 gigatonnes per year), 1983-2003 (73.8 ± 40.5 gigatonnes per year), and 2003-2010 (186.4 ± 18.9 gigatonnes per year). Furthermore, using two surface mass balance models we partition the mass balance into a term for surface mass balance (that is, total precipitation minus total sublimation minus runoff) and a dynamic term. We find that many areas currently undergoing change are identical to those that experienced considerable thinning throughout the twentieth century. We also reveal that the surface mass balance term shows a considerable decrease since 2003, whereas the dynamic term is constant over the past 110 years. Overall, our observation-based findings show that during the twentieth century the GIS contributed at least 25.0 ± 9.4 millimetres of global-mean sea level rise. Our result will help to close the twentieth-century sea level budget, which remains crucial for evaluating the reliability of models used to predict global sea level rise.

Greenland ice sheet mass balance: a review.

Over the past quarter of a century the Arctic has warmed more than any other region on Earth, causing a profound impact on the Greenland ice sheet (GrIS) and its contribution to the rise in global sea level. The loss of ice can be partitioned into processes related to surface mass balance and to ice discharge, which are forced by internal or external (atmospheric/oceanic/basal) fluctuations. Regardless of the measurement method, observations over the last two decades show an increase in ice loss rate, associated with speeding up of glaciers and enhanced melting. However, both ice discharge and melt-induced mass losses exhibit rapid short-term fluctuations that, when extrapolated into the future, could yield erroneous long-term trends. In this paper we review the GrIS mass loss over more than a century by combining satellite altimetry, airborne altimetry, interferometry, aerial photographs and gravimetry data sets together with modelling studies. We revisit the mass loss of different sectors and show that they manifest quite different sensitivities to atmospheric and oceanic forcing. In addition, we discuss recent progress in constructing coupled ice-ocean-atmosphere models required to project realistic future sea-level changes.

Fifty thousand years of Arctic vegetation and megafaunal diet.

Although it is generally agreed that the Arctic flora is among the youngest and least diverse on Earth, the processes that shaped it are poorly understood. Here we present 50 thousand years (kyr) of Arctic vegetation history, derived from the first large-scale ancient DNA metabarcoding study of circumpolar plant diversity. For this interval we also explore nematode diversity as a proxy for modelling vegetation cover and soil quality, and diets of herbivorous megafaunal mammals, many of which became extinct around 10 kyr bp (before present). For much of the period investigated, Arctic vegetation consisted of dry steppe-tundra dominated by forbs (non-graminoid herbaceous vascular plants). During the Last Glacial Maximum (25-15 kyr bp), diversity declined markedly, although forbs remained dominant. Much changed after 10 kyr bp, with the appearance of moist tundra dominated by woody plants and graminoids. Our analyses indicate that both graminoids and forbs would have featured in megafaunal diets. As such, our findings question the predominance of a Late Quaternary graminoid-dominated Arctic mammoth steppe.

Molecular- and pollen-based vegetation analysis in lake sediments from central Scandinavia.

Plant and animal biodiversity can be studied by obtaining DNA directly from the environment. This new approach in combination with the use of generic barcoding primers (metabarcoding) has been suggested as complementary or alternative to traditional biodiversity monitoring in ancient soil sediments. However, the extent to which metabarcoding truly reflects plant composition remains unclear, as does its power to identify species with no pollen or macrofossil evidence. Here, we compared pollen-based and metabarcoding approaches to explore the Holocene plant composition around two lakes in central Scandinavia. At one site, we also compared barcoding results with those obtained in earlier studies with species-specific primers. The pollen analyses revealed a larger number of taxa (46), of which the majority (78%) was not identified by metabarcoding. The metabarcoding identified 14 taxa (MTUs), but allowed identification to a lower taxonomical level. The combined analyses identified 52 taxa. The barcoding primers may favour amplification of certain taxa, as they did not detect taxa previously identified with species-specific primers. Taphonomy and selectiveness of the primers are likely the major factors influencing these results. We conclude that metabarcoding from lake sediments provides a complementary, but not an alternative, tool to pollen analysis for investigating past flora. In the absence of other fossil evidence, metabarcoding gives a local and important signal from the vegetation, but the resulting assemblages show limited capacity to detect all taxa, regardless of their abundance around the lake. We suggest that metabarcoding is followed by pollen analysis and the use of species-specific primers to provide the most comprehensive signal from the environment.