S U M M A R YFrom the dispersion of Love waves, we infer models of shear wave velocity structure underneath Bucharest (Romania) at depths down to 2 km that can contribute to seismic hazard estimation. Waves from eight regional events recorded during 10 months with a network of 34 seismic broad-band stations of the URban Seismology (URS) experiment are used. Although these events provide poor azimuthal coverage the data reliably constrain a shear wave velocity model with an interface between the Neogene and the Cretaceous sediments that is dipping northwards towards the Carpathian mountains. Array processing techniques that account for non-uniform wave propagation are used to estimate the dispersion of structural phase velocity. From this, we infer subsurface structure at three different latitudes. The Neogene sediments are represented by a gradient layer with no significant lateral variation. Shear wave velocity increases from approximately 400 m s −1 near the surface to 1 km s −1 at 1 km depth and 5 km in the south, and to 1.35 km s −1 at 1.5 km depth and 5 km in the north from the centre of Bucharest, respectively. For the half-space representing the Cretaceous sediments, we obtain shear wave velocities of 2.7-2.9 km s −1 . The results are consistent with results from boreholes and shallow seismics for the near-surface structure and results from receiver function studies and crustal refraction seismic studies for the deeper structure. The details of the Neogene layer comprising a vertical gradient fill a gap in existing models of the subsurface structure of Bucharest and can contribute to modelling of seismic hazard for the city. Since the signalto-noise ratio restricted useful data to the frequency range from 90 to 290 mHz, the inversion could not constrain the near-surface velocity independently. Due to strong trade-off between near-surface velocity and depth of half-space, the latter had to be introduced as a priori data from previous studies.