Appendix A. Whole-rock major, trace and rare earth element analyses of the Ulubey volcanics, Appendix B. Analytical standards and detection limits for analyses

Collisional and post-collisional volcanic rocks in the Ulubey (Ordu) area at the western edge of the Eastern Pontide Tertiary Volcanic Province (EPTVP) in NE Turkey are divided into four suites; Middle Eocene (49.4-44.6 Ma) aged Andesite-Trachyandesite (AT), Trachyandesite-Trachydacite-Rhyolite (TTR), Trachydacite-Dacite (TD) suites, and Middle Miocene (15.1 Ma) aged Trachybasalt (TB) suite. Local stratigraphy in the Ulubey area starts with shallow marine environment sediments of the Paleocene-Eocene time and then continues extensively with sub-aerial andesitic to rhyolitic and rare basaltic volcanism during Eocene and Miocene time, respectively. Petrographically, the volcanic rocks are composed primarily of andesites/trachyandesites, with minor trachydacites/rhyolites, basalts/trachybasalts and pyroclastics, and show porphyric, hyalo-microlitic porphyric and rarely glomeroporphyric, intersertal, intergranular, fluidal and sieve textures. The Ulubey (Ordu) volcanic rocks indicate magma evolution from tholeiitic-alkaline to calc-alkaline with medium-K contents. Primitive mantle normalized trace element and chondrite normalized rare earth element (REE) patterns show that the volcanic rocks have moderate light rare earth element (LREE)/heavy rare earth element (HREE) ratios relative to E-Type MORB and depletion in Nb, Ta and Ti. High Th/Yb ratios indicate parental magma(s) derived from an enriched source formed by mixing of slab and asthenospheric melts previously modified by fluids and sediments from a subduction zone. All of the volcanic rocks share similar incompatible element ratios (e.g., La/Sm, Zr/Nb, La/Nb) and chondrite-normalized REE patterns, indicating that the basic to acidic rocks originated from the same source. The volcanic rocks were produced by the slab dehydration-induced melting of an existing metasomatized mantle source, and the fluids from the slab dehydration introduced significant large ion lithophile element (LILE) and LREE to the source, masking its inherent HFSE-enriched characteristics. The initial 87Sr/86Sr (0.7044-0.7050) and eNd (-0.3 to +3.4) ratios of the volcanics suggest that they originated from an enriched lithospheric mantle source with low Sm/Nd ratios. Integration of the geochemical, petrological and isotopical with regional and local geological data suggest that the Tertiary volcanic rocks from the Ulubey (Ordu) area were derived from an enriched mantle, which had been previously metasomatized by fluids derived from subducted slab during Eocene to Miocene in collisional and post-collisional extension-related geodynamic setting following Late Mesozoic continental collision between the Eurasian plate and the Tauride-Anatolide platform.

Metallic Fluence Monitors and J-value standards

Utilizing the neutron-irradiation parameter J is one of the major uncertainties in 40Ar/39Ar dating. The associated error of the individual J-value for a sample of unknown age depends on the accuracy of the age of the geological standards, the fast-neutron fluence distribution in the reactor and the distances between standards and samples during irradiation. While it is generally assumed that rotating irradiation evens out radial neutron fluence gradients, we observed axial and radial variations of the J-values in sample irradiations in the rotating channels of two reactors. To quantify them, we included three-dimensionally distributed metallic fast- (Ni) and thermal- (Co) neutron fluence monitors in three irradiations and geological age standards in three more. Two irradiations were carried out under Cd-shielding in the FRG1 reactor in Geesthacht, Germany, and four without Cd-shielding in the LVR-15 reactor in Rez, Czech Republic. The 58Ni(nf,p)58Co activation reaction and ?-spectrometry of the 811 keV peak associated with the subsequent decay of 58Co to 58Fe allow to calculate the fast-neutron fluence. The fast-neutron fluences at known positions in the irradiation container correlate with the J-values determined by mass-spectrometric 40Ar/39Ar measurements of the geological age standards. Ra-dial neutron fluence gradients are up to 1.8 %/cm in FRG1 and up to 2.2 %/cm in LVR-15; the corre-sponding axial gradients are up to 5.9 and 2.1 %/cm. We conclude that sample rotation might not al-ways suffice to meet the needs of high-precision dating and gradient monitoring can be crucial.

In-situ Iron isotope measurements and mineral compositions from ODP Hole 206-1256D and IODP Hole 335-U1256D

In-situ Fe isotope measurements have been carried out to estimate the impact of the hydrothermal metamorphic overprint on the Fe isotopic composition of Fe-Ti-oxides and Fe-sulfides of the different lithologies of the drilled rocks from IODP Hole 1256D (eastern equatorial Pacific; 15 Ma crust formed at the East Pacific Rise). Most igneous rocks normally have a very restricted range in their 56Fe/54Fe ratio. In contrast, Fe isotope compositions of hot fluids (> 300 °C) from mid-ocean-ridge spreading centers define a narrow range that is shifted to lower delta 56Fe values by 0.2 per mil - 0.5 per mil as compared to igneous rocks. Therefore, it is expected that mineral phases that contain large amounts of Fe are especially affected by the interaction with a fluid that fractionates Fe isotopes during exsolution/precipitation of those minerals. We have used a femtosecond UV-Laser ablation system to determine mineral 56Fe/54Fe ratios of selected samples with a precision of < 0.1 per mil (2 sigma level) at micrometer-scale. We have found significant variations of the delta 56Fe (IRMM-014) values in the minerals between different samples as well as within samples and mineral grains. The overall observed scale of delta 56Fe (magnetite) in 1256D rocks ranges from - 0.12 to + 0.64 per mil, and of delta 56Fe (ilmenite) from - 0.77 to + 0.01 per mil. Pyrite in the lowermost sheeted dike section is clearly distinguishable from the other investigated lithological units, having positive delta 56Fe values between + 0.29 and + 0.56 per mil, whereas pyrite in the other samples has generally negative delta 56Fe values from - 1.10 to - 0.59 permil. One key observation is that the temperature dependent inter-mineral fractionations of Fe isotopes between magnetite and ilmenite are systematically shifted towards higher values when compared to theoretically expected values, while synthesized, well equilibrated magnetite-ilmenite pairs are compatible with the theoretical predictions. Theoretical considerations including beta-factors of different aqueous Fe-chlorides and Rayleigh-type fractionations in the presence of a hydrous, chlorine-bearing fluid can explain this observation. The disagreement between observed and theoretical equilibrium fractionation, the fact that magnetite, in contrast to ilmenite shows a slight downhole trend in the delta 56Fe values, and the observation of small scale heterogeneities within single mineral grains imply that a general re-equilibration of the magnetite-ilmenite pairs is overprinted by kinetic fractionation effects, caused by the interaction of magnetite/ilmenite with hydrothermal fluids penetrating the upper oceanic crust during cooling, or incomplete re-equilibration at low temperatures. Furthermore, the observation of significant small-scale variations in the 56Fe/54Fe ratios of single minerals in this study highlights the importance of high spatial-resolution-analyses of stable isotope ratios for further investigations.

Water temperature reconstructions based on paired Uk'37 and Tex86 records, and paired surface and subsurface Mg/Ca records

This dataset contains the collection of available published paired Uk'37 and Tex86 records spanning multi-millennial to multi-million year time scales, as well as a collection of Mg/Ca-derived temperatures measured in parallel on surface and subsurface dwelling foraminifera, both used in the analyses of Ho and Laepple, Nature Geoscience 2016. As the signal-to-noise ratios of proxy-derived Holocene temperatures are relatively low, we selected records that contain at least the last deglaciation (oldest sample >18kyr BP).