Spatio-temporal Analysis of COVID-19: A Global Study

authors:

avatar Sajjad Rahimi Pordanjani ORCID 1 , 2 , avatar Maryam Mohammadian ORCID 3 , avatar Somayeh Derakhshan ORCID 3 , avatar Fatemeh Hadavandsiri ORCID 4 , avatar Seyed Saeed Hashemi Nazari ORCID 5 , 6 , * , avatar Mohammad Hossein Panahi ORCID 6 , 3 , **

Social Determinants of Health Research Center, Semnan University of Medical Sciences, Semnan, Iran
Assistant Professor of Epidemiology, Department of Community Medicine, School of Medicine, Semnan University of Medical Sciences, Semnan, Iran
Department of Epidemiology, School of Public Health and Safety, Shahid Beheshti University of Medical Sciences, Tehran, Iran
Cancer Research Center, Shahid Beheshti University of Medical Sciences, Tehran, Iran
Department of Epidemiology, Prevention of Cardiovascular Disease Research Center, School of Public Health and Safety, Shahid Beheshti University of Medical Sciences, Tehran, Iran
Department of Epidemiology, Safety Promotion and Injury Prevention Research Center, School of Public Health and Safety, Shahid Beheshti University of Medical Sciences, Tehran, Iran
Corresponding Authors:

how to cite: Rahimi Pordanjani S, Mohammadian M, Derakhshan S, Hadavandsiri F, Hashemi Nazari S S, et al. Spatio-temporal Analysis of COVID-19: A Global Study. Middle East J Rehabil Health Stud. 2023;10(2):e131344. https://doi.org/10.5812/mejrh-131344.

Abstract

Background:

COVID-19 is an international public health emergency in the world.

Objectives:

The aim of the present study is to determine the geographic pattern and temporal trend of Coronavirus disease 2019 incidence, fatality, and recovery rates worldwide.

Methods:

The present ecological study is a mixed exploratory study. The study population included Patients with COVID-19, recovered individuals, and deaths from COVID-19 from October 1, 2019, until June 30, 2021, worldwide. Descriptive analysis included the calculation cumulative incidence rate (CIR), case fatality rate (CFR), and case recovery rate (CRR) of COVID-19. Global Moran's I and Anselin Local Moran's I tests were used for spatial analysis. The joinpoint regression analysis was used to examine the time trend by ArcGIS, Joinpoint, and SPSS software.

Results:

The average cumulative incidence rate was 1077 in 106 individuals; also, the average case recovery rate and average case fatality rate were %72.81 and %3.21, respectively. Global Moran's I index measured for CIR was 0.159. The results of Anselin's local Moran's I, high-high cluster, consists of some countries in South America and in southern and Western Europe and central and western Asia. The temporal trend of changes in the incidence rate and CRR of COVID-19 were incremental, and the average annual percentage change from October 2019 to June 2021 increased by 44.4% and 3.2%, respectively (P < 0.001), but CFR decreased by -0.3% and was not significant (P > 0.05).

Conclusions:

As regards the specific spatial pattern of fatality and recovery rate of COVID-19, it seems essential to consider spatial conditions and environmental factors which are related to the incidence and fatality of COVID-19 in different regions, as well as the necessity of upgrading the care system in high-risk areas, in order to have better management and control of the pandemic and optimal function in early diagnosis, proper treatment, and high vaccination coverage.

1. Background

The world faced an international public health emergency in late December 2019 that started in Wuhan, China (1) numerous cases appeared rapidly in Asian countries and spread to other continents, which led to on March 11, 2020, the World Health Organization declared it a global pandemic, the main cause of which was severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) (2). Simultaneously with the wide spread of COVID-19 in the world, the need to study this communicable disease from a geographical point of view to answer different needs and questions was highlighted (2). This is due to the fact that COVID-19 disease has unknown features (2), Many climatic and spatial variables can effect on changes in the pattern of COVID-19 occurrence (3). Geography defines spatial, territorial, locational, and landscape aspects (4).

One of the most important research methods in this aspect of the pandemic period is spatial analysis (5). Spatial analysis that considers geographical development includes a set of statistical methods and trends (6, 7) that, by examining the incidence and prevalence of diseases, prepares them in the form of geographical maps, and the result is an explanation of dispersion (8). Thus, disease mapping has become necessary in public health (9).

Disease monitoring provides the basis for conducting research programs on the cause of the disease and evaluating early public health interventions (6). The incidence, prevalence, and deaths of diseases change over time (10). Temporal statistics are used to identify time trends. The temporal trend is determined by the annual percent change (APC) (11). Since the disease has spread rapidly in many regions, the pandemic has become a challenge for health decision-makers and Health service managers (12). So understanding of the spatiotemporal dynamics of COVID-19 is crucial for the prevention of incident cases. Because it helps policymakers reduce the prevalence rate and spread of disease (3, 13). According to the climate variation and socioeconomic differences in different countries, spatial analysis can determine disease distribution worldwide (14, 15) and also help comprehend patterns of incidence, mortality, morbidity, and burden of COVID-19 (16).

2. Objectives

In this regard, in the present study, we aim to understand the geographic pattern and temporal trend of the 2019 Coronavirus disease incidence worldwide to strengthen healthcare services and actions in high-risk regions. Also, this study provides information about morbidity, mortality, and recovery of COVID-19 disease from October 1, 2019, until June 30, 2021.

3. Methods

3.1. Design and Type of Study

This ecological study was designed as a mixed exploratory study. In exploratory studies, investigators discover and describe patterns of disease. The mixed study indicates that both temporal trends (Time trend) and spatial patterns (multiple-group study) are investigated (6, 17).

3.2. Target and Study Population

The target population consists of all patients infected with COVID-19 worldwide. The study population includes countries where information on morbidity, mortality, and recovery was available from October 1, 2019, until June 30, 2021.

3.3. Data Collection and Data Source

In this study, three main types of data sources were used. The data on the population at risk of COVID-19 were obtained from the World Bank website (18), and for obtaining the number of COVID-19 patients from October 2019 until June 2021, each country used CDC and WORDOMETER websites (19). Another source contained patients’ geographical coordinates data obtained from the georeferenced layer of Google Maps in each country.

3.4. Statistical Analysis

In this study, descriptive analysis of data was performed to determine the CIR, CFR, and CRR of COVID-19 in each county from October 1, 2019, until June 30, 2021. Data analysis was done with ArcGIS, Joinpoint, and SPSS software at a 0.05 significance level.

3.5. Cumulative Incidence Rate, Case Fatality Rate, and Case Recovery Rate Measurements

In order to compute the cumulative incidence rate of COVID-19, the number of newly diagnosed cases during a specific period was divided into the total population at risk during the beginning of the period in 106 individuals (20). CIRs measure the risk of developing the disease (21). Case fatality rate (CFR) and case recovery rate (CRR) were specified as the proportion (percent) of people who died and recovered from COVID-19 among all individuals diagnosed with COVID-19 over a certain period of time at the global level, respectively. The CFR is a measure of the severity of mortality or virulence of a disease and is sometimes used to predict the prognosis of the disease. The CRR reflects the quality of care and the strength of the health system (17, 22).

3.6. Temporal Trend

With the Joinpoint regression program version 4.7.0.0, the cumulative number of new cases, CFR, and CRR were calculated monthly and annually in order to evaluate the temporal trend of COVID-19. In this study, CIR, CFR, and CRR were dependent variables, and months and countries were independent variables. In order to compute the average annual percentage change (AAPC) and the annual percentage change (APC), we employed the natural log-linear model according to the Poisson distribution of dependent variables. The Grid search method was used to calculate the joint points where the incidence rate changes drastically (10). The Monte Carlo method estimated the P-value with 4499 iterations. The precision of models was examined by the sum of squared errors (SSE) and Bayesian information criterion (BIC). Regression line slope changes and the estimated APC-AAPC values were obtained for the annual trend of COVID-19 as well as; in order to evaluate the seasonal trend of COVID-19; the incidence rate was calculated by month. Finally, monthly slope alterations of trend, joinpoint times, and amounts of APC and AAPC were computed by joinpoint regression analysis.

3.7. Spatial Pattern Analysis and Cluster Identification

A spatial pattern for COVID-19 was determined by mapping the CIR, CFR, and CRR of each country using ArcGIS software version 10.8.

Global Moran's I was used to designate spatial autocorrelation of COVID-19 in the world. The value of global Moran's I is variable from -1 to 1 (23). The values near +1 show clustered distribution, values near -1 show dispersion, and zero indicates the random distribution of COVID-19 in the area under study (23).

In the next step, Anselin Local Moran's I was used to define the low-risk and high-risk cluster locations of COVID-19. This index focuses on the central polygon and compares it with the surrounding areas. Anselin Local Moran's I split polygons into five categories, containing low-high, high-low, low-low, high-high, and not significant (10). According to this category, in the H-H cluster, firstly, the incidence of disease in these areas is higher than the average incidence of disease in the entire study area, and secondly, the incidence of disease in this area and the surrounding neighbors is high, which is called a hot spot (10). The L-L cluster occurs when the incidence of disease in that cluster is lower than the world average, and the incidence of disease in the central polygon and surrounding countries is low, which is called cold spots[a1]. This index also identifies outliers, which are identified as H-L and L-H (10).

In order to indicate more precise clusters of COVID-19 in the world, an optimized hot spot analysis (OHSA) was employed. This method, according to the Getis-Ord Gi*, with the optimization of distance band for determining hot and cold spots (24). Using incremental spatial correlation, OHSA automatically selects the optimal bandwidth value at which spatial autocorrelation and clustering intensity are maximum and then uses it as an analytical scale for the Getis-Ord Gi * statistic (24). This will show the output's best possible results to identify hot and cold spots. OHSA demonstrates results at 90%, 95%, and 99% confidence levels (24).

This study is approved under the ethical approval code of IR.SBMU.RETECH.REC.1398.876.

The webpage of the ethical approval code is: (ethics.research.ac.ir/ProposalCertificateEn.php?id=125693)

4. Results

4.1. Descriptive

In the present study, the average cumulative incidence rate was 1077 per 106 individuals. Also, the average case recovery rate and average case fatality rate were respectively %72.81 and %3.21 from the outbreak of the epidemic in October 2019 until the end of June 2021. The table in the Appendix 1 shows the CIR-CFR-CRR values of all countries in the world every month.

4.2. Spatial Pattern and Cluster Identification

Spatial patterns of CIR, CFR, and CRR, which were in association with COVID-19, have been declared in Figures 1-3 from the onset of the epidemic in October 2019 until the end of June 2021, and then the zoning of the disease was drawn in GIS environment.

Spatial distribution of cumulative incidence rate (CIR) of COVID-19 from October 2019 to June 2021 worldwide
Spatial distribution of cumulative incidence rate (CIR) of COVID-19 from October 2019 to June 2021 worldwide
Spatial distribution of case recovery rate (CRR) of COVID-19 from October 2019 to June 2021 worldwide
Spatial distribution of case recovery rate (CRR) of COVID-19 from October 2019 to June 2021 worldwide
Case fatality rate (CFR) of COVID-19 from October 2019 to June 2021 worldwide
Case fatality rate (CFR) of COVID-19 from October 2019 to June 2021 worldwide

Overall, given the CIR, some parts of North America and the next level, the entire continent of South America, whole parts of Europe apart from Eastern Europe, and a few areas of Eastern and Southern Asia reported the highest CIR, respectively. On the other hand, the lowest CIR is recorded in the entire continent of Australia and Africa except for South Africa and a few northern regions of this continent, some of the continents of Asia and Northern America. The list of countries with the highest CIR consists of Seychelles (CIR = 21614.284), Montenegro (CIR = 20506.762), Andorra (CIR = 19615.608), Gibraltar (CIR = 16330.770), Bahrain (CIR = 16138.812) whereas Micronesia (CIR = 0.182), Vanuatu (CIR = 1.302), Samoa (CIR = 1.511), Western Sahara (CIR = 1.674), Tanzania (CIR = 2.288) are five cities with the least CIR (Figure 1). Considering countries that are more affected by the epidemic and reported the highest total case, CIR could be mentioned as follows: USA (CIR = 13151.648), UK (CIR = 11144.975), India (CIR = 11144.975), Brazil (CIR = 10024.574) and Russia (CIR = 5054.579911). As for Figure 1, the USA has first place, India, Brazil, and the UK is in second place, and Russia is in third place in CIR.

Regarding Figure 2, the highest case recovery rate was relevant to Europe and West America, while areas of Middle and Eastern Africa and the entire continent of Australia, and some parts of northern America in the next level had the least CRR. By and large, CRR has been reported to be above 77% in about 80% of countries. For example, CRR was 100% in countries with less than 50 cases. Nevertheless, CRR was reported to be more than 90% even in countries in which the incidence rate passed 200,000. Having said that, CRR was too low in some countries, such as Martinique, with the lowest CRR worldwide and approximately 41 thousand incidences. New Caledonia had the slightest CRR after Martinique, with 0.99% of CRR and about 6 thousand incidences. Other countries with less than 5% of CRR include Guadeloupe and Burundi. CRR values for countries with the highest reported total cases are as follows: USA (CRR = 75.882), India (CRR = 97.766), Brazil (CRR = 95.360), UK (CRR = 80.688), and Russia (CRR = 89.130).

In view of CFR, large areas of Africa, Asia, and South America showed the highest rate, while some areas of North America, Europe, and Australia showed the least. The highest CFR at the time of study is assigned to Vanuatu, which had one death from 4 infected cases and Yemen had the second place with 9 thousand incidences (CFR = 18.88%). The third place of CFR pertains to Peru, with an incidence rate of more than 2 million (CFR = 9.175%); in addition, Mexico and Sudan had the fourth and fifth place with approximately 7.5% of CFR, respectively. The lowest CFR respectively pertains to Laos (CFR = 0.075%), Singapore (CFR = 0.084%), Bhutan (CFR = 0.115%), Faroe Is (CFR = 0.182%), and Burundi (CFR = 0.232%) (Figure 3).

4.3. Anselin's Local Moran’s I Results

The results of Global Moran’s I analysis show that the value of this index was measured at 0.159, which indicates the intensity of autocorrelation and the very high tendency of COVID-19 disease for clustering. The obtained Z score = 12.11 and P-value < 0.001 demonstrate that the disease spatial autocorrelation rate is significant even at the confidence level of 0.99. The figure in the Appendix 2 provides information about the gained results of analyzed Global Moran’s I.

4.4. Optimize Hot Spot Analysis Results

The optimal value of bandwidth is estimated to be 2411.332 kilometers based on the OHSA method. This indicator identified hot spots in countries that are mostly located in the south and west of Europe and west of Asia. These countries include Belgium, Israel, France, Netherlands, United Kingdom, Spain, [a1] Portugal, Turkey, Greece, Poland, and Tunisia in Africa (99% of confidence level), Armenia, Azerbaijan, Georgia, Argentina, Syria, Sweden (95% of confidence level) and Iraq, Jordan, Libya, Morocco, Paraguay, Uruguay, Finland (90% of confidence level).

Moreover, cold spots are recorded in countries of Western Africa and Middle Africa such as Zaira, Congo, Gabon, Chad, Cameroon, Nigeria (99% of confidence level), and other countries like Mauritania, Mali, Niger, Sudan, Kenya, Angola, Zambia (95% of confidence level) and New Guinea, Laos, Vietnam, Cambodia, Thailand, Malaysia, Brunei, Taiwan, Botswana (90% of confidence level). With an overall look at the zoning of clusters, a majority of hot spots are located in Europe, southern and western parts of it in particular, central and western parts of Asia, and to a lesser extent, some countries of South America and North Africa. On the other hand, cold spots are located in East, West, and Central Africa and some parts of South East Asia (Figure 4).

As regards the results of Anselin's local Moran’s I, the High-High cluster consists of some countries such as Brazil, French Guiana, Colombia, Uruguay, and Argentine in South America and Sweden, Ireland, United Kingdom, France, Spain, Turkey, Tunisia, Poland, Belarus, Estonia, Latvia Ukraine, Switzerland, Italy, Romania, Bulgaria, Greece, Georgia, Armenia in southern and western Europe and central and west Asia. It declares the spatial clusters of high values and high risk of incidence in these regions. Indeed, the incidence rate in these countries is more than the global average value; in addition, they are surrounded by countries with a high outbreak.

On the other hand, Australia was identified as a Low-Low cluster which reveals the spatial cluster of low values and low risk of this disease in these regions. Firstly, the incidence rate is less than the global average. Secondly, they are surrounded by countries with low outbreak (Figure 5).

Hot spots (clusters of high values) and cold spots (clusters of low values) of COVID-19 in the world using optimized hot spot analysis (OHSA)
Hot spots (clusters of high values) and cold spots (clusters of low values) of COVID-19 in the world using optimized hot spot analysis (OHSA)
Hot spots and cold spots of COVID-19 in the world using Anselin local Moran's I index
Hot spots and cold spots of COVID-19 in the world using Anselin local Moran's I index

4.5. Time Trend Results

Time trend results of new cases, case recovery rate, and case fatality rate were analyzed by joinpoint regression on a monthly basis and are shown in Figure 6 and Table 1.

Table 1.

More Details of Temporal Analysis with Joinpoint Regression Models Fitted to New Cases, Case Recovery Rate, and Case Fatality Rate of COVID-19

Number of Graph-relatedThe Time Frame of the TrendAPCP-ValueAAPCP-Value
Point Estimate95% CIPoint Estimate95% CI
7(A)10.2019 - 12.2019861.63370.6, 1865.2< 0.00144.427.3,63.8< 0.001
12.2019 - 3.202073.25-15.2,254.00.1
3.2020 - 8.202026.100.6,58.1< 0.001
8.2020 - 6.20210.10-5.3,5.81.0
7(B)10.2019 - 12.2019-30.08-35.5, -24.2< 0.0013.21.5,4.9< 0.001
12.2019 - 3.202039.0628.3,50.7< 0.001
3.2020 - 6.20205.16-3.0,14.00.2
6.2020 - 6.20211.6912,22< 0.001
7(C)10.2019 - 1.202049.6345.1,54.3< 0.001-0.3-1.3,0.80.6
1.2020 - 5.2020-16.10-18.6, -13.5< 0.001
5.2020 - 8.2020-12.73-17.9,0.72< 0.001
8.2020 - 6.2021-1.52-2.0, -1.1< 0.001

Figure 6A and Table 1 illustrate the time trend of new cases in each month and year. As for the gained results, differentiation of disease incidence had been on the rise between October 2019 (10.2019) and June 2021 (6.2021), and the average (AAPC) increased by 44.4% per year (P < 0.001). Although a joinpoint was observed in 12.2019 and new cases changed substantially in this time period, the frequency of new cases increased dramatically between 10.2019 and 12.2019 (APC = 861.63%), notwithstanding the severity of disease incidence decreased from 12.2019 to 3.2020 (APC = 0.1%). Furthermore, there was a new significant joinpoint in 8.2020, although the incidence of new cases of the disease has reached stability and is increasing at the same rate in the period between 6.2021 and 8.2020 (APC = 0.1%).

Temporal analysis with joinpoint regression models fitted to new cases of COVID-19 in the world during COVID-19 During Oct 2019 to Jun 2021. (A) Temporal analysis with joinpoint regression models fitted to the case recovery rate of COVID-19 in the world during during Oct 2019 to Jun 2021. (B) Temporal analysis with joinpoint regression models fitted to case fatality rate of COVID-19 in the world during Oct 2019 to Jun 2021. (C) *Indicates that the annual percent change (APC) is significantly different from zero at the alpha = 0.05 level
Temporal analysis with joinpoint regression models fitted to new cases of COVID-19 in the world during COVID-19 During Oct 2019 to Jun 2021. (A) Temporal analysis with joinpoint regression models fitted to the case recovery rate of COVID-19 in the world during during Oct 2019 to Jun 2021. (B) Temporal analysis with joinpoint regression models fitted to case fatality rate of COVID-19 in the world during Oct 2019 to Jun 2021. (C) *Indicates that the annual percent change (APC) is significantly different from zero at the alpha = 0.05 level

The time trend of the case recovery rate in each month and year is demonstrated in Table 1 and Figure 6C. CRR had risen substantially from the onset of the outbreak until June 2021 (6.2021), and the average (AAPC) increased by 3.2% per year (P < 0.001). But a joinpoint was observed in 12.2019, and the CRR shifted significantly in this time period, so this progression experienced a considerable drop between 10.2019 and 12.2019 (APC = -30.08%), whereas the APC increased from 12.2019 to 3.2020 (APC = 5.16%). Thus, there was another joinpoint in 6.2020, albeit CRR increasing progression seems to be more stable between 6.2020 and 6.2021. It soared dramatically, but even so, the velocity of this increase has plunged down (APC = -1.69%).

Table 1 and Figure 6C declare the time trend of the case fatality rate in each month and year. CFR has decreased since the onset of the epidemic, and this falling progression continued until June 2021 (6.2021). It averagely dropped by -0.3% per year (AAPC), which was not considerable (P > 0.05). Although a joinpoint was observed in 1.2020 and CFR varied substantially as a consequence of this joinpoint, this indicator had a significant rise between 10.2019 and 1.2020 (APC = 49.63%). Despite that, it fell considerably between 1.2020 and 5.2020 (APC = -16.10%) and continued decreasing with a lower speed in the time period between 5.2020 and 8.2020 (APC = -12.7%). Hence the second joinpoint occurred on 5.2020, and ultimately, in the period of 8.2020 - 6.2021, it significantly reached its lowest growth rate (APC = -1.5%). The third joinpoint arose on 8.2020.

5. Discussion

From the pandemic's beginning in October 2019 until June 2021, the average case fatality rate was estimated at 0.031. Although CFR only increased in the first four months of the pandemic, it decreased from May 2020 until June 2021. This falling process may be due to countries' access to diagnostic tests, which has increased recently. This theory is consistent with a study in Illinois which states that raw CFR decreases in all regions and age groups during the development of diagnostic tests (25). This study also showed an increase in CRR. There has been no specific medicine for this disease, even though various combinations of drugs and therapies are utilized, most of which are supportive and symptomatic. As for the pneumonia-like symptoms of COVID-19, early diagnosis and controlling these symptoms results in valuable findings and a better-soaring improvement rate (13, 15, 26).

The results of this study revealed that the highest CIR is recorded in America and Europe. Perhaps the reason is that more diagnostic tests are taken in these developed areas (27).

Global Moran's I = 0.159 and obtained Z score and P-value showed that spatial autocorrelation rate and tendency to spatial clustering of COVID-19 disease is very high. This finding develops the assumption of the possible association of environmental risk factors with COVID-19, such as temperature, humidity, hours of sunshine, topographic and built environment, etc. (28). Several studies have shown that temperature and humidity are the main factors that affect SARS-CoV-2. Studies are compromised on the inverse correlation between temperature and humidity in global SARS-CoV-2 transmission. COVID-19 has sessional variations; however, vague results about the effect of humidity on SARS-CoV-2 transmission were obtained with reference to other meteorological factors such as air pollution, wind speed, etc., analyses performed in local areas (29).

We can mention the results of a study that was held in 16 African countries in view of recording cold spots in Africa. It stated that higher temperatures and humidity decrease the new daily cases and fatality rate of the COVID-19 outbreak. Therefore, these factors minimize the daily incidence and mortality from COVID-19 epidemics in African countries (30).

With an overall look at the zoning of clusters, it can be concluded that most hot spots are located in Europe, Asia, and to a lesser extent, some countries in America and Africa. On the other hand, cold spots are located in Africa and some South East Asia. Finding hotspots in developed areas could be associated with high population density, a higher proportion of the elderly population (31), high population movement via road, rail (9). Despite that, recording hot spots in regions with lower socioeconomic status, such as some parts of Asia, indicates poor access to health facilities and a higher risk for disease incidence as a consequence (15, 32).

Another result of the present study is a significantly increasing progression of disease incidence between Oct 2019 (10.2019) and June 2021 (6.2021). It had the highest growth rate in the first months through the last ten months of evaluation demonstrated the lowest growth rate. The growth rate was some have stable in this period. The high growth rate at the beginning of the pandemic outbreak could be due to a dramatic increase in patients in Wuhan, China, rapid transmission to other countries, and failure to discover suitable methods to diagnose and treat the disease (33).

The drop in the speed of disease rate may have pertained to the observation of health protocols, personal hygiene, and the use of protective equipment such as face mask and shield, gown, gloves, glasses, and N95 respirator, especially in hospital environments, proper ventilation, hand washing with soap and water or alcoholic disinfectant after contact with any contaminated surface (34), a social distancing that has affected the pattern and rate of population movement, purposeful implementation of local quarantines (35), reduction of gatherings, social events and group meetings (34), use of efficient diagnostic tools such as rapid and large-scale tests which play an essential role in patient management and reducing the rate of a disease outbreak (36) and most importantly, vaccination which is essential for infection prevention, controlling disease prevalence and reducing its recurrence rate, etc. (37).

The average case recovery rate had been plummeting in the first three months of the epidemic outbreak significantly and, after that, soared until the end of the period of study. Moreover, CFR had been increasing for the first five months of the epidemic, though it started a decreasing progression which continued until June 2021 (6.2021). This increasing process of recovery rate and reduction of CFR may be due to the strengthening of the immune system and promotion of self-care by striving for a healthy lifestyle and consuming adequate nutrients (38, 39), paying efforts to accelerate the development of treatment for COVID-19 by reusing antiviral drugs which are previously produced and approved by the FDA or currently assessed for other viral infections and combining these drugs with new technologies. These treatments aim to reduce the SARS-CoV-2 infection burden in patients, decreasing the severe symptoms and slowing down the COVID-19 improvement (40-42). For the most part, timely diagnosis and early treatment of infected cases are the first steps to winning this battle (13).

5.1. Conclusions

As regards the specific spatial pattern of fatality and recovery rate of COVID-19, it seems essential to consider spatial conditions and environmental factors which are related to the incidence and fatality of COVID-19 in different regions, as do the necessity of upgrading the care system in high-risk areas, in order to have better management and control of the pandemic and optimal function in early diagnosis, proper treatment, and high vaccination coverage.

Acknowledgements

References

  • 1.

    Rolling updates on coronavirus disease (COVID-19). World Health Organization; 2020, [cited 4/24/2020]. Available from: https://www.who.int/emergencies/diseases/novel-coronavirus-2019/events-as-they-happen.

  • 2.

    Rahimi Pordanjani S, Hasanpour A, Askarpour H, Bastam D, Rafiee M, Khazaei Z, et al. Aspects of Epidemiology, Pathology, Virology, Immunology, Transmission, Prevention, Prognosis, Diagnosis, and Treatment of COVID-19 Pandemic: A Narrative Review. Int J Prev Med. 2021;12:38. [PubMed ID: 34249287]. [PubMed Central ID: PMC8218815]. https://doi.org/10.4103/ijpvm.IJPVM_469_20.

  • 3.

    Franch-Pardo I, Napoletano BM, Rosete-Verges F, Billa L. Spatial analysis and GIS in the study of COVID-19. A review. Sci Total Environ. 2020;739:140033. [PubMed ID: 32534320]. [PubMed Central ID: PMC7832930]. https://doi.org/10.1016/j.scitotenv.2020.140033.

  • 4.

    Pattison WD. The Four Traditions of Geography. J Geogr. 1964;63(5):211-6. https://doi.org/10.1080/00221346408985265.

  • 5.

    Kavousi A, Salehi M. [Spatial analysis of gastrointestinal Cancer incidence rate in Iran using Poisson Kriging]. Iran J Off Stat Stud. 2012;23(1):129-40. Persian.

  • 6.

    Pordanjani SR, Kavousi A, Mirbagheri B, Shahsavani A, Etemad K. Spatial analysis and geoclimatic factors associated with the incidence of acute lymphoblastic leukemia in Iran during 2006-2014: An environmental epidemiological study. Environ Res. 2021;202:111662. [PubMed ID: 34273372]. https://doi.org/10.1016/j.envres.2021.111662.

  • 7.

    Olfatifar M, Hosseini SM, Shokri P, Khodakarim S, Khadembashi N, Rahimi Pordanjani S. How to improve the human brucellosis surveillance system in Kurdistan Province, Iran: reduce the delay in the diagnosis time. Epidemiol Health. 2020;42. e2020058. [PubMed ID: 32777881]. [PubMed Central ID: PMC7871156]. https://doi.org/10.4178/epih.e2020058.

  • 8.

    Azizpour F, Riahi V. Spatial analysis of corona outbreak in rural area of Damavand city. J Spat Anal Environ Hazarts. 2021;7(4):21-40.

  • 9.

    Rex FE, Borges CAS, Kafer PS. Spatial analysis of the COVID-19 distribution pattern in Sao Paulo State, Brazil. Cien Saude Colet. 2020;25(9):3377-84. [PubMed ID: 32876254]. https://doi.org/10.1590/1413-81232020259.17082020.

  • 10.

    Rahimi Pordanjani S, Kavousi A, Mirbagheri B, Shahsavani A, Etemad K. Geographical Pathology of Acute Lymphoblastic Leukemia in Iran with Evaluation of Incidence Trends of This Disease Using Joinpoint Regression Analysis. Arch Iran Med. 2021;24(3):224-32. [PubMed ID: 33878881]. https://doi.org/10.34172/aim.2021.34.

  • 11.

    Kulldorff M, Athas WF, Feurer EJ, Miller BA, Key CR. Evaluating cluster alarms: a space-time scan statistic and brain cancer in Los Alamos, New Mexico. Am J Public Health. 1998;88(9):1377-80. [PubMed ID: 9736881]. [PubMed Central ID: PMC1509064]. https://doi.org/10.2105/ajph.88.9.1377.

  • 12.

    Jesus JG, Sacchi C, Candido DDS, Claro IM, Sales FCS, Manuli ER, et al. Importation and early local transmission of COVID-19 in Brazil, 2020. Rev Inst Med Trop Sao Paulo. 2020;62. e30. [PubMed ID: 32401959]. [PubMed Central ID: PMC7232955]. https://doi.org/10.1590/s1678-9946202062030.

  • 13.

    Rafiee M, Parsaei F, Rahimi Pordanjani S, Amiri V, Sabour S. A Review on Applicable and Available Paraclinical Methods for Diagnosis of Coronavirus Disease-19. Arch Iran Med. 2020;23(11):794-800. [PubMed ID: 33220699]. https://doi.org/10.34172/aim.2020.106.

  • 14.

    Ramirez-Aldana R, Gomez-Verjan JC, Bello-Chavolla OY. Spatial analysis of COVID-19 spread in Iran: Insights into geographical and structural transmission determinants at a province level. PLoS Negl Trop Dis. 2020;14(11). e0008875. [PubMed ID: 33206644]. [PubMed Central ID: PMC7710062]. https://doi.org/10.1371/journal.pntd.0008875.

  • 15.

    Khazaei Z, Mazaheri E, Hasanpour-Dehkordi A, Rahimi Pordanjani S, Naghibzadeh-Tahami A, Naemi H, et al. COVID-19 Pandemic in the World and its Relation to Human Development Index: A Global Study. Arch Clin Infect Dis. 2020;15(5). https://doi.org/10.5812/archcid.103093.

  • 16.

    Gross B, Zheng Z, Liu S, Chen X, Sela A, Li J, et al. Spatio-temporal propagation of COVID-19 pandemics. EPL (Europhysics Letters). 2020;131(5):58003. https://doi.org/10.1209/0295-5075/131/58003.

  • 17.

    Rothman KJ, Greenland S, Lash TL. Modern epidemiology. Philadelphia, USA: Wolters Kluwer Health/Lippincott Williams & Wilkins; 2008.

  • 18.

    COVID-19. World Bank; 2021. Available from: https://data.worldbank.org/.

  • 19.

    COVID-19 Coronavirus Pandemic. Worldometers; 2021. Available from: https://www.worldometers.info/coronavirus/.

  • 20.

    Lesson 3: Measures of Risk. Centers for Disease Control and Prevention; 2012. Available from: https://www.cdc.gov/csels/dsepd/ss1978/lesson3/section2.html.

  • 21.

    Krug A. Cumulative Incidence. Encyclopedia Britannica; 2016, [cited 8/1/2016]. Available from: https://www.britannica.com/science/cumulative-incidence.

  • 22.

    Battegay M, Kuehl R, Tschudin-Sutter S, Hirsch HH, Widmer AF, Neher RA. 2019-novel Coronavirus (2019-nCoV): estimating the case fatality rate - a word of caution. Swiss Med Wkly. 2020;150:w20203. [PubMed ID: 32031234]. https://doi.org/10.4414/smw.2020.20203.

  • 23.

    Pordanjani SR, Kavousi A, Mirbagheri B, Shahsavani A, Etemad K. Identification of high-risk and low-risk clusters and estimation of the relative risk of acute lymphoblastic leukemia in provinces of Iran during 2006-2014 period: A geo-epidemiological study. J Res Med Sci. 2021;26:18. [PubMed ID: 34084197]. [PubMed Central ID: PMC8106411]. https://doi.org/10.4103/jrms.JRMS_662_20.

  • 24.

    Rahimi Pordanjani S, Kavousi A, Mirbagheri B, Shahsavani A, Etemad K. Temporal trend and spatial distribution of acute lymphoblastic leukemia in Iranian children during 2006-2014: a mixed ecological study. Epidemiol Health. 2020;42. e2020057. [PubMed ID: 32777882]. [PubMed Central ID: PMC7871159]. https://doi.org/10.4178/epih.e2020057.

  • 25.

    Holden TM, Richardson RAK, Arevalo P, Duffus WA, Runge M, Whitney E, et al. Geographic and demographic heterogeneity of SARS-CoV-2 diagnostic testing in Illinois, USA, March to December 2020. BMC Public Health. 2021;21(1):1-13. [PubMed ID: 34107947]. [PubMed Central ID: PMC8189821]. https://doi.org/10.1186/s12889-021-11177-x.

  • 26.

    Di Castelnuovo A, Bonaccio M, Costanzo S, Gialluisi A, Antinori A, Berselli N, et al. Common cardiovascular risk factors and in-hospital mortality in 3,894 patients with COVID-19: survival analysis and machine learning-based findings from the multicentre Italian CORIST Study. Nutr Metab Cardiovasc Dis. 2020;30(11):1899-913. [PubMed ID: 32912793]. [PubMed Central ID: PMC7833278]. https://doi.org/10.1016/j.numecd.2020.07.031.

  • 27.

    Fani M, Teimoori A, Ghafari S. Comparison of the COVID-2019 (SARS-CoV-2) pathogenesis with SARS-CoV and MERS-CoV infections. Future Virology. 2020;15(5):317-23. https://doi.org/10.2217/fvl-2020-0050.

  • 28.

    Fatima M, O'Keefe KJ, Wei W, Arshad S, Gruebner O. Geospatial Analysis of COVID-19: A Scoping Review. Int J Environ Res Public Health. 2021;18(5). [PubMed ID: 33673545]. [PubMed Central ID: PMC7956835]. https://doi.org/10.3390/ijerph18052336.

  • 29.

    Byun WS, Heo SW, Jo G, Kim JW, Kim S, Lee S, et al. Is coronavirus disease (COVID-19) seasonal? A critical analysis of empirical and epidemiological studies at global and local scales. Environ Res. 2021;196:110972. [PubMed ID: 33705770]. [PubMed Central ID: PMC7941024]. https://doi.org/10.1016/j.envres.2021.110972.

  • 30.

    Meo SA, Abukhalaf AA, Alomar AA, Aljudi TW, Bajri HM, Sami W, et al. Impact of weather conditions on incidence and mortality of COVID-19 pandemic in Africa. Eur Rev Med Pharmacol Sci. 2020;24(18):9753-9. [PubMed ID: 33015822]. https://doi.org/10.26355/eurrev_202009_23069.

  • 31.

    Zhang CH, Schwartz GG. Spatial Disparities in Coronavirus Incidence and Mortality in the United States: An Ecological Analysis as of May 2020. J Rural Health. 2020;36(3):433-45. [PubMed ID: 32543763]. [PubMed Central ID: PMC7323165]. https://doi.org/10.1111/jrh.12476.

  • 32.

    Macharia PM, Joseph NK, Okiro EA. A vulnerability index for COVID-19: spatial analysis at the subnational level in Kenya. BMJ Glob Health. 2020;5(8). [PubMed ID: 32839197]. [PubMed Central ID: PMC7447114]. https://doi.org/10.1136/bmjgh-2020-003014.

  • 33.

    Al Hasan SM, Saulam J, Kanda K, Hirao T. The novel coronavirus disease (COVID-19) outbreak trends in mainland China: a joinpoint regression analysis of the outbreak data from January 10 to February 11, 2020. Bull World Health Organ. 2020;17. https://doi.org/10.2471/blt.20.253153.

  • 34.

    Shoaib MH, Ahmed FR, Sikandar M, Yousuf RI, Saleem MT. A Journey From SARS-CoV-2 to COVID-19 and Beyond: A Comprehensive Insight of Epidemiology, Diagnosis, Pathogenesis, and Overview of the Progress into Its Therapeutic Management. Front Pharmacol. 2021;12:576448. [PubMed ID: 33732150]. [PubMed Central ID: PMC7957225]. https://doi.org/10.3389/fphar.2021.576448.

  • 35.

    Sartorius B, Lawson AB, Pullan RL. Author Correction: Modelling and predicting the spatio‑temporal spread of COVID‑19, associated deaths and impact of key risk factors in England. Sci Rep. 2021;11(1):1-11. [PubMed ID: 34465862]. [PubMed Central ID: PMC8406384]. https://doi.org/10.1038/s41598-021-97282-8.

  • 36.

    Safiabadi Tali SH, LeBlanc JJ, Sadiq Z, Oyewunmi OD, Camargo C, Nikpour B, et al. Tools and Techniques for Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2)/COVID-19 Detection. Clin Microbiol Rev. 2021;34(3). [PubMed ID: 33980687]. [PubMed Central ID: PMC8142517]. https://doi.org/10.1128/CMR.00228-20.

  • 37.

    Liu C, Zhou Q, Li Y, Garner LV, Watkins SP, Carter LJ, et al. Research and Development on Therapeutic Agents and Vaccines for COVID-19 and Related Human Coronavirus Diseases. ACS Cent Sci. 2020;6(3):315-31. [PubMed ID: 32226821]. [PubMed Central ID: PMC7094090]. https://doi.org/10.1021/acscentsci.0c00272.

  • 38.

    Simpson RJ, Kunz H, Agha N, Graff R. Exercise and the Regulation of Immune Functions. Prog Mol Biol Transl Sci. 2015;135:355-80. [PubMed ID: 26477922]. https://doi.org/10.1016/bs.pmbts.2015.08.001.

  • 39.

    Yoshikawa TT, High KP. Nutritional strategies to boost immunity and prevent infection in elderly individuals. Clin Infect Dis. 2001;33(11):1892-900. [PubMed ID: 11692301]. https://doi.org/10.1086/324509.

  • 40.

    Saha RP, Sharma AR, Singh MK, Samanta S, Bhakta S, Mandal S, et al. Repurposing Drugs, Ongoing Vaccine, and New Therapeutic Development Initiatives Against COVID-19. Front Pharmacol. 2020;11:1258. [PubMed ID: 32973505]. [PubMed Central ID: PMC7466451]. https://doi.org/10.3389/fphar.2020.01258.

  • 41.

    Saul S, Einav S. Old Drugs for a New Virus: Repurposed Approaches for Combating COVID-19. ACS Infect Dis. 2020;6(9):2304-18. [PubMed ID: 32687696]. [PubMed Central ID: PMC7416635]. https://doi.org/10.1021/acsinfecdis.0c00343.

  • 42.

    Binesh M, Pahlevanian A, Rahimi Pordanjani S, Ahmadizadeh Z. Relationship Between COVID-19-related Factors and Self-management Behaviors in People with Type-2 Diabetes: A Cross-sectional Study. Middle East J Rehabil Health Stud. 2022;9(3). https://doi.org/10.5812/mejrh-129017.