Creating and applying SIR modified compartmental model for calculation of COVID-19 lockdown efficiency

Highlights • SIR compartmental model of COVID-19 spread may help to evaluate efficiency of lockdown measures.• Epidemiological data collected in fifteen european countries, are taken as input parameters.• Herd immunity level and the time of its formation are considered as factors indicating efficiency of lockdowns.• Lockdown and no-lockdown mode of containment lead to roughly similar results. Lockdowns do not stop COVID-19 spread.• Rationale for lockdowns is avoidance of healthcare system overburdening.• Reliance on lockdowns as the only administrative measures that may contain a pandemic, is a misstep that should be avoided in the future.


Introduction
Efficiency and necessity of lockdown measures implemented on a scale of the whole world, cause much controversy [1][2][3] . In the current paper, we analyse whether total lockdowns are helpful in stopping spread of Coronavirus disease -2019 (COVID- 19) and future similar global diseases, by means of investigating herd immunity formation to Severe acute respiratory syndrome-related coronavirus-2, the viral causative agent of COVID-19 disease (SARS-CoV-2). In the current absence of a vaccine, herd immunity remains the only way to stabilise human population reaction to the novel viral pathogen.
It was repeatedly emphasised that creating of so-called herd (population, block, natural) immunity is important for slowing down the rate of COVID-19 spread in population and, actually, for stopping the pandemic [4][5][6][7][8][9][10] . However, as of 25 June 2020, not many official investigations of SARS-CoV-2 dissemination in the whole ecosystems and populations (e.g. random mass testing, representative sample screening, closed ecosystem studies, etc.) have means of social control and political compelling/blackmailing in many parts of the world. Excepting a very limited number of countries, the almost global and universal governmental response to COVID-19 pandemic evidently put the world on the verge of the Orwell's 1984 scenario. Now, in September 2020, when media and many politicians do not stop to speculate about second and consecutive waves of SARS-CoV-2 pandemic, lockdown measure are looming on the horizon again. In UK, in the small city of Aberdeen in Scotland, exceeding the acceptable SARS-CoV-2 infection limits defined by the government merely by fifty humans, plunged the city in the full lockdown of indefinite length again. The Brazilian population vast in size continues to suffer from incessant lockdowns. In USA, each state administration regards COVID-19 differently and the country remains severely split in terms of the SARS-CoV-2 containment strategy. In Russia, all governmental media do not cease to intimidate people by a possibility that soon everybody will be locked again and every citizen will pay for his/her relaxed summer rest (basically, normal mode of life without coercive deprivation of personal freedoms). Is the time not come when we should critically evaluate with the help of mathematical models backed by experimental epidemiological information collected thus far, to which extent the lockdown measures and strategies should be used again, if used at all?
Currently, a lot of mathematical models of SARS-CoV-2 spread have been elaborated, both theoretical [16,17] and computational, many of them already described in Special Issues of Chaos, Solitons and Fractals journal "Modeling and forecasting of epidemic spreading: the case of Covid-19 and beyond" [18,19] . By dint of our model, we hope to contribute to our mutual understanding of the lockdown measures efficacy.

Methodology
We propose a Susceptible-Infected-Recovered (SIR) modified compartmental model. Using epidemiological data as input parameters for the model, we calculate times and levels of herd immunity formation for different modes of containment (lockdown and no-lockdown). Finally, we evaluate the efficiency of lockdown measures.

Privately collected primary epidemiological data
COVID-19 primary statistical data sets were collected since 2 March to 21 April 2020. In May-August 2020, additional data have been collected that changed our understanding of SARS-CoV-2 spread substantially. The Acknowledgement section contain a number of names of people who assisted in collecting the data and expressed their explicit consent for their contribution to be noted.
The primary information and its sources are summarised in Table 1 for the initial period (March-April 2020).

Modelling the COVID-19 spread in population
We assume that the dissemination of COVID-19 may be explained by a continuous-time Markov process model [35,36] . This process can be generally described by a Susceptible-Infected-Recovered (SIR) compartmental model [37,38] . Some authors used Susceptible-Exposed-Infected-Recovered (SEIR) model for their calculations [39][40][41] . We chose to use a common SIR compartmental model altered by introducing several significant modifications to the corresponding Markov process structure.
Regarding lockdown measures imposed on citizens, we used two models: a) common; b) Swedish. In the former, full lockdown measures were implemented during the surge of COVID-19 disease (only Home -Pharmacy -Grocery Store -Walk with an Animal (e.g. dog) -Hospital possibilities), in the latter no lockdowns were in force. Most of countries followed the lockdown way. Very few countries implemented no lockdowns, e.g. Sweden, Iceland, Belarus, Japan, Taiwan, Republic of Korea. The Swedish model is more complex, and the corresponding Markov process includes 13 knots (in our modelling). The knots (states) of the Markov processes of people relocation in a community, are shown in Fig. 1 . Fig. 2 demonstrates our SIR modified compartmental model as a continuous-time Markov process modelling the infection transmission, with the knots (states) corresponding to different possible situations of a human in regard to the disease. We have no information about the total number of COVID-19 tests made on a scale of a whole country, if the information is not disclosed publicly. * * As of 21 March 2020. * * * In Lombardy only. * * * * In-patient (hospital only) PHE statistics, as of 21 April 2020.

Table. 2
Estimation of herd immunity formation in Europe (95% confidence interval) during the "first wave" of SARS-CoV-2 according to our SIR modified modelling with the privately collected input data on different European countries, for different possible basic reproduction numbers r 0 .

Mode of containment
Full-lockdown No-lockdown Basic reproduction number r 0 r 0 = 1.6 r 0 = 5.6 r 0 = 1.6 Negative time values stand for the time before the beginning of the infection process in a given community. * Scenarios:. a world disease peak ("first wave" of SARS-CoV-2) has not yet been passed on 11 March 2020 (announcement of the pandemic by World Health Organisation (WHO));. b there are two peaks a and c ; a current situation is their overlapping;. c disease peak has been already passed unnoticed before the pandemic was announced by WHO.

Distribution and relocation of population in space and time in two different modes of containment
Let us consider a human who may be in any of the five places in a full lockdown mode of life. Then the probabilities of his/her stay in these five places are p 0 (home), p 1 (pharmacy), p 2 (grocery store), p 3 (street saunter with an animal, most commonly a dog), p 4 (hospital / any other medical institution, in inpatient / outpatient modes of treatment). It is obvious that First, let us write a system of Kolmogorov linear differential Eqs.
(2)-(6) for the movement of a human in Model 2 (in a country with a full lockdown mode). For the considerations of simplicity, we assume that a hospitalisation may be made only from his/her home, and intensities are not functions of t : (3) For the second model (Swedish), we have fourteen Kolmogorov equations. Taking into account that Markov process in the Swedish model is ergodic, in comparison with the full lockdown Markov process, we receive the following system of Kolmogorov equations: Typical values of λ ij and μ ij are taken from the similar models of human relocation, as described in the works [42][43][44][45] .
To estimate the number of the immune people and, therefore, the level of herd immunity created thus far, we have to observe how Kolmogorov equation systems 2-6 and 7, and the corresponding master equations behave (how the solutions p i ( t ) change) on varying the epidemiologic and demographic parameters. The solutions, i.e. probability functions p i ( t ) are complex functions of t , and they have these parameters as variables: where r 0 is basic reproduction number, c drop contagiousness coefficient of droplet transmission, c surf contagiousness coefficient of touching surface transmission, N the total number of humans in a community, N imm the number of people in the community who already have the immunity to SARS-CoV-2.
Using coefficients of transmission is a rather new concept. In fact, these coefficients represent probabilities of a nonimmune person of becoming infected through a direct contact with a person already infected or a contaminated surface: c =

N in f ect ed a f t er contact N tot al cont acts
, and, therefore, 0 ≤ c drop. + c surf. ≤ 1. In a more broad sense, 0 ≤ c drop. + c surf. ≤ 2, but the value 2 for the sum is not achievable even for the most contagious diseases known to humanity thus far (e.g. varicella for droplet transmission or Ebola haemorrhagic fever for body liquids transmission). Basic reproduction number r 0 is estimated differently by different research groups. The difference is tremendous, from 1.6 [36] or 2.6 [46] to 5.6 as a median with 6.7 as the maximal value [47] and even almost 15 for the virus spread estimations on Diamond Princess cruise ship [20] . In its report of 17 and 24 April 2020, Robert Koch Institute made evaluations of r 0 as low as 0.5-0.6 [29] . Such enormous difference in evaluating basic reproduction number may result from the fact that various research groups studied different samplings incomparable with each other in terms of closeness and rate of human contacts. We will use r 0 = 1.6 and 5.6 as reference points.
The sum c drop + c surf was initially estimated in 0.3-0.4 range (nearly 3 or 4 of 10 persons directly contacting an infected individ-ual of surfaces with a full virus titre, will be infected by SARS-CoV-2 virus) [48][49][50] , where c drop may be 0.1-0.2 [50] . Values of N for different communities are taken from various demographic sources publicly available over the Internet.
is required to be assessed in our model.

The model of the virus transmission
To perform estimations of N imm ( t ), an additional Kolmogorov equation system should be included in our analysis. Within our SIR modified compartmental model, this system of linear differential equations describes intensities and directions of flows in the transmission model ( Fig. 2 ). Analysing the behaviour of Markov process that describes the virus spread, along with the behaviour of systems 2-6 and 7 of space-time population distribution, enabled us to evaluate the level and time of herd immunity formation in the countries studied.
A non-infected person may contact with: a) an ambiguous person (he/she may be either infected or non-infected) (knot 1); b) a proven SARS-CoV-2 carrier (knot 2); c) a contaminated surface (knot 3); d) an ambiguous surface (a virus amount may be present or may not) (knot 4). The knots 5-7 correspond to the virus carriers proven by an RT-PCR or RT-LAMP test: asymptomatic OR those in the incubation period (knot 5); mild symptomatic (knot 6) and severe symptomatic with atypical pneumonia clinical picture (knot 7). Finally, the outcomes are either decease (knot 8) or recovery (knots 9-10). amongst the recovered patients, there may be the contagious ones (knot 9) and completely recovered ones without active virions in their bodies, and therefore, without any SARS-CoV-2 contagiousness (knot 10). The Markov process is non-ergodic, non-stationary. We are interested in the characteristics of the final state of this process. Potentially, with relatively large time values, the Markov process described will reach a stationary state of the maximal herd immunity formation (one global spike of disease) or several quasi-stationary states of local herd immunity plateaus formation (two or more local spikes of disease), with the last plateau being the global maximum of herd immunity achievable in population. The stationary state of the Markov process discussed will not change over time any longer; the quasistationary states will not change over some time. These states correspond with herd immunity global and local plateaus. We shall find the times and values of these plateaus by finding the inflection points of herd immunity dependence on time, that coincide with extremum points of herd immunity first time derivative.
The system of Kolmogorov equations for this SIR modified disease transmission compartmental model is the following (9)(10)(11)(12)(13)(14)(15)(16)(17)(18)(19): (19) where probability changes (first derivatives) in the left parts of the equations are calculated by means of separate probabilities p i and current intensities λ ij. A separate probability p i stands for a description that in moment t a human being will be in state S i (one of the described states). A major assumption of the model that oversimplifies it is that a said human cannot be in two states simultaneously. Plus signs are for straight direction (from initial events in the disease to further events, i.e. classical progression of the infection) (black arrows in Fig. 2 ) and minus signs for backward events that return a human to previous states (blue arrows in Fig. 2 ). The exact structure of possible connections and, therefore, currents, were taken from the book of Karin VanMeter and Robert Hubert [35] . Besides, the condition (1) is met for this system, as m m = 8 where m is the number of Markov process knots, in our case 11, and we consider that the number of deaths is small.
Several obvious border conditions and definitions may be further explicated, such as p 0 (0) = N (0) / N < 1; the inequality assumes that at the initial moment of time there is a portion of persons that already have the immunity to SARS-CoV-2.
we may take T → ∞ for the simplicity), is the final herd immunity level; while p 10 ( t ) + p 9 ( t ) may be not equal to N immune ( t ) / N for any moment of time. For any t the intermediary herd immunity level may be measured as Likewise, p 8 ( T ) = N d ( T ) / N, is total population fatality rate (TPFR). For any t we use instead. It is population fatality rate (PFR) that shows the percentage of the deceased in the total population.
for any t , is population infection rate (PIR); and is the infection fatality rate (IFR). We composed the system (9-19) just for one person. Solving the system numerically for each person for each state in (2-6) and 7 with a subsequent averaging, and observing how the solutions would change on varying epidemiological parameters, is not a best algorithm. In a community there are N people, and this number is constant (the number of deaths is small in regard to the total population size). Therefore, we have a multinomial distribution. We may assume that the time of a contact with a carrier is negligibly small in comparison with other times (e.g. the time of selfisolated life or a stay in a hospital). Then, according to Sanov theorem about large deviations for a multinomial distribution [51] further explicated by Borovkov [52] , we receive the equation: Indeed, (26) can be obtained for our case in such a way. Infection process development may be described as a continuous Ehrenfests' chain [44,53,54] . Let ξ ( t ) stand for these Ehrenfest chain within our virus spread system. ν m and ν m + 1 are parts (proportions) of Markov chains η i that in time moment t are in states m and m + 1, i.e. they can be defined as Ehrenfest frequencies: Then, as explained by A. Gasnikov in details in his work [54] , let us express the probability of these frequencies being probabilities p 1 and p 2 (later we may transfer similarly to all our probabilities p m ): Indeed, Ehrenfest chain ξ ( t ) with N + 1 knots may be expressed via N independent or dependant Markov chains η i with only one border condition P ( η 0 (t) = 0 ) : (in conditional designation) Hence we have Switching in (29) to the limit at t → ∞ (we assume that the virus spread may be as prolonged as we prefer and not restricted by time), we receive After that, we proceed to the limit at N → ∞ (the human population size where SARS-CoV-2 is spreading may be as large as we prefer). Using Stirling formula for factorial approximation n ! = n n e −n √ It is obvious that the right part may be re-written using logarithmic functions of probabilities ψ( p i ) = -p 1 ln p 1 -p 2 ln p 2 and exponent, and then we have: We investigate our system of SARS-CoV-2 dissemination for stability -precisely what we do with studying the behaviour of its Lyapunov functions. Here the stability may be defined as Sanov [51] , Borovkov [52] and Gasnikov [54] define it, viz. a stable state of virus spread system ( p 1 , p 2, …, p m ) in whose little neighbourhood stationary measure (probability of frequencies ν m receiving values m , m + 1, m + 2, …) is being concentrated, independently of t . It is simpler to consider stationary states with t → ∞ , but non-stationary states may be obtained as easily. Maximising It is not difficult to see that we received exactly our system of Eqs. (9 -19 ) with Lyapunov functions ψ( p i ), Q.E.D. The same procedure may be carried out first for considering limit at N → ∞ and only after that at t → ∞ . In such a case, the completely similar result may be received using the theorem formulated by Thomas G. Kurtz [53] . Now, having assured ourselves that system (36) is in fact system (9)(10)(11)(12)(13)(14)(15)(16)(17)(18)(19), we can analyse the behaviour of our multinomial distribution (26) without solving system (9)(10)(11)(12)(13)(14)(15)(16)(17)(18)(19) for each state in (2-6) and 7. We can vary r 0, c drop + c surf and compare the outcome. A proprietary C# algorithm was written for estimating levels and times of herd immunity formation in the countries studies. We observed the stability of the multinomial distribution of virus spread in the population (26) by analysing the behaviour of Lyapunov functions of master equation corresponding to the process (9)(10)(11)(12)(13)(14)(15)(16)(17)(18)(19) in the non-negative orthant R n + (all probabilities and herd immunity levels are non-negative real numbers from 0 to 1; times are non-negative real numbers). Finding the attractors (if any) by means of analysing Lyapunov functions behaviour in the orthant R n + would allow us to determine the most probable values of the projected herd immunity and the times of its formation.

Software packages
Microsoft Visual Studio 2010, OriginLab OriginPro 8.1 and PTC Mathcad 6.0 were used.

Time of preparing the paper
The paper was prepared in its initial state in April 2020 and was re-considered several times in May, June and August with serious corrections of the model proposed. The final revision was completed by 24 September 2020.

Change of population infection rate over time
In addition to the mean population infection rate (PIR) level, its change over time is of importance. In Fig. 3 , the time dependence of instantaneous PIR is shown for United Kingdom, Russia, the Netherlands (lockdowns) and Sweden (no lockdown). Instantaneous PIR is calculated as P I R inst. = N new pos. cases N new tests made , i.e. the ratio of daily increase of positive cases to daily increase of the tests performed.
We observed mixed dependency on time here. By 21 April 2020, symbate dynamics of instantaneous PIR and daily increment of positive cases is observed only for Russia. For UK, there is a maximum of instantaneous PIR that does not coincide with the maximum of confirmed positive cases. For the Netherlands, the approximation is very close to a constant line. For Sweden, we may localise a decay, low plateau and growth. On the whole, for the four countries concerned, instantaneous PIR as a function of time does not correlate with confirmed cases dynamics.
Asymbate dynamics of instantaneous PIR and confirmed positive cases at the initial stages of the pandemic may be generally explained by three factors. First, the samplings in at least one case (confirmed positive cases measurements or mass screening procedures) may be not fully representative on a scale of a population. Second, the change in new tests performed on a daily basis may grow with different speed than the infection dissemination, e.g. faster constant growth of number of tests made than the disease spread may result in an instantaneous PIR decay on growing confirmed cases dynamics. Third, several spikes of disease may be present of which we are currently detecting only one. Along with relatively low values of PIRs, it enables us to suggest that three scenarios for European countries studied are possible, regarding real (not observed) disease spread extremums: a disease peak has not yet been passed on 21 April 2020; b there may be several disease peaks, and different European countries may be situated in different time places regarding the peaks (e.g. between them); a current situation may be an overlapping of different peaks; c disease peak has been already passed unnoticed by the world community, by the moment of the pandemic announcement (early March 2020).
Taking into account these three scenarios, creating herd immunity may become an even more important factor for population adaptation to COVID-19. It was so in April and continues to be in August 2020.

Herd immunity formation under lockdowns and without lockdowns
To estimate times and values of maximums of herd immunity formation in the countries studied, PC numerical solving the Kolmogorov equation systems described in Materials and Methods section with modelling the Lyapunov functions behaviour for the corresponding master equations (finding global attractors) has been performed for the three scenarios described above, for fulllockdown and no-lockdown modes of SARS-CoV-2 containment.
On maximising herd immunity first derivative we receive the inflection points for herd immunity curves. The schematics for Sweden is shown in Fig. 4 . For the interval 0-120 days from Day Zero, confirmed positive cases dynamics (official WHO data) may be approximated by an extreme function  [53][54][55][56] . The extreme function provided above was chosen of the list of the most widely used functions by criterion of achieving the best approximation results (the lowest χ dim.

and highest R adj.
2 ), where the parameters were calculated as follows: m is the numbers of points; n -1 is the number of degrees of freedom; x i are experimental epidemiological data; C i are approximated values; and σ 2 is variance.
The adjusted coefficient of determination where R is Pearson correlation coefficient between experimental data and the approximation and k is the number of explanatory  terms (descriptors). Of course, another approximation could be used instead of the proposed one.
Using private data on Sweden PIR, taking into account its nolockdown mode of containment, r 0 = 5.6 and maximising herd immunity first time derivative by Lyapunov function analysis in the positive quadrant, give us a double exponential function dh (t) dt (green curve in Fig. 4 ): The coefficients: D = 3417.46 days -1 t z = 62 days p = 1.6 q = 150.0 days u = 1 v = 50.50 days Finally, integrating dh (t) dt over t , we receive a sigmoidal function of h ( t ) (black curve in the inset in Fig. 4 ), with its maximum value of approximately 540,0 0 0 people (based on the epidemiological data collected for Sweden in June 2020). That means the overall herd immunity rate. Taking into consideration the population of Sweden (currently 10,319,600 people), it gave 20.80 per cent of population in April and 5.20 per cent in June 2020. Thus, we see a dramatic correction of the model to be discussed below.

Iteration-based correcting and refining the model
SARS-CoV-2 is an evolving process that change our understanding on a constant basis. Obviously, SARS-CoV-2 has not yet reached limits of its dissemination. In our truly global world without any fixed cultural boundaries that might help to slow down tourism and migrant currents, spreading a new virus will continue till a stable co-existence with homo sapiens population has been attained [57][58][59][60][61] . Recently, Oleg Donskikh [62] and Alexandre Gnes [63] demonstrated the influence of cultural implications on epidemiological situation. As Wolfgang Sassin evidently showed in his modelling human evolutionary processes, in the global open world with population more than 1.2-1.5 billion people, even a local viral pathogen outbreak with 100-150 humans infected may hit the whole world in two-three months and the process of its dissemination will not be stopped by any restrictive means until a new evolutionary niche with co-existence of homo and the virus has been created [59,60,64,65] . Regarding SARS-CoV-2, this niche is not created yet. New data on SARS-CoV-2 situation are appearing incessantly that help to improve the model by the possibility to input the most novel and relevant epidemiological information.
Four major iterations of model correction/refinement have been made (mid-March prediction of April, mid-April prediction of May, mid-May prediction of June, and mid-June prediction of further   SARS-CoV-2 spread). It is very important that the newest epidemiological information on the pandemic changed the model dramatically since March 2020, when Europe had entered the first SARS-CoV-2 wave. Additionally, different type of information has emerged since that time, viz. antibody testing programmes have been launched in many European countries, that are auxiliary sources of our knowledge about the virus spread. Fig. 5 demonstrates that now, much later after SARS-CoV-2 initial outbreak in China, PIR diminished significantly, reaching a lower plateau independently of the country or mode of containment (full lockdown / partial lockdown / no lockdown) concerned. Of six states the data on which are presented in Fig. 5 , UK and Russia chose full lockdown mode, the Netherlands partial lockdown while Sweden, Iceland and Belarus no lockdown at all. Surprisingly, since May 2020, a clear trend of PIR decrease was observed everywhere. In most of countries, PIR can be estimated now as some 2-4 per cent, in Sweden and Belarus 6-7 per cent, whereas initial March and April numbers were substantially higher. It may be possibly explained by improving representativeness of the virus (RT-PCR) and antibody testing programmes. At the beginning of the pandemic, only people who were seriously ill with COVID-19 got into the focus of the tests, whether official, government-sponsored or commercial. This fact unfolds a likely origin of highly overstated epidemiological figures for UK and Belgium in Match-April (PIR of nearly 35-40 per cent). Later milder symptomatic patients were included in the tests too. But only since May and especially in sum-mer 2020, a broad and much more representative cohort of European population commenced to be tested [66][67][68][69][70] . Using the May and June data as input parameters for the model enabled us to augment its quality and predicting force. June iteration based on June European epidemiological situation resulted in the herd immunity level prediction for August approximately four times less than the April prediction.

Statistical treatment
Different predicting force of the model has been observed for different types of ambience: closed communities, semi-open premises and open space ( Fig. 6 ). In Fig. 6 , the June modelling iteration is provided that was performed at the time when the paper was returned for revision. It was based on the data collected by us during surveillance procedures in UK, Russia, the Netherlands and Sweden. As modelling error estimations show ( χ dim.
2 and R adj.
2 ), the model allowed to receive the best coincidence of predicted herd immunity with the experimental data for semi-open premises (multi-storey living houses, condominiums, fenced living districts, large offices, trading floors, warehouse premises, etc.), where lowest χ dim. 2 and highest R adj. 2 were achieved within apartments and small houses) gave less exact results ( Fig. 6 , panels C and D). A biphasic structure of the virus spread predicted by the model was in a satisfactory agreement with the experimental epidemiological data. Modelling spreading SARS-CoV-2 in closed communities (family environments within apartments and small houses) gave less exact results ( Fig. 6 , panels A and B). Predicting the virus dissemination in open spaces was the least precise ( Fig. 6 , panels A and B).
In Fig. 7 , the comparison of the model precision is given for the four main iterations of correcting the model. Precision coefficient is calculated as the mean of precision factors for each point i of n : where C x i is the largest by module of { C mod i; C real i }; C mod i is a modelled value in point i ; and C real i is an observed value in this point. In Fig. 7 , one can see that the model changed essentially since the beginning of the pandemic. In March 2020, PIR = 19.19% was taken as the maximum of herd immunity that may be achieved without any lockdown and restrictive measures. In August 2020, it is already becoming clear that value is the maximum only for closed communities ( Fig. 6 , panel A). The fact is not very surprising, as such places with interlinked ventilation as Diamond Princess cruise ship may be considered a completely closed ecosystem [20][21][22][23][71][72][73][74][75] . An opener ambience leads to much less PIRs ( Fig. 6 , panels C and E). Current improving the model does not mean, however, that we achieved a desired level of precision. A possible surge of SARS-CoV-2 infected humans' figures in autumn 2020 and winter 2020-2021 may introduce further significant corrections.

Conclusions
The analysis of both epidemiological data and simulation results indicates that the initially anticipated herd immunity level for SARS-CoV-2 of nearly two thirds of a population or even higher, is hardly ever achievable. The real herd immunity for the current virus is three-six times as less. Therefore, COVID-19 contagiousness is not so high as it was initially thought in January-March 2020. Almost universal and worldwide implementation of lockdown measures and complete switching off the economies, as it has been done in Germany, France, Russia or UK, may be reckless. Despite some governmental and administrative assumptions that only strict quarantine might lead to diminishing SARS-CoV-2 spread, our study does not confirm it. Neither modelling, nor PIR statistical data on the European countries collected and studied by us so far, may corroborate that full-lockdown modes are any better than the Swedish no-lockdown mode in terms of the virus dissemination ( Fig. 8 , Table 2). Fig. 8 demonstrates the June prediction iteration with 95% and 99% CI. According to the model, only at the very beginning of the epidemic, the lockdown mode may demonstrate better results in containment of the virus than no-lockdown one (pink line is lower than blue line). After 7-8 weeks of the virus spread, both modes start to give similar outcomes. This simulation is corroborated sufficiently by the experimental data obtained. Countries with full lockdowns and no lockdowns reached very similar results in PIR and herd immunity, viz. 3-5 per cent predicted by the model after the June iteration ( Fig. 8 ) and 3-7 per cent of the total population observed in reality ( Fig. 5 ). The basic reproduction numbers are also similar [75,76] .
Why so? First, during lockdowns people were still allowed to go outside, e.g. to grocery stores and saunter with pets. Meetings that caused reciprocal infection in elevators, on streets and in close premises, were imminent. Second, ventilation system in large multi-storey buildings disseminated the virus further. Third, public transport was still on, where cross-infection possibility was high.
Was any sense in lockdowns at all? We think yes. They helped to avoid healthcare system collapse in many countries. However, their near and distant harm to health of different groups of population, economies, business and world supply chains is still to be assessed in the future.
We suppose that a real hazard of COVID-19 lockdowns is associated with the common governmental belief backed by several World Health Organisation's announcements that it is the lockdowns that saved humanity from excessive mortality connected with COVID-19. Our research proves that it is a very dangerous misbelief with far-reaching consequences. It is SARS-CoV-2 relatively low contagiousness and case fatality rate that led to avoidance of millions of deaths, not lockdowns. We agree with the Editors that the health consequences of the pandemic are devastating [77] . However, non-evidence-based reliance of governments just on total lockdown as a universal measure of the pandemic containment may be much more devastating in the future, in case of possible consecutive waves of SARS-CoV-2 or any other viral pathogen. Governments should not use simplest lockdown ways of "containment" the virus with an obvious threat to the democracy and human rights in the whole world instead of elaborating complex and effective epidemiological and social strategies.
Why may the belief in lockdowns as the only means of pandemic containment be harmful? In our truly global world with global tourism and transport routes, an emergence and possible immediate worldwide spread of a new viral pathogen, much more dangerous than SARS-CoV-2, may be only a matter of time. In case this virus may have a mortality of Ebola, contagiousness of varicella and worldwide spread, governmental and administrative reliance just on lockdowns may cost hundred millions of human lives.
In the end, lockdowns may be ineffective because they do not and potentially cannot stop a respiratory virus spread [78,79] . Instead, effective healthcare and public policy containment measures must be depended upon, such as strictest epidemiological surveillance in airports, international railways and bus stations; managing tourist routes at the very beginning of early epidemic warnings; applying antiviral cleansing and sanitation of publicly accessible territories, especially in urban environment; habitual wearing individual protective units in places of mass congestion of people and during seasonal respiratory disease surges (e.g., according to Japanese positive experience); and preparedness of medical care systems for global pandemics and other "Black Swan" emergencies. All these measures can be much more effective than lockdowns.

Funding
This work has been funded by Government program of basic research in Koltzov Institute of Developmental Biology, Russian Academy of Sciences in 2020, No. 0108-2019-0 0 02.

Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper..