(Supplement 1) Selected chemical composition of Icelandic river waters and bedrock age

In this study we map the spatial distribution of selected dissolved constituents in Icelandic river waters using GIS methods to study and interpret the connection between river chemistry, bedrock, hydrology, vegetation and aquatic ecology. Five parameters were selected: alkalinity, SiO2, Mo, F and the dissolved inorganic nitrogen and dissolved inorganic phosphorus mole ratio (DIN/DIP). The highest concentrations were found in rivers draining young rocks within the volcanic rift zone and especially those draining active central volcanoes. However, several catchments on the margins of the rift zone also had high values for these parameters, due to geothermal influence or wetlands within their catchment area. The DIN/DIP mole ratio was higher than 16 in rivers draining old rocks, but lowest in rivers within the volcanic rift zone. Thus primary production in the rivers is limited by fixed dissolved nitrogen within the rift zone, but dissolved phosphorus in the old Tertiary catchments. Nitrogen fixation within the rift zone can be enhanced by high dissolved molybdenum concentrations in the vicinity of volcanoes. The river catchments in this study were subdivided into several hydrological categories. Importantly, the variation in the hydrology of the catchments cannot alone explain the variation in dissolved constituents. The presence or absence of central volcanoes, young reactive rocks, geothermal systems and wetlands is important for the chemistry of the river waters. We used too many categories within several of the river catchments to be able to determine a statistically significant connection between the chem¬istry of the river waters and the hydrological categories. More data are needed from rivers draining one single hydrological category. The spatial dissolved constituent distribution clearly revealed the difference between the two extremes, the young rocks of the volcanic rift zone and the old Tertiary terrain.

(Table 1) Age control points and dating methods for sediment cores of the Sea of Okhotsk

Marine- and terrestrial-derived biomarkers (alkenones, brassicasterol, dinosterol, and long-chain n-alkanes), as well as carbonate, biogenic opal, and ice-rafted debris (IRD), were measured in two sediment cores in the Sea of Okhotsk, which is located in the northwestern Pacific rim and characterized by high primary productivity. Down-core profiles of phytoplankton markers suggest that primary productivity abruptly increased during the global Meltwater Pulse events 1A (about 14 ka) and 1B (about 11 ka) and stayed high in the Holocene. Spatial and temporal distributions of the phytoplankton productivity were found to be consistent with changes in the reconstructed sea ice distribution on the basis of the IRD. This demonstrates that the progress and retreat of sea ice regulated primary productivity in the Sea of Okhotsk with minimum productivity during the glacial period. The mass accumulation rates of alkenones, CaCO3, and biogenic opal indicate that the dominant phytoplankton species during deglaciation was the coccolithophorid, Emiliania huxleyi, which was replaced by diatoms in the late Holocene. Such a phytoplankton succession was probably caused by an increase in silicate supply to the euphotic layer, possibly associated with a change in surface hydrography and/or linked to enhanced upwelling of North Pacific Deep Water.

(Table 2) Age determination of brock-red sandy mud beds and detrital carbonate beds in the Scotian margin

Piston cores from the continental margin off Nova Scotia show up to four discrete intervals of "brick-red sandy mud", which are up to 20 cm thick. The ages of these intervals are bracketed by several radiocarbon dates, and three fall in the range 12.5-14.1 ka (radiocarbon years with -0.4 kyr reservoir correction). The youngest dates from ~10.4 ka, placing it within the Younger Dryas. The distribution of the beds and their petrographic character indicate a source in the Gulf of Saint Lawrence. The grain size of these beds suggests that they comprise a coarse component transported by ice rafting that diminishes distally and a fine component that represents suspension fallout from a surface plume and resulting nepheloid layers. Graded brick-red beds in some cores were probably redeposited from turbidity currents. The lowermost bed on the Laurentian Fan and East Scotian Rise is immediately overlain by a carbonate-rich interval that can be identified all around the margin of the Grand Banks. This interval is correlated with detrital carbonate bed DC-1 in the Labrador Sea and Heinrich event H1 in the North Atlantic. The sequential occurrence of the two beds suggests that they may be a response to the same trigger, probably sea level rise, but that the Gulf of Saint Lawrence source was more easily destabilized.

(Table 1) Downcore age and bulk geochemistry of sediment core GS7202-35

The distribution of redox-sensitive metals in sediments is potentially a proxy for past ocean ventilation and productivity, but deconvolving these two major controls has proved difficult to date. Here we present a 740 kyr long record of trace element concentrations from an archived sediment core collected at ~15°S on the western flank of the East Pacific Rise (EPR) on 1.1 Myr old crust and underlying the largest known hydrothermal plume in the world ocean. The downcore trace element distribution is controlled by a variable diagenetic overprint of the inferred primary hydrothermal plume input. Two main diagenetic processes are operating at this site: redox cycling of transition metals and ferrihydrite to goethite transition during aging. The depth of oxidation in these sediments is controlled by fluctuations in the relative balance of bottom water oxygen and electron donor input (organic matter and hydrothermal sulfides). These fluctuations induce apparent variations in the accumulation of redox-sensitive species with time. Subsurface U and P peaks in glacial age sediments, in this and other published data sets along the southern EPR, indicate that basin-wide changes in deep ocean ventilation, in particular at glacial-interglacial terminations II, III, IV, and V, alter the depth of the oxidation front in the sediments. These basin-wide changes in the deep Pacific have significant implications for carbon partitioning in the ocean-atmosphere system, and the distribution of redox-sensitive metals in ridge crest sediment can be used to reconstruct past ocean conditions at abyssal depths in the absence of alternative proxy records.

Age determination of ODP Site 165-1002

An expanded Cariaco Basin 14C chronology is tied to 230Th-dated Hulu Cave speleothem records in order to provide detailed marine-based 14C calibration for the past 50,000 years. The revised, high-resolution Cariaco 14C calibration record agrees well with data from 230Th-dated fossil corals back to 33 ka, with continued agreement despite increased scatter back to 50 ka, suggesting that the record provides accurate calibration back to the limits of radiocarbon dating. The calibration data document highly elevated Delta14C during the Glacial period. Carbon cycle box model simulations show that the majority of observed Delta14C change can be explained by increased 14C production. However, from 45 to 15 ka, Delta14C remains anomalously high, indicating that the distribution of radiocarbon between surface and deep ocean reservoirs was different than it is today. Additional observations of the magnitude, spatial extent and timing of deep ocean Delta14C shifts are critical for a complete understanding of observed Glacial Delta14C variability.

Sedimentology, age models and stable isotope ratios of sediment cores from the equatorial Pacific

Based on detailed reconstructions of global distribution patterns, both paleoproductivity and the benthic d13C record of CO2, which is dissolved in the deep ocean, strongly differed between the Last Glacial Maximum and the Holocene. With the onset of Termination I about 15,000 years ago, the new (export) production of low- and mid-latitude upwelling cells started to decline by more than 2-4 Gt carbon/year. This reduction is regarded as a main factor leading to both the simultaneous rise in atmospheric CO2 as recorded in ice cores and, with a slight delay of more than 1000 years, to a large-scale gradual CO2 depletion of the deep ocean by about 650 Gt C. This estimate is based on an average increase in benthic d13C by 0.4-0.5 per mil. The decrease in new production also matches a clear 13C depletion of organic matter, possibly recording an end of extreme nutrient utilization in upwelling cells. As shown by Sarnthein et al., [1987], the productivity reversal appears to be triggered by a rapid reduction in the strength of meridional trades, which in turn was linked via a shrinking extent of sea ice to a massive increase in high-latitude insolation, i.e., to orbital forcing as primary cause.

Age model and palaeoclimate records over the past 22,000 years of sediment core GeoB10053-7

The Australian-Indonesian monsoon is an important component of the climate system in the tropical Indo-Pacific region. However, its past variability, relation with northern and southern high-latitude climate and connection to the other Asian monsoon systems are poorly understood. Here we present high-resolution records of monsoon-controlled austral winter upwelling during the past 22,000 years, based on planktic foraminiferal oxygen isotopes and faunal composition in a sedimentary archive collected offshore southern Java. We show that glacial-interglacial variations in the Australian-Indonesian winter monsoon were in phase with the Indian summer monsoon system, consistent with their modern linkage through cross-equatorial surface winds. Likewise, millennial-scale variability of upwelling shares similar sign and timing with upwelling variability in the Arabian Sea. On the basis of element composition and grain-size distribution as precipitation-sensitive proxies in the same archive, we infer that (austral) summer monsoon rainfall was highest during the Bølling-Allerød period and the past 2,500 years. Our results indicate drier conditions during Heinrich Stadial 1 due to a southward shift of summer rainfall and a relatively weak Hadley cell south of the Equator. We suggest that the Australian-Indonesian summer and winter monsoon variability were closely linked to summer insolation and abrupt climate changes in the northern hemisphere.

Grain size distribution and sea level changes in Gåsfjärden fjord

We analysed long-term variations in grain-size distribution in sediments from Gåsfjärden, a fjord-like inlet on the south-west Baltic Sea, and explored potential drivers of the recorded changes in sediment grain-size data. Over the last 5.4 thousand years (ka), the relative sea level decreased 17 m in the study region, caused by isostatic land uplift. As a consequence, Gåsfjärden has been transformed from an open coastal setting into a semi-closed inlet surrounded on the east by numerous small islands. To quantitatively estimate the morphological changes in Gåsfjärden over the last 5.4 ka and to further link the changes to our grain-size data, a digital elevation model (DEM)-based openness index was calculated. In the period between 5.4 and 4.4 ka BP, the inlet was characterised by the largest openness index. During this interval, the highest sand contents (~0.4 %) and silt/clay ratios (~0. 3) in the sediment sequence were recorded, indicating relatively high bottom water energy. After 4.4 ka BP, the average sand content was halved to ~0.2 % and the silt/clay ratios showed a significant decreasing trend over the last 4 ka. These changes are found to be associated with the gradual embayment of Gåsfjärden as represented in the openness index. The silt/clay ratios exhibited a delayed and slower change compared with the sand contents, which further suggest that finer particles are less sensitive to changes in hydrodynamic energy. Our DEM-based coastal openness index has proved to be a useful tool for interpreting the sedimentary grain-size record.

(Table 1) Distribution of upper Campanian-lower Maestrichtian calcareous nannofossils at DSDP Hole 71-511

Calcareous nannofossils from upper Campanian-lower Maestrichtian Deep Sea Drilling Project Leg 71 Cores 511-23 and 511-24 are described and correlated with assemblages of similar age from piston and drill cores on the Falkland Plateau, South Atlantic Ocean. The Leg 71 cores partially fill a drilling gap of at least 20 meters left within a thick (50 m) carbonate section first drilled by DSDP Leg 36 at Site 327. Cores 511-23 and 511-24 both fall within the upper portion of the Biscutum coronum Zone of Wind and demonstrate an overlap in the range of Monomarginatus quaternarius with the ranges of M. pectinatus, Misceomarginatus pleniporus, and Biscutum coronum across the Campanian/ Maestrichtian boundary. Resolution of the sequence of highest occurrence datums for the latter species must await the recovery of a more complete section. Comparison of the Site 511 assemblages with those from Mas Orcadas Core 07-75-44 to the north confirms earlier speculation that the Falkland Plateau served as an important boundary between major water masses during the Late Cretaceous.

Distribution of heavy minerals in sediment cores from the Laptev Sea continental margin and the central Arctic Ocean

Focussing on heavy-mineral associations in the Laptev-Sea continental margin area and the eastern Arctic Ocean, 129 surface sediment samples, two short and four long gravity cores have been studied. By means of the accessory components, heavy-mineral associations of surface sediment samples from the Laptev-See continental slope allowed the distinction into two different mineralogical provinces, each influenced by fluvial input of the Siberian river Systems. Transport pathways via sea ice from the shallow shelf areas into the Arctic Ocean up to the final ablation areas of the Fram Strait can be reconstructed by heavy-mineral data of surface sediments from the central Arctic Ocean. The shallow shelf of the Laptev Sea seems to be the most important source area for terrigenous material, as indicated by the abundant occurence of amphiboles and clinopyroxenes. Underneath the mixing Zone of the two dominating surface circulation Systems, the Beaufort- Gyre and Transpolar-Drift system, the imprint of the Amerasian shelf regions up to the Fram Strait is detectable because of a characteristical heavy-mineral association dominated by detrital carbonate and opaque minerals. Based On heavy-mineral characteristics of the potential circum-Arctic source areas, sea-ice drift, origin and distribution of ice-rafted material can be reconstructed during the past climatic cycles. Different factors controlling the transport of terrigenous material into the Arctic Ocean. The entrainment of particulate matter is triggered by the sea level, which flooded during highs and lows different regions resulting in the incorporation of sediment from different source areas into the sea ice. Additionally, the fluvial input even at low stands of sea level is responsible for the delivery of material of distinct sources for entrainment into the sea ice. Glacials and interglacials of climate cycles of the last 780 000 years left a characteristical signal in the central Arctic Ocean sediments caused by the ice- rafted material from different sources in the circum-Arctic regions and its change through time. Changes in the heavy-mineral association from an amphibole-dominated into a garnet-epidote-assemblage can be related to climate-related changes in source areas and directions of geostrophic winds, the dominating drive of the sea-ice drift. During Marine Isotope Stage (MIS) 6, the central Arctic Ocean is marked by an heavy-mineral signal, which occurs in recent sediments of the eastern Kara Sea. Its characteristics are high amounts of epidote, garnet and apatite. On the other hand, during the Same time interval a continuous record of Laptev Sea sediments is documented with high contents of amphiboles on the Lomonosov Ridge near the Laptev Sea continental slope. A nearly similar Pattern was detected in MIS 5 and 4. Small-scale glaciations in the Putorana-mountains and the Anabar-shield may have caused changes in the drainage area of the rivers and therefore a change in fluvial input. During MIS 3, the heavy-mineral association of central Arctic sediments show similar patterns than the Holocene mineral assemblage which consists of amphiboles, ortho- and clinopyroxenes with a Laptev Sea source. These minerals are indicating a stable Transpolar-Drift system similar to recent conditions. An extended influence of the Beaufort Gyre is only recognized, when sediment material from the Amerasian shelf areas reached the core location PS2757-718 during Termination Ib. Based On heavy-mineral data from Laptev-Sea continental slope Core PS2458-4 the paleo-sea-ice drift in the Laptev Sea during 14.000 years was reconstructed. During Holocene sea-level rise, the bathymetrically deeper parts of the Western shelf were flooded first. At the beginning of the Atlantic stage, nearly the entire shelf was marine influenced by fully marine conditions and the recent surface circulation was established.

Planktic foraminiferal distribution and stable isotope ratios of sediment core MSM05/5_712-1 from the Arctic Ocean

The Arctic is responding more rapidly to global warming than most other areas on our planet. Northward flowing Atlantic Water is the major means of heat advection towards the Arctic and 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 ~150 years. Here, we present a multidecadal scale record of ocean temperature variations during the last 2000 years, derived from marine sediments off Western Svalbard (79°N). We find that early-21st-century temperatures of Atlantic Water entering the Arctic Ocean are unprecedented over the past 2000 years and are presumably linked to the Arctic Amplification of global warming.