Al-Shomrani, M.M.Yusuf, A.2026-03-152026-03-1520262473-698810.3934/math.20261672-s2.0-105029938452https://doi.org/10.3934/math.2026167https://hdl.handle.net/20.500.14517/8949The 2022–2023 global monkeypox (Mpox) virus outbreak, which affected more than 100 nations, has highlighted the importance of effective public health measures and predictive modeling techniques. In this work, we formulated a mathematical model (MD) of Mpox to assess the impact of early therapeutic education on seasonal variation. To guarantee epidemiological viability, we demonstrated the boundedness and positivity of the model. The next-generation matrix method was employed to determine the basic reproduction number and evaluate the stability of the disease-free equilibrium at the local and global levels. We assessed the existence of an arbitrary equilibrium. The center manifold theorem was employed to assess the backward bifurcation, which shows that imperfect vaccination causes the backward bifurcation. An evaluation of the existence of an endemic equilibrium point was also conducted. The global stability of an endemic equilibrium point was established through the application of a nonlinear Lyapunov function of the Goh-Voltterra type. Data fitting was done to show the model’s fitness and estimate the data’s parameters. Sensitivity analysis was ascertained. Seasonal variation was assessed. In a numerical simulation, we assessed five controls, which indicated that early therapeutic education is the most effective way to control the re-emergence of Mpox disease in the human (HUM) population. © 2026 the Author(s), licensee AIMS Press.eninfo:eu-repo/semantics/openAccessEducationMathematical ModelMonkeypoxParameter EstimationSeasonal VariationTherapeuticArticle