<p>This study investigates the Love wave propagation in a multilayered geological structure consisting of an isotropic homogeneous layer overlying a fluid-saturated porous medium under an initial stress, which rests on a non-homogeneous elastic half-space with a parabolic interfacial irregularity. The governing equations of motion for each medium are formed using the theory of elasticity by Biot. Analytical solutions are obtained using Fourier transformation techniques. The dispersion relation for the Love waves is derived using the perturbation method and Willis’ integral formula, revealing the coupled influence of the material properties and the geometric parameters on wave characteristics. Computational simulations, conducted via MATLAB, demonstrate how the inhomogeneity parameter, anisotropy factor, porosity, initial stress and interface irregularity depth collectively affect phase velocity dispersion. The numerical results indicate that the phase velocity decreases with increasing wave number, inhomogeneity, porosity and anisotropy. A notable reversal in the phase velocity trend is observed around a critical wave number (~1000), which appears to depend on the depth-to-height ratio of the irregularity, in the case of the variation in inhomogeneity parameter (from 0 to 3). The results provide new insights into the Love wave behaviour in complex media, particularly highlighting the significant impact of small-scale interface imperfections that are often neglected in conventional models. These findings advance the theoretical understanding of seismic wave propagation in stratified porous systems and have practical implications in seismic hazard assessment and geotechnical engineering applications.</p>