S U M M A R YIn recent years, there has been increasing interest in the study of seismic noise interferometry as it can provide a complementary approach to active source or earthquake-based methods for imaging and continuous monitoring the shallow structure of the Earth. This meaningful information is extracted from wavefields propagating between those receiver positions at which seismic noise was recorded. Until recently, noise-based imaging relied mostly on Rayleigh waves. However, considering similar wavelengths, a combined use of Rayleigh and Love wave tomography can succeed in retrieving velocity heterogeneities at depth due to their different sensitivity kernels. Here, we present a novel one-step algorithm for simultaneously inverting Rayleigh and Love wave dispersion data aiming at identifying and describing complex 3-D velocity structures. The algorithm may help to accurately and efficiently map the shear wave velocities and the Poisson ratio of the surficial soil layers. In the high-frequency range, the scattered part of the correlation functions stabilizes sufficiently fast to provide a reliable estimate of the velocity structure not only for imaging purposes but also allows for changes in the medium properties to be monitored. Such monitoring can be achieved with a high spatial resolution in 3-D and with a time resolution as small as a few hours. In this paper, we describe a recent array experiment in a volcanic environment in Solfatara (Italy) and we show that this novel approach has identified strong velocity variations at the interface between liquids and gas-dominated reservoirs, allowing localizing a region which is highly dynamic due to the interaction between the deep convection and its surroundings.