At plate boundaries, interseismic surface strain can be interpreted in terms of fault coupling (i.e., the extent of kinematic locking of a fault) and allows to infer where a fault is prone to release elastic energy through slip, either seismically or aseismically (e.g., Avouac, 2015;Bürgmann, 2018). However, such interpretation does not account for temporal variations in slip rate, though frequently observed in various seismo-tectonic contexts and at different spatial and temporal scales (Jolivet & Frank, 2020, and references therein). Long-term (>1 year) variations of interseismic coupling have been detected in Sumatra, Parkfield and along the Cascadia subduction zone (Barbot et al., 2013;Materna et al., 2019;Tsang et al., 2015). On shorter time scale (<1 year), slow slip events (SSEs) have been observed in areas previously thought to be either locked or slipping continuously, both in subduction zones (e.g.,