Phase III pilot trial for population suppression.

1 Introduction

The present pilot trial for population suppression of Anopheles arabiensis is structured over 3 distinct areas, each divided into 3 sectors.

Sector number 3, in the west bank area was selected for the sit treatment, were releases of sterile males were performed. The other 8 sectors were used for control.

A total of 94 sampling units of 1-hectare in size and each belonging to one of four land-class-land-use (LCLU) classes were selected and distributed across sectors. 3 adult traps were installed in each sampling unit.

A series of larvae surveys have been conducted in breeding sites within sectors.

A series of swarm surveys have been conducted in the sit sector following some of the releases of sterile males.

See the descriptive report for more details on the experimental set up and a description of the collected data.

A series of 18 ground releases of about 3k - 20k marked sterile males were performed from 11 points spread over the sit sector (Table 1.1, Figure 1.1).

Table 1.1: Release events of marked sterile males in the sit sector.
Release number Release date Number of sterile males released
1 2014-05-11 12300
2 2014-05-15 9700
3 2014-05-27 8400
4 2014-08-06 8100
5 2014-08-19 14744
6 2015-01-16 11505
7 2015-02-13 7450
8 2015-03-27 11480
9 2015-05-22 16600
10 2016-02-11 10000
11 2016-02-26 6900
12 2016-04-28 9000
13 2016-07-21 13700
14 2016-10-16 15000
15 2016-12-16 2900
16 2016-12-30 20000
17 2017-01-23 14000
18 2017-02-15 3000
Number of individuals released by date.

Figure 1.1: Number of individuals released by date.

The purpose of the present analysis is to quantify the suppression rate of the mosquito density in the field resulting from the releases of sterile males, and to elaborate on the effectiveness of the SIT as a strategy to control the malaria vector Anopheles arabiensis.

1.1 Sterile pressure

The sterile pressure is a calculated variable that estimates roughly the number of sterile individuals alive in the population at a given date in a logarithmic scale, based on a hypothesised daily survival probability. It takes into account the number of sterile individuals released prior to the target date, and the time elapsed since each release. Let \(\pi\) be the daily survival probability, the number of sterile individuals alive at time \(t\) is: \[\begin{equation} z_\pi(t) = \sum_{i:\; t(i) < t} R_i \exp(t \, \log\pi) \end{equation}\] where \(R_i\) are the released number of sterile individuals in release \(i\), for which the release times \(t(i)\) are prior to the target time.

The sterile pressure \(P(t) = \log(z_\pi(t))\) is then calculated at each observation time. We have used \(\pi = 0.9\), which yields the temporal estimates of sterile population and pressure shown in Figure 1.2. Other values have been tested, without improved association.

Sterile population and corresponding pressure over time estimated as the log-number of sterile individuals in the population, assuming a daily survival probability of 0.9.

Figure 1.2: Sterile population and corresponding pressure over time estimated as the log-number of sterile individuals in the population, assuming a daily survival probability of 0.9.

2 Exploratory analysis of the field survey data

Given the experimental design, where a series of releases of sterile males have been performed in a sit sector, we can appreciate the potential impact of the intervention in terms of:

  1. Difference in the population density in the sit sector, with respect to the other sectors.

    This requires the distinction between the variation that would be naturally expected between sectors and the excess variation that can be attributed to the impact of the intervention.

  2. Variation in the population density within the sit sector, as a function of the time since the last release of sterile males.

    Indeed, we can expect the population density to progressively drop following a release down to a minimum value from which it slowly recovers up to normal values after some time, if left alone.

Furthermore, variations in the population densities should reflect to different extents in adult captures, larvae surveys and swarmings.

In the present section we explore these two effects on the three types of surveys that have been conducted during the experiment.

2.1 Adult surveys

Figure 2.1 shows the number of adult Anopheles arabiensis captured in each trap, at each collection day, by sex and sector. Most often, the traps were empty, since they were actually resting sites and don’t accumulate captures. Only the 3% of non-zero values are represented in the figure.

Number of individuals catched in adult traps, by survey date, sex, trap and sector, in a logarithmic scale. Only non-zero catches are represented, as the vast majority of surveys are 0. Sterile males were released in sector 3, highlighted in orange, at the intervention times represented with light vertical lines.

Figure 2.1: Number of individuals catched in adult traps, by survey date, sex, trap and sector, in a logarithmic scale. Only non-zero catches are represented, as the vast majority of surveys are 0. Sterile males were released in sector 3, highlighted in orange, at the intervention times represented with light vertical lines.

Still, most of the raw outcomes are clustered on the lower-end of the scale, with only a few extreme observations that are actually visible. Thus, we need to summarise these data from various angles in order to describe the patterns properly.

Consider the capture rates (CR), measured as the average number of adult Anopheles captured per trap in a single survey. Figure 2.2 shows the capture rate across traps and sexes by sector and date.

There seems to be a initial drop in the average catches in the first months of the study (late 2014) at the sit sector, reaching almost zero by 2015 and staying at a low level with some sporadic peaks in 2016 and 2017.

However, the individual variability is very large, with a vast majority of zeroes and a few larger values for some specific surveys. This produces very noisy averages with large standard deviations when large values occur.

Furthermore, the peaks and drops are not obviously associated with the release times.

Capture rate by survey date and sector, across sexes and traps. The sit sector is highlighted in orange with a band of ±1 SD. Vertical lines represent release times.

Figure 2.2: Capture rate by survey date and sector, across sexes and traps. The sit sector is highlighted in orange with a band of ±1 SD. Vertical lines represent release times.

Figure 2.3 summarises further the data, aggregating the observations by year and displaying capture rates by year and sector.

Here again, the initial decline in the capture rate at the sit sector is apparent. Yet, the variations are still very high, and consistent with natural variation in the control groups.

Capture rate by year and sector, across sexes, traps and surveys. The sit sector is highlighted in orange with a band of ±1 SD. Vertical lines represent release times. Numbers at the right-hand side label the control sectors with highest capture rates.

Figure 2.3: Capture rate by year and sector, across sexes, traps and surveys. The sit sector is highlighted in orange with a band of ±1 SD. Vertical lines represent release times. Numbers at the right-hand side label the control sectors with highest capture rates.

Proportion of zeroes in catches by sector (points) and sex, across traps and time. Proportions in the sit sector are highlighted in orange.

Figure 2.4: Proportion of zeroes in catches by sector (points) and sex, across traps and time. Proportions in the sit sector are highlighted in orange.

Empirical distribution of the number non-zero catches by sector and sex, across traps and time, in a logarithmic scale. Distributions from the sit sector are highlighted in orange.

Figure 2.5: Empirical distribution of the number non-zero catches by sector and sex, across traps and time, in a logarithmic scale. Distributions from the sit sector are highlighted in orange.

From figures 2.4 and 2.5, there does not seem to be evidence of any decrease in the abundance of the wild population of Anopheles arabiensis in the sit sector, beyond the natural variation across sectors.

However, the impact of the releases of sterile males could be limited in time, and thus not noticeable when outcomes are averaged across time.

We explore this question next, by looking at results as a function of time since the last release.

Proportion of zeroes in catches by sector, and time since last release across traps and sexes. Proportions in the sit sector are highlighted in orange.

Figure 2.6: Proportion of zeroes in catches by sector, and time since last release across traps and sexes. Proportions in the sit sector are highlighted in orange.

There might be a trend, but the pattern is also consistent with random variation.

Distribution of non-zero catches by role and time since last release, across traps, control sectors and sexes. Distributions from the sit sector are highlighted in orange.

Figure 2.7: Distribution of non-zero catches by role and time since last release, across traps, control sectors and sexes. Distributions from the sit sector are highlighted in orange.