Additional file 1 of Microtremor array surveys and development of the velocity model in the Hakodate Plain, Hokkaido, Japan

Additional file 1: Table S1. Search ranges of parameters in the GA search for each observation site. Fig. S1. Array configurations at HAK site in Hakodate city. Red, green, and blue symbols and lines correspond to L array, S array, and N array, respectively. The circles indicate the locations of sensors. The black inverted triangle represents K-NET HKD160 strong motion station. Fig. S2. Array configurations at MHR site in Mihara, Hakodate city. Red, green, and blue symbols and lines correspond to L array, S array, and N array, respectively. The circles indicate the locations of sensors. The black inverted triangle represents JMA D66 strong motion station. Fig. S3. Array configurations at KMI site in Kamiiso, Hokuto city. Red, green, and blue symbols and lines correspond to L array, S array, and N array, respectively. The circles indicate the locations of sensors. The black inverted triangle represents HKDP08 strong motion station. Fig. S4. Array configurations at OON site in Ono, Hokuto city. Red, green, and blue symbols and lines correspond to L array, S array, and N array, respectively. The circles indicate the locations of sensors. The black inverted triangle represents HKDP09 strong motion station. Fig. S5. Array configurations at NNE site in Nanae town. Red, green, and blue symbols and lines correspond to L array, S array, and N array, respectively. The circles indicate the locations of sensors. The black square represents the location of the hot-spring well (H480-017). Fig. S6. Photo of the huddle test on September 25, 2017, at Omori Park in Hakodate city. Eight blue cylindrical-shaped sensors are LE-3D/5s manufactured by Lennartz-Electronic (one was for spare, which was not used in this field work). Data loggers DATAMARK LS-8800 connected to LE-3D/5s are stored in orange boxes. Four duralumin boxes are portable all-in-one accelerographs SMAR-6A3P (three accelerometers and data logger LS-8800). Fig. S7. Results of the huddle test for LE-3D/5s sensors. (Top) Velocity waveforms of the vertical component observed by seven LE-3D/5s sensors (0.1–10 Hz). The origin time of the time axis is at 10:08:18.72 on September 25, 2017 (JST). (Bottom Left) Ensemble average of the power spectral densities of five selected time-windows of the vertical component for all seven sensors. (Bottom Right) Ensemble average of the coherences and phase differences of the vertical components for all possible combinations of the seven sensors. Fig. S8. Results of the huddle test for SMAR-6A3P seismographs. (Top) Acceleration waveforms of the vertical component observed by four SMAR-6A3P seismographs (0.1–10 Hz). The origin time of the time axis is at 10:04:56.16 on September 25, 2017 (JST). (Bottom Left) Ensemble average of the power spectral densities of five selected time-windows of the vertical component for all four sensors. (Bottom Right) Ensemble average of the coherences and phase differences of the vertical components for all possible combinations of the four sensors. Fig. S9. Comparison of the S-wave velocity model estimated by the microtremor array observation with the geological profile of nearby hot spring well (Geological Survey of Hokkaido 1995, 2004). The locations of the hot spring wells were shown in Fig. 1(a). H: Holocene, Fm: Fumizuki Formation, Tm: Tomikawa Formation, Tv: Togeshita pyroclastic rocks, Mh: Mohejigawa Formation, Td: Todagawa Formation, Hk: Hekirijigawa Formation, K: Kamiiso Group. Fig. S10. Goodness-of-Fit (GOF) scores evaluated for five frequency bands (0.05–0.1 Hz, 0.1–0.2 Hz, 0.2–0.5 Hz, 0.5–1 Hz, and 1–2 Hz) based on the GOF criterion proposed by Anderson (2004). The square, circle, and triangle correspond to NS, EW, and UD component, respectively. The number above each panel is the average of the GOF scores of all frequency bands and components. Red: This study, purple: J-SHIS V2 model.