Next Article in Journal
Haze Occurrence Caused by High Gas-to-Particle Conversion in Moisture Air under Low Pollutant Emission in a Megacity of China
Previous Article in Journal
Association between Body Mass Index with Sugar-Sweetened and Dairy Beverages Consumption in Children from the Mexico–USA Border
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

How Seasonality and Control Measures Jointly Determine the Multistage Waves of the COVID-19 Epidemic: A Modelling Study and Implications

1
State Key Laboratory of Organic Geochemistry, Guangzhou Institute of Geochemistry, Chinese Academy of Sciences, Guangzhou 510640, China
2
CAS Center for Excellence in Deep Earth Science, Guangzhou 510640, China
3
University of Chinese Academy of Sciences, Beijing 100049, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2022, 19(11), 6404; https://doi.org/10.3390/ijerph19116404
Submission received: 26 March 2022 / Revised: 15 May 2022 / Accepted: 20 May 2022 / Published: 25 May 2022

Abstract

:
The current novel Coronavirus Disease 2019 (COVID-19) is a multistage epidemic consisting of multiple rounds of alternating outbreak and containment periods that cannot be modeled with a conventional single-stage Suspected-Exposed-Infectious-Removed (SEIR) model. Seasonality and control measures could be the two most important driving factors of the multistage epidemic. Our goal is to formulate and incorporate the influences of seasonality and control measures into an epidemic model and interpret how these two factors interact to shape the multistage epidemic curves. New confirmed cases will be collected daily from seven Northern Hemisphere countries and five Southern Hemisphere countries from March 2020 to March 2021 to fit and validate the modified model. Results show that COVID-19 is a seasonal epidemic and that epidemic curves can be clearly distinguished in the two hemispheres. Different levels of control measures between different countries during different seasonal periods have different influences on epidemic transmission. Seasonality alone cannot cause the baseline reproduction number R0 to fall below one and control measures must be taken. A superposition of a high level of seasonality and a low level of control measures can lead to a dramatically rapid increase in reported cases.

1. Introduction

The current novel Coronavirus Disease 2019 (COVID-19) is a multistage epidemic consisting of multiple rounds of alternating outbreak and containment phases. In a conventional single-stage epidemic, the number of daily new confirmed cases increases over time, gradually decreases after reaching its maximum, and eventually declines to zero despite sporadic relapses [1,2]. However, in a multistage epidemic, the number increases, reaches the maximum, and declines but does not return to zero. Instead, a new round of outbreak occurs, repeating the above steps [3]. By 20 March 2022, 470,1098,898 confirmed cases and 6,098,267 deaths have been reported, and the number continues to rise steadily [4]. So far, most countries have experienced multiple rounds of large-scale outbreaks, but the driving factors of the multistage epidemic and its specific influences remain to be elucidated.
Control measures may be an important factor of the multistage epidemic [5,6], which interrupt transmission chains by reducing contact between people [7], and different levels of control measures may have a different impact on the transmission of COVID-19 [8,9,10]. Mild control measures include wearing masks, maintaining social distance, and routine temperature measurements, while strict control measures include closing public facilities (e.g., schools and restaurants), border closures, and traffic restrictions [10]. In China’s most severe case during the first round of the outbreak in 2020, the entire country was locked down and people were told to stay home unless they went out for daily needs, which effectively prevented the transmission of COVID-19 [11,12,13]. In the case of the United States, control measures could not be fully implemented during both the outbreak and containment phases, resulting in a significantly higher number of new cases confirmed daily compared with other countries. Therefore, the formulation and evaluation of the extent and impact of control measures are important for the interpretation of the multistage patterns.
Despite the control measures, seasonality could be another important factor influencing the multistage epidemic. It is still controversial whether the transmission of COVID-19 follows the similar seasonal patterns as other seasonal epidemics [14,15,16,17]. Epidemic curves show that outbreaks are concentrated in two hemispheres in different time periods [4]. In Northern Hemisphere countries, the number of new cases confirmed dramatic daily increases between November and January, whereas in Southern Hemisphere countries, a distinct outbreak period is observed between May and August. A number of studies have attempted to establish the relationship between the occurrence of COVID-19 and seasonally varying meteorological factors. Confirmed cases were found to be concentrated in air temperatures of 5 to 15 °C and absolute humidity of 3 to 10 g/m3 [18]. Experimental data and regional analyses show that coronavirus has a lower survival rate in a high temperature environment, suggesting that the virus may be less infectious in summer [19]. Lower wind speed significantly correlates with higher COVID-19 cases by affecting the transmission of aerosols [20,21], which are considered important virus carriers [22]. Despite the evidence for seasonality, there is debate that COVID-19 is not seasonal. Some studies find no relationship between temperature and the occurrence of COVID-19 [23,24,25] or state that the relationship exists only at the global level, but not at the regional level [26]. Given the recent outbreak in the winter of 2022, there is growing evidence that COVID-19 is a seasonal epidemic, but for policy decisions, the influences of seasonality have yet to be formulated and verified in model simulations.
Many studies have been conducted to examine the influences of control measures or seasonality separately, but few studies have been conducted to examine the joint influences of these two factors. In addition, most modeling studies for COVID-19 epidemic simulation focus on the first round of the outbreak and are not able to simulate the multistage epidemic. In this paper, we attempt to formulate the influences of seasonality and control measures on the transmission of the COVID-19 epidemic. We modify the most commonly used Suspected-Exposed-Infectious-Removed (SEIR) epidemic model by incorporating the influences of these two factors. We simulate the epidemic curves of twelve countries from two hemispheres from March 2020 to March 2021 using our modified model and interpret how these two factors interact and ultimately influence the transmission of the COVID-19 epidemic.

2. Materials and Methods

Seven Northern Hemisphere countries (United States, United Kingdom, Italy, France, Germany, Spain, and Turkey) and five Southern Hemisphere countries (Brazil, Chile, Australia, South Africa, Argentina) were selected and daily new confirmed cases from these countries are collected daily [4]. The selected countries, where considerable confirmed cases have been reported and where obvious multistage patterns are observed, are representative cases for the study of seasonality between two hemispheres. Countries near the equator, whose climates do not vary as much with the seasons, are not included [27]. The study period for each country begins when the number of new daily confirmed cases exceeds 100, to avoid the high uncertainty in reported data due to limited testing capabilities in the very early phase of the epidemic, and includes a full year for an entire seasonal cycle (approximately from March 2020 to March 2021). The influences of new variants and vaccination are excluded from the study period, and only seasonality and control measures are formulated.
The exact timing of implementation and mitigation of each round of control measures are gathered from official news reports or announcements (Table 1) to formulate the influences of control measures. An outbreak period is defined as the period from the request for control measures to the mitigation of control measures, while a containment period is inversely defined as the period from the mitigation of control measures to the request for a new round of control measures. In the selected countries, at least two rounds of outbreaks and weakening are observed.
The SEIR model divides the population into four subpopulations, the suspect population (S), the exposed population (E), the infectious population (I), and removed population (R), whose relationship can be formulated by the following differential equations:
d S ( t ) d t = γ R 0 ( t ) N S ( t ) I ( t ) ,
d E ( t ) d t = γ R 0 ( t ) N S ( t ) I ( t ) α E ( t ) ,
d I ( t ) d t = α E ( t ) γ I ( t ) ,
where R0 is the baseline reproduction number (defined as the average number of secondary infections caused by a single infected individual in a fully susceptible population [28]), N is the total population, γ is the removal rate (fixed at 0.34), and α is the infection rate (fixed at 0.19). The model is fitted with daily data of new confirmed cases and therefore the equation of the removed population is not presented here.
In a conventional single-stage SEIR model, R0 is usually fixed as a fixed value. However, in our multistage model, we replace R0 in Equations (1) and (2) with an adjusted, time-varying R 0 ( t ) , which is defined as follows:
R 0 ( t ) = R 0 ( t ) q ( t ) ,
where t denotes the days since the study period, q(t) is a time-varying quarantine index, and R0(t) is a time-varying and periodic parameter defined by the following periodic function to formulate the seasonality of COVID-19:
R 0 ( t ) = A 0 c o s   ( 2 π 365 t + ω ) + A 1 ,
with A0 being the amplitude of seasonal fluctuation, A1 the constant term, ω the initial phase, and the length of the periodic cycle is set as 365 days.
We construct a time-varying quarantine index q ( t ) , which falls between (0, 1], to formulate the influences of control measures. The quarantine index is a parameter to give a rough but quantitative estimation about the strict levels of control measures implemented in the whole country. In other words, to what extend can R 0 ( t ) be brought down by control measures. In a mathematical sense, when q ( t ) is equal to zero, R 0 ( t ) is zero and the epidemic is cut off completely. When q ( t ) is equal to one, R 0 ( t ) is equal to R 0 ( t ) and the epidemic spreads without any restrictions. In reality, q ( t ) cannot be one except the period without any control measures in the very early epidemic stage. In addition, q ( t ) cannot be zero because so far the transmission of COVID-19 cannot be cut off completely in any country. Thus, a lower bound q o and an upper bound q m , with 0 <   q o   <   q m   1, are set to constrain the values of q ( t ) , and the variations of q ( t ) during outbreak and containment periods are shown in Equations (6) and (7):
outbreak   period   q ( t ) = ( q m q o ) e ( ε o ( t t o ) ) + q o ,
containment   period   q ( t ) = ( q o q m ) e ( ε m ( t t m ) ) + q m ,
with subscripts o and m denoting the outbreak period and containment period, respectively, to and tm being the start dates of outbreak period and containment period. Parameters ε o and ε m are applied to controlling the the ascending and descending speed of exponent functions in Equations (6) and (7). In other words, higher values of ε o mean that the epidemic curves climb rapidly during outbreak periods, while lower values of ε m mean that the epidemic curves decline slowly during containment periods. For brevity, the values of q o , q m , ε o and ε m are assumed to share the same value in each epidemic stage. It should be noted that q m is set as one in the earliest stage as there were no control measures. In the modified SEIR model, in total seven parameters ( q o , q m , ω , A 0 , A 1 , ε o , and ε m ) should be fitted. Genetic algorithm is applied in this study to approximating the best parameters combination [29]. To accelerate the solution, coarse lower and upper bounds of these seven parameters are explored upfront, as shown in Table 2.

3. Results

Our modified SEIR model performs well in fitting the multistage COVID-19 epidemic curves with R2 ranging from 0.781 to 0.969 (Figure 1). The best fitting performance, with R2 greater than 0.94, is obtained in Germany, the United Kingdom, Italy, and Australia, where two or three rounds of alternating outbreak and downturn periods can be simulated. Our model manages to capture the main trend of the epidemic curves during each period and neglect the smaller noise (e.g., Brazil). In some cases, despite high fitness, some outbreaks cannot be captured. For example, the second round of the outbreak period in the United States and Argentina is ignored.
The results of the seven parameters with best fitting performance are listed in Table 3. According to the definition of the periodic function (Equation (5)), R 0 ( t ) reaches the maximum value at A 0   +   A 1 (t = 365- ω ) and the minimum value at A 1 A 0 (t = 182- ω ). The maximum value ranges from 3.53 (Brazil) to 6.42 (Germany) while the minimum value ranges from 1.02 (Turkey) to 4.76 (United States). The lowest minimum value of 1.02 in Turkey indicates that R 0 ( t ) cannot fall below one with only seasonal fluctuation and control measures must be implemented. The maximum values of seasonally fluctuating R 0 ( t ) of selected countries are projected in Figure 2. Countries from two hemispheres can be clearly distinguished, with R 0 ( t ) reaching maximum values in November and January in the Northern Hemisphere and March and June in the Southern Hemisphere. The value of q o ranges from 0.09 (Australia) to 0.32 (Brazil), indicating that control measures can reduce the capability to transmit during the outbreak period to at least one-third. The value of q m falls between 0.21 (United States) and 0.79 (Turkey), implying that the transmission capability can increase to at most four-fifths during the containment period. The major turning points of the epidemic curves are near the initial dates of the outbreak or containment periods when control measures are implemented or contained. In most countries, implementation of control measures can result in a significant decrease in new daily confirmed cases.
The relations of R 0 ( t ) , R 0 ( t ) , and q(t) are shown in Figure 3. At the earliest stage of the epidemic, when no control measures are implemented and COVID-19 spreads unimpeded, q ( t ) is equal to one and R 0 ( t ) is equal to R 0 ( t ) . During the outbreak periods, control measures are implemented and q ( t ) falls exponentially from maximum to minimum (from q m to q o ) while R 0 ( t ) also decreases, and the outbreak is prevented. An outbreak period is followed by a containment period, in which control measures are mitigated, with q resuming to q m and R 0 ( t ) rebounding as well. Outbreak and containment periods increase alternately, with q ( t ) varying between q m and q o and R 0 ( t ) goes up and down around one.

4. Discussion

The dates when the periodic R 0 ( t ) reach maximum values in two hemispheres are concentrated in two distinct periods, demonstrating that the transmission ability of COVID-19 indeed fluctuates seasonally. The amplitudes of R 0 ( t ) vary among countries, which may be attributed to the atmospheric pollution [30,31], demographic structure [20,32], healthcare condition, and other environmental or social factors [33]. In general, our model performs better in the northern hemisphere cases than the southern hemisphere cases. In some southern hemisphere cases (i.e., Argentina, Chile), the anti-seasonal outbreaks in January, 2021 could not be explained by seasonality and further efforts should be put on.
The formulation of the influences of control measures is significant for explaining the different forms of epidemic curves among selected countries. We define a time-varying quarantine index q ( t ) and assume that q rises or falls in an exponential form after the implementation or mitigation of control measures, and stays stable when it reaches maximum or minimum values. Our early attempts show that an exponential function has better performance in formulating the variation processes of quarantine index compared with a polynomial function or sinusoidal function. The upper bound q m is assumed to be greater than the lower bound q o in simulation, representing a lower level of control measures in the containment period than that of the outbreak period. The values ε o and ε m control the speed at which the quarantine index reaches its maximum or minimum values, in other words, how fast the restriction or mitigation policies take effect. In practice, the upper bound of ε o and ε m is set to three, which means that q reaches its maximum or minimum values within one day; thus, a value greater than three is meaningless. The quarantine index relies on the exact dates when the policies are claimed, and it does not perform well if these dates cannot be collected. For instance, control measures in the United States are claimed state-to-state rather than nationwide, and thus the second round of outbreak is not caught by this model. In addition, the quarantine index may be greater than one in some extreme cases with the contact rate of the populations far beyond the usual (i.e., the large-scale concerts or sports events), which are not considered in this study.
Seasonality and control measures are the two most significant driving factors that influence the transmission ability of COVID-19, and a single factor itself is not able to interpret the different multi-stage epidemic curves. The epidemic curve of the United Kingdom is an ideal case to reveal the combined effects of these two factors (Figure 4). At the earliest stage, the epidemic spread without any restrictions and R 0 ( t ) was equal to R 0 ( t ) . The first round of control measures was implemented on 23 March 2020, after which q ( t ) dropped from one to q o and R 0 ( t ) dropped as well. The control measures were mitigated on 11 May 2020, and the daily confirmed cases maintained a low level during May and August. With the seasonality force strengthening, the curve started to rise up in September, 2020 and the second round of outbreak occurred in October, 2020. With control measures being tightened again on 5 November 2020, R 0 ( t ) fell down to less than one and a short-term decline was recorded. However, when the cases were still at a high level, control measures were mitigated immediately on 2 December 2020. Under the superposition of high values of q ( t ) and R 0 ( t ) during November and January, the cases rebound rapidly and soon the third round of outbreak occurred. It’s an interesting phenomenon that the first round of mitigation resulted in a period with cases at a low level, while the second round of mitigation resulted in a rapid rebound, and this phenomenon can be observed in most northern hemisphere countries (i.e., United Kingdom; France; Italy; Spain). The major difference between these two rounds of mitigation demonstrates the existence of seasonality.
It is well known that the epidemic enters an expansion stage when the initial baseline reproduction number R0 is greater than one, and that a large-scale outbreak may subsequently occur. However, few studies have provided a reasonable and quantitative interpretation for the increase in R0. We replace the original R0 with our time-varying adjusted R 0 ( t ) in the modified model, and our study suggests that the increase in R 0 ( t ) is due to the strengthening of seasonality or to the weakening of control measures. The former is a gradual and periodic process, while the latter is a policy-decided and abrupt process, and a superposition of these two factors will lead to a large-scale outbreak.
In previous studies, a quarantine rate is often used to separate the quarantine population from the infectious population and to formulate the influences of control measures [34]. However, the quarantine rate is usually assigned manually and discretely [35]. A study on the Ebola epidemic formulated the influences of human intervention by an exponential function, but was only suitable for a single-stage epidemic [36]. Our model succeeds in establishing the continuous variation functions of the quarantine index during each period, and we multiply the periodic R 0 ( t ) by q ( t ) and keep the rest of the model unchanged instead of adding an additional subpopulation to the SEIR model. In the simulation, the fitted values of the parameters can vary widely and form different combinations but achieve similarly high R2 values, so all further interpretations for these parameters should be taken with caution. For example, q o is 0.17 in both the United Kingdom and the United States, but this does not mean that the control measures of these two countries are at the same level.
Many previous studies have indicated that the asymptomatic group and the undetected group play an important role in the transmission of the COVID-19 epidemic [37,38]. For simplicity, these two groups are not included in our model, but satisfactory adjustment performance can still be achieved. Moreover, new variants of COVID-19 (i.e., Beta; Delta; Omicron) are spreading throughout the world and dramatically changing the epidemic evolution patterns [39,40,41], and more data should be collected to verify the seasonality and intensity of their transmission ability. So far, the basic reproductive number of the newest variant Omicron has not been clearly elucidated, and it should be treated separately from other strains in the model simulation. Vaccination can effectively prevent transmission by reducing the suspected population. However, due to production limitation and validity of vaccination in new variants, herd immunity is difficult to achieve and control measures are still needed. In this paper, the influences of new variants and vaccination are excluded by shortening the study period before March 2021. If the study period is extended, our model does not perform well mainly because of these two factors.
It has been more than two years since COVID-19 appeared and swept across the world. Since there is no indication that the epidemic will disappear in a short period of time, it is of great importance to uncover the driving factors of the multistage COVID-19 epidemic, which will help policy makers to control the epidemic and avoid a large-scale outbreak. Some countries scaled down control measures only when daily confirmed cases increased to unacceptable levels, resulting in a significant number of deaths. Therefore, in terms of public health, we propose real-time calculation and monitoring of the R 0 ( t ) proposed in this paper. Once R 0 ( t ) breaks through an upward value and remains at a fairly high level for several days, control measures should be restricted in advance regardless of the strengthening of seasonality or the weakening of control measures to avoid a superposition of these two factors, since we cannot change the fluctuation of seasonality.

5. Conclusions

The aim of this study is to interpret how seasonality and control measures jointly determine the multistage waves of the COVID-19 epidemic. We formulate and integrate these two factors into our modified SEIR model and test our model with data from twelve countries in two hemispheres. Our model achieves high fitness and captures the multistage patterns well, and the influences of seasonality and control measures can be well formulated and interpreted. The multilevel epidemic is shaped by the joint influence of control measures and seasonality. If only seasonality plays a role, the basic reproductive number R0 cannot fall below one and control measures must be taken. An overlay of strong seasonality and a low level of control measures can lead to a dramatically rapid increase in reported cases.
Our study also has some limitations. First, some anti-seasonal outbreaks cannot be captured by our model. Second, the use of the quarantine index relies on the exact time of the measures, and does not work well if these time points cannot be collected. Finally, other long-term factors such as new variants, vaccination, reinfection and demographics are not accounted for in our study.

Author Contributions

Conceptualization, Y.Z. and Y.W.; methodology, Y.Z.; investigation, Y.Z. and Y.W.; data curation, Y.Z.; writing—original draft preparation, Y.Z.; writing—review and editing, Y.W.; visualization, Y.Z.; funding acquisition, Y.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by The Natural Science Foundation of China, grant number U1901215 and 42007205, and Natural Science Fund of Guangdong Province, grant number 2021A1515011375. This is contribution No. IS-3198 from GIGCAS.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fukutome, A.; Watashi, K.; Kawakami, N.; Ishikawa, H. Mathematical modeling of severe acute respiratory syndrome nosocomial transmission in Japan: The dynamics of incident cases and prevalent cases. Microbiol. Immunol. 2007, 51, 823–832. [Google Scholar] [CrossRef] [Green Version]
  2. Oraby, T.; Tyshenko, M.G.; Balkhy, H.H.; Tasnif, Y.; Quiroz-Gaspar, A.; Mohamed, Z.; Araya, A.; Elsaadany, S.; Al-Mazroa, E.; Alhelail, M.A.; et al. Analysis of the Healthcare MERS-CoV Outbreak in King Abdulaziz Medical Center, Riyadh, Saudi Arabia, June-August 2015 Using a SEIR Ward Transmission Model. Int. J. Environ. Res. Public Health 2020, 17, 2936. [Google Scholar] [CrossRef]
  3. Khedher, N.B.; Kolsi, L.; Alsaif, H. A multi-stage SEIR model to predict the potential of a new COVID-19 wave in KSA after lifting all travel restrictions. Alex. Eng. J. 2021, 60, 3965–3974. [Google Scholar] [CrossRef]
  4. COVID-19 Coronavirus Pandemic. Available online: https://www.worldometers.info/coronavirus (accessed on 20 March 2022).
  5. Lopez, L.; Rodo, X. A modified SEIR model to predict the COVID-19 outbreak in Spain and Italy: Simulating control scenarios and multi-scale epidemics. Results Phys. 2021, 21, 103746. [Google Scholar] [CrossRef]
  6. Liu, X.; Huang, J.; Li, C.; Zhao, Y.; Wang, D.; Huang, Z.; Yang, K. The role of seasonality in the spread of COVID-19 pandemic. Environ. Res. 2021, 195, 110874. [Google Scholar] [CrossRef]
  7. Wang, H.; Xu, K.; Li, Z.; Pang, K.; He, H. Improved Epidemic Dynamics Model and Its Prediction for COVID-19 in Italy. Appl. Sci. 2020, 10, 4930. [Google Scholar] [CrossRef]
  8. Zhang, R.; Li, Y.; Zhang, A.; Wang, Y.; Molina, M.J. Identifying airborne transmission as the dominant route for the spread of COVID-19. Proc. Natl. Acad. Sci. USA 2020, 117, 14857–14863. [Google Scholar] [CrossRef]
  9. Ruan, L.; Wen, M.; Zeng, Q.; Chen, C.; Huang, S.; Yang, S.; Yang, J.; Wang, J.; Hu, Y.; Ding, S.; et al. New Measures for the Coronavirus Disease 2019 Response: A Lesson from the Wenzhou Experience. Clin. Infect. Dis. 2020, 71, 866–869. [Google Scholar] [CrossRef] [Green Version]
  10. Nistal, R.; de la Sen, M.; Gabirondo, J.; Alonso-Quesada, S.; Garrido, A.J.; Garrido, I. A Study on COVID-19 Incidence in Europe through Two SEIR Epidemic Models Which Consider Mixed Contagions from Asymptomatic and Symptomatic Individuals. Appl. Sci. 2021, 14, 6266. [Google Scholar] [CrossRef]
  11. Hellewell, J.; Abbott, S.; Gimma, A.; Bosse, N.I.; Jarvis, C.I.; Russell, T.W.; Munday, J.D.; Kucharski, A.J.; Edmunds, W.J.; Sun, F.; et al. Feasibility of controlling COVID-19 outbreaks by isolation of cases and contacts. Lancet Glob. Health 2020, 8, e488–e496. [Google Scholar] [CrossRef] [Green Version]
  12. Pan, A.; Liu, L.; Wang, C.; Guo, H.; Hao, X.; Wang, Q.; Huang, J.; He, N.; Yu, H.; Lin, X.; et al. Association of Public Health Interventions with the Epidemiology of the COVID-19 Outbreak in Wuhan, China. JAMA 2020, 323, 1915–1923. [Google Scholar] [CrossRef] [Green Version]
  13. Prem, K.; Liu, Y.; Russell, T.W.; Kucharski, A.J.; Eggo, R.M.; Davies, N.; Flasche, S.; Clifford, S.; Pearson, C.A.; Munday, J.D.; et al. The effect of control strategies to reduce social mixing on outcomes of the COVID-19 epidemic in Wuhan, China: A modelling study. Lancet Public Health 2020, 5, e261–e270. [Google Scholar] [CrossRef] [Green Version]
  14. Batabyal, S. COVID-19: Perturbation dynamics resulting chaos to stable with seasonality transmission. Chaos Solitons Fractals 2021, 145, 110772. [Google Scholar] [CrossRef]
  15. Neher, R.A.; Dyrdak, R.; Druelle, V.; Hodcroft, E.B.; Albert, J. Potential impact of seasonal forcing on a SARS-CoV-2 pandemic. Swiss Med. Wkly. 2020, 150, w20224. [Google Scholar] [CrossRef] [Green Version]
  16. Steel, J.; Palese, P.; Lowen, A.C. Transmission of a 2009 pandemic influenza virus shows a sensitivity to temperature and humidity similar to that of an H3N2 seasonal strain. J. Virol. 2011, 85, 1400–1402. [Google Scholar] [CrossRef] [Green Version]
  17. Eccles, R. An explanation for the seasonality of acute upper respiratory tract viral infections. Acta Otolaryngol. 2002, 122, 183–191. [Google Scholar] [CrossRef]
  18. Huang, Z.; Huang, J.; Gu, Q.; Du, P.; Liang, H.; Dong, Q. Optimal temperature zone for the dispersal of COVID-19. Sci. Total Environ. 2020, 736, 139487. [Google Scholar] [CrossRef]
  19. Yao, M.; Zhang, L.; Ma, J.; Zhou, L. On airborne transmission and control of SARS-Cov-2. Sci. Total Environ. 2020, 731, 139178. [Google Scholar] [CrossRef]
  20. Coskun, H.; Yildirim, N.; Gunduz, S. The spread of COVID-19 virus through population density and wind in Turkey cities. Sci. Total Environ. 2021, 751, 141663. [Google Scholar] [CrossRef]
  21. Rendana, M. Impact of the wind conditions on COVID-19 pandemic: A new insight for direction of the spread of the virus. Urban Clim. 2020, 34, 100680. [Google Scholar] [CrossRef]
  22. Van Doremalen, N.; Bushmaker, T.; Morris, D.H.; Holbrook, M.G.; Gamble, A.; Williamson, B.N.; Tamin, A.; Harcourt, J.L.; Thornburg, N.J.; Gerber, S.I.; et al. Aerosol and Surface Stability of SARS-CoV-2 as Compared with SARS-CoV-1. N. Engl. J. Med. 2020, 382, 1564–1567. [Google Scholar] [CrossRef]
  23. Xie, J.; Zhu, Y. Association between ambient temperature and COVID-19 infection in 122 cities from China. Sci. Total Environ. 2020, 724, 138201. [Google Scholar] [CrossRef]
  24. Yao, Y.; Pan, J.; Liu, Z.; Meng, X.; Wang, W.; Kan, H.; Wang, W. No association of COVID-19 transmission with temperature or UV radiation in Chinese cities. Eur. Respir. J. 2020, 55. [Google Scholar] [CrossRef] [Green Version]
  25. Jüni, P.; Rothenbühler, M.; Bobos, P.; Thorpe, K.E.; Da Costa, B.R.; Fisman, D.N.; Slutsky, A.S.; Gesink, D. Impact of climate and public health interventions on the COVID-19 pandemic: A prospective cohort study. Can. Med. Assoc. J. 2020, 192, E566–E573. [Google Scholar] [CrossRef]
  26. Meyer, A.; Sadler, R.; Faverjon, C.; Cameron, A.R.; Bannister-Tyrrell, M. Evidence That Higher Temperatures Are Associated with a Marginally Lower Incidence of COVID-19 Cases. Front. Public Health 2020, 8, 367. [Google Scholar] [CrossRef]
  27. Lowen, A.C.; Steel, J.; Mubareka, S.; Palese, P. High temperature (30 degrees C) blocks aerosol but not contact transmission of influenza virus. J. Virol. 2008, 82, 5650–5652. [Google Scholar] [CrossRef] [Green Version]
  28. Kissler, S.M.; Tedijanto, C.; Goldstein, E.; Grad, Y.H.; Lipsitch, M. Projecting the transmission dynamics of SARS-CoV-2 through the postpandemic period. Science 2020, 368, 860–868. [Google Scholar] [CrossRef]
  29. Yarsky, P. Using a genetic algorithm to fit parameters of a COVID-19 SEIR model for US states. Math. Comput. Simul. 2021, 185, 687–695. [Google Scholar] [CrossRef]
  30. Bashir, M.F.; Benghoul, M.; Numan, U.; Shakoor, A.; Komal, B.; Bashir, M.A.; Bashir, M.; Tan, D. Environmental pollution and COVID-19 outbreak: Insights from Germany. Air Qual. Atmos. Health 2020, 13, 1385–1394. [Google Scholar] [CrossRef]
  31. Coccia, M. Factors determining the diffusion of COVID-19 and suggested strategy to prevent future accelerated viral infectivity similar to COVID. Sci. Total Environ. 2020, 729, 138474. [Google Scholar] [CrossRef]
  32. Mollalo, A.; Vahedi, B.; Rivera, K.M. GIS-based spatial modeling of COVID-19 incidence rate in the continental United States. Sci. Total Environ. 2020, 728, 138884. [Google Scholar] [CrossRef] [PubMed]
  33. Lian, X.; Huang, J.; Zhang, L.; Liu, C.; Liu, X.; Wang, L. Environmental Indicator for COVID-19 Non-Pharmaceutical Interventions. Geophys. Res. Lett. 2021, 149, 2020GL090344. [Google Scholar] [CrossRef] [PubMed]
  34. Li, W.; Gong, J.; Zhou, J.; Zhang, L.; Wang, D.; Li, J.; Shi, C.; Fan, H. An evaluation of COVID-19 transmission control in Wenzhou using a modified SEIR model. Epidemiol. Infect. 2021, 149, e2. [Google Scholar] [CrossRef] [PubMed]
  35. Carcione, J.M.; Santos, J.E.; Bagaini, C.; Ba, J. A Simulation of a COVID-19 Epidemic Based on a Deterministic SEIR Model. Front. Public Health 2020, 8, 230. [Google Scholar] [CrossRef] [PubMed]
  36. Chowell, G.; Hengartner, N.W.; Castillo-Chavez, C.; Fenimore, P.W.; Hyman, J.M. The basic reproductive number of Ebola and the effects of public health measures: The cases of Congo and Uganda. J. Theor. Biol. 2004, 229, 119–126. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Sun, D.; Long, X.; Liu, J. Modeling the COVID-19 Epidemic with Multi-Population and Control Strategies in the United States. Front. Public Health 2021, 9, 751940. [Google Scholar] [CrossRef]
  38. Liu, Z.; Magal, P.; Seydi, O.; Webb, G. Understanding Unreported Cases in the COVID-19 Epidemic Outbreak in Wuhan, China, and the Importance of Major Public Health Interventions. Biology 2020, 9, 50. [Google Scholar] [CrossRef] [Green Version]
  39. Yadav, P.D.; Sapkal, G.N.; Ella, R.; Sahay, R.R.; Nyayanit, D.A.; Patil, D.Y.; Deshpande, G.; Shete, A.M.; Gupta, N.; Mohan, V.K.; et al. Neutralization of Beta and Delta variant with sera of COVID-19 recovered cases and vaccinees of inactivated COVID-19 vaccine BBV152/Covaxin. J. Travel Med. 2021, 28, taab104. [Google Scholar] [CrossRef]
  40. Tareq, A.M.; Emran, T.B.; Dhama, K.; Dhawan, M.; Tallei, T.E. Impact of SARS-CoV-2 delta variant (B.1.617.2) in surging second wave of COVID-19 and efficacy of vaccines in tackling the ongoing pandemic. Hum. Vaccines Immunother. 2021, 17, 4126–4127. [Google Scholar] [CrossRef]
  41. Zhao, H.; Lu, L.; Peng, Z.; Chen, L.L.; Meng, X.; Zhang, C.; Ip, J.D.; Chan, W.M.; Chu, A.W.; Chan, K.H.; et al. SARS-CoV-2 Omicron variant shows less efficient replication and fusion activity when compared with delta variant in TMPRSS2-expressed cells. Emerg. Microbes Infect. 2022, 11, 277–283. [Google Scholar] [CrossRef]
Figure 1. The curves of fitted and actual daily new confirmed cases for the selected countries. (a) United States; (b) France; (c) Brazil; (d) United Kingdom; (e) Italy; (f) Turkey; (g) Germany; (h) Australia; (i) Spain; (j) Argentina; (k) South-Africa; (l) Chile.
Figure 1. The curves of fitted and actual daily new confirmed cases for the selected countries. (a) United States; (b) France; (c) Brazil; (d) United Kingdom; (e) Italy; (f) Turkey; (g) Germany; (h) Australia; (i) Spain; (j) Argentina; (k) South-Africa; (l) Chile.
Ijerph 19 06404 g001
Figure 2. The maximum values of R 0 of the selected countries.
Figure 2. The maximum values of R 0 of the selected countries.
Ijerph 19 06404 g002
Figure 3. The curves of R 0 ( t ) , R 0 ( t ) and q(t) of each selected country. (a) United States; (b) France; (c) Brazil; (d) United Kingdom; (e) Italy; (f) Turkey; (g) Germany; (h) Australia; (i) Spain; (j) Argentina; (k) South-Africa; (l) Chile.
Figure 3. The curves of R 0 ( t ) , R 0 ( t ) and q(t) of each selected country. (a) United States; (b) France; (c) Brazil; (d) United Kingdom; (e) Italy; (f) Turkey; (g) Germany; (h) Australia; (i) Spain; (j) Argentina; (k) South-Africa; (l) Chile.
Ijerph 19 06404 g003
Figure 4. The co-influences of seasonality and control measures on the multi-stage COVID-19 epidemic in the case of the United Kingdom.
Figure 4. The co-influences of seasonality and control measures on the multi-stage COVID-19 epidemic in the case of the United Kingdom.
Ijerph 19 06404 g004
Table 1. Start dates of outbreak and containment periods for the selected countries in this study.
Table 1. Start dates of outbreak and containment periods for the selected countries in this study.
CountriesOutbreak
Period 1
Containment Period 1Outbreak
Period 2
Containment Period 2Outbreak
Period 3
Containment Period 3
United States21 March 20201 May 2020NANA7 February 2021NA
France16 March 202011 May 20203 October 202011 November 20203 April 2020NA
Brazil25 March 20201 June 202024 November 2020NANANA
United Kingdom23 March 202011 May 20205 November 20202 December 20204 January 202117 May 2021
Italy9 March 20204 May 202025 October 202010 January 20213 April 20212 June 2021
Germany16 March 202020 April 202012 December 2020NANANA
Turkey1 April 202012 May 20208 November 202025 January 202129 April 202117 May 2021
Australia2 March 202027 April 20202 August 202013 September 20201 January 202129 January 2021
Spain13 March 202013 April 202025 October 202023 November 20208 January 20219 May 2021
Argentina20 March 202016 May 20201 July 202018 July 202026 October 20201 December 2020
South-Africa26 March 20201 June 202012 July 202017 August 202029 December 2020NA
Chile18 March 20207 August 20203 January 202123 March 2021NANA
Table 2. Upper and lower bounds of fitted parameters for the model.
Table 2. Upper and lower bounds of fitted parameters for the model.
ParametersUpper and Lower Bounds
q o (0, 0.5]
q m [0.1, 0.9]
ω [1, 365]
A 0 [0.1, 0.9]
A 1 [0.8, 2.0]
ε o (0, 3]
ε m (0, 3]
Table 3. Results of fitted parameters for the selected countries of this study.
Table 3. Results of fitted parameters for the selected countries of this study.
Countries ω A0A1Maximum Value
(A1 + A0)
Minimum Value
(A1A0)
q o q m ε o ε m
United States860.565.325.884.760.170.210.253.00
France1211.073.584.652.510.130.470.050.02
Brazil140.842.693.531.850.320.500.070.70
United Kingdom980.504.144.643.640.170.311.023.00
Italy1081.244.025.262.780.160.310.120.22
Germany1101.255.176.423.920.160.230.480.96
Turkey711.822.854.671.020.140.790.080.02
Australia2741.104.885.993.780.090.230.470.11
Spain660.633.273.902.640.230.400.540.05
Argentina3390.744.605.343.860.230.270.233.00
South-Africa3471.804.316.112.510.150.400.300.19
Chile3391.272.603.871.330.310.690.093.00
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zheng, Y.; Wang, Y. How Seasonality and Control Measures Jointly Determine the Multistage Waves of the COVID-19 Epidemic: A Modelling Study and Implications. Int. J. Environ. Res. Public Health 2022, 19, 6404. https://doi.org/10.3390/ijerph19116404

AMA Style

Zheng Y, Wang Y. How Seasonality and Control Measures Jointly Determine the Multistage Waves of the COVID-19 Epidemic: A Modelling Study and Implications. International Journal of Environmental Research and Public Health. 2022; 19(11):6404. https://doi.org/10.3390/ijerph19116404

Chicago/Turabian Style

Zheng, Yangcheng, and Yunpeng Wang. 2022. "How Seasonality and Control Measures Jointly Determine the Multistage Waves of the COVID-19 Epidemic: A Modelling Study and Implications" International Journal of Environmental Research and Public Health 19, no. 11: 6404. https://doi.org/10.3390/ijerph19116404

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop