Mathematical Model of COVID-19 Transmission Dynamics in South Korea: The Impacts of Travel Restrictions, Social Distancing, and Early Detection

: The novel coronavirus disease (COVID-19) poses a severe threat to public health ofﬁcials all around the world. The early COVID-19 outbreak in South Korea displayed signiﬁcant spatial heterogeneity. The number of conﬁrmed cases increased rapidly in the Daegu and Gyeongbuk (epicenter), whereas the spread was much slower in the rest of Korea. A two-patch mathematical model with a mobility matrix is developed to capture this signiﬁcant spatial heterogeneity of COVID-19 outbreaks from 18 February to 24 March 2020. The mobility matrix is taken from the movement data provided by the Korea Transport Institute (KOTI). Some of the essential patch-speciﬁc parameters are estimated through cumulative conﬁrmed cases, including the transmission rates and the basic reproduction numbers (local and global). Our simulations show that travel restrictions between the epicenter and the rest of Korea effectively prevented massive outbreaks in the rest of Korea. Furthermore, we explore the effectiveness of several additional strategies for the mitigation and suppression of Covid-19 spread in Korea, such as implementing social distancing and early diagnostic interventions.


Introduction
The first case of the novel coronavirus, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), was reported by Chinese health officials in Wuhan City in December 2019. Coronaviruses are a family of viruses that can cause illnesses, such as severe acute respiratory syndrome and the Middle East respiratory syndrome [1]. The clinical spectrum of COVID-19 infection is broad, ranging from no symptoms to severe pneumonia. Approximately half of the COVID-19 patients (40-50%) in one study did not show any symptoms [2]. Other patients developed fever, body aches, nausea, or diarrhea [3] typically 2 to 14 days after exposure to the virus. During the initial phase of COVID-19 in China (10-23 January 2020), only 14% of total infections were confirmed. The remaining 86% were not identified or quarantined, contributing to a community spread in China. Later, on 23 January, the Chinese government implemented a total lockdown of Wuhan City, which prevented further community spread [4]. Owing to this strict intervention, the number of new COVID-19 cases in China dropped to the single-digits in early March. Other countries, however, started to report an increasing number of confirmed cases. COVID-19 spread from Wuhan to the rest of China and neighboring Asian countries such as Thailand, Japan, Korea, Singapore, and Hong Kong. As of 31 March, the number of confirmed cases worldwide had exceeded 820,000, and the number of deaths had reached more than 40,000 in more than 110 countries. The World Health Organization declared COVID-19 a pandemic on 11 March 2020 ( [5,6]).
In South Korea, the first confirmed case of COVID-19 was reported on 20 January, and the number of cases steadily increased to 30 by 17 February. Beyond the initial sporadic outbreaks from overseas inflows, including China, there was a large-scale community spread around specific regions. Beginning with the 31st confirmed case on 18 February, however, an unprecedented rapid increase of confirmed cases was recorded in the area of Daegu and North Gyeongsang Province (Gyeongbuk). It is still unclear how or from whom the 31st case became infected. Daegu is Korea's 4th largest city and is surrounded by Gyeongbuk (see Figure 1). In the process of tracing the contacts of the 31st case on 18 February, it was revealed that she attended worship services at the Sincheonji Daegu Church on 9 February and 16 February, raising the possibility of community spread through the church service. Indeed, an inspection of more than 1000 Sincheonji Daegu members on 20 February resulted in the identification of more than 111 cases. Epidemiological surveillance of these 111 newly identified cases revealed that most cases presented to Cheongdo Daenam Hospital in Gyeongbuk. Health authorities later imposed cohort isolation on the hospital after a large cluster of COVID-19 cases was reported. Although most of the confirmed cases were concentrated in Daegu or Gyeongbuk, COVID-19 spread to nearly every part of the country [7]. Across the nation, 79.7% of the cases have been proven to be linked to cluster transmission, and 62.5% have been linked to Shincheonji [8]. As of 31 March 2020, 9786 patients were confirmed in Korea, 7984 of whom (81%) were confirmed in Daegu and Gyeongbuk [8]. The spread of COVID-19 in Korea displayed a spatial heterogeneity, a clustered rapid spread in Daegu and Gyeongbuk, and a steady slow spread in the rest of the country.
Mathematical modeling can provide useful insights for mitigating emerging infectious diseases such as recent SARS-CoV, MERS-CoV, COVID-19 outbreaks. Particularly, the mathematical modeling of the spread of diseases becomes more important as the movement of people and other living organisms increases due to globalization [9,10]. Multi-patch models have been developed to understand the global spread of transmission dynamics in various emerging and reemerging infectious diseases [11][12][13][14]. Since the COVID-19 outbreak began, many studies have been carried out to investigate the basic or control reproductive numbers, whereas others predicted the outbreak peak time and evolution of the epidemic to guide intervention policies ( [15][16][17][18]). The national and global spread of COVID-19 was explored, accounting for the effect of the metropolitan-wide quarantine of Wuhan and surrounding cities [19]. A stochastic transmission model was combined with COVID-19 data from Wuhan City to estimate the temporal evolution of transmission [20]. The study also predicted the probability of how newly introduced cases would generate outbreaks in other areas. Time-varying basic reproduction number has been investigated in Italy and forecasted the COVID-19 spread in the region [21]. The critical role of a latency period is highlighted under three different modeling of latency periods in SEIR-based models [22]. Some studies investigated the effectiveness of COVID-19 mitigation strategies, such as travel restrictions ( [23,24]) and social distancing ( [25][26][27][28]). These nonpharmaceutical measures reduced transmissibility by a maximum of 34% without resorting to a strict lockdown strategy [29].
In this study, we develop a two-patch model of COVID-19 transmission dynamics, representing a hot spot (patch 1: Daegu and Gyeongbuk) and a slow-spreading area (patch 2: the rest of Korea). In each patch, we employed a mathematical model with five epidemiological compartments (susceptible, exposed, infected, quarantined, and recovered), which is based on the previous work [14]. The two patches are connected by a mobility matrix that captures the number of people traveling in South Korea (provided by the Korea Transport Institute (KOTI) in 2016) [30]. It includes origin and destination data on various modes of transportation, including airplanes, cars, buses, ships, and trains. The incubation and recovery rates were assumed to be the same in the two areas and other parameters such as transmission rate were assumed to be different in the two patches and were estimated using cumulative incidence data from 18 February to 24 March 2020. Using the estimated model parameters, we carried out a patch-specific sensitivity analysis. Furthermore, we investigated the impacts of various intervention strategies on the patch-specific transmission dynamics, including limiting traveling between the two regions, implementing social distancing, and early diagnosis. This paper is organized as follows. In Section 2, we present a data description followed by a single-patch SEIQR model and a two-patch SEIQR model with a mobility matrix and discuss the basic reproduction numbers (the local reproduction number for a single-patch model and the global reproduction number for a two-patch model). We have carried out parameter estimation and sensitivity analysis in Section 3. Numerical simulations under various mitigation scenarios have been investigated in Section 4. The paper concludes with a discussion of results and conclusions in Sections 5 and 6. Following Section 6, we include the Appendix A containing the mathematical work, showing the derivation of the basic reproduction number for a two-patch SEIQR model.

Characteristics of the Early COVID-19 Outbreaks in South Korea
Daily confirmed cases in South Korea were publicly available from the Korea Centers for Disease Control and Prevention (KCDC) [31,32]. Due to the 31st case, a rapid increase occurred, especially in the Daegu/Gyeongbuk area, and Figure 2 displays the daily confirmed cases in Daegu/Gyeongbuk and the rest of Korea. The spatial characteristics of cumulative incidence are displayed in Figure 1. Note that there were hot spots in Daegu and Gyeongbuk, which were around 80% of the total cases, implying a high level of spatial heterogeneity in the number of confirmed cases of Korea.
The timeline of administrative measures in Daegu and Gyeongbuk is displayed in Figure 2 and Table 1. The Korean government and Dague and Gyeongbuk implemented administrative countermeasures in response to COVID-19 outbreaks (as shown in Table 1) [31,32]. On 17 February, the first drive-through testing facility was installed, where symptom checks and sample collection could be performed in less than 10 min. The Korean government declared Daegu and Gyeongbuk as "special care zones" on 21 February. The Korean government upgraded the infectious disease crisis alert to the severe level on 23 February and requested citizens of Daegu to refrain from voluntarily going out and restricting their movement for at least 2 weeks. All schools were closed, and employees at many companies were working at home [7]. Prompt detection followed by isolation was advised to mitigate the spread of the infectious disease. Moreover, many medical professionals and administrative officials from the private sector volunteered to serve in Daegu. The citizens of Daegu begin to implement the "328 Dague movement", including the minimization of movement, as well as social distancing.

A Single-Patch Model: SEIQR Model
We first considered a single-patch SEIQR model. The birth and natural death of the population were not considered. The population was divided into five groups according to disease state. First, at time t, S(t) is the number of susceptible individuals and E(t) is the number of exposed individuals. Next, I(t) denotes the number of infectious individuals, including asymptotic infection but not undetected (or not confirmed) individuals. Finally, Q(t) represents the number of confirmed and therefore, quarantined individuals, while R(t) represents the number of recovered individuals. From the schematic diagram depicted in Figure 3, we have the following system of five ordinary differential equations: where Here, β represents the transmission rate, and κ is the rate of exposed people (E(t)) to become infectious (I(t)), p represents a proportion from E(t) to I(t) (asymptomatic and undetected but infectious individuals), α stands for the rate of quarantined, and γ represents a recovery rate (either from I(t) or Q(t)). The definitions and values of all model parameters are listed in Table 2. The dynamic behavior of the model can be classified by the basic reproductive number. This threshold condition determines whether an infectious disease will spread in a susceptible population when the disease is introduced into the population [33]. The threshold is calculated by using the next-generation method as where ρ is defined as the spectral radius of the next-generation matrix FV −1 . Here, F is the rate of appearance of new infections, and V is the transfer of individuals out. Then, we find the Jacobian matrix of F (x) and V (x), and denote F = [∂F i /∂x j ] and V = [∂V i /∂x j ](i, j = 1, 2, · · · , 5) evaluated at the disease free equilibrium point E DFE = (N, 0, 0, 0, 0). More details are given in the Appendix A. We can explicitly calculate the basic reproduction number, R 0 (or the local reproduction number for a single-patch model) as We concluded that E DFE is locally asymptotically stable if, and only if, R 0 < 1. When R 0 > 1, the matrix has a positive real eigenvalue, and this means that the disease-free equilibrium is unstable.

Two-Patch Transmission Model
In this subsection, we extended the single-patch model to the two-patch transmission model by incorporating the mobility matrix between the two patches. Other state variables and parameters are the same as the single-patch model given in the previous subsection. The mobility matrix is defined by M = (m ij ) 2×2 , as proposed in [13,14]. The mobility matrix captures the movement of individuals between patch 1 and patch 2 and we set m ij to be the number of people traveling from a patch i to a patch j per unit time (day) and m ji to be the number of people traveling from a patch j to a patch i per unit time (day). In this study, we assumed that the number of traveling between two patches is symmetric (m ij = m ji ). Moreover, we assume that the quarantined individuals do not move to the other patch, the number of susceptible people entering patch i is m ji S j /(N j − Q j ), and the number of susceptible people leaving patch i is m ij S i /(N i − Q i ). It is also assumed that the traveling occurs instantaneously, so no infection is transmitted during the traveling process. As shown in the schematic diagram of the two-patch model, Figure 4, the following system of differential equations describes the resulting two-patch transmission model: where . It is assumed that p is the same in the both patches p = p 1 = p 2 ) and that κ (or γ) is the same in the both patches (κ = κ 1 = κ 2 and γ = γ 1 = γ 2 ). The descriptions and values of all model parameters are listed in Table 2.   For the two-patch model, the (global) basic reproduction number, R 0 , is computed by the next generation matrix approach [33] as where ρ is defined as the spectral radius of the next-generation matrix FV −1 . Here, F is the rate of appearance of new infections, and V is the transfer of individuals out. Then, we find the Jacobian matrix of F (x) and V (x), and denote F = [∂F i /∂x j ] and V = [∂V i /∂x j ] (i, j = 1, 2, · · · , 4) evaluated at the disease free equilibrium point, E DFE = (N 1 , 0, 0, 0, 0, N 2 , 0, 0, 0, 0). Since the basic reproduction number is obtained by evaluating at the disease equilibrium under the assumption of the initial stage (i.e., S i (t) ≡ N i (t)), thus, we can replace N i (t) − Q i (t) with N i (t) (Q i (t) ≡ 0 at the the disease equilibrium point). The dynamics of infected compartments for a two-patch model can thus be described by where m 12 = m 12 /N 1 and m 21 = m 21 /N 2 . Indeed, we obtain the (global) basic reproduction number, R 0 of the two-patch model as where D 1 = κ 1 + m 12 , D 2 = γ + α 1 + m 12 , D 3 = κ 2 + m 21 , D 4 = γ + α 2 + m 21 , and The details of the above derivations are given in Appendix A.

Estimated Parameter
As shown in Figure 2, the daily confirmed cases of COVID-19 in the Daegu and Gyeongbuk area show special features; the number of confirmed cases demonstrated a slow initial increase at the beginning, shortly after the 31st confirmed case, a rapid increase of cases. After the peak of confirmed cases, the number of cumulative incidences remained steady (almost no increase). To capture this special feature, the entire time interval ([0, 40] days) was divided into the following 3 subintervals: a small increase interval [t 0 , t 1 ), a rapid increase interval [t 1 , t 2 ), and a slower increase interval [t 2 , t end ). Therefore, we estimated three different transmission rates, a small β 1,1 during [t 0 , t 1 ), a rapid increase with a transmission rate β 1,2 during time [t 1 , t 2 ), and a slower increase of incidence after peak with a transmission rate β 1,3 during [t 1 , t end ) ( β 1,3 < β 1,1 < β 1,2 ).
We set the initial susceptible population to be 5.12 million and 46.52 million in the Daegu and Gyeongbuk and the rest of Korea, respectively (ie, S 1 (t 0 ) = 5.12 × 10 6 & S 2 (t 0 ) = 46.52 × 10 6 ) as the population size of the two patches. The initial quarantine numbers of Q 1,2 (t 0 ) were taken from the cumulative incidence data on 18 February 2020. We assumed the number of exposed E 1,2 (t 0 ) and recovered R 1,1 (t 0 ) to be zero. The number of initially infected individuals (I 1,2 (t 0 )) was not measurable and thus was estimated using the method below.
Using the report by Workman [36], we set the percentage of asymptomatic individuals to be 33.3% (i.e., p = 0.333). To estimate a typical number of traveling in Korea, we used data from the 2016 Korea Transport Data Base, a national organization that produces, manages, and provides various transportation data for Korea (KOTI) [30]. The remaining 10 parameters and 2 initial numbers of infected (I 1,2 (t 0 )) were estimated. An inverse least-squares formulation was used to find a parameter set (θ) that minimized the least-squares error between cumulative incidence of the model outputs (C 1,2 ) and actual cumulative incidence data (D 1,2 ) from 18 February to 24 March 2020.
The least-squares error J is defined as where θ represents a set of parameters, t j time points (with j = 1, · · · , N), and D 1,2 actual cumulative incidence data (a subscript denotes a patch; C 1 , D 1 for patch 1 and C 2 , D 2 for patch 2). The optimal parameter set is given byθ = arg min θ J(θ).
We used the Nelder-Mead simplex method [38], which aims to minimize the least-squares error through a direct search using MATLAB's built-in optimization function fminsearch. The optimal parameter values for the two-patch model are given in Table 2, and the fitted curves generated by the estimated parameters are shown in Figure 5. Cumulative incidence of the two-patch model (solid blue curve) and COVID-19 data (red circles) in each patch are in good agreement.

Sensitivity Analysis
To assess the variation of the model outputs in response to small perturbations in the parameter values, we performed a global sensitivity analysis that described how a model output changed in response to varying the parameters over a subspace of the parameter space. A global sensitivity method developed by [39] is a variance decomposition method. First, input parameters are varied, which causes model output to be varied. Next, the variance of the model output is quantified. Then, a partial variance is calculated by partitioning of model output variance to determine what fraction of the model output variance is caused by the variation of each model parameter input.
The model outputs we considered were C 1 , cumulative incidences in the Daegu and Gyeongbuk, and C 2 , cumulative incidences in the rest of Korea except for Daegu and Gyeongbuk. The sensitivity indices computed by the method described in [39] are shown in Figure 6. The left panel for patch 1 shows various sensitivity analysis; the model output C 1 was highly sensitive to the parameters α 1 , t 1 , β 1,2 , p 1 , t 2 , where α 1 indicates a detection rate from I 1 to Q 1 , and β 1,2 represents a transmission rate during (t 1 , t 2 ), while p 1 is a proportion from E 1 to I 1 .
Next, the right panel for patch 2 shows various sensitivity analysis; the model output C 2 was highly sensitive to the parameters β 2 , t 1 , β 1,2 , p 2 , α 1 , α 2 , where β 2 indicates a transmission rate of patch 2, t 1 represents the time when a rapid increase of incidences was observed in the Daegu/Gyeongbuk area, β 1,2 is a transmission rate during (t 1 , t 2 ) in the Daegu/Gyeongbuk area, p 2 indicates a proportion from E 2 to I 2 , and α 1,2 represents a detection rate from I 1,2 to Q 1,2 , respectively. Cumulative incidence in patch 1, C 1 is sensitive to m 21 (mobility from patch 2 to patch 1) while cumulative incidence in patch 2, C 2 is sensitive to m 12 (mobility from patch 1 to patch 2).

The Impact of Mobility on the COVID-19 Transmission Dynamics
First, we investigate the impacts of movement (travel or mobility) on the two-patch transmission dynamics: Daegu and Gyeongbuk (the epicenter) are represented by patch 1 and the rest of Korea is patch 2. Other parameters are fixed as the baseline parameter values given in Table 2. The average number of travelers between the two areas was obtained from the Korea Transport DataBase in 2016 [30] (m 12 and m 21 are given in Table 2). In 2016, the average number of people traveling was 1.9 million per day, however, the average daily nationwide level of traffic for the first three months of 2020 was 9.7% lower than the same period in 2019. After the alert level of COVID-19 raised to the highest level on 22 February 2020, traffic in the 4th week of February was down by 26.1% compared with the one of 2019 [40].
We assume that a reduction (0.5-30%) of the average number of mobility to mimic the travel reduction due to the COVID-19 outbreak. Now, we present the impact of mobility on the global R 0 given in (5)  . The increase was more dramatic when the transmission rate in the epicenter was higher (right). A previous study also reported similar simulation results [14], arguing that this is a result of the dilution effect (people move more than COVID-19 spread from the epicenter to other areas). Next, we present the impacts of mobility on daily and cumulative incidence in the two patches as migration rate changes from 1% to 10% (see Figure 8). The top panel shows the baseline results under m 12 = m 21 = 1%. The peak size of daily incidences and cumulative confirmed cases in patch 1 was highest (519 cases and 7821 cases, respectively). The number of confirmed cases in the rest of Korea remained relatively small (peak, 67 cases; cumulative incidences, 3497 cases) as the actual COVID outbreak with the highest spatial heterogeneity.
The middle panel shows the results under m 12 = 10%, m 21 = 1%, which implies that more people were traveling out of the Daegu and Gyeongbuk to the rest of Korea. We observed a decrease in the confirmed cases (approximately a 50% reduction) in the Daegu and Gyeongbuk, whilst a sharp increase in cases was expected (∼400%) in the rest of Korea ( Figure 8; the first panel vs. the second panel).
Interestingly, increases in traveling between the two areas (epicenter and the rest of Korea) resulted in a reduced number of confirmed cases in the epicenter, as shown in the bottom panel under m 12 = 10%, m 21 = 10%. It turns out that travel restrictions from Daegu and Gyeongbuk to other areas prevented the spread of COVID-19 across the whole country.

The Impact of Social Distancing and Early Detection on the COVID-19 Transmission Dynamics
In this subsection, we explore the impacts of two other interventions (social distancing and early detection) on two-patch COVID-19 transmission dynamics. Assuming minimal levels of travel between the two areas (epicenter and the rest of Korea), we evaluated the impacts of social distancing and early detection on the patch-specific peak size and cumulative incidence. Social distancing was simply modeled by reducing the transmission rate, σβ i , i = 1, 2 with σ = 0.5, where β i was an estimated value and given in Table 2. Figure 9 shows the effects of σβ i on patch-specific incidence of COVID-19 under two different combinations of σβ i , i = 1, 2. The top panel shows the results under the baseline scenarios as given in Table 2, and the bottom panel shows the results assuming 50% reduced transmission rates. Halved transmission rates (50% reduction) led to an approximate 67% reduction in the spread in the two areas, respectively. This is consistent with the results of the global sensitivity analysis (patch-specific cumulative incidence was the most sensitive to the transmission rates in patch 2, as shown in Figure 6). Next, we investigate the impact of early detection on two-patch COVID-19 transmission dynamics. What if asymptomatic patients were efficiently traced, identified, and quarantined? To address this question, we simulated community spread increasing the estimated model parameters α 1,2 . The early detection was modeled by incorporating a re-scaling constant, η in the rate, ηα i , i = 1, 2 with η = [1,500], where the estimated α 1 was 0.031 in the epicenter while the value was 0.35 in the rest of Korea. Figure 10 shows the effects of ηα i on the patch-specific incidence of COVID-19 under two different combinations of ηα i , i = 1, 2. The top panel shows the results under η = 1 and the baseline scenarios given in Table 2. The bottom panel shows the results when the early detection rate was increased. Five hundred-fold increases of α 1,2 resulted in an approximate 70% reduced spread in each area. Taken together, the identification of asymptomatic individuals by extensive testing contributed to controlling massive outbreaks. This is consistent with the results of the global sensitivity analysis (patch-specific cumulative incidence was the most sensitive to the early detection rates in patch 1, as shown in Figure 6).

Discussion
In this work, we developed a single-patch and a two-patch model for COVID-19 transmission dynamics, focusing on the largest outbreak in Daegu and Gyeongbuk and the rest of South Korea. We investigate the impacts of three different intervention strategies on the COVID-19 transmission dynamics. These include travel restrictions, social distancing, and the efficient early detection of confirmed cases followed by quarantine.
First, travel restrictions are particularly useful in the early stage of the COVID-19 outbreak confined to a certain area. However, travel restrictions may be less effective once the outbreak is more widespread [23]. Our simulations showed that travel restrictions delayed the spread of the disease into the rest of Korea from the epicenter and reduced the number of total incidences (see the leftmost panel of Figure 11, 12,000 cases at 1% vs 16,000 cases at 9%). An increased travel movement between Daegu and Gyeongbuk and the rest of Korea resulted in a reduced number of confirmed cases in the epicenter but not in the rest of Korea. The travel ban in Wuhan and the international travel restrictions did not show noticeable differences in the epidemic in Wuhan but delayed ∼3 days for other locations in mainland China [24]. Another study demonstrated that internal mobility restrictions are effective only if prohibitively high (50% of travel) restrictions are applied [26]. In another infection study, Ref. [41] suggested that short-term mobility between heterogeneous patches does not always contribute to overall increases in the prevalence of tuberculosis. In [14], the strategy of the migration restriction may result in locally bigger outbreaks rather than extinguishing the pandemic spread, even though the spread is delayed. Our results are consistent with other studies showing that travel restrictions may not be particularly effective to reduce the cumulative incidence.
Second, social distancing measures aim to reduce the frequency of contact, and increase the physical distance between individuals, reducing the risks of person-to-person transmission. A large number of published studies examining social distancing conducted during the 2009 influenza pandemic concluded that social distancing controlled the epidemic [26]. Data-driven models of COVID-19 spread in China [27] have shown that strict social distancing for 30 days had significant benefits in limiting community spread Wuhan and Hubei. Choi et al suggested that social distancing is crucial in suppressing the spread of the COVID-19 in Korea [28]. Furthermore, the transmission of viruses was lower with the physical distancing of 1 meter or more, compared with a distance of less than 1 m [25]. Our simulations also showed that strict social distancing led to a significant reduction in total incidences. Approximately 2000 fewer cases were expected as transmission rates were decreased by 10% (see the middle panel of Figure 11). From 12 April to 18 April, after four weeks of strict social distancing was enforced, the number of confirmed cases decreased dramatically [8]. In the situation of a lack of therapeutics and without vaccines, social distancing is one of the basic protective measures against COVID-19.
Lastly, one of the notable things about COVID-19 is the high rate of asymptomatic but infectious individuals. In Korea, approximately 33.3% of the confirmed cases were asymptomatic [36]. In Italy, 42.5% of the confirmed infections detected were asymptomatic [42]. Thus, identification of asymptomatic virus carriers by aggressive testing and contact tracing could potentially reduce the size of the epidemic. Our simulations varying the detection rate α suggested that approximately 2000 confirmed cases (approximately 20% of the total incidence) could be reduced by doubling detection rates (see the rightmost panel of Figure 11). Another study also demonstrated that the identification of undocumented COVID-19 cases suppressed the spread of COVID-19 [4]. Early detection is also a necessary method for preventing transmission. To mitigate the spread of the disease, it is important to detect the outbreak of the patient early, conduct quarantine and treatment, and identify contacts through effective epidemiological surveillance systems and investigations.

Conclusions
We have investigated the COVID-19 transmission dynamics using a two-patch SEIQR system with a mobility matrix. The two patches represent the epicenter (Daegu/Gyeongbuk) and the rest of Korea to capture two distinct characteristics of COVID-19 outbreaks with the highest spatial heterogeneity. Some of the essential epidemic parameters including patch-specific transmission rates and the local and global R 0 have been estimated using empirical cumulative incidence data of COVID-19 from 18 February to 24 March 2020. We have investigated the impacts of several interventions implemented by the Koran government including travel restrictions, social distancing, and early diagnostic interventions. Our results indicate that the two-patch dynamics are strongly dependent on the following three key factors: mobility, transmission rates, early diagnostic interventions. Overall, implementing all possible interventions together would be the most effective in controlling the COVID-19 outbreaks, however, that would increase the higher economic burden and social costs as well. The present study focused only on the first two months of COVID-19 outbreaks, and we would like to extend our work to a multi-city model, which can capture the entire course of COVID-19 transmission in South Korea for our future study.

Conflicts of Interest:
The authors declare no conflict of interest.