This paper formulates an SEIRSHM epidemic model with general birth rate, media report and limited medical resources. Firstly, the well-posedness of the solutions and the extinction of the disease are discussed. Then, the existence of the endemic equilibrium is discussed and we find when R∗ > 1 and R 0 = 1, there exhibits a backward bifurcation, if R∗ < 1 and R 0 = 1, there exhibits a forward bifurcation. Finally, numerical simulations are carried out to illustrate the main results and show that media report and limited medical resources have a great impact on disease transmission.

Due to the popularity of the Internet and the extensive use of new media, after the occurrence of infectious diseases, the spread of social media information greatly affects the group's opinion and cognition and even the health behaviors they take, thus affecting the spread of infectious diseases. Therefore, this paper studies the event evolution from multiple dimensions.

To address this gap, we developed a three-layer model framework of major infectious disease event evolution based on subject modeling. This framework integrates three key factors-health transmission, perspective interaction, and risk perception-to analyze group perspective evolution, behavioral change, and virus transmission processes. The model's effectiveness was evaluated through simulation and sensitivity analysis. In addition, we conducted an empirical analysis by constructing a social media health transmission effect index system to identify the critical factors affecting health transmission.

Simulation results reveal that among the three factors, health transmission has the most significant impact on the evolution of group perspectives during infectious disease events. Moreover, the dynamics of public viewpoint evolution influence individual decisions regarding the adoption of non-pharmacological interventions, which are shown to effectively reduce both the transmission rate of the virus and the peak number of infections.

The findings of this study enhance our understanding of the complex mechanisms and evolutionary pathways in infectious disease events. By integrating multiple dimensions of event evolution, the proposed model offers valuable insights for the design of effective countermeasures and strategies in emergency management and response to infectious disease outbreaks.

Human behaviors play a crucial role in the intertwined dynamics of information diffusion and epidemic spreading. In turn, the spread of information and epidemics also affects individual behavioral changes. Understanding how information and epidemics propagate when individuals make coordinated decisions is essential to designing practicable control policies. To delve into this, we modify the unaware-aware-unaware susceptible-infected-susceptible (UAU-SIS) model by a two-strategy game-theory dilemma and explore how individual protective behaviors drive the interaction between information diffusion and epidemic spreading. Our theoretical analysis reveals that at the onset of an epidemic, individuals will not take any preventive measures, with the epidemic threshold being determined primarily by the topological structure of the epidemic layer. Extensive simulations help us explore the emergence of protective behaviors. A key finding is the existence of a crucial protection threshold, beyond which aware individuals begin to adopt preventive measures. Furthermore, our findings suggest that the high recovery rate or cost associated with contracting the disease, coupled with the poor failure rate of preventive measures or the low forgotten rate, leads to a significant number of aware individuals participating in self-protection, curbing the spread of epidemics. Moreover, even though individual protective decisions evolve in the information layer, the topology of the epidemic layer profoundly impacts both information diffusion and epidemic spreading. This work offers a new insight into the intertwined processes between information diffusion and epidemic spreading driven by human behaviors, which could help decision-makers gain some viable approaches to intervening in diseases.

This study aims at exploring a general and adaptive control strategy to confront the rapid evolution of an emerging infectious disease ('Disease X'), drawing lessons from the management of COVID-19 in China. We employ a dynamic model incorporating age structures and vaccination statuses, which is calibrated using epidemic data. We therefore estimate the cumulative infection rate (CIR) during the first epidemic wave of Omicron variant after China relaxed its zero-COVID policy to be 82.9% (95% CI: 82.3%, 83.5%), with a case fatality rate (CFR) of 0.25% (95% CI: 0.248%, 0.253%). We further show that if the zero-COVID policy had been eased in January 2022, the CIR and CFR would have decreased to 81.64% and 0.205%, respectively, due to a higher level of immunity from vaccination. However, if we ease the zero-COVID policy during the circulation of Delta variant from June 2021, the CIR would decrease to 74.06% while the CFR would significantly increase to 1.065%. Therefore, in the face of a 'Disease X', the adaptive strategies should be guided by multiple factors, the 'zero-COVID-like' policy could be a feasible and effective way for the control of a variant with relative low transmissibility. However, we should ease the strategy as the virus matures into a new variant with much higher transmissibility, particularly when the population is at a high level of immunity.

This study aims to examine the development trend of COVID-19 in China and propose a model to assess the impacts of various prevention and control measures in combating the COVID-19 pandemic. Using COVID-19 cases reported by the National Health Commission of China from January 2, 2020, to January 2, 2022, we established a Susceptible-Exposed-Infected-Asymptomatic-Quarantined-Vaccinated-Hospitalized-Removed (SEIAQVHR) model to calculate the COVID-19 transmission rate and Rt effective reproduction number, and assess prevention and control measures. Additionally, we built a stochastic model to explore the development of the COVID-19 epidemic. We modeled the incidence trends in five outbreaks between 2020 and 2022. Some important features of the COVID-19 epidemic are mirrored in the estimates based on our SEIAQVHR model. Our model indicates that an infected index case entering the community has a 50%-60% chance to cause a COVID-19 outbreak. Wearing masks and getting vaccinated were the most effective measures among all the prevention and control measures. Specifically targeting asymptomatic individuals had no significant impact on the spread of COVID-19. By adjusting prevention and control parameters, we suggest that increasing the rates of effective vaccination and mask-wearing can significantly reduce COVID-19 cases in China. Our stochastic model analysis provides a useful tool for understanding the COVID-19 epidemic in China.

A well-connected transportation network unites localities but also accelerates the transmission of infectious diseases. Subways-an important aspect of daily travel in big cities-are high-risk sites for the transmission of urban epidemics. Intensive research examining the transmission mechanisms of infectious diseases in subways is necessary to ascertain the risk of disease transmission encountered by commuters.

In this study, we improve the susceptible-exposed-infected-recovered (SEIR) model and propose the susceptible-exposed-infected-asymptomatic infected (SEIA) model. First, we added asymptomatic patients to the improved model as a parameter to explore the role of asymptomatic patients in the transmission of infectious diseases in a subway. The numbers of boarding and alighting passengers were added to the model as two time-varying parameters to simulate the exchange of passengers at each station.

The improved model could simulate the transmission of infectious diseases in subways and identify the key factors of transmission. We then produced an example of the transmission of coronavirus disease (COVID-19) in a subway using real subway passenger data substituted into the model for the calculations.

We ascertained that the number of exposed people continuously increased with the operation of the subway. Asymptomatic patients had a greater impact on the transmission of infectious diseases than infected people in the course of transmission. The SEIA model constructed in this study accurately determined the spread of infectious diseases in a subway and may also be applicable to studies on the transmission of infectious diseases in other urban public transport systems.

Control and prevention strategies are indispensable tools for managing the spread of infectious diseases. This paper examined biological models for the post-vaccination stage of a viral outbreak that integrate two important mitigation tools: social distancing, aimed at reducing the disease transmission rate, and vaccination, which boosts the immune system. Five different scenarios of epidemic progression were considered: (ⅰ) the "no control" scenario, reflecting the natural evolution of a disease without any safety measures in place, (ⅱ) the "reconstructed" scenario, representing real-world data and interventions, (ⅲ) the "social distancing control" scenario covering a broad set of behavioral changes, (ⅳ) the "vaccine control" scenario demonstrating the impact of vaccination on epidemic spread, and (ⅴ) the "both controls concurrently" scenario incorporating social distancing and vaccine controls simultaneously. By comparing these scenarios, we provided a comprehensive analysis of various intervention strategies, offering valuable insights into disease dynamics. Our innovative approach to modeling the cost of control gave rise to a robust computational algorithm for solving optimal control problems associated with different public health regulations. Numerical results were supported by real data for the Delta variant of the COVID-19 pandemic in the United States.

In this paper we develop a new multi-objective simulated annealing (MOSA) algorithm to generate optimal testing protocols for infectious diseases, using the COVID-19 pandemic as our context. A SEIR (susceptible-exposed-infected-recovered) epidemiological model is embedded as the computational platform for our MOSA algorithm to optimize testing protocols for screening across three joint objectives: minimum cost of test materials, minimum total infections over the testing horizon, and minimum number of false negatives over the horizon. We demonstrate the application of this optimization tool to recommend screening protocols for K-12 school districts in the U.S. State of North Carolina. Our approach is scalable by population coverage and can be employed at the level of individual school districts or regional collections of districts, individual schools or collections of schools across a district, business sites, or nursing homes, among other congregate settings where individuals may be screened prior to gaining entry to the site. The algorithm can be solved two ways, generating either independent optimal protocols across individual testing locations, or a common protocol covering all locations in the collection of testing sites. Our findings can be used to inform policy decisions to guide the development of effective testing strategies for controlling the spread of COVID-19 or other pandemic diseases in a wide range of congregate settings across various geographic regions.

We propose a versatile model with a flexible choice of control for an early-pandemic outbreak prevention when vaccine/drug is not yet available. At that stage, control is often limited to non-medical interventions like social distancing and other behavioral changes. For the SIR optimal control problem, we show that the running cost of control satisfying mild, practically justified conditions generates an optimal strategy, u(t), t ∈ [0, T], that is sustainable up until some moment τ ∈ [0, T). However, for any t ∈ [τ, T], the function u(t) will decline as t approaches T, which may cause the number of newly infected people to increase. So, the window from 0 to τ is the time for public health officials to prepare alternative mitigation measures, such as vaccines, testing, antiviral medications, and others. In addition to theoretical study, we develop a fast and stable computational method for solving the proposed optimal control problem. The efficiency of the new method is illustrated with numerical examples of optimal control trajectories for various cost functions and weights. Simulation results provide a comprehensive demonstration of the effects of control on the epidemic spread and mitigation expenses, which can serve as invaluable references for public health officials.

Compartmental models have emerged as robust computational frameworks and have yielded remarkable success in the fight against COVID-19. This study proposes a vaccination-based compartmental model for COVID-19 transmission dynamics. The model reflects the specific stages of COVID-19 infection and integrates a vaccination strategy, allowing for a comprehensive analysis of how vaccination rates influence the disease spread. We fit this model to daily confirmed COVID-19 cases in Tennessee, United States of America (USA), from June 4 to November 26, 2021, in a Bayesian inference approach using the Hamiltonian Monte Carlo (HMC) algorithm. First, excluding vaccination dynamics from the model, we estimated key epidemiological parameters like infection, recovery, and disease-induced death rates. This analysis yielded a basic reproduction number (R 0 ) of 1.5. Second, we incorporated vaccination dynamics and estimated the vaccination rate for three vaccines: 0.0051 per day for both Pfizer and Moderna and 0.0059 per day for Janssen. The fitted curves show reductions in the epidemic peak for all three vaccines. Pfizer and Moderna vaccines bring the peak down from 8,029 infected cases to 5,616 infected cases, while the Janssen vaccine reduces it, to 6,493 infected cases. Simulations of the model by varying the vaccination rate and vaccine efficacy were performed. A highly effective vaccine (95% efficacy) with a daily vaccination rate of 0.006 halved COVID-19 infections, reducing cases from 8,029 to around 4,000. The results also show that the model's prediction accuracy for new observations improves with the number of observed data used to train the model.

Accurate prediction of epidemics is pivotal for making well-informed decisions for the control of infectious diseases, but addressing heterogeneity in the system poses a challenge. In this study, we propose a novel modelling framework integrating the spatio-temporal heterogeneity of susceptible individuals into homogeneous models, by introducing a continuous recruitment process for the susceptibles. A neural network approximates the recruitment rate to develop a Universal Differential Equations (UDE) model. Simultaneously, we pre-set a specific form for the recruitment rate and develop a mechanistic model. Data from a COVID Omicron variant outbreak in Shanghai are used to train the UDE model using deep learning methods and to calibrate the mechanistic model using MCMC methods. Subsequently, we project the attack rate and peak of new infections for the first Omicron wave in China after the adjustment of the dynamic zero-COVID policy. Our projections indicate an attack rate and a peak of new infections of 80.06% and 3.17% of the population, respectively, compared with the homogeneous model's projections of 99.97% and 32.78%, thus providing an 18.6% improvement in the prediction accuracy based on the actual data. Our simulations demonstrate that heterogeneity in the susceptibles decreases herd immunity for ~37.36% of the population and prolongs the outbreak period from ~30 days to ~70 days, also aligning with the real case. We consider that this study lays the groundwork for the development of a new class of models and new insights for modelling heterogeneity.

Parameter identification involves the estimation of undisclosed parameters within a system based on observed data and mathematical models. In this investigation, we employ DAISY to meticulously examine the structural identifiability of parameters of a within-host SARS-CoV-2 epidemic model, taking into account an array of observable datasets. Furthermore, Monte Carlo simulations are performed to offer a comprehensive practical analysis of model parameters. Lastly, sensitivity analysis is employed to ascertain that decreasing the replication rate of the SARS-CoV-2 virus and curbing the infectious period are the most efficacious measures in alleviating the dissemination of COVID-19 amongst hosts.

The complex interactions were performed among non-pharmaceutical interventions, vaccinations, and hosts for all epidemics in mainland China during the spread of COVID-19. Specially, the small-scale epidemic in the city described by SVEIR model was less found in the current studies. The SVEIR model with control was established to analyze the dynamical and epidemiological features of two epidemics in Jinzhou City led by Omicron variants before and after Twenty Measures. In this study, the total population (N) of Jinzhou City was divided into five compartments: the susceptible (S), the vaccinated (V), the exposed (E), the infected (I), and the recovered (R). By surveillance data and the SVEIR model, three methods (maximum likelihood method, exponential growth rate method, next generation matrix method) were governed to estimate basic reproduction number, and the results showed that an increasing tendency of basic reproduction number from Omicron BA.5.2 to Omicron BA.2.12.1. Meanwhile, the effective reproduction number for two epidemics were investigated by surveillance data, and the results showed that Jinzhou wave 1 reached the peak on November 1 and was controlled 7 days later, and that Jinzhou wave 2 reached the peak on November 28 and was controlled 5 days later. Moreover, the impacts of non-pharmaceutical interventions (awareness delay, peak delay, control intensity) were discussed extensively, the variations of infection scales for Omicron variant and EG.5 variant were also discussed. Furthermore, the investigations on peaks and infection scales for two epidemics in dynamic zero-COVID policy were operated by the SVEIR model with control. The investigations on public medical requirements of Jinzhou City and Liaoning Province were analyzed by using SVEIR model without control, which provided a possible perspective on variant evolution in the future.

Accurate forecasting and analysis of emerging pandemics play a crucial role in effective public health management and decision-making. Traditional approaches primarily rely on epidemiological data, overlooking other valuable sources of information that could act as sensors or indicators of pandemic patterns. In this paper, we propose a novel framework, MGLEP, that integrates temporal graph neural networks and multi-modal data for learning and forecasting. We incorporate big data sources, including social media content, by utilizing specific pre-trained language models and discovering the underlying graph structure among users. This integration provides rich indicators of pandemic dynamics through learning with temporal graph neural networks. Extensive experiments demonstrate the effectiveness of our framework in pandemic forecasting and analysis, outperforming baseline methods across different areas, pandemic situations, and prediction horizons. The fusion of temporal graph learning and multi-modal data enables a comprehensive understanding of the pandemic landscape with less time lag, cheap cost, and more potential information indicators.

Differences in transmissibility of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) variants of concern (VOCs) in different districts are hard to assess. To address this, our study focused on calculating the Real-time reproduction number (R t ) for these variants in different regions.

According to the criteria defined by the World Health Organization (WHO), the global landscape was categorized into six distinct regions. In each region, the predominant SARS-CoV-2 variant was first identified based on the proportion of variant sequencing analysis results. Then, using serial interval (SI) parameters, we calculated R t for the relevant Variant of Concern (VOC) in each region. This approach enabled us to compare the R t values of the same variant across different regions and analyze the transmissibility of each region's variant in relation to the overall situation in that region.

The progression of VOC for SARS-CoV-2 shows regional variations. However, a common sequence of evolution is observed: Wild-type → Alpha → Beta → Delta → Omicron. Moreover, an increasing trend is discerned within diverse regions where the shift in R t of distinct VOC corresponds with the overarching R t route of SARS-CoV-2 in specific regions.

As the COVID-19 pandemic advances, regional epidemiological trends are aligning, likely due to similar virus mutations and shared public health strategies, suggesting opportunities for standardized global responses.

To effectively combat emerging infectious diseases like COVID-19, it is crucial to adopt strict prevention and control measures promptly to effectively contain the spread of the epidemic. In this paper, we propose a transmission model to investigate the influence of two control strategies: reducing contact numbers and improving medical resources. We examine these strategies in terms of constant control and time-varying control. Through sensitivity analysis on two reproduction numbers of the model with constant control, we demonstrate that reducing contact numbers is more effective than improving medical resources. Furthermore, these two constant controls significantly influence the peak values and timing of infections. Specifically, intensifying control measures can reduce peak values, albeit at the expense of delaying the peak time. In the model with time-varying control, we initially explore the corresponding optimal control problem and derive the characteristic expression of optimal control. Subsequently, we utilize real data from January 10th to April 12th, 2020, in Wuhan city as a case study to perform parameter estimation by using our proposed improved algorithm. Our findings illustrate that implementing optimal control measures can effectively reduce infections and deaths, and shorten the duration of the epidemic. Then, we numerically explore that implementing control measures promptly and increasing intensity to reduce contact numbers can make actual control be more closer to optimized control. Finally, we utilize the real data from October 31st to November 18th, 2021, in Hebei province as a second case study to validate the feasibility of our proposed suggestions.

In this paper, a new stochastic epidemic model is established and the dynamical behavior of its solutions is studied for this model. A deterministic epidemic model (ordinary differential equation) is first proposed by considering the isolation mechanism, and the transmission probability function is determined by a Wells-Riley model method to analyze the transmission in the quarantine. For this deterministic model, the basic reproduction number R0 is computed and it is used to determine the existence of disease-free and positive equilibria. The linearized stability of the equilibria is also discussed by analyzing the distribution of eigenvalues of the linear system. Following that, a corresponding stochastic epidemic model is further established by introducing stochastic disturbance. Then, the extinction result of the model is derived also with the help of the basic reproduction number R0s. Furthermore, by applying the theory of Markov semigroups, it is proved that the densities of the distributions of the solutions can converge to an invariant density or sweeping under certain conditions. At last, some numerical simulations are provided and discussed to illustrate the practicability of the model and the obtained theoretical results.

The design of optimized non-pharmaceutical interventions (NPIs) is critical to the effective control of emergent outbreaks of infectious diseases such as SARS, A/H1N1 and COVID-19 and to ensure that numbers of hospitalized cases do not exceed the carrying capacity of medical resources. To address this issue, we formulated a classic SIR model to include a close contact tracing strategy and structured prevention and control interruptions (SPCIs). The impact of the timing of SPCIs on the maximum number of non-isolated infected individuals and on the duration of an infectious disease outside quarantined areas (i.e. implementing a dynamic zero-case policy) were analyzed numerically and theoretically. These analyses revealed that to minimize the maximum number of non-isolated infected individuals, the optimal time to initiate SPCIs is when they can control the peak value of a second rebound of the epidemic to be equal to the first peak value. More individuals may be infected at the peak of the second wave with a stronger intervention during SPCIs. The longer the duration of the intervention and the stronger the contact tracing intensity during SPCIs, the more effective they are in shortening the duration of an infectious disease outside quarantined areas. The dynamic evolution of the number of isolated and non-isolated individuals, including two peaks and long tail patterns, have been confirmed by various real data sets of multiple-wave COVID-19 epidemics in China. Our results provide important theoretical support for the adjustment of NPI strategies in relation to a given carrying capacity of medical resources.

During the global COVID-19 pandemic, densely populated megacities engaged in active international exchanges have faced the most severe impacts from both the disease and the associated infodemic. This study examines the factors influencing public participation behavior on government microblogs in these megacities during the pandemic. It guides megacities in disseminating epidemic information, promoting knowledge on epidemic prevention, managing public opinion, and addressing related matters.

Utilizing the elaboration likelihood model's central and peripheral routes, drawing on an empirical analysis of 6,677 epidemic-related microblogs from seven Chinese megacities, this study analyses the influence mechanisms influencing public participation behavior and reveals the regulatory role of confirmed case numbers. Meanwhile,a qualitative comparative analysis examines and discusses diferent confgurations of ixn fuential factors.

The study reveals that microblog content richness demonstrates a U-shaped impact on public participation behavior. Conversely, content interaction, content length, and the number of fans positively impact participation, while update frequency has a negative impact. Additionally, the number of new confrmed cases positively regulates the impact of microblog content and publisher characteristics on public participation behavior. Public participation behavior also varies based on publishing time and content semantic features. This study further revealed the different confgurations of influential factors by QCA method.

This study reveals the impact mechanism of the microblog content and publisher characteristics on public participation behavior. It also demonstrates the regulatory role of newly confrmed cases in the way content and publishers' characteristics influence public participation behavior. This study is of great significance for the operation of government microblogs, the release of emergency information, and the promotion of public participation.

In recent years, advanced regularization techniques have emerged as a powerful tool aimed at stable estimation of infectious disease parameters that are crucial for future projections, prevention, and control. Unlike other system parameters, i.e., incubation and recovery rates, the case reporting rate, Ψ, and the time-dependent effective reproduction number, R e ( t ) , are directly influenced by a large number of factors making it impossible to pre-estimate these parameters in any meaningful way. In this study, we propose a novel iteratively-regularized trust-region optimization algorithm, combined with SuSvIuIvRD compartmental model, for stable reconstruction of Ψ and R e ( t ) from reported epidemic data on vaccination percentages, incidence cases, and daily deaths. The innovative regularization procedure exploits (and takes full advantage of) a unique structure of the Jacobian and Hessian approximation for the nonlinear observation operator. The proposed inversion method is thoroughly tested with synthetic and real SARS-CoV-2 Delta variant data for different regions in the United States of America from July 9, 2021, to November 25, 2021. Our study shows that case reporting rate during the Delta wave of COVID-19 pandemic in the US is between 12% and 37%, with most states being in the range from 15% to 25%. This confirms earlier accounts on considerable under-reporting of COVID-19 cases due to the impact of "silent spreaders" and the limitations of testing.

Non-pharmaceutical interventions (NPIs) and vaccines have been widely used to manage the COVID-19 pandemic. However, uncertainty persists regarding the effectiveness of these interventions due to data quality issues, methodological challenges, and differing contextual factors. Accurate estimation of their effects is crucial for future epidemic preparedness.

To address this, we developed a population-based mechanistic model that includes the impact of NPIs and vaccines on SARS-CoV-2 transmission and hospitalization rates. Our statistical approach estimated all parameters in one step, accurately propagating uncertainty. We fitted the model to comprehensive epidemiological data in France from March 2020 to October 2021. With the same model, we simulated scenarios of vaccine rollout.

The first lockdown was the most effective, reducing transmission by 84 % (95 % confidence interval (CI) 83-85). Subsequent lockdowns had diminished effectiveness (reduction of 74 % (69-77) and 11 % (9-18), respectively). A 6 pm curfew was more effective than one at 8 pm (68 % (66-69) vs. 48 % (45-49) reduction), while school closures reduced transmission by 15 % (12-18). In a scenario without vaccines before November 2021, we predicted 159,000 or 168 % (95 % prediction interval (PI) 70-315) more deaths and 1,488,000 or 300 % (133-492) more hospitalizations. If a vaccine had been available after 100 days, over 71,000 deaths (16,507-204,249) and 384,000 (88,579-1,020,386) hospitalizations could have been averted.

Our results highlight the substantial impact of NPIs, including lockdowns and curfews, in controlling the COVID-19 pandemic. We also demonstrate the value of the 100 days objective of the Coalition for Epidemic Preparedness Innovations (CEPI) initiative for vaccine availability.

Since the global outbreak of COVID-19, the virus has continuously mutated and can survive in the air for long periods of time. This paper establishes and analyzes a model of COVID-19 with self-protection and quarantine measures affected by viruses in the environment to investigate the influence of viruses in the environment on the spread of the outbreak, as well as to develop a rational prevention and control measure to control the spread of the outbreak. The basic reproduction number was calculated and Lyapunov functions were constructed to discuss the stability of the model equilibrium points. The disease-free equilibrium point was proven to be globally asymptotically stable when $ R_0 < 1 $, and the endemic equilibrium point was globally asymptotically stable when $ R_0 > 1 $. The model was fitted using data from COVID-19 cases in Chongqing between November 1 to November 25, 2022. Based on the numerical analysis, the following conclusion was obtained: clearing the virus in the environment and strengthening the isolation measures for infected people can control the epidemic to a certain extent, but enhancing the self-protection of individuals can be more effective in reducing the risk of being infected and controlling the transmission of the epidemic, which is more conducive to the practical application.

This study aimed to develop a modified susceptible-exposed-infected-recovered (SEIR) model to evaluate monkeypox epidemics in the United States and explore more optimized prevention and control measures. To further assess the impact of public health measures on the transmission of monkeypox, different intervention scenarios were developed based on the classic SEIR model, considering reducing contact, enhancing vaccination, diagnosis delay, and environmental transmission risk, respectively. We evaluated the impact of different measures by simulating their spread in different scenarios. During the simulation period, 8709 people were infected with monkeypox. The simulation analysis showed that: (1) the most effective measures to control monkeypox transmission during the early stage of the epidemic were reducing contact and enhancing vaccination, with cumulative infections at 51.20% and 41.90% of baseline levels, respectively; (2) shortening diagnosis time would delay the peak time of the epidemic by 96 days; and (3) the risk of environmental transmission of monkeypox virus was relatively low. This study indirectly proved the effectiveness of the prevention and control measures, such as reducing contact, enhancing vaccination, shortening diagnosis time, and low risk of environmental transmission, which also provided an important reference and containment experience for nonepidemic countries.

We develop a multi-group and multi-patch model to study the effects of population dispersal on the spatial spread of vector-borne diseases across a heterogeneous environment. The movement of host and/or vector is described by Lagrangian approach in which the origin or identity of each individual stays unchanged regardless of movement. The basic reproduction number [Formula: see text] of the model is defined and the strong connectivity of the host-vector network is succinctly characterized by the residence times matrices of hosts and vectors. Furthermore, the definition and criterion of the strong connectivity of general infectious disease networks are given and applied to establish the global stability of the disease-free equilibrium. The global dynamics of the model system are shown to be entirely determined by its basic reproduction number. We then obtain several biologically meaningful upper and lower bounds on the basic reproduction number which are independent or dependent of the residence times matrices. In particular, the heterogeneous mixing of hosts and vectors in a homogeneous environment always increases the basic reproduction number. There is a substantial difference on the upper bound of [Formula: see text] between Lagrangian and Eulerian modeling approaches. When only host movement between two patches is concerned, the subdivision of hosts (more host groups) can lead to a larger basic reproduction number. In addition, we numerically investigate the dependence of the basic reproduction number and the total number of infected hosts on the residence times matrix of hosts, and compare the impact of different vector control strategies on disease transmission.

Several studies on COVID-19 pandemic have shown that the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) originating from human stool are detected in raw sewage for several days, leading to potential health risks for workers due to the production of bioaerosols and droplets during wastewater treatment process. In this study, data of SARS-CoV-2 concentrations in wastewater were gathered from literatures, and a quantitative microbial risk assessment with Monte Carlo simulation was used to estimate the daily probability of infection risk through exposure to viable infectious viral airborne particles of the workers during four seasons and under six environmental conditions. Inhalation of bioaerosols and direct ingestion of wastewater droplets were selected as exposure pathways. Spearman rank correlation coefficients were used for sensitivity analysis to identify the variables with the greatest influence on the infection risk probability. It was found that the daily probability of infection risk decreased with temperature (T) and relative humidity (RH) increase. The probability of direct droplet ingestion exposure pathway was higher than that of the bioaerosol inhalation pathway. The sensitivity analysis indicated that the most sensitive variable for both exposure pathways was the concentration of SARS-CoV-2 in stool. So, appropriate aeration systems, covering facilities, and effective ventilation are suggested to implement in wastewater treatment plants (WWTPs) to reduce emission concentration. Further to this, the exposure time (t) had a larger variance contribution than T and RH for the bioaerosol inhalation pathway. Implementing measures such as adding more work shifts, mandating personal protective equipment for all workers, and implementing coverage for treatment processes can significantly reduce the risk of infection among workers at WWTPs. These measures are particularly effective during environmental conditions with low temperatures and humidity levels.

China managed coronavirus disease 2019 (COVID-19) with measures against Class B infectious diseases, instead of Class A infectious diseases, in a major shift of its epidemic response policies. We aimed to generate robust information on the transmission dynamics of novel coronavirus infection in Shanxi, a province located in northern China, after the implementation of the "Class B infectious disease Class B management" policy.

We consolidated infection data in Shanxi province from December 6, 2022 to January 14, 2023 through a network questionnaire survey and sentinel surveillance. A dynamics model of the SEIQHCVR was developed to track the infection curves and effective reproduction number (R t ).

Our model was effective in estimating the trends of novel coronavirus infection, with the coefficient of determination (R 2 ) above 90% in infections, inpatients, and critically ill patients. The number of infections in Shanxi province as well as in urban and rural areas peaked on December 20, 2022, with the peak of inpatients and critically ill patients occurring 2 to 3 weeks after the peak of infections. By the end of January 2023, 87.72% of the Shanxi residents were predicted to be infected, and the outbreak subsequently subsided. A small wave of COVID-19 infections may re-emerge at the end of April. In less than a month, the R t values of positive infections, inpatients and critically ill patients were all below 1.0.

The outbreak in Shanxi province is currently at a low prevalence level. In the face of possible future waves of infection, there is a strong need to strengthen surveillance and early warning.

Rational allocation of limited vaccine resources is one of the key issues in the prevention and control of emerging infectious diseases. An age-structured infectious disease model with limited vaccine resources is proposed to explore the optimal vaccination ages. The effective reproduction number [Formula: see text] of the epidemic disease is computed. It is shown that the reproduction number is the threshold value for eradicating disease in the sense that the disease-free steady state is globally stable if [Formula: see text] and there exists a unique endemic equilibrium if [Formula: see text]. The effective reproduction number is used as an objective to minimize the disease spread risk. Using the epidemic data from the early spread of Wuhan, China and demographic data of Wuhan, we figure out the strategies to distribute the vaccine to the age groups to achieve the optimal vaccination effects. These analyses are helpful to the design of vaccination schedules for emerging infectious diseases.

The outbreak of the COVID-19 pandemic sparked numerous studies on policy options for managing public health emergencies, especially regarding how to choose the intensity of prevention and control to maintain a balance between economic development and disease prevention.

We constructed a cost-benefit model of COVID-19 pandemic prevention and control policies based on an epidemic transmission model. On this basis, numerical simulations were performed for different economies to analyse the dynamic evolution of prevention and control policies. These economies include areas with high control costs, as seen in high-income economies, and areas with relatively low control costs, exhibited in upper-middle-income economies.

The simulation results indicate that, at the outset of the COVID-19 pandemic, both high-and low-cost economies tended to enforce intensive interventions. However, as the virus evolved, particularly in circumstances with relatively rates of reproduction, short incubation periods, short spans of infection and low mortality rates, high-cost economies became inclined to ease restrictions, while low-cost economies took the opposite approach. However, the consideration of additional costs incurred by the non-infected population means that a low-cost economy is likely to lift restrictions as well.

This study concludes that variations in prevention and control policies among nations with varying income levels stem from variances in virus transmission characteristics, economic development, and control costs. This study can help researchers and policymakers better understand the differences in policy choice among various economies as well as the changing trends of dynamic policy choices, thus providing a certain reference value for the policy direction of global public health emergencies.

In 2022, Omicron outbreaks occurred at multiple sites in China. It is of great importance to track the incidence trends and transmission dynamics of coronavirus disease 2019 (COVID-19) to guide further interventions.

Given the population size, economic level and transport level similarities, two groups of outbreaks (Shanghai vs. Chengdu and Sanya vs. Beihai) were selected for analysis. We developed the SEAIQRD, ARIMA, and LSTM models to seek optimal modeling techniques for waves associated with the Omicron variant regarding data predictive performance and mechanism transmission dynamics, respectively. In addition, we quantitatively modeled the impacts of different combinations of more stringent interventions on the course of the epidemic through scenario analyses.

The best-performing LSTM model showed better prediction accuracy than the best-performing SEAIQRD and ARIMA models in most cases studied. The SEAIQRD model had an absolute advantage in exploring the transmission dynamics of the outbreaks. Regardless of the time to inflection point or the time to Rt curve below 1.0, Shanghai was later than Chengdu (day 46 vs. day 12/day 54 vs. day 14), and Sanya was later than Beihai (day 16 vs. day 12/day 20 vs. day 16). Regardless of the number of peak cases or the cumulative number of infections, Shanghai was higher than Chengdu (34,350 vs. 188/623,870 vs. 2,181), and Sanya was higher than Beihai (1,105 vs. 203/16,289 vs. 3,184). Scenario analyses suggested that upgrading control level in advance, while increasing the index decline rate and quarantine rate, were of great significance for shortening the time to peak and Rt below 1.0, as well as reducing the number of peak cases and final affected population.

The LSTM model has great potential for predicting the prevalence of Omicron outbreaks, whereas the SEAIQRD model is highly effective in revealing their internal transmission mechanisms. We recommended the use of joint interventions to contain the spread of the virus.

Contact tracing is commonly recommended to control outbreaks of COVID-19, but its effectiveness is unclear. Following PRISMA guidelines, we searched four databases using a range of terms related to contact tracing effectiveness for COVID-19. We found 343 papers; 32 were included. All were observational or modelling studies. Observational studies (n = 14) provided consistent, very-low certainty evidence that contact tracing (alone or in combination with other interventions) was associated with better control of COVID-19 (e.g. in Hong Kong, only 1084 cases and four deaths were recorded in the first 4.5 months of the pandemic). Modelling studies (n = 18) provided consistent, high-certainty evidence that under assumptions of prompt and thorough tracing with effective quarantines, contact tracing could stop the spread of COVID-19 (e.g. by reducing the reproduction number from 2.2 to 0.57). A cautious interpretation indicates that to stop the spread of COVID-19, public health practitioners have 2-3 days from the time a new case develops symptoms to isolate the case and quarantine at least 80% of its contacts.

A local COVID-19 outbreak with two community clusters occurred in a large industrial city, Shaoxing, China, in December 2021 after serial interventions were imposed. We aimed to understand the reason by analysing the characteristics of the outbreak and evaluating the effects of phase-adjusted interventions. Publicly available data from 7 December 2021 to 25 January 2022 were collected to analyse the epidemiological characteristics of this outbreak. The incubation period was estimated using Hamiltonian Monte Carlo method. A well-fitted extended susceptible-exposed-infectious-recovered model was used to simulate the impact of different interventions under various combination of scenarios. There were 387 SARS-CoV-2-infected cases identified, and 8.3% of them were initially diagnosed as asymptomatic cases. The estimated incubation period was 5.4 (95% CI 5.2-5.7) days for all patients. Strengthened measures of comprehensive quarantine based on tracing led to less infections and a shorter duration of epidemic. With a same period of incubation, comprehensive quarantine was more effective in containing the transmission than other interventions. Our findings reveal an important role of tracing and comprehensive quarantine in blocking community spread when a cluster occurred. Regions with tense resources can adopt home quarantine as a relatively affordable and low-impact intervention measure compared with centralized quarantine.

SARS COV-2 (Covid-19) has imposed a monumental socio-economic burden worldwide, and its impact still lingers. We propose a deterministic model to describe the transmission dynamics of Covid-19, emphasizing the effects of vaccination on the prevailing epidemic. The proposed model incorporates current information on Covid-19, such as reinfection, waning of immunity derived from the vaccine, and infectiousness of the pre-symptomatic individuals into the disease dynamics. Moreover, the model analysis reveals that it exhibits the phenomenon of backward bifurcation, thus suggesting that driving the model reproduction number below unity may not suffice to drive the epidemic toward extinction. The model is fitted to real-life data to estimate values for some of the unknown parameters. In addition, the model epidemic threshold and equilibria are determined while the criteria for the stability of each equilibrium solution are established using the Metzler approach. A sensitivity analysis of the model is performed based on the Latin Hypercube Sampling (LHS) and Partial Rank Correlation Coefficients (PRCCs) approaches to illustrate the impact of the various model parameters and explore the dependency of control reproduction number on its constituents parameters, which invariably gives insight on what needs to be done to contain the pandemic effectively. The foregoing notwithstanding, the contour plots of the control reproduction number concerning some of the salient parameters indicate that increasing vaccination coverage and decreasing vaccine waning rate would remarkably reduce the value of the reproduction number below unity, thus facilitating the possible elimination of the disease from the population. Finally, the model is solved numerically and simulated for different scenarios of disease outbreaks with the findings discussed.

There are many plausible reasons for recurrent outbreaks of emerging infectious diseases. In this paper, we develop a mathematical model to illustrate how population behavioural adaption and adaptation implementation delay, in response to the perceived infection risk, can lead to recurrent outbreak patterns. We consider the early phase of an infection outbreak when herd immunity is not reached, pathogen mutation is not considered, and seasonality is ruled out as a major contributor. We derive a transmission dynamics model coupled with the renewal equation for the disease transmission effective contacts (contact rate per unit time multiplied by the transmission probability per contact). The model incorporates two critical parameters: the population behavioural adaptation flexibility index and the behavioural change implementation delay. We show that when the behavioural change implementation delay reaches a critical value, the number of infections starts to oscillate in an equilibrium that is determined by the population behavioural adaptation flexibility. We also show that the numbers of infections at the subsequent peaks can exceed that of the first peak. This was an oblique observation globally during the early phase of the COVID-19 pandemic before variants of concern emerged, and it was an observed phenomena with the Omicron variant induced wave in areas where early interventions were successful in preventing the large outbreaks. Our model and analyses can provide partially explanation for these observations.

This study proposes a non-linear mathematical model for analysing the effect of COVID-19 dynamics on the student population in higher education institutions. The theory of positivity and boundedness of solution is used to investigate the well-posedness of the model. The disease-free equilibrium solution is examined analytically. The next-generation operator method calculates the basic reproduction number ( R 0 ) . Sensitivity analyses are carried out to determine the relative importance of the model parameters in spreading COVID-19. In light of the sensitivity analysis results, the model is further extended to an optimal control problem by introducing four time-dependent control variables: personal protective measures, quarantine (or self-isolation), treatment, and management measures to mitigate the community spread of COVID-19 in the population. Simulations evaluate the effects of different combinations of the control variables in minimizing COVID-19 infection. Moreover, a cost-effectiveness analysis is conducted to ascertain the most effective and least expensive strategy for preventing and controlling the spread of COVID-19 with limited resources in the student population.

Time series analysis is a valuable tool in epidemiology that complements the classical epidemiological models in two different ways: Prediction and forecast. Prediction is related to explaining past and current data based on various internal and external influences that may or may not have a causative role. Forecasting is an exploration of the possible future values based on the predictive ability of the model and hypothesized future values of the external and/or internal influences. The time series analysis approach has the advantage of being easier to use (in the cases of more straightforward and linear models such as Auto-Regressive Integrated Moving Average). Still, it is limited in forecasting time, unlike the classical models such as Susceptible-Exposed-Infectious-Removed. Its applicability in forecasting comes from its better accuracy for short-term prediction. In its basic form, it does not assume much theoretical knowledge of the mechanisms of spreading and mutating pathogens or the reaction of people and regulatory structures (governments, companies, etc.). Instead, it estimates from the data directly. Its predictive ability allows testing hypotheses for different factors that positively or negatively contribute to the pandemic spread; be it school closures, emerging variants, etc. It can be used in mortality or hospital risk estimation from new cases, seroprevalence studies, assessing properties of emerging variants, and estimating excess mortality and its relationship with a pandemic.

Objective The aim of this study was to investigate change in the stigma that emerged during the COVID-19 pandemic over time and the factors responsible for the change. Methods Individuals with COVID-19 who presented to Ankara Medicalpark and VM Medicalpark Hospitals' Internal Diseases and Chest Diseases polyclinic between May 2021 and April 2022 were examined. The volunteers were divided into two groups: those who had COVID-19 within the first six months of the pandemic (group 1) and those who had it in the second six months (group 2). The questionnaire assessing stigma consisted of 29 propositions that participants could mark whether they agreed with them or not. Results The median age of the volunteers was 38 years. Eighty-eight (69.3%) had the disease in the first six months of the pandemic and 39 (30.7%) in the second six months. Moreover, 76.1% of the participants in the first group and 94.9% of those in the second group did not agree with the statement "I thought COVID-19 was a punishment for me" (p=0.011). Further, 56.8% of the participants in the first group and 97.4% of those in the second group stated that they did not agree with the statement "Employers may terminate the employment of employees who they find out have contracted COVID-19" (p<0.001). 80.7% of the participants in the first group and 38.5% of those in the second group agreed with the statement "There was social discrimination against people who caught COVID-19" (p<0.001).  Conclusions At the beginning of the pandemic, the participants had concerns about losing their status and jobs, but this anxiety decreased over time. Stigma in the first six months of the pandemic was greater than that in the second six months, and discrimination related to stigma decreased with recognition of the disease and the increase in experience.

During the COVID-19 pandemic, control measures, especially massive contact tracing following prompt quarantine and isolation, play an important role in mitigating the disease spread, and quantifying the dynamic contact rate and quarantine rate and estimate their impacts remain challenging. To precisely quantify the intensity of interventions, we develop the mechanism of physics-informed neural network (PINN) to propose the extended transmission-dynamics-informed neural network (TDINN) algorithm by combining scattered observational data with deep learning and epidemic models. The TDINN algorithm can not only avoid assuming the specific rate functions in advance but also make neural networks follow the rules of epidemic systems in the process of learning. We show that the proposed algorithm can fit the multi-source epidemic data in Xi'an, Guangzhou and Yangzhou cities well, and moreover reconstruct the epidemic development trend in Hainan and Xinjiang with incomplete reported data. We inferred the temporal evolution patterns of contact/quarantine rates, selected the best combination from the family of functions to accurately simulate the contact/quarantine time series learned by TDINN algorithm, and consequently reconstructed the epidemic process. The selected rate functions based on the time series inferred by deep learning have epidemiologically reasonable meanings. In addition, the proposed TDINN algorithm has also been verified by COVID-19 epidemic data with multiple waves in Liaoning province and shows good performance. We find the significant fluctuations in estimated contact/quarantine rates, and a feedback loop between the strengthening/relaxation of intervention strategies and the recurrence of the outbreaks. Moreover, the findings show that there is diversity in the shape of the temporal evolution curves of the inferred contact/quarantine rates in the considered regions, which indicates variation in the intensity of control strategies adopted in various regions.

Since the end of 2019 the COVID-19 repeatedly surges with most countries/territories experiencing multiple waves, and mechanism-based epidemic models played important roles in understanding the transmission mechanism of multiple epidemic waves. However, capturing temporal changes of the transmissibility of COVID-19 during the multiple waves keeps ill-posed problem for traditional mechanism-based epidemic compartment models, because that the transmission rate is usually assumed to be specific piecewise functions and more parameters are added to the model once multiple epidemic waves involved, which poses a huge challenge to parameter estimation. Meanwhile, data-driven deep neural networks fail to discover the driving factors of repeated outbreaks and lack interpretability. In this study, aiming at developing a data-driven method to project time-dependent parameters but also merging the advantage of mechanism-based models, we propose a transmission dynamics informed neural network (TDINN) by encoding the SEIRD compartment model into deep neural networks. We show that the proposed TDINN algorithm performs very well when fitting the COVID-19 epidemic data with multiple waves, where the epidemics in the United States, Italy, South Africa, and Kenya, and several outbreaks the Omicron variant in China are taken as examples. In addition, the numerical simulation shows that the trained TDINN can also perform as a predictive model to capture the future development of COVID-19 epidemic. We find that the transmission rate inferred by the TDINN frequently fluctuates, and a feedback loop between the epidemic shifting and the changes of transmissibility drives the occurrence of multiple waves. We observe a long response delay to the implementation of control interventions in the four countries, while the decline of the transmission rate in the outbreaks in China usually happens once the implementation of control interventions. The further simulation show that 17 days' delay of the response to the implementation of control interventions lead to a roughly four-fold increase in daily reported cases in one epidemic wave in Italy, which suggest that a rapid response to policies that strengthen control interventions can be effective in flattening the epidemic curve or avoiding subsequent epidemic waves. We observe that the transmission rate in the outbreaks in China is already decreasing before enhancing control interventions, providing the evidence that the increasing of the epidemics can drive self-conscious behavioural changes to protect against infections.

This study aimed to predict the transmission trajectory of the 2019 Corona Virus Disease (COVID-19) and analyze the impact of preventive measures on the spread of the epidemic. Considering that tracking a long-term epidemic trajectory requires explanatory modeling with more complexities than short-term predictions, an improved Susceptible-Exposed-Infected-Removed (SEIR) transmission dynamic model is established. The model depends on defining various parameters that describe both the virus and the population under study. However, it is likely that several of these parameters will exhibit significant variations among different states. Therefore, regression algorithms and heuristic algorithms were developed to effectively adapt the population-dependent parameters and ensure accurate fitting of the SEIR model to data for any specific state. In this study, we consider the second outbreak of COVID-19 in Italy as a case study, which occurred in August 2020. We divide the epidemic data from February to September of the same year into two distinct stages for analysis. The numerical results demonstrate that the improved SEIR model effectively simulates and predicts the transmission trajectories of the Italian epidemic during both periods before and after the second outbreak. By analyzing the impact of anti-epidemic measures on the spread of the disease, our findings emphasize the significance of implementing anti-epidemic preventive measures in COVID-19 modeling.

The COVID-19 infection cases have surged globally, causing devastations to both the society and economy. A key factor contributing to the sustained spreading is the presence of a large number of asymptomatic or hidden spreaders, who mix among the susceptible population without being detected or quarantined. Due to the continuous emergence of new virus variants, even if vaccines have been widely used, the detection of asymptomatic infected persons is still important in the epidemic control. Based on the unique characteristics of COVID-19 spreading dynamics, here we propose a theoretical framework capturing the transition probabilities among different infectious states in a network, and extend it to an efficient algorithm to identify asymptotic individuals. We find that using pure physical spreading equations, the hidden spreaders of COVID-19 can be identified with remarkable accuracy, even with incomplete information of the contract-tracing networks. Furthermore, our framework can be useful for other epidemic diseases that also feature asymptomatic spreading.

Non-pharmaceutical interventions (NPIs) implemented on China-bound travel have successfully mitigated cross-regional transmission of COVID-19 but made the country face ripple effects. Thus, adjusting these interventions to reduce interruptions to individuals' daily life while minimizing transmission risk was urgent.

An improved Susceptible-Infected-Recovered (SIR) model was built to evaluate the Delta variant's epidemiological characteristics and the impact of NPIs. To explore the risk associated with inbound travelers and the occurrence of domestic traceable outbreaks, we developed an association parameter that combined inbound traveler counts with a time-varying initial value. In addition, multiple time-varying functions were used to model changes in the implementation of NPIs. Related parameters of functions were run by the MCSS method with 1,000 iterations to derive the probability distribution. Initial values, estimated parameters, and corresponding 95% CI were obtained. Reported existing symptomatic, suspected, and asymptomatic case counts were used as the training datasets. Reported cumulative recovered individual data were used to verify the reliability of relevant parameters. Lastly, we used the value of the ratio (Bias2/Variance) to verify the stability of the mathematical model, and the effects of the NPIs on the infected cases to analyze the sensitivity of input parameters.

The quantitative findings indicated that this improved model was highly compatible with publicly reported data collected from July 21 to August 30, 2021. The number of inbound travelers was associated with the occurrence of domestic outbreaks. A proportional relationship between the Delta variant incubation period and PCR test validity period was found. The model also predicted that restoration of pre-pandemic travel schedules while adhering to NPIs requirements would cause shortages in health resources. The maximum demand for hospital beds would reach 25,000/day, the volume of PCR tests would be 8,000/day, and the number of isolation rooms would reach 800,000/day within 30 days.

With the pandemic approaching the end, reexamining it carefully helps better address future outbreaks. This predictive model has provided scientific evidence for NPIs' effectiveness and quantifiable evidence of health resource allocation. It could guide the design of future epidemic prevention and control policies, and provide strategic recommendations on scarce health resource allocation.

Media coverage can greatly impact the spread of infectious diseases. Taking into consideration the impacts of media coverage, we propose an SEIR model with a media coverage mediated nonlinear infection force. For this novel disease model, we identify the basic reproduction number using the next generation matrix method and establish the global threshold results: If the basic reproduction number $ \mathcal{R}_{0} < 1 $, then the disease-free equilibrium $ P_{0} $ is stable, and the disease dies out. If $ \mathcal{R}_{0} > 1 $, then the endemic equilibrium $ P^{*} $ is stable, and the disease persists. Sensitivity analysis indicates that the basic reproduction number $ \mathcal{R}_{0} $ is most sensitive to the population recruitment rate $ \Lambda $ and the disease transmission rate $ \beta _{1} $.

On September 28, 2022, the first case of Omicron subvariant BF.7 was discovered among coronavirus disease 2019 (COVID-19) infections in Hohhot, China, and then the epidemic broke out on a large scale during the National Day holiday. It is imminently necessary to construct a mathematical model to investigate the transmission dynamics of COVID-19 in Hohhot.

In this study, we first investigated the epidemiological characteristics of COVID-19 cases in Hohhot, including the spatiotemporal distribution and sociodemographic distribution. Then, we proposed a time-varying Susceptible-Quarantined Susceptible-Exposed-Quarantined Exposed-Infected-Asymptomatic-Hospitalized-Removed (SQEIAHR) model to derive the epidemic curves. The next-generation matrix method was used to calculate the effective reproduction number (Re). Finally, we explored the effects of higher stringency measures on the development of the epidemic through scenario analysis.

Of the 4,889 positive infected cases, the vast majority were asymptomatic and mild, mainly concentrated in central areas such as Xincheng District. People in the 30-59 age group primarily were affected by the current outbreak, accounting for 53.74%, but females and males were almost equally affected (1.03:1). Community screening (35.70%) and centralized isolation screening (26.28%) were the main ways to identify positive infected cases. Our model predicted the peak of the epidemic on October 6, 2022, the dynamic zero-COVID date on October 15, 2022, a number of peak cases of 629, and a cumulative number of infections of 4,963 (95% confidential interval (95%CI): 4,692 ~ 5,267), all four of which were highly consistent with the actual situation in Hohhot. Early in the outbreak, the basic reproduction number (R0) was approximately 7.01 (95%CI: 6.93 ~ 7.09), and then Re declined sharply to below 1.0 on October 6, 2022. Scenario analysis of higher stringency measures showed the importance of decreasing the transmission rate and increasing the quarantine rate to shorten the time to peak, dynamic zero-COVID and an Re below 1.0, as well as to reduce the number of peak cases and final affected population.

Our model was effective in predicting the epidemic trends of COVID-19, and the implementation of a more stringent combination of measures was indispensable in containing the spread of the virus.

By 31 May 2022, original/Alpha, Delta and Omicron strains induced 101 outbreaks of COVID-19 in mainland China. Most outbreaks were cleared by combining non-pharmaceutical interventions (NPIs) with vaccines, but continuous virus variations challenged the dynamic zero-case policy (DZCP), posing questions of what are the prerequisites and threshold levels for success? And what are the independent effects of vaccination in each outbreak? Using a modified classic infectious disease dynamic model and an iterative relationship for new infections per day, the effectiveness of vaccines and NPIs was deduced, from which the independent effectiveness of vaccines was derived. There was a negative correlation between vaccination coverage rates and virus transmission. For the Delta strain, a 61.8% increase in the vaccination rate (VR) reduced the control reproduction number (CRN) by about 27%. For the Omicron strain, a 20.43% increase in VR, including booster shots, reduced the CRN by 42.16%. The implementation speed of NPIs against the original/Alpha strain was faster than the virus's transmission speed, and vaccines significantly accelerated the DZCP against the Delta strain. The CRN ([Formula: see text]) during the exponential growth phase and the peak time and intensity of NPIs were key factors affecting a comprehensive theoretical threshold condition for DZCP success, illustrated by contour diagrams for the CRN under different conditions. The DZCP maintained the [Formula: see text] of 101 outbreaks below the safe threshold level, but the strength of NPIs was close to saturation especially for Omicron, and there was little room for improvement. Only by curbing the rise in the early stage and shortening the exponential growth period could clearing be achieved quickly. Strengthening China's vaccine immune barrier can improve China's ability to prevent and control epidemics and provide greater scope for the selection and adjustment of NPIs. Otherwise, there will be rapid rises in infection rates and an extremely high peak and huge pressure on the healthcare system, and a potential increase in excess mortality.