This study investigates the threshold behavior of a population-varying stochastic susceptible-infectious-recovered-susceptible (SIRS) model driven by a logarithmic Ornstein-Uhlenbeck process. Introducing the logarithmic Ornstein-Uhlenbeck process to account for random environmental fluctuations enhances the biological significance of the model. By applying the Itô stochastic integral, we construct a suitable Lyapunov function, proving the existence and uniqueness of the global positive solution of the model, thereby ensuring biological feasibility. Then, a critical threshold parameter $ R_0^s $ is derived: If $ R_0^s > 1 $, the system admits a unique invariant probability measure; if $ R_0^s < 1 $, the infection dies out almost surely around the disease-free equilibrium. Furthermore, near the quasi-equilibrium point, the invariant probability density admits a local Gaussian approximation and converges weakly to a normal distribution as the environmental noise tends to zero. Numerical simulations in MATLAB illustrate the theoretical results, further reveal the sensitivity of the stochastic threshold to the key parameters, and confirm that population variation influences the threshold structure and long-term infection level.
Citation: Yanan Zhao, Chang Liu. Threshold behavior in a stochastic SIRS epidemic model with a logarithmic Ornstein-Uhlenbeck process[J]. Electronic Research Archive, 2026, 34(6): 4290-4324. doi: 10.3934/era.2026191
This study investigates the threshold behavior of a population-varying stochastic susceptible-infectious-recovered-susceptible (SIRS) model driven by a logarithmic Ornstein-Uhlenbeck process. Introducing the logarithmic Ornstein-Uhlenbeck process to account for random environmental fluctuations enhances the biological significance of the model. By applying the Itô stochastic integral, we construct a suitable Lyapunov function, proving the existence and uniqueness of the global positive solution of the model, thereby ensuring biological feasibility. Then, a critical threshold parameter $ R_0^s $ is derived: If $ R_0^s > 1 $, the system admits a unique invariant probability measure; if $ R_0^s < 1 $, the infection dies out almost surely around the disease-free equilibrium. Furthermore, near the quasi-equilibrium point, the invariant probability density admits a local Gaussian approximation and converges weakly to a normal distribution as the environmental noise tends to zero. Numerical simulations in MATLAB illustrate the theoretical results, further reveal the sensitivity of the stochastic threshold to the key parameters, and confirm that population variation influences the threshold structure and long-term infection level.
| [1] |
W. O. Kermack, A. G. McKendrick, A contribution to the mathematical theory of epidemics, Proc. R. Soc. London, Ser. A, 115 (1927), 700–721. https://doi.org/10.1098/rspa.1927.0118 doi: 10.1098/rspa.1927.0118
|
| [2] |
S. Busenberg, P. Van den Driessche, Analysis of a disease transmission model in a population with varying size, J. Math. Biol., 28 (1990), 257–270. https://doi.org/10.1007/BF00178776 doi: 10.1007/BF00178776
|
| [3] |
R. Rifhat, Z. Teng, C. Wang, Extinction and persistence of a stochastic SIRV epidemic model with nonlinear incidence rate, Adv. Differ. Equations, 2021 (2021), 200. https://doi.org/10.1186/s13662-021-03347-3 doi: 10.1186/s13662-021-03347-3
|
| [4] |
X. Bao, X. Han, Dynamic behavior and numerical simulation of stochastic epidemic models with time delay, Discrete Contin. Dyn. Syst. Ser. S, 18 (2025), 1923–1937. https://doi.org/10.3934/dcdss.2024209 doi: 10.3934/dcdss.2024209
|
| [5] |
J. Ge, W. Zuo, D. Jiang, Stationary distribution and density function analysis of a stochastic epidemic HBV model, Math. Comput. Simul., 191 (2022), 232–255. https://doi.org/10.1016/j.matcom.2021.08.003 doi: 10.1016/j.matcom.2021.08.003
|
| [6] |
S. Hussain, O. Tunç, G. ur Rahman, H. Khan, E. Nadia, Mathematical analysis of stochastic epidemic model of MERS-corona & application of ergodic theory, Math. Comput. Simul., 207 (2023), 130–150. https://doi.org/10.1016/j.matcom.2022.12.023 doi: 10.1016/j.matcom.2022.12.023
|
| [7] |
Z. Shi, D. Jiang, Dynamical behaviors of a stochastic HTLV-I infection model with general infection form and Ornstein–Uhlenbeck process, Chaos, Solitons Fractals, 165 (2022), 112789. https://doi.org/10.1016/j.chaos.2022.112789 doi: 10.1016/j.chaos.2022.112789
|
| [8] |
Y. Zhou, D. Jiang, Dynamical behavior of a stochastic SIQR epidemic model with Ornstein–Uhlenbeck process and standard incidence rate after dimensionality reduction, Commun. Nonlinear Sci. Numer. Simul., 116 (2023), 106878. https://doi.org/10.1016/j.cnsns.2022.106878 doi: 10.1016/j.cnsns.2022.106878
|
| [9] |
J. Shang, W. Li, Dynamical behaviors of a stochastic SIRV epidemic model with the Ornstein–Uhlenbeck process, Adv. Contin. Discrete Models, 2024 (2024), 9. https://doi.org/10.1186/s13662-024-03807-6 doi: 10.1186/s13662-024-03807-6
|
| [10] |
W. Wang, Y. Cai, Z. Ding, Z. Gui, A stochastic differential equation SIS epidemic model incorporating Ornstein–Uhlenbeck process, Physica A, 509 (2018), 921–936. https://doi.org/10.1016/j.physa.2018.06.099 doi: 10.1016/j.physa.2018.06.099
|
| [11] |
K. Mamis, M. Farazmand, Stochastic compartmental models of COVID-19 pandemic must have temporally correlated uncertainties, Proc. R. Soc. A, 479 (2023), 20220568. https://doi.org/10.1098/rspa.2022.0568 doi: 10.1098/rspa.2022.0568
|
| [12] |
B. Zhou, D. Jiang, B. Han, T. Hayat, Threshold dynamics and density function of a stochastic epidemic model with media coverage and mean-reverting Ornstein–Uhlenbeck process, Math. Comput. Simul., 196 (2022), 15–44. https://doi.org/10.1016/j.matcom.2022.01.014 doi: 10.1016/j.matcom.2022.01.014
|
| [13] |
Z. Shi, D. Jiang, Stochastic modeling of SIS epidemics with logarithmic Ornstein–Uhlenbeck process and generalized nonlinear incidence, Math. Biosci., 365 (2023), 109083. https://doi.org/10.1016/j.mbs.2023.109083 doi: 10.1016/j.mbs.2023.109083
|
| [14] |
B. Han, D. Jiang, Complete characterization of dynamical behavior of stochastic epidemic model motivated by black-Karasinski process: COVID-19 infection as a case, J. Franklin Inst., 360 (2023), 14841–14877. https://doi.org/10.1016/j.jfranklin.2023.10.007 doi: 10.1016/j.jfranklin.2023.10.007
|
| [15] |
M. Gao, Y. Jiang, D. Jiang, Threshold dynamics of a stochastic SIRS epidemic model with transfer from infected individuals to susceptible individuals and log-normal Ornstein-Uhlenbeck process, Electron. Res. Arch., 33 (2025), 3037–3064. https://doi.org/10.3934/era.2025133 doi: 10.3934/era.2025133
|
| [16] | R. Khasminskii, Stochastic Stability of Differential Equations, Springer Science & Business Media, 2012. https://doi.org/10.1007/978-3-642-23280-0 |
| [17] |
Q. Liu, D. Jiang, T. Hayat, A. Alsaedi, B. Ahmad, A stochastic SIRS epidemic model with logistic growth and general nonlinear incidence rate, Physica A, 551 (2020), 124152. https://doi.org/10.1016/j.physa.2020.124152 doi: 10.1016/j.physa.2020.124152
|
| [18] |
R. Rudnicki, Long-time behaviour of a stochastic prey-predator model, Stochastic Processes Appl., 108 (2003), 93–107. https://doi.org/10.1016/S0304-4149(03)00090-5 doi: 10.1016/S0304-4149(03)00090-5
|
| [19] |
B. Zhou, X. Zhang, D. Jiang, Dynamics and density function analysis of a stochastic SVI epidemic model with half saturated incidence rate, Chaos, Solitons Fractals, 137 (2020), 109865. https://doi.org/10.1016/j.chaos.2020.109865 doi: 10.1016/j.chaos.2020.109865
|
| [20] |
B. Zhou, H. Wang, T. Wang, D. Jiang, Stochastic generalized Kolmogorov systems with small diffusion: I. Explicit approximations for invariant probability density function, J. Differ. Equations, 382 (2024), 141–210. https://doi.org/10.1016/j.jde.2023.10.057 doi: 10.1016/j.jde.2023.10.057
|
| [21] | X. Mao, Stochastic Differential Equations and Applications, Elsevier, 2007. |
| [22] |
Y. Zhao, D. Jiang, X. Mao, A. Gray, The threshold of a stochastic SIRS epidemic model in a population with varying size, Discrete Contin. Dyn. Syst. - Ser. B, 20 (2015), 1277–1295. https://doi.org/10.3934/dcdsb.2015.20.1277 doi: 10.3934/dcdsb.2015.20.1277
|
| [23] | D. R. Bell, The Malliavin calculus, Courier Corporation, 2012. |
| [24] |
N. H. Du, D. H. Nguyen, G. G. Yin, Conditions for permanence and ergodicity of certain stochastic predator-prey models, J. Appl. Probab., 53 (2016), 187–202. https://doi.org/10.1017/jpr.2015.18 doi: 10.1017/jpr.2015.18
|
| [25] |
N. Etemadi, An elementary proof of the strong law of large numbers, Z. Wahrscheinlichkeitstheorie verw Gebiete, 55 (1981), 119–122. https://doi.org/10.1007/BF01013465 doi: 10.1007/BF01013465
|
| [26] | B. Øksendal, Stochastic Differential Equations: An Introduction with Applications, Springer Science & Business Media, 2013. https://doi.org/10.1007/978-3-642-14394-6 |
| [27] | Z. Ma, Y. Zhou, C. Li, Qualitative and Stability Methods for Ordinary Differential Equations, Science Press, Beijing, 2015. |
| [28] |
D. J. Higham, An algorithmic introduction to numerical simulation of stochastic differential equations, SIAM Rev., 43 (2001), 525–546. https://doi.org/10.1137/S0036144500378302 doi: 10.1137/S0036144500378302
|
| [29] |
J. B. Plotkin, J. Dushoff, S. A. Levin, Hemagglutinin sequence clusters and the antigenic evolution of influenza A virus, Proc. Natl. Acad. Sci. U.S.A., 99 (2002), 6263–6268. https://doi.org/10.1073/pnas.082110799 doi: 10.1073/pnas.082110799
|
| [30] |
X. B. Zhang, X. H. Zhang, The threshold of a deterministic and a stochastic SIQS epidemic model with varying total population size, Appl. Math. Modell., 91 (2021), 749–767. https://doi.org/10.1016/j.apm.2020.09.050 doi: 10.1016/j.apm.2020.09.050
|