In this paper, we present MUSIQUE, an algorithm designed to analyze seismic signals recorded by arrays of three-component seismic sensors. Using the original MUSIC algorithm [1], azimuth and slowness (or velocity) of incident waves are estimated. The quaternion-MUSIC algorithm [2, 3] is then used to characterize the polarization properties of the waves. In this way, Love and Rayleigh waves are distinguished and their respective properties retrieved. This allows the characterization of the seismic wave field and the soil structure underneath the seismic array, which are important parameters for the estimation of seismic hazards