Comparison of iron isotope variations in modern and Ordovician siliceous Fe oxyhydroxide deposits

Formation pathways of ancient siliceous iron formations and related Fe isotopic fractionation are still not completely understood. Investigating these processes, however, is difficult as good modern analogues to ancient iron formations are scarce. Modern siliceous Fe oxyhydroxide deposits are found at marine hydrothermal vent sites, where they precipitate from diffuse, low temperature fluids along faults and fissures on the seafloor. These deposits exhibit textural and chemical features that are similar to some Phanerozoic iron formations, raising the question as to whether the latter could have precipitated from diffuse hydrothermal fluids rather than from hydrothermal plumes. In this study, we present the first data on modern Fe oxyhydroxide deposits from the Jan Mayen hydrothermal vent fields, Norwegian-Greenland Sea. The samples we investigated exhibited very low δ56Fe values between -2.09‰ and -0.66‰. Due to various degrees of partial oxidation, the Fe oxyhydroxides are with one exception either indistinguishable from low-temperature hydrothermal fluids from which they precipitated (-1.84‰ and -1.53‰ in δ56Fe) or are enriched in the heavy Fe isotopes. In addition, we investigated Fe isotope variations in Ordovician jasper beds from the Løkken ophiolite complex, Norway, which have been interpreted to represent diagenetic products of siliceous ferrihydrite precursors that precipitated in a hydrothermal plume, in order to compare different formation pathways of Fe oxyhydroxide deposits. Iron isotopes in the jasper samples have higher δ56Fe values (-0.38‰ to +0.89‰) relative to modern, high-temperature hydrothermal vent fluids (ca. -0.40‰ on average), supporting the fallout model. However, formation of the Ordovician jaspers by diffuse venting cannot be excluded, due to lithological differences of the subsurface of the two investigated vent systems. Our study shows that reliable interpretation of Fe isotope variations in modern and ancient marine Fe oxyhydroxide deposits depends on comprehensive knowledge of the geological context. Furthermore, we demonstrate that very negative δ56Fe values in such samples might not be the result of microbial dissimilatory iron reduction, but could be caused instead by inorganic reactions.

Magmatic–hydrothermal molybdenum isotope fractionation and its relevance to the igneous crustal signature

We analysed the Mo isotope composition of a comprehensive series of molybdenite samples from the porphyry- type Questa deposit (NM, USA), as well as one rhyolite and one granite sample, directly associated with the Mo mineralization. The δ98Mo of the molybdenites ranges between −0.48‰ and +0.40‰, with a median at −0.05‰. The median Mo isotope composition increases from early magmatic (−0.29‰) to hydrothermal (−0.05‰) breccia mineralization (median bulk breccia = −0.17‰) to late stockwork veining (+0.22‰). Moreover, variations of up to 0.34‰ are found between different molybdenite crystals within an individual hand specimen. The rhyolite sample with 0.12 μg g−1 Mo has δ98Mo = −0.57‰ and is lighter than all molybde- nites from the Questa deposit, interpreted to represent the igneous leftover after aqueous ore fluid exsolution. We recognize three Mo isotope fractionation processes that occur between about 700 and 350 °C, affecting the Mo iso- tope composition of magmatic–hydrothermal molybdenites. Δ1Mo: Minerals preferentially incorporate light Mo isotopes during progressive fractional crystallization in subvolcanic magma reservoirs, leaving behind a melt enriched in heavy Mo isotopes. Δ2Mo: Magmatic–hydrothermal fluids preferentially incorporate heavy Mo iso- topes upon fluid exsolution. Δ3Mo: Light Mo isotopes get preferentially incorporated in molybdenite during crys- tallization from an aqueous fluid, leaving behind a hydrothermal fluid that gets heavier with progressive molybdenite crystallization. The sum of all three fractionation processes produces molybdenites that record heavier δ98Mo compositions than their source magmas. This implies that the mean δ98Mo of molybdenites published so far (~0.4‰) likely represents a maximum value for the Mo isotope composition of Phanerozoic igneous upper crust.

Exhumation rates in the Archean from pressure-ime paths: Examplefrom the Skjoldungen Orogen (SE Greenland)

The discussions on the orogenic evolution during Earth's history converge to the question of a different thermal structure in the Archean compared to the Phanerozoic and the applicability of the plate tectonic paradigm. However, geothermal structures are transient in orogens and are difficult to translate into large-scale tectonics and exhumation rates. Therefore, we propose depth–time data in the Archean Skjoldungen Orogen (SE Greenland, North Atlantic Craton) that allow for reconstruction of an exhumation rate independent of geothermal gradients. The resulting exhumation rate of ca. 0.4 km/Ma is similar to exhumation rates during erosion-controlled processes in modern orogens. These exhumation rates can only be established by erosion time constants similar to modern orogens. The occurrence of erosion-controlled exhumation is best explained by a stiff foreland promoting localized deformation in the orogen. Therefore, a switch from magmatic-dominated processes to localized deformation is proposed in the Skjoldungen Orogen area. This is supported by a change in magma composition and volume, from widespread granodiorite to localized alkaline intrusions. In addition, the involved metasedimentary rocks include detrital zircons of the only 50 Ma older foreland, which also correspond to erosion and tectonics as in modern orogens, i.e. flysh-type sediments. Relatively fast exhumation rates and the structural-magmatic evolution of the Neoarchean Skjoldungen Orogen thus indicate modern-style tectonic processes where stiff Mesoarchean continental crust forms a foreland to a collisional orogen instead of typical accretionary tectonics of weak island arc-like terranes in granite-greenstone terranes.

The Interplay between Melting, Refertilization and Carbonatite Metasomatism in Off-Cratonic Lithospheric Mantle under Zealandia: an Integrated Major, Trace and Platinum Group Element Study

It is widely accepted that stabilization of the continental crust requires the presence of sub-continental lithospheric mantle. However, the degree of melt depletion required to stabilize the lithosphere and whether widespread refertilization is a significant process remain unresolved. Here, major and trace element, including platinum group elements (PGE), characterization of 40 mantle xenoliths from 13 localities is used to constrain the melt depletion, refertilization and metasomatic history of lithospheric mantle underneath the micro-continent Zealandia. Our previously published Re–Os isotopic data for a subset of these xenoliths indicate Phanerozoic to Paleoproterozoic ages and, reinterpreted with the new major and trace element data presented here, demonstrate that a large volume (>2 million km3) of lithospheric mantle with an age of 1·99 ± 0·21 Ga is present below the much younger crust of Zealandia. A peritectic melting model using moderately incompatible trace elements (e.g. Yb) in bulk-rocks demonstrates that these peridotites experienced a significant range of degrees of partial melting, between 3 and 28%. During subsolidus equilibration clinopyroxene gains significant rare earth elements (REE), which then leads to the underestimation of the degree of partial melting by ≤12% in fertile xenoliths. A new approach taking into account the effects of subsolidus re-equilibration on clinopyroxene composition effectively removes discrepancies in the calculated degree of melting and provides consistent estimates of between 4 and 29%. The estimated amount of melting is independent of the Re–Os model ages of the samples. The PGE patterns record simple melt depletion histories and the retention of primary base metal sulfides in the majority of the xenoliths. A rapid decrease in Pt/IrN observed at c. 1·0 wt % Al2O3 is a direct result of the exhaustion of sulfide in the mantle residue at c. 20–25% partial melting and the inability of Pt to form a stable alloy phase. Major elements preserve evidence for refertilization by a basaltic component that resulted in the formation of secondary clinopyroxene and low-forsterite olivine. The majority of xenoliths show the effects of cryptic metasomatic overprinting, ranging from minor to strong light REE enrichments in bulk-rocks (La/YbN = 0·16–15·9). Metasomatism is heterogeneous, with samples varying from those with weak REE enrichment and notable positive Sr and U–Th anomalies and negative Nb–Ta anomalies in clinopyroxene to those that have extremely high concentrations of REE, Th–U and Nb. Chemical compositions are consistent with a carbonatitic component contributing to the metasomatism of the lithosphere under Zealandia. Notably, the intense metasomatism of the samples did not affect the PGE budget of the peridotites as this was controlled by residual sulfides.

Cretaceous stable isotopic record and sea surface temperture estimation for ODP Site 207-1259

The Turonian (93.5 to 89.3 million years ago) was one of the warmest periods of the Phanerozoic eon, with tropical sea surface temperatures over 35°C. High-amplitude sea-level changes and positive d18O excursions in marine limestones suggest that glaciation events may have punctuated this episode of extreme warmth. New d18O data from the tropical Atlantic show synchronous shifts ~91.2 million years ago for both the surface and deep ocean that are consistent with an approximately 200,000-year period of glaciation, with ice sheets of about half the size of the modern Antarctic ice cap. Even the prevailing supergreenhouse climate was not a barrier to the formation of large ice sheets, calling into question the common assumption that the poles were always ice-free during past periods of intense global warming.

Biostratigraphy and Sr isotopic composition of sediments from DSDP Legs 12, 14, 15, 17, 20-22, and 30

We present the data used to construct the Cenozoic and Cretaceous portion of the Phanerozoic curve of seawater 87Sr/86Sr that had been given in summary form by W.H. Burke co-workers. All Cenozoic samples (128) and 22 Cretaceous samples are foram-nannofossil oozes and limestones from DSDP cores distributed among 13 sites in the Atlantic, Pacific and Indian Oceans, and the Caribbean Sea. Non-DSDP Cretaceous samples (126) include limestone, anhydrite and phosphate samples from North America, Europe and Asia. Determination of the 87Sr/86Sr value of seawater at particular times in the past is based on comparison of ratios derived from coeval marine samples from widely separated geographic areas. These samples are characterized by a wide variety of diagenetic and burial histories. The large size and cosmopolitan nature of the data set decreases the likelihood that, among coeval data, systematic error has been introduced by a similar pattern of diagenetic alteration of the ratios. There is good clustering of data points throughout the Cenozoic and Cretaceous curve. The consistency of data is illustrated by Cenozoic and Cretaceous data plots that include a separate symbol for each DSDP site and non-DSDP sample location. More than 98% of the data points are enclosed by upper and lower lines that define a narrow band. For any given time, the correct seawater ratio probably lies within this band. A line drawn within the band represents our estimate of the actual seawater ratio as a function of time. The general configuration of the Cenozoic and Cretaceous curve appears to be strongly influenced by the history of plate interactions and sea-floor spreading. Specific rises and falls in the 87Sr/86Sr of seawater, however, may be caused by a variety of factors such as variation in lithologic composition of the crust exposed to weathering, configuration and topographic relief of continents, volcanic activity, rate of sea-floor spreading, extent of continental inundation by epeiric seas, and variations in both climate and paleooceanographic conditions. Many or all of these factors are probably related to global tectonic processes, yet their combined effect on the temporal variation of seawater 87Sr/86Sr can complicate a direct platetectonic interpretation for portions of the seawater curve.

Trace element contents in zircons from oceanic crust

We present newly acquired trace element compositions for more than 300 zircon grains in 36 gabbros formed at the slow-spreading Mid-Atlantic and Southwest Indian Ridges. Rare earth element patterns for zircon from modern oceanic crust completely overlap with those for zircon crystallized in continental granitoids. However, plots of U versus Yb and U/Yb versus Hf or Y discriminate zircons crystallized in oceanic crust from continental zircon, and provide a relatively robust method for distinguishing zircons from these environments. Approximately 80% of the modern ocean crust zircons are distinct from the field defined by more than 1700 continental zircons from Archean and Phanerozoic samples. These discrimination diagrams provide a new tool for fingerprinting ocean crust zircons derived from reservoirs like that of modern mid-ocean ridge basalt (MORB) in both modern and ancient detrital zircon populations. Hadean detrital zircons previously reported from the Acasta Gneiss, Canada, and the Narryer Gneiss terrane, Western Australia, plot in the continental granitoid field, supporting hypotheses that at least some Hadean detrital zircons crystallized in continental crust forming magmas and not from a reservoir like modern MORB.

Compositions of paleosols from the Salarevo Formation, Sukhona River valley (Severnaya Dvina basin)

Paleosols crop out in the Sukhona River valley as several members up to 10 m thick embedded into the Salarevo Formation sediments. Principal characteristics of the paleosols include a dense network of root channels, indications of eluvial gley alteration, redistribution and formation of secondary carbonates represented by several generations, and formation of block-prismatic soil structure with specific clayey films at structural jointing faces. The paleosols are divided into a number of genetically interrelated horizons (from top to bottom): presumably organogenic accumulation (AElg), eluvial gley horizon (Elg), illuvial horizons (B1 and B2), illuvial gley horizon (Bg), and transitional horizons (ElBg and BElg). The paleosols formed under conditions of a semiarid climate with sharp seasonal or secular and multisecular oscillations of atmospheric precipitation. Such soils point to specific ecological environments existed in the northern semiarid belt of the Earth before the greatest (in Phanerozoic) biospheric crisis at the Permian-Triassic boundary.

Accumulation and composition of terrigenous material on Ceara rise over the past 55 Ma

Ceara Rise, located east the Amazon River mouth, is covered with a thick blanket of pelagic carbonate and hemipelagic terrigenous sediment. The terrigenous component has been extracted from 57 bulk sediment samples at Ocean Drilling Program (ODP) Sites 925 and 929 on Ceara Rise to obtain a Cenozoic record of riverine discharge from northern South America. From the early Eocene to early Miocene (55-20 Ma), terrigenous accumulation was dominated by moderate amounts of generally large-grained, gray to green sediment especially depleted in elements that are enriched in post-Archaean shale (e.g. Cs, Th, Yb). However, pulsed inputs of relatively small-grained, gray to green terrigenous sediment less depleted in the above elements occurred in the late Eocene and Oligocene. The accumulation of terrigenous sediment decreased significantly until 16.5 Ma. In the middle Miocene (16.5-13 Ma), terrigenous accumulation was dominated by small amounts of small-grained, tan sediment notably depleted in Na and heavy rare earth elements. The accumulation rate of terrigenous sediment increased markedly from the latest Miocene (10 Ma) to the present day, a change characterized by deposition of gray-green sediment enriched in elements that are enriched in post-Archaean shale. Observed changes in terrigenous sediment at Ceara Rise record tectonism and erosion in northern South America. The Brazil and Guyana shields supplied sediment to the eastern South American margin until the middle Miocene (20-16.5 Ma) when a period of thrusting, shortening and uplift changed the source region, probably first to highly weathered and proximal Phanerozoic sediments. By the late Miocene (9 Ma), there was a transcontinental connection between the Andes and eastern South America. Weathering products derived from the Andes have increasingly dominated terrigenous deposition at Ceara Rise since the Late Miocene and especially since the late Pliocene.

Seawater carbonate chemistry and carbon and oxygen isotopes during experiments with planktonic foraminifera Orbulina universa and Globigerina bulloides, 1997

Stable oxygen and carbon isotope measurements on biogenic calcite and aragonite have become standard tools for reconstructing past oceanographic and climatic change. In aquatic organisms, 18O/16O ratios in the shell carbonate are a function of the ratio in the sea water and the calcification temperature (Epstein et al., 1953). In contrast, 13C/12C ratios are controlled by the ratio of dissolved inorganic carbon in sea water and physiological processes such as respiration and symbiont photosynthesis (Spero et al., 1991, doi:10.1029/91PA02022). These geochemical proxies have been used with analyses of foraminifera shells to reconstruct global ice volumes (Shackleton and Opdyke, 1973, doi:10.1016/0033-5894(73)90052-5), surface and deep ocean temperatures (Broecker, 1986, doi:10.1016/0033-5894(86)90087-6; Labeyrie et al., 1987, doi:10.1038/327477a0), ocean circulation changes (Duplessy et al., 1988, doi:10.1029/PA003i003p00343) and glacial-interglacial exchange between the terrestrial and oceanic carbon pools (Sackleton, 1977). Here, we report experimental measurements on living symbiotic and non-symbiotic plankton foraminifera (Orbulina universa and Globigerina bulloides respectively) showing that the 13C/12C and 18O/16O ratios of the calcite shells decrease with increasing seawater [CO3 2-]. Because glacial-period oceans had higher pH and [CO3 2-] than today (Sanyal et al., 1995, doi:10.1038/373234a0), these new relationships confound the standard interpretation of glacial foraminiferal stable-isotope data. In particular, the hypothesis that the glacial-interglacial shift in the 13C/12C ratio was due to a transfer of terrestrial carbon into the ocean(Shackleton ,1977) can be explained alternatively by an increase in ocean alkalinity (Lea et al., 1996). A carbonate-concentration effect could also help explain some of the extreme stable-isotope variations during the Proterozoic and Phanerozoic aeons (Kaufman et al., 1993, doi:10.1016/0012-821X(93)90254-7).

An insight into the breakup of Gondwana: Identifying events through low-temperature thermochronology from the basement rocks of Madagascar

Fission track analysis was applied to the Precambrian suites of Madagascar in order to identify the lower-temperature cooling histories and their relationships to the Phanerozoic events that affected the island. Apatite ages range from 431 to 68 Ma, and zircon ages range from 452 to 238 Ma. Thermochronologically, the island can be divided into a southern, central, and northern region each with a subdivision on an east-west basis. The southern region is sharply separated from the central region by strongly contrasting apparent apatite ages over the northwest-southeast striking Ranotsara Shear Zone (RSZ). The change in apparent ages over the RSZ is indicative of later reactivation along younger brittle faults. The central region has the oldest ages of the island and has a diffuse contact to the third region northward. Along the entire western margin of the Precambrian basement initial Paleozoic exhumation was followed by heating (burial by sediments) during Jurassic and Cretaceous times. A decrease in ages along the eastern margin from 119 to 68 Ma coincides with the predicted positions of the Marion hot spot after effects of erosion are considered. On the other hand, these ages may represent progressive opening of the margin in a southward direction together with associated denudation of the rift shoulder. The eastern part of the central region has remained very stable since at least Devonian times, undergoing only long-term very slow exhumation at rates of 1–5 m/Myr.

(Table S1) Boron isotope and trace element data and reconstructed carbonate system parameters, ODP Holes 122-761B and 154-926A

The middle Miocene Climatic Optimum (17-15 Ma; MCO) is a period of global warmth and relatively high CO2 and is thought to be associated with a significant retreat of the Antarctic Ice Sheet (AIS). We present here a new planktic foraminiferal d11B record from 16.6 to 11.8 Ma from two deep ocean sites currently in equilibrium with the atmosphere with respect to CO2. These new data demonstrate that the evolution of global climate during the middle Miocene (as reflected by changes in the cyrosphere) was well correlated to variations in the concentration of atmospheric CO2. What is more, within our sampling resolution (~1 sample per 300 kyr) there is no evidence of hysteresis in the response of ice volume to CO2 forcing during the middle Miocene, contrary to what is understood about the Antarctic Ice Sheet from ice sheet modelling studies. In agreement with previous data, we show that absolute levels of CO2 during the MCO were relatively modest (350-400 ppm) and levels either side of the MCO are similar or lower than the pre-industrial (200-260 ppm). These new data imply the presence of either a very dynamic AIS at relatively low CO2 during the middle Miocene or the advance and retreat of significant northern hemisphere ice. Recent drilling on the Antarctic margin and shore based studies indicate significant retreat and advance beyond the modern limits of the AIS did occur during the middle Miocene, but the complete loss of the AIS was unlikely. Consequently, it seems that ice volume and climate variations during the middle Miocene probably involved a more dynamic AIS than the modern but also some component of land-based ice in the northern hemisphere.

Geological Time Scale 2012. Chapter 10. Oxygen isotope stratigraphy, Appendix 1

Variations in the 18O/16O ratios of marine fossils and microfossils record changes in seawater 18O/16O and temperature and form the basis for global correlation. Relying on previous compilations and new data, this chapter presents oxygen isotope curves for Phanerozoic foraminifera, mollusks, brachiopods, and conodonts, and for Precambrian limestones, dolostones, and cherts. Periodic oxygen-isotopic variations in deep-sea foraminifera define marine isotope stages that, when combined with biostratigraphy and astronomical tuning, provide a late Cenozoic chronostratigraphy with a resolution of several thousand years. Oxygen isotope events of early Cenozoic, Mesozoic, and Paleozoic age serve as chemostratigraphic markers for regional and global correlation. Precambrian oxygen isotope stratigraphy, however, is hampered by the lack of unaltered authigenic marine sediments.

Susceptibility measurements at Gubbio sections

The complex interplay between extraterrestrial events and earth-bound processes that triggered one of the greatest biological crises of the Phanerozoic requires a high resolution timescale. Detailed magnetic susceptibility measurements at the Contessa Highway and Bottaccione sections (Italy) span the Cretaceous-Paleogene boundary and reveal clear orbital signatures in the sedimentary record. Identification of precession and 405 kyr eccentricity cycles allows an estimate of 324+/-40 kyr for the duration of the Maastrichtian part of Chron C29r. We present in the same high resolution time frame sites in Spain and the North and South Atlantic and bio-horizons, biotic changes, stable isotopic excursions and the decrease in Osmium isotopes recorded in these sections. The onset of 187 Os/ 188 Os decrease coincides with the d13 C negative excursion K-PgE1, thus suggesting a first pulse in Deccan volcanism at 66.64 Ma. The K-PgE3 d13 C negative excursion is possibly the expression of a second pulse at 66.26 Ma. Late Maastrichtian d13 C negative excursions are of low intensity and span durations of one to two eccentricity cycles, whereas early Danian excursions are brief (about 30 kyr) and acute. In Biotic response to late Maastrichtian perturbations occurred with a delay of ca. 200 kyr after the beginning of K-PgE1 shortly before K-PgE3. The biotic perturbation could be thus either a delayed response to K-PgE1, or a direct response to K-PgE3, and possibly, a threshold response to the stepwise buildup of CO2 atmospheric injections. No delay is evident in response to early Danian hyperthermal events. These differences suggest that short-lived, volcanically-derived environmental perturbations were buffered within the stable late Maastrichtian oceanic realm whereas they were amplified by the more sensitive and highly disturbed early Danian oceanic ecosystem.