Search

The Neotropical region exhibits the greatest worldwide diversity and the diversification history of several clades is related to the puzzling geomorphologic and climatic history of this region. The freshwater Amazon ecoregion contains the main hydrographic basins of the Neotropical region that are highly dendritic and ecologically diverse. It contains a rich and endemic fish fauna, including one of its most iconic and economically important representatives, the bony‐tongue Arapaima gigas (Teleostei, Osteoglossiformes). Here, we evaluated the projected distribution of the genus in different historical periods (Present, Last Glacial Maximum, Last Interglacial Maximum and Near Future) and interpreted these results in light of the genomic diversity and modeled historical demography. For that, we combined species distribution models, population genetic analysis using SNPs and deep learning model selection. We analyzed a representative sample of the genus from the two basins where it naturally occurs, four localities in the Amazon (Am) and three in the Tocantins‐Araguaia (To‐Ar) basin, as well as individuals from three fish farms. We inferred a potentially smaller distribution in the glacial period, with a possible refuge in central Am. Our genetic data agrees with this result, suggesting a higher level of genetic diversity in the Am basin, compared to that observed in To‐Ar. Our deep learning model comparison indicated that the To‐Ar basin was colonized by the population from the Am basin. Considering a global warming scenario in the near future, A. gigas could reach an even larger range, especially if anthropogenic related dispersal occurs, potentially invading new areas and impacting their communities. [ABSTRACT FROM AUTHOR]

Environmental conditions were particularly severe during the Last Glacial Maximum, altering the distribution of the Southern Hemisphere biota, particularly at higher latitudes. The copepod Boeckella poppei is the only macroscopic continental invertebrate species known to be distributed today across the three main biogeographic regions in Antarctica as well as in southern South America. Signy Island (South Orkney Islands) is a unique location for the study of Antarctic freshwater ecosystems due to its location and geographic isolation; it contains 17 lakes in several low altitude catchments. We conducted phylogeographic and demographic analyses using the cox1 gene on 84 individuals of B. poppei from seven lakes across Signy Island. We recorded low levels of genetic diversity and a strong genetic differentiation signal between the eastern and western valleys within the island. Phylogeographic structure and demographic inference analyses suggested at least one asymmetrical dispersal event from west to east. Demographic inference detected a strong signal of population growth during the deglaciation process, which may have followed either (1) a strong genetic bottleneck due to a reduction in population size during the last glacial period, or (2) a founder effect associated with postglacial recolonization of Signy Island from elsewhere. The genetic architecture of this island's populations of B. poppei shows that historical events, rather than continuous dispersal events, likely played a major role in the species' current distribution. Finally, our study considers possible mechanisms for dispersal and colonization success of the most dominant species in the Antarctic freshwater community. [ABSTRACT FROM AUTHOR]

Many formerly glaciated valleys in the western United States preserve detailed glacial features that span the penultimate glaciation through the last deglaciation; however, numerical age control is limited in many of these systems. We report 35 new cosmogenic 10Be surface exposure ages of moraine boulders in the Sawatch Range, Colorado. Eight ages suggest Bull Lake moraines in Lake Creek (range: 132–120 ka, n = 4) and Clear Creek (range: 187–133 ka, n = 4) valleys may correlate with Marine Isotope Stage 6. In Lake Creek valley, 22 10Be ages from Pinedale end moraines average 20.6 ± 0.6 ka, and 5 10Be ages from a recessional moraine average 15.6 ± 0.7 ka, indicating that glaciers occupied two extended positions at ~21–20 and ~16 ka. The glacial extent dated to ~16 ka was nearly as great as that of the earlier glacial phase, suggesting that climate conditions in the Colorado Rocky Mountains at this time were similar to those of the last glacial maximum. Combining these moraine ages with seven previously published 10Be ages from cirque and valley-bottom bedrock reveals that the Lake Creek paleoglacier lost 82% of its full glacial length in ~1.5 ka and was completely deglaciated by ~14 ka. [ABSTRACT FROM AUTHOR]

Historical factors (colonization scenarios, demographic oscillations) and contemporary processes (population connectivity, current population size) largely contribute to shaping species' present‐day genetic diversity and structure. In this study, we use a combination of mitochondrial and nuclear DNA markers to understand the role of Quaternary climatic oscillations and present‐day gene flow dynamics in determining the genetic diversity and structure of the newt Calotriton asper (Al. Dugès, 1852), endemic to the Pyrenees. Mitochondrial DNA did not show a clear phylogeographic pattern and presented low levels of variation. In contrast, microsatellites revealed five major genetic lineages with admixture patterns at their boundaries. Approximate Bayesian computation analyses and linear models indicated that the five lineages likely underwent separate evolutionary histories and can be tracked back to distinct glacial refugia. Lineage differentiation started around the Last Glacial Maximum at three focal areas (western, central and eastern Pyrenees) and extended through the end of the Last Glacial Period in the central Pyrenees, where it led to the formation of two more lineages. Our data revealed no evidence of recent dispersal between lineages, whereas borders likely represent zones of secondary contact following expansion from multiple refugia. Finally, we did not find genetic evidence of sex‐biased dispersal. This work highlights the importance of integrating past evolutionary processes and present‐day gene flow and dispersal dynamics, together with multilocus approaches, to gain insights into what shaped the current genetic attributes of amphibians living in montane habitats. [ABSTRACT FROM AUTHOR]

Indian Summer Monsoon is an integrated component of the global climate system. The spatial movement of ITCZ at seasonal and orbital time scales is revealed in the ensemble of terrestrial and marine records covering the last 50,000 years. We deduced an evolutionary shift in the precipitation pattern near western North America, Mediterranean and South East Asia based on the estimated water isotopic composition from the d 18Ocarbonate and ?47 thermometer. Record revealed three stages of climate transition: mid-Holocene optima, Younger Dryas and Last glacial maxima. Estimated mean arid air temperature was 9-21°C and 12-35°C during the glacial and interglacial periods respectively. The June summer solar insolation at 65°N is captured in the temperature record linked with ice volume, atmospheric CO2 levels and sea surface temperature; factors influencing the monsoon precipitation near mid-latitude region worldwide. [ABSTRACT FROM AUTHOR]

In this study, we reconstruct the style and timing of the late Quaternary glaciation pattern in a monsoon‐influenced precipitation‐deprived region (Dhauliganga basin) of the Central Himalaya, India. We present detailed glacial geomorphological mapping, aided by topographic maps, satellite data, Google Earth images, global positioning system survey, and fieldwork (2012–2016), supported by optically stimulated luminescence dating. Our estimations for the magnitude and timing of glaciations in the Dunagiri Glacier valley suggest five sets of glacial advances. The oldest and most extensive glaciation in the valley is Dunagiri Glacial Stage‐I (DGS‐I), represented by a discontinuous moraine ridge that extends down to ~10 km (~2700 m a.s.l.) from the present‐day snout and formed 13±3 ka ago. The DGS‐II stage is present in the form of a small moraine ridge (~3000 m a.s.l.) that formed 10±2 ka ago. The DGS‐III and DGS‐IV glacial stages are 5–6 and 3.5 ka old, respectively. The youngest stage (DGS‐V) of glaciation belongs to the neoglaciation episode ~1 ka old. DGS‐I was the most extensive glaciation in the valley and corresponds to the global Last Glacial Maximum. The glaciation pattern has been correlated with the estimates of present‐day glacial sensitivity to climate and explained as forcing factors sensitive to temperature and/or precipitation changes in the region. The observations based on prevailing climate correlations aided with available glacial chronologies underline the fact that the glacial valleys situated in precipitation‐deprived regions have preserved a delayed response of glaciations to regional climatic events. Estimates for the timing and pattern of glaciation in the present study suggest a different glacial response for modern climatic changes than previously assumed. We suggest that individual glacier advances differ regionally and, therefore, should not be interpreted as being reflective of a change in the whole geographical belt. [ABSTRACT FROM AUTHOR]

Recent studies have revealed that areas with sediments pre‐dating the Last Glacial Maximum (LGM) are more common in southern Scandinavia than previously thought. To describe and date these sediments, four deep core drillings were made in the province of Småland and one in the province of Skåne, southern Sweden, and an existing sediment core drilled in the 1960s from western Skåne was re‐examined and re‐dated, all resulting in a total of 22 AMS14C ages and 19 OSL ages. With some variation, the Småland cores revealed a stratigraphy consisting of two tills separated by proglacial glacifluvial and glaciolacustrine sediments, dated to between 46 and 25 ka BP (Marine Isotope Stage (MIS) 3 to 2). We suggest that these sediments were deposited during the ice retreat of the MIS 3 Ristinge glaciation over southern Scandinavia and the subsequent re‐advance over the area in late MIS 3 and the transition to MIS 2, when the Fennoscandian Ice Sheet advanced towards its LGM marginal positions. This contrasts with the glacial history of southernmost Skåne, where an additional MIS 3 glacial advance through the Baltic Basin is suggested by the occurrence of a till named Allarp Till. New ages from sediments under‐ and overlying this till suggest that it was deposited by the same Baltic ice advance that deposited the Klintholm till in Denmark (Klintholm Advance, c. 34–30 ka). [ABSTRACT FROM AUTHOR]

The Kellerjoch forms a small isolated massif at the northernmost rim of the central Eastern Alps of Tyrol and shows a number of geomorphological features of glacial and periglacial origin. Mapping yields evidence of two local glaciations postdating the Last Glacial Maximum. Using a simple glaciological approach the palaeoglaciers related to these events were reconstructed. The older glaciation yields an equilibrium line altitude (ELA) ranging from 1660 m for the maximum extent to 1800 m a.s.l. for the innermost moraine. For the younger glaciation, ELAs were reconstructed at 1905 m and 1980 m (depending on the reconstruction) for the Kellerjoch palaeoglacier 2, as well as 1870 m and 2060 m a.s.l. for the Proxen palaeoglacier and the Gart palaeoglacier, respectively. A comparison with published data from the Eastern Alps shows that the older glaciation in the Kellerjoch region likely corresponds to the Gschnitz stadial. Low basal shear stresses of the glacier tongues point towards a cold and dry climate, similar to the reconstruction for the Gschnitz type locality at Trins. The younger glaciation cannot unambiguously be assigned to a specific Late Glacial ice advance, but a Younger Dryas age is a distinct possibility. [ABSTRACT FROM AUTHOR]

Anthropogenic climate change raises growing concerns about its potential catastrophic impacts on both ecosystems and human societies. Yet, several studies on damages induced on the economy by unmitigated global warming have proposed a much less worrying picture of the future, with only a few points decrease in the world GDP per capita by the end of the century, even for a global warming above +4 °C. We consider here two different empirically estimated functions linking GDP growth or GDP level to temperature at the country level and apply them to a global cooling of -4 °C in 2100, corresponding to a return to glacial conditions. We show that the alleged impact on global average GDP per capita runs from -1.8 %, if temperature impacts GDP level, to +36 %, if the impacts are rather on GDP growth. These results are then compared to the hypothetical environmental conditions faced by humanity, taking the last glacial maximum as a reference. The modeled impacts on the world's GDP appear strongly underestimated given the magnitude of climate and ecological changes recorded for that period. After discussing the weaknesses of the aggregated statistical approach to estimate economic damages, we conclude that, if these functions cannot reasonably be trusted for such a large cooling, they should also not be considered as providing relevant information on potential damages in the case of a warming of similar magnitude, as projected in the case of unabated greenhouse gas emissions. [ABSTRACT FROM AUTHOR]

The population dynamics of freshwater organisms are expected to be related to the connectivity among comparable streams, ponds, or rivers in a patchy habitat. Here, we investigated the population dynamics of the giant water bug, Belostoma angustum Lauck 1964 (Hemiptera: Belostomatidae), in a fine-scale spatial sampling, and evaluated which gene flow model previously described for freshwater organisms could explain the genetic–morphological variation in this species. For these purposes, we evaluated genetic and morphological variations, as well as the demographic history of this freshwater insect. Our genetic analyses showed a lack of geographical structure within B. angustum populations across the evaluated range, concordant with widespread gene flow model. Our findings of the demographic history of B. angustum suggest recent and rapid expansion beginning during the late Pleistocene after the Last Glacial Maximum. Likewise, we did not find geographically structured morphological variation in B. angustum , except for body size. The lack of structure of genetic–morphological variation in B. angustum could be explained by a stepping ponds system resulting in the widespread gene flow detected among populations of this species. The warmer and wetter climatic conditions after the last glacial period may have favored the demographic expansion of B. angustum populations due to the increasing of potential freshwater habitats and food resources. This favorable habitat probably allowed the stepping ponds dispersal mode resulting in the verified geographically unstructured genetic–morphological variation. [ABSTRACT FROM AUTHOR]

The Late Pleistocene landscape in northern Eurasia and North America was inhabited by a specific megafaunal complex, which largely disappeared during the Pleistocene/Holocene transition. Vegetation changes are considered as one of the factors responsible for these extinctions, but the structure and composition of the Pleistocene vegetation are still poorly known. Here we complement previous studies by comparing the taxonomic composition of the plant remains found in the gastrointestinal tracts of the frozen carcasses of Pleistocene megaherbivores with the species composition of the current Siberian vegetation. We compiled a dataset of palaeobotanical records from frozen individuals of Pleistocene megaherbivores found in northern Siberia and Beringia and dated to the period from more than 50 kyr BP to 9 kyr BP. We also compiled a dataset of vegetation plots from several regions in Siberia. We analysed the similarity in taxonomic composition of plants between these two datasets using a novel method that accounts for variable taxonomic resolution in palaeobotanical data. For most megaherbivore individuals, plant remains in their gastrointestinal tracts corresponded to tundra, forest and mire vegetation, while they showed low similarity to steppe. This pattern was relatively constant over time, showing no remarkable differences between the Last Glacial Maximum and the periods before and afterwards. This suggests that during the Upper Pleistocene, a mosaic of mesic and wet vegetation types such as tundra with patches of forests and mires was common in northern Siberia and Beringia. In contrast, the steppe was rare to absent in the landscape or underused by the megaherbivores as a pasture since they found enough food in the widespread mesic and wet habitats with more productive vegetation. [ABSTRACT FROM AUTHOR]

The dust cycle plays an important role in the long-term evolution of the climate and environment. In this paper, an improved climate model including aerosol processes was used to carry out a set of sensitivity experiments and comparative analyses of the effects of high-latitude ice-sheet extent and abnormal dust erosion, as well as Earth's orbital parameters and atmospheric greenhouse gas content, on dust activities during the last glacial maximum. The comparative analysis found that incorporating the abnormal surface erosion factor alone could increase dust emissions by 2.77-fold and 3.77-fold of the present-day global and Asian dust emissions, respectively. The high-latitude ice-sheet factor caused global dust emissions to increase by 1.25-fold that of the present day. Sensitivity experiments showed that increased surface erosion in Asia during the last glacial maximum made the greatest contribution to the increased dust emissions in Asia, followed by the high-latitude ice-sheet factor, while the contributions of the greenhouse gas content and orbital parameters were relatively weak. Strong dust emissions during the glacial period were therefore not only dependent on the development of the high-latitude ice sheets but were strongly associated with the underlying surface characteristics of local dust source regions. [ABSTRACT FROM AUTHOR]

The centre–periphery hypothesis (CPH) predicts a decrease in population performance from the centre of the species range towards the edge, hindering further species expansion. To overcome ecological limitation, local adaptation of peripheral populations is assumed necessary to extend niche space and thus to potentially facilitate species' range expansion. However, adaptive changes do not necessarily correspond to current ecological marginality. Instead, population history may provide a fuller context for understanding patterns of local adaptation within the species range. Here, we test whether local adaptation to current conditions or population history explains current species distribution in an annual heterocarpic species, Atriplex tatarica, in Europe.The ecological marginality and population history were estimated using species distribution modelling (SDM) and 13 microsatellite markers. To test for local adaptation, we evaluated the population performance of both currently central and currently peripheral populations by examining germination and growth traits under three temperature regimes which correspond to climates within, and also beyond, the species' current range. In addition, we compared Bayesian estimates of population differentiation in neutral genetic markers and quantitative traits to test whether phenotypic differentiation evolved by local adaptation and/or by genetic drift.We found population performances were not higher under temperature regimes corresponding to their own habitats, suggesting that patterns of local adaptation in A. tatarica do not correspond to current centre–periphery gradient. In contrast, population genetic structure indicated a strong influence of population history on variation in the germination phenotype. Species distribution modelling suggested a more cold–tolerant genetic lineage was located at the species' range periphery during the last glacial period. Meanwhile, a less cold‐tolerant lineage occurred in more suitable areas of the Balkan Peninsula and southeastern Europe. Thus, the strong adaptive divergence of evolutionary lineages may reflect a past spatial arrangement of populations, suggesting evolutionary processes were ongoing before the species expansion in the Holocene.Synthesis. There is continuing debate concerning the speed of trait evolution during range expansion with climate change. Our results indicate that even though organisms expand rapidly, the accompanying evolutionary changes in phenotype may be much slower. [ABSTRACT FROM AUTHOR]

The well‐known vicariance and dispersal models dominate in understanding the allopatric pattern for related species and presume the simultaneous occurrence of speciation and biogeographic events. However, the formation of allopatry could postdate the species divergence. We examined this hypothesis using DNA sequence data from three chloroplast fragments and five nuclear loci of Dipelta floribunda Maxim. and D. yunnanensis Franch, two shrub species with the circum Sichuan Basin distribution, combining the climatic niche modeling approach. The best‐fit model supported by the approximate Bayesian computation analysis indicated that D. floribunda and D. yunnanensis diverged during the mid‐Pleistocene period, consistent with the largest glacial period in the Qinghai–Tibet Plateau. The historically interspecific gene flow was identified, but seemed to have ceased after the last interglacial period, when the range of D. floribunda moved northward from the south of the Sichuan Basin. Furthermore, populations of D. floribunda had expanded obviously in the north of the Sichuan Basin after the last glacial maximum (LGM). Relatively, the range of D. yunnanensis expanded before the LGM, and reduced during the post‐LGM especially in the north of the Sichuan Basin, reflecting the asynchronous responses of related species to contemporary climate changes. Our results suggested that complex topography should be considered in understanding distributional patterns, even for closely related species and their demographic responses. [ABSTRACT FROM AUTHOR]

Orographic and climatic oscillations have played crucial roles in shaping the nucleotide diversity and evolutionary history of many species across the Northern Hemisphere. In this study, based on 10 nuclear loci and a chloroplast DNA marker, we analyzed the nucleotide polymorphisms and demographic history of the endangered conifer species Pinus bungeana in Northwest China and investigated the phylogenetic relationships between P. bungeana and two related species, that is, Pinus gerardiana and Pinus squamata. We found that P. bungeana exhibited an extremely low level of nucleotide diversity (πsil = 0.00159). Demographic simulations based on DIYABC analysis showed that P. bungeana underwent demographic expansion and contraction during the Miocene. According to ecological niche modeling, we found that this species survived in situ during the glacial period and was not restricted to southern refugia. We speculate that P. bungeana may have experienced widespread population shrinkage from the Last Interglacial to the Last Glacial Maximum due to geological or climatic events. Isolation‐with‐migration analysis revealed that the divergence (~2.4–4.2 Ma) among P. bungeana and its related species was significantly associated with the Qinghai–Tibetan Plateau uplift events in the mid‐to‐late Tertiary period. Species tree analyses suggested that these three related Pinus species formed a monophyletic clade with high bootstrap support. These results suggest that the Miocene–Pliocene and Pleistocene geological and climatic fluctuations might have profoundly affected the nucleotide diversity and demography of this psychrotolerant conifer species in western China. [ABSTRACT FROM AUTHOR]

Knowledge of the glacial chronologies for the Last Glacial Maximum (LGM) helps in understanding the interactions between climate, topography and glacier development. In this sense, the investigation of the Lake Orta moraine amphitheatre (Alpine foreland, northern Italy) allowed spatial and temporal reconstruction of the Orta Glacier. The end‐moraine system was investigated by means of geomorphological field surveys, analysis of 13 rock samples for cosmogenic 10Be and 36Cl concentrations, and remote sensing analysis. The dating results indicate that the age of the outer moraine belt is concordant with the LGM culmination at 26.5–23 ka, as found in other amphitheatres in the Alps. This new age estimate of the outermost moraines shows that the maximum extent of the Orta Glacier during the LGM was significantly bigger than recently suggested. A younger stabilization phase of the glacier front at about 19 ka indicates that the onset of the withdrawal of glaciers from the lower Alpine valleys started later. Provenance analysis of the boulders shows that the greatest contribution of ice to the Orta Glacier came from the Anzasca Valley rather than the major Ossola Valley. This reflects the closeness (about 45 km) to the foreland of the high‐elevated accumulation area of the Monte Rosa massif (4634 m a.s.l.), whose eastern glacier seems to have reached the lower valley faster than the trunk Toce Glacier. This fact underlines the key role played by high‐elevation accumulation areas that are located close to the foreland in controlling the path and geometry of major glaciers in the Alps. [ABSTRACT FROM AUTHOR]

Subaqueous calcite deposited on the walls of Devils Hole 2 cave (Nevada, USA) represents a unique archive for geochemical variations within the regional aquifer. Here, we present a 475,000-year record of initial 234U/238U activity ratios in delta notation (δ234U0). Results show a range in values from 1851-1616‰. Variations in δ234U0 coincide with interglacial-glacial cycles over the past 475,000 years. Maximum δ234U0 values correspond to the last five glacial intervals, during which southwest Nevada experienced cool, pluvial conditions. Minimum δ234U0 values correspond to interglacial intervals, during which this region experienced warm, arid conditions. We propose that an elevated water table during glacial periods inundated previously dry bedrock and basin sediments, thereby leaching excess 234U accumulated in these materials. We interpret Devils Hole 2 cave δ234U0 as a proxy for water-rock interactions in this regional aquifer, which is ultimately tied to the surface moisture conditions at recharge zones. The mechanism proposed here serves as a testable hypothesis and possible analogue for future subaqueous speleothem studies in similar hydrogeologic settings. Due to its unprecedented duration, the Devils Hole 2 cave δ234U0 record provides the first paleo-moisture record in southwest Nevada for marine isotope stages 10-12. In addition, high-precision δ234U measurements of modern groundwaters sampled from Devils Hole 2 cave are presented. [ABSTRACT FROM AUTHOR]

The article presents a research on transient glacial incision in the Patagonian Andes. Topics discussed include decline in local erosion rates may represent a return toward a balance between rock uplift and erosion; spatially averaged erosion rate of the source region after the onset of alpine glaciations; and negative feedback operates at the scale of the glacial network.

Larix laricina (eastern larch, tamarack) is a transcontinental North American conifer with a prominent disjunction in the Yukon isolating the Alaskan distribution from the rest of its range. We investigate whether in situ persistence during the last glacial maximum (LGM) or long‐distance postglacial migration from south of the ice sheets resulted in the modern‐day Alaskan distribution. We analyzed variation in three chloroplast DNA regions of 840 trees from a total of 69 populations (24 new sampling sites situated on both sides of the Yukon range disjunction pooled with 45 populations from a published source) and conducted ensemble species distribution modeling (SDM) throughout Canada and United States to hindcast the potential range of L. laricina during the LGM. We uncovered the genetic signature of a long‐term isolation of larch populations in Alaska, identifying three endemic chlorotypes and low levels of genetic diversity. Range‐wide analysis across North America revealed the presence of a distinct Alaskan lineage. Postglacial gene flow across the Yukon divide was unidirectional, from Alaska toward previously glaciated Canadian regions, and with no evidence of immigration into Alaska. Hindcast SDM indicates one of the broadest areas of past climate suitability for L. laricina existed in central Alaska, suggesting possible in situ persistence of larch in Alaska during the LGM. Our results provide the first unambiguous evidence for the long‐term isolation of L. laricina in Alaska that extends beyond the last glacial period and into the present interglacial period. The lack of gene flow into Alaska along with the overall probability of larch occurrence in Alaska being currently lower than during the LGM suggests that modern‐day Alaskan larch populations are isolated climate relicts of broader glacial distributions, and so are particularly vulnerable to current warming trends. [ABSTRACT FROM AUTHOR]

We present triple oxygen isotope data from speleothems obtained by an O 2 -CO 2 Pt-catalyzed oxygen-isotope equilibration method. The high precision (9 per meg or better, 1σ SD) of our new speleothem Δ17O (carbonate 17O anomaly) data is sufficient to resolve subtle hydroclimatic signals. In addition, we determined triple oxygen isotope fractionation factors through two sets of modern paired carbonate-dripwater samples collected at a temperature of 17 ± 1 °C, which is prerequisite to the calculation of triple oxygen isotope compositions of parent meteoric waters at cave sites from speleothem proxy data. Based on this calibration, we back calculated triple oxygen isotope compositions of parent waters across well-characterized climate transitions using speleothem proxy data from three regions. Resulting oxygen-isotope data closely track the Global Meteoric Water Line (GMWL), providing a preliminary validation of the method. Our speleothem Δ17O data indicate a 21 per meg difference between Marine Isotope Stage 5d and 5e in samples from Central Asia and a 15 per meg difference between the Middle Holocene and Last Glacial Maximum in samples from the eastern Mediterranean, suggesting a shift in moisture source and/or fractionation history. Unexpectedly, there were no measurable Δ17O differences between glacial and interglacial samples from both the South American (western Amazon) and Asian (southern China) monsoon domains, implying consistent moisture-source conditions across glacial and interglacial cycles, at least in terms of relative humidity. Remarkably, Δ17O values from the western Amazonian samples are significantly higher (∼20 per meg) than those from Asian monsoon regions, suggesting lower relative humidity along moisture trajectories in the western Amazon during the Middle Holocene and Last Glacial period. Similarly, Δ17O values of the eastern Mediterranean samples are significantly higher (19–55 per meg) than those of coeval samples from Central Asia, implying different hydrological environments or moisture sources despite being in the same westerly circulation domain. Speleothem Δ17O data may thus provide new and important constraints for understanding regional and global hydroclimate dynamics. [ABSTRACT FROM AUTHOR]

Human activity and adaptation to the environmental changes on the Tibetan Plateau are under debate. Lhasa is the political and economic center of the Tibetan Autonomous Region, and was founded in 633 CE. The potential influence of changes in the landscape on the city's development remains unclear. In this paper, Optically Stimulated Luminescence and radiocarbon dating techniques were applied to determine the ages at which landforms have changed since the Last Glacial Maximum. Combining these results with archaeological evidence taken from two sites (Changguogougou and Qugong) in the Lhasa area, the relation between river migration and human adaptation was discussed. The results showed that the two terraces (T1 and T2) of the Lhasa River are fill and fill-cut terraces, and that these were formed after ∼0.9 ka and 0.08 cal ka BP, respectively. Land below an altitude of ∼3650 m above sea level was lower than the contemporary river level before ∼0.9 ka. With the down cutting of the Lhasa River, more land became available for habitation, and the flood risk to the city was also reduced. The sediments of the lower land around the Changguogou area appear to be composed of fluvial-lacustrine facies during the Last Glacial Maximum to Late Glacial period. Humans occupied the Changguogou area between ∼3.255 and 3.059 cal ka BP when there was weaker aeolian activity, and abandoned this area when aeolian activity increased. [ABSTRACT FROM AUTHOR]

Late Paleolithic sites are widely distributed in southern China (the Lingnan region), which mainly includes Guangxi, Guangdong and Hainan Provinces. Although the pebble-tool industry continued to exist in this region during the late Pleistocene and Holocene, small flake-tool assemblages also occurred and coexisted with the pebble-tool industry. The pebble-tool industries are characterized by unifacially flaked cobble tools, dominated by choppers, with a lack of handaxes. The small flake-tool industries have the characteristics of retouched tools made on small fakes, a striking platform of cores without preparation, and scrapers dominating the assemblage. Technologically and typologically, pebble tools of the Late Palaeolithic are similar to those of the Early Palaeolithic in this region and belong to the same cultural tradition. The small flake tools, however, are very different from the cobble tools. They belong to a new technology which originated from northern China during the period of the Last Glacial Maximum. The Late Palaeolithic industries in Southern China are similar to those in mainland Southeast Asia. [ABSTRACT FROM AUTHOR]

Background: The biodiversity and distributions of terrestrial snails at local and regional scales are influenced by their low vagility and microhabitat specificity. The accessibility of large-bodied species and their characteristically high levels of genetic polymorphism make them excellent ecological and evolutionary models for studies on the phylogeography, phylogenetics, and conservation of organisms in fragmented populations. This study aims to elucidate the biodiversity, systematics, and distributions of genetic lineages within the genus Oreohelix at the northern and western periphery of their range. Results: We found four mitochondrial clades, three of which are putative subspecies of Oreohelix subrudis. One clade was geographically widespread, occurring within numerous sites in Cypress Hills and in the Rocky Mountains, a second was geographically restricted to the Rocky Mountains in Alberta, and a third was restricted to the Cypress Hills region. A fourth clade was the small-bodied species, O. cooperi. ITS2 sequence and screening data revealed three genetic clusters, of which one was O. cooperi. Cluster 1 contained most individuals in COI clade X and some from clade B and cluster 2 was predominantly made up of individuals from COI clades B and B′ and a few from clade X. ITS2 alleles were shared in a narrow contact zone between two COI clades, suggestive of hybridization between the two. Conclusions: A sky island known as Cypress Hills, in southeastern Alberta, Canada, is a biodiversity hotspot for terrestrial land snails in the genus Oreohelix. The observed phylogeographic patterns likely reflect reproductive isolation during the Last Glacial Maximum, followed by secondary contact due to passive, long-range dispersal resulting from low vagility, local adaptation, and complex glacial history. [ABSTRACT FROM AUTHOR]

We describe a revised understanding of the extent and dynamics of the Cordilleran Ice Sheet on the continental shelf of northern British Columbia, Canada. During the Local Last Glacial Maximum ice streams occupied two of the three shelf‐crossing troughs in Queen Charlotte Sound (Goose Island Trough and Mitchell's Trough). A 25‐km‐wide outlet glacier – the Hecate Glacier – flowed south in Hecate Strait, parallel to the modern coast. It reached a grounding line at the head of Moresby Trough, beyond which an ice shelf may have extended to the edge of the continental shelf. The southern part of the independent Haida Gwaii Ice Cap formed piedmont lobes dissected by tunnel valleys. An emergent area east of the ice cap corresponds with the 'Hecate Refugium'. Subsequently grounded ice retreated from the shelf troughs. The Hecate Glacier retreated incrementally towards the north. Mainland ice then stabilized along a north–south margin near modern coasts, marked by submarine moraines in the Chatham Sound region and elsewhere. This margin was probably contemporaneous with the margin at Mt. Buxton in the south of the study area, dated at c. 17.6 to 16.6 ka (Darvill et al. 2018). Off southern Haida Gwaii, following a retreat, a glacier margin was established east of Juan Perez Sound, proximal to a sandur plain graded to a water level of −150 m. Proglacial lakes to the north of here were the likely source for an outburst flood that created gravel bedforms and an area of sea‐floor scour. In a final phase, diminished mainland ice was confined to fjords, and retreat inland is marked by a succession of submarine moraines, beginning with a moraine at the Douglas Channel sill. The study has implications for the migration pathways for human migration into the Americas following deglaciation. [ABSTRACT FROM AUTHOR]

Several regions of tropical America show imprints of past glacial activity. These relict landforms can support the understanding of past climate conditions, such as during the Last Glacial Maximum (LGM), and the implications that these paleoclimatic conditions could have had on landscape change. Here, we present and analyze glacial morphologies for the Chirripó National Park in Costa Rica based on aerial imagery (1:25,000), detailed Digital Elevation Models, geomorphic mapping, as well as geomorphic assessments in the field to determine and validate landforms. This study adds valuable insights into the reconstruction of the maximum expansion of tropical glaciation during the LGM in Costa Rica and into tropical America glacial landscapes in general. [ABSTRACT FROM AUTHOR]

In order to improve our understanding of the glacier development in central Switzerland after the Last Glacial Maximum (LGM) we have investigated moraine sequences of the Meiental ("Meien Valley"). After detailed field mapping, 34 rock samples from different moraines were dated with cosmogenic 10Be. The results indicate that the investigated moraines were deposited mainly in multi-phased advances during the Younger Dryas and the Early Holocene. This is in contrast to previous studies suggesting that these moraines represent equivalents of the older Egesen and Daun stadials. Our research shows that factors like the moraine stratigraphy, geomorphologic characteristics and Equilibrium Line Altitudes (ELAs) can be highly variable at the local scale, which can make a stadial correlation of moraines between different valleys or regions within the Alps without age constraints ambiguous and unreliable. Against this background the interregional use of ΔELAs (ELA depressions) for the parallelization of paleo-glacier extents throughout the Alps according to traditional concepts of the Lateglacial deglaciation seems to produce erroneous results. Geomorphologic observations and relative age relationships in moraine successions based on stratigraphic principles remain the basis for this kind of studies. But as we show here it is necessary to base local-to-regional comparisons of moraine sequences additionally on exposure ages to gain a more reliable understanding of the Alpine glacier development after the LGM. • Marked multi-phased glacier advances in the Meiental were dated to the YD. • Several Early Holocene moraines prove separate glacier advance/standstill phases. • YD age of supposed Daun equivalent moraines challenges the Daun as separate stadial. • The ELA approach for the parallelization of paleo glaciers has evident limitations. • Exposure dating results increasingly update the concept of the Alpine deglaciation. [ABSTRACT FROM AUTHOR]

The Upper Gunnison Basin (UGB), Colorado, is a montane region characterized by unusual physiography and topographic isolation. Excavations of three caves in the UGB provide one of the most diverse records of high-elevation late Quaternary vertebrates in North America. The localities, Haystack Cave (2450 m above sea level [m asl]), Cement Creek Cave (2860 m asl), and Signature Cave (3055 m asl), together provide a near-continuous record of vertebrate communities that extends from before the last glacial maximum to the present. These communities largely represent a sagebrush (Artemisia spp.) steppe-tundra environment that prevailed throughout the UGB in the late Pleistocene. At least five taxa of extinct large mammals disappear from the UGB by the Early Holocene; one small mammal (the short-faced skunk Brachyprotoma cf. B. brevimala) also became extinct. The fossil record further indicates that only four small extant mammals (Sorex preblei , Dicrostonyx sp., Lemmiscus curtatus , and Urocitellus elegans) were extirpated from the UGB by the Early Holocene, in part because of community restructuring and loss of open habitats with expansion of forests to higher elevations. An analysis of taxonomic richness and evenness at Cement Creek Cave indicates high resilience in the small mammal community despite major climate shifts over the past 40,000+ yr. [ABSTRACT FROM AUTHOR]

We present a 3D reconstruction of ice thickness distribution across the New Zealand Southern Alps at the Last Glacial Maximum (LGM, c. 30–18 ka). To achieve this, we used a perfect plasticity model which could easily be applied to other regions, hereafter termed REVOLTA (Reconstruction of Volume and Topography Automation). REVOLTA is driven by a Digital Elevation Model (DEM), which was modified to best represent LGM bed topography. Specifically, we removed contemporary ice, integrated offshore bathymetry and removed contemporary lakes. A review of valley in-fill sediments, uplift and denudation was also undertaken. Down-valley ice extents were constrained to an updated geo-database of LGM ice limits, whilst the model was tuned to best-fit known vertical limits from geomorphological and geochronological dating studies. We estimate a total LGM ice volume of 6,800 km3, characterised predominantly by valley style glaciation but with an ice cap across Fiordland. With a contemporary ice volume of approximately 50 km3, this represents a loss of 99.25% since the LGM. Using the newly created ice surface, equilibrium line altitudes (ELAs) for each glacier were reconstructed, revealing an average ELA depression of approximately 950 m from present. Analysis of the spatial variation of glacier-specific ELAs and their depression relative to today shows that whilst an east-west ELA gradient existed during the LGM it was less pronounced than at present. The reduced ELA gradient is attributed to an overall weakening of westerlies, a conclusion consistent with those derived from the latest independent climate models. • LGM ice thickness reconstructed on a glacier-by-glacier basis, driven by geomorphological evidence. • During the LGM the Southern Alps contained approximately 6,800 km3 ice, of which 99.25% has subsequently been lost. • Predominant valley style of glaciation across Southern Alps but ice cap across Fiordland. • The spatial patterns of precipitation and temperature were markedly different at the LGM compared to present. • Westerly circulation was apparently reduced during the LGM when compared to present. [ABSTRACT FROM AUTHOR]

Ravine slopes at the recently discovered Nesseltalgraben site in southeastern Germany provide a unique last glacial sediment record for the Northern Calcareous Alps. The 21 m-long profile is dominated by fine-grained lacustrine-palustrine sediments overlain by several metres of glacifluvial gravels and lodgement tills of the Last Glacial Maximum and underlain by a diamicton. The age model includes 29 radiocarbon analyses and one paleomagnetic anomaly (Laschamp event) together providing a modelled age range from c. 59 to 29.6 ka cal BP, i.e. Marine Isotope Stage (MIS) 3. Apart from a description of the lithofacies, X-ray-fluorescence (XRF) scanning and elemental analyses provide high-resolution records of the geochemical composition. Multivariate analyses of XRF data separate Ca from other major elements. Carbonate contents, represented by Ca and total inorganic carbon, reach maxima in repeatedly occurring calcareous silty to sandy layers and are related to glacigenic origin. These indicate repeated inner-alpine local glaciations during stadials of MIS 3. Scanning electron microscopy and XRF data confirm the detrital origin of these layers. In contrast, organic matter and elements more resistant to chemical weathering (Si, Ti, Zr) accumulated during interstadials and concurrently elevated Rb/Sr ratios indicate intensified weathering. The high-frequency proxy variations determined for Nesseltalgraben reflect interstadial-stadial climate variability comparable with oxygen-isotope records from Greenland ice cores and Alpine speleothems. Thus, Nesseltalgraben is among the very few independently dated sediment records from continental Europe covering the entire MIS 3 and reflecting the full Dansgaard-Oeschger climate variability. • A 21 m sediment profile from the northern Alps provides a complete MIS 3 record. • Repeated inner Alpine glaciations are indicated by calcareous silt layers. • Geochemical variations reflect interstadials and stadials. • The Nesseltalgraben record resembles Dansgaard-Oeschger climate variability. [ABSTRACT FROM AUTHOR]

We present a new speleothem trace element and stable isotope record, which extends previous paleoclimate evidence from Cuban speleothems to the last 96 ka. Multiple proxies were used to reconstruct Caribbean hydroclimate and vegetation variability on orbital to millennial timescales. In particular, our proxies indicate a shift to more abundant C4 (compared to C3) vegetation and/or reduced soil activity during the Last Glacial Maximum. Our results confirm the close link of Caribbean precipitation to North Atlantic climate variability, insolation and the strength of the AMOC. The associated variations in Caribbean SSTs and the shift of the ITCZ led to substantial changes of the hydrological cycle. During Heinrich stadials and the Younger Dryas, climatic conditions in Cuba were comparably cool and/or dry. In contrast, warm Dansgaard/Oeschger cycles were accompanied by wetter conditions. This is in agreement with previous evidence for the Holocene, when North Atlantic cold (warm) events coincided with drier (wetter) conditions in Cuba. However, our record suggests that this connection was less dominant during MIS 4 to 2, displayed by a weaker response of Cuban precipitation proxies to North Atlantic climatic shifts than during MIS 5 in particular. We hypothesize that this weaker connection was mediated by the southerly shifted ITCZ inducing a weakening of the northern branch of the Caribbean Low Level Jet which usually transports moisture into the northern Caribbean. This is strikingly evident during the deglacial, when driest and/or coolest conditions prevailed in western Cuba from HS 1 into the Younger Dryas, including the Bølling/Allerød interstadial. • We present a speleothem stable isotope and trace element record from Cuba covering the last 96 ka. • Our results confirm the link of Caribbean precipitation to NA climate, insolation and the AMOC. • Last glacial NA cold (warm) (inter)-stadials coincided with drier (wetter) conditions in Cuba. • During the LGM, our data suggest a shift to more abundant C4 vegetation and/or reduced soil activity. • We observe a weaker sensitivity to NA climate during MIS 4 to 2 compared to MIS5 and the Holocene. [ABSTRACT FROM AUTHOR]

For many species, climate oscillations drove cycles of population contraction during cool glacial periods followed by expansion during interglacials. Some groups, however, show evidence of uniform and synchronous expansion, while others display differences in the timing and extent of demographic change. We compared demographic histories inferred from genetic data across marine turtle species to identify responses to postglacial warming shared across taxa and to examine drivers of past demographic change at the global scale. Using coalescent simulations and approximate Bayesian computation (ABC), we estimated demographic parameters, including the likelihood of past population expansion, from a mitochondrial data set encompassing 23 previously identified lineages from all seven marine turtle species. For lineages with a high posterior probability of expansion, we conducted a hierarchical ABC analysis to estimate the proportion of lineages expanding synchronously and the timing of synchronous expansion. We used Bayesian model averaging to identify variables associated with expansion and genetic diversity. Approximately 60% of extant marine turtle lineages showed evidence of expansion, with the rest mainly exhibiting patterns of genetic diversity most consistent with population stability. For lineages showing expansion, there was a strong signal of synchronous expansion after the Last Glacial Maximum. Expansion and genetic diversity were best explained by ocean basin and the degree of endemism for a given lineage. Geographic differences in sensitivity to climate change have implications for prioritizing conservation actions in marine turtles as well as for identifying areas of past demographic stability and potential resilience to future climate change for broadly distributed taxa. [ABSTRACT FROM AUTHOR]

• Abrupt mid-Holocene ice sheet thinning in western Weddell Sea Embayment since LGM. • At least 385 m of ice surface lowering revealed by in situ 14C exposure dating. • Glacial-geological data neither preclude nor mandate late Holocene ice thinning. The glacial history of the westernmost Weddell Sea sector of Antarctica since the Last Glacial Maximum is virtually unknown, and yet it has been identified as critical for improving reliability of glacio-isostatic adjustment models that are required to correct satellite-derived estimates of ice sheet mass balance. Better knowledge of the glacial history of this region is also important for validating ice sheet models that are used to predict future contribution of the Antarctic ice sheet to sea level rise. Here we present a new Holocene deglacial chronology from a site on the Lassiter Coast of the Antarctic Peninsula, which is situated in the western Weddell Sea sector. Samples from 12 erratic cobbles and 18 bedrock surfaces from a series of presently-exposed ridges were analysed for cosmogenic 10Be exposure dating, and a smaller suite of 7 bedrock samples for in situ 14C dating. The resulting 10Be ages are predominantly in the range 80–690 ka, whereas bedrock yielded much younger in situ 14C ages, in the range 6.0–7.5 ka for samples collected from 138–385 m above the modern ice surface. From these we infer that the ice sheet experienced a period of abrupt thinning over a short time interval (no more than 2700 years) in the mid-Holocene, resulting in lowering of its surface by at least 250 m. Any late Holocene change in ice sheet thickness — such as re-advance, postulated by several modelling studies — must lie below the present ice sheet surface. The substantial difference in exposure ages derived from 10Be and 14C dating for the same samples additionally implies ubiquitous 10Be inheritance acquired during ice-free periods prior to the last deglaciation, an interpretation that is consistent with our glacial-geomorphological field observations for former cold-based ice cover. The results of this study provide evidence for an episode of abrupt ice sheet surface lowering in the mid-Holocene, similar in rate, timing and magnitude to at least two other locations in Antarctica. [ABSTRACT FROM AUTHOR]

Mongolian glaciers have been the subject of relatively little research, resulting in less geochronological constraint than other parts of Central Asia. The Khentey Mountains (latitude 47–51°N, longitude 105–112°E) are a typical landlocked mountain range exhibiting clear geomorphic evidence of late Quaternary glaciation. Yet, compared to western parts of Mongolia such as the Mongolian Altay, Gobi Altay, Khangay, and Khovsgol, glacial history of the Khentey Mountains is poorly understood. To address this, and permit comparison of the Khentey glacier–climate record with other alpine regions in Mongolia, we performed geomorphological mapping and cosmogenic 10Be surface‐exposure dating in two glaciated regions of the Khentey Mountains: Yestii and Khagiin Khar. Specifically, we measured 10Be in 34 samples collected from five moraine sequences, which, together with morphostratigraphy, correspond to four main glacial stages: (i) The My1 terminal moraine sequence for Yestii (21.0±4.9 ka) and the Mk1 moraine for Khagiin Khar (19.6±2.6 ka), both of which represent the Last Glacial Maximum; (ii) the Lateglacial Mk2 moraine, dated to 16.0±3.5 ka; (iii) the Mk3 moraine, dated to either 17.6±7.0 ka (Lateglacial) or 12.1±1.1 ka (Younger Dryas); and (iv) the currently undated Mk4 moraine (~2200 m a.s.l.), to which we assign a Holocene age. Our results suggest that the timing of maximum glacier extent in Mongolia was regionally variable. In the Khentey Mountains, maximum glaciation occurred during Marine Isotope Stage (MIS) 2, whereas the maximum occurred during MIS 3 in Khangay and Khovsgol and during MIS 4 in the Altay. The MIS 2 glacial maximum in the Khentey Mountains coincided with the global sea level minimum during the Last Glacial Maximum, and was followed by at least three glacial re‐advances during the Lateglacial to possibly the Holocene. [ABSTRACT FROM AUTHOR]

Geomorphic evidence of past glaciation, such as U-shaped valleys, aretes, glacial lakes, and moraines, is preserved in the highland surrounding Cerro Chirripó in the Cordillera de Talamanca, Costa Rica. Previous work to establish a glacial chronology has focused on relative age dating of moraines and on radiocarbon dating of basal lake sediments to infer the timing of deglaciation. We used cosmogenic 36Cl surface exposure dating to constrain the ages of moraines within two formerly glaciated valleys, the Morrenas and Talari valleys. Forty-nine boulder samples were processed and measured from four moraine complexes in the Morrenas Valley and two moraine complexes in the Talari Valley. The exposure ages of these samples indicate a major glacial event occurred in this area from ~25 to 23 ka, broadly synchronous with the global last glacial maximum. Our results also indicate periods of glacial retreats and standstills from the deglacial period to the Early Holocene (~16–10 ka) before the complete disappearance of glaciers in this highland. These findings provide important insights into the glacial chronology and paleoclimate of tropical America. [ABSTRACT FROM AUTHOR]

The study reports new progress in the glacial chronology from new sites of the Chandrabhaga basin (CBB), in the transition between the monsoon-dominated Pir-Panjal and westerly-dominated Trans-Himalaya, in Lahaul Himalaya, northern India. In Lahaul, asynchronous timing and extent of late Quaternary glaciation have been reported (Owen et al., 1996, 1997, 2001; Eugster et al., 2016). The established oldest and most extensive Chandra glacial stage remains undated and speculative. We reconstrain the late Quaternary glacial history of the CBB using geomorphological mapping, Infrared Stimulated Luminescence (IRSL), and 10Be Terrestrial Cosmogenic Nuclide (TCN) surface exposure dating. Four local glacial stages, the Chandra, Batal, Dali, and Mayar, are dated to 152 ± 37 ka, 43 ± 3 ka, 21 ± 2 ka, and Little Ice Age (LIA), respectively. The Chandra glacial stage occurred during the late Marine Isotope Stage 6 (MIS 6; 186 ± 30 – 113 ± 10 ka), characterized by thick lacustrine deposits, rock benches, and ice-polished bedrock along the trunk valleys. We provide the first account of the numerical dating for the much-awaited timing of the oldest Chandra glacial stage (Local Last Glacial Maxima; LLGM) in the region (Owen et al., 1996, 1997, 2001). The LLGM in the CBB is much older than previously thought, consistent with some north-western Himalayan glacial records (Owen et al., 2006; Dortch et al., 2010; Hedrick et al., 2011; Orr et al., 2017, 2018). The Batal glacial stage occurred during the early MIS 3 and is dated from 47 ± 2 to 41 ± 2 ka, represented by the thick lacustrine deposits and impressive terraces. The Dali glacial stage occurred during the Global Last Glacial Maxima (GLGM), characterized by a series of latero-frontal recession moraines, and is dated from 33 ± 3 to 13 ± 1 ka. The Mayar advance is characterized by confined tributary valley glaciers, possibly linked to the LIA advance reported in the Kulti valley (Saha et al., 2019). The equilibrium line altitude (ELA) rose by ∼486 m from the Chandra stage to the present day. Reconstructed ELA altitudes and glacial landforms demonstrate that the glaciations have become progressively less extensive through time for this part of the western Himalayas. The change of glaciation pattern during the later part of the Quaternary is like other regions of the orogen, such pattern of glaciation reflects regional climatic forcing. • We present novel 31 k-feldspar IRSL and 12 cosmogenic 10Be ages for Chandrabhaga basin, Lahaul Himalaya, northern India. • We reconstrain the chronology for previously reported glacial stages of Owen et al. (1996, 1997, 2001). • We assign previously reported but undated oldest Chandra glacial stage with absolute dates. • Findings suggest that the Local Last Glacial Maxima in this region is much older than previously thought. • Our findings strongly question the previously reported extensions and rate of deglaciation during Global Last Glacial Maxima by Eugster et al. (2016). • Each glaciation becomes progressively less extensive in the Chandrabhaga basin over time. [ABSTRACT FROM AUTHOR]

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 CO 2 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. • Fluvial‒aeolian interplay during the Last Glacial Period in the Mu Us Desert is studied. • Regional fluvial‒aeolian interplay during the LGP is correlated with millennial-scale global climatic events. • Regional fluvial‒lacustrine environments mainly developed during the LGMa and B/A. • Fluvial‒lacustrine and aeolian processes during the LGP were mutually reinforcing and inhibiting. • High-northern-latitude forcing modulated millennial-scale fluvial‒aeolian interplay in the MUD during the LGP. [ABSTRACT FROM AUTHOR]

The area of Prados del Cervunal (PC) is an intra-morainic topographic depression located at 1800 m asl in the divide or interfluve between Garganta de Gredos and Garganta del Pinar valleys (Central Gredos; Iberian Central System, ICS). Both valleys, along with the adjacent Hoya Nevada, were occupied by glaciers during the Upper Pleistocene, leading to the development of the Prados del Cervunal moraine complexes studied in this work. Using cartographic methods and morphostratigraphic analysis, the three main glacial formations established in the Regional Chrono-Evolutionary Pattern for the ICS, Peripheral Deposits (PD), Principal Moraine (PM) and Internal Deposits (ID), have been identified and mapped in this area. The chronology of these formations has been implemented by Cosmic Ray Exposure (CRE) techniques using 10Be (new data) and 26Cl (previous data, recalculated in this work) in samples from morainic boulders. With these data, the following chrono-evolutionary sequence has been established: (stage 1) local-Maximum Ice Extent (MIE), dated in 25.0 ± 1.4 ka and corresponding to the maximum age obtained in these paleoglaciers; (stage 2) period of oscillations around the MIE, corresponding to the development of the PD Formation between ∼25 ka and ∼21 ka; (stage 3) period of readvance and stabilisation, dated after ∼21 ka (average age obtained for the PD moraines attached to PM moraines) and previous to ∼18 (minimum age obtained for a main crest of the PM formation); and (stage 4) onset of deglaciation dated around to ∼18 ka (average of ages obtained for the first main crest of the ID formation). During the stages of maximum ice expansion, these three glaciers formed an Ice field whose tongues were interconnected on the PC flat by an ice transfluence system (stages 1 and 2, Plateau Glacier Period). In later stages, the ice masses were partitioned, giving rise to valley glaciers and large moraines forming morainic complexes like those of PC (stages 2, 3 and 4, Valley Glaciers Period). The local MIE and onset of deglaciation stages in this area show a good fit with the ages stablished to global level for the global Last Glacial Maximum (LGM) and the onset of the Last Glacial Termination (Termination I). They also show good correlation at local (with other areas of the ICS), peninsular (with other Iberian mountains) and continental (some areas of the Alps and mountains of Central Europe) level. Finally, this evolutionary sequence and its correlations allowed us to adjust and validate some units of the Regional Chrono-Evolutionary Pattern model and propose the Gredos-Pinar-Cabeza Nevada glacial system as benchmark for the glaciation of Marine Isotope Stage (MIS) 2 in the Iberian Peninsula. • Glacial Geomorphology in Gredos-Pinar-Hoya Nevada paleoglaciers (Sierra de Gredos) was established. • New chronological data using 10Be-CRE and previous data using 36Cl-CRE allowed precise dating of the morainic sequence. • Previous morphology, aspect and debris supply as control factors of the layout and the dynamic of ice masses. • The local-MIE and the subsequent major readvance stage occurred within the LGM period. • Model from MIS 2 glaciation in the area: MIE and readvance post-MIE intra-LGM immediately followed by Termination-I. [ABSTRACT FROM AUTHOR]

Greenland ice cores are high-resolution archives of atmospheric dust loadings, and a better understanding of the driving factors for Greenland ice core dust concentration changes across different timescales provides valuable insights into changes in atmospheric circulation pattern and benefits accurate prediction of the future climate change. In this review, based on the dust record of the NGRIP ice core, we systematically characterize the Greenland ice core dust concentration changes during the last glacial, and explore the major driving factors. Based on the previously published and own isotope geochemical data, we estimated the contributions of provenances for the Greenland ice core dust using Bayesian mixing model. The results show that the dust was primarily derived from the Asian sources (71.5%) during the last glacial. Dust contributions from Africa include 10% (6.9%–12.0%) for North Africa and 6.9% (4.6%–8.3%) for West Africa; East Central & Eastern European sources yield a contribution of 11.7% (8.3%–14.2%). Subsequently, we revealed the contributions of multiple influencing factors to the Greenland dust concentration changes on orbital scale using the Lindeman-Merenda-Gold method. Overall, the findings highlight a crucial role of the aeolian activity in sources in controlling the Greenland dust concentration changes during the last glacial (63%). A greater effect of the atmospheric transport efficiency is also found (31%); by contrast, the spilt of jet stream contributes little (5%). We further argue that the variations in contribution of intensity of aeolian activity in sources through MIS2–5 may be related to glacial grinding by high-altitude mountain glaciers, with more significant influences exerted on variations in the Greenland dust concentration during MIS4. Inversely, more contributions from the atmospheric transport efficiency (wet deposition and reduced atmospheric residence time of dust in transport pathways) are observed during MIS2 and MIS3 than those during MIS4 and maybe MIS5. The results about changes of the contributing proportions of intensities of the different source areas throughout MIS2–5 reveal the major roles of sediment availability, displacement of the westerlies, and the Northern Hemisphere ice volume. The variations in contribution of the split of jet stream through MIS2–5 also demonstrate its minor role in influencing the Greenland dust concentration changes. In addition, it is suggested that the 21-ka precessional period found in the Greenland dust records possibly reflects the integrated response of paleoclimate changes to Milankovitch forcing (northern summer insolation). As for variations in the last glacial maximum (LGM) dust concentration, the enhanced dust emissions induced by changes of the North Atlantic jet stream over the Central and Eastern Europe and the Tarim Basin promoted the sharp dust concentration increases during the early-LGM. Such a scenario was intimately related to more frequent Rossby-wave breakings under negative NAO/AO phases. Furthermore, loss of dust aerosol due to wet deposition en route and ice accumulation rates could also have affected the Greenland dust concentration changes during the LGM. The findings are thought to enable interpretation of the observed decoupled relationship between the CLP loess accumulation rates and the Greenland dust concentrations during the LGM. • Systematically characterizing the Greenland ice core dust flux changes. • Provenances of the Greenland ice core dust were identified in detail. • Drivers of the last glacial dust content changes were quantitatively assessed. • Providing a full explanation of the sharply changing dust contents during the LGM. [ABSTRACT FROM AUTHOR]

The NE Tibetan Plateau (NETP) is one of the most ecologically vulnerable areas, with large expansion and retreat of desert area since the Last Glacial Maximum. However, the long-term history of desertification in this area and its mechanisms are still unknown due to lack of long-term and well-dated records. Here we present the results of sand particle content (a proxy of local materials) and detrital-zircon ages from a high-resolved loess-paleosol section spanning the past 420 ka in the Menyuan (MY) Basin lying in the NETP. Our results show that the provenance of loess units (formed in glacial periods) was strongly influenced by local desertified materials that were mainly produced by alpine glaciers surrounding the basin, which was significantly different from the paleosol units (formed in interglacial periods). The strong signals of ∼100 kyr and ∼ 20 kyr cycles in the content of local materials suggest that the desertification of the MY Basin was controlled by both the late Quaternary ∼100-kyr ice age and the ∼20-kyr precession cycles, which is obviously different from the periodicities of the classic Luochuan loess section in the mainland of the Chinese Loess Plateau (CLP). The strong ∼20-kyr cycle in the MY Basin suggest that the alpine glaciation activities in the NETP were driven by local summer insolation. These results allow us to characterize the long-term natural variability of desertification in the NETP, as well as the mechanisms for the spatial and glacial-interglacial variations of loess provenance on the CLP. • The provenance of MY loess varied over glacial-interglacial cycles. • Periodic desertification of the NETP was driven by orbital cycles. • Alpine glaciation activities in the NETP were driven by local summer insolation. • NETP is an important dust source for loess on the CLP during glacial periods. [ABSTRACT FROM AUTHOR]

The Southern Hemisphere Westerlies (SHW) and their linkages with key ocean and atmospheric processes have the potential to drive abrupt climate change. The westerlies migrate seasonally and are thought to have shifted during past climate events, such as the last glacial termination. However, the timing, magnitude, direction, and mechanisms behind such shifts remain a topic of ongoing study. Here, we contribute to the understanding of past temporal and spatial changes in the SHW by mapping their fluctuations in the South Atlantic region from the last glacial maximum (LGM) to present. We use lake sediment proxies, including plant wax isotopes extracted from two tarns on Mt. Usborne, East Falkland (51.7 °S) to infer changes in paleoclimate and in the mean annual position of the SHW over the last ∼23,000 years. Together, the proxies indicate that the position of the SHW during the LGM lay north of the Falkland Islands. We interpret our plant wax isotopic record as showing a southward migration of the mean annual position of the wind belt beginning just before ∼21 ka, accompanied by warming at Mt. Usborne at ∼16.5 ka. The late-glacial climate was variable at the field site, and the mean annual position of the SHW may have fluctuated around the latitude of the Falkland Islands. The early Holocene was characterized by relatively warm, dry conditions with high evaporation rates, and the SHW may have lain south of the islands. A brief southern excursion at 7–6 ka was followed by gradual northward migration of the SHW for the remainder of the Holocene to the pre-industrial position. This northward migration was paired with increasingly wet, cool conditions in the Falkland Islands, and the islands appear to be experiencing some of the wettest conditions of the last 23,000 years during the late Holocene. • Lake sediment proxies from East Falkland track the SHW over the last 23,000 years. • SHW moved south at ∼21 ka before reversing the trend at the start of the Holocene. • The SHW fluctuated near the Falkland Islands during the late glacial period. • The early Holocene was warm with high evaporation, and a southern SHW position. • The SHW migrated north in the second half of the Holocene. [ABSTRACT FROM AUTHOR]

Aim: The evolutionary process of an organism provides valuable data toward an understanding of the Earth evolution history. To investigate the relationship between the uplift of the Qinghai–Tibet Plateau (QTP) and mammalian evolution since the late Cenozoic, the geographic distribution of genetic variations in the Tibetan hamster Cricetulus kamensis was investigated using phylogeographical methods. In particular, population divergence, demographic history, genetic variation, and the prediction of species distribution area were investigated. Location: The Qinghai–Tibet Plateau. Methods: A total of 53 specimens, representing 13 geographic populations, were collected from the QTP. The phylogeographical pattern and demographic history of C. kamensis were analyzed, and the probable factors in the QTP uplift and the Quaternary glacial periods were inferred from one nuclear and four mitochondrial genes. Furthermore, the species distribution model (SDM) was used to predict changes in potentially suitable habitats since the last Interglacial. Results: Phylogenetic analysis demonstrated that two major genetic differentiations of the C. kamensis population occurred during the Early Pleistocene that were influenced by the Qing‐Zang tectonic movement from the Middle Pliocene to the Early Pleistocene. Genetic distance between two major clades indicated low genetic divergence. Demographic history analysis showed that the C. kamensis population was affected by the Quaternary glacial period. SDM analysis indicated that C. kamensis was endemic to the QTP and the suitable habitat was affected by climate change, especially during the Last Glacial Maximum (LGM). Main conclusion: Our results indicated that the QTP uplift led to the population divergence of C. kamensis, and vicariance well accounted for the geographic distribution of genetic variation in C. kamensis as a result of genetic divergence and lack of gene flow. The genetic distance shows that C. alticola may be a subspecies of C. kamensis. Demographic history analysis suggests that the QTP was affected by the last glacial period. SDM analysis supports that almost the entire QTP is covered by a huge ice sheet during the LGM. [ABSTRACT FROM AUTHOR]

The Indo-Pacific Warm Pool is the largest source of heat and rainfall on Earth, but the mechanisms driving long-term climate variability in this region remain uncertain. Some studies suggest that precessional variations in insolation exert a primary role; alternatively, sea level change may be more important, with exposure of the Sunda and Sahul shelves leading to weakened Walker circulation across the Indian Ocean during glacial periods. The limited number of paleoclimate records from the warm pool that extend beyond the Last Glacial Maximum makes it difficult to test these competing hypotheses. Here, we present a 450,000-year long reconstruction of sea surface temperatures, sub-surface temperatures, and regional vegetation from southern Sumatra. We show that the dominant signal in each record is the 100 ky glacial cycle. We find that both surface and sub-surface temperatures cool during glacial periods, but the sub-surface cools more, reflecting a shoaling of the thermocline. Southern Sumatra and western Java consistently exhibit increased C 4 grass coverage and an enhanced dry season during glacial periods. The observed changes in the thermocline and rainfall seasonality are consistent with the shelf exposure mechanism, suggesting that sea level changes exert a prominent role in warm pool climate over glacial-interglacial timescales. • Tropical eastern Indian Ocean thermocline shoaled during last five glacial periods. • C 4 grasses expand over exposed Sunda Shelf during glacial periods. • Enhanced rainfall seasonality supports glacial drying over the Sunda Shelf. • The 100 ky glacial cycle dominates precession-driven insolation in our records. • Sea level exerts long-term control on coupled climate in the Indo-Pacific. [ABSTRACT FROM AUTHOR]

Correlative ecological niche models are increasingly used to estimate potential distributions during the Last Glacial Maximum (LGM) for biogeographical research. In the case of presence‐background/pseudoabsences techniques, cold environments that are poorly represented in existing geography can complicate the process of model calibration and transfer into more extreme cold environments that were very common during the LGM (non‐analog conditions). This may lead to biologically unrealistic estimations. Using one cold‐adapted North American mammal, we explore a real scenario to better understand the effect of restricting the range of environmental conditions over which niche models are calibrated and then transferred to LGM conditions. We performed two sets of experiments in Maxent: 1) we calibrated models in the context of only present‐day climate conditions, which is the most common practice, and compared predictions under LGM conditions based on two extrapolation methods (clamping versus unconstrained); 2) we calibrated single models using both present‐day and LGM conditions as part of the same background in order to include more extreme environments in the model calibration. Our experiments led to dramatically different estimates of species' potential distributions, showing notable differences with respect to latitudinal and elevational shifts during the LGM. Models calibrated using present‐day climates yielded biologically unrealistic estimations, suggesting that species survived in the glaciers during the LGM. Even more unrealistic estimations were achieved when clamping was enforced as the method to extrapolate. Models calibrated in the context of both modern and past climates reduced the required degree of extrapolation and allowed more realistic potential distributions, suggesting that the species avoided extremely cold conditions during the LGM. This study alerts to the possibility of obtaining implausible potential distributions during the LGM due to restricted background datasets and offers recommendations that should promote better strategies to estimate distributional changes during glaciations. [ABSTRACT FROM AUTHOR]

Proglacial lakes can affect the stability of mountain glaciers and can partly disengage glacier behaviour from climatic perturbations. However, their role in controlling the onset and progression of deglaciation from the Last Glacial Maximum (LGM) remains poorly understood. This lack of understanding is partly because the evidence required to consistently and robustly identify the location and evolution of ice-contact lakes is not standardised. In this paper we therefore firstly present a new set of criteria for identifying the landform and sedimentary evidence that defines and characterises ice-marginal lakes. Secondly, we then apply these key criteria with the aid of high-resolution topographic mapping to produce the first holistic definition and assessment of major proglacial lake landforms and sediments pertaining to the end of the LGM across South Island, New Zealand. The major findings of this assessment can be grouped to include that: (i) The localised constraints to proglacial lake extent were topography, glacier size and meltwater/sediment fluxes, (ii) Lake damming was initiated by outwash fan-heads that interrupted water and sediment flows down-valley, and (iii) New Zealand LGM lakes were unequivocally in contact with a calving ice margin. These findings will be useful for reconstructing ice dynamics and landscape evolution in this region. • Key criteria to identify, characterise and interpret landforms and sediments associated with ice-contact proglacial lakes. • Exploitation of high-resolution LiDAR topographic data for proglacial lakes across South Island, New Zealand. • Synthesis analysis of the 3D geometry, landforms and sediments, and interpretation of the evolution of New Zealand LGM lakes. [ABSTRACT FROM AUTHOR]

Here we summarize the motivation and issues surrounding the responses of ice sheets and sea level to past climate warming as part of the PALeo constraints on SEA level rise (PALSEA) working group. Papers in this special issue of Quaternary Science Reviews focus on the timescale of glaciations during the late Pliocene, the magnitude of ice-sheet fluctuations and volume leading up to and during the last glacial maximum, the timing and persistence of ice-sheet impacts on deglacial and future relative sea-level change, and relative sea-level change during peak interglacial climate. A more dynamic cryosphere is noted under both late Pliocene and last glacial cycle climate conditions, while relative sea-level changes during the last deglaciation appear to correspond closely with individual ice-sheet deglaciation. Lastly, relative sea-level change during peak interglacial conditions may have fluctuated by as much as a meter, although the sources of such variability (Greenland, Antarctica or elsewhere) remain elusive. [ABSTRACT FROM AUTHOR]

Calculating cosmogenic-nuclide surface-exposure ages is critically dependent on a knowledge of the altitude of the sample site. Changes in altitude have occurred through time as a result of glacial isostatic adjustment (GIA), potentially altering local nuclide production rates and, therefore, surface-exposure ages. Here we assess the impact of GIA on surface-exposure dating by calculating global time-dependent production rates since the Last Glacial Maximum using surface elevations that were corrected and uncorrected for GIA. We find that the magnitude of the GIA effect is spatially and temporally variable. Nuclide production could be reduced by up to 50% in the interior of large ice masses (in North America, Scandinavia and West Antarctica) at times of maximum glacial isostatic depression. Although smaller, the effect is still significant at ice sheet margins, where nuclide production is reduced by >5% and potentially >10%, making exposure ages older in those areas. Away from the ice sheet margins, land surfaces can be isostatically elevated, which can increase nuclide production by >5% and, therefore, make exposure ages younger. Areas that were more recently exposed or that are distal to large ice masses will generally be less affected. Importantly, we find that the effect at the primary 10Be production calibration sites is <1%. Applying a GIA correction to surface-exposure data may help resolve mismatches between some chronologies, but not necessarily in all regions, implying that additional factors may need to be considered. Past atmospheric changes can amplify or reduce the impact of GIA on nuclide production, and the combined effects should be fully accounted for in the future. These time-dependent influences on surface-exposure dating have potentially important implications for interpreting chronologies and for using the data to constrain ice sheet models. • Glacial isostatic adjustment (GIA) can make exposure age estimates older or younger. • Nuclide production can be reduced by up to 50% from isostatic depression. • The GIA effect at the primary 10Be production rate calibration sites is negligible. • Correcting existing exposure datasets can make ages up to ∼2,500 years (15%) older. • Atmospheric changes may amplify or reduce the impact of GIA on exposure ages. [ABSTRACT FROM AUTHOR]

Graphical abstract Highlights • First glacial and de-glacial chronology from the Sikkim higher Himalaya. • Glaciation dates back to the MIS 2 and started at ∼28 ± 2.8 ka. • We ascribe the glaciation to lower temperature and increase moisture. Abstract The present study investigates the relics of the most extensive glacier advances and their retreat in the transitional climatic zone of the Chopta valley, northern Sikkim Himalaya with an aim to understand the forcing factors. The optically stimulated luminescence (OSL) dating of the stratigraphically oldest glacier advance and its retreat in the Chopta valley revealed that the glacier advanced at ∼28 ± 2.8 ka and sustained through the global Last Glacial Maximum (gLGM; 26.5 ka and 19 ka; Mix et al., 2001). The deglaciation is manifested by recessional moraine (ridge) dated at ∼18 ± 2.2 ka. Comparing the results of the present study with the existing glacial chronologies from the Himalayan region, it has been observed that irrespective of the geographical locations; the intensified mid latitude-westerlies during the gLGM and the associated decrease in winter temperature was the major driver of glacier advance during the MIS-2. [ABSTRACT FROM AUTHOR]

The only glaciers existing today in the Iberian Peninsula are small features located in the Pyrenees, though their number and extension has undergone significant changes over the Late Quaternary. The wide range of glacial landforms and deposits distributed across different Iberian ranges suggests the occurrence of several past periods with larger glacial systems. The objective of this research is to summarize the current knowledge on the spatial and temporal patterns of glacial activity in the Iberian mountains during the Late Quaternary. To this purpose, the chronological framework was divided in six periods: glaciations prior to the Last Glacial Cycle (Middle Pleistocene), Last Glacial Cycle (Late Pleistocene), Termination-1, Holocene, Little Ice Age (LIA) and present-day. The data were geographically divided considering the mountain systems where glacial evidence exists: Pyrenees, Cantabrian Range, NW ranges, Central Range, Iberian Range and Sierra Nevada. During Quaternary cold stages, ice accumulated in the head valleys of these mountain ranges and glaciers flowed down-valleys. In all cases, glaciers remained confined within the mountain systems and did not reach the surrounding lowlands. Depending on the combination of temperatures and moisture conditions, more or less ice was stored. In some ranges, there is evidence of Middle Pleistocene glaciations, one potentially correlating with marine isotope stage (MIS) 12 and another correlating with MIS 6 with glaciation dated to ca. 130-170 ka. However, most of the glacial records correspond to the Last Glacial Cycle and subsequent Termination. The maximum glacial expansion of this last Pleistocene glaciation stage occurred well before the global Last Glacial Maximum (LGM) between 30 and 60 ka in the Cantabrian Mountains and Pyrenees, at ca. 30 ka in Sierra Nevada and NW ranges, and (almost) synchronously to the LGM in the Central Range and Iberian Range. A massive glacial retreat occurred in all ranges at 19-20 ka, but the long-term deglaciation process was interrupted by cold intervals, such as the Oldest and Younger Dryas, which favoured glacial expansion in the highest mountains. Temperature increase recorded during the Holocene conditioned the melting of glaciers, which only reappeared in the highest massifs during the coldest periods, such as the LIA. However, post-LIA warming led to glacier disappearance in the Cantabrian Mountains, Sierra Nevada and most massifs of the Pyrenees, together with an accelerated shrinkage of the small glaciers still existing in this range at elevations near 3000 m. [ABSTRACT FROM AUTHOR]

The Islands of Four Mountains island group of the Aleutian island arc is remote and difficult to access. Consequently, little fieldwork has focused on geomorphic processes and their relationship to island morphology, climatic change, and human settlement. We investigated glacial, fluvial, and slope processes on the morphologically different Carlisle, Cleveland, and Tana volcanoes. The islands were extensively glaciated at the last glacial maximum (LGM), and there is evidence for a Neoglacial advance. On the highly dissected Tana volcano, a large basin is likely the result of a pre-LGM sector collapse and subsequent glacial erosion into weak hydrothermally altered rock. Valley and moraine morphology is also influenced by hydrothermal alteration. On both Tana and Carlisle, there are sediment fans composed dominantly of thick debris flow deposits mantled by ~3 m of layered tephras, fine-grained alluvium, and anthropogenic deposits. Debris flow deposition was favored during the unstable paraglacial landscape of the early Holocene–latest Pleistocene. The earliest direct archaeological evidence for settlement is 3.8 cal ka BP, but soil geochemical evidence suggests that the islands were inhabited by 7.3 cal ka BP. This discrepancy in the archaeological record may be explained by geomorphic processes including coastal erosion and unstable geomorphic surfaces. [ABSTRACT FROM AUTHOR]