#### Abstract

We proposed in this study a deterministic mathematical model of malaria transmission with climate variation factor. In the first place, fundamental properties of the model, such as positivity of solution and boundedness of the biological feasibility of the model, were proved whenever all initial data of the states were nonnegative. The next-generation matrix method is used to compute a basic reproduction number with respect to the disease-free equilibrium point. The Jacobian matrix and the Lyapunov function are used to check the local and global stability of disease-free equilibriums. If the basic reproduction number is less than one, the model’s disease-free equilibrium points are both locally and globally asymptotically stable; otherwise, an endemic equilibrium occurs. The results of the sensitivity analysis of the basic reproduction numbers were obtained, and its biological interpretation was provided. The existence of bifurcation was discussed, and the model exhibits forward and backward bifurcations with respect to the first and second basic reproduction numbers, respectively. Secondly, using the maximum principle of Pontryagin, the optimal malaria reduction strategies are described with three control measures, namely, treated bed nets, infected human treatment, and indoor residual spraying. Finally, based on numerical simulations of the optimality system, the combination of treatment and indoor spraying is the most efficient and least expensive strategy for malaria eradication.

#### 1. Introduction

Malaria is a serious vector-borne disease caused by a parasite called *Plasmodium* and transmitted between humans (hosts) and mosquitoes (vectors) through bites from infected female *Anopheles* mosquitoes [1]. It can spread through blood transfusions, needle sharing, or congenital causes. The parasites then reproduce in the human liver and bloodstream. Malaria is often associated with poverty, which has a negative impact on the country’s economic development. According to the most recent global malaria report, published in December 2021, an estimated 241 million cases and 627000 deaths were reported worldwide, with the African region accounting for 94% of all cases in 2020 [2]. The use of insecticide-treated bed nets, antimalaria drug treatment, and indoor residual spraying are the most effective methods of preventing malaria transmission dynamics [3]. Most mosquito distributions worldwide have been influenced by climatic factors such as temperature, rainfall, and malaria transmission increases with temperatures ranging from to , which create conducive conditions for mosquito breeding rates [4].

Ross [5] initiated a mathematical model on the dynamics of malaria spread. He developed SISSI models for human and malaria populations. According to Ross, reducing the number of mosquitoes to below a certain threshold is sufficient for malaria control. Following the Ross model, several models were proposed by scholars who extended the Ross model by considering various factors such as the exposure time in humans and mosquitoes [6–8] and climate variability on mosquito death rate, contact rate, and birth rate [9–11].

A group of researchers studied an optimal malaria transmission control model to determine the role of control measures in disease transmission. For example, Okosun et al. [12] presented the malaria transmission dynamics with optimal control and cost-effectiveness analysis. The authors used a combination of three control measures: treated bed nets, antimalarial drug treatment of infected humans, and indoor residual spraying. They conclude that the combination of drug treatment of infected humans and indoor insecticide spraying is the most effective and least expensive way to prevent malaria spread. Makinde et al. [13] proposed a malaria model for disease spread with optimal control and cost-effectiveness analysis based on three control measures. Using the incremental cost-effectiveness ratio methods, a cost-effectiveness analysis was performed to support the results of the optimal control problem. The authors used numerical simulations to supplement the theoretical results. Finally, they proposed that a combination of treated bed nets, drug treatment, and indoor residual spraying is the most effective and less expensive intervention strategy. Okosun et al. [14] presented the SEIRS-SEI malaria transmission model with optimal control and cost-effectiveness analysis using combinations of three control measures, including treated bed nets, infected human treatment, and indoor residual spraying. The authors presented the use of infected human treatment with drugs and indoor residual spraying as the most cost-effective controls to prevent malaria spread. Otieno et al. [15] presented a malaria transmission model with an optimal control problem, employing four time-dependent control measures. Pontryagin’s minimum principle was described in order to derive the necessary conditions for optimal malaria control. Their findings concluded that a combination of treated bed nets, treatment, and spraying is the most effective strategy for mitigating disease. Temesgen et al. [16] presented SIRS for the human population and SI malaria model for malaria transmission by incorporating three continuous controls. The authors suggested that combinations of treatment of infected humans and indoor residual spraying are the best strategies for the prevention of the disease.

However, none of these models take into account the impact of climate variability on malaria epidemics with optimal control and cost-effectiveness analysis with a logistic growth of climate variation with respect to mosquito breeding and malaria infection. In this paper, the malaria transmission model [17] is extended to the optimal control problem in the presence of climate variability in mosquito breeding rate and malaria infection, and further analysis on the least cost strategy is considered using the concept of cost-effectiveness analysis.

This paper is arranged as follows: In Section 2, we propose a malaria transmission model that shows the role of climate variability in malaria epidemics. In Section 3, the model analysis of the model is illustrated. In Section 4, we describe the analysis of sensitivity. In Section 5, the optimal control of the malaria transmission model is analytically analyzed using the maximum principle of Pontryagin. In Section 6, we show the simulation of theoretical results. In Section 7, the analysis of cost-effectiveness is depicted. The conclusion of the work is discussed in Section 8.

#### 2. Model Description and Formulation

In this section, we developed a malaria transmission dynamics model in which the human population is represented by the SIRS model and the population of mosquitoes is described by the SI model due to their short lifespan. The total population of human at time denoted by is divided into three subpopulations based on their disease status; susceptible humans *S _{h}(t)* are individuals who are vulnerable to risk of developing an infection from the mosquito. Infected humans

*I*are individuals who show the symptom of the disease and can transmit the disease to mosquito. Recovered humans

_{h}(t)*R*are individuals who got temporary immunity, so that the total human populations is given by . Individuals who are born or migrate are assumed to be susceptible human populations and considered as the recruitment rate denoted by . Susceptible human populations become infected if they come into contact with infected mosquitoes at the rate which depends on climate change, represents the human contact rate with vector when there is no climate variation, and denotes the human incremental contact rate due to climate variation. Humans leave the total population with a death rate , and due to malaria disease, the induced death rate reduces the population of human. Infected human recovered through treatment rate using antimalarial drugs. The recovered human population whose immunity is not permanent can lose it and become susceptible to reinfection at the rate . The total vector population given by at time is subgrouped into susceptible mosquito and infected mosquito . Hence, the total vector population is given by . The vector population is recruited at the rate which depends on climate change, whereas denotes the vector birth rate when there is no climate variation and is the vector incremental birth rate due to variation of climate. A mosquito population gets an infection when it contacts with infected human at the rate which depends on climate variation, represents the vector contact rate with a human when there is no climate change, and denotes incremental contact rate due to variations in climate. The vector natural death rate is . We assume that mosquitoes do not recover from malaria disease. Thus, mosquitoes do not die due to disease infection. The growth rate of the temperature follows a logistic function, is the temperature-dependent rate of precipitation, and represents the maximum temperature for the mosquito to be most active, whereas the minimum temperature for the mosquito to be less active is denoted by . Also, we have considered that the parameters stated in Table 1 are positive. Figure 1 illustrates the transmission dynamics of the malaria diagram.

_{h}(t)Depending on the flowchart diagram in Figure 1, we write a model that governs the transmission dynamics of malaria disease using a system of ordinary differential equations.wherewith initial conditions

Based on the fact that the average temperature at the Earth’s surface rises, more evaporation occurs, which in turn increases overall precipitation (rainfall). Therefore, a warming climate is expected to increase precipitation (rainfall) in many areas. This implies that rainfall pattern is an increasing function of temperature. Note that is the temperature-dependent rate of precipitation and is the amount of precipitation at . In essence,

#### 3. Mathematical Analysis of the Model

##### 3.1. Invariant Region

The invariant region is used to determine where the model’s solution is constrained. The proposed model (1) has two kinds of populations. Firstly, the total population of humans is given by and differentiating with respect to time. Then, by adding the first three equations from system (1), we obtain

This implies that equation (5) becomes

By solving equation (6), we obtain that

Thus, the invariant region of system (1) for the human population is a given by

Secondly, the total number of the mosquito populations from system (1) is given as

Hence, the differentiating with respect to time is given by equation

By solving equation (10), we obtain . Hence, the invariant region of system (1) for the mosquito population is given by

Consequently, the dynamics of system (1) is studied in the biological feasible region of the formwhich is a positive invariant set under the flow induced by the solution set of system (1).

##### 3.2. Positivity of the Solution

For model (1), we will show that all solutions of system with positive initial data will remain positive for all times .

Theorem 1. *If and are nonnegative, then the solutions and of system (1) are nonnegative for .*

*Proof. *The first equation from system (1) is given byIntegrating equation (13) with respect to time and using the method of separation variables with applying the initial conditions, we obtainClearly, equation (13) is nonnegative at all time . The other state variables , and are nonnegative for all time by the same procedure. As a result, the malaria transmission model stated in system (1) is both epidemiologically significant and mathematically well-posed in equation (12).

##### 3.3. Basic Reproduction Numbers

In the first place, we calculate the disease-free equilibrium (DFE) points of system (1) before calculating an expression for basic reproduction numbers. To do so, we equate all equations of model (1) to zero with and . In this case, we obtain two malaria-free equilibrium points of model (1), which are given by or where

One advantage of determining the DFE of system (1) is to help us to calculate an expression for the basic reproduction number, which is defined as the average amount of secondary infections caused by a primary infectious in the given period [18]. It can be obtained using the approach of the next-generation matrix and that is the dominant eigenvalue of the next-generation matrix. For model (1), to obtain and we rewrite model (1), beginning with newly infective classes of humans and mosquitoes, which are given by

The right-hand side (RHS) of equation (16) can be written in the form of , where

The partial derivatives of and in equation (17) evaluated at DFE of system (1) produce two matrices and , respectively,

As a result, the basic reproduction number , where is the dominant eigenvalue of . Thus, the first basic reproduction number calculated at the first disease-free equilibrium points is given by

Similarly, the second basic reproduction number calculated at the second disease-free equilibrium points is given bywhere , and .

Generally, the basic reproduction number at maximum temperature and rainfall variation in terms of is given bywhere is the first basic reproduction number at for the mosquito to be less active to breed.

##### 3.4. Local Stability of Disease-Free Equilibrium

Theorem 2. *The disease-free equilibrium point of system (1) is locally asymptotically stable in if .*

*Proof. *We begin with finding that the Jacobian matrix of system (1) calculated at DFE is given byFrom equation (22), we can obtain the characteristic equation of the following form:whereHence, from equation (23), we haveAnother equation from equation (23) will be a second-degree polynomialBy using the Routh–Hurwitz criteria [19], equation (26) has a negative real root if and . Hence, we can observe that since it is the sum of nonnegative parameters and the value of at is given byIt is worth noting from equation (27) that if . Since , so the DFE is locally asymptotically stable if .

##### 3.5. Global Stability of Disease-Free Equilibrium

Theorem 3. *If , , then the disease-free equilibrium point(s) of system (1) is globally asymptotically stable in .*

*Proof. *To discuss the global behavior of DFE of system (1), we use a technique implemented by the Lyapunov theorem [20] by first constructing a suitable Lyapunov functionBy differentiating equation (28) with respect to time (*t*), the result obtained becomesConsequently, we observe that if and if and only if . So, the dominant bounded invariant set is the singleton which is DFE in . Therefore, by the well-known LaSalle’s invariant principle [21], every solution that begins in the domain bounded invariant set approaches DFE as time tends to infinity when . Therefore, the DFE is globally asymptotically stable in if .

##### 3.6. Malaria Present Equilibrium

Malaria endemic equilibrium point is a situation where malaria exists in the human population regardless of time and other factors. The endemic equilibrium point can be obtained by setting RHS of all model equations of system (1) equal to zero. From model (1), the malaria present equilibrium point at () = () is given by

From equation (30), the endemic equilibrium easily satisfies the following polynomial and is computed from the equation: where

Hence, and whenever Solving for in equation (31), we obtain . Thus, the model has no nonnegative malaria present equilibrium whenever . This guarantees that backward bifurcation does not exist in the model if .

Similarly, the endemic equilibrium point at () = (), and solving for as expression of parameters, we obtainwhere , and .

From equation (33), the endemic equilibrium satisfies the following polynomial and is computed from the equationwherewhere , and .

From equation (35), we see that if , then it immediately implies that and DFE exists for both and . However, does not immediately imply as the value of will be larger than unity that shows that whereas present DFE. Thus, may become an endemic situation or exhibit backward bifurcation while undergoes forward bifurcation only.

#### 4. Sensitivity Analysis

Malaria control and eradication strategies should focus on key parameters that have a significant impact on the basic reproduction number. The purpose of carrying out the sensitivity analysis of the basic model parameter is to identify a parameter that affects the basic reproduction number. Essentially, the robustness of system predictions to parameter values can be expressed using the normalized sensitivity index of the basic reproduction number to the given basic parameters because the values of those parameters can increase or decrease a basic reproduction number and vice versa. This method relies on our knowledge of the parameters that have a large influence on the basic reproduction number () in order to design the best disease control strategies. We used the technique described in [13, 22, 23] to perform the sensitivity analysis.

*Definition 1. *(See [13, 22, 24, 25]): The forward sensitivity index of , which is differentiable with respect to a given basic parameter , is defined asThe sensitivity index of of model (1) with respect to parameter , for instance, is obtained asUsing the same approach with respect to the rest of the parameters, are computed, and the sensitivity indices are presented in Table 2.

By the same procedure, the sensitivity index of of the model (1) with respect to is given asSimilarly, with respect to other basic parameters, are computed, and the sensitivity indices are described in Table 3.

##### 4.1. Interpretation of the Sensitivity Indices

We described sensitivity indices of basic reproduction number () with respect to eight basic parameters in Table 2. The output showed that the parameters with a positive sensitivity index increased the value of *R*_{01} as their values increased, while the rest of the parameters remained constant. And if the values of the parameters with negative indices are increased while the values of the other parameters remain constant, the value of decreases. Similarly, the sensitivity indices of with respect to eleven basic parameters are shown in Table 3. The parameters with positive sensitivity indices have a high impact on malaria transmission in the community as their values increase. The basic parameters with negative sensitivity indices increase the malaria disease if their values decrease while the other parameters remain constant. For example, shows that decreasing (increasing) the rate of the human recruitment by 10% decreases (increases) the basic reproduction number by 5%; similarly, indicates that decreasing (increasing) the rate of the mosquito death by 10% increases (decreases) the basic reproduction number by 10%.

#### 5. Extension of the Model into Optimal Control

According to the necessity and severity of a disease, many intervention strategies are applied by public health officials to control the disease [25]. In this study, we extended the malaria transmission model (1) to an optimal control problem to determine control strategy decisions involving a mathematical model of biological situations [30]. Using this method, we hope to demonstrate the best malaria prevention strategy. After incorporating the controls into the malaria transmission model (1), the obtained state equations are as follows: where .

The control functions represent that the use of a treated bed net to protect against mosquitoes is , the treatment of infected humans with antimalarial drugs is , and the use of indoor residual spraying to kill mosquitoes is . The objective functional that we developed for the optimal control model (39) is given aswhere denoted the terminal time, and are weight constants for the infected human and mosquito, respectively, and , , and are weight constants for each control, respectively. The expression represents the cost function that corresponds to the controls and is quadratic in the other pieces of literature [18, 31–34, 36–38]. The goal of the objective functional (40) is to reduce the total number of infected humans , infected mosquitoes , and control costs . The primary goal is to compute a triple optimal control , and such thatwhere : such that , and are Lebesgue measurable on with is the control set. The obtained Hamiltonian (H) function of the optimal control problem that consists of equations (39) and (40) is represented as

From equation (42), the minimize Hamiltonian function with respect to is given bywhere , and are adjoint variables. Next, to obtain the costate variables by using Pontryagin’s maximum principle [32], with the existence result of [35], the following theorem is stated.

Theorem 4. *Given optimal controls and a solution of the corresponding state system that minimize J() over subject to equation (39), then adjoint variables , and hold the adjoint system**With transversality conditions,**Furthermore, the optimal controls are represented by*

*Proof. *To obtain the form of the costate equations, we compute the derivative of the Hamiltonian function (H) equation (42) with respect to , and , respectively. Then, the adjoint or costate equation obtained is given byWith transversality conditions,To obtain the control value, we compute the partial derivative of Hamiltonian given byObviously, after a partial derivative of Hamiltonian with respect to the controls, the result becomesMoreover, solving for controls variables from equation (50), we obtainRearranging the solution of (51) with the boundary condition of each control, we gotNext, we will see the simulation of the optimality system to identify an optimal strategy that is most optimal to minimize the spread of malaria transmission.

#### 6. Numerical Simulation

In this study, we solved the optimality system and used the forward-backward sweep to solve the state and adjoint systems in order to obtain the optimal strategy. We used the forward fourth-order Runge–Kutta to solve the state equations (39) due to the initial value of state variables. Similarly, the adjoint equations are solved using backward fourth-order Runge–Kutta due to the transversality condition (45) having the solution of state functions and the value of optimal controls. The controls are then updated using a convex combination of the previous controls and the value from the optimality conditions (46). This situation will continue until two consecutive iterations are close enough to each other [30]. The initial conditions that we used for numerical simulation of the optimality system are *T*(0) = 16°C, and *R*(0) = 11 mm, as well as the parameter values from Table 4. We used the following weight constant values for the states and controls: , and . We used the following four strategies with different combinations of two controls at a time and three controls at a time to determine the impact of each control on malaria reduction.

##### 6.1. Strategy A: Combination of Use of Treated Bed Nets () and Treatment of Infected Humans ()

We optimized the objective function (39) using the treated bed net and the treatment of infected humans with antimalarial drugs when the value of the indoor residual spraying was set to zero. Also, as shown in Figure 2(a), if there are controls, the number of infected humans decreases and then tends to zero, whereas the number of infected humans increases if there are no controls. Similarly, in Figure 2(b), we can see that the infected mosquito decreases as the control strategy is used, whereas the infected mosquito increases in the uncontrolled case. The control profiles shown in Figure 2(c) with this strategy suggest that controls on treated bed net kept maximum level (100%) until the end of the implementation while infected human treatment using antimalarial drugs retained its maximum bound for 6 days then gradually tended to a minimum after the 82nd day.

**(a)**

**(b)**

**(c)**

##### 6.2. Strategy B: Combination of Use of Treated Bed Net () and Indoor Spraying ()

This strategy combined the treated bed net as a control and indoor residual spraying as a control to reduce total infected populations and costs associated with the treatment of infected humans . We can see in Figure 3(a) that the amount of infected human with controls decreases and reaches its lowest point. In contrast, when there are no controls, the number of infected humans increases to a certain point. According to Figure 3(b), an infected mosquito decreases in the control strategy and decreases. To its minimum point, an infected mosquito increases in the uncontrolled case. According to the control profiles depicted in Figure 3(c), control on treated bed net is maintained at its maximum(100%) for 171 days, while indoor residual spraying is maintained at its maximum for 8 days before dropping to its lowest value after the 120th day.

**(a)**

**(b)**

**(c)**

##### 6.3. Strategy C: Combination of Treatment of Infected Humans () and Indoor Spraying ()

In this study, we used a combination treatment of infected humans as control and indoor residual spraying as control to reduce total infected populations and save money. Figure 4(a) shows that when using this control strategy, the number of infected humans decreases to its lowest value and becomes smaller when controls are used. In Figure 4(b), the infected mosquito decreases in the occurrence of control strategy and then drops to its minimum values, whereas in the uncontrolled case, an infected mosquito increase is observed. With this approach, the control profiles in Figure 4(c) conclude that controls on the treatment of infected human maintain upper level (100 percent) for 9 days, while the insecticides spray maintains maximum level for 7 days and then declines to lower bound at the end of the 80th day.

**(a)**

**(b)**

**(c)**

##### 6.4. Strategy D: Use of Treated Bed Net (), Treatment (), and Indoor Spraying ()

We used a combination of three continuous controls on the treated bed net ; treatment of infected human ; and indoor residual spraying to reduce the objective function (39). By implementing those three control strategies, we discovered that, in Figure 5(a), the amount of infected human decreases and decreases to a minimum level, whereas the amount of infected human increases to a certain point if no control is used. In Figure 5(b), the infected population of mosquitoes is decreasing and minimized to its lower bound when a control strategy is used, whereas in the absence of a control strategy, an increase in infected mosquitoes is observed. Using this strategy, the control profile in Figure 5(c) depicted that controls on treated bed net and treatment were kept at 100% coverage for 5 and 8 days, respectively. Then, treated bed net and treatment tend to their lowest bound after 80 days, whereas indoor spraying retains its highest values (100 percent) for 162 days before dropping to the lowest level on the 100th day.

**(a)**

**(b)**

**(c)**

Figure 6 shows the existence of bifurcation for malaria transmission model. Bifurcation occurs. This implies that, in equation (13), if , this automatically implies that and DFE exists for both and . However, do not automatically describe , , because the value of may be larger than unity that shows that DFE () may depict backward bifurcation while DFE () only depicts forward bifurcation.

#### 7. Cost-Effectiveness Analysis

In this study, we must determine the most effective and least expensive method of reducing disease transmission. We employed the incremental cost-effectiveness ratio method to arrive at this strategy (ICER). Using this technique, we incrementally compare more than one competing intervention; for example, one intervention could be compared to a second, less effective alternative. This approach was defined as the ratio of the difference in averted costs between two strategies to the difference in the total number of infections saved [23]. We calculated the total cost avoided and total infections saved from the numerical simulation of the optimal control problem, and the control strategy is ordered in increasing order based on the total infections saved, as shown in Table 5. The amount of total infection saved is calculated by subtracting the total number of humans infected with malaria using control from the total number of humans infected with malaria without control, while the cost averted of each strategy was obtained by using the cost function represented by , , and over the time [17, 28]. The amount of the total infection saved and the total cost of all strategies with their ICER is given in Table 6.

The value of the cost-effectiveness ratio (ICER) is computed from the total number of saved populations and total cost of averted for each strategy given in Table 5, which is used to compare the differences between the two strategies as obtained and given by

With the above results, the number of infections saved with ICER for four different strategies is shown in Table 6.

Table 6 compares the interventions B and C. ICER(C) is less than ICER(B) as shown in the table.

It implies that strategy B is costly and has a low likelihood of saving people. As a result, C saves more people than B. Then, from the competing strategies, we have removed B. Then, as shown in Table 7, calculate the ICER for the remaining strategies C, A, and D.

The ICER(A) is greater than the ICER(C) from the competing intervention strategies listed in Table 7. This demonstrates that ICER(C) strategy outperforms ICER(A). As a result, ICER(A) is less effective and more expensive than ICER(C). As a result, we removed strategy A from the group of competing strategies and calculated the ICER as shown in Table 8.

Table 8 with intervention strategies C and D shows that ICER(C) is less than ICER(D). This indicates that strategy C strongly outperforms strategy D. As a result, the C strategy has the lowest total cost and is the most optimal. Based on the findings of the analysis, we believe that intervention C, which consists of treating infected humans and spraying indoors, is the best optimal and least expensive strategy for limiting disease spread.

#### 8. Conclusion

The impact of climate variability on malaria epidemics is described in this paper using deterministic mathematical modeling. The model analysis revealed qualitatively that the model’s solution is both bounded and positive in the fixed domain. The method of the next-generation matrix is used to calculate the basic reproduction number with respect to the disease-free equilibrium. The Jacobian matrix and the Lyapunov method are used to demonstrate the local stability and global stability of the disease-free equilibrium. Furthermore, if the basic reproduction number is less than one, the disease-free equilibrium is locally and globally asymptotically stable; otherwise, a positive endemic equilibrium occurs. The model’s sensitivity has been described in detail, and the model demonstrates forward and backward bifurcation. According to the analytical findings, the best optimal way to prevent malaria epidemics is to reduce human-mosquito contact, increase mosquito death rates, and increase the treated rate of an infected human. Furthermore, we extended the model to an optimal control problem with three continuous controls, such as personal prevention with a treated bed net, drug treatment of infected humans, and indoor residual spraying for mosquito minimizing strategy. The maximum principle of Pontryagin is applied to obtain the necessary conditions for optimal control, and the cost-effectiveness analysis is described for all combinations of the controls considered in the study. Based on the numerical analysis, we propose that the combination of treatment and indoor residual spraying is the best strategy for effectively reducing malaria.

#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

#### Conflicts of Interest

The authors have no particular conflicts of interest within the manuscript.