Differential impacts of contact tracing and lockdowns on outbreak size in COVID-19 model applied to China ========================================================================================================= * Cameron J. Browne * Hayriye Gulbudak * Joshua C. Macdonald ## Abstract Rapid growth of the COVID-19 epidemic in China induced extensive efforts of contact tracing and social-distancing/lockdowns, which quickly contained the outbreak and has been replicated to varying degrees around the world. We construct a novel infectious disease model incorporating these distinct quarantine measures (contact tracing and self-quarantine) as reactionary interventions dependent on current infection levels. Derivation of the final outbreak size leads to a simple inverse proportionality relationship with self-quarantine rate, revealing a fundamental principle of exponentially increasing cumulative cases when delaying mass quarantine or lockdown measures beyond a critical time period. In contrast, contact tracing results in a proportional reduction in reproduction number, flattening the epidemic curve but only having sizable impact on final size when a large proportion of contacts are “perfectly” traced. We fit the mathematical model to data from China on reported cases and quarantined contacts, finding that lockdowns had an overwhelming influence on outbreak size and duration, whereas contact tracing played a role in reducing peak number of infected. Sensitivity analysis and simulations under different re-opening scenarios illustrate the differential effects that responsive contact tracing and lockdowns can have on current and second wave outbreaks. ## 1 Introduction The current COVID-19 pandemic began in Wuhan, China, where infections grew rapidly and spread throughout the country in late December 2019 and January 2020. In order to contain the virus, drastic measures, such as travel restrictions alongside extensive lockdowns and contact tracing efforts, were implemented. The overall success of these control strategies in suppressing the outbreak in China has been recognized in several studies (*1, 2*). An important question is which intervention had the largest impact, or in more detail, quantifying the effect of each intervention on case reduction. The problem is relevant not only for retrospective analysis, as all countries including China face the task of controlling ongoing or possible second wave outbreaks of COVID-19. The strategies currently available for the fight against COVID-19 are often classified as non-pharamacuetical interventions (NPIs), since consensus vaccines or treatments have not been found to date. The effectiveness and aims of NPIs may vary by country and type of intervention. While the goal of large-scale lockdowns and social distancing is often characterized as “flattening the curve”, whereas successful contact tracing may suppress outbreaks, a more nuanced picture of their potential impact on epidemic trajectories is necessary. A few studies have quantified impact of travel restrictions (*3, 4*) and lockdowns inducing large-scale changes in contact patterns or depletion of susceptible individuals (*5, 6*), showing the efficacy of these interventions in China. Yet, the precise qualitative and quantitative effect of brute force interventions such as lockdowns (or widespread social distancing), versus the more targeted strategy of contact tracing, on the outbreak shape is less explored. Traditionally the influence of control strategies on outbreaks has been theoretically investigated in compartmental ordinary differential equation models of the susceptible-infected-recovered (SIR) type. Analysis yields the herd immunity (or critical vaccination) threshold for suppressing an outbreak by proportionally reducing the reproduction number, ℛ, below one, along with a nonlinear relationship between ℛ and final outbreak size when ℛ is above one. Furthermore, inference of parameters by fitting the model to data can help to determine the effect of interventions. However both the analytical and parameter estimation approaches are challenged by the dynamic nature of control strategies as public health authorities and individuals react to an evolving outbreak. While the early phase of COVID-19 can be characterized by exponential growth, case saturation occurred much earlier than would be predicted by the basic SIR model due to the comprehensive control measures that have been deployed. In particular, stringent lockdown with broad (self- and contact tracing) quarantine interventions reduced the pool of susceptible individuals, effective contact rate and secondary transmissions. Several models have utilized time-dependent transmission or isolation rates to capture the dynamics (*3, 7*), and recent work has also considered removal of susceptible individuals at a constant rate (*6*). Here we develop a generalized SIR-type model incorporating a total (government mandated and individual) social distancing rate, along with contact tracing, both *depending on overall infection rate*, in order to fit an observed reactionary public health system and derive novel formulae for outbreak size. In order to quantify the impacts of contact tracing and comprehensive social distancing (self-quarantine or lockdowns), we simultaneously utilize case and quarantined contact data from China to estimate parameters in our model. Furthermore, through computational and theoretical analysis of the model, we can explore the sensitivity of distinct epidemic measures (e.g. outbreak size, peak number of infected, timing and extent of self-quarantine) to interpretable control parameters. These investigations allow us to dissect how combinations of NPIs, such as contact tracing and lockdowns, may influence sequential outbreaks through loosening and tightening of control measures. The emergent picture is of distinct qualitative impacts of contact tracing and lockdowns on the outbreak, variable in scope and timing, and dependent on underlying disease parameters. A better understanding of these differential effects can help shape or suppress the epidemic curve of COVID-19 in a sustainable and acceptable manner to societies. ## 2 Model with Social Distancing and Contact Tracing We formulate a SEIR model (Fig.1 and generalized equations are given in SI), which modifies a detailed differential equation system of contact tracing during outbreaks (*8*). The model variables include: susceptible (*S*), exposed (*E*) and infectious (*I*) individuals; social-distanced (or self-quarantined) susceptible (*S**q*), exposed (*E**q*) and infectious (*I**q*) individuals; contact traced susceptible (*S**c*), exposed (*E**c*), and infectious (*I**c*) individuals; and the decoupled compartments of (safely) isolated reported cases (*R*) with a subset of currently quarantined contact-traced cases (*R**c*). The full system of equations, along with a table of variables and parameters, are given in the *SI*. Here we highlight a few key features of the model. Parameters *β, β**c*, and *β**q* represent transmission rates of reported un-quarantined, contact traced quarantined and social-distanced infected individuals, respectively, where *β**c* and *β**q* reflect reductions in transmission due to contact tracing and social distancing which are generally imperfect (e.g. tracing individual after they become infectious, looseness in following stay-at-home orders). A critical control parameter is the total rate of susceptible transition to (contact traced or self-) quarantine state, *ψλS*, with ![Graphic][1], depending on force of infection *λ* = *βI* + *β**q**I**q* + *β**c**I**c*, the proportion of contacts traced *ϕ*, the probability of transmission upon contact *p* and the self-quarantine (social distancing or lockdown) factor *σ*. The dependence on force of infection reflects mechanism of contact tracing (*8*), along with the responsive nature of broader social distancing/quarantine measures to current transmission. Other important parameters include *α**q*, the rate of return to susceptible from social distancing, and *ν**q* the susceptibility of social-distanced individuals measuring the looseness of the social distancing measures. A simplified version of the model assuming *perfect indefinite* quarantine and *perfect* contact tracing is given by the following system: ![Formula][2] where *R* is decoupled, and the additional decoupled compartments of self-quarantined and contact traced (fit to data) are detailed in the SI. In Section 4, we fit both the simplified model and full model above to total case and quarantined contact data of China. ## 3 Reproduction Number and Outbreak Size The *basic reproduction number, ℛ*, of the “perfect quarantine” model (1) is derived as: ![Formula][3] where *S*(0) is the initial susceptible population size. We also consider (time-dependent) *effective reproduction number ℛ**e*, governed by *S*(*t*), calculated continuously during the outbreak. Note that *R**e* for the full model (Fig. 1) depends upon the quarantined populations *S**q*(*t*) and *S**c*(*t*), along with their imperfect susceptibility and transmission parameters (*ν**q*, *ν**c* and *β**q*, *β**c*) and is formulated in *SI*. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/06/15/2020.06.10.20127860/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2020/06/15/2020.06.10.20127860/F1) Figure 1: Full compartmental model of COVID-19 with reactive contact tracing and self-quarantine. We extend the basic SEIR model describing infection of susceptibles (*S*) at rate *β*, pre-infectious compartment (*E*) of average length *τ* and infectious compartment (*I*) of avg. period *T*. Contact tracing is incorporated by tracking a proportion *ϕ* of reported case contacts whom are either susceptible (*S**c*) with reduction in susceptibility *ν**c* for average time of 1*/α**c*, or infected (*E**c*, *I**c*) and have reduced reproduction number *β**c**T**c*. Self-quarantine (social distancing) of susceptibles (*S**q*) is implemented at rate *σ* proportional to “force of infection” *λ* = *βI* + *β**q* *I**q* + *β**c**I**c* with reduction in susceptibility *ν**q* for average time of 1*/α**q*, and can be infected (*E**q*, *I**q*) with reduced reproduction number *β**q* *T**q*. In the simplified system (1), we consider “indefinite” self-quarantine (*α**q* = 0) where contact traced susceptibles always transition (*θ**c* = 1), along with perfect contact tracing and self-quarantine (*ν**c* = *ν**q* = *β**c* = *β**q* = 0). Next we present novel theoretical results of relations between ℛ and final outbreak size in terms of contact tracing and quarantine/social distancing. In order to obtain the results on final outbreak size, we assume an indefinite quarantine period (*α**q* = 0), as in (1). In this way, the outbreak size can represent magnitude of first, second or subsequent waves dependent on contact tracing and self-quarantine parameters. Define the final proportion of susceptible individuals ![Graphic][4] and the *final (cumulative) epidemic size* ![Graphic][5]. Under simplifying conditions, we derive the final size of *U**∞* and *𝒞**∞* dependent on the (assumed) uniform susceptibility of self-quarantined and contact traced individuals, *ν* (detailed in *SI*). In the best case scenario where self-quarantine or contact tracing perfectly prevents susceptible infection (*ν* = 0), we obtain the exact formulae: ![Formula][6] where *ψ* depends on self-quarantine factor (*σ*) and ℛ depends on contact tracing proportion (*ϕ*). Note that each formula can account for arbitrary initial conditions in order to quantify how ℛ*e* and quarantine measures affect outbreak size beginning at any stage (see *SI* for formal derivations). In the above formulae (3), note the classical relation between final susceptible proportion *U**∞* and ℛ. Then, since self-quarantine factor (*σ*) has almost all weight in the total quarantine parameter (*ψ* = *σ* + *ϕ*(1 *− p*)*/p*), *𝒞**∞* has a simple inverse proportionality relationship with *σ*. For instance, given that *σ* is relative to force of infection, to reduce the outbreak size by 90% (compared to no lockdown case), the authorities should implement strict quarantine at approximately 9 times the rate of infection. The rates can be translated to time periods for more interpretability since the doubling time of cumulative incidence depends directly on the force of infection. For this example, in each cumulative incidence doubling period, “perfectly quarantined” individuals should increase 10-fold. Authorities should do much better than keeping pace with new infections and strive for very large values of *σ*, which as we will see from fitting results, was instrumental for China rapidly curbing their epidemic. ## 4 Data Fitting & Efficacy of Quarantine Measures in China We utilize data on total reported cases and quarantined contacts in mainland China published in daily reports by NHC (*9*). Although there are certain issues with the reported case data (*10*), qualitative results on how contact tracing and social distancing/lockdown measures affected outbreak size were robust when fitting raw or smoothed data (see Supplementary). We utilize data on total reported cases and quarantined contacts in mainland China published in daily reports by NHC (*9*). Despite inconsistencies in the reported case data (*10*), qualitative results on how contact tracing and social distancing/lockdown measures affected outbreak size were robust when fitting raw or smoothed data (see *SI*). We fit both our full model ((1) in *SI* and Fig.1) and simplest model (1) simultaneously to (cumulative) reported case data and (daily number of) quarantined contacts. Overall the models can fit the two datasets well with several parameter sets. We constrained the estimations by fixing the incubation period (time to infectiousness, *τ* = 3 days (*11*)), and infectious period (time to isolation, *T* = 4.64 days (*12*)). Note for the full model, we constrain the reduction in transmission due to contact tracing based on a large study of cases and their contacts in Shenzhen, China (*12*). With the above assumptions, fixing the baseline reproduction number without any control (ℛ0,*b* = *βT*) produced similar results compared to when *β* was a fitted parameter. The value of ℛ0,*b* chosen was 6, in line with other studies (*6,7*). Overall the extra parameters in the full model or adding unreported cases only slightly reduced error from the simplified model fit, however the additional detail in the full model allow us to vary more features of the social distancing and contact tracing interventions (see *SI* and Fig. 3(b)). Although the proportion of contact traced cases varied when fitting the models under different assumptions of ℛ0,*b* or incorporating a proportion of unreported cases into model, a consistent pattern emerged on the impact of the contact tracing probability *ϕ* on the outbreak. Despite impacting ℛ, larger estimates of *ϕ* tend to correlate with larger baseline ℛ0,*b* values, which diminishes any effect on outbreak size. Based on these observations, we utilize the parameter fitting from the simplified model (1) with ℛ0,*b* = 6, displayed in Fig. 2. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/06/15/2020.06.10.20127860/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2020/06/15/2020.06.10.20127860/F2) Figure 2: Modelrecapitulatescaseandcontact tracing data, estimates large rate of self-quarantine (via lockdown) rapidly contains outbreak. (a) Cumulative reported cases and daily number of quarantined contacts (inserted figure) simultaneously fit to model (1). (b) The corresponding daily reported cases in model with data and inferred subset of contact traced infected. (c) Average and daily impact of contact tracing on ℛ*e* from estimated parameters utilizing statistical method based on generation time distribution. Observe the reduction in ℛ*e* withorwithoutcontacttracingduetosocialdistancing/lockdowns is the major factor to rapidly contain the outbreak. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/06/15/2020.06.10.20127860/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2020/06/15/2020.06.10.20127860/F3) Figure 3: Sensitivity analysis of contact tracing and self-quarantine. All together the figures illustrate how contact tracing (*ϕ*) has less impact on outbreak size (*𝒞**∞*) than self-quarantine (*σ*) due to underlying concavity of their nonlinear relationships derived in final size formula (3). Fitted parameter values for China show that rapidly implemented lockdowns contained the outbreak, where a delay of 2 weeks in 95% of susceptible population being self-quarantined would result in a 10-fold increase in cumulative cases. Although contact tracing had little effect (2%) on final outbreak size, it “flattened the curve”, reducing peak infected by 34%. (a) Contact tracing (CT) proportion *ϕ* versus outbreak size *𝒞**∞* (nonlinear relationship) and reproduction number ℛ (linear relationship) for 3 levels of self-quarantine (SQ) factor *σ* differing by order of magnitude of 10. (b) *𝒞**∞* as a (nonlinear) function of *σ* for 3 levels of SQ susceptibility reduction *ν**q*. (c) Total (final) self-quarantined individuals, along with time until 50%, 75%, 95% of population is self-quarantined as SQ factor *σ* varies. (d) Epidemic curve trajectories for 3 levels of CT proportion *ϕ* marked in (a). (e,f) Contour maps of outbreak size and peak infected (with time to peak inserted) depending on CT proportion *ϕ* and SQ rate *σ*. The best fit value of ℛ was found to be 3.74 (CI 3.37-3.80), where the estimated value of the proportion of traced contacts/reported case (*ϕ*) is 0.38 (CI 0.35-0.46). The self-quarantine factor (relative to force of infection) was consistently estimated to be high (*σ* = 1240 in displayed fit). The initial amount of infected individuals, *I*, was estimated at 778 on January 21, 2020, when the dataset begins. Note that we fit the model starting at this date, close to the time when major lockdowns began (e.g. a *cordon sanitaire* implemented in Wuhan on Jan. 23). However our force-of-infection (*λ*) dependent rate formulation of mass self-quarantine (*σλ*) actually allows for a very similar epidemic trajectory when initiating the model a month earlier with one infected individual (*I* = 1) and all other parameters the same as our fit starting from Jan. 21 (Fig.S1). Complete parameter values and uncertainty analysis are presented in *SI*. An alternative direct calculation for (daily) ℛ*e* utilizes the daily case data and estimates of the serial interval (generation time) distribution (*12, 13*). Here, by incorporating quarantined contact data, we also infer the efficacy of contact tracing (*8*). Although missing information, in particular the amount of infected quarantined contacts, hinders our ability to directly evaluate contact tracing impact on a daily basis, we utilize the predicted relative transmission and incidence of contact traced individuals from our model fit to assess ℛ*e* with and without contact tracing. Despite noisiness inherent in the daily case data and reproduction number, the fitted compartmental model ℛ*e* captures the general trend, and indicates the strict population-wide lockdowns were the main quarantine measure (as opposed to contact tracing) which rapidly contained the outbreak in China (Figs. 3(b),3(c)). To determine the effect of the main control parameters on the final outbreak size, we perform sensitivity analysis (Fig. 3). While the inferred contact tracing level for China was not found to significantly reduce the final outbreak size, our results suggest a larger effect on reducing peak infection levels by flattening the curve. In particular, by varying contact tracing proportion *ϕ*, we observe the total reduction of 34% in peak infected size as compared to the 2% impact on cumulative outbreak size (Fig. 3(d)). In general, we observe that the time to peak increases with *ϕ*, reflecting the curve flattening, however this time period eventually decreases for sufficiently large values of *ϕ* as contact tracing effectively suppresses the outbreak (Fig. S7). In addition, with sufficiently large contact tracing coverage, outbreak size can be significantly reduced when there is less stringent lockdown (less total quarantined and more time to enact quarantine). Yet even in this case, some level of broader social distancing measures is almost certainly needed in combination with contact tracing. The self-quarantine (lockdown) factor (*σ*) has a large impact on outbreak size even though it does not affect ℛ. As predicted by our derived inverse proportionality relationship (3), there are escalating costs as *σ* decreases, i.e. as self-quarantine action lessens relative to ongoing infection rate (Fig 3(b)). For example, if the estimated time for 50% of initial susceptible population of China to be self-quarantined (*∼* 2 weeks from Jan. 21) had been delayed by just one week, then the total number of cases would be approximately 10 times larger (Fig. 3(c)). In the full model, the additional parameters *ν**q*, *β**q* (measuring looseness of the lockdown) also impacts the outbreak size (Fig. 3(b)). Going from *ν**q*, *β**q**/β ≈* 0 (as predicted for China) to *ν**q* = *β**q**/β* = 0.05 to *ν**q* = *β**q**/β* = 0.15, the outbreak size would increase by a factor of 1.14 and 1.45, respectively. In addition, the estimated rate of return from self-quarantine (*α**q*) for China was estimated to be very small, emphasizing the strictness of the lockdown. ## 5 Quarantine Interventions for COVID-19 2nd Wave A major question is how public health authorities should guide loosening of broad lockdown measures after initial containment of COVID-19, while optimally responding to any subsequent outbreaks induced by the relaxations. Here we analyze how the scale and rate of different reactive contact-based interventions affect 2nd wave outbreaks under two different scenarios of loosening, namely *Instantaneous Return of Several Sectors* (IRSS) or via *Gradual Return of Self-Quarantined* (GRSQ). The goal is to attain qualitative insights on the timing and allocation of control strategies for shrinking, flattening or delaying the subsequent outbreak curves. By varying self-quarantine factor *σ*, contact tracing probability *ϕ* and looseness of social distancing *ν* under the distinct relaxation policies in our model parameterized to data from China, we observe potential consequences of different strategies. Simulating the instantaneous return of 80% of self-quarantined individuals (IRSS strategy), with no change in parameters (and crucially the same “reactive lockdown” factor *σ*) we observe that the cumulative number of infected cases for the 2nd wave (outbreak size) and peak infected was 75% and 58%, respectively, of the 1st wave. Furthermore, a similar number of individuals as during the first wave lockdown re-enter self-quarantine about 6 weeks after relaxation (see Fig.4(a)). When the contact tracing efforts are enhanced after lockdown (to *ϕ* = .65), outbreak size and peak infected are 54% and 16%, respectively, of the 1st wave, and the curve is flattened, i.e. the peak outbreak size shrunk and the time to peak outbreak size increased. Finally if contact tracing is doubled to *ϕ* = 0.75, the 2nd wave outbreak size and peak infected are 25% and 3%, respectively, of the 1st wave. In addition, the number of individuals re-entering self-quarantine was reduced, revealing that contact tracing can be an effective tool for managing the epidemic with a less stringent lockdown. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/06/15/2020.06.10.20127860/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2020/06/15/2020.06.10.20127860/F4) Figure 4: Combined impact of contact tracing and responsive social distancing (self-quarantine) measures under distinct re-opening strategies. Model simulations predict 2nd wave outbreak shape after instantaneous re-normalization of several sectors (IRSS) by releasing 80% of self-quarantined (a,b,c) versus gradual return of self-quarantined (GRSQ) to normalcy (d,e,f). Increased contact tracing levels (*ϕ*) can flatten or suppress (when highly effective) subsequent outbreak in (a,d). Social distancing measures responsive to new incidence shrink the outbreak size dependent on rate (*σ*) and looseness (*ν**q*) of self-quarantine in (b,e) and (c,f), respectively. In addition to the differential effects on daily infected counts of contact tracing and responsive self-quarantine (flattening versus crushing), observe the much longer delay in 2nd wave under the GRSQ re-opening policy, which can buy time for effective treatments or vaccines. In the case of GRSQ strategy, after containing initial outbreak with lockdown, we increase the return to “normalcy” rate to *α* = 0.01, where half the social-distanced return to normalcy in the approximate half-life time given by *t*1*/*2 = ln 2*/α* = 72 *days*. Assuming other parameters remain constant (including the reactive SQ factor *σ*) the second peak, emerging with a 100 day delay, reduced to 42% of the first wave, however the number of infected individuals settle into a rather large quasi-equilibrium resulting in more cumulative cases (see Fig.4(d)). Here there is a balance of force of infection induced self-quarantine (*σλ*) and reversion of individuals to their normal contact behavior (*α*), leading to an insufficient amount of population social distancing for reducing cases below a certain level. On the other hand, after loosening the lockdown, when the contact tracing efforts are enhanced or doubled, the peak size significantly diminished (27% or 0.3% of 1st wave), along with the number of self-quarantined. Importantly, for about 6 months (or the whole year in the case of doubling *ϕ*), the number of infected cases stayed significantly low. This suggests that gradual release of self-quarantined individuals with increasing contact tracing efforts can be utilized as a strategy to gain time until vaccination, while reinstating societal interactions in a carefully measured stepwise fashion. Responsive re-implementation of lockdown (or social distancing) measures is crucial for reducing any second wave outbreak. Reduction in SQ factor *σ* by 1*/*2 (or 1*/*4), as predicted simply by the inverse proportionality in the derived final size formula (3), results in twice (or four times) more cumulative cases for the 2nd wave, and the simulations show the same relations between peak size (see Fig.4(b)). Although the number of self-quarantined individuals eventually become the same with the different SQ rates, the delay in implementing large-scale self-quarantine (in response to incidence) makes significant differences in the final (and peak) outbreak size. For the simulations presented in Fig. 4(b), a delay of just 9 days from the baseline parameter case results in twice as many infections, and a delay of 18 days induces four times the infected individuals. Compared to instantaneous release, the increased quarantine exit rate (*α*) under gradual return resulted in larger (but delayed) peak and total outbreak size inversely proportional to declines in SQ factor *σ* (see Fig.4(e)). Finally, varying the looseness of the quarantine (measured by uniform susceptibility and infectivity values *ν**q* = *β**q**/β*) from perfect quarantine to 25% (or to 50%) looseness, leads to approximately 1.3 times (or to 2 times) more total and peak infections during the outbreak. Different from the rate of SQ, the proportionality relations are nonlinear, thus a slight looseness in quarantine can still offer an effective intervention, but the cases will increase at a growing rate as the measures become less strict (see Fig.4(c),4(f)). ## 6 Discussion In this study, we compare how two distinct types of contact-based control strategies, contact tracing and large-scale lockdowns/self-quarantine or social distancing, impact the characteristics of single or sequential COVID-19 outbreaks. We find that contact tracing generally is less effective in decreasing outbreak size for rapidly spreading pathogens (high baseline reproduction number ℛ0,*b*), unless the tracing is very efficient. On the other hand, widespread lockdowns/social distancing interventions can lower outbreak size inversely proportional to an increase in the rate of self-quarantine. Our analysis indicates that China benefited from the heavy influence of lockdowns by rapidly containing the quickly growing COVID-19 cases, and, despite massive efforts, contact tracing was less influential in bringing down the epidemic. Despite the difference in the targeted nature of contact tracing versus the more indiscriminate lockdown measures, we contend there is a similar reactive quality to both control strategies. Contact tracing reacts to reported cases by tracking and (to varying degrees) quarantining individuals whom have been contacted. Mass social distancing or self-quarantine reflects a natural response by both governments and individuals which intensifies as cases build, a phenomenon that has been labeled as “exponential whiplash” (*14*). These features motivate us to construct a COVID-19 model with both contact tracing (mechanistically) and self-quarantine (phenomenologically) dependent on force of infection. In contrast to another model which assumes a linear rate of self-quarantine (*6*), the nonlinear social distancing rate captures a contagion-like behavioral response to infected cases, and allows us to derive novel formulae for final outbreak size. Furthermore the model provides a good simultaneous fit to both cumulative reported cases and daily quarantined contact data from China. An important distinction between contact tracing and lockdowns is their mode of action, namely preventing onward secondary infections by early tracking of likely infected cases in the former and large-scale depletion (or shielding) of susceptible individuals for the latter. This contrast determines how they affect the major epidemiological quantities of reproduction number and outbreak size in our “transmission-reactive” formulation. In particular, contact tracing proportionally reduces ℛ, akin to vaccination, leading to a nonlinear relationship with final outbreak size, which decreases substantially only as ℛ approaches one. The responsive self-quarantine factor does not affect ℛ, and we derive a simple inverse proportionality with outbreak size. This can be translated to a time of action for quarantine measures, analytically demonstrating the escalating impacts of delaying implementation of responsive lockdowns beyond a critical time period, which has been observed in other studies via simulation (*15, 16*). Even though similar levels of self-quarantine would eventually be reached in our model as incidence grows, the cost of delays can result in a large excess of cases. Although we find that the extensive lockdowns and social distancing was a much larger factor in controlling COVID-19 *outbreak size* in China, our sensitivity analysis shows that contact tracing did dampen and delay *peak number of infected* despite its more limited impact on the cumulative count. In this way, contact tracing can flatten the incidence curve, easing the strain on limited hospital resources. A combination of expediently enacted contact-based interventions may be the best strategy, where effective contact tracing and responsive social distancing measures can synergistically and efficiently suppress an outbreak. However COVID-19 has proved to be a particular challenge and large-scale lockdowns have been a needed antidote for controlling outbreaks in several countries. The drastic self-quarantine orders can also reduce case numbers to a more manageable level and hopefully allow for effective contact tracing in the event of incidence occurrence after easing restrictions. The capacity to respond to the continuing threat of COVID-19 will be vital for minimization of sequential epidemic waves. We investigated control measures under an instantaneous normalization of contact for a large portion (or several sectors) of the population versus a more gradual release of self-quarantined individuals back into social interactions. Our results show that increased contact tracing efforts can alter the second outbreak shape, either reducing and spreading out the number of infected or completely suppressing cases for highly efficient tracing. Social distancing or lockdown measures responsive to incidence can effectively compress the second peaks, with the timing being critical again. Either measure will depend upon sufficient case detection and reporting, highlighting the importance of testing. Furthermore, indefinite or reoccurring strict lockdowns are likely to impart too high of an economic cost, and our model shows that looser restrictions and contact tracing can still reduce a second wave to manageable levels. Additionally, the strategy of gradual release of quarantined sectors can substantially delay the second wave, possibly buying time for effective treatments or vaccines to be developed. There are factors not considered in our current study which may have played an important role in determining the COVID-19 epidemic in China. For example, we neglect a more fine-grained regional structure within China for our model and data fitting. However, a major novelty of this work was to incorporate data on the quarantined contacts, which was compiled solely for the whole of China. Obtaining provincial quarantine records may allow for simultaneous fitting of the heterogeneous spread of the virus in different regions of China. Furthermore, more detailed contact tracing data quantifying the proportion of reported cases whom were traced can allow for superior accuracy in estimating efficacy of contact tracing, which can add confidence to our conclusion that lockdowns had substantially larger influence in controlling the COVID-19 outbreak data. Nevertheless, the analytical and qualitative results here illustrate the differential effects that reactive contact tracing and lockdowns or mass self-quarantine have on outbreak shape. This knowledge and further investigation may offer insights for the public health response to COVID-19 outbreaks. ## Data Availability Data publicly available or upon request. ## Supporting Information (SI) SI Appendix Figs. S1 to S7 Tables S1 to S3 ## Acknowledgments CJB, HG, and JCM are supported by a U.S. National Science Foundation RAPID grant (DMS-2028728). HG was also supported by a grant from the Simons Foundation/SFARI(638193). CJB is partially supported by an NSF grant (DMS-1815095). * Received June 10, 2020. * Revision received June 15, 2020. * Accepted June 15, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## References 1. 1.WHO, Report of the WHO-China Joint Mission on Coronavirus Disease 2019 (COVID-19), [https://www.who.int/docs/default-source/coronaviruse/who-china-joint-mission-on-covid-19-final-report.pdf](https://www.who.int/docs/default-source/coronaviruse/who-china-joint-mission-on-covid-19-final-report.pdf) (2020). 2. 2. M. U. Kraemer, et al., Science 368, 493 (2020). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjgvNjQ5MC80OTMiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wNi8xNS8yMDIwLjA2LjEwLjIwMTI3ODYwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 3. 3. S. Lai, et al., Nature pp. 1–7 (2020). 4. 4. H. Tian, et al., Science 368, 638 (2020). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjgvNjQ5MS82MzgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wNi8xNS8yMDIwLjA2LjEwLjIwMTI3ODYwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 5. 5. J. Zhang, et al., Science (2020). 6. 6. B. F. Maier, D. Brockmann, Science 368, 742 (2020). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjgvNjQ5Mi83NDIiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wNi8xNS8yMDIwLjA2LjEwLjIwMTI3ODYwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 7. 7. B. Tang, et al., *Available at SSRN 3537099* (2020). 8. 8. C. Browne, H. Gulbudak, G. Webb, Journal of theoretical biology 384, 33 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jtbi.2015.08.004&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26297316&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F15%2F2020.06.10.20127860.atom) 9. 9.NHC, National Health Commission of the People’s Republic of China, [http://en.nhc.gov.cn](http://en.nhc.gov.cn) (2020). 10. 10. T. K. Tsang, et al., The Lancet Public Health (2020). 11. 11. X. He, et al., Nature medicine 26, 672 (2020). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F15%2F2020.06.10.20127860.atom) 12. 12. Q. Bi, et al., The Lancet Infectious Diseases (2020). 13. 13. L. Ferretti, et al., Science 368 (2020). 14. 14.The Economist (2020). 15. 15. J. Dehning, et al., Science (2020). 16. 16. S. Pei, S. Kandula, J. Shaman, medRxiv (2020). [1]: /embed/inline-graphic-1.gif [2]: /embed/graphic-1.gif [3]: /embed/graphic-2.gif [4]: /embed/inline-graphic-2.gif [5]: /embed/inline-graphic-3.gif [6]: /embed/graphic-4.gif