Snow depths of volcanic events and accumulation rates of firn cores from the IPY traverse, Antarctica

Volcanic signatures in ice-core records provide an excellent means to date the cores and obtain information about accumulation rates. From several ice cores it is thus possible to extract a spatio-temporal accumulation pattern. We show records of electrical conductivity and sulfur from 13 firn cores from the Norwegian-USA scientific traverse during the International Polar Year 2007-2009 (IPY) through East Antarctica. Major volcanic eruptions are identified and used to assess century-scale accumulation changes. The largest changes seem to occur in the most recent decades with accumulation over the period 1963-2007/08 being up to 25% different from the long-term record. There is no clear overall trend, some sites show an increase in accumulation over the period 1963 to present while others show a decrease. Almost all of the sites above 3200 m above sea level (asl) suggest a decrease. These sites also show a significantly lower accumulation value than large-scale assessments both for the period 1963 to present and for the long-term mean at the respective drill sites. The spatial accumulation distribution is influenced mainly by elevation and distance to the ocean (continentality), as expected. Ground-penetrating radar data around the drill sites show a spatial variability within 10-20% over several tens of kilometers, indicating that our drill sites are well representative for the area around them. Our results are important for large-scale assessments of Antarctic mass balance and model validation.

Chronology for the EDML ice core from Dronning Maud Land, Antarctica, over the last 150 000 years

A chronology called EDML1 has been developed for the EPICA ice core from Dronning Maud Land (EDML). EDML1 is closely interlinked with EDC3, the new chronology for the EPICA ice core from Dome-C (EDC) through a stratigraphic match between EDML and EDC that consists of 322 volcanic match points over the last 128 ka. The EDC3 chronology comprises a glaciological model at EDC, which is constrained and later selectively tuned using primary dating information from EDC as well as from EDML, the latter being transferred using the tight stratigraphic link between the two cores. Finally, EDML1 was built by exporting EDC3 to EDML. For ages younger than 41 ka BP the new synchronized time scale EDML1/EDC3 is based on dated volcanic events and on a match to the Greenlandic ice core chronology GICC05 via 10Be and methane. The internal consistency between EDML1 and EDC3 is estimated to be typically ~6 years and always less than 450 years over the last 128 ka (always less than 130 years over the last 60 ka), which reflects an unprecedented synchrony of time scales. EDML1 ends at 150 ka BP (2417 m depth) because the match between EDML and EDC becomes ambiguous further down. This hints at a complex ice flow history for the deepest 350 m of the EDML ice core.

Deuterium analysis of ice core GISP2

The Greenland Ice Sheet Project 2 (GISP2) core can enhance our understanding of the relationship between parameters measured in the ice in central Greenland and variability in the ocean, atmosphere, and cryosphere of the North Atlantic Ocean and adjacent land masses. Seasonal (summer, winter) to annual responses of dD and deuterium excess isotopic signals in the GISP2 core to the seesaw in winter temperatures between West Greenland and northern Europe from A.D. 1840 to 1970 are investigated. This seesaw represents extreme modes of the North Atlantic Oscillation, which also influences sea surface temperatures (SSTs), atmospheric pressures, geostrophic wind strength, and sea ice extents beyond the winter season. Temperature excursions inferred from the dD record during seesaw/extreme NAO mode years move in the same direction as the West Greenland side of the seesaw. Symmetry with the West Greenland side of the seesaw suggests a possible mechanism for damping in the ice core record of the lowest decadal temperatures experienced in Europe from A.D. 1500 to 1700. Seasonal and annual deuterium excess excursions during seesaw years show negative correlation with dD. This suggests an isotopic response to a SST/ land temperature seesaw. The isotopic record from GISP2 may therefore give information on both ice sheet and sea surface temperature variability. Cross-plots of dD and d show a tendency for data to be grouped according to the prevailing mode of the seesaw, but do not provide unambiguous identification of individual seesaw years. A combination of ice core and tree ring data sets may allow more confident identification of GA and GB (extreme NAO mode) years prior to 1840.

Match points used for synchronization of NGRIP, GRIP, and GISP2

We here present a synchronization of the NGRIP, GRIP, and GISP2 ice cores based mainly on volcanic events over the period 14.9-32.45 ka b2k (before AD 2000), corresponding to Marine Isotope Stage 2 (MIS 2) and the end of MIS 3. The matching provides a basis for applying the recent NGRIP-based Greenland Ice Core Chronology 2005 (GICC05) time scale to the GRIP and GISP2 ice cores, thereby making it possible to compare the synchronized palaeoclimate profiles of the cores in detail and to identify relative accumulation differences between the cores. Based on the matching, a period of anomalous high accumulation rates in the GISP2 ice core is detected within the period 16.5-18.3 ka b2k. The d18O and [Ca2+] profiles of the three cores are presented on the common GICC05 time scale and generally show excellent agreement across the stadial-interstadial transitions and across the two characteristic dust events in Greenland Stadial 3. However, large differences between the d18O and [Ca2+] profiles of the three cores are seen in the same period as the 7-9% increase in the GISP2 accumulation rate. We conclude that changes of the atmospheric circulation are likely to have occurred in this period, altering the spatial gradients in Greenland and resulting in larger variations between the records.

Synchronization of the NGRIP, GRIP, and GISP2 ice cores

We here present a synchronization of the NGRIP, GRIP, and GISP2 ice cores based mainly on volcanic events over the period 14.9-32.45 ka b2k (before AD 2000), corresponding to Marine Isotope Stage 2 (MIS 2) and the end of MIS 3. The matching provides a basis for applying the recent NGRIP-based Greenland Ice Core Chronology 2005 (GICC05) time scale to the GRIP and GISP2 ice cores, thereby making it possible to compare the synchronized palaeoclimate profiles of the cores in detail and to identify relative accumulation differences between the cores. Based on the matching, a period of anomalous high accumulation rates in the GISP2 ice core is detected within the period 16.5-18.3 ka b2k. The d18O and [Ca2+] profiles of the three cores are presented on the common GICC05 time scale and generally show excellent agreement across the stadial-interstadial transitions and across the two characteristic dust events in Greenland Stadial 3. However, large differences between the d18O and [Ca2+] profiles of the three cores are seen in the same period as the 7-9% increase in the GISP2 accumulation rate. We conclude that changes of the atmospheric circulation are likely to have occurred in this period, altering the spatial gradients in Greenland and resulting in larger variations between the records.

Miocene Mediterranean Nd, Sr, C and O isotope composition

Isotopic compositions of marine sediments and fossils have been investigated from northern basins of the Mediterranean to help constrain local oceanographic and climatic changes adjacent to the uplifting Alps. Stable C and O isotope compositions of benthic and planktonic foraminifera from the Umbria-Marche region (UMC) have an offset characteristic for their habitats and the changes in composition mimic global changes, suggesting that the regional conditions of climate and the carbon cycle were controlled by global changes. The radiogenic isotope composition of these fossil assemblages allows recognition of three distinct periods. In the first period, from 25 to 19 Ma, high epsilon-Nd values and low 87Sr/86Sr of sediments and fossils support intense tectonism and volcanism, related to the opening of the western Mediterranean. In the second period, from 19 to 13 Ma the 87Sr/86Sr ratio of Mediterranean (UMC) deviate from the global ocean, which is compatible with rapid uplift of the hinterland and intense influx of Sr from Mesozoic carbonates of the western Apennines. This local control on the seawater was driven by a humid and warm climate and indicates restricted exchange of water with the global ocean. Generally, the epsilon-Nd values of the fossils are very similar to those of Indian Ocean water, with brief periods of a decrease in the epsilon-Nd values coinciding with volcanic events and maybe sea level variation at 15.2 Ma. In the third period, from 13 to 10 Ma the fossils have 87Sr/86Sr similar to those of Miocene seawater while their epsilon-Nd values change considerably with time. This indicates fluctuating influence of the Atlantic versus the Paratethys and/or locally evolved seawater in the Mediterranean driven by global sea level changes. Other investigated localities near the Alps and from the ODP 900 site are compatible with this oceanographic interpretation. However, in the late early Miocene, enhanced local control, reflecting erosion of old crustal silicate rocks near the Alps, results in higher 87Sr/86Sr.

Investigations on seven sediment cores from the New Ireland Basin, east of Papua New Guinea

Sediment cores were recovered from the New Ireland Basin, east of Papua New Guinea, in order to investigate the late Quaternary eruptive history of the Tabar-Lihir-Tanga-Feni (TLTF) volcanic chain. Foraminifera d18O profiles were matched to the low-latitude oxygen isotope record to date the cores, which extend back to the early part of d18O Stage 9 (333 ka). Sedimentation rates decrease from >10 cm/1000 yr in cores near New Ireland to ~2 cm/1000 yr further offshore. The cores contain 36 discrete ash beds, mostly 1-8 cm thick and interpreted as either fallout or distal turbidite deposits. Most beds have compositionally homogeneous glass shard populations, indicating that they represent single volcanic events. Shards from all ash beds have the subduction-related pattern of strong enrichment in the large-ion lithophile elements relative to MORB, but three distinct compositional groups are apparent: Group A beds are shoshonitic and characterised by >1300 ppm Sr, high Ce/Yb and high Nb/Yb relative to MORB, Group B beds form a high-K series with MORB-like Nb/Yb but high Ce/Yb and well-developed negative Eu anomalies, whereas Group C beds are transitional between the low-K and medium-K series and characterised by flat chondrite-normalised REE patterns with low Nb/Yb relative to MORB. A comparison with published data from the TLTF chain, the New Britain volcanic arc and backarc including Rabaul, and Bagana on Bougainville demonstrates that only Group A beds share the distinctive phenocryst assemblage and shoshonitic geochemistry of the TLTF lavas. The crystal- and lithic-rich character of the Group A beds point to a nearby source, and their high Sr, Ce/Yb and Nb/Yb match those of Tanga and Feni lavas. A youthful stratocone on the eastern side of Babase Island in the Feni group is the most probable source. Group A beds younger than 20 ka are more fractionated than the older Group A beds, and record the progressive development of a shallow level magma chamber beneath the cone. In contrast, Group B beds represent glass-rich fallout from voluminous eruptions at Rabaul, whereas Group C beds represent distal glass-rich fallout from elsewhere along the volcanic front of the New Britain arc.

Talos Dome Ice Core Deuterium Data

During the 1996 Programma Nazionale di Ricerche in Antartide-International Trans-Antarctic Scientific Expedition traverse, two firn cores were retrieved from the Talos Dome area (East Antarctica) at elevations of 2316 m (TD, 89 m long) and 2246 m (ST556, 19 m long). Cores were dated by using seasonal variations in non-sea-salt (nss) SO42- concentrations coupled with the recognition of tritium marker level (1965-1966) and nss SO42- spikes due to the most important volcanic events in the past (Pinatubo 1991, Agung 1963, Krakatoa 1883, Tambora 1815, Kuwae 1452, Unknown 1259). The number of annual layers recognized in the TD and ST556 cores was 779 and 97, respectively. The dD record obtained from the TD core has been compared with other East Antarctic isotope ice core records (Dome C EPICA, South Pole, Taylor Dome). These records suggest cooler climate conditions between the middle of 16th and the beginning of 19th centuries, which might be related to the Little Ice Age (LIA) cold period. Because of the high degree of geographical variability, the strongest LIA cooling was not temporally synchronous over East Antarctica, and the analyzed records do not provide a coherent picture for East Antarctica. The accumulation rate record presented for the TD core shows a decrease during part of the LIA followed by an increment of about 11% in accumulation during the 20th century. At the ST556 site, the accumulation rate observed during the 20th century was quite stable.

Major and trace element analyses of ODP Leg 120 volcanic rocks (Table 1)

Widespread Lower Cretaceous magmatism occurred along the Indian-Australian/Antarctic margins, and in the juvenile Indian Ocean, during the rifting of eastern Gondwana. The formation of this magmatic province probably began around 120-130 Ma with the eruption of basalts on the Naturaliste Plateau and at Bunbury, western Australia. On the northeast margin of India, activity began around 117 Ma with the Rajmahal continental basalts and associated lamprophyre intrusions. The formation of the Kerguelen Plateau in the Indian Ocean began no later than 114 Ma. Ultramafic lamprophyres (alnoites) were emplaced in the Prince Charles Mountains near the Antarctic continental margin at ~ 110 Ma. These events are considered to be related to a major mantle plume, the remnant of which is situated beneath the region of Kerguelen and Heard islands at the present day. Geochemical data are presented for each of these volcanic suites and are indicative of complex interactions between asthenosphere-derived magmas and the continental lithosphere. Kerguelen Plateau basalts have Sr and Nd isotopic compositions lying outside the field for Indian Ocean mid-ocean ridge basalts (MORB) but, with the exception of Site 738 at the southern end of the plateau, within the range of more recent hotspot basalts from Kerguelen and Heard Islands. However, a number of the plateau tholeiites are characterized by lower 206Pb/204Pb ratios than are basalts from Kerguelen Island, and many also have anomalously high La/Nb ratios. These features suggest that the source of the Kerguelen Plateau basalts suffered contamination by components derived from the Gondwana continental lithosphere. An extreme expression of this lithospheric signature is shown by a tholeiite from Site 738, suggesting that the southernmost part of the Kerguelen Plateau may be underlain by continental crust. The Rajmahal tholeiites mostly fall into two distinct geochemical groups. Some Group I tholeiites have Sr and Nd isotopic compositions and incompatible element abundances, similar to Kerguelen Plateau tholeiites from Sites 749 and 750, indicating that the Kerguelen-Heard mantle plume may have directly furnished Rajmahal volcanism. However, their elevated 207Pb/204Pb ratios indicate that these magmas did not totally escape contamination by continental lithosphere. In contrast to the Group I tholeiites, significant contamination is suggested for Group II Rajmahal tholeiites, on the basis of incompatible element abundances and isotopic compositions. The Naturaliste Plateau and the Bunbury Basalt samples show varying degrees of enrichment in incompatible elements over normal MORB. The Naturaliste Plateau samples (and Bunbury Basalt) have high La/Nb ratios, a feature not inconsistent with the notion that the plateau may consist of stretched continental lithosphere, near the ocean-continent divide.

Chemistry of smectite and volcanic glass of ODP Site 159-959

The Ocean Drilling Program (ODP) Site 959 was drilled in the northern border of the Côte d'Ivoire-Ghana Ridge at a water depth of 2100 m. Pleistocene total thickness does not exceed 20 m. Winnowing processes resulted in a low accumulation rate and notable stratigraphic hiatuses. During the Late Pleistocene, bottom circulation was very active and controlled laminae deposition (contourites) which increased the concentration of glauconitic infillings of foraminifera, and of volcanic glass and blue-green grains more rarely, with one or several subordinate ferromagnesian silicates. Volcanic glass generally was X-ray amorphous and schematically classified as basic to intermediate (44-60% SiO2). Opal-A or opal-CT suggested the beginning of the palagonitisation process, and previous smectitic deposits may have been eroded mechanically. The blue-green grains presented two main types of mineralogic composition: (1) neoformed K, Fe-smectite associated with zeolite (like phillipsite) and unequal amounts of quartz and anorthite; (2) feldspathic grains dominated by albite but including quartz, volcanic glass and smectites as accessory components. They were more or less associated with the volcanic glass. On the basis of their chemical composition, the genetic relationship between the blue-green grains and the volcanic glass seemed to be obvious although some heterogeneous grains seemed to be primary ignimbrite and not the result of glass weathering. The most reasonable origin of these pyroclastic ejecta would be explosive events from the Cameroon Volcanic Ridge, especially from the Sao Thome and Principe Islands and Mount Cameroon area. This is supported both by grain geochemistry and the time of volcanic activity, i.e. Pleistocene. After westward wind transport (some 1200 km) and ash fall-out, the subsequent winnowing by bottom currents controlled the concentration of the volcanic grains previously disseminated inside the hemipelagic sediment. Palagonitisation, and especially phillipsite formation, may result from a relatively rapid reaction during burial diagenesis (<1 m.y.), in deep-sea deposits at relatively low sedimentation rate. However, it cannot be excluded that the weathering had begun widely on the Cameroon Ridge before the explosive event.

Geochemistry of Pliocene-Pleistocene volcanic sands of ODP Site 136-842

Several distinct, thin (2-7 cm), volcanic sand layers ("ashes") were recovered in the upper portions of Holes 842A and 842B. These holes were drilled 320 km west of the island of Hawaii on the outer side of the arch that surrounds the southern end of the Hawaiian chain. These layers are Pliocene to Pleistocene in age, graded, and contain fresh glass and mineral fragments (mainly olivine, plagioclase, and clinopyroxene) and tests of Pleistocene to Eocene radiolarians. The glass fragments are weakly vesicular and blocky to platy in shape. The glass and olivine fragments from individual layers have large ranges in composition (i.e, larger than expected for a single eruption). These features are inconsistent with an explosive eruption origin for the sands. The only other viable mechanism for transporting these sands hundreds of kilometers from their probable source, the Hawaiian Islands, is turbidity currents. These currents were probably related to several of the giant debris slides that were identified from Gloria sidescan images around the islands. These currents would have run over the ~500-m-high Hawaiian Arch on their way to Site 842. This indicates that the turbidity currents were at least 325 m thick. Paleomagnetic and biostratigraphic data allow the ages of the sands to be constrained and, thus, related to particular Hawaiian debris flows. These correlations were checked by comparing the compositions of the glasses from the sands with those of glasses and rocks from islands with debris flows directed toward Site 842. Good correlations were found for the 110-ka slide from Mauna Loa and the ~1.4-Ma slide from Lanai. The correlation with Kauai is poor, probably because the data base for that volcano is small. The low to moderate sulfur content of the sand glasses indicates that they were derived from moderately to strongly degassed lavas (shallow marine or subaerially erupted), which correlates well with the location of the landslide scars on the flanks of the Hawaiian volcanoes. The glass sands may have been formed by brecciation during the landslide events or spallation and granulation as lava erupted into shallow water.

Biostratigraphy of selected taxa and radiolarian events in ODP Hole 183-1138A sediments

Well-preserved radiolarian assemblages of late middle Miocene to early Pliocene age are found in Ocean Drilling Program (ODP) Hole 1138A (Cores 183-1138A-12R to 20R), which was rotary drilled into the Central Kerguelen Plateau. The faunas are typical for Antarctic assemblages of this time interval, and the site appears to have been south of the Polar Front during the time period studied. Despite only moderate drilling recovery of the section, most late middle to early Pliocene radiolarian zones are present, although at the sample resolution used, subzones could not be identified. A significant discontinuity in the section is present at the boundary between lithologic Units I and II (between Cores 183-1138A-12R and 13R), corresponding to an interval from at least 4.6 to 6.1 Ma. Mixed late Miocene-early Pliocene assemblages are seen in the base of Core 183-1138A-12R (Sample 183-1138A-12R-3, 20 cm), and the overlying basal Pliocene Tau Zone appears to be absent. It cannot be determined if the discontinuity is due to incomplete recovery of the section and drilling disturbance or if it reflects a primary sedimentary structure - a hiatus or interval of condensed sedimentation.

Geochemistry of tephra from the Taupo Volcanic Zone

The Taupo Volcanic Zone (TVZ), central North Island, New Zealand, is the most frequently active Quaternary rhyolitic system in the world. Silicic tephras recovered from Ocean Drilling Programme Site 1123 (41°47.16'S, 171°29.94'W; 3290 m water depth) in the southwest Pacific Ocean provide a well-dated record of explosive TVZ volcanism since ~1.65 Ma. We present major, minor and trace element data for 70 Quaternary tephra layers from Site 1123 determined by electron probe microanalysis (1314 analyses) and laser ablation inductively coupled plasma mass spectrometry (654 analyses). Trace element data allow for the discrimination of different tephras with similar major element chemistries and the establishment of isochronous tie-lines between three sediment cores (1123A, 1123B and 1123C) recovered from Site 1123. These tephra tie-lines are used to evaluate the stratigraphy and orbitally tuned stable isotope age model of the Site 1123 composite record. Trace element fingerprinting of tephras identifies ~4.5 m and ~7.9 m thick sections of repeated sediments in 1123A (49.0-53.5 mbsf [metres below seafloor]) and 1123C (48.1-56.0 mbsf), respectively. These previously unrecognised repeated sections have resulted in significant errors in the Site 1123 composite stratigraphy and age model for the interval 1.15-1.38 Ma and can explain the poor correspondence between d18O profiles for Site 1123 and Site 849 (equatorial Pacific) during this interval. The revised composite stratigraphy for Site 1123 shows that the 70 tephra layers, when correlated between cores, correspond to ~37-38 individual eruptive events (tephras), 7 of which can be correlated to onshore TVZ deposits. The frequency of large-volume TVZ-derived silicic eruptions, as recorded by the deposition of tephras at Site 1123, has not been uniform through time. Rather it has been typified by short periods (25-50 ka) of intense activity bracketed by longer periods (100-130 ka) of quiescence. The most active period (at least 1 event per 7 ka) occurred between ~1.53 and 1.66 Ma, corresponding to the first ~130 ka of TVZ rhyolitic magmatism. Since 1.2 Ma, ~80% of tephras preserved at Site 1123 and the more proximal Site 1124 were erupted and deposited during glacial periods. This feature may reflect either enhanced atmospheric transport of volcanic ash to these sites (up to 1000 km from source) during glacial conditions or, more speculatively, that these events are triggered by changes in crustal stress accumulation associated with large amplitude sea-level changes. Only 8 of the ~37-38 Site 1123 tephra units (~20%) can be found in all three cores, and 22 tephra units (~60%) are only present in one of the three cores. Whether a tephra is preserved in all three cores does not have any direct relationship to eruptive volume. Instead it is postulated that tephra preservation at Site 1123 is 'patchy' and influenced by the vigorous nature of their deposition to the deep ocean floor as vertical density currents. At this site, at least 5 cores would need to have been drilled within a proximity of 10's to 100's of metres of each other to yield a >99% chance of recovering all the silicic tephras deposited on the ocean surface above it in the past 1.65 Ma.