Advertisement
Research Article| Volume 102, P1-9, January 2021

Spatial variability in reproduction number and doubling time across two waves of the COVID-19 pandemic in South Korea, February to July, 2020

Open AccessPublished:October 07, 2020DOI:https://doi.org/10.1016/j.ijid.2020.10.007

      Highlights

      • South Korea has experienced two spatially heterogeneous waves of COVID-19.
      • Seoul and Gyeonggi Province experienced two waves of COVID-19 — in March and June 2020.
      • In densely populated Seoul and nearby areas, reproduction numbers exceeded 3.0.
      • The easing of social distancing measures resulted in the second wave.

      Abstract

      Objectives

      In South Korea, 13 745 cases of coronavirus disease (COVID-19) had been reported as of 19 July, 2020. To examine spatiotemporal changes in the transmission potential, we aimed to present regional estimates of the doubling time and reproduction number (Rt) for COVID-19 in the country.

      Methods

      Daily series of confirmed COVID-19 cases in the most affected regions were extracted from publicly available sources. We employed established mathematical and statistical methods to investigate the time-varying reproduction numbers and doubling time for COVID-19 in Korea.

      Results

      At the regional level, Seoul and Gyeonggi Province experienced the first peak of COVID-19 in early March, followed by a second wave in early June, withRt exceeding 3.0 and mean doubling time ranging from 3.6 to 10.1 days. As of 19 July, 2020, Gyeongbuk Province and Daegu had yet to experience a second wave of the disease. During the first wave, mean Rt for these areas reached 3.5–4.4, and doubling time ranged from 2.8 to 4.6 days.

      Conclusions

      Our findings support the effectiveness of control measures against COVID-19 in Korea. However, the easing of restrictions that had been imposed by the government in May 2020 facilitated a second wave in the greater Seoul area.

      Keywords

      Introduction

      Since the first COVID-19 cases were reported in Wuhan, Hubei Province, China in December 2019, more than 24.7 million cases of coronavirus disease (
      • Coronavirus
      Coronavirus: South Korea confirms second wave of infections.
      ), including more than 830 000 related deaths, had been reported worldwide (WHO) as of August 30, 2020. In South Korea, the novel coronavirus was first diagnosed in a 36-year-old Chinese woman who entered the country on January 20, 2020. South Korea has since experienced two heterogeneous waves of the disease, with a total of 13 745 cases, including 295 deaths, as of July 19, 2020 (
      • KCDC
      The Updates of COVID-19 in Republic of Korea.
      ).
      During the early phase of the COVID-19 outbreak in South Korea, public health authorities primarily conducted strict contact tracing and isolation of confirmed cases, as well as applying quarantine for those suspected of infection with the novel coronavirus (
      • Covid-19 National Emergency Response Center E
      • Case Management Team KCfDC, Prevention
      Contact transmission of COVID-19 in South Korea: novel investigation techniques for tracing contacts.
      ). As the number of COVID-19 cases continued to increase, Korean public health authorities set the alert to the highest level (Level 4) on February 23, and mandated the population to report any symptoms related to COVID-19 for further screening and testing. In addition, the country rapidly adopted a ‘test, trace, isolate, and treat’ strategy that has been deemed effective in eliminating localized outbreaks of the novel coronavirus (
      • KCDC
      The Updates of COVID-19 in Republic of Korea.
      ). However, the total number of confirmed cases in South Korea spiked from 31 cases on February 18 to 433 on February 22. According to the Korea Centers for Disease Control and Prevention (KCDC), this sudden jump was mainly attributed to a super-spreader (the 31 st case) who had participated in a religious gathering of attendees of the Shincheonji Church of Jesus in Daegu (
      • KCDC
      The Updates of COVID-19 in Republic of Korea.
      ). These superspreading events occurred in the Daegu and Gyeongbuk provincial regions, leading to more than 5210 secondary COVID-19 cases in Korea (
      • KCDC
      The Updates of COVID-19 in Republic of Korea.
      ,
      • Ryall J.
      Coronavirus: Surge in South Korea virus cases linked to church ‘super-spreader’.
      ). These events facilitated sustained transmission chains, with 38% of the cases in the country associated with the church cluster in Daegu (
      • Shim E.
      • Tariq A.
      • Choi W.
      • Lee Y.
      • Gerardo C.
      Transmission potential and severity of COVID-19 in South Korea.
      ).
      On March 8, the KCDC announced that 79.4% of all cases had epidemiological links, while the remaining 20.6% cases were either sporadic cases or under investigation (
      • KCDC
      The Updates of COVID-19 in Republic of Korea.
      ). Case clusters started to accumulate from churches in the Seoul capital area, and on March 17, 79 church attendees developed COVID-19 after a service at the River of Grace Community Church. In spite of social distancing orders put forward by the government, some churches continued to conduct services, which led to new clusters of infection. For instance, the Manmin Central Church in Seoul was involved in one of the clusters, with 41 infections linked to a gathering in early March; SaengMyeongSu Church in Gyeonggi Province was another cluster linked to 50 cases (
      • Park C.K.
      Coronavirus cluster emerges at another South Korean church, as others press ahead with Sunday services.
      ).
      As SARS-CoV-2 infection spread rapidly outside Korea, the number of imported cases started to increase, resulting in 476 imported cases out of 9661 total cases (4.9%) as of March 30. Consequently, as of 1 April, the KCDC implemented self-quarantine measures for travellers from Europe and the USA (
      • KCDC
      The Updates of COVID-19 in Republic of Korea.
      ). In addition, incoming travellers with symptoms but negative test results for coronavirus, as well as asymptomatic short-term visitors, were ordered to follow a 2-week quarantine in the government facilities (
      • KCDC
      The Updates of COVID-19 in Republic of Korea.
      ).
      Such control measures undertaken by South Korea have been deemed successful in limiting the spread of the outbreak, without locking down entire cities (
      • Normile D.
      Coronavirus cases have dropped sharply in South Korea. What’s the secret to its success?.
      ). Therefore, after a sustained period of low incidence with fewer than 20 cases per day (April 16 to May 5), the government eased its strict nationwide social distancing guidelines on May 6, with a phased reopening of schools starting mid-May, 2020. However, a new cluster linked to nightclubs in Itaewon emerged in central Seoul in early May, resulting in a resurgence of cases that led to a second wave of COVID-19 in the greater areas of Seoul. As of May 29, the number of cases linked to this cluster had reached 266 (
      • KCDC
      The Updates of COVID-19 in Republic of Korea.
      ). Accordingly, the Seoul city government ordered all clubs, bars, and other nightlife establishments in the city to close indefinitely (
      • KCDC
      The Updates of COVID-19 in Republic of Korea.
      ). Simultaneously, another cluster emerged from an e-commerce warehouse in the Gyeonggi Province, resulting in 108 cases as of May 30.
      In the last week of May, around 40–80 daily new cases of COVID-19 were being reported (
      • KCDC
      The Updates of COVID-19 in Republic of Korea.
      ). Following this spike in the number of new COVID-19 infections — the first in nearly 2 months — public health authorities reimplemented strict lockdown measures in Seoul, along with school closures once more across the nation. In June, it was announced that the strict social distancing campaign would be extended indefinitely as a preventive measure in Seoul, Incheon, and Gyeonggi Province; however, phased reopening of schools was initiated on May 20. It was reported by the KCDC that a holiday weekend in early May triggered a new wave of infections focused in the greater Seoul area, the so-called second wave of COVID-19 in South Korea (2020). In Seoul, the average number of daily new cases reported from June 4 to June 17 was 43 (
      • KCDC
      The Updates of COVID-19 in Republic of Korea.
      ). This was followed by sporadic clusters of infection across the country in July, most of them associated with religious facilities and door-to-door salespeople, especially in the densely populated Seoul region and adjacent areas. As a result, the government banned churches from organizing small gatherings other than regular worship services from July 10 (
      • KCDC
      The Updates of COVID-19 in Republic of Korea.
      ). As of September 23, 23 216 cases of COVID-19 had been reported in South Korea, comprizing 13.4% imported cases, 59.7% cases linked to local clusters, 14.5% unlinked local cases, and 12.4% cases under investigation (
      • KCDC
      The Updates of COVID-19 in Republic of Korea.
      ).
      To estimate the regional and temporal variability in the reproduction number for COVID-19 in South Korea, including the second wave concentrated in the greater Seoul areas, we analysed the spatiotemporal progression of the epidemic in the country from mid-February to mid-July, 2020. Our focus was on estimating and interpreting the doubling time and effective reproduction number Rt, a metric that quantifies the time-dependent transmission potential of the disease, incorporating the effect of control measures, susceptible depletion, and behavioural changes. This key epidemiological parameter, Rt, represents the average number of secondary cases generated per case whenever conditions persist as they were at time t. Epidemic doubling times refer to the sequence of intervals at which the cumulative incidence doubles (
      • Lee W.
      • Hwang S.S.
      • Song I.
      • Park C.
      • Kim H.
      • Song I.K.
      • et al.
      COVID-19 in South Korea: epidemiological and spatiotemporal patterns of the spread and the role of aggressive diagnostic tests in the early phase.
      ,
      • Muniz-Rodriguez K.
      • Chowell G.
      • Cheung C.H.
      • Jia D.
      • Lai P.Y.
      • Lee Y.
      • et al.
      Doubling time of the COVID-19 epidemic by Chinese province.
      ). Therefore, an increase in the doubling time implies a decline in disease transmission. In this report, we estimated the doubling time and the effective reproduction number in relation to two epidemic waves of the COVID-19 epidemic in South Korea by employing the time series of cases by date of symptom onset for the four most affected Korean regions: Seoul, Gyeonggi Province, Gyeongbuk Province, and Daegu. We also discuss the spatiotemporal variability in the reproduction number in terms of the public health policies that were put in place by the Korean government.

      Methods

      Data

      Daily series of confirmed local COVID-19 cases in South Korea were collected from 20 January to 19 July; these were published by national and local public health authorities, including city or provincial departments of public health in South Korea (
      • KCDC
      The updates on COVID-19 in Korea as of 25 February.
      ). Our analysis focused on the regions with the highest caseloads, including Seoul, Gyeonggi Province, Gyeongbuk Province, and Daegu (Figure 1).
      Figure 1
      Figure 1Map showing the location of Seoul, Gyeonggi Province, Gyeongbuk Province, and Daegu.

      Imputing the date of onset

      For a more accurate estimation of epidemic growth rates, an epidemic curve should be analyzed according to the date of symptom onset rather than the date of reporting, because reporting delays can fluctuate substantially over the course of an epidemic. Reporting delays distort the incidence pattern of epidemics, misrepresenting the outbreak trajectory, and thus possibly affecting the estimation of reproduction number (
      • Tariq A.
      • Roosa K.
      • Mizumoto K.
      • Chowell G.
      Assessing reporting delays and the effective reproduction number: the Ebola epidemic in DRC, May 2018–January 2019.
      ). A prior study suggested that obtaining knowledge about reporting parameters, such as delay patterns and structure, improves the estimating of reproduction numbers (
      • Azmon A.
      • Faes C.
      • Hens N.
      On the estimation of the reproduction number based on misreported epidemic data.
      ). However, for the COVID-19 data in Korea, the date of symptom onset was only available for 732 cases reported in Gyeonggi Province, which yielded a mean of 4.5 days and standard deviation of 4.4 days for the distribution of delays from symptom onset to reporting of cases. Therefore, we utilized the empirical distribution of these 732 reporting delays from the onset of symptoms to reporting to impute the missing dates of onset for the remaining cases (
      • Shim E.
      • Tariq A.
      • Choi W.
      • Lee Y.
      • Chowell G.
      Transmission potential and severity of COVID-19 in South Korea.
      ). Specifically, we reconstructed 300 epidemic curves according to the date of symptom onset, from which we derived the mean incidence curve of local case incidence (
      • Shim E.
      • Tariq A.
      • Choi W.
      • Lee Y.
      • Chowell G.
      Transmission potential and severity of COVID-19 in South Korea.
      ,
      • Tariq A.
      • Roosa K.
      • Mizumoto K.
      • Chowell G.
      Assessing reporting delays and the effective reproduction number: the Ebola epidemic in DRC, May 2018–January 2019.
      ). For the calculation of Rt, the estimated mean incidence curve, based on the date of symptom onset, was used for the regions of interest (i.e., Seoul, Gyeonggi Province, Gyeongbuk Province, and Daegu) (Figure 2). Using the reconstructed mean incidence curve for local case incidence, we removed the first and last three data points to adjust for the reporting delays in our real-time analysis. We assumed that the first wave ended when the mean incidence became less than 0.2 individuals per day. Similarly, we assumed that the second wave began when the mean incidence of local cases became greater than 0.5 individuals per day. Slight variations to these thresholds did not affect our results.
      Figure 2
      Figure 2Timeline of confirmed cases of COVID-19 in Seoul, Gyeonggi Province, Gyeongbuk Province, and Daegu. The daily numbers of COVID-19 cases by date of report and by date of symptom onset are shown. The empirical distribution of reporting delays from the onset to diagnosis for 732 cases was used to impute the dates of onset for the remainder of the cases with missing data.

      Calculation of the doubling time

      We analyzed the number of times COVID-19 cumulative incidence doubled and the evolution of the doubling times in the four most affected areas in Korea (i.e. Seoul, Gyeonggi Province, Gyeongbuk Province, and Daegu) from January 20 to July19. Using regional-level daily cumulative incidence data, we calculated the times at which cumulative incidence doubled, denoted by tdi. Specifically, we assume that:
      2C(tdi)=C(tdi+1)


      where td0=0, C(td0) = C0 (i=0,1,2,3,,nd), and C(tdi) denotes the cumulative number of cases at time tdi (
      • Muniz-Rodriguez K.
      • Chowell G.
      • Cheung C.H.
      • Jia D.
      • Lai P.Y.
      • Lee Y.
      • et al.
      Doubling time of the COVID-19 epidemic by Chinese province.
      ). Here, nd is defined as the total number of times cumulative incidence doubles. Specifically, the sequence of doubling times are described as dj=Δtdi=tdi-tdi-1 where j=1,2,3,,nd. In addition, we used parametric bootstrapping with a Poisson error structure around the harmonic mean of doubling times to obtain the 95% confidence interval (
      • Chowell G.
      • Ammon C.E.
      • Hengartner N.W.
      • Hyman J.M.
      Transmission dynamics of the great influenza pandemic of 1918 in Geneva, Switzerland: Assessing the effects of hypothetical interventions.
      ,
      • Chowell G.
      • Shim E.
      • Brauer F.
      • Diaz-Duenas P.
      • Hyman J.M.
      • Castillo-Chavez C.
      Modelling the transmission dynamics of acute haemorrhagic conjunctivitis: application to the 2003 outbreak in Mexico.
      ).

      Calculation of Rt

      We assume that Rt can be estimated by the ratio of the number of new infections generated at time step t (It) to the total infectiousness of infected individuals at time t, given by s=1tIt-sws (
      • Chong K.C.
      • Zee B.C.Y.
      • Wang M.H.
      Approximate Bayesian algorithm to estimate the basic reproduction number in an influenza pandemic using arrival times of imported cases.
      ,
      • Fraser C.
      Estimating individual and household reproduction numbers in an emerging epidemic.
      ). Here, ws denotes the infectivity profile of the infected individual, which is dependent on the time since infection (s) but independent of calendar time (t) (
      • He X.
      • Lau E.H.Y.
      • Wu P.
      • Deng X.
      • Wang J.
      • Hao X.
      • et al.
      Temporal dynamics in viral shedding and transmissibility of COVID-19.
      ,
      • Wallinga J.
      • Teunis P.
      Different epidemic curves for severe acute respiratory syndrome reveal similar impacts of control measures.
      ). Specifically, ws is defined as a probability distribution describing the average infectiousness profile after infection. Individual biological factors such as pathogen shedding or symptom severity can affect the distribution ws. For example, an individual would be most infectious at time s when ws is the largest. Thus, s=1tIt-sws indicates the sum of infection incidence up to time step t–1, weighted by the infectivity function ws. Steady values of Rt above 1 indicate sustained disease transmission, whereas values less than 1 indicate that the number of new cases is expected to follow a declining trend.
      The infectivity profile, ws, can be approximated by the distribution of the generation time; however, times of infection are rarely observed, making it difficult to measure the distribution of the generation time (
      • Fraser C.
      Estimating individual and household reproduction numbers in an emerging epidemic.
      ). Therefore, the timing of symptom onset is often used to estimate the distribution of the serial interval (SI) instead, which is defined as the time interval between symptom onset in two successive cases in a chain of transmission (
      • Cori A.
      • Ferguson N.M.
      • Fraser C.
      • Cauchemez S.
      A new framework and software to estimate time-varying reproduction numbers during epidemics.
      ). Specifically, the infectiousness of a patient is a function of the time since infection and is proportional to ws if we set the timing of infection in the primary case as the time zero of ws and assume that the generation interval equals the SI. The SI was assumed to follow a gamma distribution with a mean of 4.8 days and a standard deviation of 2.3 days (
      • Nishiura H.
      • Linton N.M.
      • Akhmetzhanov A.R.
      Serial interval of novel coronavirus (COVID-19) infections.
      ). Analytical estimates of Rt were obtained within a Bayesian framework using the EpiEstim R package in R language, version 3.6.3 (R Foundation for Statistical Computing, Vienna, Austria) (
      • Cori A.
      • Ferguson N.M.
      • Fraser C.
      • Cauchemez S.
      A new framework and software to estimate time-varying reproduction numbers during epidemics.
      ). Rt was estimated at 7-day intervals, and the median and 95% credible interval (CrI) were reported.

      Results

      City of Seoul

      As of July 19, Seoul had reported a total of 1474 cases (10.7% of the total reported in South Korea), including 323 imported cases and 10 deaths, yielding an incidence rate estimated at 151 cases per million. In Seoul, the first peak based on the estimated dates of symptom onset occurred during the second week of March (March 8–14), with 18 new cases reported each day as the number of new cases linked to a Guro-gu call centre kept rising. Based on the estimated dates of symptom onset, the 7-day moving average of daily cases reached 19 cases on March 9 (Figure 3), whereas the highest value of Rt was estimated at Rt ∼ 2.9 (95% CrI: 1.6–4.7) on February 19, which continued to stay above 1 until March 6 (Figure 3).
      Figure 3
      Figure 3The epidemic trajectory of COVID-19 in Seoul as of July 19, 2020. Upper panel: The epidemic curve shows the daily number of new cases by the imputed date of symptom onset. The dates of symptom onset for cases with missing data were imputed based on the empirical distribution of delay from the onset of symptoms to reporting. Lower panel: Real-time estimates of the time-varying reproduction number (Rt) in Seoul. The solid line indicates the daily estimated Rt and the grey area indicates the 95% credible interval for Rt. The dotted line indicates the epidemic threshold of Rt = 1.
      After its first peak in February, the number of daily new cases by date of symptom onset in Seoul gradually declined, dropping below five on April 1 and remaining under five new cases per day for about a month (Figure 3). However, in early May, despite a steady decline in imported cases, locally transmitted infections surged throughout the Seoul metropolitan area, with case clusters traced to clubs, churches, and sports facilities. Therefore, Rt increased, reaching 3.0 (95% CrI: 1.6–5.0) on May 4. During the first wave, the doubling time was estimated to be 7.5 (95% CI: 7.0–8.2) days in Seoul (Table 1).
      Table 1Regional variations in doubling times in days for COVID-19 cumulative incidence, and their 95% CI, from January 20 to July 19, 2020: Seoul, Gyeonggi Province, Gyeongbuk Province, and Daegu.
      RegionMean doubling time (95% CI)
      First waveSecond wave
      Seoul7.5 (7.0–8.2)6.0 (5.4–6.7)
      Gyeonggi Province4.6 (4.2–5.5)7.5 (5.4–8.6)
      Gyeongbuk Province3.6 (3.5–4.0)10.1 (4.6–14.5)
      Daegu2.8 (2.5–4.0)10.0 (7.1–13.4)
      The number of cases continued to increase thereafter, and in the first week of June, the average daily number of confirmed COVID-19 cases in the capital surpassed the previous high point recorded in the middle of March. The major clusters in Seoul were linked to nightclubs (139 cases), the Guro-gu call centre (99 cases), Manmin Central Church (41 cases), Richway (97 cases), Yangcheon-gu table tennis club (41 cases), and Newly Planted Church in the Seoul Metropolitan Region (37 cases), as of June 18. On June 14, the average Rt in the capital dropped below 1 (95% CrI: 0.8–1.2), implying that the spread of the virus had slowed down substantially in the city (Figure 3). During the second wave, the doubling time in Seoul decreased to 6.0 (95% CI: 5.4–6.7) days, indicating faster transmission compared with that during the first wave (Table 1). As of July 15, the Rt in Seoul was estimated at 0.9 (95% CrI: 0.7–1.2), straddling the epidemic threshold of 1.0, and suggesting potential for further transmission of the virus.

      Gyeonggi Province

      Gyeonggi Province (literally meaning the ‘province surrounding Seoul’) is located in the western central region of Korea and is the most populous province in South Korea, with a population of 13.5 million people. In Gyeonggi Province, the daily number of new cases by date of symptom onset during the last weeks of February averaged 6.3 (Figure 4). Accordingly, the first peak of Rt occurred on February 22, reaching 8.9 (95% CrI: 4.8–14.2), with an estimated doubling time of 4.6 (95% CI: 4.2–5.5) days (Table 1). In the second week of March, South Korea recorded continuous drops in the number of daily new infections as large-scale testing of the followers of a religious sect in the south-eastern city of Daegu, the epicentre of COVID-19, was nearing its end; thereafter, the number of cases in Gyeonggi Province gradually decreased.
      Figure 4
      Figure 4The epidemic trajectory of COVID-19 in Gyeonggi Province, as of July 19, 2020. Upper panel: The epidemic curve shows the daily number of new cases by the imputed date of symptom onset. The dates of symptom onset for cases with missing data were imputed based on the empirical distribution of delay from the onset of symptoms to reporting. Lower panel: Real-time estimates of the time-varying reproduction number (Rt) in Gyeonggi Province. The solid line indicates the daily estimated Rt and the grey area indicates the 95% credible interval for Rt. The dotted line indicates the epidemic threshold of Rt = 1.
      However, clusters of infections in Gyeonggi Province raised concerns about further community spread, with a resurgence of cases in the province occurring in late May and resulting in the highest peak in early June. From June 1–13, an average of 14 new cases were reported each day in Gyeonggi Province. The second peak of Rt in the region occurred on May 12, with an estimated Rt value of 4.8 (95% CrI: 3.0–7.0) and the doubling time estimated at 7.5 (95% CI: 5.4–8.6) days (Table 1). Following its second peak, Rt gradually decreased (Figure 4); however, a series of sporadic clusters continued to occur. Major clusters in Gyeonggi Province included Grace River Church (67), Coupang warehouse (67), nightclubs (59), Richway (58), Uijeongbu St Mary’s Hospital (50), Guro-gu call centre/Bucheon SaengMyeongSu Church (50), door-to-door sales in the Seoul Metropolitan Region (32), and Yangcheon-gu sports facility (28). As of July 19, the number of local cases in Gyeonggi Province was 1027 (10.4% of the total reported cases in South Korea), including 29 deaths, with an Rt estimated at 0.8 (Figure 4). The incidence rate in the province was estimated at 108 per million.

      Gyeongbuk Province

      The first case in the Sincheonji cult cluster (the largest COVID-19 cluster in South Korea) appeared on February 18, resulting in sustained transmission chains, with 39% of the cases associated with the church cluster in Gyeongbuk Province. Consequently, the virus alert level was raised to ‘red’ (the highest level) on February 23, and the health authorities focused on halting the spread of the virus in Daegu and Gyeongbuk Provinces. Figure 5 shows that the peak of the epidemic occurred in the first week of March (with a reproduction number greater than 1 until March 9) (Figure 5). The doubling time in Gyeongbuk Province reached as low as 3.6 (96% CI: 3.5–4.0) days (Table 1). As of July 18, the number of cases in Gyeongbuk Province was 1393, including 54 deaths. Among these cases, 566 were related to the Shincheonji cluster. The incidence rate in Gyeongbuk Province was 523 per million, accounting for 10.2% of all confirmed cases in South Korea (
      • KCDC
      The Updates of COVID-19 in Republic of Korea.
      ). The major clusters in Gyeongbuk Province were linked to Cheongdo Daenam Hospital (119 cases), Bonghwa Pureun Nursing Home (68 cases), Gyeongsan Seo Convalescent Hospital (66 cases), pilgrimage to Israel (41 cases), Yecheon-gun (40 cases), and Gumi Elim Church (11 cases).
      Figure 5
      Figure 5The epidemic trajectory of COVID-19 in Gyeongbuk Province, as of July 19, 2020. Upper panel: The epidemic curve shows the daily number of new cases by the imputed date of symptom onset. The dates of symptom onset for cases with missing data were imputed based on the empirical distribution of delay from the onset of symptoms to reporting. Lower panel: Real-time estimates of the time-varying reproduction number (Rt) in Gyeongbuk Province. The solid line indicates the daily estimated Rt and the grey area indicates the 95% credible interval for Rt. The dotted line indicates the epidemic threshold of Rt = 1.

      City of Daegu

      The epicentre of the South Korean COVID-19 outbreak has been identified in Daegu, a city of 2.5 million people, approximately 150 miles south-east of Seoul. The rapid spread of COVID-19 in Daegu was attributed to a superspreading event in a religious group called Shincheonji, resulting in an explosive outbreak of 4511 infections in the city of Daegu, resulting in a relatively short doubling time of 2.8 (95% CI: 2.5–4.0) days (Table 1 and Figure 6). Other major clusters in Daegu included the second Mi-Ju Hospital (196 cases), Hansarang Convalescent Hospital (124 cases), Daesil Convalescent Hospital (101 cases), and Fatima Hospital (39 cases). Daegu was the most severely affected area in South Korea, with 6932 cumulative cases as of July 19, accounting for 51.0% of all confirmed cases in Korea. According to our model, the number of new cases, based on the onset of symptoms, was estimated to be the highest on February 27, with the number gradually decreasing thereafter. Accordingly, the estimated Rt was above 2 until February 27 and dropped to below 1 on March 5, although recent sporadic infections have caused Rt to fluctuate around 1 (Figure 6).
      Figure 6
      Figure 6The epidemic trajectory of coronavirus disease 2019 in Daegu, as of July 19, 2020. Upper panel: The epidemic curve shows the daily number of new cases by the imputed date of symptom onset. The dates of symptom onset for cases with missing data were imputed based on the empirical distribution of delay from the onset of symptoms to reporting. Lower panel: Real-time estimates of the time-varying reproduction number (Rt) in Daegu. The solid line indicates the daily estimated Rt and the grey area indicates the 95% credible interval for Rt. The dotted line indicates the epidemic threshold of Rt = 1.

      Discussion

      Estimates of the transmission potential of COVID-19 in Korea have displayed substantial spatiotemporal variation. Indeed, several factors influence the value of the reproduction number, including the transmissibility of an infectious agent, individual susceptibility, individual contact rates, and control measures (
      • Anderson R.
      • May R.
      Infectious Diseases of Humans: Dynamics and Control.
      ). Our results indicated that the effective reproduction number for COVID-19 declined to low levels after the first wave and straddled the epidemic threshold of 1.0 in March and April, suggesting that social distancing measures had a significant effect on mitigating the spread of the novel coronavirus. Estimates of early national Rt values for South Korea retrieved from other studies — 2.9 (95% CrI 2.0–3.9) in February (
      • Ryu S.
      • Ali S.T.
      • Jang C.
      • Kim B.
      • Cowling B.J.
      Effect of Nonpharmaceutical interventions on transmission of severe acute respiratory syndrome Coronavirus 2, South Korea, 2020.
      ) and 2.6 (95% CI: 2.3–2.9) in March — are in good agreement with our Rt estimates (
      • Zhuang Z.
      • Zhao S.
      • Lin Q.
      • Cao P.
      • Lou Y.
      • Yang L.
      • et al.
      Preliminary estimates of the reproduction number of the coronavirus disease (COVID-19) outbreak in Republic of Korea and Italy by 5 March 2020.
      ).
      Our results suggest that South Korea has experienced two spatially heterogeneous waves of the novel coronavirus. At the regional level, Seoul and Gyeonggi Province have experienced two waves, whereas Daegu and Gyeongbuk Provinces are yet to experience a second wave. The highest epidemic peak occurred in Daegu and Gyeongbuk Province in late February and early March, with Rt estimated at 4.4 (95% CrI: 2.6–6.6) and 3.5 (95% CrI: 0.9–7.3), respectively. During these epidemic peaks, the doubling times were estimated at 2.8 (95% CI 2.5–4.0) days and 3.6 (95% CI 3.5–4.0) days, respectively, which is similar to a prior estimate of doubling time of 3.8 (95% CI: 3.4–4.2) days (
      • Lee W.
      • Hwang S.S.
      • Song I.
      • Park C.
      • Kim H.
      • Song I.K.
      • et al.
      COVID-19 in South Korea: epidemiological and spatiotemporal patterns of the spread and the role of aggressive diagnostic tests in the early phase.
      ). Gyeonggi Province and Seoul experienced their first wave in late February and early March, respectively. However, sporadic clusters of infections appeared in Seoul and near Gyeonggi Province, immediately after the government eased its strict nationwide social distancing guidelines on May 6. This resurgence of infections in Seoul and Gyeonggi Province (i.e. the province surrounding Seoul), after a sustained period with fewer than five cases per day in each region, led to a second epidemic wave with sub-exponential growth dynamics. In Seoul, the mean doubling time decreased from 7.5 (95% CI: 7.0–8.2) days during the first wave to 6.0 (95% CI: 5.4–6.7) days during the second wave, indicating faster transmission during the case resurgences. Accordingly, our findings revealed sustained local transmission in Seoul and Gyeonggi Province, with the estimated reproduction number estimated to be above 1 until the end of May. In late May, the country implemented 2 weeks of strict social distancing measures, incorporating stringent virus prevention guidelines for the metropolitan area. These measures included the shutting down of public facilities, and regulating bars and karaoke rooms. In the second week of June, South Korea decided to indefinitely extend a period of strict social distancing measures, as nearly all locally transmitted cases were in the metropolitan area.
      Although Korea has had a relatively low number of reported cases compared with other countries such as the USA and China, it is believed that South Korea is currently experiencing yet another resurgence of the virus (
      • WHO
      Coronavirus disease (COVID-2019) situation reports.
      ). Originally, South Korean authorities predicted a resurgence of the virus in the fall or winter; however, this possible second wave started in and around Seoul, which, with 51.6 million inhabitants, accounts for about half of the entire population of the country. Secondary waves of the disease can result from multiple factors, including easing of travel restrictions and the resumption of social activities, especially in the high-population-density areas of Seoul and Gyeonggi Province. Furthermore, a substantial proportion of COVID-19 cases are asymptomatic (
      • Mizumoto K.
      • Kagaya K.
      • Zarebski A.
      • Chowell G.
      Estimating the asymptomatic proportion of coronavirus disease 2019 (COVID-19) cases on board the Diamond Princess cruise ship, Yokohama, Japan, 2020.
      ); thus, they are not detected by surveillance systems, resulting in underestimation of the epidemic growth curve. It was also recently reported that individuals aged 20–39 years in South Korea drove the COVID-19 epidemic throughout society, with multiple rebounds, and an increase in infection among the elderly was significantly associated with an elevated transmission risk among young adults (
      • Yu X.
      • Duan J.
      • Jiang Y.
      • Zhang H.
      Distinctive trajectories of the COVID-19 epidemic by age and gender: a retrospective modeling of the epidemic in South Korea.
      ).
      Our study shows some limitations, including the lack of dates of symptom onset for all cases, relying on a statistical reconstruction of the epidemic curve by dates of symptom onset from a previous study (
      • Shim E.
      • Tariq A.
      • Choi W.
      • Lee Y.
      • Chowell G.
      Transmission potential and severity of COVID-19 in South Korea.
      ). Overall, using the most up-to-date epidemiological data from South Korea, our study highlights the effectiveness of strong control interventions in South Korea, and emphasizes the need to maintain firm social distancing and contact tracing efforts to mitigate the risk of additional waves of the disease.

      Contributions

      ES retrieved, managed, and conceptualized the analysis of the data. ES, GC, and AT analyzed the data. ES and GC wrote the first draft of the paper. All authors contributed to the writing of the paper.

      Financial support

      This work was supported by a National Research Foundation of Korea (NRF) grant to ES, funded by the Korean government (MSIT)[No. 2018R1C1B6001723]. For GC and AT, this work was supported by RAPID NSF No. 2026797.

      Conflicts of interest statement

      None.

      Ethical approval

      Not required

      References

        • Coronavirus
        Coronavirus: South Korea confirms second wave of infections.
        BBC, 2020
        • Anderson R.
        • May R.
        Infectious Diseases of Humans: Dynamics and Control.
        Oxford University Press, New York1991 (See p. 311, eq. (12.23) and discussion)
        • Azmon A.
        • Faes C.
        • Hens N.
        On the estimation of the reproduction number based on misreported epidemic data.
        Stat Med. 2014; 33: 1176-1192
        • Chong K.C.
        • Zee B.C.Y.
        • Wang M.H.
        Approximate Bayesian algorithm to estimate the basic reproduction number in an influenza pandemic using arrival times of imported cases.
        Travel Med Infect Dis. 2018; 23: 80-86
        • Chowell G.
        • Ammon C.E.
        • Hengartner N.W.
        • Hyman J.M.
        Transmission dynamics of the great influenza pandemic of 1918 in Geneva, Switzerland: Assessing the effects of hypothetical interventions.
        J Theor Biol. 2006; 241: 193-204
        • Chowell G.
        • Shim E.
        • Brauer F.
        • Diaz-Duenas P.
        • Hyman J.M.
        • Castillo-Chavez C.
        Modelling the transmission dynamics of acute haemorrhagic conjunctivitis: application to the 2003 outbreak in Mexico.
        Stat Med. 2006; 25: 1840-1857
        • Cori A.
        • Ferguson N.M.
        • Fraser C.
        • Cauchemez S.
        A new framework and software to estimate time-varying reproduction numbers during epidemics.
        Am J Epidemiol. 2013; 178: 1505-1512
        • Covid-19 National Emergency Response Center E
        • Case Management Team KCfDC, Prevention
        Contact transmission of COVID-19 in South Korea: novel investigation techniques for tracing contacts.
        Osong Public Health Res Perspect. 2020; 11: 60-63
        • Fraser C.
        Estimating individual and household reproduction numbers in an emerging epidemic.
        PLoS One. 2007; 2: e758
        • Fraser C.
        Estimating individual and household reproduction numbers in an emerging epidemic.
        PLoS One. 2007; 2: e758
        • He X.
        • Lau E.H.Y.
        • Wu P.
        • Deng X.
        • Wang J.
        • Hao X.
        • et al.
        Temporal dynamics in viral shedding and transmissibility of COVID-19.
        Nat Med. 2020; 26: 672-675
        • KCDC
        The Updates of COVID-19 in Republic of Korea.
        Centers for Disease Control and Prevention, Korea2020
        • KCDC
        The updates on COVID-19 in Korea as of 25 February.
        Korea Centers for Disease Control and Prevention, Seoul, Korea2020
        • Lee W.
        • Hwang S.S.
        • Song I.
        • Park C.
        • Kim H.
        • Song I.K.
        • et al.
        COVID-19 in South Korea: epidemiological and spatiotemporal patterns of the spread and the role of aggressive diagnostic tests in the early phase.
        Int J Epidemiol. 2020; (Online ahead of print)
        • Mizumoto K.
        • Kagaya K.
        • Zarebski A.
        • Chowell G.
        Estimating the asymptomatic proportion of coronavirus disease 2019 (COVID-19) cases on board the Diamond Princess cruise ship, Yokohama, Japan, 2020.
        Euro Surveillance: Bull Maladies Transmissibles = Eur Commun Dis Bull. 2020; 25
        • Muniz-Rodriguez K.
        • Chowell G.
        • Cheung C.H.
        • Jia D.
        • Lai P.Y.
        • Lee Y.
        • et al.
        Doubling time of the COVID-19 epidemic by Chinese province.
        medRxiv. 2020; (preprint)
        • Nishiura H.
        • Linton N.M.
        • Akhmetzhanov A.R.
        Serial interval of novel coronavirus (COVID-19) infections.
        Int J Infect Dis. 2020; 93: 284-286
        • Normile D.
        Coronavirus cases have dropped sharply in South Korea. What’s the secret to its success?.
        Science. 2020; (Magazine article)
        • Park C.K.
        Coronavirus cluster emerges at another South Korean church, as others press ahead with Sunday services.
        South China Morning Post. 2020;
        • Ryall J.
        Coronavirus: Surge in South Korea virus cases linked to church ‘super-spreader’.
        The Telegraph: Telegraph Media Group Limited. 2020;
        • Ryu S.
        • Ali S.T.
        • Jang C.
        • Kim B.
        • Cowling B.J.
        Effect of Nonpharmaceutical interventions on transmission of severe acute respiratory syndrome Coronavirus 2, South Korea, 2020.
        Emerg Infect Dis. 2020; 26: 2406-2410
        • Shim E.
        • Tariq A.
        • Choi W.
        • Lee Y.
        • Chowell G.
        Transmission potential and severity of COVID-19 in South Korea.
        Int J Infect Dis. 2020; 93: 339-344
        • Shim E.
        • Tariq A.
        • Choi W.
        • Lee Y.
        • Gerardo C.
        Transmission potential and severity of COVID-19 in South Korea.
        Int J Infect Dis. 2020; 93: 339-344
        • Tariq A.
        • Roosa K.
        • Mizumoto K.
        • Chowell G.
        Assessing reporting delays and the effective reproduction number: the Ebola epidemic in DRC, May 2018–January 2019.
        Epidemics. 2019; 26: 128-133
        • Wallinga J.
        • Teunis P.
        Different epidemic curves for severe acute respiratory syndrome reveal similar impacts of control measures.
        Am J Epidemiol. 2004; 160: 509-516
        • WHO
        Coronavirus disease (COVID-2019) situation reports.
        2020
        • Yu X.
        • Duan J.
        • Jiang Y.
        • Zhang H.
        Distinctive trajectories of the COVID-19 epidemic by age and gender: a retrospective modeling of the epidemic in South Korea.
        Int J Infect Dis. 2020; 98: 200-205
        • Zhuang Z.
        • Zhao S.
        • Lin Q.
        • Cao P.
        • Lou Y.
        • Yang L.
        • et al.
        Preliminary estimates of the reproduction number of the coronavirus disease (COVID-19) outbreak in Republic of Korea and Italy by 5 March 2020.
        Int J Infect Dis. 2020; 95: 308-310