Search

Three seasonally resolved ice core records covering the 20th century were extracted in 1994, 2004, and 2012 at a nearly identical location from the Col du Dôme (4250 m above sea level, m a.s.l.; Mont Blanc, French Alps) drill site. Here, we complete and combine chemical records of major ions and radiometric measurements of 3H and 210Pb obtained from these three cores with a 3D ice flow model of the Col du Dôme glacier to investigate in detail the origin of discontinuities observed in the depth–age relation of the ice cores drilled in 2004 and 2012. Taking advantage of the granitic bedrock at Col du Dôme, which makes the ice core 210Pb records sensitive to the presence of upstream crevasses, and the fact that the depth–age disturbances are observed at depths for which absolute time markers are available, we draw an overall picture of a dynamic crevasse formation. This can explain the non-disturbed depth–age relation of the ice core drilled in 1994 and the perturbations observed in those drilled in 2004 and 2012. Since crevasses are common at high-Alpine glacier sites, our study points to the important need for rigorous investigations of the depth–age scale and glaciological conditions upstream of drill sites before interpreting high-alpine ice core records in terms of atmospheric changes.

Carbon monoxide (CO) is a naturally occurring atmospheric trace gas, a regulated pollutant, and one of the main components determining the oxidative capacity of the atmosphere. Evaluating climate–chemistry models under different conditions than today and constraining past CO sources requires a reliable record of atmospheric CO mixing ratios ([CO]) that includes data since preindustrial times. Here, we report the first continuous record of atmospheric [CO] for Southern Hemisphere (SH) high latitudes over the past 3 millennia. Our continuous record is a composite of three high-resolution Antarctic ice core gas records and firn air measurements from seven Antarctic locations. The ice core gas [CO] records were measured by continuous flow analysis (CFA), using an optical feedback cavity-enhanced absorption spectrometer (OF-CEAS), achieving excellent external precision (2.8–8.8 ppb; 2σ) and consistently low blanks (ranging from 4.1±1.2 to 7.4±1.4 ppb), thus enabling paleo-atmospheric interpretations. Six new firn air [CO] Antarctic datasets collected between 1993 and 2016 CE at the DE08-2, DSSW19K, DSSW20K, South Pole, Aurora Basin North (ABN), and Lock-In sites (and one previously published firn CO dataset at Berkner) were used to reconstruct the atmospheric history of CO from ∼1897 CE, using inverse modeling that incorporates the influence of gas transport in firn. Excellent consistency was observed between the youngest ice core gas [CO] and the [CO] from the base of the firn and between the recent firn [CO] and atmospheric [CO] measurements at Mawson station (eastern Antarctica), yielding a consistent and contiguous record of CO across these different archives. Our Antarctic [CO] record is relatively stable from −835 to 1500 CE, with mixing ratios within a 30–45 ppb range (2σ). There is a ∼5 ppb decrease in [CO] to a minimum at around 1700 CE during the Little Ice Age. CO mixing ratios then increase over time to reach a maximum of ∼54 ppb by ∼1985 CE. Most of the industrial period [CO] growth occurred between about 1940 to 1985 CE, after which there was an overall [CO] decrease, as observed in Greenland firn air and later at atmospheric monitoring sites and attributed partly to reduced CO emissions from combustion sources. Our Antarctic ice core gas CO observations differ from previously published records in two key aspects. First, our mixing ratios are significantly lower than reported previously, suggesting that previous studies underestimated blank contributions. Second, our new CO record does not show a maximum in the late 1800s. The absence of a [CO] peak around the turn of the century argues against there being a peak in Southern Hemisphere biomass burning at this time, which is in agreement with (i) other paleofire proxies such as ethane or acetylene and (ii) conclusions reached by paleofire modeling. The combined ice core and firn air [CO] history, spanning −835 to 1992 CE, extended to the present by the Mawson atmospheric record, provides a useful benchmark for future atmospheric chemistry modeling studies.

Changes in sea ice conditions and atmospheric circulation over the Southern Ocean play an important role in modulating Antarctic climate. However, observations of both sea ice and wind conditions are limited in Antarctica and the Southern Ocean, both temporally and spatially, prior to the satellite era (1970 onwards). Ice core chemistry data can be used to reconstruct changes over annual, decadal, and millennial timescales. To facilitate sea ice and wind reconstructions, the CLIVASH2k (CLimate Variability in Antarctica and the Southern Hemisphere over the past 2000 years) working group has compiled a database of two species, sodium [Na+] and sulfate [SO42-], commonly measured ionic species. The database (https://doi.org/10.5285/9E0ED16E-F2AB-4372-8DF3-FDE7E388C9A7; Thomas et al., 2022) comprises records from 105 Antarctic ice cores, containing records with a maximum age duration of 2000 years. An initial filter has been applied, based on evaluation against sea ice concentration, geopotential height (500 hPa), and surface wind fields to identify sites suitable for reconstructing past sea ice conditions, wind strength, or atmospheric circulation.

The temperature of the Earth is one of the most important climate parameters. Proxy records of past climate changes, in particular temperature, represent a fundamental tool for exploring internal climate processes and natural climate forcings. Despite the excellent information provided by ice core records in Antarctica, the temperature variability of the past 2000 years is difficult to evaluate from the low-accumulation sites in the Antarctic continent interior. Here we present the results from the Aurora Basin North (ABN) ice core (71∘ S, 111∘ E, 2690 m a.s.l.) in the lower part of the East Antarctic plateau, where accumulation is substantially higher than other ice core drilling sites on the plateau, and provide unprecedented insight into East Antarctic past temperature variability. We reconstructed the temperature of the last 2000 years using two independent methods: the widely used water stable isotopes (δ18O) and by inverse modelling of borehole temperature and past temperature gradients estimated from the inert gas stable isotopes (δ40Ar and δ15N). This second reconstruction is based on three independent measurement types: borehole temperature, firn thickness, and firn temperature gradient. The δ18O temperature reconstruction supports stable temperature conditions within 1 ∘C over the past 2000 years, in agreement with other ice core δ18O records in the region. However, the gas and borehole temperature reconstruction suggests that surface conditions 2 ∘C cooler than average prevailed in the 1000–1400 CE period and supports a 20th century warming of 1 ∘C. A precipitation hiatus during cold periods could explain why water isotope temperature reconstruction underestimates the temperature changes. Both reconstructions arguably record climate in their own way, with a focus on atmospheric and hydrologic cycles for water isotopes, as opposed to surface temperature for gas isotopes and boreholes. This study demonstrates the importance of using a variety of sources for comprehensive paleoclimate reconstructions.

Cosmic rays entering the Earth's atmosphere produce showers of secondary particles such as protons, neutrons, and muons. The interaction of these particles with oxygen-16 (16O) in minerals such as ice and quartz can produce carbon-14 (14C). In glacial ice, 14C is also incorporated through trapping of 14C-containing atmospheric gases (14CO2, 14CO, and 14CH4). Understanding the production rates of in situ cosmogenic 14C is important to deconvolve the in situ cosmogenic and atmospheric 14C signals in ice, both of which contain valuable paleoenvironmental information. Unfortunately, the in situ 14C production rates by muons (which are the dominant production mechanism at depths of >6 m solid ice equivalent) are uncertain. In this study, we use measurements of in situ 14C in ancient ice (>50 ka) from the Taylor Glacier, an ablation site in Antarctica, in combination with a 2D ice flow model to better constrain the compound-specific rates of 14C production by muons and the partitioning of in situ 14C between CO2, CO, and CH4. Our measurements show that 33.7 % (±11.4 %; 95 % confidence interval) of the produced cosmogenic 14C forms 14CO and 66.1 % (±11.5 %; 95 % confidence interval) of the produced cosmogenic 14C forms 14CO2. 14CH4 represents a very small fraction (%) of the total. Assuming that the majority of in situ muogenic 14C in ice forms 14CO2, 14CO, and 14CH4, we also calculated muogenic 14C production rates that are lower by factors of 5.7 (3.6–13.9; 95 % confidence interval) and 3.7 (2.0–11.9; 95 % confidence interval) for negative muon capture and fast muon interactions, respectively, when compared to values determined in quartz from laboratory studies (Heisinger et al., 2002a, b) and in a natural setting (Lupker et al., 2015). This apparent discrepancy in muogenic 14C production rates in ice and quartz currently lacks a good explanation and requires further investigation.

The injection of sulfur into the stratosphere by volcanic eruptions is the dominant driver of natural climate variability on interannual to multidecadal timescales. Based on a set of continuous sulfate and sulfur records from a suite of ice cores from Greenland and Antarctica, the HolVol v.1.0 database includes estimates of the magnitudes and approximate source latitudes of major volcanic stratospheric sulfur injection (VSSI) events for the Holocene (from 9500 BCE or 11 500 years BP to 1900 CE), constituting an extension of the previous record by 7000 years. The database incorporates new-generation ice-core aerosol records with a sub-annual temporal resolution and a demonstrated sub-decadal dating accuracy and precision. By tightly aligning and stacking the ice-core records on the WD2014 chronology from Antarctica, we resolve long-standing inconsistencies in the dating of ancient volcanic eruptions that arise from biased (i.e., dated too old) ice-core chronologies over the Holocene for Greenland. We reconstruct a total of 850 volcanic eruptions with injections in excess of 1 teragram of sulfur (Tg S); of these eruptions, 329 (39 %) are located in the low latitudes with bipolar sulfate deposition, 426 (50 %) are located in the Northern Hemisphere extratropics (NHET) and 88 (10 %) are located in the Southern Hemisphere extratropics (SHET). The spatial distribution of the reconstructed eruption locations is in agreement with prior reconstructions for the past 2500 years. In total, these eruptions injected 7410 Tg S into the stratosphere: 70 % from tropical eruptions and 25 % from NH extratropical eruptions. A long-term latitudinally and monthly resolved stratospheric aerosol optical depth (SAOD) time series is reconstructed from the HolVol VSSI estimates, representing the first Holocene-scale reconstruction constrained by Greenland and Antarctica ice cores. These new long-term reconstructions of past VSSI and SAOD variability confirm evidence from regional volcanic eruption chronologies (e.g., from Iceland) in showing that the Early Holocene (9500–7000 BCE) experienced a higher number of volcanic eruptions (+16 %) and cumulative VSSI (+86 %) compared with the past 2500 years. This increase coincides with the rapid retreat of ice sheets during deglaciation, providing context for potential future increases in volcanic activity in regions under projected glacier melting in the 21st century. The reconstructed VSSI and SAOD data are available at https://doi.org/10.1594/PANGAEA.928646 (Sigl et al., 2021).

Carbon monoxide (CO) is a regulated pollutant and one of the key components determining the oxidizing capacity of the atmosphere. Obtaining a reliable record of atmospheric CO mixing ratios ([CO]) since preindustrial times is necessary to evaluate climate–chemistry models under conditions different from today and to constrain past CO sources. We present high-resolution measurements of CO mixing ratios from ice cores drilled at five different sites on the Greenland ice sheet that experience a range of snow accumulation rates, mean surface temperatures, and different chemical compositions. An optical-feedback cavity-enhanced absorption spectrometer (OF-CEAS) was coupled with continuous melter systems and operated during four analytical campaigns conducted between 2013 and 2019. Overall, continuous flow analysis (CFA) of CO was carried out on over 700 m of ice. The CFA-based CO measurements exhibit excellent external precision (ranging from 3.3 to 6.6 ppbv, 1σ) and achieve consistently low blanks (ranging from 4.1±1.2 to 12.6±4.4 ppbv), enabling paleoatmospheric interpretations. However, the five CO records all exhibit variability that is too large and rapid to reflect past atmospheric mixing ratio changes. Complementary tests conducted on discrete ice samples demonstrate that these variations are not artifacts of the analytical method (i.e., production of CO from organics in the ice during melting) but are very likely related to in situ CO production within the ice before analysis. Evaluation of the signal resolution and co-investigation of high-resolution records of CO and total organic carbon (TOC) suggest that past atmospheric CO variations can be extracted from the records' baselines with accumulation rates higher than 20 cm w.e.yr-1 (water equivalent per year). Consistent baseline CO records from four Greenland sites are combined to produce a multisite average ice core reconstruction of past atmospheric CO for the Northern Hemisphere high latitudes, covering the period from 1700 to 1957 CE. Such a reconstruction should be taken as an upper bound of past atmospheric CO abundance. From 1700 to 1875 CE, the record reveals stable or slightly increasing values in the 100–115 ppbv range. From 1875 to 1957 CE, the record indicates a monotonic increase from 114±4 to 147±6 ppbv. The ice core multisite CO record exhibits an excellent overlap with the atmospheric CO record from Greenland firn air which spans the 1950–2010 CE time period. The combined ice core and firn air CO history, spanning 1700–2010 CE, provides useful constraints for future model studies of atmospheric changes since the preindustrial period.

Volcanic fallout in polar ice sheets provides important opportunities to date and correlate ice-core records as well as to investigate the environmental impacts of eruptions. Only the geochemical characterization of volcanic ash (tephra) embedded in the ice strata can confirm the source of the eruption, however, and is a requisite if historical eruption ages are to be used as valid chronological checks on annual ice layer counting. Here we report the investigation of ash particles in a Greenland ice core that are associated with a volcanic sulfuric acid layer previously attributed to the 79 CE eruption of Vesuvius. Major and trace element composition of the particles indicates that the tephra does not derive from Vesuvius but most likely originates from an unidentified eruption in the Aleutian arc. Using ash dispersal modeling, we find that only an eruption large enough to include stratospheric injection is likely to account for the sizable (24–85 µm) ash particles observed in the Greenland ice at this time. Despite its likely explosivity, this event does not appear to have triggered significant climate perturbations, unlike some other large extratropical eruptions. In light of a recent re-evaluation of the Greenland ice-core chronologies, our findings further challenge the previous assignation of this volcanic event to 79 CE. We highlight the need for the revised Common Era ice-core chronology to be formally accepted by the wider ice-core and climate modeling communities in order to ensure robust age linkages to precisely dated historical and paleoclimate proxy records.

Volcanic eruptions are a key source of climatic variability, and reconstructing their past impact can improve our understanding of the operation of the climate system and increase the accuracy of future climate projections. Two annually resolved and independently dated palaeoarchives – tree rings and polar ice cores – can be used in tandem to assess the timing, strength and climatic impact of volcanic eruptions over the past ∼ 2500 years. The quantification of post-volcanic climate responses, however, has at times been hampered by differences between simulated and observed temperature responses that raised questions regarding the robustness of the chronologies of both archives. While many chronological mismatches have been resolved, the precise timing and climatic impact of two major sulfate-emitting volcanic eruptions during the 1450s CE, including the largest atmospheric sulfate-loading event in the last 700 years, have not been constrained. Here we explore this issue through a combination of tephrochronological evidence and high-resolution ice-core chemistry measurements from a Greenland ice core, the TUNU2013 record. We identify tephra from the historically dated 1477 CE eruption of the Icelandic Veiðivötn–Bárðarbunga volcanic system in direct association with a notable sulfate peak in TUNU2013 attributed to this event, confirming that this peak can be used as a reliable and precise time marker. Using seasonal cycles in several chemical elements and 1477 CE as a fixed chronological point shows that ages of 1453 CE and 1458 CE can be attributed, with high precision, to the start of two other notable sulfate peaks. This confirms the accuracy of a recent Greenland ice-core chronology over the middle to late 15th century and corroborates the findings of recent volcanic reconstructions from Greenland and Antarctica. Overall, this implies that large-scale Northern Hemisphere climatic cooling affecting tree-ring growth in 1453 CE was caused by a Northern Hemisphere volcanic eruption in 1452 or early 1453 CE, and then a Southern Hemisphere eruption, previously assumed to have triggered the cooling, occurred later in 1457 or 1458 CE. The direct attribution of the 1477 CE sulfate peak to the eruption of Veiðivötn, one of the most explosive from Iceland in the last 1200 years, also provides the opportunity to assess the eruption's climatic impact. A tree-ring-based reconstruction of Northern Hemisphere summer temperatures shows a cooling in the aftermath of the eruption of −0.35 ∘C relative to a 1961–1990 CE reference period and −0.1 ∘C relative to the 30-year period around the event, as well as a relatively weak and spatially incoherent climatic response in comparison to the less explosive but longer-lasting Icelandic Eldgjá 939 CE and Laki 1783 CE eruptions. In addition, the Veiðivötn 1477 CE eruption occurred around the inception of the Little Ice Age and could be used as a chronostratigraphic marker to constrain the phasing and spatial variability of climate changes over this transition if it can be traced in more regional palaeoclimatic archives.

Water levels in streams and aquifers often exhibit daily cycles during rainless periods, reflecting daytime extraction of shallow groundwater by evapotranspiration (ET) and, during snowmelt, daytime additions of meltwater. These cycles can aid in understanding the mechanisms that couple solar forcing of ET and snowmelt to changes in streamflow. Here we analyze 3 years of 30 min solar flux, sap flow, stream stage, and groundwater level measurements at Sagehen Creek and Independence Creek, two snow-dominated headwater catchments in California's Sierra Nevada mountains. Despite their sharply contrasting geological settings (most of the Independence basin is glacially scoured granodiorite, whereas Sagehen is underlain by hundreds of meters of volcanic and volcaniclastic deposits that host an extensive groundwater aquifer), both streams respond similarly to snowmelt and ET forcing. During snow-free summer periods, daily cycles in solar flux are tightly correlated with variations in sap flow, and with the rates of water level rise and fall in streams and riparian aquifers. During these periods, stream stages and riparian groundwater levels decline during the day and rebound at night. These cycles are reversed during snowmelt, with stream stages and riparian groundwater levels rising during the day in response to snowmelt inputs and falling at night as the riparian aquifer drains. Streamflow and groundwater maxima and minima (during snowmelt- and ET-dominated periods, respectively) lag the midday peak in solar flux by several hours. A simple conceptual model explains this lag: streamflows depend on riparian aquifer water levels, which integrate snowmelt inputs and ET losses over time, and thus will be phase-shifted relative to the peaks in snowmelt and evapotranspiration rates. Thus, although the lag between solar forcing and water level cycles is often interpreted as a travel-time lag, our analysis shows that it is mostly a dynamical phase lag, at least in small catchments. Furthermore, although daily cycles in streamflow have often been used to estimate ET fluxes, our simple conceptual model demonstrates that this is infeasible unless the response time of the riparian aquifer can be determined. As the snowmelt season progresses, snowmelt forcing of groundwater and streamflow weakens and evapotranspiration forcing strengthens. The relative dominance of snowmelt vs. ET can be quantified by the diel cycle index, which measures the correlation between the solar flux and the rate of rise or fall in streamflow or groundwater. When the snowpack melts out at an individual location, the local groundwater shifts abruptly from snowmelt-dominated cycles to ET-dominated cycles. Melt-out and the corresponding shift in the diel cycle index occur earlier at lower altitudes and on south-facing slopes, and streamflow integrates these transitions over the drainage network. Thus the diel cycle index in streamflow shifts gradually, beginning when the snowpack melts out near the gauging station and ending, months later, when the snowpack melts out at the top of the basin and the entire drainage network becomes dominated by ET cycles. During this long transition, snowmelt signals generated in the upper basin are gradually overprinted by ET signals generated lower down in the basin. The gradual springtime transition in the diel cycle index is mirrored in sequences of Landsat images showing the springtime retreat of the snowpack to higher elevations and the corresponding advance of photosynthetic activity across the basin. Trends in the catchment-averaged MODIS enhanced vegetation index (EVI2) also correlate closely with the late springtime shift from snowmelt to ET cycles and with the autumn shift back toward snowmelt cycles. Seasonal changes in streamflow cycles therefore reflect catchment-scale shifts in snowpack and vegetation activity that can be seen from Earth orbit. The data and analyses presented here illustrate how streams can act as mirrors of the landscape, integrating physical and ecohydrological signals across their contributing drainage networks.

The last glacial period is characterized by a number of millennial climate events that have been identified in both Greenland and Antarctic ice cores and that are abrupt in Greenland climate records. The mechanisms governing this climate variability remain a puzzle that requires a precise synchronization of ice cores from the two hemispheres to be resolved. Previously, Greenland and Antarctic ice cores have been synchronized primarily via their common records of gas concentrations or isotopes from the trapped air and via cosmogenic isotopes measured on the ice. In this work, we apply ice core volcanic proxies and annual layer counting to identify large volcanic eruptions that have left a signature in both Greenland and Antarctica. Generally, no tephra is associated with those eruptions in the ice cores, so the source of the eruptions cannot be identified. Instead, we identify and match sequences of volcanic eruptions with bipolar distribution of sulfate, i.e. unique patterns of volcanic events separated by the same number of years at the two poles. Using this approach, we pinpoint 82 large bipolar volcanic eruptions throughout the second half of the last glacial period (12–60 ka). This improved ice core synchronization is applied to determine the bipolar phasing of abrupt climate change events at decadal-scale precision. In response to Greenland abrupt climatic transitions, we find a response in the Antarctic water isotope signals (δ18O and deuterium excess) that is both more immediate and more abrupt than that found with previous gas-based interpolar synchronizations, providing additional support for our volcanic framework. On average, the Antarctic bipolar seesaw climate response lags the midpoint of Greenland abrupt δ18O transitions by 122±24 years. The time difference between Antarctic signals in deuterium excess and δ18O, which likewise informs the time needed to propagate the signal as described by the theory of the bipolar seesaw but is less sensitive to synchronization errors, suggests an Antarctic δ18O lag behind Greenland of 152±37 years. These estimates are shorter than the 200 years suggested by earlier gas-based synchronizations. As before, we find variations in the timing and duration between the response at different sites and for different events suggesting an interaction of oceanic and atmospheric teleconnection patterns as well as internal climate variability.

Antarctica is well known to be highly susceptible to atmospheric and oceanic warming. However, due to the lack of long-term and in situ meteorological observations, little is known about the magnitude of the warming and the meteorological conditions in the intersection region between the Antarctic Peninsula (AP), the West Antarctic Ice Sheet (WAIS) and the East Antarctic Ice Sheet (EAIS). Here we present new stable water isotope data (δ18O, δD, d excess) and accumulation rates from firn cores in the Union Glacier (UG) region, located in the Ellsworth Mountains at the northern edge of the WAIS. The firn core stable oxygen isotopes and the d excess exhibit no statistically significant trend for the period 1980–2014, suggesting that regional changes in near-surface air temperature and moisture source variability have been small during the last 35 years. Backward trajectory modelling revealed the Weddell Sea sector, Coats Land and Dronning Maud Land (DML) to be the main moisture source regions for the study site throughout the year. We found that mean annual δ18O (δD) values in the UG region are negatively correlated with sea ice concentrations (SICs) in the northern Weddell Sea but not influenced by large-scale modes of climate variability such as the Southern Annular Mode (SAM) and the El Niño–Southern Oscillation (ENSO). Only mean annual d-excess values show a weak positive correlation with the SAM. On average annual snow accumulation in the UG region amounts to 0.245 m w.e. a−1 in 1980–2014 and has slightly decreased during this period. It is only weakly related to sea ice conditions in the Weddell Sea sector and not correlated with SAM and ENSO. We conclude that neither the rapid warming nor the large increases in snow accumulation observed on the AP and in West Antarctica during the last decades have extended inland to the Ellsworth Mountains. Hence, the UG region, although located at the northern edge of the WAIS and relatively close to the AP, exhibits rather stable climate characteristics similar to those observed in East Antarctica.

The South Pole Ice Core (SPICEcore) was drilled in 2014–2016 to provide a detailed multi-proxy archive of paleoclimate conditions in East Antarctica during the Holocene and late Pleistocene. Interpretation of these records requires an accurate depth–age relationship. Here, we present the SPICEcore (SP19) timescale for the age of the ice of SPICEcore. SP19 is synchronized to the WD2014 chronology from the West Antarctic Ice Sheet Divide (WAIS Divide) ice core using stratigraphic matching of 251 volcanic events. These events indicate an age of 54 302±519 BP (years before 1950) at the bottom of SPICEcore. Annual layers identified in sodium and magnesium ions to 11 341 BP were used to interpolate between stratigraphic volcanic tie points, yielding an annually resolved chronology through the Holocene. Estimated timescale uncertainty during the Holocene is less than 18 years relative to WD2014, with the exception of the interval between 1800 to 3100 BP when uncertainty estimates reach ±25 years due to widely spaced volcanic tie points. Prior to the Holocene, uncertainties remain within 124 years relative to WD2014. Results show an average Holocene accumulation rate of 7.4 cm yr−1 (water equivalent). The time variability of accumulation rate is consistent with expectations for steady-state ice flow through the modern spatial pattern of accumulation rate. Time variations in nitrate concentration, nitrate seasonal amplitude and δ15N of N2 in turn are as expected for the accumulation rate variations. The highly variable yet well-constrained Holocene accumulation history at the site can help improve scientific understanding of deposition-sensitive climate proxies such as δ15N of N2 and photolyzed chemical compounds.

New ice cores retrieved from the Taylor Glacier (Antarctica) blue ice area contain ice and air spanning the Marine Isotope Stage (MIS) 5–4 transition, a period of global cooling and ice sheet expansion. We determine chronologies for the ice and air bubbles in the new ice cores by visually matching variations in gas- and ice-phase tracers to preexisting ice core records. The chronologies reveal an ice age–gas age difference (Δage) approaching 10 ka during MIS 4, implying very low snow accumulation in the Taylor Glacier accumulation zone. A revised chronology for the analogous section of the Taylor Dome ice core (84 to 55 ka), located to the south of the Taylor Glacier accumulation zone, shows that Δage did not exceed 3 ka. The difference in Δage between the two records during MIS 4 is similar in magnitude but opposite in direction to what is observed at the Last Glacial Maximum. This relationship implies that a spatial gradient in snow accumulation existed across the Taylor Dome region during MIS 4 that was oriented in the opposite direction of the accumulation gradient during the Last Glacial Maximum.

Abstract Changes in the composition of dust trapped in ice provide evidence of past atmospheric circulation and earth surface conditions. Investigations of dust provenance in Antarctic ice during glacial and interglacial periods indicate that South America is the primary dust source during both climate regimes. Here, we present results from a new ice core dust archive extracted from the Taylor Glacier in coastal East Antarctica during the deglacial transition from Marine Isotope Stage 6 to 5e. Radiogenic strontium and neodymium isotopes indicate that last interglacial dust is young and volcanic, in contrast to the observed preindustrial and Holocene (Marine Isotope Stage 1) dust composition. The dust composition differences from the last interglacial and current interglacial period at the site require a profound difference in atmospheric transport and environmental conditions. We consider several potential causes for enhanced transport of volcanic material to the site, including increased availability of volcanic material and large‐scale atmospheric circulation changes.

In this study, vanillic acid was measured in the Tunu ice core from northeastern Greenland in samples covering the past 1700 years. Vanillic acid is an aerosol-borne aromatic methoxy acid, produced by the combustion of lignin during biomass burning. Air mass trajectory analysis indicates that North American boreal forests are likely the major source region for biomass burning aerosols deposited to the ice core site. Vanillic acid levels in the Tunu ice core range from

Long-range transport of aerosol from lower latitudes to the high Arctic may be a significant contributor to climate forcing in the Arctic. To identify the sources of key contaminants entering the Canadian High Arctic an intensive campaign of snow sampling was completed at Alert, Nunavut, from September 2014 to June 2015. Fresh snow samples collected every few days were analyzed for black carbon, major ions, and metals, and this rich data set provided an opportunity for a temporally refined source apportionment of snow composition via positive matrix factorization (PMF) in conjunction with FLEXPART (FLEXible PARTicle dispersion model) potential emission sensitivity analysis. Seven source factors were identified: sea salt, crustal metals, black carbon, carboxylic acids, nitrate, non-crustal metals, and sulfate. The sea salt and crustal factors showed good agreement with expected composition and primarily northern sources. High loadings of V and Se onto Factor 2, crustal metals, was consistent with expected elemental ratios, implying these metals were not primarily anthropogenic in origin. Factor 3, black carbon, was an acidic factor dominated by black carbon but with some sulfate contribution over the winter-haze season. The lack of K+ associated with this factor, a Eurasian source, and limited known forest fire events coincident with this factor's peak suggested a predominantly anthropogenic combustion source. Factor 4, carboxylic acids, was dominated by formate and acetate with a moderate correlation to available sunlight and an oceanic and North American source. A robust identification of this factor was not possible; however, atmospheric photochemical reactions, ocean microlayer reaction, and biomass burning were explored as potential contributors. Factor 5, nitrate, was an acidic factor dominated by NO3−, with a likely Eurasian source and mid-winter peak. The isolation of NO3− on a separate factor may reflect its complex atmospheric processing, though the associated source region suggests possibly anthropogenic precursors. Factor 6, non-crustal metals, showed heightened loadings of Sb, Pb, and As, and correlation with other metals traditionally associated with industrial activities. Similar to Factor 3 and 5, this factor appeared to be largely Eurasian in origin. Factor 7, sulfate, was dominated by SO42− and MS with a fall peak and high acidity. Coincident volcanic activity and northern source regions may suggest a processed SO2 source of this factor.

Growing evidence suggests that the sea ice surface is an important source of sea salt aerosol and this has significant implications for polar climate and atmospheric chemistry. It also suggests the potential to use ice core sea salt records as proxies for past sea ice extent. To explore this possibility in the Arctic region, we use a chemical transport model to track the emission, transport, and deposition of sea salt from both the open ocean and the sea ice, allowing us to assess the relative importance of each. Our results confirm the importance of sea ice sea salt (SISS) to the winter Arctic aerosol burden. For the first time, we explicitly simulate the sea salt concentrations of Greenland snow, achieving values within a factor of two of Greenland ice core records. Our simulations suggest that SISS contributes to the winter maxima in sea salt characteristic of ice cores across Greenland. However, a north–south gradient in the contribution of SISS relative to open-ocean sea salt (OOSS) exists across Greenland, with 50 % of winter sea salt being SISS at northern sites such as NEEM (77° N), while only 10 % of winter sea salt is SISS at southern locations such as ACT10C (66° N). Our model shows some skill at reproducing the inter-annual variability in sea salt concentrations for 1991–1999, particularly at Summit where up to 62 % of the variability is explained. Future work will involve constraining what is driving this inter-annual variability and operating the model under different palaeoclimatic conditions.

Englacial ice contains a significant reservoir of organic material (OM), preserving a chronological record of materials from Earth's past. Here, we investigate if OM composition surveys in ice core research can provide paleoecological information on the dynamic nature of our Earth through time. Temporal trends in OM composition from the early Holocene extending back to the Last Glacial Maximum (LGM) of the West Antarctic Ice Sheet Divide (WD) ice core were measured by fluorescence spectroscopy. Multivariate parallel factor (PARAFAC) analysis is widely used to isolate the chemical components that best describe the observed variation across three-dimensional fluorescence spectroscopy (excitation–emission matrices; EEMs) assays. Fluorescent OM markers identified by PARAFAC modeling of the EEMs from the LGM (27.0–18.0 kyr BP; before present 1950) through the last deglaciation (LD; 18.0–11.5 kyr BP), to the mid-Holocene (11.5–6.0 kyr BP) provided evidence of different types of fluorescent OM composition and origin in the WD ice core over 21.0 kyr. Low excitation–emission wavelength fluorescent PARAFAC component one (C1), associated with chemical species similar to simple lignin phenols was the greatest contributor throughout the ice core, suggesting a strong signature of terrestrial OM in all climate periods. The component two (C2) OM marker, encompassed distinct variability in the ice core describing chemical species similar to tannin- and phenylalanine-like material. Component three (C3), associated with humic-like terrestrial material further resistant to biodegradation, was only characteristic of the Holocene, suggesting that more complex organic polymers such as lignins or tannins may be an ecological marker of warmer climates. We suggest that fluorescent OM markers observed during the LGM were the result of greater continental dust loading of lignin precursor (monolignol) material in a drier climate, with lower marine influences when sea ice extent was higher and continents had more expansive tundra cover. As the climate warmed, the record of OM markers in the WD ice core changed, reflecting shifts in carbon productivity as a result of global ecosystem response.

Rapidly rising temperatures and loss of snow and ice cover have demonstrated the unique vulnerability of the high Arctic to climate change. There are major uncertainties in modelling the chemical depositional and scavenging processes of Arctic snow. To that end, fresh snow samples collected on average every 4 days at Alert, Nunavut, from September 2014 to June 2015 were analyzed for black carbon, major ions, and metals, and their concentrations and fluxes were reported. Comparison with simultaneous measurements of atmospheric aerosol mass loadings yields effective deposition velocities that encompass all processes by which the atmospheric species are transferred to the snow. It is inferred from these values that dry deposition is the dominant removal mechanism for several compounds over the winter while wet deposition increased in importance in the fall and spring, possibly due to enhanced scavenging by mixed-phase clouds. Black carbon aerosol was the least efficiently deposited species to the snow.

Wildfires and their emissions have significant impacts on ecosystems, climate, atmospheric chemistry, and carbon cycling. Well-dated proxy records are needed to study the long-term climatic controls on biomass burning and the associated climate feedbacks. There is a particular lack of information about long-term biomass burning variations in Siberia, the largest forested area in the Northern Hemisphere. In this study we report analyses of aromatic acids (vanillic and para-hydroxybenzoic acids) over the past 2600 years in the Eurasian Arctic Akademii Nauk ice core. These compounds are aerosol-borne, semi-volatile organic compounds derived from lignin combustion. The analyses were made using ion chromatography with electrospray mass spectrometric detection. The levels of these aromatic acids ranged from below the detection limit (0.01 to 0.05 ppb; 1 ppb = 1000 ng L−1) to about 1 ppb, with roughly 30 % of the samples above the detection limit. In the preindustrial late Holocene, highly elevated aromatic acid levels are observed during three distinct periods (650–300 BCE, 340–660 CE, and 1460–1660 CE). The timing of the two most recent periods coincides with the episodic pulsing of ice-rafted debris in the North Atlantic known as Bond events and a weakened Asian monsoon, suggesting a link between fires and large-scale climate variability on millennial timescales. Aromatic acid levels also are elevated during the onset of the industrial period from 1780 to 1860 CE, but with a different ratio of vanillic and para-hydroxybenzoic acid than is observed during the preindustrial period. This study provides the first millennial-scale record of aromatic acids. This study clearly demonstrates that coherent aromatic acid signals are recorded in polar ice cores that can be used as proxies for past trends in biomass burning.

Reconstruction of past changes in Arctic sea ice extent may be critical for understanding its future evolution. Methanesulfonate (MSA) and bromine concentrations preserved in ice cores have both been proposed as indicators of past sea ice conditions. In this study, two ice cores from central and north-eastern Greenland were analysed at sub-annual resolution for MSA (CH3SO3H) and bromine, covering the time period 1750–2010. We examine correlations between ice core MSA and the HadISST1 ICE sea ice dataset and consult back trajectories to infer the likely source regions. A strong correlation between the low-frequency MSA and bromine records during pre-industrial times indicates that both chemical species are likely linked to processes occurring on or near sea ice in the same source regions. The positive correlation between ice core MSA and bromine persists until the mid-20th century, when the acidity of Greenland ice begins to increase markedly due to increased fossil fuel emissions. After that time, MSA levels decrease as a result of declining sea ice extent but bromine levels increase. We consider several possible explanations and ultimately suggest that increased acidity, specifically nitric acid, of snow on sea ice stimulates the release of reactive Br from sea ice, resulting in increased transport and deposition on the Greenland ice sheet.

Contemporary climate warming over the Arctic is accelerating mass loss from the Greenland Ice Sheet through increasing surface melt, emphasizing the need to closely monitor its surface mass balance in order to improve sea-level rise predictions. Snow accumulation is the largest component of the ice sheet's surface mass balance, but in situ observations thereof are inherently sparse and models are difficult to evaluate at large scales. Here, we quantify recent Greenland accumulation rates using ultra-wideband (2–6.5 GHz) airborne snow radar data collected as part of NASA's Operation IceBridge between 2009 and 2012. We use a semiautomated method to trace the observed radiostratigraphy and then derive annual net accumulation rates for 2009–2012. The uncertainty in these radar-derived accumulation rates is on average 14 %. A comparison of the radar-derived accumulation rates and contemporaneous ice cores shows that snow radar captures both the annual and long-term mean accumulation rate accurately. A comparison with outputs from a regional climate model (MAR) shows that this model matches radar-derived accumulation rates in the ice sheet interior but produces higher values over southeastern Greenland. Our results demonstrate that snow radar can efficiently and accurately map patterns of snow accumulation across an ice sheet and that it is valuable for evaluating the accuracy of surface mass balance models.

Advances in trace gas analysis allow localised, non-atmospheric features to be resolved in ice cores, superimposed on the coherent atmospheric signal. These high-frequency signals could not have survived the low-pass filter effect that gas diffusion in the firn exerts on the atmospheric history and therefore do not result from changes in the atmospheric composition at the ice sheet surface. Using continuous methane (CH4) records obtained from five polar ice cores, we characterise these non-atmospheric signals and explore their origin. Isolated samples, enriched in CH4 in the Tunu13 (Greenland) record are linked to the presence of melt layers. Melting can enrich the methane concentration due to a solubility effect, but we find that an additional in situ process is required to generate the full magnitude of these anomalies. Furthermore, in all the ice cores studied there is evidence of reproducible, decimetre-scale CH4 variability. Through a series of tests, we demonstrate that this is an artifact of layered bubble trapping in a heterogeneous-density firn column; we use the term “trapping signal” for this phenomenon. The peak-to-peak amplitude of the trapping signal is typically 5 ppb, but may exceed 40 ppb. Signal magnitude increases with atmospheric CH4 growth rate and seasonal density contrast, and decreases with accumulation rate. Significant annual periodicity is present in the CH4 variability of two Greenland ice cores, suggesting that layered gas trapping at these sites is controlled by regular, seasonal variations in the physical properties of the firn. Future analytical campaigns should anticipate high-frequency artifacts at high-melt ice core sites or during time periods with high atmospheric CH4 growth rate in order to avoid misinterpretation of such features as past changes in atmospheric composition.

We present the WD2014 chronology for the upper part (0–2850 m; 31.2 ka BP) of the West Antarctic Ice Sheet (WAIS) Divide (WD) ice core. The chronology is based on counting of annual layers observed in the chemical, dust and electrical conductivity records. These layers are caused by seasonal changes in the source, transport, and deposition of aerosols. The measurements were interpreted manually and with the aid of two automated methods. We validated the chronology by comparing to two high-accuracy, absolutely dated chronologies. For the Holocene, the cosmogenic isotope records of 10Be from WAIS Divide and 14C for IntCal13 demonstrated that WD2014 was consistently accurate to better than 0.5 % of the age. For the glacial period, comparisons to the Hulu Cave chronology demonstrated that WD2014 had an accuracy of better than 1 % of the age at three abrupt climate change events between 27 and 31 ka. WD2014 has consistently younger ages than Greenland ice core chronologies during most of the Holocene. For the Younger Dryas–Preboreal transition (11.595 ka; 24 years younger) and the Bølling–Allerød Warming (14.621 ka; 7 years younger), WD2014 ages are within the combined uncertainties of the timescales. Given its high accuracy, WD2014 can become a reference chronology for the Southern Hemisphere, with synchronization to other chronologies feasible using high-quality proxies of volcanism, solar activity, atmospheric mineral dust, and atmospheric methane concentrations.

The role of sea ice in the Earth climate system is still under debate, although it is known to influence albedo, ocean circulation, and atmosphere–ocean heat and gas exchange. Here we present a reconstruction of 1950 to 1998 AD sea ice in the Laptev Sea based on the Akademii Nauk ice core (Severnaya Zemlya, Russian Arctic). The chemistry of halogens bromine (Br) and iodine (I) is strongly active and influenced by sea ice dynamics, in terms of physical, chemical and biological process. Bromine reacts on the sea ice surface in autocatalyzing "bromine explosion" events, causing an enrichment of the Br / Na ratio and hence a bromine excess (Brexc) in snow compared to that in seawater. Iodine is suggested to be emitted from algal communities growing under sea ice. The results suggest a connection between Brexc and spring sea ice area, as well as a connection between iodine concentration and summer sea ice area. The correlation coefficients obtained between Brexc and spring sea ice (r = 0.44) as well as between iodine and summer sea ice (r = 0.50) for the Laptev Sea suggest that these two halogens could become good candidates for extended reconstructions of past sea ice changes in the Arctic.

Observed changes in the surface elevation of the Greenland Ice Sheet are caused by ice dynamics, basal elevation change, basal melt, surface mass balance (SMB) variability, and by compaction of the overlying firn. The last two contributions are quantified here using a firn model that includes compaction, meltwater percolation, and refreezing. The model is forced with surface mass fluxes and temperature from a regional climate model for the period 1960–2014. The model results agree with observations of surface density, density profiles from 62 firn cores, and altimetric observations from regions where ice-dynamical surface height changes are likely small. In areas with strong surface melt, the firn model overestimates density. We find that the firn layer in the high interior is generally thickening slowly (1–5 cm yr−1). In the percolation and ablation areas, firn and SMB processes account for a surface elevation lowering of up to 20–50 cm yr−1. Most of this firn-induced marginal thinning is caused by an increase in melt since the mid-1990s and partly compensated by an increase in the accumulation of fresh snow around most of the ice sheet. The total firn and ice volume change between 1980 and 2014 is estimated at −3295 ± 1030 km3 due to firn and SMB changes, corresponding to an ice-sheet average thinning of 1.96 ± 0.61 m. Most of this volume decrease occurred after 1995. The computed changes in surface elevation can be used to partition altimetrically observed volume change into surface mass balance and ice-dynamically related mass changes.

The West Antarctic Ice Sheet Divide (WAIS Divide, WD) ice core is a newly drilled, high-accumulation deep ice core that provides Antarctic climate records of the past ∼68 ka at unprecedented temporal resolution. The upper 2850 m (back to 31.2 ka BP) have been dated using annual-layer counting. Here we present a chronology for the deep part of the core (67.8–31.2 ka BP), which is based on stratigraphic matching to annual-layer-counted Greenland ice cores using globally well-mixed atmospheric methane. We calculate the WD gas age–ice age difference (Δage) using a combination of firn densification modeling, ice-flow modeling, and a data set of δ15N-N2, a proxy for past firn column thickness. The largest Δage at WD occurs during the Last Glacial Maximum, and is 525 ± 120 years. Internally consistent solutions can be found only when assuming little to no influence of impurity content on densification rates, contrary to a recently proposed hypothesis. We synchronize the WD chronology to a linearly scaled version of the layer-counted Greenland Ice Core Chronology (GICC05), which brings the age of Dansgaard–Oeschger (DO) events into agreement with the U/Th absolutely dated Hulu Cave speleothem record. The small Δage at WD provides valuable opportunities to investigate the timing of atmospheric greenhouse gas variations relative to Antarctic climate, as well as the interhemispheric phasing of the "bipolar seesaw".

Biomass burning is a major source of greenhouse gases and influences regional to global climate. Pre-industrial fire-history records from black carbon, charcoal and other proxies provide baseline estimates of biomass burning at local to global scales spanning millennia, and are thus useful to examine the role of fire in the carbon cycle and climate system. Here we use the specific biomarker levoglucosan together with black carbon and ammonium concentrations from the North Greenland Eemian (NEEM) ice cores (77.49° N, 51.2° W; 2480 m a.s.l) over the past 2000 years to infer changes in boreal fire activity. Increases in boreal fire activity over the periods 1000–1300 CE and decreases during 700–900 CE coincide with high-latitude NH temperature changes. Levoglucosan concentrations in the NEEM ice cores peak between 1500 and 1700 CE, and most levoglucosan spikes coincide with the most extensive central and northern Asian droughts of the past millennium. Many of these multi-annual droughts are caused by Asian monsoon failures, thus suggesting a connection between low- and high-latitude climate processes. North America is a primary source of biomass burning aerosols due to its relative proximity to the Greenland Ice Cap. During major fire events, however, isotopic analyses of dust, back trajectories and links with levoglucosan peaks and regional drought reconstructions suggest that Siberia is also an important source of pyrogenic aerosols to Greenland.

In Antarctica, uncertainties in mass input and output translate directly into uncertainty in glacier mass balance and thus in sea level impact. While remotely sensed observations of ice velocity and thickness over the major outlet glaciers have improved our understanding of ice loss to the ocean, snow accumulation over the vast Antarctic interior remains largely unmeasured. Here, we show that an airborne radar system, combined with ice-core glaciochemical analysis, provide the means necessary to measure the accumulation rate at the catchment-scale along the Amundsen Sea coast of West Antarctica. We used along-track radar-derived accumulation to generate a 1985–2009 average accumulation grid that resolves moderate- to large-scale features (>25 km) over the Pine Island–Thwaites glacier drainage system. Comparisons with estimates from atmospheric models and gridded climatologies generally show our results as having less accumulation in the lower-elevation coastal zone but greater accumulation in the interior. Ice discharge, measured over discrete time intervals between 1994 and 2012, combined with our catchment-wide accumulation rates provide an 18-year mass balance history for the sector. While Thwaites Glacier lost the most ice in the mid-1990s, Pine Island Glacier's losses increased substantially by 2006, overtaking Thwaites as the largest regional contributor to sea-level rise. The trend of increasing discharge for both glaciers, however, appears to have leveled off since 2008.

The 17O excess (Δ17O = δ17O−0.52 × δ18O) of sulfate and nitrate reflects the relative importance of their different production pathways in the atmosphere. A new record of sulfate and nitrate Δ17O spanning the last 2400 years from the West Antarctic Ice Sheet Divide ice core project shows significant changes in both sulfate and nitrate Δ17O in the most recent 200 years, indicating changes in their formation pathways. The sulfate Δ17O record exhibits a 1.1 ‰ increase in the early 19th century from (2.4 ± 0.2) ‰ to (3.5 ± 0.2) ‰, which suggests that an additional 12–18% of sulfate formation occurs via aqueous-phase production by O3, relative to that in the gas phase. Nitrate Δ17O gradually decreases over the whole record, with a more rapid decrease between the mid-19th century and the present day of 5.6 ‰, indicating an increasing importance of RO2 in NOx cycling between the mid-19th century and the present day in the mid- to high-latitude Southern Hemisphere. The former has implications for the climate impacts of sulfate aerosol, while the latter has implications for the tropospheric O3 production rate in remote low-NOx environments. Using other ice core observations, we rule out drivers for these changes other than variability in extratropical oxidant (OH, O3, RO2, H2O2, and reactive halogens) concentrations. However, assuming OH, H2O2, and O3 are the main oxidants contributing to sulfate formation, Monte Carlo box model simulations require a large (≥ 260%) increase in the O3 / OH mole fraction ratio over the Southern Ocean in the early 19th century to match the sulfate Δ17O record. This unlikely scenario points to a~deficiency in our understanding of sulfur chemistry and suggests other oxidants may play an important role in sulfate formation in the mid- to high-latitude marine boundary layer. The observed decrease in nitrate Δ17O since the mid-19th century is most likely due to an increased importance of RO2 over O3 in NOx cycling and can be explained by a 60–90% decrease in the O3 / RO2 mole fraction ratio in the extratropical Southern Hemisphere NOx-source regions.

We present high-resolution measurements of carbon monoxide (CO) concentrations from a shallow ice core of the North Greenland Eemian Ice Drilling project (NEEM-2011-S1). An optical-feedback cavity-enhanced absorption spectrometer (OF-CEAS) coupled to a continuous melter system performed continuous, online analysis during a four-week measurement campaign. This analytical setup generated stable measurements of CO concentrations with an external precision of 7.8 ppbv (1σ), based on repeated analyses of equivalent ice core sections. However, this first application of this measurement technique suffered from a poorly constrained procedural blank of 48 ± 25 ppbv and poor accuracy because an absolute calibration was not possible. The NEEM-2011-S1 CO record spans 1800 yr and the long-term trends within the most recent section of this record (i.e., post 1700 AD) resemble the existing discrete CO measurements from the Eurocore ice core. However, the CO concentration is highly variable (75–1327 ppbv range) throughout the ice core with high frequency (annual scale), high amplitude spikes characterizing the record. These CO signals are too abrupt and rapid to reflect atmospheric variability and their prevalence largely prevents interpretation of the record in terms of atmospheric CO variation. The abrupt CO spikes are likely the result of in situ production occurring within the ice itself, although the unlikely possibility of CO production driven by non-photolytic, fast kinetic processes within the continuous melter system cannot be excluded. We observe that 68% of the CO spikes are observed in ice layers enriched with pyrogenic aerosols. Such aerosols, originating from boreal biomass burning emissions, contain organic compounds, which may be oxidized or photodissociated to produce CO within the ice. However, the NEEM-2011-S1 record displays an increase of ~0.05 ppbv yr−1 in baseline CO level prior to 1700 AD (129 m depth) and the concentration remains elevated, even for ice layers depleted in dissolved organic carbon (DOC). Thus, the processes driving the likely in situ production of CO within the NEEM ice may involve multiple, complex chemical pathways not all related to past fire history and require further investigation.

As part of the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP), we evaluate the historical black carbon (BC) aerosols simulated by 8 ACCMIP models against observations including 12 ice core records, long-term surface mass concentrations, and recent Arctic BC snowpack measurements. We also estimate BC albedo forcing by performing additional simulations using offline models with prescribed meteorology from 1996–2000. We evaluate the vertical profile of BC snow concentrations from these offline simulations using the recent BC snowpack measurements. Despite using the same BC emissions, the global BC burden differs by approximately a factor of 3 among models due to differences in aerosol removal parameterizations and simulated meteorology: 34 Gg to 103 Gg in 1850 and 82 Gg to 315 Gg in 2000. However, the global BC burden from preindustrial to present-day increases by 2.5–3 times with little variation among models, roughly matching the 2.5-fold increase in total BC emissions during the same period. We find a large divergence among models at both Northern Hemisphere (NH) and Southern Hemisphere (SH) high latitude regions for BC burden and at SH high latitude regions for deposition fluxes. The ACCMIP simulations match the observed BC surface mass concentrations well in Europe and North America except at Ispra. However, the models fail to predict the Arctic BC seasonality due to severe underestimations during winter and spring. The simulated vertically resolved BC snow concentrations are, on average, within a factor of 2–3 of the BC snowpack measurements except for Greenland and the Arctic Ocean. For the ice core evaluation, models tend to adequately capture both the observed temporal trends and the magnitudes at Greenland sites. However, models fail to predict the decreasing trend of BC depositions/ice core concentrations from the 1950s to the 1970s in most Tibetan Plateau ice cores. The distinct temporal trend at the Tibetan Plateau ice cores indicates a strong influence from Western Europe, but the modeled BC increases in that period are consistent with the emission changes in Eastern Europe, the Middle East, South and East Asia. At the Alps site, the simulated BC suggests a strong influence from Europe, which agrees with the Alps ice core observations. At Zuoqiupu on the Tibetan Plateau, models successfully simulate the higher BC concentrations observed during the non-monsoon season compared to the monsoon season but overpredict BC in both seasons. Despite a large divergence in BC deposition at two Antarctic ice core sites, some models with a BC lifetime of less than 7 days are able to capture the observed concentrations. In 2000 relative to 1850, globally and annually averaged BC surface albedo forcing from the offline simulations ranges from 0.014 to 0.019 W m−2 among the ACCMIP models. Comparing offline and online BC albedo forcings computed by some of the same models, we find that the global annual mean can vary by up to a factor of two because of different aerosol models or different BC-snow parameterizations and snow cover. The spatial distributions of the offline BC albedo forcing in 2000 show especially high BC forcing (i.e., over 0.1 W m−2) over Manchuria, Karakoram, and most of the Former USSR. Models predict the highest global annual mean BC forcing in 1980 rather than 2000, mostly driven by the high fossil fuel and biofuel emissions in the Former USSR in 1980.

When contaminated by absorbing particles, such as refractory black carbon (rBC) and continental dust, snow's albedo decreases and thus its absorption of solar radiation increases, thereby hastening snowmelt. For this reason, an understanding of rBC's affect on snow albedo, melt processes, and radiation balance is critical for water management, especially in a changing climate. Measurements of rBC in a sequence of snow pits and surface snow samples in the eastern Sierra Nevada of California during the snow accumulation and ablation seasons of 2009 show that concentrations of rBC were enhanced sevenfold in surface snow (~25 ng g–1) compared to bulk values in the snowpack (~3 ng g–1). Unlike major ions, which were preferentially released during the initial melt, rBC and continental dust were retained in the snow, enhancing concentrations well into late spring, until a final flush occurred during the ablation period. We estimate a combined rBC and continental dust surface radiative forcing of 20 to 40 W m−2 during April and May, with dust likely contributing a greater share of the forcing.

We present an automated process for determining the annual layer chronology of an ice-core with a strong annual signal, utilising the hydrogen peroxide record from an Antarctic Peninsula ice-core as a test signal on which to count annual cycles and explain the methods. The signal is de-trended and normalised before being split into sections with a deterministic cycle count and those that need more attention. Possible reconstructions for the uncertain sections are determined which could be used as a visual aid for manual counting, and a simple method for assigning probability measures to each reconstruction is discussed. The robustness of this process is explored by applying it to versions of two different chemistry signals from the same stretch of the NGRIP (North Greenland Ice Core Project) ice-core, which shows more variation in annual layer thickness, with and without thinning to mimic poorer quality data. An adapted version of these methods is applied to the more challenging non-sea-salt sulphur signal from the same Antarctic Peninsula core from which the hydrogen peroxide signal was taken. These methods could readily be adapted for use on much longer datasets, thereby reducing manual effort and providing a robust automated layer-counting methodology.

Refractory black carbon aerosols (rBC) emitted by biomass burning (fires) and fossil fuel combustion, affect global climate and atmospheric chemistry. In the Southern Hemisphere (SH), rBC is transported in the atmosphere from low- and mid-latitudes to Antarctica and deposited to the polar ice sheet preserving a history of emissions and atmospheric transport. Here, we present two high-resolution Antarctic rBC ice core records drilled from the West Antarctic Ice Sheet divide and Law Dome on the periphery of the East Antarctic ice sheet. Separated by ~3500 km, the records span calendar years 1850–2001 and reflect the rBC distribution over the Indian and Pacific ocean sectors of the Southern Ocean. Concentrations of rBC in the ice cores displayed significant variability at annual to decadal time scales, notably in ENSO-QBO and AAO frequency bands. The delay observed between rBC and ENSO variability suggested that ENSO does not directly affect rBC transport, but rather continental hydrology, subsequent fire regimes, and aerosol emissions. From 1850 to 1950, the two ice core records were uncorrelated but were highly correlated from 1950 to 2002 (cross-correlation coefficient at annual resolution: r = 0.54, p < 0.01) due to a common decrease in rBC variability. The decrease in ice-core rBC from the 1950s to late 1980s displays similarities with inventories of SH rBC grass fires and biofuel emissions, which show reduced emission estimates over that period.

Refractory black carbon aerosols (rBC) from biomass burning and fossil fuel combustion are deposited to the Antarctic ice sheet and preserve a history of emissions and long-range transport from low- and mid-latitudes. Antarctic ice core rBC records may thus provide information with respect to past combustion aerosol emissions and atmospheric circulation. Here, we present six East Antarctic ice core records of rBC concentrations and fluxes covering the last two centuries with approximately annual resolution (cal. yr. 1800 to 2000). The ice cores were drilled in disparate regions of the high East Antarctic ice sheet, at different elevations and net snow accumulation rates. Annual rBC concentrations were log-normally distributed and geometric means of annual concentrations ranged from 0.10 to 0.18 μg kg−1. Average rBC fluxes were determined over the time periods 1800 to 2000 and 1963 to 2000 and ranged from 3.4 to 15.5 μg m−2 a−1 and 3.6 to 21.8 μg m−2 a−1, respectively. Geometric mean concentrations spanning 1800 to 2000 increased linearly with elevation at a rate of 0.025 μg kg−1/500 m. Spectral analysis of the records revealed significant decadal-scale variability, which at several sites was comparable to decadal ENSO variability.

Desert dust perturbs climate by directly and indirectly interacting with incoming solar and outgoing long wave radiation, thereby changing precipitation and temperature, in addition to modifying ocean and land biogeochemistry. While we know that desert dust is sensitive to perturbations in climate and human land use, previous studies have been unable to determine whether humans were increasing or decreasing desert dust in the global average. Here we present observational estimates of desert dust based on paleodata proxies showing a doubling of desert dust during the 20th century over much, but not all the globe. Large uncertainties remain in estimates of desert dust variability over 20th century due to limited data. Using these observational estimates of desert dust change in combination with ocean, atmosphere and land models, we calculate the net radiative effect of these observed changes (top of atmosphere) over the 20th century to be −0.14 ± 0.11 W/m2 (1990–1999 vs. 1905–1914). The estimated radiative change due to dust is especially strong between the heavily loaded 1980–1989 and the less heavily loaded 1955–1964 time periods (−0.57 ± 0.46 W/m2), which model simulations suggest may have reduced the rate of temperature increase between these time periods by 0.11 °C. Model simulations also indicate strong regional shifts in precipitation and temperature from desert dust changes, causing 6 ppm (12 PgC) reduction in model carbon uptake by the terrestrial biosphere over the 20th century. Desert dust carries iron, an important micronutrient for ocean biogeochemistry that can modulate ocean carbon storage; here we show that dust deposition trends increase ocean productivity by an estimated 6% over the 20th century, drawing down an additional 4 ppm (8 PgC) of carbon dioxide into the oceans. Thus, perturbations to desert dust over the 20th century inferred from observations are potentially important for climate and biogeochemistry, and our understanding of these changes and their impacts should continue to be refined.

We present and discuss a new dataset of gridded emissions covering the historical period (1850–2000) in decadal increments at a horizontal resolution of 0.5° in latitude and longitude. The primary purpose of this inventory is to provide consistent gridded emissions of reactive gases and aerosols for use in chemistry model simulations needed by climate models for the Climate Model Intercomparison Program #5 (CMIP5) in support of the Intergovernmental Panel on Climate Change (IPCC) Fifth Assessment report (AR5). Our best estimate for the year 2000 inventory represents a combination of existing regional and global inventories to capture the best information available at this point; 40 regions and 12 sectors are used to combine the various sources. The historical reconstruction of each emitted compound, for each region and sector, is then forced to agree with our 2000 estimate, ensuring continuity between past and 2000 emissions. Simulations from two chemistry-climate models are used to test the ability of the emission dataset described here to capture long-term changes in atmospheric ozone, carbon monoxide and aerosol distributions. The simulated long-term change in the Northern mid-latitudes surface and mid-troposphere ozone is not quite as rapid as observed. However, stations outside this latitude band show much better agreement in both present-day and long-term trend. The model simulations indicate that the concentration of carbon monoxide is underestimated at the Mace Head station; however, the long-term trend over the limited observational period seems to be reasonably well captured. The simulated sulfate and black carbon deposition over Greenland is in very good agreement with the ice-core observations spanning the simulation period. Finally, aerosol optical depth and additional aerosol diagnostics are shown to be in good agreement with previously published estimates and observations.

No abstract available.