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.

Crete fault scarp point cloud data

This paper assesses the along strike variation of active bedrock fault scarps using long range terrestrial laser scanning (t-LiDAR) data in order to determine the distribution behaviour of scarp height and the subsequently calculate long term throw-rates. Five faults on Cretewhich display spectacular limestone fault scarps have been studied using high resolution digital elevation model (HRDEM) data. We scanned several hundred square metres of the fault system including the footwall, fault scarp and hanging wall of the investigated fault segment. The vertical displacement and the dip of the scarp were extracted every metre along the strike of the detected fault segment based on the processed HRDEM. The scarp variability was analysed by using statistical and morphological methods. The analysis was done in a geographical information system (GIS) environment. Results show a normal distribution for the scanned fault scarp's vertical displacement. Based on these facts, the mean value of height was chosen to define the authentic vertical displacement. Consequently the scarp can be divided into above, below and within the range of mean (within one standard deviation) and quantify the modifications of vertical displacement. Therefore, the fault segment can be subdivided into areas which are influenced by external modification like erosion and sedimentation processes. Moreover, to describe and measure the variability of vertical displacement along strike the fault, the semi-variance was calculated with the variogram method. This method is used to determine how much influence the external processes have had on the vertical displacement. By combining of morphological and statistical results, the fault can be subdivided into areas with high external influences and areas with authentic fault scarps, which have little or no external influences. This subdivision is necessary for long term throw-rate calculations, because without this differentiation the calculated rates would be misleading and the activity of a fault would be incorrectly assessed with significant implications for seismic hazard assessment since fault slip rate data govern the earthquake recurrence. Furthermore, by using this workflow areas with minimal external influences can be determined, not only for throw-rate calculations, but also for determining samples sites for absolute dating techniques such as cosmogenic nuclide dating. The main outcomes of this study include: i) there is no direct correlation between the fault's mean vertical displacement and dip (R² less than 0.31); ii) without subdividing the scanned scarp into areas with differing amounts of external influences, the along strike variability of vertical displacement is ±35%; iii) when the scanned scarp is subdivided the variation of the vertical displacement of the authentic scarp (exposed by earthquakes only) is in a range of ±6% (the varies depending on the fault from 7 to 12%); iv) the calculation of the long term throw-rate (since 13 ka) for four scarps in Crete using the authentic vertical displacement is 0.35 ± 0.04 mm/yr at Kastelli 1, 0.31 ± 0.01 mm/yr at Kastelli 2, 0.85 ± 0.06 mm/yr at the Asomatos fault (Sellia) and 0.55 ± 0.05 mm/yr at the Lastros fault.

Planktic and benthic foraminifreral isotope data for ODP Site 207-1258

The mid-Cretaceous is widely considered the archetypal ice-free greenhouse interval in Earth history, with a thermal maximum around Cenomanian-Turonian boundary time (ca. 90 Ma). However, contemporaneous glaciations have been hypothesized based on sequence stratigraphic evidence for rapid sea-level oscillation and oxygen isotope excursions in records generated from carbonates of questionable preservation and/or of low resolution. We present new oxygen isotope records for the mid-Cenomanian Demerara Rise that are of much higher resolution than previously available, taken from both planktic and benthic foraminifers, and utilizing only extremely well preserved glassy foraminifers. Our records show no evidence of glaciation, calling into question the hypothesized ice sheets and rendering the origin of inferred rapid sea-level oscillations enigmatic. Simple mass-balance calculations demonstrate that this Cretaceous sea-level paradox is unlikely to be explained by hidden ice sheets existing below the limit of d18O detection.

Homogenization of Portuguese long-term temperature data series: Lisbon, Coimbra and Porto

Three long-term temperature data series measured in Portugal were studied to detect and correct non-climatic homogeneity breaks and are now available for future studies of climate variability. Series of monthly minimum (Tmin) and maximum (Tmax) temperatures measured in the three Portuguese meteorological stations of Lisbon (from 1856 to 2008), Coimbra (from 1865 to 2005) and Porto (from 1888 to 2001) were studied to detect and correct non-climatic homogeneity breaks. These series together with monthly series of average temperature (Taver) and temperature range (DTR) derived from them were tested in order to detect homogeneity breaks, using, firstly, metadata, secondly, a visual analysis and, thirdly, four widely used homogeneity tests: von Neumann ratio test, Buishand test, standard normal homogeneity test and Pettitt test. The homogeneity tests were used in absolute (using temperature series themselves) and relative (using sea-surface temperature anomalies series obtained from HadISST2 close to the Portuguese coast or already corrected temperature series as reference series) modes. We considered the Tmin, Tmax and DTR series as most informative for the detection of homogeneity breaks due to the fact that Tmin and Tmax could respond differently to changes in position of a thermometer or other changes in the instrument's environment; Taver series have been used, mainly, as control. The homogeneity tests show strong inhomogeneity of the original data series, which could have both internal climatic and non-climatic origins. Homogeneity breaks which have been identified by the last three mentioned homogeneity tests were compared with available metadata containing data, such as instrument changes, changes in station location and environment, observing procedures, etc. Significant homogeneity breaks (significance 95% or more) that coincide with known dates of instrumental changes have been corrected using standard procedures. It was also noted that some significant homogeneity breaks, which could not be connected to the known dates of any changes in the park of instruments or stations location and environment, could be caused by large volcanic eruptions. The corrected series were again tested for homogeneity: the corrected series were considered free of non-climatic breaks when the tests of most of monthly series showed no significant (significance 95% or more) homogeneity breaks that coincide with dates of known instrument changes. Corrected series are now available in the frame of ERA-CLIM FP7 project for future studies of climate variability.

Detection histories for eight species of Amazona parrots in Venezuela during the NeoMaps bird surveys in 2010

Documenting changes in distribution is necessary for understanding species' response to environmental changes, but data on species distributions are heterogeneous in accuracy and resolution. Combining different data sources and methodological approaches can fill gaps in knowledge about the dynamic processes driving changes in species-rich, but data-poor regions. We combined recent bird survey data from the Neotropical Biodiversity Mapping Initiative (NeoMaps) with historical distribution records to estimate potential changes in the distribution of eight species of Amazon parrots in Venezuela. Using environmental covariates and presence-only data from museum collections and the literature, we first used maximum likelihood to fit a species distribution model (SDM) estimating a historical maximum probability of occurrence for each species. We then used recent, NeoMaps survey data to build single-season occupancy models (OM) with the same environmental covariates, as well as with time- and effort-dependent detectability, resulting in estimates of the current probability of occurrence. We finally calculated the disagreement between predictions as a matrix of probability of change in the state of occurrence. Our results suggested negative changes for the only restricted, threatened species, Amazona barbadensis, which has been independently confirmed with field studies. Two of the three remaining widespread species that were detected, Amazona amazonica, Amazona ochrocephala, also had a high probability of negative changes in northern Venezuela, but results were not conclusive for Amazona farinosa. The four remaining species were undetected in recent field surveys; three of these were most probably absent from the survey locations (Amazona autumnalis, Amazona mercenaria and Amazona festiva), while a fourth (Amazona dufresniana) requires more intensive targeted sampling to estimate its current status. Our approach is unique in taking full advantage of available, but limited data, and in detecting a high probability of change even for rare and patchily-distributed species. However, it is presently limited to species meeting the strong assumptions required for maximum-likelihood estimation with presence-only data, including very high detectability and representative sampling of its historical distribution.

Icebergs in the Amundsen Sea (Antarctica) determined from ENVISAT ASAR WSM data using object-based image analysis approach, link to ESRI GDB files

ENVISAT ASAR WSM images with pixel size 150 × 150 m, acquired in different meteorological, oceanographic and sea ice conditions were used to determined icebergs in the Amundsen Sea (Antarctica). An object-based method for automatic iceberg detection from SAR data has been developed and applied. The object identification is based on spectral and spatial parameters on 5 scale levels, and was verified with manual classification in four polygon areas, chosen to represent varying environmental conditions. The algorithm works comparatively well in freezing temperatures and strong wind conditions, prevailing in the Amundsen Sea during the year. The detection rate was 96% which corresponds to 94% of the area (counting icebergs larger than 0.03 km**2), for all seasons. The presented algorithm tends to generate errors in the form of false alarms, mainly caused by the presence of ice floes, rather than misses. This affects the reliability since false alarms were manually corrected post analysis.

Collection of data on particulate and dissolved soil nitrogen in the Jena Experiment (Main Experiment, time series since 2002)

This data set contains four time series of particulate and dissolved soil nitrogen measurements from the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below). In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. 1. Total nitrogen from solid phase: Stratified soil sampling was performed every two years since before sowing in April 2002 and was repeated in April 2004, 2006 and 2008 to a depth of 30 cm segmented to a depth resolution of 5 cm giving six depth subsamples per core. In 2002 five samples per plot were taken and analyzed independently. Averaged values per depth layer are reported. In later years, three samples per plot were taken, pooled in the field, and measured as a combined sample. Sampling locations were less than 30 cm apart from sampling locations in other years. All soil samples were passed through a sieve with a mesh size of 2 mm in 2002. In later years samples were further sieved to 1 mm. No additional mineral particles were removed by this procedure. Total nitrogen concentration was analyzed on ball-milled subsamples (time 4 min, frequency 30 s-1) by an elemental analyzer at 1150°C (Elementaranalysator vario Max CN; Elementar Analysensysteme GmbH, Hanau, Germany). 2. Total nitrogen from solid phase (high intensity sampling): In block 2 of the Jena Experiment, soil samples were taken to a depth of 1m (segmented to a depth resolution of 5 cm giving 20 depth subsamples per core) with three replicates per block ever 5 years starting before sowing in April 2002. Samples were processed as for the more frequent sampling but were always analyzed independently and never pooled. 3. Mineral nitrogen from KCl extractions: Five soil cores (diameter 0.01 m) were taken at a depth of 0 to 0.15 m (and between 2002 and 2004 also at a depth of 0.15 to 0.3 m) of the mineral soil from each of the experimental plots at various times over the years. In addition also plots of the management experiment, that altered mowing frequency and fertilized subplots (see further details below) were sampled in some later years. Samples of the soil cores per plot (subplots in case of the management experiment) were pooled during each sampling campaign. NO3-N and NH4-N concentrations were determined by extraction of soil samples with 1 M KCl solution and were measured in the soil extract with a Continuous Flow Analyzer (CFA, 2003-2005: Skalar, Breda, Netherlands; 2006-2007: AutoAnalyzer, Seal, Burgess Hill, United Kingdom). 4. Dissolved nitrogen in soil solution: Glass suction plates with a diameter of 12 cm, 1 cm thickness and a pore size of 1-1.6 µm (UMS GmbH, Munich, Germany) were installed in April 2002 in depths of 10, 20, 30 and 60 cm to collect soil solution. The sampling bottles were continuously evacuated to a negative pressure between 50 and 350 mbar, such that the suction pressure was about 50 mbar above the actual soil water tension. Thus, only the soil leachate was collected. Cumulative soil solution was sampled biweekly and analyzed for nitrate (NO3-), ammonium (NH4+) and total dissolved nitrogen concentrations with a continuous flow analyzer (CFA, Skalar, Breda, The Netherlands). Nitrate was analyzed photometrically after reduction to NO2- and reaction with sulfanilamide and naphthylethylenediamine-dihydrochloride to an azo-dye. Our NO3- concentrations contained an unknown contribution of NO2- that is expected to be small. Simultaneously to the NO3- analysis, NH4+ was determined photometrically as 5-aminosalicylate after a modified Berthelot reaction. The detection limits of NO3- and NH4+ were 0.02 and 0.03 mg N L-1, respectively. Total dissolved N in soil solution was analyzed by oxidation with K2S2O8 followed by reduction to NO2- as described above for NO3-. Dissolved organic N (DON) concentrations in soil solution were calculated as the difference between TDN and the sum of mineral N (NO3- + NH4+).

Gridded bathymetry from multibeam echosounder ELAC BottomChart MkII data of the cruise POS317/3 (2004)

Bathymetry based on data recorded during POS317-3 between 19.09.2004 and 13.10.2004 in the Black Sea. This cruise concentrated on bathymetric mapping and mapping of gas seeps by hydro-acoustic detection of gas flares in the water column and the quantification of microbial turnover of gassy sediments and microbial mats. The major objective during POS317-3 was the characterization and identification of microorganisms involved in the anaerobic methane oxidation in the sediment and in microbial mats. As part of these investigations characteristic organic molecules will be identified, which can be used as biomarkers for anaerobic methane oxidizing microorganisms.

Detection of a variable intracellular acid-labile carbon pool in Thalassiosira weissflogii (Heterokontophyta) and Emiliania huxleyi (Haptophyta) in response to changes in the seawater carbon system

Accumulation of an intracellular pool of carbon (C(i) pool) is one strategy by which marine algae overcome the low abundance of dissolved CO2 (CO2 (aq) ) in modern seawater. To identify the environmental conditions under which algae accumulate an acid-labile C(i) pool, we applied a (14) C pulse-chase method, used originally in dinoflagellates, to two new classes of algae, coccolithophorids and diatoms. This method measures the carbon accumulation inside the cells without altering the medium carbon chemistry or culture cell density. We found that the diatom Thalassiosira weissflogii [(Grunow) G. Fryxell & Hasle] and a calcifying strain of the coccolithophorid Emiliania huxleyi [(Lohmann) W. W. Hay & H. P. Mohler] develop significant acid-labile C(i) pools. C(i) pools are measureable in cells cultured in media with 2-30 µmol/l CO2 (aq), corresponding to a medium pH of 8.6-7.9. The absolute C(i) pool was greater for the larger celled diatoms. For both algal classes, the C(i) pool became a negligible contributor to photosynthesis once CO2 (aq) exceeded 30 µmol/l. Combining the (14) C pulse-chase method and (14) C disequilibrium method enabled us to assess whether E. huxleyi and T. weissflogii exhibited thresholds for foregoing accumulation of DIC or reduced the reliance on bicarbonate uptake with increasing CO2 (aq) . We showed that the C(i) pool decreases with higher CO2 :HCO3 (-) uptake rates.

European mink (Mustela lutreola) and American mink (Neovison vison) detection histories in Northern Spain between 2000 and 2011

Understanding changes over time in the distribution of interacting native and invasive species that may be symptomatic of competitive exclusion is critical to identify the need for and effectiveness of management interventions. Occupancy models greatly increase the robustness of inference that can be made from presence/absence data when species are imperfectly detected, and recent novel developments allow for the quantification of the strength of interaction between pairs of species. We used a two-species multi-season occupancy model to quantify the impact of the invasive American mink on the native European mink in Spain through the analysis of their co-occurrence pattern over twelve years (2000 - 2011) in the entire Spanish range of European mink distribution, where both species were detected by live trapping but American mink were culled. We detected a negative temporal trend in the rate of occupancy of European mink and a simultaneous positive trend in the occupancy of American mink. The species co-occurred less often than expected and the native mink was more likely to become extinct from sites occupied by the invasive species. Removal of American mink resulted in a high probability of local extinction where it co-occurred with the endemic mink, but the overall increase in the probability of occupancy over the last decade indicates that the ongoing management is failing to halt its spread. More intensive culling effort where both species co-exist as well as in adjacent areas where the invasive American mink is found at high densities is required in order to stop thedecline of European mink.