Palaeoenvironment and palaeoclimate of the Middle Miocene lake in the Steinheim basin, SW Germany: A reconstruction from C, O, and Sr isotopes of fossil remains

A differentiated reconstruction of palaeolimnologic, -environmental, and -climatic conditions is presented for the Middle Miocene long-term freshwater lake (14.3 to 13.5 Ma) of the Steinheim basin, on the basis of a combined C, 0, and Sr isotope study of sympatric skeletal fossils of aquatic and terrestrial organisms from the lake sediments. The oxygen isotope composition for lake water of the Steinheim basin (delta O-18(H2O) = +2.0 +/- 0.4 parts per thousand VSMOW, n = 6) was reconstructed from measurements of delta O-18(PO4) of aquatic turtle bones. The drinking water calculated from the enamel of large mammals (proboscideans, rhinocerotids, equids, cervids, suids) has delta O-18(H2O) values (delta(OH2O)-O-18 = -5.9 +/- 1.7 parts per thousand VSMOW, n = 31) typical for Middle Miocene meteoric water of the area. This delta O-18(H2O) value corresponds to a mean annual air temperature (MAT) of 18.8 +/- 3.8 degrees C, calculated using a modem-day delta(OH2O)-O-18-MAT relation. Hence, large mammals did not use the lake water as principal drinking water. In contrast, small mammals, especially the then abundant pika Prolagus oeningensis drank from O-18-enriched water sources (delta O-18(H2O) = +2.7 +/- 2.3 parts per thousand VSMOW, n = 7), such as the lake water. Differences in Sr and 0 isotopic compositions between large and small mammal teeth indicate different home ranges and drinking behaviour and support migration of some large mammals between the Swabian Alb plateau and the nearby Molasse basin, while small mammals ingested their food and water locally. Changes in the lake level, water chemistry, and temperature were inferred using isotopic compositions of ostracod and gastropod shells from a composite lake sediment profile. Calcitic ostracod valves (Ilyocypris binocularis; delta O-18 = +1.7 +/- 1.2 parts per thousand VPDB, delta C-18 = -0.5 +/- 0.9 parts per thousand, VPDB, n = 68) and aragonitic, gastropod shells (Gyraulus spp.; delta O-18 = +2.0 +/- 13 parts per thousand VPDB, delta C-13 = -1.1 +/- 1.3 parts per thousand VPDB, n = 89) have delta O-18 and delta C-13 values similar to or even higher than those of marine, carbonates. delta C-13 values:of the biogenic carbonates parallel lake level fluctuations while delta O-18 values scatter around +2 +/- 2 parts per thousand and reflect the short term variability of meteoric water inflow vs. longer term evaporation. Sr-87/Sr-86 ratios of aragonitic Gyraulus spp. gastropod shells parallel the lake level fluctuations, reflecting variable inputs of groundwater and surface waters. Using a water delta O-18(H2O) value of +2.0 parts per thousand VSMOW, water temperatures calculated from skeletal tissue delta O-18 values of ostracods are 16.7 +/- 5.0 degrees C, gastropods 20.6 +/- 5.6 degrees C, otoliths 21.8 +/- 1.4 degrees C, and fish teeth 17.0 +/- 2.7 degrees C. The calculated MAT (similar to 19 degrees C), lake water temperatures (similar to 17 to 22 degrees C) and the O-18-enriched water compositions are indicative of warm-temperate climatic conditions, possibly with a high humidity during this period. Vegetation in the area surrounding the basin was largely of the C-3-type, as indicated by carbon isotopic compositions of tooth enamel from large mammals (delta C-13 = -11.1 +/- 1.1 parts per thousand VPDB, n = 40). (c) 2006 Elsevier B.V. All rights reserved.

Assessing sedimentary evolution by means of Sr-isotope ratios : 3 case studies on the Caribbean plate : (Cretaceous: Nicoya Peninsula, Costa Rica, Tertiary: Hess Rise, and La Désirade, Guadeloupe, France)

The understanding of sedimentary evolution is intimately related to the knowledge of the exact ages of the sediments. When working on carbonate sediments, age dating is commonly based on paleontological observations and established biozonations, which may prove to be relatively imprecise. Dating by means of strontium isotope ratios in marine bioclasts is the probably best method in order to precisely date carbonate successions, provided that the sample reflects original marine geochemical characteristics. This requires a precise study of the samples including its petrography, SEM and cathodoluminescence observations, stable carbon and oxygen isotope geochemistry and finally the strontium isotope measurement itself. On the Nicoya Peninsula (Northwestern Costa Rica) sediments from the Piedras Blancas Formation, Nambi Formation and Quebrada Pavas Formation were dated by the means of strontium isotope ratios measured in Upper Cretaceous Inoceramus shell fragments. Results have shown average 87Sr/86Sr values of 0.707654 (middle late Campanian) for the Piedras Blancas Formation, 0.707322 (Turonian-Coniacian) for the Nambi Formation and 0.707721 (late Campanian-Maastrichtian) for the Quebrada Pavas Formation. Abundant detrital components in the studied formations constitute a difficulty to strontium isotope dating. In fact, the fossil bearing sediments can easily contaminate the target fossil with strontium mobilized form basalts during diagenesis and thus the obtained strontium isotope ratios may be influenced significantly and so will the obtained ages. The new and more precise age assignments allow for more precision in the chronostratigraphic chart of the sedimentary and tectonic evolution of the Nicoya Peninsula, providing a better insight on the evolution of this region. Meteor Cruise M81 dredged shallow water carbonates from the Hess Rise and Hess Escarpment during March 2010. Several of these shallow water carbonates contain abundant Larger Foraminifera that indicates an Eocene-Oligocene age. In this study the strontium isotope values ranging from 0.707847 to 0.708238 can be interpreted as a Rupelian to Chattian age of these sediments. These platform sediments are placed on seamounts, now located at depths reaching 1600 m. Observation of sedimentologic characteristics of these sediments has helped to resolve apparent discrepancies between fossil and strontium isotope ages. Hence, it is possible to show that the subsidence was active during early Miocene times. On La Désirade (Guadeloupe France), the Neogene to Quaternary carbonate cover has been dated by microfossils and some U/Th-ages. Disagreements subsisted in the paleontological ages of the formations. Strontium isotope ratios ranging from 0.709047 to 0.709076 showed the Limestone Table of La Désirade to range from an Early Pliocene to Late Pliocene/early Pleistocene age. A very late Miocene age (87Sr/86Sr =0.709013) can be determined to the Detrital Offshore Limestone. The flat volcanic basement had to be eroded by wave-action during a long-term stable relative sea-level. Sediments of the Table Limestone on La Désirade show both low-stand and high-stand facies that encroach on the igneous basement, implying deposition during a major phase of subsidence creating accommodation space. Subsidence is followed by tectonic uplift documented by fringing reefs and beach rocks that young from the top of the Table Limestone (180 m) towards the present coastline. Strontium isotope ratios from two different fringing reefs (0.707172 and 0.709145) and from a beach rock (0.709163) allow tentative dating, (125ky, ~ 400ky, 945ky) and indicate an uplift rate of about 5cm/ky for this time period of La Désirade Island. The documented subsidence and uplift history calls for a new model of tectonic evolution of the area.

Characterization of modern and fossil mineral dust transported to high altitude in the Western Alps: Saharan sources and transport patterns

Mineral dust aerosols recently collected at the high-altitude Jungfraujoch research station (46 degrees 33'51 `' N, 7 degrees 59'06 `' E; 3580 m a.s.l.) were compared to mineral dust deposited at the Colle Gnifetti glacier (45 degrees 52'50 `' N, 7 degrees 52'33 `' E; 4455 m a.s.l.) over the last millennium. Radiogenic isotope signatures and backward trajectories analyses indicate that major dust sources are situated in the north-central to north-western part of the Saharan desert. Less radiogenic Sr isotopic compositions of PM10 aerosols and of mineral particles deposited during periods of low dust transfer likely result from the enhancement of the background chemically-weathered Saharan source. Saharan dust mobilization and transport were relatively reduced during the second part of the Little Ice Age (ca. 1690-1870) except within the greatest Saharan dust event deposited around 1770. After ca. 1870, sustained dust deposition suggests that increased mineral dust transport over the Alps during the last century could be due to stronger spring/summer North Atlantic southwesterlies and drier winters in North Africa. On the other hand, increasing carbonaceous particle emissions from fossil fuel combustion combined to a higher lead enrichment factor point to concomitant anthropogenic sources of particulate pollutants reaching high-altitude European glaciers during the last century.

PyRate: a new program to estimate speciation and extinction rates from incomplete fossil data

Despite the advancement of phylogenetic methods to estimate speciation and extinction rates, their power can be limited under variable rates, in particular for clades with high extinction rates and small number of extant species. Fossil data can provide a powerful alternative source of information to investigate diversification processes. Here, we present PyRate, a computer program to estimate speciation and extinction rates and their temporal dynamics from fossil occurrence data. The rates are inferred in a Bayesian framework and are comparable to those estimated from phylogenetic trees. We describe how PyRate can be used to explore different models of diversification. In addition to the diversification rates, it provides estimates of the parameters of the preservation process (fossilization and sampling) and the times of speciation and extinction of each species in the data set. Moreover, we develop a new birth-death model to correlate the variation of speciation/extinction rates with changes of a continuous trait. Finally, we demonstrate the use of Bayes factors for model selection and show how the posterior estimates of a PyRate analysis can be used to generate calibration densities for Bayesian molecular clock analysis. PyRate is an open-source command-line Python program available at http://sourceforge.net/projects/pyrate/.

Bayesian estimation of speciation and extinction from incomplete fossil occurrence data.

The temporal dynamics of species diversity are shaped by variations in the rates of speciation and extinction, and there is a long history of inferring these rates using first and last appearances of taxa in the fossil record. Understanding diversity dynamics critically depends on unbiased estimates of the unobserved times of speciation and extinction for all lineages, but the inference of these parameters is challenging due to the complex nature of the available data. Here, we present a new probabilistic framework to jointly estimate species-specific times of speciation and extinction and the rates of the underlying birth-death process based on the fossil record. The rates are allowed to vary through time independently of each other, and the probability of preservation and sampling is explicitly incorporated in the model to estimate the true lifespan of each lineage. We implement a Bayesian algorithm to assess the presence of rate shifts by exploring alternative diversification models. Tests on a range of simulated data sets reveal the accuracy and robustness of our approach against violations of the underlying assumptions and various degrees of data incompleteness. Finally, we demonstrate the application of our method with the diversification of the mammal family Rhinocerotidae and reveal a complex history of repeated and independent temporal shifts of both speciation and extinction rates, leading to the expansion and subsequent decline of the group. The estimated parameters of the birth-death process implemented here are directly comparable with those obtained from dated molecular phylogenies. Thus, our model represents a step towards integrating phylogenetic and fossil information to infer macroevolutionary processes.

Effects of season, sex and body size on the feeding ecology of turtle-headed seasnakes (Emydocephalus annulatus) on IndoPacific inshore coral reefs

In terrestrial snakes, many cases of intraspecific shifts in dietary habits as a function of predator sex and body size are driven by gape-limitation - and hence, are most common in species that feed on relatively large prey, and exhibit a wide body-size range. Our data on seasnakes reveal an alternative mechanism for intraspecific niche partitioning, based on sex-specific seasonal anorexia induced by reproductive activities. Turtle-headed seasnakes (Emydocephalus annulatus) on coral reefs in the New Caledonian Lagoon feed entirely on the eggs of demersal-spawning fishes. DNA sequence data (cytochrome b gene) on eggs that we palpated from stomachs of 37 snakes showed that despite this ontogenetic-stage specialization, the prey come from a taxonomically diverse array of species including damselfish (41% of samples, at least 5 species), blennies (41%, 4 species) and gobies (19%, 5 species). The composition of snake diets shifted seasonally (with damselfish dominating in winter but not summer), presumably reflecting seasonality of fish reproduction. That seasonal shift affects male and female snakes differently, because reproduction is incompatible with foraging. Adult female seasnakes ceased feeding when they became heavily distended with developing embryos in late summer, and males ceased feeding while they were mate-searching in winter. The sex divergence in foraging habits may be amplified by sexual size dimorphism; females grow larger than males, and larger snakes (of both sexes) feed more on damselfish (which often lay their eggs in exposed sites) than on blennies and gobies (whose eggs are hidden within narrow crevices). Specific features of reproductive biology of coral-reef fish (seasonality and nest type) have generated intraspecific niche partitioning in these seasnakes, by mechanisms different from those that apply to terrestrial snakes.

Nd and Sr isotope compositions in modern and fossil bones - Proxies for vertebrate provenance and taphonomy

Rare earth elements (REE), while not essential for the physiologic functions of animals, are ingested and incorporated in ppb concentrations in bones and teeth. Nd isotope compositions of modern bones of animals from isotopically distinct habitats demonstrate that the (143)Nd/(144)Nd of the apatite can be used as a fingerprint for bedrock geology or ambient water mass. This potentially allows the provenance and migration of extant vertebrates to be traced, similar to the use of Sr isotopes. Although REE may be enriched by up to 5 orders of magnitude during diagenesis and recrystallization of bone apatite, in vivo (143)Nd/(144)Nd may be preserved in the inner cortex of fossil bones or enamel. However, tracking the provenance of ancient or extinct vertebrates is possible only for well-preserved archeological and paleontological skeletal remains with in vivo-like Nd contents at the ppb-level. Intra-bone and -tooth REE analysis can be used to screen for appropriate areas. Large intra-bone Nd concentration gradients of 10(1)-10(3) are often measured. Nd concentrations in the inner bone cortex increase over timescales of millions of years, while bone rims may be enriched over millenial timescales. Nevertheless, epsilon(Nd) values are often similar within one epsilon(Nd) unit within a single bone. Larger intra-bone differences in specimens may either reflect a partial preservation of in vivo values or changing epsilon(Nd) values of the diagenetic fluid during fossilization. However, most fossil specimens and the outer rims of bones will record taphonomic (143)Nd/(144)Nd incorporated post mortem during diagenesis. Unlike REE patterns, (143)Nd/(144)Nd are not biased by fractionation processes during REE-uptake into the apatite crystal lattice, hence the epsilon(Nd) value is an important tracer for taphonomy and reworking. Bones and teeth from autochthonous fossil assemblages have small variations of +/- 1 epsilon(Nd) unit only. In contrast, fossil bones and teeth from over 20 different marine and terrestrial fossil sites have a total range of epsilon(Nd) values from -13.0 to 4.9 (n = 80), often matching the composition of the embedding sediment. This implies that the surrounding sediment is the source of Nd in the fossil bones and that the specimens of this study seem not to have been reworked. Differences in epsilon(Nd) values between skeletal remains and embedding sediment may either indicate reworking of fossils and/or a REE-uptake from a diagenetic fluid with non-sediment derived epsilon(Nd) values. The latter often applies to fossil shark teeth, which may preserve paleo-seawater values. Complementary to epsilon(Nd) values, (87)Sr/(86)Sr can help to further constrain the fossil provenance and reworking. (C) 2011 Elsevier Ltd. All rights reserved.

Accessory phases as a tool to investigate magmatic underplating in polymetamorphic metasediments : an example from a fossil lower crust : Ivrea Verbano zne, NW Italy

The Ivrea and the Strona-Ceneri zones, NW italy and S Switzerland, offer the possibility to study the continental crust of the Southern Alps. Because of its high metamorphic degree and the abundant Permo- Carboniferous mafic intrusions, the Ivrea Zone is classically interpreted an exposed section trough the Permian lower crust. The present work is focused here on metasedimentary slices (septa) intercalated within Permian gabbro (mafic complex). In particular I studied the evolution of accessory phases such as rutile and zircon and the chemistry of the metasediments. The septa build an irregular and discontinuous band that cut obliquely the mafic complex from its deepest part (N) to its roof (S). The chemistry of the metasediments evolves along the band and the chemical evolution can be compared with that observed in the country-rock surrounding the mafic intrusion to the NE and overprinted by a main regional metamorphic event. This suggests that the degree of chemical depletion of the septa was mainly established during the same regional metamorphic event. Moreover it suggests that incorporation of the septa within the gabbro did not modify their original stratigraphie distribution within the crust. It implies that the mafic complex has been emplaced following a dynamic substantially different from the classic model of « gabbro glacier » (Quick et al., 1992; Quick et al., 1994). It is more likely that it has been emplaced by repeated injections of sills at different depths during a protracted period of time. Zircon trace elements and U-Pb ages suggest that regional metamorphism occurred 330-320Ma, the first sills in the deepest part of the Mafic Complex are injected at ~300Ma, the mafic magmas reached higher levels in the crust at 285Ma and the magmatic activity continued locally until 275Ma. The ages of detrital cores in zircons fix the maximal sedimentation age at ~370Ma, this age corresponds therefore with the maximal age of the incorporation of the Ivrea zone within the lower crust. I propose that the Ivrea zone has been accreted to the lower crust during the Hercynian orogeny sensu lato. The analysis of detrital ages suggests that the source terrains for the Ivrea zone and those for the Strona-Ceneri zone have a completely different Palaeozoic history. The systematic analysis of rutile in partially molten metasediments of the Ivrea zone reveals the occurrence of two generations. The two generations are characterized by a different chemistry and textural distribution. A first generation is formed during pro-grade metamorphism in the restitic counterpart. The second generation is formed in the melts during cooling at the same time that part of the first generation re-equilibrate. Re-equilibration of the first generation seems to be spatially controlled by the presence of fluids. Locally the second generation forms overgrowths on the first generation. Considered the different diffusivity of U and Pb in rutile, U heterogeneities have important implication for U-Pb dating of rutile. ID-TIMS and LA-ICPMS dating coupled with a careful textural investigation (SEM) suggest that rutile grains are characterized by multiple path along which Pb diffusion can occur: volume diffusion is an important process, but intragrain and subgrain boundaries provide additional high diffusivity pathways for Pb escape and reduce drastically the effective diffusion length. -- La zone d'Ivrea et la zone de Strona-Ceneri, en Italie nord-occidentale et Suisse méridionale, offrent la possibilité d'étudier la croûte continentale des Alpes du Sud. En raison du haut degré métamorphique et l'abondance d'intrusions mafiques d'âge Permo-Carbonifère [complexe mafique), la zone d'Ivrea est interprétée classiquement comme de la croûte inférieure permienne. Ce travail ce concentre sur des bandes metasédimentaires (septa) incorporées dans les magmas mafiques lors de l'intrusion. Les septa forment une bande irrégulière qui coupe obliquement le complexe mafique du bas (N) vers le haut (S). La chimie des septa évolue du bas vers le haut et l'évolution chimique se rapproche de l'évolution observé dans la roche encaissante l'intrusion affecté par un événement métamorphique régionale. Cette relation suggère que le degré d'appauvrissement chimique des septa a été établit principalement lors de l'événement métamorphique régional. De plus l'incorporation dans les gabbros n'a pas perturbée la distribution stratigraphique originelle des septa. Ces deux observations impliquent que le métamorphisme dans la roche encaissante précède la mise en place du gabbro et que cette dernière ne se fait pas selon le modèle classique (« gabbro glacier » de Quick et al., 1992, 1994), mais se fait plutôt par injections répétées de sills a différentes profondeurs. Les âges U-Pb et les éléments traces des zircons suggèrent que le métamorphisme régionale a eu lieu 330-320Ma, alors que les premiers sills dans la partie profonde du Mafic Complex s'injectent à ~300Ma, le magmatisme mafique atteigne des niveaux supérieurs à 285Ma et continue localement jusqu'à 270Ma. Les âges des coeurs détritiques des zircons permettent de fixer l'âge maximale de sédimentation à ~370Ma ce qui correspond donc à l'âge maximale de l'incorporation de la zone d'Ivrea dans la croûte inférieur. L'analyse systématique des rutiles, nous a permit de montrer l'existence de plusieurs générations qui ont une répartition texturale et une chimie différente. Une génération se forme lors de l'événement UHT dans les restites, une autre génération se forme dans les liquides lors du refroidissement, au même temps qu'une partie de la première génération se rééquilibre au niveau du Zr. Localement la deuxième génération peut former des surcroissances autour de la première génération. Dans ces cas, des fortes différences en uranium entre les deux générations ont des importantes implications pour la datation U-Pb sur rutile. Classiquement les ratios Pb/U dans le rutile sont interprétés comme indiquant l'âges du refroidissement du minéral sous une température à la quelle la diffusion du Pb dans le minéral n'est plus détectable et la diffusion à plus hautes températures est assumée se faire par «volume diffusion» dans le grain (Mezger et al., 1989). Par des datations ID-TIMS (sur grain entier) et LA-ICPMS (in-situ) et une analyse texturale (MEB) approfondie nous montrons que cette supposition est trop simpliste et que le rutile est repartie en sous-domaines. Chacun de ces domaines a ça propre longueur ou chemin de diffusion spécifique. Nous proposons donc une nouvelle approche plus cohérente pour l'interprétation des âges U-Pb sur rutile.