Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

# The network limits of infectious disease control via occupation-based targeting

## Abstract

Policymakers commonly employ non-pharmaceutical interventions to reduce the scale and severity of pandemics. Of non-pharmaceutical interventions, physical distancing policies—designed to reduce person-to-person pathogenic spread – have risen to recent prominence. In particular, stay-at-home policies of the sort widely implemented around the globe in response to the COVID-19 pandemic have proven to be markedly effective at slowing pandemic growth. However, such blunt policy instruments, while effective, produce numerous unintended consequences, including potentially dramatic reductions in economic productivity. In this study, we develop methods to investigate the potential to simultaneously contain pandemic spread while also minimizing economic disruptions. We do so by incorporating both occupational and contact network information contained within an urban environment, information that is commonly excluded from typical pandemic control policy design. The results of our methods suggest that large gains in both economic productivity and pandemic control might be had by the incorporation and consideration of simple-to-measure characteristics of the occupational contact network. We find evidence that more sophisticated, and more privacy invasive, measures of this network do not drastically increase performance.

## Introduction

Containment of infectious disease requires implementation of strategies to reduce the spread of the disease1. Until a robust vaccine is widely available, as was the case for many months of the COVID-19 epidemic, the leading Non-Pharmaceutical Intervention (NPI) is physical distancing including closure of businesses2: a reduction of the number of people in physical proximity.

Since the European Industrial Revolution workers have increasingly migrated to towns and cities to convene in centralised workplaces3. A large literature of work relates growth, creativity and other metrics of economic performance directly to the increasing returns and economies of scale made possible by the density and proximity of workers in regions, cities and workplaces4,5,6. Therefore, physical distancing is unavoidably disruptive to these workers, workplaces and economies and is likely to cause long term shifts in working patterns. The highly infectious nature of COVID-19 along with the disease’s long incubation period7,8 have shown the cost of such disruption even more starkly than previous epidemics. As a result, working hour losses globally could be as high as 12% in the third quarter of 20209 and a 6.6% drop in global GDP in 202010. For this among other reasons, scholars have noted that the COVID-19 epidemic has the potential to affect the future of work substantially over the long term11,12.

More precisely, physical distancing policies affect workplaces directly and workers indirectly. While some work environments are conducive to physical distancing (e.g. construction sites), others are inherently not (e.g. gyms & nightclubs)13. As a result, physical distancing requirements have produced heterogeneous effects on workers that are primarily determined by workers’ ability to physically distance in their workplace and/or to work from home (WFH)14.

Thus policymakers have the difficult task of moving beyond coarse economic interventions e.g. furloughing all workers, to more targeted interventions that balance the human toll of infection with the economic cost of interventions. Stated precisely, workers should be partitioned into groups who may respectively (i) remain working in their workplaces (ii) work effectively from home or (iii) be furloughed, such that the spread of the epidemic is minimised while also minimising economic disruption associated with the intervention.

This optimisation task is currently impossible on a practical level for many reasons. These include the lack of sufficiently detailed information on human mobility and social interactions, the complex network dynamics of disease spread15, the interconnected nature of modern economies16,17,18, and divergent opinions on the proper balance between human and economic loss19.

Epidemics typically pose differing levels of risk to individuals based on demographic categories, such as age groups. This may be due to characteristics of the specific disease or due to differing physical contact patterns encoded in contact matrices20. This heterogeneity in the spread between members of a population also contributes to the overall complexity of accurately modeling the disease spread.

It is notoriously difficult to measure, at scale, physical proximity as it relates to airborne or droplet mediated transmission of infectious disease21,22. Deriving a network of contacts is even more challenging compared to a simple count of the number of contacts. This is due to privacy concerns and limitations of measurement accuracy as individuals need to be uniquely identified throughout the measurement23. Typically contact interactions are measured through self-reported surveys and networks measured using sensors24,25,26,27,28. Previous work has constructed contact matrices between subpopulations stratified by age and other demographics29,30. Full contact networks have been measured in large scale field studies in schools31,32, dormitories33, hospitals34, and conferences35 as well as being approximated from the use of location-based services36.

Yet, demographics represent only one set of individual characteristics that likely play a role in modulating differential epidemic spread. An individual’s work – their occupation – encapsulates many determinants of their epidemic risk. Commuting patterns, the extent of contact with others at their work, and the ease with which non-pharmaceutical interventions such as wearing of personal protective equipment and remote working can be adopted are all occupation-specific factors that shape both individual risk and the propensity for individuals in any given occupation to shape epidemic spread on human networks more generally.

While studies have measured contact networks in workplaces37, to date we are not aware of any in-depth study on contact networks stratified by occupation, and as such the true structure of such a network is unknown.

In this work, we examine the role of occupation in an epidemic spread in an urban environment. We also evaluate the efficacy of occupation-based disease control measures within a simulation of epidemic dynamics using the Susceptible-Exposed-Infected-Recovered (SEIR) framework. More specifically, we investigate the effects of several network-based interventions and compare them to the outcomes produced by more coarse heuristics. A focus on occupational interventions, coupled with detailed data on the distribution of the workforce across occupations, wage and workplace proximity allows us to simulate the economic impact of particular containment strategies alongside each intervention’s epidemiological impact.

Our methods enable us to approximate the degree of physical contact between individuals within and across occupations without explicit contact matrices (such matrices to our knowledge do not exist across the full empirical network due to the cost and methodological complexity of measuring such contacts). Our approach is general and applicable to any infectious disease, however, we parameterise our simulations with estimations of the epidemiological characteristics of COVID-19 for illustrative purposes.

Our work complements a rapidly growing body of work analysing the effects of physical distancing and furloughing policies on human mobility and behaviour in the COVID-19 epidemic13,38,39,40,41,42,43,44,45,46 with several focusing on the economic aspects in particular47,48.

Using these methods, we make three marked contributions to the science of epidemic control. First, we contribute a technique for constructing occupation-based epidemic simulations using publicly available data. Second, we compare several occupation-based containment policies based on their epidemiological and economic costs. Third, having identified contact degree as an effective heuristic for containment, we extend to investigate the fundamental limits of the controllability of epidemics in networks. Our findings have implications, not only for the containment of future epidemics, but also for the future of work and the interaction between disease risk at work, automation risk and the economics of labour.

## Data and methods

We present a general method for simulating occupation-based epidemic control policy interventions. However, we consider data from New York City as a paradigmatic example49. We make use of the Occupational Information Network (O*NET) which is a public repository of occupational data in the US and statistical economic information collected by the Bureau of Labour Statistics (BLS).

### O*NET: occupations and work characteristics

We make use of the O*NET data on “work characteristics” to derive a composite measure of proximity from five distinct work characteristics that are likely to be correlated with the degree of in-person contact a worker is required to have. These are Exposed to disease or infections, Performing for or working directly with the public, Communicating with persons outside of the organisation, Deal with external customers, Physical proximity. We construct a composite measure of proximity from these five dimensions using a Principal Component Analysis (PCA). We inspect the occupations that are assigned high and low degree to manually validate this method (see SI). The projection on the first principal component (PC) is listed in Fig. 1 for the top/bottom five jobs. We note that this first order approximate measure excludes certain dimensions that are likely to be important in preventing disease transmission in work environments such as the ability to incorporate protective measures into the workplace.

### Contact data: proportions of home/work/transit

In agreement with previous empirical studies of contact networks50, we define three categories of contacts among the adult working population: home, work and transport. The mean total contact degree has been found to be 75 in a recent study of proximity via Bluetooth in New York City51. Fixing the mean home degree at two based on the mean household size52, and considering an equal work and transportation mean degree, we set the ratio of the mean contact degree to be 4:73:73 for home:work:transport, respectively (see SI). We fit the empirical data found in30 to derive a log-normal distribution for each category rescaling the mean to fit the above ratios. For each worker node, we sample each of the three distributions to determine total degree. The number of work contacts for each worker depends on her occupation, while the number of home and transportation links are independent of occupation. Consider a worker in a high proximity occupation such as a Retail Salesperson, links are assigned as follows:

• The proximity score of the occupation is taken as the projection of the job on the first principal component derived from the PCA of the proximity variables. For the case of a retail Salesperson this is 0.55 in the range [0-1].

• This specific PCA value is mapped to a log-normal distribution (mean of 36.5) of work degree by mapping the percentiles of the distributions yielding a work degree of 49.

• A transport degree is assigned to the worker node by independently sampling a log-normal distribution with a mean of 36.5 and a standard deviation of 25.55 (see SI).

• A home degree is assigned to the worker by independently sampling a log-normal distribution with a mean of 2 and a standard deviation of 1.4.

Links from each category are functionally equivalent in terms of disease transmission in our simulations. However each layer is independent, and different policy interventions retain or remove links from each layer as described later.

### Work from home data

Some workers can work from home without disruption to productivity and thus may be effectively removed from the occupational contact network without loss of economic contribution. We consider a binary index for O*NET jobs derived by Dingel and Neiman14.

### Essentialness data

Since one motivation for NPI selectively furloughs workers according to perceptions of the essentialness of their contributions, we assign a measure of ‘essentialness’ using the data of del Rio-Chanona et al (2020)47,53.

### Wages data

Another NPI we consider here to balance between the loss of economic productivity and the epidemic expansion is the furlough of workers prioritising those on a low income. We use wage data by occupation for New York, based on the Occupational Employment Statistics survey by US Bureau Labor Statistics54. For most of the occupations, the dataset includes the average annual wages while for hourly-paid occupations like actors the dataset includes the mean hourly wages. For both cases, we estimated the average daily wage per occupation.

### Construction of contact network

We construct the contact network as follows.

1. 1.

We define N=200,000 representative workers each with a specified O*NET occupation as nodes in our contact network. The proportion of workers in each occupation matches workforce data for New York City. The total number of workers is the minimum number that allows for at least one worker of each occupation.

2. 2.

Each worker is assigned a contact degree for home, work and transport. The work contact degree is determined by the proximity score for the worker occupation. Home and transport contact degrees are drawn from fixed distributions independent of occupation (see SI).

3. 3.

The nodes are connected to form the contact network using the configuration model55 (see Fig. 2).

The synthetic networks produced by this approach have some resemblance to a configuration model. However, the configuration model typically samples from a degree distribution to assign stubs to nodes, whereas in our case the node degree sequence is known and determined by the distribution of occupations in the workforce.

Each node has as many half links as the number of home, work and transportation contacts of the worker represented by the node. Pairs of half links for each category of contacts are then connected randomly until all half-links are paired up generating the final contact network. We conduct epidemic simulations in one realisation of the contact network constructed by the configuration model. In SI, we repeat the analysis in a second realisation of the contact network drawn from the same configuration model. We make the assumption that the contact network is fixed throughout the epidemic period. While this is unrealistic, the changing patterns of physical colocation are prohibitively complex to model as to be out of scope of this paper.

### Epidemic simulation

Our simulation of worker based epidemic spread and containment proceeds as follows.

1. 1.

We construct an SEIR model56 in which each worker may be in one of the four epidemic compartments. At the initial time step all nodes are considered as susceptible. An infection is seeded with one randomly chosen infected node.

2. 2.

Once the epidemic has reached 200 nodes (i.e. 0.1% of the synthetic population), we implement one of the intervention strategies described below. Epidemics that terminate before reaching 200 nodes are discarded.

3. 3.

The simulation continues in daily time steps until the epidemic terminates.

4. 4.

For each epidemic simulation, we calculate the severity of the epidemic and the economic cost of furloughing workers and workers who are unproductive due to illness.

### Epidemic parameters

We emphasise that our methodology is designed for infectious disease epidemics in general. In this case we parameterise our SEIR simulations using the epidemic parameters within the reported ranges for the COVID-19 virus (see Table 1).

We make use of a standard SEIR model in which each worker is assigned to an epidemic state. Specifically, during the epidemic process, each node is in one of four states: susceptible (S), exposed (E), infected (I), and recovered (R). A susceptible node will become exposed for a certain period after being infected, and will transit into the infected state with a certain probability. At each time step, an infected node can recover with a fixed recovery rate. Both exposed and infected nodes can transmit the disease to a susceptible neighbor in the network with the same infection rate.

With a fixed basic reproduction number, $$R_0$$ (that is, the average number of secondary infections an infected individual causes before the individual recovers) and given the average degree, $$\langle d \rangle$$, of our synthetic contact network, we derive an estimated transmission rate per contact (see SI)

\begin{aligned} p_{inf} = 1 - (1 -( R_0 \langle d \rangle )^{(1/(\langle t_{recovery} \rangle + \langle t_{incubation} \rangle ))}) \end{aligned}

with the mean-field approximation assuming that each node has the same connectedness equal to the mean degree and recovers in $$\langle t_{recovery} \rangle + \langle t_{incubation} \rangle$$ days if the individual is infected.

Each epidemic simulation is seeded with one randomly infected node. Once the epidemic grows to more than 200 infected nodes, we implement a policy intervention into the network. If the epidemic dies out before reaching 200 nodes then the simulation is discarded, as we are here interested in modelling only those critical scenarios where human or policy interventions may be important factors in affecting epidemic outcomes.

### Intervention strategies

We seek to find the dually optimal strategy – both in infection control and economic terms – that a policy maker might use to contain the epidemic. This requires leveraging occupation-specific information to selectively remove worker nodes from the contact network, i.e. furloughing workers. This must be done in such a way that both the severity of the epidemic, as measured by the peak of the number of infected workers, and the economic cost of furloughing workers is jointly minimised. In our simulation study, we furlough workers based on the intervention strategies shown in Table 2.

### Economic productivity

Approximation of the economic impact of a worker unable to work either due to infection or furlough is a challenging measurement task66. We investigated two metrics: (i) the average wage of the furloughed/infected worker and (ii) the contribution to macro-economic productivity derived from the occupational share of the larger industry-level productivity. Due to the coarseness of the best available data (and the subsequent risk of substantial measurement error) for (ii) (see SI) we proceed with measure (i).

The total cost of lost productivity of infected workers is calculated as the sum of the daily wages over the period that each worker is infected (the daily wage of each worker is based on his/her occupation). The total cost of workers who are furloughed is calculated as the sum of their daily wages over the furloughed period (the furloughed period is the difference between the length of the epidemic period and the length of the intervention period).

In the former case (wage-based calculation), the wage of a worker is considered proportional to her economic contribution and her loss of income as a loss in disposable income that can be spent to stimulate supply. We acknowledge that wage is an imperfect measure of the economic contribution of a job. However we conclude from the best quality data on industry level economic contributions that estimation is not possible (see SI). Results in the main paper correspond to (i).

### Ethics declaration

All analyses including human data followed relevant ethical guidelines.

## Results

Figure 3 presents a comparison between the ‘no intervention’ and ‘furloughing 10% based on degree’ strategies, respectively. Under an uninhibited outbreak (strategy -1) the peak of infected nodes occurs at a mean of 129 days and with a mean 19% of nodes infected which provides reasonable correspondence to observed outbreak dynamics. This changes to a mean of 164 days and 6% infected respectively under the strategy in which workers in occupations that can work from home have transport and work ties cut. Due to the difficulty in meaningfully representing the full contact network, we have instead represented the full network aggregated to  ~ 700 nodes representing occupations while still capturing both (i) the number of workers in each occupation (the node size) and (ii) the proportion that are infected (the proportion of the node colour).

Figure 4 compares each strategy as a function of the stringency of the strategy (percentage of workers furloughed). In particular, the economic cost of furloughing workers calculated from wage (left) the economic cost of infected workers (middle) and the size of the peak of infected workers (right). We consider these distinct measures of epidemic severity as independent measures that policy makers would like to jointly minimise. The strategies are ordered in terms of total cost of furloughing across the range of stringency; most to least expensive, top to bottom. Generally, the measured economic loss due to infection is an order of magnitude less than the loss due to furloughing, since the infection period is ~ 14 days while the furlough period is the duration of the entire epidemic process (from the day that the intervention is triggered until there are no infected or exposed individuals); it is cheaper to allow a worker to be infected than to pay her wages for the duration of the epidemic.

We acknowledge that epidemic control is a careful balance of economic and human losses that is determined by some unknown and subjective multiplier of the pure cost of wages of infected workers due to secondary effects of infection e.g. lasting after effects of the disease, deaths caused by the disease, loss of quality years of life, psychological effects and so on. Likewise, furlough has negative secondary effects on individual workers66,67. Further, we recognise that the implementation of these strategies would be a logistic challenge that might favour a lower level of stringency. The range of severity of mobility restrictions has varied vastly between countries, from a laissez-faire approach (Sweden) to a ‘COVID-zero’ strategy warranting total lockdown (New Zealand). Our rationale is to thoroughly explore the behaviour across the phase space. For this same reason we do not attempt to directly compare losses on these three distinct dimensions or to prescribe the extent to which workers should be furloughed. Instead we seek to compare the general behaviour of various interventions in terms of the complex dynamics of the contact network, labour market and epidemic dynamics. We also note that the economic cost includes both the case of the state bearing the cost of covering the salaries of furloughed workers i.e. a strong social safety net and also the case that the worker herself must use savings, family support etc to replace lost wages. We consider the performance of each strategy across the entire range of severity of the strategy implementation i.e. the percentage of workers who are furloughed.

The loss due to furloughing shows a peak, at around 50% of removed nodes, for the worst performing strategies (left column rows a-c). This is related to the duration of the epidemic process (the size of the plot markers) as workers are furloughed until the epidemic process terminates. When an intermediate number of workers are removed, the epidemic is still able to spread but takes longer to fully terminate. Note that this is not attributable to herd immunity as the proportion of recovered nodes remains well below the 70-85% range. The better performing strategies (rows d-f) remove nodes more effectively, causing the epidemic to terminate sooner. In contrast to the other strategies, removing more than ~ 30% of nodes based on controllability, centrality and degree strategies (d-f) does not cause the epidemic period to increase significantly. This leads to a roughly linear increase in furlough costs as more workers are furloughed for the comparable duration.

In contrast, the loss due to infected workers (middle column) decreases monotonically; the total number of infected workers decreases as network edges are cut. Likewise the size of the peak of the infection (right column), representing the strain on health facilities, decreases monotonically as the ‘curve is flattened’.

It should be noted that the varying behaviour of the strategies is generally not trivially attributable to the breaking apart of the underlying worker network, particularly as furloughing cuts only work and transport links preserving home contacts analogously to a home quarantine (see SI). For each strategy, the size of the Largest Connected Component decreases at roughly the same rate and contains no less than 80% of nodes across all strategies and the range of severity. Our results are broadly unchanged when home links are also severed i.e. a centralised quarantine (see SI).

We find that removing workers according to essentialness underperforms even random node removal. This is consistent with only a moderate correlation ($$\rho = 0.19$$, see SI) between our proximity score and essentialness score; non-essential workers do not necessarily have a high contact degree and so their removal still allows efficient epidemic spread. Likewise removing based on wage alone performs poorly, consistent with a weak correlation between wage and proximity ($$\rho = 0.05$$, see SI). For example, construction or mechanical occupations have low proximity to others. As a result, the epidemic can spread quite effectively until a large proportion of workers are furloughed (n% $$\sim$$ 100%).

Removing workers based on contact degree shows a drastic increase in performance (38.6% of economic cost of the worst performing strategy; removing nodes according to essentialness) when considering the cost of infected workers. Notably, the best performing centrality metric is almost indistinguishable from node degree (38.6% vs 38.6%). Despite centrality metrics incorporating full knowledge of the network structure, they were unable to halt the spread of the epidemic markedly more effectively than simple knowledge of a workers local environment. Likewise for the controllability based strategy (38.6% vs. 38.0%). Figure 5 evaluates the overall performance of each strategy across the range of stringency (n%). We repeat our analysis using an alternate random network instantiation created using the same methodology as described above. We find our results are robust to these alternative specifications (see SI).

Control nodes of a network are defined as the nodes which, if their state can be manipulated, allows for the state of the entire network to be driven to a desired configuration. Under the model of switchboard controllability68 (considered the most appropriate model for epidemics), we find that 47.9% of nodes in our network are identified as control nodes. Such a large proportion implies that the contact degree has a low capacity to drive controllability, which is consistent with previous studies of social networks68 in comparison to more explicitly hierarchical networks such as those belonging to organisations or neuronal structures. Consequently, we find that a controllability based strategy offers no significant improvements in performance over a degree based strategy. This is despite incorporating full knowledge of the network structure.

These empirical results illustrate the difficulty of fully controlling an epidemic ex-ante, even with the advantage of full information of the contact network and justify the use of manual38,59,69 and digital contact tracing70,71 for containment as a local and dynamic strategy (as opposed to static and global one). Correspondingly, we find that removing control nodes preferentially is not a strategy that performs well, which suggests that the dynamics of epidemics might not amenable to network controllability68.

Given that our contact network is not empirical, it is possible that our results arise as an artifact of how the network was constructed from occupational characteristics. We emphasise that to our knowledge there does not exist a publicly available contact network larger than $$\sim 10^2$$ nodes, and none of any size that considers occupational subpopulations, that could be employed directly or bootstrapped against. Nevertheless we compare our generated network to a number of smaller empirical networks from different contexts including workplaces, schools and conferences (see SI). We find congruence between most metrics excluding those that were not explicitly coded in our configuration model such as transitivity and assortativity. This includes the high density of control nodes, in the range 40-50%.

## Discussion

In this paper we have described a method to construct epidemic contact matrices between detailed occupational sub-populations using publicly available data. We use these contact matrices as input to SEIR simulations in an urban environment and compare various epidemic containment strategies in terms of their economic and human cost. We emphasise that our goal in this paper is not to predict the spread of COVID-19, but rather to compare various occupation-based containment strategies more generally using the COVID-19 epidemic as a motivating case. We find that the heuristic of worker node removal according to network degree (the number of physical contacts that a worker has) performs approximately the same as more complex metrics based on complete network structure or other occupational characteristics. More broadly we note that epidemic contact networks exhibit low levels of controllability.

Our findings demonstrate that the structure of the contact network heavily influences disease dynamics in non-trivial ways. For example, furloughing a small proportion of workers can lead to pruning of the network in such a way that the epidemic persists for a long time, albeit at low levels, leading to a long and costly furlough. Intuitive strategies such as furloughing workers based on the essentialness of their job, by wage or at random all perform poorly on this basis. In contrast, network-based metrics such as degree and centrality are able to reduce the peak of the infection (flattening the curve) and also reduce the epidemic period.

Our findings have been shown to apply in New York City; a large, dense and wealthy city with a diverse workforce. We would expect that other dense megapolises with widespread use of public transport to share similar mixing patterns. However the distribution of occupations are likely to be very different depending on urban economic specialisations. In general, larger cities support more diverse job markets, and vica versa for smaller cities. Therefore we might expect more homogeneity in the contact degree across workers in smaller urban areas and rural areas, this would lead to lower efficacy in degree based strategies. Likewise intra-household mixing is likely to play a larger role when cultural and economic constraints lead to larger family units.

Our methods exist in a relative paucity of sufficiently rich data to simulate real-world dynamics with high fidelity. Empirical contact networks, whether self-reported or passively sensed, are only freely available on a small scale. Networks with any demographic information are rarer still, while contact networks with occupational information are not to our knowledge available at all. In a workplace the correspondence between physical contact and social ties is likely to be weaker as some jobs, particularly service jobs, require a great deal of face to face contact. Consequently social ties are likely not to be a trivial proxy for physical contact ties, as observed previously36. Thus social networks are not a viable proxy for worker contact networks. In general we find that occupational attributes relevant to epidemic containment (proximity, essentialness) are only weakly related to wage. As a result there is a non-trivial relationship between occupational characteristics and epidemic dynamics that make the dual optimisation of economic and human costs challenging. Finally we note, we have made use of the best available data on occupations linked to wages; this necessarily excludes informal work and work in kind e.g. as a housewife or househusband.

While we have successfully attributed much of our results to node degree, we expect that our results would differ if our method were repeated in a real contact network reflecting transitivity and hierarchical organisational structure. While the role of empirical structure in network spreading is non-trivial, it is likely that the global spread would be slower72 than observed in our configuration model. As such we should consider the rate and depth of epidemic spread reported here to be a lower bound.

A more precise critique is that our configuration model connects edge stubs between nodes at random, once a degree has been assigned based on occupation. It might be expected that the edge structure in an empirical contact network due to organisational structure or other heterogeneities might lead to very different behaviour than those we have observed. Shuffling the edges of our contact networks i.e. an Erdos-Renyi network with equivalent average edge density, we confirm that considerable structure present in our modified configuration model is destroyed. This can be seen by comparing the distribution of various centrality measures (see SI). Likewise our benchmark empirical contact networks exhibit comparable behaviour. The finding that the function of a network is to large extent determined by the node degree sequence, and less by the precise connections is consistent with findings on the controllability of a large corpus of networks68.

An additional caveat relates to the inherent difficulty of assessing the full economic effect of workers being unable to work. While such epidemic containment policies are generally considered on the basis of cities or states, industries are dynamically linked through national economies and global supply chains. The full economic and social cost of epidemics and associated non-pharmaceutical interventions is so complex to model as to be out of scope of this single paper.

It is also certainly true that social networks are not static73. However there is no consensus on how social ties change during lockdowns and certainly not on how this is related to occupations which is the main focus of our paper. We believe that the abrupt, disruptive and heterogeneous nature of the COVID containment measures are likely to be difficult to model in a consistent way. For example, previous work has shown that the abrupt and short trauma of a natural disaster can lead to an increase in future social ties among a student population74. However the prolonged, and for the majority of the population of high income countries, relatively less traumatic nature of the COVID pandemic it is unclear how social networks and contact networks might each be affected. It could be argued that since many social contacts are at least partially maintained at a distance e.g. through phone calls, social media, video chat; that these contacts would not be greatly disrupted. Yet it is important to distinguish between social network ties (friendships) and contact network ties (physical interactions); our interest here is not in social networks so much as contact networks; these are not trivially related36 Likewise, various occupations work together in workplaces or within supply chains to produce goods at different regional, and even global, scales. Thus, the productivity of a worker is never truly fully independent of another. Practically this means that furloughing one worker could lead to another being effectively furloughed and in this case our estimate of economic cost would represent an underestimate.

Subject to the limitations noted above, our findings may have important implications for the practical implementation of epidemic control. Public health services have a limited number of tools available for epidemic containment, primarily furloughing of workers, contact tracing of identified infected individuals and pre-emptive surveillance of citizens. The precise implementations of these strategies require careful trade-offs between economic, social and privacy costs.

Our findings make the case for better data on contact patterns broken down by occupation, as well as thorough classification of occupations and workplaces in terms of their ability to incorporate mitigating strategies e.g. Personal Protective Equipment and frequent testing. However we demonstrate that a policy of furloughing workers based on total number of contacts is far more effective at minimising epidemic spread and the cost of furloughing workers as compared to other heuristics such as furloughing based upon essentialness or furloughing a proportion of the workforce at random. However we find that there are diminishing returns to increased knowledge of the contact network beyond degree. Heuristics such as centrality or control nodes are not able to significantly outperform a simple count of contact degree. This is despite a huge increase in both the amount of invasive personal data required to reconstruct a full contact network and the complexity of its collection. Practically speaking, worker degree could be straightforwardly estimated using Bluetooth proximity to other phones, or GPS co-location. Each contact need not be identified beyond a unique untraceable ID to uniquely identify repeated contacts and no further information need be shared between individual devices or with a centralised authority.

The COVID-19 epidemic has caused many profound societal changes that are unlikely to be reversed even once the disease abates due to mass vaccination. This includes scrutinisation of the nature of work in light of vast changes in demand across sectors, variable infection across occupations, the large-scale adoption of remote working and challenging deeply ingrained understandings of workplaces.

The ability to successfully automate the skills of workers, often known as skills based technological change, depends not only on overcoming the underlying physical or engineering challenges. Moreover the automating technology must be economical and must be socially and politically acceptable in order to successfully replace human labour. The changes put in motion by the COVID-19 pandemic are likely to drastically alter these considerations, in some cases leading to ‘automation forcing’11.

Presently ‘goods producing’ jobs that can be performed in isolation e.g. construction workers or truck drivers are considered to be at high risk of automation. In contrast ‘service producing’ jobs with a high degree of personal contact e.g. fitness instructors or nannies, are at relatively lower risk. Consequently an increase in employment share has been observed in the latter category75.

We have demonstrated that the nature and number of connections between workers; the structure of the contact network, has vastly more potential to effectively stop epidemic spread than a single attribute of each occupation such as wage or essentialness. Looking to the future, we might reasonably expect that research and development of automation technologies will refocus to target the skills present in occupations with a high degree of contact with others. This could in turn lead to a change in these prevailing dynamics within low wage jobs.

## References

1. 1.

Fowler, J. H., Hill, S. J., Levin, R. & Obradovich, N. Stay-at-home orders associate with subsequent decreases in covid-19 cases and fatalities in the united states. PLoS ONE 16, 1–15. https://doi.org/10.1371/journal.pone.0248849 (2021).

2. 2.

Brauner, J. M. et al. Inferring the effectiveness of government interventions against COVID-19. Science 371, (2020).

3. 3.

Glaeser, E. Viewpoint: Triumph of the City. J. Transp. Land Use 5, (2012).

4. 4.

Bettencourt, L. M. A., Lobo, J., Helbing, D., Kuhnert, C. & West, G. B. Growth, innovation, scaling, and the pace of life in cities. Proc. Natl. Acad. Sci. 104, 7301–7306 (2007).

5. 5.

Neffke, F. M. H. The value of complementary co-workers. Sci. Adv. 5, eaax3370 (2019).

6. 6.

Claudel, M., Massaro, E., Santi, P., Murray, F. & Ratti, C. An exploration of collaborative scientific production at MIT through spatial organization and institutional affiliation. PLoS ONE 12, e0179334 (2017).

7. 7.

Hu, B., Guo, H., Zhou, P. & Shi, Z.-L. Characteristics of SARS-CoV-2 and COVID-19. Nat. Rev. Microbiol. 19, 141–154 (2020).

8. 8.

Zhang, J. et al. Evolving epidemiology and transmission dynamics of coronavirus disease 2019 outside Hubei province, China: a descriptive and modelling study. Lancet Infect. Dis. 20, 793–802 (2020).

9. 9.

International Labour Organization. ILO monitor: COVID-19 and the world of work. https://www.ilo.org/wcmsp5/groups/public/@dgreports/@dcomm/documents/briefingnote/wcms_749399.pdf (2020). Accessed 17 Oct 2020.

10. 10.

The Economist. What is the economic cost of covid-19? https://www.economist.com/finance-and-economics/2021/01/09/what-is-the-economic-cost-of-covid-19 (2021). Accessed 18 Jan 2021.

11. 11.

Autor, D. & Reynolds, E. The nature of work after the COVID crisis: Too few Low-Wage jobs. https://www.hamiltonproject.org/assets/files/AutorReynolds_LO_FINAL.pdf (2020). Accessed 10 May 2021.

12. 12.

Kniffin, K. M. et al. COVID-19 and the workplace: Implications, issues, and insights for future research and action. Am. Psychol. 76, 63–77 (2020).

13. 13.

Chang, S. et al. Mobility network models of COVID-19 explain inequities and inform reopening. Nature 589, 82–87 (2021).

14. 14.

Dingel, J. I. & Neiman, B. How many jobs can be done at home?. J. Public Econ. 189, 104235 (2020).

15. 15.

Pastor-Satorras, R., Castellano, C., Van Mieghem, P. & Vespignani, A. Epidemic processes in complex networks. Rev. Modern Phys. 87, 925–979 (2015).

16. 16.

Farmer, J. D. Economics needs to treat the economy as a complex system. https://www.ineteconomics.org/uploads/papers/farmer_berlinpaper.pdf (2012). Accessed 17 Oct 2020.

17. 17.

Hidalgo, C. A. & Hausmann, R. The building blocks of economic complexity. Proc. Natl. Acad. Sci. 106, 10570–10575 (2009).

18. 18.

Schweitzer, F., Fagiolo, G., Sornette, D., Vega-Redondo, F. & White, D. R. Economic networks: What do we know and what do we need to know?. Advs. Complex Syst. 12, 407–422 (2009).

19. 19.

Trust and the Covid-19 Pandemic. https://www.edelman.com/sites/g/files/aatuss191/files/2020-05/2020%20Edelman%20Trust%20Barometer%20Spring%20Update.pdf (2020). Accessed 18 Jan 2021.

20. 20.

Chowell, G., Tariq, A. & Hyman, J. M. A novel sub-epidemic modeling framework for short-term forecasting epidemic waves. BMC Med. 17, (2019).

21. 21.

Read, J. M., Edmunds, W. J., Riley, S., Lessler, J. & Cummings, D. A. T. Close encounters of the infectious kind: methods to measure social mixing behaviour. Epidemiol. Infect. 140, 2117–2130 (2012).

22. 22.

Stopczynski, A., Pentland, A. s. & Lehmann, S. How physical proximity shapes complex social networks. Sci. Rep. 8 (2018).

23. 23.

Stopczynski, A., Pietri, R., Pentland, A., Lazer, D. & Lehmann, S. Privacy in Sensor-Driven human data collection: a guide for practitioners (2014).

24. 24.

Grijalva, C. G. et al. A household-based study of contact networks relevant for the spread of infectious diseases in the highlands of Peru. PLoS ONE 10, e0118457 (2015).

25. 25.

Béraud, G. et al. The French connection: the first large population-based contact survey in France relevant for the spread of infectious diseases. PLoS ONE 10, e0133203 (2015).

26. 26.

Ajelli, M. & Litvinova, M. Estimating contact patterns relevant to the spread of infectious diseases in Russia. J. Theor. Biol. 419, 1–7 (2017).

27. 27.

van Hoek, A. J. et al. The social life of infants in the context of infectious disease transmission; social contacts and mixing patterns of the very young. PLoS ONE 8, e76180 (2013).

28. 28.

Sun, L., Axhausen, K. W., Lee, D.-H. & Huang, X. Understanding metropolitan patterns of daily encounters. Proc. Natl. Acad. Sci. USA 110, 13774–13779 (2013).

29. 29.

Prem, K., Cook, A. R. & Jit, M. Projecting social contact matrices in 152 countries using contact surveys and demographic data. PLoS Comput. Biol. 13, e1005697 (2017).

30. 30.

Mossong, J. et al. Social contacts and mixing patterns relevant to the spread of infectious diseases. PLoS Med. 5, e74 (2008).

31. 31.

Salathé, M. et al. A high-resolution human contact network for infectious disease transmission. Proc. Natl. Acad. Sci. USA 107, 22020–22025 (2010).

32. 32.

Stehlé, J. et al. High-resolution measurements of face-to-face contact patterns in a primary school. PLoS ONE 6, e23176 (2011).

33. 33.

Aharony, N., Pan, W., Ip, C., Khayal, I. & Pentland, A. Social fMRI: investigating and shaping social mechanisms in the real world. Pervasive Mob. Comput. 7, 643–659 (2011).

34. 34.

Vanhems, P. et al. Estimating potential infection transmission routes in hospital wards using wearable proximity sensors. PLoS ONE 8, e73970 (2013).

35. 35.

Cattuto, C. et al. Dynamics of Person-to-Person interactions from distributed RFID sensor networks. PLoS ONE 5, e11596 (2010).

36. 36.

Coviello, L., Franceschetti, M., García-Herranz, M. & Rahwan, I. Predicting and containing epidemic risk using on-line friendship networks. PLoS ONE 14, e0211765 (2019).

37. 37.

Potter, G. E., Smieszek, T. & Sailer, K. Modeling workplace contact networks: the effects of organizational structure, architecture, and reporting errors on epidemic predictions. Netw. Sci. (Camb. Univ. Press) 3, 298–325 (2015).

38. 38.

Aleta, A. et al. Modelling the impact of testing, contact tracing and household quarantine on second waves of COVID-19. Nat. Hum. Behav. 4, 964–971 (2020).

39. 39.

Block, P. et al. Social network-based distancing strategies to flatten the COVID-19 curve in a post-lockdown world. Nat. Hum. Behav. 4, 588–596 (2020).

40. 40.

Benzell, S. G., Collis, A. & Nicolaides, C. Rationing social contact during the COVID-19 pandemic: Transmission risk and social benefits of US locations. Proc. Natl. Acad. Sci. USA 117, 14642–14644 (2020).

41. 41.

Oliver, N. et al. Mobile phone data for informing public health actions across the COVID-19 pandemic life cycle. Sci. Adv. 6, eabc0764 (2020).

42. 42.

Chinazzi, M. et al. The effect of travel restrictions on the spread of the 2019 novel coronavirus (COVID-19) outbreak. Science 368, 395–400 (2020).

43. 43.

Maier, B. F. & Brockmann, D. Effective containment explains subexponential growth in recent confirmed COVID-19 cases in China. Science 368, 742–746 (2020).

44. 44.

Kraemer, M. U. G. et al. The effect of human mobility and control measures on the COVID-19 epidemic in China. Science 368, 493–497 (2020).

45. 45.

Holtz, D. et al. Interdependence and the cost of uncoordinated responses to COVID-19. Proc. Natl. Acad. Sci. USA 117, 19837–19843 (2020).

46. 46.

Mistry, D. et al. Inferring high-resolution human mixing patterns for disease modeling. Nat. Commun. 12, 323 (2021).

47. 47.

del Rio-Chanona, R. M., Mealy, P., Pichler, A., Lafond, F. & Farmer, J. D. Supply and demand shocks in the COVID-19 pandemic: An industry and occupation perspective. Oxford Rev. Econ. Policy 36, S94–S137 (2020).

48. 48.

Acemoglu, D., Chernozhukov, V., Werning, I. & Whinston, M. Optimal targeted lockdowns in a Multi-Group SIR model (2020). National Bureau of Economic Research Working Paper 27102.

49. 49.

U.S. Bureau of Labor Statistics. Occupational employment and wage statistics. https://www.bls.gov/oes/current/oes_35620.htm. Accessed 10 May 2021.

50. 50.

Social contact data. http://www.socialcontactdata.org/data/. Accessed 17 Oct 2020.

51. 51.

Bakker, M., Berke, A., Groh, M., Pentland, A. & Moro, E. Effect of social distancing measures in the New York City metropolitan area. https://connection.mit.edu/sites/default/files/publication-pdfs/Effect%20of%20social%20distance%20measures%20in%20social%20distancing%20in%20the%20NY%20area.pdf (2020). Boston: Massachusetts Institute of Technology, Accessed 10 May 2021.

52. 52.

Ellen, I. G. & O’Flaherty, B. Social programs and household size: evidence from New York city. Popul. Res. Policy Rev. 26, 387–409 (2007).

53. 53.

del Rio Chanona, R. M., Mealy, P., Pichler, A., Lafond, F. & Farmer, J. D. Data accompany the paper “Supply and demand shocks in the COVID-19 pandemic: An industry and occupation perspective”. https://zenodo.org/record/3751068#.YThykC0RrfY (2020). Accessed 16 Oct 2020.

54. 54.

55. 55.

Newman, M. Networks: An Introduction (Oxford University Press, 2010).

56. 56.

Brauer, F., van den Driessche, P. & Wu, J. Mathematical Epidemiology (Springer, Berlin, 2008).

57. 57.

Anderson, R. M., Heesterbeek, H., Klinkenberg, D. & Hollingsworth, T. D. How will country-based mitigation measures influence the course of the COVID-19 epidemic?. Lancet 395, 931–934 (2020).

58. 58.

Imai, N. et al. Transmissibility of 2019-nCoV. https://www.imperial.ac.uk/media/imperial-college/medicine/sph/ide/gida-fellowships/Imperial-College-COVID19-transmissibility-25-01-2020.pdf (2020). Technical Report 3, Imperial College London, UK, Accessed 10 May 2021.

59. 59.

Hellewell, J. et al. Feasibility of controlling COVID-19 outbreaks by isolation of cases and contacts. Lancet Glob. Health 8, e488–e496 (2020).

60. 60.

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 (2020). Accessed 14 Oct 2020.

61. 61.

Li, Q. et al. Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia. N. Engl. J. Med. 382, 1199–1207 (2020).

62. 62.

Lauer, S. A. et al. The incubation period of coronavirus disease 2019 (COVID-19) from publicly reported confirmed cases: estimation and application. Ann. Intern. Med. 172, 577–582 (2020).

63. 63.

Linton, N. M. et al. Incubation period and other epidemiological characteristics of 2019 novel coronavirus infections with right truncation: A statistical analysis of publicly available case data. J. Clin. Med. Res. 9, (2020).

64. 64.

Nepusz, T. & Vicsek, T. Controlling edge dynamics in complex networks. Nat. Phys. 8, 568–573 (2012).

65. 65.

Nepusz, T. netctrl: controllability of complex networks with node and edge dynamics. https://github.com/ntamas/netctrl.

66. 66.

Baranik, L. E., Cheung, J. H., Sinclair, R. R. & Lance, C. E. What happens when employees are furloughed? a resource loss perspective. J. Career Dev. 46, 381–394 (2019).

67. 67.

Lee, S. & Sanders, R. M. Fridays are furlough days: the impact of furlough policy and strategies for human resource management during a severe economic recession. Rev. Public Personnel Adm. 33, 299–311 (2013).

68. 68.

Liu, Y.-Y., Slotine, J.-J. & Barabási, A.-L. Controllability of complex networks. Nature 473, 167–173 (2011).

69. 69.

Huerta, R. & Tsimring, L. S. Contact tracing and epidemics control in social networks. Phys. Rev. E 66, 056115 (2002).

70. 70.

Farrahi, K., Emonet, R. & Cebrian, M. Epidemic contact tracing via communication traces. PLoS ONE 9, e95133 (2014).

71. 71.

Cebrian, M. The past, present and future of digital contact tracing. Nat. Electron. 4, 2–4 (2021).

72. 72.

Karsai, M. et al. Small but slow world: how network topology and burstiness slow down spreading. Phys. Rev. E 83, 025102 (2011).

73. 73.

Miritello, G., Lara, R., Cebrian, M. & Moro, E. Limited communication capacity unveils strategies for human interaction. Sci. Rep. 3, 1–7 (2013).

74. 74.

Phan, T. Q. & Airoldi, E. M. A natural experiment of social network formation and dynamics. Proc. Natl. Acad. Sci. 112, 6595–6600 (2015). https://www.pnas.org/content/112/21/6595.

75. 75.

Autor, D. H. & Dorn, D. The growth of low-skill service jobs and the polarization of the US labor market. Am. Econ. Rev.Am. Econ. Rev. 103, 1553–1597 (2013).

## Acknowledgements

We thank Dr. Weihua Li for contributions made to earlier stages of this project. The authors were unable to contact Dr Li to receive approval for submission of the completed manuscript. We were therefore unfortunately precluded by journal submission policies from including Dr Li as an author. We also acknowledge Jurgen Rossbach for help with editing figures. This research made use of the Rocket High Performance Computing service at Newcastle University.

## Funding

Open Access funding enabled and organized by Projekt DEAL.

## Author information

Authors

### Contributions

D.A.: literature search, formal analysis, methodology, software, validation, visualisation, writing - original draft N.O.: conceptualisation, formal analysis, methodology, writing - review & editing N.P.: literature search, methodology, investigation, writing - review & editing M.C.: conceptualisation, methodology, supervision, writing - original draft, writing - review & editing A.R.: conceptualisation, investigation, methodology, project administration, supervision, validation, visualisation, writing - original draft, writing - review & editing.

### Corresponding authors

Correspondence to Manuel Cebrian or Alex Rutherford.

## Ethics declarations

### Competing interests

The authors declare no competing interests.

### Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Avraam, D., Obradovich, N., Pescetelli, N. et al. The network limits of infectious disease control via occupation-based targeting. Sci Rep 11, 22855 (2021). https://doi.org/10.1038/s41598-021-02226-x

• Accepted:

• Published:

• DOI: https://doi.org/10.1038/s41598-021-02226-x