Relationship between seismic P-wave velocity and the composition of anhydrous igneous and meta-igneous rocks

Author Posting. © American Geophysical Union 2003. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry Geophysics Geosystems 4 (2003): 1041, doi:10.1029/2002GC000393.This study presents a new approach to quantitatively assess the relationship between the composition and seismic P-wave velocity of anhydrous igneous and meta-igneous rocks. We perform thermodynamic calculations of the equilibrating phase assemblages predicted for all igneous composition space at various pressure and temperature conditions. Seismic velocities for each assemblage are then estimated from mixing theory using laboratory measurements of the elastic parameters for pure mineral phases. The resultant velocities are used to derive a direct relationship between Vp and major element composition valid to ±0.13 km/s for pressure and temperature conditions along a normal crustal geotherm in the depth range of 5–50 km and equilibration pressures ≤12 kbar. Finally, we use the calculated velocities to invert for major element chemistry as a function of P-wave velocity assuming only the in situ temperature and pressure conditions are known. Compiling typical velocity-depth profiles for the middle and lower continental and oceanic crust, we calculate compositional bounds for each of these geologic environments. We find that the acceptable compositional range for the middle (15–30 km) and lower continental (≥35 km) crust is broad, ranging from basaltic to dacitic compositions, and conclude that P-wave velocity measurements alone are insufficient to provide fundamental constraints on the composition of the middle and lower continental crust. However, because major oxides are correlated in igneous rocks, joint constraints on Vp and individual oxides can narrow the range of acceptable crustal compositions. In the case of the lower oceanic crust (≥2 km), observed velocities are 0.2–0.3 km/s lower than velocities calculated based on the average bulk composition of gabbros in drill cores and exposed ophiolite sequences. We attribute this discrepancy to a combination of residual porosity at crustal depths less than ∼10 km and hydrous alteration phases in the lower crust, and suggest caution when inferring mantle melting parameters from observed velocities in the lower oceanic crust.This research was supported by National Science Foundation Grants OCE- 9819666, EAR-9910899, and EAR-0087706 (P.B. Kelemen)

Highly depleted cratonic mantle in West Greenland extending into diamond stability field in the Proterozoic

This study presents electron microprobe data for dunite xenoliths from a lamprophyre dyke located on the island of Qeqertaa, West Greenland. The minimum age of this dyke is Palaeoproterozoic and it experienced amphibolite facies metamorphism and deformation during that era. The samples consist of nearly 200 xenoliths with a size range of 0.5-8 cm. These dunite xenoliths have olivine Mg#, that range from 80.3 to 94.6 (n = 579) with a mean of 92.6. Orthopyroxene is found in three xenoliths and garnet in five others. The latter suggests the depth of the Qeqertaa xenolith suite to be near the diamond stability-field, which is substantiated by the finding of diamonds in bulk samples of the Qeqertaa dyke. This further indicates the presence of a lithospheric mantle domain dominated by high-Mg# dunite to this depth in Palaeoproterozoic time. Cr-rich spinel, in the 0.1–0.2 mm size range, is found within and between olivine grains in individual xenoliths. These Cr-spinels yield Fe–Mg exchange temperatures of 400–600 °C. However, the presence of intermediate spinel compositions spanning the lower temperature solvus suggests that equilibration temperatures were > 550 °C. Fe3 +#, expressed as 100 × Fe3 +/(Fe3 + + Al + Cr), is shown to be a useful parameter in order to screen for altered spinel (Fe3 +# > 10) with disturbed Mg# and Cr#. The screened spinel data (Fe3 +# < 10) show a distinctly different trend in terms of spinel Cr# versus Mg#, compared to unmetamorphosed xenoliths in Tertiary lavas and dikes from Ubekendt Ejland and Wiedemann Fjord, respectively, also located within the North Atlantic craton. This difference likely reflects amphibolite facies metamorphic resetting of the Qeqertaa xenolith suite by Fe–Mg exchange. Given the similarity of the Qeqertaa xenolith suite with the Ubekendt and Wiedemann suites, in terms of their olivine Mg# and spinel Cr# distribution, high-Mg# dunite is likely to be an important component of the subcontinental lithospheric mantle beneath the North Atlantic craton and appears to have spanned a vertical distance of at least 150 km in this region, even during the Palaeoproterozoic

Peridotite enclaves hosted by Mesoarchaean TTG-suite orthogneisses in the Fiskefjord region of southern West Greenland

AbstractThis study presents bulk-rock major, trace, and platinum-group element data, as well as mineral chemistry for peridotites which form large enclaves (up to 500×1000m) within Mesoarchaean orthogneisses of the Akia terrane in the Fiskefjord region, southern West Greenland. The largest peridotite body, known as Seqi, contains highly fosteritic olivine with a median Mg# of 92.6 and hosts extensive layers of chromitite, which can be traced for tens of metres with thicknesses of up to 30cm. Thinner (<100m thick), but extensive (up to 2000m long) peridotite sheets are associated with coarse norite and orthopyroxenite with distinct cumulate textures in the Amikoq complex, located a few tens of kilometres south of Seqi. Intercalated amphibolites of tholeiitic basaltic composition show complementary geochemical evolution to the peridotites, consistent with igneous crystal fractionation trends. The U-shaped trace element patterns of the peridotites may either reflect the parental melt composition from which these olivine-rich rocks were derived, or alternatively this feature may be the result of melt-rock interaction. Overall, we interpret the Fiskefjord region peridotites to have formed as ultramafic cumulates derived from Archaean high-Mg, low Ca/Al magmas, although their geodynamic setting remains to be established

Coexisting serpentine and quartz from carbonate-bearing serpentinized peridotite in the Samail Ophiolite, Oman

Tectonically exposed mantle peridotite in the Oman Ophiolite is variably serpentinized and carbonated. Networks of young carbonate veins are prevalent in highly serpentinized peridotite, particularly near low-temperature alkaline springs emanating from the peridotite. An unusual feature in some samples is the coexistence of serpentine and quartz, which is not commonly observed in serpentinites. This assemblage is unstable with respect to serpentine + talc or talc + quartz under most conditions. Serpentine in the carbonated serpentinites in this study is more iron rich than in most serpentinites reported in previous studies, and samples with co-existing quartz contain the most iron-rich serpentines. Calculations of thermodynamic equilibria in the MgO–SiO2–H2O–CO2 system suggest that serpentine + quartz may be a stable assemblage at low temperatures (e.g., <~15–50 °C) and is stabilized to higher temperatures by preferential cation substitutions in serpentine over talc. Based on these calculations, serpentine + quartz assemblages could result from serpentinization at near-surface temperatures. Clumped isotope thermometry of carbonate veins yields temperatures within error of the observed temperatures in Oman groundwater for all samples analyzed, while the δ^18O of water calculated to be in equilibrium with carbonate precipitated at those temperatures is within error of the observed isotopic composition of Oman groundwater for the majority of samples analyzed. As groundwater geochemistry suggests that carbonate precipitation and serpentinization occur concomitantly, this indicates that both hydration and carbonation of peridotite are able to produce extensive alteration at the relatively low temperatures of the near-surface weathering environment

A detailed geochemical study of island arc crust : the Talkeetna Arc section, south–central Alaska

Author Posting. © The Author, 2006. This is the author's version of the work. It is posted here by permission of Oxford University Press for personal use, not for redistribution. The definitive version was published in Journal of Petrology 47 (2006): 1051-1093, doi:10.1093/petrology/egl002.The Early to Middle Jurassic Talkeetna Arc section exposed in the Chugach Mountains of south central Alaska is 5-18 km wide and extends for over 150 km. This accreted island arc includes exposures of upper mantle to volcanic upper crust. The section comprises six lithologic units, in order of decreasing depth: (1) residual upper mantle harzburgite (with lesser proportions of dunite); (2) pyroxenite; (3) basal gabbronorite; (4) lower crustal gabbronorite; (5) mid-crustal plutonic rocks; and (6) volcanic rocks. The pyroxenites overlie residual mantle peridotite, with some interfingering of the two along the contact. The basal gabbronorite overlies pyroxenite, again with some interfingering of the two different units along their contact. Lower crustal gabbronorite (≤10 km thick) includes abundant rocks with well developed modal layering. The mid-crustal plutonic rocks include a heterogeneous assemblage of gabbroic rocks, dioritic to tonalitic rocks (30-40% area), and concentrations of mafic dikes and chilled mafic inclusions. The volcanic rocks (~7 km thick) range from basalt to rhyolite. Many of the evolved volcanic compositions are a result of fractional crystallisation processes whose cumulate products are directly observable in the lower crustal gabbronorites. For example, Ti and Eu enrichments in lower crustal gabbronorites are mirrored by Ti and Eu depletions in evolved volcanics. In addition, calculated parental liquids from ion microprobe analyses of clinopyroxene in lower crustal gabbronorites indicate that the clinopyroxenes crystallised in equilibrium with liquids whose compositions were the same as the compositions of volcanic rocks. The compositional variation of the main series of volcanic and chilled mafic rocks can be modeled through fractionation of observed phase compositions and phase proportions in lower crustal gabbronorite (i.e. cumulates). Primary, mantle-derived melts in the Talkeetna Arc underwent fractionation of pyroxenite at the base of the crust. Our calculations suggest that more than 25 wt % of the primary melts crystallised as pyroxenites at the base of the crust. The discrepancy between the observed proportion of pyroxenites (less than 5% of the arc section) and the proportion required by crystal fractionation modeling (more than 25%) may be best understood as the result of gravitational instability, with dense ultramafic cumulates, probably together with dense garnet granulites, foundering into the underlying mantle during the time when the Talkeetna Arc was magmatically active, or in the initial phases of slow cooling (and sub-solidus garnet growth) immediately after the cessation of arc activity.This study was supported by National Science Foundation Grant EAR-9910899

A Mg isotopic perspective on the mobility of magnesium during serpentinization and carbonation of the Oman Ophiolite

Author Posting. © American Geophysical Union, 2021. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Solid Earth 126(2), (2021): e2020JB020237, https://doi.org/10.1029/2020JB020237.Alteration of mantle peridotite in the Samail ophiolite forms secondary minerals, mainly serpentine and Mg‐rich carbonates. Magnesium accounts for ∼25 – 30% of peridotite mass and its mobility can be used to trace this alteration. We report the first set of Mg isotope measurements from peridotites and their alteration products in Oman. Partially serpentinized peridotites have Mg isotope ratios that are indistinguishable from estimates for the average mantle and bulk silicate earth (δ26Mg = −0.25 ± 0.04‰). However, more extensively altered peridotite samples show large shifts in Mg isotopic composition. The range of δ26Mg values for our suite of alteration products from the mantle section is ∼4.5‰ (from −3.39‰ to 1.19‰), or >60% of the total range of terrestrial variability in δ26Mg values. Serpentine veins are typically enriched in 26Mg (max δ26Mg value = 0.96‰) whereas Mg‐carbonate veins are associated with low 26Mg/24Mg ratios (magnesite δ26Mg = −3.3‰, dolomite δ26Mg = −1.91‰). Our preferred explanation for the range in δ26Mg values involves coprecipitation of serpentine and carbonates at water‐to‐rock ratios >103. The coincidence of alteration products characterized by δ26Mg values that are both lower and higher than bulk silicate Earth and the finite 14C ages of the carbonates suggest that both serpentinization and carbonation are ongoing in Oman. Rates of calcite precipitation in travertines inferred from Δ26Mgcal‐fl suggest that travertine formation in Oman sequesters a total of 106–107 kg CO2/yr, consistent with previous estimates.This work was supported through the Sloan Foundation–Deep Carbon Observatory (Grant 2014‐3‐01, Kelemen PI), the U.S.‐National Science Foundation (NSF‐EAR‐1516300, Kelemen lead PI).2021-06-0

Timing of magnetite growth associated with peridotite-hosted carbonate veins in the SE Samail ophiolite, Wadi Fins, Oman

Author Posting. © American Geophysical Union, 2020. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Solid Earth 125(5), (2020): e2019JB018632, doi:10.1029/2019JB018632.Carbonate‐altered peridotite are common in continental and oceanic settings and it has been suggested that peridotite‐hosted carbonate represent a significant component of the carbon‐cycle and provide an important link in the CO2 dynamics between the atmosphere, hydrosphere, and lithosphere. The ability to constrain the timing of carbonate and accessory phase growth is key to interpreting the mechanisms that contribute to carbonate alteration, veining, and mineralization in ultramafic rocks. Here we examine a mantle section of the Samail ophiolite exposed in Wadi Fins in southeastern Oman where the peridotite is unconformably overlain by Late Cretaceous‐Paleogene limestone and crosscut by an extensive network of carbonate veins and fracture‐controlled alteration. Three previous 87Sr/86Sr measurements on carbonate vein material in the peridotite produce results consistent with vein formation involving Cretaceous to Eocene seawater (de Obeso & Kelemen, 2018, https://doi.org/10.1098/rsta.2018.0433). We employ (U‐Th)/He chronometry to constrain the timing of hydrothermal magnetite in the calcite veins in the peridotite. Magnetite (U‐Th)/He ages of crystal sizes ranging from 1 cm to 200 μm record Miocene growth at 15 ± 4 Ma, which may indicate (1) fluid–rock interaction and carbonate precipitation in the Miocene, or (2) magnetite (re)crystallization within pre‐existing veins. Taken together with published Sr‐isotope values, these results suggest that carbonate veining at Wadi Fins started as early as the Cretaceous, and continued in the Miocene associated with magnetite growth. The timing of hydrothermal magnetite growth is coeval with Neogene shortening and faulting in southern Oman, which points to a tectonic driver for vein (re)opening and fluid‐rock alteration.This research was supported by a National Science Foundation (NSF) Graduate Research Fellowship to E.H.G. Cooperdock, the UTChron Laboratory at The University of Texas at Austin, the Chevron (Gulf) Centennial Professorship to D.F. Stockli, and by a Sloan Foundation grant awarded to P.B. Kelemen. We are grateful to Desmond Patterson for assistance and training with He measurements and data reduction, to Jessie Maisano for technical support with the X‐Ray Computed Tomography. These data and images were produced at the High‐Resolution X‐ray Computed Tomography Facility of the University of Texas at Austin. EHGC is grateful to Jaime Barnes, Richard Ketcham, Frieder Klein and Othmar Müntener for helpful comments on an earlier version of this manuscript. Thank you to Fin Stuart and Uwe Ring for their helpful reviews, and Stephen Parman for feedback and editorial handling of the manuscript. The (U‐Th)/He data in this manuscript are available in the GeoChron repository (https://www.geochron.org) and sample IGSNs are in the SESAR database (http://www.geosamples.org).2020-10-0

Diapirs as the source of the sediment signature in arc lavas

Author Posting. © The Author(s), 2011. This is the author's version of the work. It is posted here by permission of Nature Publishing Group for personal use, not for redistribution. The definitive version was published in Nature Geoscience 4 (2011): 641-646, doi:10.1038/ngeo1214Many arc lavas show evidence for the involvement of subducted sediment in the melting process. There is debate whether this “sediment melt” signature forms at relatively low temperature near the fluid-saturated solidus or at higher temperature beyond the breakdown of trace-element-rich accessory minerals. We present new geochemical data from high- to ultrahigh-pressure rocks that underwent subduction and show no significant depletion of key trace elements in the sediment melt component until peak metamorphic temperatures exceeded ~1050ºC from 2.7 to 5 GPa. These temperatures are higher than for the top of the subducting plate at similar pressures based on thermal models. To address this discrepancy, we use instability calculations for a non-Newtonian buoyant layer in a viscous half-space to show that, in typical subduction zones, solid-state sediment diapirs initiate at temperatures between 500–850ºC. Based on these calculations, we propose that the sediment melt component in arc magmas is produced by high degrees of dehydration melting in buoyant diapirs of metasediment that detach from the slab and rise into the hot mantle wedge. Efficient recycling of sediments into the wedge by this mechanism will alter volatile fluxes into the deep mantle compared to estimates based solely on devolatilization of the slab.Funding for this work was provided by NSF and WHOI’s Deep Ocean Exploration Institute