This paper estimates seabed shear-wave velocity profiles and their uncertainties using interface-wave dispersion curves extracted from data generated by a shear-wave source. The shear-wave source generated a seismic signature over a frequency range between 2 and 60 Hz and was polarized in both in-line and cross-line orientations. Low-frequency Scholte- and Love-waves were recorded. Dispersion curves of the Scholte- and Love-waves for the fundamental mode and higher-order modes are extracted by three time-frequency analysis methods. Both the vertically and horizontally polarized shear-wave velocity profiles in the sediment are estimated by the Scholte- and Love-wave dispersion curves, respectively. A Bayesian approach is utilized for the inversion. Differential evolution, a global search algorithm is applied to estimate the most-probable shear-velocity models. Marginal posterior probability profiles are computed by Metropolis-Hastings sampling. The estimated vertically and horizontally polarized shear-wave velocity profiles fit well with the core and in situ measurements.