2,795
Views
36
CrossRef citations to date
0
Altmetric
Original Articles

Modelling releases of sterile mosquitoes with different strategies

&
Pages 1-14 | Received 14 Apr 2014, Accepted 10 Oct 2014, Published online: 07 Nov 2014

Figures & data

Figure 1. With parameters a=20 and k=0.3 as given in Example 3.2, four different values b=1.3,2.5,5,7, are used. With the same initial value w=1, chaotic behaviour is exhibited in the upper left figure, and then stable 4-cycle, 2-cycle, and 1-cycle are presented in the upper right figure, the lower left and right figures, respectively.

Figure 1. With parameters a=20 and k=0.3 as given in Example 3.2, four different values b=1.3,2.5,5,7, are used. With the same initial value w=1, chaotic behaviour is exhibited in the upper left figure, and then stable 4-cycle, 2-cycle, and 1-cycle are presented in the upper right figure, the lower left and right figures, respectively.

Figure 2. This is a schematic diagram to show the existence of positive fixed points. The figures on the left and right correspond to k1 and k<1, respectively. The intersection between the curve of ((1+w)2/w)+b and the curve of a(1+w)ekw gives a positive fixed point of Equation (Equation19). Keep the curve of a(1+w)ekw fixed. As parameter b increases gradually from zero to exceeding bc, the curve of ((1+w)2/w)+b moves up gradually. Accordingly, there exist two, one, or no intersections, and hence Equation (Equation19) has two, one, or no positive fixed points.

Figure 2. This is a schematic diagram to show the existence of positive fixed points. The figures on the left and right correspond to k≥1 and k<1, respectively. The intersection between the curve of ((1+w)2/w)+b and the curve of a(1+w)e−kw gives a positive fixed point of Equation (Equation19(19) wn+1=awn1+wn+B(wn)wne−kwn=awn(1+wn)(1+wn)2+bwnwne−kwn.(19) ). Keep the curve of a(1+w)e−kw fixed. As parameter b increases gradually from zero to exceeding bc, the curve of ((1+w)2/w)+b moves up gradually. Accordingly, there exist two, one, or no intersections, and hence Equation (Equation19(19) wn+1=awn1+wn+B(wn)wne−kwn=awn(1+wn)(1+wn)2+bwnwne−kwn.(19) ) has two, one, or no positive fixed points.

Figure 3. The model equation for the upper figure (a) is (Equation8) with the constant releases. The lower two figures (b) on the left and (c) on the right are based on equation (Equation14) where the number of releases is proportional to the wild mosquito population size, and Equation (Equation19) where the number of releases is proportional to the wild mosquito population size plus the saturation, respectively. We fix the same parameters a=20 and k=0.3, and use b as the bifurcation parameter. For bs<b<bc, there exists a unique positive fixed point. As b decreases, backward period-doubling bifurcations occur for all of the model equations. Notice, in all of the figures, bs corresponds to the b point where there exists two fixed points on the left but only one on the right. The threshold value bc corresponds to the b value where there is a positive fixed on the left and no positive fixed point on the right.

Figure 3. The model equation for the upper figure (a) is (Equation8(8) wn+1=awnwn+bwne−kwn.(8) ) with the constant releases. The lower two figures (b) on the left and (c) on the right are based on equation (Equation14(14) wn+1=awn1+(1+b)wnwne−kwn.(14) ) where the number of releases is proportional to the wild mosquito population size, and Equation (Equation19(19) wn+1=awn1+wn+B(wn)wne−kwn=awn(1+wn)(1+wn)2+bwnwne−kwn.(19) ) where the number of releases is proportional to the wild mosquito population size plus the saturation, respectively. We fix the same parameters a=20 and k=0.3, and use b as the bifurcation parameter. For bs<b<bc, there exists a unique positive fixed point. As b decreases, backward period-doubling bifurcations occur for all of the model equations. Notice, in all of the figures, bs corresponds to the b point where there exists two fixed points on the left but only one on the right. The threshold value bc corresponds to the b value where there is a positive fixed on the left and no positive fixed point on the right.

Table 1. Summary table for the threshold release values.