In this paper, a deterministic and stochastic fractional-order model of the tri-trophic food chain model incorporating harvesting is proposed and analysed. The interaction between prey, middle predator and top predator population is investigated. In order to clarify the characteristics of the proposed model, the analysis of existence, uniqueness, non-negativity and boundedness of the solutions of the proposed model are examined. Some sufficient conditions that ensure the local and global stability of equilibrium points are obtained. By using stability analysis of the fractional-order system, it is proved that if the basic reproduction number R 0 < 1, the predator free equilibrium point E 1 is globally asymptotically stable. The occurrence of local bifurcation near the equilibrium points is investigated with the help of Sotomayor's theorem. Some numerical examples are given to illustrate the theoretical findings. The impact of harvesting on prey and the middle predator is studied. We conclude that harvesting parameters can control the dynamics of the middle predator. A numerical approximation method is developed for the proposed stochastic fractional-order model.