Specific conductivity and hydrogen ions measured on two ice core (NEEM and NGRIP)

A stratigraphy-based chronology for the North Greenland Eemian Ice Drilling (NEEM) ice core has been derived by transferring the annual layer counted Greenland Ice Core Chronology 2005 (GICC05) and its model extension (GICC05modelext) from the NGRIP core to the NEEM core using 787 match points of mainly volcanic origin identified in the electrical conductivity measurement (ECM) and dielectrical profiling (DEP) records. Tephra horizons found in both the NEEM and NGRIP ice cores are used to test the matching based on ECM and DEP and provide five additional horizons used for the timescale transfer. A thinning function reflecting the accumulated strain along the core has been determined using a Dansgaard-Johnsen flow model and an isotope-dependent accumulation rate parameterization. Flow parameters are determined from Monte Carlo analysis constrained by the observed depth-age horizons. In order to construct a chronology for the gas phase, the ice age-gas age difference (Delta age) has been reconstructed using a coupled firn densification-heat diffusion model. Temperature and accumulation inputs to the Delta age model, initially derived from the water isotope proxies, have been adjusted to optimize the fit to timing constraints from d15N of nitrogen and high-resolution methane data during the abrupt onset of Greenland interstadials. The ice and gas chronologies and the corresponding thinning function represent the first chronology for the NEEM core, named GICC05modelext-NEEM-1. Based on both the flow and firn modelling results, the accumulation history for the NEEM site has been reconstructed. Together, the timescale and accumulation reconstruction provide the necessary basis for further analysis of the records from NEEM.

Constraints on the magnitude and patterns of ocean cooling at the Last Glacial Maximum - Supplementary material

Observation-based reconstructions of sea surface temperature from relatively stable periods in the past, such as the Last Glacial Maximum, represent an important means of constraining climate sensitivity and evaluating model simulations. The first quantitative global reconstruction of sea surface temperatures during the Last Glacial Maximum was developed by the Climate Long-Range Investigation, Mapping and Prediction (CLIMAP) project in the 1970s and 1980s. Since that time, several shortcomings of that earlier effort have become apparent. Here we present an updated synthesis of sea surface temperatures during the Last Glacial Maximum, rigorously defined as the period between 23 and 19 thousand years before present, from the Multiproxy Approach for the Reconstruction of the Glacial Ocean Surface (MARGO) project. We integrate microfossil and geochemical reconstructions of surface temperatures and include assessments of the reliability of individual records. Our reconstruction reveals the presence of large longitudinal gradients in sea surface temperature in all of the ocean basins, in contrast to the simulations of the Last Glacial Maximum climate available at present.

Geochemistry of reduction haloes in ODP Holes 135-834A and 135-837A

The sediment column overlying basement in the Lau Basin consists of a sequence of volcaniclastic turbidites interbedded with hemipelagic clayey nannofossil mixed sediments, overlain in turn by a sequence of hemipelagic clayey nannofossil oozes containing sporadic calcareous turbidites. The clayey nannofossil oozes and mixed sediments are pervasively stained by hydrothermally derived iron and manganese oxyhydroxides. Sharply defined, lighter colored bands occur in the hemipelagic sediments, immediately beneath some (but by no means all) volcaniclastic and calcareous turbidites. These are identified as reduction haloes, of a type previously identified in quite different turbidite/pelagic sequences. The haloes are attributed to the burial of labile surficial Corg by turbidites, followed by the remineralization of this Corg with Mn and Fe oxyhydroxides as electron acceptors. The resultant characteristic Mn and Fe concentration/depth profiles are described, and a model is proposed for their development. The color alteration of the halo is ascribed to the removal of Mn oxyhydroxides, because, although the Fe content fluctuates through the haloes, this does not appear to affect their color. Other elements (Co, Cu, and Ni) are also at low concentration levels in the haloes like Mn, consistent with remobilization and migration out of the halo section, although the profile shapes are not identical with those of Mn. The behavior of V is distinctive in that it appears to have migrated into the haloes to be enriched there. Haloes are unlikely to form if turbidite emplacement is erosive and removes the near-surface layer, which generally is the most fluid part of the sediment and contains the highest levels of reactive Corg to drive the reduction process. Conversely, the presence of a halo implies that emplacement of the overlying turbidite did not significantly erode the pre-existing sediment/water interface.

Digital elevation model (DEM) of the ice sheet in western Dronning Maud Land, Antarctica

In this paper, a new digital elevation model (DEM) is derived for the ice sheet in western Dronning Maud Land, Antarctica. It is based on differential interferometric synthetic aperture radar (SAR) from the European Remote Sensing 1/2 (ERS-1/2) satellites, in combination with ICESat's Geoscience Laser Altimeter System (GLAS). A DEM mosaic is compiled out of 116 scenes from the ERS-1 ice phase in 1994 and the ERS-1/2 tandem mission between 1996 and 1997 with the GLAS data acquired in 2003 that served as ground control. Using three different SAR processors, uncertainties in phase stability and baseline model, resulting in height errors of up to 20 m, are exemplified. Atmospheric influences at the same order of magnitude are demonstrated, and corresponding scenes are excluded. For validation of the DEM mosaic, covering an area of about 130,000 km**2 on a 50-m grid, independent ICESat heights (2004-2007), ground-based kinematic GPS (2005), and airborne laser scanner data (ALS, 2007) are used. Excluding small areas with low phase coherence, the DEM differs in mean and standard deviation by 0.5 +/- 10.1, 1.1 +/- 6.4, and 3.1 +/- 4.0 m from ICESat, GPS, and ALS, respectively. The excluded data points may deviate by more than 50 m. In order to suppress the spatially variable noise below a 5-m threshold, 18% of the DEM area is selectively averaged to a final product at varying horizontal spatial resolution. Apart from mountainous areas, the new DEM outperforms other currently available DEMs and may serve as a benchmark for future elevation models such as from the TanDEM-X mission to spatially monitor ice sheet elevation.

Geochemistry and age determination of zircons in Cretaceous to Quaternary sediments

In central Antarctica, drainage today and earlier back to the Paleozoic radiates from the Gamburtsev Subglacial Mountains (GSM). Proximal to the GSM past the Permian-Triassic fluvial sandstones in the Prince Charles Mountains (PCM) are Cretaceous, Eocene, and Pleistocene sediment in Prydz Bay (ODP741, 1166, and 1167) and pre-Holocene sediment in AM04 beneath the Amery Ice Shelf. We analysed detrital zircons for U-Pb ages, Hf-isotope compositions, and trace elements to determine the age, rock type, source of the host magma, and "crustal" model age (T(C)DM). These samples, together with others downslope from the GSM and the Vostok Subglacial Highlands (VSH), define major clusters of detrital zircons interpreted as coming from (1) 700 to 460 Ma mafic granitoids and alkaline rock, epsilon-Hf 9 to -28, signifying derivation 2.5 to 1.3 Ga from fertile and recycled crust, and (2) 1200-900 Ma mafic granitoids and alkaline rock, epsilon-Hf 11 to -28, signifying derivation 1.8 to 1.3 Ga from fertile and recycled crust. Minor clusters extend to 3350 Ma. Similar detrital zircons in Permian-Triassic, Ordovician, Cambrian, and Neoproterozoic sandstones located along the PaleoPacific margin of East Antarctica and southeast Australia further downslope from central Antarctica reflect the upslope GSM-VSH nucleus of the central Antarctic provenance as a complex of 1200-900 Ma (Grenville) mafic granitoids and alkaline rocks and older rocks embedded in 700-460 Ma (Pan-Gondwanaland) fold belts. The wider central Antarctic provenance (CAP) is tentatively divided into a central sector with negative ?Hf in its 1200-900 Ma rocks bounded on either side by positive epsilon-Hf. The high ground of the GSM-VSH in the Permian and later to the present day is attributed to crustal shortening by far-field stress during the 320 Ma mid-Carboniferous collision of Gondwanaland and Laurussia. Earlier uplifts in the ~500 Ma Cambrian possibly followed the 700-500 Ma assembly of Gondwanaland, and in the Neoproterozoic the 1000-900 Ma collisional events in the Eastern Ghats-Rayner Province at the end of the 1300-1000 Ma assembly of Rodinia.

Lithogeochemical and isotope (Rb-Sr, Sm-Nd, Pb-Pb whole rock and Lu-Hf zircon) composition of samples from the Shackleton Range, East Antarctica

Three distinct, spatially separated crustal terranes have been recognised in the Shackleton Range, East Antarctica: the Southern, Eastern and Northern Terranes. Mafic gneisses from the Southern Terrane provide geochemical evidence for a within-plate, probably back-arc origin of their protoliths. A plume-distal ridge origin in an incipient ocean basin is the favoured interpretation for the emplacement site of these rocks at c. 1850 Ma, which, together with a few ocean island basalts, were subsequently incorporated into an accretionary continental arc/supra-subduction zone tectonic setting. Magmatic underplating resulted in partial melting of the lower crust, which caused high-temperature granulite-facies metamorphism in the Southern Terrane at c. 1710-1680 Ma. Mafic and felsic gneisses there are characterised by isotopically depleted, positive Nd and Hf initials and model ages between 2100 and 2000 Ma. They may be explained as juvenile additions to the crust towards the end of the Palaeoproterozoic. These juvenile rocks occur in a narrow, c. 150 km long E-W trending belt, inferred to trace a suture that is associated with a large Palaeoproterozoic accretionary orogenic system. The Southern Terrane contains many features that are similar to the Australo-Antarctic Mawson Continent and may be its furthermost extension into East Antarctica. The Eastern Terrane is characterised by metagranitoids that formed in a continental volcanic arc setting during a late Mesoproterozoic orogeny at c. 1060 Ma. Subsequently, the rocks experienced high-temperature metamorphism during Pan-African collisional tectonics at 600 Ma. Isotopically depleted zircon grains yielded Hf model ages of 1600-1400 Ma, which are identical to Nd model ages obtained from juvenile metagranitoids. Most likely, these rocks trace the suture related to the amalgamation of the Indo-Antarctic and West Gondwana continental blocks at ~600 Ma. The Eastern Terrane is interpreted as the southernmost extension of the Pan-African Mozambique/Maud Belt in East Antarctica and, based on Hf isotope data, may also represent a link to the Ellsworth-Whitmore Mountains block in West Antarctica and the Namaqua-Natal Province of southern Africa. Geochemical evidence indicates that the majority of the protoliths of the mafic gneisses in the Northern Terrane formed as oceanic island basalts in a within-plate setting. Subsequently the rocks were incorporated into a subduction zone environment and, finally, accreted to a continental margin during Pan-African collisional tectonics. Felsic gneisses there provide evidence for a within-plate and volcanic arc/collisional origin. Emplacement of granitoids occurred at c. 530 Ma and high-temperature, high-pressure metamorphism took place at 510-500 Ma. Enriched Hf and Nd initials and Palaeoproterozoic model ages for most samples indicate that no juvenile material was added to the crust of the Northern Terrane during the Pan-African Orogeny but recycling of older crust or mixing of crustal components of different age must have occurred. Isotopically depleted mafic gneisses, which are spatially associated with eclogite-facies pyroxenites, yielded late Mesoproterozoic Nd model ages. These rocks occur in a narrow, at least 100 km long, E-W trending belt that separates alkaline ocean island metabasalts and within-plate metagranitoids from volcanic arc metabasalts and volcanic arc/syn-collisional metagranitoids in the Northern Terrane. This belt is interpreted to trace the late Neoproterozoic/early Cambrian Pan-African collisional suture between the Australo-Antarctic and the combined Indo-Antarctic/West Gondwana continental blocks that formed during the final amalgamation of Gondwana.

Helium isotope ratios and sedimentation rate models of ODP holes

A continuous age model for the brief climate excursion at the Paleocene-Eocene boundary has been constructed by assuming a constant flux of extraterrestrial 3He (3He[ET]) to the seafloor. 3He[ET] measurements from ODP Site 690 provide quantitative evidence for the rapid onset (and short duration (<120 kyr) of global warming and of the associated disturbance to the Earth's surficial carbon budget at this time. These observations support astronomically calibrated age models indicating extremely rapid release of isotopically light carbon, possibly from seafloor methane hydrate, as the proximal cause of the event. However, the 3He[ET] technique indicates a previously unrecognized and extreme increase in sedimentation rate coincident with the return of climate proxies to pre-event values. The 3He[ET]-based age model thus suggests a far more rapid recovery from the climatic perturbation than previously proposed or predicted on the basis of the modern carbon cycle, and so may indicate additional or accelerated mechanisms of carbon removal from the ocean-atmosphere system during this period. 3He[ET] was also measured at ODP Site 1051 to test the validity of the Site 690 chronology. Comparison of these data sets seems to require removal of several tens of kyr of sediment within the climatic excursion at Site 1051, an observation consistent with sediment structures and previous age modeling efforts. The Site 1051 age model shows a ~30 kyr period in which climate proxies return toward pre-event values, after which they remain invariant for ~80 kyr. If this rise represents the recovery interval identified at Site 690, then the 3HeET-based age models of the two sites are in good agreement. However, alternative interpretations are possible, and work on less disrupted sites is required to evaluate the reliability of the proposed new chronology of the climate excursion. Regardless of these details, this work shows that the 3HeET technique can provide useful independent evidence for the development and testing of astronomically calibrated age models.

A probabilistic model estimating oil spill clean-up costs - a case study for the Gulf of Finland

Existing models estimating oil spill costs at sea are based on data from the past, and they usually lack a systematic approach. This make them passive, and limits their ability to forecast the effect of the changes in the oil combating fleet or location of a spill on the oil spill costs. In this paper we make an attempt towards the development of a probabilistic and systematic model estimating the costs of clean-up operations for the Gulf of Finland. For this purpose we utilize expert knowledge along with the available data and information from literature. Then, the obtained information is combined into a framework with the use of a Bayesian Belief Networks. Due to lack of data, we validate the model by comparing its results with existing models, with which we found good agreement. We anticipate that the presented model can contribute to the cost-effective oil-combating fleet optimization for the Gulf of Finland. It can also facilitate the accident consequences estimation in the framework of formal safety assessment (FSA).

Monthly Bonaire coral Sr/Ca and oxygen isotope data from 118 kyr ago (coral BON-5-D)

The end of the last interglacial period, ~118 kyr ago, was characterized by substantial ocean circulation and climate perturbations resulting from instabilities of polar ice sheets. These perturbations are crucial for a better understanding of future climate change. The seasonal temperature changes of the tropical ocean, however, which play an important role in seasonal climate extremes such as hurricanes, floods and droughts at the present day, are not well known for this period that led into the last glacial. Here we present a monthly resolved snapshot of reconstructed sea surface temperature in the tropical North Atlantic Ocean for 117.7±0.8 kyr ago, using coral Sr/Ca and d18O records. We find that temperature seasonality was similar to today, which is consistent with the orbital insolation forcing. Our coral and climate model results suggest that temperature seasonality of the tropical surface ocean is controlled mainly by orbital insolation changes during interglacials.

Sea-surface temperature reconstruction of sediments from the North Pacific and North Atlantic

Holocene climate variability is investigated in the North Pacific and North Atlantic realms, using alkenone-derived sea-surface temperature (SST) records as well as a millennial scale simulation with a coupled atmosphere-ocean general circulation model (AOGCM). The alkenone SST data indicate a temperature increase over almost the entire North Pacific from 7 cal kyr BP to the present. A dipole pattern with a continuous cooling in the northeastern Atlantic and a warming in the eastern Mediterranean Sea and the northern Red Sea is detected in the North Atlantic realm. Similarly, SST variations are opposite in sign between the northeastern Pacific and the northeastern Atlantic. A 2300 year long AOGCM climate simulation reveals a similar SST seesaw between the northeastern Pacific and the northeastern Atlantic on centennial time scales. Our analysis of the alkenone SST data and the model results suggests fundamental inter-oceanic teleconnections during the Holocene.

Stable isotope analysis and biomarkers of carbonates from the Columbia River in southwestern Washington, USA

Exotic limestone masses with silicified fossils, enclosed within deep-water marine siliciclastic sediments of the Early to Middle Miocene Astoria Formation, are exposed along the north shore of the Columbia River in southwestern Washington, USA. Samples from four localities were studied to clarify the origin and diagenesis of these limestone deposits. The bioturbated and reworked limestones contain a faunal assemblage resembling that of modern and Cenozoic deep-water methane-seeps. Five phases make up the paragenetic sequence: (1) micrite and microspar; (2) fibrous, banded and botryoidal aragonite cement, partially replaced by silica or recrystallized to calcite; (3) yellow calcite; (4) quartz replacing carbonate phases and quartz cement; and (5) equant calcite spar and pseudospar. Layers of pyrite frequently separate different carbonate phases and generations, indicating periods of corrosion. Negative d13Ccarbonate values as low as -37.6 per mill V-PDB reveal an uptake of methane-derived carbon. In other cases, d13Ccarbonate values as high as 7.1 per mill point to a residual, 13C-enriched carbon pool affected by methanogenesis. Lipid biomarkers include 13C-depleted, archaeal 2,6,10,15,19-pentamethylicosane (PMI; d13C: -128 per mill), crocetane and phytane, as well as various iso- and anteiso-carbon chains, most likely derived from sulphate-reducing bacteria. The biomarker inventory proves that the majority of the carbonates formed as a consequence of sulphate-dependent anaerobic oxidation of methane. Silicification of fossils and early diagenetic carbonate cements as well as the precipitation of quartz cement - also observed in other methane-seep limestones enclosed in sediments with abundant diatoms or radiolarians - is a consequence of a preceding increase of alkalinity due to anaerobic oxidation of methane, inducing the dissolution of silica skeletons. Once anaerobic oxidation of methane has ceased, the pH drops again and silica phases can precipitate.

10Be data for sediments from DSDP Leg 56 and ODP Leg 170

Sediment accretion and subduction at convergent margins play an important role in the nature of hazardous interplate seismicity (the seismogenic zone) and the subduction recycling of volatiles and continentally derived materials to the Earth's mantle. Identifying and quantifying sediment accretion, essential for a complete mass balance across the margin, can be difficult. Seismic images do not define the processes by which a prism was built, and cored sediments may show disturbed magnetostratigraphy and sparse biostratigraphy. This contribution reports the first use of cosmogenic 10Be depth profiles to define the origin and structural evolution of forearc sedimentary prisms. Biostratigraphy and 10Be model ages generally are in good agreement for sediments drilled at Deep Sea Drilling Project Site 434 in the Japan forearc, and support an origin by imbricate thrusting for the upper section. Forearc sediments from Ocean Drilling Program Site 1040 in Costa Rica lack good fossil or paleomagnetic age control above the decollement. Low and homogeneous 10Be concentrations show that the prism sediments are older than 3-4 Ma, and that the prism is either a paleoaccretionary prism or it formed largely from slump deposits of apron sediments. Low 10Be in Costa Rican lavas and the absence of frontal accretion imply deeper sediment underplating or subduction erosion.

Oxygen and silica flux from 9 multicorer profiles

Pore water profiles from 24 stations in the South Atlantic (located in the Guinea, Angola, Cape, Guyana, and Argentine basins) show good correlations of oxygen and silicon, suggesting microbially mediated dissolution of biogenic silica. We used simple analytical transport and reaction models to show the tight coupling of the reconstructed process kinetics of aerobic respiration and silicon regeneration. A generic transport and reaction model successfully reproduced the majority of Si pore water profiles from aerobic respiration rates, confirming that the dissolution of biogenic silica (BSi) occurs proportionally to O2 consumption. Possibly limited to well-oxygenated sediments poor in BSi, benthic Si fluxes can be inferred from O2 uptake with satisfactory accuracy. Compared to aerobic respiration kinetics, the solubility of BSi emerged as a less influential parameter for silicon regeneration. Understanding the role of bacteria for silicon regeneration requires further investigations, some of which are outlined. The proposed aerobic respiration control of benthic silicon cycling is suitable for benthic-pelagic models. The empirical relation of BSi dissolution to aerobic respiration can be used for regionalization assessments and estimates of the silicon budget to increase the understanding of global primary and export production patterns.