Branched glycerol dialkyl glycerol tetraether (GDGT) distributions observed in a sediment core from Lake McKenzie were utilized to quantitatively reconstruct the pattern of mean annual air temperature (MAAT) from coastal subtropical eastern Australia between 37 and 18.3 cal ka BP and 14.0 cal ka BP to present. Both the reconstructed trend and amplitude of MAAT changes from the top of the sediment core were nearly identical to a local instrumental MAAT record from Fraser Island, providing confidence that in this sediment core branched GDGTs could be used to produce a quantitative record of past MAAT. The reconstructed trend of MAAT during 37 to 18.3 cal ka BP and timing of the Last Glacial Maximum (LGM) in the Lake McKenzie record were in agreement with previously published nearby marine climate records. The amplitude of lower-than-present MAAT during the LGM potentially provides information on the latitude of separation of the Tasman Front from the East Australian current in the subtropical western Pacific. The Lake McKenzie record shows an earlier onset of near modern day warm temperatures in the early Holocene compared to marine records and the presence of a warmer than present day period during the mid-Holocene.

Sediments containing terrestrial and aquatic mollusks and ostracodes from the Ziegler Reservoir fossil site (2705 m elevation) near Snowmass Village, Colorado, span ~130–87 ka (MIS 5e through 5b). The southeastern area of the site where taxa were recovered was a relatively fresh, shallow, well-vegetated wetland during MIS 5e through 5c time, approximately 2 m deep, with a total dissolved solids value of ~200–1000 mg L− 1. The wetland was seasonally or annually variable and groundwater discharged along the margins of the bounding moraine. Groundwater likely contributed solutes to the system and may have contributed 18O-enriched water. Based on stable isotopes from ostracode calcite (δ18OOST and δ13COST), seasonal evaporation occurred and the dissolved inorganic carbon pool was unexpectedly enriched in 13C. The mollusk and ostracode faunas changed little across the MIS 5e/5d/5c boundaries, whereas a distinct change in the ostracode fauna occurred between the deposition of Unit 11 and Unit 13, which corresponds in time to the MIS 5c/5b boundary, indicating some combination of increased surface and/or groundwater flow, a decrease in water temperature, and a freshening and a possible deepening of the wetland.

The Northern Prince Gustav Ice Stream located in Prince Gustav Channel, drained the northeastern portion of the Antarctic Peninsula Ice Sheet during the last glacial maximum. Here we present a chronology of its retreat based on in situ produced cosmogenic 10Be from erratic boulders at Cape Lachman, northern James Ross Island. Schmidt hammer testing was adopted to assess the weathering state of erratic boulders in order to better interpret excess cosmogenic 10Be from cumulative periods of pre-exposure or earlier release from the glacier. The weighted mean exposure age of five boulders based on Schmidt hammer data is 12.9 ± 1.2 ka representing the beginning of the deglaciation of lower-lying areas (< 60 m a.s.l.) of the northern James Ross Island, when Northern Prince Gustav Ice Stream split from the remaining James Ross Island ice cover. This age represents the minimum age of the transition from grounded ice stream to floating ice shelf in the middle continental shelf areas of the northern Prince Gustav Channel. The remaining ice cover located at higher elevations of northern James Ross Island retreated during the early Holocene due to gradual decay of terrestrial ice and increase of equilibrium line altitude. Schmidt hammer R-values are inversely correlated with 10Be exposure ages and could be used as a proxy for exposure history of individual granite boulders in this region and favour the hypothesis of earlier release of boulders with excessive 10Be concentrations from glacier directly at this site. These data provide evidences for an earlier deglaciation of northern James Ross Island when compared with other recently presented cosmogenic nuclide based deglaciation chronologies, but this timing coincides with rapid increase of atmospheric temperature in this marginal part of Antarctica.

Marine Oxygen Isotope Stage (MIS) 2, with its profound environmental and climatic changes from before the last glacial maximum (LGM) to the last deglaciation, is an ideal period for understanding the evolution of the East Asian summer monsoon (EASM) and Indian summer monsoon (ISM), two Asian monsoon sub-systems. With 875 stable oxygen isotope ratios and 43 230Th dates from stalagmites in Sanxing Cave, southwestern China, we construct and interpret a new, replicated, Asian summer monsoon (ASM) record covering 30.9–9.7 ka with decadal resolution. δ18O records from this site and other reported Chinese caves display similar long-term orbitally dominated trends and synchronous millennial-scale strong and weak monsoonal events associated with climate changes in high northern latitudes. Interestingly, Sanxing δ18O and Arabian Sea records show a weakening ISM from 22 to 17 ka, while the Hulu and Qingtian records from East and Central China express a 3-ka intensifying EASM from 20 to 17 ka. This decoupling between EASM and ISM may be due to different sensitivities of the two ASM sub-systems in response to internal feedback mechanisms associated with the complex geographical or land-ocean configurations.

In North America, terrestrial records of biodiversity and climate change that span Marine Oxygen Isotope Stage (MIS) 5 are rare. Where found, they provide insight into how the coupling of the ocean–atmosphere system is manifested in biotic and environmental records and how the biosphere responds to climate change. In 2010–2011, construction at Ziegler Reservoir near Snowmass Village, Colorado (USA) revealed a nearly continuous, lacustrine/wetland sedimentary sequence that preserved evidence of past plant communities between ~140 and 55 ka, including all of MIS 5. At an elevation of 2705 m, the Ziegler Reservoir fossil site also contained thousands of well-preserved bones of late Pleistocene megafauna, including mastodons, mammoths, ground sloths, horses, camels, deer, bison, black bear, coyotes, and bighorn sheep. In addition, the site contained more than 26,000 bones from at least 30 species of small animals including salamanders, otters, muskrats, minks, rabbits, beavers, frogs, lizards, snakes, fish, and birds. The combination of macro- and micro-vertebrates, invertebrates, terrestrial and aquatic plant macrofossils, a detailed pollen record, and a robust, directly dated stratigraphic framework shows that high-elevation ecosystems in the Rocky Mountains of Colorado are climatically sensitive and varied dramatically throughout MIS 5.

The ‘Mystery Interval’ (MI, 17.5−14.5 ka) was the first stage of the last deglaciation, a key interval for understanding mechanisms of glacial–interglacial cycles. To elucidate possible causes of the MI, here we present three high-resolution, precisely dated oxygen-isotope records of stalagmites from Qingtian and Hulu Caves in China, reflecting changes in the East Asian summer monsoon (EASM) then. Based on well-established chronologies using precise 230Th dates and annual-band counting results, the two-cave δ18O profiles of ~7-yr resolution match well at decadal timescales. Both of the two-cave records document an abrupt weakening (2‰ of δ18O rise within 20 yr) in the EASM at ~16.1 ka, coinciding with the transition of the two-phased MI reconstructed from New Mexico's Lake Estancia. Our results indicate that the maximum southward displacement of the Intertropical Convergence Zone and associated southward shift of polar jet stream may generate this two-phase feature of the MI during that time. We also discover a linear relationship among decreasing EASM intensity, rising atmospheric CO2 and weakening Atlantic Meridional Overturning Circulation between the MI and Younger Dryas episodes, suggesting a strong coupling of atmospheric/oceanic circulations in response to the millennial-scale forcing, which in turn regulates global climate changes and carbon cycles.

Relatively few radiometrically dated records are available for the central Mediterranean spanning the marine oxygen isotope stage 6–5 (MIS 6–5) transition and the first part of the Last Interglacial. Two flowstone cores from Tana che Urla Cave (TCU, central Italy), constrained by 19 U/Th ages, preserve an interval of continuous speleothem deposition between ca. 159 and 121 ka. A multiproxy record (δ18O, δ13C, growth rate and petrographic changes) obtained from this flowstone preserves significant regional-scale hydrological changes through the glacial/interglacial transition and multi-centennial variability (interpreted as alternations between wetter and drier periods) within both glacial and interglacial stages. The glacial stage shows a wetter period between ca. 154 and 152 ka, while the early to middle Last Interglacial period shows several drying events at ca. 129, 126 and 122 ka, which can be placed in the wider context of climatic instability emerging from North Atlantic marine and NW European terrestrial records. The TCU record also provides important insights into the evolution of local environmental conditions (i.e. soil development) in response to regional and global-scale climate events.

Geological maps of South Carolina, covering >6800 km2, confirm the existence of eight preserved Pleistocene shorelines above current sea level: Marietta (+ 42.6 m), Wicomico (+ 27.4 m), Penholoway (+ 21.3 m), Ladson (+ 17.4 m), Ten Mile Hill (+ 10.7 m), Pamlico (+ 6.7 m), Princess Anne (+ 5.2 m), and Silver Bluff (+ 3 m). Current geochronologic data suggest that these eight shorelines correlate with Marine Oxygen Isotope Stages (MIS) as follows: Marietta—older than MIS 77; Wicomico—MIS 55–45; Penholoway—MIS 19 or 17; Ladson—MIS 11; Ten Mile Hill—MIS 7; Pamlico—MIS 5; Princess Anne—MIS 5; and Silver Bluff—MIS 5 or 3. Except for the MIS 5e Pamlico, and possibly the MIS 11 Ladson, the South Carolina elevations are higher than predicted by isotope proxy-based reconstructions. The <4 m of total relief from the Pamlico to the Silver Bluff shoreline in South Carolina, for which other reconstructions suggest an expected relief of ~ 80 m, illustrates the lack of match. Our results suggest that processes affecting either post-depositional changes in shoreline elevations or the creation of proxy sea-level estimates must be considered before using paleo sea-level position on continental margins.

Quaternary post-Barreiras sediments are widespread along Brazil's passive margin. These deposits are well exposed in the onshore Paraíba Basin, which is one of the rift basins formed during the Pangean continental breakup. In this area, the post-Barreiras sediments consist of sandstones with abundant soft-sediment deformation structures related to seismicity contemporaneous with deposition. The trace fossils Thalassinoides and Psilonichnus are found up to 38 m above modern sea level in sandstones dated between 60.0 (± 1.4) and 15.1 (± 1.8) ka. The integration of ichnological and sedimentary facies suggests nearshore paleoenvironments. Such deposits could not be related to eustatic sea-level rise, as this time coincides with the last glaciation. Hence, an uplift of 0.63 mm/yr, or 1.97 mm/yr if sea level was 80 m lower in the last glaciation, would have been required to ascend the post-Barreiras sediments several meters above the present-day sea level during the last 60 ka. This would suggest that the post-rift stage of the South American eastern passive margin may have experienced tectonic reactivation more intense than generally recognized. Although more complete data are still needed, the information presented herein may play an important role in studies aiming to decipher the Quaternary evolution of this passive margin.

Accurate sea-level reconstruction is critical in understanding the drivers of coastal evolution. Inliers of shallow marine limestone and aeolianite are exposed as zeugen (carbonate-capped erosional remnants) on the southern coast of the Arabian/Persian Gulf. These have generally been accepted as evidence of a eustatically driven, last-interglacial relative sea-level highstand preceded by a penultimate glacial-age lowstand. Instead, recent optically stimulated luminescence (OSL) dating suggests a last glacial age for these deposits, requiring >100 m of uplift since the last glacial maximum in order to keep pace with eustatic sea-level rise and implying the need for a wholesale revision of tectonic, stratigraphic and sea-level histories of the Gulf. These two hypotheses have radically different implications for regional neotectonics and land–sea distribution histories. Here we test these hypotheses using OSL dating of the zeugen formations. These new ages are remarkably consistent with earlier interpretations of the formations being last interglacial or older in age, showing that tectonic movements are negligible and eustatic sea-level variations are responsible for local sea-level changes in the Gulf. The cause of the large age differences between recent studies is unclear, although it appears related to large differences in the measured accumulated dose in different OSL samples.

Thirty-six new and previously published radiocarbon dates constrain the relative sea-level history of Arviat on the west coast of Hudson Bay. As a result of glacial isostatic adjustment (GIA) following deglaciation, sea level fell rapidly from a high-stand of nearly 170 m elevation just after 8000 cal yr BP to 60 m elevation by the mid Holocene (~ 5200 cal yr BP). The rate of sea-level fall decreased in the mid and late Holocene, with sea level falling 30 m since 3000 cal yr BP. Several late Holocene sea-level measurements are interpreted to originate from the upper end of the tidal range and place tight constraints on sea level. A preliminary measurement of present-day vertical land motion obtained by repeat Global Positioning System (GPS) occupations indicates ongoing crustal uplift at Arviat of 9.3 ± 1.5 mm/yr, in close agreement with the crustal uplift rate inferred from the inferred sea-level curve. Predictions of numerical GIA models indicate that the new sea-level curve is best fit by a Laurentide Ice Sheet reconstruction with a last glacial maximum peak thickness of ~ 3.4 km. This is a 30–35% thickness reduction of the ICE-5G ice-sheet history west of Hudson Bay.

Comparisons of historical and modern assemblages of mammals can yield important insights into patterns and processes of environmental change. Here, we present the first analyses of small mammal assemblages present in northern Patagonia during the last glacial maximum (LGM). Using remains obtained from owl pellets excavated from an archeological cave site (Arroyo Corral I, levels VII–V, carbon dates of 22,400–21,530 cal yr BP), we generate estimates of the minimum number of individuals for all species detected; these estimates, in turn are used to determine relative species abundances. Comparisons of these data with similar analyses of small mammal remains obtained from a second archeological site (ACoII, levels IV–V, carbon dates of 10,010–9220 cal yr BP) as well as from modern owl pellets reveal pronounced changes in relative species abundance since the LGM. In particular, Euneomys chinchilloides and Ctenomys sociabilis – the predominant species during the LGM – declined markedly, suggesting a change from open, bare habitat punctuated by patches of wet meadows and shrubs to the more densely vegetated mosaic of ecotone habitats found in this region today. These data provide important new insights into the environmental changes that have occurred in northern Patagonia over the last 20,000 years.

This work presents the first reconstruction of late Pleistocene glacier fluctuations on Uturuncu volcano, in the Southern Tropical Andes. Cosmogenic 3He dating of glacial landforms provides constraints on ancient glacier position between 65 and 14 ka. Despite important scatter in the exposure ages on the oldest moraines, probably resulting from pre-exposure, these 3He data constrain the timing of the moraine deposits and subsequent glacier recessions: the Uturuncu glacier may have reached its maximum extent much before the global LGM, maybe as early as 65 ka, with an equilibrium line altitude (ELA) at 5280 m. Then, the glacier remained close to its maximum position, with a main stillstand identified around 40 ka, and another one between 35 and 17 ka, followed by a limited recession at 17 ka. Then, another glacial stillstand is identified upstream during the late glacial period, probably between 16 and 14 ka, with an ELA standing at 5350 m. This stillstand is synchronous with the paleolake Tauca highstand. This result indicates that this regionally wet and cold episode, during the Heinrich 1 event, also impacted the Southern Altiplano. The ELA rose above 5450 m after 14 ka, synchronously with the Bolling–Allerod.

Although there is evidence from molecular studies for the existence of central European last glacial refugia for temperate species, there is still a great lack of direct fossil records to confirm this theory. Here we bring such evidence in the form of fossil shells from twenty strictly forest land snail species, which were recorded in radiocarbon-dated late glacial or older mollusc assemblages of nine non-interrupted mollusc successions situated in the Western Carpathians, and one in the Bohemian Massif. We proposed that molluscs survived the last glacial period in central Europe in isolated small patches of broadleaf forest, which we unequivocally demonstrate for two sites of last glacial maximum age.

Little is known concerning climate changes in the Eastern Baltic region during the last interglacial–glacial cycle and in particular, climate changes during the Weichselian. In this study, a quantitative reconstruction of the mean January and July temperature for the Medininkai-117 site in Lithuania is presented. The reconstruction is based on pollen and plant macrofossils from this site, which reveal that the vegetation was characteristic of many northern Europe sites during the Eemian and Early Weichselian. Gradual evolution of the vegetation suggests that relatively uniform climate conditions existed during the Eemian. Our reconstructions support the view of a relatively stable Eemian, with short cooling phases of low amplitude. A strong increase in temperature was apparent during the beginning of the interglacial and decrease during the transition to the Weichselian. Reconstructed July temperatures of the Eemian interglacial were approximately 2 °C higher than today (18.5 °C; today: 16.2 °C) and were similar to today for January (−5.2 °C; today: −5.1 °C). July temperatures during the Early Weichselian were only ~2°C lower than during the Eemian, whereas the January temperatures gradually decreased. Winter temperatures were relatively high (above − 10 °C) during the Early Weichselian.

Dansgaard–Oeschger (D–O) cycles had far-reaching effects on Northern Hemisphere and tropical climate systems during the last glacial period, yet the climatic response to D–O cycles in western North America is controversial, especially prior to 55 ka. We document changes in precipitation along the western slope of the central Sierra Nevada during early Marine Oxygen Isotope Stages (MIS) 3 and 4 (55–67 ka) from a U-series dated speleothem record from McLean's Cave. The timing of our multi-proxy geochemical dataset is coeval with D–O interstadials (15–18) and stadials, including Heinrich Event 6. The McLean's Cave stalagmite indicates warmer and drier conditions during Greenland interstadials (GISs 15–18), signified by elevated δ18O, δ13C, reflectance, and trace element concentrations, and less radiogenic 87Sr/86Sr. Our record extends evidence of a strong linkage between high-latitude warming and reduced precipitation in western North America to early MIS 3 and MIS 4. This record shows that the linkage persists in diverse global climate states, and documents the nature of the climatic response in central California to Heinrich Event 6.

Sediment-based reconstructions of late-Quaternary lake levels provide direct evidence of hydrologic responses to climate change, but many studies only provide approximate lake-elevation curves. Here, we demonstrate a new method for producing quantitative time series of lake elevation based on the facies and elevations of multiple cores collected from a lake's margin. The approach determines the facies represented in each core using diagnostic data, such as sand content, and then compares the results across cores to determine the elevation of the littoral zone over time. By applying the approach computationally, decisions are made systematically and iteratively using different facies classification schemes to evaluate the associated uncertainty. After evaluating our assumptions using ground-penetrating radar (GPR), we quantify past lake-elevation changes, precipitation minus evapotranspiration (ΔP−ET), and uncertainty in both at Lake of the Woods and Little Windy Hill Pond, Wyoming. The well-correlated (r = 0.802 ± 0.002) reconstructions indicate that water levels at both lakes fell at >11,300, 8000–5500, and 4700–1600 cal yr BP when ΔP − ET decreased to −50 to −250 mm/yr. Differences between the reconstructions are typically small (10 ± 24 mm/yr since 7000 cal yr BP), and the similarity indicates that our reconstruction method can produce statistically comparable paleohydrologic datasets across networks of sites.

Quantitative analysis of macroecological patterns for late Pleistocene assemblages can be useful for disentangling the causes of late Quaternary extinctions (LQE). However, previous analyses have usually assumed linear relationships between macroecological traits, such as body size and range size/range shift, that may have led to erroneous interpretations. Here, we analyzed mammalian datasets to show how macroecological patterns support climate change as an important driver of the LQE, which is contrary to previous analyses that did not account for more complex relationships among traits. We employed quantile regression methods that allow a detailed and fine-tuned quantitative analysis of complex macroecological patterns revealed as polygonal relationships (i.e., constraint envelopes). We showed that these triangular-shaped envelopes that describe the macroecological relationship between body size and geographical range shift reflect nonrandom extinction processes under which the large-bodied species are more prone to extinction during events of severe habitat loss, such as glacial/interglacial transitions. Hence, we provide both a theoretical background and methodological framework to better understand how climate change induces body size-biased species sorting and shapes complex macroecological patterns.

Several Gigantopithecus faunas associated with taxonomically undetermined hominoid fossils and/or stone artifacts are known from southern China. These faunas are particularly important for the study of the evolution of humans and other mammals in Asia. However, the geochronology of the Gigantopithecus faunas remains uncertain. In order to solve this problem, a program of geochronological studies of Gigantopithecus faunas in Guangxi Province was recently initiated. Chuifeng Cave is the first studied site, which yielded 92 Gigantopithecus blacki teeth associated with numerous other mammalian fossils. We carried out combined ESR/U-series dating of fossil teeth and sediment paleomagnetic studies. Our ESR results suggest that the lower layers at this cave can be dated to 1.92 ± 0.14 Ma and the upper layers can be dated to older than 1.38 ± 0.17 Ma. Correlation of the recognized magnetozones to the geomagnetic polarity timescale was achieved by combining magnetostratigraphic, biostratigraphic and ESR data. The combined chronologies establish an Olduvai subchron (1.945–1.778 Ma) for the lowermost Chuifeng Cave sediments. We also analyzed the enamel δ13C values of the Gigantopithecus faunas. Our results show that southern China was dominated by C3 plants during the early Pleistocene and that the Gigantopithecus faunas lived in a woodland-forest ecosystem.

The Dead Sea fault (DSF) is one of the most active plate boundaries in the world. Understanding the Quaternary history and sediments of the DSF requires investigation into the Neogene development of this plate boundary. DSF lateral motion preceded significant extension and rift morphology by ~10 Ma. Sediments of the Sedom Formation, dated here between 5.0 ± 0.5 Ma and 6.2−2.1+inf Ma, yielded extremely low 10Be concentrations and 26Al is absent. These reflect the antiquity of the sediments, deposited in the Sedom Lagoon, which evolved in a subdued landscape and was connected to the Mediterranean Sea. The base of the overlying Amora Formation, deposited in the terminal Amora Lake which developed under increasing relief that promoted escarpment incision, was dated at 3.3−0.8+0.9 Ma. Burial ages of fluvial sediments within caves (3.4 ± 0.2 Ma and 3.6 ± 0.4 Ma) represent the timing of initial incision. Initial DSF topography coincides with the earliest Red Sea MORB's and the East Anatolian fault initiation. These suggest a change in the relative Arabian–African plate motion. This change introduced the rifting component to the DSF followed by a significant subsidence, margin uplift, and a reorganization of relief and drainage pattern in the region resulting in the topographic framework observed today.