The transfer of material from the downgoing plate to the overriding plate at depth exerts a first-order control on the mechanical and thermal evolution of subduction zones. However, the geometry, size, and temporal organization of deep tectonic slices formed during basal accretion remain poorly constrained, due to the limited resolution of geophysical imaging and the rarity of geological analogues preserving deep accretionary architectures formed during continental subduction.

Here, we investigate the spatial extent and stacking dynamics of tectonic slices by reconstructing the architecture of the deep paleo-accretionary wedge through the study of the exhumed Phyllite-Quartzite (PQ) nappe, derived from microcontinental protoliths, which was formed by basal accretion during the Oligo-Miocene along the Hellenic subduction zone. We carried out a multidisciplinary study of this now-exhumed PQ nappe, which crops out discontinuously from Crete to the Peloponnese (Greece). Preserved in a fore-arc position and weakly overprinted by later tectonic events, this natural laboratory provides direct access to deep accretionary processes.

An integrated petro-structural study conducted across southeastern Peloponnese and Kythira combines detailed mapping, structural analysis, petrological observations, Raman spectroscopy of carbonaceous material, and thermobarometric modeling. This approach allows us to distinguish several tectono-metamorphic sub-units within the PQ nappe stack, each recording a distinct P-T evolution that constrains the depth of basal accretion for successive episodes. Hypotheses of lateral continuity between these sub-units provide first-order constraints on their present-day spatial extent and on the minimum size of individual accretionary slices.

In southeastern Peloponnese, two HP-LT sub-units are identified within the PQ nappe stack, while at least two equivalent sub-units are recognized on Kythira. These sub-units record a systematic increase in peak temperature from the base to the top of the HP-LT nappe stack, consistent with successive episodes of basal accretion. Reconstructed P-T conditions indicate that basal accretion occurred at depths of ~50-60 km along the subduction interface. Based on spatial correlations between structurally equivalent HP–LT sub-units exposed in neighboring regions along strike, we infer a minimum present-day lateral continuity of individual accretionary slices. On this basis, deep tectonic slices formed during basal accretion are inferred to currently extend over several tens of kilometers in the trench-perpendicular direction and up to a hundred kilometers along strike.

This study provides new quantitative constraints on the depth, lateral extent, and dynamics of tectonic underplating, with direct relevance for the Hellenic margin, where such processes may still be active, and for active subduction zones worldwide.

How to cite: Bouhot, M., Menant, A., Ganino, C., Angiboust, S., Oncken, O., Deldicque, D., Jolivet, L., and Skarpelis, N.: Spatial extent of deep slab slicing events:  Insights from the Phyllite-Quartzite paleo-accretionary wedge (Hellenic subduction zone), EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-6595, https://doi.org/10.5194/egusphere-egu26-6595, 2026.

Subducting bathymetric reliefs, such as seamounts, modify the slip behavior of megathrusts, thereby potentially dictating seismic segmentation, rupture dynamics, and the structural evolution of the subduction channel and upper plate. While geodetic data often suggest that the megathrust near subducting seamounts is weakly coupled and dominated by aseismic creep or microseismicity, several "seamount earthquakes" have been documented. The role of subducting topography in governing fault coupling, rupture dynamics, and the spatial distribution of rupture remains poorly understood.

Laboratory seismotectonic experiments provide an effective means of simulating earthquake cycles and observing fault slip behaviour with high spatiotemporal precision, thereby overcoming the limitations of sparse onshore and missing offshore geodetic networks, as well as short historical records. In our experiments, a topographic high with seamount geometry was subducted along a 15° dipping, velocity-weakening seismogenic zone accompanied by hundreds of analogue earthquake cycles. The model upper plate is a wedge composed of an elastoplastic granular material that can respond to seismic cycles and seamount-induced stresses. We constrained the interface slip distribution by combining analogue geodetic slip inversion of surface displacement with direct monitoring of the interface via side-view imaging.

The results reveal that during the early stages of seamount subduction, when the seamount has partially subducted beneath the upper plate, along-strike rupture propagation is arrested at the seamount, which acts as a barrier, producing partial ruptures. Progressively, as the main portion or the entire seamount becomes subducted, another consistent spatial pattern emerges: coseismic slip concentrates at the leading downdip edge of the seamount, while the center and updip regions remain largely aseismic, with minor shallow slip reflecting slope instabilities triggered by upper-plate extensional structures. This pattern aligns well with interseismic high-coupling patches, which can also extend to the deep flank of the seamount.

Our findings indicate that, while subducting seamounts inhibit earthquake nucleation and broadly arrest rupture propagation, they still allow slip to extend onto the seamount-bearing interface. This explains why the deeper flank of a subducting seamount or ridge remains seismically active. A series of earthquakes (1996 Mw 6.7 and Mw 6.8; 2024 Mw 7.1; 2025 Mw 6.8) systematically occurred around the downdip edge of a Kyushu-Palau Ridge. Similar rupture behavior has been documented for a series of Mw ~ 7 events in the southern Japan Trench and for the two Mw > 8 events in central Nankai. This spatial pattern is further supported by geological evidence of pseudotachylytes, which are only localized on the downdip side of the exhumed fossil seamount.

Beyond slip kinematics, our experiments demonstrate that subducting seamounts perturb the megathrust stress field, leading to heterogeneous stress accumulation along dip, consistent with previous numerical mechanical-hydrological modeling studies. This suggests that seamount-induced coupling enhances upper-plate deformation and long-term structural features, including forearc uplift, fault reactivation, and localized fracturing. The short- and long-term upper-plate deformation patterns provide a key means of identifying subducted topographic features and assessing their impact on earthquake and tsunami hazards.

How to cite: Tan, H., Kosari, E., Rosenau, M., Gao, X., and Oncken, O.: Spatial Patterns of Megathrust Seismogenic Behavior Modulated by a Subducting Seamount, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-9476, https://doi.org/10.5194/egusphere-egu26-9476, 2026.

Current studies on seamount subduction propose contrasting effects on megathrust behavior. Some suggest that subducting seamounts increase normal stress and promote large earthquakes, while others argue that seamounts fracture the upper plate, enhancing microseismicity and aseismic creep that may inhibit major ruptures. Field- and microscale observations from the Azuero Accretionary Complex, an exhumed intra-oceanic accretionary complex in Panama, provide new constraints on the deformation processes associated with seamount subduction.

Coastal exposures on the Azuero Peninsula expose the contact between the autochthonous Azuero Plateau and the allochthonous Azuero Accretionary Complex. The Azuero Plateau forms part of the Caribbean Large Igneous Province and consists mainly of massive to pillowed oceanic plateau basalts with minor Upper Cretaceous chert. In contrast, the accretionary complex comprises massive to pillowed basalts and volcanic breccias with ocean island affinity, locally interbedded with Paleogene carbonates. A ~3 km wide deformation zone, the Azuero Mélange, separates these units, and is inferred to be a deformed portion of the plateau based on new field observations and geochemical data.

The accreted seamount lithologies show pervasive faulting, cataclasites, abundant zeolite veins, and a chlorite-rich shear zone located ~60 m below the mélange. These rocks lack evidence for large displacement through-going faults. In contrast, plateau-derived rocks record both brittle faulting and ductile deformation. Ductile strain localized within the Azuero Mélange, where clay-rich cataclasites accommodated deformation through cataclastic flow and dissolution–precipitation creep. At the structural base of the mélange, a ~10 m thick shear zone composed of foliated cataclasites with basalt and limestone clasts within a clay-rich matrix is observed, interpreted to be sheared seamount lithologies.  

Fluids and alteration played a crucial role in localizing strain within the upper plate and the décollement, enhancing mechanical weakening and diffusive mass transfer. Cataclasis increased permeability, enabling fluid infiltration and the formation of mechanically weak phases such as clays and chlorite. These processes promoted strain localization and facilitated deformation by cataclastic flow and dissolution-precipitation creep. The interplay between these alteration and deformation processes likely favored aseismic creep during seamount subduction.

How to cite: Maglalang, E. J., Fagereng, A., Buchs, D., and Toffol, G.: Deformation and alteration during seamount subduction: Insights from an exhumed intra-oceanic accretionary complex, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-14561, https://doi.org/10.5194/egusphere-egu26-14561, 2026.

Understanding the structure of the incoming Pacific Plate is essential for reconstructing the tectonic evolution of the Japan Trench. International Ocean Discovery Program Expedition 405, “JTRACK,” drilled through the frontal prism and plate boundary fault zone of the Japan Trench and into the underlying Pacific Plate at Site-C0019, located close to the hypocentre of the 2011 Mw 9.0 Tohoku‑Oki earthquake. Here, we present new petrographic and geochemical results from five basaltic samples recovered from below the plate boundary fault zone (930–946 mbsf), which show the first direct geological evidence for a subducted volcanic body beneath this segment of the margin.

Core and image logs from Expedition 405 indicate that the drilled interval comprises a repeated sequence of sheeted dykes, massive basalts, and pillow lavas, with at least one interbedded sedimentary horizon separating the successions. This architecture is incompatible with simple ocean-plate stratigraphy and instead indicates a later phase of volcanic activity interacting with pre-existing crust. Of the five samples analysed, one was taken from below the intercalated sedimentary horizon and four from above it. The four samples from above the sediments are systematically more evolved than the sample below. In the samples above, plagioclase shows anorthitic cores overprinted by albitised rims, pyroxenes are Mg-rich, and ilmenite is abundant. Several samples contain K-rich clays, disseminated sulphides, and Zr-rich domains, including possible baddeleyite. These features suggest interaction with sedimentary cover and progressive evolution of magma chemistry during the later stages of emplacement, alongside post-emplacement hydrothermal alteration. Whole-rock major and trace element data show enrichment in incompatible elements (K, Pb, Cs, Rb, Ba, Th, U) relative to typical MORB. K-metasomatism discrimination diagrams indicate that these signatures are not produced by alteration, supporting a primary magmatic origin. Together, the mineralogical and chemical characteristics point to a volcanic body formed off-axis on the Pacific Plate rather than at a spreading ridge.

Integration with seismic interpretations suggests that this volcanic body, likely a seamount chain, was faulted prior to subduction and now lies partially perched on a horst, forming a structural high directly beneath the plate boundary fault zone. These results provide the first physical confirmation of a subducted, faulted seamount chain beneath the Japan Trench, offering new constraints on the structure and evolution of the incoming Pacific Plate.

How to cite: Gough, A., Jurado, M.-J., Ishikawa, T., Fukuchi, R., Webb, M., Nakamura, Y., Yamaguchi, A., Conin, M., Nicholson, U., Gürer, D., Rasbury, T., Fulton, P., Kirkpatrick, J., Kodaira, S., Regalla, C., Ujiie, K., Eguchi, N., Maeda, L., Okutsu, N., and Toczko, S. and the Expedition 405 Scientists: Petrological and Geochemical Evidence for a Subducted Off-Axis Seamount Chain Beneath the Japan Trench Plate Boundary Fault Zone, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-9894, https://doi.org/10.5194/egusphere-egu26-9894, 2026.

Accurately assessing strain accumulation and release in subduction zones is contingent upon robust detection and characterization of locking and slip along the megathrust. However, the distribution of slip on shallow, offshore plate boundaries is not well-resolved with onshore GNSS networks. At the Hikurangi Subduction Zone offshore Aotearoa-New Zealand, extensive investment has been made into seafloor geodetic techniques such as seafloor pressure and GNSS-acoustic, which have significantly improved observation and characterization of offshore SSEs. Despite their utility, oceanographic noise limits the ability of these seafloor techniques to detect SSEs. Formation pore pressure changes (as a proxy for volumetric strain) detected in borehole observatories have an enhanced signal-to-noise ratio and can reliably resolve deformation at the 10s of nanostrain-level, providing an improved view of shallow crustal deformation offshore.

Here, we report on a suite of SSEs observed in two IODP borehole observatories in the northern Hikurangi Subduction Zone between 2018 and 2023 and model their slip distribution and magnitude. During this time, five SSEs were clearly recorded in the borehole pore pressure data. Four of these occurred spontaneously, and the borehole pressure changes correlate with surface displacement observed at onshore GNSS stations. In contrast, in early 2021, the M_w 7.2 East Cape earthquake triggered a near-trench SSE that was only captured by the observatories. We jointly invert changes in pore pressure with onshore GNSS displacements and seafloor pressure (when available) for slip distribution along a 2D transect for each of the events. Our inversions incorporate realistic elastic properties constrained by high-resolution seismic velocity models and logging-while-drilling data, which is crucial for accurately resolving slip distribution and magnitude. We find large differences in slip initiation and evolution characteristics during the 2021 triggered SSE compared to the spontaneous events. We also find that, in total, SSEs accommodate most (>80%) of the plate convergence budget along the shallow (<10 km) megathrust. The 2021 triggered event was particularly important for filling in a slip deficit near the trench. Our results have implications for the role of SSEs in accommodating the megathrust strain budget near the trench at subduction zones.

How to cite: Carlson, M., Wallace, L., Saffer, D., and Williams, C.: Slow Slip Accommodates the Full Plate Convergence Budget at the Northern Hikurangi Subduction Zone, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-8339, https://doi.org/10.5194/egusphere-egu26-8339, 2026.

Sulawesi and Borneo are tectonically complex islands in South East Asia that were assembled from numerous crustal fragments throughout the Cenozoic. Past P wave tomography studies of this region have used land stations and hence image seismic velocity structure primarily beneath the two islands. By adding data from Ocean Bottom Seismometers deployed between 2019 and 2020 in the Makassar Strait, which lies between Sulawesi and Borneo, and incorporating core-going converted P phases with steeper ray paths (PKP and PKIKP), we illuminate the mantle beneath this offshore region to depths of around 800 km. The new tomographic velocity model robustly images a high-velocity north-west dipping tabular anomaly at depths between ~300 and ~660 km beneath the strait, which is interpreted as an aseismic ancient slab. By combining a slab age-depth relationship with a plate tectonic reconstruction, we uncover the palaeosubduction boundary responsible for this slab, thereby providing conclusive evidence for the previously hypothesised north-westward subduction of the Celebes Sea beneath Northern Borneo around ~15 Ma (plus/minus a few million years). Following slab break-off and northward plate migration of Sulawesi, sinking of this northwest Celebes Sea slab may have contributed to the initiation of subduction of the Celebes Sea southwards beneath Northern Sulawesi, which today is confidently imaged by a Benioff zone in addition to seismic tomography.

How to cite: Rawlinson, N., Li, Y., Pilia, S., Kesumastuti, L., Lü, C., Widiyantoro, S., and Hao, T.: Cenozoic Subduction Polarity Reversal Within the Celebes Sea Inferred from Teleseismic Tomography, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-13387, https://doi.org/10.5194/egusphere-egu26-13387, 2026.

Slip behavior on the shallowest part of the subduction zone is heavily influenced by the lithostratigraphy of the incoming plate. Coseismic slip that propagates toward the seafloor is a major tsunami hazard as seen in the 2011 Tohoku-Oki earthquake. In Nankai Trough, scientific drilling combined with seismic facies mapping has allowed scientists to characterize the lithostratigraphy of the Shikoku Basin which in turn have been linked with the geotechnical properties, structural architecture, and seismicity. Seismic reflection surveys have been conducted here for more than 25 years covering over 730 x 150 km and mapping the sediment input using the traditional visual interpretation at this scale is inefficient. Taking advantage of this large seismic dataset combined with scientific drilling to map the fault properties can help improve hazard assessments. This study presents a framework for a semi-automatic classification using 2D seismic reflection profiles acquired by different surveys. Three trench-parallel lines within the accretionary wedge acquired by three different surveys were used for the case study. The amplitudes of the western and eastern profiles were first scaled based on the seafloor reflection of the central profile. Features were extracted using a rectangular window measuring 500 m above the top of the oceanic basement with varying widths measuring 1 km, 3 km, and 5 km. A total of 121 features divided into three groups were extracted. The statistical group (13) describes the strength of reflections, the spectral group (30) describes the presence or absence of internal reflections, and textural group (78) describes the continuity of reflections. The three principal components of each group were extracted and altogether subjected to K-Means clustering with 6 clusters. The 5 km window showed the most comparable classification with visual interpretation and the consistent classification in the overlap between profiles indicate a satisfactory performance of our method. Comparing the classification with previous drilling, Cluster 0 located in the overlap between the central and eastern profile is associated with turbidites occurring in basement lows. Cluster 1 is classified as noise. Cluster 2 in the western and central profile are likely siliciclastic turbidites from the Kyushu Fan. Cluster 3 and 4 in the eastern profile appears to have no drilling analog. Cluster 5 in the central profile is associated with hemipelagic mudstones. These initial results appear promising and will be tested in larger datasets and other subduction zones in the future.

How to cite: Flores, P. C., Fujie, G., Shiraishi, K., Nakamura, Y., Kimura, G., Su, J., Agata, R., and Kodaira, S.: An unsupervised classification scheme for seismic facies mapping of sediment input in Nankai Trough using reflection amplitudes from 2D profiles, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-11054, https://doi.org/10.5194/egusphere-egu26-11054, 2026.

Numerous seismic surveys have been conducted to uncover the geophysical structure and fluid distribution in subduction zones, since seismic velocities are primarily influenced by pore fluids within rocks. To link seismic velocity and the pore fluid state, laboratory measurements of elastic wave velocity using rocks samples from subduction zone environments have been conducted, revealing the role of microscopic pores and cracks in controlling seismic velocity. However, natural geological systems are heterogeneous and contain defects of different scales at each scale of observation. Therefore, the velocities observed at in situ scales should be affected not only by microscopic pore structures but also by larger-scale defects, such as fractures and faults. Such large-scale defects should also play a role in fluid drainage system, since permeability of rocks depends strongly on the dimensions of conduits. In this study, we compare laboratory-measured ultrasonic velocity measured on core samples from the Susaki area in the Shimanto accretionary complex, SW Japan, with sonic velocity measured by borehole logging experiments. P-wave velocities were measured at a frequency of 1 MHz under dry conditions at 5 cm intervals along core sections spanning a total length of 128 m. The measured values were then converted to velocities under wet conditions using an effective medium model, enabling comparison with sonic velocities acquired under groundwater-saturated conditions. Results show that P-wave velocity decreases from the laboratory (~6 km/s) to the borehole scales (~5 km/s). This scale-variant effect can be explained by an effective medium model whereby mesoscale porosity that is undetectable at the ultrasonic wavelength is introduced into the matrix phase with microscale porosity. Assuming typical apertures for micro- and mesoscale fractures, we estimate that the effective permeability can increase to 10^–12–10^–11 m² with increasing in the mesoscale porosity and decreasing P-wave velocity down to 4–5 km/s. These results indicate that seismic velocity anomalies and related seismic activity are associated with the presence of mesoscale fractures in subduction zones.

How to cite: Akamatsu, Y., Okuda, H., Kitamura, M., and Sawai, M.: Scale-dependent seismic velocity and permeability in subduction zones caused by mesoscale fractures, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-4431, https://doi.org/10.5194/egusphere-egu26-4431, 2026.

Understanding how deformation is distributed within accretionary prisms is essential to constrain their structural evolution, internal organisation and seismo-tsunamigenic potential. Researchers would typically use seismic images to characterise large-scale subsurface structure, but accretionary prisms often appear internally chaotic-to-transparent in such data, lacking coherent reflectors. This is likely due to strong lateral heterogeneity and stratal disruption around the scale of the seismic resolution, resulting from intense deformation of accreted sediments and remobilised mass transport deposits.

IODP Expedition 405 “JTRACK” (October-December 2024) drilled the frontal prism of the Japan Trench subduction zone at Site C0019, targeting the decollement that hosted extremely large slip (>50 m) to the trench during the 2011 Mw 9.1 Tohoku-Oki earthquake. The expedition involved continuous coring and logging-while-drilling (LWD) through the prism, resulting in an improved time-depth tie with regional 2-D seismic profiles. At this site the seismic images show a highly chaotic internal prism, which introduces large uncertainties in structural interpretation and inhibits the upscaling and lateral extrapolation of borehole-scale observations from Site C0019.

Here we present an approach to characterise this heterogeneous internal structure by inverting for scale-independent geostatistical parameters (vertical and lateral correlation lengths, dip) from the combined seismic image and LWD sonic velocity data. We use a Bayesian Markov Chain Monte Carlo approach that gives spatially varying, probabilistic estimates of the geostatistical parameters. The lateral correlation length, specifically, can be used as a proxy for the degree of stratal disruption/deformation, as shorter correlation lengths reflect more intense deformation. While the aspect ratio of the correlation lengths is generally well-constrained for seismic data, the estimate of vertical correlation length from the borehole logs is necessary to resolve the other individual geostatistical parameters away from the borehole.

Initial results reveal that the degree of deformation varies significantly within the frontal prism, with the lateral correlation length varying between around 10-50 m. The degree of deformation appears to be compartmentalised by steeply dipping structures that could correspond to fault zones also observed in cores from C0019. These distinct changes in lateral correlation length correspond to lithological units interpreted from core observations. They also coincide with significant changes in vertical correlation length from the sonic log, anisotropy of seismic velocity from core samples and bedding orientation from borehole images. The results demonstrate that seismically-derived geostatistical parameters can delineate internal compartmentalisation of the frontal prism, providing a framework for tectonic and mass transport deposit interpretation and for the extrapolation of core-scale observations. Future work will extend this analysis using parallel profiles along the Japan Trench axis, allowing for mapping of frontal prism internal deformation in three dimensions away from existing drill sites.

How to cite: Ford, J., Nakamura, Y., Jurado, M.-J., Schaible, K., Nicholson, U., Doan, M.-L., Pei, P., Hamada, Y., Miyakawa, A., Conin, M., Fulton, P., Kirkpatrick, J., Kodaira, S., Regalla, C., Ujiie, K., Eguchi, N., Maeda, L., Okutsu, N., and Toczko, S. and the Expedition 405 Scientists: Geostatistical characterisation of internal deformation of the Japan Trench frontal prism using seismic and logging-while-drilling data (Site C0019), EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-12829, https://doi.org/10.5194/egusphere-egu26-12829, 2026.

Making a link between hydraulic properties and lithology is important for understanding fluid flow and pore pressure evolution in subduction zone sediments. However, it's a difficult challenge to get the in-situ downhole hydraulic responses due to the complex lithological variations and strong vertical heterogeneity.  IODP Expedition 405 - JTRACK, drilled Site C0026 in the Japan Trench, penetrating a thick sedimentary sequence by hemipelagic mud, pelagic clay, chert and basalt intervals. Continuous logging-while-drilling (LWD) and coring operations provide a unique opportunity to examine lithology-related hydraulic properties in the sedimentary structures.

In this study, we integrate the interpretated downhole annular pressure (DHAP) from LWD time series dataset and other physical properties, including gamma ray, resistivity, sonic velocity, and caliper logs, together with interpretated lithogical logging units. Comparing forward-modelled DHAP with measured DHAP data, the results indicates that the in-situ fluid pressure evolution is highly correlated with lithological variability. In the hemipelagic mud at shallow depth, apparent inflow is largely influenced by borehole enlargement, where caliper increases strongly affect flow modelling and make it difficult to identify the real formation inflow. In the pelagic clay interval, little to no fluid inflow or loss is observed, indicating the absence of significant overpressure or under pressure. Localized fluid loss is identified within chert layers, consistent with their brittle behaviour and increased caliper values, suggesting fracture-controlled fluid escape. In contrast, a clear and active inflow is observed in the basalt interval, where fractures identified from borehole image logs provide efficient pathways for fluid flow to the borehole.

These observations indicate that hydraulic properties at Site C0026 are strongly controlled by lithological contrasts, leading to vertically variations fluid flow behavior. This study highlights the importance of integrating with borehole logs lithological information to constrain fluid transport processes in subduction zone sedimentary sequences.

How to cite: Pei, P. and Doan, M.-L. and the Expedition 405 Scientists: Lithology-related hydraulic properties of subduction zone sediments at Japan Trench, IODP Exp.405 Site C0026, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-13415, https://doi.org/10.5194/egusphere-egu26-13415, 2026.

The rheological and compositional evolution of a subducting plate critically depends on its thermal structure. The temperature evolution of the subducting plate depends on its history prior to subduction, on its interaction with the overriding plate, and on how it interacts with the ambient mantle. Many processes that are associated with subduction such as deep seismicity and fluid release, which are responsible for arc-volcanism, can be understood through the temperature evolution of the slab. Studying the temperature evolution of the subducting slab, however, is not straightforward because of the lack of direct observations and a complete subduction record. Thus, to assist the interpretation of the available data it is necessary to use forward geodynamic modelling (FWG).

FWG can either study a full dynamic system or study the evolution of the slab using a kinematic model, in which the slab geometry and velocity are fully prescribed as boundary conditions. Kinematic models are more suitable to study specific subduction zones, as the dynamic models would require high computational cost to fit the available data. However, kinematic models suffer from several assumptions that oversimplify the complexity of the subduction process. For example, kinematic models are frequently performed with constant convergence velocity, age of subducting plate, and thermal properties. 

In this contribution, we aim to test whether the pressure- and temperature-dependent thermal properties improve our ability to interpret the natural data, and if the improvements are worth the additional complexity.  

How to cite: Piccolo, A., Craig, T. J., Van Zelst, I., and Thieulot, C.: The effect of pressure-temperature dependent material properties on thermal evolution of the slab, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-14028, https://doi.org/10.5194/egusphere-egu26-14028, 2026.

The Sierra del Convento mélange in southeastern Cuba is one of a limited number of jadeite jade occurrences worldwide. This exhumed, pervasively serpentinized fossil subduction interface hosts tectonic blocks up to several tens of meters in scale that were reworked during intense metasomatism. The flat-lying mélange spans a 300 m thick vertical exposure, progressing from a block-rich lower portion toward a block free, serpentinite-dominated upper region. We conducted structurally-controlled sampling of blocks and matrix to produce a spatially-resolved geochemical and microstructural profile through the thickness of the Convento mélange. Pale green, nearly pure jadeitite constitutes the dominant block population in the southern part of this mélange, and may contain up to 30% epidote and/or white mica by volume. Minor metasedimentary and metamafic block populations, recording variable degrees of HP metamorphism, coexist alongside the jade blocks. We report a newly identified metre-scale zoning within the jadeitite bodies, consistent with that documented in other jade localities. Green jadeitite locally occurs surrounding an older core of zoisitite containing relict jadeite crystals rimmed by omphacite. Green jade is cross-cut by fractures infilled by dark-colored brecciated jade, which is in turn rimmed by a late, pyroxene-free rind composed mainly of weakly foliated phengite + albite ± epidote. The latter facies occupies a similar position to chloritite blackwalls previously described from the Convento jade occurrence. All jade varieties except for these latest phengite-albite rinds and chloritite blackwalls are cross-cut by fractures infilled by jadeite and omphacite. Parts of the main jade bodies exhibit prismatic radial and comb jadeite microstructures, consistent with descriptions of P-type jade, which precipitate directly into open fractures from hydrothermal fluids. However, Convento jade contains paragonite with up to 2 wt.% K₂O and jadeite-omphacite exsolution domains brecciated into jigsaw-like fragments recemented by jadeite and/or omphacite. These observations are consistent with at least part of the jade in the Sierra del Convento mélange representing near-total high-temperature metasomatic replacement of high-pressure anatectic trondhjemite protoliths, which originated as partial melts of garnet amphibolite at ~15 kbar, resembling an R-type (replacement) jade paragenesis. To overcome the considerable ambiguity inherent to geochronology datasets from HP igneous and metasomatic rocks from this locality (105 to 115 Ma, U-Pb zircon), we are conducting detailed multi-mineral and multi-system geochronology, including further U-Pb on zircon, titanite, and apatite, Ar/Ar on white mica, and Rb-Sr on white mica. This multi-chronometric approach will establish relative and numerical chronology for the diverse jade facies of the Convento occurrence, resolving timescales of the multiple associated fluidization events within the subduction channel.

How to cite: Ducharme, T., Angiboust, S., Cambeses, A., Peverelli, V., Raimondo, T., Núñez-Cambra, K., Blanco-Quintero, I. F., Cárdenas-Párraga, J., and Garcia-Casco, A.: Timescales of metasomatism in a hot subduction channel: a microstructural and (radio-)isotopic study of Sierra del Convento jade, EGU General Assembly 2026, Vienna, Austria, 3–8 May 2026, EGU26-16928, https://doi.org/10.5194/egusphere-egu26-16928, 2026.