Introduction
Explosive volcanic eruptions, one of the natural forcings causing short-term climatic variations, inject large amounts of ash particles and gases into the atmosphere (Reference HofmannHofmann, 1987; Reference McCormick, Thompson and TrepteMcCormick and others, 1995). The ash particles can block sunlight and darken the skies visibly, resulting in reduced solar heating (Reference Cole-DaiCole-Dai, 2010). However, such effects are typically short-lived and geographically limited since the ash is rapidly removed from the local atmosphere (Reference RobockRobock, 1981, Reference Robock2000). Of the abundant gaseous emissions, sulphur compounds (mainly SO2) are the most climatologically important component. The SO2 is subsequently converted into the chemically stable H 2SO4-H2O or sulphuric acid aerosols. The net effect of aerosols is the reduction of energy receipt near the surface; therefore the most significant climatic impact of volcanic eruptions is the cooling at the surface and in the lower troposphere (Reference Cole-DaiCole-Dai, 2010). One of the best-known historical examples of the climatic impact of volcanic eruptions is ‘the year without a summer’ (1816) following the great ad 1815 Tambora eruption on an Indonesian island (Reference Cole-Dai, Ferris, Lanciki, Savarino, Baroni and ThiemensCole-Dai and others, 2009).
Ice cores from the polar regions provide perhaps the best means to evaluate the impact of past volcanism on global climate (Reference RobockRobock, 2000). With the detection and measurement of volcanic acids (Reference KarlöfKarlöf and others, 2000; Reference UdistiUdisti and others, 2000) or sulphur compounds (Reference Delmas, Kirchner, Palais and PetitDelmas and others, 1992; Reference Cole-Dai, Mosley-Thompson and ThompsonCole-Dai and others, 1997) in polar ice cores, the history of volcanic eruptions can be recovered. Many ice cores from Greenland (Reference HammerHammer, 1980; Reference ZielinskiZielinski and others, 1994, Reference Zielinski, Mayewski, Meeker, Whitlow and Twickler1996, Reference Zielinski1997;Reference ClausenClausen and others, 1997) and Antarctica (Reference Moore, Narita and MaenoMoore and others, 1991;Reference Delmas, Kirchner, Palais and PetitDelmas and others, 1992; Reference Cole-Dai, Mosley-Thompson, Wight and ThompsonCole-Dai and others, 2000; Reference Palmer, Van Ommen, Curran, Morgan, Souney and MayewskiPalmer and others, 2001; Reference Traversi, Becagli, Castellano, Migliori, Severi and UdistiTraversi and others, 2002; Reference Traufetter, Oerter, Fischer, Weller and MillerTraufetter and others, 2004; Reference CastellanoCastellano and others, 2005; Reference KurbatovKurbatov and others, 2006; Reference ZhouZhou and others, 2006; Reference RenRen and others, 2010) have been used to reconstruct the history of volcanic eruptions. These records improve the overall quality of the chronological record of global volcanism (Reference KurbatovKurbatov and others, 2006), which in turn assists with assessing the role of volcanism in radiative forcing and climate change. Furthermore, the volcanic records provide an effective way to evaluate and improve the climate models (Reference Gao, Robock and AmmannGao and others, 2008; Reference Schneider, Ammann, Otto-Bliesner and KaufmanSchneider and others, 2009), which are the primary tools to predict future climate change. However, the timing and magnitude of a particular volcanic signal is quite variable from site to site where the ice core is drilled (Reference Dai, Mosley-Thompson and ThompsonDai and others, 1991; Reference GaoGao and others, 2006), due to variations in the atmospheric transport of volcanic substances, and their deposition and preservation on large ice sheets. Therefore, the current ice-core volcanic records of various length and quality need to be augmented with records of ice cores from additional polar locations to improve our understanding of the volcanism-climate system.
Dome Argus (Dome A), located along the main glacio- logical dividing line of the East Antarctic plateau, has the highest altitude in East Antarctica (Fig. 1). Preliminary investigation has shown that the annual mean temperature (measured at 10 m below the surface) at Dome A is -58.5°C, the lowest annual mean temperature ever recorded on the surface of the Earth (Reference Hou, Li, Xiao and RenHou and others, 2007). The average snow accumulation rate during the past several decades (AD 1966-2004) at Dome A is 23mmw.e. a-1, which is similar to that at other Antarctic inland sites (Reference Hou, Li, Xiao and RenHou and others, 2007). In addition, the detailed radar survey and Antarctic climate history indicated that the subglacial Gamburtsev mountains at Dome A are probably older than 34 × 106 years and were the main centre for ice-sheet growth (Reference SunSun and others, 2009). The glaciology research suggests that the Dome A region holds high potential for ‘oldest ice’ cores (Reference Xiao, Li, Hou, Allison, Bian and RenXiao and others, 2008), and it has attracted attention from ice-core researchers. During the 21st Chinese Antarctic Research Expedition (CHINARE 21) in the 2004/05 austral summer, a 109.91 m shallow ice core (DA2005 ice core) was recovered at a site ~300m from the summit of Dome A (80°22′ S, 77°22′ E;4092.5 m a.s.l. (Reference Zhang, Wang, Zhou and ShenZhang and others, 2007)) (Fig. 1). This is the first ice core retrieved from the Dome A summit region.
Location of ice-core sites in Antarctica referred to in the text.
We present a new regionally representative record of volcanic eruptions over the last 2840 years from the DA2005 ice core using a methodology similar to that used in previous studies (Reference Cole-Dai, Mosley-Thompson, Wight and ThompsonCole-Dai and others, 2000;Reference Budner, Cole-Dai, Robock and OppenheimerBudner and Cole-Dai, 2003). This DA2005 record is constructed using sulphate measurement of nearly 8000 samples from the top 100.42 m of the DA2005 core. Comparisons with several previously published Antarctic ice-core volcanic records are made to assess the quality of this new record. The mean accumulation rates at the Dome A region in the last three millennia are also reported.
Ice-Core Sampling and Analysis
The DA2005 ice core, drilled with an electromechanical drill, started at ~0.4m from the 2005 snow surface and reached 109.91 m depth. The bulk density of each of the 80 cm long snow/ice cylinders was measured in the field. The cylinders were then wrapped in clean plastic sheets and shipped frozen to the Polar Research Institute of China in Shanghai.
One-half (cross section) of the DA2005 core was transported to the Ice Core and Environmental Chemistry Laboratory at South Dakota State University, USA. The top 100.42 m was analyzed for major chemical impurities, and the bottom 9.49 m of the core is reserved for other analysis. The traditional discrete sampling method with stringent contamination control procedures (Reference Cole-Dai, Thompson and Mosley-ThompsonCole-Dai and others, 1995) was used to sample the porous top 1.76 m of the DA2005 ice core, with an average depth resolution of 55 mm per sample. The samples were analyzed by ion chromatography for the concentrations of major chemical impurities (Na+, NH4 +, K+, Mg2+, Ca2+, Cl-, NO3 -, SO4 2-). The technique of continuous flow analysis coupled with ion chromatography (CFA-IC) was used to analyze the core from 1.76 to 100.42 m. The CFA-IC system, as described by Reference Cole-Dai, Budner and FerrisCole-Dai and others (2006), consists of an ice-core melter, eight ion chromatographs (ICs; four for cation measurement and four for anion measurement) and an interface that distributes meltwater to the ICs.
The CFA-IC system was set up to perform one analysis per minute of all the ions in the continuous meltwater stream from the melter. In order to achieve high temporal resolution while supplying the ICs with sufficient meltwater, the ice samples were melted at relatively slow melt rates (10-22 mm min-1). Altogether, 7927 samples were analyzed to 100.42 m depth with both discrete sampling and the CFA-IC system, and the average depth resolution was 13 mm per sample. With the bulk density data of each snow/ice cylinder, a third-order polynomial was fitted to the density- depth profile for the DA2005 ice core and was used to calculate the density of each sample. Then the snow/ice depth was converted to the water equivalent depth using the calculated density.
Results
Ice-core dating and error estimates
Reference Hou, Li, Xiao and RenHou and others (2007), using field density measurement and the β-activity horizon in snow from the 1960s atmospheric nuclear tests, estimated the mean accumulation rate at Dome A during the period 1966-2004 to be ~23mmw.e.a-1. Such a low accumulation rate, one of the lowest in Antarctica, along with the sampling resolution (13 mm or ~8mmw.e. per sample) for the chemical measurement, indicated that the DA2005 core could not be dated by counting annual layers.
Another common method to date an ice core is to use a constant or average annual accumulation rate to calculate the age of each snow layer or depth in the core, after all depths are converted to water equivalent. Because of the possible variation in accumulation rate, a mean accumulation rate for a period longer than the 40 years was needed to date the DA2005 core using this method. This was accomplished by identifying prominent volcanic signals in the DA2005 non-sea-salt sulphate (nssSO4 2-) profile (Fig. 2). Using the 23mmw.e. a-1 accumulation rate by Reference Hou, Li, Xiao and RenHou and others (2007) as a guide, we found a very large nssSO4 2- signal at a depth of 17.111 m w.e. which is likely the fallout from a massive eruption by an unknown volcano in 1259 (Reference Delmas, Kirchner, Palais and PetitDelmas and others, 1992; Reference Cole-Dai, Mosley-Thompson, Wight and ThompsonCole-Dai and others, 2000). To calculate a mean annual accumulation rate, and because the DA2005 core did not begin at the 2005 snow surface, another time-stratigraphic marker was needed. Another large nssSO4 2- signal was found at 0.756 m w.e. and was identified as the fallout in 1964 of the 1963 Agung (Indonesia) eruption (Reference Delmas, Kirchner, Palais and PetitDelmas and others, 1992; Reference Cole-Dai, Mosley-Thompson, Wight and ThompsonCole-Dai and others, 2000). The computed annual accumulation rate from the depths of these two time-stratigraphic markers, 23.2 mm w.e. a-1, is similar to that determined by Reference Hou, Li, Xiao and RenHou and others (2007) for 1966-2004 and indicates that, at Dome A, annual accumulation rate averaged over a period of at least 40 years is relatively constant. This mean annual accumulation rate was used to estimate that snow at the top of the DA2005 core (~0.4 m from the 2005 surface) was deposited in 1998.
Continuous profile of non-sea-salt sulphate concentrations in the DA2005 ice core as a function of snow depth. The solid horizontal line indicates the nonvolcanic background, and the dashed line represents the detection threshold (background + 2σ).
We found the prominent nssSO4 2- signals of several other well-known volcanic eruptions. The depth, year of eruption and year of expected appearance (year of eruption plus 1) of each signal, and the age of the snow layer calculated by the constant accumulation rate method are given in Table 1. The good agreement between the year of expected appearance and the calculated age suggests that the computed mean annual accumulation (23.2 mm w.e.) is temporally representative and that the variability in annual snow accumulation rate, averaged over a relatively long time period, is likely quite small at Dome A. Therefore, the 23.2 mmw.e.a-1 accumulation rate is used to date the entire DA2005 core (top axis in Fig. 2). According to this timescale, the 109.91 m core covers the last 3186 years before present ( bp; present = end of AD 1998), and the analyzed 100.42 m part corresponds to the last 2840 years, from 840 BC to AD 1998.
Well-documented volcanic eruptions in the last two millennia and their calculated dates in the DA2005 record. All dates are calendar years. Calculated dates refer to the event years computed using the mean accumulation rate (23.2 mm w.e. a-1), and the AD 1259 Unknown event is used as the time reference for calculation. The Difference column represents the difference between the expected appearance date and the calculated date of a volcanic event in the core. Pluses denote the calculated date is later than the expected date; minus denotes the opposite
The data in Table 1 also provide a measure of the uncertainty of the dating. Because dating is achieved with an averaged accumulation rate, the error in the calculated age of a snow layer results when the actual accumulation rate deviates from the average rate. The differences between the expected date and the calculated date shown in the last column in Table 1 are used to represent the dating uncertainty of the DA2005 core. The largest deviation, at 11 years for the 1453 Kuwae (Vanuatu) eruption (Reference GaoGao and others, 2006), appears to suggest that the dating uncertainty for DA2005 is significantly smaller than those reported for ice cores from other Antarctic locations where the snow accumulation rates are low (e.g. Plateau Remote (Reference Cole-Dai, Mosley-Thompson, Wight and ThompsonCole-Dai and others, 2000), DT401 (Reference RenRen and others, 2010)). We offer no explanation or suggestion for this apparent difference, except to note that the estimated annual accumulation rate (23.2 mmw.e.a-1) is the lowest among these locations.
There are no known volcanic markers beyond that of Taupo, New Zealand (below a depth of 42.155 m w.e.). Assuming similar accumulation rate and similarly small variability of the rate for the deeper part of the core, dating errors are not expected to be significantly larger than those indicated in the last column of Table 1. The annual layers thin at depth as a result of ice flow. Therefore, the 23.2 mm layer thickness is expected to decrease with depth. However, the thinning is estimated to be ~0.4mma_1, or ~ 1% of the average layer thickness, at the bottom of the shallow core (110 m), relative to the ice-sheet thickness (~3000 m) at the summit of Dome A (Reference CuiCui and others, 2010). The error of age determination due to thinning (~10 years at the bottom of the core) is therefore smaller than the uncertainty caused by variations in the average accumulation rate. Therefore, no correction due to thinning is made to the timescale.
Criteria for the detection of volcanic signals
Explosive volcanic eruptions are not the only source of sulphate in Antarctic snow. The presence of SO4 2- in Antarctic snow is also linked to marine (sea-salt) inputs and dimethylsulphide (DMS) from marine biogenic emissions (Reference Prospero, Savoie, Saltzman and LarsonProspero and others, 1991; Reference Legrand and MayewskiLegrand and Mayewski, 1997). Sea-salt sulphate and the oxidation products of DMS constitute the nonvolcanic or background sulphate. The background sulphate concentration varies temporally in an ice core, and the sulphate from volcanic eruptions is superimposed on this variable background. To detect volcanic signals in an ice core, a threshold must be established to distinguish volcanic sulphate from the background. In this work, the volcanic threshold is estimated using a method similar to that described by Reference Cole-Dai, Mosley-Thompson and ThompsonCole-Dai and others (1997). However, two details are slightly different from the earlier method. First, instead of the sulphate concentrations, the nssSO4 2- concentrations, as calculated from the measured total SO4 2- and Na+ concentrations (Karkas and others, 2005), were used. Second, a minimum duration of 1 year (Reference Delmas, Kirchner, Palais and PetitDelmas and others, 1992) was required for elevated nssSO4 2- concentration to qualify for a volcanic event in this work.
The calculated DA2005 background nssSO4 2- concentration is 95.0µgkg-1 (solid line in Fig. 2), with a standard deviation of 28.3 µg kg-1. The threshold of 151.6 µgkg-1 (the background plus two standard deviations) is indicated by the dashed line in Figure 2. Altogether, 78 volcanic events are found with this threshold and the 1 year duration criterion. A complete list of the events is shown in Table 2. Usually, the appearance of an eruption in Antarctic snow lags the date of a low-latitude eruption by 1 or 2 years (Reference Cole-Dai and Mosley-ThompsonCole-Dai and Mosley-Thompson, 1999; Reference Legrand and WagenbachLegrand and Wagenbach, 1999). Therefore, the actual eruption years may be 1 or 2 years earlier than the dates in Table 2. The events and their associated dates and duration are designated as the DA2005 volcanic record and numbered in chronological order. In the following discussion, the volcanic events are referred to by their numbers.
Volcanic events found in the DA2005 ice core. The date for a volcanic event is assigned to the year of appearance of the sulphate peak. Negative event dates represent years BC
Volcanic fluxes
Volcanic sulphate mass flux, f, of a sample is calculated by first subtracting background nssSO4 2- from the sample nssSO42- concentration and then multiplying by the sample length in water equivalent (Reference Cole-Dai and Mosley-ThompsonCole-Dai and Mosley-Thompson, 1999). The total flux for a volcanic event is the sum of the volcanic flux of all samples associated with that event. All the volcanic fluxes are shown in Table 2.
The rate of deposition, i.e. flux, depends on several local factors at the ice-core site, such as surface irregularity, elevation, temperature, wind redistribution and relative contribution of wet-dry deposition, so the volcanic flux of a particular eruption is quite variable from site to site (Reference Clausen and HammerClausen and Hammer, 1988; Reference GaoGao and others, 2006). In this work, the normalized flux (volcanic flux normalized against that of the ad 1815 Tambora eruption, as described by Reference Cole-Dai, Mosley-Thompson and ThompsonCole-Dai and others (1997)) which may minimize the location- specific effects, is used to compare the magnitude of a volcanic event found in different ice cores. The volcanic events are categorized into three groups according to their f/f T values (Table 2), where f T is the volcanic sulphate mass flux of the AD 1815 Tambora eruption. Large eruptions (L) are those with f/f T ≥ 1, moderate eruptions (M) are those with 0.5 ≤ f/f T < 1, and small eruptions (S) are those with f/f T < 0.5. As seen in Table 2, a total of 12 large events, possibly explosive eruptions with global climatic implications, are recorded in the DA2005 ice core. It is worth noting that volcanic flux of an eruption in polar snow is related to the latitude location of the source volcano (Reference Langway, Clausen and HammerLangway and others, 1988). A small or moderate eruption in the Antarctic and sub-Antarctic region may result in a large event in the DA2005 ice core. However, since only a few volcanoes are known to be active in the Antarctic and sub-Antarctic regions, few such events are expected and they may be differentiated by comparison with other Antarctic ice-core records (Reference Cole-Dai, Mosley-Thompson, Wight and ThompsonCole-Dai and others, 2000).
Discussion
Comparing volcanic events recorded in ice cores from different sites helps to improve ice-core dating and also to remove spurious sulphate signals arising from atmospheric and glaciological effects which may be locally important, but unrelated to volcanic aerosols (Reference Cole-Dai, Mosley-Thompson, Wight and ThompsonCole-Dai and others, 2000). Depending on the availability of existing well-dated volcanic records for comparison, the following discussion of the DA2005 record is divided into three time periods: (1)the last 1000 years ( ad 1000-1998), (2) the period AD 1-1000 and (3) the period 840-1 bc. Comparisons with several Antarctic ice-core volcanic records are made to corroborate the volcanic events in the last 1000 years of the DA2005 record. For the period AD 1-1000, in which the Taupo eruption at 181 is the only known stratigraphic horizon, the DA2005 record is compared with four Antarctic records covering this time period. Comparisons are also made with Greenland records to identify possible low-latitude eruptions. The lack of well-dated stratigraphic horizons during the period 840-1 bc makes it difficult to compare specific events between the DA2005 record and other Antarctic and Greenland ice-core records. Therefore, the discussion of this period focuses on the largest event (DA69).
The last 1000 years ( ad 1000-1998)
Twenty-eight volcanic events are detected in the DA2005 ice core during the period AD 1000-1998. Table 3 lists all volcanic events in the last 1000 years recorded in the DA2005 ice core, and those in six other Antarctic ice cores covering this time period: the PR core from Plateau Remote (Reference Cole-Dai, Mosley-Thompson, Wight and ThompsonCole-Dai and others, 2000), the EDC96 core from Dome C (Reference CastellanoCastellano and others, 2005), the DT401 core from the East Antarctica plateau (Reference RenRen and others, 2010), the SP2001 core from South Pole (Reference Budner, Cole-Dai, Robock and OppenheimerBudner and Cole-Dai, 2003), the DML05 core from Dronning Maud Land (Reference Traufetter, Oerter, Fischer, Weller and MillerTraufetter and others, 2004) and the NBY89 core from Byrd Station (Reference Langway, Osada, Clausen, Hammer, Shoji and MitaniLangway and others, 1994). As seen in Table 3, most of the volcanic events found in the other Antarctic ice-core records are also detected in the DA2005 record. These include the well-known volcanic events in the last millennium: Agung (1963);Krakatau, Indonesia (1883); Cosiguina, Nicaragua (1835); Tambora (1815); an unknown eruption (1809); Unknown (1693); Mount Parker, Philippines (1641);Deception Island, Antarctica (1641);Huaynaputina, Peru (1600); Kuwae (1453); and Unknown (1259). A signal with its nssSO4 2- concentration above the threshold was dated at 1994 which is around the expected appearance date of the 1991 Pinatubo (Philippines) eruption. However, this signal does not satisfy the 1 year duration criterion for the volcanic event and is not included in the list. The 1641 eruption of Mount Parker (VEI = 6) and a contemporaneous sub-Antarctic volcanic eruption on Deception Island are found as a single sulphate event (DA12), similar to event PR9 in the PR record. A doublet, identified as Krakatau (1883; VEI = 6) and Tarawera, New Zealand (1886; VEI =5), was found in the DML05 record (Trauffetter and others, 2004), but only one signal is detected in the DA2005 record. As Reference Cole-Dai, Mosley-Thompson, Wight and ThompsonCole-Dai and others (2000) stated, it is not unusual for two volcanic signals within a few years to appear as a continuous event in ice cores from low-accumulation sites.
Volcanic events during the last 1000 years found in the DA2005, DML05, EDC96, PR, SP2001, DT401 and NBY89 ice cores. Dates are eruption years ( ad) given in each core. Events in each core are numbered sequentially
The DA2005 record also contains a number (four) of large and moderately large eruptions in the 13th century, with the 1259 event (f/f T = 3.55) being the most outstanding (Table 2). An event that appears around 1230 in several other records (DML05-27, EDC96-17, PR18, SP2001-24, DT401-16 and NBY89-23) as a relatively large event is not present in the DA2005 record. Reference Cole-Dai, Mosley-Thompson, Wight and ThompsonCole-Dai and others (2000) found by bipolar comparison that the 1230, 1259 (DA23) and 1287 (DA20) events are likely from large low-latitude eruptions. The 1230 event may not be detectable at Dome A or other locations due to local features, such as surface snow redistribution, accumulation rates, and frequency of snowfalls.
In addition to the above well-dated volcanic events, six other events (DA8, DA10, DA13, DA18, DA19 and DA24) are also in agreement with the other Antarctic ice-core records. The small event DA10, dated around 1678 in the DA2005 record, was assigned to the 1673 Gamkonora (Indonesia) eruption with VEI = 5 (Reference Traufetter, Oerter, Fischer, Weller and MillerTraufetter and others, 2004; Reference CastellanoCastellano and others, 2005). Event DA19 is detected in all the above records and assigned to the 1325 ± 75 Cerro Bravo (Colombia) eruption by Reference Traufetter, Oerter, Fischer, Weller and MillerTraufetter and others (2004). The moderate event DA24 (1194) and the small event DA8 (1764) may be signals of eruptions in the mid- to high latitudes of the Southern Hemisphere, since corresponding events were not found in either the Dye 3 or the Greenland Icecore Project (GRIP) record (Reference ClausenClausen and others, 1997). Events dA13 (1623) and DA18 (1388) are small and likely from volcanoes in the mid- to high southern latitudes, as corresponding signals are not found in any Greenland cores (Reference ZielinskiZielinski and others, 1994; Reference ClausenClausen and others, 1997). Four other small events (DA2, DA11, DA16 and DA25) have a corresponding peak in only one of the other six Antarctic records. Event DA2 (1982) is contemporaneous with a similarly small event in the DML05 record (DML05-2), and they may be the result of the El Chichon (Mexico) eruption (Reference Traufetter, Oerter, Fischer, Weller and MillerTraufetter and others, 2004). Event DA11 (1658) confirms that the small signal PR8 (1653) in the PR record (Reference Cole-Dai, Mosley-Thompson, Wight and ThompsonCole-Dai and others, 2000) is likely from a very minor volcanic eruption. And event DA16 (1511) may suggest that the small event SP2001-14 (1508) in the SP2001 record (Budner and Cole-Dai, 2003) is not a false positive detection but from a minor eruption.
Four events (DA1, DA26, DA27 and DA28) in the DA2005 record are not found in any of the six other Antarctic records. The small event DA1 is likely a spurious signal due to contamination from discrete sampling. As seen in Table 3, almost no contemporaneous events are found among the seven records during the period AD 1000-1100 when events DA26, DA27 and DA28 were recorded. This may be due to the dating errors of each core during this time period and may also indicate these Antarctic ice cores record different volcanic eruptions in the mid- to high southern latitudes. Also the possibility cannot be excluded that some of these signals are spurious.
Several events found in the DML05, EDC96, SP2001, DT401 and NBY89 records are not detected in the DA2005 record. Most of these events are detected in only one of the above records, with no corresponding signals in the other records. For example, no volcanic signal around 1969 is found in the EDC96, PR, SP2001, DT401 and NBY89 records to support event DML05-3 in the DML05 record. Some of these signals are so small that they may be spurious signals or easily missed in ice cores drilled at sites with quite low accumulation rates such as Dome A. Although no event around 1172 is recorded in DA2005, an event was detected in the DML05, EDC96, SP2001 and NBY89 records and was also recorded in another Antarctic ice core, the SP78 core (Reference Langway, Osada, Clausen, Hammer and ShojiLangway and others, 1995). In the Greenland Ice Sheet Project 2 (GISP2; Reference ZielinskiZielinski and others, 1994), GRIP and Dye 3 (Reference ClausenClausen and others, 1997) records, an event was found at 1175, 1179 and 1180, respectively, and was attributed to an Icelandic eruption. Therefore, the event around 1172 may represent coincident eruptions in the mid- to high latitudes of both hemispheres.
AD 1-1000 (2000-1000 years BP)
The PR, DT401 and EDC96 records mentioned above also cover this time period, while the DML05 record (dated by counting annual layers) only presents the volcanic history as old as AD 186. Figure 3 compares the DA2005 record for the first millennium with those in the PR, DT401, EDC96 and DML05 records. Dates and event numbers of contemporaneous events found in the DA2005 record and at least two of the other four Antarctic records are listed in Table 4. We also found several events in the DA2005 record with potential counterparts in only one of the other four records. These events are not included in Table 4. Unlike the numerous well-dated volcanic eruptions in the last 1000 years, the only well-known event during this 1000year period is the Taupo eruption at 181. Thus, a detailed discussion of all the contemporaneous volcanic events is not warranted for this 1000 year period and we summarize below the most prominent events.
Comparison of volcanic profiles (volcanic flux versus age) for the first millennium AD (AD 1–1000) from sulphate measurement in the PR, DT401, EDC96, DML05 and DA2005 ice cores. Shaded area indicates period not covered by ice core.
Contemporaneous events during the first millennium AD (AD 1-1000) found in DA2005 and at least two of the DML05, EDC96, PR and DT401 ice cores (all from East Antarctica)
As seen in Table 4, ten events in the DA2005 record (a total of 29 events) have corresponding signals in at least two of the other four records. The largest age difference of corresponding events is 30 years between EDC96-29 and DA51. However, all differences are well within the dating errors of these cores. Of the above ten events, five are large (DA38, DA40, DA41, DA51 and DA52) but none is of the 1259 Unknown event magnitude, and another four (DA45, DA48, DA49 and DA56) are moderate events. It is worth noting that there is good agreement among volcanic signatures in the period AD 180-700 at three different sites (Dome A, Dome C and DML) far away from each other (Fig. 1). Indeed, nine of fourteen events in the DML05 record and eight of nine in the EDC96 record are contemporaneous to the major events in the DA2005 record during the period AD 180-700.
The Taupo eruption at AD 186± 10 (Reference Wilson, Ambraseys, Bradley and WalkerWilson and others, 1980) is regarded as one of the most significant eruptions (VEI=6+) in the Southern Hemisphere in the first millennium (Reference Cole-Dai, Mosley-Thompson, Wight and ThompsonCole-Dai and others, 2000). It has been detected in several Antarctic and Greenland ice cores (Reference ZielinskiZielinski and others, 1994; Reference Cole-Dai, Mosley-Thompson, Wight and ThompsonCole-Dai and others, 2000; Reference Traufetter, Oerter, Fischer, Weller and MillerTraufetter and others, 2004; Reference RenRen and others, 2010). As seen in Figure 3, the Taupo event was a very large signal in both the PR and DT401 cores. However, in the DML05 core (Reference Traufetter, Oerter, Fischer, Weller and MillerTraufetter and others, 2004) and in an ice core from Siple Dome, West Antarctica (SDMA; Reference KurbatovKurbatov and others, 2006), Taupo was either a moderate or small signal. The large signal dated at 181 (event DA52) in DA2005 may correspond to the Taupo eruption, with f/f T = 1.27 compared to 3.10 in the PR record (Reference Cole-Dai, Mosley-Thompson, Wight and ThompsonCole-Dai and others, 2000). According to Reference LarsenLarsen and others (2008), three contemporaneous events (dated at 674/675, 567/568 and 533/534, respectively) in three Greenland ice- core (Dye 3, GRIPand NorthGRIP) records correspond to the DML deposits at 685 (DML05-36), 578 (DML05-39) and 542 (DML05-40), respectively. With the DML05 record as the reference, DA38 (698), DA40 (588) and DA41 (551) may be the corresponding events to the three events in the Greenland records. Hence, events DA38, DA40 and DA41 are most likely low-latitude eruptions. The moderate event DA56 dated at 68 appears to correspond to a large event (PR32 dated at 74) in the PR record (Reference Cole-Dai, Mosley-Thompson, Wight and ThompsonCole-Dai and others, 2000) and a small event (DT401-28 dated at 82) in the DT401 record (Reference RenRen and others, 2010). In the GISP2 (Reference ZielinskiZielinski and others, 1994) and GRIP (Reference ClausenClausen and others, 1997) records, an event was also found at AD 77 and 79, respectively. However, the event in the above three Antarctic records and that in the two Greenland records are unlikely from the same eruption since the AD 77 event may have resulted from the AD 79 Mount Vesuvius (Italy) eruption in the Northern Hemisphere (Reference ZielinskiZielinski and others, 1994).
These comparisons with the other Antarctic volcanic records suggest that the DA2005 record during this 1000year period is reliable. More ice-core-based volcanic records are needed to clarify the source and magnitude of the contemporaneous events found in DA2005 and previous Antarctic ice cores.
840-1 BC (2840-2000 years BP)
A total of 21 events are detected in the DA2005 ice core during the period 840-1 BC, compared to 8 in the PR core (Reference Cole-Dai, Mosley-Thompson, Wight and ThompsonCole-Dai and others, 2000) and 5 in the EDC96 core (Reference CastellanoCastellano and others, 2005). Eruption frequency (2.8 per century) for the most recent 2000 years in DA2005 is also higher than those in the Antarctic ice cores from sites with low accumulation rates (e.g. the PR, EDC96 and DT401 cores). These differences may be due to a number of glaciological and atmospheric factors variable from site to site. Another influencing factor may be the temporal resolution of the sulphate measurement. Lower temporal resolution (larger sample length) could dilute the sulphate concentration of a volcanic event such that the event is below the detection threshold. Therefore, some small events may not be detected in ice cores with low-resolution analysis. For example, the DA2005 and DT401 cores are from locations only 120 km from each other in East Antarctica (Fig. 1), which may share similar glaciological and atmospheric features. The fewer volcanic signatures in DT401 may also be due to the lower temporal resolution of the analysis of the DT401 core (one measurement every 30 or 35 mm in DT401 compared to 13 mm in DA2005).
The largest volcanic event in the DA2005 record (DA69) is dated at 425 bc. The maximum sulphate concentration of DA69 (1009.8 mgkg-1) is similar to that of the 1259 Unknown eruption (1011.0 mgkg-1), but the volcanic flux (156.74 kg km-2) of DA69 is more than twice that of the 1259 event (63.27 kg km-2). The duration of DA69 (~12.3 years) is significantly longer than the typical duration of 1-3 years expected for volcanic events in Antarctic ice cores. The long duration may be attributed to the combined effect of a long atmospheric residence time of its volcanic aerosols and postdeposition modification (Reference Cole-Dai, Mosley-Thompson, Wight and ThompsonCole-Dai and others, 2000). No such volcanic signals were found in the PR and DT401 records during this time period. However, in the EDC96 record, a large signal was detected at 384 bc (Reference CastellanoCastellano and others, 2005). Also in the SDMA record a very large event was found at 325 BC (Reference KurbatovKurbatov and others, 2006). It is not possible, given that the differences in the ages are within the dating uncertainties of these cores, to determine whether these three signals are from the same eruption.
Conclusions
A total of 78 volcanic eruptions are identified in the top 100.42 m of the DA2005 ice core covering the period 840 BC-AD 1998. Of these, 12 are probable large events with fluxes larger than that of the ad 1815 Tambora eruption (VEI = 7). The largest event is dated at 425 bc, with its sulphate flux almost nine times that of the Tambora eruption. The mean eruption frequency for the entire 2840 years (2.7 events per century) in DA2005 is higher than those in previous Antarctic ice cores recovered from sites with low accumulation rates.
Comparisons with previous Antarctic ice-core volcanic records suggest that the DA2005 record is reliable in the most recent 2000 years. A series of volcanic events in the first millennium AD (AD 1-1000) in the DA2005 record are also found in several other Antarctic ice-core records. Three events, DA38, DA40 and DA41, during the period AD 1-1000 are likely from low-latitude eruptions since corresponding signals were found in several Greenland ice-core volcanic records. However, better-dated Antarctic ice-core volcanic records, ideally from sites with high accumulation rates, are needed to corroborate the DA2005 events in the older (840-1 BC) period.
The mean accumulation rate between the 1963 Agung and 1259 Unknown eruptions is 23.2 mm w.e. a-1, which is quite similar to the value between 1966 and 2004 measured by Reference Hou, Li, Xiao and RenHou and others (2007). Dating results from the 700year mean accumulation rate show that there is good agreement between the year of expected appearance and the calculated age of all the well-dated volcanic events in the last two millennia. This suggests that there is neither an indication of a change nor a trend in the accumulation rate apparent in the last 2000 years, which may indicate that no drastic change in deposition has occurred at Dome A within this time period.
Acknowledgements
We thank all the members of the CHINARE 21 Inland Traverse (2004/05) team. Special thanks go to Shugui Hou and Cunde Xiao for their contribution in collecting the ice- core samples. This research was financially supported by the Natural Science Foundation of China (40906098, 40773074, 40806074) and Ministry of Science and Technology of China (2006BAB18B01). Ice-core analysis at South Dakota State University was supported in part by the US National Science Foundation.
