Search

Abstract Subducting oceanic plates carry large amounts of carbon into the Earth’s interior. The subducted carbon is mobilized by fluid and encounters ultramafic rocks in the mantle wedge, resulting in changes to the mineral assemblage and mechanical properties of the mantle. Here, we use thermodynamic modeling of interactions between carbon-bearing multi-component fluids and mantle rocks to investigate the down-dip variation in mineral assemblage in the forearc mantle along subduction megathrusts. We found that fluids rich in aqueous carbon are preferentially generated in a warm subduction zone (e.g., Nankai, SW Japan), causing a change in mineral assemblage from serpentine-rich at the mantle wedge corner to talc + carbonate-rich at greater depths. The transition caused by the infiltration of aqueous carbon may influence the depth of the boundary between the seismogenic and aseismic zones, and the down-dip limit of episodic tremor and slip.

Abstract The slab–mantle interface in subduction zones is one of the geological boundaries with the most significant chemical potential gradients, which leads to fluid-mediated metasomatic reactions and chemical transport. As subducting sediment and basaltic crust often contain silica in various forms, the Si-metasomatism of mantle rocks is thought to occur along the subduction zone interface. However, growing evidence from the geochemistry of altered rocks and thermodynamic modelling has revealed the presence of multi-component fluids at the slab interface. Here, we review the laboratory experiments, geochemical models, and natural observations that improve our understanding of mass transport and metasomatic reactions at the crust–mantle interface, focusing on the relative mobility of Mg and Si. Hydrothermal experiments using analogues for the boundary between mantle (olivine) and crust (quartz or plagioclase) under vapor-saturated pressures indicate that Si is preferentially transported from crust to mantle, whereas Mg is immobile. This result is consistent with the distribution of talc rocks in oceanic lithosphere. On the other hand, at the contact between ultramafic (e.g., serpentinite) and crustal (pelitic schist or basaltic rocks) rocks in high-pressure metamorphic terranes, a large volume of chlorite rocks form in the crustal rocks, and the volume of chlorite often exceeds talc in serpentinites. Geochemical modeling reveals that in the shallow part of a subduction zone, the dissolved Si content of fluids in equilibrium with pelitic schist (C Si,crust) is significantly higher than the dissolved Mg content of fluids in equilibrium with mantle peridotite (C Mg,mantle); however, C Mg,mantle becomes dominant at depth, resulting in the Mg-metasomatism of crustal rocks to form chlorite rocks. This Mg-metasomatism is more widespread in warmer subduction zones (e.g., the Nankai and Cascadia subduction zones) than in colder subduction zones (e.g., in Northeast Japan). In addition, the infiltration of CO2-bearing fluid can form talc (along with carbonates) in ultramafic rocks without Si-metasomatism. Variations in the relative mobility of Si and Mg at the subduction zone interface produce variations in the overall solid volume change of mantle (expansion or contraction), the types of sheet silicates (talc versus chlorite), and the fluid budget (dehydration or hydration) during metasomatic reactions, which affects the pore fluid pressure, frictional strength of the subduction megathrust, and the location of seismicity around the mantle wedge corner.

Veins of aragonite in deep ocean trenches can record discharges of CO2-rich paleoseawater and suggest that exposed serpentinized forearcs can act as carbon reservoirs, according to residence time calculations based on geochemical analyses of mantle rocks from offshore Japan.

Carbonation of a serpentinized mantle wedge is associated with volume reduction and fluid over-pressurization which drives brittle fracture and seismicity, according to chemical analyses of carbonate veins from the Sanbagawa belt Japan and thermodynamic calculations.

Heterogeneous reactions are chemical reactions that occur at the interfaces of multiple phases, and often show a nonlinear dynamical behavior due to the effect of the time-variant surface area with complex reaction mechanisms. It is important to specify the kinetics of heterogeneous reactions in order to elucidate the microscopic elementary processes and predict the macroscopic future evolution of the system. In this study, we propose a data-driven method based on a sparse modeling algorithm and sequential Monte Carlo algorithm for simultaneously extracting substantial reaction terms and surface models from a number of candidates by using partial observation data. We introduce a sparse modeling approach with non-uniform sparsity levels in order to accurately estimate rate constants, and the sequential Monte Carlo algorithm is employed to estimate time courses of multi-dimensional hidden variables. The results estimated using the proposed method show that the rate constants of dissolution and precipitation reactions that are typical examples of surface heterogeneous reactions, necessary surface models, and reaction terms underlying observable data were successfully estimated from only observable temporal changes in the concentration of the dissolved intermediate products.

The article presents a study which explores the geological records of transient fluid drainage into the shallow mantle wedge. It mentions about constraining the duration and velocity of fluid flow through a shallow mantle wedge based on the analyses of vein networks consisting of high-temperature serpentine in hydrated ultramafic rocks from the Oman ophiolite.

Talc is widely distributed over the Earth's surface and is predicted to be formed in various tectonic settings. Talc is a very soft and anisotropic sheet silicate showing very low friction behavior. Therefore, the formation of talc is expected to weaken the strength of talc-bearing rocks and may be associated with the initiation of subduction, and with a decrease in the coupling coefficient resulting in aseismic movements along faults and shear zones within subduction zones. For these reasons, understanding the crystallographic preferred orientation (CPO) of talc is important to quantify the anisotropy and physical properties of the host rock. However, it is difficult to measure a significant number of talc crystal orientations and to evaluate the accuracy of the measurements using electron-backscattered diffraction (EBSD). Therefore, talc CPO has not been reported, and there is uncertainty regarding the estimation of the strength of deformed talc-bearing rocks. Using methods developed for antigorite, we report the first successful EBSD measurements of talc CPO from a talc schist formed due to Si-metasomatism of ultramafic rocks by subduction zone fluids. We used a combination of W-SEM and FE-SEM measurements to examine domains of various grain sizes of talc. In addition, we used TEM measurements to evaluate the accuracy of the EBSD measurements and discuss the results of talc CPO analysis. Talc CPO in the present study shows a strong concentration of the pole to the (001) plane normal to the foliation. The strongest concentration of the [100] direction is parallel to the lineation. The talc schist produces similar S-wave splitting and P- and S-wave anisotropy as antigorite schist in deeper domains, thus identifying talc-rich layers in subduction zones may require a combination of geophysical surveys, seismic observations, and anisotropy modeling. The presence of strong talc CPO in rocks comprising the slab-mantle interface boundary may promote spatial expansion of the slip area during earthquakes along the base of the mantle wedge. [ABSTRACT FROM AUTHOR]