In a 2007 experiment conducted in the northern North Sea, observations of a low-frequency seismo-acoustic wave field with a linear horizontal array of vector sensors located on the seafloor revealed a strong, narrow peak around 38 Hz in the power spectra and a presence of multi-mode horizontally and vertically polarized interface waves with phase speeds between 45 and 350 m/s. Dispersion curves of the interface waves exhibit piece-wise linear dependences between the logarithm of phase speed and logarithm of frequency with distinct slopes at large and small phase speeds, which suggests a seabed with a power-law shear speed dependence in two distinct sediment layers. The power spectrum peak is interpreted as a manifestation of a seismo-acoustic resonance. A simple geoacoustic model with a few free parameters is derived that quantitatively reproduces the key features of the observations. This article's approach to the inverse problem is guided by a theoretical analysis of interface wave dispersion and resonance reflection of compressional waves in soft marine sediments containing two or more layers of different composition. Combining data from various channels of the vector sensors is critical for separating waves of different polarizations and helps to identify various arrivals, check consistency of inversions, and evaluate sediment density.