Plate boundaries at shallow depth are complex and highly deformed zones showing structures from both distributed and localized deformation (Collot et al., 2011;Sibson, 2013). Depending on the heterogeneity of the rock assemblage (i.e., competent vs. incompetent materials and/or clay contents) within the shear zone and the strain rates, shearing can be localized along discrete faults during seismic slip, while aseismic slip produces a scaly fabric associated with distributed deformation (Fagereng & Sibson, 2010;Kirkpatrick et al., 2015). Deciphering the stress and strain distribution across plate boundary shear zones is critical to understand the physical processes involved in the nucleation and behavior of megathrust faults.As magnetic minerals are sensitive to stress regime and fluid-rock interactions, the investigation of the magnetic properties of rocks can be an effective tool in studying faulting processes at intraplate shear zones (Yang et al., 2020). Anisotropy of magnetic susceptibility (AMS) provides insights into the preferred orientation of mineral grains and qualitative relationships between petrofabrics and deformation intensity (