Oral traditions describing details of ancient volcanic eruptions and their effects survive throughout the inhabited world. Many such eruptions, especially those having catastrophic environmental and societal consequences, proved sufficiently memorable to form the basis of enduring oral traditions. Using global databases, we identified 2306 such eruptions from 477 inhabited locations that occurred before the start of the Common Era (CE) and are therefore likely to have been the subject of oral traditions. Of these, we selected 20 events (‘remembered’ Holocene eruptions) for which there are extant oral (-derived) traditions that demonstrate how such traditions can reveal details of past volcanism that often are undetectable by retrodictive geoscientific enquiry. We also selected 20 events (‘forgotten’ Holocene eruptions) about which no oral traditions are known and discuss the possible reasons for this. Such oral traditions, while often challenging for conventionally trained geoscientists to interpret, are valuable yet largely overlooked sources of information about the nature and effects of Holocene volcanism that can usefully complement geoscientific enquiry. In particular, we identified locations where memories of such volcanism appear ‘forgotten’ in the hope that scientists might focus their attention on revealing, identifying, and analyzing local traditions.

The fluvial capture of endorheic basins represents a milestone in basin chronology, implying a profound disequilibrium that triggers critical geomorphological, sedimentological, paleogeographic, and even paleoecological transformations. The primary goal of many geomorphological studies is to determine the timing of endorheic-to-exorheic transitions with the objective of unveiling the dynamics that follow the capture event. The age of the Guadix-Baza Basin capture in the Central Betic Cordillera (S Spain) remains a subject of controversy, with proposed estimates ranging from 17 to 600 ka. In this study, we present new 234U/230Th and optically stimulated luminescence ages from exorheic deposits exposed within the basin's main fluvial valley, the Guadiana Menor River. We acquired the oldest numerical age recorded to date for a postcapture deposit within the basin. This age corresponds to a travertine platform formed 240.8 ± 25 ka on a surface level that was already incised into the glacis surface at approximately 250 m. Using these data, we estimate that basin capture took place earlier than ca. 240 ka, plus the time required for the river to incise 250 m to the position of the travertine. Furthermore, the proximity of the Matuyama-Brunhes reversal (781 ka) to the top of the endorheic succession and the ages of the paleontological sites (> ca. 750 ka) throughout the basin suggest that the capture could have occurred earlier than the oldest previously proposed age of 600 ka.

This study explores patterning in δ18O values of tooth enamel in contemporary African herbivores from mainly C3-dominated ecosystems. Evapotranspiration causes plants to lose H216O to a greater extent than H218O, leaving leaves enriched in 18O. In eastern Africa, ES species (evaporation-sensitive species: those obtaining water from food) tend to have more positive δ18Oenamel values than EI species (evaporation-insensitive species: those heavily dependent on drinking water); the magnitude of the difference increases with increasing aridity. We find the same pattern applies in the winter and year-round rainfall region of southern Africa, allowing us to use δ18Oenamel in fossil animals to examine paleo-aridity. We apply this approach to infer aridity at Quaternary fossil assemblages from present-day winter and year-round rainfall zones, including Elandsfontein (ca. 1–0.6 Ma), Hoedjiespunt (ca. 300–130 ka), and Nelson Bay Cave (23.5–3 ka). This analysis suggests that (1) at various times during the Pleistocene, Elandsfontein and Hoedjiespunt environments were wetter than last glacial maximum (LGM) to Holocene environments at Nelson Bay Cave (year-round rainfall zone); and (2) considered alongside other evidence from the year-round rainfall zone, wetter conditions across the Pleistocene–Holocene transition at Nelson Bay Cave suggests that climate changes at near-coastal sites may be out of phase with the adjacent interior.

Geotechnical drilling for a tunnel between Port Moody and Burnaby, BC, Canada, uncovered a buried fjord. Its sedimentary fill has a thickness of at least 130 m and extends more than 37 m below present mean sea level. Recovered sediments record cyclical growth and decay of successive Cordilleran ice sheets. The oldest sediments comprise 58 m of almost stoneless silt conformably overlying ice-proximal sediments and till, which in turn overlie bedrock. These sediments may predate Marine Isotope Stage (MIS) 4. Glacial sediments assigned to MIS 4 overlie this basal succession and, in turn, are overlain by MIS 3 interstadial sediments and sediments from two MIS 2 glacial advances. Indicators of relative sea-level elevations that bracket glacial deposits of MIS 4 and 2 indicate the cyclic existence of moat-like isostatic depressions in the front of expanding ice sheets. Compared with present sea level, these depressions were at least 160 m during the onsets of MIS 4 and MIS 2. Assuming a maximum eustatic drawdown of 120 m during MIS 2, isostatic depression may have exceeded 200 m during retreat of glacial ice from the Evergreen tunnel area. This is consistent with region-specific low mantle viscosity and rapid Cordilleran Ice Sheet buildup and wasting.

Based on calibrated radiocarbon ages of terrestrial gastropod shells (Succineidae, Discus, Stenotrema, Webbhelix), the chronology of Peoria Silt (loess) deposition in the Central Lowlands is updated. These taxa provide reliable ages (within ~0.2 ka), based on historical shell dating, shell-organic age comparisons, and stratigraphic consistency. A compilation of 53 new and 36 published Peoria Silt shell ages (calibrated), from 12 localities, date from 30.0 to 17.4 ka. Proximal (fossiliferous) loess from 10 sections had mean loess accumulation rates of 0.6–2.2 mm/yr. Study sites along the upper Mississippi, Illinois, to mid-Mississippi, and Ohio-Wabash Valleys suggest Peoria loess accumulated from ~27 to 15 ka, ~29 to 18 ka, and ~30 to 18 ka, respectively. The cessation age for Peoria Silt, based on surface extrapolations, is ~1–6 ka earlier than some prior Illinois estimates, even assuming slower loess accumulation in the modern solum. Younger loess in northwestern Illinois likely reflects, in part, Superior and Des Moines Lobe glacial-meltwater sediment, and Iowan Erosion Surface inputs to the upper Mississippi Valley, after the Lake Michigan Lobe receded. Furthermore, stronger winds, drier conditions, and reduced vegetation cover in valley deflation areas may have favored higher accumulation rates and later loess deposition in northwestern relative to southeastern areas.

High-resolution sedimentological and geochronological investigations of paleochannel systems in the Ried Central d'Alsace (northeastern France) allow for the reconstruction of the late glacial and Holocene fluvial evolution of this section of the Upper Rhine alluvial plain. During the Oldest Dryas, the landscape featured a dominant braided Rhine system and, to a lesser extent, a braided Fecht system. The shift to the Bølling-Allerød saw a narrowing of the Rhine's active channel belt, the development of a complex channel pattern, and the genesis of the Ill River. The river channel patterns remained unchanged during the Younger Dryas. In the Early Holocene, the Rhine's active belt narrowed further, and the Rhine and Ill Rivers developed braided-anastomosing and anastomosing channel patterns, respectively. Throughout the Holocene, both rivers maintained their channel patterns while migrating east and west across the alluvial plain, respectively. In the late glacial, fluvial dynamics in this section of the Upper Rhine plain were primarily influenced by climate-related environmental and hydrogeomorphological changes. Conversely, during the Holocene, the evolution of the fluvial hydrosystems was driven by a complex interaction of climatic and non-climatic factors, including human activity at the catchment scale, alluvial plain architecture, and local neotectonics.

Reconstructing Pleistocene glaciation timing and extent is vital for understanding paleoclimate. Whereas late Pleistocene glaciation has been studied extensively in western North American mountain ranges, the glacial history of the western Elk Range in Colorado remains understudied, particularly in the East River watershed, a site of intense scientific focus. Here we use cosmogenic nuclide exposure and depth–profile dating methods to determine the timing of glaciation in the East River watershed. We use glacier modeling to reconstruct paleoglacier extents and quantify past climate conditions. Our findings indicate that the East River glacier retreated from its maximum position approximately 17–18 ka, moving to recessional positions between 13 and 15 ka, before experiencing more substantial retreat to high-elevation cirques around 13 ka. Glacier modeling suggests that the maximum ice extents at 17–18 ka could have been sustained by temperature depressions of approximately −6.5°C compared to modern conditions, assuming consistent precipitation. Additionally, the ice position at 13–15 ka could have been supported by temperature depressions of around −4.0°C. These results offer insights into the deglaciation timeline in the East River watershed and broader western Elk Range as well as paleoclimate conditions during the late Pleistocene, which may aid future research on critical zone evolution in the East River watershed.

The role of human hunting behavior versus climate change in the mass extinction of megafauna during the Late Quaternary is much debated. To move beyond monocausal arguments, we treat human–megafauna–environment relationships as social–ecological systems from a complex adaptive systems perspective, to create an agent-based model that tests how human hunting may interact with environmental stress and animal life history to affect the probability of extinction. Using the extinction of Syncerus antiquus in South Africa at 12–10 ka as a loose inspirational case study, we parameterized a set of experiments to identify cross-feedbacks among environmental dynamics, prey life history, and human hunting pressure that affect extinction probability in a non-linear way. An important anthropogenic boundary condition emerges when hunting strategies interrupt prey animal breeding cycles. This effect is amplified in patchy, highly seasonal environments to increase the chances of extinction. This modeling approach to human behavior and biodiversity loss helps us understand how these types of cross-feedback effects and boundary conditions emerge as system components interact and change. We argue that this approach can help translate archaeological data and insight about past extinction for use in understanding and combating the current mass extinction crisis.

We present new stable oxygen and carbon isotope composite records (δ18O, δ13C) of speleothems from Sandkraal Cave 1 (SK1) on the South African south coast for the time interval between 104 and 18 ka (with a hiatus between 48 and 41 ka). Statistical comparisons using kernel-based correlation analyses and semblance analyses based on continuous wavelet transforms inform the relationships of the new speleothem records to other proxies and their changes through time. Between 105 and ~70 ka, changes of speleothem δ18O values at SK1 are likely related to rainfall seasonality. Variations of δ13C values are associated with changes of vegetation density, prior carbonate precipitation (PCP), CO2 degassing in the cave, and possibly variations of the abundance of C3 and C4 grasses in the vegetation. The relationships of δ18O with other proxies shift between ~70 and 48 ka (Marine Isotope Stages 4–3) so that both stable isotope records now reflect CO2 degassing, evaporation, and PCP. Similar relationships also continue after the hiatus for the deposition phase between 42 and 18 ka. Our findings support modeling results suggesting drier conditions in the study area when the Southern Hemisphere westerlies are shifted north and the paleo–Agulhas Plain is exposed.

Mining operations in the Enmynveem valley, northeastern Siberia, exposed a well-preserved right hind leg of Mammuthus primigenius (woolly mammoth), dated to ca. 37,500 cal yr BP. The leg had a fracture that crosscut the midsections of the tibia and fibula. Additional skeletal and soft tissue remains, including two mummified adults (Berezovka, ca. 47,200 cal yr BP; Bolshoi Lyakhovsky, ca. 37,000 cal yr BP), document the presence of mammoths in interior mountain valleys and across both northern and southern coasts of far northeastern Siberia during Marine Isotope Stage (MIS) 3. A mosaic of herb-dominated tundra communities characterized the vegetation of the Enmynveem site during late to middle MIS 3 and MIS 2 (ca. 37,000–17,000 cal yr BP). Shrubs were limited to Salix during the late Pleistocene, whereas Betula also may have been present in sheltered sites during MIS 3. Herb communities remained dominant during the late Pleistocene–Early Holocene transition, although shrub Betula increased during this interval. By ca. 10,200 cal yr BP, the vegetation was Betula–Alnus shrub tundra. Larix and Pinus pumila were established in the valley by ca. 8700 cal yr BP and ca. 5700 cal yr BP, respectively.

Recent research has shown the potential of speleothem δ13C to record a range of environmental processes. Here, we report on 230Th-dated stalagmite δ13C records for southwest Sulawesi, Indonesia, over the last 40,000 yr to investigate the relationship between tropical vegetation productivity and atmospheric methane concentrations. We demonstrate that the Sulawesi stalagmite δ13C record is driven by changes in vegetation productivity and soil respiration and explore the link between soil respiration and tropical methane emissions using HadCM3 and the Sheffield Dynamic Global Vegetation Model. The model indicates that changes in soil respiration are primarily driven by changes in temperature and CO2, in line with our interpretation of stalagmite δ13C. In turn, modelled methane emissions are driven by soil respiration, providing a mechanism that links methane to stalagmite δ13C. This relationship is particularly strong during the last glaciation, indicating a key role for the tropics in controlling atmospheric methane when emissions from high-latitude boreal wetlands were suppressed. With further investigation, the link between δ13C in stalagmites and tropical methane could provide a low-latitude proxy complementary to polar ice core records to improve our understanding of the glacial–interglacial methane budget.

Lamalunga Cave (Altamura, Southern Italy) is renowned for the discovery in 1993 of an excellently preserved Neanderthal skeleton. Given the importance of the findings and the potential use of Lamalunga speleothems for paleoclimate reconstructions, a detailed monitoring program was undertaken to investigate the connections between microclimate parameters and speleothem growth. The cave air temperature is characterized by annual sinusoidal cycles with increasing phase shift and decreasing thermal amplitude from ± 2.1°C to ± 0.04°C as a function of increasing rock overburden, and daily to sub-daily cycles induced by surface air pressure (SAP) variations characterized by strong 24-hour and 12-hour solar harmonic frequencies, with thermal amplitude decreasing from 0.0054°C to 0.0021°C in the deeper parts of the cave. The ventilation regime is mainly controlled by SAP fluctuations. Fast SAP rises can trigger “emptying events” during which most of the cave air is replaced and the CO2 concentration falls towards near-atmospheric values. The steady and gentle ventilation created by SAP fluctuations also influences the growth and morphology of calcite coralloids, with larger popcorn assemblages concentrated in the northern galleries and delicate branching morphologies prevailing along the South Gallery where stable temperature and relative humidity conditions occur. The study is a seminal example of how high-precision measurements of cave air temperature can provide a wealth of information on cave ventilation and thermal regime and provide valuable support for robust paleoclimate reconstruction from speleothems.

The depositional history of the Bonneville Salt Flats, a perennial saline pan in Utah's Bonneville basin, has poor temporal constraints, and the climatic and geomorphic conditions that led to saline pan formation there are poorly understood. We explore the late Pleistocene to Holocene depositional record of Bonneville Salt Flats cores. Our data challenge the assumption that the saline pan formed from the desiccation of Lake Bonneville, the largest late Pleistocene lake in the Great Basin, which covered this area from 30 to 13 cal ka BP. We test two hypotheses: whether climatic transitions from (1) wet to arid or (2) arid to wet led to saline pan deposition. We describe the depositional record with radiocarbon dating, sedimentological structures, mineralogy, diatom, ostracode, and portable X-ray fluorescence spectrometer measurements. Gypsum and carbonate strontium isotope ratio measurements reflect changes in water sources. Three shallow saline lake to desiccation cycles occurred from >45 and >28 cal ka BP. Deflation removed Lake Bonneville sediments between 13 and 8.3 cal ka BP. Gypsum deposition spanned 8.3 to 5.4 cal ka BP, while the oldest halite interval formed from 5.4 to 3.5 cal ka BP during a wetter period. These findings offer valuable insights for sedimentologists, archaeologists, geomorphologists, and land managers.

Predicting coastal change depends upon our knowledge of postglacial relative sea-level variability, partly controlled by glacio-isostatic responses to ice-sheet melting. Here, we reconstruct the postglacial relative sea-level changes along the Caribbean and Pacific coasts of northwestern South America by numerically solving the sea-level equation with two scenarios of mantle viscosity: global standard average and high viscosity. Our results with the standard model (applicable to the Pacific coast) agree with earlier studies by indicating a mid-Northgrippian high stand of ~2 m. The high-viscosity simulation (relevant to the Caribbean coast) shows that the transition from far- to intermediate-field influence of the Laurentide Ice Sheet occurs between Manzanillo del Mar and the Gulf of Morrosquillo. South of this location, the Colombian Caribbean coast has exhibited a still stand with a nearly constant Holocene relative sea level. By analyzing our simulations considering sea-level indicators, we argue that tectonics is more prominent than previously assumed, especially along the Caribbean coast. This influence prevents a simplified view of regional relative sea-level changes on the northwestern South American coast.

Loess paleoenvironmental reconstructions on regional to supra-regional scales have recently gained much attention. Geochemistry comparisons in relation to reference datasets, such as the Upper Continental Crust (UCC) data, have furthered our understanding of the climatic and geomorphological conditions under which terrestrial sites have developed. However, UCC data differs from loess, thereby obscuring important features, and the existing “average loess” datasets also are not sufficient for modern investigations.

In this study, we examine the youngest Polish loess (L1LL1 = MIS 2, ca. 26–15 ka) for its suitability as a new, loess-focused reference dataset. Eighty-nine samples from seven sites were analyzed, using inductively coupled plasma spectrometry. The loess had assumedly been homogenized during transportation and/or sedimentary recycling (LaN/SmN = 3.34–4.06, median 3.78; Eu/Eu* = 0.46–0.66, median 0.55; GdN/YbN = 1.08–1.49, median 1.26), and weakly affected by pre- or post-depositional weathering (CIA = 53.64–69.12, median 57.69). The statistically significant differences between sites in elemental medians were mostly conditioned by variations in grain size and in the “fresh” to “re-deposited” sediment ratio. Nonetheless, the overall geochemical composition homogeneity provided a basis for the estimation of Polish Median Loess (PML) data, as determined for 41 chemical elements. When used, PML data highlight differences between loess regions in Europe, thereby providing a tool for cross-continental comparisons.

Connecting changes in erosion and vegetation is necessary for predicting topographic and ecologic change in thawing permafrost landscapes. Formerly periglacial landscapes serve as potential analogs for understanding modern permafrost landscape change, yet compared to paleoenvironmental records at these sites, less is known about concurrent geomorphic processes, particularly their rates and relationships to climate change. Here, we target sediments preserved in a central Appalachian peat bog to reconstruct sedimentation across the last deglacial warming. We use ground-penetrating radar and geochemistry of cored bog sediments to quantify sedimentation timing, style, and provenance. Using 14C dating of sedimentary and geochemical shifts, we connect depositional changes to global climate and local vegetation change. We show that deglacial warming promoted deep soil disturbances via solifluction at ca. 14 ka. In contrast, relatively wetter conditions from ca. 10–9 ka promoted shallow disturbance of hillslopes via slopewash, which corresponds to a time of vegetation change. Our results highlight climate-modulated erosion depth and processes in periglacial and post-periglacial landscapes. The existence of similar erosion and vegetation records preserved regionally implies these dynamics were pervasive across unglaciated Appalachian highlands, aiding in reconstructing erosion responses to warming at a resolution with implications for predicting high-latitude landscape responses to disturbance.

Cave carbonate mineral deposits (speleothems) contain trace elements that are intensively investigated for their significance as palaeoclimate and environmental proxies. However, chlorine, which is abundant in marine and meteoric waters, has been overlooked as a potential palaeo-proxy, while cosmogenic 36Cl could, in principle, provide a solar irradiance proxy. Here, total Cl concentrations analysed from various speleothems were low (3–14 mg/kg), with variations linked to crystal fabrics. High-resolution synchrotron radiation micro X-ray fluorescence (μ-XRF) trace element mapping showed Cl often associated with Na, Si, and Al. We propose that speleothems incorporate Cl in two fractions: (1) water soluble (e.g., fluid inclusions) and (2) water insoluble and strongly bound (e.g., associated with detrital particulates). However, disparities indicated that alternate unidentified mechanisms for Cl incorporation were present, raising important questions regarding incorporation of many trace elements into speleothems. Our first measurements of 36Cl/Cl ratios in speleothems required large samples due to low Cl concentrations, limiting the potential of 36Cl as a solar irradiance proxy. Critically, our findings highlight a knowledge gap into how Cl and other trace elements are incorporated into speleothems, how the incorporation mechanisms and final elemental concentrations are related to speleothem fabrics, and the significance this may have for how trace elements in speleothems are interpreted as palaeoclimate proxies.

This study presents the results of a detailed paleoclimate investigation on stalagmite YL-1 (Yelini Cave, Türkiye). YL-1 grew between 117.13 (+0.57/−0.44) ka and 114.87 (+1.63/−2.89) ka within Greenland Stadial 26, indicating a positive moisture balance during the stadial conditions in this semi-arid region. Rainfall is significantly affected by sub-cloud and surface evaporation and decreasing net effective winter precipitation is recorded by high isotope values. Enriched δ18O and δ13C at 116.65 (+0.51/−0.39) ka are interpreted as a drought event that took place ca. 400 years before the end of the MIS 5e. This event, which was reported simultaneously in marine and terrestrial archives in the Northern Hemisphere, is a result of decreased cyclone activity linked to weakening of the Atlantic meridional overturning circulation. 87Sr/86Sr values of YL-1 are close to the host-rock values. Decreased 87Sr/86Sr ratio at 116.67 (+0.54/−0.38) ka reflects the intensified water–rock interaction due to lower precipitation. Along with prior calcite precipitation effect, this is also observed by increased Mg/Ca and Sr/Ca, while low P, Cu, Be, Y, and Zr concentrations indicate a lowered amount of soil-derived colloidal material. The MIS 5e/5d transition is marked by reduced insolation and enriched δ18O at 116.24 (+0.53/−0.86) ka. The Greenland Interstadial 25 phase at 115.87 (+0.83/−1.71) ka is represented by more negative δ18O and δ13C.