Explanation of supplementary files.

1. Data S1 (e.g., 1997-2000.out)
  The estimation results of the strain rate fields for the 19 periods in Figure 2. We set Wt=12 for Phase 1 (1997-2000 to 2008-2011) and Phase 3 (2016-2019 to 2018-2021), and Wt=4 for Phase 2 (2012-2015 to 2015-2018). The first column of data represents the longitude of the grid point, the second does the latitude, the third does the optimized smoothing distance D [km], the fourth does the dilatational strain rate [nstrain/yr], the fifth does the maximum shear strain rate [nstrain/yr], and the sixth does the azimuth of the maximum contraction of strain rate tau_Hmax [degree]. 


2. Data S2 (e.g., 2005-2011diff.out)
  The comparative results of the strain rate and stress fields for the periods of 2005-2011 in Phase 1 and 2016-2019 in Phase 3. The first column of data represents the longitude of the grid point, the second does the latitude, the third does the optimized smoothing distance D [km], the fourth does the dilatational strain rate [nstrain/yr], the fifth does the maximum shear strain rate [nstrain/yr], the sixth does the azimuth of the maximum contraction of strain rate tau_Hmax [degree], the seventh does the azimuth of the maximum compressive stress S_Hmax [degree], the eighth does the number of earthquakes, and the ninth does the absolute difference between the azimuths of tau_Hmax and S_Hmax [degree].