Next Article in Journal
Fock Quantization of a Klein–Gordon Field in the Interior Geometry of a Nonrotating Black Hole
Next Article in Special Issue
Optimal Control Theory and Calculus of Variations in Mathematical Models of Chemotherapy of Malignant Tumors
Previous Article in Journal
Mixed-Type Hypergeometric Bernoulli–Gegenbauer Polynomials
Previous Article in Special Issue
Necessary Conditions for the Optimality and Sustainability of Solutions in Infinite-Horizon Optimal Control Problems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Control of a Two-Patch Dengue Epidemic under Limited Resources

by
Edwin Barrios-Rivera
1,*,
Olga Vasilieva
1 and
Mikhail Svinin
2
1
Department of Mathematics, Universidad del Valle, Cali 760032, Colombia
2
Department of Information Science and Engineering, Ritsumeikan University, Kusatsu 525-8577, Japan
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(18), 3921; https://doi.org/10.3390/math11183921
Submission received: 14 August 2023 / Revised: 9 September 2023 / Accepted: 11 September 2023 / Published: 15 September 2023

Abstract

:
Despite the ongoing preventive measures of vector control, dengue fever still presents outbreaks, and daily commuting of people also facilitates its propagation. To contain the disease spread after an outbreak has already occurred, the local healthcare authorities are compelled to perform insecticide spraying as a corrective measure of vector control, thus trying to avoid massive human infections. Several issues concerned with the practical implementation of such corrective measures can be solved from a mathematical standpoint, and the purpose of this study is to contribute to this strand of research. Using as a basis a two-patch dengue transmission ODE model, we designed the patch-dependent optimal strategies for the insecticide spraying with the optimal control approach. We also analyzed the response of the optimal strategies to three alternative modes of budget cuts under different intensities of daily commuting. Our approach illustrated that trying “to save money” by reducing the budget for corrective control is completely unwise, and the anticipated “savings” will actually turn into considerable additional public spending for treating human infections, which could have been averted by a timely corrective intervention.

1. Introduction

Dengue is a viral infection caused by the dengue virus (DENV1-DENV4), which is primarily transmitted to human individuals through the bite of infected Aedes aegypti and Aedes albopictus mosquitoes. Since only females need to ingest human blood to mature their eggs and male mosquitoes do not bite people, only the Aedes females are considered vectors of dengue and other mosquito-borne infections.
Tropical countries featuring the presence and abundance of Aedes mosquitoes usually carry the burden of dengue and other mosquito-borne viral infections in the human population. The local healthcare authorities conduct preventive control actions, such as eliminating mosquito breeding sites, stagnant water and rainwater treatment, environmental management, and educational campaigns targeting local communities, regardless of the presence or absence of disease outbreaks. However, these preventive measures cannot eliminate virus circulation, and disease outbreaks still occur.
To reduce the spread of dengue fever after an outbreak has already occurred, the local healthcare authorities need to implement the so-called corrective measure of vector control in order to curb the epidemics and contain the disease’s spread. Among the corrective actions, the most common is spreading insecticides in urban areas, excluding the interior of the houses. Notably, the costs of such actions are not usually anticipated in the yearly budget of local healthcare entities, and these actions are typically dependent on the available stock of the insecticide and human power needed to implement the spraying.
Among numerous questions regarding the practical implementation of corrective insecticide-based vector control measures, several issues can be solved from a mathematical standpoint, and the purpose of this study is to contribute to this strand of research. To make our study more realistic, we introduce a corrective control action modeling of the insecticide spraying into a two-patch model of dengue transmission that mimics the 2013 dengue outbreak in the city of Cali, Colombia, and its suburban areas [1]. This model has been thoroughly analyzed in the mentioned work, and its parameters have been estimated using the data gathered by the Municipal Secretariat of Public Health of Cali, Colombia.
The two-patch model in [1] was a good starting point because it considers daily commuting from suburban areas to Cali, using the Lagrange approach and residence times [2]. It is worth noting that the impact of daily commuting on the disease transmission was captured in other works (see, e.g., [3,4,5,6,7]); however, the design of patch-dependent strategies for corrective interventions based on insecticide spraying has not yet been addressed from the standpoint of dynamic optimization. In the present study, we derive the patch-dependent strategies for insecticide spraying using the optimal control framework, interpret their practical implementation, and estimate their underlying costs. We also show the resilience and robustness of the proposed strategies with respect to the intensity of people commuting. This result bears a remarkable value for the practical implementation of corrective control measures. Namely, healthcare entities may design the patch-dependent strategies for insecticide spraying basing only on the current state of the disease in both patches, and without estimating the actual number of daily commuters.
Following the approach developed in [8] for a one-patch model, we adjust the patch-dependent strategies for insecticide spraying to the budget restrictions imposed by healthcare entities. Generally speaking, the budget reductions can be made either in both patches or in only one of them. Here, it is pertinent to point out that the healthcare authorities usually prioritize corrective vector control in the big cities, thus neglecting the suburban areas. However, as our results show, those common practices are less cost-effective than their alternatives. In fact, they lead to a significant increase in the total public expenditures, i.e., the expenses for insecticide spraying and disease-related costs in the human population.
Lastly, we would like to emphasize that the present study is mainly focused on applying the optimal control methodology in the practical context of corrective vector control rather than exhibiting any new findings in the optimal control theory. From this standpoint, our results are expected to contribute to the practical aspects of the corrective disease control, whether featuring resource limitations or not, and the latter constitutes the primary novelty of our study. If the dengue outbreak has already begun, insecticide-based mosquito control is the most accessible and cost-effective measure to reduce the number of human infections and prevent further spread of the disease. In this context, our approach clearly illustrates that trying “to save money” at the beginning of dengue epidemics by reducing the budget for corrective control is entirely unwise and that such anticipated “savings” will turn into considerable additional public spending for treating human infections that could have been averted by a timely insecticide spraying at a much lower cost.
The paper is organized as follows. Section 2 provides the comprehensive description of the two-patch dengue transmission model originally introduced in [1]. In Section 3, we introduce two exogenous variables mimicking the insecticide spraying in both patches, formulate the optimal control problem of dynamic optimization, and justify the existence of its solution.
The main results are presented in Section 4, which consists of four subsections. Section 4.1 provides realistic data referring to the optimal control model that further allows us to estimate the total public expenditures related to the costs of corrective control strategies and the disease-related outlays. Section 4.2 provides numerical solutions of the optimal control problem with different commuting intensities when no budget limits are imposed. From these solutions, we estimate the costs necessary for implementing the optimal insecticide spraying in both patches. The situation when only a share of resources needed for implementing the corrective control strategies is available is thoroughly analyzed in Section 4.3, where we review three alternative modes of budget cuts combined with different intensities of daily commuting. The impact of the budget cut distributions between two patches on the disease propagation is examined in Section 4.4, and the modes of budget cuts corresponding to public expenditures combined with different intensities of daily commuting are estimated. From the latter, we conclude that the uniform budget reduction for both patches is the best way to act in resource-limited settings. Finally, Section 5 summarizes the key outcomes of the present work.

2. Two-Patch Dengue Transmission Model

Let us begin by revisiting the formulation of the two-patch dengue transmission model of the SEIR(S)-SEI type originally presented in [1]. This model was derived based on the following standard assumptions:
  • The populations of human residents in both patches, N h 1 and N h 2 , as well as the populations of vectors inhabiting both patches, N v 1 and N v 2 , remain essentially invariant in time.
  • All populations N h i , N v i , i = 1 , 2 are homogeneous regarding susceptibility, exposure, and attraction.
  • There is no dengue-induced mortality, neither for human individuals nor for vectors.
  • Once DENV-infected, the mosquitoes never recover and stay infectious during their lifetime.
  • There is no vertical DENV transmission, so all the newborns (humans and vectors) are DENV-free.
  • Human individuals can move between the patches, while mosquitoes remain in the same patch during their lifetime.
The total populations of human individuals N h i , i = 1 , 2 residing in each patch are subdivided into four disjoint classes: susceptible (DENV-free people), exposed (DENV-carrying people who are not infectious yet), infectious (DENV-carrying people who are capable of transmitting the virus to female mosquitoes during the bite), and individuals recovered from the disease. The fractions of the corresponding human groups regarding the total populations N h i , i = 1 , 2 residing in each patch are denoted with S h i ,   E h i ,   I h i , and R h i , respectively. In turn, the total populations of vectors N v i , i = 1 , 2 inhabiting each patch are subdivided into three disjoint classes: susceptible (DENV-free female mosquitoes), exposed (DENV-carrying female mosquitoes which are not infectious yet), and infectious (DENV-carrying female mosquitoes which are capable of transmitting the virus to human individuals during the bite). The fractions of the corresponding groups of vectors regarding the total populations N v i , i = 1 , 2 inhabiting each patch are denoted with S v i ,   E v i , and I v i , respectively. Notably, all the population fractions or proportions introduced above vary in time and satisfy the following relationships valid for all t 0 :
S h i + E h i + I h i + R h i = 1 , S v i + E v i + I v i = 1 , i = 1 , 2 .
Numerous studies have confirmed that dengue disease caused by any of the four dengue serotypes DENV1-DENV4 might result either in asymptomatic infection or range from mild febrile illness with flu-like symptoms to severe dengue or dengue shock syndrome (DSS). There are no exact estimations of the proportion of asymptomatic infections among all dengue morbidity, but they may be even higher in density than the symptomatic cases [9,10,11]. Therefore, we may assume that some infectious human individuals from the classes I h i , i = 1 , 2 may exhibit no symptoms and be unaware of being contagious. Such people may commute between the two patches and contribute to the disease spreading.
To model the demographic changes while keeping the populations of humans and vectors essentially invariant over time, we assume that both human populations N h i , i = 1 , 2 and both vector populations N v i , i = 1 , 2 have equal birth and death rates. Namely, the human birth and death rate is denoted by μ h , while the birth and death rate for mosquitoes is denoted by μ v . For convenience, Table 1 summarizes the descriptions of all baseline parameters related to disease transmission dynamics.
The exposed or latent human and vector individuals from the classes E h i and E v i , i = 1 , 2 become infectious after the virus incubation period and then pass to the classes I h i and I v i , i = 1 , 2 , respectively. We suppose the intrinsic and extrinsic virus incubation (in humans and mosquitoes) lasts for 1 / κ h and 1 / κ v days, respectively. Thus, the conversion rates from latency to infectiousness are denoted by κ h and κ v for humans and vectors, respectively.
We also assume that the recovery from the disease takes 1 / γ days, and the recovered individuals become susceptible to heterologous DENV serotypes after 1 / τ days if various DENV serotypes circulate simultaneously in the environment. Notably, no reinfection is expected under the strong dominance of only one homologous DENV serotype, and in such a case, it is supposed that τ = 0 .
It is worthwhile to point out that the susceptible vectors inhabiting Patch 1 (resp. Patch 2) may acquire the virus after biting either an infectious human resident from class I h 1 (resp. class I h 2 ) or an asymptomatic visitor to Patch 1 (resp. Patch 2) belonging to the class I h 2 (resp. class I h 1 ). Similarly, the susceptible human individuals residing in Patch 1 (resp. Patch 2) may acquire the virus after being bitten by either an infectious vector from their home patch, belonging to the class I v 1 (resp. I v 2 ), or by an infectious vector from the other patch I v 2 (resp. I v 1 ) while visiting Patch 2 (resp. Patch 1). Thus, it becomes clear that people commuting between the two patches may affect the virus transmission, even though all vectors remain in their home patches.
The commuting of people between two patches is modeled using the so-called Lagrangian approach [2,3], under which the movement of individuals between patches is mimicked in terms of their residence times, i.e., the shares of time spent at each patch. Let
Q = q 11 q 12 q 21 q 22
denote the residence time matrix [3] that couples both patches and whose elements q i j express the shares or fractions of time a person residing in Patch i spends, on average, in Patch j. Notably, the elements of Q satisfy the following relationships:
0 q i j 1 , j = 1 2 q i j = 1 , i , j = 1 , 2 .
Using the elements of (2) and following the idea exposed in [2], we introduce another matrix
P = P 11 P 12 P 21 P 22 = q 11 N h 1 q 11 N h 1 + q 21 N h 2 q 12 N h 1 q 12 N h 1 + q 22 N h 2 q 21 N h 2 q 11 N h 1 + q 21 N h 2 q 22 N h 2 q 12 N h 1 + q 22 N h 2 ,
which is referred to as the matrix of effective residence times, since its components are defined in terms of the elements of the residence times matrix Q (see (2)), weighted by the respective shares or proportions of the effective population sizes at each patch. Notably, the elements of P satisfy the following relationships:
0 P i j 1 , i = 1 2 P i j = 1 , i , j = 1 , 2
and they express the proportions of the residents from Patch i effectively present in Patch j. For example, P 12 is the fraction of residents of Patch 1 who visit Patch 2 (that is, q 12 N h 1 ) relative to the effective population size of Patch 2 (which is equal to q 12 N h 1 + q 22 N h 2 ).
To model the disease incidence, we recall from other works (see, e.g., [1,2,3]) that the per capita “human-to-vector” virus transmission rate β v is usually assumed to be proportional to the mosquito biting rate α and the probability p v [ 0 , 1 ] of the virus transmission, that is, β v = α p v . Here, the mosquito biting rate refers to the number of times a female mosquito needs to have a blood meal by biting people on average daily. As the parameters α and p v are unaffected by the people’s movements, the same “human-to-vector” virus transmission rate β v can be assumed for both patches.
Usually, for a female mosquito, it suffices to bite one person for having her blood meal. On the other hand, the same person may receive bites from more than one female mosquito. Therefore, the “vector-to-human” virus transmission rate β h depends not only on the mosquito biting rate α and the probability p h [ 0 , 1 ] of the virus transmission but also on the so-called vectorial density, i.e., the average number of female mosquitoes per one human individual. Even though all vectors stay in their patches and both vector populations N v i remain essentially invariant, the number of human individuals effectively present at each patch varies when people are commuting between the patches. Therefore, the “vector-to-human” virus transmission rate β h depends on the so-called effective vectorial density at each patch [1]. The latter can be defined using the components of the residence time matrix Q in the following way:
Patch 1 : N v 1 q 11 N h 1 + q 21 N h 2 , Patch 2 : N v 2 q 12 N h 1 + q 22 N h 2
Thus, the effective vectorial density expresses an average number of female mosquitoes in Patch j per one human individual effectively present in Patch j for j = 1 , 2 . Using the quantities defined by (4), the patch-dependent “vector-to-human” contact rates β h i , i = 1 , 2 can be written as follows:
β h 1 = α p h N v 1 q 11 N h 1 + q 21 N h 2 , β h 2 = α p h N v 2 q 12 N h 1 + q 22 N h 2 .
Finally, the normalized two-patch dengue transmission model can be written as a dynamical system of fourteen ordinary differential equations:
S ˙ v 1 ( t ) = μ v β v S v 1 ( t ) j = 1 2 P j 1 I h j ( t ) μ v S v 1 ( t ) E ˙ v 1 ( t ) = β v S v 1 ( t ) j = 1 2 P j 1 I h j ( t ) ( μ v + κ v ) E v 1 ( t ) I ˙ v 1 ( t ) = κ v E v 1 ( t ) μ v I v 1 ( t ) S ˙ h 1 ( t ) = μ h S h 1 ( t ) j = 1 2 β h j q 1 j I v j ( t ) μ h S h 1 ( t ) + τ R h 1 ( t ) E ˙ h 1 ( t ) = S h 1 ( t ) j = 1 2 β h j q 1 j I v j ( t ) ( μ h + κ h ) E h 1 ( t ) I ˙ h 1 ( t ) = κ h E h 1 ( t ) ( γ + μ h ) I h 1 ( t ) R ˙ h 1 ( t ) = γ I h 1 ( t ) ( μ h + τ ) R h 1 ( t ) S ˙ v 2 ( t ) = μ v β v S v 2 ( t ) j = 1 2 P j 2 I h j ( t ) μ v S v 2 ( t ) E ˙ v 2 ( t ) = β v S v 2 ( t ) j = 1 2 P j 2 I h j ( t ) ( μ v + κ v ) E v 2 ( t ) I ˙ v 2 ( t ) = κ v E v 2 ( t ) μ v I v 2 ( t ) S ˙ h 2 ( t ) = μ h S h 2 ( t ) j = 1 2 β h j q 2 j I v j ( t ) μ h S h 2 ( t ) + τ R h 2 ( t ) E ˙ h 2 ( t ) = S h 2 ( t ) j = 1 2 β h j q 2 j I v j ( t ) ( μ h + κ h ) E h 2 ( t ) I ˙ h 2 ( t ) = κ h E h 2 ( t ) ( γ + μ h ) I h 2 ( t ) R ˙ h 2 ( t ) = γ I h 2 ( t ) ( μ h + τ ) R h 2 ( t )
where dots in the left-hand sides denote the derivative with respect to time, that is,
X ˙ ( t ) : = d X d t .
The above-formulated two-patch dengue transmission model was thoroughly studied in [1] under different human mobility scenarios and without vector control measures. In the sequel, this model will serve as a basis for testing different strategies of vector control intervention applied in both patches. The following section encompasses the formulation of dengue control problem from the standpoint of optimal control.

3. Optimal Control Approach

As millions of people worldwide are exposed to the risk of vector-borne infections, dengue control is a collective effort involving individuals, communities, healthcare systems, and government authorities. With several types of vaccines against all DENV serotypes still under development [12], the disease control efforts continue targeting the reduction in the local mosquito density. Vector control actions conducted by local healthcare authorities in dengue-endemic areas can be divided into two majors groups:
  • Preventive actions include eliminating mosquito breeding sites, stagnant water and rainwater treatment, environmental management, and educational campaigns targeting local communities. These are routine interventions carried out regardless of the presence or absence of disease outbreaks.
  • Corrective actions, such as insecticide spraying, are conducted as part of an emergency response to curb the disease outbreak or a cluster of identified cases in a particular area.
In the present study, we will address the second group of control actions from the standpoint of dynamic optimization and accounting for people commuting between two patches. Taking as an example the epidemiological situation of dengue in the city of Cali and its suburbs (Republic of Colombia), we will try to model the current practices of local healthcare authorities for implementing corrective control actions seeking eventual suppression of dengue cases in the city. As the corrective measures do not always have a sufficient budget, we will also address the flexibility and robustness of corrective actions under budget constraints.
To model the corrective control actions of insecticide spraying applied in both patches, we introduce two exogenous variables u i ( t ) [ 0 , 1 ] ,   i = 1 , 2 ,   t [ 0 , T ] where T > 0 denotes the time horizon. These variables u 1 ( t ) and u 2 ( t ) will directly affect the mosquito mortality rate, μ v , in both patches. However, as the model presented in Section 2 is normalized, the variables u i ( t ) ,   i = 1 , 2 should also be scaled accordingly while bearing a meaningful relation to the amount of insecticide sprayed.
To derive a meaningful definition of u i ( t ) ,   i = 1 , 2 , it is worthwhile to recall that there are generally accepted norms and guidelines for daily insecticide spraying regarding the volume of insecticide used. Specific regulations may vary based on the type of insecticide, target mosquito species, and local environmental policies. One of the key indicators is the so-called insecticide application rate. It refers to the amount of a diluted insecticide (volume) applied daily per unit area. Let v i ( t ) [ 0 , D ] ,   i = 1 , 2 represent the total volume of a diluted insecticide destined for application in Patch i on day t, and let D > 0 denote the maximal volume of insecticide that is allowed for daily application over a unit area at each patch. Thus, the exogenous variables u i ( t ) ,   i = 1 , 2 can be viewed as
u i ( t ) = v i ( t ) area of Patch i × D , i = 1 , 2 .
Under the above definition, u 1 ( t ) and u 2 ( t ) become dimensionless and acquire meaningful interpretations. For instance, u 1 ( t ) = 1 ,   u 2 ( t ) = 0.5 implies that, on day t, the maximal allowed volume of a diluted insecticide should be uniformly sprayed in Patch 1, while only 50 % of the maximal allowed volume should be uniformly spread in Patch 2.
It is worthwhile to point out that the maximum volume of insecticide allowed for daily application over a unit area (denoted by D) may vary depending on the insecticide dilution. The concentration of the insecticide also defines the insecticide’s level of toxicity or lethality for target mosquito species and usually ranges between 20 % and 90 % [13]. In the sequel, the insecticide lethality will be denoted by η ( 0 , 1 ) , while it is also understood that smaller values of η will require a larger volume v i ( t ) ,   i = 1 , 2 of diluted insecticide and the constant D may take larger values (cf. (6)). In accordance with definition (6), we can determine the so-called set of admissible controls as U : = [ 0 , 1 ] × [ 0 , 1 ] .
As insecticide spraying is a corrective measure of vector control, healthcare authorities may consider applying this type of intervention when there is a significant increase in dengue cases among the human population. Thus, the ultimate purpose of insecticide spraying consists in containing dengue epidemics. Still, the effect of this control action is not immediate and can be observed within 15–30 days after the insecticide spraying begins. It is worth pointing out that the duration and frequency of insecticide-spraying campaigns depend on the environmental norms (expressed by the quantity D) and other factors, such as the total stock of insecticides, available personnel, equipment, etc. Hence, this type of intervention is usually limited in terms of costs and time and may continue only for several days at most.
Therefore, the global challenge for the healthcare authorities consists in choosing an optimal strategy for insecticide spraying in both patches, u 1 * ( t ) , u 2 * ( t ) U ,   t [ 0 , T ] , that minimizes the total costs associated with human infections during the dengue outbreak together with the overall costs of control interventions while meeting the budget constraint whenever such a constraint is imposed. From the mathematical standpoint, this goal can be expressed by the following objective functional:
J u 1 ( t ) , u 2 ( t ) = 0 T W 1 N h 1 I h 1 ( t ) + W 2 N h 2 I h 2 ( t ) + W 3 A 1 u 1 ( t ) + W 4 A 2 u 2 ( t ) d t ,
where N h i , i = 1 , 2 refer to the total human populations residing in both patches and A i , i = 1 , 2 denote their areas. Here, the weight coefficients W 1 , W 2 express the medical and societal cost related to one human infection, while W 3 , W 4 stand for the cost of insecticide spraying per unit of area.
Thus, the dynamic optimization problem can be formulated as minimizing the objective functional (7), that is,
min ( u 1 , u 2 ) U J u 1 ( t ) , u 2 ( t )
subject to
S ˙ v i ( t ) = μ v β v S v i ( t ) j = 1 2 P j i I h j ( t ) μ v + η u i ( t ) S v i ( t )
E ˙ v i ( t ) = β v S v i ( t ) j = 1 2 P j i I h j ( t ) μ v + κ v + η u i ( t ) E v i ( t )
I ˙ v i ( t ) = κ v E v i ( t ) μ v + η u i ( t ) I v i ( t )
S ˙ h i ( t ) = μ h S h i ( t ) j = 1 2 β h j q i j I v j ( t ) μ h S h i ( t ) + τ R h i ( t )
E ˙ h i ( t ) = S h i ( t ) j = 1 2 β h j q i j I v j ( t ) ( μ h + κ h ) E h i ( t )
I ˙ h i ( t ) = κ h E h i ( t ) ( γ + μ h ) I h i ( t )
R ˙ h i ( t ) = γ I h i ( t ) ( μ h + τ ) R h i ( t )
S v i ( 0 ) , E v i ( 0 ) , I v i ( 0 ) , S h i ( 0 ) , E h i ( 0 ) , I h i ( 0 ) , R h i ( 0 ) ( 1 , 0 , 0 , 1 , 0 , 0 , 0 ) , i = 1 , 2
and ( u 1 , u 2 ) U implies that
0 u i ( t ) 1 , t [ 0 , T ] , i = 1 , 2 .
It is important to highlight that the initial conditions (9h) are nonnegative, have constant values between 0 and 1, and satisfy the following relationships:
S h i ( 0 ) + E h i ( 0 ) + I h i ( 0 ) + R h i ( 0 ) = 1 , S v i ( 0 ) + E v i ( 0 ) + I v i ( 0 ) = 1 , i = 1 , 2 .
Numerical values given in (9h) correspond to the initial proportions of the vector and human populations in both patches before the control intervention begins. Furthermore, the control functions (10) are assumed to be piecewise continuous.
The necessary condition for the solvability of the optimal control problem (8)–(10) is the solvability of the dynamical system (9) with any admissible controls (10). The following statement gives the validity of this condition.
Proposition 1. 
For each pair of piecewise continuous real functions given by (10), the initial value problem (9) has a unique solution defined for all t 0 and T > 0 , while all the components of this solution remain within the range [ 0 , 1 ] .
Proof. 
Let X = X 1 , X 2 R + 14 denote the vector of the state variables of the dynamical system (9), where X i = S v i , E v i , I v i , S h i , E h i , I h i , R h i ,   i = 1 , 2 . The controlled system (9) can then be written in the form
X ˙ ( t ) = F X ( t ) , u 1 ( t ) , u 2 ( t ) ,
where F X ( t ) , u 1 ( t ) , u 2 ( t ) is the vector field corresponding to the right-hand sides of Equation (9). Let X ( t ; u 1 , u 2 ) denote the solution to (9) corresponding to ( u 1 , u 2 ) U .
In the absence of control intervention, that is, when u 1 = u 2 = 0 , it was shown in [1] that the dynamical system (9) has a unique nonnegative solution X ( t ; 0 , 0 ) engendered by any initial conditions complying with (11). Furthermore, this solution continuously depends on the initial conditions, can be extended for all t 0 , and its components satisfy the relationships (1) while staying within the range [ 0 , 1 ] .
Moreover, as F X ( t ) , u 1 ( t ) , u 2 ( t ) remains continuous and locally Lipschitz in X for each pair ( u 1 , u 2 ) U , a unique solution to (9) exists for all u i ( t ) , i = 1 , 2 defined by (10) for any T > 0 . More details regarding this standard result can be consulted in [14,15] or similar textbooks. Notably, the right-hand sides of the dynamical system (9) are non-increasing with respect to the exogenous variables u i , i = 1 , 2 . Therefore, the growth rate F X ( t ) , u 1 ( t ) , u 2 ( t ) with u i ( 0 , 1 ] , i = 1 , 2 is smaller than F X ( t ) , 0 , 0 . Thus, all the components of the solution to (9) with ( u 1 , u 2 ) U are also bounded and remain within the range [ 0 , 1 ] .  □
Once the solvability of the controlled dynamical system (9) is established for all admissible controls ( u 1 , u 2 ) U , we proceed to formulate the following statement for ensuring the solvability of the optimal control problem (8)–(10).
Proposition 2. 
There exists a pair of optimal controls u 1 * ( t ) , u 2 * ( t ) satisfying (10), and the corresponding solutions S v i * ( t ) , E v i * ( t ) , I v i * ( t ) , S h i * ( t ) , E h i * ( t ) , I h i * ( t ) , R h i * ( t ) ,   i = 1 , 2 of the initial value problem (9), which minimize the functional (7).
Proof. 
To formally prove the existence of the optimal controls, we make use of the standard existence result given [16] and verify that all the hypotheses established are satisfied for the optimal control problem (8)–(10). In this context, we point out the validity of the following statements:
(i)
In view of Proposition 1, the state trajectories X ( t ; u 1 , u 2 ) of the controlled dynamical system (9) are bounded for all admissible controls u 1 , u 2 U and for all t 0 .
(ii)
The set U = [ 0 , 1 ] × [ 0 , 1 ] of admissible controls is compact and convex in R + 2 .
(iii)
The right-hand sides of the dynamical system (9) are linear with respect to u 1 and u 2 .
(iv)
The set of all possible initial states X ( 0 ) and the set of all possible terminal states X ( T ) , 0 < T < are both compact because each state variable remains within the range [ 0 , 1 ] for all t 0 .
(v)
As the integrand of the objective functional (7) is linear with respect to u 1 and u 2 , it is also convex in u 1 ,   u 2 and is bounded from below.
Conditions (i)–(v) indicate that the sufficient conditions for the existence of optimal controls formulated in [16] (see Theorem 4.1 and Corollary 4.1, pp. 68–69) are satisfied. This ensures the existence of u 1 * : [ 0 , T ] [ 0 , 1 ] and u 2 * : [ 0 , T ] [ 0 , 1 ] that minimize functional (7) over any finite time interval t [ 0 , T ] , 0 < T < .  □
The optimal control problem (8)–(10) has fourteen state variables plus two control variables. Due to its high dimensionality, there are no prospects of obtaining an exact solution in the closed form; therefore, the optimal control problem (8)–(10) can only be solved numerically. Before we proceed to present its numerical solutions (Section 4), let us expose some arguments concerning the form of the optimal strategies u 1 * ( t ) , u 2 * ( t ) and their interpretations.
Knowing from Proposition 2 that the optimal solution exists and observing that the right-hand sides of the dynamical system (9), as well as the integrand of the objective functional (7), are linear with respect to u 1 , u 2 , one may expect an optimal solution in the bang-bang form. Notably, the bang-bang form of optimal controls fully agrees with the common practices of insecticide-spreading campaigns in many countries affected by dengue and other mosquito-borne diseases [8].
By formally applying the solution procedure based on the Pontryagin minimum principle (see, e.g., [17] or similar textbooks), we define the Hamiltonian corresponding to the optimal control problem (8)–(10) that has the following form:
H ( Λ , X , u 1 , u 2 ) = W 1 N h 1 I h 1 + W 2 N h 2 I h 2 + W 3 A 1 u 1 + W 4 A 2 u 2 + Λ ( t ) , F X , u 1 , u 2 ,
where · , · stands for scalar product and Λ = Λ 1 , Λ 2 R 14 ,   Λ i = ( λ 1 i , , λ 7 i ) ,   i = 1 , 2 denotes the vector of adjoint or co-state variables corresponding to the state variables X = X 1 , X 2 ,   X i = S v i , E v i , I v i , S h i , E h i , I h i , R h i ,   i = 1 , 2 . Here, the co-state vector function Λ ( t ) is absolutely continuous [17] and satisfies the adjoint ODE system with transversality condition
Λ ˙ ( t ) = H X , Λ ( T ) = 0 .
The detailed form of the adjoint system (13) is given in Appendix A.
In view of Proposition 2, the optimal control problem (8)–(10) has an optimal solution, say u 1 * ( t ) , u 2 * ( t ) U , which is defined for all t [ 0 , T ] ,   0 < T < . This solution, together with the corresponding optimal state X * ( t ) = X t ; u 1 * ( t ) , u 2 * ( t ) and the underlying co-state Λ ( t ) satisfying the adjoint ODE system (13), must comply with the necessary condition of optimality in the form of the Pontryagin minimum principle. In other words, the Hamiltonian (12) must satisfy the following condition almost at each t [ 0 , T ] along the optimal path X * ( t ) :
H Λ , X * , u 1 * , u 2 * H Λ , X * , u 1 , u 2 , ( u 1 , u 2 ) U .
As the Hamiltonian is linear with respect to u 1 and u 2 , we have that its partial derivatives
H u i = W 2 + i A i η λ 1 i S v i + λ 2 i E v i + λ 3 i I v i , i = 1 , 2
are independent of u 1 and u 2 and, generally speaking, the Hamiltonian attains its minimum with respect to u i ,   i = 1 , 2 at the extreme points of the interval [ 0 , 1 ] . Then we can derive that
u i * ( t ) = 0 if φ i ( t ) > 0 1 if φ i ( t ) < 0 i = 1 , 2
where
φ i ( t ) : = W 2 + i A i η λ 1 i ( t ) S v i * ( t ) + λ 2 i ( t ) E v i * ( t ) + λ 3 i ( t ) I v i * ( t ) , i = 1 , 2 ,
denote the switching functions which represent the right-hand sides of the relationships in (14) evaluated at each t [ 0 , T ] . Thus, the optimal control u i * ( t ) ,   i = 1 , 2 is piecewise constant and switches between the “on” and “off” modes expressed by u i * ( t ) = 1 and u i * ( t ) = 0 , respectively. The number of switches equals the number of roots of the switching function φ i ( t ) . However, if there exists an interval [ t 1 , t 2 ] [ 0 , T ] ,   t 1 < t 2 within which it holds that φ i ( t ) 0 ,   t [ t 1 , t 2 ] , there is a possibility that the optimal control is interior, that is, 0 < u i * ( t ) < 1 ,   t [ t 1 , t 2 ] .
Optimal controls of the type (15), referred to as “bang-bang” control and defined by the underlying switching functions (16), have a straightforward interpretation linking the optimal strategy’s costs and benefits. Indeed, in the switching functions defined by (16), the positive term W 2 + i A i refers to the cost of the strategy u i * ( t ) at each t [ 0 , T ] , while the negative term of the form η stands for the benefit rendered by this strategy. If the cost exceeds the benefit for some t ^ [ 0 , T ] (that is, φ i ( t ^ ) > 0 ), it is optimal to assign u i * ( t ^ ) = 0 and take no action. Alternatively, if the cost drops below the benefit for some t ^ [ 0 , T ] (that is, φ i ( t ^ ) < 0 ), it is optimal to assign u i * ( t ^ ) = 1 and use all available resources.
Thus, if φ ( t ^ ) = 0 for some t ^ [ 0 , T ] , it implies that the cost equals the benefit and a change in the strategy is expected immediately after t ^ , that is, at t = t ^ + Δ t . If the switching function φ i ( t ) merely changes its sign at t ^ , meaning that φ i ( t ^ ) = 0 while φ i ( t ± Δ t ) 0 with Δ t arbitrarily small, the change in the strategy u i * ( t ) at t = t ^ is a jump from 1 to 0 or vice versa. However, if the condition φ ( t ^ + Δ t ) = 0 holds for all t [ t ^ , t ^ + Δ t ] ,   Δ t > 0 then it is said that u i ( t ) has a singular arc on the segment [ t ^ , t ^ + Δ t ] so that u i ( t ) ( 0 , 1 ) when t [ t ^ , t ^ + Δ t ] . Given the high dimensionality of our model, further analysis of the possible appearance of singular arcs seems rather hard to perform. Still, more details regarding the methodology for singular arc detection in lower-dimensional models can be consulted in [18] and similar textbooks.
The foregoing rationale will help us interpret the optimal solutions found numerically in the next section.

4. Numerical Results and Discussion

In this section, we will present the numerical results of the optimal control problem (8) and (9) obtained using the GPOPS-II software package (note that GPOPS-II Manual with basic descriptions can be downloaded from http://www.gpops2.com/ (accessed on 12 August 2023)) developed by PR Optimization Research LLC for the MATLAB platform [19,20]. This toolbox can solve different kinds of optimal control problems and implements the variable-order Gaussian quadrature collocation methods numerically. The GPOPS-II solver does not use numerical integration and employs the Legendre–Gauss–Radau quadrature orthogonal collocation technique to convert the continuous-time optimal control problem into a large sparse nonlinear programming problem [21].

4.1. Parameter Values and Initial Data

To make our numerical results more meaningful, we use as an example the propagation of dengue in the city of Cali, Colombia (Patch 1) and its nearby localities (Patch 2) because this geographic area is considered hyper-endemic with regards to dengue fever [22], and the values of the key parameters of the model (9) have been estimated in [1] using the data of the 2013 dengue epidemic outbreak. The values of these parameters are given in Table 1.
The total sizes of the human and vector populations, N h i and N v i , i = 1 , 2 , will also be assumed to be the same as in [1], that is, N h 1 = 2,319,655 (Patch 1, the city of Cali) and N h 2 = 527,091 (Patch 2, the suburbs). Here, we take into account that Cali has three major satellite towns, Palmira (57%), Jamundi (23%), and Yumbo (20%), many of whose residents commute daily to Cali; however, there are no reliable data regarding the estimated share of such commuters. Using the estimations of the corresponding mosquito populations obtained in [1], we also have N v 1 = 4.36 × 10 6 and N v 2 = 1.49 × 10 6 together with the initial conditions (9h) that correspond to the beginning of the 2013 dengue outbreak:
S v 1 ( 0 ) = 1 E v 1 ( 0 ) I v 1 ( 0 ) , E v 1 ( 0 ) = 3.2721 × 10 4 , I v 1 ( 0 ) = 4.2936 × 10 4 , S v 2 ( 0 ) = 1 E v 2 ( 0 ) I v 2 ( 0 ) , E v 2 ( 0 ) = 3.0832 × 10 5 , I v 2 ( 0 ) = 4.0458 × 10 5 ,
E h 1 ( 0 ) = 1.0907 × 10 4 , I h 1 ( 0 ) = 1.4312 × 10 4 , R h 1 = 1 × 10 4 , E h 2 ( 0 ) = 1.0907 × 10 4 , I h 2 ( 0 ) = 1.4312 × 10 4 , R h 2 = 1 × 10 5 ,
S h 1 ( 0 ) = 1 E h 1 ( 0 ) I h 1 ( 0 ) R h 1 ( 0 ) , S h 2 ( 0 ) = 1 E h 2 ( 0 ) I h 2 ( 0 ) R h 2 ( 0 ) .
Note that the above initial conditions are already normalized to the population sizes N h i , N v i , i = 1 , 2 . They reflect the beginning of the dengue outbreak in Cali (Patch 1), where the fraction of recovered people is larger than in the suburbs (Patch 2). At the same time, the average mosquito density (i.e., the mean number of vectors per one human resident) is assumed to be lower in the city than in the suburbs. This assumption stays in line with other works (see, e.g., [1,7,23]) and captures four important factors:
  • Compared to the metropolitan areas, the suburbs have wider green zones that harbor mosquitoes and provide breeding sites close to the houses.
  • There are many high-story apartment blocks in the city compared to the suburbs; therefore, the human population density is higher in the city than in the suburbs, while the average vectorial density is lower.
  • Suburban areas also have problems with sanitation and intermittent water supplies, so the residents store potable water in big tanks, which provide additional breeding sites for mosquitoes.
  • In major Colombian cities located in dengue-endemic areas (such as Cali), the local health authorities implement routine vector control actions. These measures include a periodic larvicide treatment of rainwater catch basins in residential and commercial areas that are regarded as the main breeding sites of Aedes aegypti and other mosquito species. In the suburbs, however, such actions are either sporadic or absent.
The entries q i j , i , j = 1 , 2 of the time residence matrix Q will be further defined by different scenarios considered in the sequel. Notably, the residence times q i j will also determine the values of P i j and β h j in the dynamical system (9).
Let us now specify some plausible values for the entries of the objective functional (7). First, we have to establish the length of the observation period T > 0 . In this context, it is worthwhile to recall that, in dengue-endemic zones, corrective insecticide-spraying campaigns are usually conducted in response to threats of dengue outbreak. Their effect is not immediate and can be perceived in 2–4 weeks. Such interventions may last for a few days or weeks, focusing on specific areas with high mosquito density or reported clusters of disease cases [24]. Therefore, the time horizon will be defined as t [ 0 , T ] with T = 30 days and t = 0 marking the beginning of insecticide spraying.
Second, as the values of weight coefficients W i , i = 1 , 2 refer to the medical and societal costs of one human infection, we can assume that W 1 = W 2 , meaning that such expenses are the same in both patches. In [25], the authors have assessed that one ambulatory dengue infection in Colombia costs US$327, while one hospitalization with dengue infection leads to a higher cost of US$771. On the other hand, another study [26] estimated that 53% of dengue infections in Colombia require hospitalization, while the rest (47%) are cured with ambulatory treatment. Using these data, we can determine the mean cost of one dengue infection as
W 1 = W 2 = US $ 771 × 53 % ( inpatient ) + US $ 327 × 47 % ( outpatient ) = US $ 562.32 .
Using the above definition, the total cumulative cost of all dengue infections in both patches over T = 30 days is expressed by the first two terms of the objective functional (7).
Third, to calculate the total cost of the control intervention measures (defined by the last two terms of the functional (7)), we have to estimate the areas of both patches, A i , i = 1 , 2 , and establish the values of weight coefficients W 3 and W 4 .
According to [27], the urban area A 1 of Cali is 120 km 2 , and the overall urban area of the suburbs, A 2 , is obtained by summing up the urban areas of Palmira (23 km 2 ), Jamundi (42 km 2 ), and Yumbo (26 km 2 ) taken from the local sources [28,29,30]. Thus, we have that A 1 = 120 km 2 and A 2 = 91 km 2 .
The weights W 3 and W 4 express the costs related to insecticide spraying in both patches per unit area (that is, per 1 km 2 ), and their values may depend on different factors. It has already been mentioned in Section 3 that, by the definition of the control variables u i , i = 1 , 2 (see (6)), the total volume v i ( t ) of a diluted insecticide destined for application in Patch i on day t, as well as the maximal volume of insecticide allowed for daily application over a unit area at each patch (expressed by the constant D) depend on the insecticide lethality η .
Several studies addressing the costs of corrective vector control measures [31,32,33,34] point out that, besides the direct cost of the insecticide itself, this type of intervention should also account for other indirect and operational expenses (portable fogging equipment, trucks for spatial spraying, protective clothing and salaries for the personnel, etc.). In the absence of accurate information regarding the total costs of insecticide spraying in Cali, Colombia, we may assume the control strategies u i ( t ) ,   i = 1 , 2 induce costs proportional to the area treated with insecticide and use the data provided by the local contractor called SANICONTROL S.A.S. This local company (URL: https://www.sanicontrol.com/ (accessed on 12 August 2023)) is usually contracted by the local healthcare authorities to implement the corrective vector control measures in Cali and neighboring towns when an outbreak of the disease is detected. This company usually employs low-lethality insecticides bearing residual effects of about one day to comply with the local ambiental norms and provides price quoting per area [35]. Namely, the space spraying outside the houses and other buildings on a street-by-street basis of 1 km 2 of urban area costs around US$297–303 (depending on the variable exchange rate) when using an insecticide with a 40% lethality (meaning that η = 0.4 ). This is the total cost, including the insecticide and other operational and indirect expenses. Therefore, we may assume that W 3 = W 4 = US $ 300 in the objective functional (7).
On the other hand, more and more former residents of big cities are becoming accustomed to daily commuting after they move to nearby towns and suburbs, seeking to avoid overpopulation, heavy traffic, higher leasing costs, etc. What can be the impacts of daily commuting on the propagation of infectious diseases that are regarded as persistent and endemic? Does daily commuting increase or reduce the overall morbidity in the human population?
We will focus on the more realistic case regarding the people commuting between the two patches. As the suburban areas are mainly residential, they do not offer enough employment opportunities to their inhabitants. Therefore, a considerable share of the suburban residents are employed in the city and commute to work. Moreover, in the last decade, many former residents of Cali have moved out to the nearby towns (Jamundi, Palmira, and Yumbo) while keeping their jobs in the city to avoid overpopulation, pollution, higher leasing costs, heavy traffic, etc. This implies that the people’s flow from the suburbs (Patch 2) to the city (Patch 1) is significant, while the inverse flow is negligible. Thus, we will assume that the residents of Cali (Patch 1) remain ( q 12 = 0 ). We will also consider different intensities of the people flow from Patch 2 to Patch 1 ( 0 < q 21 < 0.5 and 0.5 < q 21 < 1 ) and compare the results with the uncoupled case ( q 11 = q 22 = 1 ). So, we have three basic scenarios, described as follows:
  • No commuting ⇒ uncoupled case denoted by S 0
  • Commuting  weak coupling ( S 1 ) strong coupling ( S 2 )
Table 2 provides the settings for the elements q i j ,   i , j = 1 , 2 of the residence time matrix Q for the cases S 1 and S 2 .

4.2. Optimal Solutions without Budget Constraints

To solve the optimal control problem (8) and (9), we assign the values of its parameters defined in the preceding subsection and assume that there is no constraint on insecticide use, meaning that the healthcare authorities dispose of sufficient budget to be spent on this type of corrective intervention. Furthermore, the healthcare authorities do not have any information regarding the directions and intensities of people commuting between the two patches. Therefore, their decisions are usually made without considering it, i.e., for the uncoupled case S 0 .
The results of the numerical solution of the optimal control problem (8) and (9) by the GPOPS-II solver are presented in Figure 1a, where both optimal controls u 1 * ( t ) and u 2 * ( t ) , drawn in blue and red colors, respectively, are of the bang-bang type. Their interpretations are rather straightforward. As u 1 * ( t ) drops from 1 to 0 at t ^ 1 = 6.1192 , it implies that the insecticide spraying in Patch 1 should go on for a bit more than 6 consecutive days at the maximum allowed spraying rate. In turn, u 2 * ( t ) drops from 1 to 0 at t ^ 2 = 4.3024 , meaning that, in practical terms, the insecticide spraying in Patch 2 should be performed for a little more than 4 consecutive days, also at the maximum allowed spraying rate.
Notably, numerical solutions of the optimal control problem (8) and (9) with the other two residence time configurations S 1 and S 2 (see Table 2) render optimal controls u 1 * ( t ) and u 2 * ( t ) that are almost identical to the ones presented in Figure 1a. The only difference is that the action of u 1 * ( t ) is reduced by up to 0.2 of a day (about four hours), while the action of u 2 * ( t ) is extended by a similar time. This outcome can be explained by the variability of effective vectorial densities (4) that express an average number of female mosquitoes in Patch i per one human individual effectively present in the same Patch i for i = 1 , 2 . As more suburban residents (Patch 2) commute to the city (Patch 1) than vice versa, the suburban area bears higher effective vectorial density, and a slightly longer intervention may be required.
As the healthcare authorities do not actually possess information regarding the directions and intensities of people’s commuting, the optimal controls derived for the uncoupled case S 0 (see Figure 1a) can be adopted in practice as the corrective measure for dengue reduction due to the resilience and robustness exhibited for other commuting scenarios S 1 , S 2 bearing different intensities.
The impact of the optimal controls u 1 * ( t ) and u 2 * ( t ) can be perceived by calculating the total number of new human infections that appear during the observation period [ 0 , T ] with control intervention versus the baseline case where u 1 ( t ) = u 2 ( t ) = 0 , t [ 0 , T ] . The corresponding epidemiological metric is referred to as cumulative incidence and is defined as
C h ( t ; u 1 , u 2 ) = C h 1 ( t ; u 1 , u 2 ) + C h 2 ( t ; u 1 , u 2 ) = N h 1 0 t κ h E h 1 ( ν ; u 1 , u 2 ) d ν + N h 2 0 t κ h E h 2 ( ν ; u 1 , u 2 ) d ν = 0 t κ h N h 1 E h 1 ( ν ; u 1 , u 2 ) + N h 2 E h 2 ( ν ; u 1 , u 2 ) d ν .
Here, E h i ,   i = 1 , 2 are solutions to the dynamical system (9) with specified u i ( t ) ,   i = 1 , 2 and C h i ( t ; u 1 , u 2 ) ,   i = 1 , 2 and also C h ( t ; u 1 , u 2 ) are non-decreasing positive real functions. Note that (18) only includes new infections that appear within the period [ 0 , t ] without accounting for infections that started before t = 0 . Figure 1b presents the graphs of C h ( t ; 0 , 0 ) (black curve) and C h ( t ; u 1 * , u 2 * ) (blue curve) corresponding to the baseline case and the optimal control intervention, respectively. Furthermore, in the absence of control intervention ( u 1 ( t ) = u 2 ( t ) = 0   , t [ 0 , T ] ), we have
C h 1 ( T ; 0 , 0 ) = 2965 , C h 2 ( T , 0 , 0 ) = 960 , C h ( T ; 0 , 0 ) = 3925 human infections
with an estimated total cost of US$2,207,106 using the average unit cost of US$562.32 per one infection as in (17) and no additional costs of the control intervention. On the other hand, when optimal controls u 1 * , u 2 * given in Figure 1a are in action, we have
C h 1 ( T ; u 1 * , u 2 * ) = 666 , C h 2 ( T , u 1 * , u 2 * ) = 257 , C h ( T ; u 1 * , u 2 * ) = 923 human infections .
Here, the total dengue-related public spending consists of two parts. The first part corresponds to the societal cost of 923 dengue infections in humans, and its estimation is US$519,021, assuming the average unit cost of US$562.32 per infection. The second part corresponds to the budget required for the control intervention in both patches. Let B 1 and B 2 denote the budgets in US$ required for implementing the optimal control strategies u 1 * ( t ) and u 1 * ( t ) , respectively. Then the total budget is
B : = B 1 + B 2 = 0 t ^ 1 W 3 A 1 u 1 * ( t ) d t + 0 t ^ 2 W 4 A 2 u 2 * ( t ) d t .
Using the above formula with W 3 , W 4 , A 1 , and A 2 defined in Section 4.1, we obtain the following estimations:
B 1 = 220 , 290 , B 2 = 117 , 450 , B = 337 , 740 .
Thus, by employing the optimal control strategies u 1 * , u 2 * given in Figure 1a, it is possible to cut the dengue-related public spending from US$2,207,106 to US$856,761, meaning more than a 2.5-fold reduction compared to the no-intervention case. However, this favorable result can only be guaranteed if the healthcare authorities have solvency for implementing the insecticide spraying in both patches. Yet, in the practical settings, the healthcare authorities may not possess sufficient funds to implement the optimal control strategy u 1 * , u 2 * .
Following the approach presented in [8], we address this issue in the next section, exploring different forms of budget reduction.

4.3. Optimal Solutions with Budget Constraints

In the previous subsection, we obtained that the necessary budget for implementing the optimal control strategy in both patches is B = US$337,740 (cf. (20)). Let us now revise the situation when only a share σ ( 0 , 1 ) of B is available for the insecticide spraying in both patches. When the local healthcare authorities face a lack of budget, they have different ways to deal with it. For instance, they may reduce uniformly by σ -share the budgets B 1 and B 2 . Alternatively, they may keep the budget B 1 necessary for implementing the optimal control u 1 * ( t ) while drastically reducing the budget B 2 , or vice versa. These three principal options can be formulated as follows:
Case 1. 
Uniform or unbiased budget cut in both patches σ B = σ B 1 + σ B 2 can be modeled by adding the following constraints to the optimal control problem (8) and (9):
0 T W 3 A 1 u 1 ( t ) d t = σ B 1 , 0 T W 4 A 2 u 2 ( t ) d t = σ B 2 .
Case 2. 
Keeping the required budget for Patch 1 with a drastic budget cut in Patch 2 only: σ B = B 1 + σ B 2 ( 1 σ ) B 1 . This situation is viable only if σ B 2 ( 1 σ ) B 1 > 0 and can be modeled by adding the following constraints to the optimal control problem (8) and (9):
0 T W 3 A 1 u 1 ( t ) d t = B 1 , 0 T W 4 A 2 u 2 ( t ) d t = σ B 2 ( 1 σ ) B 1 .
Case 3. 
Keeping the required budget for Patch 2 with a drastic budget cut in Patch 1 only: σ B = B 2 + σ B 1 ( 1 σ ) B 2 . This situation is viable only if σ B 1 ( 1 σ ) B 2 > 0 and can be modeled by adding the following constraints to the optimal control problem (8) and (9):
0 T W 3 A 1 u 1 ( t ) d t = σ B 1 ( 1 σ ) B 2 , 0 T W 4 A 2 u 2 ( t ) d t = B 2 .
It is worthwhile to point out that the most common way of dealing with budget cuts in practice is described by Case 2, while the rarest one is defined by Case 3. In the sequel, we can also validate what type of budget cut distribution will show better results by assessing the overall number of human infections in both patches.
Now we proceed to solve numerically the optimal control problem (8) and (9) with additional constraints (21)–(23) and under different residence time configurations (cf. Table 2). We have chosen σ 0.7 , 0.8 , 0.9 to comply with the conditions mentioned in Cases 2 and 3. Thus, we will consider the budget cuts of 10 % , 20 % , and 30 % with respect to to the required total budget B = US$337,740.
Numerical solutions obtained by GPOPS-II have once again culminated in optimal controls u 1 * ( t ) , u 2 * ( t ) that are almost identical for all three residence time configurations S 0 , S 1 , S 2 . The left column of Figure 2 presents the optimal controls u 1 * ( t ) , u 2 * ( t ) corresponding to the weak coupling ( S 1 ), while its right column displays the effect of these controls on the total cumulative incidence in both patches. The rows of Figure 2 correspond to Cases 1, 2, and 3 described above. The outcomes for the other two scenarios (uncoupled S 0 and with strong coupling S 2 ) look very much like the ones presented in Figure 2, and they are omitted here to save space. Nonetheless, the graphs of C h i ( t ; u 1 * , u 2 * ) ,   i = 1 , 2 are presented in Appendix B for all commuting scenarios S 0 , S 1 , S 2 , and the impact of commuting intensity on the number of new disease cases is briefly discussed.
By revising the left column of Figure 2, we notice that the periods of insecticide spraying ( u i * ( t ) = 1 ,   i = 1 , 2 ) are shortened proportionally to the budget reduction in all Cases 1–3. This feature was also detected in [8], where a more simple one-patch dengue transmission model of the Ross–MacDonald type was studied under resource limitation for insecticide spraying. On the other hand, the right column of Figure 2 brings forward some interesting insights. For a 10% budget cut (i.e., σ = 0.9 ), we perceive no difference in C h ( T ; u 1 * , u 2 * ) between Cases 1, 2, and 3. However, when the budget cuts become more significant (20% or 30%), it becomes visible that Case 1 (uniform or unbiased distribution of the budget reduction) renders better results than the two biased distributions (Cases 2 and 3). Notably, the most common way of dealing with budget cuts in practice (Case 2, middle chart in the right column of Figure 2) turns out to be the worst choice for dengue control.
In the next section, we bring forward more arguments about the cost of insecticide spraying (with and without budget cuts) and the expenses needed to cure human infections that will still arise after the corrective intervention measures.

4.4. Discussion

Even though dengue-related costs, both direct and indirect, vary significantly across different countries, regions, and socioeconomic contexts, the burden of dengue is often higher in resource-limited settings. Here, the direct costs principally refer to the implementation of vector control measures, both preventive and corrective, and the medical expenses assumed by the local public health units (such as diagnosing and treating dengue cases, including hospitalization, medications, laboratory tests, and specialized care for severe cases). On the other hand, the indirect costs encompass productivity losses by sickness leaves, work absenteeism, and other social consequences such as stress and emotional distress for infected individuals and their families that negatively affect people’s well-being.
In this section, we intend to estimate the direct costs for the three commuting scenarios S 0 , S 1 , S 2 (cf. Table 2) and under three different distributions of budget cuts, described in Section 4.3 as Cases 1, 2, and 3. We will only account for the costs of corrective control measures (insecticide spraying) conducted as an emergency response to suppress the disease outbreak while disregarding the costs of the preventive control measures, which are routinely applied regardless of the disease status. It is worthwhile to point out that the monetary assessment of the indirect costs is a rather challenging task that goes beyond the scope of our study.
The total direct cost related to the dengue outbreak during the observation period of T = [ 0 , 30 ] days will then include the following two components:
  • Spending for the insecticide spraying in both patches, expressed by σ B , where B is defined in (19) and σ 0.7 , 0.8 , 0.9 , 1 stands for the budget reduction factor.
  • Spending for the medical treatment of infected individuals residing in both patches, estimated using the Formulas (17), (18) as W 1 × C h ( T ; u 1 * , u 2 * ) , where u 1 * ( t ) and u 2 * ( t ) are solutions to the optimal control problem (8) and (9) with budget constraints (21)–(23) obtained for different commuting scenarios S 0 , S 1 , S 2 described in Table 2.
Figure 3 displays the total costs corresponding to Cases 1, 2, and 3 (by rows) under different commuting scenarios S 0 , S 1 , and S 2 (by columns) in the form of bar diagrams. Each bar sums up the two types of spending described above. A cursory glance at the bar diagrams in this figure reveals that cutting the resources for corrective control measures will always lead to a noticeable increase in the total cost. In other words, whatever is “saved” on the insecticide spraying will not suffice to cover the additional medical costs, and the total public spending will increase. Therefore, the best option is not to reduce the funds required for the corrective control measures because their reduction will always induce additional expenses related to illness costs.
Furthermore, we arrive at the following conclusion by comparing the three rows in Figure 3 corresponding to different budget cut distributions (Cases 1–3). If the healthcare authorities only have a fraction σ ( 0 , 1 ) of the budget B = B 1 + B 2 required for the corrective control measures, the uniform or unbiased budget reduction seems the most meaningful way to implement them. Moreover, Figure 3 also reveals that the habitual modus operandi for budget reduction (expressed by Case 2) is the most unsuitable option.
Finally, comparing the columns in Figure 3 that correspond to different commuting scenarios S 0 , S 1 , and S 2 , we observe no significant differences between them. In other words, the intensity of people commuting does not have a critical impact on the disease spread, at least when the necessary resources for corrective measures are fully available, meaning no budget reduction occurs. The numbers given in Appendix B (see Table A1, Table A2 and Table A3) indicate that under no budget reduction, a minimal number of expected new infections (merely one or two) can be caused by commuting. Nonetheless, this number increases progressively with the budget reduction.

5. Conclusions

In this paper, we have studied a realistic example of a two-patch dengue transmission model consisting of one big city (Patch 1) and its suburbs (Patch 2) whose residents daily commute to the city for work. It was also assumed that, even though both patches are located in the same dengue-endemic area, the city features lower vectorial density than the suburbs. An epidemic outbreak may occur when the routine preventive measures of vector control do not yield satisfactory results. In such a situation, the local healthcare authorities are compelled to implement urgent insecticide spraying, thus seeking to suppress the number of new infections in humans as quickly as possible. This corrective intervention does not replace the routine preventive measure of vector management because it may harm other non-target living organisms, and it is usually intended only for short-term application.
To design the patch-dependent optimal strategies for corrective control, we formulated an optimal control problem where the objective functional expresses the direct dengue-related costs to be minimized, namely, the costs destined for the insecticide spraying and the public spending for treating dengue infections in humans. Using numerical solutions to the formulated optimal control problem, we showed that the patch-dependent optimal strategies for insecticide spraying have almost no dependency on the intensity of people commuting but their effects significantly depend on the commuting intensity.
Regarding the societal disease-related costs, our approach illustrated that trying “to save money” by reducing the budget for corrective control is entirely unwise. Namely, the anticipated savings on insecticide spraying will further turn into considerable additional public spending for treating human infections, and such extra spending is avoidable by timely corrective intervention.
Furthermore, we have revised a situation when only a share of resources needed to implement the optimal control strategies is available and designed the optimal strategies for both patches under budget constraints. We also considered three different ways of distributing the reduced budget between the two patches. In practice, when there are not enough resources for insecticide spraying in the city and the suburbs, the local healthcare authorities first secure the necessary funds for the control measures in the city, and the rest are used in the suburbs. Our simulation has shown that this habitual modus operandi is the most unwise, for it yields the highest total costs. On the other hand, we have also corroborated that a uniform budget reduction for both patches is the wisest way of acting in resource-limited settings.
Summarizing the outcomes of the present work, we have provided mathematical and computational arguments to support the following recommendations:
  • The patch-dependent strategies for corrective vector control can be designed without knowing the intensity of people commuting between the patches. This finding simplifies the design of such strategies because the commuting intensity exhibits variability and uncertainties, while its estimation is complicated and requires additional efforts together with underlying societal costs.
  • Saving on corrective vector control measures during the ongoing epidemics is unwise. Such savings will always entail extra spending on disease treatment and enlarge the overall disease-related public expenditures.
  • In resource-limited settings, a uniform budget reduction for both patches is the wisest way to distribute the funding because it will produce smaller overall disease-related public expenditures than biased distribution scenarios.
In conclusion, we hope that the outcomes of this study will help in designing patch-dependent strategies for corrective intervention and also send a message to the public healthcare authorities in dengue-endemic countries regarding the necessary adjustments in resource-limited settings.

Author Contributions

Conceptualization, E.B.-R. and O.V.; methodology, E.B.-R., O.V. and M.S.; Software, E.B.-R. and M.S.; validation, O.V.; formal analysis, E.B.-R. and O.V; investigation, E.B.-R.; data curation, E.B.-R.; writing—original draft preparation, E.B.-R., O.V. and M.S.; writing—review and editing, O.V. and M.S.; visualization, E.B.-R.; supervision, O.V. and M.S.; project administration, O.V.; funding acquisition, O.V. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Fund for Science, Technology, and Innovation (Autonomous Heritage Fund Francisco José de Caldas) by way of the Research Program No. 1106-852-69523 (Principal Investigator: Hector J. Martinez-Romero), Contract: CT FP 80740-439-2020 (Colombian Ministry of Science, Technology, and Innovation—Minciencias), Grant ID: CI-71241 (Universidad del Valle, Colombia).

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Acknowledgments

The authors would like to thank Carmen Alicia Ramírez Bernate (Universidad del Valle, Colombia) for helping to obtain the real price quotes for the insecticide spraying from SANICONTROL S.A.S. This added considerable practical value to our study and made it more realistic.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ODEOrdinary Differential Equation(s)
DENVdengue virus
SEIR(S)Susceptible–Exposed–Infectious–Recovered (Susceptible)
SEISusceptible–Exposed–Infectious
GPOPSGeneral Purpose Optimal Control Software
MATLABMatrix Laboratory (matrix programming language)
S.A.S.Simplified Stock Company (Spanish: Sociedad por Acciones Simplificadas)

Appendix A

The adjoint ODE system (13) has the following form:
λ ˙ 1 i = H S v i = ( μ v + η u i ) λ 1 i + β v j = 1 2 P j i I h j ( λ 1 i λ 2 i ) λ ˙ 2 i = H E v i = ( μ v + κ v + η u i ) λ 2 i κ v λ 3 i λ ˙ 3 i = H I v i = ( μ v + η u i ) λ 3 i + β β h i q 1 i ( λ 41 λ 51 ) S h 1 + β h i q 2 i ( λ 42 λ 52 ) h 2 λ ˙ 4 i = H S h i = μ h λ 4 i + j = 1 2 β h j q i j I v j ( λ 4 i λ 5 i ) λ ˙ 5 i = H E h i = ( μ h + κ h ) λ 5 i κ h λ 6 i λ ˙ 6 i = H I h i = β v P i 1 S v 1 ( λ 11 λ 21 ) + β v P i 2 S v 2 ( λ 12 λ 22 ) + ( γ + μ h ) λ 6 i γ λ 7 i W i N h i λ ˙ 7 i = H R h i = τ λ 4 i + ( μ h + τ ) λ 7 i
with i = 1 , 2 .

Appendix B

Figure 2, presented in Section 4.3, only shows the graphs of the total cumulative incidence C h ( t ; u 1 * , u 2 * ) (left column) for the weak coupling scenario S 1 . To showcase the impact of people’s commuting intensity on the disease incidence, we present in the sequel more detailed charts and numbers corresponding to three commuting scenarios S 0 , S 1 , S 2 , where the cumulative incidences are separated into C h i ( t ; u 1 * , u 2 * ) ,   i = 1 , 2 and where the three distribution modes of budget reduction are accounted for.
As we can see, the black curves in all charts of Figure A1, Figure A2 and Figure A3, corresponding to the no-budget-cut situation, are essentially the same. Let us now revise consecutively Cases 1, 2, and 3 referring to different budget cut distributions while paying more attention to the effects of people’s commuting on the expected number of new infections in both patches. These numbers are presented in Table A1, Table A2 and Table A3.
Figure A1. Uncoupled scenario S 0 : cumulative incidences C h i ( t ; u 1 * , u 2 * ) ,   i = 1 , 2 ,   t [ 0 , T ] in Patch 1 (left column) and Patch 2 (right column).
Figure A1. Uncoupled scenario S 0 : cumulative incidences C h i ( t ; u 1 * , u 2 * ) ,   i = 1 , 2 ,   t [ 0 , T ] in Patch 1 (left column) and Patch 2 (right column).
Mathematics 11 03921 g0a1
First, we contrast the charts from the first rows in Figure A1, Figure A2 and Figure A3 (Case 1: uniform or unbiased budget cut distribution). Here, we observe almost no differences for the moderate budget reduction of 10% ( σ = 0.9 ). However, when the budget cuts become more significant (20% or 30%), it becomes visible that the number of new infections in Patch 1 decreases as the intensity of commuting increases from S 0 to S 1 , and further S 2 (cf. Figure A1, Figure A2 and Figure A3, left charts in the first rows). At the same time, the number of infections in Patch 2 increases with the intensity of coupling (see exact numbers in the upper rows of Table A1, Table A2 and Table A3).
Table A1. Uncoupled scenario S 0 : total number of new infections in Patch 1 and Patch 2 expected over 30 days.
Table A1. Uncoupled scenario S 0 : total number of new infections in Patch 1 and Patch 2 expected over 30 days.
Scenario S 0 New Infections in 30 Days
Patch 1Patch 2Both
Case 1100%666257923
90%708288996
80%7613281089
70%8293781207
Case 2100%666257923
90%6663431009
80%6664941160
70%6667631429
Case 3100%666257923
90%7472571004
80%8712571128
70%10602571317
Now we contrast the middle rows in Figure A1, Figure A2 and Figure A3 (Case 2: keeping the required budget for Patch 1 with a drastic budget cut in Patch 2 only). Even though the left-hand charts look almost the same, we can still spot the effect of people commuting; the total number of new infections in Patch 1 slightly decreases as the intensity of commuting increases from S 0 to S 1 , and further S 2 . This tendency is also confirmed by the numbers in the middle lows of Table A1, Table A2 and Table A3.
Figure A2. Weak coupling scenario S 1 : cumulative incidences C h i ( t ; u 1 * , u 2 * ) ,   i = 1 , 2 ,   t [ 0 , T ] in Patch 1 (left column) and Patch 2 (right column).
Figure A2. Weak coupling scenario S 1 : cumulative incidences C h i ( t ; u 1 * , u 2 * ) ,   i = 1 , 2 ,   t [ 0 , T ] in Patch 1 (left column) and Patch 2 (right column).
Mathematics 11 03921 g0a2
On the other hand, the situation in Patch 2 is the opposite. By comparing the right-hand charts of the middle rows in Figure A1, Figure A2 and Figure A3, we observe that more intense commuting results in a greater number of new infections among the residents of Patch 2, while the drastic budget reductions only exacerbate this tendency. This conclusion is illustrated by the number presented in the middle rows of Table A1, Table A2 and Table A3.
Table A2. Weak coupling scenario S 1 : total number of new infections in Patch 1 and Patch 2 expected over 30 days.
Table A2. Weak coupling scenario S 1 : total number of new infections in Patch 1 and Patch 2 expected over 30 days.
Scenario S 1 New Infections in 30 Days
Patch 1Patch 2Both
Case 1100%632292924
90%671328999
80%7203721092
70%7834291212
Case 2100%632292924
90%6333801013
80%6335321165
70%6348031437
Case 3100%632292924
90%7062991005
80%8193091128
70%9923251317
Figure A3. Strong coupling scenario S 2 : cumulative incidences C h i ( t ; u 1 * , u 2 * ) ,   i = 1 , 2 ,   t [ 0 , T ] in Patch 1 (left column) and Patch 2 (right column).
Figure A3. Strong coupling scenario S 2 : cumulative incidences C h i ( t ; u 1 * , u 2 * ) ,   i = 1 , 2 ,   t [ 0 , T ] in Patch 1 (left column) and Patch 2 (right column).
Mathematics 11 03921 g0a3
Table A3. Strong coupling scenario S 2 : total number of new infections in Patch 1 and Patch 2 expected over 30 days.
Table A3. Strong coupling scenario S 2 : total number of new infections in Patch 1 and Patch 2 expected over 30 days.
Scenario S 2 New Infections in 30 Days
Patch 1Patch 2Both
Case 1100%617308925
90%6553451000
80%7033921095
70%7634511214
Case 2100%617308925
90%6183961014
80%6195491168
70%6208211441
Case 3100%617308925
90%6883181006
80%7973331130
70%9633561319
Finally, by comparing the charts in the last rows of Figure A1, Figure A2 and Figure A3 (Case 3: keeping the required budget for Patch 2 with a drastic budget cut in Patch 1 only), we observe tendencies similar to Case 2. Namely, as the intensity of commuting increases from S 0 to S 1 , and further S 2 , the number of infections decreases in Patch 1 and increases in Patch 2. Moreover, the role of people commuting from Patch 2 to Patch 1 becomes more visible. Even though there is no budget reduction for corrective intervention in Patch 2, the suburban residents commuting to Patch 1 are more likely to be infected in the city where the intervention measures are less efficient due to budget cuts. This situation is reflected in the right charts of the last rows in Figure A2 and Figure A3 corresponding to Case 3, and also in the lower rows of Table A1, Table A2 and Table A3. Thus, the budget cuts in Patch 1 contribute to the disease propagation not only in the city itself but also in the suburban areas through the people commuting.

References

  1. Barrios, E.; Lee, S.; Vasilieva, O. Assessing the effects of daily commuting in two-patch dengue dynamics: A case study of Cali, Colombia. J. Theor. Biol. 2018, 453, 14–39. [Google Scholar] [CrossRef] [PubMed]
  2. Bichara, D.; Castillo-Chavez, C. Vector-borne diseases models with residence times—A Lagrangian perspective. Math. Biosci. 2016, 281, 128–138. [Google Scholar] [CrossRef] [PubMed]
  3. Lee, S.; Castillo-Chavez, C. The role of residence times in two-patch dengue transmission dynamics and optimal strategies. J. Theor. Biol. 2015, 374, 152–164. [Google Scholar] [CrossRef] [PubMed]
  4. Kim, J.; Lee, H.; Lee, C.; Lee, S. Assessment of optimal strategies in a two-patch dengue transmission model with seasonality. PLoS ONE 2017, 12, e0173673. [Google Scholar] [CrossRef]
  5. Bichara, D.; Iggidr, A. Multi-patch and multi-group epidemic models: A new framework. J. Math. Biol. 2018, 77, 107–134. [Google Scholar] [CrossRef]
  6. Mishra, A.; Gakkhar, S. Non-linear dynamics of two-patch model incorporating secondary dengue infection. Int. J. Appl. Comput. Math. 2018, 4, 1–22. [Google Scholar] [CrossRef]
  7. Lasluisa, D.; Barrios, E.; Vasilieva, O. Optimal strategies for dengue prevention and control during daily commuting between two residential areas. Processes 2019, 7, 197. [Google Scholar] [CrossRef]
  8. Sepulveda-Salcedo, L.S.; Vasilieva, O.; Svinin, M. Optimal control of dengue epidemic outbreaks under limited resources. Stud. Appl. Math. 2020, 144, 185–212. [Google Scholar] [CrossRef]
  9. Jan, R.; Khan, M.; Gomez-Aguilar, J. Asymptomatic carriers in transmission dynamics of dengue with control interventions. Optim. Control. Appl. Methods 2020, 41, 430–447. [Google Scholar] [CrossRef]
  10. Rodriguez-Barraquer, I.; Solomon, S.; Kuganantham, P.; Srikrishnan, A.; Vasudevan, C.; Iqbal, S.; Balakrishnan, P.; Solomon, S.; Mehta, S.; Cummings, D. The hidden burden of dengue and Chikungunya in Chennai, India. PLoS Neglected Trop. Dis. 2015, 9, e0003906. [Google Scholar] [CrossRef]
  11. Thai, K.; Nishiura, H.; Hoang, P.; Tran, N.; Phan, G.; Le, H.; Tran, B.; Nguyen, N.; de Vries, P. Age-specificity of clinical dengue during primary and secondary infections. PLoS Neglected Trop. Dis. 2011, 5, e1180. [Google Scholar] [CrossRef] [PubMed]
  12. Deng, S.Q.; Yang, X.; Wei, Y.; Chen, J.T.; Wang, X.J.; Peng, H.J. A review on dengue vaccine development. Vaccines 2020, 8, 63. [Google Scholar] [CrossRef] [PubMed]
  13. Al-Zahrani, M.; Gharsan, F.; Al-Ghamd, K.; Mahyoub, J.; Alghamdi, T. Toxicity of two groups of pesticides against the mosquito Aedes aegypti. GSC Biol. Pharm. Sci. 2020, 13, 148–155. [Google Scholar] [CrossRef]
  14. Coddington, E.; Levinson, N. Theory of Ordinary Differential Equations; McGraw-Hill: New York, NY, USA, 1987. [Google Scholar]
  15. Walter, W. Ordinary differential equations. In Graduate Texts in Mathematics; Springer: Berlin/Heidelberg, Germany, 1998; Volume 182. [Google Scholar]
  16. Fleming, W.; Rishel, R. Deterministic and Stochastic Optimal Control; Springer: New York, NY, USA, 1975. [Google Scholar]
  17. Lenhart, S.; Workman, J.T. Optimal Control Applied to Biological Models; Chapman & Hall/CRC: Boca Raton, FL, USA, 2007. [Google Scholar]
  18. Harmand, J.; Lobry, C.; Rapaport, A.; Sari, T. Optimal Control in Bioprocesses: Pontryagin’s Maximum Principle in Practice; Wiley: London, UK, 2019. [Google Scholar]
  19. Rao, A.; Benson, D.; Darby, C.; Patterson, M.; Francolin, C.; Sanders, I.; Huntington, G. Algorithm 902: GPOPS, a MATLAB software for solving multiple-phase optimal control problems using the Gauss pseudospectral method. ACM Trans. Math. Softw. 2010, 37, 22. [Google Scholar] [CrossRef]
  20. Patterson, M.; Rao, A. GPOPS-II: A MATLAB software for solving multiple-phase optimal control problems using hp-adaptive Gaussian quadrature collocation methods and sparse nonlinear programming. ACM Trans. Math. Softw. 2014, 41, 1. [Google Scholar] [CrossRef]
  21. Garg, D.; Patterson, M.; Francolin, C.; Darby, C.; Huntington, G.; Hager, W.; Rao, A. Direct trajectory optimization and costate estimation of finite-horizon and infinite-horizon optimal control problems using a Radau pseudospectral method. Comput. Optim. Appl. 2011, 49, 335–358. [Google Scholar] [CrossRef]
  22. Desjardins, M.; Eastin, M.; Paul, R.; Casas, I.; Delmelle, E. Space–Time Conditional Autoregressive Modeling to Estimate Neighborhood-Level Risks for Dengue Fever in Cali, Colombia. Am. J. Trop. Med. Hyg. 2020, 103, 2040. [Google Scholar] [CrossRef]
  23. Padilla, J.; Rojas, D.; Gómez, R. Dengue en Colombia: Epidemiología de la Reemergencia a la Hiperendemia; Dengue in Colombia: Epidemiology of Hyperendemic Re-emergence; Guías de Impresión Ltda.: Bogota, Colombia, 2012. [Google Scholar]
  24. Esu, E.; Lenhart, A.; Smith, L.; Horstick, O. Effectiveness of peridomestic space spraying with insecticide on dengue transmission; Systematic review. Trop. Med. Int. Health 2010, 15, 619–631. [Google Scholar] [CrossRef]
  25. Rodríguez, R.; Rueda-Gallardo, J.; Avella-Niño, M. The burden of dengue illness and its economics costs in the Americas: A review on the most affected countries. In Dengue Fever—A Resilient Threat in the Face of Innovation; Falcón-Lezama, J., Betancourt-Clavioto, M., Tapis-onyer, R., Eds.; IntechOpen: London, UK, 2019; Chapter 2. [Google Scholar]
  26. Hernández, L.; Durán, D.; Buitrago, D.; Garnica, C.; Gómez, L.; Bados, D.; Bernal, M.; Páez, L. Epidemiology and geo-referencing of the dengue fever in a hospital of second level in Colombia, 2010–2014. J. Infect. Public Health 2018, 11, 558–565. [Google Scholar] [CrossRef]
  27. Escobar-Morales, G. (Ed.) Cali en Cifras 2013; Departamento Administrativo de Planeación, Alcaldía de Santiago de Cali: Cali, Colombia, 2014. [Google Scholar]
  28. Agenda Ambiental del Municipio de Palmira; Alcaldía de Palmira; Ministerio del Medio Ambiente: Bogota, Colombia, 2002. Available online: http://documentacion.ideam.gov.co/openbiblio/bvirtual/005574/home/Pal.pdf (accessed on 12 August 2023).
  29. Plan Municipal de Gestión del Riesgo de Desastres, Alcaldía de Jamundí. 2012. Available online: https://1library.co/document/q7wg71oz-municipal-gestion-riesgo-desastres-municipio-jamundi-valle-cauca.html (accessed on 12 August 2023).
  30. Lesmes-Duque, M.C.; Peña, P.L.; Abreu, C.P. Análisis de Situación de Salud con el Modelo de los Determinantes Sociales de Salud. 2018. Available online: https://www.valledelcauca.gov.co/loader.php?lServicio=Tools2&lTipo=viewpdf&id=29565 (accessed on 12 August 2023).
  31. Orellano, P.W.; Pedroni, E. Cost-benefit analysis of vector control in areas of potential dengue transmission. Pan Am. J. Public Health 2008, 24, 113–119. [Google Scholar] [CrossRef]
  32. Baly, A.; Toledo, M.; Rodriguez, K.; Benitez, J.; Rodriguez, M.; Boelaert, M.; Vanlerberghe, V.; van der Stuyft, P. Costs of dengue prevention and incremental cost of dengue outbreak control in Guantanamo, Cuba. Trop. Med. Int. Health 2012, 17, 123–132. [Google Scholar] [CrossRef]
  33. Packierisamy, P.R.; Ng, C.W.; Dahlui, M.; Venugopalan, B.; Halasa, Y.A.; Shepard, D.S. The cost of dengue vector control activities in Malaysia by different service providers. Asia Pac. J. Public Health 2015, 27, 73S–78S. [Google Scholar] [CrossRef] [PubMed]
  34. Thompson, R.; Martin Del Campo, J.; Constenla, D. A review of the economic evidence of Aedes-borne arboviruses and Aedes-borne arboviral disease prevention and control strategies. Expert Rev. Vaccines 2020, 19, 143–162. [Google Scholar] [CrossRef] [PubMed]
  35. Ramírez-Bernate, C.; Martínez-Romero, H.; Erazo-Borja, D. Control strategies in the spatial population dynamics of Aedes aegypti vector using sterile mosquitoes and insecticides. Univ. Sci. 2022, 27, 206–232. [Google Scholar] [CrossRef]
Figure 1. (a) Optimal controls u 1 * ( t ) and u 2 * ( t ) of bang-bang type for the uncoupled case S 0 ; (b) total cumulative incidence in both patches.
Figure 1. (a) Optimal controls u 1 * ( t ) and u 2 * ( t ) of bang-bang type for the uncoupled case S 0 ; (b) total cumulative incidence in both patches.
Mathematics 11 03921 g001
Figure 2. Weak coupling configuration S 1 : optimal controls u 1 * ( t ) and u 2 * ( t ) of bang-bang type under budget cuts (left column); total cumulative incidence C h ( t ; u 1 * , u 2 * ) ,   t [ 0 , T ] in both patches (right column).
Figure 2. Weak coupling configuration S 1 : optimal controls u 1 * ( t ) and u 2 * ( t ) of bang-bang type under budget cuts (left column); total cumulative incidence C h ( t ; u 1 * , u 2 * ) ,   t [ 0 , T ] in both patches (right column).
Mathematics 11 03921 g002
Figure 3. Total costs corresponding to Cases 1, 2, and 3 (by rows) under different commuting scenarios S 0 , S 1 , and S 2 (by columns).
Figure 3. Total costs corresponding to Cases 1, 2, and 3 (by rows) under different commuting scenarios S 0 , S 1 , and S 2 (by columns).
Mathematics 11 03921 g003
Table 1. Description of the model’s parameters with values borrowed from [1] (epidemic scenario).
Table 1. Description of the model’s parameters with values borrowed from [1] (epidemic scenario).
ParameterDescriptionUnitsAssumed Value
κ v Latent-to-infectious conversion rate in vectorsday 1 1 / 8
κ h Latent-to-infectious conversion rate in humansday 1 1 / 4
γ Rate of human recoveryday 1 1 / 6
μ v Per capita vector birth/death rateday 1 1 / 30
μ h Per capita of human births/deathsday 1 1 / ( 75 × 365 )
α Mosquito biting rateday 1 0.4517
τ Reinfection rate for human individualsday 1 0.0262
p v Human-to-vector transmission probability of the pathogen 0.2378
p h Vector-to-human transmission probability of the pathogen 0.219
Table 2. Configurations of the elements q i j ,   i , j = 1 , 2 of the residence time matrix Q assumed for the numerical simulations.
Table 2. Configurations of the elements q i j ,   i , j = 1 , 2 of the residence time matrix Q assumed for the numerical simulations.
ScenariosCouplingResidence Times
S 0 None q 11 = 1 , q 12 = 0
q 21 = 0 , q 22 = 1
S 1 Weak q 11 = 1 , q 12 = 0
q 21 = 0.4 , q 22 = 0.6
S 2 Strong q 11 = 1 , q 12 = 0
q 21 = 0.6 , q 22 = 0.4
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Barrios-Rivera, E.; Vasilieva, O.; Svinin, M. Optimal Control of a Two-Patch Dengue Epidemic under Limited Resources. Mathematics 2023, 11, 3921. https://doi.org/10.3390/math11183921

AMA Style

Barrios-Rivera E, Vasilieva O, Svinin M. Optimal Control of a Two-Patch Dengue Epidemic under Limited Resources. Mathematics. 2023; 11(18):3921. https://doi.org/10.3390/math11183921

Chicago/Turabian Style

Barrios-Rivera, Edwin, Olga Vasilieva, and Mikhail Svinin. 2023. "Optimal Control of a Two-Patch Dengue Epidemic under Limited Resources" Mathematics 11, no. 18: 3921. https://doi.org/10.3390/math11183921

APA Style

Barrios-Rivera, E., Vasilieva, O., & Svinin, M. (2023). Optimal Control of a Two-Patch Dengue Epidemic under Limited Resources. Mathematics, 11(18), 3921. https://doi.org/10.3390/math11183921

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop