ALBERT

Description: The Argentine margin contains important sedimentological, paleontological and chemical records of regional and local tectonic evolution, sea level, climate evolution and ocean circulation since the opening of the South Atlantic in the Late Jurassic–Early Cretaceous as well as the present-day results of post-depositional chemical and biological alteration. Despite its important location, which underlies the exchange of southern- and northern-sourced water masses, the Argentine margin has not been investigated in detail using scientific drilling techniques, perhaps because the margin has the reputation of being erosional. However, a number of papers published since 2009 have reported new high-resolution and/or multichannel seismic surveys, often combined with multi-beam bathymetric data, which show the common occurrence of layered sediments and prominent sediment drifts on the Argentine and adjacent Uruguayan margins. There has also been significant progress in studying the climatic records in surficial and near-surface sediments recovered in sediment cores from the Argentine margin. Encouraged by these recent results, our 3.5-day IODP (International Ocean Discovery Program) workshop in Buenos Aires (8–11 September 2015) focused on opportunities for scientific drilling on the Atlantic margin of Argentina, which lies beneath a key portion of the global ocean conveyor belt of thermohaline circulation. Significant opportunities exist to study the tectonic evolution, paleoceanography and stratigraphy, sedimentology, and biosphere and geochemistry of this margin.

Description: The current generation of marine biogeochemical modules in Earth system models (ESMs) considers mainly the effect of marine biota on the carbon cycle. We propose to also implement other biologically driven mechanisms in ESMs so that more climate-relevant feedbacks are captured. We classify these mechanisms in three categories according to their functional role in the Earth system: (1) "biogeochemical pumps", which affect the carbon cycling; (2) "biological gas and particle shuttles", which affect the atmospheric composition; and (3) "biogeophysical mechanisms", which affect the thermal, optical, and mechanical properties of the ocean. To resolve mechanisms from all three classes, we find it sufficient to include five functional groups: bulk phyto- and zooplankton, calcifiers, and coastal gas and surface mat producers. We strongly suggest to account for a larger mechanism diversity in ESMs in the future to improve the quality of climate projections.

Description: EGU2010-13373 The frequency of volcanic activity varies on a wide rangeof spatial and temporal scales, from 〈1 yr. periodicities in single volcanic systems to periodicities of 106 yrs. in global volcanism. The causes of these periodicities are poorly understood although the long-term global variations are likely linked to plate-tectonic processes. Here we present evidence for temporal changes in eruption frequencies at an intermediate time scale (104 yrs.) using the Pleistocene to recent records of widespread tephras of sub-Plinian to Plinian, and occasionally co-ignimbrite origin, along the Pacific Ring of Fire, which accounts for about half of the global length of 44,000 km of active subduction. Eruptions at arc volcanoes tend to be highly explosive and the well-preserved tephra records from the ocean floor can be assumed to be representative of how eruption frequencies varied with time. Volcanic activity along the Pacific Ring of Fire evolved through alternating phases of high and low frequency; although there is modulation by local and regional geologic conditions, these variations have a statistically significant periodicity of 43 ka that overlaps with the temporal variation in the obliquity of the Earth’s rotation axis, an orbital parameter that also exerts a strong control on global climate changes. This may suggest that the frequency of volcanic activity is controlled by effects of global climate changes. However, the strongest physical effects of climate change occur at 100 ka periods which are not seen in the volcanic record. We therefore propose that the frequency of volcanic activity is directly influenced by minute changes in the tidal forces induced by the varying obliquity resulting in long-period gravitational disturbances acting on the upper mantle.

Description: EGU2010-10518 Seafloor compliance is the transfer function between pressure and vertical displacement at the seafloor Infragravity waves in the oceanic layer have long periods in the range of 30 – 500 s and obey a simple frequencywavenumber relation. Seafloor compliance from infragravity waves can be analyzed with single station recordings to determinate sub-seafloor shear wave velocities. Previous studies in the Pacific Ocean have demonstrated that reliable near-surface shear wave profiles can be derived from infragravity wave compliance. However, these studies indicate that, beside the water depth the compliance measurements are limited by instrument sensitivity, calibration uncertainties and possibly other effects. In this work seafloor compliance and infragravity waves are observed at two different locations in the Atlantic Ocean: the Logatchev hydrothermal field at the Mid Atlantic Ridge and the Azores (Sao Miguel Island). The data was acquired with the broadband ocean compliance station developed at the University of Hamburg as well as ocean station from the German instrument pool for amphibian seismology (DEPAS) equipped with broadband seismometers and pressure sensors. Vertical velocity and pressure data were used to calculate power spectral densities and normalized compliance along two profiles (one in each location). Power spectral densities show a dominant peak at low frequencies (0.01-0.035Hz) limited by the expected cut-off frequency, which is dependent on the water depth at each station. The peak has been interpreted as a strong infragravity wave with values between 10-14 and 10-11 (m/s2)2/Hz and 104 and 106 (Pa2)2/Hz for acceleration and pressure respectively. The results show compliance values between 10-10 and 10-8 1/Pa and its estimations take into account the coherence between seismic and pressure signals in order to confirm that the seismic signals in the infragravity waves are caused by pressure sources. Shear wave velocity models, with depth resolution from 200 to 2500 m for the deep water stations, were derived from compliance. Preliminary results indicate shear wave velocity increasing from 200 to 3500 m/s.

Description: EGU2010-9841 Active mud volcanoes, where changing salinities of pore fluids, large temperature gradients and occurrences of free gas are frequently observed, should potentially exhibit significant variability in their internal resistivity structure. This is due to the fact that the bulk resistivity is mainly determined by the porosity of sediments and the electrical resistivity of the pore filling contained therein. The resistivity variations may be derived from controlled source electromagnetic (CSEM) measurements. CSEM systems consist of an electric dipole transmitter producing a time varying source field and electric dipole receivers, which measure the earth´s response to this signal. For a RWE Dea funded investigation of fluid and gas leakages at the North Alex Mud Volcano (NAMV) - a comparatively small target with an area of about 1km2 - we have developed a new high resolution CSEM system. The system consists of several autonomous electric dipole receivers and a lightweight electric dipole transmitter, which can be mounted on a small remotely operated underwater vehicle (ROV). The use of a ROV allows for a precise placement of the transmitter, which is a necessary prerequisite for the investigation of such a small target. Furthermore, electromagnetic signals may be transmitted from different directions with respect to the stationary receivers, allowing for a 3D-style tomographic experiment. In this experiment, ten receivers were deployed over the surface of NAMV at a total of 16 receiver locations. During three successful dives with a Cherokee ROV (Ghent University, Belgium), the transmitter was deployed at a total of 80 locations. Here we present first quantitative results consisting of apparent resistivity estimations from the CSEM time domain data for each transmitter-receiver pair. The apparent resistivity map shows that the NAMV indeed has a heterogeneous resistivity structure with apparent resistivities varying by at least a factor of two: low apparent resistivities (~ 0.8Ωm) are found towards the center of the MV, whereas higher apparent resistivities (~ 1.6Ωm) prevail away from the center. In a second step, we interpret the time-domain data based on 1D inversions. Good data fits can be achieved by models containing 2-3 layers. Generally, the models indicate low resistivities at the surface, which can be associated with penetrating salt water and/or high temperatures. Toward greater depths, increasing resistivities presumably are due to a combination of compaction of sediments (i.e. reduced pore space), an increased presence of fresh water and possible occurrences of free gas. For some 1D models, the increase in resistivity exceeds a factor of 10 or more and layer interfaces are indicated down to depths of up to 70m. The derived resistivity variations observed at the NAMV will be interpreted in conjunction with temperature (Feseker, this session), fluid flow (Brückmann et al., this session) and seismic data (Bialas et al., this session) acquired. Temperature variations measured in the upper few meters are related to fluid flow, where high temperatures are indicative of upwelling fluids of low salinity and low temperature of either a downward flow of saline fluids or no flow activity. This type of surface measurement constitutes an integrative fluid flow gauge, which we can resolve vertically with our resistivity models. Seismic data yield a background structure to our resistivity model. New analysis of seismic data shows that seismic activity may also be linked to fluid flow activity, which we aim to match with resistivity variations and oscillations, which were observed in the electric and magnetic fields (Lefeldt et al., this session).

Description: Recently Hathorne et al. (2009) documented large intratest trace element (TE) variations in planktonic foraminifera from a single sediment trap sample that could not be explained by variations in water column properties. The laser ablation ICP-MS depth profiles of trace elements through the test walls revealed strong positive correlations between Li, Mg, Mn and Ba resulting from the mixing of a lower TE outer calcite with a higher TE inner calcite. In contrast Sr/Ca ratios remained relatively constant throughout the test wall. These intratest TE variations likely result from biomineralization processes and therefore should be explained by any valid biomineralization model. However, changes in calcite precipitation rate, crystal structure, or the chemical composition of the internal calcification reservoir could not, by themselves, fully account for the pattern of cation intratest variability. Here I expand on this work and investigate if a model of coral biomineralization by Sinclair and Risk (2006) can be adapted to explain the pattern of intratest TE variability in foraminifera. It is clear that the low Mg calcite secreting foraminifera must reduce the Mg/Ca ratio of the calcifying solution by at least a factor of 10 (e.g. Hathorne et al., 2009) and it has been suggested this is achieved by active removal of Mg from the calcification reservoir, although the actual mechanisms remain debatable (e.g. Bentov and Erez, 2006). However, a recent study of the calcification of a low Mg calcite species in the laboratory found a large shortcoming in the amount of Ca potentially provided by seawater transported to the site of calcification in vacuoles compared to a conservative estimate of the amount required to form the new calcite wall (de Nooijer et al., 2009a). This suggests active Ca transport to the site of calcification is required to provide sufficient Ca. If Ca specific, this Ca addition would effectively dilute the TE content (including Mg) of the calcification reservoir to varying degrees and potentially cause the positive TE correlations seen across the test wall. Sinclair and Risk (2006) used this dilution model to successfully explain some TE correlations in coral skeletons. This model can be effectively adapted to foraminifera as it accounts for recent observations of foraminiferal calcification including the transport of seawater by liquid endocytosis to the calcification site and an elevated pH at the site of calcification (Bentov et al., 2009; de Nooijer et al., 2009a, 2009b). This model therefore provides a powerful tool with which to integrate constraints from experimental observation with those from micro-analytical measurements to improve the accuracy, precision and scale of the palaepalaeoceanographic application of foraminiferal geochemistry. Bentov and Erez (2006) Geochem. Geophys. Gepsyst. 7, Q01P08. Bentov et al. (2009) PNAS 106, 21500. de Nooijer et al. (2009a) Biogeosciences 6, 2669. de Nooijer et al. (2009b) PNAS 106, 15374. Hathorne et al. (2009) Paleoceanography 24, PA4204. Sinclair and Risk (2006) Geochim. Cosmochim. Acta 70, 3855.

Description: We present Mg/Ca analyses performed via a Flow Through sequential dissolution device connected to an ICP-OES on the planktonic foraminifer Globorotalia inflata. The aim of the study is to explore the possibility to reconstruct the thermal gradient in the water column by separating non-crusted and crusted calcite phases in the tests of G. inflata using the difference between their Mg/Ca ratios as a measure of the thermal gradient. An important assumption is that the non-crusted part of the tests is calcified in shallow, warmer water than the crusted part. For analyses a range of different preparation steps were used to determine the ideal way of separating the phases. Foraminifer tests were (not) cleaned, (not) crushed, and (not) pulverized before online analysis with the FT device. To analyze samples with a FT device the foraminifer tests are placed on a filter with a mesh of 0.45 μm preventing clay minerals to wash through. A sequential dissolution protocol first rinses the samples with buffered Seralpur water before QD HNO3 is added in small steps to create a ramp of increasing acid strength. As acid is kept constant at each concentration for several minutes, dissolution of a specific calcite phase can take place. Initial results show that it is most effective to slightly crush the tests without applying standard cleaning procedures, but rather analyze them without cleaning. Samples were selected from the South Atlantic (core tops and specific downcore samples) and the Mediteterranean (plankton tows). All samples were chosen based on previous work on them to provide comparison with routinely analysed Mg/Ca ratios. The South Atlantic samples have been analyzed extensively as bulk samples separated in difference size fractions and crusted vs. non-crusted (Groeneveld and Chiessi). The Mediterranean samples were not only analyzed as bulk samples but also by Laser Ablation ICP-MS (von Raden et al.). Results show that bulk analyses are reliably reproduced by the FT method, especially for samples which are dominated by crusted calcite. Samples which were uncrusted often gave much higher Mg/Ca ratios than the bulk analyses. These higher Mg/Ca ratios mainly occur in the plankton tow samples and were also identified with Laser Ablation ICP-MS. A possible reason for this could be the presence of a high Mg amorphous calcite layer on the outside of foraminifer tests which have not completed their calcification yet as was recently also pointed out in several other studies. Identification of the crusted and uncrusted phases, and therewith a thermal gradient, seems to give the expected differences but a more rigorous statistical treatment is needed to pinpoint singular dissolution phases.

Description: Earthquake history shows that the Sunda subduction zone of the Indonesian margin produces great earthquakes offshore Sumatra, whereas earthquakes of comparable magnitude are lacking offshore Java and the Lesser Sunda islands. Morphological structures from multibeam bathymetric data across the forearc relate with the extent of the seismogenic zone (SZ). Off Java and the Lesser Sunda islands the Indo-Australian plate subducts almost normal underneath the oceanic plate of the Indonesian archipelago. Landward of the trench, the outer wedge of the slope break is ~50 km uniformly wide with uniform bathymetric gradients. The slope of the outer wedge is locally cut by one/two steeper ridges of ~5 km extent. The sharp slope break corresponds to the updip limit of the SZ, which is also associated with the seawardmost part of the outer arc high. Landward of the slope break we find narrow, uniform outer arc ridges. The landward termination of these ridges coincides with the downdip limit of the SZ. The intersection of the shallow upper plate mantle with the subduction thrust fault marks the downdip limit of the SZ beneath the forearc. Off Sumatra the Indo-Australian plate subducts obliquely underneath the continental part of the Indonesian Sunda margin. Landward of the trench, the outer wedge varies, being mostly ~70 km wide, in some areas narrowing to 50 km width. The lower slope bathymetric gradients are steep. The outer wedge slope is made up of several steeper ridges of ~5 km extent. The slope break is only locally sharp, and corresponds to the updip limit of the SZ. The outer arc ridges off Sumatra are, in comparison with the forearc structures off Java and the Lesser Sunda islands, wider and partly elevated above sea level forming the Mentawai forearc islands. The downdip limit of the SZ coincides with the intersection of a deeper upper plate mantle with the subduction thrust fault beneath the forearc. Sunda Strait marks a transition zone between the Sumatra and Java margins. Seafloor morphology enables the identification of the seismogenic zone (SZ) across the entire Sunda margin. The SZ is uniformly wide for the Sumatra margin and narrows off Sunda Strait. Sunda Strait is the transition between the Sumatra margin and the uniformly narrow extent of the SZ of the Java/Lesser Sunda margin. Comparing the Java and Lesser Sunda islands with the Sumatra margin we find the differences along the Sunda margin, especially the wider extent of the SZ off Sumatra, producing larger earthquakes, to result from the combination of various causes: The sediment income on the oceanic incoming plate and the subduction direction; we attribute a major role to the continental/oceanic upper plate nature of Sumatra/Java influencing the composition and deformation style along the forearc and subduction fault. Off Sumatra the SZ is up to more than twice as wide as off Java/Lesser Sunda islands, enlarging the unstable regime off Sumatra and thus the risk of sudden stress release in a great earthquake.

Description: EGU2011-407 The spatial and temporal distribution of sea ice in the subpolar North Atlantic is mainly controlled by the advection of warm Atlantic Water via the Norwegian and West Spitsbergen Current in eastern Fram Strait. Simultaneously, polar water and sea ice from the Arctic Ocean is transported southward by the East Greenland Current. Hence, variations in the strength of this oceanic circulation regime may either stimulate or reduce the sea ice extent. Based on organic geochemical studies of a high-resolution sediment core from eastern Fram Strait we provide new evidence for the highly variable character of the sea ice conditions in this area. The combination of the sea ice proxy IP25 (Belt et al., 2007) with phytoplankton derived biomarkers (e.g. brassicasterol, dinosterol; Volkman 2006) enables a reliable reconstruction of sea surface and sea ice conditions, respectively (Müller et al., 2009; 2010). By means of these biomarkers, we trace gradually increasing sea ice occurrences from the Mid to the Late Holocene – consistent with the neoglacial cooling trend. Throughout the past ca. 3,000 years (BP) we observe a significant short-term variability in the biomarker records, which points to rapid advances and retreats of the sea ice cover at the continental margin of West Spitsbergen. The co-occurrence of IP25 and phytoplankton markers, however, suggests that the primary productivity benefits from these sea ice surges. As such, higher amounts of open-water phytoplankton biomarkers together with peak abundances of IP25 indicate recurring periods of enhanced ice-edge phytoplankton blooms at the core site. To what extent a seesawing of temperate Atlantic Water may account for these sea ice fluctuations requires further investigation. Concurrent variations in Siberian river discharge (Stein et al., 2004) or Norwegian glacier extents (Nesje et al., 2001), however, strengthen that these fluctuations may be assigned to variations in the North Atlantic/Arctic Oscillation (NAO/AO) and (hence) a weakened/accelerated Atlantic Water input and Arctic sea ice export.

Description: The upwelling area in the eastern equatorial Pacific off Peru is one of the most pronounced oxygen minimum zones (OMZs) of the modern ocean. Modeling scenarios predict an expansion of the OMZs in the course of global change in the coming decades. As a consequence, the Peruvian continental margin represents a key locality for studies on biogeochemical dynamics in the future ocean. We present pore water and sediment data for redox-sensitive metals (Fe, Mn, V, Mo, and U) that have been collected along a transect across the Peruvian margin at 11°S. The results are used to evaluate the behavior of trace metals in a wide range of biogeochemical and hydrodynamic settings. In the core of the OMZ, where permanently anoxic conditions prevail, redox sensitive metals exhibit diagenetic behaviors largely consistent with previous studies. Vanadium and Mo are released from Fe oxihydroxides and subsequently recycled through diffusion across the benthic boundary or trapped through formation of authigenic V phases and sequestration of Mo by authigenic pyrite. Some U is delivered through diffusion across the benthic boundary, reduction and precipitation of UO2 and incorporation into phosphorites. The utmost part of the buried U, however, is delivered in particulate form, most likely as bioauthigenic U which cannot be recycled in the suboxic waters overlying the anoxic sediments. In contrast to sediments in the core of the OMZ, sediments on the shelf experience frequent oxygenation episodes related to the passage of internal waves and the regular recurrence of El Niño events. These oxygenation episodes lead to the re-oxidation and remobilization of authigenic U and V. In contrast to that, the authigenic accumulation of Mo is favored by the occasional occurrence of slightly oxidizing conditions. This is most likely due to enhanced formation of sulfur intermediates necessary for pyrite formation and the increased stability of pyrite, the major Mo sink, under oxidizing conditions, compared to authigenic V and U phases. Redox oscillations in the Peruvian OMZ thus lead to a discrimination of U against Mo, a mechanism that should be considered in the interpretation of U/Mo systematics in paleo redox studies. Overall our results provide valuable constraints on how trace metal inventories of marginal sediments may respond to expanding shelf anoxia and to short term perturbations of sediment redox conditions.

Description: EGU2010-2934 The exploration of the arctic seas require an integrated approach applying different infrastructures. In Fall 2009 German and Russian scientists performed a geo marine cruise off Kamchatka and in the western Bering Sea within the frame of the KALMAR-Project. Two main research subjects formed the scientific backbone of the cruise: The first objective focuses on the geodynamic and volcanological magmatic development of the Kuril-Kamchatka island arc system and the Kamchatka Aleutean Islands Triple-Junction. Very little is known about the composition of the mantle and the oceanic crust as well as of the seamounts including their ages. The best studied site is the Volcanologist’s Massif located between the Bering- and the Alpha Fracture Zone (Tsvetkov 1990, Volynets et al. 1992, Yogodzinsky et al. 1994), which structurally belongs to the Komandorsky Basin. The oldest rocks of the Volcanologist´s Massif show very similar trace element and isotope signatures like those rocks cropping out in the volcanoes on Kamchatka in the prolongation of the Alpha Fracture Zone (Portnyagin et al. 2005a), indicating similar conditions of magma formation. The top of the Volcanologist´s Massif is characterized by the young (〈 0.5 Ma) and hydrothermally active Piip volcano, which consists of special magnesium rich andesites ("Piip type"). Another hot site was the Meiji-Seamount which is the northernmost Seamount of the hotspot spur of the Hawaii-Emperor-Seamount chain, having an age of probably 〉 85 Ma. The only existing basement rocks from this seamount were gained during DSDP Leg 19. These are basalts with MORB like trace element and isotope signatures (Keller et al. 2000, Regelous et al. 2003). These data indicate that the Hawaii-Hotspot was at a MOR in Cretaceous time and that large volumes of depleted mantle material played a´role during the magma formation. The second objective focuses on paleo-oceanographic investigations concentrating on the sediments along the eastern continental slope of Kamchatka, in the Komandorsky Basin, and on the Shirshov Ridge in order to explore paleoclimate archives to better understand the subpolar water mass transfer and the oceanographic and climatic development in the subarctic NW-Pacific. Comparisons of Late Pleistocene and Holocene temperature changes within the near surface water masses between the NW-Pacific and the N-Atlantic resulted in a new hypothesis, the "Atlantic-Pacific seesaw" (Kiefer et al. 2001, Kim et al. 2004, Kiefer and Kienast, 2005). This Atlantic-Pacific pattern of opposite temperature variations dominates the last 60ka on millennial timescales. Modelling results of Saenko et al. (2004) support the hypothesis of the "Atlantic-Pacific seesaw" and they postulate a mechanistic connection between the two regions driven by salinity variations, which couples both regions through the thermohaline circulation. A different model relates the Holocene Atlantic-Pacific dipole to the atmospheric tele-connection between the Arctic Oscillation/N-Atlantic Oscillation and the Pacific N-American Oscillation (Kim et al. 2004). http://kalmar.ifm-geomar.de

Description: EGU2011-12780 A temporary passive seismic network of 31 broad-band stations was deployed in the region around Talca and Constitución between 35°S to 36°S latitude and 71°W to 72.5°W longitude. The network was operated between March and October 2008. Thus, we recorded data prior the magnitude Mw=8.8 earthquake of 27 February 2010 at a latitude of the major slip and surface uplift. The experiment was conducted to address fundamental questions on deformation processes, crustal and mantle structures, and fluid flow. We present results of a teleseismic P receiver function study that covers the coastal region and reaches to the Andes. The aim is to determine the structure and thickness of the continental crust and constrain the state of hydration of the mantle wedge. The P-wave receiver function technique requires large teleseismic earthquakes from different distances and backazimuths. A few percent of the incident P-wave energy from a teleseismic event will be converted into S-wave (Ps) at significant and relatively sharp discontinuities beneath the station. A small converted S phase is produced that arrives at the station within the P wave coda directly after the direct P-wave. The converted Ps phase and their crustal multiples contain information about crustal properties, such as Moho depth and the crustal vp/vs ratio. We use teleseismic events with magnitudes mb 〉 5.5 at epicentral distances between 30° and 95° to examine P-to-S converted seismic phases. Our preliminary results provide new information about the thickness of the continental crust beneath the coastal region in Central Chile. At most of the stations we observed significant energy from P to S converted waves between 4 and 5 s after the direct P-wave within a positive phase interpreted as the Moho, occurring at 35 to 40 km. The great Maule earthquake of 27 February 2010 nucleated up-dip of the continental Moho. The rupture of this earthquake seems to have propagated down-dip of the Moho. The Moho reflection show a positive polarity, indicating that the mantle is either dry or only moderately hydrated. We observed converted energy from an intracrustal boundary at around 2 s that disappears near the coast. Further, positive polarity peaks occur that are possibly caused by the down going plate.

Description: EGU2011-8738 At present, the Arctic is responding faster to global warming than most other areas on earth, as indicated by rising air temperatures, melting glaciers and ice sheets and a decline of the sea ice cover. As part of the meridional overturning circulation which connects all ocean basins and influences global climate, northward flowing Atlantic Water is the major means of heat and salt advection towards the Arctic where it strongly affects the sea ice distribution. Records of its natural variability are critical for the understanding of feedback mechanisms and the future of the Arctic climate system, but continuous historical records reach back only ca. 150 years. To reconstruct the history of temperature variations in the Fram Strait Branch of the Atlantic Current we analyzed a marine sediment core from the western Svalbard margin. In multidecadal resolution the Atlantic Water temperature record derived from planktic foraminifer associations and Mg/Ca measurements shows variations corresponding to the well-known climatic periods of the last millennium (Medieval Climate Anomaly, Little Ice Age, Modern/Industrial Period). We find that prior to the beginning of atmospheric CO2 rise at ca. 1850 A.D. average summer temperatures in the uppermost Atlantic Water entering the Arctic Ocean were in the range of 3-4.5°C. Within the 20th century, however, temperatures rose by ca. 2°C and eventually reached the modern level of ca. 6°C. Such values are unprecedented in the 1000 years before and are presumably linked to the Arctic Amplification of global warming. Taking into account the ongoing rise of global temperatures, further warming of inflowing Atlantic Water is expected to have a profound influence on sea ice and air temperatures in the Arctic.

Description: EGU2011-1847-3 Lake Van is a lake by volume of 607 km3 and a maximum depth of 450 meters in a tectonically active zone in eastern Anatolia, Turkey. In summer 2010, Lake Van was the target of a deep drilling campaign (PaleoVan) in the frame of ICDP (International Continental Scientific Drilling Program). Two sites were drilled based on reflection seismic data collected during a seismic campaign in 2004. Here we present a first joint interpretation of the seismic and drilling data. Interpretation of seismic reflection data from lake reveals three physiographic provinces: a lacustrine shelf, a lacustrine slope, and a deep, relatively flat lake basin. The most prominent features of the lacustrine shelf and slope are prograding deltaic sequences, numerous unconformities, submerged channels, as well as closely spaced U- and/or V-shaped depressions, reflecting the variable lake level history of Lake Van. The seismic units of the shelf are dominantly composed of low-to-good continuity, variable amplitude reflections interpreted as fluvial deposits. The lake consists of three prominent basins (Tatvan, Deveboynu, and Northern Basins), separated by basement highs or ridges (Ahlat Ridge). The seismic units corresponding to these basins mainly consist of low to very high amplitude, well-stratified reflection patterns. Chaotic reflections are seen in parts of these basins. The Deveboynu Basin consists mainly of chaotic reflections. The Tatvan and Northern Basins are characterized by an alternating succession of well-stratified and chaotic reflecting layers. The chaotic seismic facies are interpreted as slump and slide deposits, which are probably the result of quick lake level fluctuations and/or earthquakes. The moderateto high amplitude, well-stratified facies seen in the deep parts of the basins away from the terrigenous sediment sources are interpreted as lacustrine deposits and tephra layers. The total sediment thickness in the deep parts of the lake is over 400 m. Prominent clinoforms indicate the initial flooding of Lake Van about 500 ka ago. The acoustic basement and the sediments lying on top of the basement in the southern part of the lake are disrupted by various intrusions and extrusions suggesting active volcanism. Synthetic seismograms calculated based on core logging, wire-line logging and check shot data will allow the correlation between seismic and drill data. This approach will allow extrapolating the stratigraphy from the wells to 3D-space by using the seismic data.

Description: EGU2011-4235 The Arctic is undergoing rapid environmental and economic transformations. Recent climate warming, which is simplifying access to oil and gas resources, enabling trans Arctic shipping, and shifting the distribution of harvestable resources, has brought the Arctic Ocean to the top of national and international political agendas. Scientific knowledge of the present status of the Arctic Ocean and the process-based understanding of the mechanics of change are urgently needed to make useful predictions of future conditions throughout the Arctic region. These are required to plan for the consequences of climate change. A step towards improving our capacity to predict future Arctic change was undertaken with the Second International Conference on Arctic Research Planning (ICARP II) meetings in 2005 and 2006, which brought together scientists, policymakers, research managers, Arctic residents, and other stakeholders interested in the future of the Arctic region. The Arctic in Rapid Transition (ART) Initiative developed out of the synthesis of the several resulting ICARP II science plans specific to the marine environment. This process started in October 2008 and has been driven by early career scientists. The ART Initiative is an integrative, international, multi-disciplinary, long-term pan-Arctic network to study changes and feedbacks with respect to physical characteristics and biogeochemical cycles in the Arctic Ocean in a state of rapid transition and its impact on the biological production. The first ART workshop was held in Fairbanks, Alaska, in November 2009 with 58 participants from 9 countries. Workshop discussions and reports were used to develop a science plan that integrates, updates, and develops priorities for Arctic Marine Science over the next decade. The science plan was accepted and approved by the International Arctic Science Committee (IASC) Marine Group, the former Arctic Ocean Science Board. The second ART workshop was held in Winnipeg, Canada, in October 2010 with 20 participants from 7 countries to develop the implementation plan. Our focus within the ART Initiative will be to bridge gaps in knowledge not only across disciplinary boundaries (e.g., biology, geochemistry, geology, meteorology, physical oceanography), but also across geographic (e.g., international boundaries, shelves, margins, and the central Arctic Ocean) and temporal boundaries (e.g., alaeo/geologic records, current process observations, and future modeling studies). This approach of the ART Initiative will provide a means to better understand and predict change, particularly the consequences for biological productivity, and ultimate responses in the Arctic Ocean system. More information about the ART Initiative can be found at http://aosb.arcticportal.org/art.html.

Description: EGU2011-13199 Images of crustal construction provide a key to understand the interplay of magmatism and tectonism while oceanic crust is build up. Bathymetric data show that the crustal construction is highly variable. Areas that are dominated by magmatic processes are adjacent to areas that are highly tectonised and where mantle rocks were found. The Mid-Atlantic Ridge at 22°N shows this high variability along the ridge axis, within the TAMMAR segment, and from segment to segment. However, this strong variability occurs also off-axis, spreading parallel, representing different times in the same area of the ridge. A fracture zone, with limited magma supply, has been replaced by a segment centre with a high magmatic budget. Roughly 4.5 million years ago, the growing magmatic active TAMMAR segment, propagated into the fracture zone, started the migration of the ridge offset to the south, and stopped the formation of core complexes. We present data from seismic refraction and wide-angle reflection profiles that surveyed the crustal structure across the ridge crest of the TAMMAR segment. These yield the crustal structure at the segment centre as a function of melt supply. The results suggest that crust is ~8 km thick near the ridge and decreases in thickness with offset to the ridge axis. Seismic layer 3 shows profound changes in thickness and becomes rapidly one kilometre thicker approx. 5 million years ago. This correlates with gravimetric data and the observed “Bull’s eye” anomaly in that region. Our observations support a temporal change from thick lithosphere with oceanic core complex formation to thin lithosphere with focussed mantle upwelling and segment growing. The formation of ‘thick-crust’ volcanic centre seems to have coincided with the onset of propagation 4.5 million years ago.

Description: EGU2011-6081 Natural stable isotopes are a powerful tool in marine sciences to investigate biological processes, such as present and past nutrient utilization. In this study we present the first dissolved silicon isotope data in the upwelling area off Peru, where one of the world’s largest Oxygen Minimum Zones (OMZ) is located. Silicon is the most important component required for phytoplankton (diatom) growth, which dominates primary productivity in this region. Stable Si isotopes are fractionated during diatom growth in that the lighter Si isotopes are preferentially incorporated into diatoms with a fractionation factor of -1.1 promille. The Si isotope composition of dissolved silicic acid of the corresponding surface waters is therefore left isotopically heavier. The Si isotope composition, 30Si/28Si, is expressed as δ30Si values, which stand forh deviations from a given standard (NBS28). Investigation of the dissolved seawater Si isotope composition thus provides a measure for the utilization and, combined with information on the Si isotope composition of the water masses upwelling off Peru, it is a measure for the supply pathways of Si to the coastal upwelling centres. Surface waters on the shelf off Peru are mainly fed by the Equatorial Undercurrent, which mainly consists of waters originating from the western and Central Pacific and which has a characteristic δ30Si of +1.5 promille. In areas and during phases of intense upwelling the fractionation of Si isotopes was observed to be weaker due to upwelling-driven supply of less fractionated Si (δ30Si = 1.7 promille, from water depths of around 100-150 m, whereas under weak upwelling conditions fractionation is higher (δ30Si ~3 promille due to a more complete utilization of the available dissolved silicate. The distribution of dissolved δ30Si correlates strongly with particulate biogenic silicate (opal) concentrations in that highest opal concentrations in the surface waters show the lowest δ30Si values thus strongest upwelling intensity. The most extreme δ30Si values in surface waters (δ30Si = 4.5 promille are observed offshore where silicic acid concentrations are nearly zero. Furthermore we compare the δ30Si data with the dissolved nitrogen isotope distribution, which in addition to nitrate utilization is mainly influenced by denitrification and annamox processes in the OMZ. Combined silicon and nitrogen isotope compositions can thus help to disentangle different fractionation processes within the nitrogen cycle.

Description: EGU2011-2455 The current interglacial has gone through a variety of warmer and colder periods. Consistent with the decreasing solar insolation during the Holocene, warmest conditions have occurred particularly within its earliest phase. We studied high-resolution sediment sequences from the Western Svalbard margin covering the last ca 10,000 years in order to reconstruct the variations of Atlantic Water advection to the Arctic, the sea ice extent, and the structure of the water column on the Westspitsbergen continental margin. The Fram Strait, often referred to as the Arctic Gateway, is the only deep-water passage for Atlantic-derived water masses to enter the Arctic Ocean. Northward advection of relatively warm and saline Atlantic Water masses keeps the eastern part of the Fram Strait ice-free all year. It therefore plays a crucial role for the heat budget of the Arctic. A multiproxy data set including geochemical, micropaleontological, and sedimentological parameters was established with centennial to multidecadal time resolution. Records of foraminiferal oxygen and carbon isotopes, planktic foraminifer assemblages, and the amount of ice rafted debris clearly reveal distinct variations between climatically warmer and colder intervals throughout this period. Planktic foraminifer assemblages reveal warmest conditions for the early Holocene period (ca 10-8 ka). A second warming pulse is detected between 5 and 6 ka. In the second half of the Holocene, increased IRD contents are indicative of a significant cooling trend. Despite of the decreasing solar insolation planktic foraminiferal assemblages suggest a return of slightly strengthened Atlantic Water advection around 3 to 2 ka and a strong warming event in the present, anthropogenically influenced period.

Description: We determined the isotopic composition of neodymium (Nd) and lead (Pb) of past seawater to reconstruct water mass exchange and erosional input between the Arctic Ocean and the Norwegian-Greenland Seas over the past 5 Ma. For this purpose, sediments of ODP site 911 (leg 151) located at 900 m water depth on the Yermak Plateau in the Fram Strait were used. The paleo-seawater variability of Nd and Pb isotopes was extracted from the sea water-derived metal oxide coatings on the sediment particles following the leaching method of Gutjahr et al. (2007). All radiogenic isotope data were acquired by Multi-Collector (MC) ICP-MS. The site 911 stratigraphy of Knies et al. (2009) was applied. Surface sediment Sr and Nd isotope data, as well as downcore Sr isotope data obtained on the same leaches are close to seawater and confirm the seawater origin of the Nd and Pb isotope signatures. The deep water Nd isotope time series extracted from site 911 was in general more radiogenic ("Nd = -7.5 to -10) than present day deep water ("Nd = -9.8 to -11.8) in the area of the Fram Strait (Andersson et al., 2008) and does not show a systematic trend with time. In contrast, the radiogenic isotope composition of Pb evolved from 206Pb/204Pb ratios around 18.7 to more radiogenic values around 19.2 between 2 Ma and today. The data indicate that mixing of water masses from the Arctic Ocean and the Norwegian-Greenland Seas has controlled the Nd isotope signatures of deep waters on the Yermak Plateau over the past 5 Ma. Prior to 1.7 Ma the Nd isotope signatures on the Yermak Plateau were less radiogenic than waters from the same depth in the central Arctic Ocean (Haley et al., 2008) pointing to a greater influence from the Norwegian-Greenland Seas. After 1.7 Ma the central Arctic and Yermak Plateau data have varied around similar values indicating water mass mixing overall similar to today. In contrast, the Pb isotope composition of deep waters in the Fram Strait appears to have been dominated by weathering inputs from glacially weathering old continental landmasses, such as Greenland or parts of Svalbard since 2 Ma. A similar control over the Pb isotope evolution of seawater since the onset of Northern Hemisphere Glaciation was recorded by ferromanganese crusts that grew from North Atlantic DeepWater in the western North Atlantic. References: Gutjahr, M., Frank, M., Stirling, C.H., Klemm, V., van de Flierdt, T. and Halliday, A.N. (2007): Reliable extraction of a deepwater trace metal isotope signal from Fe-Mn oxyhydroxide coatings of marine sediments.- Chemical Geology 242, 351-370 Haley B. A., M. Frank, R.F. Spielhagen and A. Eisenhauer (2008): Influence of brine formation on Arctic Ocean circulation over the past 15 million years. Nature Geoscience 1, 68–72 Andersson, P.S., Porcelli, D., Frank, M., Björk, G., Dahlqvist, R. and Gustafsson, Ö. (2008): Neodymium isotopes in seawater from the Barents Sea and Fram Strait Arctic- Atlantic gateways.- Geochim. Cosmochim. Acta 72, 2854-2867 Knies, J., J. Matthiessen, C. Vogt, J.S. Laberg, B.O. Hjelstuen, M.Smelror, E. Larsen, K. Andreassen, T. Eidvin and T.O. Vorren (2009): The Plio-Pleistocene glaciation of the Barents Sea–Svalbard region: a new model based on revised chronostratigraphy - Quaternary Science Reviews 28, 9-10, 812-829

Description: In March 2007 the sea floor drill rig MeBo (short for “Meeresboden-Bohrgerät”, ‘sea floor drill rig’ in German) returned from a 17-day scientific cruise with the new German research vessel Maria S. Merian. Four sites between 350 m and 1700 m water depth were sampled at the continental slope off Morocco by push coring and rotary drilling. Up to 41.5-m-long sediment cores were recovered from Miocene, Pliocene, and Pleistocene marls. MeBo bridges the gap between conventional sampling methods from standard multipurpose research vessels (gravity corer, piston corer, dredges) and drill ships. Most bigger research vessels will be able to support deployment of the MeBo. Since the drill system can be easily transported within 20-ft containers, worldwide operation from vessels of opportunity is possible. With the MeBo a new system is available for marine geosciences that allows the recovery of high quality samples from soft sediments and hard rock from the deep sea without relying on the services of expensive drilling vessels.

Description: Heinrich events are among the dominant modes of glacial climate variability. During these events, massive iceberg armadas were released by the Laurentide Ice Sheet, sailed across the Atlantic, and caused large-scale climate changes. We study these events in a fully coupled complex ice sheet–climate model with synchronous coupling between ice sheets and oceans. The ice discharges occur as internal variability of the model with a recurrence period of 5kyr, an event duration of 1–1.5kyr, and a peak discharge rate of about 50mSv, roughly consistent with reconstructions. The climate response shows a two-stage behavior, with freshwater release effects dominating the surge phase and ice-sheet elevation effects dominating in the post-surge phase. As a direct response to the freshwater discharge during the surge phase, the deepwater formation in the North Atlantic decreases and the North Atlantic deepwater cell weakens by 3.5Sv. With the reduced oceanic heat transport, the surface temperatures across the North Atlantic decrease, and the associated reduction in evaporation causes a drying in Europe. The ice discharge lowers the surface elevation in the Hudson Bay area and thus leads to increased precipitation and accelerated ice sheet regrowth in the post-surge phase. Furthermore, the jet stream widens to the north and becomes more zonal. This contributes to a weakening of the subpolar gyre, and a continued cooling over Europe even after the ice discharge. This two-stage behavior can explain previously contradicting model results and understandings of Heinrich Events.

Description: In this study, we propose a statistical method to validate sea-level reconstructions using geological records known as sea-level indicators (SLIs). SLIs are often the only available data to retrace late-glacial relative sea level (RSL). Determining the RSL from SLI height is not straight forward, the elevation at which an SLI was found usually does not represent the past RSL. In contrast, it has to be related to past RSL by investigating sample’s type, habitat and deposition conditions. For instance, water distribution at which a specific specimen is found today can be related to the indicator's depositional height range. Furthermore, the precision of dating varies between geological samples, and, in case of radiocarbon dating, the age has to be calibrated using a non-linear calibration curve. To avoid an a-priori assumption like normal-distributed uncertainties, we define likelihood functions which take into account the indicative meaning’s available error information and calibration statistics represented by joint probabilities. For this conceptional study, we restrict ourselves to one type of indicators, shallow-water shells, which are usually considered as low-grade samples giving only a lower limit of former sea level, as the depth range in which they live spreads over several tens of meters, and does not follow a normal distribution. The presented method is aimed to serve as a strategy for glacial isostatic adjustment reconstructions, in this case for the German Paleo-Climate Modelling Initiative PalMod (https://www.palmod.de/en) and by extending it to other SLI types.

Description: Dynamic vegetation models simulate global vegetation in terms of fractional coverages of a few plant functional types (PFTs). Although these models often share the same concept, they differ with respect to the number and kind of PFTs, complicating the comparability of simulated vegetation distributions. Pollen-based reconstructions are initially only available in form of time-series of individual taxa that are not distinguished in the models. Thus, to evaluate simulated vegetation distributions, the modelling results and pollen-based reconstructions have to be converted into a comparable format. The classical approach is the method of biomisation, but hitherto, PFT-based biomisation methods were only available for individual models. We introduce and evaluate a simple, universally applicable technique to harmonize PFT-distributions by assigning them into nine mega-biomes that follow the definitions commonly used for vegetation reconstructions. The method works well for all state-of the art dynamic vegetation models, independent of the spatial resolution or the complexity of the models. Large biome belts (such as tropical forest) are well represented, but regionally confined biomes (warm-mixed forest, Savanna) are only partly captured. Overall, the PFT-based biomisation is able to keep up with the conventional biomisation approach of forcing biome models (here: BIOME1) with the background climate states. The new method has, however, the advantage that it allows a more direct comparison and evaluation of the vegetation distributions simulated by Earth System Models. Thereby, the new method provides a powerful tool for the evaluation of Earth System Models in general.

Description: We have developed a new module to calculate soil organic carbon (SOC) accumulation in perennially frozen ground in the land surface model JSBACH. Running this offline version of MPI-ESM we have modelled permafrost carbon accumulation and release from the Last Glacial Maximum (LGM) to the Pre-industrial (PI). Our simulated near-surface PI permafrost extent of 16.9Miokm2 is close to observational evidence. Glacial boundary conditions, especially ice sheet coverage, result in profoundly different spatial patterns of glacial permafrost extent. Deglacial warming leads to large-scale changes in soil temperatures, manifested in permafrost disappearance in southerly regions, and permafrost aggregation in formerly glaciated grid cells. In contrast to the large spatial shift in simulated permafrost occurrence, we infer an only moderate increase of total LGM permafrost area (18.3Miokm2) – together with pronounced changes in the depth of seasonal thaw. Reconstructions suggest a larger spread of glacial permafrost towards more southerly regions, but with a highly uncertain extent of non-continuous permafrost. Compared to a control simulation without describing the transport of SOC into perennially frozen ground, the implementation of our newly developed module for simulating permafrost SOC accumulation leads to a doubling of simulated LGM permafrost SOC storage (amounting to a total of ~150PgC). Despite LGM temperatures favouring a larger permafrost extent, simulated cold glacial temperatures – together with low precipitation and low CO2 levels – limit vegetation productivity and therefore prevent a larger glacial SOC build-up in our model. Changes in physical and biogeochemical boundary conditions during deglacial warming lead to an increase in mineral SOC storage towards the Holocene (168PgC at PI), which is below observational estimates (575PgC in continuous and discontinuous permafrost). Additional model experiments clarified the sensitivity of simulated SOC storage to model parameters, affecting long-term soil carbon respiration rates and simulated active layer depths. Rather than a steady increase in carbon release from the LGM to PI as a consequence of deglacial permafrost degradation, our results suggest alternating phases of soil carbon accumulation and loss as an effect of dynamic changes in permafrost extent, active layer depths, soil litter input, and heterotrophic respiration.

Description: Climate reconstructions based on proxy records recovered from marine sediments, such as alkenone records or geochemical parameters measured on foraminifera, play an important role in our understanding of the climate system. They provide information about the state of the ocean ranging back hundreds to millions of years and form the backbone of paleo-oceanography. However, there are many sources of uncertainty associated with the signal recovered from sediment archived proxies. These include seasonal or depth habitat biases in the recorded signal, a frequency dependent reduction in the amplitude of the recorded signal due to bioturbation of the sediment, aliasing of high frequency climate variation onto a nominally annual, decadal or centennial resolution signal, and additional sample processing and measurement error introduced when the proxy signal is recovered. Here we present a forward model for sediment archived proxies that jointly models the above processes, so that the magnitude of their separate and combined effects can be investigated. Applications include the interpretation and analysis of uncertainty in existing proxy records, parameter sensitivity analysis to optimize future studies, and the generation of pseudo-proxy records that can be used to test reconstruction methods. We provide examples, such as the simulation of individual foraminifera records, that demonstrate the usefulness of the forward model for paleoclimate studies. The model is implemented as a user-friendly R package, sedproxy, the use of which we hope will contribute to a better understanding of both the limitations and potential of marine sediment proxies to inform about past climate.

Description: EGU2011-12864 The Woodlark Basin east of Papua New Guinea represents one of the few places on Earth where a spreading axis propagates into continental crust. This special tectonic setting allows insights into the evolution of magma composition as continental extension and break-up changes to the formation of ocean crust. We report here geochemical results on samples collected in 2009 from the four segments closest to the continental breakup, from segment 1 which abuts the detachment fault responsible for continental extension on Moresby Seamount in the West, to segment 4, representing mature oceanic crust in the East. A total of 208 glass samples have been analyzed for their major (EMPA) and trace element (LA-ICPMS) compositions. The data show strong E-W variations. Samples ranging from tholeiitic basalt and basaltic andesite to andesite and rhyolite are found on Segment 1. They have generally high alkali values and a wide range of trace element contents and ratios. Segments 2 to 4 magmas in contrast only comprise tholeiitic basalt with lower alkali contents and a more restricted range of trace element chemistry. The geochemical differences between the segments cannot be attributed to differentiation processes alone, and different sources are required. High Ba/La, (La/Sm)N, Rb/Sr, and Th/La on Segment 1 suggest a derivation from an enriched mantle source, while low Nd/Pb and Nb/U suggest that some of the enrichment may also reflect the influence of continental crust during magma genesis. Whether this continental signature is present in the form of recycled material in the mantle or as rafted continental blocks in the axial region is at present unclear. In contrast to rocks from segment 1, trace element compositions of volcanic glasses from segments 2 to 4 show a stronger MORB signature, presumably reflecting more mature spreading in this part of the basin. The influence of continental material appears to be minimal, suggesting that uncontaminated asthenosphere quickly flows into the rift and/or that continental blocks are not retained in the axial region for long time periods following the rifting-spreading transition.

Description: The Global Ocean Data Analysis Project (GLODAP) is a synthesis effort providing regular compilations of surface-to-bottom ocean biogeochemical bottle data, with an emphasis on seawater inorganic carbon chemistry and related variables determined through chemical analysis of seawater samples. GLODAPv2.2022 is an update of the previous version, GLODAPv2.2021 (Lauvset et al., 2021). The major changes are as follows: data from 96 new cruises were added, data coverage was extended until 2021, and for the first time we performed secondary quality control on all sulphur hexafluoride (SF6) data. In addition, a number of changes were made to data included in GLODAPv2.2021. These changes affect specifically the SF6 data, which are now subjected to secondary quality control, and carbon data measured onboard the RV Knorr in the Indian Ocean in 1994–1995 which are now adjusted using CRM measurements made at the time. GLODAPv2.2022 includes measurements from almost 1.4 million water samples from the global oceans collected on 1085 cruises. The data for the now 13 GLODAP core variables (salinity, oxygen, nitrate, silicate, phosphate, dissolved inorganic carbon, total alkalinity, pH, CFC-11, CFC-12, CFC-113, CCl4, and SF6) have undergone extensive quality control with a focus on systematic evaluation of bias. The data are available in two formats: (i) as submitted by the data originator but converted to World Ocean Circulation Experiment (WOCE) exchange format and (ii) as a merged data product with adjustments applied to minimize bias. For the present annual update, adjustments for the 96 new cruises were derived by comparing those data with the data from the 989 quality controlled cruises in the GLODAPv2.2021 data product using crossover analysis. SF6 data from all cruises were evaluated by comparison with CFC-12 data measured on the same cruises. For nutrients and ocean carbon dioxide (CO2) chemistry comparisons to estimates based on empirical algorithms provided additional context for adjustment decisions. The adjustments that we applied are intended to remove potential biases from errors related to measurement, calibration, and data handling practices without removing known or likely time trends or variations in the variables evaluated. The compiled and adjusted data product is believed to be consistent to better than 0.005 in salinity, 1 % in oxygen, 2 % in nitrate, 2 % in silicate, 2 % in phosphate, 4 μmol kg-1 in dissolved inorganic carbon, 4 μmol kg-1 in total alkalinity, 0.01–0.02 in pH (depending on region), and 5 % in the halogenated transient tracers. The other variables included in the compilation, such as isotopic tracers and discrete CO2 fugacity (fCO2), were not subjected to bias comparison or adjustments.

Description: Perturbations in stratospheric aerosol due to explosive volcanic eruptions are a primary contributor to natural climate variability. Observations of stratospheric aerosol are available for the past decades, and information from ice cores has been used to derive estimates of stratospheric sulfur injections and aerosol optical depth over the Holocene (approximately 10 000 BP to present) and into the last glacial period, extending back to 60 000 BP. Tephra records of past volcanism, compared to ice cores, are less complete but extend much further into the past. To support model studies of the potential impacts of explosive volcanism on climate variability across timescales, we present here an ensemble reconstruction of volcanic stratospheric sulfur injection (VSSI) over the last 140 000 years that is based primarily on terrestrial and marine tephra records. VSSI values are computed as a simple function of eruption magnitude based on VSSI estimates from ice cores and satellite observations for identified eruptions. To correct for the incompleteness of the tephra record, we include stochastically generated synthetic eruptions assuming a constant background eruption frequency from the ice core Holocene record. While the reconstruction often differs from ice core estimates for specific eruptions due to uncertainties in the data used and reconstruction method, it shows good agreement with an ice-core-based VSSI reconstruction in terms of millennial-scale cumulative VSSI variations over the Holocene. The PalVol reconstruction provides a new basis to test the contributions of forced vs. unforced natural variability to the spectrum of climate and the mechanisms leading to abrupt transitions in the palaeoclimate record with low- to high-complexity climate models. The PalVol volcanic forcing reconstruction is available at https://doi.org/10.26050/WDCC/PalVolv1 (Toohey and Schindlbeck-Belo, 2023).