Skip to main content
Download PDF
Download ePub

Spatiotemporal and meteorological relationships in dengue transmission in the Dominican Republic, 2015–2019

Tropical Medicine and Health volume 51, Article number: 32 (2023) Cite this article

Abstract

Dengue has broadened its global distribution substantially in the past two decades, and many endemic areas are experiencing increases in incidence. The Dominican Republic recently experienced its two largest outbreaks to date with 16,836 reported cases in 2015 and 20,123 reported cases in 2019. With continued increases in dengue transmission, developing tools to better prepare healthcare systems and mosquito control agencies is of critical importance. Before such tools can be developed, however, we must first better understand potential drivers of dengue transmission. To that end, we focus in this paper on determining relationships between climate variables and dengue transmission with an emphasis on eight provinces and the capital city of the Dominican Republic in the period 2015–2019. We present summary statistics for dengue cases, temperature, precipitation, and relative humidity in this period, and we conduct an analysis of correlated lags between climate variables and dengue cases as well as correlated lags among dengue cases in each of the nine locations. We find that the southwestern province of Barahona had the largest dengue incidence in both 2015 and 2019. Among all climate variables considered, lags between relative humidity variables and dengue cases were the most frequently correlated. We found that most locations had significant correlations with cases in other locations at lags of zero weeks. These results can be used to improve predictive models of dengue transmission in the country.

Introduction

Global incidence of dengue fever has increased substantially in recent decades, with the range of dengue expanding from only nine countries before 1970 to at least 129 countries today [1,2,3,4]. In addition to its rapid global spread, dengue outbreaks in endemic regions are resulting in increasingly larger numbers of cases and contributing to a growing burden on public health systems. Today, it is estimated that over 390 million people are at risk of contracting dengue [5]. Dengue is primarily distributed across regions of the world with tropical and subtropical climates, although in the past two decades, dengue cases have occurred with greater frequency in temperate zones as well [6,7,8,9]. In 2021, 1,254,648 cases and 436 deaths were reported in the Americas [10]. According to the Pan American Health Organization (PAHO) the countries in the Caribbean reporting the most cases of dengue between 2014 and 2021 are the Dominican Republic, Martinique, Guadeloupe, French Guiana, and Cuba, with the Dominican Republic having 60% more cases in that time as the country reporting the second highest number [10]. The Dominican Republic also reports the highest number of severe dengue cases and deaths in the Caribbean [10]. In 2019, the Dominican Republic experienced its largest outbreak to date with a 1145% increase in cases from 2018 [11]. The cumulative incidence in 2019 was 194.85 cases per 100,000 people, which is a 142% increase from the average incidence between 2005 and 2014 [12, 13].

With outbreaks becoming increasingly severe in the Dominican Republic and other regions, it is more imperative than ever to understand drivers of epidemic dengue. Potential drivers of global spread of dengue include increases in urbanization, more frequent global travel, and changes in temperature and precipitation [14,15,16]. At local scales, transmission of dengue can also be a function of socioeconomic and demographic characteristics, connectivity to other regions, human behavior, volume of tourism, and rates of migration [14,15,16,17,18]. Many of these variables play an important role in developing and sustaining an environment that is suitable for the vectors of dengue, Aedes aegypti and Aedes albopictus, which in turn amplifies risk of transmission [19,20,21].

Because there is an inherent delay between human cases of dengue resulting from the intermediate mosquito host and a serial interval of 15–17 days, early detection of new dengue outbreaks can be complicated [22, 23]. In the last decade, efforts have been made to improve early detection of dengue outbreaks by improving surveillance and warning [22, 24,25,26]. These early warning systems are mathematical and statistical models that integrate data to provide predictions for changes in dengue transmission that may indicate outbreaks. Chief among these data are climate variables such as temperature, precipitation, or humidity which are all positively correlated with Aedes mosquito populations and dengue transmission [22]. However, before such early warning systems can be developed, relationships between climate variables and dengue cases must be explored to determine which climate variables are most important to local and regional dengue transmission. Recently, Freitas et al. thoroughly analyzed the 2019 dengue outbreak in the Dominican Republic and found that a model incorporating temperature and rainfall with delays of 2–5 weeks provided good predictions of dengue transmission [27]. This work marks an important first step in developing better predictive models for the country.

Herein, we aim to build upon this work by analyzing dengue activity in the Dominican Republic between 2015 and 2019 and exploring relationships between climate variables and dengue cases. We present descriptive analysis of each data set used in the study along with analysis of correlations in lags between variables. We further investigate correlations in lags between provinces in the Dominican Republic to understand potential movement of dengue throughout the country. The work presented herein provides a foundation on which statistical and mathematical models can be constructed to further study drivers of previous outbreaks and to predict future outbreaks.

Materials and methods

Study site

This study was conducted in the Dominican Republic, a Caribbean country that occupies the eastern two-thirds of the Island of Hispaniola. The 2019 estimate of the population size of the Dominican Republic is 10,448,499 people [27]. The country is divided geopolitically into 31 provinces and Distrito Nacional (the capital city) [28].

The Dominican Republic has perhaps the most diverse climate of all of the Caribbean islands because of the presence of high mountains and abundant coastal regions [29]. Much of the country, however, has a tropical climate with mean annual temperatures ranging between 22 and 31 °C [29,30,31,32]. Rainy seasons vary geographically with the northern part of the country experiencing heavier rain from November to January and much of the rest of the country having its rainy season May–November. Mean annual precipitation in the country ranges between 400 mm in the southwest to more than 2200 mm in the mountain regions [29]. The majority of southern coastal regions experience a mean rainfall of around 1000 mm while northern coasts typically have a higher mean annual rainfall of 1600 mm or more [29].

In this work, we focus on nine provinces throughout the country. These nine were chosen because they are the provinces for which we were able to obtain both meteorological and epidemiological data between 2015 and 2019. The nine provinces included in this study are Barahona, La Altagracia, La Romana, Monte Cristi, Puerto Plata, Samaná, Santiago, Santo Domingo, and Distrito Nacional (Fig. 1). The provinces cover each of the three major regions of the country: North (Monte Cristi, Puerto Plata, Santiago, Samaná), South (Barahona), and East/Southeast (Santo Domingo, Distrito Nacional, La Romana, La Altagracia). The nine provinces include 6,295,775 people (2019 estimate), representing 60.78% of the total population of the country.

Fig. 1
figure 1

Provinces of the Dominican Republic. Provinces at the focus of this study are highlighted in green

Full size image

Data collection

Dengue cases

The number of weekly reported cases for the period January 2015 to December 2019 was provided by Sistema Nacional de Vigilancia Epidemiológica de la Dirección General de Epidemiológica (Ministerio de Salud Pública). The epidemiological week was defined as Sunday to Saturday. Cases include suspected and laboratory-confirmed cases aggregated at the province level according to surveillance definitions [33]. Dengue Incidence Rate (DIR) was calculated using the number of new cases, divided by the local population each year, multiplied by 100,000 inhabitants. Figure 2a shows dengue incidence for the five years across the nine provinces included in the study.

Fig. 2
figure 2

Time series of epidemiological and meteorological metrics across the five years and nine provinces at the focus of this study. Provinces are arranged by latitude (northernmost to southernmost). a Dengue incidence per 100,000 people; b total weekly precipitation (mm3; log scale); c average weekly relative humidity (%); (d–f) average, minimum, and maximum weekly temperature, respectively (°C)

Full size image

Meteorological data

Meteorological data were obtained by supplementing official national data (from ONAMET) with data provided by the U.S. National Aeronautics and Space Administration (NASA). Previous studies have confirmed this approach for collecting data, especially where there are gaps in reliable data [24, 34]. We calculated summary values (minimum, maximum, mean, sum) of meteorological parameters by epidemiologic weeks. Figure 2b–d show average weekly temperature, total weekly precipitation, and average weekly relative humidity for the 5 years across the nine provinces included in the study.

Population data

We obtained population data from Oficina Nacional de Estadísticas (ONE) [27]. This data includes the total population and the population density for each province (Table 1). Distrito Nacional has the highest population density, and Santo Domingo province has the highest population.

Table 1 Demographic characteristics from provinces of the Dominican Republic included in this study, 2019
Full size table

Statistical analysis

Summary statistics

We conducted a preliminary statistical analysis of epidemiological data to highlight changes in dengue cases across provinces and across years. We calculated summary statistics of dengue cases for the entire country each year. We focus much of our analysis of cases on the years 2015 and 2019, when large epidemics took place. Our initial findings support subsequent analysis at the province level to assess the association between explanatory variables and the distribution of the disease over space and time. We calculate dengue incidence in each province in 2015 and 2019 as well as descriptive statistics such as mean, maximum, minimum, and standard deviation of meteorological variables for each province. The results were obtained by calculations in Microsoft Excel, through the Excel Data Analysis Tool. Maps with spatial data were generated by QGIS 3.16.6. These maps improve our understanding of dengue transmission in different regions and how transmission evolves spatially over time [35].

Correlated lags analysis

We compared time series of dengue case data across provinces with climate data by conducting an analysis of correlations of lags between data sets. We assume a unidirectional relationship between dengue cases and meteorological variables. We first implemented a standard prewhitening approach to remove any effects of autocorrelation within data sets [36]. We fit each time series data set to seasonal autoregressive integrated moving average (SARIMA) models. We calculated cross-correlation functions between residuals of time series for different weekly summary data for these variables with lags up to 10 weeks. We chose this cutoff for lags because of the relatively short time scale of our data (5 years) and because it is a biologically reasonable time for weather events to impact mosquito development and disease transmission given the development time and generation time of mosquitoes as well as the serial interval of dengue [23, 37, 38]. We tested for significance of correlated lags with a two-tailed t-test to test the null hypothesis that the correlation was equivalent to 0. We report the lags with the highest correlation along with p-values at the 0.10, 0.05, and 0.01 confidence levels.

We also conducted a correlated lag analysis among provinces. We calculated correlations between lags in residuals of time series of cases in each province. We determined the significance of these correlations with a two-tailed t-test to test the null hypothesis that the correlation was equivalent to 0. We report the lags with the highest correlation along with p-values at the 0.10, 0.05, and 0.01 confidence levels. All correlation analyses were conducted in R 4.1.0 [39]. Fitting of SARIMA models was conducted using the auto.arima() function in the forecast package [40, 41]. Cross-correlation functions were calculated with the ccf() function in the R base installation [42].

Results

Dengue cases: spatiotemporal analysis

We first calculated descriptive statistics for dengue cases in the country each year (total cases, mean cases per week, standard deviation of cases per week, minimum number of cases per week, maximum number of cases per week, and the dengue incidence rate per year). Table 2 presents these results. In 2015 and 2019 there were major outbreaks with dengue incidence rates of 168.69 and 194.27 cases per 100,000 residents, respectively.

Table 2 Descriptive statistics for dengue cases each year, 2015–2020
Full size table

The two largest outbreaks each started during late spring (April–May) and continued for approximately one year (Fig. 3). The 2015 outbreak peaked in October, while the 2019 outbreak peaked in August. Central provinces experienced the highest incidence during the 2015 outbreaks, while in 2019, provinces in the north and south experienced the highest incidence (Fig. 4). The increase in incidence in the northern and southern provinces could be related to socioeconomic factors or differences in climate variables between the 2 years. Barahona in the southwest, along with Hermanas Mirabel, Sánchez Ramirez, and San José de Ocoa in the center and Hato Mayor in the east all experienced similarly high incidence in both 2015 and 2019. Puerto Plata in the north along with Distrito Nacional in the southeast and many other coastal provinces experienced similarly lower incidence in both outbreaks. Puerto Plata and Distrito Nacional are both popular destinations for international travel and thus are likely to employ more aggressive mosquito control and dengue prevention practices. Table 3 shows incidence calculations for the nine provinces of focus for this study along with the national incidence for both 2015 and 2019.

Fig. 3
figure 3

National dengue incidence by week from 2015 to 2019 in the Dominican Republic

Full size image
Fig. 4
figure 4

Spatial description of dengue incidence in 2015 (a) and 2019 (b)

Full size image
Table 3 Dengue incidence per 100,000 people by province for the 2015 and 2019 outbreaks
Full size table

Climate variables and dengue cases

The majority of the Dominican Republic has a tropical climate with hot temperatures all year and the warmest months being May to October. There is a rainy season between late April and October, while the northern coast, exposed to the trade winds, is rainy throughout the year. On the southern coast, there is a considerable amount of precipitation because it is not protected by mountains. As a Caribbean country, the rains occur mainly as short showers and thunderstorms which are sometimes intense and often concentrated in short periods of time. Table 4 summarizes the climate statistics for the nine provinces studied in the two outbreaks. In both years, the values for all climate variables did not vary so much.

Table 4 Summary statistics of climate variables
Full size table

Although average values and ranges of climate variables are useful for pointing out variations across years, it is important to consider the temporal variation in climate variables and how they might relate to dengue outbreaks. As an example, we show in Fig. 5 temporal variation in climate variables and dengue cases in 2015 and 2019 in Distrito Nacional. Behavior across provinces were similar and are excluded here for brevity. Temperature and relative humidity are relatively stable throughout the year, although temperatures increase from the beginning of each year until epidemiological weeks 30–35. The cumulative precipitation per week is rather variable and could potentially have more influence on dengue cases. We investigate this relationship, along with relationships between weekly variation in other climate variables and dengue transmission, in the next section.

Fig. 5
figure 5

Climate variables and dengue cases in 2015 (a–c) and 2019 (d–f). Variables included are (a, d) temperature (mean temperature is given by the solid curve; °C); (b, e) precipitation (mm3); and (c, f) relative humidity (%)

Full size image

Cross-correlation analysis

Correlations in lags between dengue cases and climate variables

Table 5 contains correlations between lags in climate variables and dengue cases in the 9 provinces. Although we tested 14 variables, we present only 9 here. Additional file 1: Table S1 in supporting information shows values for the remaining variables we tested. Maximum and average weekly relative humidity were the variables most often significantly correlated with dengue cases at the α = 0.05 or stronger confidence level (6 and 5 of 9 provinces, respectively), followed by weekly minimum temperature (4 of 9 provinces), relative humidity range, mean daily temperature, and maximum weekly temperature (3 of 9 provinces).

Table 5 Correlations in lags between dengue cases and climate variables
Full size table

In Barahona, La Altagracia, La Romana, Monte Cristi, Puerto Plata, and Santo Domingo, lags between dengue cases and average relative humidity were significantly correlated at the α = 0.05 or stronger confidence level with lags of 4–10 weeks. For most provinces, these correlations were negative, suggesting that, for example, decreasing average relative humidity may be associated with increases in dengue cases. Lags of 3–10 weeks between dengue cases and maximum relative humidity were significantly correlated at the α = 0.05 or stronger confidence level in Barahona, La Romana, Monte Cristi, Santo Domingo, and Distrito Nacional. For all provinces except Distrito Nacional, these correlations were also negative. For lags with the weekly range of relative humidity, lags of 5–10 weeks were significantly positively correlated at the α = 0.05 or stronger confidence level for Barahona, La Romana, and Puerto Plata provinces, suggesting higher ranges in weekly relatively humidity could be associated with dengue cases. For provinces for which we found significant positive correlations between dengue cases and relative humidity variables, rates of increases in cases lagged behind those of other provinces, suggesting that perhaps timing of cases in provinces relative to one another played a role in these relationships. We explore this further later in this section.

In general, correlations between lags in dengue cases and temperature variables were significant less often than those we found for relative humidity; however, lags between cases and minimum weekly temperature were significantly correlated at the α = 0.05 or stronger level for lags of 1–7 weeks for the provinces of Barahona, Monte Cristi, Samaná, and Santo Domingo. In all four of these provinces, correlations were positive, suggesting that higher minimum weekly temperatures were associated with more dengue cases. Lags of 2–5 weeks between maximum weekly temperature and dengue cases were also strongly correlated at the α = 0.05 or stronger level for La Romana, Monte Cristi, and Puerto Plata provinces, although correlations in La Romana and Puerto Plata were negative whereas the correlation in Monte Cristi was positive. For mean daily temperature, lags of 1–5 weeks between this variable and dengue cases were significantly correlated at the α = 0.05 or stronger confidence level. Notably, lags of 1 week between mean daily temperature and dengue cases were significant in Santo Domingo and Distrito Nacional; however, the direction of these correlations differed. Again, it is possible that these differences in the directions of relationships are a result of lags between the timing of outbreaks in different provinces.

No correlations of lags with total precipitation and average precipitation were significant at the α = 0.05 confidence level. This result is surprising because we would expect dengue transmission in tropical climates to be positively correlated with precipitation given the important role of water in the mosquito’s life cycle [43].

The provinces with the most significantly correlated lags in climate variables at α = 0.05 or stronger confidence level were Santo Domingo (7 variables); Puerto Plata (6); Barahona, La Romana, and Distrito Nacional (5); and Monte Cristi (4). For all other provinces included in the study, only 2 or fewer climate variables were significantly correlated.

Correlations in lags between cases in different provinces

Almost all of the largest correlations in lags between cases in different provinces were significant (p < 0.05), and most were highly significant (p < 0.01). Only three of 72 tested correlations were not significant at the α = 0.05 or stronger confidence level: cases in Barahona as a predictor of cases in Distrito Nacional, La Altagracia as a predictor for Santiago, and Santo Domingo as a predictor for Santiago. In most cases, correlations between provinces with a lag of τ = 0 were the strongest (20 of 72 or 27.78%), suggesting that cases were largely synced across provinces (Table 6). Correlations with lags of 2 weeks (16.67%), 1 week (12.5%), and 4 weeks (11.11%) were also among some of the strongest, indicating that cases across the country closely followed cases in other parts of the country.

Table 6 Correlations in lags (weeks) between dengue cases in each province
Full size table

Cases in Puerto Plata and Santo Domingo were generally in sync with cases in other regions, with cases in both provinces with a lag of τ = 0 having strong positive correlations with cases in four of the eight other provinces. In fact, cases in Puerto Plata had strong positive correlations with cases in all other provinces with a lag of 0–3 weeks. The only exception is Samaná (τ = − 8). For Santo Domingo, cases were positively correlated with cases in other provinces with a lag of 0 or 4 weeks, but cases were negatively correlated with cases in Barahona (τ = − 9) and Distrito Nacional (τ = − 1). It is possible that the strong positive correlations with short-term lags between cases in Puerto Plata and Santo Domingo are due either population size, high rates of tourism, or some combination thereof.

Cases in Barahona, La Romana, Samaná, and Santiago also were positively correlated with all other provinces with lags of 0–10, 0–7, 0–10, and 0–8 weeks, respectively. In the remaining provinces, correlations were mostly positive but cases in each of these five provinces were negatively correlated with cases in one other province. These correlations were with cases at lags of 5–9 weeks, indicating that perhaps decreases in cases in some provinces preceded increases elsewhere (or vice-versa) by several weeks. Moreso than in other provinces, cases in Samaná, Monte Cristi, and La Romana were correlated with cases in other provinces with lags greater than one week, indicating that cases in these provinces may often occur 2 or more weeks before cases in other provinces. Notably, cases in Distrito Nacional were correlated with cases in all other provinces with lags of 4–9 weeks, suggesting that increases in transmission in the capital district may precede outbreaks elsewhere in the country.

Discussion

Herein we characterized dengue incidence at the province level in the Dominican Republic between 2015 and 2019, a period in which the country and the Caribbean region experienced two important epidemics. We focused our study on nine provinces that included all major geographic regions of the country that represented different climate patterns. In our study, we observed different potential drivers of dengue activity in different regions of the country. We anticipate that this study will be a foundation upon which early warning systems and models aimed at predicting dengue activity may be built.

We noted that both major outbreaks (2015 and 2019) occurred after the 30th epidemiological week, which corresponds approximately to late July. This, together with the fluctuations noted even in the years in which no epidemic occurred, indicate a seasonal pattern of dengue transmission. When comparing the epidemiology of dengue in the Dominican Republic with the epidemiology throughout the region of the Americas, a similar behavior was observed for 2015 and 2019, the latter being the year with the highest number of cases recorded in the history of dengue in the Americas [28, 44]. The reduction in the number of cases between 2016 and 2018 could be explained in part by the vector control actions implemented by the Ministry of Health, the adaptation of the pathogen, reduction of susceptible population, or partial immunity to dengue conferred by the wave of Zika virus that moved through the region between 2015 and 2016 [45, 46]. Another possible factor that could have an impact on dengue transmission is the El Niño Southern Oscillation (ENSO). During 2015 and at the beginning of 2019 there were warm periods related to ocean–atmosphere temperature [47]. These results are in line with previous studies that found evidence that ENSO is associated with dengue outbreaks [48, 49]. Considering the number of dengue cases that occurred in 2015 and 2019, the ENSO phenomenon is an important factor to be considered in any early warning system. Even with the short period of time in our analysis (5 years), the major outbreaks corresponded to the only years that were classified as warm periods in relationships to ocean–atmosphere features.

We found that the southwestern province of Barahona had the largest dengue incidence in both 2015 (273.9 per 100,000 people) and 2019 (456.2). Furthermore, dengue activity in Barahona preceded dengue activity in six other provinces by 1–10 weeks, and among all the provinces studied here, cases in Barahona were most often negatively correlated at longer time lags (6–9 weeks) with cases elsewhere, suggesting that outbreaks beginning in Barahona may be starting to die out as outbreaks are still growing in other provinces. It is possible that new cases are introduced to the Dominican Republic in this region through immigration from Haiti or via tourism. It is also possible that individuals who have dengue must travel to other provinces for medical care given that this province has only one public hospital [50]. This could lead to movement of cases into other provinces and throughout the country. Cases in Distrito Nacional preceded cases in all other provinces. Distrito Nacional is the most densely populated among the provinces we considered, and as the capital district, it is important to the national economy and tourism. It is highly connected to other provinces in the country through highways which facilitate movement of people and potentially mosquitoes via movement of tires and other containers that serve as breeding sites [51,52,53]. This result is consistent with other studies showing the importance of population-dense urban areas in regional and national transmission of dengue [54, 55]. These initial findings from the present work suggest that provinces such as Barahona and Distrito Nacional will be important in the development of future predictive models and early warning systems.

In fact, our analyses of lags in cases between provinces could help determine how cases spread spatially in the country by identifying “source” provinces where dengue cases begin (such as Barahona and Distrito Nacional) and “sink” provinces where dengue cases later appear. For example, cases in La Romana and Santo Domingo often trailed cases in other parts of the country. Both provinces are home to several popular tourist attractions and may benefit from increased surveillance and vector control [56]. It is possible that as cases are reported elsewhere in the country, control efforts delay significant amounts of transmission in these provinces. Despite these lags found, most of the lags between provinces were zero, indicating that outbreaks occurred throughout most of the region at about the same time. This could be explained by human movement inside the country. A similar study in the country with Zika virus also suggested that the human mobility and the infrastructure level of each region could influence the transmission of diseases that had Aedes aegypti as a vector [57].

We analyzed climate variables that could contribute to the dengue transmission cycle by their impacts on the Ae. aegypti life cycle. When all five years of dengue case and climate data are considered, lags between relative humidity variables and dengue cases were most highly correlated, indicating relative humidity as a good predictor of dengue transmission throughout the region. This is consistent with previous studies showing strong correlations between dengue cases and average relative humidity, whether those correlations are positive [58, 59] or negative [60]. Our results showing different directional relationships between dengue and relative humidity (i.e., positive correlations for some provinces and negative correlations for others) emphasize the importance of considering effects of local climate variables on dengue cases and highlight a need for more thorough data collection and analyses of these relationships.

Lags with temperature variables, too, were significantly correlated with cases in many provinces. This result is supported by work showing that temperature influences dengue transmission through its impacts on both the vector life cycle and the virus [15, 19, 61,62,63]. Surprisingly, lags between precipitation and dengue cases were not found to be significantly correlated with cases in any of the provinces. In studies of other tropical regions, precipitation and humidity are often found to be positively correlated with arbovirus activity [29, 58, 59]. It is possible that this is because the Dominican Republic’s unique topography interferes with weather patterns and results in having rainy seasons at different times of the year in different regions [29]. The landscape of the Dominican Republic is composed of chains and valleys with active faults which decreases precipitation from northeast to southwest in winter and spring and increases aridity in western areas where fewer mountain chains are observed. The mountain chains serve as a barrier to trade winds, affecting the humidity of such areas. Additionally, during short periods of time in April and in November, there is a subtropical atmosphere which is influenced by anomalies in sea surface temperatures [29]. While cases could be impacted locally by changes in precipitation, this may not correspond to times at which dengue transmission is occurring elsewhere, which may lead to impacts on correlation. These results are in line with the ones achieved in [64], showing that temperature and humidity have impacts on the transmission chain.

Conclusions

The short period of data included in this study is insufficient for making strong characterizations of relationships. However, in this work we developed a better understanding of which variables have been most strongly associated with dengue cases in this time frame, which includes two important outbreaks. Both outbreaks occurred in the middle of the year, which indicate a need for more heightened vigilance during this period of the year. Additionally, humidity and temperature are the climate variables with the highest correlations with the number of cases. These findings will help inform future work for building predictive models that incorporate climate and spatiotemporal data to characterize province risk and refine public health responses. This initial analysis may also provide the foundation for models based on ARIMA, SARIMA, SARIMAX, or other frameworks for predictive models [65,66,67]. A reliable warning system built on such models and adapted to the intrinsic characteristics of each province (namely climate, demographics, and landforms) could help in the monitoring both the vector and arbovirus transmission [68]. This in turn could lead to a faster intervention of health and entomological authorities, thus decreasing dengue incidence.

This study contributes an important analysis of recent dengue transmission on which more complex spatiotemporal analyses can be conducted. For instance, Distrito Nacional has the most correlations in dengue cases with lags greater than 0 weeks, and cases in Barahona with longer time lags are highly correlated with other regions. It is possible these areas could provide early warning of nationwide outbreaks. The general characterizations of climate and dengue activity along with the correlated lags analysis across the nine provinces included here provide a foundation upon which future studies may build to investigate more intricate relationships between dengue and climate, human movement, and human activity.

Availability of data and materials

The datasets analyzed during the current study are available https://github.com/marobert-biomath/dominican_republic_dengue.

References

  1. WHO (World Health Organization). Dengue and severe dengue. https://www.who.int/news-room/fact-sheets/detail/dengue-and-severe-dengue. Accessed 10 Mar 2022.

  2. Brady OJ, Gething PW, Bhatt S, Messina JP, Brownstein JS, Hoen AG, et al. Refining the global spatial limits of dengue virus transmission by evidence-based consensus. PLoS Negl Trop Dis. 2012;6: e1760.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Simo FBN, Bigna JJ, Kenmoe S, Ndangang MS, Temfack E, Moundipa PF, et al. Dengue virus infection in people residing in Africa: a systematic review and meta-analysis of prevalence studies. Sci Rep. 2019;9:13626.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Lun X, Wang Y, Zhao C, Wu H, Zhu C, Ma D, et al. Characteristics and temporal-spatial analysis of overseas imported dengue fever cases in outbreak Provinces of China, 2005–201. 2021.

  5. Bhatt S, Gething PW, Brady OJ, Messina JP, Farlow AW, Moyes CL, et al. The global distribution and burden of dengue. Nature. 2013;496:504–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Wilder-Smith A. The expanding geographic range of dengue in Australia. Med J Aust. 2021;215:171.

    Article  PubMed  Google Scholar 

  7. Robert MA, Tinunin DT, Benitez EM, Ludueña-Almeida FF, Romero M, Stewart-Ibarra AM, et al. Arbovirus emergence in the temperate city of Córdoba, Argentina, 2009–2018. Scientific Data. 2019;6:276.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Stephenson C, Coker E, Wisely S, Liang S, Dinglasan RR, Lednicky JA. Imported dengue case numbers and local climatic patterns are associated with dengue virus transmission in Florida, USA. Insects. 2022;13:163.

    Article  PubMed  PubMed Central  Google Scholar 

  9. López MS, Jordan DI, Blatter E, Walker E, Gómez AA, Müller GV, et al. Dengue emergence in the temperate Argentinian province of Santa Fe, 2009–2020. Sci Data. 2021;8:134.

    Article  PubMed  PubMed Central  Google Scholar 

  10. PAHO PAHO, WHO WHO. PAHO/WHO Data—Dengue cases|PAHO/WHO. Pan American Health Organization/World Health Organization. https://www3.paho.org/data/index.php/en/mnu-topics/indicadores-dengue-en/dengue-nacional-en/252-dengue-pais-ano-en.html. Accessed 27 Mar 2022.

  11. ACAPS. Dominican Republic Dengue Fever. ACAPS.org. https://www.acaps.org/sites/acaps/files/products/files/20190916_acaps_start_dengue_fever_dominican_republic__0.pdf. Accessed 27 Mar 2022.

  12. Ministerio de Salud Pública. Boletín Epidemiológico Semanal 52-2019 [Internet. Dirección General de Epidemiología. https://digepi.gob.do/docs/Boletines%20epidemiologicos/Boletines%20semanales/2019/Boletin%20Semanal%2052-2019.pdf. Accessed 12 Mar 2022.

  13. Ávila-Agüero ML, Camacho-Badilla K, Brea-Del-Castillo J, Cerezo L, Dueñas L, Luque M, et al. Epidemiología del dengue en Centroamérica y República Dominicana. Rev Chilena Infectol. 2019;36:698–706.

    Article  PubMed  Google Scholar 

  14. Young PR. Arboviruses: a family on the move. In: Hilgenfeld R, Vasudevan SG, editors. Dengue and Zika: control and antiviral treatment strategies. Singapore: Springer; 2018. p. 1–10.

    Google Scholar 

  15. Mordecai EA, Cohen JM, Evans MV, Gudapati P, Johnson LR, Lippi CA, et al. Detecting the impact of temperature on transmission of Zika, dengue, and chikungunya using mechanistic models. PLoS Neglect Trop D. 2017;11: e0005568.

    Article  Google Scholar 

  16. Bowman LR, Tejeda GS, Coelho GE, Sulaiman LH, Gill BS, McCall PJ, et al. Alarm variables for dengue outbreaks: a multi-centre study in Asia and Latin America. PLoS ONE. 2016;11: e0157971.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Young PR, Ng LFP, Hall RA, Smith DW, Johansen CA. 14. Arbovirus Infections. In: Farrar J, Hotez PJ, Junghanss T, Kang G, Lalloo D, White NJ, editors. Manson’s tropical infectious diseases (Twenty-third Edition). London: W.B. Saunders; 2014. p. 129- 161.e3.

    Chapter  Google Scholar 

  18. Velazquez-Castro J, Anzo-Hernandez A, Bonilla-Capilla B, Soto-Bajo M, Fraguela-Collar A. Vector-borne diesease risk indexes in spatially structured populations. PLOS Negl Trop Dis. 2018;12:e0006234.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Mordecai EA, Caldwell JM, Grossman MK, Lippi CA, Johnson LR, Neira M, et al. Thermal biology of mosquito-borne disease. Ecol Lett. 2019;22:1690–708.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Brady OJ, Hay SI. The global expansion of dengue: how aedes aegypti mosquitoes enabled the first pandemic arbovirus. Annu Rev Entomol. 2020;65:191–208.

    Article  CAS  PubMed  Google Scholar 

  21. Estallo EL, Sippy R, Stewart-Ibarra AM, Grech MG, Benitez EM, Ludueña-Almeida FF, et al. A decade of arbovirus emergence in the temperate southern cone of South America: dengue, Aedes aegypti and climate dynamics in Córdoba, Argentina. Ecology. 2020;13:255.

    Google Scholar 

  22. Kuhn K, Campbell-Lendrum D, Haines A, Corvalan C, Anker M. Using climate to predict infectious disease epidemics. 2005.

  23. Aldstadt J, Yoon I-K, Tannitisupawong D, Jarman RG, Thomas SJ, Gibbons RV, et al. Space-time analysis of hospitalised dengue patients in rural Thailand reveals important temporal intervals in the pattern of dengue virus transmission. Tropical Med Int Health. 2012;17:1076–85.

    Article  Google Scholar 

  24. Lowe R, Bailey TC, Stephenson DB, Graham RJ, Coelho CAS, Sá Carvalho M, et al. Spatio-temporal modelling of climate-sensitive disease risk: towards an early warning system for dengue in Brazil. Comput Geosci. 2011;37:371–81.

    Article  Google Scholar 

  25. Lowe R, Bailey TC, Stephenson DB, Jupp TE, Graham RJ, Barcellos C, et al. The development of an early warning system for climate-sensitive disease risk with a focus on dengue epidemics in Southeast Brazil. Stat Med. 2013;32:864–83.

    Article  PubMed  Google Scholar 

  26. Racloz V, Ramsey R, Tong S, Hu W. Surveillance of dengue fever virus: a review of epidemiological models and early warning systems. PLoS Negl Trop Dis. 2012;6.

  27. Estadística (ONE) ON de. Datos y Estadísticas. Oficina Nacional de Estadística (ONE). https://www.one.gob.do/datos-y-estadisticas/. Accessed 30 Jan 2022.

  28. Pan American Health Organization (PAHO). Health in the Americas+, 2017 Edition. Summary: Regional Outlook and Country Profiles. https://iris.paho.org/handle/10665.2/34321. Accessed 12 Mar 2022.

  29. Izzo M, Rosskopf CM, Aucelli PPC, Maratea A, Méndez R, Pérez C, et al. A new climatic map of the Dominican Republic based on the thornthwaite classification. Phys Geogr. 2010;31:455–72.

    Article  Google Scholar 

  30. Beck HE, Zimmermann NE, McVicar TR, Vergopolan N, Berg A, Wood EF. Present and future Köppen-Geiger climate classification maps at 1-km resolution. Scientific Data. 2018;5: 180214.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Peel MC, Finlayson BL, McMahon TA. Updated world map of the Köppen-Geiger climate classification. Hydrol Earth Syst Sci. 2007;11:1633–44.

    Article  Google Scholar 

  32. Cano-Ortiz A, Musarella CM, Pinar Fuentes JC, Pinto Gomes CJ, Cano E. Forests and landscapes of Dominican Republic. Br J Appl Sci Technol. 2015;9:231.

    Article  Google Scholar 

  33. Ministerio de Salud Pública. Protocolo de atencion para el manejo del dengue. 2017.

  34. Chang AY, Parrales ME, Jimenez J, Sobieszczyk ME, Hammer SM, Copenhaver DJ, et al. Combining Google Earth and GIS mapping technologies in a dengue surveillance system for developing countries. Int J Health Geogr. 2009;8:49.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Pigott DM, Howes RE, Wiebe A, Battle KE, Golding N, Gething PW, et al. Prioritising infectious disease mapping. PLoS Negl Trop Dis. 2015;9: e0003756.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Fuenzalida H, Rosenblüth B. Prewhitening of climatological time series. J Clim. 1990;3:382–93.

    Article  Google Scholar 

  37. Siraj AS, Oidtman RJ, Huber JH, Kraemer MUG, Brady OJ, Johansson MA, et al. Temperature modulates dengue virus epidemic growth rates through its effects on reproduction numbers and generation intervals. PLoS Negl Trop Dis. 2017;11: e0005797.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Sowilem MM, Kamal HA, Khater EI. Life table characteristics of Aedes aegypti (Diptera:Culicidae) from Saudi Arabia. Trop Biomed. 2013;30:301–14.

    PubMed  Google Scholar 

  39. Team RC. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2018.

  40. Hyndman R, Khandakar Y. Automatic time series forecasting: the forecast package for R. J Stat Softw. 2008;26:1–22.

    Google Scholar 

  41. Hyndman R, Athanasopoulos G, Bergmeir C, Caceres G, Chhay L, O’Hara-Wild M, et al. forecast: Forecasting functions for time series and linear models. 2023.

  42. ccf function - RDocumentation. https://www.rdocumentation.org/packages/tseries/versions/0.1-2/topics/ccf. Accessed 12 Apr 2023.

  43. Tran B-L, Tseng W-C, Chen C-C, Liao S-Y. Estimating the threshold effects of climate on dengue: a case study of Taiwan. Int J Environ Res Public Health. 2020;17:1392.

    Article  PubMed  PubMed Central  Google Scholar 

  44. PAHO PAHO. PAHO/WHO Data—Dengue. https://www3.paho.org/data/index.php/en/mnu-topics/indicadores-dengue-en.html. Accessed 27 Mar 2022.

  45. Perez F, Llau A, Gutierrez G, Bezerra H, Coelho G, Ault S, et al. The decline of dengue in the Americas in 2017: discussion of multiple hypotheses. Trop Med Int Health. 2019;24:442–53.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Borchering RK, Huang AT, Mier-y-Teran-Romero L, Rojas DP, Rodriguez-Barraquer I, Katzelnick LC, et al. Impacts of Zika emergence in Latin America on endemic dengue transmission. Nat Commun. 2019;10:5730.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Center NCP. NOAA’s Climate Prediction Center. https://origin.cpc.ncep.noaa.gov/products/analysis_monitoring/ensostuff/ONI_v5.php. Accessed 12 Apr 2023.

  48. Vincenti-Gonzalez MF, Tami A, Lizarazo EF, Grillet ME. ENSO-driven climate variability promotes periodic major outbreaks of dengue in Venezuela. Sci Rep. 2018;8:5727.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Ferreira HDS, Nóbrega RS, da Brito PVS, Farias JP, Amorim JH, Moreira EBM, et al. Impacts of El Niño Southern Oscillation on the dengue transmission dynamics in the Metropolitan Region of Recife, Brazil. Rev Soc Bras Med Trop. 2022;55:e0671.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Hospitales. Servicio Nacional de Salud. https://sns.gob.do/descarga-documentos/hospitales/. Accessed 22 Sep 2022.

  51. Bennett KL, Gómez Martínez C, Almanza A, Rovira JR, McMillan WO, Enriquez V, et al. High infestation of invasive Aedes mosquitoes in used tires along the local transport network of Panama. Parasites Vectors. 2019;12:264.

    Article  PubMed  PubMed Central  Google Scholar 

  52. González MA, Rodríguez-Sosa MA, Vásquez-Bautista YE, Rosariodel EC, Durán-Tiburcio JC, Alarcón-Elbal PM. A survey of tire-breeding mosquitoes (Diptera: Culicidae) in the Dominican Republic: considerations about a pressing issue. Biomedica. 2020;40:507–15.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Mohammadi A, Mostafavi E, Zaim M, Enayati A, Basseri HR, Mirolyaei A, et al. Imported tires; a potential source for the entry of Aedes invasive mosquitoes to Iran. Travel Med Infect Dis. 2022;49: 102389.

    Article  PubMed  Google Scholar 

  54. Khalid B, Bueh C, Ghaffar A. Assessing the factors of dengue transmission in urban environments of Pakistan. Atmosphere. 2021;12:773.

    Article  Google Scholar 

  55. Barcellos C, Lowe R. Expansion of the dengue transmission area in Brazil: the role of climate and cities. Tropical Med Int Health. 2014;19:159–68.

    Article  Google Scholar 

  56. de Playa Dorada, Asociacion de Hoteles. Overview and Hotspots Analysis of the Tourism Value Chain in Dominican Republic. 2019.

  57. Kingston R, Routledge I, Bhatt S, Bowman LR. Novel epidemic metrics to communicate outbreak risk at the municipality level: dengue and Zika in the Dominican Republic. Viruses. 2022;14:162.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Monintja TCN, Arsin AA, Amiruddin R, Syafar M. Analysis of temperature and humidity on dengue hemorrhagic fever in Manado Municipality. Gac Sanit. 2021;35:S330–3.

    Article  PubMed  Google Scholar 

  59. Sumi A, Telan EFO, Chagan-Yasutan H, Piolo MB, Hattori T, Kobayashi N. Effect of temperature, relative humidity and rainfall on dengue fever and leptospirosis infections in Manila, the Philippines. Epidemiol Infect. 2017;145:78–86.

    Article  CAS  PubMed  Google Scholar 

  60. Sirisena P, Noordeen F, Kurukulasuriya H, Romesh TA, Fernando L. Effect of climatic factors and population density on the distribution of dengue in Sri Lanka: a GIS based evaluation for prediction of outbreaks. PLoS ONE. 2017;12: e0166806.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Brady OJ, Johansson MA, Guerra CA, Bhatt S, Golding N, Pigott DM, et al. Modelling adult Aedes aegypti and Aedes albopictus survival at different temperatures in laboratory and field settings. Parasites Vectors. 2013;6:351.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Tjaden NB, Thomas SM, Fischer D, Beierkuhnlein C. Extrinsic incubation period of dengue: knowledge, backlog, and applications of temperature dependence. PLoS Negl Trop Dis. 2013;7: e2207.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Freitas A, Rodrigues HS, Martins N, Iutis A, Robert MA, Herrera D, et al. Multiplicative mixed-effects modelling of dengue incidence: an analysis of the 2019 outbreak in the Dominican Republic. Axioms. 2023;12:150.

    Article  Google Scholar 

  64. Petrone ME, Earnest R, Lourenço J, Kraemer MUG, Paulino-Ramirez R, Grubaugh ND, et al. Asynchronicity of endemic and emerging mosquito-borne disease outbreaks in the Dominican Republic. Nat Commun. 2021;12:151.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Navarro Valencia V, Díaz Y, Pascale JM, Boni MF, Sanchez-Galan JE. Assessing the effect of climate variables on the incidence of dengue cases in the metropolitan region of Panama City. Int J Environ Res Public Health. 2021;18:12108.

    Article  PubMed  PubMed Central  Google Scholar 

  66. Adeola AM, Botai JO, Rautenbach H, Adisa OM, Ncongwane KP, Botai CM, et al. Climatic variables and malaria morbidity in mutale local municipality, South Africa: a 19-year data analysis. Int J Environ Res Public Health. 2017;14:1360.

    Article  PubMed  PubMed Central  Google Scholar 

  67. Gharbi M, Quenel P, Gustave J, Cassadou S, Ruche GL, Girdary L, et al. Time series analysis of dengue incidence in Guadeloupe, French West Indies: forecasting models using climate variables as predictors. BMC Infect Dis. 2011;11:166.

    Article  PubMed  PubMed Central  Google Scholar 

  68. Ndii MZ, Anggriani N, Messakh JJ, Djahi BS. Estimating the reproduction number and designing the integrated strategies against dengue. Results Phys. 2021;27: 104473.

    Article  Google Scholar 

Download references

Acknowledgements

The authors would like to thank Albert Rodriguez, Heyliana Marte, and Pedro Vegas for their contributions to this project.

Funding

This project was supported by the Fund for Innovation and Scientific and Technological Development—Ministry of Higher Education, Science and Technology of the Dominican Republic. MAR was supported in part by a Burroughs Wellcome Fund Climate and Health Interdisciplinary Award (BWF #1022621).

Author information

Author notes

  1. Michael A. Robert and Helena Sofia Rodrigues contributed equally to this work

Authors and Affiliations

  1. Department of Mathematics and Center for Emerging, Zoonotic, and Arthropod-Borne Pathogens (CeZAP), Virginia Polytechnic Institute and State University, Blacksburg, VA, USA

    Michael A. Robert

  2. Escola Superior de Ciências Empresariais, Instituto Politécnico de Viana do Castelo, Valença, Portugal

    Helena Sofia Rodrigues

  3. Centro de Investigação e Desenvolvimento em Matemática e Aplicações, Universidade de Aveiro, Aveiro, Portugal

    Helena Sofia Rodrigues

  4. Centro de Investigación en Salud Dr. Hugo Mendoza, Hospital Pediátrico Dr. Hugo Mendoza, Santo Domingo, Dominican Republic

    Demian Herrera

  5. Departamento de Medicina Preventiva y Salud Pública y Microbiología, Facultad de Medicina, Universidad Autónoma de Madrid, Madrid, Spain

    Juan de Mata Donado Campos & Javier Del Águila Mejía

  6. Instituto de Investigación Sanitaria del Hospital Universitario La Paz (IdiPAZ), Universidad Autónoma de Madrid, Madrid, Spain

    Juan de Mata Donado Campos

  7. Consorcio de Investigación Biomédica en Red de Epidemiología y Salud Pública (CIBERESP), Instituto de Salud Carlos III, Calle Monforte de Lemos 3-5, 28029, Madrid, Spain

    Juan de Mata Donado Campos

  8. Departamento de Informática y Automática, Escuela Técnica Superior de Ingeniería Informática, Universidad Nacional de Educación a Distancia, Madrid, Spain

    Fernando Morilla

  9. Instituto Tecnológico de Santo Domingo (INTEC), Santo Domingo, Dominican Republic

    María Elena Guardado & Manuel Colomé-Hidalgo

  10. Dirección General de Epidemiología, Ministerio de Salud, Santo Domingo, Dominican Republic

    Ronald Skewes

Authors
  1. Michael A. Robert
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Helena Sofia Rodrigues
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Demian Herrera
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Juan de Mata Donado Campos
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Fernando Morilla
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Javier Del Águila Mejía
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. María Elena Guardado
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Ronald Skewes
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Manuel Colomé-Hidalgo
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

MAR, HSR, DH, and MC-H conceived and designed the study. MAR and HSR conducted statistical analyses. MAR and HSR drafted the manuscript. All authors contributed to interpretation of the data and revisions of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Michael A. Robert.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

: Table S1. Correlations in lags between dengue cases and climate variables not included in the main text. Lags are given as the number of weeks prior to dengue cases. Lags are listed with correlations in parentheses. Stars indicate confidence levels for testing significance: *** p<.01, ** p< .05, *p <.10.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Robert, M.A., Rodrigues, H.S., Herrera, D. et al. Spatiotemporal and meteorological relationships in dengue transmission in the Dominican Republic, 2015–2019. Trop Med Health 51, 32 (2023). https://doi.org/10.1186/s41182-023-00517-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s41182-023-00517-9

Keywords

Tropical Medicine and Health

ISSN: 1349-4147

Contact us