Investigating solitary waves is critical for predicting the tsunami waves features in the coastal region. The Smooth Particle Hydrodynamic method (SPH) is the Lagrangian meshless method that is commonly used to simulate nonlinear waves and free surface problems. In this study, the SPH method is used to investigate the propagation of the breaking and non-breaking solitary waves over the linear and nonlinear sharp slopes sea bed. The presented SPH model is validated in comparing the experimental and another numerical model. The results show that the presented method prepares powerful tools to simulate the solitary wave propagation over the variable sea bed. Then, the transmitted solitary wave from the linear and nonlinear sharp slopes of the sea bed were compared. The results show that the solitary wave energy loss for waves propagating through the linear slopes is more than that of non-linear slopes. Moreover, the solitary wave transmission coefficient increases with increases in wave height in the linear and nonlinear slope. Also, the transmitted solitary wave height through the nonlinear slope is slightly more than the linear slope. Finally, the breaking Solitary waves when passing through the variable depth in the near shore with nonlinear slopes is investigated by comparing the wave free surface, streamlines, and orbital velocities before and after wave breaking. The results show that the solitary wave horizontal orbital velocity values after the wave breaking are increasing in comparison with before breaking.