Dynamic analysis and optimal control of Zika virus transmission with immigration

: In this paper, a type of Zika virus model with immigration is considered. Additionally based on the risk of infected immigrants, we propose a control measure of screening for immigrants and a three-measure control model of combined mosquito prevention and killing. The existence and stability of the equilibrium in the Zika virus model are analyzed. The necessary conditions for the existence of the optimal solution are given using Pontryagin’s maximum principle. We focused on testing screening of the immigrating population to ensure a reduction in the transmission of the virus. Models have demonstrated that in combination with routine mosquito control measures and the appropriate use of mosquitoicides, the transmission of Zika virus in the population can be e ﬀ ectively reduced.


Introduction
Zika virus belongs to flavivirus genus, flavivirus family, which is an arbovirus transmitted by mosquitoes, mostly transmitted by Aedes aegypti, Aedes albopictus and Culex mosquitoes, etc. Zika virus was first detected in rhesus monkeys in the Zika jungle in Uganda in 1947 [1], and then in 1952, first found in populations in Uganda [2]. During 2015 and 2016, Brazil reported more than 4,000 cases of microcephalic teratoma in pregnant women infected with Zika virus [3,4], a 20 fold increase compared to previous years. By early 2016, Zika virus has spread to South America, Asia, Oceania and other regions [5][6][7]. According to the World Health Organization, 86 countries and territories have experienced outbreaks of the Zika virus since the outbreak began.
Human immigration is a very normal phenomenon. The movement or travel of a group of people, either from an endemic area to a healthy area or from a healthy area to an endemic area. Therefore, population movement or immigration is the main cause of the rapid spread of infectious diseases. Fred Brauer et al. [8] considered the immigration of infected individuals based on the traditional SI and SIS models of infectious diseases and assumed a constant rate of population immigration. Their results suggest that isolation of migrating infected individuals is necessary. Molalegn Ayana and Purnachandra Rao Koya [9] considered the impact of having immigrants infected with Zika virus through a model and proposed that infected immigrants facilitate the spread of Zika virus. Traoré et al. [10] analyzed a vector-borne infectious disease model that takes into account vector and human immigration. Their results also indicate that human migration is a problem that cannot be ignored and may cause the spread of disease in non-infected areas. So what control measures are more appropriate for this situation? Should we take immigration testing or require pre-migration vaccination?
There are many scholars who have focused their research on optimal cost-effectiveness and cost analysis. Kouidere et al. [11] studied a mathematical model that describes the transmission dynamics of African swine fever virus (ASFV) between pigs and ticks. To reduce the number of infected pigs and ticks, several strategies are proposed, and Pontryagin's Maximum principle and cost analysis are used to find the solution of optimal control. In 2022, they proposed a mathematical model in another paper [12] to describe the spread of COVID-19 in Peru, and characterized the optimal control through Pontryagin's Maximum principle. Using an optimal control model, Abdulfatai and Armin [13] concluded that by comparing treatment of symptomatic infected individuals and indoor residual spraying is the most cost-effective strategy. Miyaoka et al. [14] developed a response-diffusion model of Zika transmission, suggesting that the best control strategy is to immunize susceptible populations with vaccination as the control variable. Bonyah Ebenezer et al. [15] proposed and analyzed a SEIR-Zika epidemic model and established an optimal control model. They only proposed the optimal control strategy and did not mention the cost. Similarly, there are a number of studies that consider only optimal control strategies [16][17][18][19]. Screening controls for in-migrants are largely absent from these control strategies. As of now, there is some wait time for a Zika virus vaccine to become available, but testing for Zika virus is currently available. The major contributions of this work are as follows: A Zika virus model with immigrants is proposed to explain the risk of virus transmission by immigrants. We also propose appropriate screening measures for immigrants to find the optimal control scheme by building an optimal control model as well as cost analysis, pointing out the feasibility of appropriate screening in Zika virus control.
The structure of this paper is as follows. In the next section, the Zika virus model with immigration is proposed. In Section 3, we discussed the stability of the equilibriums point. In Section 4, combining with the actual situation, the control strategy is proposed and the optimal control model is established. The optimal control analysis is carried out. Cost-effectiveness analysis is given in Section 5. Then, the conclusion and discussion will be made in Section 6.

Establishment of the model
Considering that mosquitoes move in a relatively small distance, we overlook the immigration of mosquitoes. We divide the human population into four sub-classes, namely susceptible humans S H , exposed humans E H , infected humans I H , and recovered humans R H . By this virtue, the total human population can be represented by: Similarly, N M (t) is the total number of mosquitoes, which can be rationally partitioned into susceptible mosquitoes S M , exposed mosquitoes E M , infected mosquitoes I M . Hence, One problem for the model proposed by Molalegn Ayana and Purnachandra Rao Koya [9] is that if there are no infected people in the current infected population, then no one in the immigrant population is infected, which is clearly inappropriate. Imagine a city in which there are no infected people, but if people infected by the virus move in, then that must be one reason why the virus would be present in that city. Therefore, based on the idea of [10], the following model developed: where Λ H denotes the recruitment of humans, which also means the rate at the new individuals enter human population. Of these new individuals, we hypothesize that a fraction of 1 − P E − P I − P R are susceptibles, P E are exposed, P I are infected and P R are are recovered. Λ M represents susceptible mosquitoes recruitment, β H signifies the direct transmission rate of the disease from infectious mosquitoes to susceptible humans, β M denotes the rate of transmission from infected humans to susceptible mosquitoes, β H ρ represents the rate of transmission from infected humans to susceptible humans, α H stands for the rate of exposed humans moving into infectious class, δ M is the rate flow from E M to I M . Natural mortality levels associated with each subpopulation of humans and mosquitoes are denoted by µ H and µ M , respectively, r is the natural recovery rate, and η denotes the treatment rate. All the parameters here are positive. Figure 1 shows the compartment model, which clearly shows the construction process of the model. Table 1 defines all the parameters of the model.  The growth rate of immigration 1 − P E − P I − P R The proportion of susceptible immigrant P E , P I , P R The proportion of exposed persons infected persons and recovered persons β H The rate of transmission of the virus from mosquitoes to humans ρ The rate of transmission of the virus from person to person µ H The natural mortality rate of people α H The transfer rate from exposed to infected persons r Human natural recovery rate η Recovery rate due to treatment Λ M Growth rate of susceptible mosquitoes β M The rate of transmission of the virus from humans to mosquitoes µ M Mosquito natural mortality rate δ M The transfer rate of exposed mosquitoes to infected mosquitoes We can get that Ω is the positive invariant set of the system (2.1).

The disease-free equilibrium and the basic reproduction number
When P E = 0, P I = 0, that means that no infected or exposed person moves in, the model can be written as (3.1) The disease-free equilibrium point of system . By the next generation operator method, F and V are respectively The basic reproduction number can be obtained by ρ(FV −1 ), that is , .

The stability of the disease-free equilibrium
The Jacobi matrix of system (3.1) at E 0 is The Jacobi matrix shows that the three eigenvalues are −µ H , −µ H , −µ M . The remaining four eigenvalues can be determined by the following characteristic equation where .
From Lemma 1, we know that 2R 1 +R 2 < 1 if and only if R 0 < 1 holds. Therefore all the order principal minor determinants are positive, when K i , i = 1, 2, 3, 4. Thus, when R 0 < 1, according to Hurwitz's criterion, system (3.1) is locally asymptotically stable at the equilibrium point E 0 ; otherwise, it is unstable. Next we give the global asymption stability of DFE.
Theorem 2. The disease-free equilibrium point E 0 of system (3.1) is globally asymptotically stable if R 0 < 1; otherwise, it is unstable.
Proof. Define the Lyapunov function The time derivative of the Lyapunov function is The disease-free equilibrium point of system We get Therefore the largest compact invariant set in Ω is the singleton set E 0 . According to LaSalle's invariance principle [20], the disease-free equilibrium point E 0 of system (3.1) is globally asymptotically stable if R 0 < 1; otherwise, it is unstable.

The existence of the endemic equilibrium
Obviously, when P E 0, P I 0, there is no disease-free equilibrium. Let each equation of system (2.1) equal to 0, we can get Let the endemic equilibrium of system (3.2) be E 1 = (S H1 , E H1 , I H1 , R H1 , S M1 , E M1 , I M1 ). Denoting It can be derived from Eq (3.3) that . (3.5) According to Eqs (3.4) and (3.5), we assume that The derivatives of Eqs (3.6) and (3.7) are taken as follows Continue with the second order derivative as follows Because of f 2 (I H ) < 0, f 1 (I H ) > 0, f (I H ) > 0 holds. That admits that f (I H ) increases monotonically as I H falling in (0, I 1 ). Let I H = 0, then I) If P E , P I > 0, the system has a unique endemic equilibrium point. II) If P E , P I = 0, R 0 > 1, the system has a unique endemic equilibrium point [15].

The stability of the endemic equilibrium
In order to simplify the system, system (2.1) can be written as the following equivalent system (3.9) Theorem 4. For system (3.9): I) If P E , P I > 0, the only endemic equilibrium of the system is globally asymptotically stable. II) If P E , P I = 0, R 0 > 1, the only endemic equilibrium of the system is globally asymptotically stable.
Proof. For conclusion II), Theorem 5.2 in Reference [15] has been proved, now we only prove conclusion I).
Define the Lyapunov function as The total derivative of the Lyapunov function with respect to the system (3.9) is (3.11) Calculated separately, we can get

Optimal control analysis
Optimal control theory is used to determine the method to achieve the minimum cost and maximum performance under various assumptions [21,22]. In this section, to reduce infection vectors and control the spread of diseases, based on the previous analysis and combined with reality, several control strategies were proposed.
The two most common ways to control Zika virus are the efforts on preventing Zika infections (e.g. using mosquito nets, condoms, and so on) and the efforts on eliminating mosquitoes by insecticides. For the risk of someone carrying the virus in the immigrant population, we propose a measure of virus screening.
It is important to note that the control variables are primarily acting on the parameters corresponding to each control measure. For example, measures to prevent mosquito bites will change the rate of mosquito bite infection, so the reduction in infection rate is expressed in the form ( Similarly, if a mosquito killer is to be used, then it is to some extent increasing the mortality rate of mosquitoes, so expressed in the form (1 + u 3 )µ M . The implementation of virus screening will reduce the proportion of exposure and infection in the immigrant population. So, we denote it by the form (1 − u 2 )P E , (1 − u 2 )P I . Therefore, our optimal control model is given by where u 1 , u 2 , u 3 are described in Table 2. Efforts on preventing zika infections through mosquito nets, condoms, and so on. u 2 Efforts on reducing the rate of infected populations in immigration people by screening of migrant populations. u 3 Efforts on harvesting through spray insecticide.
In general, if the basic reproduction number exists, the control scheme can be proposed by limiting the basic reproduction number to less than 1 and backtracking the values of the control parameters. However, in this model, if P I 0, P E 0, there is no the basic reproduction number, so the above method cannot be used to study the control measures. Therefore, in this paper, the transmission model is first qualitatively analyzed to reveal the dynamics behavior of continuous virus transmission. Then, the Pontryagin maximum principle [23] will be used to find the optimal control scheme to achieve virus control in a limited time.
The objective function is defined as follows where ω E H , ω I H , ω E M , ω I M , a 1 , a 2 , a 3 represents the weight coefficients of the control variables, which are designed to maintain a balance among the items of the integration function so that no dominant individual term emerges. t f is the terminal moment when the control policy is implemented. Our goal is to find a set of control parameters (u * 1 , u * 2 , u * 3 ) satisfying It should be noticed that the right side of model (4.2) is bounded and the objective function is convex. A sufficient condition for the existence of optimal control indicates that the model has optimal control. The Lagrange function is

Define the Hamilitonian function as
From Pontryagin extreme value principle, the control set u(t) should satisfy the following necessary conditions: Theorem 5. There is a set of u 1 , u 2 , u 3 , so that J u * 1 , u * 2 , u * 3 = min J (u 1 , u 2 , u 3 ). The optimal control expressions are as follows u * 1 = max min 1, Proof. From the extremum condition of the control equation, we have ∂H Then, we can obtain the solution of the equations as follows Therefore, the optimal control solution can be expressed as

The fitting of parameters
Since the duration of Zika virus infection is about one week, we assumed a natural recovery rate of 0.862 per week in humans. Then, based on the weekly infection data of the Brazilian Zika virus in the first 20 weeks of 2015 in Reference [25], we used the least square method to fit the remaining parameters on the DEDiscover software, and the fitting results are shown in Table 3 and Figure 4. The sum of squares of the fitting residuals is 0.116.

Optimal control results
In order to find out the optimal control, we used the forward-backward Runge-Kutta method [26] to solve the optimal solution. Select the balance weight coefficient as ω E H = 50, ω I H = 30, ω E M = 0.2, ω I M = 0.1, a 1 = 30, a 2 = 20, a 3 = 50.
The results showed that the number of exposed people, infected people, exposed mosquitoes and infected mosquitoes were all on the rise before the control. Among them, the number of infected mosquitoes even reached a peak of 30,000 before the control, and the number of infected people was as high as 15,000 before the control. However, the number of infected and exposed mosquitoes and people decreased significantly under the control measures u 1 , u 2 , and u 3 . After control, the number of infected and exposed people tended to zero (see Figure 5(b),(d)) and the number of infected mosquitoes stabilized below 20 (see Figure 5(f),(h)). The spread of Zika virus was well controlled. Figure 6 shows the time-varying control profile.
(a) Exposed people before control (b) Exposed people after control (c) Infectious people before control (d) Infectious people after control (e) Exposed mosquito before control (f) Exposed mosquito after control (g) Infectious mosquito before control (h) Infectious mosquito after control Figure 5. Figure 2(a)-(h) shows the comparison of exposed people, infected people, exposed mosquitoes and infected mosquitoes before and after control.

Cost-effectiveness analysis
Next, we analyze the cost-effectiveness of investments in Zika virus prevention and control. Assume that the percentage of testing is the number of existing infections as a percentage of the immigrating population. The total cases averted (TCA) though control and the total costs (TC) associated with intervention are given by the following formula screening is estimated to be 10.3 dollars per person per week based on online quotes for commercially available Zika virus test reagents. Based on the data simulation results, we can eliminate the virus by using u 1 and u 3 . If only u 3 works, the number of patients will reduce, but the virus will not disappear.
Here, incremental cost-benefit ratio (ICER) is used to analyze the cost-effectiveness results of the three schemes. The formula is as follows .
ICER represents the incremental cost per unit of incremental health outcome, which is a classic method to analyze cost-effectiveness [26]. It can be seen from Table 4 that the ICER value of control option 2 is 326.0178819117, which is significantly higher than that of options 1 and 3. Therefore, control option 2 requires higher cost and lower efficiency, so option 2 is excluded from the alternative plan. Comparing options 1 and 3, you will find the ICER value of option 1 is higher than that of option 3. Therefore, it can be concluded that option 3 has the lowest ICER value, so it is the lowest cost and the highest efficiency option. From Figure 7, it is clear that option 3 can achieve the effect of option 1. From the perspective of environmental sustainability, option 3 is less harmful to the environment and human beings.

Conclusions
This paper focuses on the optimal control and cost-effective analysis of the Zika virus model with migration. Non-negativity and boundedness of the model are also shown. When no infected population enters, a disease-free equilibrium point exists in the system, and the disease-free equilibrium point was globally asymptotically stable. When cases migrate, there is no disease-free equilibrium in the model, only an epidemic equilibrium with global asymptotic stability. Based on this model, we propose three control measures u 1 , which means reducing mosquito bites through mosquito nets, and mosquito repellants. The factor u 2 stands for reducing the likelihood of infection and exposure among immigrants through testing and screening immigrants. The parameter u 3 represents the reduction of the mosquito population by spraying insecticides to prevent the spread of the Zika virus. The necessary conditions for the existence of an optimal solution are given using Pontriagin's maximum principle. Based on Control option 1, the cost-effectiveness of three control strategies was compared and analyzed. According to the incremental cost-effectiveness ratio results (see Table 3), the ICER value of Control option 2 is 326.0178819117, which is much higher than options 1 and 3. This means that Control option 2 is the least efficient and most expensive. On the other hand, Control option 3 with the lowest ICER value of -64.3479200904 exhibits the lowest cost and the highest return. Especially from the perspective of human health and environmental protection, the combination of these three options is capable of reducing our dependence on pesticides. Indeed, although we are concerned with screening, if immunization against the Zika virus is possible and the immigrant population is fully vaccinated and protected, this factor could also be reflected in the control item. Therefore, both vaccination and screening are indispensable control measures to control the Zika virus.
In Zika virus transmission, the environment also has an impact on the transmission of the virus, such as temperature, humidity, etc. Subsequently, there will be corresponding measures adjusted in virus control, which will be considered in the next study.

Use of AI tools declaration
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.