key: cord-0837746-lhv8u6ef
authors: Schimit, P.H.T.
title: A model based on cellular automata to estimate the social isolation impact on COVID-19 spreading in Brazil
date: 2020-11-09
journal: Comput Methods Programs Biomed
DOI: 10.1016/j.cmpb.2020.105832
sha: e89544f096e984d7fa31405beddb1675d1b7dd90
doc_id: 837746
cord_uid: lhv8u6ef

Background and objective Many countries around the world experienced a high increase in the number of COVID-19 cases after a few weeks of the first case, and along with it, excessive pressure on the healthcare systems. While medicines, drugs, and vaccines against the COVID-19 are being developed, social isolation has become the most used method for controlling the virus spreading. With the social isolation, authorities aimed to slow down the spreading, avoiding saturation of the healthcare system, and allowing that all critical COVID-19 cases could be appropriately treated. By tuning the proposed model to fit Brazil’s initial COVID-19 data, the objectives of the paper are to analyze the impact of the social isolation features on the population dynamics; simulate the number of deaths due to COVID-19 and due to the lack of healthcare infrastructure; study combinations of the features for the healthcare system does not collapse; and analyze healthcare system responses for the crisis. Methods In this paper, a Susceptible-Exposed-Infected-Removed model is described in terms of probabilistic cellular automata and ordinary differential equations for the transmission of COVID-19, flexible enough for simulating different scenarios of social isolation according to the following features: the start day for the social isolation after the first death, the period for the social isolation campaign, and the percentage of the population committed to the campaign. Results Results showed that efforts in the social isolation campaign must be concentrated both on the isolation percentage and campaign duration to delay the healthcare system failure. For the hospital situation in Brazil at the beginning of the pandemic outbreak, a rate of 200 purchases per day of intensive care units and mechanical ventilators is the minimum rate to prevent the collapse of the healthcare system. Conclusions By using the model for different scenarios, it is possible to estimate the impact of social isolation campaign adhesion. For instance, if the social isolation percentage increased from 40% to 50% in Brazil, the purchase rate of 150 intensive care units and mechanical ventilators per day would be enough to prevent the healthcare system to collapse. Moreover, results showed that a premature relaxation of the social isolation campaign can lead to subsequent waves of contamination.

range of symptoms, with some cases presenting no symptoms, and others presenting fever and severe respiratory problems [60] , complicating health advice at the onset of symptoms, for instance [40] ; and the third is the colossal number of cases requiring medical treatment at the same time, collapsing many healthcare systems around the globe [15, 48, 51] .

Governments, scientists, and specialists have proposed protocols for dealing with the situation, like washing hands carefully, a proper sneezing and coughing etiquette, self-isolation, social distancing, quarantine, and complete lockdown [19] . Due to the rapid escalation of cases, countries quickly moved to the lockdown to stop the COVID-19 spreading [48, 51] . The term lockdown has been used for describing different scenarios, but it means the "banning of any non-essential public gatherings, closure of educational and public/cultural institutions, ordering people to stay home apart from exercise and essential tasks" [15] . Therefore, according to the level of decreased social contact, a term has been used. Here, we will consider different levels of social isolation, which starts on self-isolation, until a lockdown. However, in order to use a quantity of isolation, we will refer to any isolation by social isolation, and the percentage is the decreased level of social contact.

Regarding what we know so far about COVID-19, Susceptible-Exposed-Infected-Removed (SEIR) models have been used for predictions and analysis [19] [20] [21] 28, 53, 61] . In this model, susceptible individuals become exposed when contacting infected, and after an incubation period, exposed individuals become infected. When infected, individuals may die due to the disease or get cured, changing to removed. For the exposed stage ( E I ), researchers are dealing with uncertainties regarding the start of the disease. Some studies show that the incubation period for COVID-19 is 4-5 days [15, 25] , there are cases of transmission during incubation [62] , and evidence of an early peak of SARS-CoV-2 replication (day 5) than SARS (days 7-10) [54] . One issue related to regular SEIR models is the absence of the asymptomatic state, and the use of one single infected state for hospitalized, non-hospitalized, and hospitalized in Intensive Care Units (UCI). Cellular automata have been used for investigating the spread of contagious diseases, and some of the advantages are the flexibility of local rules and configuration of state transitions [59] . Instead of using deterministic state transitions, the Probabilistic Cellular Automata (PCA) is more suitable for epidemiological studies [1, 2, 45] . Accordingly, it has been used for studying migratory movements on the persistence of contagious disease [5, 16] , the impact of the time delay in the spreading of a disease based on SEIR model [44] , the adaptation of cellular automata for using real population density maps [18] , and disease infection in groups of individuals [37] . Also, for large homogeneous and well-mixed populations, Ordinary Differential Equations (ODE) can be interpreted as a mean-field approximation of the PCA model [4, 32, 42] .

Therefore, here we propose a model based on Probabilistic Cellular Automata and Ordinary Differential Equations, where the states are: Susceptible ( S), Exposed ( E I ), Infected Symptomatic ( I S ), Infected Asymptomatic ( I A ), Infected non-hospitalized ( I S1 ), Infected hospitalized ( I S2 ), Infected in ICU ( I S3 ), and Recovered/Removed ( R ). A S-individual, when exposed to E I , I A , I S , I S1 , I S2 , and I S3 states (not necessarily all and at the same time), has a probability of being infected, changing the state to E I . For E I -individuals, after a period of incubation, may either change to symptomatic ( I S ) or asymptomatic ( I A ). For I A -individuals, there is a probability of getting recovered, becoming R . After a period of symptoms, an I Sindividual has a probability of not needing hospital, being an I S1individual, and an I S1 -individual has a probability of getting recovered, becoming R . If an I S needs a hospital, the state is updated to I S2 (which has its own probabilities of cure and death), and if a hospitalized case needs an ICU or mechanical ventilator, the state is updated to I S3 (also, with its own probabilities of cure and death). All individuals may die due to other causes, with a Sindividual replacing them, to keep the total number of individuals constant.

It is important to emphasize that the main idea in proposing a model with eight states and fifteen parameters (to be presented) is to make the model flexible and easy to input data considering the COVID-19 uncertainties. Also, COVID-19 requires different approaches for the infected state due to different outcomes for the disease case-by-case. Other works also consider parameters for externalities of the epidemiological model (such as the parameters to be presented here which are related to the healthcare system dynamic) [3, 28, 39, 61] .

By using available data of COVID-19, and considering the Brazilian situation on hospital beds, ICUs and mechanical ventilator at the beginning of pandemic outbreak, we propose a SEIR model adjusted to the disease requirements. The objective of this paper is to explore social isolation features such as the start of the social isolation campaign, the percentage of reduced contacts and the duration of the campaign on the spread of the COVID-19. Also, we will explore continuous and periodic isolation operations in order to understand the dynamics of the disease and the impact on the health system for a medium to long-range period.

The paper is organized as follows: the next section contains the proposed SEIR model based on PCA and ODE approaches; the third section has the results of the simulations, and in the fourth section, the results are discussed.

The PCA model consists of a lattice with side n and a total size of n × n = N cells with periodic boundary conditions (forming a toroidal surface to eliminate edge effects), where each cell represents an individual in one of the eight states: S, E i , I A , I S , I S1 , I S2 , I S3 , or R . Cellular automata models usually rely on local interactions to represent an infectious process, and here we expand the neighbourhood, as presented on [42] , where, at the beginning of each time step, individuals make C connections to neighbours inside a Moore radius r. Therefore, the probability of connection between an individual centred on a square of side 2 r + 1 , and another individual at layer i is given by q i,r = 2(r + 1 − i ) /r(r + 1) .

The layer is composed by the 8 l cells at the border of the square of side 2 l + 1 , with l ≤ r. For example, if r = 4 , the probability of interacting to one cell of layers 1, 2, 3, and 4 are q 1 , 4 = 0 . 4 , q 2 , 4 = 0 . 3 , q 3 , 4 = 0 . 2 , and q 4 , 4 = 0 . 1 , respectively, and after randomly choosing the layer, a random cell of the layer is also randomly chosen. The random contact network formed by these connections may have a "high" clustering coefficient and "small" average shortest path length, characteristics of social networks [42] .

Consider the example shown in Fig. 1 with r = 4 . The first layer i = 1 is blue, the second layer i = 2 is grey, the third layer i = 3 is red, and the fourth layer i = 4 is yellow. The white cells are examples of C = 5 connections for a maximum radius of r = 4 . After randomly picking up the layer according to q i,r , a random cell in that layer is randomly chosen.

Due to COVID-19 uncertainties, parameters used for the states transitions are either based on probabilities, or periods before changing the state. Therefore, at each time step a S-individual infectiousness of I S due to limited contact with their close relatives and health care personnel. Although states I IS2 and I IS3 are related to individuals at hospitals beds and ICUs, they are responsible for spreading the disease due to health care overwhelmed situation [8, 56] . We consider a lower infectivity for E I -individuals due to conclusion about the infectiousness starting 12 hours before the onset of symptoms for symptomatic, and 4.6 days when asymptomatic [15] . Also, there are some cases of infection from a recent exposed individual [62] . Therefore, E I -individuals are mildly infective for a period T E I . After the T E I period, P I S is the probability of the E I individual becomes symptomatic I S , and 1 − P I S to become asymptomatic. When asymptomatic, I A -individuals have an average period of T I A days of infectivity, and after this period, they become recovered ( R -individual). For I S -individuals, there is an average delay period of T I S days before going to their final stages of the disease. After this period, a fraction P hosp needs hospitalization, going to I S2 state, and 1 − P hosp becomes I S1 . An I S1 -individual has an average period of T CIS1 days for recovering from the disease. For I S2 -individuals, the average period for a hospital stay is T HIS 2 , and after this period, either I S2 -individual is transferred for an ICU with probability P ICU , die with probability P d , or get cured. Finally, I S3 -individuals have an average hospital stay of T HIS 3 days, and after this period, there is a probability P d−IS3 of dying due to disease complications and 1 − P d−IS3 of getting cured. Finally, all individuals may die due to natural causes with probability P n per time step. Therefore, one time step consists of generating C connections per individual; for S-individuals, the transition S → E is tested according to P i (v i ) ; for the other states the transitions are tested conforming to the respective periods and probabilities; and, after all individuals calculations, their states are updated synchronously at the end of the time step.

An approach for the PCA model can be described in terms of the following ODE if we consider that individuals are homogeneously distributed over space: 

where a is the infection rate constant; b a , b S1 , b S2 , and b S3 are the recovering rate constant for I A , I S1 , I S2 , and I S3 individuals, respectively; d S2 and d S3 are the death rate constant due to disease for I S2 and I S3 individuals; e is the death rate constant due to other causes; α A and α S are the symptomatic and asymptomatic constant rates from exposed individuals; and β S1 , β S2 , and β S3 are constant rates for symptomatic individuals related to non-hospitalization, hospitalization and hospitalized cases requiring ICUs for I S , I S1 , and I S2 , respectively. Fig. 2 contains the diagram of the transitions with the associated variables.

Note that

therefore the total number of individuals is constant, and

We assume that for the period analyzed, the number of deaths is equal to the number of births.

The basic reproduction number ( R 0 ) can be derived from the ODE by using the Next-Generation Matrix (NGM), which is an alternative for models with finite categories of individuals. Therefore, we decompose the Jacobian matrix from the system of ODEs as T + , where T is the transmission part, containing the production of new infected, and is the transition part, containing changes in the states [11] , only for the infected states. With these matrices, we calculate the dominant eigenvalue (spectral radius) ρ of the matrix −T −1 . For the ODE (1) state with E I = I S = I A = I S1 = I S2 = I S3 = 0 , we have: 

Then, the spectral radius ρ(−T −1 ) returns:

with:

Here, we will consider the effective basic reproduction number, R t , which is similar to R 0 by assuming that susceptibility is not complete in the population, and can be used for a partially recovered population. It is given by:

where S(t) is the normalized concentration of susceptible individual at time step t in the population.

The constant parameters from Equation (1) can be obtained from PCA simulation, since the ODE is a mean-field for the PCA model [42] . Therefore,

is the increase per time step of exposed individuals due to contamination process; I S (t) E I → I S is the increase per time step of symptomatic from exposed individuals; I S (t) E I → I A is the increase per time step of asymptomatic from exposed individuals; I S1 (t) I S → I S1 is the increase per time step of nonhospitalized from symptomatic individuals; I S1 (t) I S → I S2 is the increase per time step of hospitalized from symptomatic individuals; I S1 (t) I S → I S1 is the increase per time step of ICU patients from hospitalized individuals; S(t ) I S2 → S is the increase per time step of susceptible individuals due to the death caused by the disease of hospitalized individuals;

time step of susceptible individuals due to the death for other causes of exposed, symptomatic, asymptomatic, non-hospitalized, hospitalized, and ICU patients, respectively. In the next sections, the results are presented.

In this section, we present the results of the simulations with the proposed model. The results are divided into five situations, starting with a comparison between the PCA and ODE approaches, followed by a simulation of the model using real data from Brazil without social distancing control and hospital beds and UCIs unlimited. Then, an analysis of the social isolation features contrasting with the condition with absolutely no isolation between individuals is presented. The fourth situation is an investigation of the minimum percentage of isolation in the population in order to the healthcare system can deal with the ill individuals who require hospital beds, ICUs and mechanical ventilators. Finally, the periodic characteristics of the social isolation campaigns and individuals perceptions of the disease are explored. 3 shows the temporal evolution of PCA ( Fig. 3 a) and ODE ( Fig. 3 b) simulations. For the PCA simulation, the effective basic reproduction number ( R t ) is also plotted. Note the good agreement between the two approaches when we compare the steady states of the simulation (where the ODE parameters are based on PCA simulation). The initial transient is not similar as a result of data from the permanent regime of PCA being used to simulate the ODE, and ODE being a mean-field approach for PCA. Other combinations of parameters have also been tested showing similar results.

Concerning the COVID-19 parameters, we will use the following data: an average period of T E i = 5 days [15] for E I -individuals, probabilities of P I S = 0 . 5 , and P I A = 1 − P I S = 0 . 5 [15, 52] for an E Iindividual becoming symptomatic or asymptomatic, respectively. Once in one of these states, the average period is T I S = 4 days [15] for symptomatic, and T I A = 14 days [15, 40] for asymptomatic. Symptomatic individuals soon need hospital, and they are probably the most tested individuals. Therefore, T I S < T I A , because after this period either I S -individuals need a hospital with probability P hosp = 0 . 2 [27] becoming a I S2 -individual, or they are considered a non-hospitalized case, becoming an I S1 -individual.

I S1 -individual remains in this state for an average period of T CIS1 = 14 days [7] , and then it is cured ( R -individual). A hospitalized I S2 -individual stays at the hospital bed for an average period of T HIS 2 = 8 days [15] . After this period, either there is a probability of P ICU = 0 . 33 [8] to require an ICU treatment, or die due to disease with probability P d = 0 . 0165 [26] , or get cured ( R -individual). Finally, an ICU patient remains in intensive care for an average period of T HIS 3 = 10 days [15] with a mortality of P d−IS3 = 0 . 5 [15] after this period (or 1 − P d−IS3 = 0 . 5 of getting cured). All individuals may die due to other causes with probability P n = 0 . 0 0 0 036 [50] (considering the data from Brazilian population).

The last parameter for the model is k, related to the infectivity of the disease. This parameter will be tuned by using COVID-19 ac- cumulated deaths for the first twenty-five days of the outbreak in Brazil after the first known death, which happened in March 6 th , 2020. In this period, there is no, if any, effect of initial isolation actions (which started somehow at March 21 th , 2020). This parameter has been adjusted manually by considering the accumulated number of death in the period from March 6 th , 2020 to March 30 th , 2020 [31] . Note in Fig. 4 that the adjustment with k = 0 . 15 is conservative when comparing with real data [31] . Table 1 contains a summary of all variables and values considered here. By using these values, the PCA simulation is shown in Fig. 5 for the initial conditions. It is worth noting that the model used here does not consider limits of hospital and ICU/mechanical ventilators beds. Therefore, at the peak of the hospital and ICU beds demand, it would be required 228,243 hospital beds and 93,692 ICU/mechanical ventilators. The total number of deaths until the fade-out of any type of infected individuals in the population is 3,887,423. Also, after the initial peak of infections with R t ≈ 7 , the effective basic reproduction number remains constant with R t ≈ 1 . 2 while the disease is active. These values of R t have been computed by using Eqs. 2, 3 , and 4 .

However, Brazil is far from offering 10 0,0 0 0 ICU/mechanical ventilators for this crisis. Data acquired at the beginning of the pandemic showed that there were around 426,380 hospital beds, 55,100 ICU beds, and 10311 mechanical ventilators available in Brazil [6, 38] . Although the recommended occupation rate of ICU beds in Brazil is 85% [49] , usually the rate is higher than 90%. Therefore, by considering that the medical staff can handle the ICU and mechanical ventilators management during this period, we consider that I S3 may require either ICU or mechanical ventilator. Also, by considering an occupation rate of 80% on hospital beds, 85% on ICU beds, and 50% on mechanical ventilators, there will be available to COVID-19 patients a total of 85,276 hospital beds and 13,421 ICU/mechanical ventilators. From here, we will refer to the management of ICU and mechanical ventilators only by ICU. Therefore, these limits on patients treatment are set on the model on the state transitions I S → I S2 , when symptomatic individuals require a hospital, and I S2 → I S3 , when hospitalized individuals require an ICU. We consider that if there are beds available, the transition is done; otherwise, the individual dies. In this case, transitions I S → S and I S2 → I S are counted. In this restriction of hospital and ICU beds situation, there would be a total of 13,591,960 deaths, with 2,051,427 deaths due to ICU limit, and 11,540,533 due to hospital beds limit. Fig. 6 shows the temporal evolution for this situation.

The social isolation is modelled by considering three parameters: the beginning of the social isolation after the first confirmed death ( Q S ), the duration of the social isolation ( Q P ), and the percentage of social contact reduction ( Q φ ). By using these three features of a social isolation campaign, it is possible to understand their effects on the number of deaths due to healthcare system collapse and the time for the healthcare system to collapse. Therefore, the same simulation configuration presented in the previous section to generate Fig. 6 is used here with the following modifications: After Q S time steps of the first death, the social 1 , 2 , . . . 12 , resulting in 648 cases. If the first death occurs in time step t d , each simulation runs for t d + Q S + 180 time steps. In this period, the system does not achieve a permanent regime, but it is enough for making some considerations.

The first data to be analysed is the number of days for the healthcare system to collapse, when the deaths due to limits of hospital and ICU beds start to happen. This result is shown in Fig. 7 , where the number of days for the healthcare system to collapse is plotted in function of the duration of social isolation ( Q P , in weeks) and the isolation percentage ( Q φ ), for the six values of the isolation campaign start ( Q S ). The figure colour is related to the gradient g(x, y ) of the surface f (x, y ) , given by g(x, y ) = (δ f /δx ) 2 + (δ f /δy ) 2 .

Note that there are only small shifts between the surfaces of different values of Q S . Also, the gradient of the figure is useful to check that combining the effort s on Q P and Q φ is a better option than invest only on either the period or percentage of social isolation. Thus, the isolation campaign should start with at least Q P 6 and Q φ 0 . 6 .

A second analysis consists of the saved lives in function of the social isolation parameters. By considering the same simulations of the previous analysis, the number of saved lives when comparing with simulation without any social isolation ( Fig. 6 ) is shown in Fig. 8 . The period used for comparison is the 180 days follow- ing the beginning of social isolation campaign. The deaths difference due to unavailable hospital and ICU beds is plotted in function of Q P (in weeks), and Q φ , for the six values of Q S . The colour of the figure is related to the gradient of the surface. As in the previous analysis, the isolation campaign should start with at least Q P 6 and Q φ 0 . 6 , considering that any extra effort from this point could have more significant gains than increasing effort s on low values of Q P and Q φ . In this section, the number of hospital and ICU beds, and mechanical ventilator were fixed during the whole simulation. Of course, once the COVID-19 outbreak takes place, government and healthcare association act to increase these numbers. Also, for the cases simulated here, independent of the social isolation features, the COVID-19 is not entirely eliminated of the population. All results showed that the COVID-19 returns to rapidly spread once the social isolation campaign is down. These points are better explored in the next sub-sections and in the supplementary file.

One of the primary responses of healthcare systems around the globe to the pandemics crisis is the purchase of ICU and mechanical ventilators due to the high increase of the demand. Therefore, the minimum daily purchase rate of ICU beds to the healthcare system does not collapse is investigated. We consider a constant daily purchase rate ( τ ) during a simulation from the beginning of social isolation campaign ( Q S ), and the minimum rate is taken. The Fig. 9 contain the results of the simulations for Q S = Note that the beginning of the campaign does not influence the ICU purchase rate. On the other hand, an increase in the social isolation from Q φ = 0 . 4 to Q φ = 0 . 6 may reduce the daily constant purchase rate of ICU by around 100 unities.

A supplementary file is provided with different formats of social isolation campaigns and their effects on the time evolution of the disease.

The objective of this paper was not to propose a breakthrough model which would change the current analysis of COVID-19 whose results should be strictly taken into consideration for public authorities. The objective was to propose a model based on reliable models used for past epidemics to analyse how social isolation can interfere with the COVID-19 dynamics. We used initial data from Brazil, but the same model could be used for any territory or country. It is clear that if the whole population is confined at home, the disease slows down. However, it is not clear that this is the only variable to consider for decisions concerning the COVID-19 future. As discussed in [22] , numerous mathematical models have been suggested since the beginning of the pandemic, and many of them have limitations (usually with an optimistic view) for forecasting the COVID-19 spreading with precision. Not to mention the uncertainty regarding the SARS-Cov-2 or the disease statistics, for instance.

What mathematical models can do is to bring light on COVID-19 issues. Here, we attempted to cover how the social isolation features may change the disease dynamics. Therefore, by proposing a model based on probabilistic cellular automata and ordinary differential equations for a modified Susceptible-Exposed-Infected-Removed case, we added the social isolation variables, as the beginning and the period of the campaign, and the percentage of social contact reduction, for evaluating the impact on the population dynamics. The healthcare system has also been considered, by classifying the infected individuals according to the disease stage and adding healthcare infrastructure data, like hospital beds, intensive care units, and mechanical ventilator numbers. Cellular automata have been used for studying many diseases, such as dengue [10, 34] , Chagas disease [46] , foot and mouth disease in feral pigs [12] , bubonic plague [23] , and hepatitis B [57] , with some of these models considering some heterogeneities in the spatial distributions of the population. Also, in [41] , ordinary differential equations were a good approach for heterogeneous population models based on a wide range of complex random networks. Here, probabilistic cellular automata with a randomly varying neighbourhood could replicate the first days of the COVID-19 spreading. Fig. 10 contains the expansion of the data presented in Fig. 4 (that is, the simulation runs without social isolation), and the divergent trend is visible few days after some social distancing procedures, which started in March 24 th , 2020 (also using data from [31] ).

It is worth mentioning that the flexibility for implementing a model in cellular automata has its drawback. In a regular PC with a 4GHz processor with 16 GB RAM, the simulations of this paper would take two months to be completed. Some efforts have been made to run the model in a Graphics Processing Unit (GPU), which reduced the time to one week and a half.

The first result of the paper is the model itself. Although not usual, the classification of infected states according to the infection stage, with the representation of each stage by a state in the model may help to describe the COVID-19 specification fully. In a study for the Italian case [17] , the infected state has also been split into five states: infected, diagnosed, ailing, recognized and threatened, similar to the states E, I A , I S , I S1 , I S2 , and I S3 presented here. The comparison between the PCA and ODE models showed a good agreement between the approaches.

The second result helps to understand the impact of the social isolation on the Brazilian healthcare system, and it was presented in Fig. 7 . The number of days for the healthcare system to collapse was shown in terms of the isolation percentage and the duration of the campaign. For a short period of social isolation with a low percentage of reduced contacts, the system would collapse in 50 days after the first registered death. The company Inloco , which estimates the social isolation level based on cell phones location data from networks providers [29] , showed that this estimation is around 40% since March, 24 th , 2020.

The estimative was that the hospital and ICU beds will be full in around 80 days after the first death. Since the first registered death occurred in March, 6 th , 2020, the prediction of the paper would be May 25 th , 2020. The peak of hospitals crisis was around June, extending to early August, depending on the region. The long-term results presented here reduced the importance of the day when the isolation campaign starts. Of course that in a short-term, the evolution of infected states is delayed, as also reported in [60] .

Certainly, there would be responses from the Brazilian healthcare system. Therefore, we simulated the minimum purchase rate of intensive care units to consider all patients that could require such treatment. Considering the social isolation features in Brazil for the beginning of the pandemic outbreak, the rate was around 200 ICUs/day ( Q φ = 0 . 4 and Q S = 18 ) according to Fig. 9 . Also, a small increase in social isolation percentage to Q φ = 0 . 5 would reduce the rate to around 150 ICUs/day. It is hard to estimate the exact number of ICU and mechanical ventilators purchased for the healthcare system in Brazil, but a data from Brazilian Ministry Of Health updated in April, 20 th , 2020, showed that 14,100 mechanical ventilators had been purchased, to be fully available for use in July 20 th , 2020. Considering the beginning of the social isolation (in some states) as of March 24 th , 2020, the average purchase rate per day is around 120. This number does not consider private healthcare or purchases made directly by states.

Finally, some social isolation campaign configurations were presented in the supplementary file. The conclusion is that if we consider one data from the epidemic evolution or one data from the healthcare system for the decision about starting and stopping a campaign, there will be a long disease persistence in the population, with oscillations of the infected states over time.Such oscillatory behaviour for public health issues has also been found in a vaccination campaign, for instance [43] . Further waves of COVID-19 propagation have been concerned specialists due to premature relaxation of social isolation campaigns worldwide [58] . Furthermore, campaigns that hardly slow the virus propagation make the disease endemic in population for many years.

Overall, delaying the SARS-Cov-2 spreading seems to be the most effective way to avoid the healthcare system to collapse and save lives. There are some variables to delay the spreading, and here we attempted to investigate some of them. Splitting the campaign decisions into different variables may help to deal with some side effects of social distancing, as reduced sleep quality of selfisolated individuals [56] , extreme social isolation of people living with HIV [30] or living in long-term care facilities [13] , for instance. The next steps of this work could be in the following directions: use a population density map to build the spatial cellular automata; consider a network of cities, to study how the disease spread over the country, and; consider groups of individuals [14, 33, 35, 37] as an approach for the social contact.

Ethics approval was not required for this study.

PHTS is supported by grants #307194/2019-1 and #402874/2016-1 of Conselho Nacional deDesenvolvimento Científico e Tecnológico (CNPq) and grant #2017/12671-8, São Paulo Research Foundation (FAPESP)

The author declares that there are no conflicts of interest regarding the publication of this paper.

On modeling epidemics. including latency, incubation and variable susceptibility

On modeling hepatitis b transmission using cellular automata

Accounting for symptomatic and asymptomatic in a SEIR-type model of COVID-19

Techniques of spatially explicit individual-based models: Construction, simulation, and mean-field analysis

Modeling the impact of immigration and climatic conditions on the epidemic spreading based on cellular automata approach

Cenário dos hospitais no brasil

Gastrointestinal Manifestations of SARS-CoV-2 Infection and Virus Load in Fecal Samples from the Hong Kong Cohort and Systematic Review and Metaanalysis

Clinical Characteristics of 54 medical staff with COVID-19: A retrospective study in a single center in Wuhan

Coronavirus disease (covid-2019) situation reports

Modeling the dynamic transmission of dengue fever: Investigating disease persistence

The construction of nextgeneration matrices for compartmental epidemic models

Simulating the spatial dynamics of foot and mouth disease outbreaks in feral pigs and livestock in Queensland, Australia, using a susceptible-infected-recovered cellular automata model

Breaking Social Isolation Amidst COVID-19: A Viewpoint on Improving Access to Technology in Long-Term Care Facilities

Modelling disease outbreaks in realistic urban social networks

Impact of non-pharmaceutical interventions (NPIs) to reduce COVID-19 mortality and healthcare demand

The impact of imported cases on the persistence of contagious diseases

Modelling the COVID-19 epidemic and implementation of population-wide interventions in Italy

Epidemiological modeling with a population density map-based cellular automata simulation system

The effectiveness of the quarantine of Wuhan city against the Corona Virus Disease 2019 (COVID-19): well-mixed SEIR model analysis

Spatial-temporal distribution of COVID-19 in China and its prediction: A data-driven modeling analysis

A Simulation on Potential Secondary Spread of Novel Coronavirus in an Exported Country Using a Stochastic Epidemic SEIR Model

Predictive Mathematical Models of the COVID-19 Pandemic: Underlying Principles and Value of Projections

Bubonic plague: A metapopulation model of a zoonosis

An Epidemiological Study on COVID-19: A Rapidly Spreading Disease

Epidemiologic characteristics of early cases with 2019 novel coronavirus (2019-nCoV) disease in Korea

Correcting under-reported COVID-19 case numbers: estimating the true scale of the pandemic

Articles First-wave COVID-19 transmissibility and severity in China outside Hubei after control measures, and secondwave scenario planning : a modelling impact assessment

COVID-19 pandemics modeling with modified determinist SEIR, social distancing, and age stratification. The effect of vertical confinement and release in Brazil

Physical Distancing in COVID-19 May Exacerbate Experiences of Social Isolation among People Living with HIV

Ministério da saúde -covid-19

Big cities : Shelters for contagious diseases

Social contacts and mixing patterns relevant to the spread of infectious diseases

Dengue fever spreading based on probabilistic cellular automata with two lattices

Dynamics of dengue epidemics in urban contexts

Projeções e estimativas da população do brasil e das unidades da federação

Disease spreading on populations structured by groups

Recursos físicos -hospitalar -leitos de internação -brasil, tabnet, datasus

Why is it difficult to accurately predict the COVID-19 epidemic?

Effect of delay in diagnosis on transmission of COVID-19

Disease spreading in complex networks: A numerical study with Principal Component Analysis

On the basic reproduction number and the topological properties of the contact network: An epidemiological study in mainly locally connected cellular automata

A vaccination game based on public health actions and personal decisions

Impact of time delay on the dynamics of SEIR epidemic model using cellular automata

A cellular automaton model for the effects of population movement and vaccination on epidemic propagation

A cellular automata model for Chagas disease

COVID-19: Real-time dissemination of scientific information to fight a public health emergency of international concern

Covid-19: How doctors and healthcare systems are tackling coronavirus worldwide

Tábuas completas de mortalidade

Indications for healthcare surge capacity in European countries facing an exponential increase in coronavirus disease (COVID-19) cases

Estimates of the severity of COVID-19 disease

When will the battle against novel coronavirus end in Wuhan: A SEIR modeling analysis

Virological assessment of hospitalized patients with COVID-2019

The Effects of Social Support on Sleep Quality of Medical Staff Treating Patients with Coronavirus Disease

A probability cellular automaton model for hepatitis B viral infections

Beware of the second wave of COVID-19

Cellular automaton modeling of epidemics

Clinical characteristics and imaging manifestations of the 2019 novel coronavirus disease (COVID-19):A multi-center study in Wenzhou city

Modified SEIR and AI prediction of the epidemics trend of COVID-19 in China under public health interventions

A Familial Cluster of Infection Associated With the 2019 Novel Coronavirus Indicating Possible Person-to-Person Transmission During the Incubation Period

Supplementary material associated with this article can be found, in the online version, at 10.1016/j.cmpb.2020.105832