(Table S1) Diatom species richness from the Early Pleistocene to present

The extent to which the spatial distribution of marine planktonic microbes is controlled by local environmental selection or dispersal is poorly understood. Our ability to separate the effects of these two biogeographic controls is limited by the enormous environmental variability both in space and through time. To circumvent this limitation, we analyzed fossil diatom assemblages over the past ~1.5 million years from the world oceans and show that these eukaryotic microbes are not limited by dispersal. The lack of dispersal limitation in marine diatoms suggests that the biodiversity at the microbial level fundamentally differs from that of macroscopic animals and plants for which geographic isolation is a common component of speciation.

Model-based changes in global annual mean surface temperature change (Delta T_g) and radiative forcing due to land ice albedo changes (Delta R_[LI]) over the last 5 Myr, supplementary material

It is still an open question how equilibrium warming in response to increasing radiative forcing - the specific equilibrium climate sensitivity S - depends on background climate. We here present palaeodata-based evidence on the state dependency of S, by using CO2 proxy data together with a 3-D ice-sheet-model-based reconstruction of land ice albedo over the last 5 million years (Myr). We find that the land ice albedo forcing depends non-linearly on the background climate, while any non-linearity of CO2 radiative forcing depends on the CO2 data set used. This non-linearity has not, so far, been accounted for in similar approaches due to previously more simplistic approximations, in which land ice albedo radiative forcing was a linear function of sea level change. The latitudinal dependency of ice-sheet area changes is important for the non-linearity between land ice albedo and sea level. In our set-up, in which the radiative forcing of CO2 and of the land ice albedo (LI) is combined, we find a state dependence in the calculated specific equilibrium climate sensitivity, S[CO2,LI], for most of the Pleistocene (last 2.1 Myr). During Pleistocene intermediate glaciated climates and interglacial periods, S[CO2,LI] is on average ~ 45 % larger than during Pleistocene full glacial conditions. In the Pliocene part of our analysis (2.6-5 Myr BP) the CO2 data uncertainties prevent a well-supported calculation for S[CO2,LI], but our analysis suggests that during times without a large land ice area in the Northern Hemisphere (e.g. before 2.82 Myr BP), the specific equilibrium climate sensitivity, S[CO2,LI], was smaller than during interglacials of the Pleistocene. We thus find support for a previously proposed state change in the climate system with the widespread appearance of northern hemispheric ice sheets. This study points for the first time to a so far overlooked non-linearity in the land ice albedo radiative forcing, which is important for similar palaeodata-based approaches to calculate climate sensitivity. However, the implications of this study for a suggested warming under CO2 doubling are not yet entirely clear since the details of necessary corrections for other slow feedbacks are not fully known and the uncertainties that exist in the ice-sheet simulations and global temperature reconstructions are large.

Susceptibility and time-dependent magnetization of basalts at DSDP Hole 76-534A

Basalts from Hole 534A are among the oldest recovered from the ocean bottom, dating from the opening of the Atlantic 155 Ma. Upon exposure to a 1-Oe field for one week, these basalts acquire a viscous remanent magnetization (VRM), which ranges from 4 to 223% of their natural remanent magnetization (NRM). A magnetic field of similar magnitude is observed in the paleomagnetic lab of the Glomar Challenger, and it is therefore doubtful if accurate measurements of magnetic moment in such rocks can be made on board unless the paleomagnetic area is magnetically shielded. No correlation is observed between the Konigsberger ratio (beta), which is usually less than 3, and the ability to acquire a VRM. The VRM shows both a log t dependence and a Richter aftereffect. Both of these, but especially the log t dependence, will cause the susceptibility measurements (made by applying a magnetic field for a very short time) to be minimum values. The susceptibility and derived Q should therefore be used cautiously for magnetic anomaly interpretation, because they can cause the importance of the induced magnetization to be underestimated.

Modeling the sensitivity of precipitation oxygen isotopes to the Last Glacial Maximum boundary conditions: A study using Community Atmosphere Model (CAM3.0)

To understand the validity of d18O proxy records as indicators of past temperature change, a series of experiments was conducted using an atmospheric general circulation model fitted with water isotope tracers (Community Atmosphere Model version 3.0, IsoCAM). A pre-industrial simulation was performed as the control experiment, as well as a simulation with all the boundary conditions set to Last Glacial Maximum (LGM) values. Results from the pre-industrial and LGM simulations were compared to experiments in which the influence of individual boundary conditions (greenhouse gases, ice sheet albedo and topography, sea surface temperature (SST), and orbital parameters) were changed each at a time to assess their individual impact. The experiments were designed in order to analyze the spatial variations of the oxygen isotopic composition of precipitation (d18Oprecip) in response to individual climate factors. The change in topography (due to the change in land ice cover) played a significant role in reducing the surface temperature and d18Oprecip over North America. Exposed shelf areas and the ice sheet albedo reduced the Northern Hemisphere surface temperature and d18Oprecip further. A global mean cooling of 4.1 °C was simulated with combined LGM boundary conditions compared to the control simulation, which was in agreement with previous experiments using the fully coupled Community Climate System Model (CCSM3). Large reductions in d18Oprecip over the LGM ice sheets were strongly linked to the temperature decrease over them. The SST and ice sheet topography changes were responsible for most of the changes in the climate and hence the d18Oprecip distribution among the simulations.

Chemistry of magmatic samples of ODP Holes 134-832B and 134-833B

New petrographic and compositional data were reported for 143 samples of core recovered from Sites 832 and 833 during Ocean Drilling Program (ODP) Leg 134. Site 832 is located in the center and Site 833 is on the eastern edge of the North Aoba Basin, in the central part of the New Hebrides Island Arc. This basin is bounded on the east (Espiritu Santo and Malakula islands) and west (Pentecost and Maewo islands) by uplifted volcano-sedimentary ridges associated with collision of the d'Entrecasteaux Zone west of the arc. The currently active Central Belt volcanic front extends through the center of this basin and includes the shield volcanoes of Aoba, Ambrym, and Santa Maria islands. The oldest rocks recovered by drilling are the lithostratigraphic Unit VII Middle Miocene volcanic breccias in Hole 832B. Lava clasts are basaltic to andesitic, and the dominant phenocryst assemblage is plagioclase + augite + orthopyroxene + olivine. These clasts characteristically contain orthopyroxene, and show a low to medium K calc-alkaline differentiation trend. They are tentatively correlated with poorly documented Miocene calc-alkaline lavas and intrusives on adjacent Espiritu Santo Island, although this correlation demands that the measured K-Ar of 5.66 Ma for one clast is too young, due to alteration and Ar loss. Lava clasts in the Hole 832B Pliocene-Pleistocene sequence are mainly ankaramite or augite-rich basalt and basaltic andesite; two of the most evolved andesites have hornblende phenocrysts. These lavas vary from medium- to high-K compositions and are derived from a spectrum of parental magmas for which their LILE and HFSE contents show a broad inverse correlation with SiO2 contents. We hypothesize that this spectrum results from partial melting of an essentially similar mantle source, with the low-SiO2 high HFSE melts derived by lower degrees of partial melting probably at higher pressures than the high SiO2, low HFSE magmas. This same spectrum of compositions occurs on the adjacent Central Chain volcanoes of Aoba and Santa Maria, although the relatively high-HFSE series is known only from Aoba. Late Pliocene to Pleistocene lava breccias in Hole 833B contain volcanic clasts including ankaramite and augite + olivine + plagioclase-phyric basalt and rare hornblende andesite. These clasts are low-K compositions with flat REE patterns and have geochemical affinities quite different from those recovered from the central part of the basin (Hole 832B). Compositionally very similar lavas occur on Merelava volcano, 80 km north of Site 833, which sits on the edge of the juvenile Northern (Jean Charcot) Trough backarc basin that has been rifting the northern part of the New Hebrides Island Arc since 2-3 Ma. The basal sedimentary rocks in Hole 833B are intruded by a series of Middle Pliocene plagioclase + augite +/- olivine-phyric sills with characteristically high-K evolved basalt to andesite compositions, transitional to shoshonite. These are compositionally correlated with, though ~3 m.y. older than, the high-HFSE series described from Aoba. The calc-alkaline clasts in Unit VII of Hole 832B, correlated with similar lavas of Espiritu Santo Island further west, presumably were erupted before subduction polarity reversal perhaps 6-10 Ma. All other samples are younger than subduction reversal and were generated above the currently subduction slab. The preponderance in the North Aoba Basin and adjacent Central Chain islands of relatively high-K basaltic samples, some with transitional alkaline compositions, may reflect a response to collision of the d'Entrecasteaux Zone with the arc some 2-4 Ma. This may have modified the thermal structure of the subduction zone, driving magma generation processes to deeper levels than are present normally along the reminder of the New Hebrides Island Arc.

Effect of simulated ocean acidification on the acute toxicity of Cu and Cd to Tigriopus japonicus

Heavy metals pollution in marine environments has caused great damage to marine biological and ecological systems. Heavy metals accumulate in marine creatures, after which they are delivered to higher trophic levels of marine organisms through the marine food chain, which causes serious harm to marine biological systems and human health. Additionally, excess carbon dioxide in the atmosphere has caused ocean acidification. Indeed, about one third of the CO2 released into the atmosphere by anthropogenic activities since the beginning of the industrial revolution has been absorbed by the world's oceans, which play a key role in moderating climate change. Modeling has shown that, if current trends in CO2 emissions continue, the average pH of the ocean will reach 7.8 by the end of this century, corresponding to 0.5 units below the pre-industrial level, or a three-fold increase in H+ concentration. The ocean pH has not been at this level for several millions of years. Additionally, these changes are occurring at speeds 100 times greater than ever previously observed. As a result, several marine species, communities and ecosystems might not have time to acclimate or adapt to these fast changes in ocean chemistry. In addition, decreasing ocean pH has the potential to seriously affect the growth, development and reproduction reproductive processes of marine organisms, as well as threaten normal development of the marine ecosystem. Copepods are an important part of the meiofauna that play an important role in the marine ecosystem. Pollution of the marine environment can influence their growth and development, as well as the ecological processes they are involved in. Accordingly, there is important scientific value to investigation of the response of copepods to ocean acidification and heavy metals pollution. In the present study, we evaluated the effects of simulated future ocean acidification and the toxicological interaction between ocean acidity and heavy metals of Cu and Cd on T. japonicus. To accomplish this, harpacticoids were exposed to Cu and Cd concentration gradient seawater that had been equilibrated with CO2 and air to reach pH 8.0, 7.7, 7.3 and 6.5 for 96 h. Survival was not significantly suppressed under single sea water acidification, and the final survival rates were greater than 93% in both the experimental groups and the controls. The toxicity of Cu to T. japonicus was significantly affected by sea water acidification, with the 96h LC50 decreasing by nearly threefold from 1.98 to 0.64 mg/L with decreasing pH. The 96 h LC50 of Cd decreased with decreasing pH, but there was no significant difference in mortality among pH treatments. The results of the present study demonstrated that the predicted future ocean acidification has the potential to negatively affect survival of T. japonicus by exacerbating the toxicity of Cu. The calculated safe concentrations of Cu were 11.9 (pH 7.7) and 10.5 (pH 7.3) µg/L, which were below the class I value and very close to the class II level of the China National Quality Standard for Sea Water. Overall, these results indicate that the Chinese coastal sea will face a

(Table 1) Physical oceanography from several POLARSTERN cruises to the Arctic Ocean, 1993 - 2012

Freshwater in the Arctic Ocean plays an important role in the regional ocean circulation, sea ice, and global climate. From salinity observed by a variety of platforms, we are able, for the first time, to estimate a statistically reliable liquid freshwater trend from monthly gridded fields over all upper Arctic Ocean basins. From 1992 to 2012 this trend was 600±300 km**3/yr. A numerical model agrees very well with the observed freshwater changes. A decrease in salinity made up about two thirds of the freshwater trend and a thickening of the upper layer up to one third. The Arctic Ocean Oscillation index, a measure for the regional wind stress curl, correlated well with our freshwater time series. No clear relation to Arctic Oscillation or Arctic Dipole indices could be found. Following other observational studies, an increased Bering Strait freshwater import to the Arctic Ocean, a decreased Davis Strait export, and enhanced net sea ice melt could have played an important role in the freshwater trend we observed.