Last Glacial Research Articles

Chronology and erosion rate of the Pinedale glaciation, Colorado Front Range (USA), inferred from the sedimentary record of glacial Lake Devlin

Glacial and periglacial sediments and landforms record the chronology of glaciation and amount of Pleistocene erosion during colder periods that added substantially to global sediment budgets and contributed to the global CO2 cycle. The now-drained glacial Lake Devlin, dammed in a Front Range tributary valley by a glacier in the North Branch of Boulder Creek (Colorado, USA) preserves an important sedimentary archive of the ca. 32−14 ka Pinedale glaciation, recording both paleoclimate information and an integrated measure of glacial and periglacial erosion rates over a full glacial cycle. Despite rapid erosion of fine-grained deposits after the lake drained, most sediment generated during Pinedale time remains as legacy deposits in the catchment. Geomorphic evidence and dating of glaciolacustrine sediment from surface exposures demonstrate that the ca. 30 ka Pinedale glacial advance was nearly as extensive as the local Late Glacial Maximum at ca. 20 ka. Sedimentary archives dated by 14C, optically stimulated luminescence, and cosmogenic nuclides extend earlier studies (Madole et al., 1973) of pollen and magnetic susceptibility (MS) in cores from the glaciolacustrine deposits of Lake Devlin and of Pinedale climate. Records suggest short-term warming and biotic change at ca. 15 ka after ∼14 kyr of cold, dry conditions punctuated by MS peaks at ca. 26.5 ka, 20 ka, and 16.5 ka. Lake Devlin drained catastrophically after ca. 14 ka, millennia after ice had retreated upvalley from the lateral moraine that dammed the lake. Sediment production during the Pinedale was equivalent to a periglacial and glacial erosion rate of ∼70 mm kyr−1, several times higher than long-term rates in the adjacent Front Range, but much lower than rates measured where modern glaciers are eroding weak bedrock in zones of rapid rock uplift, such as SSE Alaska, USA. Data from the Lake Devlin basin contribute to contemporary discussions of how glacial erosion influences the global CO2 cycle.

High-northern-latitude forcing modulated fluvial‒aeolian interplay during the Last Glacial Period in the Mu Us Desert, northern central China

Substantial stratigraphic evidence from the Late Quaternary at orbital- and millennial‒scales indicates that the formation and evolution of aeolian‒sand landscapes results from the interaction between regional aeolian and fluvial‒lacustrine activity. In many deserts (sandy lands) across the planet, fluvial‒aeolian interplay at the glacial‒interglacial scale since the Late Quaternary has been responsive to monsoonal precipitation, and linked to atmospheric CO2 and periodic variations in insolation. However, the patterns of fluvial‒aeolian interplay and their relationship with climate change are not always consistent over a given timescale during the Late Pleistocene. The paucity of studies on fluvial‒aeolian interplay in the Mu Us Desert (MUD) during the Last Glacial Period (LGP) hinders our understanding of the relationships between the dynamics of aeolian and fluvial‒lacustrine landscapes and their drivers. Here, based on a comprehensive field survey and stratigraphic dating of a typical fluvial‒aeolian sequence in the MUD, we synthesized the chronology of fluvial‒lacustrine and aeolian strata during the LGP. Our results show that extensive development of a fluvial‒lacustrine environment in the MUD during the LG mainly occurred during the late Last Glacial Maximum (LGMa, 20.1–18.4 ka BP) and Bølling/Allerød interstadial (14.8–12.6 ka BP), while thicker aeolian sands extensively accumulated during the Heinrich Stadial 2, Heinrich Stadial 1, and Younger Dryas periods. However, the development of fluvial‒lacustrine environments showed considerable spatial heterogeneity, mainly occurring on the southern and eastern margins of the desert, while they were extremely rare in the central and western parts of the desert owing to strong wind erosion. Comparisons with other high-resolution geological records from the Asian monsoon region reveal that the fluvial‒aeolian interplay in the MUD during the LGP was highly consistent with the evolution pattern of the East Asian summer monsoon (EASM) on the millennial‒centennial scale, and indicate that the development of fluvial‒lacustrine environments was mainly controlled by an increase in monsoonal precipitation during strengthening of the EASM. Consistency with other broad changes at high latitudes in the Northern Hemisphere suggests that fluvial‒aeolian interplay during the LGP in northern Chinese deserts, as represented by the MUD, was mainly controlled by changes in high-northern-latitude ice volume and the Atlantic Meridional Overturning Circulation.

Paleohydrogeology of the Horonobe area, Northern Hokkaido, Japan: Groundwater flow conditions during glacial and postglacial periods estimated from chemical and isotopic data for fracture and pore water

Understanding the difference in groundwater flow between glacial and interglacial periods is crucial for predicting the impact of future climate changes on groundwater movement. This study assesses the difference in groundwater flow between the last glacial period (LGP) and the postglacial period (PGP) in fractured mudstones of the Horonobe area, Japan, by combining the data for stable isotopes (δD and δ18O) and Cl− concentration of fracture and pore waters with radiocarbon (14C) age. The isotopic compositions of fractures and pore waters indicate that groundwater at 28–250 m deep in a borehole closest to the recharge area comprises meteoric water, recharged under the same climates as the present. The fracture water has isotopic compositions more similar to meteoric water than the matrix pore water near the fracture. The 14C age of fracture water suggests meteoric water recharge during the PGP. At greater depths in the borehole and sampling points in other boreholes, the isotopic compositions indicate the mixing of glacial meteoric and altered connate water, with the fracture water having comparable isotopic compositions with the matrix pore water. The recharge timing of meteoric water is inferred to be the LGP or before based on 14C dating. These results suggest that the meteoric water recharged during the PGP flows at a shallow depth, whereas the meteoric water recharged during the LGP intruded to greater depths. This result is consistent with previous inferences from surface geophysical and geological surveys that the depths of local valleys during the LGP were greater by < 50 m than the present ones and enhanced the downward hydraulic gradient. Combining the chemical and isotopic compositions of groundwater with 14C age helps assess the groundwater flow during the LGP and PGP in fractured rocks.

Global clustering of recent glacier surges from radar backscatter data, 2017–2022

AbstractUsing global Sentinel-1 radar backscatter data, we systematically map the locations of glaciers with surge-type activity during 2017–22. Patterns of pronounced increases or decreases in the strongest backscatter between two winter seasons often indicate large changes in glacier crevassing, which we treat here as a sign of surge-type activity. Validations against velocity time series, terminus advances and crevassing found in optical satellite images confirm the robustness of this approach. We find 115 surge-type events globally between 2017 and 2022, around 100 of which on glaciers already know as surge-type. Our data reveal a pronounced spatial clustering in three regions, (i) Karakoram, Pamirs and Western Kunlun Shan (~50 surges), (ii) Svalbard (~25) and (iii) Yukon/Alaska (~9), with only a few other scattered surges elsewhere. This spatial clustering is significantly more pronounced than the overall global clustering of known surge-type glaciers. The 2017–22 clustering may point to climatic forcing of surge initiation.

Unraveling the climate control on debris-free glacier evolution in the Everest region (Nepal, central Himalaya) during the Holocene

Current mass balance and meteorological surveys of Mera glacier located about 30 km south of Mount Everest in Nepal show the dominant role of Asian monsoon precipitation on interannual mass balance variability while temperature controls the altitude of snow-rain threshold. As these observations on mass balance variability only explore the recent decades, studies on paleo glacial extents are useful to investigate the long-term climate forcing on glacier evolution. To do so, we investigated the long-term evolution of the debris-free Mera glacier and a neighbouring small debris-free South Khare glacier. Fifty-one 10Be CRE ages were obtained from samples collected on moraine boulders and roches moutonnées. 10Be CRE ages of the boulders span from the end of the Lateglacial (19.0–11.7 ka) to the Little Ice Age (∼0.6–0.1 ka). The oldest dated moraine in this study was observed at the base of South Khare glacier with an age of 13.6 ± 0.5 ka. The two glaciers subsequently experienced their largest Holocene extent in the Early Holocene with moraines dated to 11.0 ± 0.3 ka at the base of Mera glacier and 10.8 ± 0.5 ka at the base of South Khare glacier. We did not observe any moraine from the Mid-Holocene. During the Late Holocene several glacier advances were recorded around 2.3 ± 0.2 ka, 1.5 ka and then during the last centuries at Mera glacier and around 2.8 ± 0.6 ka, and during the Little Ice Age at South Khare glacier. To explore the links between long-term Nepalese glacier changes and climate, we used oceanic and terrestrial Indian Summer monsoon reconstructions and temperature and precipitation output from two transient global climate models TraCE and LOVECLIM. These climate data outputs were corrected by a reconstruction of the Atlantic Meridional Overturning Circulation (AMOC) over the Holocene and its associated climatic impacts. We also used sensitivity experiments from the IPSL (Institut Pierre Simon Laplace) model to discuss the possible influence of horizontal resolution, land hydrology, vegetation and runoff on changes in Asian summer monsoon. Importantly, we show this long-term Nepalese glacier pattern does not perfectly conform neither to the Indian monsoon precipitation that is documented from terrestrial and marine records nor to temperature and precipitation changes simulated by the models. While the maximum glacier extent in the Early Holocene corresponds to enhanced precipitation documented by proxies and models, the Late Holocene glacier advance remains puzzling. We claim that new paleo glacier records and improved climate simulations are necessary to get a better understanding of past glacier changes and the associated climate dynamics, which might be crucial to gain confidence in both glacier and climate future evolutions.

Controls on fracture openness and reactivation in Forsmark, Sweden

In crystalline bedrock, the open fraction of the fracture network constitutes the main pathways for fluids. Many observations point out that the state of stress influences the open fraction, likely indicating recent reactivation. But how this occurs is still unresolved. We analyse the conditions for fracture reactivation from fracture data collected in the uppermost 1 km of bedrock in Forsmark, Sweden. The open fraction is mainly correlated to the stress acting normally on the fracture but even away from critical failure, leading us to analyse the potential fluid pressure required for reactivation, {P}_{c}. We observe that 100% of the fractures are open when {P}_{c} is hydrostatic, and the ratio decreases exponentially to a plateau of ~ 17% when {P}_{c} is lithostatic and above. Exceptions are the oldest fractures, having a low open fraction independent of {P}_{c}. We suggest that these results reflect past pressure build-ups, potentially related to recent glaciations, and developing only if the preexisting open fraction is large enough.

The Yale Geochronometric Laboratory and the Rewriting of Global Environmental History.

Beginning in the nineteenth century, scientists speculated that the Pleistocene megafauna-species such as the giant ground sloth, wooly mammoth, and saber-tooth cat-perished because of rapid climate change accompanying the end of the most recent Ice Age. In the 1950s, a small network of ecologists challenged this view in collaboration with archeologists who used the new tool of radiocarbon dating. The Pleistocene overkill hypothesis imagined human hunting, not climate change, to be the primary cause of megafaunal extinction. This article situates the Pleistocene overkill hypothesis in a broader history of the emergence of historical ecology as a distinct sub-discipline of paleoecology. Tracing the work of the Yale Geochronometric Laboratory and an interdisciplinary research network that included Paul Sears, Richard Foster Flint, Edward Deevey, Kathryn Clisby, and Paul S. Martin, it reveals how both the methods and the meaning of studying fossil pollen shifted between the 1910s and 1960s. First used as a tool for fossil fuel extraction, fossil pollen became a means of envisioning climatic history, and ultimately, a means of reimagining global ecological history. First through pollen stratigraphy and then through radiocarbon dating, ecologists reconstructed past biotic communities and rethought the role of humans in these communities. By the 1980s, the discipline of historical ecology would reshape physical environments through the practice of ecological restoration.

Characterizing 4 decades of accelerated glacial mass loss in the west Nyainqentanglha Range of the Tibetan Plateau

Abstract. Glacier retreat is altering the water regime of the Tibetan Plateau (TP) as the region's climate changes, but there remain substantial gaps in our knowledge of recent glacier loss in this region due to the difficulty of making direct high-mountain observations, and this limits our ability to predict the future of this important water resource. Here, we assessed 44 years of glacier area and volume changes in the major west Nyainqentanglha Range (WNR) that supplies meltwater to the densely populated Lhasa River basin and Nam Co, the second largest endorheic lake on the TP. Between the two periods 1976–2000 and 2000–2020, we found that the glacier areal retreat rate more than doubled (from -0.54±0.21 % a−1 to -1.17±0.30 % a−1), and surface lowering also accelerated (from -0.26±0.09 to -0.37±0.15 m w.e a−1) with particularly intense melting after 2014. This acceleration is similar in both timing and magnitude to that observed for Himalayan glaciers farther south. Besides, the areal retreat rate and mass loss rate of most glaciers in the WNR were not synchronized. To understand the sensitivity of WNR glaciers to climate forcing, we examined the effects of topography, debris cover and the presence of proglacial lakes on our observed changes. We found consistently faster areal retreat but slower thinning rates on steeper slopes and an inconsistent relationship with aspect. We concluded that our observed spatial and temporal patterns of glacier change were dominated by observed local variations in temperature and precipitation, the melt-reducing role of supraglacial debris, and the increasing influence of ice-marginal lakes on glacier ablation.

Recent glacier and lake elevation changes in the western Kunlun Mountains (Tibetan Plateau) revealed by CryoSat-2 satellite altimetry data

Mountain glaciers and alpine lakes are significant indicators of regional climate change. The temperature rise has resulted in the retreat of glaciers, further influencing the ecological cycle and the water balance. The present work mainly focuses on recent glacier and lake elevation changes in the western Kunlun Mountains. The CryoSat-2 satellite altimetry with relatively dense ground track spacing can measure more elevation samples for small glaciers and lakes surrounded by complex topography. The CryoSat-2 data with various algorithms were employed to monitor the water levels of the three typical lakes and the surface elevation changes of their supplying glaciers. The glacier surface elevations showed a slightly increasing trend and remarkable spatial heterogeneity from 2010 to 2021. Meanwhile, the Bangda Lake and Aksayqin Lake displayed clear and strong increasing trends over the observation period with an annual change rate of +0.65 and +0.20 m / yr, respectively. Surprisingly, Gozha Lake exhibited a decreasing trend with the rate of −0.04 m / yr. Regional climate warming led to more melt-water from the glaciers and permafrost flowing into the three typical lakes. The increase in precipitation played a leading role in lake expansion over the western Kunlun Mountains. The present work demonstrated that the satellite altimetry data provide important supplementary data to help us perform the cryosphere research in the data-scarce northern Tibetan Plateau depopulated zone.

Geophysical–geotechnical methodology for assessing the spatial distribution of glacio‐lacustrine sediments: The case history of Lake Seracchi

AbstractProglacial lakes are distinctive features of deglaciated landscapes and often act as sediment sinks, collecting solid material from subglacial erosion or washout of deglaciated areas. The solid transport flow, strongly linked to the glaciers and periglacial landforms, may rise due to the rapid changes driven by climate warming, causing deep transformations in the basin hydrology, and even the appearance or disappearance of lakes at a decadal timescale.The goal of this study was to present a geophysical–geotechnical approach that integrates several techniques, to quantify the sediment distribution in a proglacial lake. A geophysical survey is performed with ground‐penetrating radar (GPR) installed on a boat, whereas a time‐domain reflectometer (TDR) measures the electrical conductivity and permittivity of the lakebed sediments. Unperturbed samples are collected and analyzed to measure the main geotechnical properties of the sediment: grain‐size distribution, plastic limit, and liquid limit. Such properties support the interpretation of the GPR data and the detection of spatial variations of the sediment facies.To validate the proposed methodology, field tests were carried out at Lake Seracchi, the largest lake of the Rutor glacier, Italian Alps. It formed around 1880 because of the recent glacier shrinkage, as chronicled by valuable historical documents. Its greyish waters carry a significant amount of suspended sediment recognized as glacial flour, which gradually accumulates on the bottom of the lake.The obtained bathymetry and sediment thickness maps of Lake Seracchi show the strength of the approach: from only a few manual samples, it is possible to extrapolate the geotechnical properties of interest, such as friction angle or hydraulic conductivity, to wider areas, surveyed by the geophysical techniques. This is achieved by investigating the spatial distribution of key geophysical properties linked to the geotechnical properties of interest.

Late Holocene glacier variations indicated by the δ18O of ice core enclosed gaseous oxygen in the central Tibetan Plateau

Late Holocene glacier variations indicated by the δ18O of ice core enclosed gaseous oxygen in the central Tibetan Plateau

Distribution of Devensian glacial erratics and related evidence elucidate complex ice flow changes across a former ice divide: Northern England

The Lune Gorge and the uplands of the southern Shap Fells represent a key area in developing an understanding of the dynamics of the Late Devensian glaciation (Dimlington Stadial) of northern England. Here ice masses emanating from southern Scotland, the Lake District and the Howgill Fells interacted in the area of the upper Lune valley. Glacial landforms are recorded and tills noted. The dispersal pattern of Shap granite (Sg) erratics is mapped as these clasts are an important tracer of regional ice movements and local ice dynamics. This new information is synthesised with existing literature to provide an understanding of ice dynamics in an area of complex ice flow history. In particular, the ice flow interactions over an area of the western Pennines and the southern Shap Fells have been defined.Early Dimlington ice flow in the vicinity of the upper Lune valley was easterly. Later a northerly shift in the position of a regionally significant west–east aligned ice divide led to topographic steering of southerly basal ice flows, resulting in the development of a western ice stream (Mint Sg plume) and an eastern ice stream (Lune Sg plume); both flowing south around the massif of the southern Shap Fells. At that time, southerly flow of basal ice over the highest ridges of the Grayrigg massif in the southern Shap Fells was relatively weak. Whereas the western stream extended into the Lancashire lowlands, the eastern ice stream was impeded in the Lune Gorge by ice emanating from Borrowdale which forced northern ice to the eastern side of the Lune Gorge where it was blocked by western flowing Howgill ice; the latter extending as far as just east of Kendal.During the late Dimlington, the ice masses over the Mint valley and the southern Shap Fells largely thinned and retreated to the north and west, with a surge in northerly ice movement within the upper Lune valley that did not override the western Pennines. Local ice dynamics are well-illustrated in the western margin of the Pennines (Crosby Ravensworth Fell Gaythorne Plain), where the disposition of erratics reflects local late west-to-east weak down-wasting ice flow and the presence of subglacial meltwater drainage channels. However, ice ultimately decayed in situ on Crosby Ravensworth Fell Gaythorne Plain. Similarly, after complex variable easterly and southerly ice flow during the early Dimlington, there was weak northerly ice flow in the later phase over Birkbeck Fells Common before ice thinned and retreated from Grayrigg Forest. An extensive ice stream, fed by ice emanating from an ice dispersal centre in the eastern Lake District, occupied the large trough of Borrowdale that transverses the southern Shap Fells, but its extension and recession dynamics remain enigmatic.“Glacial Theory …, in its application to the transport of blocks across Stainmoor, involves such obvious mechanical absurdities that the author considers it totally unworthy of the attention of the Society.” William Hopkins, 1842

Natural and anthropogenic halite karst subsidence in north Cheshire, UK; comparison of Rostherne Mere, Melchett Mere, Tatton Mere and their surroundings

Most of the north Cheshire (Knutsford Group) of meres (lakes) in the UK formed naturally by dissolution of Triassic halite after the Devensian glaciation. Anthropogenic brine extraction in the 19th and 20th centuries produced further subsidence that enlarged some meres and formed the new lake of Melchett Mere. The characteristic features of three meres, Rostherne, Melchett and Tatton, are compared here by historical surveys, maps, photographs and light detection and ranging (LiDAR) interpretations. These illustrate the similarities of the natural and anthropogenic subsidence features, which can be separated only by temporal evidence of their formation. Rostherne Mere and Tatton Mere are mainly natural, though deepened or made larger by anthropogenic salt dissolution; Melchett Mere is completely anthropogenic and mainly formed between 1927 and 2003. All three meres are surrounded by landslip scars related to the subsidence. Former brine pumping at Northwich, Plumley and possibly Agden is implicated in the formation of Melchett Mere and the reactivation of natural subsidence at Rostherne and Tatton meres plus The Mere along with Tabley, Pickmere and Budworth meres to the SW. The brine run linkages between these abstraction areas and the subsidence cross the route of the proposed HS2 railway.

A search for piroplasmids and spirochetes in threatened pudu (Pudu puda) and associated ticks from Southern Chile unveils a novel Babesia sp. and a variant of Borrelia chilensis.

Cervids are important hosts for ticks and although they are refractory to some tick-borne agents such as Borrelia, they do act as reservoirs for others such as Babesia. Babesia and Borrelia are commonly transmitted by Ixodes spp. associated with deer, and most of the knowledge on their biological cycles comes from northern latitudes of the globe. In this study, we performed genetic screenings to detect tick-borne agents in blood and Ixodes stilesi ticks collected from an insular population of threatened pudu (Pudu puda), a pygmy deer species that inhabits temperate rainforests of southern South America. Inferred by phylogenetic analyses for 18S rRNA, COI and cytb genes, our results unveiled a novel genospecies of Babesia (Babesia sp. pudui) genetically related to Babesia odocoilei, a species that infects Odocoileus virginianus deer in North America. Although blood of the deer was negative for Borrelia infection, multilocus sequencing typing performed in one I. stilesi tick revealed the occurrence of a novel genetic variant of Borrelia chilensis, differing 0.93% and 0.18% in flaB and pepX genes with the type of strain for the species, respectively. Such a genetic divergence could be the result of thousands of years of isolation because of recent glaciation events that separated pudus and their tick populations at Chiloé Island approximately 437,000years ago. The finding of a Babesia sp. has no precedents for wild and domestic ungulates in Chile and shows a novel piroplasmid that must be considered now on in rehabilitation centres and zoos that attend pudu deer. Further research is now necessary to confirm pathogenic roles.

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.

Ecological conservation, biodiversity, and agricultural education as integrated approaches for envisioning the future of sustainable agriculture in North America

ABSTRACT In the ancestral heritage, mankind was in mutual interdependence with nature through hunting of animals, gathering of plant materials, and animal and plant domestication. However, agricultural advancements detached mankind from that heritage. In U.S. Midwest, it is evident that technological advancements over the past half-century have threatened man’s survival for centuries to come. Midwest is the prime agricultural land area in North America whose recent glaciation, deep-rooted prairie vegetation, and moderate precipitation favored the development of highly productive rain-fed farming systems. Evolution of productionist paradigm in Midwest focused on crop and livestock output but also resulted in numerous environmental externalities. Massive use of agrochemicals, for instance, resulted in pollution of Mississippi River watershed from upperstream to downstream creating a hypoxic zone in Gulf of Mexico. In this situation, it is appropriate to ask what best alternatives for mankind could be to reconnect with nature to achieve sustainability. The study explored three approaches including ecological conservation, biodiversity, and agricultural education. Conservation requires the passion of humanity to reconnect with ecological community and become part of it rather than to manipulate and seek control over it to reap ecosystem services. In creating interrelated diversity through domestication of culturally important wild species, agroforestry, horticulture, animal husbandry, and inclusion of relevant language and culture we can define a sustainable world. Agricultural education brings in three tenets of teaching/learning, research/discovery, and extension/engagement to provide, create, and apply knowledge, respectively, which are important in communicating the science of sustainability.

Late glacial and Holocene glaciation history of North and Northeast Greenland

ABSTRACT Northeast Greenland is the place where the Greenland Ice Sheet (GrIS) experienced the largest areal changes since the Last Glacial Maximum. However, the age constraints of the last deglaciation are in some areas sparse. In this study, we use forty-seven new 10Be cosmogenic exposure ages to constrain the deglaciation of the present-day ice-free areas in Northeast Greenland. Our results show that the outer coast region was deglaciated between 12.8 ± 0.6 and 11.5 ± 0.2 ka and the region close to the present ice margin was deglaciated 2 to 4 ka later between 9.2 ± 0.3 to 8.6 ± 0.3 ka. By combining our new results with previously published data from North and Northeast Greenland, we find that the ice sheet advanced to the shelf edge between 26 and 20 cal. ka BP. The outer coast was deglaciated between 12.8 and 9.7 ka and the present ice extent was reached between 10.8 to 5.8 ka. The ice margin continued to retreat farther inland during the Middle Holocene before it readvanced toward its Little Ice Age position. The deglaciation was probably forced by a combination of increased atmospheric and ocean temperatures, but local topography also played an important role. These results add to the growing knowledge about the glaciation history of the GrIS and add useful constraints for future ice sheet models.

Enigmatic super-heavy pyrite formation: Novel mechanistic insights from the aftermath of the Sturtian Snowball Earth

It is not well understood how, in the immediate aftermath of the Sturtian Snowball Earth, marine sulfur cycling resulted in a global distribution of sedimentary pyrite with δ34S values higher than coeval seawater. Here, we analyze the quadruple sulfur isotope systematics of organic-bound sulfur (OS) from the lowermost post-Sturtian Datangpo Formation, South China, and identify two generations of OS formation, each sampling an isotopically distinct sulfate reservoir (δ34S ≈ 26‰ and 52–93‰) that differentially impacted its respective, co-occurring pyrite. Combining several lines of geochemical evidence, we argue that the first OS generation was the product of a sulfate-impoverished meltwater-influenced setting, with OS preservation being the result of resistance to acid hydrolysis. However, the second OS generation was sourced from H2S produced in sediments during early diagenesis via microbial reduction of a 34S-enriched sulfate pool derived from overlying euxinic or ferruginous seawater. This is the first ancient marine data set where all observed pyrite is more enriched in 34S than its associated OS. Our proposed origin may be applied to global superheavy pyrite (SHP) immediately after the Sturtian and is comparable to processes linked to freshwater-to-marine transitions during rising sea level in the wake of recent glaciation.

Lacustrine diatom oxygen isotopes as palaeo precipitation proxy - Holocene environmental and snowmelt variations recorded at Lake Bolshoye Shchuchye, Polar Urals, Russia

The diatom oxygen isotope composition (δ18Odiatom) from lacustrine sediments helps tracing the hydrological and climate dynamics in individual lake catchments, and is generally linked to changes in temperature and δ18Olake. Lake Bolshoye Shchuchye (67°53′N; 66°19′ E; 186 m a.s.l) is the largest and deepest freshwater reservoir in the Polar Urals, Arctic Russia. The diatom oxygen isotope interpretation is supported by modern (isotope) hydrology, local bioindicators such as chironomids, isotope mass-balance modelling and a digital elevation model of the catchment.The Bolshoye Shchuchye δ18Odiatom record generally follows a decrease in summer insolation and the northern hemisphere (NH) temperature history. However, it displays exceptional, short-term variations exceeding 5‰, especially in Mid and Late Holocene. This centennial-scale variability occurs roughly contemporaneously with and similar in frequency to Holocene NH glacier advances. However, larger Holocene glacier advances in the Lake Bolshoye Shchuchye catchment are unknown and have not left any significant imprint on the lake sediment record. As Lake Bolshoye Shchuchye is deep and voluminous, about 30–50% of its volume needs to be exchanged with isotopically different water within decades to account for these shifts in the δ18Odiatom record. A plausible source of water with light isotope composition inflow is snow, known to be transported in surplus by snow redistribution from the windward to the leeward side of the Polar Urals. Here, we propose snow melt variability and associated influx changes being the dominant mechanism responsible for the observed short-term changes in the δ18Odiatom record. This is the first time such drastic, centennial-scale hydrological changes in a catchment have been identified in Holocene lacustrine diatom oxygen isotopes, which, for Lake Bolshoye Shchuchye, are interpreted as proxy for palaeo precipitation and, on millennial timescales, for summer temperatures.

Time-varying uplift in Svalbard—an effect of glacial changes

SUMMARYWe analyse Global Navigation Satellite System (GNSS) data from Svalbard to understand how uplift rates are controlled by the elastic and viscoelastic response of the solid Earth to changes in glacier mass on annual, interannual, decadal, centennial and millennial timescales. To reveal local patterns of deformation, we filter the GNSS time-series with an enhanced common-mode filtering technique where the non-tidal loading signal is incorporated. This technique reduces the estimated uncertainties for 5-yr time-series from 0.8 to 0.3 mm yr–1. Analysis of the GNSS data with different software–GAMIT, GipsyX, and GINS–produce consistent results that all indicate large temporal variations in uplift. For example, at the Ny-Ålesund GNSS station, uplift varies between 6 and 12 mm yr–1 for different 5-yr periods, and also shows a significant increase in the last 15 yr. We show that this increase is due to climate change-related ice mass loss in Svalbard. We constrain recent glacier retreat on Svalbard using a series of digital elevation models, and then correct the GNSS-derived uplift records for the elastic signal from these ice mass changes. The residual uplift signal is relatively constant, confirming the hypothesis that current ice mass changes exert a strong influence on GNSS observations. The relatively constant record of residual uplift can be used to constrain other geophysical signals such as the viscoelastic response of the solid Earth to ice loading during the Little Ice Age and the Last Glacial Period. We review uplift results from previous viscoelastic modelling studies and show that the residual signal cannot yet be fully explained. Our new uplift results thus motivate the need for new viscoelastic modelling of the glacial isostatic adjustment process in Svalbard.