Effect of human mobility in Dengue spreading: Study cases for Caldas (CO) ========================================================================= * Carolina Ospina-Aguirre * David Soriano * Gerard Olivar-Tost * Cristian. C. Galindo-González * Jesús Gómez-Gardeñes * Gustavo Osorio ## Abstract According to the World Health Organization (WHO), dengue is the most common acute arthropod-borne viral infection in the world. The spread of dengue and other infectious diseases is closely related to human activity and mobility. In this paper we analyze the effect on the total number of dengue cases within a population after introducing mobility restrictions as a public health policy. To perform the analysis, we use a complex metapopulation in which we implement a compartmental propagation model coupled with the mobility of individuals between the patches. This model is used to investigate the spread of dengue in the municipalities of Caldas (CO). Two scenarios corresponding to different types of mobility restrictions are applied. In the first scenario, the effect of restricting mobility is analyzed in three different ways: a) limiting the access to the endemic node but allowing the movement of its inhabitants, b) restricting the diaspora of the inhabitants of the endemic node but allowing the access of outsiders, and c) a total isolation of the inhabitants of the endemic node. In this scenario, the best simulation results are obtained when endemic nodes are isolated during a dengue outbreak, obtaining a reduction of up to 22.51% of dengue cases. Finally, the second scenario simulates a total isolation of the network, *i.e*., mobility between nodes is completely limited. We have found that this control measure reduces the number of total dengue cases in the network by up to 42.67%. **Author summary** For the World Health Organization, dengue is a disease of public health concern. In recent years there is an increasing trend in the number of dengue cases despite existing prevention and control campaigns. The mobility of the population is considered an important factor in dengue dispersion. In this paper, we are interested in addressing how restrictions to human mobility might affect the incidence of dengue in a region. Our research is relevant because the model can be adapted to other regions or scales, and the mobility control measures can be taken into account for the development of public health policies in endemic regions. ## 1 Introduction According to the World Health Organization (WHO), dengue fever is the acute arthropod-borne viral infection with the highest incidence in humans in the world [1]. It is transmitted to humans mainly by the bite of a mosquito *Aedes aegypti* infected with Dengue virus (DENV). As a vector-borne disease the transmission between humans only occurs through the bite of infected mosquitoes, never from one person to another. In other words, the mosquito does not cause the disease directly, but acts as a bridge between two people, one with the virus and the other without it. In turn, the mosquito becomes infected when it feeds on the blood of a person infected with dengue and transmits the virus when it bites healthy people. Mathematical models are extremely useful to understand mechanisms that drive a healthy population toward an epidemic or endemic state, as well as for evaluating containment measures that help to suppress, or at least mitigate, the incidence of a given communicable disease. The usual mathematical approach to such tasks is the use of compartmental models. In these models, individuals in a population can be divided into classes or compartments according to their epidemiological state. For instance, in the celebrated Susceptible-Infectious-Recovered (SIR) model, individuals are divided into susceptible (healthy people who may acquire the virus), infectious (people who have acquired and can transmit the virus) and recovered (people who cannot propagate the pathogen and have acquired immunity to the virus) [2–4]. Yet simple and minimal, the SIR model has been pervasively used to analyse a plethora of viral infections such as measles [5, 6], rubella [7], malaria [8], zika [9], COVID-19 [10, 11], dengue [12–18], among others. In the case of vector-borne disease, more refined compartmental models have been introduced in which both human and vector populations are divided into several compartments. For instance in [19, 20] dynamics of human population is modeled through a SIR model whereas vectors are divided into Susceptible and Infectious. This division creates a SIR-SI model which is the compartmental dynamics adopted in our work. Alternative approaches for the study of the transmission of DENV include the Ross-Macdonald epidemic model [21] or the addition of further compartments to capture the growth dynamics of the vector population [22, 23]. In all these approaches, the main goal is to define control strategies that favor the eradication of the virus in humans population in a population. Most of control strategies rely on improving hygiene measures, the use of pesticides or, more recently, the release of Wolbachia-infected vectors in high incidence habitats. However, human behavior plays an important role in the geographical spread of pathogens [24, 25], and vector-borne diseases are not the exception. One of the most salient features of human behavior affecting the spread of vector-borne diseases is human mobility [26]. For the particular case of dengue transmission, mobility determines the degree of exposure to the disease vectors [27, 28] and, evidently, it is essential to foster local dengue outbreaks in low incidence areas through the importation of cases from distant and high incidence regions [29, 30]. A remarkable recent example of such phenomenon was the reintroduction of the dengue virus in Singapore [31]. The contribution of human mobility in the spread of dengue cases has become an essential factor given the expansion of urban environments and the increased frequency of international travels [32–34]. It has been projected that by 2030 more than half of the world’s population will reside in urban areas in the tropics, due to population growth and migration from rural areas [35]. The importance of human behavior and, in particular, human mobility in the spread of communicable diseases has motivated the formulation of different mathematical frameworks to study the contribution of mobility in the spread of vector-borne diseases [12, 36]. Most of these approaches bridge the gap between single-population models to metapopulations by incorporating the complex architecture of human flows in the form of networks [37] defined through data-driven frameworks [38–40]. Equipped with this framework, we focus on the analysis of how mobility constraints affect dengue transmission considering two control scenarios. In the first, mobility constraints are concentrated in the nodes with the highest incidence. In scenario, the entire network is connected and two types of mobility restrictions are analyzed. The first measure is to allow residents to leave but prevent outsiders from entering. In this scenario, preventing outsiders from entering and allowing residents to leave a high-incidence municipality was found to be the most effective measure, resulting in a 22.51% reduction in the number of dengue cases in the entire region. The second is to allow outsiders to enter a high-incidence municipality, but not to allow its residents to leave. And finally isolation consists of preventing residents from leaving and outsiders from entering a municipality. The second scenario analyzes the effect of confinement, *i.e*., totally restricting mobility throughout the network, preventing visitors from entering and residents from leaving a municipality. When isolation is applied, a reduction of cases of 42.67% is obtained compared to the unrestricted case. This article is organized as follows. Section 2, Materials and methods, describes the proposed compartmental model, the methodology used for estimating the model parameters and the coupling of the epidemiological model with the mobility network. Section 3 presents the results of the case study in which the model was applied in the department of Caldas - Colombia, where we analyze the effectiveness of the different mobility restrictions described above. The manuscript rounds off in section 4 with a brief discussion of the results and their implication for informing public health decision makers. ## 2 Materials and methods According to Satorras [37], we are currently witnessing a golden age in epidemic modeling: models are improving significantly thanks to the continuous addition of data, while the powerful computational resources available now make it possible to extend simulations to new limits. Throughout the years, SIR models have been used to analyse viral infections such as measles [5, 6], rubella [7], malaria [8], zika [9], COVID-19 [10, 11], dengue [12–18] and others. ### 2.1 SIR-SI compartmental model Let us denote the total size of the populations for humans and mosquitoes by *N**h* and *N**m*, respectively. Human population is divided into susceptible (*S**h*), infectious (*I**h*) and recovered (*R**h*), while mosquitoes population is divided into susceptible (*S**m*) and infectious (*I**m*) ones. To define the model, the following assumptions are made: * The size of the human population *N**h* is assumed to be constant, corresponding to the steady state on the underlying population dynamics, since the period of the disease is very small in relation with population dynamics. * Incubation period is neglected for both, humans and mosquitoes. * Deaths caused by the disease are not considered for either humans or mosquitoes. * Mosquitoes cannot recover after being infected. * There is not super-infection for either humans or mosquitoes. * Susceptible humans can only get infected trough an infected mosquito bite. * Human recruitment rate is constant. * Mosquito recruitment rate depends on environmental conditions and is constant in each region. The dynamics of dengue transmission between humans and mosquitoes, is defined by the system of equations (1). A susceptible mosquito can become infectious if it has contact with an infectious individual, according to the transmission probability of the vector *λ**m*, the bite rate *β* and the proportion of infectious humans in the node ![Graphic][1]. Therefore, the infection rate for susceptible mosquitoes is given by ![Graphic][2]. For the case of humans, the rate of infection is proportional to the biting rate *β* and the number of infectious mosquitoes *I**M*, transmission probability *λ**h* and inversely proportional to the total population *N**H*, determining the probability of a given individual being bitten. Therefore, the infection rate for susceptible humans is then expressed as ![Graphic][3]. Considering the former mechanisms, the dynamical evolution of each compartment associated to a given municipality is described by the following equations: ![Formula][4] with *N**h* = *S**h* + *I**h* + *R**h*. In table 1 we specify the definition of the parameters used in the former equations of the SIR-SI model. View this table: [Table 1.](http://medrxiv.org/content/early/2023/01/12/2023.01.10.23284416/T1) Table 1. Parameters description for the SIR-SI model. The dynamics of the infectious set in the SIR model according to *Allen* [41] can yield two different behaviors, one in which there is an epidemic peak and the other in which the number of infected individuals decreases monotonically until the end of the outbreak. In [42] the analysis developed by *Allen* is extended to the SIR-SI model (1) and it is shown that the dynamics for the infected population depends on ![Formula][5] When *ρ >* 1, a maximum occurs (epidemic peak) and then decreases to zero. When these conditions generate *ρ <* 1 the infected have a decreasing behavior until they reach zero. #### 2.1.1 Parameters estimation and initialization The parameters governing contagion dynamics and the mosquito’s life-cycle are extracted from the work published by Helmersson *et al*. [43]. In this manuscript, the authors proposed different equations capturing the reduction of vectorial capacity of *Aedes aegypti* when being exposed to either very high (*T >* 34°*C*) or very low (*T <* 12°*C*) temperatures. The behavior of the parameters as a function of temperature can be seen in the Fig 1. In particular, the average bite rate reads: ![Formula][6] the probability that a mosquito becomes infectious after biting and infectious human is: ![Formula][7] instead, in the range 12.286 ≤ *T* ≤ 32.461, the infection probability for humans can be modeled as: ![Formula][8] while the mosquito death rate goes as: ![Formula][9] for 10.54 ≤ *T* ≤ 33.41. In Fig. 1 shows the variation of these four parameters as the temperature increases from 8°*C* to 40.7°*C* ![Fig 1.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/01/12/2023.01.10.23284416/F1.medium.gif) [Fig 1.](http://medrxiv.org/content/early/2023/01/12/2023.01.10.23284416/F1) Fig 1. Temperature dependence of the parameters. (a) Bite rate *β*, (b) Rate of death of mosquitoes *δ**m*, (c) Likelihood of a person becoming infectious *λ**h*, (d) Likelihood of a mosquito becoming infectious *λ**m*. The other four parameters are defined as follows. For Λ*h*, the approach proposed by [44] is used, where it is defined as ![Graphic][10]. The recovery rate for humans *μ* is set to 0.32288 according to the study developed by Hamdan [14], while *δ**h* is estimated using known specific statistics from the geographical region under analysis. Finally, to estimate Λ*m* we start from the equation for the vital dynamics of the vector when *S**m* = 0, (*i.e*. ![Graphic][11]), where the carrying capacity (number of mosquitoes at steady state) is obtained as ![Graphic][12]. The recruitment rate, therefore, can be expressed in terms of the carrying capacity ![Graphic][13]. To set the initial conditions for the model, some considerations have to be mentioned. First of all, it has to be noted that *N**h* can be fixed based on the known population of each municipality, while *N**m* can be estimated with the Index of Domiciliary Infestation (number of houses with one or more containers positive for immature *Aedes aegypti* divided by the number of houses sampled multiplied by 100) and the average number of people per household in the study region. Considering Eqs. (1) and the parametrization described above we shown in Fig 2 the behavior of the model when simulating a region with an average temperature of 31°*C* (*β* = 0.227 *λ**m* = 1, *δ**m* = 0.0298, *λ**h* = 0.732, *μ* = 0.329) and the mosquito recruitment rate Λ*m* = 86.634. When the initial conditions are *N**m* = 148.352, *S**m* = 118.352, *I**m* = 30, *N**h* = 76963 *S**h* = 76963 *I**h* = 0, that produces a *ρ >* 1 and an epidemic peak. When the initial conditions are *N**m* = 148.352, *S**m* = 118.352, *I**m* = 0, *N**h* = 76963, *S**h* = 76928, *I**h* = 35, that results in a *ρ <* 1 and thus the number of infected decreases to zero. ![Fig 2.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/01/12/2023.01.10.23284416/F2.medium.gif) [Fig 2.](http://medrxiv.org/content/early/2023/01/12/2023.01.10.23284416/F2) Fig 2. Behavior of the human population over time in the SIR-SI model. (a,c) Susceptible, (b,e) Infected and (c,f) Recovered. When *ρ <* 1 (a, b and c), and *ρ >* 1 (d,e and f). ### 2.2 SIR-SI model with human mobility The SIR-SI can be modified to account for the impact of recurrent human mobility patterns. For this purpose, let us assume that each individual has his/her residence located inside a given patch but can move to other destination, identified as the workplace for instance. In the absence of data, we assume these movements to be governed by a fully-connected weighted matrix ϒ. We construct this matrix synthetically by assuming that 90% of the population inside each patch remain there whereas 10% of the population move to another area. We compute the flows connecting different patches following the gravity model [39], implying that the number of connections between two patches *i* and *j*, hereinafter denoted by *W**ij*, can be expressed as: ![Formula][14] Taking into account the flows distribution, the elements of the mobility matrix ϒ read: ![Formula][15] Agents’ movements between patches *i* and *j* cause a redistribution of the population across the system; therefore it is necessary to adjust the equations (1) to account for the actual number of people that are in any patch *i* at any given time. In particular, the effective population of patch *i* is defined as: ![Formula][16] which accounts for the distribution of the residential population ![Graphic][17] and the mobility patterns of the individuals of the metapopulation ϒ. Likewise, mobility also changes the effective number of infected individuals in each patch *i*, which now reads: ![Formula][18] Finally, we assume that mosquitoes stay inside their associated area. To model the spatio-temporal evolution of the disease, we define the quantity *x**i*(*t*) as the occupation of each compartment (*x* ∈ {*S**h*, *I**h*, *R**h*, *S**m*, *I**m*}) inside each patch *i*. Following the assumptions of the model, the time evolution of these quantities is given by: ![Formula][19] In Fig. 3 we show an scheme of the metapopulation model in which the SIR-SI compartmental model is integrated with human mobility flows. ![Fig 3.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/01/12/2023.01.10.23284416/F3.medium.gif) [Fig 3.](http://medrxiv.org/content/early/2023/01/12/2023.01.10.23284416/F3) Fig 3. Schematic view of the SIR-SI metapopulation model. (a) shows the compartmental dynamics of the SIR-SI model with cross-infections between vectors and humans while panel (b) shows a toy metapopulation of 3 patches connected through links whose weights are given by matrix ϒ. ## 3 Results Once presented the general formalism to tackle the dissemination of dengue fever across a metapopulation we now apply it to particular scenarios. We first address the analysis of three simple cases that illustrate how mobility between a small set of patches have profound implications in the overall extent of dengue cases. Then, we present our main case study: the spread of dengue fever in the department of Caldas (CO). We describe the main attributes of the department, their map into the metapopulation framework, and the parametrization of the model according to observed dengue cases. Finally, this section presents the results of two different mobility restriction and show their impact on the total dengue fever cases. ### 3.1 Mobility effects over dengue cases in simple metapopulations We start by developing the analysis of three very simple metapopulations to unveil the effect of mobility on the spread of dengue. To this aim, we define three case examples in which the distance between nodes is set as a control parameter so that its increase implies, according to the gravity model, the decrease of mobility. The three scenarios are sketched in Fig. 4. ![Fig 4.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/01/12/2023.01.10.23284416/F4.medium.gif) [Fig 4.](http://medrxiv.org/content/early/2023/01/12/2023.01.10.23284416/F4) Fig 4. Schematic plot of the three simple metapopulations covered as case examples. In panel (a), corresponding to case A, node 1 is a zero-incidence node, i.e., there are no mosquitoes on it, however in node 2, there are environmental conditions for vector breeding and the disease is present. Panel (b) corresponds to case B. In this case node 1 is defined as endemic and is connected to a non-endemic node. Note that in both nodes the disease is present, but in node 1 the vector breeding conditions are better, so there is a wider spread of the disease than in node 2. Finally in panel (c) we define case C for which there are three equidistant connected nodes. Node 1 has zero incidence and its population is the largest while nodes 2 and 3 have characteristics of endemic regions. However, node 2 has more inhabitants than node 3 and, therefore, the probability of infection in node 3 is larger than in the other two patches. #### 3.1.1 Case A As sketched in Fig. 4.a this case example corresponds to a patch (node 1) with zero incidence connected to an endemic patch (node 2). The model parameters used to simulate this example are shown in Table 10(a). In Fig 5a, we show the evolution of dengue cases in both patches as the distance increases. This graph shows that increasing the distance between patches, i.e. reducing mobility between patches, leads to an overall reduction in the total number of cases. Note however that the effects of increasing distance (reducing mobility) for both patches are opposite. In particular, while increasing mobility increase the exposure of individuals from non-endemic areas leading to a higher incidence, movement across patches is beneficial for the population in endemic areas, decreasing the total incidence of the virus on this population. ![Fig 5.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/01/12/2023.01.10.23284416/F5.medium.gif) [Fig 5.](http://medrxiv.org/content/early/2023/01/12/2023.01.10.23284416/F5) Fig 5. Evolution of dengue cases in patches as the distance increases when connecting a) High and a zero incidence node. b) High and medium incidence node. c) High, medium and zero incidence node. The green line shows the denge cases in the zero incidence node, the blue line shows the cases in the medium incidence node and the red line shows the cases in the high incidence node. The black line is used for the total number of cases. #### 3.12 Case B In this example, node 1 is a patch with high disease presence and node 2 a patch where there are fewer dengue cases (see Table 11(a) for the parameters used in this case). In Fig. 5b, we observe a qualitatively similar pattern to the one described before but much less pronounced due to the higher similarity among the connected patches. #### 3.1.3 Case C In this example three equidistant nodes are connected. According to the parameters chosen in Table 12, node 1 has zero incidence and the largest population, node 2 has a larger population than node 3, but the latter has a higher incidence of the disease (97.5 cases per thousand inhabitants versus 75.2 for node 2) when patches are isolated. The results of this case are shown in Fig. 5 Case C. Again geographical location, *i.e*. the increase of the distance between patches, allows for a reduction in the total number of cases. However, locally,the situation of node 3 (zero incidence) worsens when the distance is reduced, while for the medium and high incidence nodes, proximity reduces the number of dengue cases by up to 50%. The results of the former cases clearly show that mobility restrictions have a great impact on reducing the number of cases, although locally the epidemiological situation can worsen as a byproduct of this contention measures. Motivated by these results we now propose scenarios to evaluate the effect of applying mobility restrictions in a more complex metapopulation inspired in the municipalities of department of Caldas (CO). ### 3.2 Defining the metapopulation of the department of Caldas The department of Caldas is located in the central west of the Colombian Andean region. It is crossed by the Central and Western mountain ranges, which makes it a department of mountainous topography. It has an area of 7888 *km*2 and a population of 987991 inhabitants [45], while its elevation above sea level ranges between 5400 masl and 170 masl. The department has 27 municipalities, ten of which have direct communication with 6 other municipalities in neighboring departments. To apply our SIR-SI model to the department of Caldas, we define a metapopulation with 33 patches, each with a number of residents equal to that reported by the census. The connections (links) between patches are calculated as a function of the road distance between them, resulting in the ϒ matrix containing the transition rates between pairs of municipalities. This matrix does not take into consideration whether the municipalities belong to the department of Caldas or they are external municipalities that by proximity have a relationship with the previous ones. The precise transition rates between the 33 municipalities (patches) are reported in Table 5). For the numerical simulations, the parameters *β, λ**h*, *λ**m* and *δ**m* are calculated with the procedure mentioned in section 2.1.1 and shown in Table 3 for each of the analyzed municipalities. The estimations are based on the average temperature value in each of them (see Table 2). In turn, the recruitment rate for humans is calculated with the value of life expectancy (*EV H* = 74.64), defined in public health records whereas humans mortality rate is obtained from DANE statistics. View this table: [Table 2.](http://medrxiv.org/content/early/2023/01/12/2023.01.10.23284416/T2) Table 2. Average altitude and temperature of the municipalities of the department of Caldas. Source: Caldas Health Observatory Indicator Booklets. View this table: [Table 3.](http://medrxiv.org/content/early/2023/01/12/2023.01.10.23284416/T3) Table 3. Estimated values of the model parameters. With the demographic, epidemiological and mobility parameters, initial conditions were set for each municipality (presented in Table 4). To calibrate the model to capture the number of dengue cases historically occurred in the department of Caldas (CO) during some outbreaks in the last decade, periodic perturbations in the number of infected mosquitoes during two rainy seasons in a year were introduced. As a necessary simplification, we only considered one dengue serotype; this assumption implies that individuals can only be infected once. The model was adjusted by minimizing the error between simulated and real cases. In particular, the SIR-SI simulation model with the 33 municipalities connected after calibration yields 1099 cases of dengue in Caldas (see Table 6). View this table: [Table 4.](http://medrxiv.org/content/early/2023/01/12/2023.01.10.23284416/T4) Table 4. Model initial conditions View this table: [Table 5.](http://medrxiv.org/content/early/2023/01/12/2023.01.10.23284416/T5) Table 5. Matrix ϒ. Each cell represents the probability that an inhabitant of a municipality moves to another one. View this table: [Table 6.](http://medrxiv.org/content/early/2023/01/12/2023.01.10.23284416/T6) Table 6. Dengue cases full conection View this table: [Table 7.](http://medrxiv.org/content/early/2023/01/12/2023.01.10.23284416/T7) Table 7. Dengue incidence rate per 100 000 inhabitants, by municipality of origin. Source: Public Health Observatory. The agreement between real cases and those obtained after the SIR-SI model has been calibrated is shown in Fig. 6. In particular the map on the left shows the dengue cases reported by the health entities and on the right the dengue cases estimated with the model. It can be seen at the bottom that there is a correlation of 0.76 between the data estimated with the model and the cases occurred in the department of Caldas during 2015. ![Fig 6.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/01/12/2023.01.10.23284416/F6.medium.gif) [Fig 6.](http://medrxiv.org/content/early/2023/01/12/2023.01.10.23284416/F6) Fig 6. On the left, cases of dengue fever occurred in the department of Caldas during one year (2015). On the right, dengue cases predicted by the model during the same year. Below correlation between both variables. ### 3.3 Control Scenario 1: Mobility restrictions focused on high incidence nodes The purpose with this scenario is to get some insights on the contribution to the total number of dengue cases in the network, after applying mobility constraints to the five most affected cities (one at a time). We analyze the variation of dengue cases according to mobility constraints when: i) Access: No foreigners are allowed to enter the municipality, ii) Exit: no inhabitants can leave the municipality, iii) Isolation: No entry or exit allowed for residents or visitors. These control measures were applied in simulations for the five municipalities with the highest number of dengue cases during during the period 2015 − 2019 according to the public health office of Caldas [46]. These municipalities are La Dorada, Norcasia, Marmato, Chinchiná and Manizales. Fig. 7 shows the variation of the total number of dengue cases in the department of Caldas when the three mobility control measures proposed in this work are applied in the five municipalities mentioned. In Fig.7 it can be seen that the total number of dengue cases in the department varies as mobility in each municipality is reduced. ![Fig 7.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2023/01/12/2023.01.10.23284416/F7.medium.gif) [Fig 7.](http://medrxiv.org/content/early/2023/01/12/2023.01.10.23284416/F7) Fig 7. Variation of dengue cases in the department Caldas according to mobility restrictions when: i) no foreigners are allowed to enter the municipality (blue), no inhabitants can leave the municipality (green), the municipality is totally isolated (orange). #### 3.3.1 Access This control measure, when applied in each of the five municipalities with the highest burden of the disease, reduces the total number of dengue cases in the department. In Fig 7 it can be seen that the most significant reduction occurs when this control measure is applied in La Dorada and Marmato, with a decrease of 159 and 210, respectively. These values indicate that completely restricting access to this endemic municipality has a positive effect on reducing the total number of dengue cases in the department. #### 3.3.2 Exit When the exit of the inhabitants of the municipality of La Dorada is restricted, but the entry of foreigners is allowed, the number of dengue cases increase in 20 with respect to the free mobility scenario. This same behavior occurs when the control measure is applied in Norcasia and Marmato where the increase of 9 and 5 dengue cases. When the control measure is applied in Manizales, dengue cases are reduced by 247. This is due to the fact that its inhabitants are infected when traveling to endemic areas. #### 3.3.3 Isolation Isolating a municipality in this study means limiting mobility to zero, therefore its inhabitants cannot travel to other municipalities and people are not allowed to enter. Applying this control measure in La Dorada, Norcasia, Chinchiná and Manizales reduces the total number of dengue cases in the department by 139, 90, 51 and 270, respectively. The opposite occurs when applied in Marmato, where dengue cases increase by 46. ### 3.4 Control scenario 2: Confinement During the COVID-19 pandemic, governmental entities with the objective of containing the spread of the virus have opted to isolate neighborhoods, even entire countries. In this work we simulate the effect of confining the municipalities in order to analyze how this restriction influences the total number of dengue cases. Confinement consists of limiting the exit of inhabitants and entry of foreigners to each municipality. When confinement is applied, 630 cases occur, which implies a reduction of 42.67% with respect to the 1099 dengue cases that occurred when there are no restrictions. The dengue cases before and after confinement are shown in Table 8, where it can be seen that mobility affects each municipality differently. The column (variation of cases) shows positive values when the mobility restriction reduces dengue cases and negative values when there is an inverse effect. It can be observed that in the municipalities of Manzanares, Salamina, Aguadas, Risaralda, Aranzazu, La Merced and Pácora, confining population reduces dengue cases to zero. Confinement reduces dengue cases in the municipalities of La Dorada, Victoria and Riosucio, where a higher percentage of travelers go to areas with a higher incidence of dengue. A similar phenomenon occurs in Norcasia, Marmato and Palestina, but in the opposite direction, with the highest percentage of travelers from these three municipalities going to areas with lower incidence. The percentages of travelers to municipalities of zero, lower and higher incidence than the municipality of residence can be seen in Table 9. The results obtained in this scenario coincide with those obtained by *Conceição* in his study conducted in Sao Paulo [47] where there was a reduction in dengue cases as a collateral result of population confinement to hinder the advance of COVID-19 pandemic. View this table: [Table 8.](http://medrxiv.org/content/early/2023/01/12/2023.01.10.23284416/T8) Table 8. Dengue cases with and without mobility restrictions. Negative number means increase and positive number represents reduce after applying the mobility control measure (quarantine). View this table: [Table 9.](http://medrxiv.org/content/early/2023/01/12/2023.01.10.23284416/T9) Table 9. Percentage of travelers to 0, lower and higher incidence areas than their city of residence. View this table: [Table 10.](http://medrxiv.org/content/early/2023/01/12/2023.01.10.23284416/T10) Table 10. Example A. View this table: [Table 11.](http://medrxiv.org/content/early/2023/01/12/2023.01.10.23284416/T11) Table 11. Example B. View this table: [Table 12.](http://medrxiv.org/content/early/2023/01/12/2023.01.10.23284416/T12) Table 12. Example C. ## 4 Discussion Several scenarios were analyzed to provide policy makers with insight into the implementation of control measures to mitigate the spread of dengue in the department of Caldas. The study used a compartmental model at each node of a network and investigated the effects on the total number of dengue cases when applying mobility control measures such as Access: No foreigners are allowed to enter the municipality, Isolation: No residents or visitors are allowed to enter or leave, and Exit: No inhabitant is allowed to leave the municipality. The results obtained explain phenomena that have been observed in the historical data, and the coexistence of different types of trends according to local conditions and connectivity. Such is the case of the capital Manizales (zero incidence since environmental conditions are not suitable for the mosquito), which has reported cases of the disease in the health system. The inhabitants of this municipality contract the disease when moving to areas where dengue fever is present. The simulations show this phenomenon and it is observed that limiting mobility from nodes of zero incidence to nodes with the presence of the disease prevents the spread of the disease. The results obtained in this study are similar to those obtained by Jassen, who showed that bans on international air travel during the COVID-19 pandemic led to a sharp decrease in dengue cases in Australia, where the disease is predominantly imported [48]. Another phenomenon observed was that the isolation of Marmato (a node with a high incidence of the disease) generates an increase in dengue cases. Given that when human movements are limited to one’s own home and its surroundings, contact between people and vectors may increase, resulting in a higher risk of exposure and transmission of the virus, as occurred in Thailand where the isolation generated by COVID-19 generated an increase in dengue cases [49]. On the contrary, when residents of Marmato travel to nearby municipalities with a lower incidence of the disease, the probability of contagion is reduced with displacement. In the case of La Dorada, which is also a municipality with a high incidence, the increase in cases is generated by the departure of its inhabitants. This is due to the fact that most people move to areas with a higher incidence of the disease, which increases the probability of infection and, therefore, the number of infected people. This phenomenon was evidenced in Singapore where the mobility of inhabitants was limited to public spaces, such as workplaces or schools, where the density of mosquitoes was high, which led to a decrease in dengue transmission during the COVID-19 quarantine [50]. Therefore, in La Dorada or any municipality with these mobility characteristics, the control measure that generates the greatest impact on the reduction of dengue cases is to limit the exit of its inhabitants. In this study it was observed that the reduction of dengue cases after applying quarantine is not significantly higher than applying other measures, which is why it is suggested to control entities analyze the relevance of this measure when making decisions since, as has been observed in the current pandemic, isolation has a negative influence on people. It causes alterations in eating habits [51], increases the use of drugs [52], has negative psychological [53, 54], and economic [55] effects on people around the world. The analysis carried out in this study can be used, both in the department of Caldas as well as in any other region, to design mobility reduction strategies to control or mitigate a dengue outbreak. ## 5 Conclusion The model presented in this study makes it possible to analyze the effect of applying control measures to reduce dengue cases in a region. By comparing the results of the model with the actual cases, it has been established that the estimated values have a realistic correspondence. Based on the results, it is recommended that health experts and administrative authorities limit mobility to areas with high incidence of the disease during dengue outbreaks in order to prevent its spread. This model can be extrapolated to any region, if the average temperature, altitude, number of inhabitants, human vital statistics and aedic records are known. Moreover, the mobility control measures here discussed can be taken into account for the development of public health policies in endemic regions. ## Data Availability All relevant data are within the manuscript and its Supporting Information files. ## 7 Supporting information ## 6 Acknowledgments We acknowledge the Health observatory of Caldas, Territorial of Health of Caldas and National Coffee Research Center - Cenicafé for allowing access to the database and departmental information for the development of this study. Carolina Ospina acknowledges funding from Colciencias “convocatoria doctorados nacionales 647 de 2014”. We acknowledge financial support from grant PID2020-113582GB-I00 funded by MCIN/AEI/10.13039/501100011033, and from Departamento de Industria e Innovación del Gobierno de Aragón y Fondo Social Europeo through projects no. E36 17R (FENOL group). * Received January 10, 2023. * Revision received January 10, 2023. * Accepted January 12, 2023. * © 2023, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. 1.Scientific Working Group on Dengue Meeting (2006: Geneva S, for Research UBSP, in Tropical Diseases T. Report of the Scientific Working Group meeting on Dengue, Geneva, 1–5 October 2006; 2007. 2. 2.Anderson RM, May RM. Infectious diseases in humans. Oxford University Press, Oxford; 1992. 3. 3.Keeling M, Rohani P. Modeling Infectious Diseases in Humans and Animals. Princeton University Press,Princeton; 2007. 4. 4.Angulo Salazar RC, Díaz Cervo Y, Pardo Pinzón R. Índice de Pobreza Multidimensional para Colombia (IPM-Colombia) 1997-2010. Archivos de Economía. 2011;382. 5. 5.Zhou L, Wang Y, Xiao Y, Li MY. Global dynamics of a discrete age-structured SIR epidemic model with applications to measles vaccination strategies. Mathematical Biosciences. 2019;308:27–37. 6. 6.Bistritz I, Bambos N, Kahana D, Ben-Gal I, Yamin D. Controlling Contact Network Topology to Prevent Measles Outbreaks. In: 2019 IEEE Global Communications Conference (GLOBECOM); 2019. p. 1–6. 7. 7.1. Jezic G, 2. Chen-Burger YHJ, 3. Kusek M, 4. Šperka R, 5. Howlett RJ, 6. Jain LC, editors Kurahashi S. An Agent-Based Infectious Disease Model of Rubella Outbreaks. In: Jezic G, Chen-Burger YHJ, Kusek M, Šperka R, Howlett RJ, Jain LC, editors. Agents and Multi-agent Systems: Technologies and Applications 2019. Singapore: Springer Singapore; 2020. p. 237–247. 8. 8.Affandi P, Faisal. Optimal control mathemathical SIR model of malaria spread in South Kalimantan. Journal of Physics. 2018;. 9. 9.Bonyah E, Khan M, Okosun K, Islam S. A theoretical model for Zika virus transmission. PLoS ONE. 2017;12(10). 10. 10.Giordano G, Blanchini F, Bruno R, Colaneri P, Di Filippo A, Di Matteo A, et al. Modelling the COVID-19 epidemic and implementation of population-wide interventions in Italy.Nature Medicine. 2020;. 11. 11.Biswas K, Khaleque A, Sen P. Covid-19 spread: Reproduction of data and prediction using a SIR model on Euclidean network. arXiv e-prints. 2020;. 12. 12.Zhu D, Ren J, Zhu H. Spatial-temporal basic reproduction number and dynamics for a dengue disease diffusion model. Mathematical Methods in the Applied Sciences. Science of the Total Environment. 2018;41(14):5388–5403. 13. 13.Nur W, Rachman H, Abdal NM, Abdy M, Side S. IR Model Analysis for Transmission of Dengue Fever Disease with Climate Factors Using Lyapunov Function. Journal of Physics: Conference SeriesDocumento CEDE. 2018;1028. 14. 14.Hamdan NI, Kilicman A. A fractional order SIR epidemic model for dengue transmission. Chaos, Solitons Fractals. 2018;114:52–62. 15. 15.Khan MA, Khan A, Elsonbaty A, Elsadany AA. Modeling and simulation results of a fractional dengue model. The European Physical Journal Plus. 2019;134(8). 16. 16.Aguiar M, Stollenwerk N. SHAR and effective SIR models: from dengue fever toy models to a COVID-19 fully parametrized SHARUCD framework. Mathematical Modelling in Biosciences. 2020;3(1):60–89. 17. 17.Bouzid L, Belhamiti O, Belgacem FBM. Efect of human mobility on predective spatio-temporal model of denge epidemic transmission. Journal of Biological Systems. 2022;30(4):721–739. 18. 18.Chen Y, Liu T, Yu X, Zeng Q, Cai Z, Wu H, et al. An ensemble forecast system for tracking dynamics of dengue outbreaks and its validation in China. PLoS Comput Biol. 2022;18(6). 19. 19.Derouich M, Boutayeb A, Twizell E. A model of dengue fever. BioMed Central. 2003;. 20. 20.Derouich M, Boutayeb A. Dengue fever: Mathematical modelling and computer simulation. Applied Mathematics and Computation. 2006;17:528–544. 21. 21.Sepúlveda LS. Manejo óptimo y viable en modelos epidemiológicos del dengue. Universidad Autónoma de Occidente; 2015. 22. 22.Lopez Montenegro LE. Modelos y Control Optimo Poblacional del Aedes aegypti con Retardos de Tiempo. Universidad Nacional de Colombia sede Manizales; 2012. 23. 23.Reyna-Lara A, Soriano-Paños D, Arias-Castro J, Martínez H, Gómez-Gardeñes J. A metapopulation approach to identify targets for Wolbachia-based dengue control. Chaos: An Interdisciplinary Journal of Nonlinear Science. 2022;32(4):041105. 24. 24.Fenichel EP, Castillo-Chavez C, Ceddia MG, Chowell G, Gonzalez Parra PA, Hickling GJ, et al. Addaptive human behavior in epidemiological models. Proceedings of the National Academy of Sciences. 2011;108(15):6306–6311. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTA4LzE1LzYzMDYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wMS8xMi8yMDIzLjAxLjEwLjIzMjg0NDE2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 25. 25.Funk S, Gilad E, Watkins C, Jansen VAA. The spread of awareness and its impact on epidemic outbreaks. Proceedings of the National Academy of Sciences. 2009;106:6872–7. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTA2LzE2LzY4NzIiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wMS8xMi8yMDIzLjAxLjEwLjIzMjg0NDE2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 26. 26.Stoddard ST, Morrison AC, Vazquez Prokopec GM, Soldan VP, Kochel TJ, Kitron U, et al. The role of human movement in the transmission of vector-borne pathogens. PLoS Negl Trop Dis. 2009;3(7). 27. 27.Stoddard S, Forsheycde BM, Morrisona AC, az Soldanf VA, Vazquez-Prokopecb PGM, Astetec H, et al. House-to-house human movement drives dengue virus transmission. Proceedings of the National Academy of Sciences. 2013;110(3):994–999. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czo5OiIxMTAvMy85OTQiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wMS8xMi8yMDIzLjAxLjEwLjIzMjg0NDE2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 28. 28.Tao H, Liu Y, Wang K, Zhuo L. Assessing Impacts of Traffic Flows on the Spatial Distribution of Early Dengue in Guangzhou Subdistricts. In: IGARSS 2019 - 2019 IEEE International Geoscience and Remote Sensing Symposium; 2019. p. 3468–3470. 29. 29.Barrios E, Lee S, Vasilievaa O. Assessing the effects of daily commuting in two-patch dengue dynamics: A case study of Cali, Colombia. Journal of Theoretical Biology. 2018;453:14–39. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jtbi.2018.05.015&link_type=DOI) 30. 30.Kiang MV, Santillana M, Chen JT, Onnela JP, Krieger N, Engø-Monsen K, et al. Incorporating human mobility data improves forecasts of Dengue fever in Thailand. Scientific Reports. 2021;11(923). 31. 31.Ooi E, Gubler D. Dengue in Southeast Asia:epidemiological characteristics and strategic challenges in disease prevention. Cad Saude Publica. 2009;25:115–124. 32. 32.Gubler DJ. The global pandemic of dengue/dengue haemorrhagic fever: current status and prospects for the future. Ann Acad Med Singapore. 1998;27:227–234. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9663316&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F01%2F12%2F2023.01.10.23284416.atom) 33. 33.Russell RC, Currie BJ, Lindsay MD, Mackenzie JS, Ritchie SA, Whelan PI. Dengue and climate change in Australia: predictions for the future should incorporate knowledge from the past. Medical Journal of Australia. 2009;190:265–268. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19296793&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F01%2F12%2F2023.01.10.23284416.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000265400000015&link_type=ISI) 34. 34.Tian H, Sun Z, Faria NR, Yang J, Cazelles B, Huang S, et al. Increasing airline travel may facilitate co-circulation of multiple dengue virus serotypes in Asia. PLoS Negl Trop Dis. 2017;11(8). 35. 35.Seto KC, Guneralp B, Hutyra LR. Global forecasts of urban expansion to 2030 and direct impacts on biodiversity and carbon pools. Proc Natl Acad Sci USA. 2012;109(16083–16088). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTA5LzQwLzE2MDgzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDEvMTIvMjAyMy4wMS4xMC4yMzI4NDQxNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 36. 36.Bomfim R, Pei S, Shaman J, Yamana T, Makse H, Andrade Jr J, et al. Predicting dengue outbreaks at neighbourhood level using human mobility in urban areas. J R Soc. 2020;17(20200691). 37. 37.Pastor-Satorras R, Castellano C, Mieghem PV, Vespignani A. Epidemic processes in complex networks. Reviews of Modern Physics. 2015;87(925). 38. 38.Simini F, Gonzalez MC, Maritan A, Barabasi AL. A universal model for mobility and migration patterns. Nature. 2012;484(7392):96–100. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature10856&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22367540&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F01%2F12%2F2023.01.10.23284416.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000302343400041&link_type=ISI) 39. 39.Barbosa H, Barthelemy M, Ghoshal G, James CR, Lenormand M, Louail T, et al. Human mobility: Models and applications. Physics Reports. 2018;734:1–74. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.physrep.2018.01.001&link_type=DOI) 40. 40.Xia Y, Bjørnstad ON, Grenfell BT. SMeasles metapopulation dynamics: A gravity model for epidemiological coupling and dynamics. The American Naturalist. 2004;164:267–281. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/422341&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15278849&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F01%2F12%2F2023.01.10.23284416.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000223099300014&link_type=ISI) 41. 41.Allen JSL. An Introduction to Mathematical Biology. Pearson Prentice hall’; 2006. 42. 42.Ospina Aguirre C. Modelado de brotes epidémicos de dengue para la toma de decisiones en salud pública: Efecto de la movilidad en el departamento de Caldas [tesis de doctorado, Universidad Nacional de Colombia]. Repositorio, Universidad Nacional de Colombia sede Manizales; 2021. 43. 43.Liu-Helmersson J, Stenlund H, Wilder-Smith A, Rocklöv J. Vectorial Capacity of Aedes aegypti: Effects of Temperature and Implications for Global Dengue Epidemic Potential. 2014;9(3). 44. 44.Aldila D, Situngkir N, Nareswari K. Understanding resistant effect of mosquito on fumigation strategy in dengue control program. Journal of Physics: Conference Series. 2018;. 45. 45.DANE. Censo general de Colombia, 2005; 2009. Bogotá: Departamento Administrativo Nacional de Estadística;. 46. 46.Instituto nacional dS. Vigilancia rutinaria departamental. Ministerio de Salud y Protección Socia; 2015. 47. 47.Conceição GMdS, Barbosa GL, Lorenz C, Dias Bocewicz AC, Reis Santana LM, Marques CCdA, et al. Effect of social isolation in dengue cases in the state of Sao Paulo, Brazil: an analysis during the COVID-19 pandemic. Travel Med Infect Dis. 2021;44. 48. 48.CC J, JM D, FA B, MA S, AF vdH. Impact of COVID-19 Mitigation Measures on Mosquito-Borne Diseases in 2020 in Queensland, Australia. Viruses. 2021;13(6). 49. 49.Jt L, Bsl D, Lzx C, Elw C, JR K, J A. Impact of sars-cov-2 interventions on dengue transmission. PLoS Negl Trop Dis. 2020;14(10). 50. 50.Jt L, Bl D, j A, Vj L, Ar C. Decreased dengue transmission in migrant worker populations in Singapore attributable to SARS-CoV-2 quarantine measures. J Travel Med. 2021;28(2). 51. 51.Sudria M, Andreatta M, Defagó M. Los efectos de la cuarentena por coronavirus (Covid-19) en los hábitos alimentarios en Argentina. Asociación Argentina de Dietistas y Nutricionistas Dietistas; Diaeta. 2020;38(171):10–19. 52. 52.Rojas C. Cuarentena, aislamiento forzado y uso de drogas. Panamerican journal of neuropsychology. 2020;14(1). 53. 53.Ribot Reyes V, Chang Paredes N, Gonzáles Castillo A. Efectos de la COVID-19 en la salud mental de la población. Revista Habanera de Ciencias médicas. 2020;19. 54. 54.Gómez Tejeda J, Besteiro Arjona E, Hernández Pérez C, Góngora Villares Y. Impacto psicológico causado por la pandemia de COVID-19. Revista Científica Estudiantil de Cienfuegos INMEDSUR. 2020;3(2). 55. 55.Chiatchoua C, Lozano C, Macías-Durán J. Análisis de los efectos del COVID-19 en la economía mexicana. Revista Del Centro De Investigación De La Universidad La Salle. 2020;14(53):265–290. [1]: /embed/inline-graphic-1.gif [2]: /embed/inline-graphic-2.gif [3]: /embed/inline-graphic-3.gif [4]: /embed/graphic-1.gif [5]: /embed/graphic-3.gif [6]: /embed/graphic-4.gif [7]: /embed/graphic-5.gif [8]: /embed/graphic-6.gif [9]: /embed/graphic-7.gif [10]: /embed/inline-graphic-4.gif [11]: /embed/inline-graphic-5.gif [12]: /embed/inline-graphic-6.gif [13]: /embed/inline-graphic-7.gif [14]: /embed/graphic-10.gif [15]: /embed/graphic-11.gif [16]: /embed/graphic-12.gif [17]: /embed/inline-graphic-8.gif [18]: /embed/graphic-13.gif [19]: /embed/graphic-14.gif