Hidden melting signatures recorded in the Troodos ophiolite plutonic suite: evidence for widespread generation of depleted melts and intra-crustal melt aggregation

The main plutonic complex of the Troodos ophiolite, north of the Arakapas Fault Zone, has been re-examined both from field and geochemical perspectives. Ion microprobe analyses of clinopyroxene crystal cores show that the range of melt compositions added to the lower crust far exceeds that of published lavas in the main Troodos massif. This suggests that the lower crust acted as a filter into which a large range of melt compositions were added and out of which a homogenised (and generally fractionated) derivative was extracted. This crustal-level aggregation homogenised diverse melt fractions from a broad range of degrees of melting. Depleted melts with U-shaped rare earth element (REE) patterns were a significant component of the melts added to the crust, but because of their low incompatible element abundances, mixing with less depleted melts prior to eruption masked their signature in the lavas. The discovery that highly depleted melts constituted a significant component of the melts added to the Troodos crust, but not of the lavas, demonstrates that the spatial distribution of lava-types is not necessarily a good indicator of where different parental melt compositions are generated within the mantle. Compared with normal mid-ocean ridge basalts, the Troodos parental melts were (1) generally depleted in immobile incompatible trace elements, (2) less depleted in light REE (LREE) than would be expected for the concomitant depletion in middle and heavy REE, (3) enriched in Sr with respect to the LREE and (4) more oxidised. Modelling of these characteristics suggests a mantle source that had previously lost a significant melt fraction under relatively reducing conditions. This was followed by remelting under more oxidising conditions in an environment in which Sr and LREE were added to the source consistent with previous models of a supra-subduction zone setting.     

The time scales of magma mixing and mingling involving primitive melts and melt–mush interaction at mid-ocean ridges

We have studied the chemical zoning of plagioclase phenocrysts from the slow-spreading Mid-Atlantic Ridge and the intermediate-spreading rate Costa Rica Rift to obtain the time scales of magmatic processes beneath these ridges. The anorthite content, Mg, and Sr in plagioclase phenocrysts from the Mid-Atlantic Ridge can be interpreted as recording initial crystallisation from a primitive magma (~11 wt% MgO) in an open system. This was followed by crystal accumulation in a mush zone and later entrainment of crystals into the erupted magma. The initial magma crystallised plagioclase more anorthitic than those in equilibrium with any erupted basalt. Evidence that the crystals accumulated in a mush zone comes from both: (1) plagioclase rims that were in equilibrium with a Sr-poor melt requiring extreme differentiation; and (2) different crystals found in the same thin section having different histories. Diffusion modelling shows that crystal residence times in the mush were <140 years, whereas the interval between mush disaggregation and eruption was ≤1.5 years. Zoning of anorthite content and Mg in plagioclase phenocrysts from the Costa Rica Rift show that they partially or completely equilibrated with a MgO-rich melt (>11 wt%). Partial equilibration in some crystals can be modelled as starting <1 year prior to eruption but for others longer times are required for complete equilibration. This variety of times is most readily explained if the mixing occurred in a mush zone. None of the plagioclase phenocrysts from the Costa Rica Rift that we studied have Mg contents in equilibrium with their host basalt even at their rims, requiring mixing into a much more evolved magma within days of eruption. In combination these observations suggest that at both intermediate- and slow-spreading ridges: (1) the chemical environment to which crystals are exposed changes on annual to decadal time scales; (2) plagioclase crystals record the existence of melts unlike those erupted; and (3) disaggregation of crystal mush zones appears to precede eruption, providing an efficient mechanism by which evolved interstitial melt can be mixed into erupted basalts.     

Oxygen isotope evidence for short-lived high-temperature fluid flow in the lower oceanic crust at fast-spreading ridges

Millimeter-scale amphibole veins in the lower oceanic crust record fracture-controlled fluid flow at high-temperatures but the importance of this fluid flow for the thermal and chemical evolution of the lower oceanic crust is unclear. In the section of lower oceanic crust recovered at Hess Deep from ODP Hole 894G, which formed at the fast-spreading East Pacific Rise, these veins are randomly distributed with an average spacing of  1 m. We unravel the history of fluid flow through one of these veins by combining in situ O-isotope analyses of wall-rock plagioclase with major element analyses, geothermometry and diffusion modeling. Thermometry indicates vein sealing by amphibole at  720 °C over a narrow temperature interval (± 20 °C). In situ O-isotope analyses by ion microprobe, with a precision of < 0.5‰, reveal zoning of O-isotopes in plagioclase adjacent to the vein. The zoning profiles can be reproduced using a diffusion model if the duration of O-isotope exchange was ≤ 100 yr. A similar interval of fluid–rock exchange is suggested by modeling potassium depletion in plagioclase adjacent to the vein. If representative of fracture controlled fluid flow in the lower oceanic crust the limited duration of fluid flow, and its occurrence over a narrow temperature interval, suggest that high-temperature fluid flow in this porosity network does not transport significant heat.     

Petrology and geochemistry of primitive lower oceanic crust from Pito Deep: implications for the accretion of the lower crust at the Southern East Pacific Rise

A suite of samples collected from the uppermost part of the plutonic section of the oceanic crust formed at the southern East Pacific Rise and exposed at the Pito Deep has been examined. These rocks were sampled in situ by ROV and lie beneath a complete upper crustal section providing geological context. This is only the second area (after the Hess Deep) in which a substantial depth into the plutonic complex formed at the East Pacific Rise has been sampled in situ and reveals significant spatial heterogeneity in the plutonic complex. In contrast to the uppermost plutonic rocks at Hess Deep, the rocks studied here are generally primitive with olivine forsterite contents mainly between 85 and 88 and including many troctolites. The melt that the majority of the samples crystallized from was aggregated normal mid-ocean ridge basalt (MORB). Despite this high Mg# clinopyroxene is common despite model predictions that clinopyroxene should not reach the liquidus early during low-pressure crystallization of MORB. Stochastic modeling of melt crystallisation at various levels in the crust suggests that it is unlikely that a significant melt mass crystallized in the deeper crust (for example in sills) because this would lead to more evolved shallow level plutonic rocks. Similar to the upper plutonic section at Hess Deep, and in the Oman ophiolite, many samples show a steeply dipping, axis-parallel, magmatic fabric. This suggests that vertical magmatic flow is an important process in the upper part of the seismic low velocity zone beneath fast-spreading ridges. We suggest that both temporal and spatial (along-axis) variability in the magmatic and hydrothermal systems can explain the differences observed between the Hess Deep and Pito Deep plutonics. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Experimental determination of the diffusion coefficient for calcium in olivine between 900°C and 1500°C

Experiments have been carried out to determine the temperature, oxygen fugacity (fO₂) and compositional dependence of the tracer diffusion coefficient (D) of calcium in olivine. These data constrain the diffusion coefficient over the temperature range 900 to 1500°C for the three principal crystallographic axes. Well constrained linear relationships between the reciprocal of the absolute temperature and log(D) exist at any given oxygen fugacity. There is a strong dependence of the diffusion coefficient on oxygen fugacity with D ∝ fO₂^(1/3). This makes a knowledge of the T-fO₂ path followed by geological samples a prerequisite for modelling Ca diffusion in olivine. The best fitting preexponential factor (Do) and activation energy (E) to the Arrhenius equation log (D) = log [Do exp(−E/RT)] + 0.31Δ log fO₂ for Ca diffusion in olivine at a given oxygen fugacity (fO₂*) are given by:diffusion along [100]: log [Do (m²/s)] = −10.78 ± 0.43; E = 193 ± 11 kJ/moldiffusion along [010]: log [Do (m²/s)] = −10.46 ± 0.37; E = 201 ± 10 kJ/moldiffusion along [001]: log [Do (m²/s)] = −10.02 ± 0.29; E = 207 ± 8 kJ/molwhere Δ log fO₂ = log[fO₂*] − log[10⁻¹²] with fO₂* in units of bars. There is no measurable compositional dependence of the diffusion coefficient between Fo₈₃ and Fo₉₂. Diffusion in Fo₁₀₀ has a much higher activation energy than in Fe-bearing olivine and has a weaker fO₂ dependence.     

Rates of hydrothermal cooling of new oceanic upper crust derived from lithium-geospeedometry

Episodic emplacement and cooling of lavas and dikes at mid-ocean ridges leads to large fluctuations in hydrothermal fluxes and biological activity. However, the processes operating beneath the seafloor during these transient events such as permeability creation and dike cooling are poorly understood. We have developed a new approach to determine the cooling rate of the sheeted dike complex based on the extent of diffusion of lithium from plagioclase into clinopyroxene during cooling. We have calibrated this Li-geospeedometer using new high-temperature experiments to determine both the temperature dependence of the partitioning of Li between plagioclase and clinopyroxene and the diffusion coefficient for Li in clinopyroxene. Application of this method to lavas and dikes from ODP Hole 504B shows that cooling rates vary dramatically with depth in the upper oceanic crust. Extremely rapid cooling rates (> 450 °C hr^− 1) in the upper part of the sheeted dike complex are sufficient to power hydrothermal megaplume formation within the overlying water column.     

Contrasting Cooling Rates in the Lower Oceanic Crust at Fast- and Slow-spreading Ridges Revealed by Geospeedometry

Two approaches to determining the high-temperature (1000°Cto 600°C) cooling rate of the lower oceanic crust and uppermantle are presented and critically evaluated. The first isbased on the down-temperature diffusive exchange of Ca betweenolivine and clinopyroxene. The second, less well-constrained,approach is based on the down-temperature diffusive exchangeof Mg and Fe between olivine and spinel. Cooling rates basedon olivine–spinel geospeedometry are approximately anorder of magnitude faster than those from Ca-in-olivine geospeedometry.In contrast, cooling rates derived from thermochronology andremanent magnetism are approximately an order of magnitude slowerthan those derived by Ca-in-olivine geospeedometry; this isprobably because they record cooling at lower temperatures.Using the Ca-in-olivine geospeedometer, the cooling rate ofsamples from the lower oceanic crust and upper oceanic mantleformed in the Oman ophiolite and in the three main ocean basinshas been determined. Samples from the lower oceanic crust formedat fast-spreading ridges show a large decrease in cooling ratebetween the top and base of the gabbroic section, with mostof the variation occurring within the upper kilometre. Thisis consistent with vertical heat loss (within the crustal frameof reference) dominating the thermal evolution at fast-spreadingridges. Samples from Ocean Drilling Program Hole 735B, whichformed at the slow-spreading Southwest Indian Ridge, show novariation in cooling rate over 1500 m depth range and cooledsubstantially faster than rocks from the deeper portion of thegabbros in the Oman ophiolite, where the change in cooling ratewith depth is limited. These observations are consistent withheat loss from small plutons emplaced in cool lithosphere atthe slow-spreading ridge. Alternatively, they could be explainedby cooling through the Ca-in-olivine closure interval duringuplift towards the surface. KEY WORDS: geospeedometry; lower oceanic crust; Hess Deep; Pito Deep; ODP Hole 735B; ODP Leg 153