We investigated the shear wave velocity structure beneath the Damavand volcano, Iran, using analysis of Rayleigh wave empirical Green's functions obtained from the cross-correlation of 1-year seismic ambient noise and 26 local earthquakes. After the observed phase, velocity dispersion curves were inverted to obtain phase velocity tomographic maps; the synthetic phase velocity dispersion curves were then estimated for each evenly spaced geographic grid point. Estimated synthetic phase velocities were inverted using the nonlinear damped least-squares inversion method to obtain a quasi-3D shear wave velocity for the study area. Analysis of obtained quasi-3D shear wave velocity model reveals the presence of three distinct low or high anomalies in the upper 5 km. A low-velocity layer (V (S) similar to 2.8 km/s) was observed in the upper similar to 2 km of the crust that includes the sedimentary sequence of carbonate, siliciclastic, and volcanic rocks. According to our results, there is an indication of a high-velocity layer (V (S) a parts per thousand yenaEuro parts per thousand 3.0 km/s) at the depth range of 2 to 5 km, which may indicate the presence of dense, cooled magma. There is a clear indication of intrusion of the low-velocity anomaly/body into the higher velocity layer, in the depth range of approximately 3.0 to 4.5 km. The corresponding low-velocity body is interpreted as a hot magma chamber associated with the young eruption of Damavand. The chamber is located southwest of the crater, which stretches from the west to beneath the Damavand cone.