Solubility and Diffusivity of Water in Lunar Basalt, and Chemical Zonation in Olivine-Hosted Melt Inclusions

I report the solubility and diffusivity of water in lunar basalt and an iron-free basaltic analogue at 1 atm and 1350 °C. Such parameters are critical for understanding the degassing histories of lunar pyroclastic glasses. Solubility experiments have been conducted over a range of fO₂ conditions from three log units below to five log units above the iron-wüstite buffer (IW) and over a range of pH₂/pH₂O from 0.03 to 24. Quenched experimental glasses were analyzed by Fourier transform infrared spectroscopy (FTIR) and secondary ionization mass spectrometry (SIMS) and were found to contain up to ~420 ppm water. Results demonstrate that, under the conditions of our experiments: (1) hydroxyl is the only H-bearing species detected by FTIR; (2) the solubility of water is proportional to the square root of pH₂O in the furnace atmosphere and is independent of fO₂ and pH₂/pH₂O; (3) the solubility of water is very similar in both melt compositions; (4) the concentration of H₂ in our iron-free experiments is <3 ppm, even at oxygen fugacities as low as IW-2.3 and pH₂/pH₂O as high as 24; and (5) SIMS analyses of water in iron-rich glasses equilibrated under variable fO₂ conditions can be strongly influenced by matrix effects, even when the concentrations of water in the glasses are low. Our results can be used to constrain the entrapment pressure of the lunar melt inclusions of Hauri et al. (2011).

Diffusion experiments were conducted over a range of fO₂ conditions from IW-2.2 to IW+6.7 and over a range of pH₂/pH₂O from nominally zero to ~10. The water concentrations measured in our quenched experimental glasses by SIMS and FTIR vary from a few ppm to ~430 ppm. Water concentration gradients are well described by models in which the diffusivity of water (D^*_water) is assumed to be constant. The relationship between D^*_water and water concentration is well described by a modified speciation model (Ni et al. 2012) in which both molecular water and hydroxyl are allowed to diffuse. The success of this modified speciation model for describing our results suggests that we have resolved the diffusivity of hydroxyl in basaltic melt for the first time. Best-fit values of D^*_water for our experiments on lunar basalt vary within a factor of ~2 over a range of pH₂/pH₂O from 0.007 to 9.7, a range of fO₂ from IW-2.2 to IW+4.9, and a water concentration range from ~80 ppm to ~280 ppm. The relative insensitivity of our best-fit values of D^*_water to variations in pH₂ suggests that H₂ diffusion was not significant during degassing of the lunar glasses of Saal et al. (2008). D^*_water during dehydration and hydration in H₂/CO₂ gas mixtures are approximately the same, which supports an equilibrium boundary condition for these experiments. However, dehydration experiments into CO₂ and CO/CO₂ gas mixtures leave some scope for the importance of kinetics during dehydration into H-free environments. The value of D^*_water chosen by Saal et al. (2008) for modeling the diffusive degassing of the lunar volcanic glasses is within a factor of three of our measured value in our lunar basaltic melt at 1350 °C.

In Chapter 4 of this thesis, I document significant zonation in major, minor, trace, and volatile elements in naturally glassy olivine-hosted melt inclusions from the Siqueiros Fracture Zone and the Galapagos Islands. Components with a higher concentration in the host olivine than in the melt (MgO, FeO, Cr₂O₃, and MnO) are depleted at the edges of the zoned melt inclusions relative to their centers, whereas except for CaO, H₂O, and F, components with a lower concentration in the host olivine than in the melt (Al₂O₃, SiO₂, Na₂O, K₂O, TiO₂, S, and Cl) are enriched near the melt inclusion edges. This zonation is due to formation of an olivine-depleted boundary layer in the adjacent melt in response to cooling and crystallization of olivine on the walls of the melt inclusions concurrent with diffusive propagation of the boundary layer toward the inclusion center.

Concentration profiles of some components in the melt inclusions exhibit multicomponent diffusion effects such as uphill diffusion (CaO, FeO) or slowing of the diffusion of typically rapidly diffusing components (Na₂O, K₂O) by coupling to slow diffusing components such as SiO₂ and Al₂O₃. Concentrations of H2O and F decrease towards the edges of some of the Siqueiros melt inclusions, suggesting either that these components have been lost from the inclusions into the host olivine late in their cooling histories and/or that these components are exhibiting multicomponent diffusion effects.

A model has been developed of the time-dependent evolution of MgO concentration profiles in melt inclusions due to simultaneous depletion of MgO at the inclusion walls due to olivine growth and diffusion of MgO in the melt inclusions in response to this depletion. Observed concentration profiles were fit to this model to constrain their thermal histories. Cooling rates determined by a single-stage linear cooling model are 150–13,000 °C hr^-1 from the liquidus down to ~1000 °C, consistent with previously determined cooling rates for basaltic glasses; compositional trends with melt inclusion size observed in the Siqueiros melt inclusions are described well by this simple single-stage linear cooling model. Despite the overall success of the modeling of MgO concentration profiles using a single-stage cooling history, MgO concentration profiles in some melt inclusions are better fit by a two-stage cooling history with a slower-cooling first stage followed by a faster-cooling second stage; the inferred total duration of cooling from the liquidus down to ~1000 °C is 40 s to just over one hour.

Based on our observations and models, compositions of zoned melt inclusions (even if measured at the centers of the inclusions) will typically have been diffusively fractionated relative to the initially trapped melt; for such inclusions, the initial composition cannot be simply reconstructed based on olivine-addition calculations, so caution should be exercised in application of such reconstructions to correct for post-entrapment crystallization of olivine on inclusion walls. Off-center analyses of a melt inclusion can also give results significantly fractionated relative to simple olivine crystallization.

All melt inclusions from the Siqueiros and Galapagos sample suites exhibit zoning profiles, and this feature may be nearly universal in glassy, olivine-hosted inclusions. If so, zoning profiles in melt inclusions could be widely useful to constrain late-stage syneruptive processes and as natural diffusion experiments.

Investigations into the degassing and eruption mechanisms of Nyamuragira volcano, Democratic Republic of the Congo (Africa)

One of two active volcanoes in the western branch of the East African Rift, Nyamuragira (1.408ºS, 29.20ºE; 3058 m) is located in the D.R. Congo. Nyamuragira emits large amounts of SO2 (up to ~1 Mt/day) and erupts low-silica, alkalic lavas, which achieve flow rates of up to ~20 km/hr. The source of the large SO2 emissions and pre-eruptive magma conditions were unknown prior to this study, and 1994-2010 lava volumes were only recently mapped via satellite imagery, mainly due to the region’s political instability. In this study, new olivine-hosted melt inclusion volatile (H2O, CO2, S, Cl, F) and major element data from five historic Nyamuragira eruptions (1912, 1938, 1948, 1986, 2006) are presented. Melt compositions derived from the 1986 and 2006 tephra samples best represent pre-eruptive volatile compositions because these samples contain naturally glassy inclusions that underwent less post-entrapment modification than crystallized inclusions. The total amount of SO2 released from the 1986 (0.04 Mt) and 2006 (0.06 Mt) eruptions are derived using the petrologic method, whereby S contents in melt inclusions are scaled to erupted lava volumes. These amounts are significantly less than satellite-based SO2 emissions for the same eruptions (1986 = ~1 Mt; 2006 = ~2 Mt). Potential explanations for this observation are: 1) accumulation of a vapor phase within the magmatic system that is only released during eruptions, and/or 2) syn-eruptive gas release from unerupted magma. Post-1994 Nyamuragira lava volumes were not available at the beginning of this study. These flows (along with others since 1967) are mapped with Landsat MSS, TM, and ETM+, Hyperion, and ALI satellite data and combined with published flow thicknesses to derive volumes. Satellite remote sensing data was also used to evaluate Nyamuragira SO2 emissions. These results show that the most recent Nyamuragira eruptions injected SO2 into the atmosphere between 15 km (2006 eruption) and 5 km (2010 eruption). This suggests that past effusive basaltic eruptions (e.g., Laki 1783) are capable of similar plume heights that reached the upper troposphere or tropopause, allowing SO2 and resultant aerosols to remain longer in the atmosphere, travel farther around the globe, and affect global climates.

Continuous eclogite melting and variable refertilisation in upwelling heterogeneous mantle

Large-scale tectonic processes introduce a range of crustal lithologies into the Earth's mantle. These lithologies have been implicated as sources of compositional heterogeneity in mantle-derived magmas. The model being explored here assumes the presence of widely dispersed fragments of residual eclogite (derived from recycled oceanic crust), stretched and stirred by convection in the mantle. Here we show with an experimental study that these residual eclogites continuously melt during upwelling of such heterogeneous mantle and we characterize the melting reactions and compositional changes in the residue minerals. The chemical exchange between these partial melts and more refractory peridotite leads to a variably metasomatised mantle. Re-melting of these metasomatised peridotite lithologies at given pressures and temperatures results in diverse melt compositions, which may contribute to the observed heterogeneity of oceanic basalt suites. We also show that heterogeneous upwelling mantle is subject to diverse local freezing, hybridization and carbonate-carbon-silicate redox reactions along a mantle adiabat.

Element concentrations of plagioclase and clinopyroxene in poikilitic olivine gabbros from ODP Hole 153-923A (Table 4)

Ocean Drilling Program Hole 923A, located on the western flank of the Mid-Atlantic Ridge south of the Kane Fracture Zone, recovered primitive gabbros that have mineral trace element compositions inconsistent with growth from a single parental melt. Plagioclase crystals commonly show embayed anorthitic cores overgrown by more albitic rims. Ion probe analyses of plagioclase cores and rims show consistent differences in trace element ratios, indicating variation in the trace element characteristics of their respective parental melts. This requires the existence of at least two distinct melt compositions within the crust during the generation of these gabbros. Melt compositions calculated to be parental to plagioclase cores are depleted in light rare earth elements, but enriched in yttrium, compared to basalts from this region of the Mid-Atlantic Ridge, which are normal mid-ocean ridge basalt (N-MORB). Clinopyroxene trace element compositions are similar to those predicted to be in equilibrium with N-MORB. However, primitive clinopyroxene crystals are much more magnesian than those produced in one-atmosphere experiments on N-MORB, suggesting that the major element composition of the melt was unlike N-MORB. These data require that the diverse array of melt compositions generated within the mantle beneath mid-ocean ridges are not always fully homogenised during melt extraction from the mantle and that the final stage of mixing can occur efficiently within crustal magma chambers. This has implications for the process of melt extraction from the mantle and the liquid line of descent of MORB

Trace element concentrations in zircons from ODP Holes 147-894G and 153-923A

The trace element compositions of Hadean zircons have been used in two ways to argue for the existence of Hadean continental crust. One argument is based on low crystallization temperatures of Hadean zircons that have been determined using a novel geothermometer based on the Ti content of zircons in equilibrium with rutile. The second argument is based on using the trace element abundances in zircons to calculate their parental melt compositions, especially the rare earth elements. Here we demonstrate that zircons that grow from a melt formed by basalt differentiation at modern mid-ocean ridges cannot be unambiguously distinguished from Hadean zircons on either of these grounds. Thus, we conclude that the trace element compositions of Hadean zircons are permissive of models that do not include the generation of continental crust in the Hadean.

Stable oxygen isotope ratios in zircons from oceanic crust

Lower ocean crust is primarily gabbroic, although 1-2% felsic igneous rocks that are referred to collectively as plagiogranites occur locally. Recent experimental evidence suggests that plagiogranite magmas can form by hydrous partial melting of gabbro triggered by seawater-derived fluids, and thus they may indicate early, high-temperature hydrothermal fluid circulation. To explore seawater-rock interaction prior to and during the genesis of plagiogranite and other late-stage magmas, oxygen-isotope ratios preserved in igneous zircon have been measured by ion microprobe. A total of 197 zircons from 43 plagiogranite, evolved gabbro, and hydrothermally altered fault rock samples have been analyzed. Samples originate primarily from drill core acquired during Ocean Drilling Program and Integrated Ocean Drilling Program operations near the Mid-Atlantic and Southwest Indian Ridges. With the exception of rare, distinctively luminescent rims, all zircons from ocean crust record remarkably uniform d18O with an average value of 5.2 ± 0.5 per mil (2SD). The average d18O(Zrc) would be in magmatic equilibrium with unaltered MORB [d18O(WR) ~5.6-5.7 per mil], and is consistent with the previously determined value for equilibrium with the mantle. The narrow range of measured d18O values is predicted for zircon crystallization from variable parent melt compositions and temperatures in a closed system, and provides no indication of any interactions between altered rocks or seawater and the evolved parent melts. If plagiogranite forms by hydrous partial melting, the uniform mantle-like d18O(Zrc) requires melting and zircon crystallization prior to significant amounts of water-rock interactions that alter the protolith d18O. Zircons from ocean crust have been proposed as a tectonic analog for >3.9 Ga detrital zircons from the earliest (Hadean) Earth by multiple workers. However, zircons from ocean crust are readily distinguished geochemically from zircons formed in continental crustal environments. Many of the >3.9 Ga zircons have mildly elevated d18O (6.0-7.5 per mil), but such values have not been identified in any zircons from the large sample suite examined here. The difference in d18O, in combination with newly acquired lithium concentrations and published trace element data, clearly shows that the >3.9 Ga detrital zircons did not originate by processes analogous to those in modern mid-ocean ridge settings.

Mineral composition of mid-ocean-ridge basalts from ODP Leg 187 sites

The chemical compositions of olivine, plagioclase, pyroxene, and spinel in lavas collected during Ocean Drilling Program Leg 187 in the Australian Antarctic Discordance, Southeast Indian Ridge (41°-46°S, 126°-135°E) were analyzed, and modeling of the theoretical equilibrium petrogenetic conditions between olivine and melt was conducted. The cores of larger olivine phenocrysts, particularly in the isotopic Indian-type mid-ocean-ridge basalt (MORB), are not equilibrated with melt compositions and are considered to be xenocrystic. Larger plagioclase phenocrysts with compositionally reversed zonation are also xenocrystic. The compositions of primary magma were calculated using a "maximum olivine fractionation" model for primitive MORB that should fractionate only olivine. Olivine compositions equilibrated with calculated primary magma and compositions of calculated primary magma suggest that (1) isotopic Pacific-type MORB is more fractionated than Indian-type MORB, (2) Pacific-type MORB was produced by higher degrees of partial melting than Indian-type MORB, and (3) primary magma for Indian-type MORB was segregated from mantle at 10 kbar (~30 km depth), whereas that for Pacific-type MORB was segregated at 15 kbar (~45 km depth).

Whole rock and Cr-rich spinel geochemistry of ODP Legs 152 and 163 and DSDP Hole 38-338 samples

Coring during Ocean Drilling Program and Deep Sea Drilling Project Legs 163, 152, 104, 81, and 38 recovered sequences of altered basalt from North Atlantic seaward-dipping reflector sequences (SDRS) erupted during the initial rifting of Greenland from northern Europe and likely associated with excessive mantle temperatures caused by an impacting mantle plume head. Cr-rich spinel is found abundantly as inclusions and groundmass crystals within the olivine-rich lavas of Hole 917A (Leg 152) cored into the Southeast Greenland SDRS, but only rarely as inclusions within plagioclase in the lavas of the Vøring Plateau SDRS, and it is absent from other cored SDRS lavas from the Rockall Plateau and Southeast Greenland. Eruptive melt compositions determined from inferred, thermodynamically-defined, spinel-melt exchange equilibria indicate that the most primitive melts represented by Hole 917A basalts have Mg/(Mg + Fe2+) at least as high as 0.70 and approach near-primary mantle melt compositions. In contrast, Cr-rich spinels from Hole 338 (Leg 38) lavas on the Vøring Plateau SDRS give evidence for melt with Mg/(Mg + Fe2+) only as high as 0.64. This study underlines that primitive melts similar to those from Hole 917A comprise only a small fraction of the eruptive North Atlantic SDRS melts, and that most SDRS basalts were, in fact, too evolved to have precipitated Cr-rich spinel, with true melt Mg/(Mg + Fe2+) likely below 0.60. The evolved nature of the SDRS basalts implies large amounts of fractionation at the base of the crust or deep within it, consistent with seismic results that indicate an abnormally thick Layer 3 underlying the SDRS.

Modes of occurrence and geochemistry at DSDP Holes 69-505A and 69-505B

The compositions of chrome spinels of Costa Rica Rift basalts from Deep Sea Drilling Project Site 505 vary depending on their occurrences as (1) inclusions in olivine crystals, (2) inclusions in Plagioclase crystals, and (3) isolated crystals in variolitic or glassy samples. The variations are a consequence of (1) changes of melt compositions as crystallization proceeds, and (2) contrasting behavior of olivine and Plagioclase in competition with spinels for Al and Mg. Some spinels have skeletal rims compositionally less magnesian than mineral cores; however, the cores do not appear to be xenocrysts, unlike some texturally similar spinels in Mid-Atlantic Ridge basalts.

The igneous charnockite—high-K alkali-calcic I-type granite—incipient charnockite association in Trivandrum Block, southern India

The Pan-African (640 Ma) Chengannoor granite intrudes the NW margin of the Neoproterozoic high-grade metamorphic terrain of the Trivandrum Block (TB), southern India, and is spatially associated with the Cardamom hills igneous charnockite massif (CM). Geochemical features characterize the Chengannoor granite as high-K alkali-calcic I-type granite. Within the constraints imposed by the high temperature, anhydrous, K-rich nature of the magmas, comparison with recent experimental studies on various granitold source compositions, and trace- and rare-earth-element modelling, the distinctive features of the Chengannoor granite reflect a source rock of igneous charnockitic nature. A petrogenetic model is proposed whereby there was a period of basaltic underplating; the partial melting of this basaltic lower crust formed the CM charnockites. The Chengannoor granite was produced by the partial melting of the charnoenderbites from the CM, with subsequent fractionation dominated by feldspars. In a regional context, the Chengannoor I-type granite is considered as a possible heat source for the near-UHT nature of metamorphism in the northern part of the TB. This is different from previous studies, which favoured CM charnockite as the major heat source. The Occurrence of incipient charnockites (both large scale as well as small scale) adjacent to the granite as well as pegmatites (which contain CO2, CO2-H2O, F and other volatiles), suggests that the fluids expelled from the alkaline magma upon solidification generated incipient charnockites through fluid-induced lowering of water activity. Thus the granite and associated alkaline pegmatites acted as conduits for the transfer of heat and volatiles in the Achankovil Shear Zone area, causing pervasive as well as patchy charnockite formation. The transport Of CO2 by felsic melts through the southern Indian middle crust is suggested to be part of a crustal-scale fluid system that linked mantle heat and CO2 input with upward migration of crustally derived felsic melts and incipient charnockite formation, resulting in an igneous charnockite - I-type granite - incipient charnockite association.

Trace element composition of IODP Site 304-U1309 and 305-U1309 (Table A2)

IODP Site U1309 was drilled at Atlantis Massif, an oceanic core complex, at 30°N on the Mid-Atlantic Ridge (MAR). We present the results of a bulk rock geochemical study (major and trace elements) carried out on 228 samples representative of the different lithologies sampled at this location. Over 96% of Hole U1309D is made up of gabbroic rocks. Diabases and basalts cross-cut the upper part of the section; they have depleted MORB compositions similar to basalts sampled at MAR 30°N. Relics of mantle were recovered at shallow depth. Mantle peridotites show petrographic and geochemical evidence of extensive melt-rock interactions. Gabbroic rocks comprise: olivine-rich troctolites (> 70% modal olivine) and troctolites having high Mg# (82-89), high Ni (up to 2300 ppm) and depleted trace element compositions (Yb 0.06-0.8 ppm); olivine gabbros and gabbros (including gabbronorites) with Mg# of 60-86 and low trace element contents (Yb 0.125-2.5 ppm); and oxide gabbros and leucocratic dykes with low Mg# (< 50), low Ni (~65 ppm) and high trace element contents (Yb up to 26 ppm). Troctolites and gabbros are amongst the most primitive and depleted oceanic gabbroic rocks. The main geochemical characteristics of Site U1309 gabbroic rocks are consistent with a formation as a cumulate sequence after a common parental MORB melt, although (lack of systematic) downhole variations indicate that the gabbroic series were built by multiple magma injections. In detail, textural and geochemical variations in olivine-rich troctolites and gabbronorites suggest chemical interaction (assimilation?) between the parental melt and the intruded lithosphere. Site U1309 gabbroic rocks do not represent the complementary magmatic product of 30°N volcanics, although they sample the same mantle source. The bulk trace element composition of Site U1309 gabbroic rocks is similar to primitive MORB melt compositions; this implies that there was no large scale removal of melts from this gabbro section. The occurrence of such a large magmatic sequence implies that a high magmatic activity is associated with the formation of Atlantis Massif. Our results suggest that almost all melts feeding this magmatic system stays trapped into the intruded lithosphere.