 Research article
 Open Access
 Published:
Modelling the pulse populationwide nucleic acid screening in mitigating and stopping COVID19 outbreaks in China
BMC Infectious Diseases volume 23, Article number: 280 (2023)
Abstract
Background
During 20212022, mainland China experienced multiple times of local COVID19 outbreaks in several cities, including Yangzhou, Xi’an etc., and the Chinese government persistently adopted the zeroCOVID policy in combating with the local outbreaks.
Methods
We develop a mathematical model with pulse populationwide nucleic acid screening, part of the zeroCOVID policy, to reveal its role in controlling the spread of COVID19. We calibrate the model by fitting the COVID19 epidemic data of the local outbreaks in Yangzhou and Xi’an, China. Sensitivity analysis is conducted to investigate the impact of populationwide nucleic acid screening on controlling the outbreak of COVID19.
Results
Without the screening, the cumulative number of confirmed cases increases by \(77.7\%\) and \(62.2\%\) in Yangzhou and Xi’an, respectively. Meanwhile, the screening program helps to shorten the lockdown period for more than one month when we aim at controlling the cases into zero. Considering its role in mitigating the epidemics, we observe a paradox phenomenon of the screening rate in avoiding the runs on medical resource. That is, the screening will aggravate the runs on medical resource when the screening rate is small, while it helps to relieve the runs on medical resource if the screening rate is high enough. We also conclude that the screening has limited effects on mitigating the epidemics if the outbreak is in a high epidemic level or there has already been runs on medical resources. Alternatively, a smaller screening population per time with a higher screening frequency may be a better program to avoid the runs on medical resources.
Conclusions
The populationwide nucleic acid screening strategy plays an important role in quickly controlling and stopping the local outbreaks under the zeroCOVID policy. However, it has limited impacts and even increase the potential risk of the runs on medical resource for containing the large scale outbreaks.
Background
An emerging coronavirus, known as “Severe Acute Respiratory Syndrome coronavirus type 2” (SARSCoV2) caused a worldwide coronavirus disease 2019 (COVID19) pandemic [1, 2]. Up until October 3, 2022, the global number of confirmed cases exceeds 620 million [3]. In China, the nationwide epidemic wave centered on Wuhan, Hubei province was contained by April 2020 with strict and effective nonpharmaceutical interventions (NPIs). In the nearly two years, the Chinese government persistently adopted the zeroCOVID policy in combating with the COVID19 epidemics. Here, the center purpose of the zeroCOVID policy is to control the infections into zero in a short time period with a small epidemic size [4,5,6,7]. To this end, a package of strict NPIs, including the contact distancing and quarantine strategies, also the largescale populationwide nucleic acid screening program, are carried out. Therefore, although China experienced more than 100 independent outbreaks (with a starting time and end time of reported cases) in different cities, the implementation of zeroCOVID policy has helped to control the infections into zero in almost one month, consequently there is long epidemicfree period in almost cities except the period of the local outbreaks. That is, the real epidemic data from the National Health Commission of the People’s Republic of China strongly supports the fact that the zeroCOVID policy indeed achieved its goal [8]. Therefore, quantifying the role of the zeroCOVID policy would be essential not only for controlling COVID19, but also for other emerging infectious diseases.
As we mentioned above, in the zeroCOVID policy in China, the populationwide nucleic acid testing is one of the key control interventions. Here, the population wide testing aims at testing the whole population in the city with epidemic every few days (like two or three days), where “teninone” or “twentyinone” pooled nucleic acid screening test method is adapted [9, 10]. The mentioned pooled testing method refers to mix multiple (ten or twenty) samples before testing. There is no infected individuals in the sample if the testing result is negative, while the further individual sample method is needed to identify the infected individual if the result is positive [11], and it was found that up to thirty samples pooled in a pool can increase test capacity with existing test resources and detect positive samples with sufficient diagnostic accuracy [12]. Obviously, the populationwide testing can help to quickly detect the infected individuals, including the asymptomatic infections, hence, it helps to quickly isolate the infections. Therefore, the largescale populationwide screening program gradually becomes an effective intervention to rapidly contain the spread of epidemic during an outbreak response [13,14,15]. Actually, the implementation scheme of largescale populationwide testing depends on the existing equipment and test kits. Considering the limitation of test resources, Hernandez and Valentinotti studied a optimal sample selection strategy which selected 5% or 10% of the whole population with excluding the sample of individuals who have a low probability of being infected (recently tested ones). Then under fixed numbers of tests, this strategy can effectively reduce the overall infectedpersondays [16]. However, it should be mentioned that on the background of the COVID19 local outbreaks in China, the nucleic acid test is implemented every one to three days of the whole population [17, 18]. For example, during the local outbreaks happening at Yangzhou (2021/7/28 to 2021/8/26) and Xi’an (2021/12/9 to 2022/1/9), the comprehensive and routine population wide screening strategy ensures the timeliness of early case detection and interruption of the spread of COVID19, which plays critically important roles in the prevention and control strategy. However, the quantitative role of largescale populationwide nucleic acid screening in mitigating the COVID19 epidemics in China remains unclear.
Mathematical models are frequently used as the powerful tools to study the transmission dynamics of COVID19 [1, 2], evaluate the effectiveness of control interventions [19, 20] or the mass vaccination program [21, 22]. Given the populationwide nucleic acid testing, the authors used a statistical model to study the optimal number of samples combined into a pool under different infection rates [11]. It should be mentioned that the populationwide nucleic acid testings are carried out every few days, which are actually the combination of the dynamic system and discrete events. In the existing modelling frameworks, fixedmoments impulsive differential equations were widely used in modelling the intermittent prevention and control measures for containing the spread of infectious diseases, and these models assumed that measures are carried out at fixed discrete times [23,24,25,26]. Therefore, the fixedmoments impulsive differential equations can be the potential modelling frameworks to describe the pulsed largescale testing. To our best knowledge, few modelling studies quantitatively investigated the impact of pulsed population wide nucleic acid screening strategy on controlling and stopping the spread of COVID19. Qualifying the issues through mathematical models follows the scope of this study.
The main purpose of this study is to propose the epidemic model considering nucleic acid screening, focus on exploring its role to contain the COVID19 outbreaks, and providing the important guidance for the implementation of the largescale population wide screening in controlling the future waves. The rest of the paper is organized as follows. In Methods section, we firstly extend the classic SEIR model characterizing the disease progress to a deterministic impulsive dynamical model to describe the transmission of COVID19 with the pulse nucleic acid screening. The data and the fitting methods are also described in details, and we calibrate the model by fitting it to the epidemic data of the local outbreaks in Yangzhou city and Xi’an city, China. In Results section, we reveal the role of nucleic acid screening through the scenario analysis in reducing the cumulative cases and shortening the time period of lockdown. We further extend the original model to a pulse model with periodic screening, and discuss the impact of the largescale screening on mitigating the epidemic and the runs on medical resources. In Discussion and conclusion section, we give some discussions and conclusions.
Methods
Model
According to the epidemic status of COVID19 infections and the implementation of control interventions, including the close contact tracing and quarantine and the nucleic acid screening adopted in the local outbreaks in Yangzhou and Xi’an, the total population N is divided into seven compartments: susceptible (S), exposed (E), infected (I), quarantined susceptible (\(S_q\)), quarantined exposed (\(E_q\)), confirmed and isolated (H) and recovered (R). It should be noticed that in our model framework, the infected population (I) will move into the confirmed and isolated class (H) through two ways: fever clinic (called opportunistic confirmed) and population wide nucleic acid screening. As we mentioned in the introduction, almost all the infections can be detected through the nucleic acid screening and then isolated, hence the asymptomatic compartment is not involved in our model. The transmission diagram is shown in Fig. 1.
Let \(\tau _i\) denote the time when the infected individuals are confirmed and isolated through the population wide nucleic acid screening for the whole society, with \(i=1,2,\cdots\) being the times implementing screening. We assume that this part of infected individuals can be diagnosed simultaneously, and immediately move to compartment H at \(\tau _i^+\). Based on the above assumptions and previous studies, the transmission dynamics is governed by the following impulsive model:
Considering the enhanced control interventions by the government, we introduced four timedependent parameters into the system (i.e. model (1)). With the implementation of contact tracing, a proportion of q of individuals exposed to the virus is quarantined. Let \(\beta\) be the transmission probability and c be the contact rate, then the quarantined individuals can move into compartment \(E_q\) (or \(S_q\)) at a rate of \(\beta cq\) (or \((1\beta )cq\)) if they are effectively infected (or not effectively infected). While the other proportion, \(1q\), missed from the contact tracing, will move into the exposed compartment E at a rate of \(\beta c(1q)\) once effectively infected or stay in the susceptible compartment S otherwise. Moreover, as the implementation of lockdown, the contact rate c(t) should be a decreasing function of time t with the following form [27, 28]
where \(c_0\) is the contact rate before lockdown (i.e., before 2021/7/31 in Yangzhou and 2021/12/23 in Xi’an, respectively), \(c_b\) is the minimum contact rate, and \(r_c\) is the corresponding exponential decreasing rate of the contact rate. Similarly, we set the quarantined rate q(t) to be an increasing function of time t after lockdown, which takes the following form
where \(q_0\) is the quarantine rate before lockdown, \(q_b\) is the maximum quarantine rate, and \(r_q\) is the exponential increasing rate of the quarantine rate.
In addition, the rate of diagnosis of infected individuals and quarantined exposed individuals in Xi’an are set to be piecewise functions [29], and the basis of function segmentation is consistent with the time of lockdown strategy [30], which are given by
and
We note that \(c_i\) denotes the number of diagnosed individuals due to the ith time population wide pooled nucleic acid testing during the local outbreaks. Then the term \(\text{ min }\{I(\tau _i),c_i\}\) describes the newly confirmed and hospitalized individuals through the nucleic acid screening at time \(\tau _i^+\), and the term \(\text{ max }\{I(\tau _i)c_i,0\}\) denotes the number of infected individuals who still remain in compartment I at time \(\tau _i^+\).
Data
We obtained the epidemic data related to the local COVID19 outbreaks in Yangzhou city of Jiangsu Province of China during 2021/7/28 to 2021/8/26 and Xi’an city of Shaanxi Province of China during 2021/12/9 to 2022/1/9 from the Health Commission of Jiangsu Province and the Health Commission of Shaanxi Province [31, 32], respectively. The data information includes the time series of the daily opportunistically diagnosed cases, daily confirmed cases from the quarantined population, and the daily diagnosed cases through the population wide nucleic acid screening, as shown in Fig. 2. It should be mentioned that the infected individuals can be diagnosed through opportunistic diagnosis such as fever clinic and population wide nucleic acid screening.
Model fitting process
When fitting model (1) to the epidemic data of Yangzhou and Xi’an, we initially informated several parameters and initial conditions from the database and literuture [33,34,35,36]. This includes latent period (\(1/\sigma\)), the releasing rate of quarantined susceptible individuals (\(\lambda\)), initial susceptible (S(0)), quarantined susceptible (\(S_q(0)\)), quarantined exposed (\(E_q(0)\)), confirmed and isolated (H(0)), and recovered (R(0)) populations, as listed in Table 1. As mentioned in the introduction, the COVID19 epidemic pattern in China is that each outbreak is an independent outbreak with a starting time and end time of confirm cases. Incorporating the public data, we set the initial conditions of \(S_q(0)\), \(E_q(0)\), H(0) and R(0) as 0, 0, 2, 0 in Yangzhou, and as 0, 0, 1, 0 in Xi’an, respectively.
By assuming that each data point follows a Poisson distribution with a mean assumed as the observed counts, we sample 200 data of each time point from the corresponding distributions, consequently, we obtain 200 time series data of daily opportunistically confirmed cases, daily confirmed cases from quarantined population and the cumulative confirmed cases [29, 37]. Then, for each sampled time series, we use the least squared method to simultaneously fit the daily opportunistically confirmed cases, daily confirmed cases from quarantined population and the cumulative confirmed cases, where a priori distribution for each unknown parameter and initial condition is provided. Here, the ODE system is solved by the “ODE45” function while the “fmincon” function is used to search the optimal solutions of the objective function. The objective function is defined as the residual sum of squares between the real data of the time series of the opportunistic daily confirmed cases, daily confirmed cases from the quarantined population and the cumulative confirmed cases and the correspondingly predicted numbers by solving system (1). Therefore, we obtain 200 outputs of the unknown parameters and initial conditions, from which we can obtain their means. Also, we then calculate the 2.5% and 97.5% percentiles of the 200 solutions to generate the 95% confidence interval (CI) for quantile of the fitting results.
Results
Model fitting results and senario analysis
By fitting model (1) to daily opportunistically confirmed cases, daily confirmed cases from quarantined population and the cumulative confirmed cases simultaneously, we showed the fitting results of the local outbreaks in Xi’an and Yangzhou in Fig. 3. Based on the fitting results, we obtained the mean and the standard diviation of the unknown parameters and initial conditions, as listed in Table 1. Further, to test if the key parameters correlate to each other, we yield the symmetric covariance matrix of seven parameters \(\beta\), \(c_0\), \(q_0\), \(\gamma _I\), \(\delta _{I_0}\), \(\delta _{q_0}\) and \(\gamma _H\) estimated by fitting the data of Xi’an in Table 2. It follows from Table 2 that the estimated values of the initial contact rate (\(c_0\)) and the recovery rate \(\gamma _I\) are of the biggest positive correlation coefficient. It highlights the importance to use the least squared methods with a priori distribution to each parameter for including the useful information of the parameters and to fix them into a reasonable range.
In this subsection, we investigate the impact of population wide nucleic acid screening on controlling the outbreak of COVID19 and further discuss how the frequency and time of implementing population wide nucleic acid screening affect the spread of COVID19 and the case clearing time. Firstly, we compared the previous screening scheme and no screening administered in Yangzhou and Xi’an, and plotted the daily number of opportunistic confirmed cases, the confirmed cases from \(E_q\) and the cumulative number of confirmed cases under the above two cases, shown in Fig. 4. It follows from Fig. 4 that compared to the previous screening strategy, the daily confirmed number of infected cases from I through opportunistic diagnosis and quarantined exposed cases from \(E_q\), and the cumulative number of confirmed cases remarkably increase both in Yangzhou and Xi’an without population wide nucleic acid screening strategy. To be specific, the cumulative number of confirmed cases will increase by \(77.7\%\) in Yangzhou and \(62.2\%\) in Xi’an without implementing population wide nucleic acid screening, shown in Table 3. These results indicate that previous screening strategy can effectively reduce the final epidemic size of infections, which play an important role in containing the spread of COVID19.
Then we focus on discussing the effectiveness of three different screening schemes: In the first scheme, we assume that the first half screenings of the previous screening strategy are remained while the last half screenings will not be carried out for Yangzhou (i.e. the screenings before 2022/1/1 of the previous screening strategy are remained while the screenings after 2022/1/1 will not be carried out for Xi’an); In the second scheme, we assume that one cycle screening is implemented every other two days; In the third scheme, we remain the last half screenings of the previous screening strategy while the first half screenings will not be administered for Yangzhou (i.e. remaining screenings after 2022/1/1 of the previous screening strategy while the screenings before 2022/1/1 will not be administered for Xi’an). Under different screening schemes, we illustrate the impacts of various screening strategies on the COVID19 epidemic in Yangzhou and Xi’an in Fig. 5. Among the above three screening strategy, we find that the first scheme results in the most significant decrease of the cumulative number of confirmed cases while the third scheme has little effect on reducing the cumulative number of confirmed cases both in Yangzhou and Xi’an. We further compared the cumulative number of confirmed cases according to different screening schemes in Table 3. It can be observed that when reducing the times of population wide nucleic acid screening, the earlier the nucleic acid screening, the impact on controlling the spread of COVID19 is more remarkable. However, the nucleic acid screening in the late stage of the epidemic is of little help to mitigate the outbreak of COVID19.
The above analysis quantitatively showed how the population wide necleic acid screening can help to stop the spread of COVID19 when we aimed at controlling the cases into zero quickly. In the next section, we focus on discussing the role of nucleic acid screening in mitigating the COVID19 when we aimed at controlling the epidemics in a low level. Correspondingly, we assume that the lockdown strategy would be released.
Impacts of the screening on mitigating COVID19 outbreaks and the runs on medical resource
In this section, we focus on investigating how alternative population wide screening program affects the transmission dynamics of COVID19, especially, analysing the impacts of the screenings on the runs on medical resources. To this end, we assume to screening the population using nucleic acid test at a constant rate \(q_s\). That is, a constant ratio of infected population (class I) will be tested at the screening time. In addition, we set \(d_e\) as the detection efficiency rate of the screening, hence, \(d_eq_s\) denotes the effective screening rate. The periodical screening strategy is considered with a period of \(T_s\) (days), that is, every \(T_s\) days, we carry out one cycle of the population wide nucleic acid screening with the ratio \(d_eq_s\) of infected individuals being confirmed and isolated.
As we mentioned above, the main purpose in this section is to analyze the influence of screening in mitigating the COVID19 epidemics, especially in avoiding the runs on medical resources. That is to say, we would not tend to adopt such a strict control strategy by locking down the city. Correspondingly, the contact rate will not decreased to \(c_b\), instead, it can vary with the range \([c_b,c_0]\). The other parameter values are fixed as those estimated from the outbreak of Xi’an, as listed in Table 1.
Based on the above assumptions, by varying the contact rate and the screening rate, we plotted the infected population (I(t)) and confirmed and isolated population (H(t)) in Fig. 6. In Fig. 6(AB), fixing the parameters and initial conditions as the same as those related to the outbreak of Xi’an and increasing the contact rate to \(0.4c_0\) after Dec. 23, 2021, we present the impacts of the intensive screening program (one cycle every two days) on the transmission dynamics. In this circumstance, the outbreak in Xi’an is under control as the infected cases persistently decrease to a very low level in the 60 days. As we further increase the contact rate to \(0.6c_0\), if the screening rate is high enough (\(q_s=0.75\) or \(q_s=0.9\)), the outbreak is still under control. However, the screening program seems not to be able to control the outbreak if the screening rate is only 0.6 (the black curves in Fig. 6(CD)).
In Fig. 7, in addition to release the strict NPIs by setting a higher contact rate, we also released the intensive screening by fixing a lower screening frequency. Definitely, there can be a large outbreak in this case, hence we move our focus to avoid the runs on medical resource instead of control the cases into zero. With such a purpose, we should control the peak of the confirmed and isolated population, with a ratio of them being hospitalized, below a critical value. In this sense, we observed the paradox phenomenon of the screening rate in Fig. 7(AB), that is, the peak value of H(t) for \(q_s=0.3\) is higher than those for \(q_s=0\) and \(q_s=0.9\). In Fig. 7(E), we showed in more details that when the screening rate is small, the peak value of the hospitalized population increases as the screening rate increases. This indicates that the screening program may aggravate the runs on medical resources when the screening rate in a certain range. Comparing Fig. 7(AB) with (CD), we find that when the contact rate is bigger, the screening program can lead to a higher peak value of hospitalized population in a larger range of the screening rate, which can be also clearly seen from Fig. 7(E). Moreover, when decrease the detection efficiency \(d_e\) from 1 to 0.95 in Fig. 7, we find that the soluations of I(t), H(t) are very close to each other by comparing the dashed curves and the solid curves with a same color. Also, the change rule of peak value of H(t) with respect to parameter \(q_s\) seems to have a similar trend, particularly in terms of the existence of the paradox. This implies that the results are robust with a small change to parameter \(d_e\).
Figure 7 shows an important fact that there exists the potential risk of runs on medical resources to use the control measure of population wide screening when the COVID19 outbreaks are in high epidemic level. For this reason, instead of testing the whole population each time, an alternative strategy by testing a relatively small population each time but with a high screening frequency may be a better choice to balance the control of COVID19 epidemics and the runs on medical resources. Note that, due to the testing capacity, the screening ratio \(q_s\) is usually negatively correlated to the screening frequency, that is, a longer sreening period \(T_s\) means a higher screening during each screening period and vice versa. For this reason, we assume that the screening period \(T_s\) is proportional to the screening rate \(q_s\) with a constant coefficient. And we consider different combinations of \(T_s\) and \(q_s\) under this assumption, i.e. initially set \(T_s=10\) and \(q_s=0.2\), then double \(T_s\) to \(T_s=20\) and further to \(T_s=40\), correspondingly double \(q_s\) to \(q_s=0.4\) and futher to \(q_s=0.8\). The corresponding results of I(t) and H(t) are plotted in Fig. 8. It’s easy to see from Fig. 8 that the peak value of the hospitalized population can be smaller for the combination of a shorter screening period \(T_s\) and a lower screening ratio \(q_s\). This also supports that to avoid the runs on medical resource, we should try to adopt a screening program which screens a lower population each time with a higher frequency when the outbreak is in a high epidemic level already.
Discussion and conclusion
In this study, we proposed a mathematical model with pulse populationwide nucleic acid screening, which is a main part of the zeroCOVID policy (i.e. a policy to control the cases into zero in a short period). The model was calibrated by fitting the epidemic data related to the local COVID19 outbreaks in Yangzhou city of Jiangsu province and Xi’an city of Shaanxi province, China. With the proposed models, we focus on investigating two essential issues: a) the role of the largescale screening program in the zeroCOVID policy; b) the impact of the largescale screening on the transmission dynamics of COVID19 and the runs on medical resources.
To analyze the impact of the largescale population wide nucleic acid screening on reaching the goal of the zeroCOVID policy, we conducted senario analysis by deleting part times of the real screening program in Yangzhou and Xi’an. In line with the goal of the zeroCOVID policy, we compared the final size and the time that the reported cases reach a certain low level (threshold value of one case for example) in different schedules of the screening program. We showed that the cumulative confirmed cases increase by \(77.7\%\) in Yangzhou and \(62.2\%\) in Xi’an if we totally cancel the screening program during the outbreaks. In addition, cancelling the large scale screening program can postpone the time of reopening for more than one month. Considering to delete a fixed times of the screening in Yangzhou and Xi’an, we found that the screening, at the time period that the outbreak was controlled to a low level or in the decreasing phase already, has limited influence on the epidemics. In other words, screening earlier could be more costeffectiveness to balance the control of the outbreak and the costs of the testing.
Considering a periodic screening strategy, we extended our original model to a mathematical model with pulse screening at the fixed time points. The result, from the simulations by keeping the intensive screening and increasing the contact rate, shows that the intensive screening program can provide a certain space to release the highly strong NPIs. In other words, given the costeffectiveness, we may not need to lockdown the city even though the goal is to control the cases into zero in a short period. On the other hand, if we switch our goal to mitigate the epidemics instead of stopping its spread by the control interventions, we could release the highly strict interventions, including the NPIs and the intensive screening. It’s easy to image that the largescale screening will lead to a large ratio of infected population be diagnosed instantaneously, consequently, a large population need to be hospitalized when the new wave is in a high epidemic level. Therefore, there exists the potential risk of the runs on medical resource induced by the largescale screening. This is true and we observed the paradox phenomenon of the largescale screening in the sense of runs on medical resource, as shown in Fig. 7. The driver of the paradox phenomenon should be that the screening is not strong enough, corresponding to a limited screening rate, to drop the epidemic into a low level. This indicates that the populationwide nucleic acid screening can have very limited effects on controlling the epidemics of COVID19 if the outbreak is in a relatively high epidemic level or there has already been runs on medical resources. Given the runs of the medical resource, we further showed that if we try to screen a fixed number of the population, a higher frequency by screening a smaller population in each time can be a better choice to balance the control of the epidemics and the runs on the medical resources. It should be noticed that our study bases on the whole population nucleic acid tests of the COVID19 local outbreaks in China. For lowincome countries or those countries with limited test resources, optimal sample selection strategy with selecting a proportion of the whole population and excluding the samples of individuals who have a low probability of being infected (recently tested ones) should be explored, which can be further studied in the future work.
Availability of data and materials
Data and relevant code for this research work are stored in GitHub: https://github.com/Qian199210/nucleicacidscreening.git.
Abbreviations
 SARSCoV2:

Severe Acute Respiratory Syndrome coronavirus type 2
 COVID19:

Coronavirus disease 2019
 NPIs:

Nonpharmaceutical interventions
 ODE:

Ordinary differential equation
 CI:

Confidence interval
References
Kissler SM, Tedijanto C, Goldstein E, Grad YH, Lipsitch M. Projecting the transmission dynamics of SARSCoV2 through the postpandemic period. Science. 2020;368:860–8.
Leung K, Wu JT, Liu D, Leung GM. Firstwave COVID19 transmissibility and severity in China outside Hubei after control measures, and secondwave scenario planning: a modelling impact assessment. Lancet. 2020;395:1382–93.
World Health Organization. 2022 WHO Coronavirus Disease (COVID19) Dashboard. 2022. https://covid19.who.int/.
Wang ZH, Jin YQ, Jin X, Lu YF, Yu XP, Li LJ, Zhang YM. Preliminary Assessment of Chinese Strategy in Controlling Reemergent Local Outbreak of COVID19. Front Public Health. 2021;9:650672.
Lv X, Hui HW, Liu FF, Bai YL. Stability and optimal control strategies for a novel epidemic model of COVID19. Nonlinear Dyn. 2021;106:1491–507.
Xing YH, Wong GWK, Ni W, Hu XW, Xing QS. Rapid response to an outbreak in Qingdao. China N Engl J Med. 2020;383:e129.
Wang LY, Zhang Q, Liu JJ. On the dynamical model for COVID19 with vaccination and timedelay effects: A model analysis supported by Yangzhou epidemic in 2021. Appl Math Lett. 2022;125:107783.
National Health Commission of the People’s Republic of China. 2022. http://www.nhc.gov.cn/xcs/yqfkdt/gzbd_index.shtml.
Regen F, Eren N, Heuser I, HellmannRegen J. A simple approach to optimum pool size for pooled SARSCoV2 testing. Int J Infect Dis. 2020;100:324–6.
Majid F, Omer SB, Khwaja AI. Optimising SARSCoV2 pooled testing for lowresource settings. Lancet Microbe. 2020;1:101–2.
Zhou D, Zhou M. Mathematical Model and Optimization Methods of WideScale Pooled Sample Testing for COVID19. Mathematics. 2022;10(7):1183.
Lohse S, Pfuhl T, BerkóGöttel B, et al. Pooling of samples for testing for SARSCoV2 in asymptomatic people. Lancet Infect Dis. 2020;20(11):1231–2.
Li ZJ, Liu FF, Cui JZ, et al. Comprehensive largescale nucleic acidtesting strategies support China’s sustained containment of COVID19. Nat Med. 2021;27:740–2.
Hu ZL, Song C, Xu CJ, et al. Clinical characteristics of 24 asymptomatic infections with COVID19 screened among close contacts in Nanjing. China Sci China Life Sci. 2020;63:706–11.
Long QX, Tang XJ, Shi QL, et al. Clinical and immunological assessment of asymptomatic SARSCoV2 infections. Nat Med. 2020;26:1200–4.
Hernandez X, Valentinotti S. On an optimal testing strategy for workplace settings operating during the COVID19 pandemic. PLoS ONE. 2022;17(3):e0264060.
Han X, Li J, Chen Y, et al. SARSCoV2 nucleic acid testing is China’s key pillar of COVID19 containment. Lancet. 2022;399(10336):1690–1.
Zhu W, Zhu Y, Wen Z, et al. Quantitative assessment of the effects of massive nucleic acid testing in controlling a COVID19 outbreak. BMC Infect Dis. 2022;22(1):845.
Hsiang S, Allen D, AnnanPhan S, et al. The effect of largescale anticontagion policies on the COVID19 pandemic. Nature. 2020;584(7820):262–7.
Gatto M, Bertuzzo E, Mari L, et al. Spread and dynamics of the COVID19 epidemic in Italy: effects of emergency containment measures. PNAS. 2020;117(19):10484–91.
Wagner CE, SaadRoy CM, Grenfell BT. Modelling vaccination strategies for COVID19. Nat Rev Immunol. 2022;22:139–41.
Watson OJ, Barnsley J, Toor J, Hogan AB, Winskill P, Ghani AC. Global impact of the first year of COVID19 vaccination: a mathematical modelling study. Lancet Infect Dis. 2022;22(9):1293–302.
Shulgin B, Stone L, Agur Z. Pulse vaccination strategy in the SIR epidemic model. Bull Math Biol. 1998;60:1123–48.
Lu ZH, Chi XB, Chen LS. The effect of constant and pulse vaccination on SIR epidemic model with horizontal and vertical transmission. Math Comput Model. 2002;36(9–10):1039–57.
Yang YP, Xiao YN, Wu JH. Pulse HIV vaccination: Feasibility for virus eradication and optimal vaccination schedule. Bull Math Biol. 2013;75:725–51.
Tang SY, Pang WH. On the continuity of the function describing the times of meeting impulsive set and its application. Math Biosci Eng. 2017;14:1399–406.
Tang B, Bragazzi NL, Li Q, et al. An updated estimation of the risk of transmission of the novel coronavirus (2019nCov). Infect Dis Model. 2020;5:248–55.
Tang B, Xia F, Tang SY, et al. The effectiveness of quarantine and isolation determine the trend of the COVID19 epidemics in the final phase of the current outbreak in China. Int J Infect Dis. 2020;95:288–93.
Tang B, Xia F, Bragazzi NL, et al. Lessons drawn from China and South Korea for managing COVID19 epidemic: Insights from a comparative modeling study. ISA Trans. 2022;124:164–75.
Press conference on the prevention and control of COVID19 in Xi’an. 2022. http://live.cnwest.com/live/4496.html.
Jiangsu Commission of Health. 2022. http://wjw.yangzhou.gov.cn/yzwshjh/index.shtml.
Shaanxi Commission of Health. 2022. http://sxwjw.shaanxi.gov.cn/sy/wjyw/.
The data of the seventh national census. 2022. https://www.hongheiku.com/category/shijirenkou.
Tang B, Wang X, Li Q, Bragazzi NL, Tang SY, Xiao YN, Wu JH. Estimation of the Transmission Risk of the 2019nCoV and Its Implication for Public Health Interventions. J Clin Med. 2020;9:462.
Li Q, Tang B, Bragazzi NL, Xiao YN, Wu JH. Modeling the impact of mass influenza vaccination and public health interventions on COVID19 epidemics with limited detection capability. Math Biosci. 2020;325:108378.
Guan WJ, Ni ZY, Hu Y, et al. Clinical Characteristics of Coronavirus Disease 2019 in China. N Engl J Med. 2020;382:1708–20.
Chowell G, Luo R. Ensemble bootstrap methodology for forecasting dynamic growth processes using differential equations: application to epidemic outbreaks. BMC Med Res Methodol. 2021;21:34.
Acknowledgements
Not applicable.
Funding
This research was funded by the National Natural Science Foundation of China (grant number: 12101488 (BT), 12201493 (QL)) and the Natural Science Basic Research Program of Shaanxi Province (grant numbers: 2022JQ050(QL), 2022JM027(QL)). The funders played no role in study design, data collection, analysis, and the writing of the manuscript.
Author information
Authors and Affiliations
Contributions
Conceptualization, QL, BT; methodology, QL; validation and simulation, QL, BT; data curation, QL, YB; writingoriginal draft preparation, QL; writingreview and editing, YB, BT; All authors have read and agreed to the published version of the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Li, Q., Bai, Y. & Tang, B. Modelling the pulse populationwide nucleic acid screening in mitigating and stopping COVID19 outbreaks in China. BMC Infect Dis 23, 280 (2023). https://doi.org/10.1186/s12879023082651
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12879023082651
Keywords
 COVID19
 Nucleic acid screening
 Runs on medical resource
 Screening paradox
 Mathematical model