Auxiliary material and basic geochemical and benthic foraminiferal stable isotope data for the interval bracketing Eocene Thermal Maximum 2 (ETM2) at DSDP Sites 401 and 550 in the northeastern Atlantic Ocean

Ever since its discovery, Eocene Thermal Maximum 2 (ETM2; ~53.7 Ma) has been considered as one of the "little brothers" of the Paleocene-Eocene Thermal Maximum (PETM; ~56 Ma) as it displays similar characteristics including abrupt warming, ocean acidification, and biotic shifts. One of the remaining key questions is what effect these lesser climate perturbations had on ocean circulation and ventilation and, ultimately, biotic disruptions. Here we characterize ETM2 sections of the NE Atlantic (Deep Sea Drilling Project Sites 401 and 550) using multispecies benthic foraminiferal stable isotopes, grain size analysis, XRF core scanning, and carbonate content. The magnitude of the carbon isotope excursion (0.85-1.10 per mil) and bottom water warming (2-2.5°C) during ETM2 seems slightly smaller than in South Atlantic records. The comparison of the lateral d13C gradient between the North and South Atlantic reveals that a transient circulation switch took place during ETM2, a similar pattern as observed for the PETM. New grain size and published faunal data support this hypothesis by indicating a reduction in deepwater current velocity. Following ETM2, we record a distinct intensification of bottom water currents influencing Atlantic carbonate accumulation and biotic communities, while a dramatic and persistent clay reduction hints at a weakening of the regional hydrological cycle. Our findings highlight the similarities and differences between the PETM and ETM2. Moreover, the heterogeneity of hyperthermal expression emphasizes the need to specifically characterize each hyperthermal event and its background conditions to minimalize artifacts in global climate and carbonate burial models for the early Paleogene.

Helium isotope ratios and sedimentation rate models of ODP holes

A continuous age model for the brief climate excursion at the Paleocene-Eocene boundary has been constructed by assuming a constant flux of extraterrestrial 3He (3He[ET]) to the seafloor. 3He[ET] measurements from ODP Site 690 provide quantitative evidence for the rapid onset (models indicating extremely rapid release of isotopically light carbon, possibly from seafloor methane hydrate, as the proximal cause of the event. However, the 3He[ET] technique indicates a previously unrecognized and extreme increase in sedimentation rate coincident with the return of climate proxies to pre-event values. The 3He[ET]-based age model thus suggests a far more rapid recovery from the climatic perturbation than previously proposed or predicted on the basis of the modern carbon cycle, and so may indicate additional or accelerated mechanisms of carbon removal from the ocean-atmosphere system during this period. 3He[ET] was also measured at ODP Site 1051 to test the validity of the Site 690 chronology. Comparison of these data sets seems to require removal of several tens of kyr of sediment within the climatic excursion at Site 1051, an observation consistent with sediment structures and previous age modeling efforts. The Site 1051 age model shows a ~30 kyr period in which climate proxies return toward pre-event values, after which they remain invariant for ~80 kyr. If this rise represents the recovery interval identified at Site 690, then the 3HeET-based age models of the two sites are in good agreement. However, alternative interpretations are possible, and work on less disrupted sites is required to evaluate the reliability of the proposed new chronology of the climate excursion. Regardless of these details, this work shows that the 3HeET technique can provide useful independent evidence for the development and testing of astronomically calibrated age models.

Thermal reaction norms of growth in Atlantic and Pacific silverside fishes

How organisms may adapt to rising global temperatures is uncertain, but concepts can emerge from studying adaptive physiological trait variations across existing spatial climate gradients. Many ectotherms, particularly fish, have evolved increasing genetic growth capacities with latitude (i.e. countergradient variation (CnGV) in growth), which are thought to be an adaptation primarily to strong gradients in seasonality. In contrast, evolutionary responses to gradients in mean temperature are often assumed to involve an alternative mode, 'thermal adaptation'. We measured thermal growth reaction norms in Pacific silverside populations (Atherinops affinis) occurring across a weak latitudinal temperature gradient with invariant seasonality along the North American Pacific coast. Instead of thermal adaptation, we found novel evidence for CnGV in growth, suggesting that CnGV is a ubiquitous mode of reaction-norm evolution in ectotherms even in response to weak spatial and, by inference, temporal climate gradients. A novel, large-scale comparison between ecologically equivalent Pacific versus Atlantic silversides (Menidia menidia) revealed how closely growth CnGV patterns reflect their respective climate gradients. While steep growth reaction norms and increasing growth plasticity with latitude in M. menidia mimicked the strong, highly seasonal Atlantic coastal gradient, shallow reaction norms and much smaller, latitude-independent growth plasticity in A. affinis resembled the weak Pacific latitudinal temperature gradient.

Age models and summer sea surface temperature and winter sea ice concentration for the EPILOG-LGM time slice in the Pacific Southern Ocean

Sea surface temperatures and sea-ice extent are the most critical variables to evaluate the Southern Ocean paleoceanographic evolution in relation to the development of the global carbon cycle, atmospheric CO2 variability and ocean-atmosphere circulation. In contrast to the Atlantic and the Indian sectors, the Pacific sector of the Southern Ocean has been insufficiently investigated so far. To cover this gap of information we present diatom-based estimates of summer sea surface temperature (SSST) and winter sea-ice concentration (WSI) from 17 sites in the polar South Pacific to study the Last Glacial Maximum (LGM) at the EPILOG time slice (19,000-23,000 cal. years BP). Applied statistical methods are the Imbrie and Kipp Method (IKM) and the Modern Analog Technique (MAT) to estimate temperature and sea-ice concentration, respectively. Our data display a distinct LGM east-west differentiation in SSST and WSI with steeper latitudinal temperature gradients and a winter sea-ice edge located consistently north of the Pacific-Antarctic Ridge in the Ross sea sector. In the eastern sector of our study area, which is governed by the Amundsen Abyssal Plain, the estimates yield weaker latitudinal SSST gradients together with a variable extended winter sea-ice field. In this sector, sea-ice extent may have reached sporadically the area of the present Subantarctic Front at its maximum LGM expansion. This pattern points to topographic forcing as major controller of the frontal system location and sea-ice extent in the western Pacific sector whereas atmospheric conditions like the Southern Annular Mode and the ENSO affected the oceanographic conditions in the eastern Pacific sector. Although it is difficult to depict the location and the physical nature of frontal systems separating the glacial Southern Ocean water masses into different zones, we found a distinct temperature gradient in latitudes straddled by the modern Southern Subtropical Front. Considering that the glacial temperatures north of this zone are similar to the modern, we suggest that this represents the Glacial Southern Subtropical Front (GSSTF), which delimits the zone of strongest glacial SSST cooling (>4K) to its North. The southern boundary of the zone of maximum cooling is close to the glacial 4°C isotherm. This isotherm, which is in the range of SSST at the modern Antarctic Polar Front (APF), represents a circum-Antarctic feature and marks the northern edge of the glacial Antarctic Circumpolar Current (ACC). We also assume that a glacial front was established at the northern average winter sea ice edge, comparable with the modern Southern Antarctic Circumpolar Current Front (SACCF). During the glacial, this front would be located in the area of the modern APF. The northward deflection of colder than modern surface waters along the South American continent leads to a significant cooling of the glacial Humboldt Current surface waters (4-8K), which affects the temperature regimes as far north as into tropical latitudes. The glacial reduction of ACC temperatures may also result in the significant cooling in the Atlantic and Indian Southern Ocean, thus may enhance thermal differentiation of the Southern Ocean and Antarctic continental cooling. Comparison with temperature and sea ice simulations for the last glacial based on numerical simulations show that the majority of modern models overestimate summer and winter sea ice cover and that there exists few models that reproduce our temperature data rather well.

Extraterrestrial 3He estimation across the Paleocene-Eocene thermal maximum at ODP Site 208-1266

In the deep-sea, the Paleocene-Eocene Thermal Maximum (PETM) is often marked by clay-rich condensed intervals caused by dissolution of carbonate sediments, capped by a carbonate-rich interval. Constraining the duration of both the dissolution and subsequent cap-carbonate intervals is essential to computing marine carbon fluxes and thus testing hypotheses for the origin of this event. To this end, we provide new high-resolution helium isotope records spanning the Paleocene-Eocene boundary at ODP Site 1266 in the South Atlantic. The extraterrestrial 3He, 3HeET, concentrations replicate trends observed at ODP Site 690 by Farley and Eltgroth (2003, doi:10.1016/S0012-821X(03)00017-7). By assuming a constant flux of 3HeET we constrain relative changes in accumulation rates of sediment across the PETM and construct a new age model for the event. In this new chronology the zero carbonate layer represents 35 kyr, some of which reflects clay produced by dissolution of Paleocene (pre-PETM) sediments. Above this layer, carbonate concentrations increase for ~165 kyr and remain higher than in the latest Paleocene until 234 +48/-34 kyr above the base of the clay. The new chronology indicates that minimum d13C values persisted for a maximum of 134 +27/-19 kyr and the inflection point previously chosen to designate the end of the CIE recovery occurs at 217 +44/-31 kyr. This allocation of time differs from that of the cycle-based age model of Röhl et al. (2007, doi:10.1029/2007GC001784) in that it assigns more time to the clay layer followed by a more gradual recovery of carbonate-rich sedimentation. The new model also suggests a longer sustained d13C excursion followed by a more rapid recovery to pre-PETM d13C values. These differences have important implications for constraining the source(s) of carbon and mechanisms for its subsequent sequestration, favoring models that include a sustained release

Stable isotopes and XRF data over the Paleocene-Eocene Thermal Maximum from IODP Sites U1403 and U1409

During the Paleocene-Eocene Thermal Maximum (PETM) about 56 million years ago, thousands of petagrams of carbon were released into the atmosphere and ocean in just a few thousand years, followed by a gradual sequestration over approximately 200,000 years. If silicate weathering is one of the key negative feedbacks that removed this carbon, a period of seawater calcium carbonate saturation greater than pre-event levels is expected during the event's recovery phase. In marine sediments, this should be recorded as a temporary deepening of the depth below which no calcite is preserved - the calcite compensation depth (CCD). Previous and new sedimentary records from sites that were above the pre-PETM calcite compensation depth show enhanced carbonate accumulation following the PETM. A new record from an abyssal site in the North Atlantic that lay below the pre-PETM calcite compensation depth shows a period of carbonate preservation beginning about 70,000 years after the onset of the PETM, providing the first direct evidence for an over-deepening of the calcite compensation depth. This record confirms an overshoot in ocean carbonate saturation during the PETM recovery. Simulations with two earth system models support scenarios for the PETM that involve both a large initial carbon release followed by prolonged low-level emissions, consistent with the timing of CCD deepening in our record. Our findings indicate that sequestration of these carbon emissions was most likely the result of both globally enhanced calcite burial above the calcite compensation depth and, at least in the North Atlantic, by a temporary over-deepening of the calcite compensation depth.

Experimental data on photophysiology and growth rates of the invasive foraminifera Amphistegina lobifera showing high thermal tolerance from the Eastern Mediterranean and the Gulf of Aqaba

Invasive species allow an investigation of trait retention and adaptations after exposure to new habitats. Recent work on corals from the Gulf of Aqaba (GoA) shows that tolerance to high temperature persists thousands of years after invasion, without any apparent adaptive advantage. Here we test whether thermal tolerance retention also occurs in another symbiont-bearing calcifying organism. To this end, we investigate the thermal tolerance of the benthic foraminifera Amphistegina lobifera from the GoA (29° 30.14167 N 34° 55.085 E) and compare it to a recent "Lessepsian invader population" from the Eastern Mediterranean (EaM) (32° 37.386 N, 34°55.169 E). We first established that the studied populations are genetically homogenous but distinct from a population in Australia, and that they contain a similar consortium of diatom symbionts, confirming their recent common descent. Thereafter, we exposed specimens from GoA and EaM to elevated temperatures for three weeks and monitored survivorship, growth rates and photophysiology. Both populations exhibited a similar pattern of temperature tolerance. A consistent reduction of photosynthetic dark yields was observed at 34°C and reduced growth was observed at 32°C. The apparent tolerance to sustained exposure to high temperature cannot have a direct adaptive importance, as peak summer temperatures in both locations remain <32°C. Instead, it seems that in the studied foraminifera tolerance to high temperature is a conservative trait and the EaM population retained this trait since its recent invasion. Such pre-adaptation to higher temperatures confers A. lobifera a clear adaptive advantage in shallow and episodically high temperature environments in the Mediterranean under further warming.

Thermal conductivity of ODP Leg 110 holes

Thirty-four sediment and mudline temperatures were collected from six drill holes on ODP Leg 110 near the toe of the Barbados accretionary complex. When combined with thermal conductivity measurements these data delineate the complicated thermal structure on the edge of this convergent margin. Surface heat-flow values from Leg 110 (calculated from geothermal gradients forced through the bottom-water temperature at mudline) of 92 to 192 mW/m**2 are 80% to 300% higher than values predicted by standard heat flow vs. age models for oceanic crust, but are compatible with earlier surface measurements made at the same latitude. Measured heat flow tends to decrease downhole at four sites, suggesting the presence of heat sources within the sediments. These results are consistent with the flow of warm fluid through the complex along sub-horizontal, high-permeability conduits, including thrust faults, the major decollement zone, and sandy intervals. Simple calculations suggest that this flow is transient, occurring on time scales of tens to tens of thousands of years. High heat flow in the vicinity of 15°30'N and not elsewhere along the deformation front suggests that the Leg 110 drill sites may be situated over a fluid discharge zone, with dewatering more active here than elsewhere along the accretionary complex.