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.

(Table 1) Boron-isotope measurements of shallow-dwelling planktonic foraminifera of ODP Holes 144-817A, 144-872C and 143-865C

Knowledge of the evolution of atmospheric carbon dioxide concentrations throughout the Earth's history is important for a reconstruction of the links between climate and radiative forcing of the Earth's surface temperatures. Although atmospheric carbon dioxide concentrations in the early Cenozoic era (about 60 Myr ago) are widely believed to have been higher than at present, there is disagreement regarding the exact carbon dioxide levels, the timing of the decline and the mechanisms that are most important for the control of CO2 concentrations over geological timescales. Here we use the boron-isotope ratios of ancient planktonic foraminifer shells to estimate the pH of surface-layer sea water throughout the past 60 million years, which can be used to reconstruct atmospheric CO2 concentrations. We estimate CO2 concentrations of more than 2,000 p.p.m. for the late Palaeocene and earliest Eocene periods (from about 60 to 52 Myr ago), and find an erratic decline between 55 and 40 Myr ago that may have been caused by reduced CO2 outgassing from ocean ridges, volcanoes and metamorphic belts and increased carbon burial. Since the early Miocene (about 24 Myr ago), atmospheric CO2 concentrations appear to have remained below 500 p.p.m. and were more stable than before, although transient intervals of CO2 reduction may have occurred during periods of rapid cooling approximately 15 and 3 Myr ago.

Geometric parameters of shoreface profiles of Laptev Sea erosion coast (Table 1)

Field investigations of the Laptev Sea shoreface morphology were carried out (1) off erosional shores composed of unconsolidated sediments, (2) off the modern delta shores of the Lena River, and (3) off rocky shores. It was found that profiles off erosional shores had a concave shape. This shape is not well described by commonly applied power functions, a feature, which is in disagreement with the generally accepted concept of the equilibrium shape of shoreface profiles. The position of the lower shoreface boundary is determined by the elevation of the coastal lowland inundated during the last transgression (at -5 to -10 m) and may easily be recognized by a sharp, an order of magnitude decrease in the mean inclination of the sea floor. The mean shoreface inclination depends on sediment grain-size and ranges from 0.0022 to 0.033. The concave shape of the shoreface did not change substantially during the last 20-30 years, which indicates that shoreline retreat did not slow down and hence suggests continued intensive coastal erosion in the 21st century. The underwater part of the Lena River delta extends up to 35 km offshore. Its upper part is formed by a shallow and up to 18-km wide bench, which reaches depths of 2-3 m along the outer edge. The evolution of the delta was irregular. Whereas some parts of the delta are advancing rapidly (58 m/year), other parts are eroding. Comparison of measured profiles with older bathymetric data gave an opportunity to evaluate the changes of the underwater delta over past decades. Bathymetric surveys of the seabed around the delta can thus contribute towards a quantification of the sediment budget of the river-sea system. In addition, some sections of the Laptev Sea coast are composed of bedrock that has a comparatively low resistance to wave erosion. These sections may supply a considerable amount of sediment, especially if the cliffs are high. This source must therefore also be taken into account when assessing the contribution of shore erosion to the Laptev Sea sediment budget.

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.

Benthic foraminifera across the Cretaceous/Paleogene boundary in ODP Hole 198-1210A

Deep-sea benthic foraminifera show important but transient assemblage changes at the Cretaceous/Paleogene (K/Pg) boundary, when many biota suffered severe extinction. We quantitatively analyzed benthic foraminiferal assemblages from lower bathyal-upper abyssal (1500-2000 m) northwest Pacific ODP Site 1210 (Shatsky Rise) and compared the results with published data on assemblages at lower bathyal (~ 1500 m) Pacific DSDP Site 465 (Hess Rise) to gain insight in paleoecological and paleoenvironmental changes at that time. At both sites, diversity and heterogeneity rapidly decreased across the K/Pg boundary, then recovered. Species assemblages at both sites show a similar pattern of turnover from the uppermost Maastrichtian into the lowermost Danian: 1) The relative abundance of buliminids (indicative of a generally high food supply) increases towards the uppermost Cretaceous, and peaks rapidly just above the K/Pg boundary, coeval with a peak in benthic foraminiferal accumulation rate (BFAR), a proxy for food supply. 2) A peak in relative abundance of Stensioeina beccariiformis, a cosmopolitan form generally more common at the middle than at the lower bathyal sites, occurs just above the buliminid peak. 3) The relative abundance of Nuttallides truempyi, a more oligotrophic form, decreases at the boundary, then increases above the peak in Stensioeina beccariiformis. The food supply to the deep sea in the Pacific Ocean thus apparently increased rather than decreased in the earliest Danian. The low benthic diversity during a time of high food supply indicates a stressed environment. This stress might have been caused by reorganization of the planktic ecosystem: primary producer niches vacated by the mass extinction of calcifying nannoplankton may have been rapidly (<10 kyr) filled by other, possibly opportunistic, primary producers, leading to delivery of another type of food, and/or irregular food delivery through a succession of opportunistic blooms. The deep-sea benthic foraminiferal data thus are in strong disagreement with the widely accepted hypothesis that the global deep-sea floor became severely food-depleted following the K/Pg extinction due to the mass extinction of primary producers ("Strangelove Ocean Model") or to the collapse of the biotic pump ("Living Ocean Model").

Global compilation of last glacial maximum oxygen isotopic ratios for planktonic and benthic foraminifera

We use the fully coupled atmosphere-ocean three-dimensional model of intermediate complexity iLOVECLIM to simulate the climate and oxygen stable isotopic signal during the Last Glacial Maximum (LGM, 21 000 yr). By using a model that is able to explicitly simulate the sensor (d18O), results can be directly compared with data from climatic archives in the different realms. Our results indicate that iLOVECLIM reproduces well the main feature of the LGM climate in the atmospheric and oceanic components. The annual mean d18O in precipitation shows more depleted values in the northern and southern high latitudes during the LGM. The model reproduces very well the spatial gradient observed in ice core records over the Greenland ice-sheet. We observe a general pattern toward more enriched values for continental calcite d18O in the model at the LGM, in agreement with speleothem data. This can be explained by both a general atmospheric cooling in the tropical and subtropical regions and a reduction in precipitation as confirmed by reconstruction derived from pollens and plant macrofossils. Data-model comparison for sea surface temperature indicates that iLOVECLIM is capable to satisfyingly simulate the change in oceanic surface conditions between the LGM and present. Our data-model comparison for calcite d18O allows investigating the large discrepancies with respect to glacial temperatures recorded by different microfossil proxies in the North Atlantic region. The results argue for a trong mean annual cooling between the LGM and present (>6°C), supporting the foraminifera transfer function reconstruction but in disagreement with alkenones and dinocyst reconstructions. The data-model comparison also reveals that large positive calcite d18O anomaly in the Southern Ocean may be explained by an important cooling, although the driver of this pattern is unclear. We deduce a large positive d18Osw anomaly for the north Indian Ocean that contrasts with a large negative d18Osw anomaly in the China Sea between the LGM and present. This pattern may be linked to changes in the hydrological cycle over these regions. Our simulation of the deep ocean suggests that changes in d18Osw between the LGM and present are not spatially homogenous. This is supported by reconstructions derived from pore fluids in deep-sea sediments. The model underestimates the deep ocean cooling thus biasing the comparison with benthic calcite d18O data. Nonetheless, our data-model comparison support a heterogeneous cooling of few degrees (2-4°C) in the LGM Ocean.

Orbital-resolution Pliocene-Pleistocene simulations of CO2, sea level, global temperature and benthic d18O, and d11B-based proxy-CO2 data

During the past five million yrs, benthic d18O records indicate a large range of climates, from warmer than today during the Pliocene Warm Period to considerably colder during glacials. Antarctic ice cores have revealed Pleistocene glacial-interglacial CO2 variability of 60-100 ppm, while sea level fluctuations of typically 125 m are documented by proxy data. However, in the pre-ice core period, CO2 and sea level proxy data are scarce and there is disagreement between different proxies and different records of the same proxy. This hampers comprehensive understanding of the long-term relations between CO2, sea level and climate. Here, we drive a coupled climate-ice sheet model over the past five million years, inversely forced by a stacked benthic d18O record. We obtain continuous simulations of benthic d18O, sea level and CO2 that are mutually consistent. Our model shows CO2 concentrations of 300 to 470 ppm during the Early Pliocene. Furthermore, we simulate strong CO2 variability during the Pliocene and Early Pleistocene. These features are broadly supported by existing and new d11B-based proxy CO2 data, but less by alkenone-based records. The simulated concentrations and variations therein are larger than expected from global mean temperature changes. Our findings thus suggest a smaller Earth System Sensitivity than previously thought. This is explained by a more restricted role of land ice variability in the Pliocene. The largest uncertainty in our simulation arises from the mass balance formulation of East Antarctica, which governs the variability in sea level, but only modestly affects the modeled CO2 concentrations.

Ice rafted debris and reconstructed sea surface temperature of ODP Hole 104-642B

The mid-Piacenzian warm period (3.264-3.025 Ma) of the Pliocene epoch has been proposed as a possible reference for future warm climate states. However, there is significant disagreement over the magnitude of high latitude warming between data and models for this period of time, raising questions about the driving mechanisms and responsible feedbacks. We have developed a new set of orbital-resolution alkenone-based sea surface temperature (SST) and ice rafted debris (IRD) records from the Norwegian Sea spanning 3.264-3.14 Ma. The SSTs in the Norwegian Sea were 2-3?°C warmer than the Holocene average, likely caused by the radiative effect of higher atmospheric CO2 concentrations. There is notable obliquity-driven SST variability with a range of 4?°C, shown by evolutive spectra. The correlation of SST variability with the presence of IRD suggests a common climate forcing acting across the Nordic Seas region. Changes of the SST gradient between the Norwegian Sea and North Atlantic sites suggest that the subpolar gyre was at least as strong as during the Holocene, and that the northward heat transport by the North Atlantic Current was comparable.