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.

Elevation lines of modelled ice sheet extension during last glacial maximum

A numerical ice-sheet model was used to reconstruct the Late Weichselian glaciation of the Eurasian High Arctic, between Franz Josef Land and Severnaya Zemlya. An ice sheet was developed over the entire Eurasian High Arctic so that ice flow from the central Barents and Kara seas toward the northern Russian Arctic could be accounted for. An inverse approach to modeling was utilized, where ice-sheet results were forced to be compatible with geological information indicating ice-free conditions over the Taymyr Peninsula during the Late Weichselian. The model indicates complete glaciation of the Barents and Kara seas and predicts a "maximum-sized" ice sheet for the Late Weichselian Russian High Arctic. In this scenario, full-glacial conditions are characterized by a 1500-m-thick ice mass over the Barents Sea, from which ice flowed to the north and west within several bathymetric troughs as large ice streams. In contrast to this reconstruction, a "minimum" model of glaciation involves restricted glaciation in the Kara Sea, where the ice thickness is only 300 m in the south and which is free of ice in the north across Severnaya Zemlya. Our maximum reconstruction is compatible with geological information that indicates complete glaciation of the Barents Sea. However, geological data from Severnaya Zemlya suggest our minimum model is more relevant further east. This, in turn, implies a strong paleoclimatic gradient to colder and drier conditions eastward across the Eurasian Arctic during the Late Weichselian.

Reconstructed glaciation limits in Eurasia (digitized from map Fig.1)

Based on field investigations in northern Russia and interpretation of offshore seismic data, we have made a preliminary reconstruction of the maximum ice-sheet extent in the Barents and Kara Sea region during the Early/Middle Weichselian and the Late Weichselian. Our investigations indicate that the Barents and Kara ice sheets attained their maximum Weichselian positions in northern Russia prior to 50 000 yr BP, whereas the northeastern flank of the Scandinavian Ice Sheet advanced to a maximum position shortly after 17000 calendar years ago. During the Late Weichselian (25 000-10000 yr BP), much of the Russian Arctic remained ice-free. According to our reconstruction, the extent of the ice sheets in the Barents and Kara Sea region during the Late Weichselian glacial maximum was less than half that of the maximum model which, up to now, has been widely used as a boundary condition for testing and refining General Circulation Models (GCMs). Preliminary numerical-modelling experiments predict Late Weichselian ice sheets which are larger than the ice extent implied for the Kara Sea region from dated geological evidence, suggesting very low precipitation.

(Table 1) Age determination of Cape Shpindler, Yugorski Peninsula, northwest Russia

Depositional environments, stratigraphic relations, and 35 new AMS 14C dates at Cape Shpindler, Yugorski Peninsula, help constrain the late Pleistocene glacial and environmental history of the southern Kara Sea region. Fifteen- to fifty-meter-high coastal exposures reveal a complex package of shallow marine, fluvial, glacial, and postglacial deposits, and are documented here in a 19-km-long cross-section and eight vertical sections. The shallow marine (Unit A), estuarine or prodeltaic (Unit B), and fluvio-deltaic (Unit C) deposits contain an interglacial molluscan fauna, yield radiocarbon dates greater than 40 ka, and may correspond with a regional sea-level highstand during the Eemian. These units are overlain by a diamicton (Unit D), and are pervasively deformed by folds and low- to high-angle faults into a stacked glaciotectonic accretionary complex. The diamicton (Unit D) is a subglacial till, and associated massive ground ice with deformed debris bands (Unit E) appears to be relict glacier ice. Glaciotectonic structures document both southward- and northward-directed glacier movement. Above the till and associated glaciotectonic horizons lies 0- to 11-m-thick postglacial deposits of peatland, eolian, fluvial, and primarily lacustrine origin (Unit F). The postglacial deposits yield radiocarbon ages of 12.8 to 0.8 ka. Thus, at least one regional glaciation is prominently represented in the stratigraphy, and occurred probably after the Eemian but before 12.8 ka. We infer that the bulk of the glacial record corresponds with southward advance by an early Weichselian Kara Sea Ice Sheet, in agreement with other recently documented, regional records from Yamal Peninsula and the Pechora Basin. The timing and source of northward-directed glacier ice are less well constrained. Across the broad expanse of the Eurasian Arctic, Quaternary stratigraphy is still sparsely documented. The new data from Cape Shpindler fill a spatial gap in paleoenvironmental research.