## 4 Approach of the Algorithm with the Kik-Net Vertical Arrays Data
\begin{table}
\begin{tabular}{c c c c c c c c} \hline \hline  & \multirow{2}{*}{**Source: [26]**} & \multicolumn{2}{c}{**Calculated**} & \multirow{2}{*}{**Calculated**} & \multirow{2}{*}{**According to**} & \multirow{2}{*}{**Calculated According to**[39]} & \multirow{2}{*}{**Source: [43, 44]**} & \multirow{2}{*}{**Calculated by Equations (22) and (25)**} \\  & & & [38] & & & & & \\ \hline
**Depth to Top** & **V\({}_{\mathbf{\
u}}\)** & **\(\rho_{\text{bulk}}=\text{f(V,)}\)**,** & **\(\rho_{\text{bulk}}=\text{f(V,)}\)**,** & **\(\rho_{\text{avg}}=\text{f(V,)}\)**,** & **\(\rho_{\text{s}}\)**,** & **V\({}_{\mathbf{\
u}}=\text{f(P, $\rho_{\text{avg}}$)}\)** & **\(\Phi=\text{f(P, $\rho_{\text{avg}}$)}\)** & **\(\Phi=\text{f(P, $\rho_{\text{avg}}$)}\)** \\
**of Layer (m)** & **(m/s)** & **(m/s)** & **(kg/m\({}^{3}\))** & **(kg/m\({}^{3}\))** & **(kg/m\({}^{3}\))** & **(kg/m\({}^{3}\))** & **(kg/m\({}^{3}\))** & **V\({}_{\mathbf{\
u}}=\text{f(P, $\rho_{\text{avg}}$)}\)** & **\(\Phi=\text{f(P, $\rho_{\text{avg}}$)}\)** \\ \hline
0 & 220 & 90 & 1921 & 1339 & 1246 & 2.74 & 0.07 & 0.55 \\
1 & 590 & 130 & 1525 & 1475 & 1338 & 2.74 & 0.10 & 0.51 \\ \hline
5 & 1500 & 210 & 1926 & 1972 & 1468 & 2.74 & 0.14 & 0.46 \\
16 & 1500 & 300 & 1926 & 1836 & 1572 & 2.74 & 0.17 & 0.43 \\ \hline
36\({}^{*}\) & 3100 & 1400 & 2310 & 2749 & 2117 & 2.63 & 0.30 & 0.20 \\ \hline \hline \end{tabular} \({}^{*}\)Bedrock.  
\end{table} 
Table 2: Geotechnical parameters (initial and calculated) at the KSRH10 clayey site of the Kik-net network. Figure 8: Degradation curves of shear modulus \(G\)/\(G_{max}\) (\(\gamma\)) (**a**), constrained moduli above water table _M/M\({}_{max}\)_ (\(\gamma\)) (**b**) and under water table _M\({}^{*}\)/M\({}^{*}_{max}\)_ (\(\gamma\)) (**c**) and corresponding damping ratio curves \(D(\gamma)\) (**c**) used for modeling the responses for TCGH15 and KSRH10 Kik-net sites. Figure 10: Comparison of response spectra (damping 5%) of modeled and recorded accelerograms (components EW, NS, UD; spectral acceleration—SA, spectral velocity—SV, spectral displacement—SD) calculated for the surface of soil layers at Kik-net site KSRH10. Figure 9: Comparison of response spectra (damping 5%) of modeled and recorded accelerograms (components EW, NS, UD; spectral acceleration—SA, spectral velocity—SV, spectral displacement—SD) calculated for the surface of soil layers at Kik-net site TCGH15. Figure 11 shows the residuals of response spectra of modeled and recorded accelerators dramatically for the surface of soil layers at Kik-net sites TCGH15 and KSRH10. The obtained residuals for the TCGH15 site turned out to be on average smaller than for the KSRH10 site. However, this does not apply to the vertical components. Moreover, the residuals for the vertical component for the KSRH10 site are smaller than for the horizontal components, especially in the spectral range \(T<0.1\) s. For the TCGH15 site, the residuals of the horizontal and vertical components are mostly similar. It should also be noted that, for spectral displacements, the residuals noticeably increase in the spectral range \(T>3\) for both sites.