Palladium Isotopic Evidence for Nucleosynthetic and Cosmogenic Isotope Anomalies in IVB Iron Meteorites

The origin of ubiquitous nucleosynthetic isotope anomalies in meteorites may represent spatial and/or temporal heterogeneity in the sources that supplied material to the nascent solar nebula, or enhancement by chemical processing. For elements beyond the Fe peak, deficits in s-process isotopes have been reported in some (e.g., Mo, Ru, W) but not all refractory elements studied (e.g., Os) that, among the iron meteorites, are most pronounced in IVB iron meteorites. Palladium is a non-refractory element in the same mass region as Mo and Ru. In this study, we report the first precise Pd isotopic abundances from IVB irons to test the mechanisms proposed for the origin of isotope anomalies. First, this study determined the existence of a cosmogenic neutron dosimeter from the reaction 103Rh(n, beta-)104Pd in the form of excess 104Pd, correlated with excess 192Pt, in IVB irons. Second, all IVB irons show a deficit of the s-process only isotope 104Pd (\varepsilon 104Pd = -0.48 ± 0.24), an excess of the r-only isotope 110Pd (\varepsilon 110Pd = +0.46 ± 0.12), and no resolvable anomaly in the p-process 102Pd (\varepsilon 102Pd = +1 ± 1). The magnitude of the Pd isotope anomaly is about half that predicted from a uniform depletion of the s-process yields from the correlated isotope anomalies of refractory Mo and Ru. The discrepancy is best understood as the result of nebular processing of the less refractory Pd, implying that all the observed nucleosynthetic anomalies in meteorites are likely to be isotopic relicts. The Mo-Ru-Pd isotope systematics do not support enhanced rates of the 22Ne(alpha,n)25Mg neutron source for the solar system s-process.

Correlated cosmogenic W and Os isotopic variations in Carbo and implications for Hf-W chronology

An obstacle for establishing the chronology of iron meteorite formation using 182Hf-182W systematics (t1/2 = 8.9 Myr) is to find proper neutron fluence monitors to correct for cosmic ray modification of W isotopic composition. Recent studies showed that siderophile elements such as Pt and Os could serve such a purpose. To test and calibrate these neutron dosimeters, the isotopic compositions of W and Os were measured in a slab of the IID iron meteorite Carbo. This slab has a well-characterized noble gas depth profile reflecting different degrees of shielding to cosmic rays. The results show that W and Os isotopic ratios correlate with distance from the pre-atmospheric center. Negative correlations, barely resolved within error, were found between epsilo190Os-epsilo189Os and epsilo186Os-epsilo189Os with slopes of -0.64 ± 0.45 and -1.8(+1.9/-2.1), respectively. These Os isotope correlations broadly agree with model predictions for capture of secondary neutrons produced by cosmic ray irradiation and results reported previously for other groups of iron meteorites. Correlations were also found between epsilo182W-epsilo189Os (slope = 1.02 ± 0.37) and epsilo182W-epsilo190Os (slope = -1.38 ± 0.58). Intercepts of these two correlations yield pre-exposure epsilo182W values of -3.32 ± 0.51 and -3.62 ± 0.23, respectively (weighted average epsilo182W = -3.57 ± 0.21). This value relies on a large extrapolation leading to a large uncertainty but gives a metal-silicate segregation age of -0.5 ± 2.4 Myr after formation of the solar system. Combining the iron meteorite measurements with simulations of cosmogenic effects in iron meteorites, equations are presented to calculate and correct for cosmogenic effects on 182W using Os isotopes.

Environmental correlates of species rank − abundance distributions in global drylands

Theoretical models predict lognormal species abundance distributions (SADs) in stable and productive environments, with log-series SADs in less stable, dispersal driven communities. We studied patterns of relative species abundances of perennial vascular plants in global dryland communities to: (i) assess the influence of climatic and soil characteristics on the observed SADs, (ii) infer how environmental variability influences relative abundances, and (iii) evaluate how colonisation dynamics and environmental filters shape abundance distributions. We fitted lognormal and log-series SADs to 91 sites containing at least 15 species of perennial vascular plants. The dependence of species relative abundances on soil and climate variables was assessed using general linear models. Irrespective of habitat type and latitude, the majority of the SADs (70.3%) were best described by a lognormal distribution. Lognormal SADs were associated with low annual precipitation, higher aridity, high soil carbon content, and higher variability of climate variables and soil nitrate. Our results do not corroborate models predicting the prevalence of log-series SADs in dryland communities. As lognormal SADs were particularly associated with sites with drier conditions and a higher environmental variability, we reject models linking lognormality to environmental stability and high productivity conditions. Instead our results point to the prevalence of lognormal SADs in heterogeneous environments, allowing for more evenly distributed plant communities, or in stressful ecosystems, which are generally shaped by strong habitat filters and limited colonisation. This suggests that drylands may be resilient to environmental changes because the many species with intermediate relative abundances could take over ecosystem functioning if the environment becomes suboptimal for dominant species.

Polytomies and bayesian phylogenetic inference

Bayesian phylogenetic analyses are now very popular in systematics and molecular evolution because they allow the use of much more realistic models than currently possible with maximum likelihood methods. There are, however, a growing number of examples in which large Bayesian posterior clade probabilities are associated with very short edge lengths and low values for non-Bayesian measures of support such as nonparametric bootstrapping. For the four-taxon case when the true tree is the star phylogeny, Bayesian analyses become increasingly unpredictable in their preference for one of the three possible resolved tree topologies as data set size increases. This leads to the prediction that hard (or near-hard) polytomies in nature will cause unpredictable behavior in Bayesian analyses, with arbitrary resolutions of the polytomy receiving very high posterior probabilities in some cases. We present a simple solution to this problem involving a reversible-jump Markov chain Monte Carlo (MCMC) algorithm that allows exploration of all of tree space, including unresolved tree topologies with one or more polytomies. The reversible-jump MCMC approach allows prior distributions to place some weight on less-resolved tree topologies, which eliminates misleadingly high posteriors associated with arbitrary resolutions of hard polytomies. Fortunately, assigning some prior probability to polytomous tree topologies does not appear to come with a significant cost in terms of the ability to assess the level of support for edges that do exist in the true tree. Methods are discussed for applying arbitrary prior distributions to tree topologies of varying resolution, and an empirical example showing evidence of polytomies is analyzed and discussed.

Polytomies and bayesian phylogenetic inference

Bayesian phylogenetic analyses are now very popular in systematics and molecular evolution because they allow the use of much more realistic models than currently possible with maximum likelihood methods. There are, however, a growing number of examples in which large Bayesian posterior clade probabilities are associated with very short edge lengths and low values for non-Bayesian measures of support such as nonparametric bootstrapping. For the four-taxon case when the true tree is the star phylogeny, Bayesian analyses become increasingly unpredictable in their preference for one of the three possible resolved tree topologies as data set size increases. This leads to the prediction that hard (or near-hard) polytomies in nature will cause unpredictable behavior in Bayesian analyses, with arbitrary resolutions of the polytomy receiving very high posterior probabilities in some cases. We present a simple solution to this problem involving a reversible-jump Markov chain Monte Carlo (MCMC) algorithm that allows exploration of all of tree space, including unresolved tree topologies with one or more polytomies. The reversible-jump MCMC approach allows prior distributions to place some weight on less-resolved tree topologies, which eliminates misleadingly high posteriors associated with arbitrary resolutions of hard polytomies. Fortunately, assigning some prior probability to polytomous tree topologies does not appear to come with a significant cost in terms of the ability to assess the level of support for edges that do exist in the true tree. Methods are discussed for applying arbitrary prior distributions to tree topologies of varying resolution, and an empirical example showing evidence of polytomies is analyzed and discussed.

Zircons, isotope ratios and element analyses from tonalite and oxide gabbro of the mid-ocean ridge from ODP Hole 176-735B

A limiting factor in the accuracy and precision of U/Pb zircon dates is accurate correction for initial disequilibrium in the 238U and 235U decay chains. The longest-lived-and therefore most abundant-intermediate daughter product in the 235U isotopic decay chain is 231Pa (T1/2 = 32.71 ka), and the partitioning behavior of Pa in zircon is not well constrained. Here we report high-precision thermal ionization mass spectrometry (TIMS) U-Pb zircon data from two samples from Ocean Drilling Program (ODP) Hole 735B, which show evidence for incorporation of excess 231Pa during zircon crystallization. The most precise analyses from the two samples have consistent Th-corrected 206Pb/238U dates with weighted means of 11.9325 ± 0.0039 Ma (n = 9) and 11.920 ± 0.011 Ma (n = 4), but distinctly older 207Pb/235U dates that vary from 12.330 ± 0.048 Ma to 12.140 ± 0.044 Ma and 12.03 ± 0.24 to 12.40 ± 0.27 Ma, respectively. If the excess 207Pb is due to variable initial excess 231Pa, calculated initial (231Pa)/(235U) activity ratios for the two samples range from 5.6 ± 1.0 to 9.6 ± 1.1 and 3.5 ± 5.2 to 11.4 ± 5.8. The data from the more precisely dated sample yields estimated DPazircon/DUzircon from 2.2-3.8 and 5.6-9.6, assuming (231Pa)/(235U) of the melt equal to the global average of recently erupted mid-ocean ridge basaltic glasses or secular equilibrium, respectively. High precision ID-TIMS analyses from nine additional samples from Hole 735B and nearby Hole 1105A suggest similar partitioning. The lower range of DPazircon/DUzircon is consistent with ion microprobe measurements of 231Pa in zircons from Holocene and Pleistocene rhyolitic eruptions (Schmitt (2007; doi:10.2138/am.2007.2449) and Schmitt (2011; doi:10.1146/annurev-earth-040610-133330)). The data suggest that 231Pa is preferentially incorporated during zircon crystallization over a range of magmatic compositions, and excess initial 231Pa may be more common in zircons than acknowledged. The degree of initial disequilibrium in the 235U decay chain suggested by the data from this study, and other recent high precision datasets, leads to resolvable discordance in high precision dates of Cenozoic to Mesozoic zircons. Minor discordance in zircons of this age may therefore reflect initial excess 231Pa and does not require either inheritance or Pb loss.

Rb, Sr, Sm and Nd concentrations and isotopic composition of basalts from DSDP Leg 74 (Table 1)

Basement intersected in DSDP holes 525A, 528 and 527 on the Walvis Ridge consists of submarine basalt flows and pillows with minor intercalated sediments. These holes are situated on the crest and mid and lower northwest flank of a NNW-SSE-trending ridge block which would have closely paralleled the paleo mid-ocean ridge (Rabinowitz and LaBrecque, 1979 doi:10.1029/JB084iB11p05973, Moore et al. (1983 doi:10.1130/0016-7606(1983)94<907:TWRTDS>2.0.CO;2). The basalts were erupted approximately 70 m.y. ago, an age equivalent to that of immediately adjacent oceanic crust in the Angola Basin and coraistent with formation at the paleo mid-ocean ridge (Moore et al., 1983). The basalt types vary from aphyric quartz tholeiites on the ridge crest to highly plagioclase phyric olivine tholeiites on the ridge flank. These show systematic differences in incompatible trace element and isotopic composition. Many element and isotope ratio pairs form systematic trends with the ridge crest basalts at one end and the highly phyric ridge flank basalts at the other. The low 143Nd/144Nd (0.51238), 206Pb/204Pb (17.54), 207Pb/204Pb (15.47), 208Pb/204Pb (38.14) and high 87Sr/86Sr (0.70512) ratios of the ridge crest basalts suggest derivation from an old Nd/Sm-, Rb/Sr- and Pb/U-enriched mantle source. This isotopic signature is similar to that of alkaline basalts on Tristan da Cunha but offset to significantly lower Nd and Pb isotopic ratios. The isotopic ratio trends may be extrapolated beyond the ridge flank basalts with higher 143Nd/144Nd (0.51270), 206Pb/204Pb (18.32), 207Pb/204Pb (15.52), 208Pb/204Pb (38.77) and lower 87Sr/86Sr (0.70417) ratios in the direction of increasingly Nd/Sm-, Rb/Sr- and Pb/U-depleted source compositions. These isotopic correlations are equally consistent with mixing of depleted and enriched end member melts or partial melting of an inhomogeneous, variably enriched mantle source. However, observed Zr-Ba-Nb-Y interelement relationships are inconsistent with any simple two-component model of magma mixing, as might result from the rise of a lower mantle plume through the upper mantle. Incompatible element and Pb isotopic systematics also preclude extensive involvement of depleted (N-type) MORB material or its mantle sources. In our preferred petrogenetic model the Walvis Ridge basalts were derived by partial melting of mantle similar to an enriched (E-type) MORB source which had become heterogeneous on a small scale due to the introduction of small-volume melts and metasomatic fluids.

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.

(Table 1) U and Th isotope analyses of pyrite and anhydrite of ODP Site 158-957

Hydrothermal circulation at oceanic spreading ridges causes sea water to penetrate to depths of 2 to 3 km in the oceanic crust where it is heated to ~400 °C before venting at spectacular 'black smokers'. These hydrothermal systems exert a strong influence on ocean chemistry (Edmond et al., 1979, doi:10.1016/0012-821X(79)90061-X), yet their structure, longevity and magnitude remain largely unresolved (Elderfield and Schultz., 1996, doi:10.1146/annurev.earth.24.1.191). The active Transatlantic Geotraverse (TAG) deposit, at 26° N on the Mid-Atlantic Ridge, is one of the largest, oldest and most intensively studied of the massive sulphide mounds that accumulate beneath black-smoker fields. Here we report ages of sulphides and anhydrites from the recently drilled (Humphris et al., 1995, doi:10.1038/377713a0) TAG substrate structures -determined from 234U-230Th systematics analysed by thermal ionization mass spectrometry. The new precise ages combined with existing data (Lalou et al., 1993, doi:10.1029/92JB01898; 1998, doi:10.2973/odp.proc.sr.158.214.1998) show that the oldest material (11,000 to 37,000 years old) forms a layer across the centre of the deposit with younger material (2,300-7,800 years old) both above and below. This stratigraphy confirms that much of the sulphide and anhydrite are precipitated within the mound by mixing of entrained sea water with hydrothermal fluid (James and Elderfield, 1996, doi:10.1130/0091-7613(1996)024<1147:COOFFA>2.3.CO;2). The age distribution is consistent with episodic activity of the hydrothermal system recurring at intervals of up to 2,000 years.

Chemical analyses of ODP Hole 105-647A basalts (Table 2)

Basalt samples recovered from the lowermost 37 m of Leg 105 Hole 647A in the Labrador Sea are fine- to medium grained, have microphenocrysts of clinopyroxene, and show little evidence of alteration. Chemically, these rocks are low potassium (0.01-0.09 wt% K20), olivine- to quartz-normative tholeiites that are also depleted in other incompatible elements. In terms of many of the incompatible trace elements, the Labrador Sea samples are similar both to iV-type midocean ridge basalts (MORBs) and to the terrestrial Paleocene volcanic rocks in the Davis Strait region of Baffin Island and West Greenland. However, significant differences are found in their strontium and neodymium isotope systematics. Hole 647A samples are more depleted in epsilon-Nd (+9.3) and are anomalously rich in 87Sr/86Sr (0.7040) relative to the Davis Strait basalts (epsilon-Nd +2.54 to + 8.97; mean 87Sr/86Sr, 0.7034). We conclude that the Hole 647A and Davis Strait basalts may have been derived from a similar depleted mantle source composition. In addition, the Davis Strait magmas were generated from mantle of more than one composition. We also suggest that there is no geochemical evidence from the Hole 647A samples to support or to refute the existence of foundered continental crust in the Labrador Sea.