Shallow Cores Research Articles

Ice-core data used for the construction of the Greenland Ice-Core Chronology 2005 and 2021 (GICC05 and GICC21)

Abstract. We here describe, document, and make available a wide range of data sets used for annual-layer identification in ice cores from DYE-3, GRIP, NGRIP, NEEM, and EGRIP. The data stem from detailed measurements performed both on the main deep cores and shallow cores over more than 40 years using many different setups developed by research groups in several countries and comprise both discrete measurements from cut ice samples and continuous-flow analysis data. The data series were used for counting annual layers 60 000 years back in time during the construction of the Greenland Ice-Core Chronology 2005 (GICC05) and/or the revised GICC21, which currently only reaches 3800 years back. Now that the underlying data are made available (listed in Table 1) we also release the individual annual-layer positions of the GICC05 timescale which are based on these data sets. We hope that the release of the data sets will stimulate further studies of the past climate taking advantage of these highly resolved data series covering a large part of the interior of the Greenland ice sheet.

Internal Architectural Patterns of Bar Fingers Within Digitate Shallow-Water Delta: Insights from the Shallow Core, GPR and Delft3D Simulation Data of the Ganjiang Delta, China

Abstract Digitate shallow-water deltas are commonly found in modern lakes and bays, as well as within cratonic petroliferous basins. They develop one or multiple sinuous finger-like sands (i.e., bar fingers), including high-RSI (sinuosity ratio of distributary channel and bar finger ≥1) and low-RSI (RSI < 1) types. Bar fingers consist of four types of subenvironments, that is, distributary channels, point bars, mouth bars, and levees. However, the internal architecture within the above subenvironments is still unclear. This paper documents the internal architecture of a digitate delta based on the integration of shallow-core and ground-penetrating radar data from the Ganjiang Delta, China, coupled with Delft3D simulations. Our results show that multiple convex-up muddy-silty accretion beds are developed in mouth bars, which top lap the side of the distributary channels or point bars and down lap the bottom of the mouth bar. The accretion beds have low dip angles (<2°), which is slightly higher for the upper accretion beds. Point bars, unique to the high-RSI bar finger, develop multiple inclined silty drapes, which top lap the top of the point bar. The cohesive levee and backwater effect impede the migration of the distributary channel, resulting in silty drapes with high-dip angles (can be >10°) compared with those in the supplying river. This dip angle exhibits a negative relationship with downstream distance and a positive exponential relationship with lateral migration distance. Silty drapes become dense along the migration direction of the distributary channel. The levee develops multiple horizontal muddy accretion beds. The high-RSI bar finger develops a large number (>3) of accretion beds in mouth bars with high dip angles, and a large number of accretion beds in thick levees, compared with the low-RSI bar finger. The results of this paper provide insights into the prediction and development of cratonic digitate shallow-water delta reservoirs.

Chemical bonding in americium oxides probed by X-ray spectroscopy

The electronic structure and the chemical state in Am binary oxides and Am-doped UO_2 were studied by means of X-ray absorption spectroscopy at shallow Am core (4d and 5d) edges. In particular, the Am 5f states were probed and the nature of their bonding to the oxygen states was analyzed. The interpretation of the experimental data was supported by the Anderson impurity model (AIM) calculations which took into account the full multiplet structure due to the interaction between 5f electrons as well as the interaction with the core hole. The sensitivity of the branching ratio of the Am 4d_{3/2} and 4d_{5/2} X-ray absorption lines to the chemical state of Am was shown using Am binary oxides as reference systems. The observed ratio for Am-doped UO_2 suggests that at least at low Am concentrations, americium is in the Am(III) state in the UO_2 lattice. To confirm the validity of the applied AIM approach, the analysis of the Am 4f X-ray photoelectron spectra of AmO_2 and Am_2O_3 was also performed which revealed a good agreement between experiment and calculations. As a whole, AmO_2 can be classified as the charge-transfer compound with the 5f occupancy (n_f) equal to 5.73 electrons, while Am_2O_3 is rather a Mott–Hubbard system with n_f = 6.05.

Impact of organic carbon reworking upon GDGT temperature proxies during the Paleocene-Eocene Thermal Maximum

Glycerol dialkyl glycerol tetraethers (GDGTs) have been widely applied to coastal marine sediments to reconstruct past temperature variability. However, coastal environments are characterised by variability in the source, age and/or thermal maturity of different organic carbon (OC) pools and may bias various GDGT-based proxies. Here we analyse TEX86 and MBT5MEvalues within a shallow marine sediment core (South Dover Bridge, Maryland) from the Paleocene-Eocene Thermal Maximum (PETM; 56 million years ago (Ma)) to explore how changes in OC reworking influence GDGT-derived sea surface and terrestrial temperature estimates, respectively. We demonstrate that TEX86 values are unaffected by an increase in soil- and fossil organic carbon during the PETM. In contrast, we find large and unexpected variations in MBT5ME-derived temperature estimates (∼6 to 25 °C) during the onset and core of the PETM at some sites. This coincides with input of reworked terrestrial OC from the Cenomanian-aged Raritan Formation. However, there is also an increase in the degree of cyclisation of tetramethylated branched GDGTs, suggesting that branched GDGTs are also derived from marine in-situ production. These factors preclude terrestrial temperature reconstructions at this site. We explored whether OC reworking is problematic in other PETM-aged coastal environments. Using GDGT metrics and the Branched and Isoprenoid GDGT Machine learning Classification algorithm (BIGMaC), we demonstrate that TEX86 values are mostly unaffected by changes in OC sources. However, MBT5ME values are affected by marine and/or freshwater overprints, especially in environments with low terrestrial OC input. Taken together, this study highlights the importance of constraining the provenance of different GDGTs in marine and lacustrine environments.

An upgraded CFA – FLC – MS/MS system for the semi-continuous detection of levoglucosan in ice cores

A new Continuous Flow Analysis (CFA) system coupled with Fast Liquid Chromatography – tandem Mass Spectrometry (FLC-MS/MS) has been recently developed for determining organic markers in ice cores. In this work we present an upgrade of this innovative technique, optimized for the detection of levoglucosan in ice cores, a crucial tracer for reconstructing past fires. The upgrade involved a specific optimization of the chromatographic and mass spectrometric parameters, allowing for a higher sampling resolution (down to 1 cm) and the simultaneous collection of discrete samples, for off-line analysis of water stable isotopes and additional chemical markers. The robustness and repeatability of the method has been tested by the analysis of multiple sticks of ice cut from the same shallow alpine ice core, and running the system for several hours on different days. The results show similar and comparable trends between the ice sticks. With this upgraded system, a higher sensitivity and a lower limit of detection (LOD) was achieved compared to discrete analysis of alpine samples for levoglucosan measurements. The new LOD was as low as 66 ng L−1, a net improvement over the previous LOD of 600 ng L−1.

Environmental Changes in Antarctica Using a Shallow Ice Core from Dronning Maud Land (DML), East Antarctica

Environmental Changes in Antarctica Using a Shallow Ice Core from Dronning Maud Land (DML), East Antarctica

Export pathways of biosolid derived microplastics in soil systems – Findings from a temperate maritime climate

The environmental fate of microplastics (MPs) added to agricultural soils remains poorly understood, particularly regarding their mobility in soils. Here we investigate the potential for MP export from soil to surface waters and groundwater in two agricultural settings with a 20-year history of biosolid treatment. A third site where biosolids had never been applied served as a reference (Field R). The potential for MP export along overland and interflow pathways to surface waters was determined from MP abundances in shallow surface cores (10 cm) along ten down-slope transects (five each for Field A and B), and through MP abundances in effluent from a sub-surface land drain. The risk of vertical MP migration was assessed from 2 m cores, and from MP abundances in groundwater sampled from the core boreholes. XRF Itrax core scanning was conducted on two of the deep cores to capture high-resolution optical and 2-D radiographic imaging. Results suggest limited MP mobility at depths >35 cm, with MPs largely recovered in surface soils characterised by lower compaction. Furthermore, abundances of MPs across the surface cores were comparable, with no evidence of MP accumulations observed. Average MP abundance in the top 10 cm of soil across Field A and B was 365 ± 302 MP kg−1, with 0.3 MP l−1 and 1.6 MP l−1 recovered from the groundwater and field drainpipe water samples, respectively. MP abundances were significantly higher in fields treated with biosolids than in Field R (90 ± 32 MP kg−1 soil). Findings suggest ploughing is the most significant driver of MP mobility in upper soil layers, however the potential for overland or interflow movement cannot be excluded, particularly for fields that may be artificially drained.

New data on palaeogeographical changes during the Holocene in the Varna Lake region

The article presents the results of an integrated study of geological data based on three shallow cores from the Varna Lake area. The sedimentary sequences and the related archaeological artifacts were examined through lithological descriptions, a biostratigraphical analysis of the molluscan fauna, and an analysis of the geomorphological, archaeological and radiocarbon data. Continental, liman and marine facies conditions of sedimentation have been established by examining the sedimentary textures, stratigraphic and palaoecological value of the molluscan subfossil fauna. This complex interpretation of data through different analyses allowed the highlighting of important stages in the palaeogeographical development of the Lake of Varna and in the human adaptation to changes in the coastal landscape at the beginning of the Holocene.

RXJ0437+00: constraining dark matter with exotic gravitational lenses

ABSTRACT We present the first strong-gravitational-lensing analysis of the galaxy cluster RX J0437.1+0043 (RXJ0437; z = 0.285). Newly obtained, deep MUSE observations, Keck/MOSFIRE near-infrared spectroscopy, and Hubble Space Telescope SNAPshot imaging reveal 13 multiply imaged background galaxies, three of them (at z = 1.98, 2.97, and 6.02, respectively) in hyperbolic umbilic (H–U) lensing configurations. The H–U images are located only 20–50 kpc from the cluster centre, i.e. at distances well inside the Einstein radius where images from other lens configurations are demagnified and often unobservable. Extremely rare (only one H–U lens was known previously) these systems are able to constrain the inner slope of the mass distribution – and unlike radial arcs, the presence of H–U configurations is not biased towards shallow cores. The galaxies lensed by RXJ0437 are magnified by factors ranging from 30 to 300 and (in the case of H–U systems) stretched nearly isotropically. Taking advantage of this extreme magnification, we demonstrate how the source galaxies in H–U systems can be used to probe for small-scale (∼109 M⊙) substructures, providing additional insight into the nature of dark matter.

Visual Stratigraphy-Based Age Scale Developed for the Shallow Mount Siple Firn Core, Antarctica

Here we present a case study for using visual stratigraphy to date a shallow firn core from coastal West Antarctica. The Mount Siple ice core has the potential to reconstruct climate in this data-sparse region over recent decades. Line scanned images of the 24 m firn core were used to generate a grey-scale, which displays variability consistent with annual cycles. The resulting Mount Siple age scale spans from 1998 ± 1 to 2017 CE. This study demonstrates that the seasonal changes in the grey-scale record provide an independent method of dating firn cores. However, the presence of melt layers at this site has introduced an error of ±1 year. Visual line stratigraphy has the unique advantage over traditional annual layer counting, based on chemical or isotopic species, of being non-destructive and relatively inexpensive. Visual line stratigraphy has proved to be an effective dating method for this site.

Microplastic accumulation in endorheic river basins – The example of the Okavango Panhandle (Botswana)

The Okavango Panhandle is the main influent watercourse of the Okavango Delta, an inland sink of the entire sediment load of the Cubango-Okavango River Basin (CORB). The sources of pollution in the CORB, and other endorheic basins, are largely understudied when compared to exorheic systems and the world's oceans. We present the first study of the distribution of microplastic (MP) pollution in surface sediments of the Okavango Panhandle in Northern Botswana. MP concentrations (64 μm-5 mm size range) in sediment samples from the Panhandle range between 56.7 and 399.5 particles kg−1 (dry weight) when analysed with fluorescence microscopy. The concentrations of MP in the 20 μm to 5 mm grain size range (analysed with Raman spectroscopy) range between 1075.7 and 1756.3 particles kg−1. One shallow core (15 cm long) from an oxbow lake suggests that MP size decreases with depth while MP concentration increases downcore. Raman Spectroscopy revealed that the compositions of the MP are dominated by polyethene terephthalate (PET), polypropylene (PP), polyethene (PE), polystyrene (PS), and polyvinyl chloride (PVC). From this novel data set it was possible to estimate that 10.9–336.2 billion particles could be transported into the Okavango Delta annually, indicating that the region represents a significant sink for MP, raising concerns for the unique wetland ecosystem.

Conservation Paleobiology and Taphonomy: Differential Preservation and Time-Averaging of Echinoids and Mollusks

Fossils from surficial death assemblages and shallow cores are an important archive used in the field of conservation paleobiology. Understanding the taphonomic filters and time-averaging that affect modern biomineralized taxa is crucial for using their fossil record as a source of geohistorical data. Through comparative analysis of the live-dead patterns of echinoids and mollusks, we aim to assess multiple hypotheses regarding differences between the two taxa in terms of preservation potential and fossil record resolution. In this pilot study, six sediment samples (five surface [0-10cm] and one subsurface [30-40cm]) were collected in a shallow, subtidal habitat off the coast of Cedar Key, Florida. These samples were sieved for mollusks and echinoids. Specimens were segregated by taxa and classified into live, whole dead, and fragments. Results were consistent among surface samples, which showed that, by weight, an average of 8% of all mollusk specimens and 55% of all echinoid specimens were live collected. Fragmentation rate by weight was notable higher for echinoids than for mollusks: 77% for mollusks and >99% for echinoids. The subsurface sample lacked complete echinoid tests, but complete mollusks shells, mollusk fragments and echinoid fragments remained common. These results support the hypothesis that echinoid tests degrade more quickly than mollusk shells, making them less likely to be preserved as part of time-averaged assemblages formed under low net-accumulation rates. Echinoid fragments do not seem to follow this pattern, and instead seem to preserve similarly to mollusks. The results suggest that geohistorical records provided by echinoids and mollusks may be fundamentally distinct from one another in terms of temporal resolution and completeness of the fossil record.

Spatial Variations of Fabric and Microstructure of Blue Ice Cores at the Shear Margin of Dalk Glacier, Antarctica

The study of the fabric and microstructure of ice at the shear margin of the Antarctic ice sheet is of great significance for understanding the ice flow and its contributions to sea level rise. In this study, twenty-three one-meter-long ice cores were drilled from blue ice areas at the shear margin of the Dalk Glacier, Antarctica. The ice fabric and microstructure of these ice cores are analyzed using a G50 fabric analyzer. This study shows that the shallow ice cores in this region present a cluster fabric as a consequence of shear stress. The grain size decreases following the direction of the ice flow towards the exposed bedrock at the end of the glacier, due to the blocking and squeezing by the bedrock. The formation mechanism of the shallow ice layers is that the ice from the original accumulation area flows here, lifted by the bedrock and shaped by the summer ablation and denudation. The basal ice at the shear margin of the Dalk Glacier is strongly rubbed by the bedrock and demonstrates a cluster fabric. The analysis of stable water isotopes shows a weak negative correlation between shallow ice fabric and stable water isotopes with depth. Bedrock topography and shear stress have a greater influence on grain microstructure among different ice cores over long distances at shear margins.

Temporal variations of surface mass balance over the last 5000 years around Dome Fuji, Dronning Maud Land, East Antarctica

Abstract. We reconstructed surface mass balance (SMB) around Dome Fuji, Antarctica, over the last 5000 years using the data from 15 shallow ice cores and seven snow pits. The depth–age relationships for the ice cores were determined by synchronizing them with a layer-counted ice core from West Antarctica (WAIS Divide ice core) using volcanic signals. The reconstructed SMB records for the last 4000 years show spatial patterns that may be affected by their locations relative to the ice divides around Dome Fuji, proximity to the ocean, and wind direction. The SMB records from the individual ice cores and snow pits were stacked to reconstruct the SMB history in the Dome Fuji area. The stacked record exhibits a long-term decreasing trend at -0.037±0.005 kg m−2 per century over the last 5000 years in the preindustrial period. The decreasing trend may be the result of long-term surface cooling over East Antarctica and the Southern Ocean and sea ice expansion in the water vapor source areas. The multidecadal to centennial variations of the Dome Fuji SMB after detrending the record shows four distinct periods during the last millennium: a mostly negative period before 1300 CE, a slightly positive period from 1300 to 1450 CE, a slightly negative period from 1450 to 1850 CE with a weak maximum around 1600 CE, and a strong increase after 1850 CE. These variations are consistent with those of previously reconstructed SMB records in the East Antarctic plateau. The low accumulation rate periods tend to coincide with the combination of strong volcanic forcings and solar minima for the last 1000 years, but the correspondence is not clear for the older periods, possibly because of the lack of coincidence of volcanic and solar forcings or the deterioration of the SMB record due to a smaller number of stacked cores.

A record of volcanic eruptions over the past 2,200 years from Vostok firn cores, central East Antarctica

Introduction: The products of volcanic eruptions found in the snow, firn and ice deposits of the polar ice sheets are precious sources of information on the volcanic forcing of the climate system in the recent or remote past. On the other hand, the layers containing the traces of well-known eruptions serve as absolute age markers that help to construct the depth-age scale for the snow-firn thickness.Methods: In this study we present new records of the sulfate concentrations and electrical conductivity (ECM) from three shallow (up to 70 m depth) firn cores drilled in the vicinity of Vostok station (central East Antarctica).Results: In the non-sea-salt sulfate and ECM profiles we were able to identify 68 peaks that can be interpreted as traces of volcanic events.Discussion: 22 of these peaks can be unambiguously attributed to well-known volcanic eruptions (including Tambora 1816 CE, Huaynaputina 1601 CE, Samalas 1258 CE, Ilopango 541 CE and others), which allowed to construct a robust depth-age scale for the cores. 37 events have their counterparts in other Antarctic cores, but cannot be associated with welldated eruptions. Finally, 9 peaks do not have analogues in the other cores, i.e., they may be traces of so far unknown volcanic events. According to the newly constructed depth-age function, the deepest studied firn layers (70.20 m) are dated by 192 BCE.

Comparison of Three Methods for Distinguishing Glacier Zones Using Satellite SAR Data

Changes in glacier zones (e.g., firn, superimposed ice, ice) are good indicators of glacier response to climate change. There are few studies of glacier zone detection by SAR that are focused on more than one ice body and validated by terrestrial data. This study is unique in terms of the dataset collected—four C- and L-band quad-pol satellite SAR images, Ground Penetrating Radar data, shallow glacier cores—and the number of land ice bodies analyzed, namely, three tidewater glaciers in Svalbard and one ice cap in Iceland. The main aim is to assess how well popular methods of SAR analysis perform in distinguishing glacier zones, regardless of factors such as the morphologic differences of the ice bodies, or differences in SAR data. We test and validate three methods of glacier zone detection: (1) Gaussian Mixture Model–Expectation Maximization (GMM-EM) clustering of dual-pol backscattering coefficient (sigma0); (2) GMM-EM of quad-pol Pauli decomposition; and (3) quad-pol H/α Wishart segmentation. The main findings are that the unsupervised classification of both sigma0 and Pauli decomposition are promising methods for distinguishing glacier zones. The former performs better at detecting the firn zone on SAR images, and the latter in the superimposed ice zone. Additionally, C-band SAR data perform better than L-band at detecting firn, but the latter can potentially separate crevasses via the classification of sigma0 or Pauli decomposition. H/α Wishart segmentation resulted in inconsistent results across the tested cases and did not detect crevasses on L-band SAR data.

Microfaunal Recording of Recent Environmental Changes in the Herschel Basin, Western Arctic Ocean

ABSTRACTMicrofaunal assemblages of benthic foraminifera, ostracods, and tintinnids from two marine sediment cores retrieved from the Herschel Basin of the Canadian Beaufort Sea shelf document relationships with environmental parameters such as salinity, sea-ice cover, and turbulence. Cores YC18-HB-GC01 and PG2303-1 were collected at 18 and 32 m water depth, respectively. At these sites, sediment accumulation rates range between 0.6 and 1.7 cm yr–1 allowing a near-annual temporal resolution over the last 50 years. Multivariate analyses indicate that benthic foraminiferal assemblages respond primarily to food supply. Dissimilarities between the microfaunal assemblages of the two cores are mainly the result of bottom water salinity levels linked to water depth. High abundance of the benthic foraminiferal species Elphidium clavatum and occurrences of Elphidium bartletti point to varying, but relatively low, salinities at the shallow core site YC18-HB-GC01, which may be affected by variations in the summer halocline depth. Higher species diversity and more abundant Cassidulina reniforme and Stainforthia feylingi characterize the deeper core PG2303-1, which might reflect more stable conditions and higher bottom-water salinities throughout the studied time interval. The most important microfaunal shift of the last 50 years, observed in the shallower longer core YC18-HB-GC01, occurred at the turn of the 21st century. Prior to ∼2000 CE, the presence of Islandiella norcrossi indicates more stable and saline conditions. Since ∼2000 CE, increased abundances of Haynesina nivea and of the ciliate Tintinnopsis fimbriata suggest decreased salinity and increased turbidity. An increased abundance of Eoeponidella pulchella after ∼2000 CE suggests a concurrent increase in productivity in the last two decades. This shift is nearly synchronous with a decrease in mean summer sea-ice concentration, which can play an important role in bottom water stability on the shelf. Easterly winds can induce a reduction in the sea-ice cover, but also foster a westward spreading of the Mackenzie River plume and the upwelling of nutrient-rich Pacific waters onto the shelf. Both factors would explain the increased freshening and productivity of the Herschel Basin. The last two decades were also marked by a decrease in ostracod abundance that may relate to higher water turbidity. This study shows that combining information from benthic foraminifera, ostracods, and tintinnids provides a comprehensive insight into recent hydrographic/climatic changes in nearshore Arctic habitats, where productivity is critical for the food security of local communities.

Paleothermometric Inferences Using Elemental Mapping: An Appraisal of Ostracoda Species from Shallow Core Sediment of Bay of Bengal, India

Elemental mapping of Ostracoda valves to infer the paleothermometric fluctuations from off-Visakhapatnam, Bay of Bengal is the focus of the present paper. Two Ostracoda species such as Bairdoppilata (Bairdoppilata) alcyonicola and Actinocythereis scutigera were dominant throughout the core. The weight percentages of Mg, Sr, Ba, Fe and Mn in ostracod carapaces were estimated and the ratio with respect to Ca was analyzed to decipher the paleoclimate and oxy-redox conditions. It is observed that trace element accumulation varies with respect to different ostracod species and shell position. Thus, the elemental mapping of the ostracod carapaces by the EDS techniques proved as a primary proxy to decipher the paleotemperature fluctuations in the study area. Keywords: Paleothermometry, Elemental Mapping, Ostracoda, Visakhapatnam, Bay of Bengal

Lagena Inferocostata – An Opportunistic Index Species of Hydrocarbon Pollutants in Gokana and Khana Local Government Areas of River State, Nigeria

Thirty shallow cores ranging from 0.0 – 1.0m and composited into upper and lower section, was studied to determine relationship between petroleum hydrocarbon pollutants and benthic foraminiferal species in order to determine an index for environmental monitoring in the study area. The core samples were taken from Gokana and Khana Local Government Areas of Rivers State, Nigeria. The cores were treated with Rose Bengal solution to distinguish between then acoenose and biocoenosis benthonic foraminiferal. The cores were analyzed for sedimentology, petroleum hydrocarbon concentrations and micropaleontology. Sedimentological descriptions of the cores indicated inhibited free flow of fluids based on the rock textures which mostly consists of clays and peaty clays. Petroleum hydrocarbons analysis consisted of Total Petroleum Hydrocarbon (TPH) and Polycyclic Aromatic Hydrocarbon (PAH). The results in concentrations of the petroleum hydrocarbons indicated higher concentrations in the upper composited sections. Total Petroleum Hydrocarbon (TPH) and Polycyclic Aromatic Hydrocarbon (PAH) ranged between (20612mg/kg – 37.09mg/kg and 9.55mg/kg – <0.01mg/kg respectively) in the upper composited section and between (14731mg/kg – 14.67mg/kg and 5.71mg/kg – <0.01mg/kg respectively) in the lower composited section. Occurrence of benthonic foraminiferal in the study samples is very low to barren among very low diversity. Three benthonic foraminiferal diversities were identified in the study samples and included Lagenainferocostata, Anomalinoidesspp and Ammodiscus spp. Benthonic foraminifera showed negative correlation with TPH/PAH in the presence of carcinogenic Polycyclic Aromatic Hydrocarbon (PAH) congeners except for Lagenainferocostata which increased consistently and significantly with increasing TPH/PAH in the two composited sections. This suggested Lagenainferocostata as an Index opportunistic species for the hydrocarbon polluted study area.

Prefacing unexplored archives from Central Andean surface-to-bedrock ice cores through a multifaceted investigation of regional firn and ice core glaciochemistry

AbstractShallow firn cores, in addition to a near-basal ice core, were recovered in 2018 from the Quelccaya ice cap (5470 m a.s.l) in the Cordillera Vilcanota, Peru, and in 2017 from the Nevado Illimani glacier (6350 m a.s.l) in the Cordillera Real, Bolivia. The two sites are ~450 km apart. Despite meltwater percolation resulting from warming, particle-based trace element records (e.g. Fe, Mg, K) in the Quelccaya and Illimani shallow cores retain well-preserved signals. The firn core chronologies, established independently by annual layer counting, show a convincing overlap indicating the two records contain comparable signals and therefore capture similar regional scale climatology. Trace element records at a ~1–4 cm resolution provide past records of anthropogenic emissions, dust sources, volcanic emissions, evaporite salts and marine-sourced air masses. Using novel ultra-high-resolution (120 μm) laser technology, we identify annual layer thicknesses ranging from 0.3 to 0.8 cm in a section of 2000-year-old radiocarbon-dated near-basal ice which compared to the previous annual layer estimates suggests that Quelccaya ice cores drilled to bedrock may be older than previously suggested by depth-age models. With the information collected from this study in combination with past studies, we emphasize the importance of collecting new surface-to-bedrock ice cores from at least the Quelccaya ice cap, in particular, due to its projected disappearance as soon as the 2050s.