Search

Ice-coring auger machines can cost thousands of dollars and are generally made for vertical extraction from sea ice or glaciers. This project required specimens from ice accreted onto vertical piles at the Port of Alaska (POA) in Anchorage, Alaska, which meant the extraction was horizontal and had to be completed during the short low-tide window. In this paper, we describe the development of a method for economically harvesting and processing the accreted ice from the POA into cylindrical cores and shear beams. The ice was harvested by cutting large cubes out of the in situ mass with a chainsaw and then lifting these onto the port deck using a crane. The cubes were then transported to a laboratory walk-in freezer to be processed. Ice cores were drilled from the cubes using a custom-fabricated core bit made from metal tubing and hole saws. This simple core bit quickly produced smooth ice surfaces, provided there was a mechanism for evacuating the cut ice chips during drilling. To facilitate this, an alignment template was used to locate and drill two adjoining evacuation holes prior to coring. These evacuation holes had to be precisely located to be adjacent, parallel, and connected with the outside diameter of the forthcoming core, but were not allowed to impinge on the core. With this new method, an ice block can be converted into up to a dozen test-ready compression ice cores within a few hours. The variability in the diameter of the ice cores was found to be small, with a coefficient of variation of < 0.4%. [ABSTRACT FROM AUTHOR]

Major explosive volcanic eruptions were important triggers of abrupt climate changes during the Holocene and crucial sources of Hg to the atmosphere, yet there remains limited understanding regarding the long‐range transportation of this volcanic Hg and its imprint in natural archives. Here, we present a reconstruction of Holocene global volcanism based on the anomalies in Hg concentrations, accumulation fluxes, and Hg/C ratios in three high‐resolution peat profiles spanning Eurasia. Our reconstruction reveals that the two Tibetan peat profiles recorded 33 major explosive volcanic eruptions (with 11 eruptions being synchronously detected), which correspond with a French Pyrenees peat record and sulfate anomalies in polar ice cores. Additionally, the major explosive volcanic eruptions recorded in the TP peat profiles coincided with abrupt decreases in solar irradiance during the Holocene, suggesting these eruptions might have had a greater global climate impact. Our results suggest the atmospheric transport of volcanic Hg within the Northern Hemisphere and underscore the significant role played by major explosive volcanic eruptions in precipitating abrupt global climate and environmental changes during the Holocene. This study has implications for deciphering the configuration of volcanic eruption seasons, locations, and magnitudes during the Holocene and aligning the chronology of peat deposits with ice cores. Plain Language Summary: Major explosive volcanic eruptions, one of the principal geogenic sources of Hg to the atmosphere, are widely acknowledged as significant drivers of abrupt global environmental and climate changes during the Holocene. While the magnitudes and timing of these eruptions have been well documented in ice core records, the long‐range transportation of the volcanic materials (such as Hg and SO2) and their terrestrial dispersion remain largely unknown. In this study, we present a reconstruction of Holocene global volcanism based on the anomalies in Hg concentrations, accumulation fluxes, and Hg/C ratios in three high‐resolution peat profiles, with two collected from the Tibetan Plateau and one from the French Pyrenees. Our reconstruction reveals that the Tibetan peat profiles recorded 33 major volcanic eruptions, which correspond with a French peat record and the sulfate anomalies in polar ice cores. Additionally, we demonstrate that certain major explosive volcanic eruptions likely triggered several abrupt decreases in solar irradiance during the Holocene. Our findings suggest the convective atmospheric transport of volcanic Hg within the Northern Hemisphere extratropical region and underscore the significant role played by major explosive volcanic eruptions in precipitating abrupt global climate and environmental changes during the Holocene. Key Points: Eleven major explosive volcanic eruptions during the Holocene were synchronously recorded in two peat profiles in the Tibetan PlateauThe recorded major explosive volcanic eruptions coincide with nine abrupt decreases in solar irradiance during the HoloceneThe Hg anomalies in peat profiles have the potential to be chronology markers that were linked with volcanic eruptions [ABSTRACT FROM AUTHOR]

Establishing accurate chronological frameworks is imperative for reliably identifying lead–lag dynamics within the climate system and enabling meaningful intercomparisons across diverse paleoclimate proxy records over long time periods. Robust age models provide a solid temporal foundation for establishing correlations between paleoclimate records. One of the primary challenges in constructing reliable radiocarbon-based chronologies in the marine environment is to determine the regional marine radiocarbon reservoir age correction. Calculations of the local marine reservoir effect (ΔR) can be acquired using independent 14C dating methods, such as synchronization with other well-dated archives. The cosmogenic radionuclide 10Be offers such a synchronization tool. Its atmospheric production rate is controlled by the global changes in the cosmic ray influx, caused by variations in solar activity and geomagnetic field strength. The resulting fluctuations in the meteoric deposition of 10Be are preserved in sediments and ice cores and can thus be utilized for their synchronization. In this study, for the first time, we use the authigenic 10Be/9Be record of a Laptev Sea sediment core for the period 8–14 kyr BP and synchronize it with the 10Be records from absolutely dated ice cores. Based on the resulting absolute chronology, a benthic ΔR value of + 345 ± 60 14Cyears was estimated for the Laptev Sea, which corresponds to a marine reservoir age of 848 ± 90 14Cyears. The ΔR value was used to refine the age–depth model for core PS2458-4, establishing it as a potential reference chronology for the Laptev Sea. We also compare the calculated ΔR value with modern estimates from the literature and discuss its implications for the age–depth model. [ABSTRACT FROM AUTHOR]

Stable water isotopes recorded in Antarctic ice cores have traditionally been used to infer past surface air temperatures (SATs). During the historical period (1850 onward), observational data and good-quality ice core records overlap, yielding an opportunity to investigate key relationships between ice core stable water isotope (δ18 O) measurements and the Antarctic climate. We present a new ensemble of climate model simulations covering 1851–2004 using the UK Met Office HadCM3 general circulation model equipped with stable water isotopes. Our ensemble captures observed historical SAT and precipitation trends and weak δ18 O trends. The weak δ18 O trends mean there is no significant relationship between SAT and δ18 O over one-third of Antarctica, and also half of our considered ice core sites, though relationships are stronger when using regional averages. The strongest regional relationships occur in the West Antarctic Ice Sheet (WAIS) region. This decoupling between SAT and δ18 O occurs primarily because of the impact of autumnal sea ice loss during the simulated warming. The warming and sea ice loss are associated with (i) changes in near-coastal air mass intrusions (synoptic effects) induced by changes in the large-scale circulation and/or sea ice; (ii) direct sea-ice-driven changes in moisture pathways (especially lengths) to Antarctica; and (iii) precipitation seasonality changes, again mostly driven by sea ice changes. Consequently, when reconstructing temperatures over these timescales, changes in sea ice need to be considered, both to determine the most appropriate SAT and δ18 O relationship and to understand how uncertainties affect the inference of past temperature from ice core δ18 O measurements. [ABSTRACT FROM AUTHOR]

During the Mid-Pleistocene Transition (MPT), ca. 1200–800 000 years ago (ka), the Earth's glacial cycles changed from 41 to 100 kyr periodicity. The emergence of this longer ice age periodicity was accompanied by higher global ice volume in glacial periods and lower global ice volume in interglacial periods. Since there is no known change in external orbital forcing across the MPT, it is generally agreed that the cause of this transition is internal to the Earth system. Resolving the climate, carbon cycle, and cryosphere processes responsible for the MPT remains a major challenge in Earth and palaeoclimate science. To address this challenge, the international ice core community has prioritised recovery of an ice core record spanning the MPT interval. Here we present results from a simple generalised least-squares (GLS) model that predicts atmospheric CO2 out to 1.8 Myr. Our prediction utilises existing records of atmospheric carbon dioxide (CO2) from Antarctic ice cores spanning the past 800 kyr along with the existing LR04 benthic δ18Ocalcite stack (Lisiecki and Raymo, 2005; hereafter "benthic δ18O stack") from marine sediment cores. Our predictions assume that the relationship between CO2 and benthic δ18O over the past 800 000 years can be extended over the last 1.5 million years. The implicit null hypothesis is that there has been no fundamental change in feedbacks between atmospheric CO2 and the climate parameters represented by benthic δ18O , global ice volume, and ocean temperature. We test the GLS-model-predicted CO2 concentrations against observed blue ice CO2 concentrations, δ11B -based CO2 reconstructions from marine sediment cores, and δ13C of leaf-wax-based CO2 reconstructions (Higgins et al., 2015; Yan et al., 2019; Yamamoto et al., 2022). We show that there is no clear evidence from the existing blue ice or proxy CO2 data to reject our predictions or our associated null hypothesis. A definitive test and/or rejection of the null hypothesis may be provided following recovery and analysis of continuous oldest ice core records from Antarctica, which are still several years away. The record presented here should provide a useful comparison for the oldest ice core records and an opportunity for further constraints on the processes involved in the MPT. [ABSTRACT FROM AUTHOR]

Dansgaard–Oeschger (DO) events occurred throughout the last glacial period. Greenland ice cores show a rapid warming during each stadial to interstadial transition, alongside an abrupt loss of sea ice and a major reorganisation of the atmospheric circulation. Other records also indicate simultaneous abrupt changes to the oceanic circulation. Recently, an advanced Bayesian ramp-fitting method has been developed and used to investigate time lags between transitions in these different climate elements with a view to determining the relative order of these changes. Here, we critically review this method in both its original implementation and a new, extended implementation. Using ice core data, climate model output, and carefully synthesised data representing DO events, we demonstrate that both implementations of the method suffer from biases of up to 15 years. These biases mean that the method will tend to yield transition onsets that are too early. Further investigation of DO warming event records in climate models and ice core data reveals that the biases are on the same order of magnitude as potential timing differences between the abrupt transitions of different climate elements. Additionally, we find that higher-resolution records would not reduce these biases. We conclude that decadal-scale leads and lags between climate elements across DO events cannot be reliably detected, as we cannot exclude the possibility that they result solely from the biases we present here. [ABSTRACT FROM AUTHOR]

Ice core oxygen isotope (δ18O) records from low-latitude regions preserve high-resolution climate records in the past, yet the interpretation of these ice core δ18O records is still facing difficulty due to the uncertainty of ice core dating. Here we present a new established δ18O time series from Qiangtang (QT) No. 1 ice core retrieved from the central Tibetan Plateau. Due to the vague seasonal signals in the QT ice core, we investigated the spectral properties of δ18O record with depth and discussed the implications of significant spectral power peaks in the QT ice core. We employed a variational mode decomposition (VMD) analysis for the upper part of the QT ice core to decompose the δ18O depth series in order to separate the El Niño Southern Oscillation (ENSO) mode, a signal strongly preserved in the QT ice core δ18O record. With this approach, we established a time series of 335 years (1677–2011 CE) for the upper 50 m of the QT ice core. Subsequently, we examined the frequency of the new established δ18O time series and detected strong signals of the bidecadal and multidecadal modes of Pacific Decadal Oscillation (PDO). The PDO consists of two modes with periods of approximately 25–35 years and 50–70 years, and we found that the 50–70 years periodicity has persisted since 1700 CE, succeeded by dominance of the 25–75 years periodicity after 1900 CE. Additionally, we analyzed the δ18O series of the QT ice core during the past century and determined that the increasing frequency of El Niño events is an important factor contributing to the increase in recent ice core δ18O. [ABSTRACT FROM AUTHOR]

Pleistocene Ice Ages display abrupt Dansgaard-Oeschger (DO) climate oscillations that provide prime examples of Earth System tipping points--abrupt transition that may result in irreversible change. Greenland ice cores provide key records of DO climate variability, but gas-calibrated estimates of the temperature change magnitudes have been limited to central and northwest Greenland. Here, we present ice-core δ15N-N2 records from south (Dye 3) and coastal east Greenland (Renland) to calibrate the local water isotope thermometer and provide a Greenland-wide spatial characterization of DO event magnitude. We combine these data with existing records of δ18O, deuterium excess, and accumulation rates to create a multiproxy "fingerprint" of the DO impact on Greenland. Isotope-enabled climate models have skill in simulating the observational multiproxy DO event impact, and we use a series of idealized simulations with such models to identify regions of the North Atlantic that are critical in explaining DO variability. Our experiments imply that wintertime sea ice variation in the subpolar gyre, rather than the commonly invoked Nordic Seas, is both a sufficient and a necessary condition to explain the observed DO impacts in Greenland, whatever the distal cause. Moisture-tagging experiments support the idea that Greenland DO isotope signals may be explained almost entirely via changes in the vapor source distribution and that site temperature is not a main control on δ18O during DO transitions, contrary to the traditional interpretation. Our results provide a comprehensive, multiproxy, data-model synthesis of abrupt DO climate variability in Greenland. [ABSTRACT FROM AUTHOR]

The Mount Brown South ice core (MBS 69.111° S 86.312° E) is a new, high resolution ice core drilled in coastal East Antarctica. With mean annual accumulation estimated to be 20–30 cm ice equivalent accumulation throughout the length of the core (∼295 m), MBS represents a high resolution archive of ice core data spanning 1137 years (873–2009 CE), from an area previously underrepresented by high resolution ice core data. Here, we present a high-resolution dataset of chemistry and impurities obtained via continuous flow analysis (CFA). The dataset consists of meltwater electrolytic conductivity, sodium (Na+), ammonium (NH4+), hydrogen peroxide (H2O2), and insoluble microparticle measurements. The data are presented in three datasets: as a 1 mm depth resolution record, 3 cm averaged record, and decadal average record. The 1 mm record represents an oversampling of the true resolution, as due to smoothing effects the actual resolution is closer to 3 cm for some species. Therefore, the 3 cm resolution dataset is considered to be the minimum true resolution given the system setup. We also describe the current Copenhagen CFA system, and provide a detailed assessment of data quality, precision, and functional resolution. The 1 mm averaged, 3 cm averaged, and MBS2023 decadal averaged datasets are available at the Australian Antarctic Data Center: http://dx.doi.org/doi:10.26179/9tke-0s16 (Harlan et al., 2024). [ABSTRACT FROM AUTHOR]

Antarctic sea ice, a key component in the complex Antarctic climate system, is an important driver and indicator of the global climate. In the relatively short satellite‐observed period from 1979 to 2022 the sea ice extent has continuously increased (contrasting a major decrease in Arctic sea ice) up to a dramatic decrease between 2014 and 2017. Recent years have seen record sea ice lows in February 2022–February 2023. We use a statistical ensemble reconstruction of Antarctic sea ice to put the observed changes into the historical context of the entire 20th century. We propose a seasonal Vector Auto‐Regressive Moving Average (VARMA) model fit in a Bayesian framework using regularized horseshoe priors on the regression coefficients to create a stochastic ensemble reconstruction of monthly Antarctic Sea ice extent from 1900 to 1979. This novel model produces a set of 2,500 plausible sea ice extent reconstructions for the sea ice by sector that incorporate the autocorrelation structure of sea ice over time as well as the dependence of sea ice between the sectors. These fully observed reconstructions exhibit plausible month‐to‐month changes in reconstructed sea ice as well as plausible interactions between the sectors and the total. We reconstruct an overall higher sea ice extent earlier in the 20th century with a relatively sharp decline in the 1970s. These trends agree well with previous reconstructions of Antarctic sea ice based on ice core data, whaling locations, and climatological data, as well as early satellite observations in the reconstruction period. Plain Language Summary: Antarctic sea ice, a key component in the complex Antarctic climate system, is an important driver and indicator of the global climate. In the relatively short satellite‐observed period from 1979 to 2022 the sea ice extent has continuously increased (contrasting a major decrease in Arctic sea ice) up to a dramatic decrease between 2014 and 2017. Recent years have seen record sea ice lows in February 2022 and February 2023. We reconstruct the Antarctic sea ice extent to put the observed changes into the historical context of the entire 20th century. We propose a novel statistical model to create multiple reconstructions of monthly Antarctic Sea ice extent from 1900 to 1979. This state‐of‐the‐art model produces a set of 2,500 plausible sea ice extent reconstructions that exhibit plausible month‐to‐month changes in reconstructed sea ice as well as plausible interactions between the sectors and the total. We reconstruct an overall higher sea ice extent earlier in the 20th century with a relatively sharp decline in the 1970s. These trends agree well with previous reconstructions of Antarctic sea ice based on ice core data, whaling locations, and climatological data, as well as early satellite observations in the reconstruction period. Key Points: Introduction of statistical model for ensemble reconstruction of Antarctic sea ice extentReconstruction of unobserved 20th century Antarctic sea ice extent and imputation of missing values in satellite observed periodBayesian framework creates plausible reconstructions of Antarctic sea ice by sector displaying realistic spatio‐temporal autocorrelation [ABSTRACT FROM AUTHOR]

Ice formation is generally considered to exclude many particles and most solutes and thus be relatively pure compared to ambient waters. Because river ice forms by a combination of thermal and mechanical processes, some level of sediment entrainment in the ice column is likely, though reports of sediment in river ice are limited. We observed high and sporadic levels of silt and sand in ice of the Kuskokwim and Tanana rivers (Alaska, the United States) during routine field studies. These observations led us to make a more comprehensive survey of sediment entrainment in river ice of the Kuskokwim and Yukon rivers and several of their tributaries. We collected and subsampled 48 ice cores from 19 different river locations in March 2023, which included concurrent measurements of water turbidity, velocity, and depth. Approximately 60% of cores contained detectable levels of sediment, averaging 438 mg/L with median concentrations exceeding 1000 mg/L in three cores from the Yukon and Kuskokwim main stems. Many cores had even higher concentrations at certain intervals, with seven cores having subsamples exceeding 2000 mg/L; these were often located in the middle or lower portion of the ice column. Jumble ice, formed mechanically by frazil‐pan jamming during freeze‐up, was generally the best predictor of higher sediment entrainment, and these locations often had higher under‐ice velocities and depths. Our observation of high and widespread sediment entrainment in northern river ice, particularly in jumble‐ice fields, may have implications for sediment transport regimes, ice strength and transportation safety, and how rivers break up in the springtime. [ABSTRACT FROM AUTHOR]

Studies of the Last Interglacial (ca. 129 to 116 ka BP) provide an opportunity to study the impact of high‐latitude warm temperatures on the Earth system. To build an accurate spatio‐temporal picture of climate and environmental variability during the Last Interglacial, building robust chronologies, through which the patchwork of terrestrial, marine, and ice core archives can be correlated, is paramount. In this review, we briefly evaluate the most common approaches used to date climate and environmental archives from the Last Interglacial, as well as the chronostratigraphic tools available for direct correlation between sequences, with a focus on terrestrial archives. We then present a case study on the use of pollen biostratigraphy for correlating sequences in NW Europe, highlighting its strength as a local correlation tool, and the challenges this approach presents in comparing sequences to global records of climate and environmental change. A move towards consistently dated sequences will improve our understanding of environmental responses to Last Interglacial climate variability across different elements of the global climate system and the rates at which different elements of the climate system respond to changes in global temperatures. [ABSTRACT FROM AUTHOR]

The abrupt warming events punctuating the Termination 1 (about 11.7–18 ka Before Present, BP) were marked by sharp rises in the concentration of atmospheric methane (CH4). The role of permafrost organic carbon (OC) in these rises is still debated, with studies based on top‐down measurements of radiocarbon (14C) content of CH4 trapped in ice cores suggesting minimum contributions from old and strongly 14C‐depleted permafrost OC. However, organic matter from permafrost can exhibit a continuum of 14C ages (contemporaneous to >50 ky). Here, we investigate the large‐scale permafrost remobilization at the Younger Dryas‐Preboreal transition (ca. 11.6 ka BP) using the sedimentary record deposited at the Lena River paleo‐outlet (Arctic Ocean) to reflect permafrost destabilization in this vast drainage basin. Terrestrial OC was isolated from sediments and characterized geochemically measuring δ13C, Δ14C, and lignin phenol molecular fossils. Results indicate massive remobilization of relatively young (about 2,600 years) permafrost OC from inland Siberia after abrupt warming triggered severe active layer deepening. Methane emissions from this young fraction of permafrost OC contributed to the deglacial CH4 rise. This study stresses that underestimating permafrost complexities may affect our comprehension of the deglacial permafrost OC‐climate feedback and helps understand how modern permafrost systems may react to rapid warming events, including enhanced CH4 emissions that would amplify anthropogenic climate change. Plain Language Summary: Feedback on climate change through methane emissions from organic carbon (OC) released from thawing permafrost soils is still debated. Studies investigating rapid warming events from the last 15,000 years and measuring the radiocarbon content in methane from ancient air trapped in ice cores suggested small contributions from the oldest (10,000 years or more) and most refractory OC fraction to the deglacial methane rise. However, the younger and more reactive OC likely contributed more to methane emissions from permafrost systems. Using terrestrial material from sediments deposited at the paleo‐delta of the Lena River (Arctic Ocean), we investigated large‐scale permafrost OC remobilization across the climate transition to the Preboreal (about 11,600 years ago). Results show that relatively young (about 2,600 years) terrestrial organic matter was largely remobilized from the layer of East Siberian permafrost soils thawing seasonally that deepened fast due to abrupt warming. Using the young radiocarbon age of remobilized permafrost OC established in this study suggests that permafrost systems might have contributed more than previously thought to the deglacial methane rise. Considering more components of the heterogeneous permafrost OC pool will be crucial to better constrain past and future impacts on climate change. Key Points: A few thousand years old permafrost organic carbon (OC) was massively remobilized from East Siberian soils during the Preboreal climate transitionRemineralization of a few thousand years old OC may increase the maximum allowable deglacial methane emissions from permafrostAge complexities of permafrost OC must be considered to better constrain the permafrost carbon‐climate feedback [ABSTRACT FROM AUTHOR]

Iodine and bromine play central roles in polar atmospheric chemistry: iodine influences the atmospheric oxidative capacity and can generate cloud condensation nuclei, while bromine participates in ozone depletion reactions, known as bromine explosions. Here we present a very sensitive analytical method for Br and I speciation by coupling the ion chromatography system (IC) with an inductively coupled plasma sector field mass spectrometer (ICP-SFMS). We achieved sub-picogram-per gram (pg g−1) as limits of detection (LODs) ranging from 0.4 pg g−1 for I−, 0.8 pg g−1 for IO3−, 4 pg g−1 for Br−, and 1 pg g−1 for BrO3−, respectively. These values represent a decrease of up to 30 times compared to the LODs reported in other studies. The method was validated using deep snow samples from the Svalbard Islands, collected at the end of the polar night to quantify various oxidized compounds during their seasonal minimum. In the future, this method could prove useful in the paleoclimatic study of ice cores and snow, as well as in ice chemistry research. [ABSTRACT FROM AUTHOR]

Dust records extracted from ice cores can facilitate the reconstruction of historical atmospheric dust levels and climate change. However, interpreting dust variations in ice cores is intricate because of the compounded influence of emission, transport, and deposition processes. This study investigated dust records retrieved from the Guliya ice cap drilled in 2015 on the West Tibetan Plateau using a mean trajectory transport and deposition model. Results showed that the Guliya dust concentration has exhibited a declining trend since the 1960s (−751 μg kg−1 yr−1). Applying an attribution approach, we discovered that low dust emission (80.3%) was the main cause of the drop in dust concentration, with changes related to transportation (5.2%) and deposition (14.5%) making only minor contributions. The weakening of surface wind speed in the desert and increasing precipitation in both the desert and glacier were the primary factors driving the decrease in Guliya dust concentration. Plain Language Summary: Dust emission, transport, and deposition processes are the key factors influencing the dust record in an ice core. Consequently, relying solely on ice core records may not be sufficient for the accurate reconstruction of past desertification processes and the evolution of the climate and environment. In this study, we endeavored to separate these complex processes from the Guliya ice core records on the Tibetan Plateau by integrating a transport model, thereby quantifying the impact of these processes on ice core records over the past 50 years. We found that the dust concentration changes in the Guliya ice core were primarily attributed to changes in dust emission strength, with a secondary influence of deposition and transport changes. This study can enhance the scientific interpretation of dust records preserved in ice cores, providing novel insights for a deeper comprehension of the climate and environmental evolution of the Tibetan Plateau. Key Points: The dust concentration in the Guliya ice core has decreased over the past 50 yearsThe decreasing trend of Guliya dust is primarily controlled by the emission process in the Taklimakan DesertThe principal drivers behind the reduction in Guliya dust were decreased surface wind speed and increased precipitation [ABSTRACT FROM AUTHOR]

We provide a detailed description of the ice-sheet layer age tracer Englacial Layer Simulation Architecture (ELSA) – a model that uses a straightforward method to simulate the englacial stratification of large ice sheets – as an alternative to Eulerian or Lagrangian tracer schemes. ELSA's vertical axis is time, and individual layers of accumulation are modeled explicitly and are isochronal. ELSA is not a standalone ice-sheet model but requires unidirectional coupling to another model providing ice physics and dynamics (the "host model"). Via ELSA's layer tracing, the host model's output can be evaluated throughout the interior using ice core or radiostratigraphy data. We describe the stability and resolution dependence of this coupled modeling system using simulations of the last glacial cycle of the Greenland ice sheet using one specific host model. Key questions concern ELSA's design to maximize usability, which includes making it computationally efficient enough for ensemble runs, as well as exploring the requirements for offline forcing of ELSA with output from a range of existing ice-sheet models. ELSA is an open source and collaborative project, and this work provides the foundation for a well-documented, flexible, and easily adaptable model code to effectively force ELSA with (any) existing full ice-sheet model via a clear interface. [ABSTRACT FROM AUTHOR]

In the context of glacier retreat and increased precipitations, Arctic glacier basin slopes are subject to stress leading to visible transformations. In this work, subsurface features of a small Arctic glacier basin slopes are mapped using ground‐penetrating RADAR. In combination with surface topography data, eight transects were surveyed ranging from the areas furthest from the current glacier extent to the areas still in contact with the glacier. This allowed for a reconstitution of the successive stages ice‐cored slopes go through when glaciers retreat. It appears that slopes evolve from thick debris‐covered ice bodies connected with the glacier, to residual ice and ice/debris mixes covered in debris. At the same time, surface morphology of the slopes shifts from homogeneous ice‐cored slope gradients to more complex talus‐type slopes at the end of the process. The stages of these evolutions are in compliance with former glacier extents. The main driving factors of the slopes successive stages are the constant slope adjustments linked to debris movements, and the melting of ice cores. All these factors are exacerbated by the warmer and wetter conditions they are subject to. [ABSTRACT FROM AUTHOR]

We report high resolution measurements of the stable water isotope ratios (δ18O, δD) from the Mount Brown South ice core (MBS, 69.11° S 86.31° E). The record covers the period 873 - 2009 CE with sub-annual temporal resolution. Preliminary analyses of surface cores have shown the Mount Brown South site has relatively high annual snowfall accumulation (0.3 metres ice equivalent) with a seasonal bias toward lower snowfall during austral summer. Precipitation at the site is frequently related to intense, short term synoptic scale events from the mid-latitudes of the southern Indian Ocean. Higher snowfall regimes are associated with easterly winds, while lower snowfall regimes are associated with south-easterly winds. Isotope ratios are measured with Infra-Red Cavity Ring Down Spectroscopy, calibrated on the VSMOW/SLAP scale and reported on the MBS2023 time scale interpolated accordingly. We provide estimates for measurement precision and internal accuracy for δ18O and δD. [ABSTRACT FROM AUTHOR]

The trend in partial pressure of atmospheric CO2, P(CO2), across the 66 MYr of the Cenozoic requires elucidation and explanation. The Null Hypothesis sets sea surface temperature (SST) as the baseline driver for Cenozoic P(CO2). The crystallization and cooling of flood basalt magmas is proposed to have heated the ocean, producing the Paleocene–Eocene Thermal Maximum (PETM). Heat of fusion and heat capacity were used to calculate flood basalt magmatic Joule heating of the ocean. Each 1 million km3 of oceanic flood basaltic magma liberates ~5.4 × 1024 J, able to heat the global ocean by ~0.97 °C. Henry's Law for CO2 plus seawater (HS) was calculated using δ18O proxy-estimated Cenozoic SSTs. HS closely parallels Cenozoic SST and predicts the gas solute partition across the sea surface. The fractional change of Henry's Law constants, H n − H i H n − H 0 is proportional to ΔP(CO2)i, and H n − H i H n − H 0 × ∆ P (C O 2) + P (C O 2 ) m i n , where ΔP(CO2) = P(CO2)max − P(CO2)min, closely reconstructs the proxy estimate of Cenozoic P(CO2) and is most consistent with a 35 °C PETM ocean. Disparities are assigned to carbonate drawdown and organic carbon sedimentation. The Null Hypothesis recovers the glacial/interglacial P(CO2) over the VOSTOK 420 ka ice core record, including the rise to the Holocene. The success of the Null Hypothesis implies that P(CO2) has been a molecular spectator of the Cenozoic climate. A generalizing conclusion is that the notion of atmospheric CO2 as the predominant driver of Cenozoic global surface temperature should be set aside. [ABSTRACT FROM AUTHOR]

The Deliblato (Banat) Sand Sea, which is one of the largest areas of аeolian sand in Europe, is located near the Iron Gate, which marks the crossing of the Danube River through the biggest gorge of this river. Here, Danubian alluvium has served as the sand source for the Banat Sand Sea, which was formed primarily through southeasterly (Košava) winds. Utilizing a multi-proxy approach, the objective of this study is to gain a better understanding of the environmental dynamics of the Banat Sand Sea. To achieve this goal, we conducted an analysis of an archive representing an approximately 20-m-thick dune formation on the northern edge of this dune field. Using optically stimulated luminescence (OSL) dating, we calculated aeolian sedimentation rates and dune ages. Sand was deposited here approximately between 17 ka and 13 ka. Magnetic susceptibility, grain size, and colorimetric analyses were interpreted in terms of local paleoenvironmental conditions. Calculated sedimentation rates (SR) indicate intensive aeolian deposition during the study period that range from 483 cm/ka to 502 cm/ka. We compared our data with regional and other European archives, as well as with climatic variations recorded in the Greenland ice core North Greenland Ice Core Project (NGRIP). [ABSTRACT FROM AUTHOR]

We investigate the lock-in zone (LIZ) of the East Greenland Ice Core Project (EastGRIP) region, northeastern Greenland, in detail. We present results from the firn air-pumping campaign of the S6 borehole, forward modeling, and a novel technique for finding the lock-in depth (LID, the top of the LIZ) based on the visual stratigraphy of the EastGRIP ice core. The findings in this work help to deepen our knowledge of how atmospheric gases are trapped in ice cores. CO2 , δ15N , and CH4 data suggest that the LID lies around 58 to 61 m depth. With the pixel value intensity and bright-spot analysis based on visual stratigraphy, we can pinpoint a change in ice properties to exactly 58.3 m depth, which we define as the optical lock-in depth (OLID). This visual change in ice properties is caused by the formation of rounded and enclosed air bubbles that alter the measured refraction of the light pathways. The results for the LID and OLID agree accurately on the depth. We furthermore use the visual stratigraphy images to obtain information on the sharpness of the open- to closed-porosity transition. Combining traditional methods with the independent optical method presented here, we can now better constrain the bubble closure processes in polar firn. [ABSTRACT FROM AUTHOR]

Basal materials in ice cores hold information about paleoclimate conditions, glacial processes, and the timing of past ice-free intervals, all of which aid understanding of ice sheet stability and its contribution to sea level rise in a warming climate. Only a few cores have been drilled through ice sheets into the underlying sediment and bedrock, producing limited material for analysis. The last of three Camp Century ice cores, which the U.S. Army collected in northwestern Greenland from 1963–1966 CE, recovered about 3.5 m of subglacial material, including ice and sediment. Here, we document the scientific history of the Camp Century subglacial material. We present our recent core-cutting, sub-sampling, and processing methodology and results for this unique archive. In 1972 CE, curators at the Buffalo, New York, Ice Core Laboratory cut the original core sections into 32 segments that were each about 10 cm long. Since then, two segments were lost and are unaccounted for, two were thawed, and two were cut as pilot samples in 2019 CE. Except for the two thawed segments, the rest of the extant core has remained frozen since collection. In 2021 CE, we documented, described, and then cut each of the remaining frozen archived segments (n=26). We saved an archival half and cut the working half into eight oriented sub-samples under controlled temperature and light conditions for physical, geochemical, isotopic, sedimentological, magnetic, and biological analyses. Our approach was designed to maximize sample usage for multiproxy analysis, minimize contamination, and preserve archive material for future analyses of this legacy subglacial material. Grain size, bulk density, sedimentary features, magnetic susceptibility, and ice content, as well as pore ice pH and conductivity, suggest that the basal sediment contains five stratigraphic units. We interpret these stratigraphic units as representing different depositional environments in subglacial or ice-free conditions: from bottom to top, a diamicton with subhorizontal ice lenses (Unit 1), vertically fractured ice with dispersed fine-grained sediments (<20 % in mass) (Unit 2), a normally graded bed of pebbles to very fine sand in an icy matrix (Unit 3), bedded very fine to fine sand (Unit 4), and stratified medium to coarse sand (Unit 5). Plant macrofossils are present in all samples and are most abundant in Units 3 and 4; insect remains are present in some samples (Units 1, 3, and 5). Our approach provides a working template for future studies of ice core basal materials because it includes intentional planning of core sub-sampling, processing methodologies, and archiving strategies to optimize the collection of paleoclimate, glacial process, geochemical, geochronological, and sediment properties from archives of limited size. Our work benefited from a carefully curated and preserved archive, allowing the application of analytical techniques not available in 1966 CE. Preserving uncontaminated core material for future analyses that use currently unavailable tools and techniques is an important consideration for rare archive materials such as these from Camp Century. [ABSTRACT FROM AUTHOR]

A 10-year record of oxygen and hydrogen isotopic composition of precipitation is presented here: from 2008 to 2017, 1483 daily precipitation samples were collected year-round on a raised platform at Concordia Station, East Antarctica. Weather data were retrieved from the Italian Antarctic Meteo-Climatological Observatory automatic weather station (AWS), while ERA5 was used to estimate total precipitation. The δ –temperature relationships were moderately high for daily data (r2=0.63 and 0.64 for δ18 O and δ2 H, respectively) and stronger using monthly data (r2=0.82 for both δ18 O and δ2 H), with a slope of about 0.5 ‰ °C -1 for δ18O/TAWS (3.5 ‰ °C -1 for δ2H/TAWS), which remains consistent also using annual averages. The isotopic composition of precipitation is the input signal of the snow–ice system, and this dataset will be useful to improve the interpretation of paleoclimate records and promote a better understanding of the post-depositional processes affecting the isotopic signal in ice cores. This dataset represents a benchmark for the evaluation of isotope-enabled general circulation models. Here, the ECHAM6-wiso output was compared to experimental data, showing moderately good relationships for δ18 O and δ2 H but not for d-excess, nonetheless marking a substantial improvement from the previous release of the model. [ABSTRACT FROM AUTHOR]

The 2050 ± 50 14C yBP caldera-forming eruption of Okmok volcano, Alaska, had a global atmospheric impact with tephra deposits found in distant Arctic ice cores and a sulfate signal found in both Greenland and Antarctic ice cores. The associated global climate cooling was driven by the amount of sulfur injected into the stratosphere during the climactic phase of the eruption. This phase was dominated by pyroclastic density currents, which have complex emplacement dynamics precluding direct estimates of the sulfur stratospheric load. We simulated the dynamics of the climactic phase with the two-phase flow model MFIX-TFM under axisymmetric conditions with several combinations of mass eruption rate, jet water content, vent size, particle size and density, topography, and emission duration. Results suggest that a steady mass eruption rate of 1.2–3.9 × 1011 kg/s is consistent with field observations. Minimal stratospheric injections occur in pulses issued from the central plume initially rising above the caldera center and from successive phoenix ash-clouds caused by the encounter of the pyroclastic density currents with topography. Most of the volcanic gas is injected into the stratosphere by the buoyant liftoff of dilute parts of the currents at the end of the eruption. Overall, 58–64 wt% of the total amount of gas emitted reaches the stratosphere. A fluctuating emission rate or an efficient final liftoff due to seawater interaction is unlikely to have increased this loading. Combined with petrological estimates of the degassed S, our results suggest that the eruption injected 11–20 Tg S into the stratosphere, consistent with the subsequent climate response and Greenland ice sheet deposition. Our results also show that the combination of the source Richardson number and the mass eruption rate is able to characterize the buoyant–collapse transition at Okmok. We extended this result to 141 runs from 10 published numerical studies of eruptive jets and found that this regime diagram is able to capture the first-order layout of the buoyant–collapse transition in all studies except one. An existing multivariate criterion yields the best predictions of this regime transition. [ABSTRACT FROM AUTHOR]

To extract climatically relevant chemical signals from the deepest, oldest Antarctic ice, we must first investigate the degree to which chemical ions diffuse within solid ice. Volcanic sulfate peaks are an ideal target for such an investigation because they are high-amplitude, short-duration (∼3 years) events with a quasi-uniform structure. Here we present an analysis of the EPICA Dome C sulfate record over the last 450 kyr. We identify volcanic peaks and isolate them from the non-sea-salt sulfate background to reveal the effects of diffusion: amplitude damping and broadening of peaks in the time domain with increasing depth and age. Sulfate peak shape is also altered by the thinning of ice layers with depth that results from ice flow. Both processes must be simulated to derive effective diffusion rates. This is achieved by running a forward model to diffuse idealised sulfate peaks at different rates while also accounting for ice thinning. Our simulations suggest a median effective diffusion rate of sulfate ions of 2.4±1.7×10-7 m 2 yr -1 in Holocene ice, slightly faster than suggested by previous work. The effective diffusion rate observed in deeper ice is significantly lower, and Holocene ice shows the highest rate of the last 450 kyr. Beyond the Holocene, there is no systematic difference between the effective diffusion rates of glacial and interglacial periods despite variations in soluble ion concentrations, dust loading, and ice grain radii. Effective diffusion rates for 40 to 200 ka are relatively constant and of the order 1×10-8 m 2 yr -1. Our results suggest that the diffusion of sulfate ions within volcanic peaks is relatively fast initially, perhaps through an interconnected vein network, but slows significantly after 40 kyr. In the absence of clear evidence for a controlling influence of temperature on sulfate diffusivity with depth and age, we hypothesise that the rapid decrease in effective diffusion rate from the time of deposition to ice of 50 ka age may be due to a switch in the mechanism of diffusion resulting from the changing location of sulfate ions within the ice microstructure and/or interconnectedness of veins and grain boundaries. [ABSTRACT FROM AUTHOR]

Recent technological advancements in spectral imaging core-scanning techniques have proved to be a promising tool to study lake sediments at extremely high resolution. We used this novel analytical approach to scan core AU3 of the Pleistocene Auel maar, Western Germany. The resulting ultra-high-resolution RABD670 spectral index, a proxy for the lake's primary production, shows an almost complete succession of Greenland Interstadials of the NGRIP ice core chronology back to around 60,000 years. Using the ELSA-20 chronology and its anchor points to the NGRIP record as a stratigraphic basis, we were able to compare and fine-tune prominent climate signals occurring in both regions. This in-depth correlation yields strong evidence that the climates of Greenland and Central Europe were not only strongly coupled on timescales of stadials and interstadials but even on multidecadal scales, showing prominent climate cycles between 20 and 125 years. As climate changes in these regions were ultimately driven by variations in the North Atlantic meridional heat transport, their strong coupling becomes most apparent during cold and arid intervals. In contrast, longer-lasting warmer and more humid phases caused the activation of various regional feedback mechanisms (e.g., soil formation, forest growth), resulting in more complex patterns in the proxy records. [ABSTRACT FROM AUTHOR]

Biomass burning plays an important role in climate-forcing and atmospheric chemistry. The drivers of fire activity over the past two centuries, however, are hotly debated and fueled by poor constraints on the magnitude and trends of preindustrial fire regimes. As a powerful tracer of biomass burning, reconstructions of paleoatmospheric carbon monoxide (CO) can provide valuable information on the evolution of fire activity across the preindustrial to industrial transition. Here too, however, significant disagreements between existing CO records currently allow for opposing fire histories. In this study, we reconstruct a continuous record of Antarctic ice coreCObetween 1821 and 1995CE to overlap with direct atmospheric observations. Our record indicates that the Southern Hemisphere CO burden ([CO]) increased by 50% from a preindustrial mixing ratio of ca. 35 ppb to ca. 53 ppb by 1995 CE with more variability than allowed for by state-of-the-art chemistry-climate models, suggesting that historic CO dynamics have been not fully accounted for. Using a 6-troposphere box model, a 40 to 50% decrease in Southern Hemisphere biomass-burning emissions, coincident with unprecedented rates of early 20th century anthropogenic land-use change, is identified as a strong candidate for this mismatch. [ABSTRACT FROM AUTHOR]

Radio-echo sounding is a standard technique for imaging the englacial stratigraphy of glaciers and ice sheets. In most cases, internal reflection horizons (IRHs) represent former glacier surfaces, comprise information about past accumulation and ice deformation, and enable the linking of ice core chronologies. IRHs in the lower third of the ice column are often difficult to detect or coherently trace. In the polar ice sheets, progress in IRH detection has been made by using multistatic, phase-coherent radars, enabling focused synthetic aperture radar (SAR) processing. However, these radar systems are often not suitable for deployment on mountain glaciers. We present a proof-of-concept study for a lightweight, phase-coherent and ground-based radar system, based on the phase-sensitive radio-echo sounder (pRES). To improve the detectability of IRHs we additionally adapted a layer-optimized SAR processing scheme to this setup. We showcase the system capability at Colle Gnifetti, Switzerland and Italy, where specular reflections are now apparent down to the base of the glacier. Compared to previously deployed impulse radar systems, with the mobile pRES the age of the oldest continuously traceable IRH could be increased from 78±12 to 288±35 a. Corresponding reflection mechanisms for this glacier are linked to stratified acidic impurities which in the upper part were deposited at a higher rate due to increased industrial activity in the area. Possible improvements to the system are discussed. If successfully implemented, these may provide a new way to map the deep internal structure of Colle Gnifetti and other mountain glaciers more extensively in future deployments. [ABSTRACT FROM AUTHOR]

While the processes controlling pore closure are broadly understood, the physical mechanisms driving the associated elemental fractionation remains ambiguous. Previous studies have shown that the pore closure process leads to a depletion in small-sized molecules (e.g. H2 , O2 , Ar, Ne, He) in ice core bubbles relative to larger-sized molecules like N2. This size-dependent fractionation, identified using ice core δ(O2/N2) records, exhibits a clear anti-correlation with local summer solstice insolation, making δ(O2/N2) a valuable ice core dating tool. Mechanisms controlling this relationship are attributed to the physical properties of deep firn. In this study, we compile δ(O2/N2) records from 15 polar ice cores and show a new additional link between δ(O2/N2) and local surface temperature and/or accumulation rate. Using the Crocus snowpack model, we perform sensitivity tests to identify the response of near-surface snow properties to changes in insolation intensity, accumulation rate, and air temperature. These tests support a mechanism linked to firn grain size, such that the larger the grain size for a given density, the stronger the pore closure fractionation and, hence, the lower the δ(O2/N2) values archived in the ice. Based on both snowpack model outputs and data compilation, our findings suggest that local accumulation rate and temperature should be considered when interpreting δ(O2/N2) as a local insolation proxy. [ABSTRACT FROM AUTHOR]

Accurate estimates of water isotope diffusion lengths are crucial when reconstructing and interpreting water isotope records from ice cores. This is especially true in the deepest, oldest sections of deep ice cores, where thermally enhanced diffusive processes have acted over millennia on extremely thinned ice. Previous statistical estimation methods, used with great success in shallower, younger ice cores, falter when applied to these deep sections, as they fail to account for the statistics of the climate on millennial timescales. Here, we present a new method to estimate the diffusion length from water isotope data and apply it to the Marine Isotope Stage 19 (MIS 19) interglacial at the bottom of the EPICA Dome C (EDC, Dome Concordia) ice core. In contrast to the conventional estimator, our method uses other interglacial periods taken from further up in the ice core to estimate the structure of the variability before diffusion. Through use of a Bayesian framework, we are able to constrain our fit while propagating the uncertainty in our assumptions. We estimate a diffusion length of 31±5 cm for the MIS 19 period, which is significantly smaller than previously estimated (40–60 cm). Similar results were obtained for each interglacial used to represent the undiffused climate signal, demonstrating the robustness of our estimate. Our result suggests better preservation of the climate signal at the bottom of EDC and likely other deep ice cores, offering greater potentially recoverable temporal resolution and improved reconstructions through deconvolution. [ABSTRACT FROM AUTHOR]

We document the isotopic evolution of near-surface snow at the East Greenland Ice Core Project (EastGRIP) ice core site in northeast Greenland using a time-resolved array of 1 m deep isotope (δ18O , δD) profiles. The snow profiles were taken from May–August during the 2017–2019 summer seasons. An age–depth model was developed and applied to each profile, mitigating the impacts of stratigraphic noise on isotope signals. Significant changes in deuterium excess (d) are observed in surface snow and near-surface snow as the snow ages. Decreases in d of up to 5 ‰ occur during summer seasons after deposition during two of the three summer seasons observed. The d always experiences a 3 ‰–5 ‰ increase after aging 1 year in the snow due to a broadening of the autumn d maximum. Models of idealized scenarios coupled with prior work indicate that the summertime post-depositional changes in d (Δd) can be explained by a combination of surface sublimation, forced ventilation of the near-surface snow down to 20–30 cm , and isotope-gradient-driven diffusion throughout the column. The interannual Δd is also partly explained with isotope-gradient-driven diffusion, but other mechanisms are at work that leave a bias in the d record. Thus, d does not just carry information about source-region conditions and transport history as is commonly assumed, but also integrates local conditions into summer snow layers as the snow ages through metamorphic processes. Finally, we observe a dramatic increase in the seasonal isotope-to-temperature sensitivity, which can be explained solely by isotope-gradient-driven diffusion. Our results are dependent on the site characteristics (e.g., wind, temperature, accumulation rate, snow properties) but indicate that more process-based research is necessary to understand water isotopes as climate proxies. Recommendations for monitoring and physical modeling are given, with special attention to the d parameter. [ABSTRACT FROM AUTHOR]

Galactic cosmic rays (GCRs) interact with matter in the atmosphere and at the surface of the Earth to produce a range of cosmogenic nuclides. Measurements of cosmogenic nuclides produced in surface rocks have been used to study past land ice extent as well as to estimate erosion rates. Because the GCR flux reaching the Earth is modulated by magnetic fields (solar and Earth's), records of cosmogenic nuclides produced in the atmosphere have also been used for studies of past solar activity. Studies utilizing cosmogenic nuclides assume that the GCR flux is constant in time, but this assumption may be uncertain by 30 % or more. Here we propose that measurements of 14 C of carbon monoxide (14 CO) in ice cores at low-accumulation sites can be used as a proxy for variations in GCR flux on timescales of several thousand years. At low-accumulation ice core sites, 14 CO in ice below the firn zone originates almost entirely from in situ cosmogenic production by deep-penetrating secondary cosmic ray muons. The flux of such muons is almost insensitive to solar and geomagnetic variations and depends only on the primary GCR flux intensity. We use an empirically constrained model of in situ cosmogenic 14 CO production in ice in combination with a statistical analysis to explore the sensitivity of ice core 14 CO measurements at Dome C, Antarctica, to variations in the GCR flux over the past ≈ 7000 years. We find that Dome C 14 CO measurements would be able to detect a linear change of 6 % over 7 ka, a step increase of 6 % at 3.5 ka or a transient 100-year spike of 190 % at 3.5 ka at the 3 σ significance level. The ice core 14 CO proxy therefore appears promising for the purpose of providing a high-precision test of the assumption of GCR flux constancy over the Holocene. [ABSTRACT FROM AUTHOR]

In situ measurements of water vapour isotopic composition in polar regions has provided needed constrains of post-deposition processes involved in the archiving of the climatic signal in ice core records. During polar winter, the temperatures, and thus the specific humidity, are so low that current commercial techniques are not able to measure the vapour isotopic composition with enough precision. Here, we make use of new developments in infrared spectroscopy and combine an optical-feedback frequency-stabilised laser source (OFFS technique) using a V-shaped cavity optical feedback (VCOF) cavity and a high-finesse cavity ring-down spectroscopy (CRDS) cavity to increase the signal-to-noise ratio while measuring absorption transitions of water isotopes. We present a laboratory infrared spectrometer leveraging all these techniques dedicated to measure water vapour isotopic composition at low humidity levels. At 400 ppmv , the instrument demonstrates a precision of 0.01 ‰ and 0.1 ‰ in δ18O and d-excess, respectively, for an integration time of 2 min. This set-up yields an isotopic composition precision below 1 ‰ at water mixing ratios down to 4 ppmv , which suggests an extrapolated precision in δ18O of 1.5 ‰ at 1 ppmv. Indeed, thanks to the stabilisation of the laser by the VCOF, the instrument exhibits extremely low drift and very high signal-to-noise ratio. The instrument is not hindered by a strong isotope–humidity response which at low humidity can create extensive biases on commercial instruments. [ABSTRACT FROM AUTHOR]

We analyze the temperature time series of the EPICA Dome C ice cores in Antarctica and of the Greenland project, Summit, with durations of 800 000 and 248 000 years, respectively, with a recent mathematical tool defined through the Fourier phases of the series, known as the J -index. This data driven index can differentiate between purely random dynamics and dynamics with a deterministic component. It is sensitive to nonlinear components and robust to the presence of noise. Our J -index data analysis shows that both Greenland and Antarctica climatic fluctuations possess deterministic traits and suggests the presence of an underlying nonlinear dynamics. Furthermore, in both regions, it reveals the simultaneous occurrence of an important global event known as the "Pelukian transgression." For Antarctica, it also detects the marine isotopic stage 11. Additionally, our calculation of the time series Hurst exponents and our detrended fluctuation analysis show the presence of long-range persistent correlations for Antarctica and anti-persistent correlations for Greenland. For the latter case, our fractal dimension determinations are indicative of a more complex climatic dynamics in Greenland with respect to Antarctica. Our results are encouraging for further development of climate variability deterministic models for these regions. [ABSTRACT FROM AUTHOR]

The isotopic signature (δ18О, δ2Н, d-excess) of the glacier ice in Ushkovskii volcano (Kamchatka Peninsula) was studied. A new shallow ice core was obtained in 2022 in the Gorshkov crater. The 14 m long ice core was dated by counting the annual layers, which were also compared with known eruptions in recent years. The upper 14 m of the glacier were formed over the last 16 years (from 2006 to 2022). The values of δ18О vary from −16 to −24‰, and the values of δ2Н, from −110.5 to −177.7‰, while the mean values are −20.5 and −150.2‰, respectively. The d-excess values vary with depth from 8.7 to 21.3‰ at a mean value of 13.7‰. For winter horizons at low values of δ18О and δ2Н, an increase in d-excess is noted. Such features are associated with the origin of moisture brought to Kamchatka. The source of moisture is the Pacific Ocean, the Sea of Okhotsk, and the Sea of Japan, for which there are pronounced differences in the conditions of moisture evaporation between summer and winter. The trends of increases in the δ18О and δ2Н values from 2006 to 2022 are accompanied by a decrease in the deuterium excess, indicating an increase in summer precipitation. However, in addition to changes in the seasonal proportions of precipitation, the d-excess values of ice may reflect climatic changes in the source of moisture. [ABSTRACT FROM AUTHOR]

Carbonyl sulfide (COS) is the most abundant sulfur gas in the atmosphere with links to terrestrial and oceanic productivity. We measured COS in ice core air from an intermediate-depth ice core from the South Pole using both dry and wet extraction methods, recovering a 52 500-year record. We find evidence for COS production in the firn, altering the atmospheric signal preserved in the ice core. Mean sea salt aerosol concentrations from the same depth are a good proxy for the COS production, which disproportionately impacts the measurements from glacial period ice with high sea salt aerosol concentrations. The COS measurements are corrected using sea salt sodium (ssNa) as a proxy for the excess COS resulting from the production. The ssNa-corrected COS record displays substantially less COS in the glacial period atmosphere than the Holocene and a 2 to 4-fold COS rise during the deglaciation synchronous with the associated climate signal. The deglacial COS rise was primarily source driven. Oceanic emissions in the form of COS, carbon disulfide (CS2), and dimethylsulfide (DMS) are collectively the largest natural source of atmospheric COS. A large increase in ocean COS emissions during the deglaciation suggests enhancements in emissions of ocean sulfur gases via processes that involve ocean productivity, although we cannot quantify individual contributions from each gas. [ABSTRACT FROM AUTHOR]

Microplastic pollution is a pervasive issue, with remarkably high concentrations observed even in the most remote locations such as Arctic sea ice and snow. The reason for such large microplastic abundances in sea ice is still speculative and applies mainly to saline or freshwater conditions. In this study, we investigated seasonal ice core samples collected in March 2021 from the northern Baltic Sea (Gulf of Bothnia) for their microplastic distributions. The Baltic Sea is characterized by low salinity and can be ice-covered for up to six months annually. Microplastics were analyzed in the melted ice samples using an adsorption technique and Raman microscopy to identify their abundances, colors, shapes, and sizes to calculate their masses. Due to the strong dynamic of the ice layer and the repeated melting and freezing processes during the ice formation, no discernible trends in microplastic abundances, masses, or polymer types were observed throughout the ice core length. The average microplastic abundance (±SD) in the Baltic Sea ice was determined to be 22.3 ± 8.6 N L−1, with 64.9% of the particles exhibiting a particulate shape and 35.1% having a fibrous shape. The most prevalent polymer type was polyethylene terephthalate (PET), accounting for 44.4% of all polymers. This is likely due to the high proportion of PET fibers (93.8%). The majority of particle-shaped microplastics were identified as polyethylene (PE; 37.2%), followed by PET (17.2%), polyvinyl chloride (PVC; 15.9%), and polypropylene (PP; 15.9%). No correlations were found between microplastic concentrations and proximity to land, cities, industries, or rivers, except for PP mass concentrations and particle sizes, which correlated with distances to industries in Luleå, Sweden. [ABSTRACT FROM AUTHOR]

The obituary in the journal "Radiocarbon" commemorates the life and work of Mike Baillie, a Belfast-born scientist known for his contributions to dendrochronology and radiocarbon dating. Baillie's research focused on building oak and pine tree-ring chronologies in Ireland, which became the basis for radiocarbon calibration curves. He also explored the environmental impacts on trees and societies, suggesting links between narrow-ring events and volcanic dust or cosmic origins. Baillie's meticulous work and interdisciplinary approach left a lasting legacy in the fields of science, archaeology, and history. [Extracted from the article]

Ice-core records from the interior of the Greenland ice sheet suggest widespread thinning during the Holocene. However, the recurring underestimation of this thinning in numerical models raises concerns about both the veracity of such reconstructions and the reliability of glaciological models. Recent work suggests the 8000-year-old Northeast Greenland Ice Stream (NEGIS), including a now-extinct northern tributary, may have been an early influence on Greenland ice-sheet dynamics. Yet, the inaccurate reproduction of NEGIS-like dynamics in most models hampers investigation of whether this feature played a role in Holocene ice-sheet thinning. Here we show that grounding-line retreat in northeast Greenland triggers elevation changes at the northern summit via ice-dynamic effects modulated by the paleo NEGIS system. In our simulations, fast ice-stream flow caused by transiently imposed reduced basal shear stress following the northeast retreat explains 55% (± 18%) of the estimated ice thinning, showing that ice-stream dynamics is one of the main drivers of the NGRIP Holocene surface elevation drop. Our findings show that the ice-flow in northeast Greenland plays a large role in ice-surface elevation changes in central Greenland. Using a three-dimensional ice sheet model, this work demonstrates that more than half of the Holocene surface elevation drop recorded in ice cores from north-central Greenland can be attributed to ice dynamic changes in northeast Greenland following ice-sheet retreat during the last deglaciation. [ABSTRACT FROM AUTHOR]

Black carbon in snow (BCS) is a crucial parameter in Earth System modeling, as it influences global radiative balance. Here, simulated BCS from Coupled Model Intercomparison Project Phase 5 and 6 (CMIP5 and CMIP6) that provided BCS as a model output were evaluated. In comparison with global BCS observations, CMIP5/6 models successfully reproduced long‐term historical trends linked to human activities, but struggled capturing decadal variability caused by natural climate variability. CMIP6 models NorESM2‐MM, NorESM2‐LM, and TaiESM1 yielded the most accurate simulations of BCS concentration with modest overestimation of <50%, while the four CESM2 models underestimated concentrations by up to ∼80%. These errors effectively balanced for the CMIP6 multi‐model ensemble mean (MME), which had a relative error (RE) of −37%. However, the CMIP5 MME was less reliable due to extreme overestimation by up to 8,000% in the three MIROC models. The significant BCS concentration errors in the MIROC and CESM2 models were linked mainly to errors in handling of BC in snow processes. Conversely, marked improvements in NorESM, the only common to both CMIP5 and CMIP6, were due to improved simulation of black carbon deposition. BCS errors significantly impacted radiative forcing estimates, particularly at the poles, where model errors reached several thousandfold. CMIP6 exhibited superior results compared to CMIP5, achieving global MME RE of −33% in radiative forcing estimates. However, it's worth noting BCS output is currently limited, with only seven models available for each of CMIP5 and CMIP6 here. Additional models simulating BCS are desirable in the next CMIP generations. Plain Language Summary: Black carbon in snow (BCS) has important implications for snow albedo reduction, solar radiation absorption, and snow melt acceleration on the Earth. However, model performance in BCS simulation remains unclear in the widely used Coupled Model Intercomparison Project phases 5 and 6 (CMIP5/6). In this study, CMIP5 and CMIP6 models providing BCS information were evaluated against a global data set of BCS observations from ice cores and snow samples. Results indicate that both CMIP generations successfully reproduce long‐term trends in BCS but fail to capture the decadal variability. Marked differences exist in modeled BCS concentrations, with the best‐performing models (NorESM2‐LM, NorESM2‐MM, and TaiESM1) overestimating concentrations by <50%, while the three worst‐performing MIROC models yield extreme overestimation. Such discrepancies stem primarily from the post‐depositional treatment of black carbon in snow, which influences the estimation of climate‐warming effects of BCS. Key Points: The global pattern of black carbon in snow (BCS) is effectively simulated by CMIP5/6 models, but large errors remainPost‐depositional snow processes contribute most to the large simulated‐BCS errorsBCS‐induced radiative forcing in Antarctica is poorly simulated by CMIP5/6 models with an overestimation of ∼200,000% and 144%, respectively [ABSTRACT FROM AUTHOR]

While the simple model of the total atmospheric carbon sink effect as a linear function of concentration has provided excellent prediction results, several problems remained to be investigated and solved. The most obvious open issue is the correct treatment of land use change emissions. It turns out that the model improves by mostly neglecting these emissions after 1950. This effectively implies that land use change emissions have been constant and small since then. The key investigation starts with the observation that the total carbon sink has a short-term component that can be explained by temperature changes. The apparent paradox, why contrary to the short-term changes no temperature-caused trend can be detected, despite the fact that several contributing processes exhibit clear temperature dependency, is analyzed and explained. The result of this analysis leads to the model extension, where the total effect of absorptions and natural emissions are a linear function of concentration and temperature. This extended model not only explains current measurements but also paleo-climate data from ice core time series. [ABSTRACT FROM AUTHOR]

We report on the measurement of the index of refraction of glacial ice at radio frequencies at Summit Station, Greenland. This measurement is of particular importance for the Radio Neutrino Observatory Greenland, an experiment currently under construction at Summit Station that seeks to detect radio signals from ultra-high-energy neutrino interactions in the ice. By correlating radio reflections in the bulk ice with features in the conductivity measurements from ice cores, we determine the index of refraction as n=1.778±0.006. [ABSTRACT FROM AUTHOR]

Measurements of carbon-14-containing carbon monoxide (14 CO) in glacial ice are useful for studies of the past oxidative capacity of the atmosphere as well as for reconstructing the past cosmic ray flux. The 14 CO abundance in glacial ice represents the combination of trapped atmospheric 14 CO and in situ cosmogenic 14 CO. The systematics of in situ cosmogenic 14 CO production and retention in ice are not fully quantified, posing an obstacle to interpretation of ice core 14 CO measurements. Here we provide the first comprehensive characterization of 14 CO at an ice accumulation site (Summit, Greenland), including measurements in the ice grains of the firn matrix, firn air and bubbly ice below the firn zone. The results are interpreted with the aid of a firn gas transport model into which we implemented in situ cosmogenic 14 C. We find that almost all (≈ 99.5 %) of in situ 14 CO that is produced in the ice grains in firn is very rapidly (in <1 year) lost to the open porosity and from there mostly vented to the atmosphere. The timescale of this rapid loss is consistent with what is expected from gas diffusion through ice. The small fraction of in situ 14 CO that initially stays in the ice grains continues to slowly leak out to the open porosity at a rate of ≈ 0.6 % yr -1. Below the firn zone we observe an increase in 14 CO content with depth that is due to in situ 14 CO production by deep-penetrating muons, confirming recent estimates of 14 CO production rates in ice via the muon mechanisms and allowing for narrowing constraints on these production rates. [ABSTRACT FROM AUTHOR]

Scientific drilling in polar regions plays a crucial role in obtaining ice cores and using them to understand climate change and to study the dynamics of polar ice sheets and their impact on global environmental changes (sea level, ocean current cycle, atmospheric circulation, etc.). Mechanical rotary cutting is a widely used drilling method that drives the cutter to rotate to cut and drill through ice layers. It is necessary to conduct in-depth research on the brittle fracture behavior of ice and mechanical model and to analyze the factors and specific mechanisms (cutter's angle, rotation speed of the drill bit, and cutting depth) affecting cutting force for the rational design of ice core drill systems, improving the efficiency of ice core drilling and ensuring the drilling process runs smoothly. Therefore, in this paper, the process of ice cutting was observed, the fracture mechanics characteristics of the ice cutting process were analyzed, the formation process of ice chips was divided into three stages, and a mathematical model for the cutting force was established based on the observation results. The paper describes the damage conditions of ice failure and points out the factors and specific laws influencing cutting force. Furthermore, the cutting force generated under various experimental conditions was tested. Based on specified real-time variation curves of cutting force, the characteristics of cutting force were analyzed during the cutting and drilling process. Based on comparison to results of the average cutting force, the influence mechanism of various parameters acting on the cutting force was obtained. This proves the correctness of the mathematical model of the cutting force and provides a theoretical reference for the calculation of the cutting force during ice cutting and drilling in polar regions. [ABSTRACT FROM AUTHOR]

Terrestrial and extraterrestrial radioisotope research has been strongly dependent on the development of analytical methods which would enable to trace radioisotopes at low concentrations in subgram samples (e.g., in tree rings, ice cores, meteorites, etc.). Accelerator mass spectrometry (AMS) has become the most sensitive technique for ultralow-level analysis of long-lived radioisotopes, such as 14C, 10Be and 26Al. We review developments and applications carried out in the CENTA laboratory, and describe a recently installed fully equipped AMS line, designed for analysis of long-lived radioisotopes from tritium to curium. [ABSTRACT FROM AUTHOR]

During the Last Glacial Maximum (LGM, ∼ 23 000 to 19 000 years ago), the Patagonian Ice Sheet (PIS) covered the central chain of the Andes between ∼ 38 to 55° S. Existing paleoclimate evidence – mostly derived from glacial landforms – suggests that maximum ice sheet expansions in the Southern Hemisphere and Northern Hemisphere were not synchronized. However, large uncertainties still exist in the timing of the onset of regional deglaciation and its major drivers. Here we present an ensemble of numerical simulations of the PIS during the LGM. We assess the skill of paleoclimate model products in reproducing the range of atmospheric conditions needed to enable an ice sheet growth in concordance with geomorphological and geochronological evidence. The resulting best-fit climate product is then combined with records from southern South America offshore sediment cores and Antarctic ice cores to drive transient simulations throughout the last 70 ka using a glacial index approach. Our analysis suggests a strong dependence of the PIS geometry on near-surface air temperature forcing. Most ensemble members underestimate the ice cover in the northern part of Patagonia, while tending to expand beyond its constrained eastern boundaries. We largely attribute these discrepancies between the model-based ice geometries and geological evidence to the low resolution of paleoclimate models and their prescribed ice mask. In the southernmost sector, evidence suggests full glacial conditions during marine isotope stage 3 (MIS3, ∼ 59 400 to 27 800 years ago), followed by a warming trend towards MIS2 (∼ 27 800 to 14 700 years ago). However, in northern Patagonia, this deglacial trend is absent, indicating a relatively consistent signal throughout MIS3 and MIS2. Notably, Antarctic cores do not reflect a glacial history consistent with the geochronological observations. Therefore, investigations of the glacial history of the PIS should take into account southern midlatitude records to capture effectively its past climatic variability. [ABSTRACT FROM AUTHOR]

The continuous ice core record extends 800 000 years into the past, covering the period of 100 000-year glacial cycles but not the transition from 40 000-year glacial cycles (the mid-Pleistocene transition, 1.2–0.7 million years ago). A primary goal of the International Partnerships in Ice Core Sciences is therefore to retrieve a 1.5-million-year-old continuous ice core, increasing our understanding of this major change in the climate system and thus of fundamental climate forcings and feedbacks. However, complex glacial processes, limited bedrock data, and young basal ice in previous cores necessitate careful reconnaissance studies before extracting a full core. Ice borehole optical logging reflects the ice dust content and may be used to date ice quickly and inexpensively if a reference record is known. Here we explore the relationship between ice dust records and well-dated marine dust records from sediment cores in the southern Atlantic and Pacific oceans, which lie along paths of dust sources to Antarctica. We evaluate how representative these records are of Antarctic dust both through the existing ice core record and during the older target age range, suggesting that a newly published 1.5-million-year record from Site U1537 near South America is likely the most robust predictor of the Oldest Ice dust signal. We then assess procedures for rapid dating of potential Oldest Ice sites, noting that the ability to detect dating errors is an essential feature. We emphasize that ongoing efforts to identify, recover, date, and interpret an Oldest Ice core should use care to avoid unfounded assumptions about the 40 kyr world based on the 100 kyr world. [ABSTRACT FROM AUTHOR]

Palaeoclimate time series, reflecting the state of Earth's climate in the distant past, occasionally display very large and rapid shifts showing abrupt climate variability. The identification and characterisation of these abrupt transitions in palaeoclimate records is of particular interest as this allows for understanding of millennial climate variability and the identification of potential tipping points in the context of current climate change. Methods that are able to characterise these events in an objective and automatic way, in a single time series, or across two proxy records are therefore of particular interest. In our study the matrix profile approach is used to describe Dansgaard–Oeschger (DO) events, abrupt warmings detected in the Greenland ice core, and Northern Hemisphere marine and continental records. The results indicate that canonical events DO-19 and DO-20, occurring at around 72 and 76 ka, are the most similar events over the past 110 000 years. These transitions are characterised by matching transitions corresponding to events DO-1, DO-8, and DO-12. They are abrupt, resulting in a rapid shift to warmer conditions, followed by a gradual return to cold conditions. The joint analysis of the δ18 O and Ca 2+ time series indicates that the transition corresponding to the DO-19 event is the most similar event across the two time series. [ABSTRACT FROM AUTHOR]