Atlantic 231Pa/230Th and biogenic opal compilation of the Holocene and the LGM

The strength and geometry of the Atlantic meridional overturning circulation is tightly coupled to climate on glacial-interglacial and millennial timescales, but has proved difficult to reconstruct, particularly for the Last Glacial Maximum. Today, the return flow from the northern North Atlantic to lower latitudes associated with the Atlantic meridional overturning circulation reaches down to approximately 4,000 m. In contrast, during the Last Glacial Maximum this return flow is thought to have occurred primarily at shallower depths. Measurements of sedimentary 231Pa/230Th have been used to reconstruct the strength of circulation in the North Atlantic Ocean, but the effects of biogenic silica on 231Pa/230Th-based estimates remain controversial. Here we use measurements of 231Pa/230Th ratios and biogenic silica in Holocene-aged Atlantic sediments and simulations with a two-dimensional scavenging model to demonstrate that the geometry and strength of the Atlantic meridional overturning circulation are the primary controls of 231Pa/230Th ratios in modern Atlantic sediments. For the glacial maximum, a simulation of Atlantic overturning with a shallow, but vigorous circulation and bulk water transport at around 2,000 m depth best matched observed glacial Atlantic 231Pa/230Th values. We estimate that the transport of intermediate water during the Last Glacial Maximum was at least as strong as deep water transport today.

Comparison of glacial/interglacial carbonate concentrations in Atlantic sediments

To assess the regional effects of glaciation on sedimentation in the Atlantic Ocean we compare sediment types, distributions, and rates between Recent (core top) and last glacial maximum (LGM: ~18,000 years B.P.) stratigraphic levels. Based upon smear slides and carbonate analyses in 178 cores we find that glacial age carbonate content is generally lower than Recent. During both the Recent and LGM, carbonate content shows an east/west asymmetry with western basins exhibiting lower carbonate values. Input of ice-rafted detritus into the North Atlantic during LGM time interrupts this topographic control on carbonate distribution considerably farther south than at present; in the South Atlantic this effect is minor. Comparison of LGM and Recent sediment distributions indicates that the LGM seafloor was dominated by biogenic oozes, calcareous clays, and clays, while the Recent is dominated by biogenic oozes and marls. Coarse-grained detritus is much more prevalent in LGM sediments, derived not only from glacial input but also from fluvial and aeolian sources. Sedimentation rates, calculated from LGM to Recent sediment thickness in cores, are <4 cm/1000 yr for most of the ocean. Higher rates are typical of the continental margin off the Amazon River, the North American Basin, and a small region off west equatorial Africa.

Compilation of downward flux observations from sediment trap deployments in the Atlantic Ocean - Contribution to EURO-BASIN's Data integration

We provide a compilation of downward fluxes (total mass, POC, PON, BSiO2, CaCO3, PIC and lithogenic/terrigenous fluxes) from over 6000 sediment trap measurements distributed in the Atlantic Ocean, from 30 degree North to 49 degree South, and covering the period 1982-2011. Data from the Mediterranean Sea are also included. Data were compiled from different sources: data repositories (BCO-DMO, PANGAEA), time series sites (BATS, CARIACO), published scientific papers and/or personal communications from PI's. All sources are specifed in the data set. Data from the World Ocean Atlas 2009 were extracted to provide each flux observation with contextual environmental data, such as temperature, salinity, oxygen (concentration, AOU and percentage saturation), nitrate, phosphate and silicate.

(Table S1) Stable carbon and oxygen isotope ratios of benthic foraminifera of the western Atlantic Ocean

Oxygen and carbon isotopic data were produced on the benthic foraminiferal taxa Cibicidoides and Planulina from 25 new piston cores, gravity cores, and multicores from the Brazil margin. The cores span water depths from about 400 to 3000 m and intersect the major water masses in this region. These new data fill a critical gap in the South Atlantic Ocean and provide the motivation for updating the classic glacial western Atlantic d13C transect of Duplessy et al. (1988). The distribution of 13C of SumCO2 requires the presence of three distinct water masses in the glacial Atlantic Ocean: a shallow (~1000 m), southern source water mass with an end-member d13C value of about 0.3-0.5 per mil VPDB, a middepth (~1500 m), northern source water mass with an end-member value of about 1.5 per mil, and a deep (>2000 m), southern source water with an end-member value of less than -0.2 per mil, and perhaps as low as the -0.9 per mil values observed in the South Atlantic sector of the Southern Ocean (Ninnemann and Charles, 2002, doi:10.1016/S0012-821X(02)00708-2). The origins of the water masses are supported by the meridional gradients in benthic foraminiferal d18O. A revised glacial section of deep water d13C documents the positions and gradients among these end-member intermediate and deep water masses. The large property gradients in the presence of strong vertical mixing can only be maintained by a vigorous overturning circulation.

Last Interglacial synthesis of high-latitude temperature: temperature anomalies and associated errors for 4 time slices

The Last Interglacial (LIG, 129-116 thousand of years BP, ka) represents a test bed for climate model feedbacks in warmer-than-present high latitude regions. However, mainly because aligning different palaeoclimatic archives and from different parts of the world is not trivial, a spatio-temporal picture of LIG temperature changes is difficult to obtain. Here, we have selected 47 polar ice core and sub-polar marine sediment records and developed a strategy to align them onto the recent AICC2012 ice core chronology. We provide the first compilation of high-latitude temperature changes across the LIG associated with a coherent temporal framework built between ice core and marine sediment records. Our new data synthesis highlights non-synchronous maximum temperature changes between the two hemispheres with the Southern Ocean and Antarctica records showing an early warming compared to North Atlantic records. We also observe warmer than present-day conditions that occur for a longer time period in southern high latitudes than in northern high latitudes. Finally, the amplitude of temperature changes at high northern latitudes is larger compared to high southern latitude temperature changes recorded at the onset and the demise of the LIG. We have also compiled four data-based time slices with temperature anomalies (compared to present-day conditions) at 115 ka, 120 ka, 125 ka and 130 ka and quantitatively estimated temperature uncertainties that include relative dating errors. This provides an improved benchmark for performing more robust model-data comparison. The surface temperature simulated by two General Circulation Models (CCSM3 and HadCM3) for 130 ka and 125 ka is compared to the corresponding time slice data synthesis. This comparison shows that the models predict warmer than present conditions earlier than documented in the North Atlantic, while neither model is able to produce the reconstructed early Southern Ocean and Antarctic warming. Our results highlight the importance of producing a sequence of time slices rather than one single time slice averaging the LIG climate conditions.

(Supporting Information) 14Cage plateau tuning of local Atlantic sediment cores

Modeling studies predict that changes in radiocarbon (14C) reservoir ages of surface waters during the last deglacial episode will reflect changes in both atmospheric 14C concentration and ocean circulation including the Atlantic Meridional Overturning Circulation. Tests of these models require the availability of accurate 14C reservoir ages in well-dated late Quaternary time series. We here test two models using plateau-tuned 14C time series in multiple well-placed sediment core age-depth sequences throughout the lower latitudes of the Atlantic Ocean. 14C age plateau tuning in glacial and deglacial sequences provides accurate calendar year ages that differ by as much as 500-2500 years from those based on assumed global reservoir ages around 400 years. This study demonstrates increases in local Atlantic surface reservoir ages of up to 1000 years during the Last Glacial Maximum, ages that reflect stronger trades off Benguela and summer winds off southern Brazil. By contrast, surface water reservoir ages remained close to zero in the Cariaco Basin in the southern Caribbean due to lagoon-style isolation and persistently strong atmospheric CO2 exchange. Later, during the early deglacial (16 ka) reservoir ages decreased to a minimum of 170-420 14C years throughout the South Atlantic, likely in response to the rapid rise in atmospheric pCO2 and Antarctic temperatures occurring then. Changes in magnitude and geographic distribution of 14C reservoir ages of peak glacial and deglacial surface waters deviate from the results of Franke et al. (2008) but are generally consistent with those of the more advanced ocean circulation model of Butzin et al. (2012).

Spatially explicit estimates of stock size, structure and biomass of North Atlantic albacore tuna (Thunnus alalunga) in the North Atlantic for the period 1956-2010, compiled from the International Commission for the Conservation of Atlantic Tunas

The development of the ecosystem approach and models for the management of ocean marine resources requires easy access to standard validated datasets of historical catch data for the main exploited species. They are used to measure the impact of biomass removal by fisheries and to evaluate the models skills, while the use of standard dataset facilitates models inter-comparison. North Atlantic albacore tuna is exploited all year round by longline and in summer and autumn by surface fisheries and fishery statistics compiled by the International Commission for the Conservation of Atlantic Tunas (ICCAT). Catch and effort with geographical coordinates at monthly spatial resolution of 1° or 5° squares were extracted for this species with a careful definition of fisheries and data screening. In total, thirteen fisheries were defined for the period 1956-2010, with fishing gears longline, troll, mid-water trawl and bait fishing. However, the spatialized catch effort data available in ICCAT database represent a fraction of the entire total catch. Length frequencies of catch were also extracted according to the definition of fisheries above for the period 1956-2010 with a quarterly temporal resolution and spatial resolutions varying from 1°x 1° to 10°x 20°. The resolution used to measure the fish also varies with size-bins of 1, 2 or 5 cm (Fork Length). The screening of data allowed detecting inconsistencies with a relatively large number of samples larger than 150 cm while all studies on the growth of albacore suggest that fish rarely grow up over 130 cm. Therefore, a threshold value of 130 cm has been arbitrarily fixed and all length frequency data above this value removed from the original data set.

Distribution of planktic foraminifera in surface sediments of the Atlantic Ocean

We present a data set of 738 planktonic foraminiferal species counts from sediment surface samples of the eastern North Atlantic and the South Atlantic between 87°N and 40°S, 35°E and 60°W including published Climate: Long-Range Investigation, Mapping, and Prediction (CLIMAP) data. These species counts are linked to Levitus's [1982] modern water temperature data for the four caloric seasons, four depth ranges (0, 30, 50, and 75 m), and the combined means of those depth ranges. The relation between planktonic foraminiferal assemblages and sea surface temperature (SST) data is estimated using the newly developed SIMMAX technique, which is an acronym for a modern analog technique (MAT) with a similarity index, based on (1) the scalar product of the normalized faunal percentages and (2) a weighting procedure of the modern analog's SSTs according to the inverse geographical distances of the most similar samples. Compared to the classical CLIMAP transfer technique and conventional MAT techniques, SIMMAX provides a more confident reconstruction of paleo-SSTs (correlation coefficient is 0.994 for the caloric winter and 0.993 for caloric summer). The standard deviation of the residuals is 0.90°C for caloric winter and 0.96°C for caloric summer at 0-m water depth. The SST estimates reach optimum stability (standard deviation of the residuals is 0.88°C) at the average 0- to 75-m water depth. Our extensive database provides SST estimates over a range of -1.4 to 27.2°C for caloric winter and 0.4 to 28.6°C for caloric summer, allowing SST estimates which are especially valuable for the high-latitude Atlantic during glacial times.

(Table 3) Barium barite and excess comparison, accumulation rates and productivity from surface sediments

Since Dymond et al. (1992, doi:10.1029/92PA00181) proposed the paleoproductivity algorithm based on "Bio-Ba", which relies on a strong correlation between Ba and organic carbon fluxes in sediment traps, this proxy has been applied in many paleoproductivity studies. Barite, the main carrier of particulate barium in the water column and the phase associated with carbon export, has also been suggested as a reliable paleoproductivity proxy in some locations. We demonstrate that Ba(excess) (total barium minus the fraction associated with terrigenous material) frequently overestimates Ba(barite) (barium associated with the mineral barite), most likely due to the inclusion of barium from phases other than barite and terrigenous silicates (e.g., carbonate, organic matter, opal, Fe-Mn oxides, and hydroxides). A comparison between overlying oceanic carbon export and carbon export derived from Ba(excess) shows that the Dymond et al. (1992) algorithm frequently underestimates carbon export but is still a useful carbon export indicator if all caveats are considered before the algorithm is applied. Ba(barite) accumulation rates from a wide range of core top sediments from different oceanic settings are highly correlated to surface ocean 14C and Chlorophyll a measurements of primary production. This relationship varies by ocean basin, but with the application of the appropriate f ratio to 14C and Chlorophyll a primary production estimates, the plot of Ba(barite) accumulation and carbon export for the equatorial Pacific, Atlantic, and Southern Ocean converges to a global relationship that can be used to reconstruct paleo carbon export.

Geochemical and isotopic composition of Bransfield Basin Lavas collected during SONNE cruise SO155

Bransfield Basin is an actively extending marginal basin separating the inactive South Shetland arc from the northern Antarctic Peninsula. Rift-related volcanism is widespread throughout the central Bransfield Basin, but the wider eastern Bransfield Basin was previously unsampled. Lavas recovered from the eastern subbasin form three distinct groups: (1) Bransfield Group has moderate large-ion lithophile element (LILE) enrichment relative to normal mid-ocean ridge basalt (NMORB), (2) Gibbs Group has strong LILE enrichment and is restricted to a relic seamount interpreted as part of the South Shetland arc, and (3) fresh alkali basalt was recovered from the NE part of the basin near Spanish Rise. The subduction-related component in Bransfield and Gibbs Group lavas is a LILE-rich fluid with radiogenic Sr, Nd, and Pb isotope compositions derived predominantly from subducting sediment. These lavas can be modeled as melts from Pacific MORB source mantle contaminated by up to 5% of the subduction-related component. They further reveal that Pacific mantle, rather than South Atlantic mantle, has underlain Bransfield Basin since 3 Ma. Magma productivity decreases abruptly east of Bridgeman Rise, and lavas with the least subduction component outcrop at that end. Both the eastward decrease in subduction component and occurrence of young alkali basalts require that subduction-modified mantle generated during the lifetime of the South Shetland arc has been progressively removed from NE to SW. This is inconsistent with previous models suggesting continued slow subduction at the South Shetland Trench but instead favors models in which the South Scotia Ridge fault has propagated westward since 3 Ma generating transtension across the basin.