Fluid‐mediated mass transfer between mafic and ultramafic rocks in subduction zones

© The Author(s), 2022. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Codillo, E., Klein, F., Dragovic, B., Marschall, H., Baxter, E., Scambelluri, M., & Schwarzenbach‬, E. Fluid‐mediated mass transfer between mafic and ultramafic rocks in subduction zones. Geochemistry Geophysics Geosystems, 23, (2022): e2021GC010206, https://doi.org/10.1029/2021gc010206.Metasomatic reaction zones between mafic and ultramafic rocks exhumed from subduction zones provide a window into mass-transfer processes at high pressure. However, accurate interpretation of the rock record requires distinguishing high-pressure metasomatic processes from inherited oceanic signatures prior to subduction. We integrated constraints from bulk-rock geochemical compositions and petrophysical properties, mineral chemistry, and thermodynamic modeling to understand the formation of reaction zones between juxtaposed metagabbro and serpentinite as exemplified by the Voltri Massif (Ligurian Alps, Italy). Distinct zones of variably metasomatized metagabbro are dominated by chlorite, amphibole, clinopyroxene, epidote, rutile, ilmenite, and titanite between serpentinite and eclogitic metagabbro. Whereas the precursor serpentinite and oxide gabbro formed and were likely already in contact in an oceanic setting, the reaction zones formed by diffusional Mg-metasomatism between the two rocks from prograde to peak, to retrograde conditions in a subduction zone. Metasomatism of mafic rocks by Mg-rich fluids that previously equilibrated with serpentinite could be widespread along the subduction interface, within the subducted slab, and the mantle wedge. Furthermore, the models predict that talc formation by Si-metasomatism of serpentinite in subduction zones is limited by pressure-dependent increase in the silica activity buffered by the serpentine-talc equilibrium. Elevated activities of aqueous Ca and Al species would also favor the formation of chlorite and garnet. Accordingly, unusual conditions or processes would be required to stabilize abundant talc at high P-T conditions. Alternatively, a different set of mineral assemblages, such as serpentine- or chlorite-rich rocks, may be controlling the coupling-decoupling transition of the plate interface.M. Scambelluri acknowledges the Italian Ministry of Research MUR for granting the PRIN project n. 2017ZE49E7. This research was funded by NSF-OISE (Office of International Science & Engineering, Petrology & Geochemistry) PIRE, Award #1545903, and the WHOI Ocean Ventures Fund

Preferential formation of chlorite over talc during Si-metasomatism of ultramafic rocks in subduction zones

Talc formation via silica-metasomatism of ultramafic rocks is believed to play key roles in subduction zone processes. Yet, the conditions of talc formation remain poorly constrained. We used thermodynamic reaction-path models to assess the formation of talc at the slab-mantle interface and show that it is restricted to a limited set of pressure–temperature conditions, protolith, and fluid compositions. In contrast, our models predict that chlorite formation is ubiquitous at conditions relevant to the slab-mantle interface of subduction zones. The scarcity of talc and abundance of chlorite is evident in the rock record of exhumed subduction zone terranes. Talc formation during Si-metasomatism may thus play a more limited role in volatile cycling, strain localization, and in controlling the decoupling-coupling transition of the plate interface. Conversely, the observed and predicted ubiquity of chlorite corroborates its prominent role in slab-mantle interface processes that previous studies attributed to talc. Key Points: Limited talc formation by Si-metasomatism of ultramafic rocks in subduction zones Chlorite formation is likely pervasive at the slab-mantle interface Preferential formation of chlorite has wide-ranging chemical and physical implications for subduction zone processes Plain Language Summary: In subduction zones, talc can form during chemical reactions of mantle rocks with silica-enriched fluids at the interface between descending oceanic plates and the overriding mantle. Its formation and distribution in subduction zones are believed to affect the volatile budget, rheological properties, and the down-dip limit of the decoupling of the slab-mantle interface. Therefore, illuminating the conditions that facilitate talc formation at high pressure-temperature conditions is key in assessing its roles in fundamental subduction zone processes. Using thermodynamic reaction-path models, we show that the formation of talc at the slab-mantle interface is restricted to a limited set of environmental conditions, because its formation is highly sensitive to the compositions of the mantle rocks and reactant fluids. Contrary to common belief, talc is unlikely to form in high abundance in ultramafic rocks metasomatized by Si-rich slab-derived fluids. Rather, our models predict the ubiquitous formation of chlorite along with other silicate minerals during Si-metasomatism due to the competing effects from other dissolved components that favor their formation over talc. This study calls into question the importance of talc during Si-metasomatism in subduction zones but highlights the more predominant role of chlorite

Causes of oceanic crustal thickness oscillations along a 74-Myr Mid-Atlantic Ridge flow line

These data sets collected geophysical data: multi-beam bathymetry, gravity, magnetics, sub-bottom profile to investigate the relationships between faulting, magmatism, and sea level change.Gravity, magnetic, and bathymetry data collected along a continuous 1400-km-long spreading-parallel flow line across the Mid-Atlantic Ridge indicate significant tectonic and magmatic fluctuations in the formation of oceanic crust over a range of timescales. The transect spans from 28 Ma on the African Plate to 74 Ma on the North American plate, crossing the Mid-Atlantic Ridge at 35.8 ºN. Gravity-derived crustal thicknesses vary from 3–9 km with a standard deviation of 1 km. Spectral analysis of bathymetry and residual mantle Bouguer anomaly (RMBA) show diffuse power at >1 Myr and concurrent peaks at 390, 550, and 950 kyr. Large-scale (>10-km) mantle thermal and compositional heterogeneities, variations in upper mantle flow, and detachment faulting likely generate the >1 Myr diffuse power. The 550- and 950-kyr peaks may reflect the presence of magma solitons and/or regularly spaced ~7.7 and 13.3 km short-wavelength mantle compositional heterogeneities. The 390-kyr spectral peak corresponds to the characteristic spacing of faults along the flow line. Fault spacing also varies over longer periods (>10 Myr), which we interpret as reflecting long-lived changes in the fraction of tectonically- vs. magmatically- accommodated extensional strain. A newly discovered off-axis oceanic core complex (Kafka Dome) found at 8 Ma on the African plate further suggests extended time periods of tectonically dominated plate separation. Fault spacing negatively correlates with gravity-derived crustal thickness, supporting a strong link between magma input and fault style at mid-ocean ridges

Causes of Oceanic Crustal Thickness Oscillations Along a 74‐M Mid‐Atlantic Ridge Flow Line

International audience1. We conducted one of the longest continuous geophysical surveys along a 74-Myr 15 spreading-parallel flow line across the Mid-Atlantic Ridge. 16 2. Spectral densities of bathymetry and gravity data show concurrent peaks at 390-, 550-, 17 and 950-kyr periods, and diffuse power at >1 Myr. 18 3. A negative correlation between fault spacing and gravity-derived crustal thickness 19 suggests a link between magma input and fault style. 20 21 Abstract: 22 Gravity, magnetic, and bathymetry data collected along a continuous 1400-km-23 long spreading-parallel flow line across the Mid-Atlantic Ridge indicate significant 24 tectonic and magmatic fluctuations in the formation of oceanic crust over a range of 25 timescales. The transect spans from 28 Ma on the African Plate to 74 Ma on the North 26 American plate, crossing the Mid-Atlantic Ridge at 35.8 ºN. Gravity-derived crustal 27 thicknesses vary from 3-9 km with a standard deviation of 1 km. Spectral analysis of 28 bathymetry and residual mantle Bouguer anomaly (RMBA) show diffuse power at >1 29 Myr and concurrent peaks at 390, 550, and 950 kyr. Large-scale (>10-km) mantle thermal 30 and compositional heterogeneities, variations in upper mantle flow, and detachment 31 faulting likely generate the >1 Myr diffuse power. The 550-and 950-kyr peaks may 32 reflect the presence of magma solitons and/or regularly spaced ~7.7 and 13.3 km short-33 wavelength mantle compositional heterogeneities. The 390-kyr spectral peak corresponds 34 to the characteristic spacing of faults along the flow line. Fault spacing also varies over 35 longer periods (>10 Myr), which we interpret as reflecting long-lived changes in the 3