Abstract
Since the start of the COVID-19 pandemic on December 31st, 2019, with the World Health Organization being notified of pneumonia of unknown cause in Wuhan (China), Taiwan has successfully ended two COVID-19 community outbreaks. For 19 days, the third community outbreak has now been successfully suppressed, putting Taiwan on path to end it too around Aug. 16th based on our forecast using an exponential model. Since May 28th the 7-day average of reported confirmed infected, which peaked at 593, has been falling to 204 on June 16th and the 7-day average of reported suspected and excluded cases increased to above 25 000. Resulting in a decrease in the ratio of the 7-day average of local & unknown confirmed to suspected cases–the identified control variable–to less than one third of its peak value. The later is a hallmark of working contact tracing, which together with testing and isolation of infected are the keys to ending the community outbreak.
Introduction
The first outbreak started a little before February 16th and ended April 11th, 2020 (54 days), see Figure 1a. The second outbreak started a little before January 12th and ended February 9th, 2021 (27 days). The third outbreak started a little before April 20th and is still ongoing (56 days), but shows clear signs of having peaked on May 28th. The record of 721 confirmed locally or unknown acquired infections in a day was reported on May 22nd. Between the first and second outbreak Taiwan saw a period with 273 days without any reported local infection, excluding an isolated case of unknown source reported August 2nd and a local case reported December 22nd. Between the second and third outbreak Taiwan had 69 days without any reported local infection.
Reported number of local and unknown infected and 7-day moving average in Taiwan from the first local case 2020-01-28 until 2021-06-16 (a). The three community outbreaks are marked in blue (1st), green (2nd), and red (3rd). The cases between are not connected. The ratio of local and unknown confirmed to suspected cases versus ratio of excluded to suspected cases reported for the three outbreaks (b-d). Number of local and unknown confirmed cases versus number of suspected cases reported for the three outbreaks (e-g). The 7-day moving average of each quantity is used in b-g.
Results and Discussion
From June 2020 until the current outbreak, life in Taiwan was pre-pandemic, except for the mandatory mask and temperature checks in the public transportation introduced on March 31st1 and at large events. Large events with tens of thousands of visitors have been held, such as the Mayday New Year’s Eve pop concert at Taoyuan International Baseball Stadium on December 31st2. Both the 2nd and 3rd outbreak involved the alpha variant (B.1.1.7), which Hsu et al.3 estimated to have a 1.44-fold higher infection probability and 57% higher basic reproduction number based on household transmissions during the 1st and 2nd outbreak, in agreement with previous estimates of 43-90%4. Evidently, the measures in place to reduce transmissibility and prevent a large scale outbreak were not sufficient, so Taiwan should have had a surge in cases if community transmission existed during the two periods without local infections. In addition, the National Health Insurance Administration (NHIA) has proactively been seeking out patients with severe respiratory symptoms in its database and allowing all hospitals, clinics, and pharmacies to see the patients’ travel history obtained from the National Immigration Agency from February 18th, 2020 onwards5. The current outbreak has been traced to a cluster of infected China Airline pilots and the Novotel at Taoyuan International Airport violating the quarantine rules by housing quarantined flight crews and local guests in the same building6,7. This initial failure of the mandatory quarantine at the boarder and contact tracing during the current outbreak combined with the success of contact tracing during the 2nd outbreak and the transmissibility reduction measures remaining the same, implies that Taiwan’s success is due to border control and tracing of contacts upon suspicion. Thus we can finally with confidence say that Taiwan has ended two community outbreaks thanks to contact tracing, testing, and isolation.
The Taiwanese COVID-19 control strategy implemented by the Central Epidemic Command Center (CECC), activated on January 20th, 20205, is based on six main pilars: (a) border control with quarantine upon arrival, (b) self-health monitoring when having visited a place with known cases, (c) testing when showing symptoms and seeking medical care, (d) mandatory supervised quarantine of confirmed infected and individuals at high-risk of having been infected, (e) contact tracing, and (f) a four-level system of measures to suppress community spreading. These are devised to control the effective reproduction number, i.e. the expected number of people an infected individual will transmit the disease to while infectious. The reproduction number can be seen as the transmission risk per contact (transmissibility), times the number of contacts per day, times the number of days the person is infectious. From June 7th until May 11th no restrictions on the size of gatherings, i.e. curbing of the number of contacts, existed8,7. The current level 3 measures in place since May 19th, consist of mandatory wearing of masks at all times outside private spaces and social distancing, which reduce the transmissibility; indoor gatherings limited to five people and closure of certain businesses, partly including schools and preschools, which reduce the number of contacts; and mandatory COVID-19 testing in areas where community transmission has taken place, which reduces the number of days an infected person can transmit the disease before being quarantined9,10. As demonstrated by numerous countries, such as the United Kingdom and United States, community suppression measures, such as lockdowns, alone do not end a community outbreak. As soon as these non-pharmacological interventions (NPIs) are lifted the reproduction number will increase again, so why do we claim that Taiwan is on track to ending this community outbreak?
In short, the ratio of the 7-day moving average of confirmed infected to suspected has already fallen to less than one third of the peak value of 0.029 on May 28th, while the number of daily suspected has increased to above 25 000, and the ratio of the 7-day moving average of excluded to suspected is slightly above one, which is needed to resolve the backlog from when it was down at 0.37. In other words, the contact tracing, testing, and isolation of infected has, after initial failure and challenges, successfully been scaled up to the required volume of around 25 000 suspected per day and the number of confirmed infected is decreasing. The same pattern has been repeated during all three outbreaks: first the number of confirmed infected increase rapidly and the contact tracing, testing, and exclusion of cases lags, then focus on contact tracing brings up the number of suspected and investment into expanding testing clears the backlog, see Figure 1b-g. Thus the key to ending an outbreak is to focus on contact tracing to bring up the number of suspected and increase testing and isolation/care capacity. The focus on contact tracing and bringing up the number of suspected is key because one need to prevent every infected person from infecting a new person to end the outbreak. It makes the ratio of confirmed infected to suspected the essential metric to follow. The Lancet Commission Task Force on public health measures to suppress the pandemic also report on the effectiveness of contact tracing and quarantine measures, albeit not on its use to end outbreaks11. Since May 19th, Taiwan has implemented an SMS text message contact tracing system based on each person upon entry in each business scanning a place specific QR code that generates a unique SMS sent to 1922–the Taiwan Center for Disease Control hotline12, which helps scale the contact tracing. This data tells us that the contact tracing, testing, and isolation capacity has been scaled to the required level. Therefore, it is not a question of if the outbreak can be ended but a question of when it will be ended. It is a matter of keeping up the work until this outbreak also is brought to an end. To keep up the work requires will–will from everyone to contribute to the contact tracing, self-health monitoring, testing, and caring of each other. To have will, one need to believe that it can be done, therefore it is vital to get this message to everyone: Together we can and will also end this outbreak.
The two extensions of the level 3 measures have made some ask if vaccination is necessary to end this outbreak? Of the only three countries, with a population over 1 million, that have vaccinated over half their population, both Bahrain and Mongolia in the past week reported over 5 000 confirmed infected, while Israel reported a mere 117 infected13. Actually, both Bahrain and Mongolia have set new records in the weekly number of confirmed infected within the previous three weeks. Israel only recently, on June 1st, lifted most restrictions, including limitations on the number of people at gatherings and proof of vaccination to enter e.g. gyms, theatres, hotels, and synagogues, and have kept border controls14. In a recent simulation study, Moghadas et al.15 predict that if all unvaccinated individuals in the United States revert to pre-pandemic life on July 4th, then it would result in the worst daily incidence of COVID-19 cases despite assuming close to 70% being vaccinated, but if only the vaccinate individuals revert, then it result in a minimal increase even if they revert two weeks after the first dose. Based on the recent estimate of the basic reproduction number without NPIs in Taiwan during the second outbreak at 6.73, 85% of the Taiwanese population would need to be vaccinated to reach herd immunity. Considering the difficulties the Taiwanese government has faced in purchasing COVID-19 vaccines, it is unlikely that Taiwan could implement a population scale vaccination program and reach herd immunity any faster than other nations. A population wide vaccination campaign is likely to take a year before it could end this community outbreak, while the proven method of contact tracing, testing, and isolation is likely to end it around Aug. 16th. We forecast this end date based on fitting an exponential model to the daily ratios of the 7-day moving average of confirmed infected to suspected from the peak of 0.029 on May 28th to 0.0077 on June 16th. The calculated non-simultaneous prediction bounds on a new observation with 95% confidence puts the earliest end date on July 2nd and the latest indefinitely into the future, see Figure 2. The predicted half-time is 9.9 days (95% confidence interval: 9.1-10.8). We therefore encourage all Taiwanese to be hopeful and remain vigilant; the end around Aug. 16th is only two months away.
Evolution of the ratio of the 7-day moving mean of local & unknown confirmed to the 7-day moving mean of suspected cases and our forecast of the end of the third outbreak. The three community outbreaks are marked in blue (1st), green (2nd), and red (3rd). The cases between in grey are not connected.
Conclusion
As soon as this outbreak is ended, we hope the world would take notice and learn from the Taiwanese strategy, so unnecessary suffering and deaths can be avoided. Like Wang et al.5, we think “Taiwan is an example of how a society can respond quickly to a crisis and protect the interests of its citizens.”
Methods
To limit the influence of reporting differences between the week days and fluctuations in when cases are reported to the Taiwan CDC, we use the 7-day moving mean calculated over the past seven days. It is worth noting that our moving mean cover the mean serial interval, which was estimated to 6.2 days among the 42 infector-infectee pairs during the 1st outbreak18.
To forecast when the current third community outbreak will end, we assumed that the de-crease in the ratio of local & unknown confirmed to suspected cases can be described as an exponential growth process. More precisely, we assumed the model
where a is the initial value and the reciprocal of b is the half-time. We fitted this model to the ratio of the 7-day moving mean of local & unknown confirmed to the 7-day moving mean of suspected cases from the peak on 2021-05-28 to 2021-06-16, i.e. 19 time points, using the Curve Fitting Toolbox in Matlab R2015b by MathWorks using the ”NonlinearLeastSquares” estimator. The root mean squared error of the model fitted to the data points was 0.0011 and adjusted R-squared 0.97, indicating a good fit that explains the data well. The estimated parameters are a 0.029 (95% confidence bounds 0.028 - 0.030) and b 0.10 (95% confidence bounds 0.09 - 0.11). The later corresponds to a halving of the ratio every 9.9 days (95% confidence bounds 9.1 - 10.8 days). The previous two outbreaks have ended after the ratio has been brought down to 1.2· 10−4, see Figure 2, so we assume that the same is required to end this outbreak. Note that the length of the window of our moving mean and the baseline number of suspected affects this ratio. When only one confirmed infected remains, the nominator will be 0.14, i.e. 1/7, and the denominator then need to be around 1200, i.e. on average 1200 suspected each day for the last confirmed infected. The mean 7-day moving mean number of suspected between the 1st and 2nd outbreak was 324 and between the 2nd and 3rd outbreak it was 684. By extrapolation we forecast that this ratio would be reached on 2021-08-16 after required six halvings. This forecast depends on the contact tracing, testing, and isolation being equally effective as during the past three weeks and no new local clusters emerging from people infected abroad. The calculated non-simultaneous prediction bounds on a new observation with 95% confidence, drawn as black dashed lines in Figure 2, puts the earliest end date on July 2nd and the latest indefinitely into the future. The later in agreement with a failure to execute the contact tracing, testing, and isolation or new local clusters emerging from people infected abroad.
Funding
This study was supported by Ministry of Science and Technology, Taiwan (MOST 109-2224-E-006-003, 109-2740-B-006-001). The funding has covered a stipend to Yu-Heng Wu and general lab expenses. The funding institution has not been involved in this study.
Disclosure statement
The authors have nothing to disclose.
Declaration of competing interest
The authors have no conflicts of interest relevant to this article.
Contribution statement
TEMN did conceptualisation, funding acquisition, methodology, project administration, re-sources, supervision, and writing of the original draft. YHW did investigation, data curation, formal analysis, and software development. Both did data validation, visualisation, and review & editing of the final manuscript.
Data
All COVID-19 data were collected from the Taiwan Centers for Disease Control (Taiwan CDC), which provided an online platform for downloading datasets, including respiratory syndrome coronavirus-2 (SARS-CoV-2) infection data16. Both data downloaded from the online platform and manually extracted from their News Bulletins17 were used. All data used is for the readers convenience also included in the Appendix.
Data collection
All data is publicly available and was collected from the Taiwan Centers for Disease Control (Taiwan CDC). We collected the cumulative number of confirmed infected with SARS-CoV-2 with a know local source and unknown source, as well as the cumulative reported suspected and excluded cases, from the Taiwan CDC News Bulletins, available in Chinese at https://www.cdc.gov.tw/Bulletin/List/MmgtpeidAR5Ooai4-fgHzQ (last visited 2021-06-17). Unfortunately, these news bulletins have not been issued every day, nor at the same time of the day through out the pandemic, which means that the daily increase derived by subtracting the value in the previous report is influenced by the varying time between these reports. The way of counting the number of suspected and excluded was changed on 2020-03-06, which contain a large correction to these cumulative numbers. We therefore also downloaded the ”Daily Number of Cases Suspected SARS -CoV-2 Infection Tested” dataset from the Taiwan CDC Open Data portal https://data.cdc.gov.tw/en/dataset/daily-cases-suspected-sars-cov-2-infection_tested (last visited 2021-06-17). This dataset contains the daily number of cases suspected to have SARS-CoV-2 from three sources: (i) notifications of infectious diseases, (ii) home quarantine and inspection, and (iii) expanded monitoring, as well as the total daily number.
Data cleaning and preprocessing
We first removed the first CDC News Bulletin for the days when more than one was issued, thus keeping the later report with the larger cumulative numbers. The we corrected for the number of days in between the report for all cases with more than one day in between by evenly distributing the change in cumulative number over the days in between and ensuring that the cumulative sum equals the reported one. Due to the change in the way of counting the number of suspected and excluded on 2020-03-06, the number of suspected cases is unreasonably low before it and a jump correcting for the change occur on that day. These artefacts are clearly seen in Figure 3 as a dip in the number of suspected before and a value slightly above 10 000 on this day. These artefacts clearly influence the 7-day moving mean of the number of suspected that we need for getting a stable estimate of the ratio of confirmed infected to suspected cases. The total daily number of suspected cases in the CDC Open Data portals dataset does not contain this error, so we decided to use it instead.
Comparison of the change in number of suspected between CDC News Bulletins (red circles), the total daily number of suspected in the CDC Open Data (blue x), and the change in number of suspected between CDC News Bulletins corrected for the number of days in between (red dot), and the 7-day moving mean of the later two (blue and red lines), and the 7-day moving mean of the change in number of excluded between CDC News Bulletins corrected for the number of days in between (magenta line), and the 7-day moving mean of our estimate of the daily number of excluded based on the total daily number of suspected in the CDC Open Data (cyan line).
Relative difference in the 7-day moving mean of the change in number of suspected between CDC News Bulletins corrected for the number of days in between and the 7-day moving mean of the total daily number of suspected in the CDC Open Data (red line), and the 7-day moving mean of the change in number of excluded between CDC News Bulletins corrected for the number of days in between and the 7-day moving mean of our estimate of the daily number of excluded based on the total daily number of suspected in the CDC Open Data (magenta line).
Since the number of excluded cases only was available from in the CDC News Bulletins and contain a similar change in the way fo counting on 2020-03-06, we assumed that the ratio of excluded to suspected was correct and used the ratio to estimate the number of daily excluded based on the total number of daily suspected from the CDC Open Data dataset. To ensure that the cumulative sum of daily excluded matched the reported values in the News Bulletins after March 6th, we proportionally to the daily suspected added or removed excluded cases from the days between two reports iteratively from the later report day to the former to avoid errors from rounding to whole integers later. This resulted in an estimate of the daily total number of excluded cases that on each day with a CDC News Bulletin report match the reported number exactly. The result of our preprocessing is visualised in Figure 3, where the number of suspected cases from the corrected News Bulletin and Open Data closely match after 2020-03-06. Note that a red circle without a dot in it indicate that more than one day had passed between the consecutive News Bulletins and the change is therefore larger than expected from one day to the next. We have verified that the cumulative number of excluded always remain smaller than the cumulative number of suspected.
Since our analysis if dependent on the 7-day moving mean of suspected and excluded we also evaluated the relative difference in the 7-day moving mean daily values obtained from the CDC News Bulletin and Open Data, relative to the later. The relative difference in number of suspected remains below 41% every day after 2020-03-06, below 10% 81% of the days, and below 1% 33% of the days. The large relative differences occur on separated days with small relative differences in between and are caused by differences in on which day extreme numbers of suspected are reported, see Figure 3. This is to be expected, since the Open Data most probably contains all suspected cases on each day, while the News Bulletins contain the cases added since the last bulletin. Thus these differences may shift values by one day, but no more. We cannot tell which of these two data sources that is more reliable, but considering that the Open Data at least is corrected for the change in the way of calculating the number of suspected we decided to use it and our estimate of the total daily number of excluded based on it. This does not affect the ratio of excluded suspected cases, but it affects the ratio of local & unknown confirmed to suspected cases by up to 41% on individual days, which in practise means that specific values may be shifted in time by one day. Since our forecast of when the third outbreak will end is based on 19 time points any effects of this difference is well within the estimated confidence intervals.
Data used
The data we used after cleaning and preprocessing is included in Table 1-4.
The data used in this study.
The data used in this study (continued).
The data used in this study (continued).
The data used in this study (continued).
Acknowledgements
The authors acknowledges valuable discussions, help with presentation, and comments from Dr. Akram Ashyani, Dipesh Dhayfule, and Jose Chang.
Appendix
Here we provide details on the data collection, cleaning, and analysis to enable verification and reproduction of the results.