Müller Ice Cap sits on Umingmat Nunaat (Axel Heiberg Island), Nunavut, Canada, ~ 80°N. Its high latitude and elevation suggest it experiences relatively little melt and preserves an undisturbed paleoclimate record. Here, we present a suite of field measurements, complemented by remote sensing, that constrain the ice thickness, accumulation rate, temperature, ice-flow velocity, and surface-elevation change of Müller Ice Cap. These measurements show that some areas near the top of the ice cap are more than 600 m thick, have nearly stable surface elevation, and flow slowly, making them good candidates for an ice core. The current mean annual surface temperature is −19.6 °C, which combined with modeling of the temperature profile indicates that the ice is frozen to the bed. Modeling of the depth-age scale indicates that Pleistocene ice is likely to exist with measurable resolution (300–1000 yr m−1) 20–90 m from the bed, assuming that Müller Ice Cap survived the Holocene Climatic Optimum with substantial ice thickness (~400 m or more). These conditions suggest that an undisturbed Holocene climate record could likely be recovered from Müller Ice Cap. We suggest 91.795°W, 79.874°N as the most promising drill site.

Radiocarbon and uranium-thorium dating of microbialites and penecontemporaneous cements in a microbialite mound at Death Point at Lakeside, Utah, on the shore of Great Salt Lake, Utah, call for a revision of the Lake Bonneville hydrograph. At about 30,000 cal yr BP, the lake experienced an abrupt rise of about 20 m, then dropped back down to levels near or slightly higher than the modern average elevation of Great Salt Lake. Over the ensuing ~6000 yr the lake experienced a series of fluctuations, up to levels a few tens of meters higher than the modern average Great Salt Lake, then down again. The exact timing and amplitudes of those fluctuations are not known, but importantly, the lake did not rise to levels near the Stansbury shoreline (~80 m higher than Great Salt Lake) until after about 24,000 cal yr BP. After the Stansbury shoreline, the lake rose almost 200 m to its highest level at the Bonneville shoreline by about 17,500 cal yr BP. This interpretation is different from previously published hydrographs, many of which show a relatively steady rise to near the Stansbury shoreline between 30,000 and 25,000 cal yr BP.

We review successes and challenges from five recent subglacial bedrock drilling campaigns intended to find evidence for Antarctic Ice Sheet retreat during warm periods in the geologic past. Insights into times when the polar ice sheets were smaller than present serve as guiding information for modeling efforts that aim to predict the rate and magnitude of future sea level rise that would accompany major retreat of the Antarctic Ice Sheet. One method to provide direct evidence for the timing of deglaciations and minimum extent of prior ice sheets is to extract subglacial bedrock cores for cosmogenic nuclide analysis from beneath the modern ice sheet surface. Here we summarize the lessons learned from five field seasons tasked with obtaining bedrock cores from shallow depths (<120 m beneath ice surface) across West Antarctica since 2016. We focus our findings on drilling efforts and technology and geophysical surveys with ground-penetrating radar. Shallow subglacial drilling provides a high risk, high reward means to test for past instabilities of the Antarctic Ice Sheet, and we highlight key challenges and solutions to increase the likelihood of success for future subglacial drilling efforts in polar regions.

The Late Westphalian to Artinskian Haushi Group in the Sultanate of Oman consists of the glaciogenic Al Khlata Formation and the Gharif Formation which contains marginal marine, coastal plain, and fluvial sediments. The sequence was deposited during a global-warming event following the Permo-Carboniferous glaciation of Gondwana. Because of a varied subsidence history, these sediments range from the surface in the SE to almost 5000 m in the NW of the basin.

Mixed-layer illite-smectite (I-S) is an important constituent of the <2 µm size fraction of sandstone and shale samples in both formations at all depths. Different starting compositions lead to three distinct trends of illite layers in I-S versus temperature for different sedimentary environments and paleoclimatic conditions. The starting compositions of I-S at the surface range from an ordered I-S in the Al Khlata Formation to smectite-rich in the Upper+Middle Gharif members.

Physical, chemical and environmental factors were investigated as causes for the different starting compositions of I-S. Both formations share an identical burial history, paragenesis, thermal evolution, and source of detrital material. They differ only in environmental conditions during sedimentation. Thus, the variation in starting composition of I-S appears to be best explained by distinct weathering conditions during sedimentation of the three units. In particular, the expected low intensity of chemical weathering during glaciogenic conditions is marked by the presence of higher amounts of unstable volcanic and sedimentary rock fragments in the Al Khlata Formation.

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.

In this letter we make the case that closer integration of sediment core and passive optical remote sensing data would provide new insights into past and contemporary glacio-sedimentary processes. Sediment cores are frequently used to study past glacial processes and environments as they contain a lengthy geochemical and sedimentological record of changing conditions. In contrast, optical remote sensing imagery is used extensively to examine contemporary glacial processes, including meltwater dynamics, glacial retreat, calving, and ice accumulation. While paleoenvironmental data from sediment cores and optical remote sensing imagery are rarely used in tandem, they are complementary. Sediment core records are spatially discrete, providing long-term paleoenvironmental proxy data which require assumptions about environment-sediment linkages. Optical imagery offers precise, spatially extensive data to visualize contemporary processes often limited in their temporal extent. We suggest that methodologies which integrate optical remotely sensing with sediment core data allow direct observation of processes interpolated from sedimentological analysis and achieve a more holistic perspective on glacial processes. This integration addresses the limitations of both data sources and can achieve a stronger understanding of glacier dynamics by expanding the spatiotemporal extent of data, reducing the uncertainty of interpretations, and broadening the local analyses to regional and global scales.

The clay mineralogy of the Zhada sediments was investigated, using X-ray diffraction and scanning electron microscopy, to obtain a better understanding of climatic change and uplift of the Himalayas in the Zhada region of Tibet. The sediments of Zhada basin in the late Miocene to Pliocene consist of lacustrine and fluvial deposits >800 m thick and can be subdivided into five clay assemblage zones based on their clay-mineral composition. The upward zonation is as follows: (1) smectite-kaolinite; (2) illite-chlorite; (3) chlorite-illite-kaolinite; (4) illite-chlorite; and (5) smectite, illite, and kaolinite. The ratio of chlorite + illite to kaolinite + smectite (Ch+I/K+S) and the Kübler index indicate a warm and humid climate from 9.5 to 8.4 Ma, a cold and dry climate from 8.4 to 7.2 Ma, a warm and seasonal arid climate from 7.2 to 4.5 Ma, a cool and humid climate from 4.5 to 3.6 Ma, and a warm and seasonally humid climate from 3.6 to 3.0 Ma. Intense fluctuations in the Kübler index and in the quantities of evaporite minerals dolomite, aragonite, and gypsum, during the period 7.2–4.5 Ma suggest strong climatic fluctuations between humid and seasonally humid conditions in the Zhada basin. Rapid uplift around the Zhada basin occurred at 8.4 and 3.6 Ma, with sharp subsidence at 7.2 and 4.5 Ma. Evolution of the climate at Zhada showed a different model from that of global climate change, and tectonics-led climate change was the major contributor to climate evolution in the area.

The sedimentary series from Academician Ridge, Lake Baikal, eastern Siberia, was examined using cation exchange capacity (CEC) to estimate the amount of expandable clay minerals (ECM) and high-temperature X-ray diffraction (HT-XRD) to determine their basic classification. The comparison of the magnetic susceptibility (MS) at sub-millennial resolution and the δl8O record of a reference Atlantic core (ODP 980) was used to create an age model. The most closely studied part of the series covered the major part of the last glacial cycle (120–20 ky BP). The HT-XRD analysis is based on monitoring the course of ECM dehydration with 5°C steps between 25 and 250°C and enabled us to improve the discrimination between ECM, chlorite and micas. The CEC obtained at millennial resolution showed that the neoformation of ECM in warmer periods of the last interglacial was either insignificant or fully compensated by their dissolution or dilution. The CEC record was correlated with the main climatic stages in the period studied. Both MS and CEC records reflected the environmental changes at about millennial resolution, including climatic instabilities between 117 and 73 ky BP (late MIS5).

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.

Our limited knowledge of the climate prevailing over Europe during former glaciations is the main obstacle to reconstruct the past evolution of the ice coverage over the Alps by numerical modelling. To address this challenge, we perform a two-step modelling approach: First, a regional climate model is used to downscale the time slice simulations of a global earth system model in high resolution, leading to climate snapshots during the Last Glacial Maximum (LGM) and the Marine Isotope Stage 4 (MIS4). Second, we combine these snapshots and a climate signal proxy to build a transient climate over the last glacial period and force the Parallel Ice Sheet Model to simulate the dynamical evolution of glaciers in the Alps. The results show that the extent of modelled glaciers during the LGM agrees with several independent key geological imprints, including moraine-based maximal reconstructed glacial extents, known ice transfluences and trajectories of erratic boulders of known origin and deposition. Our results highlight the benefit of multiphysical coupled climate and glacier transient modelling over simpler approaches to help reconstruct paleo glacier fluctuations in agreement with traces they have left on the landscape.

Simulating abundances of stable water isotopologues, that is, molecules differing in their isotopic composition, within climate models allows for comparisons with proxy data and, thus, for testing hypotheses about past climate and validating climate models under varying climatic conditions. However, many models are run without explicitly simulating water isotopologues. We investigate the possibility of replacing the explicit physics-based simulation of oxygen isotopic composition in precipitation using machine learning methods. These methods estimate isotopic composition at each time step for given fields of surface temperature and precipitation amount. We implement convolutional neural networks (CNNs) based on the successful UNet architecture and test whether a spherical network architecture outperforms the naive approach of treating Earth’s latitude-longitude grid as a flat image. Conducting a case study on a last millennium run with the iHadCM3 climate model, we find that roughly 40% of the temporal variance in the isotopic composition is explained by the emulations on interannual and monthly timescale, with spatially varying emulation quality. The tested CNNs outperform simple baseline models such as random forest and pixel-wise linear regression substantially. A modified version of the standard UNet architecture for flat images yields results that are as good as the predictions by the spherical CNN. Variations in the implementation of isotopes between climate models likely contribute to an observed deterioration of emulation results when testing on data obtained from different climate models than the one used for training. Future work toward stable water-isotope emulation might focus on achieving robust climate–oxygen isotope relationships or exploring the set of possible predictor variables.

This study utilizes speleothem trace elements as climate proxies to reconstruct hydroclimate variability over approximately 350 years in the Southern Cook Islands. Stalagmites Pu17 and Pu4 from Pouatea cave were analyzed using high-resolution LA-ICP-MS for trace elements (Mg, Na, Sr, P, U, Y). By monitoring cave dripwater and conducting regression analysis, we found that Mg, Sr, and Na in Pouatea dripwater mostly originated from marine aerosols, while Sr and Ba were primarily from bedrock, with additional Ba coming from marine aerosols and weathered oceanic basalt leaching. Mg was identified as the most reliable element for hydroclimate reconstruction due to its predominantly marine aerosol origin. Infiltration, via dilution of marine aerosols and bedrock inputs, was identified as the main driver of trace element variations in Pouatea at a seasonal scale. Transfer functions were established between each trace element and effective infiltration was calculated, with Mg showing the strongest correlation. The reconstructed infiltration data were compared with climate indices, showing an overarching role of the SPCZ and ENSO in controlling rainfall in the South Pacific. This research demonstrates the potential of speleothem trace elements for paleohydroclimate reconstructions, improving understanding of rainfall variability in the climatically vulnerable South Pacific Islands over the past millennia.

Arid regions are especially vulnerable to climate change and land use. More than one-third of Earth's population relies on these ecosystems. Modern observations lack the temporal depth to determine vegetation responses to climate and human activity, but paleoecological and archaeological records can be used to investigate these relationships. Decreasing rainfall across the Late Holocene provides a case study for vegetation response to changing hydroclimate. Rock hyrax (Procavia capensis) middens preserve paleoenvironmental indicators in arid environments where traditional archives are unavailable. Pollen from modern middens collected in Dhofar, Oman, demonstrates the reliability of this archive. Pollen, stable isotope (δ13C, δ15N), and microcharcoal data from fossil middens reveal changes in vegetation, relative moisture, and fire from 4000 cal yr BP to the present. Trees limited to moister areas (e.g., Terminalia) today existed farther inland at ~3100 cal yr BP. After ~2900 cal yr BP, taxa with more xeric affiliations (e.g., Senegalia) had increased. Coprophilous fungal spores (Sporormiella) and grazing indicator pollen revealed an amplified signal of domesticate grazing at ~1000 cal yr BP. This indicates that trees associated with semiarid environments were maintained in the interior desert during ~3000–4000 yr of decreasing rainfall and that impacts of human activity intensified after the transition to a drier environment.

A high-resolution multiproxy lake sediment dataset, comprising lithology, radiography, μXRF elemental, magnetic susceptibility (MS), δ13C, and δ18O measurements since ca. AD 400 is presented in this study. Changes in lithology, radiography, magnetic susceptibility (MS), δ13C, and δ18O reflect wet/dry climate periods, whereas variability in log(Ca/K) can reflect warm/cold climate periods. Analyses of the multiproxy results allow the distinction of several climate periods, which may be associated with climatic phenomena such as changes in North Atlantic Oscillation (NAO) and/or solar activity. The influence of NAO−/NAO+ (negative/positive) is suggested to be related with the southward/northward displacement of the storm tracks resulting from the NAO−/NAO+ phases. For solar activity, the influence is explained through a direct increase in solar heating leading to calcite precipitation. The Dark Ages Cold Period (DACP, AD 450–750) reflects cold-dry climate conditions at this site, indicative of a positive North Atlantic Oscillation (NAO+) and low solar activity. The Medieval Climate Anomaly (MCA, AD 950–1250) exhibits wet-dry-wet and warm-cold-warm climate conditions. The wet/dry periods likely are associated with NAO−/NAO+, respectively, and the warm/cold period may reflect relatively high/low solar activity. The Little Ice Age (LIA, AD 1400–1850) is characterized by dry and cold climate conditions, suggesting the influence of NAO+ and low solar activity. Comparison of the results of this study with local and regional results suggests a generally similar climate pattern, which is indicative of similar climate mechanisms. The contradictions can be associated with age-related uncertainties, orographic differences, and/or other regional teleconnections.

The Buffalo National River in northwest Arkansas preserves an extensive Quaternary record of fluvial bedrock incision and aggradation across lithologies of variable resistance. In this work, we apply optically stimulated luminescence (OSL) dating to strath and fill terraces along the Buffalo River to elucidate the role of lithology and climate on the development of the two youngest terrace units (Qtm and Qty). Our OSL ages suggest a minimum strath planation age of ca. 250 ka for the Qtm terraces followed by a ca. 200 ka record of aggradation. Qtm incision likely occurred near the last glacial maximum (LGM), prior to the onset of Qty fill terrace aggradation ca. 14 ka. Our terrace ages are broadly consistent with other regional terrace records, and comparison with available paleoclimatic archives suggests that terrace aggradation and incision occurred during drier and wetter hydrological conditions, respectively. Vertical bedrock incision rates were also calculated using OSL-derived estimates of Qtm strath planation and displayed statistically significant spatial variability with bedrock lithology, ranging from ~35 mm/ka in the higher resistance reaches and ~16 mm/ka in the lower resistance reaches. In combination with observations of valley width and terrace distribution, these results suggest that vertical processes outpace lateral ones in lithologic reaches with higher resistance.

The Sichuan Basin was a part of the Yangtze Carbonate Platform (YCP) during the Cambrian–Ordovician, and marine carbonates were deposited in the basin during this interval. Although previous studies have evaluated the paleogeography, paleoclimate and paleoecology of this basin, they have primarily focused on the paleoecology and biological evolution in the basin; however, analysis of paleogeography and paleoclimate is lacking. This study integrated outcrop sedimentological and magnetic fabric data to document sedimentary differentiation and anisotropy of magnetic susceptibility (AMS) within the YCP. The aims of this study were to infer paleowind directions during each epoch of the Cambrian–Ordovician and to constrain the paleogeographic location of the YCP. The northwestern, central and southeastern sides of the YCP were characterized by high-energy deposition (e.g. sub-angular to rounded intraclasts), medium-energy deposition (e.g. sub-angular to sub-rounded intraclasts) and low-energy deposition (e.g. angular to sub-angular intraclasts), respectively. The centroid D-Kmax values for the Early, Middle and Late Cambrian were 116° ± 52°, 145° ± 57° and 159° ± 62° from the present north, respectively; corresponding values for the Early, Middle and Late Ordovician were 169° ± 70°, 139° ± 73° and 91° ± 68° from the present north, respectively. Sedimentary differentiation and AMS results indicated that the prevailing wind directions during the Early Cambrian, Middle Cambrian, Late Cambrian, Early Ordovician, Middle Ordovician and Late Ordovician were 296° ± 52°, 325° ± 57°, 339° ± 62°, 349° ± 70°, 319° ± 73° and 271° ± 68° from the present north, respectively. The present study provides evidence for the location of the YCP during the Cambrian–Ordovician via the correspondence between the paleowind directions over the YCP and the trade winds in the Northern and Southern hemispheres. The novelty of this study lies in the following aspects: (1) it integrates microfacies and AMS analyses to establish paleowind patterns; (2) it constrains the paleo-hemispheric location of the YCP during the Cambrian–Ordovician; and (3) it provides a reference for further studies of the paleoclimate and paleogeography of the YCP during the Cambrian–Ordovician.

We present stratigraphic, magnetic, diffuse reflectance spectrophotometric analyses, and chronologic data for the Chenarli loess–paleosol sequence in northeastern Iran. Eight loess units (LU) are identified, described, and constrained in time based on relative stratigraphy and luminescence dating from >130 ± 9.1 ka to ~9.8 ± 0.7 ka. Our data indicate that pedogenic magnetite/maghemite formation gives rise to magnetic enhancement in modern soil and paleosols. The top of LU 7 is demarcated by the well-developed last interglacial soil in which magnetic depletion occurred. We infer that magnetic depletion in this paleosol was produced by reducing conditions in a seasonally waterlogged soil during a warm and wet phase within Marine Isotope Stage (MIS) 5e. Units LU 6 to 1 record several episodes of dust accumulation and soil formation during the last glacial and Holocene. Increased dust accumulation rates occurred during middle-late MIS 3 and lasted into the late MIS 2, with a peak during the last glacial maximum (LU 2). These findings could be applicable to understanding magnetic enhancement/dissolution mechanism in the loess–paleosol sequences in study area. We infer paleoenvironmental changes in northeastern Iran relative to northern Iran, Eurasia, and China.

Cosmogenic nuclide measurements in glacial deposits extend our knowledge of glacier chronologies beyond the observational record. The short half-life of in situ cosmogenic 14C makes it particularly useful for studying glacier chronologies, as resulting exposure ages are less sensitive to nuclide inheritance when compared with more commonly measured, long-lived nuclides. An increasing number of laboratories using an automated process to extract carbon from quartz has led to in situ 14C measurements in Antarctic samples at an accelerating rate over the past decade, shedding light on deglaciation in Antarctica. In situ 14C has had the greatest impact in the Weddell Sea Embayment, where inferences on the thickness of ice and timing of deglaciation were limited by inheritance in other cosmogenic nuclide systems. Future subglacial measurements of the nuclide hold much potential as they can provide direct evidence of proposed Holocene thinning and subsequent re-thickening of parts of the Antarctic ice sheets.