Most of shallow water geoacoustic inversions based on modal dispersion cannot reliably estimate the deep geoacoustic parameters. Because these studies focused on the dispersions of water waves but ignored the dispersions of ground waves. Therefore, in this paper a Bayesian geoacoustic inversion is studied based on wideband modal dispersions of water waves and ground waves. Firstly, the modal dispersion curves with Airy phase components are discussed. Secondly, the Bayesian inversion theory and a novel sample-efficient inference algorithm, namely Variational Bayesian Monte Carlo, are introduced briefly. In the Bayesian inversion, the posterior probability densities of unknown parameters are inferred, which can provide the prediction closest to the observation data and the uncertainty of the prediction. Considering that the forward acoustic model is computationally intensive, the posterior analysis is carried out by using the Variational Bayesian Monte Carlo method. It is performed by finding the variational distribution closest to the target distribution and requires less computation time than the Markov chain Monto Carlo method. In the simulation study, a range-independent two-layer seabed, including the sediment layer and basement layer, is modeled, on the assumption that the water column is homogeneous. The function of shear wave in waveguide is ignored. The compressional sound speed of the sediment layer varies linearly from c_1U to c_1L between 0 and h₁, while other geoacoustic parameters are constant. By comparing the inversion results with and without the information of ground waves for different signal-to-noise ratios, it can be concluded that the deep geoacoustic parameters are more sensitive to the dispersions of ground waves. And then, a shallow-water experimental study is carried out in the Bohai Sea of China. The average water depth is about 20 m. The wideband pulse signals are recorded by a hydrophone with a sensitivity of –170 dB re 1 V/μPa. The received signals include well-defined Airy phase components, and the modal dispersion curves of water waves and ground waves are extracted accurately. The experimental results indicate that the Bayesian inversion combining water and ground wave dispersions can not only estimate the deep geoacoustic parameters reliably, but also reduce the inversion uncertainties of other model parameters, such as the shallow geoacoustic parameters, water depth, etc. The estimated source-receiver range and water sound speed are close to their measured values. The modal dispersion curves predicted by the posterior mean samples are in good consistence with those extracted from the experimental data in different ranges. In addition, the good forecast of transmission loss also demonstrates the reliability of the joint inversion.