The low-frequency vector wave field modeling method for the predic-tion of fields in layered inhomogeneous ocean and the ocean bottom en-vironment is proposed. The method is especially stable to vertically cut media multiple arbitrary thickness layers. Based on integral presentation of 2-D cylindrically symmetric source wave fields in elastic media, the method accounts for all wave types involved. Inhomogeneous media in the problem are cut in N horizontally homogeneous layers, for which 4 (N−1) equations were derived, accounting for conditions on both layers boundaries. Equation factors found by the Schmidt global matrix method are introduced in Fourier–Bessel integral expressions for local field parameters. Expressions were calculated numerically to obtain values of acoustic pressure and elastic stress, vertical and horizontal particle velocity component in liquid and elastic media, respectively. Deviations fromexact solutions expressed in modeling stability losses were obliged mainly to layers’ excessive thickness choice. A simple enough ‘‘benchmark’’ problem solution, where results for the conventional and the proposed method were compared to the exact solution, demonstrates method advan-tages. Point source of frequency 0.01–10.0 Hz, situated in a water-layer field model is presented. The influence of the bottom elastic layers’ structure on model parameters is demonstrated. The geological–acoustic bottom model depending on frequency is proposed.