Nitrogen isotopic composition and chemical composition of altered oceanic crust from ODP Legs 129 and 185

Knowledge of the subduction input flux of nitrogen (N) in altered oceanic crust (AOC) is critical in any attempt to mass-balance N across arc-trench systems on a global or individual-margin basis. We have employed sealed-tube, carrier-gas-based methods to examine the N concentrations and isotopic compositions of AOC. Analyses of 53 AOC samples recovered on DSDP/ODP legs from the North and South Pacific, the North Atlantic, and the Antarctic oceans (with larger numbers of samples from Site 801 outboard of the Mariana trench and Site 1149 outboard of the Izu trench), and 14 composites for the AOC sections at Site 801, give N concentrations of 1.3 to 18.2 ppm and d15N_air of -11.6? to +8.3?, indicating significant N enrichment probably during the early stages of hydrothermal alteration of the oceanic basalts. The N-d15N modeling for samples from Sites 801 and 1149 (n=39) shows that the secondary N may come from (1) the sedimentary N in the intercalated sediments and possibly overlying sediments via fluid-sediment/rock interaction, and (2) degassed mantle N2 in seawater via alteration-related abiotic reduction processes. For all Site 801 samples, weak correlation of N and K2O contents indicates that the siting of N in potassic alteration phases strongly depends on N availability and is possibly influenced by highly heterogeneous temperature and redox conditions during hydrothermal alteration. The upper 470-m AOC recovered by ODP Legs 129 and 185 delivers approximately 800 kg/km N annually into the Mariana margin. If the remaining less-altered oceanic crust (assuming 6.5 km, mostly dikes and gabbros) has MORB-like N of 1.5 ppm, the entire oceanic crust transfers 5100 kg/km N annually into that trench. This N input flux is twice as large as the annual N input of 2500 kg/km in seafloor sediments subducting into the same margin, demonstrating that the N input in oceanic crust, and its isotopic consequences, must be considered in any assessment of convergent margin N flux.

Noble gas concentrations of oceanic crust and sediments from the Ninety East Ridge in the Indian Ocean

We have determined the concentrations and isotopic composition of noble gases in old oceanic crust and oceanic sediments and the isotopic composition of noble gases in emanations from subduction volcanoes. Comparison with the noble gas signature of the upper mantle and a simple model allow us to conclude that at least 98% of the noble gases and water in the subducted slab returns back into the atmosphere through subduction volcanism before they can be admixed into the earth's mantle. It seems that the upper mantle is inaccessible to atmospheric noble gases due to an efficient subduction barrier for volatiles.

Geochemistry, and Nd- and Pb-isotopic composition of basalts from the upper oceanic crust of the Costa Rica Rift

The DSDP/ODP Hole 504B, drilled in the 5.9 Ma southern flank of the Costa Rica Rift, represents the deepest section through modern ocean floor basaltic basement. The hole penetrates a 570 m thick volcanic zone, a 210 m thick transition zone of volcanic rocks and dykes, and 1056 m of dykes. A representative selection of these basalt types has been investigated with respect to Nd and Pb isotopes. The epsilonNd of the basalts varies from 7.62 to 11.16. This range in the Nd-isotope composition represents about 67% of the total range reported for Pacific MORB. The Pb-isotope composition also shows significant variation, with 206Pb/204Pb varying from 17.90 to 18.82. The isotopic data show that a small volume of enriched mantle existed in the source. The large ranges in isotopic composition in a single drill hole demonstrate the importance of small-scale mantle heterogeneities in the petrogenesis of MORB. Fractional melting and extraction of small magma batches by channelled flow, and small, short-lived crustal magma reservoirs, with limited potential for mixing of the mantle derived magmas, are favored by these isotopic data.

(Table 2) Chlorine geochemistry and mineralogy of altered oceanic crust samples

Bulk chlorine concentrations and chlorine stable isotope compositions were determined for hydrothermally altered basalt (extrusive lavas and sheeted dikes) and gabbro samples (n = 50) from seven DSDP/ODP/IODP drill sites. These altered oceanic crust (AOC) samples span a range of crustal ages, tectonic settings, alteration type, and crustal depth. Bulk chlorine concentrations range from < 0.01 wt.% to 0.09 wt.%. In general, higher chlorine concentrations coincide with an increase in temperature of alteration and amphibole content. d37Cl values of whole rock AOC samples range from -1.4 to +1.8 per mil. High d37Cl values (>=~0.5 per mil) are associated with areas of higher amphibole content. This observation is consistent with theoretical calculations that estimate amphibole should be enriched in 37Cl compared to co-existing fluid. Negative to near zero d37Cl values are found in areas dominated by clay minerals. Chlorine geochemistry is a rough indicator of metamorphic grade and mineralogy. AOC is a major Cl host in the subducting oceanic lithospheric slab. Here we show that bulk chlorine concentrations are ~3 times higher than previous estimates resulting in a greater contribution of Cl to the mantle.

Thinning of the Goban Spur continental margin and formation of early oceanic crust: Constraints from forward modelling and inversion of marine magnetic anomalies

The deep seismic reflection profile Western Approaches Margin (WAM) cuts across the Goban Spur continental margin, located southwest of Ireland. This non-volcanic margin is characterized by a few tilted blocks parallel to the margin. A volcanic sill has been emplaced on the westernmost tilted block. The shape of the eastern part of this sill is known from seismic data, but neither seismic nor gravity data allow a precise determination of the extent and shape of the volcanic body at depth. Forward modelling and inversion of magnetic data constrain the shape of this volcanic sill and the location of the ocean-continent transition. The volcanic body thickens towards the ocean, and seems to be in direct contact with the oceanic crust. In the contact zone, the volcanic body and the oceanic magnetic layer display approximately the same thickness. The oceanic magnetic layer is anomalously thick immediately west of the volcanic body, and gradually thins to reach more typical values 40 km further to the west. The volcanic sill would therefore represent the very first formation of oceanic crust, just before or at the continental break-up. The ocean-continent transition is limited to a zone 15 km wide. The continental magnetic layer seems to thin gradually oceanwards, as does the continental crust, but no simple relation is observed between their respective thinnings.

Global Rates of Free Hydrogen (H2) Production by Serpentinization and other Abiogenic Processes within Young Ocean Crust

The main conclusion of this dissertation is that global H2 production within young ocean crust (<10 Mya) is higher than currently recognized, in part because current estimates of H2 production accompanying the serpentinization of peridotite may be too low (Chapter 2) and in part because a number of abiogenic H2-producing processes have heretofore gone unquantified (Chapter 3). The importance of free H2 to a range of geochemical processes makes the quantitative understanding of H2 production advanced in this dissertation pertinent to an array of open research questions across the geosciences (e.g. the origin and evolution of life and the oxidation of the Earth’s atmosphere and oceans).

The first component of this dissertation (Chapter 2) examines H2 produced within young ocean crust [e.g. near the mid-ocean ridge (MOR)] by serpentinization. In the presence of water, olivine-rich rocks (peridotites) undergo serpentinization (hydration) at temperatures of up to ~500°C but only produce H2 at temperatures up to ~350°C. A simple analytical model is presented that mechanistically ties the process to seafloor spreading and explicitly accounts for the importance of temperature in H2 formation. The model suggests that H2 production increases with the rate of seafloor spreading and the net thickness of serpentinized peridotite (S-P) in a column of lithosphere. The model is applied globally to the MOR using conservative estimates for the net thickness of lithospheric S-P, our least certain model input. Despite the large uncertainties surrounding the amount of serpentinized peridotite within oceanic crust, conservative model parameters suggest a magnitude of H2 production (~1012 moles H2/y) that is larger than the most widely cited previous estimates (~1011 although previous estimates range from 1010-1012 moles H2/y). Certain model relationships are also consistent with what has been established through field studies, for example that the highest H2 fluxes (moles H2/km2 seafloor) are produced near slower-spreading ridges (<20 mm/y). Other modeled relationships are new and represent testable predictions. Principal among these is that about half of the H2 produced globally is produced off-axis beneath faster-spreading seafloor (>20 mm/y), a region where only one measurement of H2 has been made thus far and is ripe for future investigation.

In the second part of this dissertation (Chapter 3), I construct the first budget for free H2 in young ocean crust that quantifies and compares all currently recognized H2 sources and H2 sinks. First global estimates of budget components are proposed in instances where previous estimate(s) could not be located provided that the literature on that specific budget component was not too sparse to do so. Results suggest that the nine known H2 sources, listed in order of quantitative importance, are: Crystallization (6x1012 moles H2/y or 61% of total H2 production), serpentinization (2x1012 moles H2/y or 21%), magmatic degassing (7x1011 moles H2/y or 7%), lava-seawater interaction (5x1011 moles H2/y or 5%), low-temperature alteration of basalt (5x1011 moles H2/y or 5%), high-temperature alteration of basalt (3x1010 moles H2/y or <1%), catalysis (3x108 moles H2/y or <<1%), radiolysis (2x108 moles H2/y or <<1%), and pyrite formation (3x106 moles H2/y or <<1%). Next we consider two well-known H2 sinks, H2 lost to the ocean and H2 occluded within rock minerals, and our analysis suggests that both are of similar size (both are 6x1011 moles H2/y). Budgeting results suggest a large difference between H2 sources (total production = 1x1013 moles H2/y) and H2 sinks (total losses = 1x1011 moles H2/y). Assuming this large difference represents H2 consumed by microbes (total consumption = 9x1011 moles H2/y), we explore rates of primary production by the chemosynthetic, sub-seafloor biosphere. Although the numbers presented require further examination and future modifications, the analysis suggests that the sub-seafloor H2 budget is similar to the sub-seafloor CH4 budget in the sense that globally significant quantities of both of these reduced gases are produced beneath the seafloor but never escape the seafloor due to microbial consumption.

The third and final component of this dissertation (Chapter 4) explores the self-organization of barchan sand dune fields. In nature, barchan dunes typically exist as members of larger dune fields that display striking, enigmatic structures that cannot be readily explained by examining the dynamics at the scale of single dunes, or by appealing to patterns in external forcing. To explore the possibility that observed structures emerge spontaneously as a collective result of many dunes interacting with each other, we built a numerical model that treats barchans as discrete entities that interact with one another according to simplified rules derived from theoretical and numerical work, and from field observations: Dunes exchange sand through the fluxes that leak from the downwind side of each dune and are captured on their upstream sides; when dunes become sufficiently large, small dunes are born on their downwind sides (“calving”); and when dunes collide directly enough, they merge. Results show that these relatively simple interactions provide potential explanations for a range of field-scale phenomena including isolated patches of dunes and heterogeneous arrangements of similarly sized dunes in denser fields. The results also suggest that (1) dune field characteristics depend on the sand flux fed into the upwind boundary, although (2) moving downwind, the system approaches a common attracting state in which the memory of the upwind conditions vanishes. This work supports the hypothesis that calving exerts a first order control on field-scale phenomena; it prevents individual dunes from growing without bound, as single-dune analyses suggest, and allows the formation of roughly realistic, persistent dune field patterns.

Magma Evolution in the Primitive, Intra-oceanic Tonga Arc: Petrogenesis of Basaltic Andesites at Tofua Volcano

Tofua volcano is situated midway along the Tonga oceanic arc and has undergone two phases of ignimbrite-forming activity. The eruptive products are almost entirely basaltic andesites (52 center dot 5-57 wt % SiO2) with the exception of a volumetrically minor pre-caldera dacite. The suite displays a strong tholeiitic trend with K2O <1 wt %. Phenocryst assemblages typically comprise plagioclase + clinopyroxene +/- orthopyroxene with microlites of Ti-magnetite. Olivine (Fo(83-88)) is rare and believed to be dominantly antecrystic. An increase in the extent and frequency of reverse zoning in phenocrysts, sieve-textured plagioclase and the occurrence of antecrystic phases in post-caldera lavas record a shift to dynamic conditions, allowing the interaction of magma batches that were previously distinct. Pyroxene thermobarometry suggests crystallization at 950-1200 degrees C and 0 center dot 8-1 center dot 8 kbar. Volatile measurements of glassy melt inclusions indicate a maximum H2O content of 4 center dot 16 wt % H2O, and CO2-H2O saturation curves indicate that crystallization occurred at two levels, at depths of 4-5 center dot 5 km and 1 center dot 5-2 center dot 5 km. Major and trace element models suggest that the compositions of the majority of the samples represent a differentiation trend whereby the dacite was produced by 65% fractional crystallization of the most primitive basaltic andesite. Trace element models suggest that the sub-arc mantle source is the residuum of depleted Indian mid-ocean ridge basalt mantle (IDMM-1% melt), whereas radiogenic isotope data imply addition of 0 center dot 2% average Tongan sediment melt and a fluid component derived from the subducted altered Pacific oceanic crust. A horizontal array on the U-Th equiline diagram and Ra excesses of up to 500% suggest fluid addition to the mantle wedge within the last few thousand years. Time-integrated (Ra-226/Th-230) vs Sr/Th and Ba/Th fractionation models imply differentiation timescales of up to 4500 years for the dacitic magma compositions at Tofua.

Deciphering petrological signatures of reactive melt stagnation and cooling in the oceanic mantle underneath ultraslow-spreading ridges

The global mid-ocean ridge system creates oceanic crust and lithosphere that covers more than two-thirds of the Earth. Basalts are volumetrically the most important rock type sampled at mid-ocean ridges. For this reason, our present understanding of upper mantle dynamics and the chemical evolution of the earth is strongly influenced by the study of mid-ocean ridge basalts (MORB). However, MORB are aggregates of polybarically generated small melt increments that can undergo a variety of physical and chemical processes during their ascent and consequently affect their derivative geochemical composition. Therefore, MORB do not represent “direct” windows to the underlying upper mantle. Abyssal peridotites, upper mantle rocks recovered from the ocean floor, are the residual complement to MORB melting and provide essential information on melt extraction from the upper mantle. In this study, abyssal peridotites are examined to address these overarching questions posed by previous studies of MORB: How are basaltic melts formed in the mantle, how are they extracted from the mantle and what physical and chemical processes control mantle melting? The number of studies on abyssal peridotites is small compared to those on basalts, in part because seafloor exposures of abyssal peridotites are relatively rare. For this reason, abyssal peridotite characteristics need to be considered in the context of subaerially exposed peridotites associated with ophiolites, orogenic peridotite bodies and basalt-hosted xenoliths. However, orogenic peridotite bodies are mainly associated with passive continental margins, most ophiolites are formed in supra-subduction zone settings, and peridotite xenoliths are often contaminated by their host magma. Therefore, studies of abyssal peridotites are essential to understanding the primary characteristics of the oceanic upper mantle free from the influence of continental rifting, subduction and tectonic emplacement processes. Nevertheless, numerous processes such as melt stagnation and cooling-induced, inter-mineral exchange can affect residual abyssal peridotite compositions after the cessation of melting. The aim of this study is to address these post-melting modifications of abyssal peridotites from a petrological-geochemical perspective. The samples in this study were dredged along the axis of the ultraslow-spreading Gakkel Ridge in the Arctic Ocean within the “Sparsely Magmatic Zone”, a 100 km ridge section where only mantle rocks are exposed. During two expeditions (ARK XVII-2 in 2001 and ARK XX-2 in 2004), exceptionally fresh peridotites were recovered. The boulders and cobbles collected cover a range of mantle rock compositions, with most characterized as plagioclase-free spinel peridotites or plagioclase- spinel peridotites. This thesis investigates melt stagnation and cooling processes in the upper mantle and is divided into two parts. The first part focuses on processes in the stability field of spinel peridotites (>10 kb) such as melt refertilization and cooling related trace element exchange, while the second part investigates processes in the stability field of plagioclase peridotites (< 10 kb) such as reactive melt migration and melt stagnation. The dissertation chapters are organized to follow the theoretical ascent of a mantle parcel upwelling beneath the location where the samples were collected.

Forward modelling of petrological crust-forming processes on the early Earth

Tonalite-trondhjemite-granodiorite (TTG) gneisses form up to two-thirds of the preserved Archean continental crust and there is considerable debate regarding the primary magmatic processes of the generation of these rocks. The popular theories indicate that these rocks were formed by partial melting of basaltic oceanic crust which was previously metamorphosed to garnet-amphibolite and/or eclogite facies conditions either at the base of thick oceanic crust or by subduction processes.rnThis study investigates a new aspect regarding the source rock for Archean continental crust which is inferred to have had a bulk compostion richer in magnesium (picrite) than present-day basaltic oceanic crust. This difference is supposed to originate from a higher geothermal gradient in the early Archean which may have induced higher degrees of partial melting in the mantle, which resulted in a thicker and more magnesian oceanic crust. rnThe methods used to investigate the role of a more MgO-rich source rock in the formation of TTG-like melts in the context of this new approach are mineral equilibria calculations with the software THERMOCALC and high-pressure experiments conducted from 10–20 kbar and 900–1100 °C, both combined in a forward modelling approach. Initially, P–T pseudosections for natural rock compositions with increasing MgO contents were calculated in the system NCFMASHTO (Na2O–CaO–FeO–MgO–Al2O3–SiO2–H2O–TiO2) to ascertain the metamorphic products from rocks with increasing MgO contents from a MORB up to a komatiite. A small number of previous experiments on komatiites showed the development of pyroxenite instead of eclogite and garnet-amphibolite during metamorphism and established that melts of these pyroxenites are of basaltic composition, thus again building oceanic crust instead of continental crust.rnThe P–T pseudosections calculated represent a continuous development of their metamorphic products from amphibolites and eclogites towards pyroxenites. On the basis of these calculations and the changes within the range of compositions, three picritic Models of Archean Oceanic Crust (MAOC) were established with different MgO contents (11, 13 and 15 wt%) ranging between basalt and komatiite. The thermodynamic modelling for MAOC 11, 13 and 15 at supersolidus conditions is imprecise since no appropriate melt model for metabasic rocks is currently available and the melt model for metapelitic rocks resulted in unsatisfactory calculations. The partially molten region is therfore covered by high-pressure experiments. The results of the experiments show a transition from predominantly tonalitic melts in MAOC 11 to basaltic melts in MAOC 15 and a solidus moving towards higher temperatures with increasing magnesium in the bulk composition. Tonalitic melts were generated in MAOC 11 and 13 at pressures up to 12.5 kbar in the presence of garnet, clinopyroxene, plagioclase plus/minus quartz (plus/minus orthopyroxene in the presence of quartz and at lower pressures) in the absence of amphibole but it could not be explicitly indicated whether the tonalitic melts coexisting with an eclogitic residue and rutile at 20 kbar do belong to the Archean TTG suite. Basaltic melts were generated predominantly in the presence of granulite facies residues such as amphibole plus/minus garnet, plagioclase, orthopyroxene that lack quartz in all MAOC compositions at pressures up to 15 kbar. rnThe tonalitic melts generated in MAOC 11 and 13 indicate that thicker oceanic crust with more magnesium than that of a modern basalt is also a viable source for the generation of TTG-like melts and therefore continental crust in the Archean. The experimental results are related to different geologic settings as a function of pressure. The favoured setting for the generation of early TTG-like melts at 15 kbar is the base of an oceanic crust thicker than existing today or by melting of slabs in shallow subduction zones, both without interaction of tonalic melts with the mantle. Tonalitic melts at 20 kbar may have been generated below the plagioclase stability by slab melting in deeper subduction zones that have developed with time during the progressive cooling of the Earth, but it is unlikely that those melts reached lower pressure levels without further mantle interaction.rn

Alkali contents, Sr isotope and U-Th-Pb data for ocean crust vein minerals from DSDP Holes 27-261, 61-462A, and 72-516F

Cation exchange experiments (ammonium acetate and cation resin) on celadonite-smectite vein minerals from three DSDP holes demonstrate selective removal of common Sr relative to Rb and radiogenic Sr. This technique increases the Rb/Sr ratio by factors of 2.3 to 22 without significantly altering the age of the minerals, allowing easier and more precise dating of such vein minerals. The ages determined by this technique (Site 261 - 121.4+/-1.6 m.y.; Site 462A - 105.1+/-2.8 m.y.; Site 516F - 69.9+/-2.4 m.y.) are 34, 54 and 18 m.y. younger, respectively, than the age of crust formation at the site; in the case of site 462A, the young age is clearly related to off-ridge emplacement of a massive sill/flow complex. At the other sites, either the hydrothermal circulation systems persisted longer than for normal crust (10-15 m.y.), or were reactivated by off-ridge igneous activity. Celadonites show U and Pb contents and Pb isotopic compositions little changed from their basalt precursors, while Th contents are significantly lower. Celadonites thus have unusually high alkali/U,Th ratios and low Th/U ratios. If this celadonite alteration signature is significantly imprinted on oceanic crust as a whole, it will lead to very distinctive Pb isotope signatures for any hot spot magmas which contain a component of aged subducted recycled oceanic crust. Initial Sr isotope ratios of ocean crust vein minerals (smectite, celadonite, zeolite, calcite) are intermediate between primary basalt values and contemporary sea water values and indicate formation under seawaterdominated systems with effective water/rock ratios of 20-200.

Geochemistry of abyssal oceanic magmas

Massive, nearly holocrystalline dolerites from DSDP Hole 417D contain from 0.5 to 1.5% of granophyric patches composed mainly of Na-plagioclase and quartz. These patches are compositionally similar to other crystalline silicic rocks from oceanic spreading centers and differ from rarer abyssal silicic glasses. Crystalline varieties with SiO2 > 60 wt.% generally have Na/K >10, whereas silicic glasses have Na/K in the range 3-6. While crystal fractionation readily accounts for the Na2O and K2O contents of abyssal silicic glasses, both the 417D granophyres and other crystalline abyssal silicic rocks have much lower K2O than that predicted by any reasonable crystal-liquid fractionation model. We propose that high-temperature vapor phase transport is responsible for removal of potassium during late-stage crystallization of these rocks. This allows for the formation of cogenetic silicic glassy and crystalline rocks with greatly different Na/K ratios. These observations and interpretations lead to a more confident assignment of high Na/K silicic rocks of oceanic and ophiolitic environments to a cogenetic origin with basaltic oceanic crust.