(Tables S1,2) Stable isotopes in muscle tissue of polar bears (Ursus maritimus) from various regions of the Arctic

The relative contribution of regional contamination versus dietary differences to geographic variation in polar bear (Ursus maritimus) contaminant levels is unknown. Dietary variation between Alaska, Canada, East Greenland, and Svalbard subpopulations was assessed by muscle nitrogen and carbon stable isotope (d15N, d13C) and adipose fatty acid (FA) signatures relative to their main prey (ringed seals). Western and southern Hudson Bay signatures were characterized by depleted d15N and d13C, lower proportions of C20 and C22 monounsaturated FAs and higher proportions of C18 and longer chain polyunsaturated FAs. East Greenland and Svalbard signatures were reversed relative to Hudson Bay. Alaskan and Canadian Arctic signatures were intermediate. Between-subpopulation dietary differences predominated over interannual, seasonal, sex, or age variation. Among various brominated and chlorinated contaminants, diet signatures significantly explained variation in adipose levels of polybrominated diphenyl ether (PBDE) flame retardants (14-15%) and legacy PCBs (18-21%). However, dietary influence was contaminant class-specific, since only low or nonsignificant proportions of variation in organochlorine pesticide (e.g., chlordane) levels were explained by diet. Hudson Bay diet signatures were associated with lower PCB and PBDE levels, whereas East Greenland and Svalbard signatures were associated with higher levels. Understanding diet/food web factors is important to accurately interpret contaminant trends, particularly in a changing Arctic.

(Supplement T1) Clay mineral compositions at the Aleutian Basin of the Bering Sea from IODP sites 323-U1343

The clay mineral composition at IODP Exp. 323 Site U1343 in the Bering Sea was analyzed so as to unravel their provenance over glacial-interglacial cycles for the last 2.4 Ma. Smectite was negatively correlated with the sum of illite and chlorite; therefore, their ratio [S/(I + C)] was used as an indicator of clay mineral composition changes. In general, the S/(I + C) ratio was rather similar for glacial and interglacial periods during most of the last 2.4 Ma. In addition, these results overlap with those of surface sediments in the modern East Aleutian Basin, which suggests that smectite-rich clay particles are delivered from the Aleutians by the northward Bering Slope Current (BSC) rather than from rivers in western Alaska. However, some clay mineral compositions of the glacial periods after the Mid-Pleistocene Transition (MPT: 1.25-0.7 Ma) were characterized by low illite and relatively high smectite. During this period, extensive ice-rafting might have transported the smectite-rich clay particles to Site U1343 from the glacial shelf off Alaska and/or from glacial rivers from that area.

Composition of bottom sediments from Core SO201-2-101, Shirshov Ridge, Bering Sea

Analysis of lithology, grain-size composition, clay minerals, and geochemistry of Upper Pleistocene bottom sediments from the submarine Shirshov Ridge (Bering Sea) showed that the Yukon-Tanana terrane of the Central Alaska was main source area of the sediments. Sedimentary material was transported by the Yukon River through Beringia up to the shelf break, where they were entrained by a strong north-west sea current. Lithological data revealed several pulses of ice-rafted debris deposition roughly synchronous with Heinrich events and periods of weaker bottom current intensity. Based on geochemical results we distinguished intervals of an increase in paleoproductivity and extension of the oxygen minimum zone. Our results suggest that there were three stages of deposition driven by glacioeustatic sea-level fluctuations and glacial cycles in Alaska.

Seasonality in pH and temperature for shallow zones of the Bering Sea from measurements of coral Clathromorphum nereostratum

Increasing atmospheric CO2 concentrations are potentially affecting marine ecosystems twofold, by warming and acidification. The rising amount of CO2 taken up by the ocean lowers the saturation state of calcium carbonate, complicating the formation of this key biomineral used by many marine organisms to build hard parts like skeletons or shells. Reliable time-series data of seawater pH are needed to evaluate the ongoing change and compare long-term trends and natural variability. For the high-latitude ocean, the region facing the strongest CO2 uptake, such time-series data are so far entirely lacking. Our study provides, to our knowledge, the first reconstruction of seasonal cycle and long-term trend in pH for a high-latitude ocean obtained from 2D images of stable boron isotopes from a coralline alga.

(Supplementary Table 1) Mercury isotope values for common and thick-billed murre (Uria aalge, Uria lomvia) eggs from Alaska

Atmospheric deposition of mercury to remote areas has increased threefold since pre-industrial times. Mercury deposition is particularly pronounced in the Arctic. Following deposition to surface oceans and sea ice, mercury can be converted into methylmercury, a biologically accessible form of the toxin, which biomagnifies along the marine food chain. Mass-independent fractionation of mercury isotopes accompanies the photochemical breakdown of methylmercury to less bioavailable forms in surface waters. Here we examine the isotopic composition of mercury in seabird eggs collected from colonies in the North Pacific Ocean, the Bering Sea and the western Arctic Ocean, to determine geographical variations in methylmercury breakdown at northern latitudes. We find evidence for mass-independent fractionation of mercury isotopes. The degree of mass-independent fractionation declines with latitude. Foraging behaviour and geographic variations in mercury sources and solar radiation fluxes were unable to explain the latitudinal gradient. However, mass-independent fractionation was negatively correlated with sea-ice cover. We conclude that sea-ice cover impedes the photochemical breakdown of methylmercury in surface waters, and suggest that further loss of Arctic sea ice this century will accelerate sunlight-induced breakdown of methylmercury in northern surface waters.

Laminated sediment cores in the Bering Sea

During the last glacial termination, the upper North Pacific Ocean underwent dramatic and rapid changes in oxygenation that lead to the transient intensification of oxygen minimum zones (OMZs), recorded by the widespread occurrence of laminated sediments on circum-Pacific continental margins. We present a new laminated sediment record from the mid-depth (1100 m) northern Bering Sea margin that provides insight into these deglacial OMZ maxima with exceptional, decadal-scale detail. Combined ultrahigh-resolution micro-X-ray-fluorescence (micro-XRF) data and sediment facies analysis of laminae reveal an alternation between predominantly terrigenous and diatom-dominated opal sedimentation. The diatomaceous laminae are interpreted to represent spring/summer productivity events related to the retreating sea ice margin.We identified five laminated sections in the deglacial part of our site. Lamina counts were carried out on these sections and correlated with the Bølling-Allerød and Preboreal phases in the North Greenland Ice Core (NGRIP) oxygen isotope record, indicating an annual deposition of individual lamina couplets (varves). The observed rapid decadal intensifications of anoxia, in particular within the Bølling-Allerød, are tightly coupled to short-term warm events through increases in regional export production. This dependence of laminae formation on warmer temperatures is underlined by a correlation with published Bering Sea sea surface temperature records and d18O data of planktic foraminifera from the Gulf of Alaska. The rapidity of the observed changes strongly implies a close atmospheric teleconnection between North Pacific and North Atlantic regions.We suggest that concomitant increases in export production and subsequent remineralization of organic matter in the Bering Sea, in combination with oxygen-poor waters entering the Being Sea, drove down oxygen concentrations to values below 0.1ml/l and caused laminae preservation. Calculated benthic-planktic ventilation ages show no significant variations throughout the last deglaciation, indicating that changes in formation rates or differing sources of North Pacific mid-depth waters are not prime candidates for strengthening the OMZ at our site. The age models established by our correlation procedure allow for the determination of calendar age control points for the Bølling-Allerød and the Preboreal that are independent of the initial radiocarbon-based chronology. Resulting surface reservoir ages range within 730-990 yr during the Bølling-Allerød, 800-1100 yr in the Younger Dryas, and 765-775 yr for the Preboreal.

Lithic grain and mineral composition of sand and sandstones from the North Pacific Ocean and the Bering Sea (Table 3)

Sand and sandstone compositions from different types of basins reflect provenance terranes governed by plate tectonics. One hundred and one thin sections of Upper Miocene to Holocene sand-sized material were examined from DSDP/IPOD Sites in the North Pacific Ocean and the Bering Sea. The Gazzi-Dickinson point-counting method was used to establish compositional characteristics of sands from different tectonic settings. Continental margin forearc sands from the western North America continental margin arc system are clearly different from backarc/marginal-sea sands from the Aleutian intraoceanic arc system. The forearc sands have average QFL percentages of 29-42-29, LmLvLst percentages of 32-34-34, 3 Fmwk%M and 0.82 P/F. Aleutian backarc sands have average QFL percentages of 8-22-69. LmLvLst percentages of 9-85-6, 0.5 Fmwk%M and 0.96 P/F. A trend of increasing QFL%Q and decreasing LmLvLst%Lv westward in the backarc region of the Aleutian Ridge reflects the influence of the Asiatic continental margin. Aleutian backarc sands without continental influence have average QFL percentages of 1-20-79, LmLvLst percentages of 1-98-1, 0 Fmwk%M and 0.99 P/F. Of the continental margin forearc samples, sands on the Astoria Fan (west of the Oregon-Washington trench) contain the highest LmLvLst%Lv and lowest P/F; sands from mixed transform-fault and trench settings (Delgada Fan and Gulf of Alaska samples) have slightly higher Qp/Q (0.03); and sands from the Pacific-Juan de Fuca-North America triple junction have the highest Fmwk%M. Delgada Fan and Gulf of Alaska sands have average QFL percentages of 27-38-35, LmLvLst percentages of 37-26-37, 2 Fmwk%M and 0.86 P/F. Astoria Fan sands have average QFL percentages of 35-41-24, LmLvLst percentages of 30-47-23, 3 Fmwk%M and 0.74 P/F. The triple-junction sands have average QFL percentages of 28-59-13, LmLvLst percentages of 25-26-49, 9 Fmwk%M and 0.87 P/F. The petrologic data from the modern ocean basins examined in this study can provide useful analogs for interpretation of ancient oceanic sequences. Our data suggest some refinements of, but generally substantiate, existing petrologic models relating sandstone composition to tectonic setting.

New Pliocene and Pleistocene acritarchs from the North Atlantic, Arctic Ocean and Bering Sea

Acritarchs have received limited attention in palynological studies of the Cenozoic, although they have much potential both for refining Neogene and Quaternary stratigraphy, especially in mid- and high northern latitudes, and developing palaeoceanographical reconstructions. Here we formally describe and document the stratigraphical and palaeotemperature ranges (from foraminiferal Mg/Ca) of four new acritarch species: Cymatiosphaera? aegirii sp. nov., Cymatiosphaera? fensomei sp. nov., Cymatiosphaera? icenorum sp. nov. and Lavradosphaera canalis sp. nov. In reviewing the stratigraphical distributions of all species of the genus Lavradosphaera De Schepper & Head, 2008, we demonstrate their correlation potential between the North Atlantic and Bering Sea in the Pliocene. Additionally, Lavradosphaera lucifer De Schepper & Head, 2008 and Lavradosphaera canalis sp. nov., while not themselves overlapping stratigraphically, have morphological intermediates that do partially overlap and may represent an evolutionary trend consequent upon climate cooling in the Late Pliocene. Finally, we show that the highest abundances of the acritarchs presented here were living in the eastern North Atlantic, in surface-water temperatures not very different from today.

Age models, color measurements, and calculated sea surface temperatures of sediment cores from the Bering Sea and the Sea of Okhotsk

Past changes in North Pacific sea surface temperatures and sea-ice conditions are proposed to play a crucial role in deglacial climate development and ocean circulation but are less well known than from the North Atlantic. Here, we present new alkenone-based sea surface temperature records from the subarctic northwest Pacific and its marginal seas (Bering Sea and Sea of Okhotsk) for the time interval of the last 15 kyr, indicating millennial-scale sea surface temperature fluctuations similar to short-term deglacial climate oscillations known from Greenland ice-core records. Past changes in sea-ice distribution are derived from relative percentage of specific diatom groups and qualitative assessment of the IP25 biomarker related to sea-ice diatoms. The deglacial variability in sea-ice extent matches the sea surface temperature fluctuations. These fluctuations suggest a linkage to deglacial variations in Atlantic meridional overturning circulation and a close atmospheric coupling between the North Pacific and North Atlantic. During the Holocene the subarctic North Pacific is marked by complex sea surface temperature trends, which do not support the hypothesis of a Holocene seesaw in temperature development between the North Atlantic and the North Pacific.