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.

Impact of high pCO2 and warmer temperatures on the process of silica biomineralization in the sponge Mycale grandis

Siliceous sponges have survived pre-historical mass extinction events caused by ocean acidification and recent studies suggest that siliceous sponges will continue to resist predicted increases in ocean acidity. In this study, we monitored silica biomineralization in the Hawaiian sponge Mycale grandis under predicted pCO2 and sea surface temperature scenarios for 2100. Our goal was to determine if spicule biomineralization was enhanced or repressed by ocean acidification and thermal stress by monitoring silica uptake rates during short-term (48 h) experiments and comparing biomineralized tissue ratios before and after a long-term (26 d) experiment. In the short-term experiment, we found that silica uptake rates were not impacted by high pCO2 (1050 µatm), warmer temperatures (27°C), or combined high pCO2 with warmer temperature (1119 µatm; 27°C) treatments. The long-term exposure experiments revealed no effect on survival or growth rates of M. grandis to high pCO2 (1198 µatm), warmer temperatures (25.6°C), or combined high pCO2 with warmer temperature (1225 µatm, 25.7°C) treatments, indicating that M. grandis will continue to prosper under predicted increases in pCO2 and sea surface temperature. However, ash-free dry weight to dry weight ratios, subtylostyle lengths, and silicified weight to dry weight ratios decreased under conditions of high pCO2 and combined pCO2 warmer temperature treatments. Our results show that rising ocean acidity and temperature have marginal negative effects on spicule biomineralization and will not affect sponge survival rates of M. grandis.

(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).

Sedimentology of cores from the Iberian margin

Bioturbation in marine sediments has basically two aspects of interest for palaeo-environmental studies. First, the traces left by the burrowing organisms reflect the prevailing environmental conditions at the seafloor and thus can be used to reconstruct the ecologic and palaeoceanographic situation. Traces have the advantage over other proxies of practically always being preserved in situ. Secondly, for high- resolution stratigraphy, bioturbation is a nuisance due to the stirring and mixing processes that destroy the stratigraphic record. In order to evaluate the applicability of biogenic traces as palaeoenvironmental indicators, a number of gravity cores from the Portuguese continental slope, covering the period from the last glacial to the present were investigated through X-ray radiographs. In addition, physical and chemical parameters were determined to define the environmental niche in each core interval. A number of traces could be recognized, the most important being: Thalassinoides, Planolites, Zoophycos, Chondrites, Scolicia, Palaeophycus, Phycosiphon and the generally pyritized traces Trichichnus and Mycellia. The shifts between the different ichnofabrics agree strikingly well with the variations in ocean circulation caused by the changing climate. On the upper and middle slope, variations in current intensity and oxygenation of the Mediterranean Outflow Water were responsible for shifts in the ichnofabric. Larger traces such as Planolites and Thalassinoides dominated in coarse, well oxygenated intervals, while small traces such as Chondrites and Trichichnus dominated in fine grained, poorly oxygenated intervals. In contrast, on the lower slope where calm steady sedimentation conditions prevail, changes in sedimentation rate and nutrient flux have controlled variations in the distribution of larger traces such as Planolites, Thalassinoides, and Palaeophycus. Additionally, distinct layers of abundant Chondrites correspond to Heinrich events 1, 2, and 4, and are interpreted as a response to incursions of nutrient rich, oxygen depleted Antarctic waters during phases of reduced thermohaline circulation. The results clearly show that not one single factor but a combination of several factors is necessary to explain the changes in ichnofabric. Furthermore, large variations in the extent and type of bioturbation and tiering between different settings clearly show that a more detailed knowledge of the factors governing bioturbation is necessary if we shall fully comprehend how proxy records are disturbed. A first attempt to automatize a part of the recognition and quantification of the ichnofabric was performed using the DIAna image analysis program on digitized X-ray radiographs. The results show that enhanced abundance of pyritized microburrows appears to be coupled to organic rich sediments deposited under dysoxic conditions. Coarse grained sediments inhibit the formation of pyritized burrows. However, the smallest changes in program settings controlling the grey scale threshold and the sensitivity resulted in large shifts in the number of detected burrows. Therefore, this method can only be considered to be semi-quantitative. Through AMS-^C dating of sample pairs from the Zoophycos spreiten and the surrounding host sediment, age reversals of up to 3,320 years could be demonstrated for the first time. The spreiten material is always several thousands of years younger than the surrounding host sediment. Together with detailed X-ray radiograph studies this shows that the trace maker collects the material on the seafloor, and then transports it downwards up to more than one meter in to the underlying sediment where it is deposited in distinct structures termed spreiten. This clearly shows that age reversals of several thousands of years can be expected whenever Zoophycos is unknowingly sampled. These results also render the hitherto proposed ethological models proposed for Zoophycos as largely implausible. Therefore, a combination of detritus feeding, short time caching, and hibernation possibly combined also with gardening, is suggested here as an explanation for this complicated burrow.

Impact of atmospheric pCO2, seawater temperature, and calcification rate on the delta 18O and delta 13C composition of echinoid calcite (Echinometra viridis)

The tropical echinoid Echinometra viridis was reared in controlled laboratory experiments at temperatures of approximately 20°C and 30°C to mimic winter and summer temperatures and at carbon dioxide (CO2) partial pressures of approximately 487 ppm-v and 805 ppm-v to simulate current and predicted-end-of-century levels. Spine material produced during the experimental period and dissolved inorganic carbon (DIC) of the corresponding culture solutions were then analyzed for stable oxygen (delta 18Oe, delta 18ODIC) and carbon (The tropical echinoid Echinometra viridis was reared in controlled laboratory experiments at temperatures of approximately 20°C and 30°C to mimic winter and summer temperatures and at carbon dioxide (CO2) partial pressures of approximately 487 ppm-v and 805 ppm-v to simulate current and predicted-end-of-century levels. Spine material produced during the experimental period and dissolved inorganic carbon (DIC) of the corresponding culture solutions were then analyzed for stable oxygen (delta18Oe, delta18ODIC) and carbon (delta13Ce, delta13CDIC) isotopic composition. Fractionation of oxygen stable isotopes between the echinoid spines and DIC of their corresponding culture solutions (delta18O = delta18Oe - delta18ODIC) was significantly inversely correlated with seawater temperature but not significantly correlated with atmospheric pCO2. Fractionation of carbon stable isotopes between the echinoid spines and DIC of their corresponding culture solutions (Delta delta13C = delta13Ce - delta13CDIC) was significantly positively correlated with pCO2 and significantly inversely correlated with temperature, with pCO2 functioning as the primary factor and temperature moderating the pCO2-delta13C relationship. Echinoid calcification rate was significantly inversely correlated with both delta18O and delta13C, both within treatments (i.e., pCO2 and temperature fixed) and across treatments (i.e., with effects of pCO2 and temperature controlled for through ANOVA). Therefore, calcification rate and potentially the rate of co-occurring dissolution appear to be important drivers of the kinetic isotope effects observed in the echinoid spines. Study results suggest that echinoid delta18O monitors seawater temperature, but not atmospheric pCO2, and that echinoid delta13C monitors atmospheric pCO2, with temperature moderating this relationship. These findings, coupled with echinoids' long and generally high-quality fossil record, supports prior assertions that fossil echinoid delta18O is a viable archive of paleo-seawater temperature throughout Phanerozoic time, and that delta13C merits further investigation as a potential proxy of paleo-atmospheric pCO2. However, the apparent impact of calcification rate on echinoid delta18O and delta13C suggests that paleoceanographic reconstructions derived from these proxies in fossil echinoids could be improved by incorporating the effects of growth rate.13Ce, delta13CDIC) isotopic composition. Fractionation of oxygen stable isotopes between the echinoid spines and DIC of their corresponding culture solutions (delta18O = delta18Oe - delta18ODIC) was significantly inversely correlated with seawater temperature but not significantly correlated with atmospheric pCO2. Fractionation of carbon stable isotopes between the echinoid spines and DIC of their corresponding culture solutions (delta13C = delta13Ce - delta13CDIC) was significantly positively correlated with pCO2 and significantly inversely correlated with temperature, with pCO2 functioning as the primary factor and temperature moderating the pCO2-delta13C relationship. Echinoid calcification rate was significantly inversely correlated with both delta18O and delta13C, both within treatments (i.e., pCO2 and temperature fixed) and across treatments (i.e., with effects of pCO2 and temperature controlled for through ANOVA). Therefore, calcification rate and potentially the rate of co-occurring dissolution appear to be important drivers of the kinetic isotope effects observed in the echinoid spines. Study results suggest that echinoid delta18O monitors seawater temperature, but not atmospheric pCO2, and that echinoid delta13C monitors atmospheric pCO2, with temperature moderating this relationship. These findings, coupled with echinoids' long and generally high-quality fossil record, supports prior assertions that fossil echinoid delta18O is a viable archive of paleo-seawater temperature throughout Phanerozoic time, and that delta13C merits further investigation as a potential proxy of paleo-atmospheric pCO2. However, the apparent impact of calcification rate on echinoid delta18O and delta13C suggests that paleoceanographic reconstructions derived from these proxies in fossil echinoids could be improved by incorporating the effects of growth rate.

Maintained larval growth in mussel larvae exposed to acidified undersaturated seawater

Ocean acidification (OA) is known to affect bivalve early life-stages. We tested responses of blue mussel larvae to a wide range of pH in order to identify their tolerance threshold. Our results confirmed that decreasing seawater pH and decreasing saturation state increases larval mortality rate and the percentage of abnormally developing larvae. Virtually no larvae reared at average pHT 7.16 were able to feed or reach the D-shell stage and their development appeared to be arrested at the trochophore stage. However larvae were capable of reaching the D-shell stage under milder acidification (pHT=7.35, 7.6, 7.85) including in under-saturated seawater with omega Aragonite as low as 0.54±0.01 (mean±s. e. m.), with a tipping point for normal development identified at pHT 7.765. Additionally growth rate of normally developing larvae was not affected by lower pHT despite potential increased energy costs associated with compensatory calcification in response to increased shell dissolution. Overall, our results on OA impacts on mussel larvae suggest an average pHT of 7.16 is beyond their physiological tolerance threshold and indicate a shift in energy allocation towards growth in some individuals revealing potential OA resilience.