Outputs of VSI-unstable disk and α-viscous disk (α = 5 × 10⁻⁴) simulations with embedded massive planets at the disk midplane (Z = 0). From left to right: gas density normalized by the initial value, radial velocity (υ_r), meridional velocity (υ_θ), and azimuthal velocity deviations from Keplerian rotation (υ_ϕ − υ_Kep). From top to bottom: VSI simulation with a Saturn-mass planet (m_p = 0.3m_Jup), α-viscous simulation with a Saturn-mass planet, VSI simulation with a Jupiter-mass planet (m_p = m_Jup), and α-viscous simulation with a Jupiter-mass planet. The black dashed lines in the first panel show the simulation buffer zones. The black dotted lines show the planets’ orbital radius. The white circles show the location of the planets. The diameters of the circles are equal to the planets’ Hill spheres.