The 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.
Citation: Mohammed M Al-Shomrani, Abdullahi Yusuf. Mathematical model of Mpox: assessing the impact of early therapeutic education on seasonal variation[J]. AIMS Mathematics, 2026, 11(2): 4147-4199. doi: 10.3934/math.2026167
The 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.
| [1] |
Md. A. I. Islam, M. H. M. Mubassir, A. K. Paul, S. S. Shanta, A mathematical model for understanding and controlling Mpox transmission dynamics in the USA and its implications for future epidemic management, Decoding Infection and Transmission, 2 (2024), 100031. https://doi.org/10.1016/j.dcit.2024.100031 doi: 10.1016/j.dcit.2024.100031
|
| [2] |
J. Kongson, C. Thaiprayoon, W. Sudsutad, Analysis of a mathematical model for the spreading of the Mpox virus with a constant proportional-Caputo derivative operator, AIMS Mathematics, 10 (2025), 4000–4039. https://doi.org/10.3934/math.2025187 doi: 10.3934/math.2025187
|
| [3] |
S. Anil, B. Joseph, M. Thomas, V. K. Sweety, N. Suresh, T. Waltimo, Monkeypox: a viral zoonotic disease of rising global concern, Infectious Diseases & Immunity, 4 (2024), 121–131. https://doi.org/10.1097/ID9.0000000000000124 doi: 10.1097/ID9.0000000000000124
|
| [4] |
J. Molla, I. Sekkak, A. M. Ortiz, I. Moyles, B. Nasri, Mathematical modeling of Mpox: a scoping review, One Health, 16 (2023), 100540, https://doi.org/10.1016/j.onehlt.2023.100540 doi: 10.1016/j.onehlt.2023.100540
|
| [5] |
E. M. Bunge, B. Hoet, L. Chen, F. Lienert, H. Weidenthaler, L. R. Baer, et al., The changing epidemiology of human Mpox—A potential threat? A systematic review, PLoS Negl. Trop. Dis., 16 (2022), e0010141. https://doi.org/10.1371/journal.pntd.0010141 doi: 10.1371/journal.pntd.0010141
|
| [6] | Centres for Disease Control and Prevention, About Monkeypox, 2023. Available from: https://www.cdc.gov/monkeypox/about/index.html. |
| [7] |
S. T. Rehan, H. U. Hussain, K. A. Kumar, M. Sukaina, Z. Khan, A. J. Nashwan, Global Mpox virus outbreak 2022: a bibliometric analysis, Cureus, 15 (2023), e37107. https://doi.org/10.7759/cureus.37107 doi: 10.7759/cureus.37107
|
| [8] | World Health Organization, Mpox (Mpox) outbreak 2022-Global, 2022. Available from: https://www.who.int/emergencies/situations/monkeypox-oubreak-2022/1000. |
| [9] |
A. W. Rimoin, N. Kisalu, B. Kebela-Ilunga, T. Mukaba, L. L. Wright, P. Formenty, et al., Endemic human monkey pox, Democratic Republic of Congo, 2001–2004, Emerg. Infect. Dis., 13 (2007), 934–937. https://doi.org/10.3201/eid1306.061540 doi: 10.3201/eid1306.061540
|
| [10] | Centers for Disease Control and Prevention, Morbidity and Mortality Weekly Report, Volume 52, November 7, 2003, No. 44. Available from: https://stacks.cdc.gov/view/cdc/28471/cdc_28471_DS1.pdf. |
| [11] |
O. I. Idisi, T. T. Yusuf, E. Adeniyi, A. A. Onifade, Y. T. Oyebo, A. T. Samuel, et al., A new compartmentalized epidemic model to analytically study the impact of awareness on the control and mitigation of the monkeypox disease, Healthcare Analytics, 4 (2023), 100267. https://doi.org/10.1016/j.health.2023.100267 doi: 10.1016/j.health.2023.100267
|
| [12] |
B. Cabanillas, G. Murdaca, A. Guemari, M. J. Torres, A. K. Azkur, E. Aksoy, et al., Mpox 2024 outbreak: fifty essential questions and answers, Allergy, 79 (2024), 3285–3309. https://doi.org/10.1111/all.16374 doi: 10.1111/all.16374
|
| [13] | D. N. Tanko, F. A. Abdullah, M. K. M. Ali, M. O. Adewole, J. Andrawus, Analyzing onchocerciasis transmission dynamics through mathematical models incorporating early therapeutic intervention and vector control, IAENG International Journal of Applied Mathematics, 55 (2025), 3579–3599. |
| [14] |
D. N. Tanko, F. A. Abdullah, M. K. M. Ali, M. O. Adewole, J. Andrawus, Onchocerciasis control dynamics: a focus on early treatment and vector management strategies, Journal of the Nigerian Society of Physical Sciences, 8 (2026), 2942. https://doi.org/10.46481/jnsps.2026.2942 doi: 10.46481/jnsps.2026.2942
|
| [15] |
O. J. Peter, S. Kumar, N. Kumari, F. A. Oguntolu, K. Oshinubi, R. Musa, Transmission dynamics of Mpox virus: a mathematical modelling approach, Model. Earth Syst. Environ., 8 (2022), 3423–3434. https://doi.org/10.1007/s40808-021-01313-2 doi: 10.1007/s40808-021-01313-2
|
| [16] |
B. Liu, S. Farid, S. Ullah, M. Altanji, R. Nawaz, S. W. Teklu, Mathematical assessment of Mpox disease with the impact of vaccination using a epidemiological modeling approach, Sci. Rep., 13 (2023), 13550. https://doi.org/10.1038/s41598-023-40745-x doi: 10.1038/s41598-023-40745-x
|
| [17] |
M. I. Betti, L. Farrell, J. Heffernan, A pair formation model with recovery: application to Mpox, Epidemics, 44 (2023), 100693. https://doi.org/10.1016/j.epidem.2023.100693 doi: 10.1016/j.epidem.2023.100693
|
| [18] |
T. D. Alharbi, M. R. Hasan, Modeling Mpox dynamics with human–rodent interactions and waning vaccination, AIMS Math., 10 (2025), 18660–18679. https://doi.org/10.3934/math.2025834 doi: 10.3934/math.2025834
|
| [19] |
H. Iftikhar, M. Daniyal, M. Qureshi, K. Tawaiah, R. K. Ansah, J. K. Afriyie, A hybrid forecasting technique for infection and death from the Mpox virus, Digit. Health, 9 (2023), 1–17. https://doi.org/10.1177/20552076231204748 doi: 10.1177/20552076231204748
|
| [20] |
A. Elsonbaty, W. Adel, A. Aldurayhim, A. El-Mesady, Mathematical modeling and analysis of a novel monkeypox virus spread integrating imperfect vaccination and nonlinear incidence rates, Ain Shams Eng. J., 15 (2024), 102451. https://doi.org/10.1016/j.asej.2023.102451 doi: 10.1016/j.asej.2023.102451
|
| [21] |
Y. Yang, C. Li, X. Fan, W. Long, Y. Hu, Y. Wang, et al., Effectiveness of omadacycline in a patient with chlamydia psittaci and KPC-producing gram-negative bacteria infection, Infect. Drug Resist., 18 (2025), 903–908. https://doi.org/10.2147/IDR.S505311 doi: 10.2147/IDR.S505311
|
| [22] |
Y. Wang, Q. Wang, T. W. Yang, J. M. Yin, F. Wei, H. Liu, et al., Analysis of immune and inflammatory microenvironment characteristics of noncancer end-stage liver disease, J. Interf. Cytok. Res., 43 (2023), 86–97. https://doi.org/10.1089/jir.2022.0172 doi: 10.1089/jir.2022.0172
|
| [23] |
L. Tang, Y. Chen, Q. Xiang, J. Xiang, Y. Tang, J. Li, The association between IL18, FOXP3 and IL13 genes polymorphisms and risk of allergic rhinitis: a meta-analysis, Inflamm. Res., 69 (2020), 911–923. https://doi.org/10.1007/s00011-020-01368-4 doi: 10.1007/s00011-020-01368-4
|
| [24] |
N. Ndembi, P. Mbala-Kingebeni, B. Khalifa, M. L. Mazaba, M. O. Foláyan, Are there ecological and seasonal factors influencing the resurgence of Mpox in Africa?, J. Public Health Afr., 15 (2024), a823. https://doi.org/10.4102/jphia.v15i1.823 doi: 10.4102/jphia.v15i1.823
|
| [25] |
M. Karuhanga, V. Yiga, A mathematical model for the dynamics of onchocerciasis with vector control and mass disease administration, J. Appl. Math., 2024 (2024), 6683371. https://doi.org/10.1155/2024/6683371 doi: 10.1155/2024/6683371
|
| [26] | U. E. Michael, L. O. Omenyi, E. Kafayat, E. Nwaeze, O. A. Akachukwu, G. Ozoigbo, et al., Monkeypox mathematical model with surveillance as control, Commun. Math. Biol. Neurosci., 2023 (2023), 6. |
| [27] |
S. Wang, K. Zhang, A. Liu, Flat-Lattice-CNN: a model for Chinese medical-named-entity recognition, PLoS One, 20 (2025), e0331464. https://doi.org/10.1371/journal.pone.0331464 doi: 10.1371/journal.pone.0331464
|
| [28] |
W. Wang, J.-X. Li, S.-Q. Long, Z.-N. Liu, X.-P. Li, Z.-H. Peng, et al., Emerging strategies for monkeypox: antigen and antibody applications in diagnostics, vaccines, and treatments, Military Med. Res., 12 (2025), 69. https://doi.org/10.1186/s40779-025-00660-w doi: 10.1186/s40779-025-00660-w
|
| [29] |
J. Andrawus, Y. U. Ahmad, A. A. Andrew, A. Yusuf, S. Qureshe, B. A. Denue, et al., Impact of surveillance in human-to-human transmission of Mpox virus, Eur. Phys. J. Spec. Top., 234 (2024), 483–514. https://doi.org/10.1140/epjs/s11734-024-01346-5 doi: 10.1140/epjs/s11734-024-01346-5
|
| [30] |
Y. U. Ahmad, J. Andrawus, A. Ado, Y. A. Maigoro, A. Yusuf, S. Althobaiti, et al., Mathematical modeling and analysis of human-to-human Mpox virus transmission with post-exposure vaccination, Model. Earth Syst. Environ., 10 (2024), 2711–2731. https://doi.org/10.1007/s40808-023-01920-1 doi: 10.1007/s40808-023-01920-1
|
| [31] |
G. D. Yaleu, S. Bowong, E. H. Danga, J. Kurths, Mathematical analysis of the dynamical transmission of Neisseria meningitidis serogroup A, Int. J. Comput. Math., 94 (2017), 2409–2434. http://doi.org/10.1080/00207160.2017.1283411 doi: 10.1080/00207160.2017.1283411
|
| [32] |
O. Diekmann, J. A. P. Heesterbeek, J. A. J. Metz, On the definition and the computation of the basic reproduction ratio $\mathcal{R}_c$ in models for infectious diseases in heterogeneous populations, J. Math. Biol., 28 (1990), 365–382. https://doi.org/10.1007/BF00178324 doi: 10.1007/BF00178324
|
| [33] |
P. Van De Driessche, J. Watmough, Reproduction number sub-threshold Endemic Equilibrium for compartmental models of disease transmission, Math. Biosci., 180 (2002), 29–48. https://doi.org/10.1016/S0025-5564(02)00108-6 doi: 10.1016/S0025-5564(02)00108-6
|
| [34] |
J. Andrawus, S. Abdulrahman, R. V. K. Singh, S. S. Manga, Sensitivity analysis of mathematical modeling of Ebola virus population dynamics in the presence of vaccine, Dutse Journal of Pure and Applied Sciences, 8 (2022), 40–46. https://doi.org/10.4314/dujopas.v8i2a.5 doi: 10.4314/dujopas.v8i2a.5
|
| [35] |
J. Andrawus, A. Yusuf, U. T. Mustapha, A. S. Alshomrani, D. Baleanu, Unravelling the dynamics of Ebola virus with contact tracing as control strategy, Fractals, 31 (2023), 2340159. https://doi.org/10.1142/S0218348X2340159X doi: 10.1142/S0218348X2340159X
|
| [36] |
C. Castillo-Chavez, B. Song, Dynamical models of tuberculosis and their applications, Mathematical Biosciences and Engineering, 1 (2004), 361–404. https://doi.org/10.3934/mbe.2004.1.361 doi: 10.3934/mbe.2004.1.361
|
| [37] | C. Castillo-Chavez, Z. Feng, W. Huang, On the computation of $R_0$ and its role on global stability, In: Mathematical approaches for emerging and reemerging infectious diseases: an introduction, New York: Springer, 2002,229–250. |
| [38] |
C. Castillo-Charvez, B. Song, Dynamical model of tuberculosis and their applications, Mathematical Biosciences and Engineering, 1 (2004), 361–404. https://doi.org/10.3934/mbe.2004.1.361 doi: 10.3934/mbe.2004.1.361
|
| [39] |
F. Y. Eguda, J. Andrawus, S. Babuba, The solution of a mathematical model for Dengue fever transmission using the differential method, Journal of the Nigerian Society of Physical Sciences, 1 (2019), 82–87. https://doi.org/10.46481/jnsps.2019.18 doi: 10.46481/jnsps.2019.18
|
| [40] | J. P. La Salle, The stability of dynamical systems, Philadelphia: SIAM, 1976. https://doi.org/10.1137/1.9781611970432 |
| [41] | J. Andrawus, E. Y. Felix, Mathematical analysis of a model for syphilis endemicity, International Journal of Scientific Engineering and Applied Science, 3 (2017), 48–72. |