## ABSTRACT

Since its emergence, the phenomenon of SARS-CoV-2 transmission by seemingly healthy individuals has become a major challenge in the effort to achieve control of the pandemic. Identifying the modes of transmission that drive this phenomenon is a perquisite in devising effective control measures, but to date it is still under debate. To address this problem, we have formulated a detailed mathematical model of discrete human actions (such as coughs, sneezes, and touching) and the decay of the virus in the environment. To take into account both discrete and continuous events we have extended the common modelling approach and employed a hybrid stochastic mathematical framework. This allows us to calculate higher order statistics which are crucial for the reconstruction of the observed distributions. We focused on transmission within a household, the venue with the highest risk of infection and validated the model results against the observed secondary attack rate and the serial interval distribution. Detailed analysis of the model results identified the dominant driver of pre-symptomatic transmission as the contact route via hand-face transfer and showed that wearing masks and avoiding physical contact are an effective prevention strategy. These results provide a sound scientific basis to the present recommendations of the WHO and the CDC.

## Introduction

The phenomenon of SARS-CoV-2 transmission by pre-symptomatic, otherwise seemingly healthy, individuals poses a major challenge for policy makers’ efforts to achieve control of the COVID-19 pandemic, as traditional health strategies rely on case detection through manifestation of symptoms^{1}. However, the mechanism that enables this transmission is not fully understood. Generally, respiratory viruses such as SARS-CoV-2 propagate via four modes of transmission^{2}: direct physical contact between people, indirect physical contact via intermediate object, droplets and droplet nuclei. Transmission by droplets and droplet nuclei is mediated by virus containing particles that were emitted when a person coughs, sneezes or speaks. The droplets travel less than 1.5m ^{3}, due to their size, and settle on the facial membranes of nearby individuals or on surfaces. Droplet nuclei remain suspended in the air and may infect a susceptible individual once they penetrate the respiratory tract. The commonly accepted cutoff between droplets and droplet nuclei is 5*µm* ^{2}. However, Xie et al.^{3} showed that droplets that are smaller than approximately 100*µm* evaporate to their nuclei size before reaching the ground.

The relative contribution of the different modes of transmission in indoor environments is still under debate^{4–9}. The controversy revolves about the relative importance of the droplet nuclei mode of transmission. Several studies have argued that the transmission of the SARS-CoV-2 virus is mediated primarily by close and unprotected contact (e.g., via physical contact and droplets)^{4–7}, while others have argued that breathing droplet nuclei is the main mode of transmission^{8, 9}. The close contact transmission hypothesis relies on the analysis of COVID-19 cases^{6} and the relatively low secondary attack rate (SAR, the probability of an infected person to infect a susceptible person) of 10%-16% that was observed in households^{5, 10–13}. The droplet nuclei hypothesis relies on several theoretical investigations^{8, 9}. The attempts to identify SARS-CoV-2 in air sampling taken from infection isolation rooms in hospitals and households yield conflicting results. Several studies^{14–16} found positive samples while others^{17–19} reported negative air samples (for example, Dohla et al.,^{19} reported negative air samples that were taken from households).

The aim of this study is to quantify the relative contribution of the different modes of transmission of SARS-CoV-2 to infection by pre-symptomatic individuals. We focus in this study on the transmission within a household environment, the venue with the highest risk of infection^{5, 10}. The approach taken here is an integrative detailed mechanistic modelling that describes explicitly the transfer of SARS-CoV-2 between individuals in different modes of transmission, similar to the approach used by Nicas and Sun^{20} and by Atkinson and Wein^{21} for quantifying the modes of transmission of respiratory viruses. In this work we have extended the mathematical framework of Atkinson and Wein^{21} to take into account random discrete human actions (such as coughs, sneezes and contact with objects and other people), rather than considering only the mean kinetics. This was achieved by employing a hybrid stochastic mathematical framework which allows us to explicitly calculate higher order statistics which are crucial for the reconstruction of the observed distributions. Following this, the model is validated by reconstruction of observed fundamental attributes of the pandemic, the secondary attack rate (SAR) and the serial interval distribution. Then, the model is used to assess the contribution of each of the transmission modes as well the effectiveness of different prevention measures.

### Outline of the Mathematical model

The model presented in this study incorporates the processes that influence the number of SARS-CoV-2 virus particles transferred from a pre-symptomatic infected individual, henceforth the primary, to a susceptible individual, henceforth the secondary, and the probability to become infected (Figure 1).

Contact transmission begins when the primary touches his facial membranes and, as a result, contaminates his own hands. Then, the primary transfers the virus either through direct physical contact (Figure 1, mode 1) or indirectly via small frequently touched object (fomits), such as a doorknob or a faucet, (Figure 1, mode 2) to the hands of the secondary. Eventually, the secondary places his hands into nose, mouth or eyes, which might cause an infection^{21, 22}. The droplet and droplet nuclei modes of transmission (Figure 1, modes 3 and 4, respectively), begin when the primary coughs, sneezes, or speaks and expels virus containing droplets. Droplets larger than 100*µm* settle by gravity within 1.5*m* ^{3} and contaminate large surfaces such as furniture and table tops (environmental surfaces), while smaller droplets dry out and form droplet nuclei which remain suspended in the air. As a result, the droplet nuclei may be carried over distances greater than 1.5 m by the air currents of the room^{3}. The deposition of droplets directly on the mocusa of close contacts is a rare event in workplace or household settings^{21}. Therefore, we have considered here only the contamination of environmental surfaces by the droplets after they have settled. The contaminated areas on the environmental surfaces might also contaminate the hands of the secondary individual when he touches them. The probability of infection increases with the number of SARS-CoV-2 particles that reach facial membranes of the secondary individual.

The processes described in Figure 1 consist of both discrete random events which are the actions of the individuals, such as touching each other or touching the facial membranes and continuous events which consist of environmental processes. Hence, we used a hybrid continuous and stochastic-jump framework^{23} to describe the dynamics of the transmission and infection processes using a coupled system of differential equations. The actions of the individuals are described as stochastic jump Poisson processes, while the environmental processes are described using continuum dynamics. The model explicitly tracks the dynamics of the concentration on the hands of the individuals (equation 8), the concentration on the fomites (equation 9), concentration on environmental surfaces (equation 11) and the concentration of the droplet nuclei in the air (equation 10). A complete list of the model equations and values of the corresponding parameters are provided in the Methods section. Since the actions of the individuals are represented as a stochastic process, we have conducted a Monte Carlo simulation in which multiple realizations were computed to obtain the appropriate ensemble statistics. Using a Monte-Carlo simulation the probability distributions are embedded in the input parameters directly and allow the comparison of the model results to observed distributions. Each realization begins when the primary become infected and begins an incubation period whose duration is drawn from a log-normal distribution with a mean of 5 days and standard deviation (SD) of 0.45 days^{24}. The viral load of the primary increases exponentially with time^{21} reaching a maximal level at the end of the incubation period^{25}. During that time, the primary and secondary individuals perform a series of randomized actions such as touching fomites, touching environmental surfaces, coughing, sneezing, talking, touching each other, or each touching his own face. The probability that the secondary individual will be infected is determined from his accumulated exposure over a time interval (equation 12) using the dose-response curve that was reported for SARS-CoV-1^{26} (equation 7) and assumed to be similar to SARS-CoV-2. Each realization ends when the primary develops symptoms, in accordance with the public health policy that isolates the primary at the onset of symptoms.

We define a reference simulation as a simulation which corresponds to a normal, pre-symptomatic behaviour (parameters in Table 1).

## Results

### Validation of the model

A necessary validation criteria for a model such as the one described in this study is to correctly simulate the distribution of the serial interval and the SAR. The serial interval is the time period between the symptoms’ onset in the primary and the secondary. Its distribution is closely associated with the estimation of the reproductive number and key transmission variables in epidemic models and is important for optimization of the length of the obligatory quarantine period and contact tracing strategies^{29, 30}. The serial interval distribution of COVID-19 was estimated in many countries and was usually found to be gamma distributed with mean between 4.03 to 6.3 days and standard deviation between 3 and 4.2 days (Figure 2A, shaded area)^{10, 31–33}.

The model prediction for the distribution of the serial interval and the SAR is obtained by conducting a Monte-Carlo simulation by solving equations 8 to 12 (see the Methods section) for 10000 realizations of the reference simulation (needed for convergence). Figure 2A compares the model predicted serial interval distribution (red line) with the distributions reported in the literature. As seen, the model prediction was well between the bounds of the different estimates of this distribution (Figure 2A).

As an additional validation, we compared the model prediction of the SAR to the values reported in the literature. The model predicts a SAR of 11.5% in the reference simulation, which is within the reported values ranging between 10 − 16%^{5, 10–13}. We have also analyzed the contagious period of pre-symptomatic patients by examining the cumulative SAR over time (Figure 2B). As seen, the contagious period begins approximately 30 hours before the symptoms’ onset, with increasing probability of infection as the onset of the symptoms approaches. This result is consistent with the estimation of He, et al.^{25}, that inferred from data of 77 transmission pairs (i.e., primary and secondary) a contagious period of approximately 2 days before symptoms’ onset.

Some of the parameters’ values were obtained from studies that also reported the range of variability of these values. Therefore, we have performed an extensive sensitivity analysis to check the robustness of the results (see Supplementary Information). The model’s results remain within the range of the values reported in the literature for the examined range of parameters.

### Modes of transmission in pre-syptomatic cases

Analyzing the realizations of the reference simulation, we have quantified the contribution of the different modes of transmission to the overall exposure in scenarios where the secondary was infected (Figure 3). Out of the total viral dose that was transmitted to the secondary, 64.5% (Inter quartile range, IQR: 55% − 80%) was received during direct contact events (mode 1) and 26% (IQR 13% − 32%) was received during indirect contact via fomites events (mode 2). The contribution of the large droplet route (mode 3) was negligible while the droplet nuclei transmission (mode 4), contributed 9.5% (IQR 3.6% − 12%) of of the total viral dose. Hence, according to our results, the contact mode of transmission (either direct or indirect) is the dominant mode of infection, accounting, overall, to the transfer of 90% of the viral dosage from the primary to the secondary. The main process that underlies the contact mode of transmission is the hand-face transfer. Therefore, hygienic and behavioral measures that operate on the elements that constitute the contact processes are expected to significantly reduce the risk of infection. These will be analyzed in the following section.

### Reducing the risk of infection

The fact that contact transmission is the main route of pre-symptomatic transmission, suggests that the hygienic and behavioral measures (HBMs) advised to the public should focus on reducing the contamination on the hands or somehow interrupt the hand-face transfer. We have examined five HBMs: Washing hands, cleaning fomites, avoiding physical contact (i.e., maintaining social distancing), wearing a mask and gloves. Naturally, conservative precautions measures would be an implementation of all these measures simultaneously. However, strict adherence to all these HBMs would be difficult to endure and to maintain over a long period of time. Therefore, we have tried to sort out several combinations of HBMs that should be readily implemented by the public, while significantly lowering the risk of infection. As the SAR is a proportion, it is appropriate to compare the HBMs in terms of odds ratio (*OR*), i.e., the odds that the secondary would be infected when a given combination of HBMs is taken, compared to the reference scenario in which no HBM is applied. Generally, any HBM that results in *OR* less than 1 decreases the risk of infection (i.e., provide smaller SAR than the reference)^{34}. However, in practice the lower the *OR*, the more effective the HBM combination is at lowering the risk. The values brought here are in terms of *OR* alongside with 95% confidence interval (95% CI).

Washing hands is known to remove (and also destroy) virus particles from the hands and it is the simplest measure to implement. Our simulations show that washing hands once every hour rather than 3 times a day, as in the reference simulation (Table 1), results in *OR* of 0.72 (95% CI 0.67-0.8) (Figure 4A, column H). This result is consistent with intervention studies that have shown that increased hand washing decreased respiratory illness by 20%, albeit different viruses were studied^{22}. This phenomenon seems counter intuitive, as we found that 90% of the viral dosage is transmitted through the hands and it was expected that washing it would remove the contamination. In order to understand the reason for the relatively limited effect of hand hygiene, we have examined the dynamics of the virus concentration on the hands of the secondary individual (Figure 4B). This concentration exhibits a periodic behaviour, with a period of approximately 30 to 40 minutes, that is governed by contact events on fomites and the face. Therefore, hand washing is expected to dramatically reduce the risk for infection if it occurs at at frequency higher than 40 min. Unfortunately, such frequent hand washing is unrealistic.

Cleaning the fomites more frequently reduces the virus repositories that are available. Cleaning of the fomites 10 times a day rather than twice a day, as in the reference simulation, results in *OR* of 0.84 (95% CI0.77 −0.92), similar to washing hands more frequently (Figure 4A, column F). A combined strategy that consists of frequent hand washing and cleaning fomites does not decrease the risk considerably and results in *OR* of 0.70(95% CI 0.63 − 0.76).

Wearing a surgical mask or a respirator may reduce the hand-face transfer of virus particles^{35} as well as the inhalation exposure to viral particles. Although it is difficult to asses the reduction of the transfer coefficient from hand to facial membranes due to the use of a mask or a respirator, there are measurements regarding the protection provided against airborne transmission of bacteria and viruses. Available experimental results on N95 filtering face-piece respirators and surgical masks reported a protection factor of 2 − 10 for aerosols^{36}. Hence, we have used a reduction of factor 2 in the hand-face transfer as well as in the exposure to airborne virus particles as a conservative estimate. With this value of the parameters, the simulated *OR* was 0.49 (95% CI 0.44 − 0.54) (Figure 4A, column M). Wearing gloves reduces the concentration of virus on the hands since the transfer efficiency from and to surfaces is halved with latex gloves^{37}. As a result the use of gloves resulted *OR* of 0.43 (95% CI 0.39 − 0.48) (Figure 4A, column G).

Avoiding physical contact interrupts the main route of transfer between the two individuals and leads to an *OR* of 0.37 (95% CI 0.33 − 0.41) (Figure 4A, column S). Compared to the *OR*s resulted by taking a single HBM, this is the most effective step. A combined strategy that includes wearing masks and avoiding contact results in *OR* of 0.14 (95% CI 0.12-0.16, Figure 4A, column M+S), which is lower than wearing a mask and gloves (*OR* of 0.18, 95% CI0.16 − 0.21, Figure 4A, column M+G) or frequent hand washing, cleaning fomites and avoiding physical contact (*OR* of 0.23, 95% CI 0.21-0.27, Figure 4A, column H+F+S). The above shows that wearing masks and avoiding physical contact is the most effective HBM. However, wearing masks for long periods of time is difficult. Nevertheless, cleaning fomites, washing hands, and avoiding physical contact also provides considerable reduction in the *OR*. Hence, implementing these HBMs meticulously may save people the discomfort and limitation, that are associated with having to wear a mask constantly.

## Discussion

We have analyzed the possible routes of pre-symptomatic transmission in household scenarios. Using a validated model, we were able to identify the main mode of transmission as contact associated, mostly direct contact, but also contact mediated by fomites. The principal element in this transmission is hand-face transfer. Frequent hand washing and fomite cleaning coupled with the avoidance of physical contact result in a protection similar as wearing gloves and a mask. Although the present work does not account for highly populated indoor environments (such as work or commercial spaces), the relative importance of different processes is expected to remain very similar. However, in scenarios where people are in close proximity to each other (mass transit or during medical procedures), the contribution of the droplets to the transmission is expected to be larger. Hence, such scenarios require further investigation.

The relative contribution of the airborne transmission route is currently under debate. Several theoretical investigations concluded that the contribution of the airborne transmission route is significant^{8, 9}. These studies assume that the viral load in the expelled droplets due to cough, sneeze and talking is 10^{8} − 10^{9} viruses/ml which is considerably higher than the reported median value^{38} of 7.5 · 10^{5} viruses/ml that was used in this work. In order to examine the effect of this assumption we used our model with the same maximal viral load of 10^{8} − 10^{9} viruses/ml. This resulted in a significant overestimation of the SAR. Specifically, we obtained a SAR of 70%, and to a SAR of 36% when only the airborne route was considered.

Our analysis, as with all modeling exercises, has several limitations and requires certain assumptions. The model does not account for contact patterns that prevail in households with young children and does not take into account the diurnal cycle of activity. The model parameters, such as the dose response curve, the viral shedding coefficients and transfer coefficients were chosen on the basis of knowledge of the SARS, other strains of coronavirus, or other bacteria^{26, 39}. Although the model is stable with regard to variations in these parameters, more information on the key characteristic of the disease may reduce some uncertainties.

In conclusion, our findings can provide an important tool for decision makers while advising the public of the HBMs that are necessary to impede the progression of the epidemic. As it seems, recurrent outbreaks are expected to occur, as many countries will have to establish a fine balance between posing restrictions on society and allowing citizens to lead their life as normally as possible^{40}. Under such a reality, the model presented in this study can be used to quantify the contribution of different HBM measures in order to devise guidelines that mitigate the risk of infection in scenarios of workplaces or schools, sports and cultural events, and mass transportation.

## Methods

In this section we delineate the mathematical details of the indoor transmission model that was developed in this study. The model simulates the transmission of the SARS-CoV-2 virus between two individuals that share the same indoor space, say, a room in a household or an office, using an agent based modeling approach. Specifically, it explicitly tracks the health condition of each individual, his actions, and the contamination in the indoors environment over time. The actions of the agents are randomized and therefore multiple realizations are required in order to obtain the appropriate ensemble statistics. At the beginning of the each realization, one individual, the primary, is infected but is pre-symptomatic. The other individual, the secondary, is susceptible.

The dynamics of the model are driven by the following processes:

Inoculation. Individual hands, which could be contaminated with viruses, touch the mouth and other facial membranes and exchange viruses with them.

We denote as the number of virus particles that pass from the facial membranes to the fingers of individual

*I*(*I*is either primary or susceptible) after a single contact at time*t*, where*ε*_{skin}_{−}is the fraction of viruses that transfer in a skin to skin contact,_{skin}*L*(*t*) is the current viral load of individual I at time*t*and*V*is the volume that carried from the facial membranes by the touch._{inoculation}The number of viruses that pass from the fingers to the facial membranes is denoted as

*and equals to where is the concentration of the finger and**S*is the effective surface area of the finger._{finger}Wearing gloves reduces the skin-skin transfer coefficient by a factor of and wearing masks reduces the exchange of viruses between the facial membranes and the fingers by a factor of

Coughing, sneezing, and talking. Individuals that cough, sneeze and talk emit particles within the range of 1

*µm*to 2000*µm*,^{22, 41, 42}. Infectious individuals emit particles that are loaded with viruses and thereby contaminate the nearby surfaces and the air. Particles greater than 100*µm*^{3}are referred to as large droplets. These droplets travel up to approximately 1.5m^{3}before they settle on environmental surfaces. For simplicity, we assume that the contaminated area created by the large droplets from a single cough, sneeze, and talk is a semicircle with radius of 1.5m. The deposition of droplets directly on the mocusa of close contacts is a rare event in workplace or household settings^{21}. Therefore, we have considered here only the contamination of environmental surfaces by the droplets after they have settled.Particles that are smaller than 100

*µm*evaporate before they reach the ground and remain suspended in the air as droplet nuclei^{3}. The droplet nuclei then disperse and might eventually infect a susceptible individual when he inhales them.Approximately 99% of the emitted particle volume following a sneeze, cough or talking is in droplets whose diameters are larger than 100

*µm*^{22, 41}.We assume that wearing masks completely blocks the emission of droplets and droplet nuclei.

Physical contact. Individuals also exchange viruses via physical contact with contaminated body parts. For simplicity, we consider the transfer from and to the fingers of the individuals, since the fingers (and specifically fingertips) are expected to be the most contaminated, as they are in contact with fomites and the facial mocusa.

The number of viruses that pass from the primary to the secondary is the product of the skin to skin transfer fraction

*ε*_{skin}_{−}, the concentration on the fingers and the contact surface area_{skin}*S*. Similarly, the transfer back to the primary is a product of_{finger}*ε*_{skin}_{−}and_{skin}*S*with the concentration on the fingers of the secondary ._{finger}Therefore, the overall transfer between the primary and the secondary is

Wearing gloves reduces the skin-skin transfer coefficient by a factor of

Fomite touching. Individuals touch fomite which results in virus exchanges between fomites and hands. Transfer efficiencies denoted by

*ε*_{hand}_{−}(fomite to hand) and_{fomite}*ε*_{hand}_{−}(hand to fomite), were parameterized separately based on the results in Greene, et al._{fomite}^{43}, similar to the theoretical model of Kraay, et al.^{28}.The number of viruses that transfers from the fomites to the hands is the product of the

*ε*_{hand}_{−}with the the contact surface area (_{fomite}*S*) and the viable virus concentration on the fomite (_{finger}*C*). The number of viruses that transfers back to the fomite is the product of_{fomite}*ε*_{fomite}_{−},_{hand}*S*and the concentration on the finger._{finger}Therefore, the overall transfer from fomites to the hands following a single touch is given by, where

*I*is the individual and is either primary or secondary.Wearing gloves reduces the fomite-hand and hand-fomite transfer coefficient by a factor of , as above.

Touching surfaces. The large droplets that are emitted by coughs, sneezes, and speaks settle on environmental surfaces and contaminate them. When an individual touches an environmental surface his hand gets contaminated only if he touched a contaminated area. For simplicity, we assume that the probability to touch a given contaminated area is proportional to its size (

*S*) relative to the effective surface area of the room. The effective surface area is the sum of the room’s surface area (_{contaminatedArea}*S*) and the furniture surface area (_{room}*S*). The number of virus particles that are transferred to the hand with every touch is the product of the surface area of part of the hand that touches the surface_{furniture}*S*, the viable virus concentration on the surface_{finger}*C*and the fraction that is transported to the mouth following that touch_{surface}*ε*_{surface}_{−}. Accordingly,_{hand}Where

*P*is the probability to touch a contaminated area on an environmental surface and given by the fraction of the contaminated area (_{j}*S*) to the total surface area in the room, including the furniture (_{contaminatedArea}*S*+_{room}*S*). We neglect the transfer of virus particles back to environmental surfaces._{furniture}Wearing gloves reduces the fomite-hand transfer coefficient by a factor of , as above.

Hand washing. Washing hands removes the contamination from the hands of the individual with efficiency of

*ε*. Hence, the number of viruses that is eliminated following a hand washing is_{washing}Probability to become infected. The probability to become infected is calculated from the virus dosage that was accumulated within each time interval. That is, we divide the simulation into time periods and calculate the overall exposure for each time period. The probability that a susceptible individual will become infected was inferred from the dose-response curve that was reported for SARS-CoV-1

^{26}, where*D*is the total exposure of the individual at time period_{j}*j*(see Eqn 12 below).The length of each time period was determined from a Poisson process using a timescale of 7

*h*that represents the life cycle of the virus in the body, inferred from the life cycle of the SARS-COV-1^{44}.Decay of the virus on surfaces. The viability of the virus on surfaces outside of a host body (such as fomites, environmental surfaces, and the hands) decays with time. The decay rate depends, generally, on the surface type and was inferred from measurements

^{45}.Ventilation of contaminated air. The droplet nuclei that remain suspended are carried with the air currents in the room. For simplicity, we considered the entire house as a single, well mixed compartment. In households, the rooms are usually ventilated naturally (i.e., by the wind) and not by a heating, ventilation and air conditioning (HVAC) systems common in commercial and public buildings. The natural ventilation in households exchanges the air at a rate of 0.3

*h*^{−146, 47}, which leads to a decay of the air concentration of the droplet nuclei that carry the virus.

The above nine processes consist of both discrete random events (processes 1-7) and continuous events (processes 8 and 9). Hence, we used a hybrid continuous and stochastic-jump framework to model the dynamics of the transmission and infection processes.

The model consists of a system of coupled equations that describe the dynamics of the virus on the fingers of each individual, on fomites, in each contaminated area on an environmental area, its concentration in the air, and the overall exposure of the secondary individual from the airborne and contact routes. The concentration of the virus in each contaminated area on the environmental surface is tracked individually. A list of the parameters is provided in Table 2.

The dynamics of the concentration on the fingers of the individuals is determined by inoculation (process 1), physical contact (process 3), touching fomites (process 4), environmental surfaces (process 5), hand washing (process 6), and by the decay of the virus viability on the hands (process 8). Accordingly,

Where is the virus concentration on the fingers of individual I (I is either primary or secondary). The first 5 terms on the right hand side of Eqn. 8 describe the transfer of viruses as a result of the discrete events (processes 1,3,4,5, and 6) and the last term corresponds to the decay of the virus viability on the hands where *α _{hand}* is the decay rate constant. The

*δ*equals 1 for the primary individual and −1 for the secondary individual. The discrete event times (

*t*,

_{i}*with*

_{X}*where X is inoculation, fomite, surface, washing and PhysicalContact) are determined from a Poisson distribution with rate constant*

*τ*.

_{X}The dynamics of the average concentration of virus on the fomites in the room, *C _{fomit}*, is determined by touching fomites (process 4) and the decay of the virus viability (process 8),
where

*α*is the decay rate in the fomite.

_{fomite}The dynamics of the concentration in the air is governed by the emission of droplet nuclei during the coughing, sneezing, and talking of the primary individual (process 2) and the ventilation process (process 9),
where *C _{air}* is the virus concentration in the air and

*L*is the viral load of the primary individual, and are the number of coughs, sneezes, and talks (respectively) that occur at times

*t*,

_{i}*,*

_{cough}*t*,

_{i}*and*

_{sneeze}*t*,

_{i}*, and determined from a Poisson distribution with time constants of*

_{talk}*τ*,

_{cough}*τ*, and

_{sneeze}*τ*, respectively.

_{talk}The dynamics of the concentration in the *k*-th contaminated area on an environmental surface that was created by a cough, sneeze or during a talk at time *T _{k}* is determined by the equation

Where is the concentration on the *k*-th surface, * is the probability to touch a contaminated area and **α _{surface}* is the decay rate in this surface.

The *k*-th contaminated area was created at time *T _{k}* by a cough, sneeze or during speach. The initial concentration in the contaminated area is the product of the viral load of the primary individual at time

*T*with the volume of droplets that are larger than 100

_{k}*µm*(

*V*,

_{X}*where*

_{droplets}*X*is cough, sneeze or talk).

The exposure of the secondary is calculated within the exposure time interval and calculated as the sum of the exposure in the contact and droplet nuclei mode of transmission. Following Nicas and Best^{22}, we assume that only a fraction *ε _{mouth}* of the total number of virus particles that were deposited on the mouth reach oro-and nasopharyngel target sites. Similarly, only a fraction

*ε*of the inhaled viruses deposit in the respiratory tracts. Therefore, the exposure is given by the equation:

_{breath}Where *D ^{i}* is the exposure obtained when solving the equation between the two exposure events that take place at times and and

*B*is the breathing rate.

### Numerical method

In order to solve the model’s equations 8 to 12 we chose to use the jump-adapted approximation proposed by Casella et al.^{48}. For each realization, we determine the discrete event times (*t _{i}*,

*with where X is inoculation, fomite, surface, washing, coughing, Physical Contact, sneezing and talking) are determined from a Poisson distribution with rate constant*

_{X}*τ*. Then, the events are combined and sorted in an ascending order to obtain a set of discrete times. That is, , where each

_{X}*T*is assigned to a one

_{i}*t*,

_{i}*event and is the total number of discrete events.*

_{X}The simulation is then solved implicitly between the times *T _{i}* and

*T*

_{i}_{+1}with time step of 6 seconds. At each

*T*we calculate the number of virus particles that is transferred for the corresponding event that takes place at this point in time.

_{i}The model was implemented in python version 3.6.5.

### Code availability

The computer code that was used in calculating the results of this paper is available in https://github.com/yehudarav/CoronaIndoorTransmission.git

### Parameter estimation

We have used results reported in the literature to determine the evidence based values for the parameters in the model. However, in some cases direct measurements were not available for several parameters. Therefore, their value was estimated, based on additional assumptions. In the following refer to these and we provide the details of these assumptions and their justifications.

Time scale of the exposure: Very little is known about the dynamics of the SARS-CoV-2 virus in the human body. In order to estimate the time scale of exposure, we used the results reported by Qinfen, et al.^{44} regarding the life cycle of the SARS-CoV-1 in host cells. They found that the virus assembly and maturation was first detected around 7 hours post infection. Thus, it is plausible to use this characteristic time as the exposure time interval.

Exposure to infection factor for contact and droplet nuclei mode of transmission: It is quite possible that not all the virus particles that are inhaled or reach the facial membranes cause infection. Therefore, we assumed, similar to Nicas and Best^{22}, that the fraction of the exposure via the contact route that causes infection (*ε _{mouth}*) is 0.5.

The fraction of the inhaled dose that causes infection was estimated to be roughly 0.5 of the deposition functions which provides the retention of particles in the lungs^{49} and the nasal cavity, depending on particle size^{50}. These functions were applied to the particle size distribution reported by Chen, et al.^{42}, after a correction that takes into account evaporation^{3}.

Incubation time: We have used a weighted combination of the parameters of the incubation time distributions reported by Lauer, et al.^{24} and Bi, et al.^{10}.

## Data Availability

Theoretical numerical model, all parameters are provided in the manuscript

## Author contributions statement

All authors conceived the study and wrote the main manuscript. YA developed the model and did the simulation. EF contributed to numerical analysis. ZK contributed to parameter estimation. All authors contributed to the analysis and reviewed the manuscript.

## SUPPLEMENTARY INFORMATION

Theoretical investigation of pre-symptomatic SARS-CoV-2 person-to-person transmission in households

Yehuda Arav, Ziv Klausner, Eyal Fattal

### Sensitivity analysis

Some of the parameters’ values were obtained from studies that also reported the range of these values. Therefore, we have performed an extensive sensitivity analysis to check the robustness of the results. Specifically, we have examined the sensitivity of the model’s reconstruction of the serial interval distribution and the SAR for variations in the following parameters: dose response (Figure 5A), exposure time scale (Figure 5B), decay rate on surfaces (Figure 5C), surface area of the hand (Figure 5C), surface area touched (Figure 5D) and the median viral load (Figure 5E). As seen, the model’s results remains within the range of the values reported in the literature for the examined range of parameters.

We note that since the contribution of transmission modes 3 and 4 (Indirect contact transmission via surface and droplet nuceli transmission, respectively) is very small, the model is not sensitive to the parameters that are associated with them (such as the area of the room, the decay rate in air, breathing rate and etc.).

## Acknowledgements

The authors would like to thank Dr. Shay Weiss and Dr. Yehudah Alexander for their fruitful comments. This research was not funded.