Next Article in Journal
The Incidence and Predictors of Solid- and Hematological Malignancies in Patients with Giant Cell Arteritis: A Large Real-World Database Study
Previous Article in Journal
Which National Factors Are Most Influential in the Spread of COVID-19?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling of the Long-Term Epidemic Dynamics of COVID-19 in the United States

1
Wuhan Britain-China School, No.10 Gutian Rd., Qiaokou District, Wuhan 430022, China
2
Institute of Biophysics, School of Physics, Huazhong University of Science and Technology, Wuhan 430074, China
*
Authors to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2021, 18(14), 7594; https://doi.org/10.3390/ijerph18147594
Submission received: 11 June 2021 / Revised: 4 July 2021 / Accepted: 9 July 2021 / Published: 16 July 2021
(This article belongs to the Section Global Health)

Abstract

:
Coronavirus 2019 (COVID-19) is causing a severe pandemic that has resulted in millions of confirmed cases and deaths around the world. In the absence of effective drugs for treatment, non-pharmaceutical interventions are the most effective approaches to control the disease. Although some countries have the pandemic under control, all countries around the world, including the United States (US), are still in the process of controlling COVID-19, which calls for an effective epidemic model to describe the transmission dynamics of COVID-19. Meeting this need, we have extensively investigated the transmission dynamics of COVID-19 from 22 January 2020 to 14 February 2021 for the 50 states of the United States, which revealed the general principles underlying the spread of the virus in terms of intervention measures and demographic properties. We further proposed a time-dependent epidemic model, named T-SIR, to model the long-term transmission dynamics of COVID-19 in the US. It was shown in this paper that our T-SIR model could effectively model the epidemic dynamics of COVID-19 for all 50 states, which provided insights into the transmission dynamics of COVID-19 in the US. The present study will be valuable to help understand the epidemic dynamics of COVID-19 and thus help governments determine and implement effective intervention measures or vaccine prioritization to control the pandemic.

1. Introduction

Coronavirus disease 2019 (COVID-19) is a currently ongoing pandemic that has resulted in over 187 million confirmed cases and 4.04 million deaths globally [1]. Although several countries have the pandemic under relative control, all countries around the world, including the US, are still in the process of controlling the spread of COVID-19 [2]. In the US, a total of more than 27 million confirmed cases and over 480 thousand deaths have been reported as of 14 February 2021 [1]. Unfortunately, there are still no efficient antiviral drugs for the treatment of COVID-19. Vaccines and non-pharmaceutical interventions are therefore the only two available measures to control the transmission of the disease. As of 14 February 2021, more than 10 vaccines have been permitted by at least one national regulatory agency for distribution. However, current levels of vaccination are still insufficient to control the spread of the disease. In the US, over 97,000 new cases were reported on 14 February 2021 [1]. Therefore, government decisions on policies and interventions will be essential in controlling the spread of COVID-19 in the following months. Non-pharmaceutical interventional measures including social distancing, mandatory mask-wearing, isolation, and contact tracing are still crucial in controlling the spread of the disease [3,4,5,6,7,8,9,10,11].
Given the critical importance of non-pharmaceutical interventions in the control of COVID-19, modeling the past spread and predicting the future trends of the COVID-19 epidemic will play a vital role in decision-making by health officials regarding appropriate containment measures [12,13,14,15,16,17,18,19,20,21,22]. In the past year, many studies have been carried out to model or forecast the transmission dynamics of COVID-19 in various countries, regions, or cities around the world [23,24,25,26,27,28,29]. Among these, classical compartmental models have been widely used and obtained many successful predictions. Among different compartment models, the Susceptible-Infectious-Removed (SIR) model is the simplest model and provides the most primary transmission principles of all compartmental models [30]. Therefore, the SIR model and its extended versions have been widely used to describe the epidemic pattern of infectious diseases [31,32,33,34,35,36]. One basic parameter to describe the dynamics of an epidemic is the reproduction number, R 0 , which represents the average number of secondary cases caused by one case in a large susceptible population. When R 0 > 1 , an epidemic outbreak is expected to occur in a population, when R 0 < 1 , the epidemic is under control, and when R 0 = 1 , the epidemic is in relatively stable growth [31].
Traditionally, R 0 is assumed to be constant or arbitrarily change at specific times during the spread of a disease in compartmental models. This can be true for mildly infectious diseases like the flu [32], for which no preventive interventions are normally enforced by the government. However, this would not be the case for highly infectious diseases like COVID-19, for which no effective pharmaceutical remedies are available to control the disease. As mentioned before, many non-pharmaceutical interventions like mask-wearing, social distancing, prohibition of gathering, and school closure, have been enforced to control the spread of COVID-19. Therefore, with the development of the COVID-19 epidemic, the basic reproduction number R 0 is expected to gradually reduce due to the increasing number of preventive measures implemented over the course of the epidemic [37,38,39]. As such, the basic reproduction number R 0 will be time-dependent and should be represented by a dynamic value, R t , which is defined as the average number of secondary infectious individuals generated by one infectious individual at time t. To incorporate such time-dependent dynamics, some studies have modified the traditional SIR compartmental model to track the early depleting transmission dynamics or forecast the relatively short-term trends of the COVID-19 epidemic [40,41]. Despite some successes using these models, they were only proposed for modeling the spread of COVID-19 in the early stages. These models would not be suitable for modeling the long-term transmission dynamics for an extended period of time. In addition, current models were only designed to predict the trends of epidemic situations with one outburst; however, situations with multiple outbursts cannot be predicted. Analysis of the data reported by different states in the US with high population densities shows a situation in which there is one outburst that is seemingly controlled, followed by a second outburst. The second outbursts occur at around the same time for each state. Therefore, previous models designed to predict only one outburst are unable to predict the epidemic trends in these areas. In addition, the epidemic dynamics of COVID-19 in the US are highly demography-dependent. Some states only show one apparent outbreak, while the others show two. Moreover, some outbreaks occurred in March–April 2020, but some occurred in June–August 2020, which would also be demographically dependent.
To address these issues, we propose a new time-dependent SIR model, named T-SIR, to model the transmission dynamics of COVID-19 from 22 January 2020 to 14 February 2021 for the 50 states in the US. Our T-SIR model provides a comprehensive picture of different scenarios of the long-term COVID-19 epidemic in the US. We analyzed the relationships between the epidemic parameters of our T-SIR model and the demographic properties of the states. The present study will be significant in predicting the future spread of COVID-19 at the state level. Utilizing this model, epidemic trends in certain areas can be used to help health officials make appropriate decisions on public health policies, including social distancing, wearing face masks, travel bans, lockdowns, and vaccine prioritization.

2. Methods

2.1. T-SIR Model

We have used the SIR compartmental model as the basis to build our time-dependent epidemic model for the spread of COVID-19. First, the SIR model is a classical compartmental model that was initially used to model the epidemic patterns of the flu. Despite its simplicity, the SIR model provides the most basic principles of transmission dynamics and has been widely used to model the spread of various infectious diseases [31]. Second, the SIR model is the simplest compartmental model and includes the fewest assumptions and a minimal number of parameters [40]. Therefore, the model would involve the least amount of overfitting problems possible and give a general description of the epidemic dynamics. Third, the SIR model has been used to model the transmission dynamics of COVID-19 in previous studies, which have shown the overall best fit for the early-time data of COVID-19 among different models including SIR, SEIR, and branching process models [25].
In the SIR model, individuals in a population are assigned to three compartments or categories: susceptible (S), infectious (I), and removed (R). Here, the S compartment stands for the susceptible individuals who are not immune to the virus and can become infected when they are exposed to the virus. Susceptible individuals cannot infect others because they have not carried the virus yet. The I compartment represents those individuals who are carrying the virus and infectious. The R compartment indicates those removed individuals who have either recovered or died from the disease and become immune to further infection. In this compartmental model, two important parameters determine the transmission dynamics of a disease. One is the transmission rate, β , which characterizes the probability of a susceptible individual becoming infected when the individual is exposed to an infectious individual during contact. The other is the removal rate, γ , which stands for the probability of an infectious individual becoming a removed individual during a specific time interval. For a population of N individuals, the SIR model can be described by the following set of differential equations [31]:
d S d t = β I S N , d I d t = β I S N γ I , d R d t = γ I
where the total population, N = S t + I t + R t , remains a constant during the spread of the disease.
Since the time for an infectious individual to recover is relatively stable, the removal rate, γ , is often regarded as a constant. With these two parameters, β and γ , the basic reproduction number can be defined as R 0 = β / γ . Due to the awareness-driven behavior and public health interventions in the spread of COVID-19, previous studies have shown that the transmission rate γ and corresponding reproduction number R t will gradually decrease over time. Therefore, considering the time-dependent trend in the transmission dynamics of COVID-19 and the multiple-outbreak feature of the COVID-19 epidemic in the US, we propose a dynamic epidemic model, named T-SIR, with a time-dependent transmission rate, β t , for modeling the spread of COVID-19 in the US as follows:
β ( t ) = β 1 t σ 1 γ β 2 t σ 2 t 0 < t t 1 t 1 < t t 2 t > t 2
where β 1 and β 2 stand for the initial transmission rate constants for the first and second outbreaks. The time t is in the unit of days and ranges in 1 , . Between two outbreaks, there is a relatively stable stage with a reproduction number of R t 1 , i.e., β t γ . The parameter σ 1 and σ 2 are the corresponding exponents to describe the power decreasing feature of transmission rates and normally have a value between 0 and 1. Therefore, σ can be defined as the intervention parameter reflecting the force of intervention in the epidemic control of COVID-19. When σ = 0 , our time-dependent epidemic model T-SIR returns to the original SIR model. When σ is larger, the intervention force is stronger and the transmission rate β t will reduce faster with time t . As the removal rate γ is assumed a constant in the model, the dynamic reproduction number R t = β t / γ will show a similar time-dependent relationship to the transmission rate β t . In this study, we will use our time-dependent SIR model, shown in Equations (1) and (2), to model the transmission dynamics of COVID-19 from 22 January 2020 to 14 February 2021 for the 50 states in the US.

2.2. Data Sources and Processing

2.2.1. Epidemiological Data

Similar to the previous study [12], we obtained the COVID-19 epidemic data from the COVID-19 Data Repository managed by the Center for Systems Science and Engineering (CSSE) at Johns Hopkins University (JHU) [42]. First, the time series table file, “time_series_covid19_confirmed_US.csv”, was downloaded from the JHU CCSE COVID-19 Dataset website [42]. The time series data are updated once a day and includes the time series data of the cumulative daily confirmed cases from 22 January 2020 to 14 February 2021 in the 50 US states as of our download date. The original data, which is presented daily, reported the confirmed cases at the county level in each state. We then processed the data to obtain the total number of confirmed cases on different dates for each state by summing the corresponding data for all counties in the respective state. The processed data of confirmed cases can be freely accessed at https://github.com/HuangDerek/T-SIR (accessed on 3 July 2021).

2.2.2. Geographical and Demographic Data

The most recent and accurate population data for each state were obtained from 1 July 2020 estimations from the US Census Bureau [43]. Total land area data were obtained from the US Census Bureau [44]. Population density data were obtained by calculations of the population divided by total land area for each state.

2.2.3. Mobility Data

The mobility data table of the US states during the COVID-19 pandemic, named 2020_US_Region_Mobility_Report.csv, was first downloaded from the Google COVID-19 Community Mobility Reports website at https://www.google.com/covid19/mobility/ (accessed on 14 February 2021). The original region CSV table contains six types of mobility data including “retail_and_recreation”, “grocery_and_pharmacy”, “parks”, “transit_stations”, “workplaces”, and “residential” at the country, state, and county level, which are characterized by the “percent_change_from_baseline”. We then extracted the state-level mobility data from the table file because the present study focused on the epidemic dynamics of different states. Our examinations showed that the six types of mobility data are well correlated with each other. Therefore, for simplicity, we used the “retail_and_recreation” mobility data, which is also the most complete data in the table, to represent societal mobility during the COVID-19 pandemic. The processed mobility data for the different states can be freely accessed at https://github.com/HuangDerek/T-SIR (accessed on 3 July 2021).

2.3. Fitting

To accurately model the epidemic trends, the T-SIR model parameters, β 1 , β 2 , σ 1 , σ 2 , and γ , are required to be determined by fitting the model with the cumulative data, where the corresponding dynamic reproduction numbers R 1 = β 1 / γ and R 2 = β 2 / γ . We fitted our T-SIR model using the grid search of parameters R 1 , R 2 , σ 1 , σ 2 , and γ in reasonable ranges. The ranges for R x , σ x , and γ were set to 0.1 ,   5 , 0 ,   0.3 , and 0 ,   0.3 with steps sized 0.1, 0.01, and 0.01, respectively [25]. The parameters of our T-SIR model, R 1 , R 2 , σ 1 , σ 2 , and γ , were then automatically determined through the best fit between the predicted data of T-SIR model and the reported COVID-19 data according to the mean squared error (MSE), where the three times, t 0 , t 1 , and t 2 , in Equation (2) were also automatically determined during the model fitting.

3. Results

3.1. COVID-19 Epidemics of the 50 States in the US

3.1.1. Times of COVID-19 Outbreaks

Figures S1 and S2 show the graphs of the daily numbers of confirmed cases and new cases with respect to the time period of 22 January 2020 to 14 February 2021 for each of the 50 US states, where the daily new cases were obtained by the confirmed cases of the current day minus those of the previous day. It can be seen from the figures that all states experienced at least one COVID-19 outbreak during the period. During a COVID-19 outbreak, the number of daily confirmed cases increased rapidly at an increasing rate in the early stages, which corresponds to the increasing numbers of daily new cases. After the number of daily new cases reached a maximum, the increase of daily confirmed cases became slower until it reached a relatively stable state. Therefore, a COVID-19 outbreak showed an obvious peak on the curve of daily new cases and a significant plateau on the curve of daily confirmed cases following the outbreak (Figure 1).
From Figures S1 and S2, one can also see that all the outbreaks occurred in one of the three different time windows (Figure 1). The first time window for COVID-19 outbreaks was March–April, 2020, example states for which include Connecticut, Delaware, Maryland, Massachusetts, and New Jersey (Figure 2a). The main reason for the outbreak during this time period may be due to the high infectivity of COVID-19 in the cold weather because the virus survives longer and has a higher infectious ability at lower temperatures [45]. In addition, part of the outbreak in this time window may have been caused by the early introduction of the coronavirus to these states with little to no measures enacted to control the spread because in the early stages, both the individuals and the government had little knowledge about the virus and were not aware of the high risk of the virus. The period between the first and second outbreaks where the spread of the coronavirus becomes relatively stable is the result of warmer weather in which the disease experiences a loss in the reproduction rate as well as the enacted measures to prevent the spread of the disease. The second time window for COVID-19 outbreaks is June–August 2020. Examples of states that underwent an outbreak during the period include Alabama, Arizona, California, Florida, and Texas (Figure 2a). Theoretically, an outbreak was not expected to occur during this summer period because the infectious ability of the virus decreases at a higher temperature. However, the outbreaks occurred due to human behavior more than the infectivity of the virus. This is indeed the case because several “Black Lives Matters (BLM)” protests occurred in many states during this period of time. These protests may have led to the fast spread of the COVID-19 virus due to large gatherings of people. As shown in Figure 2a,b, those states with large and extreme protests are likely to have an outbreak during this time, which confirms the important impact of protests on the COVID-19 outbreak. Additionally, summer break, which caused a higher mobility of people and lower awareness of the coronavirus, resulting in less incentive to maintain preventive measures, may have played a role. As shown in Figure 2a, many states with this outbreak have coastline access to the ocean or neighbor coastal states. As such, these states would involve more gatherings of people due to vacations, tours, beach activities, etc., which favor the spread of COVID-19. The third time window for COVID-19 outbreaks is October 2020–January 2021. Interestingly, all states show a COVID-19 outbreak during this time (Figures S1 and S2). The outbreak in this time widow could be attributed to both the infectious virus and human behavior. On one hand, the virus becomes more infectious at lower temperatures and has a higher transmission ability in cold weather. On the other hand, after a long period of fighting against the virus, individuals may tend to lower their sensitivity to the virus due to a lack of awareness and the relaxation of the public health interventions by governments due to economic reasons.

3.1.2. Properties of COVID-19 Outbreaks

The transmission dynamics of COVID-19 are also state-dependent in terms of epidemic outbreaks in the US. Roughly, the 50 states can be grouped into two categories according to the number of COVID-19 outbreaks (Figures S1 and S2). One category is those with only one COVID-19 outbreak in October 2020–January 2021 (Figure 1a). Examples of states included in this category are Alaska, Missouri, and South Dakota (Figure S1). In these states, the number of confirmed cases shows a logistic growth similar to that of the classic SIR model, which corresponds to a single peak in the curve of daily new cases. The other category is those with two outbreaks (Figure 1b,c). Examples of states included in this category are New Jersey, Arizona, and Massachusetts. The number of daily confirmed cases grows logistically then flattens out into a linear growth for a period of time and then grows logistically once again, which corresponds to two major peaks in the curves of daily new cases.
The category of states with two outbreaks can also be further divided into two subcategories. The states in these two subcategories all had a second outbreak in October 2020–January 2021 but differ in the time of the first outbreak. One subcategory had the first outbreak in March–April 2020 (Figure 1b). Examples of states included in this subcategory include Connecticut, Massachusetts, New Jersey, and New York. The other subcategory had the first outbreak in June–August 2020 (Figure 1c). Examples of states included in this subcategory are Alabama, Arizona, California, Florida, and Texas. Interestingly, a few states like Hawaii, Louisiana, and Maryland seem to show three COVID-19 outbreaks, with outbreaks occurring in March–April, June–August, and October–January 2021. Nevertheless, they may be roughly grouped into the two-outbreak category due to one outbreak having much less impact than the other two. Given such two-outbreak properties of the US states, where one outbreak can also be effectively regarded as two outbreaks (i.e., one is null and the other is real), we will be able to model the transmission dynamics of COVID-19 in the 50 states of the US using our two-outbreak epidemic model shown in Equations (1) and (2).

3.1.3. Demographic Impact on COVID-19 Outbreaks

Given that the 50 states of the US have very different populations and areas, it would be valuable to examine how such geographical and/or demographic features impact the spread of COVID-19. Table 1 highlights those two-outbreak states in ranked lists of the states according to their population densities, populations, and areas. It can be seen from the table that, overall, those states with higher populations or population densities tend to have two outbreaks. It can be understood because the higher population density causes the virus to be easier to spread due to the smaller distance between individuals and the larger population, leading to a higher possibility for the virus to spread due to the fact that there are more exposed people.
Further examination of two-outbreak states also reveals that those states with higher population densities tended to have an outbreak in March–April 2020, while states with higher populations tended to have an outbreak in June–August 2020. As shown in Table 1, the states with the top seven population densities all had two outbreaks, with the first outbreak occurring in March–April (Table 1 and Figure 2c), and the states with the top three populations all experienced two outbreaks, with the first outbreak occurring in June–August (Table 1 and Figure 2e). This kind of demography-dependent phenomenon is also understandable. As mentioned above, the outbreak in March–April is mostly due to the infectivity of the COVID-19 virus in cold weather. Therefore, it makes sense that the outbreaks in March–April tended to occur in those states with a high population density because the virus would be easier to spread due to the closer distance between individuals. In contrast, the outbreaks in June–August is mainly attributed to human behavior like BLM protests and high social mobility. For a state with a higher population, there would be more possible protests and higher social mobility due to the larger number of people, which explained the outbreaks in June–August, 2020.
Interestingly, the two-outbreak feature also seems to correlate with the total land area of the state. Namely, the outbreaks in March–April tended to occur in states with smaller areas, while the outbreaks in June–August tended to occur in states with larger areas (Table 1). However, it should be noted that such an area-dependent phenomenon on the two-outbreak feature is not truly due to the total area of the state, but rather due to an indirect effect of population density. As shown in Figure 2, small-area states tend to have a high population density (Figure 2d), and large-area states tend to have a large population (Figure 2f). In addition, some states with low population density such as Nebraska, Nevada, and Idaho, also had two outbreaks. For Nevada and Idaho, we can see that the first outbreak for these two states occurred in July. Figure 2 shows that these two states are the neighbors of large-population states like California and Washington, which had many BLM protests, explaining the increase in transmission. As for Nebraska, although it shows an increase of daily confirmed cases in March-April, the outbreak was relatively weak, which would be due to its overall low population density. These findings are expected to provide valuable guidance for COVID-19 vaccine prioritization in the US.

3.2. Modeling the COVID-19 Dynamics of the US States

3.2.1. Fitting of the T-SIR Model

Based on our two-outbreak T-SIR epidemic model in Equations (1) and (2), we have modeled the transmission dynamics of COVID-19 in the 50 US states through the fitting of the T-SIR model to the daily confirmed cases from 22 January 2020 to 14 February 2021, where the compartments I and R of the model are regarded as the confirmed cases. It can be seen from the figures that overall, our T-SIR model fits the dynamics of the daily confirmed COVID-19 cases very well and is consistent with the data from most of the states (Figures S3 and S4). The model not only fits the COVID-19 dynamics of two-outbreak states like Alabama, Connecticut, Massachusetts, New Jersey, and New York, but also describes the epidemic trends of one-outbreak states like Alaska, Colorado, Kentucky, Ohio, and West Virginia (Figure 3a–c). Similar consistency between the predicted and real data can also be observed in the fitting of daily new cases, demonstrating the robustness of our T-SIR model (Figure 3). Nevertheless, for a few states, such as Hawaii, Iowa, and Louisiana, our T-SIR model seemed to not fit the epidemic dynamics of COVID-19 very well (Figures S3 and S4), which may be due to the multiple waves of the epidemic dynamics in these states and/or inherent limitations in the reported data.

3.2.2. Implications of T-SIR Model Parameters

Through the model fitting, we can obtain the epidemic parameters of our T-SIR model, β 1 , β 2 , σ 1 , σ 2 , and γ , for the 50 states of the US, where β 1 and σ 1 are for the first outbreak of COVID-19 occurring in March–April or June–August, and β 2 and σ 2 are for the second outbreak occurring in October–January. These epidemic parameters gave a quantitative description of the COVID-19 transmission dynamics of different states and will enable the investigation of the impact of state geographic and demographic data on the spread of the disease. Such investigations will provide a deep understanding of the epidemic dynamics of COVID-19 in different states and help the government create the corresponding measures to prevent the spread of COVID-19. As such, we have investigated the relationship between the epidemic parameters and the demographic data for the 50 US states. Specifically, the epidemic parameters for fitting the model using respective state geographic and demographic data were compiled. Pairwise comparisons were made by examining their Pearson correlations. A significant correlation between two different variables corresponds to an important linear relationship between the two factors.
Figure 4a–c shows the relationships among the epidemic parameters, β 1 , σ 1 , and R 1 = β 1 / γ , for the first outbreak. High correlations can be seen among three pairs of parameters, which give a correlation coefficient of 0.954 for R 1 vs. β 1 , 0.869 for σ 1 vs. β 1 , and 0.894 for σ 1 vs. R 1 . Given the relationship of R 1 = β 1 / γ , the high correlation between the reproduction number R 1 and transmission rate β 1 suggests that the removal rate γ is relatively constant for different states (Figure 4a). This means that people infected by the virus take about the same time to recover in different states, which is consistent with the general findings for infectious diseases. It is known that the invention parameter σ is an indicator of the government’s preventive measures and individuals’ self-protection to contain the coronavirus. A higher σ corresponds to more effective preventive measures. It can be seen from Figure 4b that there is a high correlation between σ 1 and β 1 . This relationship can be understood because both individuals and the government adopt more and stricter preventive measures when the transmission rates are higher. These preventive measures may be both mandatory and voluntary, as higher caution is taken when large numbers of people are becoming infected. Figure 4c also shows a high correlation between σ 1 and R 1 . This relationship is due to σ 1 being positively related to β 1 . , and β 1 being positively correlated R 1 in the formulaic relationship R 1 = β 1 / γ , where the γ is relatively a constant. Similar relationships exist among the epidemic parameters, β 2 , σ 2 , and R 2 = β 2 / γ , for the second outbreak can also be observed, which give a correlation coefficient of 0.896 for R 2 vs. β 2 , 0.724 for σ 2 vs. β 2 , and 0.895 for σ 2 vs. R 2 (Figure 4d–f). The reasoning for these relationships is the same as that of the first outbreak.
As the intervention parameter, σ , is an indicator of the response of the government and individuals to the epidemic outbreak and thus the most important parameter of our T-SIR model, we will focus the investigation on the σ parameter. Figure 5a,b show the relationships between the σ parameters and the population densities of the 50 US states for the first and second outbreaks, respectively. It can be seen from the figure that σ 1 has a positive correlation of ρ = 0.425 with the population density for the first outbreak but shows a weakly negative correlation of ρ = 0.125 with the population density for the second outbreak. The different relationships for the two outbreaks may be understood because σ is the overall effect of both government interventions and individuals’ self-protection measures. During the first outbreak, the COVID-19 virus was new to the government and individuals. In this case, both the government and individuals in states with higher population densities tended to take more preventive measures to contain the coronavirus. For the second outbreak, the negative correlation between σ 2 and population density would be due to relaxation of preventive measures for the virus, which resulted in an apparently lower σ for states with higher population densities because the virus spread more easily in higher-density populations if no strict interventions were enforced. Compared to the case made for population density, σ shows a weaker relationship with the population (Figure 5c,d). This is understandable because it is the population density rather than the population that determines the spread of a disease in a normal epidemic.
We further investigated the relationship between σ and social mobility, where the mobility during the first outbreak is roughly taken from the average over Mach–June and that for the second outbreak is taken from the average over October–January. As shown in Figure 5e, a significant negative correlation of ρ = 0.311 between σ 1 and mobility can be seen for the first outbreak. This is explained by a higher σ indicating that the government enforced more preventive measures. Measures to control the spread of the disease including social distancing and stay-at-home isolation orders decreased people’s mobility. However, a lower correlation of ρ = 0.209 is observed for the second outbreak (Figure 5f). This may be understood because the effect of σ comes from two parts: government interventions and individual self-protection measures, where government interventions have a much higher impact on social mobility than the self-imposed measures that individuals place on themselves. With the relaxation of government interventions towards the second outbreak, σ came more from the less impactful individual self-protection measures, which explains the smaller correlation between σ 2 and mobility.
In addition, we also compared the magnitudes of the epidemic parameters, R , β , and σ , of the first and second outbreaks, which are shown in Figure 6. It can be seen from the figure that the initial transmission rates β 1 and corresponding reproduction numbers R 1 for the first outbreak have a wider range than those of the second outbreak. This is understandable because people knew little about COVID-19 at the beginning, and the government had not enforced intervention measures to prevent the transmission of the virus yet, which would lead to a high transmission rate (Figure 6a). At the same time, some states did not show an outbreak during this early stage due to little introduction of COVID-19, which yielded a low transmission rate. However, after the epidemic had gone on for a few months, both the government and individuals all had a better understanding of COVID-19 and adopted more efficient prevention measures to stop the spread of the virus in the late stages of the COVID-19 epidemic, resulting in relatively moderate transmission rates for the second outbreak (Figure 6b). Similarly, as the government and individuals adopt the corresponding preventive measures to contain the virus according to the transmission rates, the range of σ will be similar to that of β , resulting in an overall wider range of σ for the first outbreak than for the second outbreak (Figure 6c).

4. Discussion

Based on our time-dependent T-SIR model, we have modeled the transmission dynamics of COVID-19 from 22 January 2020 to 14 February 2021 in the 50 states of the US. The predicted and reported data are well consistent in both the daily confirmed cases and new cases, suggesting the robustness of our T-SIR model. Nevertheless, there may be limitations in the reported COVID-19 data because the availability of the reported data depends on various factors. First, the availability of confirmed data strongly depends on COVID-19 testing capacity. This issue was especially critical in the early stages of the epidemic when the number of infected individuals was far beyond the testing capacity for COVID-19 capacity. Second, the testing capacity is also state or city-dependent due to the differences in the hospital resources in different states or cities. Moreover, many people might not go to the hospital due to having mild or no symptoms even if they were infected, and thus these infected individuals would not be counted. In addition, compared to the confirmation of infected individuals, it is much more difficult to track recovered individuals. Therefore, it is difficult to accurately separate infected and recovered individuals in the confirmed cases. These issues in the reported data significantly limit the accurate modeling of COVID-19 epidemic dynamics, especially for the dynamics of daily new cases (Figure S4), though the basic trends can be modeled.
Despite the various limitations in the data of reported COVID-19 cases, our T-SIR model still modeled the time-dependent transmission dynamics of the virus for the 50 states of the US well, demonstrating its robustness to model the impact of non-pharmaceutical interventions on the spread of COVID-19. In other words, our T-SIR model truly reflects the unique epidemic dynamics of COVID-19 in the US. In our model, the transmission rate, β t , reduces quickly in the early stages and then changes more slowly as time passes. This is consistent with the social scenario of individuals and public health responses to an outbreak of COVID-19 in the US. During the early stages of an outbreak, the number of new cases and deaths increased rapidly, which resulted in strong pressure on society for prevention. Therefore, during this period, both the individuals and government had a tendency to adopt strict measures and/or interventions to reduce the transmission of the disease, resulting in a fast-decreasing transmission rate. When the COVID-19 epidemic entered the later stages of an outbreak after some time, the transmission rate became lower and the epidemic seemed to be under control. Correspondingly, the epidemic entered a relatively stable stage. During this stage, both the individuals and government had a tendency to relax the measures and/or interventions because of less pressure from the epidemic and economic considerations, which would result in a second epidemic outbreak, as shown in the dynamics of COVID-19 (Figures S1 and S2).
Although there is good consistency between the reported data of COVID-19 and the predicted values by our T-SIR model, it should be noted that the reported data gather only the infected people controlled by the health system, as mentioned before. These controls depend on the number of tests performed on people. For instance, people with symptoms who are not controlled for the system and the asymptomatic are not considered, although they are also infected. Therefore, the reported data are only a fraction of real infected cases. These inherent factors may need to be considered when interpreting the present T-SIR and other epidemic models.
In addition, it is known that vaccination is the ultimate measure to fight against COVID-19 due to the fact that there are no effective drugs for the treatment of the disease [47]. Starting from around December 2020, an increasing number of people around the world began to have the opportunity to receive COVID-19 vaccines. In the US, the percentage of the population that is fully vaccinated against COVID-19 increased quickly from 0.48% on 15 January 2021, to 23.47% on 15 April 2020, to 46.31% on 30 June 2021 [48]. The share of people who had received at least one dose of the COVID-19 vaccine increased from 3.17% on 15 January 2021, to 37.62% on 15 April 2020, to 54.02% on 30 June 2021 [48]. The increasing number of vaccinated people play an important role in stopping the transmission of COVID-19, which may in part explain the decreasing new cases of COVID-19 after January 2021 in all US states (Figures S1 and S2).

5. Conclusions

We have investigated the transmission dynamics of COVID-19 from 22 January 2020 to 14 February 2021 in the 50 states of the US and also proposed a time-dependent Susceptible-Infectious-Removed model, named T-SIR, to model the epidemic dynamics of the virus. It was revealed that the 50 states can be roughly divided into two groups according to the number of occurred epidemic outbreaks: the states with one outbreak and the states with two outbreaks. Those states with higher populations and/or population densities tend to have a higher chance of two COVID-19 outbreaks than the other states. Further examinations revealed that the states with high population density tend to have their first outbreak occur in March–April, which is mostly due to the high infectivity of COVID1-9 in the cold weather, while the states with a high population tended to have their first outbreak occur in June–August, which may be mostly due to protests and gatherings. All states had a COVID-19 outbreak in October–January, which would be attributed to the relaxation of both governmental intervention policies and individual self-prevention measures as well as the colder weather. We further fitted our T-SIR epidemic model to the transmission dynamics of COVID-19 in the 50 states of the US. It was shown that our T-SIR model can well describe the epidemic dynamics of all 50 states. The relationships between the corresponding epidemic parameters and demographic data were extensively investigated, providing insights into the spread of COVID-19 in different states. The present study will be valuable for not only understanding the epidemic dynamics of COVID-19 but also helping the government make effective virus-control policies such as mask wearing and vaccine prioritization in the US.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ijerph18147594/s1, Figure S1: The daily confirmed cases of COVID-19 from 22 January 2020 to 14 February 2021 for the 50 states in the US, Figure S2: The daily new cases of COVID-19 from 22 January 2020 to 14 February 2021 for the 50 states in the US, where the data were calculated from the confirmed data, Figure S3: The fitting of our T-SIR model to the daily confirmed cases of COVID-19 for the 50 states in the US, where the dashed lines indicate the predicted data of T-SIR model, Figure S4: The fitting of our T-SIR model to the daily new cases of COVID-19 for the 50 states in the US, where the dashed lines indicate the predicted data of T-SIR model. Data were smoothed using a Savitzky–Golay filter for the daily new cases.

Author Contributions

Conceptualization, D.H., S.-Y.H. and Y.X.; methodology, D.H, H.T., S.-Y.H. and Y.X.; software, D.H., Q.W. and S.-Y.H.; validation, D.H. and S.-Y.H.; formal analysis, D.H.; investigation, D.H.; resources, S.-Y.H.; data curation, D.H., H.T. and S.-Y.H.; writing—original draft preparation, D.H.; writing—review and editing, D.H., S.-Y.H. and Y.X.; visualization, D.H. and Q.W.; supervision, S.-Y.H. and Y.X.; project administration, S.-Y.H. and Y.X.; funding acquisition, S.-Y.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the startup grant of Huazhong University of Science and Technology.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used for this study are available at https://github.com/HuangDerek/T-SIR (accessed on 3 July 2021) or from the corresponding authors upon reasonable request.

Acknowledgments

The authors would like to acknowledge Chaolong Wang for the helpful discussions and for providing valuable suggestions. This work was in part supported by the startup grant of Huazhong University of Science and Technology.

Conflicts of Interest

The authors declare no conflict of interests.

References

  1. WHO Coronavirus (COVID-19) Dashboard. Available online: https://covid19.who.int/ (accessed on 15 July 2021).
  2. Coronavirus Disease (COVID-19) Situation Reports. Available online: https://www.who.int/emergencies/diseases/novel-coronavirus-2019/situation-reports (accessed on 14 February 2021).
  3. Anderson, R.M.; Heesterbeek, H.; Klinkenberg, D.; Hollingsworth, T.D. How will country-based mitigation measures influence the course of the COVID-19 epidemic? Lancet 2020, 395, 931–934. [Google Scholar] [CrossRef]
  4. Ngonghala, C.N.; Iboi, E.; Eikenberry, S.; Scotch, M.; MacIntyre, C.R.; Bonds, M.H.; Gumel, A.B. Mathematical assessment of the impact of non-pharmaceutical interventions on curtailing the 2019 novel Coronavirus. Math. Biosci. 2020, 325, 108364. [Google Scholar] [CrossRef]
  5. Teslya, A.; Pham, T.M.; Godijk, N.G.; Kretzschmar, M.E.; Bootsma, M.C.J.; Rozhnova, G. Impact of self-imposed prevention measures and short-term government-imposed social distancing on mitigating and delaying a COVID-19 epidemic: A modelling study. PLoS Med. 2020, 17, e1003166. [Google Scholar] [CrossRef]
  6. 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] [PubMed] [Green Version]
  7. Hartley, D.M.; Perencevich, E.N. Public Health Interventions for COVID-19: Emerging Evidence and Implications for an Evolving Public Health Crisis. JAMA 2020, 323, 1908–1909. [Google Scholar] [CrossRef] [Green Version]
  8. Zhang, C.; Chen, C.; Shen, W.; Tang, F.; Lei, H.; Xie, Y.; Cao, Z.; Tang, K.; Bai, J.; Xiao, L.; et al. Impact of population movement on the spread of 2019-nCoV in China. Emerg. Microbes Infect. 2020, 9, 988–990. [Google Scholar] [CrossRef] [Green Version]
  9. Cao, Z.; Tang, F.; Chen, C.; Zhang, C.; Guo, Y.; Lin, R.; Huang, Z.; Teng, Y.; Xie, T.; Xu, Y.; et al. Impact of Systematic Factors on the Outbreak Outcomes of the Novel COVID-19 Disease in China: Factor Analysis Study. J. Med. Internet. Res. 2020, 22, e23853. [Google Scholar] [CrossRef]
  10. Gilman, R.T.; Mahroof-Shaffi, S.; Harkensee, C.; Chamberlain, A.T. Modelling interventions to control COVID-19 outbreaks in a refugee camp. BMJ Glob. Health 2020, 5, e003727. [Google Scholar] [CrossRef] [PubMed]
  11. Alfonso Viguria, U.; Casamitjana, N. Early Interventions and Impact of COVID-19 in Spain. Int. J. Environ. Res. Public Health 2021, 18, 4026. [Google Scholar] [CrossRef] [PubMed]
  12. Modeling COVID-19 scenarios for the United States. Nat. Med. 2021, 27, 94–105. [CrossRef] [PubMed]
  13. Hao, X.; Cheng, S.; Wu, D.; Wu, T.; Lin, X.; Wang, C. Reconstruction of the full transmission dynamics of COVID-19 in Wuhan. Nature 2020, 584, 420–424. [Google Scholar] [CrossRef] [PubMed]
  14. Rader, B.; Scarpino, S.V.; Nande, A.; Hill, A.L.; Adlam, B.; Reiner, R.C.; Pigott, D.M.; Gutierrez, B.; Zarebski, A.E.; Shrestha, M.; et al. Crowding and the shape of COVID-19 epidemics. Nat. Med. 2020, 26, 1829–1834. [Google Scholar] [CrossRef] [PubMed]
  15. Monod, M.; Blenkinsop, A.; Xi, X.; Hebert, D.; Bershan, S.; Tietze, S.; Baguelin, M.; Bradley, V.C.; Chen, Y.; Coupland, H.; et al. Age groups that sustain resurging COVID-19 epidemics in the United States. Science 2021, 371, 6536. [Google Scholar] [CrossRef] [PubMed]
  16. Dehning, J.; Zierenberg, J.; Spitzner, F.P.; Wibral, M.; Neto, J.P.; Wilczek, M.; Priesemann, V. Inferring change points in the spread of COVID-19 reveals the effectiveness of interventions. Science 2020, 369, 6500. [Google Scholar] [CrossRef]
  17. Spassiani, I.; Sebastiani, G.; Palù, G. Spatiotemporal Analysis of COVID-19 Incidence Data. Viruses 2021, 13, 463. [Google Scholar] [CrossRef]
  18. Loeffler-Wirth, H.; Schmidt, M.; Binder, H. Covid-19 Transmission Trajectories-Monitoring the Pandemic in the Worldwide Context. Viruses 2020, 12, 777. [Google Scholar] [CrossRef]
  19. Röst, G.; Bartha, F.A.; Bogya, N.; Boldog, P.; Dénes, A.; Ferenci, T.; Horváth, K.J.; Juhász, A.; Nagy, C.; Tekeli, T.; et al. Early Phase of the COVID-19 Outbreak in Hungary and Post-Lockdown Scenarios. Viruses 2020, 12, 708. [Google Scholar] [CrossRef]
  20. Roosa, K.; Lee, Y.; Luo, R.; Kirpich, A.; Rothenberg, R.; Hyman, J.M.; Yan, P.; Chowell, G. Short-term Forecasts of the COVID-19 Epidemic in Guangdong and Zhejiang, China: February 13–23, 2020. J. Clin. Med. 2020, 9, 596. [Google Scholar] [CrossRef] [Green Version]
  21. Aguas, R.; White, L.; Hupert, N.; Shretta, R.; Pan-Ngum, W.; Celhay, O.; Moldokmatova, A.; Arifi, F.; Mirzazadeh, A.; Sharifi, H.; et al. Modelling the COVID-19 pandemic in context: An international participatory approach. BMJ Glob. Health 2020, 5, e003126. [Google Scholar] [CrossRef]
  22. Gabutti, G.; d'Anchera, E.; De Motoli, F.; Savio, M.; Stefanati, A. The Epidemiological Characteristics of the COVID-19 Pandemic in Europe: Focus on Italy. Int. J. Environ. Res. Public Health 2021, 18, 2942. [Google Scholar] [CrossRef]
  23. De Sousa, L.E.; Neto, P.H.O.; Filho, D. Kinetic Monte Carlo model for the COVID-19 epidemic: Impact of mobility restriction on a COVID-19 outbreak. Phys. Rev. 2020, 102, 032133. [Google Scholar] [CrossRef]
  24. Weitz, J.S.; Park, S.W.; Eksin, C.; Dushoff, J. Awareness-driven behavior changes can shift the shape of epidemics away from peaks and toward plateaus, shoulders, and oscillations. Proc. Natl. Acad. Sci. USA 2020, 117, 32764–32771. [Google Scholar] [CrossRef]
  25. Bertozzi, A.L.; Franco, E.; Mohler, G.; Short, M.B.; Sledge, D. The challenges of modeling and forecasting the spread of COVID-19. Proc. Natl. Acad. Sci. USA 2020, 117, 16732–16738. [Google Scholar] [CrossRef]
  26. Maugeri, A.; Barchitta, M.; Battiato, S.; Agodi, A. Estimation of Unreported Novel Coronavirus (SARS-CoV-2) Infections from Reported Deaths: A Susceptible-Exposed-Infectious-Recovered-Dead Model. J. Clin. Med. 2020, 9, 1350. [Google Scholar] [CrossRef]
  27. Reno, C.; Lenzi, J.; Navarra, A.; Barelli, E.; Gori, D.; Lanza, A.; Valentini, R.; Tang, B.; Fantini, M.P. Forecasting COVID-19-Associated Hospitalizations under Different Levels of Social Distancing in Lombardy and Emilia-Romagna, Northern Italy: Results from an Extended SEIR Compartmental Model. J. Clin. Med. 2020, 9, 1492. [Google Scholar] [CrossRef]
  28. Siraj, A.; Worku, A.; Berhane, K.; Aregawi, M.; Eshetu, M.; Mirkuzie, A.; Berhane, Y.; Siraj, D. Early estimates of COVID-19 infections in small, medium and large population clusters. BMJ Glob. Health 2020, 5, e003055. [Google Scholar] [CrossRef] [PubMed]
  29. Anđelić, N.; Baressi Šegota, S.; Lorencin, I.; Jurilj, Z.; Šušteršič, T.; Blagojević, A.; Protić, A.; Ćabov, T.; Filipović, N.; Car, Z. Estimation of COVID-19 Epidemiology Curve of the United States Using Genetic Programming Algorithm. Int. J. Environ. Res. Public Health 2021, 18, 959. [Google Scholar] [CrossRef]
  30. Chowell, G.; Sattenspiel, L.; Bansal, S.; Viboud, C. Mathematical models to characterize early epidemic growth: A review. Phys. Life Rev. 2016, 18, 66–97. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Tolles, J.; Luong, T. Modeling Epidemics with Compartmental Models. JAMA 2020, 323, 2515–2516. [Google Scholar] [CrossRef] [PubMed]
  32. Bootsma, M.C.; Ferguson, N.M. The effect of public health measures on the 1918 influenza pandemic in U.S. cities. Proc. Natl. Acad. Sci. USA 2007, 104, 7588–7593. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Tan, Z.J.; Zou, X.W.; Jin, Z.Z. Percolation with long-range correlations for epidemic spreading. Phys. Rev. 2000, 62, 8409–8412. [Google Scholar] [CrossRef]
  34. Shao, Z.G.; Tan, Z.J.; Zou, X.W.; Jin, Z.Z. Epidemics with pathogen mutation on small-world networks. Phys. Stat. Mech. Appl. 2006, 363, 561–566. [Google Scholar] [CrossRef]
  35. Leontitsis, A.; Senok, A.; Alsheikh-Ali, A.; Al Nasser, Y.; Loney, T.; Alshamsi, A. SEAHIR: A Specialized Compartmental Model for COVID-19. Int. J. Environ. Res. Public Health 2021, 18, 2667. [Google Scholar] [CrossRef]
  36. Maugeri, A.; Barchitta, M.; Battiato, S.; Agodi, A. Modeling the Novel Coronavirus (SARS-CoV-2) Outbreak in Sicily, Italy. Int. J. Environ. Res. Public Health 2020, 17, 4964. [Google Scholar] [CrossRef] [PubMed]
  37. Haushofer, J.; Metcalf, C.J.E. Which interventions work best in a pandemic? Science 2020, 368, 1063–1065. [Google Scholar] [CrossRef] [PubMed]
  38. Chinazzi, M.; Davis, J.T.; Ajelli, M.; Gioannini, C.; Litvinova, M.; Merler, S.; Pastore, Y.P.A.; Mu, K.; Rossi, L.; Sun, K.; et al. The effect of travel restrictions on the spread of the 2019 novel coronavirus (COVID-19) outbreak. Science 2020, 368, 395–400. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Cabrera, M.; Córdova-Lepe, F.; Gutiérrez-Jara, J.P.; Vogt-Geisse, K. An SIR-type epidemiological model that integrates social distancing as a dynamic law based on point prevalence and socio-behavioral factors. Sci. Rep. 2021, 11, 10170. [Google Scholar] [CrossRef]
  40. Law, K.B.; Peariasamy, K.M.; Gill, B.S.; Singh, S.; Sundram, B.M.; Rajendran, K.; Dass, S.C.; Lee, Y.L.; Goh, P.P.; Ibrahim, H.; et al. Tracking the early depleting transmission dynamics of COVID-19 with a time-varying SIR model. Sci. Rep. 2020, 10, 21721. [Google Scholar] [CrossRef]
  41. Sun, J.; Chen, X.; Zhang, Z.; Lai, S.; Zhao, B.; Liu, H.; Wang, S.; Huan, W.; Zhao, R.; Ng, M.T.A.; et al. Forecasting the long-term trend of COVID-19 epidemic using a dynamic model. Sci. Rep. 2020, 10, 21122. [Google Scholar] [CrossRef] [PubMed]
  42. JHU CSSE COVID-19 Dataset. Available online: https://github.com/CSSEGISandData/COVID-19/tree/master/csse_covid_19_data/ (accessed on 14 February 2021).
  43. Bereau, U.C. Change in Resident Population of the 50 States, the District of Columbia, and Puerto Rico: 1910 to 2020. Available online: https://www2.census.gov/programs-surveys/decennial/2020/data/apportionment/population-change-data-table.pdf (accessed on 14 February 2021).
  44. Bureau, U.C. State Area Measurements and Internal Point Coordinates. Available online: https://www.census.gov/geographies/reference-files/2010/geo/state-area.html (accessed on 14 February 2021).
  45. He, J.; Tao, H.; Yan, Y.; Huang, S.Y.; Xiao, Y. Molecular Mechanism of Evolution and Human Infection with SARS-CoV-2. Viruses 2020, 12, 428. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Savitzky, A.; Golay, M.J.E. Smoothing and Differentiation of Data by Simplified Least Squares Procedures. Anal. Chem. 1964, 36, 1627–1639. [Google Scholar] [CrossRef]
  47. Jain, V.; Schwarz, L.; Lorgelly, P. A Rapid Review of COVID-19 Vaccine Prioritization in the U.S.: Alignment between Federal Guidance and State Practice. Int. J. Environ. Res. Public Health 2021, 18, 3483. [Google Scholar] [CrossRef] [PubMed]
  48. Mathieu, E.; Ritchie, H.; Ortiz-Ospina, E.; Roser, M.; Hasell, J.; Appel, C.; Giattino, C.; Rodés-Guirao, L. A global database of COVID-19 vaccinations. Nat. Hum. Behav. 2021. [Google Scholar] [CrossRef]
Figure 1. The daily confirmed cases (upper row) and corresponding new cases (lower row) of COVID-19 for three selected states: Missouri (a), New Jersey (b), and Arizona (c), in the US that represent three categories (one outbreak; two outbreaks, with first outbreak in March–April; and two outbreaks, with first outbreak in June–August, respectively) of transmission dynamics in the US, where the dashed yellow lines indicate the epidemic outbreak peaks.
Figure 1. The daily confirmed cases (upper row) and corresponding new cases (lower row) of COVID-19 for three selected states: Missouri (a), New Jersey (b), and Arizona (c), in the US that represent three categories (one outbreak; two outbreaks, with first outbreak in March–April; and two outbreaks, with first outbreak in June–August, respectively) of transmission dynamics in the US, where the dashed yellow lines indicate the epidemic outbreak peaks.
Ijerph 18 07594 g001
Figure 2. The color-coded maps of the US states. (a) The states with two epidemic outbreaks occurring in March–April and October–January are highlighted in green, and those with two epidemic outbreaks occurring in June–August and October–January are highlighted in yellow. (b) Protests across the US by the number of cities and towns with rallies or protests in a state, the data for which were taken from https://www.usatoday.com/in-depth/graphics/2020/06/03/map-protests-wake-george-floyds-death/5310149002. (c) Population densities and (e) populations of the states with two epidemic outbreaks occurring in March–April and October–January are indicated by a red star, and those with two epidemic outbreaks occurring in June–August and October–January are indicated by a yellow star. (d) The population densities vs. areas for the 50 states in the US. (f) The populations vs. areas for the 50 states in the US.
Figure 2. The color-coded maps of the US states. (a) The states with two epidemic outbreaks occurring in March–April and October–January are highlighted in green, and those with two epidemic outbreaks occurring in June–August and October–January are highlighted in yellow. (b) Protests across the US by the number of cities and towns with rallies or protests in a state, the data for which were taken from https://www.usatoday.com/in-depth/graphics/2020/06/03/map-protests-wake-george-floyds-death/5310149002. (c) Population densities and (e) populations of the states with two epidemic outbreaks occurring in March–April and October–January are indicated by a red star, and those with two epidemic outbreaks occurring in June–August and October–January are indicated by a yellow star. (d) The population densities vs. areas for the 50 states in the US. (f) The populations vs. areas for the 50 states in the US.
Ijerph 18 07594 g002
Figure 3. The fitting of our T-SIR model to the daily confirmed cases (upper row) and corresponding new cases (lower row) of COVID-19 for three selected states: Colorado (a), New Jersey (b), and New York (c) in the US that represent three categories of transmission dynamics in the US, where the dashed lines indicate the predicted data of T-SIR model. Data were smoothed using a Savitzky–Golay filter for daily new cases [46], where the light-blue shade indicates the standard deviations of reported data.
Figure 3. The fitting of our T-SIR model to the daily confirmed cases (upper row) and corresponding new cases (lower row) of COVID-19 for three selected states: Colorado (a), New Jersey (b), and New York (c) in the US that represent three categories of transmission dynamics in the US, where the dashed lines indicate the predicted data of T-SIR model. Data were smoothed using a Savitzky–Golay filter for daily new cases [46], where the light-blue shade indicates the standard deviations of reported data.
Ijerph 18 07594 g003
Figure 4. The pairwise relationships and their correlation coefficients of three epidemic parameters, R , β , and σ , of the T-SIR model for the first outbreak (ac) and second outbreak (df), where the solid lines are the linear fittings of the data.
Figure 4. The pairwise relationships and their correlation coefficients of three epidemic parameters, R , β , and σ , of the T-SIR model for the first outbreak (ac) and second outbreak (df), where the solid lines are the linear fittings of the data.
Ijerph 18 07594 g004
Figure 5. The intervention parameter σ vs. population density (a,b), population (c,d), and social mobility (e,f) for the 50 states in the US for the first and second outbreaks, respectively. The solid lines are the linear fittings of the data and the corresponding correlation coefficients are also shown.
Figure 5. The intervention parameter σ vs. population density (a,b), population (c,d), and social mobility (e,f) for the 50 states in the US for the first and second outbreaks, respectively. The solid lines are the linear fittings of the data and the corresponding correlation coefficients are also shown.
Ijerph 18 07594 g005
Figure 6. The comparison between the three parameters of the T-SIR model for the first outbreaks (blue) and second outbreaks (green) of the 50 US states. (a) Reproduction numbers R . (b) Transmission rates β . (c) Intervention parameters σ .
Figure 6. The comparison between the three parameters of the T-SIR model for the first outbreaks (blue) and second outbreaks (green) of the 50 US states. (a) Reproduction numbers R . (b) Transmission rates β . (c) Intervention parameters σ .
Ijerph 18 07594 g006
Table 1. The ranked lists of the 50 states in the US according to their population densities, populations, and areas, where the states with two epidemic outbreaks occurring in March–April and October–January are highlighted in green, those with two epidemic outbreaks occurring in June–August and October–January are highlighted in yellow, and those with only one epidemic outbreak occurring in October–January are not highlighted.
Table 1. The ranked lists of the 50 states in the US according to their population densities, populations, and areas, where the states with two epidemic outbreaks occurring in March–April and October–January are highlighted in green, those with two epidemic outbreaks occurring in June–August and October–January are highlighted in yellow, and those with only one epidemic outbreak occurring in October–January are not highlighted.
Population DensityPopulationArea
RankStateDensity (/km2)RankStatePopulationRankStateArea (km2)
1New Jersey393.181California39,368,0781Alaska1,723,337
2Rhode Island264.222Texas29,360,7592Texas695,662
3Massachusetts252.183Florida21,733,3123California423,967
4Connecticut247.754New York19,336,7764Montana380,831
5Maryland188.475Pennsylvania12,783,2545New Mexico314,917
6Delaware153.096Illinois12,587,5306Arizona295,234
7New York136.857Ohio11,693,2177Nevada286,380
8Florida127.618Georgia10,710,0178Colorado269,601
9Pennsylvania107.179North Carolina10,600,8239Oregon254,799
10Ohio100.7210Michigan9,966,55510Wyoming253,335
11California92.8611New Jersey8,882,37111Michigan250,487
12Illinois83.9212Virginia8,590,56312Minnesota225,163
13Virginia77.5413Washington7,693,61213Utah219,882
14North Carolina76.0514Arizona7,421,40114Idaho216,443
15Indiana71.6115Massachusetts6,893,57415Kansas213,100
16Georgia69.5916Tennessee6,886,83416Nebraska200,330
17Tennessee63.0917Indiana6,754,95317South Dakota199,729
18South Carolina62.9218Missouri6,151,54818Washington184,661
19New Hampshire56.4319Maryland6,055,80219North Dakota183,108
20Hawaii49.6920Wisconsin5,832,65520Oklahoma181,037
21Kentucky42.7821Colorado5,807,71921Missouri180,540
22Texas42.2122Minnesota5,657,34222Florida170,312
23Washington41.6623South Carolina5,218,04023Wisconsin169,635
24Michigan39.7924Alabama4,921,53224Georgia153,910
25Alabama36.2525Louisiana4,645,31825Illinois149,995
26Wisconsin34.3826Kentucky4,477,25126Iowa145,746
27Louisiana34.2427Oregon4,241,50727New York141,297
28Missouri34.0728Oklahoma3,980,78328North Carolina139,391
29West Virginia28.4429Connecticut3,557,00629Arkansas137,732
30Arizona25.1430Utah3,249,87930Alabama135,767
31Minnesota25.1331Iowa3,163,56131Louisiana135,659
32Vermont25.0332Nevada3,138,25932Mississippi125,438
33Mississippi23.6533Arkansas3,030,52233Pennsylvania119,280
34Arkansas22.0034Mississippi2,966,78634Ohio116,098
35Oklahoma21.9935Kansas2,913,80535Virginia110,787
36Iowa21.7136New Mexico2,106,31936Tennessee109,153
37Colorado21.5437Nebraska1,937,55237Kentucky104,656
38Oregon16.6538Idaho1,826,91338Indiana94,326
39Utah14.7839West Virginia1,784,78739Maine91,633
40Maine14.7340Hawaii1,407,00640South Carolina82,933
41Kansas13.6741New Hampshire1,366,27541West Virginia62,756
42Nevada10.9642Maine1,350,14142Maryland32,131
43Nebraska9.6743Montana1,080,57743Hawaii28,313
44Idaho8.4444Rhode Island1,057,12544Massachusetts27,336
45New Mexico6.6945Delaware986,80945Vermont24,906
46South Dakota4.4746South Dakota892,71746New Hampshire24,214
47North Dakota4.1847North Dakota765,30947New Jersey22,591
48Montana2.8448Alaska731,15848Connecticut14,357
49Wyoming2.3049Vermont623,34749Delaware6,446
50Alaska0.4250Wyoming582,32850Rhode Island4,001
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Huang, D.; Tao, H.; Wu, Q.; Huang, S.-Y.; Xiao, Y. Modeling of the Long-Term Epidemic Dynamics of COVID-19 in the United States. Int. J. Environ. Res. Public Health 2021, 18, 7594. https://doi.org/10.3390/ijerph18147594

AMA Style

Huang D, Tao H, Wu Q, Huang S-Y, Xiao Y. Modeling of the Long-Term Epidemic Dynamics of COVID-19 in the United States. International Journal of Environmental Research and Public Health. 2021; 18(14):7594. https://doi.org/10.3390/ijerph18147594

Chicago/Turabian Style

Huang, Derek, Huanyu Tao, Qilong Wu, Sheng-You Huang, and Yi Xiao. 2021. "Modeling of the Long-Term Epidemic Dynamics of COVID-19 in the United States" International Journal of Environmental Research and Public Health 18, no. 14: 7594. https://doi.org/10.3390/ijerph18147594

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