Spatiotemporal distribution and seismic interaction of very-low-frequency earthquakes in the northern Ryukyu Trench

Slow earthquakes play a crucial role in understanding stress accumulation and release along plate interfaces in subduction zones. The northern Ryukyu Trench, where the Philippine Sea Plate subducts northwestward beneath the Eurasian Plate, experienced a major earthquake in 1911 and is currently regarded as a low-seismicity area (LSA). Understanding the seismic activity in this region, particularly the relationship between very-low-frequency earthquakes (VLFEs) and regular seismic events, is crucial for understanding subduction zone dynamics. We investigated the spatial and temporal distribution of VLFE activity in the northern Ryukyu Trench using broadband ocean-bottom seismometers deployed around Amami Island between September 2018 and June 2019. Our analysis, employing the envelope correlation method, revealed that VLFE activity is primarily concentrated northeast of Amami Island, an area characterized by low regular earthquake activity, with the distribution of VLFEs spatially segregated from that of regular earthquakes. Furthermore, we observed earthquake swarm activity at the edges of the LSA in the northern Ryukyu Trench following VLFE activity. In November 2018, intense VLFE activity northeast of Amami Island migrated northeastward, which was followed by a regular earthquake swarm at the edge of this LSA. Following VLFE activity in January 2019, additional seismic activity, including foreshocks, occurred at the edges of this LSA approximately 1 month later. The spatial segregation of VLFEs and regular earthquakes suggests that VLFE activity may be influenced by the migration of high-pressure fluids within the subducted slab. This migration appears to trigger related time-delayed seismic activity, similar to mechanisms observed in other subduction zones such as Hikurangi. Understanding these dynamics is essential for assessing the coupling state of subduction zones and associated fluid behaviors, which play a critical role in evaluating seismic hazards in LSAs.