Hostname: page-component-7c8c6479df-27gpq Total loading time: 0 Render date: 2024-03-29T10:09:57.594Z Has data issue: false hasContentIssue false

Termination V in the Vostok (Antarctica) ice core

Published online by Cambridge University Press:  08 September 2017

Makoto Suwa*
Affiliation:
Department of Geosciences, Princeton University, Guyot Hall, Washington Road, Princeton, New Jersey 08544-1003, USA E-mail: makotosuwa@alumni.duke.edu
Rights & Permissions [Opens in a new window]

Abstract

The age–depth relationship of the Vostok (Antarctica) ice core has been reconstructed in the depth interval 3300–3347 m, by comparing three gas properties in ice (CO2, CH4 and δ18Oatm) with those in the EPICA Dome C (Antarctica) core. Fourteen Vostok depths were examined in this interval, and it was found that nine samples are uniquely dated if candidate ages are restricted to the interval between 400 and 650 kyr. One of these samples is uniquely dated without restriction. The analysis supports previous reports that this section contains ice from Termination V, but that the stratigraphic order of ice is reversed. The top of the overturned layer lies between 3316 and 3319 m. At least one other stratigraphic disturbance was found between 3340 and 3343 m, as indicated by another reversal of the age–depth relationship. Finally, the oldest ice in this section is dated at ≥440 kyr, confirming the existence of ice from the cold marine isotope stage (MIS) 12 interval.

Type
Research Article
Copyright
Copyright © International Glaciological Society 2008

1. Introduction

The 3623 m long Vostok 5G core was drilled in the East Antarctic plateau between 1991 and 1997. The Vostok core contains a record of climate over at least the past ∼400 kyr (Reference PetitPetit and others, 1999). However, the stratigraphic integrity of the section below ∼3310 m has been questioned (Reference PetitPetit and others, 1999). In fact, stratigraphic disturbance of the bottom part of the ice sheet has been observed in some of the other long ice cores from polar regions, such as the GRIP and GISP2 ice cores from Summit, Greenland (Reference Chappellaz, Brook, Blunier and MalaizéChappellaz and others 1997; Reference LandaisLandais and others, 2003; Reference Suwa, von Fischer, Bender, Landais and BrookSuwa and others, 2006). More recently, Reference RaynaudRaynaud and others (2005) suggested that the section of the Vostok 5 G core between 3321.0 and 3344.9 m might contain ice from Termination V, but in reversed stratigraphic order. They reached the conclusion by visually comparing CH4 and CO2 measurements of samples from this section of the Vostok core with those of the EPICA Dome C (Antarctica) core, which now extends back to >650 kyr (Reference Siegenthaler, Stocker, Monnin, Lüthi, Schwander and StaufferSiegenthaler and others, 2005; Reference SpahniSpahni and others, 2005).

Considering that glacial terminations are periods of intense scientific interest, and the precious nature of ice samples this old (to date, only the EPICA Dome C and Dome Fuji cores include Termination V ice), it is important to reconstruct the stratigraphy of this section of the Vostok core. Therefore, this paper aims to constrain the age–depth relationship of the Vostok ice core for the section between 3300 and 3347 m using gas properties. In addition to the CH4 and CO2 records considered by Reference RaynaudRaynaud and others (2005), we incorporate measurements of δ18O of atmospheric O218Oatm) into the age reconstruction. We use the recently published ∼800 kyr long δ18Oatm record of the EPICA Dome C core for the reference (Reference DreyfusDreyfus and others, 2007) and new δ18Oatm measurements of samples from Vostok for our analysis. We note that there is no indication that the three gas properties considered in this study had been significantly modified or contaminated during ice-core extraction, core transport and gas extraction. In fact, CO2 and CH4 are within the expected range (Reference PetitPetit and others, 1999) and so are δ18O and total gas content (Reference SuwaSuwa, 2007).

2. Reconstructing the Chronology Using CH4–δ18Oatm–CO2

The method used in this study is similar to that reported by Reference Suwa, von Fischer, Bender, Landais and BrookSuwa and others (2006), who reconstructed the chronology for the bottom sections of the GISP2 and GRIP ice cores. The basis of their method is that individual samples from the stratigraphically disturbed section are dated based on their gas concentrations, which can be compared with concentration–time profiles of a stratigraphically intact ice core (i.e. Dome C). The difference between their method and the method developed here is that, in addition to the two gas properties they used, δ18Oatm and CH4, we use CO2 as a third constraint. CO2 was not useful in reconstructing the GISP2/GRIP chronologies because it cannot be measured reliably in those cores (Reference AnklinAnklin and others, 1997; Reference Smith, Wahlen, Mastroianni and TaylorSmith and others, 1997).

2.1. CH4–δ18Oatm–CO2 records of the Vostok and EPICA Dome C ice cores

For the Vostok CH4 and CO2 records, we use values reported by Reference PetitPetit and others (1999). For δ18Oatm we combine measurements by Reference PetitPetit and others (1999) and the new set of δ18Oatm measurements. The method for determining δ18Oatm for the new dataset is similar to the method described in Reference Sowers, Bender and RaynaudSowers and others (1989). We use a ∼15g sample of ice for each measurement, and apply the double melt refreeze method to extract air trapped in bubbles and clathrate hydrates in an ice sample. We then correct for the gravitational enrichment by subtracting 2 × δ15N from δ18O of the same sample to give δ18Oatm (Reference Craig, Horibe and SowersCraig and others, 1988; Reference Schwander, Oeschger and LangwaySchwander, 1989). Unfortunately, δ18Oatm (and in some cases CO2) is not measured in exactly the same sample as CH4. Therefore, we first choose 14 samples where CH4 and CO2 are measured at the same depths. Values for δ18Oatm are then interpolated for these depths. We note that δ18Oatm changes gradually in this section, which indicates that small-scale mixing is not dominant here. Figure 1 shows CH4, CO2, δ18Oatm and δDice of the Vostok ice core for the section between 3200 and 3400 m. Black solid circles indicate where original measurements are located. Red squares indicate ‘interpolated’ depths that we use for the following analysis. δDice is the temperature proxy for Vostok and Dome C.

Fig. 1. The CH4, CO2, δ18Oatm and δDice records for the section between 3200 and 3400 m of the Vostok ice core. Black solid circles indicate where original measurements are located. Red squares indicate values at interpolated depths used for the analysis.

Next, we derive the ‘reference’ line, which we assume is the true atmospheric history, based on the undisturbed record of EPICA Dome C. The ‘reference line’ in Figure 2 refers to the trajectory of δ18Oatm vs CH4 vs CO2 concentration between 400 and 650 kyr on the EPICA Dome C EDC2 timescale (EPICA Community Members, 2004). We use the published CO2 (Reference Siegenthaler, Stocker, Monnin, Lüthi, Schwander and StaufferSiegenthaler and others, 2005), CH4 (Spahni and others, 2007) and δ18Oatm (Reference DreyfusDreyfus and others, 2007) records to define the reference line.

Fig. 2. (a) An x-y-z plot of CO2 vs δ18Oatm vs CH4. The green curve traces the trajectory of CO2, CH4 and δ18Oatm points between 400 and 650 kyr in the EPICA Dome C core. Closed circles colored in shades of red and blue are Vostok data from the disturbed section of the core; the color represents δDice at the depth of sampling. The black curve connects Vostok samples in order of increasing depth. (b) As for (a) but shown only for CH4 and δ18Oatm. Gray ellipses surrounding each data point are uncertainty ranges. Black letters next to closed circles indicate the depths (m) of the Vostok sample. Orange labels next to yellow dots indicate ages (kyr) of Dome C samples.

The pooled standard deviations for EPICA Dome C and Vostok are: ±0.028‰ (1σ) (Reference DreyfusDreyfus and others, 2007) and ±0.05‰ (1σ), respectively, for δ18Oatm; ±15 ppb (1σ) (Spahni and others, 2007) and ±20 ppb (1σ) (Reference PetitPetit and others, 1999), respectively, for CH4; and ±1.5 ppm (1σ) (Reference Siegenthaler, Stocker, Monnin, Lüthi, Schwander and StaufferSiegenthaler and others, 2005) and ±3 ppm (1σ) (Reference PetitPetit and others, 1999), respectively, for CO2. These errors are added quadratically for each gas property. We adopt 2σ uncertainties, so that errors associated with δ18Oatm, CH4 and CO2 are ±0.11‰, ±50 ppb and ±6.7 ppm, respectively, in the following analysis.

2.2. Method

We first identify all compatible ages for each sample, which are our estimation of age uncertainties. This is done by seeking values for time t on the EPICA EDC2 timescale that meet the condition,

(1)

where CH4 EDC(t) is the CH4 concentration, δ18Oatm EDC(t) is the δ18Oatm value and CO2 EDC(t) is the CO2 concentration in EPICA Dome C at time t. Likewise, CH4Vtk(z), δ18Oatm Vtk(z) and CO2 Vtk(z) are CH4 concentration, δ18Oatm value and CO2 concentration at depth z in a sample from the section between 3300 and 3347 m of the Vostok core, and σ(CH4), σ18Oatm) and σ(CO2) are one standard deviation for CH4, d18Oatm and CO2, respectively, as described in section 2.1.

Second, we derive the ‘best-estimate age’ within all compatible age ranges derived according to Equation (1). The best-estimate age is defined as t which minimizes d, the distance between a sample and the reference line, calculated as

(2)

where the subscript n indicates that the term is normalized to a mean of 0 and standard deviation of 1. Normalization weights each of the three properties equally.

We first derive gas ages of samples, based on Equations (1) and (2), which are several thousand years younger than ages of ice at the same depths. This gas-age/ice-age difference (Δage) arises from the fact that the air is trapped in ice at ∼100 m below the surface of the ice sheet, and thus Δage should be taken into account when plotting proxy signals recorded in ice vs time. We compute Δage using the empirical model of Reference Herron and LangwayHerron and Langway (1980), which requires temperature and accumulation rate as input parameters. We estimate temperature from a simple linear regression of ΔT (deviation from modern Vostok temperature of −57.4°C) on δDice using values reported in Reference PetitPetit and others (1999), and we estimate accumulation rate from the equation used in Reference PetitPetit and others (1999). We understand that there are relatively large uncertainties associated with this model, especially for low-accumulation sites, but these uncertainties do not have a large impact on our results.

3. Results and Discussion

The derived ages for 14 samples examined in this study produce an excellent fit between the two ice-core records, for all four geochemical parameters. The age–depth relationship so derived (Fig. 3) allows us to draw the following five conclusions. First, 9 of 14 samples have only one compatible age range which is consistent with EPICA Dome C ages between 400 and 650 kyr. Samples from 3328, 3331 and 3334 m have multiple compatible age ranges, but all ages fall between 416 and 424 kyr. The bottommost two samples of our analysis, the samples from 3343 and 3346 m, also have two compatible age ranges. One is ∼440 kyr and the other is ∼460 and ∼450 kyr, respectively. Second, our age reconstruction indicates that the stratigraphic disturbance starts somewhere between 3316 and 3319 m. Third, four or five samples between 3319 and ∼3330 m are in reversed order. In this depth interval, one finds younger age as one goes deeper in the ice core. Fourth, it appears that there is a second stratigraphic disturbance between 3340 and 3343 m. Although the exact mechanism of development of these folds remains to be examined, layered rheological contrasts across the interglacial–glacial ice might have contributed to the process. Such folds are similar to the ones observed in the horizontal ice core retrieved from the Pâkitsoq ice-sheet margin site in Greenland (Reference Petrenko, Severinghaus, Brook, Reeh and SchaeferPetrenko and others, 2006). Lastly, the oldest sample in the section between 3300 and 3347 m is dated to at least ∼440 kyr on the EPICA Dome C EDC2 timescale.

Fig. 3. The age–depth relationship reconstructed in this study. Horizontal black bars indicate all compatible age ranges. The ‘best-estimate ages’ are plotted as solid circles. Shades of red and blue indicate dD of ice at the same depth.

Figure 4 shows the comparison of CO2, CH4, δ18Oatm and δDice records between EPICA Dome C and Vostok between 400 and 650 kyr. Black curves show the EPICA Dome C records, and blue circles are Vostok values plotted on their best-estimate ages. Our analysis supports the earlier suggestion by Reference RaynaudRaynaud and others (2005) that the temperature history during Termination V is similar at the two sites.

Fig. 4. CO2, CH4, δ18Oatm and δD records between 400 and 650 kyr. Black curves show EPICA Dome C. Blue circles indicate Vostok samples plotted on the best-estimate ages. Note that the y axis for δ18Oatm is reversed, and δDice is plotted on offset scales for EPICA Dome C and Vostok.

As stated above, we limited our age search between 400 and 650 kyr. Therefore, we cannot completely dismiss the possibility that these samples are from much older age intervals. However, we note that the EPICA Dome C record shows δ18Oatm of ∼1.5‰ during Termination V. This heavy δ18Oatm value is noteworthy as it is not reached for the rest of the published record which goes back to ∼800 kyr (Reference DreyfusDreyfus and others, 2007). Therefore, the heaviest δ18Oatm value we observed in the Vostok samples, which is ∼1.5‰, indicates that this sample is uniquely dated over the whole 800 kyr range.

Acknowledgements

I thank M. Bender, J. Severinghaus, R. Spahni and an anonymous reviewer for constructive comments on the manuscript. This work was supported by a gift from the Gary Comer Family Foundation and the Office of Polar Programs, US National Science Foundation.

References

Anklin, M. and 6 others. 1997. CO2 record between 40 and 8 kyr B.P. from the Greenland Ice Core Project ice core. J. Geophys. Res., 102(C12), 26,53926,546.CrossRefGoogle Scholar
Chappellaz, J., Brook, E., Blunier, T. and Malaizé, B.. 1997. CH4 and δ18O of O2 records from Antarctic and Greenland ice: a clue for stratigraphic disturbance in the bottom part of the Greenland Ice Core Project and the Greenland Ice Sheet Project 2 ice cores. J. Geophys. Res., 102(C12), 26,54726,557.CrossRefGoogle Scholar
Craig, H., Horibe, Y. and Sowers, T.. 1988. Gravitational separation of gases and isotopes in polar ice caps. Science, 242(4886), 16751678.CrossRefGoogle ScholarPubMed
Dreyfus, G.B. and 11 others. 2007. Anomalous flow below 2700 m in the EPICA Dome C ice core detected using δ18O of atmospheric oxygen measurements. Clim. Past Discuss., 3(1), 6393.Google Scholar
EPICA Community Members. 2004. Eight glacial cycles from an Antarctic ice core. Nature, 429(6992), 623628 Google Scholar
Herron, M.M. and Langway, Jr. 1980. Firn densification: an empirical model. J. Glaciol., 25(93), 373385.Google Scholar
Landais, A. and 12 others. 2003. A tentative reconstruction of the last interglacial and glacial inception in Greenland based on new gas measurements in the Greenland Ice Core Project (GRIP) ice core. J. Geophys. Res., 108(D18), 4563. (10.1029/2002JD003147.)Google Scholar
Petit, J.R. and 18 others. 1999. Climate and atmospheric history of the past 420,000 years from the Vostok ice core, Antarctica. Nature, 399(6735), 429436.CrossRefGoogle Scholar
Petrenko, V.V., Severinghaus, J.P., Brook, E.J., Reeh, N. and Schaefer, H.. 2006. Gas records from the West Greenland ice margin covering the Last Glacial Termination: a horizontal ice core. Quat. Sci. Rev., 25(9–10), 865875.Google Scholar
Raynaud, D. and 6 others. 2005. Palaeoclimatology: the record for marine isotopic stage 11. Nature, 436(7047), 3940.CrossRefGoogle ScholarPubMed
Schwander, J. 1989. The transformation of snow to ice and the occlusion of gases. In Oeschger, H. and Langway, Jr eds. The environmental record in glaciers and ice sheets. Chichester, etc., John Wiley and Sons, 5367.Google Scholar
Siegenthaler, U., Stocker, T.F., Monnin, E., Lüthi, D., Schwander, J. and Stauffer, B.. 2005. Stable carbon cycle–climate relationship during the late Pleistocene. Science, 310(5752), 13131317.Google Scholar
Smith, H.J., Wahlen, M., Mastroianni, D. and Taylor, K.C.. 1997. The CO2 concentration of air trapped in GISP2 ice from the Last Glacial Maximum–Holocene transition. Geophys. Res. Lett., 24(1), 14.Google Scholar
Sowers, T., Bender, M. and Raynaud, D.. 1989. Elemental and isotopic composition of occluded O2 and N2 in polar ice. J. Geophys. Res., 94(D4), 51375150.CrossRefGoogle Scholar
Spahni, R. and 10 others. 2005. Atmospheric methane and nitrous oxide of the late Pleistocene from Antarctic ice cores. Science, 310(5752), 13171321.Google Scholar
Suwa, M. 2007. Chronologies for ice cores constrained by their gas records and their implications for climate history for the past 400,000 years. (PhD thesis, Princeton University.)Google Scholar
Suwa, M., von Fischer, J.C., Bender, M.L., Landais, A. and Brook, E.J.. 2006. Chronology reconstruction for the disturbed bottom section of the GISP2 and the GRIP ice cores: implications for Termination II in Greenland. J. Geophys. Res., 111(D2), D02101. (10.1029/2005JD006032.)Google Scholar
Figure 0

Fig. 1. The CH4, CO2, δ18Oatm and δDice records for the section between 3200 and 3400 m of the Vostok ice core. Black solid circles indicate where original measurements are located. Red squares indicate values at interpolated depths used for the analysis.

Figure 1

Fig. 2. (a) An x-y-z plot of CO2 vs δ18Oatm vs CH4. The green curve traces the trajectory of CO2, CH4 and δ18Oatm points between 400 and 650 kyr in the EPICA Dome C core. Closed circles colored in shades of red and blue are Vostok data from the disturbed section of the core; the color represents δDice at the depth of sampling. The black curve connects Vostok samples in order of increasing depth. (b) As for (a) but shown only for CH4 and δ18Oatm. Gray ellipses surrounding each data point are uncertainty ranges. Black letters next to closed circles indicate the depths (m) of the Vostok sample. Orange labels next to yellow dots indicate ages (kyr) of Dome C samples.

Figure 2

Fig. 3. The age–depth relationship reconstructed in this study. Horizontal black bars indicate all compatible age ranges. The ‘best-estimate ages’ are plotted as solid circles. Shades of red and blue indicate dD of ice at the same depth.

Figure 3

Fig. 4. CO2, CH4, δ18Oatm and δD records between 400 and 650 kyr. Black curves show EPICA Dome C. Blue circles indicate Vostok samples plotted on the best-estimate ages. Note that the y axis for δ18Oatm is reversed, and δDice is plotted on offset scales for EPICA Dome C and Vostok.