Last Glacial Research Articles

The genomic natural history of the aurochs.

Now extinct, the aurochs (Bos primigenius) was a keystone species in prehistoric Eurasian and North African ecosystems, and the progenitor of cattle (Bos taurus), domesticates that have provided people with food and labour for millennia1. Here we analysed 38 ancient genomes and found 4 distinct population ancestries in the aurochs-European, Southwest Asian, North Asian and South Asian-each of which has dynamic trajectories that have responded to changes in climate and human influence. Similarly to Homo heidelbergensis, aurochsen first entered Europe around 650 thousand years ago2, but early populations left only trace ancestry, with both North Asian and European B. primigenius genomes coalescing during the most recent glaciation. North Asian and European populations then appear separated until mixing after the climate amelioration of the early Holocene. European aurochsen endured the more severe bottleneck during the Last Glacial Maximum, retreating to southern refugia before recolonizing from Iberia. Domestication involved the capture of a small number of individuals from the Southwest Asian aurochs population, followed by early and pervasive male-mediated admixture involving each ancestral strain of aurochs after domestic stocks dispersed beyond their cradle of origin.

Comparing the evolution of debris-free and debris-covered glaciers during the end of the Lateglacial and the Holocene in Dudh Koshi basin, Everest region, Nepal

Debris-covered glaciers are very frequent geomorphological features in Khumbu Himal (Nepal). Rock debris on the glacier surface play a significant role in glacier-climate relationships and glacier dynamics. These effects may cause an asynchronous evolution of debris-covered glaciers compared to debris-free glaciers at a multicentennial to millennial scale. Here, we explore this hypothesis by documenting and comparing the multi-millennial Holocene evolution of a debris-free glacier, Sabai glacier, and two debris-covered glaciers, Dig and Huuku glaciers, from adjacent catchments in Dudh Koshi basin (Everest region, Nepal). To do so, we dated rock samples collected from moraine boulders on both debris-covered and debris-free glaciers using the 10Be cosmic ray exposure (CRE) dating method. 10Be CRE ages obtained from 41 moraine boulder samples provide time constraints from ∼13.5 ka to 0.1 ka. While at Dig (debris-covered) and Sabai (debris-free) glaciers, no moraines from the Lateglacial and the Early Holocene are preserved, debris-covered Huuku glacier evidenced a large glacier extent during the Bølling-Allerød and Early Holocene with two moraines dated respectively to ∼13.5 ka and 11 ka, synchronously with most debris-free and debris-covered glaciers in this region. These two glacier advances are concomitant with enhanced monsoon precipitation supporting a qualitative relationship. The absence of debris landforms in the main valley question the nature of Huuku glacier during the Bølling-Allerød and Early Holocene, which could have been either debris-free or covered by a thin debris layer only. During the Mid Holocene, significant differences are observed in the evolution of the two glacier types. The two debris-covered glaciers recorded a significant advance at ∼4.8 ka, synchronous with that observed on other debris-covered glaciers in Khumbu valley. However, such glacier advance during the Mid Holocene was not evidenced on debris-free glaciers in the Dudh Koshi valley. Such a Mid Holocene glacier advance may have a spatial signature with frequent cases reported from both types of glaciers in the western part of High Mountain Asia, which are however infrequent in the arid and semi-arid southern and north-eastern Tibet. During the Late Holocene, both types of glaciers evolved similarly again, with moraines spanning the last two millennia, including the Little Ice Age, concomitant with enhanced monsoon precipitation.

Palaeoenvironmental Analysis of the Southern Part of the Danube–Tisza Interfluve (Hungary): The Northern Loess Wall of Katymár and the Hay Meadows and Loess Banks of Hajós

This study presents a comprehensive palaeoenvironmental reconstruction of the southern part of the Danube–Tisza interfluve in the Carpathian Basin from the Late Pleistocene to Early Holocene, addressing the region’s response to global climate forcings and local environmental factors based on multi-proxy analyses of two key protected areas: the Katymár brickyard and the hay meadows and loess banks of Hajós. By integrating radiocarbon-dated malacological, macrobotanical, pollen, phytolith, geochemical, and sedimentological analyses, it was possible to provide a picture of past climate–ecosystem interactions. The Katymár North sequence provides an extended chronology from Marine Isotope Stage 3 (MIS3, ca. 36,000 cal BP) through the Last Glacial Maximum (LGM) and into the Early Holocene, while the Hajós sequence offers high-resolution data for the LGM–Holocene transition. By the late Ice Age, humidity and surface moisture superseded temperature as primary palaeoecological limiting factors, promoting ecotone-like forest–steppe environments during cooling periods.

Revealing the Demographic History of the European Nightjar (Caprimulgus europaeus).

A species' demographic history gives important context to contemporary population genetics and a possible insight into past responses to climate change; with an individual's genome providing a window into the evolutionary history of contemporary populations. Pairwise sequentially Markovian coalescent (PSMC) analysis uses information from a single genome to derive fluctuations in effective population size change over the last ~5 million years. Here, we apply PSMC analysis to two European nightjar (Caprimulgus europaeus) genomes, sampled in Northwest and Southern Europe, with the aim of revealing the demographic history of nightjar in Europe. We successfully reconstructed effective population size over the last 5 million years. Our analysis shows that in response to global climate change, the effective population size of nightjar broadly increased under stable warm periods and decreased during cooler spans and prolonged glacial periods. PSMC analysis on the pseudo-diploid combination of the two genomes revealed fluctuations in gene flow between ancestral populations over time, with gene flow ceasing by the last-glacial period. Our results are tentatively suggestive of divergence in the European nightjar population, with timings consistent with differentiation being driven by restriction to different refugia during periods of glaciation. Finally, our results suggest that migratory behaviour in nightjar likely evolved prior to the last-glacial period, with long-distance migration seemingly persisting throughout the Pleistocene. However, further genetic structure analysis of individuals from known breeding sites across the species' contemporary range is needed to understand the extent and origins of range-wide differentiation in nightjar.

Unravelling the fragmented sediment–landform assemblage in an area of thick Quaternary sediment, western Hudson Bay Lowland, Canada

Laurentide Ice Sheet (LIS) reconstructions beyond the last glacial maximum (LGM) are difficult due to the incomplete stratigraphic record that has largely been eroded by the most recent glaciation. The Hudson Bay Lowland (HBL) is a key region to understanding the long-term evolution of the LIS, as it contains an extensive pre-LGM stratigraphic record of glacial and nonglacial events. This study uses a hybrid lithostratigraphy–allostratigraphy approach to decipher the stratigraphic record in the Kaskattama highland region of the western HBL, Manitoba, Canada. We identify five glacial (till) units and three nonglacial (sorted sediment) units that were deposited during at least three glacial–interglacial cycles, which are constrained by radiocarbon and optical ages. The provenance of till units is identified using clast-lithology and detrital hornblende 40Ar/39Ar ages and supported by matrix geochemistry. The new stratigraphic framework for this western HBL region indicates that the ice emanating from the Quebec–Labrador dome advanced into the region at the start of the last two glaciations, suggesting accelerated early growth of the Quebec–Labrador dome relative to the Keewatin Dome.

10Be exposure age data reveals last Glacial ice sheet histories in the Larsemann Hills of East Antarctica

The Larsemann Hills, near Prydz Bay is a key region to understand past glacial dynamics in East Antarctica. However, multi-proxy dating constraints for this area indicate an uncertain chronology for deglaciation history during the Last Glacial Period (LGP). In this paper, we report two 26Al exposure ages and eighteen 10Be exposure ages for bedrock samples to better understand past glacial histories. The minimum exposure ages range from 82.84 ± 4.67 ka to 15.1 ± 1.95 ka corresponding to Marine Isotope Stages (MIS) 5–2, and fourteen of these samples date from before the global Last Glacial Maximum (c. 19–26.5 ka). Analysis of 10Be and 26Al concentrations implies a complex deglaciation history, and surface lowering of bedrock may have caused systematic age underestimation. Further comparisons between two adjacent regions in the Mirror Peninsula reveal two major deglaciations involving the Lambert Glacier / Amery Ice Shelf System dated at 61.78 ka to 82.84 ka and 21.98 ka to 49.5 ka during the LGP.

Late Glacial-Holocene cirque glacier chronology on sub-Antarctic Kerguelen Archipelago (49°S) based on cosmogenic 36Cl exposure dating

Sub-Antarctic glacier chronologies can provide valuable information about the past variability of climate dynamics in the Southern Ocean region. The Kerguelen Archipelago (49°S) is advantageously located under the influence of the Southern Hemisphere's westerly wind belt, thus fluctuations of climate-sensitive glaciers on Kerguelen can provide a baseline for understanding the behavior of this atmospheric regime in response to climatic forcings. We present 17 36Cl exposure ages of moraine and erratic boulders to provide chronological constraints to paleoglacier extents of the Guynemer cirque glacier, located just north of the Cook Ice Cap. Erratic boulders show ice thinning in the Guynemer region started to occur in the Late Glacial by 13.5 ka and continued past 12.4 ka. Ice retreat was punctuated by the formation of two moraine stages, the outermost at 11.5 ± 0.4 ka followed by another at 10.4 ± 1.2 ka, which are indicative of Early Holocene glacier advances/standstills. A glacial advance occurred at 1.4 ± 0.3 ka, which corroborates other Late Holocene re-advances elsewhere on the archipelago. Finally, three undated moraine stages are found between 1.4 ka and the 1960s. The lack of moraines after 10.4 ka and through the Mid-Holocene suggests that the Guynemer glacier was significantly smaller during this extended period of the Holocene compared to its Early Holocene as well as its Late Holocene limits. The Guynemer glacier history provides unique evidence of Early Holocene moraines on Kerguelen, which have not been discovered thus far on the archipelago. Similar to glaciers in Patagonia, New Zealand and South Georgia, the Guynemer glacier was at its largest Holocene extent in the Early Holocene. However, while other southern mid-latitude glacier chronologies show progressively smaller glacial extents throughout the Mid- to Late Holocene, the Late Holocene re-advance of the Guynemer glacier, like other Kerguelen glaciers, likely exceeded its Mid-Holocene extent.

Slow Dynamics of Hydrate Systems Revealed by a Double BSR

AbstractDetermining how gas hydrate distribution evolved along continental margins in the past is essential to understanding its evolution in the future. Moreover, hydrate decomposition has been linked to several catastrophic events, including some of the largest submarine landslides on Earth and the massive release of greenhouse gases into the ocean. Offshore Romania, the presence of a second bottom‐simulating reflector (BSR) provides an opportunity to gain valuable insights into hydrate dynamics since the Last Glacial Period (LGP). We conducted transient modeling of hydrate thermodynamic stability by merging in‐situ observations with indirect assessments of sea‐bottom temperature, thermal conductivity, salinity, sedimentation rate, and sea‐level variations. We reveal a strong correlation between the BSRs and the base of the Gas Hydrate Stability Zone (GHSZ) during both the present and LGP periods. The gradual evolution of the GHSZ over the past 34 ka presented here supports a conceptual model that excludes catastrophic environmental scenarios.

Hydrological response of Andean catchments to recent glacier mass loss

Abstract. The impacts of the accelerated glacier retreat in recent decades on glacier runoff changes are still unknown in most Andean catchments, increasing uncertainties in estimating water availability. This particularly affects the outer tropics and Dry Andes, heavily impacted by prolonged droughts. Current global estimates overlook climatic and morphometric disparities, which significantly influence model parameters, among Andean glaciers. Meanwhile, local studies have used different approaches to estimate glacier runoff in a few catchments. Improving 21st-century glacier runoff projections relies on calibrating and validating models using corrected historical climate inputs and calibrated parameters across diverse glaciological zones. Here, we simulate glacier evolution and related runoff changes between the periods 2000–2009 and 2010–2019 across 786 Andean catchments (11 282 km2 of glacierized area, 11° N to 55° S) using the Open Global Glacier Model (OGGM). TerraClimate atmospheric variables were corrected using in situ data, getting a mean temperature bias by up to 2.1 °C and enhanced monthly precipitation. Glacier mass balance and volume were calibrated, where melt factor and the Glen A parameter exhibited significant alignment with varying environmental conditions. Simulation outcomes were validated against in situ data in three documented catchments (with a glacierized area > 8 %) and monitored glaciers. Our results at the Andes scale reveal an average reduction of 8.3 % in glacier volume and a decrease of 2.2 % in surface area between the periods 2000–2009 and 2010–2019. Comparing these two periods, glacier and climate variations have led to a 12 % increase in mean annual glacier melt (86.5 m3 s−1) and a decrease in rainfall on glaciers of −2 % (−7.6 m3 s−1) across the Andes, with both variables comprising the glacier runoff. We confirmed the utility of our corrected regional simulations of glacier runoff contribution at the catchment scale, where our estimations align with previous studies (e.g., Maipo 34° S, Chile) as well as provide new insights on the seasonal glaciers' largest contribution (e.g., La Paz 16° S, Bolivia) and new estimates of glacier runoff contribution (e.g., Baker 47° S, Chile).

Two-Decadal Glacier Changes in the Astak, a Tributary Catchment of the Upper Indus River in Northern Pakistan

Snow and ice melting in the Upper Indus Basin (UIB) is crucial for regional water availability for mountainous communities. We analyzed glacier changes in the Astak catchment, UIB, from 2000 to 2020 using remote sensing techniques based on optical satellite images from Landsat and ASTER digital elevation models. We used a surface feature-tracking technique to estimate glacier velocity. To assess the impact of climate variations, we examined temperature and precipitation anomalies using ERA5 Land climate data. Over the past two decades, the Astak catchment experienced a slight decrease in glacier area (−1.8 km2) and the overall specific mass balance was −0.02 ± 0.1 m w.e. a−1. The most negative mass balance of −0.09 ± 0.06 m w.e. a−1 occurred at elevations between 2810 to 3220 m a.s.l., with a lesser rate of −0.015 ± 0.12 m w.e. a−1 above 5500 m a.s.l. This variation in glacier mass balance can be attributed to temperature and precipitation gradients, as well as debris cover. Recent glacier mass loss can be linked to seasonal temperature anomalies at higher elevations during winter and autumn. Given the reliance of mountain populations on glacier melt, seasonal temperature trends can disturb water security and the well-being of dependent communities.

Internal structures of the Last Glacial Period multi-phase cold-climate aeolian dunes, Moravian Sahara, Czechia

Aeolian dunes in the Moravian Sahara developed during at least two distinct cold periods of the Last Glacial Period (LGP). Their two-phase development was previously deduced solely from the different orientations of dune crests without knowing the dunes' internal structure. Therefore, the aim of this study is to reveal the internal structure of the multi-generational dune system, developing a relationship between these dunes and bedrock. A combination of data was obtained from boreholes and ground-penetrating radar (GPR). GPR profiles were measured in places where the boreholes were drilled to obtain high-resolution combination data. Furthermore, models of aeolian sediment thickness and ground-water level were made using a spatial interpolation method (ordinary kriging with spherical semivariogram model) to interpret the GPR results and to investigate a relationship between the bedrock and dune patterns. Internal dune structures confirm the hypothesis of multi-generational development of this dune field. In addition, the aeolian sediment distribution was influenced by the bedrock topography. Finally, this study confirms that a change in LGP atmospheric circulation over Central Europe is the main controlling factor in forming the Moravian Sahara dune field.

Prokaryotic richness and diversity increased during Holocene glacier retreat and onset of an Antarctic Lake

Knowledge about biodiversity changes during transitions from glacial landscape to lake formation is limited to contemporary studies. Here, we combined analyses of lithology, chronology and geochemistry with sedimentary ancient DNA metabarcoding to assess such transition in maritime Antarctica. We inferred three paleoenvironmental stages covering the Holocene glacier retreat process. From 4900 to 3850 years before the present, we found the lowest prokaryotic richness/diversity, with bacterial taxa indicators associated to soil and terrestrial environments. From 3850 to 2650 years before the present, a higher carbon content, higher Carbon/Nitrogen variability, increased species richness/diversity, and prokaryotic taxa indicators of long-term energy starvation were detected. Finally, from 2650 to 1070 years before the present, we inferred the onset of a genuine lacustrine environment holding stable Carbon/Nitrogen ratios and the highest prokaryotic diversity, with known aquatic bacterial taxa. Our study unveils for the first time the evolution from a glacier-covered to a freshwater lake through a millennial scale.

Late Pleistocene glaciations on the sub-Antarctic Kerguelen Archipelago: new evidence from 36Cl CRE dating and comparison with other southern mid-latitude glacier records

Previous paleo-glacial studies on Kerguelen showed a singular pattern of Holocene glacier evolution on this archipelago in comparison with other southern mid-latitude glacier records. In this study, we aim to test this singularity on a longer timescale, based on 26 new in situ-produced 36Cl ages from pre-Holocene glacio-geomorphic features. Samples from moraine boulders and glacially polished bedrock were extracted at six different sites, located near the Port-aux-Français scientific station (PAF site), on Longue Island, Australia Island, on the Port-Jeanne d’Arc Peninsula (PJDA site), on the Gallieni Peninsula at Baie Larose (BLR site) and the McMurdo Island. The moraine ages indicate that glacier culminations occurred during Marine Isotopic Stage 3 (MIS 3) at 42.2 ± 4.9 ka on the PAF site, and during the global Last Glacial Maximum (gLGM) at 21.5 ± 3.2 ka on the PJDA site and at 21.4 ± 3.7 ka and 19.4 ± 2.6 on Baie Larose site. This is the first time that Late Pleistocene glacier culminations are evidenced on Kerguelen by direct moraine dating, thus allowing comparison with other moraine records from the southern mid-latitudes. While it remains speculative whether or not the MIS 3 glacial maximum at ⁓42.2 ka is in phase with other glaciers at this latitude (due to high age uncertainties), the gLGM glacial maximum is synchronous with that in other southern mid-latitude regions. 36Cl CRE ages of glacially polished bedrock surfaces sampled in different locations of the archipelago vary from ⁓39 ka to ⁓19 ka. We interpret these results as reflecting periods of deglaciation that occurred in between the two glacier culminations and right after the gLGM on Kerguelen. These ages also suggest that some places of the archipelago were free of ice at least since ⁓ 39 ka. The presence of a MIS 3 moraine at PAF site that has not been obliterated by a gLGM advance suggests that the ⁓ 42.2 ka glacier extent was at least as large as gLGM glacial maxima on the archipelago. The glacier culmination during MIS 3 being larger than that during the gLGM on the Kerguelen Archipelago matches observations in other southern mid-latitude regions. Late Pleistocene glacier culminations on Kerguelen may have been in phase with cold temperatures recorded in SST records, which suggest a cooling around Kerguelen. However, climate drivers responsible for the larger MIS 3 glacier culmination on Kerguelen still remain unclear even if we hypothesize that changes in precipitation may have superimposed on temperature changes.

Late Wolstonian and Ipswichian (MIS 6/5e) sediment fill in a limestone sinkhole, Askham Fell, northern England

ABSTRACTIn 2019 a sinkhole (doline) occurred in Late Devensian till above fissured limestone in northern England. Most sediment plugging the fissure was evacuated down into a karstic drainage system. The residual sedimentary fill comprises three main lithofacies, dated using optically stimulated luminescence to between 170.7 ± 40.0 and 56.1 ± 13.5 ka. The earliest date demonstrates fissures were present in the limestone pavement at the time of Marine Isotope Stage (MIS) 6, or shortly thereafter. The fissure filled with fine sand and silt due to surface runoff and aeolian processes probably at the MIS 6 to MIS 5e transition after Wolstonian glacial ice had retreated. The deposits then collapsed into the karst system. Further fine sand and silt deposition occurred during MIS 3; this deposit filled the central cavity surrounded by residual MIS 6/5e deposits. The sequence was capped by till as Late Devensian (MIS 2) ice transgressed the area. Solution fissures in the karst surfaces of northern England may pre‐date the Late Devensian glaciation. Moreover, fissures are repositories of pre‐Devensian sediment deposits which survived the Late Devensian glaciation and the Ipswichian interglacial. Such sites should provide information on the nature and timing of pre‐Devensian glacial–interglacial events and shed light on basal ice conditions and glaciokarst drainage behaviour.

The infill of tunnel valleys in the central North Sea: Implications for sedimentary processes, geohazards, and ice-sheet dynamics

Tunnel valleys are widespread in formerly glaciated regions such as the North Sea and record sediment transport beneath ice sheets undergoing deglaciation. However, their complex infill architecture often makes their implications for ice-sheet processes difficult to unravel. Here, we use high resolution 3D (HR3D) seismic-reflection data, improved-resolution conventional 3D seismic-reflection data, and geotechnical information from industry-acquired boreholes to image the infill architecture of buried Quaternary tunnel valleys in the North Sea in unprecedented detail. Ten cross-cutting generations of tunnel valleys are mapped beneath the seafloor of the North Sea where only seven were visible previously. Each generation of tunnel valleys potentially reflects a different glaciation, although our evidence may imply that it is possible to rapidly erode and infill multiple generations of tunnel valleys within a single glacial cycle. The infill of the oldest tunnel valley generations reflects sedimentation during relatively gradual ice-sheet retreat, with occasional episodes of overriding by re-advancing grounded ice. Tunnel valleys formed in more recent glaciations are characterised by more variable sedimentation patterns that reflect dynamic fluctuations of the ice margin, including readvances and stagnation, during valley filling and ice retreat. Numerous subglacial landforms are also imaged within the tunnel valleys; these sometimes contain shallow gas accumulations that represent geohazards for seafloor infrastructure installations. In addition, we document examples where salt diapirism has caused fluids to migrate upwards from depth through faults and into the near-surface tunnel valleys. In instances where this occurs, the relatively porous and often highly continuous subglacial landforms present within their infill may allow these fluids to spread laterally for kilometres or even escape from the seafloor; it is therefore important to consider tunnel valleys when monitoring possible CO2 leakage from carbon capture and storage efforts.

Speleothem growth and stable carbon isotopes as proxies of the presence and thermodynamical state of glaciers compared to modelled glacier evolution in the Alps

In recent years, glacier modelling proved to be an essential tool for simulating Quaternary glacier evolution in the European Alps. Yet, only sparse empirical data mostly concentrated around the Last Glacial Maximum (LGM) is available to validate these simulations. On the other hand, speleothems from the Alps are a widespread palaeoclimate archive. They provide stable carbon isotope records, which can inform about soil and vegetation conditions above a cave site but also potentially about the lack of soil during times of glacier coverage. In addition, speleothem growth in cold, high-elevation cave sites during glacials are a strong indicator of temperatures in the soil-karst-cave system above the freezing point, which is only likely to occur if the cave is covered by a warm-based glacier.Here we use existing speleothem data (growth histories and stable carbon isotopes) from Alpine caves to infer soil coverage (i.e. glacier absence) and thermodynamical states of the glaciers during the last glacial cycle and to statistically assess the compatibility with modelled glacier reconstructions. We compare data from multiple cave sites located at different elevations (870–2512 m a.s.l.) with recent glacier evolution simulations. We find a general agreement between speleothem-derived soil presence or absence and modelled glacier coverage. However, speleothem data provide evidence of surface temperatures above freezing point if covered by a glacier, which is not fully reproduced by the simulations. Our work demonstrates the unique value of speleothem-based reconstructions as proxies to assess the performance of palaeo-ice flow models in a transient manner, whereas only maximum glacier state was considered before due to lack of data.

Late Holocene glacier and climate fluctuations in the Mackenzie and Selwyn mountain ranges, northwestern Canada

Abstract. Over the last century, northwestern Canada experienced some of the highest rates of tropospheric warming globally, which caused glaciers in the region to rapidly retreat. Our study seeks to extend the record of glacier fluctuations and assess climate drivers prior to the instrumental record in the Mackenzie and Selwyn mountains of northwestern Canada. We collected 27 10Be surface exposure ages across nine cirque and valley glacier moraines to constrain the timing of their emplacement. Cirque and valley glaciers in this region reached their greatest Holocene extents in the latter half of the Little Ice Age (1600–1850 CE). Four erratic boulders, 10–250 m distal from late Holocene moraines, yielded 10Be exposure ages of 10.9–11.6 ka, demonstrating that by ca. 11 ka, alpine glaciers were no more extensive than during the last several hundred years. Estimated temperature change obtained through reconstruction of equilibrium line altitudes shows that since ca. 1850 CE, mean annual temperatures have risen 0.2–2.3 ∘C. We use our glacier chronology and the Open Global Glacier Model (OGGM) to estimate that from 1000 CE, glaciers in this region reached a maximum total volume of 34–38 km3 between 1765 and 1855 CE and had lost nearly half their ice volume by 2019 CE. OGGM was unable to produce modeled glacier lengths that match the timing or magnitude of the maximum glacier extent indicated by the 10Be chronology. However, when applied to the entire Mackenzie and Selwyn mountain region, past millennium OGGM simulations using the Max Planck Institute Earth System Model (MPI-ESM) and the Community Climate System Model 4 (CCSM4) yield late Holocene glacier volume change temporally consistent with our moraine and remote sensing record, while the Meteorological Research Institute Earth System Model 2 (MRI-ESM2) and the Model for Interdisciplinary Research on Climate (MIROC) fail to produce modeled glacier change consistent with our glacier chronology. Finally, OGGM forced by future climate projections under varying greenhouse gas emission scenarios predicts 85 % to over 97 % glacier volume loss by the end of the 21st century. The loss of glaciers from this region will have profound impacts on local ecosystems and communities that rely on meltwater from glacierized catchments.

Late Quaternary glacier advances in the Andes of Santiago, central Chile, and paleoclimatic implications

Andean mountain glaciers in central Chile are in a transitional zone between the seasonal influence of the mid-latitude westerlies and subtropical semiarid conditions to the north. Long-term glacial dynamics for these glaciers and their relationship with the paleoclimate during the late Quaternary are poorly known despite their relevancy. We estimate here the timing and extent of late Pleistocene–early Holocene glaciers in the Andes of Santiago (33°50′S) from geomorphological and geochronological analyses. Our observations evidence that a glacial stage occurred before the Last Glacial Maximum (ELGM) at the San Gabriel drift (1,300 m a.s.l.), dated as ∼46–36 ka. Glacial stages during the latest Pleistocene–early Holocene transition period, partially concomitant with the Antarctic Cold Reversal (ACR) and with the Younger Dryas (YD) chronozones, were identified and dated at La Engorda drift (2,450–2,570 m a.s.l.) at ∼15–10 ka. We propose that the San Gabriel drift represents a prolonged glacial advance driven by increased precipitation and cold conditions off central Chile during glacial times. In La Engorda drift, late glacial advances occurred associated with increased regional precipitation, in the context of a transition from humid to arid climate in central Chile, concomitantly with a general warming trend of sea surface temperatures offshore in the southeastern Pacific and with reduced austral summer insolation. The results support the sensitivity of the Andean mountain glaciers to precipitation and paleoclimate conditions, most possibly associated with periods of increased northward influence of the mid-latitude westerlies during glacial and late glacial times, in addition to the El Niño/Southern Oscillation (ENSO) impact since the mid-Holocene, driving late Quaternary glacier advances in central Chile. We estimate a maximum variation of ∼1,200 m in the position of the late Quaternary Equilibrium Line Altitude (ELA), inferred at ∼3,400 m and ∼3,600 m a.s.l. at the time of the San Gabriel and La Engorda drifts, respectively, with respect to its modern location close to 4,600 m a.s.l.

Evidence of site-specific and male-biased germline mutation rate in a wild songbird.

Germline mutations are the ultimate source of genetic variation and the raw material for organismal evolution. Despite their significance, the frequency and genomic locations of mutations, as well as potential sex bias, are yet to be widely investigated in most species. To address these gaps, we conducted whole-genome sequencing of 12 great reed warblers (Acrocephalus arundinaceus) in a pedigree spanning three generations to identify single-nucleotide de novo mutations (DNMs) and estimate the germline mutation rate. We detected 82 DNMs within the pedigree, primarily enriched at CpG sites, but otherwise randomly located along the chromosomes. Furthermore, we observed a pronounced sex bias in DNM occurrence, with male warblers exhibiting three times more mutations than females. After correction for false negatives and adjusting for callable sites, we obtained a mutation rate of 7.16 × 10-9 mutations per site per generation (m/s/g) for the autosomes, and 5.10 × 10-9 m/s/g for the Z chromosome. To demonstrate the utility of species-specific mutation rates, we applied our autosomal mutation rate in models reconstructing the demographic history of the great reed warbler. We uncovered signs of drastic population size reductions pre-dating the last glacial period (LGP), and reduced gene flow between western and eastern populations during the LGP. In conclusion, our results provide one of the few direct estimates of the mutation rate in wild songbirds and evidence for male-driven mutations in accordance with theoretical expectations.

Pre-industrial Holocene glacier variability in the tropical Andes as context for anthropogenically driven ice retreat

Pre-industrial Holocene glacier variability in the tropical Andes as context for anthropogenically driven ice retreat