Abstract:This work is devoted to the mathematical study of an optimization problem regarding control strategies of mosquito population in a heterogeneous environment. Mosquitoes are well known to be vectors of diseases, but, in some cases, they have a reduced vector capacity when carrying the endosymbiotic bacterium Wolbachia. We consider a mathematical model of a replacement strategy, consisting in rearing and releasing Wolbachia-infected mosquitoes to replace the wild population. We investigate the question of optimizing the release protocol to have the most effective replacement when the environment is heterogeneous. In other words we focus on the question: where to release, given an inhomogeneous environment, in order to maximize the replacement across the domain. To do so, we consider a simple scalar model in which we assume that the carrying capacity is space dependent. Then, we investigate the existence of an optimal release profile and prove some interesting properties. In particular, neglecting the mobility of mosquitoes and under some assumptions on the biological parameters, we characterize the optimal releasing strategy for a short time horizon, and provide a way to reduce to a one-dimensional optimization problem the case of a long time horizon. Our theoretical results are illustrated with several numerical simulations.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to optimize the mosquito population replacement strategy in a heterogeneous environment. Specifically, the researchers focus on how to most effectively replace the wild mosquito population by releasing Wolbachia - carrying mosquitoes in the case of non - uniform environments. The focus of the paper is to determine the optimal release strategy, that is, in a given non - uniform environment, where the Wolbachia - carrying mosquitoes should be released to maximize the population replacement effect in the entire area.
### Background and Motivation
Mosquitoes are vectors for many diseases, such as dengue fever, Zika virus and chikungunya. Using Wolbachia is an effective method to control the spread of these diseases because mosquitoes carrying this bacterium have a weaker ability to transmit diseases. In addition, Wolbachia has the property of being vertically transmitted to offspring and can produce cytoplasmic incompatibility (CI), which reduces the reproductive ability when male mosquitoes carrying Wolbachia mate with wild female mosquitoes. Based on these properties, researchers have developed a population replacement strategy, that is, to replace the wild mosquito population by breeding and releasing Wolbachia - carrying mosquitoes.
### Research Objectives
The main objective of the paper is to optimize the mosquito release strategy under the condition of considering environmental heterogeneity. Specifically, the researchers hope to answer the following questions:
- In a given non - uniform environment, how to determine the optimal release location to maximize the effect of population replacement?
- How to design the optimal release protocol, taking into account production limitations, to achieve the most effective population replacement?
### Mathematical Model
The researchers used a simplified mathematical model to describe the dynamics of the mosquito population. The model assumes that the carrying capacity of the environment is space - dependent and ignores the active movement of mosquitoes. The core equations of the model are as follows:
\[ \frac{\partial p}{\partial t} = f(p) + \frac{u(t, x)}{K(x)} g(p) - \frac{D \Delta K(x)}{K(x)} \psi(p) \]
where:
- \( p(t, x) \) represents the proportion of Wolbachia - carrying mosquitoes at time \( t \) and location \( x \).
- \( K(x) \) is the carrying capacity of the environment, which is a function of the spatial variable \( x \).
- \( u(t, x) \) is the release function, representing the number of Wolbachia - carrying mosquitoes released at time \( t \) and location \( x \).
- \( f(p) \) and \( g(p) \) are given functions, which respectively describe the dynamics of the mosquito population and the release effect.
- \( D \) is the diffusion coefficient, which is assumed to be zero to simplify the model.
### Optimal Control Problem
The researchers defined an optimal control problem, with the objective of minimizing the distance between the proportion of Wolbachia - carrying mosquitoes and the fully invaded state at a given final time \( T \). The specific objective function is:
\[ \min_{u_0 \in U_0, C, M} \int_{\Omega} K(x)^2 (1 - p(T, x))^2 \, dx \]
where:
- \( u_0(x) \) is the initial release function, satisfying \( 0 \leq u_0 \leq M \) and \( \int_{\Omega} u_0(x) \, dx \leq C \).
- \( K(x) \) is the carrying capacity of the environment.
- \( p(T, x) \) is the proportion of Wolbachia - carrying mosquitoes at time \( T \) and location \( x \).
### Main Results
The main results of the paper include:
1. **Short - time range**: When \( T \leq T_0 \), there is a unique optimal solution \( p^*_0 \), and it can be determined by solving a one - dimensional optimization problem.
2. **Long - time range**: When \( T > T_0 \), the optimal solution may be a unimodal function, that is, it first decreases and then increases. In this case, the optimal solution can be determined by solving a one - dimensional optimization problem.
These results provide...