We investigated the effects of the propagation path and site amplification of shallow tremors along the Nankai Trough. Using far-field S-wave propagation from intraslab earthquake data, the amplification factors at the DONET stations were 5–40 times against an inland outcrop rock site. Thick (~5 km) sedimentary layers with VS of 0.6–2 km/s beneath DONET stations have been confirmed by seismological studies. To investigate the effects of thick sedimentary layers, we synthesized seismograms of shallow tremors and intraslab earthquakes at seafloor stations. The ratios of the maximum amplitudes from the synthetic intraslab seismograms between models with and without thick sedimentary layers were 1–2. This means that the estimated large amplifications are primarily controlled by thin lower-velocity (< 0.6 km/s) sediments just below the stations. Conversely, at near-source (≤ 20 km) distances, 1-order amplifications of seismic energies for a shallow tremor source can occur due to thick sedimentary layers. Multiple S-wave reflections between the seafloor and plate interface are contaminated in tremor envelopes; consequently, seismic energy and duration are overestimated. If a shallow tremor occurs within underthrust sediments, the overestimation becomes stronger because of the invalid rigidity assumptions around the source region. After 1-order corrections of seismic energies of shallow tremors along the Nankai Trough, the scaled energies of seismic slow earthquakes were 10-10–10-9 irrespective of the region and source depth. Hence, the physical mechanisms governing seismic slow earthquakes can be the same, irrespective of the region and source depth.

We investigated the effects of the propagation path and site amplification of shallow tremors along the Nankai Trough. Using far-field S-wave propagation from intraslab earthquake data, the amplification factors at the DONET stations were 5–40 times against an inland outcrop rock site. Thick (~5 km) sedimentary layers with VS of 0.6–2 km/s beneath DONET stations have been confirmed by seismological studies. To investigate the effects of thick sedimentary layers, we synthesized seismograms of shallow tremors and intraslab earthquakes at seafloor stations. The ratios of the maximum amplitudes from the synthetic intraslab seismograms between models with and without thick sedimentary layers were 1–2. This means that the estimated large amplifications are primarily controlled by thin lower-velocity (< 0.6 km/s) sediments just below the stations. Conversely, at near-source (≤ 20 km) distances, 1-order amplifications of seismic energies for a shallow tremor source can occur due to thick sedimentary layers. Multiple S-wave reflections between the seafloor and plate interface are contaminated in tremor envelopes; consequently, seismic energy and duration are overestimated. If a shallow tremor occurs within underthrust sediments, the overestimation becomes stronger because of the invalid rigidity assumptions around the source region. After 1-order corrections of seismic energies of shallow tremors along the Nankai Trough, the scaled energies of seismic slow earthquakes were 10-10–10-9 irrespective of the region and source depth. Hence, the physical mechanisms governing seismic slow earthquakes can be the same, irrespective of the region and source depth.

We detected shallow very low frequency earthquakes (VLFEs) off the Cape Muroto and Kii Channel in the Nankai subduction zone and estimated their moment rate functions. Combining the new and previously estimated catalogs, we obtained the comprehensive catalog of shallow VLFE moment rate functions along the Nankai Trough. We defined the shallow VLFE swarms and investigated the scaling relationships of their cumulative moments, activity area, and durations in each region. Detected swarms were considered candidates for shallow slow slip events. A similar scaling relationship was observed between the cumulative moments and activity areas, irrespective of regions. It indicates similar stress drops in each region. However, the relationship between the cumulative moments and durations varied. This difference was explained by the along-strike variations in the faulting conditions of shallow slow earthquakes, such as material or hydrological properties.

Slow earthquakes follow a power-law size distribution with an exponential taper for the largest events. We investigated changes in the size distribution of shallow tectonic tremor events during two prolonged tremor episodes (>1 month) along the Nankai trough and found that the slope of the size distributions increased while the cut-off magnitudes decreased late during each episode, as tremor activity waned. Interpreting these changes with the two-dimensional probabilistic cell automaton model of slow earthquakes, we found that a decrease in event ignition probability or an increase in energy dissipation during slip can qualitatively explain the observed changes. These changes imply that a decrease in accumulated stress or pore-fluid pressure on the fault interface occurred during each tremor episode. Because the tremor source migrates during an episode, the changes in the size distribution parameters can be attributed to spatial variations or temporal changes in the source characteristics.

Fluid migration in subduction zones is a key controlling factor of slow and megathrust earthquakes at plate boundaries. During the migration, seismic velocity and heterogeneous structures in its pathways may be temporarily varied, preferably triggering slow earthquakes. Here, we show that transient changes of seismic heterogeneity occurred 0-9 months before shallow slow earthquakes in the Nankai subduction zone, Japan, using very long-term (6–10 y) records of ambient seafloor noise. The heterogeneity changes preceding to shallow slow earthquakes were observed near the margin of the source region, while concurrent changes primarily occurred in the source region. We propose that the heterogeneity changes are attributed to dynamic fluid migration, and the difference in timings reflects the pore pressure level in the corresponding source region. When fluids are supplied to a source region under relatively low pressure, fluids are leaked out from its downdip or updip side, and slow earthquakes occur not immediately but with a time delay of at most 9 months. In the high pore pressure case, slow earthquakes occur immediately with fluid migration from the source region. This study suggests that the heterogeneous seismic structure is possibly changed by fluid migration before slow earthquakes in the Nankai subduction zone.