Comment on “Low time resolution analysis of polar ice cores cannot detect impulsive nitrate events” by D.F. Smart et al.

Smart et al. (2014) suggested that the detection of nitrate spikes in polar ice cores from solar energetic particle (SEP) events could be achieved if an analytical system with sufficiently high resolution was used. Here we show that the spikes they associate with SEP events are not reliably recorded in cores from the same location, even when the resolution is clearly adequate. We explain the processes that limit the effective resolution of ice cores. Liquid conductivity data suggest that the observed spikes are associated with sodium or another nonacidic cation, making it likely that they result from deposition of sea salt or similar aerosol that has scavenged nitrate, rather than from a primary input of nitrate in the troposphere. We consider that there is no evidence at present to support the identification of any spikes in nitrate as representing SEP events. Although such events undoubtedly create nitrate in the atmosphere, we see no plausible route to using nitrate spikes to document the statistics of such events.

Millennial and sub-millennial scale climatic variations recorded in polar ice cores over the last glacial period

Since its discovery in Greenland ice cores, the millennial scale climatic variability of the last glacial period has been increasingly documented at all latitudes with studies focusing mainly on Marine Isotopic Stage 3 (MIS 3; 28–60 thousand of years before present, hereafter ka) and characterized by short Dansgaard-Oeschger (DO) events. Recent and new results obtained on the EPICA and NorthGRIP ice cores now precisely describe the rapid variations of Antarctic and Greenland temperature during MIS 5 (73.5–123 ka), a time period corresponding to relatively high sea level. The results display a succession of abrupt events associated with long Greenland InterStadial phases (GIS) enabling us to highlight a sub-millennial scale climatic variability depicted by (i) short-lived and abrupt warming events preceding some GIS (precursor-type events) and (ii) abrupt warming events at the end of some GIS (rebound-type events). The occurrence of these sub-millennial scale events is suggested to be driven by the insolation at high northern latitudes together with the internal forcing of ice sheets. Thanks to a recent NorthGRIP-EPICA Dronning Maud Land (EDML) common timescale over MIS 5, the bipolar sequence of climatic events can be established at millennial to sub-millennial timescale. This shows that for extraordinary long stadial durations the accompanying Antarctic warming amplitude cannot be described by a simple linear relationship between the two as expected from the bipolar seesaw concept. We also show that when ice sheets are extensive, Antarctica does not necessarily warm during the whole GS as the thermal bipolar seesaw model would predict, questioning the Greenland ice core temperature records as a proxy for AMOC changes throughout the glacial period.

Interhemispheric coupling, the West Antarctic Ice Sheet and warm Antarctic interglacials

Ice core evidence indicates that even though atmospheric CO2 concentrations did not exceed ~300 ppm at any point during the last 800 000 years, East Antarctica was at least ~3–4 °C warmer than preindustrial (CO2~280 ppm) in each of the last four interglacials. During the previous three interglacials, this anomalous warming was short lived (~3000 years) and apparently occurred before the completion of Northern Hemisphere deglaciation. Hereafter, we refer to these periods as "Warmer than Present Transients" (WPTs). We present a series of experiments to investigate the impact of deglacial meltwater on the Atlantic Meridional Overturning Circulation (AMOC) and Antarctic temperature. It is well known that a slowed AMOC would increase southern sea surface temperature (SST) through the bipolar seesaw and observational data suggests that the AMOC remained weak throughout the terminations preceding WPTs, strengthening rapidly at a time which coincides closely with peak Antarctic temperature. We present two 800 kyr transient simulations using the Intermediate Complexity model GENIE-1 which demonstrate that meltwater forcing generates transient southern warming that is consistent with the timing of WPTs, but is not sufficient (in this single parameterisation) to reproduce the magnitude of observed warmth. In order to investigate model and boundary condition uncertainty, we present three ensembles of transient GENIE-1 simulations across Termination II (135 000 to 124 000 BP) and three snapshot HadCM3 simulations at 130 000 BP. Only with consideration of the possible feedback of West Antarctic Ice Sheet (WAIS) retreat does it become possible to simulate the magnitude of observed warming.

Perennial Snow and Ice Variations (2000-2008) in the Arctic Circumpolar Land Area from Satellite Observations

Perennial snow and ice (PSI) extent is an important parameter of mountain environments with regard to its involvement in the hydrological cycle and the surface energy budget. We investigated interannual variations of PSI in nine mountain regions of interest (ROI) between 2000 and 2008. For that purpose, a novel MODIS data set processed at the Canada Centre for Remote Sensing at 250 m spatial resolution was utilized. The extent of PSI exhibited significant interannual variations, with coefficients of variation ranging from 5% to 81% depending on the ROI. A strong negative relationship was found between PSI and positive degree-days (threshold 0°C) during the summer months in most ROIs, with linear correlation coefficients (r) being as low as r = −0.90. In the European Alps and Scandinavia, PSI extent was significantly correlated with annual net glacier mass balances, with r = 0.91 and r = 0.85, respectively, suggesting that MODIS-derived PSI extent may be used as an indicator of net glacier mass balances. Validation of PSI extent in two land surface classifications for the years 2000 and 2005, GLC-2000 and Globcover, revealed significant discrepancies of up to 129% for both classifications. With regard to the importance of such classifications for land surface parameterizations in climate and land surface process models, this is a potential source of error to be investigated in future studies. The results presented here provide an interesting insight into variations of PSI in several ROIs and are instrumental for our understanding of sensitive mountain regions in the context of global climate change assessment.

An ice-core based history of Siberian forest fires since AD 1250

Forest fires play a key role in the global carbon cycle and thus, can affect regional and global climate. Although fires in extended areas of Russian boreal forests have a considerable influence on atmospheric greenhouse gas and soot concentrations, estimates of their impact on climate are hampered by a lack of data on the history of forest fires. Especially regions with strong continental climate are of high importance due to an intensified development of wildfires. In this study we reconstruct the fire history of Southern Siberia during the past 750 years using ice-core based nitrate, potassium, and charcoal concentration records from Belukha glacier in the continental Siberian Altai. A period of exceptionally high forest-fire activity was observed between AD 1600 and 1680, following an extremely dry period AD 1540-1600. Ice-core pollen data suggest distinct forest diebacks and the expansion of steppe in response to dry climatic conditions. Coherence with a paleoenvironmental record from the 200 km distant Siberian lake Teletskoye shows that the vegetational shift AD 1540-1680, the increase in fire activity AD 1600-1680, and the subsequent recovery of forests AD 1700 were of regional significance. Dead biomass accumulation in response to drought and high temperatures around AD 1600 probably triggered maximum forest-fire activity AD 1600-1680. The extreme dry period in the 16th century was also observed at other sites in Central Asia and is possibly associated with a persistent positive mode of the Pacific Decadal Oscillation (PDO). No significant increase in biomass burning occurred in the Altai region during the last 300 years, despite strongly increasing temperatures and human activities. Our results imply that precipitation changes controlled fire-regime and vegetation shifts in the Altai region during the past 750 years. We conclude that high sensitivity of ecosystems to occasional decadal-scale drought events may trigger unprecedented environmental reorganizations under global-warming conditions.

Recent increase in black carbon concentrations from a Mt. Everest ice core spanning 1860-2000 AD

A Mt. Everest ice core spanning 1860–2000 AD and analyzed at high resolution for black carbon (BC) using a Single Particle Soot Photometer (SP2) demonstrates strong seasonality, with peak concentrations during the winter-spring, and low concentrations during the summer monsoon season. BC concentrations from 1975–2000 relative to 1860–1975 have increased approximately threefold, indicating that BC from anthropogenic sources is being transported to high elevation regions of the Himalaya. The timing of the increase in BC is consistent with BC emission inventory data from South Asia and the Middle East, however since 1990 the ice core BC record does not indicate continually increasing BC concentrations. The Everest BC and dust records provide information about absorbing impurities that can contribute to glacier melt by reducing the albedo of snow and ice. There is no increasing trend in dust concentrations since 1860, and estimated surface radiative forcing due to BC in snow exceeds that of dust in snow. This suggests that a reduction in BC emissions may be an effective means to reduce the effect of absorbing impurities on snow albedo and melt, which affects Himalayan glaciers and the availability of water resources in major Asian rivers.

Ion fractionation in young sea ice from Kongsfjorden, Svalbard

The fractionation of major sea-water ions, or deviation in their relative concentrations from Standard Mean Ocean Water ratios, has been frequently observed in sea ice. It is generally thought to be associated with precipitation of solid salts at certain eutectic temperatures. The variability found in bulk sea-ice samples indicates that the fractionation of ions depends on the often unknown thermal history of sea ice, which affects the structure of pore networks and fate of solid salts within them. Here we investigate the distribution of ions in Arctic sea ice that is a few weeks old with a reconstructible thermal history. We separate the centrifugable (interconnected) and entrapped (likely disconnected) contributions to the ice salinity and determine their ion fractionation signatures. The results indicate that differential diffusion of ions, rather than eutectic precipitation of cryohydrates, has led to significant ion fractionation. The finding emphasizes the role of coupled diffusive–convective salt transport through complex pore networks in shaping the biogeochemistry of sea ice.

Annual layering in the NGRIP ice core during the Eemian

The Greenland NGRIP ice core continuously covers the period from present day back to 123 ka before present, which includes several thousand years of ice from the previous interglacial period, MIS 5e or the Eemian. In the glacial part of the core, annual layers can be identified from impurity records and visual stratigraphy, and stratigraphic layer counting has been performed back to 60 ka. In the deepest part of the core, however, the ice is close to the pressure melting point, the visual stratigraphy is dominated by crystal boundaries, and annual layering is not visible to the naked eye. In this study, we apply a newly developed setup for high-resolution ice core impurity analysis to produce continuous records of dust, sodium and ammonium concentrations as well as conductivity of melt water. We analyzed three 2.2 m sections of ice from the Eemian and the glacial inception. In all of the analyzed ice, annual layers can clearly be recognized, most prominently in the dust and conductivity profiles. Part of the samples is, however, contaminated in dust, most likely from drill liquid. It is interesting that the annual layering is preserved despite a very active crystal growth and grain boundary migration in the deep and warm NGRIP ice. Based on annual layer counting of the new records, we determine a mean annual layer thickness close to 11 mm for all three sections, which, to first order, confirms the modeled NGRIP time scale (ss09sea). The counting does, however, suggest a longer duration of the climatically warmest part of the NGRIP record (MIS5e) of up to 1 ka as compared to the model estimate. Our results suggest that stratigraphic layer counting is possible basically throughout the entire NGRIP ice core, provided sufficiently highly-resolved profiles become available.

Water isotopic ratios from a continuously melted ice core sample

A new technique for on-line high resolution isotopic analysis of liquid water, tailored for ice core studies is presented. We built an interface between a Wavelength Scanned Cavity Ring Down Spectrometer (WS-CRDS) purchased from Picarro Inc. and a Continuous Flow Analysis (CFA) system. The system offers the possibility to perform simultaneuous water isotopic analysis of δ18O and δD on a continuous stream of liquid water as generated from a continuously melted ice rod. Injection of sub μl amounts of liquid water is achieved by pumping sample through a fused silica capillary and instantaneously vaporizing it with 100% efficiency in a~home made oven at a temperature of 170 °C. A calibration procedure allows for proper reporting of the data on the VSMOW–SLAP scale. We apply the necessary corrections based on the assessed performance of the system regarding instrumental drifts and dependance on the water concentration in the optical cavity. The melt rates are monitored in order to assign a depth scale to the measured isotopic profiles. Application of spectral methods yields the combined uncertainty of the system at below 0.1‰ and 0.5‰ for δ18O and δD, respectively. This performance is comparable to that achieved with mass spectrometry. Dispersion of the sample in the transfer lines limits the temporal resolution of the technique. In this work we investigate and assess these dispersion effects. By using an optimal filtering method we show how the measured profiles can be corrected for the smoothing effects resulting from the sample dispersion. Considering the significant advantages the technique offers, i.e. simultaneuous measurement of δ18O and δD, potentially in combination with chemical components that are traditionally measured on CFA systems, notable reduction on analysis time and power consumption, we consider it as an alternative to traditional isotope ratio mass spectrometry with the possibility to be deployed for field ice core studies. We present data acquired in the field during the 2010 season as part of the NEEM deep ice core drilling project in North Greenland.

A sublimation technique for high-precision measurements of d13CO2 and mixing ratios of CO2 and N2O from air trapped in ice cores

In order to provide high precision stable carbon isotope ratios (δ13CO2 or δ13C of CO2) from small bubbly, partially and fully clathrated ice core samples we developed a new method based on sublimation coupled to gas chromatography-isotope ratio mass spectrometry (GC-IRMS). In a first step the trapped air is quantitatively released from ~30 g of ice and CO2 together with N2O are separated from the bulk air components and stored in a miniature glass tube. In an off-line step, the extracted sample is introduced into a helium carrier flow using a minimised tube cracker device. Prior to measurement, N2O and organic sample contaminants are gas chromatographically separated from CO2. Pulses of a CO2/N2O mixture are admitted to the tube cracker and follow the path of the sample through the system. This allows an identical treatment and comparison of sample and standard peaks. The ability of the method to reproduce δ13C from bubble and clathrate ice is verified on different ice cores. We achieve reproducibilities for bubble ice between 0.05 ‰ and 0.07 ‰ and for clathrate ice between 0.05 ‰ and 0.09 ‰ (dependent on the ice core used). A comparison of our data with measurements on bubble ice from the same ice core but using a mechanical extraction device shows no significant systematic offset. In addition to δ13C, the CO2 and N2O mixing ratios can be volumetrically derived with a precision of 2 ppmv and 8 ppbv, respectively.

Polygonal cracks in the seasonal semi-translucent CO2 ice layer in Martian polar areas

In this paper, we use morphological and numerical methods to test the hypothesis that seasonally formed fracture patterns in the Martian polar regions result from the brittle failure of seasonal CO2 slab ice. The observations by the High Resolution Imaging Science Experiment (HiRISE) of polar regions of Mars show very narrow dark elongated linear patterns that are observed during some periods of time in spring, disappear in summer and re-appear again in the following spring. They are repeatedly formed in the same areas but they do not repeat the exact pattern from year to year. This leads to the conclusion that they are cracks formed in the seasonal ice layer. Some of models of seasonal surface processes rely on the existence of a transparent form of CO2 ice, so-called slab ice. For the creation of the observed cracks the ice is required to be a continuous media, not an agglomeration of relatively separate particles like a firn. The best explanation for our observations is a slab ice with relatively high transparency in the visible wavelength range. This transparency allows a solid state green-house effect to act underneath the ice sheet raising the pressure by sublimation from below. The trapped gas creates overpressure and the ice sheet breaks at some point creating the observed cracks. We show that the times when the cracks appear are in agreement with the model calculation, providing one more piece of evidence that CO2 slab ice covers polar areas in spring.

Impact of Holocene climate changes on alpine and treeline vegetation at Sanetsch Pass, Bernese Alps, Switzerland

In order to infer reactions of treeline and alpine vegetation to climatic change, past vegetation changes are reconstructed on the basis of pollen, macrofossil and charcoal analysis. The sampled sediment cores originate from the small pond Emines, located at the Sanetsch Pass (connecting the Valais and Bern, Switzerland) at an altitude of 2288 m a.s.l. Today's treeline is at ca. 2200 m a.s.l. in the area, though due to special pass (saddle) conditions it is locally depressed to ca. 2060 m a.s.l. Our results reveal that the area around Emines was covered by treeless alpine vegetation during most of the past 12,000 years. Single individuals of Betula, Larix decidua and possibly Pinus cembra occurred during the Holocene. Major centennial to millennial-scale responses of treeline vegetation to climatic changes are evident. However, alpine vegetation composition remained rather stable between 11,500 and 6000 cal. BP, showing that Holocene climatic changes of +/− 1 °C hardly influenced the local vegetation at Emines. The rapid warming of 3–4 °C at the Late Glacial/Holocene transition (11,600 cal. BP) caused significant altitudinal displacements of alpine species that were additionally affected by the rapid upward movement of trees and shrubs. Since the beginning of the Neolithic, vegetation changes at Sanetsch Pass resulted from a combination of climate change and human impact. Anthropogenic fire increase and land-use change combined with a natural change from subcontinental to more oceanic climate during the second half of the Holocene led to the disappearance of P. cembra in the study area, but favoured the occurrence of Picea abies and Alnus viridis. The mid- to late-Holocene decline of Abies alba was primarily a consequence of human impact, since this mesic species should have benefitted from a shift to more oceanic conditions. Future alpine vegetation changes will be a function of the amplitude and rapidity of global warming as well as human land use. Our results imply that alpine vegetation at our treeline pass site was never replaced by forests since the last ice-age. This may change in the future if anticipated climate change will induce upslope migration of trees. The results of this study emphasise the necessity of climate change mitigation in order to prevent biodiversity losses as a consequence of unprecedented community and species displacement in response to climatic change.

Continuous measurements of methane mixing ratios from ice cores

This work presents a new, field-deployable technique for continuous, high-resolution measurements of methane mixing ratios from ice cores. The technique is based on a continuous flow analysis system, where ice core samples cut along the long axis of an ice core are melted continuously. The past atmospheric air contained in the ice is separated from the melt water stream via a system for continuous gas extraction. The extracted gas is dehumidified and then analyzed by a Wavelength Scanned-Cavity Ring Down Spectrometer for methane mixing ratios. We assess the performance of the new measurement technique in terms of precision (±0.8 ppbv, 1σ), accuracy (±8 ppbv), temporal (ca. 100 s), and spatial resolution (ca. 5 cm). Using a firn air transport model, we compare the resolution of the measurement technique to the resolution of the atmospheric methane signal as preserved in ice cores in Greenland. We conclude that our measurement technique can resolve all climatically relevant variations as preserved in the ice down to an ice depth of at least 1980 m (66 000 yr before present) in the North Greenland Eemian Ice Drilling ice core. Furthermore, we describe the modifications, which are necessary to make a commercially available spectrometer suitable for continuous methane mixing ratio measurements from ice cores.

Centennial mineral dust variability in high-resolution ice core data from Dome C, Antarctica

Ice core data from Antarctica provide detailed insights into the characteristics of past climate, atmospheric circulation, as well as changes in the aerosol load of the atmosphere. We present high-resolution records of soluble calcium (Ca2+), non-sea-salt soluble calcium (nssCa2+), and particulate mineral dust aerosol from the East Antarctic Plateau at a depth resolution of 1 cm, spanning the past 800 000 years. Despite the fact that all three parameters are largely dust-derived, the ratio of nssCa2+ to particulate dust is dependent on the particulate dust concentration itself. We used principal component analysis to extract the joint climatic signal and produce a common high-resolution record of dust flux. This new record is used to identify Antarctic warming events during the past eight glacial periods. The phasing of dust flux and CO2 changes during glacial-interglacial transitions reveals that iron fertilization of the Southern Ocean during the past nine glacial terminations was not the dominant factor in the deglacial rise of CO2 concentrations. Rapid changes in dust flux during glacial terminations and Antarctic warming events point to a rapid response of the southern westerly wind belt in the region of southern South American dust sources on changing climate conditions. The clear lead of these dust changes on temperature rise suggests that an atmospheric reorganization occurred in the Southern Hemisphere before the Southern Ocean warmed significantly.