The interactions between climate changes and environmental variations are complex and their cause/effect relationship is not yet fully understood. Indeed, if the large glacial/interglacial changes in the Earth climate system are attributed to orbital driving forces (see, e.g., Reference PetitPetit and others, 1999), the role played by environmental changes (e.g. sea level, ice cover, aerosol production) in fast climate changes through feedback processes is still an open question. Polar ice cores constitute a conservative and reliable environmental archive storing the gas and aerosol composition of the atmosphere for hundreds of thousands of years, and ice cores are widely used to reconstruct climatic and environmental changes over several glacial/interglacial cycles (Reference PetitPetit and others, 1999). While stable-isotopic profiles (δD or δ18O, deuterium excess) provide a reliable reconstruction of climatic conditions (temperature and relative humidity) in the deposition-site and moisture-source areas (Reference JouzelJouzel and others, 1995; Reference StenniStenni and others, 2001; Reference Vimeux, Cuffey and JouzelVimeux and others, 2002) and can be used to evaluate accumulation rates (Reference Schwander, Jouzel, Hammer, Petit, Udisti and E.WolffSchwander and others, 2001), selected chemical compounds have been proposed as actual or potential markers of environmental changes. For instance, atmospheric load of dust and sea-salt components, mainly related to atmospheric circulation processes, were inferred from the Ca2+ and Na+ content (Reference Mayewski and LegrandLegrand and others, 1988; Reference SteigSteig and others, 2000; Reference RöthlisbergerRöthlisberger and others, 2002a; Reference Ruth, Wagenbach, Bigler, Steffensen, Röthlisberger and MillerRuth and others, 2002), marine biogenic activity from methanesulphonic acid (MSA) and non-sea-salt sulphate (nssSO4 2–) (Reference McCulloch and DaviesLegrand and others, 1991, Reference Traversi, Becagli, Castellano, Largiuni, Udisti, Colacino and Giovannelli1997; Reference Hansson and SaltzmanHansson and Saltzman, 1993), continental biogenic changes from ammonium (Reference Dentener and CrutzenDentener and Crutzen, 1994; Reference McCulloch and DaviesLegrand and others, 1998), N-cycle compounds transformation and transport processes from nitrate (Reference Mayewski and LegrandMayewski and Legrand, 1990;Reference Wolff and DelmasWolff, 1995) and volcanic activity from non-biogenic sulphuric acid (Reference Delmas, Kirchner, Palais and PetitDelmas and others, 1992; Reference Cole-Dai, Mosley-Thompson, Wight and ThompsonCole-Dai and others, 2000; Reference UdistiUdisti and others, 2000; Reference ZielinskiZielinski, 2000).
In this paper, we compare the high-resolution chemical profiles of several ionic components (main and trace anions and cations) with δD and electrical conductivity measurement (ECM) profiles, all measured along the first EPICA (European Project for Ice Coring in Antarctica) Dome C ice core (EDC96). The chemical profiles constitute a continuous high-resolution record spanning the last 45 kyr, with a temporal resolution of 1–9 years depending on the depth. The EDC96 ice core, 788m long, was drilled in the field seasons from 1996/97 to 1998/99 at Dome C, central East Antarctica (75˚ 06’ S, 123˚ 24’ E; 3233ma.s.l., about 1100 km from the coastline) in the framework of the EPICA programme. Ionic components (Na+, NH4 +, K+, Mg2+, Ca2+, Cl–, NOO3 –, S4 2–, MSA) were measured by ion chromatography (IC) within the EPICA chemistry consortium, including laboratories from five European countries.
While the environmental and climatic significance of concentration changes of single components or component groups will be discussed in detail in specific papers, the aim of this paper is to give a first, preliminary overview of the pattern of all the analyzed compounds under the different climatic conditions spanning the last 45 kyr: the Holocene, the Holocene Optimum (HO), the last glacial/interglacial transition, the Last Glacial Maximum (LGM) and the Late-glacial period. Particular attention has been paid to evaluating the sensitivity of the potential environmental markers to the climatic changes that occurred during the last transition, and especially during the Antarctic Cold Reversal (ACR; a mid-deglaciation cold episode indicated by a sharp decrease in the δD record starting from about 14.0 kyr BP; see section 3), by observing trends in their concentrations and fluxes. Moreover, the paper aims to check the reliability of concentration and flux profiles to provide information on changes in source intensity and transport efficiency of atmospheric aerosol and in snow accumulation rate.
Data reported here come from IC analyses performed in different European laboratories on the entire EDC96 ice core. The core was drilled at Dome C as part of the EPICA programme aiming to reconstruct high-resolution isotopic, chemical and physical profiles for the last 1×106 years (EPICA Dome C 2001–02 Science and Drilling Teams, 2002). The first drilling attempt started during the 1996/97 campaign and stopped in 1998/99 because the drill got stuck at 788m depth. This ice core was named EDC96. A second drilling, EDC99, started in the 1999/2000 field season and reached 3100m depth, about 100m above the bedrock, during the 2002/03 campaign. The two cores were analyzed in the field for electrical properties (ECM and dielectric profiling (DEP):Reference Wolff, Wolff and BalesWolff and others, 1999; Reference Vimeux, Cuffey and JouzelUdisti and others, 2000). While EDC99 chemical and isotopic analysis is still in progress, EDC96 measurements are almost complete along the entire core. IC analysis was performed on continuous ice-core sections in the framework of the EPICA chemical consortium, involving laboratories from Denmark, France, Italy, Sweden and the United Kingdom. Ice-core sections were decontaminated following two different approaches:
1. subsamples were cut by a stainless-steel saw and packed; in the European laboratories, ice pieces were decontaminated by mechanically removing (generally with an electric plane or by manual scraper) some millimetres of the external layer;
2. sections about 110cm long (cross-section about 32 mm ×34 mm) were continuously melted as described in Reference PetitRöthlisberger and others (2000a). The inner contaminant-free part of each bar was collected in pre-cleaned Coulter Counter-type accuvettes, subdividing the melted flux into subsamples with a sequential autosampler.
Ice-core cutting was performed at Dome C, while accuvette subsampling on the last 200 m was carried out at the Alfred Wegener Institute, Bremerhaven, Germany, in October 2000.
Subsamples have a resolution of 2.5–10cm (corresponding to 1.0–4.0 and 2.2–9.0 years in the Holocene and in the glacial period, respectively), with the highest resolution in sections recording fast climatic changes (such as the transition, including the ACR).
Each laboratory received one in every five 55 cm sections, so that the accuracy of the measurements carried out by the different laboratories was continuously evaluated by comparison with the values measured on the contiguous sections. Inter-laboratory calibration rounds were periodically carried out to check the reliability of IC analysis and to tune analytical procedures and responses to obtain self-consistent profiles (Reference LittotLittot and others, 2002). Inter-laboratory comparison carried out on standards at ppb levels and on Holocene samples, where the lowest concentrations were measured, yielded errors of < 15%. Minimum errors (2.5%) were determined for sulphate, and maximum values (15%) for MSA and Mg. Very low Holocene levels of Ca, K and nitrate showed higher uncertainties (up to 50%), but these include errors in manipulation and transport of melted samples used for the inter-calibration exercise, while routine samples were decontaminated and analyzed from each laboratory just before the analysis. Glacial values were usually one to two orders of magnitude higher than the Holocene ones, so even the highest errors do not significantly affect the glacial-to-Holocene trends discussed here.
Firn sections showed high ammonium contamination by uptake of ammonia. In fact, the samples were stored in the Dome C buffer for > 1 year (two field seasons), not sealed in polyethylene bags, because the analysis processing line was set up the year after drilling started. Ammonium data related to the first 100 m of the EDC96 ice core were discarded for the Holocene mean calculation.
ECM was performed on site (Reference Wolff, Wolff and BalesWolff and others, 1999; Reference Vimeux, Cuffey and JouzelUdisti and others, 2000). Isotopic δD measurements were carried out at Saclay, France, in the framework of the EPICA isotopes consortium (Reference JouzelJouzel and others, 2001).
Since the drilling site is located on the top of a dome, EDC96 dating (EDC1) was achieved by a simple flow model considering only vertical deformations and using accumulation-rate values calculated from the relationship between δD isotopic ratio and relative humidity; model parameters were calibrated by using volcanic signatures, the end of the Greenland Younger Dryas (YD; 11.53 kyr BP) and the 10Be peak, occurring at 41 kyr BP (BP = before AD 1950), as temporal horizons (Reference Schwander, Jouzel, Hammer, Petit, Udisti and E.WolffSchwander and others, 2001). The depth at EDC96 corresponding to the end of the YD has been estimated by isotope match (EDC96/Byrd: Reference Hammer, Clausen and DansgaardHammer and others, 1994) and CH4 match (Byrd/Greenland Icecore Project (GRIP): Reference BlunierBlunier and others, 1998).
3. Results and Discussion
Figure 1 shows the concentration/age profiles for all measured chemical species and ECM, contrasted to the δD profile for the last 45 kyr. On the basis of δD trends and the EDC1 time-scale, we divided the temporal range into five periods:
1. Holocene (without Climatic Optimum): 0–9.5 kyr BP; 6–300 m;
2. HO: 9.5–11.5 kyr BP; 300–360 m;
3. Transition: 11.5–18.0 kyr BP; 360–480 m;
4. LGM: 18.0–32.0 kyr BP; 480–640 m;
5. Glacial (part): 32.0–45.0 kyr BP; 640–780 m.
Vertical lines in Figure 1 separate the different periods.
During the glacial/interglacial transition, particular attention was paid to the ACR period (see section 3.3) and to the following warming phase until the beginning of the Holocene (about 11.5 kyr BP). The mid-deglaciation cold episode (ACR) is very evident in the δD record as a sharp decrease starting from about 14.0 kyr BP The δD profile reaches a relative minimum value (end of ACR) at about 12.5 kyr BP. ACR climatic features and the temporal shift with respect to the Greenland YD have been discussed in several papers
As already pointed out by several authors (Reference Legrand, Ducroz, Wagenbach, Mulvaney and HallLegrand and others, 1988; Reference Watanabe, Kamiyama, Motoyama, Fujii, Shoji and SatowWatanabe and others, 1999; Reference SteigSteig and others, 2000), background concentration values of chemical compounds show trends anticorrelated to the temperature, which is derived from the δD isotopic profile. We mainly refer to the warm events occurring around 42.5 and 38.0 kyr BP (see δD peaks) and the following colder periods, to the LGM, to the transition and to the HO. The 38.0 kyr warm period was also found in the Byrd and Vostok ice cores and named A1 (Reference BlunierBlunier and others, 1998). Another older warm period, probably corresponding to the A2 warming event (Reference BlunierBlunier and others, 1998), is partially recorded at the end of the EDC96 ice core.
Some components, such as Ca2+, Cl–, NO3 –, Na+, Mg2+ and MSA, show high sensitivity to climatic changes (Fig. 1). Mean concentration values in the different climatic periods are reported in Figure 2. The highest concentrations were recorded in the cold periods, and especially in the LGM. Minimum background concentrations are observed in the warm stages of the glacial and in the Holocene. On the other hand, sulphate, ammonium, fluoride and, to a lesser extent, potassium, show minor changes in the different climatic periods, with little variation also during the glacial/interglacial transition. As for ammonium, the mean concentration was calculated excluding the first 100 m, affected by contamination (see section 2).
Dome C past accumulation rates were calculated by supposing them to be proportional to the derivative of the mean saturation vapour pressure at the inversion layer with respect to temperature (Reference JouzelJouzel and others, 1987), where mean surface temperature was estimated from the stable-isotope ratio in the ice (Reference Schwander, Jouzel, Hammer, Petit, Udisti and E.WolffSchwander and others, 2001). In this way, we reconstructed the accumulation-rate curve with a depth resolution of 0.55 m (corresponding to about 20 years in the Holocene and 50 years in the glacial). Since accumulation rate increases around 2.3 times from the LGM to the Holocene, we have to consider the effects of such changes on the snow content of the chemical compounds. On the basis of the EDC1 dating and taking into account the ice-layer thinning with depth, we associated to each concentration (C) the accumulation rate (A) evaluated for the time period covered by the sample resolution, calculating the net deposition fluxes (F):
Both the use of concentration and the use of flux to study the temporal profiles of climatic and environmental markers along an ice core have advantages and disadvantages. Concentration is an independent chemical measurement and thus more reliable, but changes in concentration are dependent on variation in source intensity and transport efficiency as well as in accumulation rate, especially for sites where dry deposition is important or dominant (low-accumulation-rate areas, such as Dome C: Reference LittotLegrand, 1987). Flux (net deposition mass) is calculated by multiplying concentration by accumulation rate and is potentially affected by larger errors, mainly related to the reliable reconstruction of past accumulation rate from the isotopic profile. On the other hand, flux is virtually independent of accumulation rate at Dome C, dominated by dry-deposition processes, and therefore more suitable for assessing changes in atmospheric load. We use both parameters, preferring fluxes when studying the effects of source and transport changes, and concentrations when describing the comprehensive sensitivity (driven by accumulation rate, source and transport changes) of chemical profiles to climate variations.
3.1. Mean values
Figure 2 shows mean concentration and flux values for the most relevant components measured in the EDC96 ice core in the five different climatic stages.
The mean concentration general trends confirm the profile observations. The highest mean concentrations are measured in the LGM, when both temperature and accumulation rate were lower, and minimum values are recorded in the Holocene. This pattern is shown by Na+ (Mg2+ and K+ follow identical trends) and especially by Ca2+, which exhibits the highest concentration change from high LGM values to minimum values during the HO. The parallelism between concentration and flux trends for Na+ and Ca2+ demonstrates that accumulation rate tunes only the net deposition quantities but temporal changes are driven by sources and/or transport changes. By contrast, the variations observed in sulphate concentration trend practically disappear in the flux pattern, showing the sulphate deposition processes are almost completely dominated by changes in accumulation rates. A small maximum is recorded during the HO, probably related to higher biogenic activity, as seems to be suggested by the contemporaneous increase in ammonium and MSA fluxes.
Chloride and nitrate show a more complex pattern, mainly because they are not conservative species (Reference Wolff, Wolff and BalesWolff, 1996). Indeed, post-depositional processes (diffusion, redistribution, re-emission into the atmosphere and photochemical reactions) can alter the original snow composition at low-accumulation sites (Reference Wagnon, Delmas and LegrandWagnon and others, 1999; Reference Röthlisberger, Hutterli, Sommer, Wolff and MulvaneyRöthlisberger and others, 2000a, Reference Röthlisberger, Hutterli, Sommer, Wolff and Mulvaney2002 b; Reference UdistiTraversi and others, 2000). The present Dome C accumulation rate of 2.7 cm w.e. a–1 was reduced to about 50% in the LGM (Reference Schwander, Jouzel, Hammer, Petit, Udisti and E.WolffSchwander and others, 2001), so post-depositional processes heavily affect the chloride and nitrate preservation in snow at this site. As one of the main post-depositional processes is related to reemission of the volatile acidic species (HNO3 and HCl) into the atmosphere, alkaline dust was able to fix higher concentrations of chloride and nitrate as non-volatile salts during the glacial period (Reference PetitRöthlisberger and others, 2000a, Reference Petit2002b). Therefore, chloride and nitrate concentrations and fluxes are highly dependent on changes in both accumulation rate and snow acidity. Indeed, these two components show relatively high fluxes during the HO, when snow accumulation rate was so high as to preserve larger quantities, and in the LGM, when maximum dust atmospheric load (Reference Delmonte, Petit and MaggiDelmonte and others, 2002a, Reference Delmonte, Petit and Maggib) fixed higher chloride and nitrate depositions as non-volatile salts. The better preservation of HCl deposition in the HO is also demonstrated by the Cl–/Na+ ratio which shows, in this period, a maximum mean value (2.2 w/w) in the last 45 kyr, significantly higher than the sea-water ratio (1.8 w/w). The pattern in the HO is similar to that measured in the present uppermost superficial layers at Dome C, showing the original contemporaneous deposition of NaCl and HCl was preserved by a relatively high accumulation rate (Reference Traversi, Becagli, Castellano, Largiuni, Udisti, Colacino and GiovannelliTraversi and others, 2000). In the glacial period, before the LGM, covered by the EDC96 ice core (32.0–45.0 kyr BP), the dust atmospheric load was much lower than in the LGM (Reference PetitPetit and others, 1999), and the accumulation rate similar (Reference Schwander, Jouzel, Hammer, Petit, Udisti and E.WolffSchwander and others, 2001). This environmental situation was still suitable for preserving chloride depositions, but the low snow alkalinity was probably insufficient to compensate the negative effect of low accumulation on re-emission of HNO3. Indeed, nitrate mean flux shows the lowest values in the entire core. The higher sensitivity of nitrate to post-depositional re-emission into the atmosphere is also demonstrated by the pattern of present snow deposition at Dome C: nitrate concentration (and flux) quickly decreases from high superficial values (around 150 ppb) to very low background Holocene values (around a few ppb) in the first 100 cm of snow layers (Reference Traversi, Becagli, Castellano, Largiuni, Udisti, Colacino and GiovannelliTraversi and others, 2000).
MSA is also affected by post-depositional effects, consisting of migration and re-emission processes, complicating its use as a univocal marker of past biogenic activity changes (Reference Wagnon, Delmas and LegrandWagnon and others, 1999; Reference Pasteur and MulvaneyPasteur and Mulvaney, 2000; Reference Delmas, Wagnon, Goto-Azuma, Kamiyama and WatanabeDelmas and others, 2003). Therefore, the net depositional fluxes of this component are also probably directly correlated to snow accumulation rate and alkaline dust content in the snow. Indeed, the pattern of MSA fluxes in the Holocene, HO, transition and LGM are very similar to those shown by chloride and nitrate. In the rest of the glacial, by contrast, MSA shows the highest mean flux and concentration values in all the last 45 kyr. Because the
32.0–45.0 kyr BP glacial period was characterized by low alkaline-dust atmospheric load and snow accumulation rates, we postulate that the highest MSA fluxes with respect to the LGM reflect higher phytoplanktonic activity or more efficient meridional transports from oceanic area sources of biogenic emissions or a different oxidation pathway of phytoplanktonic dimethylsulphide, producing a higher MSA/H2SO4 ratio.
3.2. Flux profiles
The flux patterns observed from multi-millennial averages can be observed in detail by continuous 45 kyr profiles. Figure 3 shows the smoothed trends of the Na+, Ca2+, Cl–, NO3 – and MSA fluxes, compared with the smoothed δD profile. The smoothing procedure was made by a weighted fitting curve. This function fits a curve to the data using a locally weighted least-squares error method. The result of this curve fit is to plot a best-fit smooth curve through the centre of the data. This is a robust fitting technique, nearly insensitive to outliers, and so able to give a reliable trend for background values. The low percentage of the smoothing (2% for δD and 3% for chemical profiles) ensures a sufficient temporal resolution for the interval, around each value, where smoothing was carried out (for the chemical dataset, 3% = 600 values). The running procedure restitutes a number of values as high as the original dataset.
To indicate the general inverse relationship and the particular features of the different profiles with respect to the isotopic profile, we use an unusual presentation, plotting the flux profiles in a reversed scale. In this way, the anticorrelation between chemical and isotopic profiles during warm and cold stages and transition appears as parallel trends, indicating particular patterns in the ACR.
In the Na+ and Ca2+ flux profiles, the inverse relationship between fluxes and isotopic temperature is very evident. All climatic changes (warm and cold stages in the glacial period, LGM, transition, ACR, HO) are recorded in the Na+ profile, indicating an anticorrelation between source intensity or transport efficiency and temperature in the different climatic periods. Ca2+ shows a still higher sensitivity for the A1 stage and, especially, for the LGM, when high atmospheric dust levels increase the depositional Ca2+ fluxes. Unlike Na+, however, Ca2+ exhibits no significant change during the ACR and the entire Holocene, demonstrating a different pattern of the dust source areas, with respect to sea spray, during the transition (see further below). Na+ sensitivity with respect to the ACR climate oscillation could be related to change in atmospheric transport processes, probably driven by redistribution of the latitudinal temperature gradient. Similar trends were observed by Reference RöthlisbergerRöthlisberger and others (2002a) using non-sea-salt Ca2+ (nssCa2+) and sea-salt Na+ (ssNa+) flux profiles obtained by continuous flow analysis (CFA; Reference PetitRöthlisberger and others, 2000b). For the ssNa+ profile, however, these authors did not discuss the situation of the ACR. Here, we preferred to use total Ca2+ and Na+ flux profiles to avoid every potential artefact in the original trends. Anyway, as the Ca2+ drops to very low fluxes about 1.0 kyr before the ACR onset (see section 3.3), the crustal contribution to the Na+ budget is only few per cent and is not sufficient to affect the observed trend.
The chloride and nitrate trends during the glacial and early transition fit the Ca2+ pattern very well, confirming the role of alkaline dust in fixing them in the snow layers as non-volatile salts. In particular, the nitrate profile is nearly superimposable on that of Ca2+ during the early transition (see, for detail, Fig. 4), in the LGM and in the rest of the glacial (with some differences in the deepest part of the ice core), showing the same sensitivity to warm and cold stages. Conversely, chloride and nitrate fluxes differ from Ca2+ in the late transition and early Holocene. Here also, nitrate has the highest sensitivity to climatic changes, showing a sharp increase in the net deposition flux in phase (Figure 3 reports fluxes as reversed scale) with temperature maxima at the ACR onset (about 14.2 kyr BP) and in the HO, probably driven by the accumulation-rate increases and not by dust content (Ca2+ fluxes are very low in the same time periods). A flux peak during the HO is also visible in the chloride profile, while no significant changes are seen in the ACR. The MSA profile is characterized by high, constant fluxes during the LGM and glacial, with low sensitivity, sometimes out of phase, to isotopic temperature changes. An inverse relationship is shown during the warm stages in the glacial, especially during A1. No significant changes are visible during the ACR, while a sharp flux increase occurs in the HO, with the same relative intensity as chloride. Finally, a progressive decreasing trend of MSA fluxes is evident during the Holocene, reaching relative stability in the last 5 kyr.
3.3 Temporal dephasing during the transition
Figure 4 shows the concentration profiles of some of the measured components during the transition, compared to δD and ECM profiles. Concentrations are plotted on a reversed scale as smoothed normalized profiles. The normalization was obtained for each component by dividing each sample concentration by the highest value measured for the same components during the LGM (when all components show maximum concentrations). Here, we preferred to report concentration instead of fluxes, to avoid any possible artefact due to an incorrect evaluation of snow accumulation rate in a temporal range where climatic changes were very fast. Besides, as already noted, concentration variations could be more suitable for evaluating the overall climate sensitivity of a snow component, including the effects of accumulation-rate changes as well as source and transport variations. Also in this case, the inverse scale allows easier comparison with the isotopic profile, better indicating changes in each profile slope.
Although ECM is mainly related to snow background acidity (Reference Hammer, Clausen and LangwayHammer and others, 1980), it also depends on the ionic content so that the smoothed ECM profile shows a different pattern with respect to the other considered parameters. For the Holocene, Reference Barnes, Wolff, Udisti, Castellano, Röthlisberger and SteffensenBarnes and others (2002) found a difference between the calculated acidity from H2SO4, HCl and HNO3 and that resulting from simple models of electrical measurements; the authors hypothesize that with a larger grain-size a greater connectivity occurs with the same amount of acid, and therefore higher conductivity results. In the glacial period the issue is more complex due to the high dust content, and further investigation is needed to fully understand the electric signals.
The majority of chemical components parallel the isotopic curve up to the ACR, showing an inverse relationship between concentration and temperature. Starting from the ACR onset, concentrations tend to reach the low Holocene values with a much slower gradient, revealing a general low sensitivity for the ACR oscillation and HO conditions. Two groups of chemical components show a quite different pattern. On the one hand, sulphate and ammonium reversed profiles follow, perfectly in phase, the δD profile up to reach the low, quite constant, Holocene values about 1.6 kyr before the ACR onset, which occurred about 14.0 kyr ago (Reference Schwander, Jouzel, Hammer, Petit, Udisti and E.WolffSchwander and others, 2001). On the other hand, nitrate and Ca2+ reversed concentration profiles precede the isotopic transition onset by about 300 years, with a dephasing very similar to that shown by ECM. After the first abrupt increase, nitrate and Ca2+ reach the low Holocene values about 1.0 kyr before the ACR onset, similarly to chloride. Whereas Ca2+ and chloride show no significant changes in the second step of the deglaciation, nitrate exhibits a concentration increase whose relative maximum is in phase with the ACR onset.
The pattern of Na+ (very similar to those shown by Mg2+ and K+) is different: the inverse normalized concentration fits the δD profile very well, following the ACR oscillation. Only MSA, of the components analyzed, shows a long delay (about 400 years) with respect to the δD increase in the deglaciation onset.
Finally, the large fluoride peak during the LGM is due to the 170 year fluoride event first identified by Reference Hammer, Clausen and LangwayHammer and others (1997) in the Byrd ice core.
Granting concentrations can reliably describe the comprehensive environmental changes that occurred during the last deglaciation, the study of the high-resolution chemical trends revealed that environmental changes in the aerosol source areas preceded or were contemporaneous (with the apparent exceptions of chloride and, especially, MSA) to the Antarctic temperature changes, as evaluated by the δD isotopic profile. This probably implies that the Antarctic site temperature (as recorded at Dome C) did not lead changes in the snow deposition composition, except for variations in the snow accumulation rates.
In particular, the Ca2+ decrease before the deglaciation onset, in phase with the acidity increase (as revealed by the ECM profile), shows that the main dust-source areas for east central Antarctica (Patagonia: Reference Basile, Grousset, Revel, Petit, Biscaye and BarkovBasile and others, 1997) experienced environmental variations in advance of central Antarctic temperature changes. Reference RöthlisbergerRöthlisberger and others (2002a) suggested that the early Ca2+ decrease could be related to changes in Patagonian climatic conditions (warmer temperature and higher humidity leading to larger vegetation growth) during the last part of the LGM (Reference McCulloch and DaviesMcCulloch and Davies, 2001). Unfortunately, the oldest record reported by these authors dates back to 17.33 kyr BP, while the isotopic temperature at Dome C increases from around 18.0 kyr BP, so that we have no definite evidence the Patagonian warming started in advance. Anyway, the Ca2+ profile trend suggests that dust-source areas affecting Dome C snow deposition started to warm a few hundred years before east central Antarctic warming and probably reached environmental conditions leading to low dust production (larger humidity and vegetation cover: Reference McCulloch and DaviesMcCulloch and Davies, 2001; Reference RöthlisbergerRöthlisberger and others, 2002a) around 1 kyr before the ACR. The dating of the Ca2+ stabilization, around 15.0 kyr BP, is congruous with the Reference McCulloch and DaviesMcCulloch and Davies (2001) observation that southern Patagonia was humid from 16.91 kyr BP. As pointed out by Reference RöthlisbergerRöthlisberger and others (2002a), it remains unknown why the re-establishment of dry, cold conditions in Patagonia, from 15.33 kyr BP (Reference McCulloch and DaviesMcCulloch and Davies, 2001), was not recorded in the Ca2+ Dome C profile. Changes in atmospheric circulation, driven by the lower Equator-to-pole temperature gradient and the increased effect of zonal circulation, with respect to meridional transport (Reference Delmonte, Petit and MaggiDelmonte and others, 2002a), could explain this pattern. Analogous explanations could help in understanding the pattern of Na+ and other chemical components showing high sensitivity to the second warming step and the HO conditions (see flux profiles in Fig. 3).
Mainly snow accumulation-rate changes could affect, on the other hand, sulphate and ammonium snow content (with some possible exception during the HO), as observed from the concentration and fluxes along the EDC96 ice core.
Finally, the question of which process has driven the acid/base equilibrium of snow during the transition cannot be reliably answered from the chemical profile. Indeed, ECM shows the highest increase in the second step of the transition, after all chemical components reached low Holocene levels.
The observation of several previous deglaciation periods, recorded in the new EDC99 ice core, will allow confirmation of the observed environment–climate relationships and a better understanding of the factors controlling snow acidity.
This work is a contribution to the ‘European Project for Ice Coring in Antarctica’ (EPICA), a joint European Science Foundation (ESF)/European Commission (EC) scientific programme, funded by the European Commission and by national contributions from Belgium, Denmark, France, Germany, Italy, the Netherlands, Norway, Sweden, Switzerland and the United Kingdom. This is EPICA publication No. 75. The research was partially supported by Ente per le Nuove Tecnologie, l’Energia e L’Ambiente (ENEA) through a cooperation agreement with the Universities of Milan– Bicocca and Venice, in the framework of the ‘Glaciology’ and ‘Chemical Contamination’ sections of Programma Nazionale di Ricerche in Antartide (PNRA).