The present study is dedicated to investigate the SH body-as well as Love-waves propagation effects in porous media with uncertain porosity and permeability. A unified formulation of the governing equations for one-dimensional (1-D) wave propagation in anisotropic porous layered media is presented deterministically. The uncertainties around the above two cited parameters are taken into account by random fields with the help of Monte Carlo Simulations (MCS). Random samples of the porosity and the permeability are generated according to the normal and lognormal distribution functions, respectively, with a mean value and a coefficient of variation for each one of the two parameters. After performing several thousands of samples, the mathematical expectation (mean) of the solution of the wave propagation equations in terms of amplification functions for SH waves and in terms of dispersion equation for Love-waves are obtained. The limits of the Love wave velocity in a porous soil layer overlaying a homogeneous half-space are obtained where it is found that random variations of porosity change the zeros of the wave equation. Also, the increase of uncertainties in the porosity (high coefficient of variation) decreases the mean amplification function amplitudes and shifts the fundamental frequencies. However, no effects are observed on both Love wave dispersion and amplification function for random variations of permeability. Lastly, the present approach is applied to a case study in the Adapazari town basin so that to estimate ground motion accelerations lacked in the fast-growing during the main shock of the damaging 1999 Kocaeli earthquake.