Paramushir is the northernmost island of the Kuril Arc comprising several Holocene volcanoes, of which two are presently active: Ebeko and Chikurachki. We present the first crustal-scale three-dimensional seismic shear-wave velocity model of the northern part of Paramushir derived with the use of data of 20 temporary stations operated in 2021-2022 and one permanent station. The continuous seismic data were used to obtain the Rayleigh wave dispersion curves for periods ranging from 0.5 s to 12 s through computing ambient seismic noise cross- correlation functions. Then the 3D shear wave velocity distribution was constructed by seismic tomography. The synthetic tests demonstrate that the model has sufficient resolution down to 7-10 km depth. We observe a series of low-velocity anomalies at the depths of 4-6 km below the volcano-origin Vernadsky Ridge. They may likely represent magma reservoirs that are responsible for Holocene eruptions along the ridge. These anomalies are overlain by the high-velocity anomalies associated with the rigid cover consisting of an upper-crustal granite- granodiorite body and consolidated magma intrusions. Along the eastern flank of Ebeko we reveal a thin low- velocity anomaly at a depth of ~1 km that can be interpreted as a water-saturated layer. The contact of this layer with hot magma intrusions leads to the formation of a large amount of steam that is ejected during the frequent phreatic eruptions of Ebeko and from numerous fumaroles in the summit area. On the other hand, this layer, which is following down to the Pacific Coast, might be promising for the geothermal energy exploitation