Spatiotemporal patterns and progression of the Delta variant of COVID-19 and their health intervention linkages in Southeast Asia

Luo Wei, Liu Zhaoyin, Zhou Yuxuan, Zhao Yumin, Li Yunyue Elita, Arif Masrur, Manzhu Yu
doi: https://doi.org/10.1101/2021.12.20.21268140

Abstract

The global pandemic of COVID-19 presented an unprecedented challenge to all countries in the world, among which Southeast Asia (SEA) countries managed to maintain and mitigate the first wave of COVID-19 in 2020. However, these countries were caught in the crisis after the Delta variant was introduced to SEA, though many countries had immediately implemented non-pharmaceutical intervention (NPI) measures along with vaccination in order to contain the disease spread. To investigate the potential linkages between epidemic dynamics and public health interventions, we adopted a prospective space-time scan method to conduct spatiotemporal analysis at the district level in the seven selected countries in SEA from June 2021 to October 2021. Results reveal the spatial and temporal propagation and progression of COVID-19 risks relative to public health measures implemented by different countries. Our research benefits continuous improvements of public health strategies in preventing and containing this pandemic.

1. Introduction

Coronavirus disease 2019 (COVID-19) is a global epidemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). SARS-CoV-2 is highly contagious as it easily spreads across human, animals and the environment via contact, droplet, airborne, fomite and other transmission modes [1]. In December 2019, COVID-19 was first identified in Wuhan – a transportation and communication hub located in central China, and rapidly spread to surrounding regions in China. Although a certain level of precautions and control regulation were taken by governments, international cases managed to emerge along with significant outbreaks worldwide. As of December 12, 2021, the disease had infected over two hundred million people across all continents except from Antarctica [2].

With the outbreak of COVID-19, countries in SEA met unprecedented challenges. Health and social care system, as well as tourism, trade in essential goods and services of SEA countries were under pressure [3]. Hence, SEA countries with export-oriented economies like Cambodia, Vietnam, Singapore etc. suffered varying degrees of damage due to COVID-19 [3, 4]. Many SEA countries had made efforts to recover the economy during the late 2020 and early 2021 by progressively easing lockdown and strengthening export orders, but they were continuously hammered by the new wave of COVID-19 delta variant [5]. The Delta variant is estimated to be twice or four times more transmissible than its original virus with a Reproductive number (R0) as nearly 7 [6]. Since April 2021, the Delta variant of SARS-CoV-2 has caused an exponential increase of new cases in SEA that has become an emerging hotspot region of COVID-19 [7]. The new outbreak caused by delta variant heavily disrupted their national and international business due to interruption of supply chain of many productions (e.g., garment, auto parts, semiconductors) which had great impacts on the manufacture industries in SEA, especially for those relying on them (e.g., Vietnam. Malaysia, Thailand, Indonesia, and Philippines) [5]. Due to the economic pressure caused by the epidemic, a growing number of Southeast Asian countries were planning to live with the virus, along with adjusting their public health intervention policies. [8]. Under this circumstance, monitoring outbreaks and identifying the space-time clusters of infection have become significant for a coordinated response to the epidemic in SEA.

Spatio-temporal analysis has been widely used in research of COVID-19 propagation to illustrate the characteristics and mechanism of COVID-19 spatial propagation, which provided public health authorities with important information to help mitigate the situation [9-11]. Among diverse spatio-temporal methods, space-time scan is one of the most popular methods adopted to explore spatiotemporal clusters in a particular region around the world, such as Mainland China [12], the United States [13], Mexico [14], Spain [15], Malaysia [16], Bangladesh [17] etc. In SEA, previous studies have applied this analysis to investigate the first wave of COVID-19 cases in SEA [18, 19]. However, those studies mainly focused on each individual country without exploring propagation patterns and progression characteristics with collective public health interventions at the regional scale. Previous research has shown that regional-wide coordination could interrupt transmission of COVID-19 in an effective way [10, 17, 20-22]. In order to contain the emerging epidemic of COVID-19 and minimize its risk, countries in SEA have deployed various preventive and containment measures such as lockdowns, distancing restrictions, and mandatory tracking and trace method [23, 24].

Delta variant transmission had shown significant spatial heterogeneity in SEA because different countries adopted different interventions along with the spread. Hence, this paper aims to identify the space-time clusters of outbreaks of COVID-19 caused by SARS-CoV-2 delta variant in SEA. We utilized district-scale daily confirmed cases of seven SEA countries from June 2021 to October 2021 to identify the active and emerging clusters of the disease outbreaks and summarized the relative policies to investigate the potential linkage between government interventions and the transmission. Our work contributes to regional surveillance of COVID-19 progression in SEA and provides essential information of COVID-19 propagation from space and time perspective to public health authorities, which is beneficial to timely policy making according to the dynamic COVID-19 situation.

2. Materials and Methods

2.1. Study areas and relevant interventions

Our study focused on seven countries in SEA, namely Indonesia, Malaysia, the Philippines, Singapore, Thailand, Vietnam and Brunei, because they disclose relevant data at a district level. During the second COVID-19 outbreak caused by the Delta variant, the dynamic of interventions implemented by different countries along the way may cause the fluctuation of transmission. For example, Thailand, Singapore, Malaysia, Philippines, and Vietnam began to relax their restrictions around August, which may cause significant changes in the epidemic patterns after then (Table 1). The diverse policies will help explain the progression and transmission of the Delta variant of COVID-19 in the following analysis. The interventions are aggregated by Center for Strategic & International Studies (CSIS) [25].

View this table:
Table 1

Major public health interventions in SEA.

2.2. COVID-19 Daily cases and population

We obtained or extracted data of COVID-19 confirmed cases from the official websites of public health authorities in the seven countries and Johns Hopkins University’s Center for Systems Science and Engineering GIS dashboard (Table 2). We adopted data from 1 June 2021 to 31 October 2021, which is the approximate date of the second COVID-19 outbreak of these seven countries in SEA (Figure 1). We aggregated the data at the first administrative level except for those in Singapore and Brunei which were aggregated at country level considering the similar magnitude of area and population in each analytic unit. We obtained or extracted population data from statistical reports and yearbooks from those countries (Table 2).

View this table:
Table 2.

Collected Data and Their Sources.

Figure 1
Figure 1

Daily new cases of COVID-19 in SEA (from January 1 to October 31)

Figure 2
Figure 2

Spatial distribution of space-time clusters of COVID-19 from June 1 to August 31.

2.3. Space-Time scan statistical analysis

To explore emerging and active space-time clusters of COVID-19 cases in SEA, we conducted SaTScanTM (version 9.6) [26] to run prospective space-time scan statistic, which is normally used in geographical clusters of epidemic detection [27, 28]. Using space-time scan statistic, we can identify and map significant clusters of the Delta variant in SEA considering the uneven distribution of population size. The space-time scan statistic adopted a cylinder to detect potential space-time clusters in SEA, which can cover each possible location, size, and time period. For each cylinder, the base represents space, the height represents time, and the center represents the centroids of study units throughout SEA. The size of the cylindrical window is expanded until reaching a specific maximum spatial and temporal upper bound which is set to 10% of the population risk and 50% of the study period respectively in this study. Besides, we defined the minimum duration of each cluster to two days for surveillance of the continuously existing clusters. Moreover, the minimum number of cases in each cluster is set to 3 in order to ensure that there must be at least three cases in each cluster.

We assumed that the COVID-19 cases follow a Poisson distribution according to the population of study units in SEA. The null hypothesis H0 indicates that the model reflects infection of COVID-19 having a constant intensity μ within or outside the cylinder, which is proportional to at-risk population. Alternative hypothesis HA indicates that the observed cases are more than expected cases, which reflects an increased risk within a cylinder. Expected cases are calculated by Equation (1) [13, 29]: Embedded Image where p represents the population within a study unit, C represents the total COVID-19 cases in our study area (i.e., seven countries in SEA), and P represents the total estimated population within our study area.

A maximum likelihood ratio test is used to evaluate the null and alternative hypotheses. It can identify scanning windows with an elevated risk for COVID-19, which is defined by equation (2) [30, 31]: Embedded Image where represents the likelihood function for cylinder Z; represents the likelihood function for H0; represents the number of COVID-19 cases in a cylinder; represents the number of expected cases in cylinder Z; N represents the total number of observed cases for the seven countries in SEA across all time periods. When the likelihood ratio is greater than 1, there is an elevated risk in the cylinder, and the cylinder with the maximum likelihood ratio should be the most likely cluster.

The relative risk of COVID-19 is assumed homogeneous throughout different districts within the same cluster. To make the results more reasonable, we calculated relative risk (RR) for each study unit within a cluster to explore the spatial heterogeneity of relative risk of COVID-19 (equation (3)) [32]: Embedded Image where c represents the total number of COVID-19 in a study unit; e represents the total number of expected cases in a study unit; C represents the total number of observed cases in the seven countries of SEA. The formula indicates that RR represents the estimated risk in a study unit divided by the risk outside that unit. Specifically, if a location (cluster or study unit) has a RR of 3, the population within the location are 3 times more likely to be exposed to COVID-19 infection than its outside.

In the following sections, the results reveal significant emerging clusters of COVID-19 cases in seven countries of SEA from 1 June to 31 October. Considering that some SEA countries (e.g., Malaysia, Philippines, Singapore, Thailand) started loosening their restrictions from August along with strengthening their vaccination plan, we divided the timeline into two parts (i.e., June 1 to August 31 and June 1 to October 31) in order to identify the dynamics of clusters. Additionally, we explored the variation of RR across each district in each half-month, which is approximately equal to the incubation period of an infection [33], using a cumulative half a month prospective scanning approach. Then we compared the interventions with the discovered space-time characteristics to identify the potential linkage between political intervention and progression of the Delta variant of COVID-19.

3. Results and Discussion

3.1. Dynamics of district-level merging clusters in SEA

3.1.1. Results between June 1 and August 31, 2021

As shown in table 3, 14 significant space-time clusters were identified during June 1 and August 31 in SEA, including 7 high-risk clusters (RR > 1) and 7 low-risk clusters (RR < 1). Most of the high-risk clusters emerged between mid-July and late-August, which means the situation of COVID-19 in SEA became severe during this period. Specifically, Cluster 1 was the most likely cluster and a transnational cluster, containing 39 high-risk districts (RR>1) out of 47 districts of Malaysia and Thailand. This cluster had the highest relative risk of 5.45, which means people in this cluster were 5.45 times more likely to be exposed to COVID-19 than other regions. Similarly, cluster 2 had relative risk of 5.62 and contained 2 districts in Vietnam, namely, Binh Duong and Ho Chi Minh city. Another transnational cluster was cluster 4, with an RR of 3.50, containing 6 districts of Malaysia, Indonesia, and Brunei. Additionally, north Thailand and north Philippines also emerged as high-risk clusters from July 22 to August 31, and from August 11 to August 31 respectively. There were also 2 clusters emerging in Indonesia and they contained only 1 district, which was Jakarta, with an RR of 2.38, and Daerah Istimewa Yogyakarta (DIY), with an RR of 3.88. Three high-risk clusters in Indonesia reveal that population within a number of regions in Indonesia were more likely to be exposed to COVID-19 compared to other regions in SEA during this period. Additionally, there were 7 low-risk clusters distributing across other regions of SEA (e.g., north of Vietnam, south of Philippines, some other districts of Indonesia), which means population within these clusters were less likely to be exposed to COVID-19. Note that, cluster 10, with an RR of 0.60, contained 2 high-risk districts in Philippines (i.e., Region VII, with an RR of 1.07, and Region X, with an RR of 1.08). Fig 3 shows the distribution of each cluster. There were a number of small-scale clusters in the south of Indonesia while the largest scale cluster appeared across the south of Thailand and the north of Malaysia. From the results, the Delta variant of COVID-19 had a wider influence in Malaysia and Indonesia in the early phase, while some regions in Vietnam and Philippines had relatively high risk as well.

View this table:
Table 3

Space-time clusters of COVID-19 from June 1 to August 31.

Figure 3
Figure 3

Spatial distribution of space-time clusters of COVID-19 from June 1 to October 31

3.1.2. Results between June 1 and October 31, 2021

Eleven significant clusters were detected between June 1 and October 31, which was three less than the early phase. Among the clusters, there were only 4 high-risk clusters (RR > 1), decreasing from 7 in the previous period (Table 4). The most likely cluster was the same as one of the previous periods, which covered partial districts of Thailand and Malaysia. The RR of cluster 1, however, decreased from 5.45 to 3.91 in this period, while the number of high-risk districts increased from 39 to 45, which implies that more districts were influenced by COVID-19. Cluster 2 was also the same as one of the previous periods, containing Binh Duong and Ho Chi Minh city. Besides, cluster 7 evolved from cluster 4 in the previous period, with the exclusion of Kalimantan Timur. The RR of this cluster increased from 3.50 to 4.62, which means that people in this cluster were more likely to be infected in this period. Similarly, Cordillera Administrative Region, Region II, and Region I merged with one high-risk cluster in the north of Philippines and formed a larger cluster. The rest of the clusters were low-risk clusters with RR less than 1. Figure 3 visualizes the distribution of clusters in this period, which shows directly that some of the clusters remained between 2 continuous periods, while a number of clusters in the previous period disappeared and some new clusters appeared in this period. Especially in Indonesia, high-risk clusters in Jakarta and DIY disappeared, which emerged low-risk clusters with other districts. Overall, the space-time scan statistic results show the transmission and dispersal of the Delta variant of COVID-19 in SEA from 2 different time periods.

View this table:
Table 4

Space-time clusters of COVID-19 from June 1 to October 31.

3.2. Temporal progression of relative risk of COVID-19 in SEA

Figure 4 presents the changes of District RR of COVID-19 in SEA between two outbreak periods, 1 June to 31 August and 1 June to 31 October. Overall, this temporal change in RR manifests different space-time characteristics of the progression of COVID-19 in different districts before and after countries changed their intervention strategies in handling COVID-19 in July and August. Amidst the seven SEA countries, Indonesia was the only country that showed an overall positive trend of decreasing RR in every district, while alarming changes of increasing RR were frequently seen in other six countries especially Singapore, Philippines, Malaysia and Vietnam.

Figure 4
Figure 4

Changes in relative risk of COVID-19 (district level) between two periods (June 1 to August 31 and June 1 to October 31)

Specifically, all districts of Indonesia observed a decrease in RR to different extents. While RR in most districts of Indonesia decreased slightly between 0.05 to 0.5, the other five districts (i.e., Jakarta, Daerah Istimewa Yogyakarta, Kalimantan Utara, Kalimantan Timur, and Kepulauan Riau) manifested a rather significant decrease (≤-0.5), in which the highest difference (−1.4) between the two phases was seen in Jakarta, the capital of Indonesia. Note that Jakarta was one of the major emerging risk districts early in the second outbreak (Figure 5) and was still faced with a relatively high RR (2.79) until the end of our study period. Nevertheless, characterized as the largest city nationwide with a high-density population, Jakarta’s success in preventing exacerbation of risk impact of COVID-19 is most likely to be related to its persistent restrictions and unchanged strategy towards COVID-19. Ever since June 1, all provinces in Indonesia have implemented Micro Community Activity Restrictions (PPKM) which includes policies of mandatory work from home, guidance on online teaching and restrictions of dine-in, social gathering and inter-province as well as international traveling [34]. On the contrary, RR of all districts in the Philippines have increased between the two outbreak periods, indicating an overall deterioration in the risk impact of COVID-19. Fortunately, amidst a total of 17 districts in the Philippines, 14 of which showed minor increases (≤0.5). Other three districts, namely National Capital Region, Region II and Cordillera Administrative Region manifested an increase between 0.59 to 1.23. Meanwhile, no significant increase (>1.5) was observed, indicating that the most severe variation of RR did not occur in the Philippines. This is mainly because Philippines loosened restrictions on a wide scale, including reopening restaurants and churches, as well as replacing a larger scale of coronavirus restrictions with localized lockdowns, though it did take certain measures to intervene in the epidemic outbreak, such as the strictest level of lockdown in Metro Manila and online teaching which may have helped prevent significant increases in RR. The gradually intensified risks of districts in the Philippines should receive continuous attention.

Figure 5
Figure 5

Spatial patterns of progression of COVID-19 relative risk in SEA (June 1 to June 15, June 1 to June 30, June 1 to July 15, June 1 to July 31, June 1 to August 15, June 1 to August 31)

Other countries (Singapore, Malaysia, Brunei, Thailand, Vietnam), on the other hand, showed a variety of RR changes in different districts. Approximately half of which (39) showed increases in RR, and half decreased amid the 77 districts in Thailand. It is obvious that coastal areas in the south of the country were faced with more elevated RR than inland areas in the central, eastern and northern districts. RR in almost all districts in Vietnam changed slightly, ranging from -0.08 to 0.28 except for 3 connected cities (i.e., Dong Nai:0.63, Ho Chi Minh: 0.72 and Binh Duong:1.73). Particularly, although a series of regulations and restrictions including lockdown, curfew, social distancing measures etc., had been specifically implemented in Ho Chi Minh, the city still manifested overall worsen symptoms. This may be related to the high contagiousness of Delta variant and overloaded health-care system in Ho Chi Minh City [35]. Also, Ho Chi Minh City has a larger population base and more developed economy than other districts which resulted in more infections via contact with a relatively larger number of crowds. It has been suggested that population density and contact intensity are the main drivers for the propagation and amplification of this virus [36]. As for Malaysia, 10 out of 16 districts showed an increase in RR and took up a major proportion of the country. Although Labuan and Kuala Lumpur observed the most obvious decrease among SEA districts, it should be noted that they were also the most severely threatened districts by COVID-19 risks early in the study period (RR>10) and remain the relatively highest RR amidst SEA until the end of the study. This is probably because the Malaysian government gradually adopted more loose measures than the other SEA countries with its phase development of the National Recovery Plan since the beginning of our study period. Additionally, Singapore saw an increased RR of 1.43, which is probably related to its intermittent relaxation of social restrictions. RR in Brunei also increased, but it still remained less than 1.

Furthermore, we detected the progression of RR at half a month interval, and 10 intervals were used to illustrate the progression of the Delta variant from June to October (Figure 5, Figure 6). In the first half of June, elevated risks were identified in the middle south of Thailand, many states in Malaysia, and many districts in Philippines, while 44 out of 63 districts in Vietnam manifested RR=0 during this period (Figure 5). These reveal that those high-risk districts were more likely to be in danger of a new wave of epidemic. From then on, the Delta variant spread in SEA, and the capital areas of several countries have been severely affected in SEA (e.g., Bangkok Metropolis and surrounding areas, Kuala Lumpur and surrounding areas, Jakarta and surrounding areas, Manila capital areas, as well as Ho Chi Minh city). By the middle of August, Thailand had been influenced by the expansion of COVID-19, and most districts showed an elevated risk. On the contrary, RR in the Manila capital region in Philippines decreased, which was probably due to suspension of travel from Malaysia, tightened restrictions, and curfew extension in late July and early August. Additionally, most high-risk areas in the previous months kept severe in this time, although the relative risk of some districts slightly declined (e.g., Bangkok, Samut Sakhon, Jakarta, Riau, Sarawak, etc.). Note that, most of the countries tightened the restriction before this time except for Malaysia, which loosened their restrictions. This phenomenon also reflected the strong infection and danger of the Delta variant. In the next one and a half month (i.e., August 15 to September 30), the situation of Indonesia had been improving, which was inseparable with continuous extension of activity restrictions in Indonesia (e.g., extended PPKM for the ninth time). The north of Philippines, however, showed increased relative risk in mid-September, which may be due to the relaxing restriction of Manila capital region on August 19 and mid-September (Figure 6). In the next month, relative risk of Singapore increased from 0.44 to 1.53, which was probably due to gradually loosened restrictions. Besides, the situation across north of Malaysia and south of Thailand was getting worse as well. Up to October 31, the situation of Thailand, Malaysia, Singapore, capital areas in Vietnam, and capital areas in Philippines was still alarming, which should be noticed with further surveillance by obtaining daily case reports in the future.

Figure 6
Figure 6

Spatial patterns of progression of COVID-19 relative risk in SEA (June 1 to September 15, June 1 to September 30, June 1 to October 15, June 1 to October 31)

4. Discussion

Principal Findings

In this study, we utilized prospective space-time scan statistics to detect the emerging and existing space-time clusters of COVID-19 in SEA. We found that most districts in Malaysia and Philippines, the capital and its surrounding areas in Thailand, Vietnam, and Indonesia exhibited high risk of COVID-19 transmission in the early phase (June to August 2021) during our study period. Space-time clusters and relative risk of districts changed along with the dynamics of policies introduced by each country after August. Indonesia successfully mitigated the risk of this epidemic by implementing continuous restrictions, while a number of regions in Malaysia, Singapore, Thailand, and Philippines remained at a relatively high risk of COVID-19 transmission due to diverse degrees of relaxation. Although propagation of COVID-19 was influenced by diverse factors (e.g., social economy, population, resource allocation, etc.), continuous strict restrictions were beneficial on epidemic control, especially for those developing regions with weak public health systems and relatively low vaccination rates. To our knowledge, this is the first attempt to explore the space-time progression of delta variant outbreak in SEA, as well as summarize the potential linkage between the epidemic dynamics and diverse public health policies. Spatio-temporal statistics is an effective method to identify high-risk clusters and regions [37], which can help explore spatial temporal patterns of epidemic progression and better understand the dynamics and characteristics of the Delta variant in SEA. Furthermore, by combining the results with relevant measures that reflect the appropriateness of public health interventions, governments can also timely improve the interventions. In the meantime, health authorities can timely adjust their interventions by comparing the RR with current measures, which may reflect the appropriateness of these interventions. In short, the prospective space-time scan can continuously be used to monitor the dynamics of COVID-19 with the most updated data, and timely adjust the potential gaps in domestic public health intervention to prevent further deterioration of the pandemic.

Implications and Recommendations

Public health interventions played an important role in epidemic containment, in which social restriction policies effectively mitigate the propagation of COVID-19 [29]. Restrictions of mass gathering and travel, keeping social distance, and reducing human mobility are beneficial to control COVID-19 because these measures could reduce the probability of exposure to virus infection [38, 39]. Our study discovered the potential linkage between the dynamics of outbreaks and interventions. This indicates that while continuously implemented strict restrictions can prevent exacerbation of the pandemic, temporary or continuous relaxation may result in an acceleration of the epidemic propagation. Appropriate restriction policies are key to preventing the pandemic because high transmission of the COVID-19 variant can easily lead to worse situations. Also, if the number of community cases exceeded imported cases, border restriction would be less valuable than domestic measures. In this case, authorities should emphasize more on domestic intervention in order to reduce community transmission [40]. Nevertheless, it is well known that those intervention measures against COVID-19 come with economic costs from other perspectives (i.e., labor market etc.). When making intervention regulations, economic and social justifications for public health policies are one of the priorities that governments would consider [41]. Hence, it is one of the challenges for all countries to live with COVID-19 and weigh the balance between epidemic development and social economic loss [42]. Given the social and economic status, most countries in SEA have decided to change their strategies from the elimination of cases to living with COVID-19 since August 2021 despite the new outbreak of Delta variant [8]. A concurrent trend observed in SEA is that all countries except for Indonesia have been gradually loosening social restrictions, allowing international communications whilst boosting vaccination attempting to achieve group immunity among citizens. For instance, the Singapore Ministry of Health believes that with the assistance of a high vaccination rate, it is assured that the growth of COVID-19 cases, though cannot be eliminated, will be controlled at a certain level. By doing this, economy and social norms could be restored without causing another disease outbreak or a breakdown in the hospital system [43]. Vaccination is increasingly essential to protect the crowd from the exacerbating threat of morbidity and mortality of this Delta variant and future variant [44]. Effective vaccination is highly beneficial to avert infection, while wide-scale vaccination is proven to be able to successfully control transmission of COVID-19 so far [45-47]. However, it should be noted that vaccinations cannot be depended on entirely. The previous experiment found that effectiveness of available vaccines against the Delta variant (B.1.617.2) showed a reduction compared to previous virus variants [48], implying that current vaccination is highly likely to become ineffective in defending against many mutations in the future [49]. On the other hand, considering the fragile health systems in SEA, contact tracing, quick isolation and strict restrictions are still essential to avoid any miscalculation and prevent another outbreak [50, 51].

International coordination including vaccine allocation, food, finance, and equitable access is essential for the world to respond to the pandemic [52]. The pandemic did not only cause the breakdown of developing countries, but also developed countries, which indicates that countries should enhance transparency, share information, and strengthen cooperation to protect international health and security [53]. As for Southeast Asia, due to strong socioeconomic connectivity among all countries, it is vital for countries in SEA to adopt collaborative policies to overcome the challenges of COVID-19. Sustainable development in SEA during the pandemic era requires multi-sectoral cooperation including tourism, economy, health etc., which can help effectively coordinate border restrictions, international trade, or allocation of medical supplies. To ensure cooperation among countries in SEA, regional surveillance and information sharing should be indispensable for each country, which can support more precise control and prevention of COVID-19. Although some countries (Singapore, Vietnam) have made efforts to donate medical equipment and provide aid to their neighbors [3], we hope that there will be more interactive assistance among countries in SEA in the future, especially along the continuous emergence of COVID-19 variants. Furthermore, multilateral collaboration with countries or regions (e.g., European Union) outside SEA is considered as an effective and valuable way to help respond to the pandemic, although there will still be a long way to go in the future [3]. Overall, SEA should focus on regional coordination with long-term plans to solve a series of threats caused by the epidemic, considering the long-term challenges caused by the COVID-19 pandemic.

Limitations and Future Work

Despite the insights from our study, there are notable limitations in the COVID-19 data. There are only seven out of twelve countries provided with data at the primary administrative district level, so we were not able to explore the complete propagation process in SEA. In addition, if high spatial resolution data become available (i.e., city, county, even block or subzone), more specific and detailed patterns could be revealed.

As a matter of fact, many previous studies were also faced with a shortage or loss of available data (i.e., insufficient pediatrics data) [54-56]. Moreover, insufficient knowledge of the data or dynamics would lead to invalidity and unreliability of responses to COVID-19 [57]. Therefore, we strongly suggest that public health authorities should disclose more representative and reliable data [58, 59]. Apart from this, COVID-19 transmission and its impacts have shown environmental inequality in terms of household income, education level, age, gender etc. [60, 61]. The potential correlation between environmental inequality and COVID-19 should be further studied using diverse data, which can provide significant insights to resource allocation and regional prevention.

Data Availability

All data produced in the present study are available upon reasonable request to the authors

Reference

  1. 1.↵
    Organization WH. Transmission of SARS-CoV-2: implications for infection prevention precautions: scientific brief, 09 July 2020. World Health Organization, 2020;. Available from: https://apps.who.int/iris/handle/10665/333114. License: CC BY-NC-SA 3.0 IGO
  2. 2.↵
    Dong E, Du H, Gardner L. An interactive web-based dashboard to track COVID-19 in real time. Lancet Infect Dis. 2020;20(5):533–4. doi: 10.1016/S1473-3099(20)30120-1.
    OpenUrlCrossRefPubMed
  3. 3.↵
    Amul GG, Ang M, Kraybill D, Ong SE, Yoong J. Responses to COVID-19 in Southeast Asia: Diverse paths and ongoing challenges. Asian Econ. Policy Rev. 2021;. doi: 10.1111/aepr.12362
    OpenUrlCrossRef
  4. 4.↵
    Ford M, Ward K. COVID-19 in Southeast Asia: Implications for workers and unions. J Ind Relat. 2021:00221856211000097. doi: 10.1177/00221856211000097
    OpenUrlCrossRef
  5. 5.↵
    ASEAN Economies Hit Hard by Escalating Covid Delta Waves. Ihs Markit. 2021 Sept 03. Available from: https://ihsmarkit.com/research-analysis/asean-economies-hit-hard-by-escalating-covid-delta-waves-sep21.html/ [accessed 2021-12-17]
  6. 6.↵
    Nath A, Aditya S. A Comparative Study on Covid-19 Corona Virus Variants: Alpha, Beta, Gamma, Delta. Researchgate.net. doi: 10.22214/ijraset.2021.37743
    OpenUrlCrossRef
  7. 7.↵
    Chookajorn T, Kochakarn T, Wilasang C, Kotanan N, Modchang C. Southeast Asia is an emerging hotspot for COVID-19. Nat Med. 2021;27(9):1495–6. doi: 10.1038/s41591-021-01471-x
    OpenUrlCrossRef
  8. 8.↵
    Zainuddin A. Southeast Asia Nations Plan To ‘Live With Virus’ As Delta Surges. Thediplomat. 2021 Sept 07. Available from: https://thediplomat.com/2021/09/southeast-asia-nationss-plan-to-live-with-virus-as-delta-surges/ [accessed 2021-12-17]
  9. 9.↵
    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. doi: 10.1016/j.scitotenv.2020.140033
    OpenUrlCrossRefPubMed
  10. 10.↵
    Yin L, Zhang H, Li Y, Liu K, Chen T, Luo W, et al. A data driven agent-based model that recommends non-pharmaceutical interventions to suppress Coronavirus disease 2019 resurgence in megacities. J R Soc Interface. 2021;18(181):20210112. doi: 10.1098/rsif.2021.0112
    OpenUrlCrossRef
  11. 11.↵
    Lai S, Ruktanonchai NW, Zhou L, Prosper O, Luo W, Floyd JR, et al. Effect of non-pharmaceutical interventions to contain COVID-19 in China. Nature. 2020;585(7825):410–3. doi: 10.1038/s41586-020-2293-x
    OpenUrlCrossRefPubMed
  12. 12.↵
    Xu M, Cao C, Zhang X, Lin H, Yao Z, Zhong S, et al. Fine-scale space-time cluster detection of COVID-19 in Mainland China using retrospective analysis. Int J Environ Res Public Health. 2021;18(7):3583. doi: 10.3390/ijerph18073583
    OpenUrlCrossRef
  13. 13.↵
    Desjardins MR, Hohl A, Delmelle EM. Rapid surveillance of COVID-19 in the United States using a prospective space-time scan statistic: Detecting and evaluating emerging clusters. Appl Geogr. 2020;118:102202. doi: 10.1016/j.apgeog.2020.102202
    OpenUrlCrossRefPubMed
  14. 14.↵
    Benita F, Gasca-Sanchez F. The main factors influencing COVID-19 spread and deaths in Mexico: A comparison between phases I and II. Appl Geogr. 2021;134:102523. doi: 10.1016/j.apgeog.2021.102523
    OpenUrlCrossRef
  15. 15.↵
    Rosillo N, Del-Águila-Mejía J, Rojas-Benedicto A, Guerrero-Vadillo M, Peñuelas M, Mazagatos C, et al. Real time surveillance of COVID-19 space and time clusters during the summer 2020 in Spain. BMC Public Health. 2021;21(1):1–11. doi: 10.1186/s12889-021-10961-z
    OpenUrlCrossRefPubMed
  16. 16.↵
    Ullah S, Nor NHM, Daud H, Zainuddin N, Gandapur MSJ, Ali I, et al. Spatial cluster analysis of COVID-19 in Malaysia (Mar-Sep, 2020). Geospat Health. 2021;16(1). doi: 10.4081/gh.2021.961
    OpenUrlCrossRef
  17. 17.↵
    Masrur A, Yu M, Luo W, Dewan A. Space-time patterns, change, and propagation of COVID-19 risk relative to the intervention scenarios in Bangladesh. Int J Environ Res Public Health. 2020;17(16):5911. doi: 10.3390/ijerph17165911
    OpenUrlCrossRef
  18. 18.↵
    Zhu M, Kleepbua J, Guan Z, Chew SP, Tan JW, Shen J, et al., editors. Early Spatiotemporal Patterns and Population Characteristics of the COVID-19 Pandemic in Southeast Asia. Healthcare; 2021: Multidisciplinary Digital Publishing Institute. doi: 10.3390/healthcare9091220
    OpenUrlCrossRef
  19. 19.↵
    Salman AM, Ahmed I, Mohd MH, Jamiluddin MS, Dheyab MA. Scenario analysis of COVID-19 transmission dynamics in Malaysia with the possibility of reinfection and limited medical resources scenarios. Comput Biol Med. 2021;133:104372. doi: 10.1016/j.compbiomed.2021.104372
    OpenUrlCrossRef
  20. 20.↵
    Bouchnita A, Jebrane A. A hybrid multi-scale model of COVID-19 transmission dynamics to assess the potential of non-pharmaceutical interventions. Chaos Solit Fractals. 2020;138:109941. doi: 10.1016/j.chaos.2020.109941
    OpenUrlCrossRef
  21. 21.
    Nian F, Shi Y, Cao J. COVID-19 Propagation Model Based on Economic Development and Interventions. Wirel Pers Commun. 2021:1–11. doi: 10.1007/s11277-021-08998-9
    OpenUrlCrossRef
  22. 22.↵
    Ruktanonchai NW, Floyd J, Lai S, Ruktanonchai CW, Sadilek A, Rente-Lourenco P, et al. Assessing the impact of coordinated COVID-19 exit strategies across Europe. Science. 2020;369(6510):1465–70. doi: 10.1126/science.abc5096
    OpenUrlAbstract/FREE Full Text
  23. 23.↵
    Khairulbahri M. Lessons learned from three Southeast Asian countries during the COVID-19 pandemic. J Policy Model. 2021;. doi: 10.1016/j.jpolmod.2021.09.002
    OpenUrlCrossRef
  24. 24.↵
    Arnakim L, Kibtiah T, editors. Response of ASEAN member states to the spread of COVID-19 in Southeast Asia. IOP Conference Series: Environ Earth Sci; 2021: IOP Publishing. doi: 10.1088/1755-1315/729/1/012100
    OpenUrlCrossRef
  25. 25.↵
    Southeast Asia Covid-19 Tracker. Center For Strategic & International Studies. 2021 Dec 15. Available from: https://www.csis.org/programs/southeast-asia-program/projects/southeast-asia-covid-19-tracker/ [accessed 2021-12-17]
  26. 26.↵
    Kulldorff M. SaTScanTM User Guide. 2021;. Available from: https://www.satscan.org/cgi-bin/satscan/register.pl/SaTScan_Users_Guide.pdf?todo=process_userguide_download/ [accessed 2021-12-17]
  27. 27.↵
    Kulldorff M. Prospective time periodic geographical disease surveillance using a scan statistic. J R Stat Soc Ser A Stat Soc. 2001;164(1):61–72. doi: 10.1111/1467-985X.00186
    OpenUrlCrossRef
  28. 28.↵
    Kulldorff M. A spatial scan statistic. Commun Stat Theory Methods. 1997;26(6):1481–96. doi: 10.1080/03610929708831995
    OpenUrlCrossRef
  29. 29.↵
    Hohl A, Delmelle EM, Desjardins MR, Lan Y. Daily surveillance of COVID-19 using the prospective space-time scan statistic in the United States. Spat Spatiotemporal Epidemiol. 2020;34:100354. doi: 10.1016/j.sste.2020.100354
    OpenUrlCrossRef
  30. 30.↵
    Desjardins M, Whiteman A, Casas I, Delmelle E. Space-time clusters and co-occurrence of chikungunya and dengue fever in Colombia from 2015 to 2016. Acta Trop. 2018;185:77–85. doi: 10.1016/j.actatropica.2018.04.023
    OpenUrlCrossRef
  31. 31.↵
    Martines MR, Ferreira RV, Toppa RH, Assunção L, Desjardins MR, Delmelle EM. Detecting space–time clusters of COVID-19 in Brazil: mortality, inequality, socioeconomic vulnerability, and the relative risk of the disease in Brazilian municipalities. J Geogr Syst 2021;23(1):7–36. doi: 10.1007/s10109-020-00344-0
    OpenUrlCrossRef
  32. 32.↵
    Kan Z, Kwan M-P, Wong MS, Huang J, Liu D. Identifying the space-time patterns of COVID-19 risk and their associations with different built environment features in Hong Kong. Sci Total Environ 2021;772:145379. doi: 10.1016/j.scitotenv.2021.145379
    OpenUrlCrossRef
  33. 33.↵
    Yang L, Dai J, Zhao J, Wang Y, Deng P, Wang J. Estimation of incubation period and serial interval of COVID-19: analysis of 178 cases and 131 transmission chains in Hubei province, China. Epidemiol Infect 2020;148. doi: 10.1017/S0950268820001338
    OpenUrlCrossRefPubMed
  34. 34.↵
    Pemerintah Terapkan PPKM Mikro Di Seluruh Provinsi Mulai 1 Juni. Cabinet Secretariat of the Republic of Indonesia. 2021 May 25. Available from: https://diskominfo.anambaskab.go.id/berita/pemerintah-terapkan-ppkm-mikro-di-seluruh-provinsi-mulai-1-juni/ [accessed 2021-12-17]
  35. 35.↵
    What the Delta variant did to South-East Asia. The Economist. 2021 Sept 11. Available from: https://www.economist.com/asia/what-the-delta-variant-did-to-south-east-asia/21804360/ [accessed 2021-12-17]
  36. 36.↵
    Rocklöv J, Sjödin H. High population densities catalyse the spread of COVID-19. J Travel Med 2020;27(3):taaa038. doi: 10.1093/jtm/taaa038
    OpenUrlCrossRefPubMed
  37. 37.↵
    Lian M, Warner RD, Alexander JL, Dixon KR. Using geographic information systems and spatial and space-time scan statistics for a population-based risk analysis of the 2002 equine West Nile epidemic in six contiguous regions of Texas. Int J Health Geogr 2007;6(1):1–10. doi: 10.1186/1476-072X-6-42
    OpenUrlCrossRefPubMed
  38. 38.↵
    Yin H, Sun T, Yao L, Jiao Y, Ma L, Lin L, et al. Association between population density and infection rate suggests the importance of social distancing and travel restriction in reducing the COVID-19 pandemic. Environ Sci Pollut Res 2021:1–7. doi: 10.1007/s11356-021-12364-4
    OpenUrlCrossRef
  39. 39.↵
    Zhang S, Wang Z, Chang R, Wang H, Xu C, Yu X, et al. COVID-19 containment: China provides important lessons for global response. Front Med 2020:1. doi: 10.1007/s11684-020-0766-9
    OpenUrlCrossRef
  40. 40.↵
    Jayant M. Delta variant makes COVID-19 border closures unsustainable in Southeast Asia. CNA. 2021. Commentary. 2021 Oct 12. Available from: https://www.channelnewsasia.com/commentary/border-closures-delta-variant-vaccination-rates-vaccinated-travel-lanes-2234886/ [accessed 2021-12-17]
  41. 41.↵
    Boettke P, Powell B. The political economy of the COVID-19 pandemic. South Econ J 2021;87(4):1090–106. doi: 10.1002/soej.12488
    OpenUrlCrossRef
  42. 42.↵
    Zhu R, Anselin L, Batty M, Kwan M-P, Chen M, Luo W, et al. The effects of different travel modes and travel destinations on COVID-19 transmission in global cities. Sci Bull. 2021;. doi: 10.1016/j.scib.2021.11.023
    OpenUrlCrossRef
  43. 43.↵
    PREPARING FOR OUR TRANSITION TOWARDS COVID RESILIENCE. Ministry of Health Singapore. 2021 Aug 06. Available from: https://www.moh.gov.sg/news-highlights/details/preparing-for-our-transition-towards-covid-resilience/ [accessed 2021-12-17]
  44. 44.↵
    Corey L, Mascola JR, Fauci AS, Collins FS. A strategic approach to COVID-19 vaccine R&D. Science. 2020;368(6494):948–50. doi: 10.1126/science.abc5312
    OpenUrlAbstract/FREE Full Text
  45. 45.↵
    Fowlkes A, Gaglani M, Groover K, Thiese MS, Tyner H, Ellingson K, et al. Effectiveness of COVID-19 vaccines in preventing SARS-CoV-2 infection among frontline workers before and during B. 1.617. 2 (Delta) variant predominance—eight US locations, December 2020–August 2021. Morb Mortal Wkly Rep 2021;70(34):1167. doi: 10.15585/mmwr.mm7034e4
    OpenUrlCrossRefPubMed
  46. 46.
    Leshem E, Wilder-Smith A. COVID-19 vaccine impact in Israel and a way out of the pandemic. The Lancet. 2021;397(10287):1783–5. doi: 10.1016/S0140-6736(21)01018-7
    OpenUrlCrossRef
  47. 47.↵
    Makhoul M, Chemaitelly H, Ayoub HH, Seedat S, Abu-Raddad LJ. Epidemiological differences in the impact of COVID-19 vaccination in the United States and China. Vaccines 2021;9(3):223. doi: 10.3390/vaccines9030223
    OpenUrlCrossRef
  48. 48.↵
    Keegan L, Truelove SA, Lessler J. Progress of the Delta variant and erosion of vaccine effectiveness, a warning from Utah. medRxiv. 2021;. doi: 10.1101/2021.08.09.21261554
    OpenUrlAbstract/FREE Full Text
  49. 49.↵
    Bewick S. Viral Variants and Vaccinations: If We Can Change the COVID-19 Vaccine, Should We? medRxiv. 2021;. doi: 10.1101/2021.01.05.21249255
    OpenUrlAbstract/FREE Full Text
  50. 50.↵
    Albouy-Llaty M, Martin C, Benamouzig D, Bothorel E, Munier G, Simonin C, et al. Positioning Digital Tracing Applications in the Management of the COVID-19 Pandemic in France. J Med Internet Res. 2021;23(10):e27301. doi: 10.2196/27301
    OpenUrlCrossRef
  51. 51.↵
    Wong CK, Wong JY, Tang EH, Au CH, Lau KT, Wai AK. Impact of national containment measures on decelerating the increase in daily new cases of COVID-19 in 54 countries and 4 epicenters of the pandemic: comparative observational study. J Med Internet Res. 2020;22(7):e19904. doi: 10.2196/19904
    OpenUrlCrossRef
  52. 52.↵
    The L. COVID-19: where do we go from here? The Lancet. 2021;398(10318):2207. doi: 10.1016/S0140-6736(21)02790-2
    OpenUrlCrossRef
  53. 53.↵
    Singh S, McNab C, Olson RM, Bristol N, Nolan C, Bergstrøm E, et al. How an outbreak became a pandemic: a chronological analysis of crucial junctures and international obligations in the early months of the COVID-19 pandemic. The Lancet. 2021;398(10316):2109–24. doi: 10.1016/S0140-6736(21)01897-3
    OpenUrlCrossRef
  54. 54.↵
    Leon-Abarca JA. Obesity and immunodeficiencies are the main pre-existing conditions associated with mild to moderate COVID-19 in children. Pediatr Obes 2020;. doi: 10.1111/ijpo.12713
    OpenUrlCrossRef
  55. 55.
    Noh JY, Yoon JG, Seong H, Choi WS, Sohn JW, Cheong HJ, et al. Asymptomatic infection and atypical manifestations of COVID-19: comparison of viral shedding duration. Journal infect 2020;81(5):816. doi: 10.1016/j.jinf.2020.05.035
    OpenUrlCrossRef
  56. 56.↵
    Woolf SH, Chapman DA, Sabo RT, Zimmerman EB. Excess deaths from COVID-19 and other causes in the US, March 1, 2020, to January 2, 2021. JAMA. 2021;325(17):1786–9. doi: 10.1001/jama.2021.5199
    OpenUrlCrossRefPubMed
  57. 57.↵
    Hirabayashi A, Kajihara T, Yahara K, Shibayama K, Sugai M. Impact of the COVID-19 pandemic on the surveillance of antimicrobial resistance. J Hosp Infect 2021;117:147–56. doi: 10.1016/j.jhin.2021.09.011
    OpenUrlCrossRef
  58. 58.↵
    Lange KW. The prevention of COVID-19 and the need for reliable data. Mov Nutr Health Dis. 2020;4. doi: 10.5283/mnhd.25
    OpenUrlCrossRef
  59. 59.↵
    Dron L, Dillman A, Zoratti MJ, Haggstrom J, Mills EJ, Park JJ. Clinical Trial Data Sharing for COVID-19–Related Research. J Med Internet Res. 2021;23(3):e26718. doi:10.2196/26718
    OpenUrlCrossRef
  60. 60.↵
    Mollalo A, Vahedi B, Rivera KM. GIS-based spatial modeling of COVID-19 incidence rate in the continental United States. Sci Total Environ 2020;728:138884. doi: 10.1016/j.scitotenv.2020.138884
    OpenUrlCrossRef
  61. 61.↵
    Spotswood EN, Benjamin M, Stoneburner L, Wheeler MM, Beller EE, Balk D, et al. Nature inequity and higher COVID-19 case rates in less-green neighbourhoods in the United States. Nat Sustain. 2021:1–7. doi: 10.1038/s41893-021-00781-9
    OpenUrlCrossRef
Back to top
Posted December 21, 2021.
Download PDF
Data/Code
Email
Spatiotemporal patterns and progression of the Delta variant of COVID-19 and their health intervention linkages in Southeast Asia
Luo Wei, Liu Zhaoyin, Zhou Yuxuan, Zhao Yumin, Li Yunyue Elita, Arif Masrur, Manzhu Yu
medRxiv 2021.12.20.21268140; doi: https://doi.org/10.1101/2021.12.20.21268140
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area