Threshold behavior in a stochastic SIRS epidemic model with a logarithmic Ornstein-Uhlenbeck process

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.