Compilation of Mg/Ca ratios of Neogloboquadrina pachyderma (sinistral) and SST estimates from the northern North Atlantic

A map of estimated calcification temperatures of the planktic foraminifer Neogloboquadrina pachyderma sinistral (TNps) for the Nordic Seas and the northern North Atlantic for the Last Glacial Maximum was produced from oxygen isotopes with support of Mg/Ca ratios. To arrive at the reconstruction, several constraints concerning the plausible salinity and ?18O-fields were employed. The reconstruction indicates inflow of temperate waters in a wedge along the eastern border of the Nordic Seas and at least seasonally ice-free waters. The reconstruction from oxygen isotopes shows similarities with Mg/Ca based paleotemperatures in the southern and southeastern sector, while unrealistically high Mg/Ca values in the central Nordic Seas prevent the application of the method in this area. The oxygen isotope based reconstruction shows some agreement with temperature reconstructions based on the modern analogue technique, but with somewhat lower temperatures and a stronger internal gradient inside the Nordic Seas. All told, our results suggest a much more ice-free and dynamic high latitude ocean than the CLIMAP reconstruction.

Theoretical model simulations for variations in thermal reaction norms

Thermal reaction norms for growth rates of six Emiliania huxleyi isolates originating from the central Atlantic (Azores, Portugal) and five isolates from the coastal North Atlantic (Bergen, Norway) were assessed. We used the template mode of variation model to decompose variations in growth rates into modes of biological interest: vertical shift, horizontal shift, and generalist-specialist variation. In line with the actual habitat conditions, isolates from Bergen (Bergen population) grew well at lower temperatures, and isolates from the Azores (Azores population) performed better at higher temperatures. The optimum growth temperature of the Azores population was significantly higher than that of the Bergen population. Neutral genetic differentiation was found between populations by microsatellite analysis. These findings indicate that E. huxleyi populations are adapted to local temperature regimes. Next to between-population variation, we also found variation within populations. Genotype-by-environment interactions resulted in the most pronounced phenotypic differences when isolates were exposed to temperatures outside the range they naturally encounter. Variation in thermal reaction norms between and within populations emphasizes the importance of using more than one isolate when studying the consequences of global change on marine phytoplankton. Phenotypic plasticity and standing genetic variation will be important in determining the potential of natural E. huxleyi populations to cope with global climate change.

Rare earth element concentrations and geochemistry of carbonate minerals from the Logatchev Hydrothermal Field and the Gakkel Ridge

The ultramafic-hosted Logatchev Hydrothermal Field (LHF) at 15°N on the Mid-Atlantic Ridge and the Arctic Gakkel Ridge (GR) feature carbonate precipitates (aragonite, calcite, and dolomite) in voids and fractures within different types of host rocks. We present chemical and Sr isotopic compositions of these different carbonates to examine the conditions that led to their formation. Our data reveal that different processes have led to the precipitation of carbonates in the various settings. Seawater-like 87Sr/86Sr ratios for aragonite in serpentinites (0.70909 to 0.70917) from the LHF are similar to those of aragonite from the GR (0.70912 to 0.70917) and indicate aragonite precipitation from seawater at ambient conditions at both sites. Aragonite veins in sulfide breccias from LHF also have seawater-like Sr isotope compositions (0.70909 to 0.70915), however, their rare earth element (REE) patterns show a clear positive europium (Eu) anomaly indicative of a small (< 1%) hydrothermal contribution. In contrast to aragonite, dolomite from the LHF has precipitated at much higher temperatures (~100 °C), and yet its 87Sr/86Sr ratios (0.70896 to 0.70907) are only slightly lower than those of aragonite. Even higher temperatures are calculated for the precipitation of deformed calcite veins in serpentine-talc fault schists form north of the LHF. These calcites show unradiogenic 87Sr/86Sr ratios (0.70460 to 0.70499) indicative of precipitation from evolved hydrothermal fluids. A simple mixing model based on Sr mass balance and enthalpy conservation indicates strongly variable conditions of fluid mixing and heat transfers involved in carbonate formation. Dolomite precipitated from a mixture of 97% seawater and 3% hydrothermal fluid that should have had a temperature of approximately 14 °C assuming that no heat was transferred. The much higher apparent precipitation temperatures based on oxygen isotopes (~ 100 °C) may be indicative of conductive heating, probably of seawater prior to mixing. The hydrothermal calcite in the fault schist has precipitated from a mixture of 67% hydrothermal fluid and 33% seawater, which should have had an isenthalpic mixing temperature of ~ 250 °C. The significantly lower temperatures calculated from oxygen isotopes are likely due to conductive cooling of hydrothermal fluid discharging along faults. Rare earth element patterns corroborate the results of the mixing model, since the hydrothermal calcite, which formed from waters with the greatest hydrothermal contribution, has REE patterns that closely resemble those of vent fluids from the LHF. Our results demonstrate, for the first time, that (1) precipitation from pure seawater, (2) conductive heating of seawater, and (3) conductive cooling of hydrothermal fluids in the sub-seafloor all can lead to carbonate precipitation within a single ultramafic-hosted hydrothermal system.

Seawater carbonate chemistry and biological processes of Mytilus edulis during experiments, 2011

Ocean acidification (OA) resulting from anthropogenic emissions of carbon dioxide (CO2) has already lowered and is predicted to further lower surface ocean pH. There is a particular need to study effects of OA on organisms living in cold-water environments due to the higher solubility of CO2 at lower temperatures. Mussel larvae (Mytilus edulis) and shrimp larvae (Pandalus borealis) were kept under an ocean acidification scenario predicted for the year 2100 (pH 7.6) and compared against identical batches of organisms held under the current oceanic pH of 8.1, which acted as a control. The temperature was held at a constant 10°C in the mussel experiment and at 5°C in the shrimp experiment. There was no marked effect on fertilization success, development time, or abnormality to the D-shell stage, or on feeding of mussel larvae in the low-pH (pH 7.6) treatment. Mytilus edulis larvae were still able to develop a shell in seawater undersaturated with respect to aragonite (a mineral form of CaCO3), but the size of low-pH larvae was significantly smaller than in the control. After 2 mo of exposure the mussels were 28% smaller in the pH 7.6 treatment than in the control. The experiment with Pandalus borealis larvae ran from 1 through 35 days post hatch. Survival of shrimp larvae was not reduced after 5 wk of exposure to pH 7.6, but a significant delay in zoeal progression (development time) was observed.

Multi-proxy record of planktonic and benthic foraminiferal d18O, and Mg/Ca-derived surface and lower thermocline temperatures in sediment core SO18471 (GIK18471), Timor Sea

Paleostudies of the Indonesian Throughflow (ITF) are largely based on temperature and salinity reconstructions of its near surface component, whereas the variability of its lower thermocline flow has rarely been investigated. We present a multi-proxy record of planktonic and benthic foraminiferal d18O, Mg/Ca-derived surface and lower thermocline temperatures, X-ray fluorescence (XRF)-derived runoff and sediment winnowing for the past 130 ka in marine sediment core SO18471. Core SO18471, retrieved from a water depth of 485 m at the southern edge of the Timor Strait close to the Sahul Shelf, sits in a strategic position to reconstruct variations in both the ITF surface and lower thermocline flow as well as to investigate hydrological changes related to monsoon variability and shelf dynamics over time. Sediment winnowing demonstrates that the ITF thermocline flow intensified during MIS 5d-a and MIS 1. In contrast during MIS 5e, winnowing was reduced and terrigenous input increased suggesting intensification of the local wet monsoon and a weaker ITF. Lower thermocline warming during globally cold periods (MIS 4 - MIS 2) appears to be related to a weaker and contracted thermocline ITF and advection of warm and salty Indian Ocean waters.

(Table 1) Oxygen isotope temperatures from hydrothermally altered breccias at DSDP Hole 60-453

Gabbro-metabasalt polymict breccias cored in Deep Sea Drilling Project Hole 453 are cemented in part by hydrothermal alteration to lower greenschist facies (chlorite-epidote-actinolite) mineral assemblages. Temperature estimates for this alteration, based on oxygen isotope determinations of secondary minerals, are nearly 100°C at the top of the breccias and over 200°C in a zone of intense alteration near the base.

Cenomanian/Turonian sea surface temperatures of the equatorial Atlantic reconstructed from geochemical data

Oceanic anoxic event 2 (OAE-2) occurring during the Cenomanian/Turonian (C/T) transition is evident from a globally recognized positive stable carbon isotopic excursion and is thought to represent one of the most extreme carbon cycle perturbations of the last 100 Myr. However, the impact of this major perturbation on and interaction with global climate remains unclear. Here we report new high-resolution records of sea surface temperature (SST) based on TEX86 and d 18O of excellently preserved planktic foraminifera and stable organic carbon isotopes across the C/T transition from black shales located offshore Suriname/French Guiana (Demerara Rise, Ocean Drilling Program Leg 207 Site 1260) and offshore Senegal (Cape Verde Basin, Deep Sea Drilling Project Leg 41 Site 367). At Site 1260, where both SST proxy records can be determined, a good match between conservative SST estimates from TEX86 and d 18O is observed. We find that late Cenomanian SSTs in the equatorial Atlantic Ocean (33°C) were substantially warmer than today (27°-29°C) and that the onset of OAE-2 coincided with a rapid shift to an even warmer (35°-36°C) regime. Within the early stages of the OAE a marked (4°C) cooling to temperatures lower than pre-OAE conditions is observed. However, well before the termination of OAE-2 the warm regime was reestablished and persisted into the Turonian. Our findings corroborate the view that the C/T transition represents the onset of the interval of peak Cretaceous warmth. More importantly, they are consistent with the hypotheses that mid-Cretaceous warmth can be attributed to high levels of atmospheric carbon dioxide (CO2) and that major OAEs were capable of triggering global cooling through the negative feedback effect of organic carbon-burial-led CO2 sequestration. Evidently, however, the factors that gave rise to the observed shift to a warmer climate regime at the onset of OAE-2 were sufficiently powerful that they were only briefly counterbalanced by the high rates of carbon burial attained during even the most extreme interval of organic carbon burial in the last 100 Myr.

Sea surface temperatures during Marine Isotope Stage (MIS) 11

The Marine Isotope Stage (MIS) 11 (424-374 ka) was characterized by a protracted deglaciation and an unusually long climatic optimum. It remains unclear to what degree the climate development during this interglacial reflects the unusually weak orbital forcing or greenhouse gas trends. Previously, arguments about the duration and timing of the MIS11 climatic optimum and about the pace of the deglacial warming were based on a small number of key records, which appear to show regional differences. In order to obtain a global signal of climate evolution during MIS11, we compiled a database of 78 sea surface temperature (SST) records from 57 sites spanning MIS11, aligned these individually on the basis of benthic (N = 28) or planktonic (N = 31) stable oxygen isotope curves to a common time frame and subjected 48 of them to an empirical orthogonal function (EOF) analysis. The analysis revealed a high commonality among all records, with the principal SST trend explaining almost 49% of the variability. This trend indicates that on the global scale, the surface ocean underwent rapid deglacial warming during Termination V, in pace with carbon dioxide rise, followed by a broad SST optimum centered at ~410 kyr. The second EOF, which explained ~18% of the variability, revealed the existence of a different SST trend, characterized by a delayed onset of the temperature optimum during MIS11 at ~398 kyr, followed by a prolonged warm period lasting beyond 380 kyr. This trend is most consistently manifested in the mid-latitude North Atlantic and Mediterranean Sea and is here attributed to the strength of the Atlantic meridional overturning circulation. A sensitivity analysis indicates that these results are robust to record selection and to age-model uncertainties of up to 3-6 kyr, but more sensitive to SST seasonal attribution and SST uncertainties >1 °C. In order to validate the CCSM3 (Community Climate System Model, version 3) predictive potential, the annual and seasonal SST anomalies recorded in a total of 74 proxy records were compared with runs for three time slices representing orbital configuration extremes during the peak interglacial of MIS11. The modeled SST anomalies are characterized by a significantly lower variance compared to the reconstructions. Nevertheless, significant correlations between proxy and model data are found in comparisons on the seasonal basis, indicating that the model captures part of the long-term variability induced by astronomical forcing, which appears to have left a detectable signature in SST trends.

Mean air temperatures of Greenland

Melt rate and surface temperature on the Greenland ice sheet are parameterized in terms of snow accumulation, mean annual air temperatur and mean July air temperature. Melt rates are calculated using positive degree-days, and firn warming (i.e. the positive deviation of the temperature at 10-15 m depth from the mean annual air temperature) is estimated from the calculated amount of refrozen melt water in the firn. A comparison between observed and calculated melt rates shows that the parameterization provides a reasonable estimate of the present ablation rates in West Greenland between 61°N and 76°N. The average equilibrium line elevation is estimated to be about 1150 m and 1000 m for West and East Greenland respectively, which is several hundred meter lower than previous estimates. However, the total annual ablation from the ice sheet is found to be about 280 km**3 of water per year which agrees well with most other estimates. The melt-rate model predicts significant melting and consequently significant firn warming even at the highest elevations of the South Greenland ice sheet, whereas a large region of central Greenland north of 70° N experiences little or no summer melting. This agrees with the distribution of the dry snow facics as given by BENSON (1962).

Tab. 1: Januar temperatures on location of the high northern hemisphere (Lustrenwerte)

Während sich Mitteleuropa bereits seit 1920 in einer Phase sinkender Januartemperatur befindet, trat dieser Klima-Umschwung in höheren Breiten immer später ein, so waren in Königsberg die beiden Januar Lustren 1916-1920 und 1921-1925 gleich warm, und in Upsala und Haparanda waren die Jahrfünfte mildester Januartemperatur in beiden Fällen 1931-35 - eine Erscheinung von übrigens weltweitem Charakter, denn auch weiteste Teile der U.S.A. hatten in dieser Zeit die mildesten Januarmonate.

Alkenone-derived UK'37 data and sea surface temperatures (SST) of sediment cores from the Gulf of Taranto (Italy)

Auxiliary data include one file with alkenone-derived UK'37 data and sea surface temperatures (SST). On these data Figs. 7 and 8 of the manuscript are based. The SST are derived from UK'37 by using the transfer function: SST = 29.876 UK'37 - 1.334 of Conte et al. (2006). The data are against the ages (in A.D.) of samples derived from cores GT91-1 (39[deg]59'23"N, 17[deg]45'25"E), GT89-3 and GT90-3 (both 39[deg]45'43"N, 17[deg]53'55"E ). Also included are composite records for UK'37 and SST. For creating the composite records, GT-89-3 was taken as reference core. In the overlapping period the GT89-3 data seem in general lower than the GT91-1 data. To accommodate for this in the composite record, the average difference (0.0343 UK'37 units; equivalent to 1.023 [deg]C) was subtracted from the GT91-1 record. Hereafter, for each depth in the overlapping interval the respective values (UK'37 or SST) of GT89-3 and GT91-1 were averaged. We have also averaged with 16 additional alkenone measurements, from 1793 to 1851, performed in the GT90-3 core.

(Table 1) Viscosity coefficients at various temperatures of samples from ODP Hole 176-735B

Ocean Drilling Program (ODP) Hole 735B was drilled to a depth of 1.5 km in a tectonic window of gabbroic lower oceanic crust created at the Southwest Indian Ridge. The gabbros have a very stable natural remanent magnetization (NRM) of reversed polarity with most unblocking temperatures slightly below the Curie temperature of magnetite. The NRM includes a drilling-induced overprint but its intensity decays strongly towards the interior of the drill core. The demagnetization data yield no or only a very small secondary magnetization component acquired during the present Brunhes chron or an earlier normal chron, suggesting cooling through most of the blocking temperature range during chron C5r and a strong resistance against the acquisition of thermoviscous magnetization. A novel furnace has been designed to measure magnetizations and their time dependences at high temperatures (up to 580 deg C) inside a commercial SQUID magnetometer. Magnetic viscosity experiments have been conducted on the gabbros at temperatures up to 550 deg C to determine the time and temperature stability of remanent magnetization. Viscosities are generally small and increase little with temperature below the main blocking temperature, where the increase becomes almost an order of magnitude. Extrapolations to geological times infer viscous acquisitions that would be 5-25% of a thermoremanence in 100 kyr and at temperatures of 200-500 deg C. At ocean bottom temperature the predicted magnetization of one sample acquired in the present Brunhes chron should be 10% of the NRM. However, this is not recognized during NRM demagnetization and partial thermoremanent magnetization (pTRM) acquisitions at 250 deg C are also much smaller than predicted. It thus appears that the NRMs are generally magnetically harder than magnetizations acquired after heating to 570 deg C in the laboratory. Susceptibility changes during heating are small (<5%) indicating a seemingly stable magneto-mineralogy, but conspicuous minima occur after heating to 520 deg C. Also, quasi paleointensity experiments reveal characteristic patterns in the NRM/pTRM ratios and also large increases in pTRM capacity after heating to 570 deg C. Moreover, anhysteretic remanent magnetization acquisition in the low field range (<=10 mT) is strongly enhanced after heating by factors up to three. The alteration of the magneto-mineralogy is interpreted to result from the annealing of defects in magnetite that originate from tectonically induced strain. The oceanic gabbros of Hole 735B are thus ideal source layer material for marine magnetic anomalies, and secondary thermoviscous acquisition, as a possible cause for anomalous skewness, is essentially absent.

Calculated sea level records and deep ocean temperatures from DSDP and ODP sites

During the mid-Pleistocene transition the dominant 41 ka periodicity of glacial cycles transitioned to a quasi-100 ka periodicity for reasons not yet known. This study investigates the potential role of deep ocean hydrography by examining oxygen isotope ratios in benthic foraminifera. Oxygen isotope records from the Atlantic, Pacific and Indian Ocean basins are separated into their ice volume and local temperature/hydrography components using a piece-wise linear transfer function and a temperature calibration. Although our method has certain limitations, the deep ocean hydrography reconstructions show that glacial deep ocean temperatures approached freezing point as the mid-Pleistocene transition progressed. Further analysis suggests that water mass reorganisation could have been responsible for these temperature changes, leading to such stable conditions in the deep ocean that some obliquity cycles were skipped until precessional forcing triggered deglaciation, creating the apparent quasi-100 ka pattern. This study supports previous work that suggests multiples of obliquity cycles dominate the quasi-100 ka glacial cycles with precession components driving deglaciations.