# Estimates of the population size and dispersal range of Anopheles arabiensis in Northern KwaZulu-Natal, South Africa: implications for a planned pilot programme to release sterile male mosquitoes

Maria L. Kaiser, Oliver R. Wood, David Damiens, Basil D. Brooke, Lizette L.

Koekemoer, Givemore Munhenga

► To cite this version:

Maria L. Kaiser, Oliver R. Wood, David Damiens, Basil D. Brooke, Lizette L. Koekemoer, et al. Estimates of the population size and dispersal range of Anopheles arabiensis in Northern KwaZulu-Natal, South Africa: implications for a planned pilot programme to release sterile male mosquitoes. Parasites & Vectors, 2021, 14 (1), pp.205. 10.1186/s13071-021-04674-w. http://hal-04975211

HAL Id: hal-04975211

https://hal.science/hal-04975211v1

Submitted on 4 Mar 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
RESEARCH

Open Access

<figure><img src="image_2.png" /></figure>

Check for
updates

# Estimates of the population size and dispersal range of *Anopheles arabiensis* in Northern KwaZulu-Natal, South Africa: implications for a planned pilot programme to release sterile male mosquitoes

Maria L. Kaiser¹,²*, Oliver R. Wood¹,², David Damiens³,⁴, Basil D. Brooke¹,², Lizette L. Koekemoer¹,² and Givemore Munhenga¹,²

**Abstract**

**Background:** *Anopheles arabiensis* is a major malaria vector, recently implicated as contributing to ongoing residual malaria transmission in South Africa, which feeds and rests both indoors and outdoors. This species is, therefore, not effectively targeted using core malaria vector control interventions alone. Additionally, increasing resistance to available insecticides necessitates investigations into complementary non-insecticide-based vector control methods for outdoor-resting mosquitoes. The feasibility of the sterile insect technique (SIT) as a complementary vector control intervention is being investigated in South Africa. Successful implementation of an SIT programme largely depends on inundating a target insect population with sterilized laboratory-bred males. Therefore, knowledge of the native population size and dispersal ability of released sterile laboratory-reared males is critical. In this study, we estimated the male *An. arabiensis* population size and the dispersal of released males in an area targeted for a pilot sterile male release programme.

**Methods:** Three separate releases were performed within a 2-year period. Approximately 5000–15,000 laboratory-reared male *An. arabiensis* (KWAG) were produced and marked for mark-release-recapture experiments. To recapture released mosquitoes, cloth tubes were deployed in widening concentric circles. The average dispersal distance of released males was calculated and the wild male *An. arabiensis* population size was estimated using two Lincoln index formulae. The natural population was sampled concurrently and *Anopheles* species diversity examined.

**Results:** The *Anopheles gambiae* complex and *An. funestus* group species made up the majority of wild collections along with other anophelines. The *An. arabiensis* population size was estimated to be between 550 and 9500 males per hectare depending on time of year, weather conditions and method used. Average dispersal distance of marked males ranged from 58 to 86 m. Marked males were found in swarms with wild males, indicating that laboratory-reared males are able to locate and participate in mating swarms.

<footnote>*Correspondence: mariak@nicd.ac.za
¹Centre for Emerging Zoonotic & Parasitic Diseases, National Institute for Communicable Diseases of the National Health Laboratory Service, Johannesburg, South Africa
Full list of author information is available at the end of the article</footnote>
**Conclusions:** It was logistically feasible to conduct mark-release-recapture studies at the current scale. The population size estimates obtained may provide a guideline for the initial number of males to use for a pending SIT pilot trial. It is promising for future SIT trials that laboratory-reared marked males participated in natural swarms, appearing at the right place at the right time.

**Keywords:** Malaria vector control, Malaria elimination, Male mosquitoes, Population size, Over-flooding ratio, Sterile insect technique

## Background

South Africa's malarious areas include the northeastern border regions of three provinces: KwaZulu-Natal (KZN), Mpumalanga and Limpopo [1]. These provinces border Eswatini (formally Swaziland), Mozambique, Zimbabwe and Botswana. Of these countries, Mozambique is an especially high-risk area for malaria [1] and poses challenges for malaria elimination initiatives in South Africa. In KZN, the highest malaria incidence occurs in the northeastern part of the province, particularly near the border with Mozambique [2].

The major malaria vectors in South Africa are *Anopheles funestus* and *An. arabiensis* [3], with the former now nearly eliminated by the provincial indoor residual spraying (IRS) programmes and the latter recently implicated in ongoing residual transmission [4]. Potential secondary vectors include *An. merus*, *An. coustani*, *An. marshallii*, *An. pharoesnis*, *An. rivulorum* and *An. leesoni* [5]. *Anopheles vaneedeni* and *An. parensis* have recently been implicated as secondary vectors in South Africa [6, 7].

In the almost complete absence of *An. funestus*, *An. arabiensis*, a species that displays variable feeding and resting behaviour, is currently implicated as the major vector in the ongoing residual malaria transmission in South Africa [4]. *Anopheles arabiensis* will feed off animals (predominantly cattle) as well as humans and may rest outdoors as well as indoors [8–10]. This poses a problem for South Africa's malaria elimination agenda [11] because the provincial malaria control programmes rely heavily on IRS and, to a lesser extent, larval source management, for vector control [3, 11, 12]. Outdoor-resting vector populations are ineffectively targeted by IRS, requiring the adoption of additional methods for their control.

Additional vector control tools that can target outdoor-feeding and resting vectors and are not insecticide based are urgently required. Environmentally friendly and species-specific, the sterile insect technique (SIT) is one such tool under consideration. SIT, a concept conceived by Knipling in the 1930s [13], has previously been considered for use against malaria vector mosquitoes, but was abandoned due to various factors, including public perceptions and political instability [14]. The technique has been used successfully against several insect pests,

including the New World screwworm fly [15], pink boll-worm [16], fruit flies [17], tsetse fly [18] and codling moth [19]. SIT involves the mass production and release of sterilized laboratory-reared male insects that have retained sufficient fitness to compete with wild males for wild females. Inundative releases of sterile males have been shown to induce sufficient sterility in wild females to effectively suppress, and in some instances eliminate, targeted insect populations [15–19].

A project to assess the feasibility of using the SIT for the control of *An. arabiensis* in South Africa has been underway since 2011 [20]. A suitable field site for pilot field releases has been established at Mamfene, northern KwaZulu-Natal Province [4]. Regular vector surveillance and community engagement has been ongoing at this site since 2014 [4, 21], providing baseline data for annual and seasonal fluctuations in malaria vector density and composition. An *An. arabiensis* genetic sexing strain has been developed by introgressing a genetic sexing mechanism into the genetic background of the target population at Mamfene [22, 23]. The optimal irradiation dose for sterilizing males has been established, and the competitiveness of wild versus laboratory-reared sterilized males determined [23, 24].

The Mamfene field site is well suited for an SIT intervention as the malaria incidence in KwaZulu-Natal is currently less than one case per 1000 residents [11]. The site is located in a malaria hotspot where annual vector control using IRS is conducted (National Department of Health). The surrounding areas are not sprayed unless a malaria case is detected, in which case spraying is conducted using foci clearing protocols [11].

Before a pilot SIT trial can be implemented, it is necessary to estimate the *An. arabiensis* population size at the target site in order to be able to estimate the number of sterilized males required to inundate the target population at the required ratio (over-flooding ratio). One of the most common methods used to estimate population size is the mark-release-recapture technique (MRR) [25]. Using MRR, population size can be estimated by releasing a known number of marked individuals into the target area. Individuals of that species are subsequently captured and, based on the proportion of marked individuals captured compared to wild (unmarked)
individuals captured, an estimate of the population size can be calculated. MRR can also be used to measure dispersal [25, 26]. Many MRR studies on mosquitoes have been conducted [25–30], mostly on females. However, some recent studies have focused on males [27, 28, 30].

The aim of this study was to estimate the *An. arabiensis* population size at the Mamfene intervention site, as well as gain insight into dispersal ability of laboratory-reared males. In addition, the logistics and feasibility of producing laboratory-reared males in Johannesburg and then transporting and releasing them in Mamfene, KwaZulu-Natal Province were determined.

# Methods

## Field experimental site

All field activities were conducted at Mamfene in the uMkhanyakude district of KwaZulu-Natal Province (Fig. 1a). The experimental area is divided into three sites, two control sites and one target or intervention

site (Fig. 1b). Mamfene is predominantly rural with subsistence farming at the household level, and most of the surrounding areas are occupied by sugarcane plantations. The local community is predominantly housed in small westernized buildings with a few isolated traditional mud structures present. Cattle and goats are kept by some members of the community and many of them are allowed to graze freely in the fields adjacent to the households, and along the roads. During the rainy season, puddles are formed, creating potential oviposition sites for gravid *An. arabiensis* which tend to oviposit in small temporary water bodies dependent on rainfall [10, 31]. There is a marshy area located between control site one and the target site that also provides suitable mosquito larval sites. A stream runs along the southern side of control site two. Rainfall is highest from September to April, coinciding with the malaria season. Average temperature is 22.3 °C, with June being the

<figure><img src="image_3.png" /><figcaption>**Fig. 1** a Location of the study site in KwaZulu-Natal Province, South Africa. b View of the Jozini Municipality, showing the location of Mamfene. c Study site showing the target site and the two control sites. Control site one is adjacent to the target site and control site two is approximately 8 km away. The area of each section is: target site, 0.48 km²; control site one, 0.70 km²; control site two, 0.26 km²</figcaption></figure>
coolest and February the hottest month, respectively. Average rainfall is 569 mm per year [32].

## Biological material

The mosquitoes were drawn from KWAG, an *An. arabiensis* colony that originated from Mamfene, KwaZulu-Natal Province and established in 2005 [33]. This colony has been reared through many generations in the laboratory since its establishment; however, generation numbers were not recorded. KWAG is maintained in the Hugh Patterson insectary at the National Institute for Communicable Diseases (NICD) in Johannesburg.

## Production of KWAG males for release

The Hugh Patterson insectary has two rearing rooms. The room reserved for KWAG larval stages was maintained at (mean ± standard deviation [SD]) 32±3 °C and 80±5% relative humidity (RH) to reduce larval developmental time. Larvae were fed a mixture of finely ground dog biscuits (Beeno®, Kempton Park, South Africa) and brewer's yeast (Vital®; Vital Health Foods, Cape Town, South Africa) [34]. A second room was used to house the adults and was maintained at the standard (25 °C) rearing temperature [34]. All adult mosquitoes were fed 10% sugar solution *ad libitum* and adult females were offered blood meals at least twice per week using an artificial membrane feeding system (Hemotek®; Hemotek Ltd., Blackburn, UK).

To obtain males for MRR, pupae were separated from larvae in batches using the vortex induction technique described in the International Atomic Energy Agency (IAEA) guidelines [35]. The number of pupae collected was established using volumetric estimation in batches of no greater than 5000 pupae [35]. After separation from larvae, pupae were placed into BugDorm® cages (Bug-Dorm-1: 30 × 30 × 30 cm; BugDorm, MegaView Science Co., Ltd, Taichung, Taiwan; hereafter referred to as “cage”) in batches of 1500 pupae per cage, provided with 10% sugar solution and allowed to emerge overnight. One to two days after adult emergence males were separated from females based on antennal morphology and transferred into new cages at a density of 1000 males per cage. These were maintained on 10% sugar solution *ad libitum* until marking.

## Marking of adult males prior to release

Males were marked with fluorescent dust (Day Glo series; Day-Glo Color Corp., Cleveland, OH, USA) 2 to 3 days prior to release. The following dust colours were used: pink, yellow and orange (see Additional file 1: Table S1). Prior to marking, the sugar water was removed from each cage, and the number of dead males was recorded and the dead males removed. The mosquitoes were

immobilized by placing each cage into a refrigerator for 2–5 min. Immobilized mosquitoes were tipped into a pre-dusted plastic container using a large funnel. Plastic containers (1 l) were prepared for dusting 1000 males by adding 0.05 g of fluorescent dust, weighed using an analytical balance (model Kern ABT 120-5DM, United Scientific Equipment Pte Ltd, Singapore; model Adam PW124, Labex, 88 17th Avenue, Edenvale, Johannesburg, 1609) according to the methods developed by Culbert et al. [36], but modified for larger numbers. Any escaped mosquitoes or those not successfully transferred into the dusting container were counted, recorded and excluded from the calculation of final male release number. The 10% sugar solution was provided *ad libitum* immediately following marking. This method successfully marked all of the males. In a similar separate experiment, no significant effect on longevity on marked versus unmarked males was observed [36]. For the mark-release-recapture experiments with more than one release location, two dust colours were used, one for each site.

## Cloth tube construction

Previously designed sentinel style resting sites, named cloth tubes (Fig. 2), were used for collections [24]. Cloth tubes with a diameter of 18–20 cm and length of 50–54 cm were constructed using a tubular frame made from plastic garden mesh fencing, a black fabric skin and a pair of netting sleeves capping each end. A couple of hand-stitched tacks between the frame and fabric were used to keep the frame and skin together. Laminated labels were stapled to the hem for trap identification, and nails were driven through the hem to anchor the traps to the ground when deployed. The cloth tubes for the pilot trial were capped with a net at one end, with the other end left open; however, the design was subsequently improved, as described above, so that it was possible to close both ends during collection. This latter design prevented the escape of any mosquitoes inside and also aided in distributing and collecting more tubes at once (Fig. 2a). Cloth tubes were placed either horizontally or at an angle of approximately 45° to the ground when deployed (Fig. 2b, c).

## Transport of males from the insectary to the field laboratory and release sites

The NICD insectary is located approximately 500 km from the field site by road. Marked males were transported in BugDorm-1® cages placed in the enclosed canopy at the back of a field vehicle (pick-up truck), at ambient temperature with access to 10% sugar solution. To maintain humidity and a cooler temperature during transportation, damp towels were draped over the cages. This method is used routinely by the laboratory
<figure><img src="image_4.png" /><figcaption>**Fig. 2** a Cloth tubes used during mark-release-recapture experiments strung onto a pole for collection and deployment of multiple tubes at once.</figcaption></figure>

<caption>b, c Typical placement in the field. Note the netting sleeves (one closed and one open) used to close the tube for more rapid collection</caption>

to maintain mosquitoes in good condition during transport. On arrival at the field laboratory, marked males were left inside the transport vehicles overnight. These males were provided with fresh 10% sugar solution the following morning and moved into the field laboratory until release. During this holding phase, all cages were visually scanned for females aided by a volunteer who placed an arm onto the mesh side of the cage. Any females found were recorded and removed. The number of females found was generally low, ≤ 5 per cage (0.5%), with some cages free of females. Males were transported in their cages to the release sites in the field vehicles. For mark-release-recapture one (see following section), with two release sites, two field vehicles were used, one for each release site (and therefore one vehicle for each dust colour).

**Release location selection and the setting up of cloth tubes in the field**

Three releases, referred to as mark-release-recapture pilot trial (MRR<sup>PT</sup>), the mark-release-recapture one (MRRI) experiment and the mark-release-recapture two (MRRII) experiment, were conducted between November 2016 and October 2018 (Additional file 1: Table S1).

The release site (household) used in MRR<sup>PT</sup>, MRRI and MRRII was selected based on its central position in the target site (Figs. 3–5). This site was named C (centre), with co-ordinates 27°24′19.68″S; 32°12′34.37″E. For MRRI, an additional release point at the edge of the target site, approximately 550 m from the central release site, was selected and named E (edge), with co-ordinates 27°24′7.70″S; 32°12′49.64″E (Fig. 4). Prior to setting up the tubes, permission was sought from household owners. The density of the cloth tubes deployed to the field
<figure><img src="image_5.png" /></figure>

increased with each release in an attempt to improve the chances of capturing and recapturing mosquitoes (Table 1). In MRR<sup>PT</sup>, the cloth tubes were placed at a maximum distance of 400 m (Fig. 3; Table 1) from the release site; however, in subsequent experiments they were placed up to a maximum of 150 m from the release site as marked males were only captured in tubes up to 100 m from the release site during MRR<sup>PT</sup> (Figs. 4, 5; Table 1). The position of each cloth tube was mapped using a hand-held GPS device (Garmin eTrex 20×; Garmin Ltd., Schaffhausen, Switzerland) (Figs. 3–5).

## Male releases

The males were transported (as described previously) from the field laboratory to the release locations in Bug-Dorm-1 cages. The duration of the trip from the field laboratory to release sites is approximately 25 min, and arrival at the release locations was timed to be at least 30 min before the planned releases. All releases were conducted outdoors between 10 and 20 min before sunset using 2- to 5-day-old marked males.

To release males, the cage tops were removed and the cages tapped to encourage the marked males to fly out. Any dead males, or those incapable of flying out of the cages, were retained. Cages were re-closed after releases and the remaining males later counted to calculate the actual number of males released. In the MRRI experiment, males were simultaneously released at two sites (Fig. 4) by two independent release teams.

## Mosquito collection from cloth tubes

All collections occurred daily, beginning at dawn just as it became light, and were complete within 1–2 h. During MRR<sup>PT</sup>, mosquitoes were collected directly from the cloth tubes (which were open at one end and capped at the other end) in the field by illuminating the trap using a torch light and collecting any mosquitoes using a manual aspirator. Collected mosquitoes were immediately placed into polystyrene cups labelled with the tube name and date of collection, and subsequently transported to the field laboratory for morphological species identification using standard keys [8, 9]. During subsequent releases (MRRI and II), entire cloth tubes were collected and
<figure><img src="image_6.png" /></figure>

replaced straightaway with a new ones labelled as previ-
ously described. Cloth tubes were placed into the field
with one sleeve closed and the other end open to allow
mosquito entry. Upon collection, the open sleeve was
closed for transport to the laboratory.

Swarm searches

Daily swarm searches were performed each evening fol-
lowing releases. Searches started approximately half an
hour before sunset and ended when it became too dark
to see. Swarms were detected by observing the hori-
zon towards the sun in potential swarming areas, as
described in Diabaté et al. [37]. Potential swarming sites
were identified above small bushes, walking path inter-
sections, above mounds or depressions in the ground and
near the leafy border of tall trees. Such contrasting areas
are considered to be potential swarm markers [38, 39].
Samples were collected from identified swarms (Fig. 4)
using a sweep net (except directly after releases). Sweep
nets were constructed specifically for collecting swarms
using a plastic hula-hoop (diameter: approx. 70 cm) fit-
ted with a large conical net (length: 1.5 m) and a square

lashed to a notched handle (length: 1.8–2 m) made of
thick bamboo (Additional file 2). Swarm catches were
immediately inspected (in sweep nets) by torch light for
the presence of mosquitoes. If mosquitoes were present,
these were transferred to polystyrene cups using an aspi-
rator and kept for morphological identification [8, 9] in
the field laboratory the following day. The date, time,
swarm height, swarm size and location of each swarm
were recorded.

Mosquito identification and processing in the laboratory

All collected mosquito specimens were transported to
the field laboratory daily for initial processing. Each cloth
tube taken to the laboratory was checked for the presence
of mosquitoes by looking through one end of the tube
towards the light while gently tapping and rotating the
tube. If mosquitoes were found, these were carefully aspi-
rated out. The cloth tubes were then fully opened onto a
sheet of white paper to recover any dead mosquitoes. All
anopheline mosquitoes found in each tube were placed
into disposable cups and appropriately labelled with the
tube name and collection date. Anopheline mosquitoes
<figure><img src="image_7.png" /><figcaption>**Table 1** Number of cloth tubes by distance from release site for mark-release-recapture pilot trial and mark-release-recapture one and two experiments</figcaption></figure>

<table><thead><tr><th rowspan="2">Distance of cloth tube from the release point (m)</th><th colspan="3">MRR experiment (release location) and number of cloth tubes per radius</th></tr><tr><th>MRR<sup>PT</sup> (C)</th><th>MRI (C and E)</th><th>MRII (C)</th></tr></thead><tbody><tr><td>25</td><td>–</td><td>4</td><td>4</td></tr><tr><td>50</td><td>5</td><td>6</td><td>8</td></tr><tr><td>75</td><td>–</td><td>8</td><td>13</td></tr><tr><td>100</td><td>8</td><td>11</td><td>18</td></tr><tr><td>125</td><td>–</td><td>14</td><td>24</td></tr><tr><td>150</td><td>–</td><td>17<sup>a</sup></td><td>29</td></tr><tr><td>200</td><td>11</td><td>–</td><td>–</td></tr><tr><td>400</td><td>12</td><td>–</td><td>–</td></tr><tr><td>Total cloth tubes</td><td>36</td><td>60<sup>a</sup></td><td>96</td></tr></tbody></table>

C, Central release point; E, edge release point; MRR<sup>PT</sup>, mark-release-recapture pilot trial; MRRI, MRRII, mark-release-recapture one and two experiments, respectively

collected alive were anaesthetised using ethyl acetate or dry ice before morphological identification [8, 9]. All mosquitoes were scanned for fluorescent dust markings under the microscope with the aid of a UV torch (9 LED UV Scorpion flashlight; Zartek, Johannesburg, South Africa), although markings were readily visible by eye. After this initial processing, all specimens were stored in individual silica tubes and transported to Johannesburg by road at the completion of the field work for further analyses. Individual specimen records were maintained on tube identification, date, species (by morphology) and whether males captured were marked (including colour) or unmarked (wild males). Once in the Johannesburg laboratory, *An. gambiae* complex and *An. funestus* group specimens were identified to species using the PCR-based methods of Scott et al. [40] and Koekemoer et al. [41], respectively.

## Dispersal distance

The average dispersal of released males for each experiment was calculated by adding the distance of each trap (from the release point) in which each marked male was

<footnote><sup>a</sup> For MRRI, E had an additional trap at 150 m (two traps were mistakenly labelled 150.14 so one became 150.14B), bringing the total number of traps at E to 61</footnote>
collected and then dividing this by the total number of males recaptured for that experiment. Maximum dispersal distance of marked males was also determined.

## Population size estimation and calculation of daily survival probability and average life expectancy

The number of *An. arabiensis* males was estimated using two modified Lincoln indexes for low recapture rates. Both indices were used to give a minimum and a maximum number for release. This was done out of concern about releasing too few sterile males as this would reduce the efficacy of the potential SIT trial.

The first index takes into account daily survival rate ([42] as referenced in [43]) according to the following equation: $P = [as^t(n - r + 1)/(r + 1)$, where P is the estimated population size, a is the number of marked males released, s is the estimated probability of daily survival, t is the sampling day post release, n is the total number of marked and unmarked males captured and r is the number of marked males recaptured ([42], as referenced in [43]). Daily survival probability (s in equation) is calculated by the regression of the total number of marked males transformed by $log(x + 1)$ in all traps on the day post release. This was calculated for 5 days after release. Survival probability was estimated using Statistix version 8.0 (Analytical Software, Tallahassee, FL, USA) from the results of the antilog of the slope of the regression line [44, 45]. Average life expectancy was calculated by the formula: $1/ - \ln(s)$ (see [46]).

The second index is the standard Lincoln index used to estimate populations from MRR experiments when r, the number of recaptured individuals, is < 20 [43]. This index assumes constant daily survival and is described by the equation: $P = \frac{a(n+1)}{r+1}$

, where a is the number of samples originally marked, n is the total number of recaptures and r is the number of recaptured marked individuals [43].

The estimated number of males is expressed as the total population estimate and as the number of males per hectare.

## Weather data

All weather data were obtained from the Makhathini Research Station (27°23′42.45″S; 32°10′48.48″E), located approximately 3.2 km from release location C. Statistical analyses were performed using Statistix version 8.0. Mean temperature, relative humidity and mean daily wind speed on the day of release and subsequent mosquito collection days were calculated. Mean temperature and relative humidity for the release and collection periods were compared for each experiment using one-way analysis of variance (ANOVA); if significant, a Bonferroni post-hoc test was conducted to separate the means.

Average wind direction and wind speed were also determined for each release. In addition, the rainfall pattern was investigated for the 4-week period prior to MRRII.

## Results

### Numbers released, percentage recapture, survival probability, average life expectancy and population size estimation

The number of marked males released, the number recaptured, the daily survival probability and the estimated population size for each experiment are given in Table 2. The population sizes for the two release locations used in MRRI were calculated separately. The largest population size estimate was obtained during MRRI at E and the smallest was obtained during MRRII at C. The percentage recapture of marked males from deployed cloth tubes was 0.09, 0.11, 0.06 and 0.04% for MRR<sup>PT</sup>, MRRI at C, MRRI at E and MRRII, respectively. The average life expectancy for each experiment, using the data from cloth tubes only (including the 1 clay pot used in MRR<sup>PT</sup>), was 2.94 days for MRR<sup>PT</sup>, 2.59 days for MRRI C, 4.48 days for MRRI E and 2.59 days for MRRII. The average *An. arabiensis* male population sizes, as calculated for the first 2 days of collections, were 539 (MRR<sup>PT</sup>), 783 (MRRI at C), 6213 (MRRI at E) and 401 (MRRII) males per hectare using the index that takes daily survival into account [42]. Using the simple index [43], the population size averages for the first 2 days after releases were: 2819 (MRR<sup>PT</sup>) 1936 (MRRI at C), 9643 (MRRI at E) and 2193 (MRRII) males per hectare (see Table 2 for more details).

### Swarms

Released males were observed to form their own swarms on the day of release near the release location during MRR<sup>PT</sup> and MRR II. Figure 4 indicates swarms that were found during MRRI. Of these, 55.55% (n=5 swarms) contained marked males (Table 3). Notably, a single marked male was recaptured 8 days after release while checking for swarms during routine surveillance work (Table 3). Although swarms were found during MRRII, these were *An. marshallii* complex with a few *An. cossutani*. Table 4 shows the number of marked males released per site in MRRI and the number of *An. arabiensis* males captured in swarms. Due to the limited amount of data obtained from swarms in this series of experiments, no population size estimates were made from the swarms.

## Species composition, distribution and dispersal
### Pilot trial (MRR<sup>PT</sup>)

During MRR<sup>PT</sup> 49 mosquitoes were collected from 12 of the 36 cloth tubes. Of these 80% (n=39) were *An. gambiae* complex, 10% (n=5) were *An. funestus* species group and 8% (n=4) were other anophelines. The
<table><caption>**Table 2** Estimated total and per hectare population size of *Anopheles arabiensis* males in Mamfene, KwaZulu-Natal, South Africa</caption><thead><tr><th>Experiment/Day</th><th>Number released (day 0)</th><th>Number of wild males</th><th>Number of marked males</th><th>Survival probability per day (s<sup>1</sup>)<sup>a</sup></th><th>Estimate 1 (per ha)<sup>b</sup></th><th>Estimate 2 (per ha)<sup>b</sup></th></tr></thead><tbody><tr><td>MRR<sup>PT</sup> C/1</td><td rowspan="6">7871</td><td>0</td><td>3</td><td>0.71</td><td>1397 (445)</td><td>7871 (2505)</td></tr><tr><td>MRR<sup>PT</sup> C/2</td><td>1</td><td>3</td><td>0.5</td><td>1984 (632)</td><td>9839 (3132)</td></tr><tr><td>MRR<sup>PT</sup> C/3</td><td>0</td><td>0</td><td>0.36</td><td>2817 (897)</td><td>7871 (2505)</td></tr><tr><td>MRR<sup>PT</sup> C/4</td><td>0</td><td>1</td><td>0.25</td><td>1000 (318)</td><td>7871 (2505)</td></tr><tr><td>MRR<sup>PT</sup> C/5</td><td>2</td><td>0</td><td>0.18</td><td>4260 (1356)</td><td>23,613 (7516)</td></tr><tr><td>MRRI C/1</td><td rowspan="6">5262</td><td>6</td><td>4</td><td>0.68</td><td>5009 (709)</td><td>11,576 (1638)</td></tr><tr><td>MRRI C/2</td><td>4</td><td>1</td><td>0.46</td><td>6051 (856)</td><td>15,786 (2233)</td></tr><tr><td>MRRI C/3</td><td>4</td><td>1</td><td>0.31</td><td>4078 (577)</td><td>15,786 (2233)</td></tr><tr><td>MRRI C/4</td><td>4</td><td>0</td><td>0.21</td><td>5525 (782)</td><td>26,310 (3722)</td></tr><tr><td>MRRI C/5</td><td>2</td><td>1<sup>c</sup></td><td>0.15</td><td>2368 (335)</td><td>15,786 (2233)</td></tr><tr><td>MRRI E/1</td><td rowspan="6">5453</td><td>9</td><td>2</td><td>0.8</td><td>14,541 (2057)</td><td>21,812 (3086)</td></tr><tr><td>MRRI E/2</td><td>20</td><td>0</td><td>0.64</td><td>73,288 (10,368)</td><td>114,513 (16,200)</td></tr><tr><td>MRRI E/3</td><td>30</td><td>1</td><td>0.51</td><td>43,275 (6122)</td><td>87,248 (12,343)</td></tr><tr><td>MRRI E/4</td><td>8</td><td>0</td><td>0.41</td><td>20,102 (2844)</td><td>49,077 (6943)</td></tr><tr><td>MRRI E/5</td><td>13</td><td>0</td><td>0.33</td><td>25,040 (3542)</td><td>76,342 (10,800)</td></tr><tr><td>MRRII C/1</td><td rowspan="6">~ 15,500</td><td>0</td><td>4</td><td>0.68</td><td>2108 (298)</td><td>15,500 (2193)</td></tr><tr><td>MRRII C/2</td><td>0</td><td>1</td><td>0.46</td><td>3565 (504)</td><td>15,500 (2193)</td></tr><tr><td>MRRII C/3</td><td>0</td><td>1</td><td>0.31</td><td>2403 (340)</td><td>15,500 (2193)</td></tr><tr><td>MRRII C/4</td><td>0</td><td>0</td><td>0.21</td><td>3255 (460)</td><td>15,500 (2193)</td></tr><tr><td>MRRII C/5</td><td>1</td><td>0</td><td>0.15</td><td>4650 (658)</td><td>31,000 (4386)</td></tr></tbody></table>

<table><caption>**Table 3** *Anopheles* male mosquitoes captured in swarms (ordered by day post release) at the central and edge locations during MRRI in Mamfene, Jozini, KwaZulu-Natal Province, South Africa</caption><thead><tr><th>Day</th><th>Location (GPS co-ordinates)</th><th>Marker</th><th>Anopheles species present</th><th>Number of marked to wild An. arabiensis males</th></tr></thead><tbody><tr><td>1</td><td>1+ km from E (27°23′52.16″S; 32°13′24.84 E)</td><td>Ground depression</td><td>An. arabiensis</td><td>0:12</td></tr><tr><td>2</td><td>150 m from E in Marsh (27°24′4.38″S; 32°12′45.55″E)</td><td>Vegetation on edge of marshy area</td><td>Other species</td><td>NA</td></tr><tr><td>2</td><td>200 m from E in Marsh (27°24′5.19″S; 32°12′42.59″E)</td><td>Mound in marsh</td><td>An. pharoensis</td><td>NA</td></tr><tr><td>2</td><td>125 m from E (27°24′6.31″S; 32°12′45.66″E)</td><td>Mound with log</td><td>An. arabiensis</td><td>1or:4</td></tr><tr><td>3</td><td>125 m from E (27°24′6.31″S; 32°12′45.66″E)</td><td>Mound with log</td><td>An. pharoensis, An arabiensis</td><td>0:1</td></tr><tr><td>3</td><td>95 m from E (27°24′5.01″S; 32°12′47.88″E)</td><td>Small bush</td><td>An. arabiensis</td><td>0:14</td></tr><tr><td>3</td><td>55 m from C (27°24′18.68″S; 32°12′32.88″E)</td><td>Bare ground</td><td>An. arabiensis</td><td>1ye:1</td></tr><tr><td>3</td><td>60 m from C (27°24′17.87″S; 32°12′33.84″E)</td><td>Between two bushes (~ 9 m to each bush)</td><td>An. arabiensis</td><td>3ye:9</td></tr><tr><td>3</td><td>215 m from C (27°24′15.95″S; 32°12′27.89″E)</td><td>Clearing between two bushes</td><td>An. arabiensis, An. parensis</td><td>1ye:9</td></tr><tr><td>4</td><td>95 m from E (27°24′5.01″S; 32°12′47.88″E)</td><td>Small bush</td><td>An. arabiensis</td><td>0:2</td></tr><tr><td>8</td><td>210 m from C (27°24′17.26″S; 32°12′27.35″E)</td><td>Large bush near CP 9B</td><td>An. arabiensis</td><td>1ye:0</td></tr></tbody></table>

CP, Clay pot; NA, not applicable (no *An. arabiensis* in swarm); or, orange marked male; ye, yellow marked male

<footnote><sup>a</sup> Survival probability as calculated for estimate 1</footnote>

<footnote><sup>b</sup> Estimate 1 was obtained using the procedure of [42], and Estimate 2 using [43]</footnote>

<footnote><sup>c</sup> An orange marked male was captured at C which had travelled over 600 m to reach the trap at E. This male was excluded from the population size estimation because it was beyond the area used for the calculation, reducing the total captured to 26 and the marked recaptures to 6 for C during MRRI</footnote>
<table><caption>**Table 4** Total *Anopheles arabiensis* males and marked released males captured from swarms only during MRRI, Mamfene, Jozini, KwaZulu-Natal Province, South Africa</caption><thead><tr><th>Release point</th><th>Number of marked males released</th><th>Total <i>An. arabiensis</i> males in swarms</th><th>Total marked males recaptured (swarms)</th></tr></thead><tbody><tr><td>C</td><td>5262</td><td>24</td><td>5</td></tr><tr><td>E</td><td>5453</td><td>22</td><td>1</td></tr></tbody></table>

DNA of one anopheline female, identified to the *An. gambiae* complex by morphology, did not amplify on multiple attempts (Table 5). Seven of the 16 *An. arabiensis* males (44%) were recaptured marked males. Marked males were only detected until the fourth day of 7 days of collections.

Analysis of *An. arabiensis* caught per cloth tube and distance from the release point showed that these were not equally productive in terms of numbers caught (Table 6). Resting site 200.8 was actually the clay pot that is used for routine surveillance which was positioned within 1 m of the location chosen for the cloth tube. The clay pot was, therefore, used instead. It should be noted that mosquitoes collected beyond 100 m from the release point were not used for further analyses as no marked males were recaptured beyond 100 m during MRR<sup>PT</sup>. Marked males (7 in total) were found only in the cloth tubes placed at 50 m (two marked males) and 100 m (5 marked males) (Table 6), making the average dispersal distance of marked males 85.7 m.

### MRRI

In MRRI, many more mosquitoes were captured at E (n=197) than at C (n=54) (Table 6). *Anopheles arabiensis* made up the majority of the anophelines collected, comprising 64% (n=126) and 79% (n=43) of collections at E and C, respectively. *Anopheles parensis* (of the *An. funestus* group) was the second most abundant species,

at 28.5% (n=56) and 9.3% (n=5) at E and C, respectively. *Anopheles quadriannulatus*, *An. merus* and *An. rivulorum* also contributed to the collections (Table 7). Three of 83 (3.6%) *An. arabiensis* males at E and seven of 27 (26%) at C were marked males (Table 8). Marked males were recaptured up to 125 m from release point C (one orange-marked male and one yellow-marked male), but only up to 75 m from release point E (Table 7). The average dispersal distance for released males recaptured from cloth tubes was 70.83 m for males released at C and 58.33 m for males released at E; maximum dispersal distances were 125 and 75 m, respectively. However, one orange-marked male released at E was recaptured at C on day 5 of recaptures. This male had reached a cloth tube located on the 125 m radius from C, approximately 600 m from its original release location (Fig. 4). If this male is included in the average dispersal calculation for males released at E, the average dispersal distance increases to about 194 m.

### MRRII experiment

Sixteen anophelines were captured over five consecutive collection days from six of the 96 cloth tubes (Table 9), of which 50% were *An. arabiensis* (n=8). The remainder of the wild specimens captured were morphologically identified as *An. marshallii* complex with one *An. funestus* group female. Six of the seven (85.72%) *An. arabiensis* males captured were marked males.

### Weather data

The mean daily temperature differed significantly between experiments (ANOVA: $F_{(2, 15)}=18.7 P<0.001$). Bonferroni post-hoc analysis (critical T=2.694) showed that the mean (± SD) temperature was significantly lower during MRRII (19.20 ± 2.17 °C) than during MRR<sup>PT</sup> (24.78 ± 2.30 °C) and MRRI (25.38 ± 1.20 °C). The means (± SD) of the daily mean minimum temperatures followed a similar trend, being significantly

<table><caption>**Table 5** Summary of anopheline mosquitoes by gender and species sampled from cloth tubes during MRR<sup>PT</sup>, Mamfene, Jozini, KwaZulu-Natal Province, South Africa</caption><thead><tr><th>Species</th><th>Number of males</th><th>Number of females</th><th>Total (%)</th></tr></thead><tbody><tr><td>Marked <i>An. arabiensis</i></td><td>7</td><td>NA</td><td>7 (14.3)</td></tr><tr><td>Wild <i>An. arabiensis</i></td><td>9</td><td>21</td><td>30 (61.2)</td></tr><tr><td>Other <i>An. gambiae</i> complex</td><td>0</td><td>2 (<i>An. merus</i>)</td><td>2 (4.1)</td></tr><tr><td>Total <i>An. gambiae</i> complex</td><td>16</td><td>23</td><td>39 (79.6)</td></tr><tr><td><i>An. funestus</i> group</td><td>1 (<i>An. rivulorum</i>)</td><td>4 (<i>An. parensis</i>)</td><td>5 (10.2)</td></tr><tr><td>Other anopheline species</td><td>1</td><td>3</td><td>4 (8.2)</td></tr><tr><td>Failed amplification</td><td>0</td><td>1</td><td>1 (2.0)</td></tr><tr><td>Total collected</td><td>18</td><td>31</td><td>49 (100)</td></tr></tbody></table>

NA, Not applicable (there were no marked *An. arabiensis* females)
<table><caption>**Table 6** *Anopheles arabiensis* mosquitoes captured per cloth tube per day during MRR<sup>PT</sup>, Mamfene, Jozini, KwaZulu-Natal Province, South Africa</caption><thead><tr><th rowspan="2">Collection day</th><th colspan="11">Cloth tube number<sup>a</sup></th><th rowspan="2">Totals</th></tr><tr><th>50.1</th><th>50.5</th><th>100.2</th><th>100.4</th><th>100.6</th><th>100.7</th><th>200.5</th><th>200.8<sup>b</sup></th><th>200.9</th><th>400.2</th><th>400.7</th><th>400.12</th></tr></thead><tbody><tr><td colspan="14"><strong>Day 1</strong></td></tr><tr><td>Marked ♂</td><td>2</td><td>-</td><td>-</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>3</td></tr><tr><td>Wild ♂</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>3</td><td>-</td><td>-</td><td>-</td><td>-</td><td>3</td></tr><tr><td>Wild ♀</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td><td>3</td><td>-</td><td>-</td><td>-</td><td>-</td><td>4</td></tr><tr><td colspan="14"><strong>Day 2</strong></td></tr><tr><td>Marked ♂</td><td>-</td><td>-</td><td>-</td><td>-</td><td>2</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>3</td></tr><tr><td>Wild ♂</td><td>-</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>2</td></tr><tr><td>Wild ♀</td><td>-</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>2</td><td>-</td><td>-</td><td>-</td><td>-</td><td>3</td></tr><tr><td colspan="14"><strong>Day 3</strong></td></tr><tr><td>Marked ♂</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>0</td></tr><tr><td>Wild ♂</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>0</td></tr><tr><td>Wild ♀</td><td>-</td><td>-</td><td>-</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td><td>1</td><td>-</td><td>-</td><td>3</td></tr><tr><td colspan="14"><strong>Day 4</strong></td></tr><tr><td>Marked ♂</td><td>-</td><td>-</td><td>-</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td></tr><tr><td>Wild ♂</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td></tr><tr><td>Wild ♀</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td><td>-</td><td>1</td><td>1</td><td>-</td><td>3</td></tr><tr><td colspan="14"><strong>Day 5</strong></td></tr><tr><td>Marked ♂</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>0</td></tr><tr><td>Wild ♂</td><td>-</td><td>-</td><td>-</td><td>2</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>2</td></tr><tr><td>Wild ♀</td><td>-</td><td>-</td><td>-</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td></tr><tr><td colspan="14"><strong>Day 6</strong></td></tr><tr><td>Marked ♂</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>0</td></tr><tr><td>Wild ♂</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td></tr><tr><td>Wild ♀</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>2</td><td>-</td><td>-</td><td>-</td><td>1</td><td>3</td></tr><tr><td colspan="14"><strong>Day 7</strong></td></tr><tr><td>Marked ♂</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>0</td></tr><tr><td>Wild ♂</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>0</td></tr><tr><td>Wild ♀</td><td>-</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>2</td><td>-</td><td>-</td><td>3</td></tr><tr><td colspan="14"><strong>Totals</strong></td></tr><tr><td>Marked ♂</td><td>2</td><td>0</td><td>0</td><td>2</td><td>2</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>7</td></tr><tr><td>Wild ♂</td><td>0</td><td>1</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>6</td><td>0</td><td>0</td><td>0</td><td>0</td><td>9</td></tr><tr><td>Wild ♀</td><td>0</td><td>2</td><td>1</td><td>1</td><td>0</td><td>0</td><td>1</td><td>8</td><td>1</td><td>4</td><td>1</td><td>1</td><td>20</td></tr></tbody></table>

Only the cloth tubes from which *An. arabiensis* mosquitoes were collected are shown

lower (ANOVA: $F_{(2, 15)}=34.8, P < 0.001$) during MRRII (12.75 ± 2.11 °C) than during MRR<sup>PT</sup> (19.58 ± 1.87 °C) and MRRI (19.80 ± 0.58 °C) following Bonferroni correction (critical $T=2.694$). Humidity did not differ significantly between the experiments (ANOVA: $F_{(2, 15)}=0.73, P>0.50$), with means (± SD) of 71.67 ± 10.61, 74.67 ± 4.27 and 76.33 ± 2.33% RH for MRR<sup>PT</sup>, MRRI and MRRII, respectively. Wind speed was significantly lower during MRRI (mean = 1.08 m/s, SD = 0.098 m/s) than during MRR<sup>PT</sup> (mean = 2.23 m/s, SD = 0.716 m/s) and

MRRII (mean = 2.30 m/s, SD = 0.596 m/s) (ANOVA: $F_{(2, 15)}=9.61, P=0.002$) after Bonferroni adjustment (critical $T=2.694$); however, there was no significant difference between MRR<sup>PT</sup> and MRRII (Additional file 1: Table S2). No rain was recorded during MRR<sup>PT</sup>. During MRRI there was very light precipitation during the morning of the fourth collection day (3 March 2018). During MRRII, rainfall occurred 4 weeks prior to release as well as during the week of the release on the following days: 3 October (4.2 mm), 4 October (8.2 mm), 13 October (4.6 mm),

<footnote><sup>a</sup> Cloth tube numbers are given with the first number indicating the distance from the release point and the second number indicating the cloth tube number for that distance. For example, 50.5 is the fifth cloth tube at a distance of 50 m from the release point</footnote>

<footnote><sup>b</sup> 200.8 was not a cloth tube but a clay pot. All other resting sites were cloth tubes</footnote>
<table><caption>**Table 7** Anopheline mosquito species sorted by gender caught from cloth tubes at the central and edge locations during MRRI, Mamfene, Jozini, KwaZulu-Natal Province, South Africa</caption><thead><tr><th rowspan="2">Species identification and gender</th><th colspan="2">Release point and total captured, N (%)</th></tr><tr><th>E</th><th>C</th></tr></thead><tbody><tr><td>Total <i>An. arabiensis</i> males</td><td>83 (42.13)</td><td>27 (50.00)</td></tr><tr><td>Of which were marked <i>An. arabiensis</i> males</td><td>3 (3.61)</td><td>7 [1<sup>a</sup>] (25.92)</td></tr><tr><td><i>An. arabiensis</i> females</td><td>43 (21.83)</td><td>16 (29.63)</td></tr><tr><td><i>An. parensis</i> males</td><td>33 (16.75)</td><td>3 (5.56)</td></tr><tr><td><i>An. parensis</i> females</td><td>23 (11.68)</td><td>2 (3.70)</td></tr><tr><td><i>An. merus</i> males</td><td>1 (0.51)</td><td>0</td></tr><tr><td><i>An. merus</i> females</td><td>1 (0.51)</td><td>2 (3.70)</td></tr><tr><td><i>An. rivulorum</i> males</td><td>0</td><td>0</td></tr><tr><td><i>An. rivulorum</i> females</td><td>1 (0.51)</td><td>2 (3.70)</td></tr><tr><td><i>An. quadriannulatus</i> males</td><td>1 (0.51)</td><td>0</td></tr><tr><td><i>An. quadriannulatus</i> females</td><td>0</td><td>1 (1.85)</td></tr><tr><td>No ID male GC</td><td>3 (1.52)</td><td>0</td></tr><tr><td>No ID female GC</td><td>1 (0.51)</td><td>0</td></tr><tr><td>No ID male FG</td><td>1 (0.51)</td><td>0</td></tr><tr><td>No ID female FG</td><td>2 (1.02)</td><td>1 (1.85)</td></tr><tr><td>Escaped female</td><td>4 (2.03)</td><td>0</td></tr><tr><td>Total</td><td>197 (100.00)</td><td>54 (100.00)</td></tr></tbody></table>

FG, funestus group by morphology, GC, gambiae complex by morphology; N, number of individuals captured,

16 October (6.6 mm), 17 October (5.6 mm), 20 October (3.8 mm) (during the day and evening following the release day). For a summary of the experiments, releases, captures, recaptures and weather data, refer to Additional file 1: Table S1.

## Discussion and conclusion

During this series of MRR experiments the following species were collected from cloth tubes: *An. arabiensis*; *An. merus*; *An. quadriannulatus*; *An. parensis*, *An. rivulorum*; and *An. marshallii* complex. The majority of wild anopheline specimens captured during MRR<sup>PT</sup> (71.4%) and MRRI (63.4% at E., 76.6% at C) were *An. arabiensis*, while only 20% of wild specimens were *An. arabiensis* in MRRII, with the remainder being *An. marshallii* complex and one *An. funestus* group specimen. Sampling in the same area in 2014 and 2015 confirmed the perennial presence of *An. arabiensis*, with the following species collected: *An. arabiensis*; *An. merus*; *An. quadriannulatus*; *An. parensis*; *An. vaneedeni*; *An. leesoni*; and *An. rivulorum* [4]. These results correspond with the species captured during the present series of MRR experiments, except that no *An. vaneedeni* or *An. leesoni* were detected in our experiments. However, the percentages differed markedly in MRRI at E, where a much higher proportion

of *An. parensis* (28.43%) was captured than reported in Dandalo et al. [4] and in the other two MRR experiments. In addition, during MRRII, more *An. marshallii* than any other species were captured. This was the first record of this species since the beginning of the SIT project. The *An. marshallii* complex has been implicated as a potential secondary malaria vector [47] and has become more common in recent field collections. The species diversity we observed is to be expected as the cloth tubes were located outdoors, providing more resting surfaces and possibly being attractive to many different mosquito species. In the experiments performed by Dandalo et al. [4], females were captured more often than males; this was also true in our study for MRR<sup>PT</sup> but not for MRRI and MRRII. Male-specific traps or resting sites have not been successfully developed yet (but see [48]); however, clay pots and the cloth tubes used in this study have previously been found to be effective for the surveillance of both male and female anophelines [4].

The marked male recapture rates in the present series of MRR experiments were low, ranging from 0.04 to 0.11% for recaptures from cloth tubes. Epopa et al. [27] obtained similarly low recapture rates (0.03–0.08%) using clay pots in Burkina Faso during a series of four MRR experiments. In Epopa et al. [27], clay pots were the least effective method at capturing males (recaptures were performed using three methods: indoor pyrethroid spray catches; clay pots and swarms), with extremely low recapture rates (0.03–0.08%) and these data were excluded from their analyses.

Interestingly, we noted that one clay pot used in place of a cloth tube in MRR<sup>PT</sup> captured by far the most mosquitoes. This warrants further investigation into resting site or collection device attractiveness. For example: Are clay pots more attractive than cloth tubes? Does the length of time the device has remained in the field for have an effect on collection efficacy? Does it depend on previous use by mosquitoes and the deposition of pheromones that mark the device as a preferred resting site? Is it simply a positional effect? This clay pot had been in place for an extended period of time, having been used for routine surveillance in the area, suggesting that the length of time a trap or resting site has been in the field may have an influence on the number of mosquitoes captured.

The high number of mosquitoes captured at E compared to C during MRRI indicates that the distribution of mosquitoes throughout the target area is heterogeneous. There is possibly a source of mosquitoes near E, in the marshy area between the intended target site and the nearby control site. Depending on resource availability, it may be better to conduct future pilot trials over the area of the intended target site as well as the nearby

<footnote>a One male captured at C was an orange-marked male (which was released at E)</footnote>
<table><caption>**Table 8** *Anopheles arabiensis* caught per tube distance from each release location during MRRI, Mamfene, Jozini, KwaZulu-Natal Province, South Africa</caption><thead><tr><th>Distance:</th><th colspan="2">25 m</th><th colspan="2">50 m</th><th colspan="2">75 m</th><th colspan="2">100 m</th><th colspan="2">125 m</th><th colspan="2">150 m</th><th colspan="2">Total</th></tr><tr><th>Location:</th><th>C</th><th>E</th><th>C</th><th>E</th><th>C</th><th>E</th><th>C</th><th>E</th><th>C</th><th>E</th><th>C</th><th>E</th><th>C</th><th>E</th></tr></thead><tbody><tr><td>Collection day 1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Marked</td><td>2</td><td>-</td><td>-</td><td>2</td><td>1</td><td>-</td><td>-</td><td>-</td><td>1</td><td>-</td><td>-</td><td>-</td><td>4</td><td>2</td></tr><tr><td>Wild ♂</td><td>2</td><td>3</td><td>1</td><td>1</td><td>2</td><td>1</td><td>-</td><td>2</td><td>-</td><td>2</td><td>1</td><td>-</td><td>6</td><td>9</td></tr><tr><td>Wild ♀</td><td>1</td><td>2</td><td>-</td><td>1</td><td>-</td><td>1</td><td>-</td><td>4</td><td>-</td><td>1</td><td>2</td><td>1</td><td>3</td><td>10</td></tr><tr><td>Collection day 2</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Marked</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td><td>-</td></tr><tr><td>Wild ♂</td><td>1</td><td>1</td><td>1</td><td>7</td><td>-</td><td>3</td><td>1</td><td>4</td><td>-</td><td>1</td><td>1</td><td>4</td><td>4</td><td>20</td></tr><tr><td>Wild ♀</td><td>3</td><td>2</td><td>-</td><td>3</td><td>-</td><td>4</td><td>2</td><td>3</td><td>-</td><td>2</td><td>-</td><td>2</td><td>5</td><td>15</td></tr><tr><td>Collection day 3</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Marked</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td><td>1</td></tr><tr><td>Wild ♂</td><td>1</td><td>5</td><td>1</td><td>9</td><td>2</td><td>2</td><td>-</td><td>4</td><td>-</td><td>7</td><td>-</td><td>3</td><td>4</td><td>30</td></tr><tr><td>Wild ♀</td><td>-</td><td>-</td><td>3</td><td>3</td><td>1</td><td>2</td><td>-</td><td>2</td><td>1</td><td>2</td><td>-</td><td>1</td><td>5</td><td>10</td></tr><tr><td>Collection day 4</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Marked</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Wild ♂</td><td>-</td><td>-</td><td>-</td><td>2</td><td>-</td><td>1</td><td>2</td><td>1</td><td>1</td><td>1</td><td>1</td><td>3</td><td>4</td><td>8</td></tr><tr><td>Wild ♀</td><td>-</td><td>-</td><td>-</td><td>1</td><td>-</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td><td>2</td><td>1</td><td>4</td></tr><tr><td>Collection day 5</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Marked</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1<sup>a</sup></td><td>-</td><td>-</td><td>-</td><td>1</td><td>-</td></tr><tr><td>Wild ♂</td><td>2</td><td>-</td><td>-</td><td>3</td><td>-</td><td>-</td><td>-</td><td>2</td><td>-</td><td>4</td><td>-</td><td>4</td><td>2</td><td>13</td></tr><tr><td>Wild ♀</td><td>1</td><td>1</td><td>1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td><td>2</td><td>2</td></tr><tr><td>Totals</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>43</td><td>124</td></tr><tr><td>Marked</td><td>2</td><td>0</td><td>0</td><td>2</td><td>2</td><td>1</td><td>1</td><td>0</td><td>2<sup>a</sup></td><td>0</td><td>0</td><td>0</td><td>7<sup>a</sup></td><td>3</td></tr><tr><td>Wild ♂</td><td>6</td><td>9</td><td>3</td><td>22</td><td>4</td><td>7</td><td>3</td><td>13</td><td>1</td><td>15</td><td>3</td><td>14</td><td>20</td><td>80</td></tr><tr><td>Wild ♀</td><td>5</td><td>5</td><td>4</td><td>8</td><td>1</td><td>7</td><td>2</td><td>9</td><td>1</td><td>5</td><td>3</td><td>7</td><td>16</td><td>41</td></tr></tbody></table>

control site and then use the second control site as the sole control. Alternatively, the target site may be changed to the distant control site. Careful consideration will be required to ensure that the impact of sterile male releases on the population size is measurable. In addition, more estimates, spaced throughout the year, to estimate seasonal population size fluctuations will be beneficial to optimize the timing of releases (as done in [28] and [49]). This is ideally when the population is at its smallest but still receptive to reproduction. Research indicates that reproduction may take place throughout the year for *An. arabiensis* in the study area [50].

The cold conditions (more than 3 °C below the usual October average), along with rainfall during MRRII, probably contributed to the low numbers of mosquitoes captured during MRRII. Paaijmans et al. [51] found that rains caused loss of larvae from water bodies due to the flushing, ejection and killing effects of rainfall. In addition, Charlwood and Edoh [52] suggest that *An. arabiensis* may be more severely affected by rainfall than *An*.

*gambiae*, possibly because *An. arabiensis* immatures are typically found in small temporary larval sites [9, 10] that are prone to desiccation and flooding.

The average wind speed was significantly less during MRRI than during the other two experiments, which may account for the higher swarming activity observed during MRRI than during MRRPT and MRRII. Swarm searches carried out before the MRR experiments were conducted experienced very little success, particularly on windier days. In addition, studies have indicated that flight and swarming activity are negatively affected by wind. Gillies and Wilkes [53] performed upwind tunnel tests using animal bait (CO₂ and attractive odour source) and found a maximum flight speed for various mosquito species of 1.4–1.8 m/s. At 2 m/s, mosquitoes were no longer captured, instead preferring to settle. Service [54] suggests that wind speeds of ≥ 3 km/h (0.83 m/s) significantly reduce host-seeking flights (female mosquitoes). Munhenga et al. [55] also described a strong correlation between wind speed and mosquito collection in their

<footnote>- = No *An. arabiensis* mosquitoes captured</footnote>

<footnote>a One marked male was captured at C on day 5 was an orange-marked male that had been released at E</footnote>
<table><caption>**Table 9** *Anopheles arabiensis* mosquitoes captured in cloth tubes during MRRII, Mamfene, Jozini, KwaZulu-Natal Province, South Africa</caption><thead><tr><th rowspan="2">Collection day</th><th colspan="6">Trap distance</th><th rowspan="2">Total</th></tr><tr><th>25</th><th>50</th><th>75</th><th>100</th><th>125</th><th>150</th></tr></thead><tbody><tr><td>1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Marked ♂</td><td>2</td><td>-</td><td>2</td><td>-</td><td>-</td><td>-</td><td>4</td></tr><tr><td>Wild ♂</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Wild ♀</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>2</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Marked ♂</td><td>-</td><td>-</td><td>1</td><td>-</td><td>-</td><td>-</td><td>1</td></tr><tr><td>Wild ♂</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>0</td></tr><tr><td>Wild ♀</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>3</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Marked ♂</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>1</td><td>1</td></tr><tr><td>Wild ♂</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Wild ♀</td><td>-</td><td>-</td><td>1</td><td>-</td><td>-</td><td>-</td><td>1</td></tr><tr><td>4</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Marked ♂</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Wild ♂</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Wild ♀</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>5</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Marked ♂</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Wild ♂</td><td>-</td><td>-</td><td>-</td><td>1</td><td>-</td><td>-</td><td>1</td></tr><tr><td>Wild ♀</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Totals</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Marked ♂</td><td>2</td><td>-</td><td>3</td><td>-</td><td>-</td><td>1</td><td>6</td></tr><tr><td>Wild ♂</td><td>-</td><td>-</td><td>-</td><td>1</td><td>-</td><td>-</td><td>1</td></tr><tr><td>Wild ♀</td><td>-</td><td>-</td><td>1</td><td>-</td><td>-</td><td>-</td><td>1</td></tr></tbody></table>

study, with the proportion of *An. gambiae* complex mosquitoes captured decreasing as wind speed increased. Almost no mosquitoes were collected at wind speeds > 3.1 m/s. At higher wind speeds mosquito flight direction would not be controlled by the mosquito, but it would be subject to the direction of the wind. In this situation, it is likely minimal active dispersal takes place [54]. Therefore, despite the presence of attractive resting sites, the mosquitoes may not be able to direct their flight to select such sites and are, instead, swept along by the wind. This would also likely greatly impact swarming behavior of males. Indeed, in a field cage experiment by Achinko et al. [56], there were no swarms detected on rainy and windy days and females did not leave their resting areas, supporting our interpretation.

In previous work related to the South African SIT project, Munhenga et al. [24] determined that sterilized *An. arabiensis* males needed to be released at a minimum ratio of 3:1 to their wild counterparts to be competitive for females. Using the population sizes calculated from the MRR trials in this study and employing a minimum sterilized to wild male release ratio of 6:1, the numbers required per week for inundation

would be approximately 3300 to 37,000 males per hectare using the Lincoln index, which takes daily survival into account [42]. These numbers are based on estimates from MRR<sup>PT</sup> (November/December) and from MRRI. MRRII was not included in this analysis as no wild males were captured during the first 4 days of collections. The fact that MRRII gave a very low population size estimate may be an indication of how low the population can get during unfavourable periods. Periods when the population is reduced, such as just after winter or before the rainy season in South Africa, may be good times to initiate sterile male releases. Targeting the population after other suppression efforts (such as the application of IRS and larviciding of water bodies) at potential points of low population size will increase the inundative ratio and enhance the likelihood of successful target population suppression using the sterile insect technique.

The average distance travelled by marked males captured in cloth tubes was between 58 and 86 m in this series of experiments. Mean distance travelled, as calculated by Ageep et al. [30] for *An. arabiensis*, was 162 m using males captured from swarms. Mean net dispersal, as determined for *An. coluzzii* by Epopa et al. [27], ranged
from 40 to 549 m. Our values fall within the lower ranges reported in these studies. However, mosquitoes are likely able to travel considerably further than the average values obtained, as indicated by the single male recaptured more than 600 m from its release location (and see [57] and [58]). As the cloth tubes in which marked males were recaptured were located at distances of ≤ 150 m from the release points and our recapture rate was low, the dispersal distances are not likely to be representative of actual dispersal, and further studies will be needed to get more accurate results.

The success of the SIT feasibility project hinges on the effective integration of the sterilized laboratory males into the local population in order to compete with wild males for wild females. Laboratory males that can locate and participate in swarms must be at the right place at the right time. During MRRI, we found swarms of wild males that contained 11–50% marked males for swarm sizes of ≥ 2 males. This result is similar to results reported by Ageep et al. [30] who, during their MRR study, recaptured 17.2–42.3% and 26.7–97% marked males from a single large swarm and two small swarms, respectively. Similarly, Epopa et al. [27] captured 76% of marked males from swarms. These results indicate that swarms are probably the ideal locations for obtaining representative samples of males and can be used as traps for population size studies if they are present and accessible.

That laboratory-reared males could locate and partici-
pate in swarms with wild males suggests that they have
retained natural mating behaviours despite a long colo-
nization period. The colonized strain originated from the
area of the study, which may have contributed to their
participation in swarms. Marked laboratory-reared males
have also been captured in known wild male swarms in
Sudan [30, 59]. Epopa et al. [27] used *An. gambiae* (s.l)
males sourced from the wild, as well as colonized *An.
coluzzii* in separate MRR experiments, and found no sig-
nificant differences in the proportions of recaptures from
swarms when these collections were compared. Whether
the mating success of laboratory-reared males in the field
is equal to that of wild males remains unknown and may
be addressed in future experiments. Due to the difficulty
in finding swarms in Mamfene, swarms have not been a
reliable method for sampling mosquitoes in the region.
Several swarms have been observed, but they are not
consistently present.

The marking, transport and release of laboratory-
reared males from Johannesburg was feasible when per-
formed at the scale reported in this study (5000–15,000
males). Challenges for upscaling the number of mos-
quitoes needed for the MRR include: sex-separation of
mosquitoes, which is laborious and time consuming;
marking large numbers of males using a small fridge

for knockdown; and transportation to the field site (live
adults in BugDorm-1® cages occupy a significant amount
of space). Potential ways to mitigate these challenges may
include: the development of a sexing system (a genetic
sexing strain [22] is available); knockdown in a walk-in
cold room for dusting [36] or the use of alternative mark-
ing methods, such as Rhodamine B dye [60]; and trans-
port of marked males in compact containers at reduced
temperatures [61].

We conclude that the population size estimates
obtained in this study can be used as a guide to deter-
mine the initial number of males to be used for an SIT
pilot trial. Based on the lower population size estimates
and a sterile: wild male inundative ratio of 6:1, a minimum
of approximately 165,000 sterile males will need to be
released per week to cover an area of 50 hectares (3300
sterile males per hectare) in the study site. This may
mean that SIT releases should be timed to take place just
before the rainy season (before the population size begins
to increase), or following a drive to suppress the target
population using other methods, such as after IRS cam-
paigns. Due to the heterogeneous nature of the study site
and the close proximity to a potential larval site, adjust-
ment of the target site selected may be required prior to
conducting SIT pilot trials.

Abbreviations

C: Central release area; E: Edge release area; IRS: Indoor residual spraying; KZN: KwaZulu-Natal; MRR: Mark-release-recapture; RH: Relative humidity; SD: Standard deviation; SIT: Sterile insect technique.

**Supplementary Information**

The online version contains supplementary material available at https://doi.org/10.1186/s13071-021-04674-w.

**Additional file 1: Table S1.** Mamfene, KwaZulu-Natal MRR experiments summary information. **Table S2.** Wind direction and speed averages for the pilot trial (MRR<sup>P1</sup>), MRRI and MRRII.

**Additional file 2.** Sweep net design. **A** Sweep net. Net length: ~ 1.5 m, hoop diameter: ~ 70 cm, handle: ~ 2 m. **B** Mid notch and lashing. **C** Top notch and lashing.

Acknowledgements

The following people are sincerely thanked for providing assistance with fieldwork and/or the preparation and processing of mosquitoes: Jabulani Zikhali, Dumsani Dlamini, Leanne Lobb, Nelius Venter, Erica Erlank, Michael Samuel, Ashley Burke, Jacek Zawada, Power Tshikae, Nonhlanhla Ntoyi, Eunice Jamesboy, Nondumiso Mabaso, Fortunate Moletsane, Malibongwe Zulu. Bheki Qwabe and Zuziwe Nyawo are thanked for the use of the KZN malaria offices and insectary during the field work. The DST team are thanked for their assistance during the MRRII fieldtrip. Hamidou Maïga is thanked for providing the initial training on swarm finding. Hanano Yamada, Jeremie Gilles and Mark Benedict are sincerely thanked for providing advice and support. Theresa Mazarire is thanked for ably producing the maps for Figs. 1, 3, 4 and 5. The Makhathini Research Station is thanked for providing the weather data. Thank you to two anonymous reviewers for helpful comments and suggestions to improve the manuscript. Finally, a big thank you to the Mamfene community.
## Authors' contributions

MLK designed, co-ordinated and performed the experiments, identified mosquitoes by morphology and drafted the manuscript. ORW produced the mosquitoes, assisted in performing experiments, morphologically identified mosquitoes and assisted with drafting of the manuscript. DD assisted in experimental design, conducting experiments and drafting of the manuscript. LLK and GM conceived the study and assisted with experimental design and drafting of the manuscript. BDB assisted with the drafting of the manuscript. All authors read and approved the final manuscript.

## Funding

This project was funded by the International Atomic Energy Agency through their technical cooperation programme (SAF5014) and Research Contract No. 17904; 17953; 19099, the Department of Science and Technology through their Health Innovation Scheme. The study was also partly funded by a Global Diseases Detection/CDC Grant (U19GH000622-01 MAL01) and the National Research Foundation (NRF) through the Medical Faculty Research Endowment Fund (00140184612015121105RMVIRL015) and a DST/NRF South African Research Chairs Initiative grant to Prof Maureen Coetzee and NRF grants (Grant numbers 119765 and 107428) awarded to Dr. Givemore Munhenga. The ANTI-VeC Pump-Priming Grant AV/PP0027/1 from University of Glasgow is also acknowledged for support.

## Availability of data and materials

The dataset generated by this study is available from the corresponding author upon reasonable request.

## Declarations

### Ethics approval and consent to participate

This study was approved by the KwaZulu-Natal Department of Health Research Committee (Ref. HRKM337/14) and the University of the Witwatersrand Human Research Ethics Committee (Clearance Certificate No. M141023). Oral informed consent was obtained from all household owners, after being provided with all necessary information regarding goals, objectives, benefits as well as risks associated with the study. Permission was obtained from community leaders prior to conducting the experiments.

### Consent for publication

Not applicable.

### Competing interests

The authors declare no competing interests.

### Author details

¹ Centre for Emerging Zoonotic & Parasitic Diseases, National Institute for Communicable Diseases of the National Health Laboratory Service, Johannesburg, South Africa. ² Wits Research Institute for Malaria, School of Pathology, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa. ³ UMR IRD 224, Maladies Infectieuses et Vecteurs-Ecologie-Génétique, Evolution et Contrôle (MIVEGEC), Institut de Recherche Pour Le Développement (IRD) CNRS 5290-Université de Montpellier, Montpellier, France. ⁴ IRD Réunion/GIP CYROI (Recherche Santé Bio-innovation), Sainte Clotilde, Réunion Island, France.

Received: 29 December 2020 Accepted: 7 March 2021

Published online: 19 April 2021

## References

1. World Health Organization. World malaria report 2019. Geneva: World Health Organization. 2019.

2. Malaria risk map for South Africa. http://www.nicd.ac.za/wp-content/uploads/2018/12/south_africa_malaria_risk_dec2018_final.pdf. Accessed 27 Nov 2019.

3. Brooke B, Koekemoer L, Kruger P, Urbach J, Misiani E, Coetzee M. Malaria vector control in South Africa. S Afr Med J. 2013;103(10 Suppl 2):784-8.

4. Dandalo LC, Brooke BD, Munhenga G, Lobb LN, Zikhali J, Ngxongo SP, et al. Population dynamics and *Plasmodium falciparum* (Heamosporida: Plasmodiidae) infectivity rates for the malaria vector *Anopheles arabiensis* (Diptera: Culicidae) at Mamfene, KwaZulu-Natal. S Afr J Med Entomol. 2017;54(6):1758-66.

5. Christian R, Dahan-Moss D, Braack L, Munhenga G, Kaiser M, Lobb L, et al. Malaria vector surveillance report, South Africa, January-December, 2017. NICD Public Health Surveillance Bull. 2018;16(1):29-35.

6. Burke A, Dandalo L, Munhenga G, Dahan-Moss D, Mbokazi F, Ngxongo S, et al. A new malaria vector mosquito in South Africa. Sci Rep. 2017;7:43779.

7. Burke A, Dahan-Moss D, Duncan F, Qwabe B, Coetzee M, Koekemoer L, et al. *Anopheles parensis* contributes to residual malaria transmission in South Africa. Malar J. 2019;18:257.

8. Gillies MT, De Meillon B. The Anophelinae of Africa south of the Sahara. 2nd ed. Johannesburg: South African Institute for Medical Research; 1968.

9. Gillies MT, Coetzee M. A supplement to the anophelinae of Africa south of the Sahara. Johannesburg: The South African Institute for Medical Research. 1987.

10. Sinka ME, Bangs MJ, Manguin S, Coetzee M, Mbogo CM, Hemingway J, et al. The dominant *Anopheles* vectors of human malaria in Africa, Europe and the Middle East: occurrence data, distribution maps and bionomic précis. Parasites Vectors. 2010;3:117.

11. South Africa National Department of Health. Malaria elimination strategic plan for South Africa 2019-2023. Pretoria: South Africa National Department of Health; 2019.

12. Coetzee M, Kruger P, Hunt RH, Durrheim DN, Urbach J, Hansford CF. Malaria in South Africa: 110 years of learning to control the disease. S Afr Med J. 2013;103(10 Suppl 2):770-8.

13. Knipliing EF. Possibilities of insect control or eradication through the use of sexually sterile males. J Econ Entomol. 1985;48(4):459-562.

14. Benedict MQ, Robinson AS. The first releases of transgenic mosquitoes: an argument for the sterile insect technique. Trends Parasitol. 2003;19(8):349-55.

15. Wyss JH. Screwworm eradication in the Americas. Ann N Y Acad Sci. 2000;916:186-93.

16. Bloem KA, Bloem S, Carpenter JE. Impact of moth suppression/eradication programmes using the sterile insect technique or inherited sterility. In: Dyck VA, Henrichs J, Robinson AS, editors. Sterile insect technique. Principles and practice in area-wide integrated pest management. Springer: Dordrecht; 2005. p. 677-700.

17. Enkerlin W, Gutiérrez-Ruelas JM, Cortes AV, Roldan EC, Midgarden D, Lira E, et al. Area freedom in Mexico from Mediterranean Fruit Fly (Diptera: Tephritidae): a review of over 30 years of a successful containment program using an integrated area-wide SIT approach. Fla Entomol. 2015;98(2):665-81.

18. Vreysen MJB, Saleh KM, Ali MY, Abdulla AM, Zhu Z-R, Juma KG, et al. *Glossina austeni* (Diptera: Glossinidae) eradicated on the island of Unguja, Zanzibar, using the sterile insect technique. J Econ Entomol. 2000;93(1):123-35.

19. Bloem S, McCuskey A, Fugger R, Arthur S, Wood S, Carpenter J. Suppression of the codling moth *Cydia pomonella* in British Columbia, Canada using an area-wide integrated approach with an SIT component. In: Vreysen MJB, Robinson AS, Hendrichs J, editors. Area-wide control of insect pests. Dordrecht: Springer; 2007. p. 591-601.

20. Munhenga G, Brooke BD, Chirwa TF, Hunt RH, Coetzee M, Govender D, et al. Evaluating the potential of the sterile insect technique for malaria control: relative fitness and mating compatibility between laboratory colonized and a wild population of *Anopheles arabiensis* from the Kruger National Park, South Africa. Parasites Vectors. 2011;4:208.

21. Manana PH, Kuonza L, Musekiwa A, Mpangane HD, Koekemoer LL. Knowledge attitude and practices on malaria transmission in Mamfene, KwaZulu-Natal Province, South Africa 2015. BMC Public Health. 2017;18:41.

22. Dandalo LC, Munhenga G, Kaiser ML, Koekemoer LL. Development of a genetic sexing strain of *Anopheles arabiensis* for KwaZulu-Natal, South Africa. Med Vet Entomol. 2018;32(1):61-9.

23. Dandalo LC. Characterization of a local genetic sexing strain as well as a wild population of *Anopheles arabiensis* from KwaZulu Natal, South Africa. 2017. PhD thesis. University of the Witwatersrand.

24. Munhenga G, Brooke BD, Gilles JRL, Slabbert K, Kemp A, Dandalo LC, et al. Mating competitiveness of sterile genetic sexing strain males (GAMA) under laboratory and semi-field conditions: steps towards the use of the
Sterile Insect Technique to control the major malaria vector *Anopheles arabiensis* in South Africa. *Parasites Vectors*. 2016;9:122.

25. Guerra CA, Reiner RC, Perkins TA, Lindsay SW, Midega JT, Brady OJ, et al. A global assembly of adult female mosquito mark-release-recapture data to inform the control of mosquito borne pathogens. *Parasites Vectors*. 2014;7:276.

26. Morris CD, Larson VL, Lounibos LP. Measuring mosquito dispersal for control programs. J Am Mosq Control Assoc. 1991;7:608–15.

27. Epopa PS, Millogo AA, Collins CM, North A, Tripet F, Benedict MQ, et al. The use of sequential mark-release-recapture experiments to estimate population size, survival and dispersal of male mosquitoes of the *Anopheles gambiae* complex in Bana, a west African humid savannah village. *Parasites Vectors*. 2017;10:376.

28. Le Goff G, Damiens D, Ruttee A-H, Payet L, Lebon C, Dehecq J-S, et al. Field evaluation of seasonal trends in relative population sizes and dispersal patterns of *Aedes albopictus* males in support of the design of a sterile male release strategy. *Parasites Vectors*. 2019;12:81.

29. Cianci D, Van den Broek J, Caputo B, Marini F, Della Torre A, Heesterbeek H, et al. Estimating mosquito population size from mark-release-recapture data. J Med Entomol. 2013;50(3):533–42.

30. Ageep TB, Damiens D, Alsharif B, Ahmed A, Salih EHO, Ahmed FTA, et al. Participation of irradiated *Anopheles arabiensis* males in swarms following field release in Sudan. Malar J. 2014;13:484.

31. Gouagna LC, Rakotondranary M, Boyer S, Lempérière G, Dehecq J-S, Fontenille D. Abiotic and biotic factors associated with the presence of *Anopheles arabiensis* immatures and their abundance in naturally occurring and man-made aquatic habitats. *Parasit Vectors*. 2012;5:96.

32. Climate-data.org. https://en.climate-data.org/africa/southAfrica/kwazu-lu-natal/jozini. Accessed 2 May 2019.

33. Mouatcho JC, Munhenga G, Hargreaves K, Brooke BD, Coetzee M, Koekemoer LL. Pyrethroid resistance in a major African malaria vector *Anopheles arabiensis* from Mamfene, northern KwaZulu-Natal South Africa. S Afr J Sci. 2009;1005:127–31.

34. Hunt RH, Brooke BD, Pillay C, Koekemoer LL, Coetzee M. Laboratory selection for and characteristics of pyrethroid resistance in the malaria vector *Anopheles funestus*. Med Vet Entomol. 2005;19(1):271–5.

35. Food and Agriculture Organization of the UN/International Atomic Energy Agency (FAO/IAEA). Guidelines for standardized mass rearing for *Anopheles* mosquitoes. Version 1.0. 2017. http://www.naweb.iaea.org/nafa/ipc/public/manuals-ipc.html. Accessed 20 July 2020.

36. Culbert NJ, Kaiser M, Venter N, Vreysen MJB, Gilles JRL, Bouyer J. A standardized method of marking male mosquitoes with fluorescent dust. *Parasites Vectors*. 2020;13:192.

37. Diabaté A, Dao A, Yaro AS, Adamou A, Gonzalez R, Manoukis NC, et al. Spatial swarm segregation and reproductive isolation between the molecular forms of *Anopheles gambiae*. Proc R Soc B. 2009;276:4215–22.

38. Jones JA. The swarming and mating flight of diptera. Ann Rev Entomol. 1969;14:271–98.

39. Sawadogo SP, Costantini C, Pennetier C, Diabaté A, Gibson G, Dabiré RK. Differences in timing of mating swarms in sympatric populations of *Anopheles coluzzii* and *Anopheles gambiae* s.s. (formerly *An. gambiae* M and S molecular forms) in Burkina Faso, West Africa. *Parasites Vectors*. 2013;6:275.

40. Scott JA, Brogdon WG, Collins FH. Identification of single specimens of the *Anopheles gambiae* complex by the polymerase chain reaction. Am J Trop Med Hyg. 1993;49(4):520–9.

41. Koekemoer LL, Kamau L, Hunt RH, Coetzee M. A cocktail polymerase chain reaction assay to identify members of the *Anopheles funestus* (Diptera: Culicidae) group. Am J Trop Med Hyg. 2002;66(6):804–11.

42. Yasuno M, Rajagopalan PK. Population estimation of *Culex fatigans* in Dehli villages. WHO/VBC 1973;73:431:18. Geneva: World Health Organization.

43. Service MW. Chapter 9 Mark-recapture techniques and adult dispersal. In *mosquito ecology: field sampling methods*. London: Applied Science Publishers LTD; 1976. p. 424–37.

44. Charlwood JD, Graves PM, Birley MH. Capture-recapture studies with mosquitoes of the group of *Anopheles punctulatus* Donitz (Diptera: Culicidae) from Papua New Guinea. Bull Entomol Res. 1986;76:211–27.

45. Walker ED, Copeland R, Paulson SL, Munstermann LE. Adult survivorship, population density, and body size of *Aedes triseriatus* and *Aedes hender-soni* (Diptera: Culicidae). J Med Entomol. 1987;24(4):485–93.

46. Niebylski ML, Craig GB Jr. Dispersal and survival of *Aedes albopictus* at a scrap tire yard in Missouri. J Am Mosq Contr. 1994;10(3):339–43.

47. Bafort JM. *Anopheles marshallii* s.l., a secondary vector of malaria in Africa. Trans R Soc Trop Med Hyg. 1985;79:566.

48. Diabate A, Tripet F. Targeting male mosquito mating behavior for malaria control. *Parasites Vectors*. 2015;8:347.

49. Gouagna L-C, Dehecq J-S, Fontenille D, Dumont Y, Boyer S. Seasonal variation in size estimates of *Aedes albopictus* population based on standard mark-release-recapture experiments in an urban area on Reunion Island. *Acta Trop*. 2015;143:89–96.

50. Burke AM, Brooke BD, Duncan FD. Metabolic rate does not vary with seasonal change in *Anopheles arabiensis* adults in South Africa. J Insect Physiol. 2019;118:103942.

51. Paaijmans KP, Wandago MO, Githeko AK, Takken W. Unexpected high losses of *Anopheles gambiae* due to rainfall. PLoS ONE. 2007;2(11):e1146.

52. Charlwood JD, Edoh D. Polymerase chain reaction used to describe larval habitat use by *Anopheles gambiae* complex (Diptera: Culicidae) in the environs of Ifakara, Tanzania. J Med Entomol. 1996;33(2):202–4.

53. Gillies MT, Wilkes TJ. Field experiments with a wind tunnel on the flight speed of some West African mosquitoes (Diptera: Culicidae). Bull Ent Res. 1981;71:65–70.

54. Service, MW. Effects of wind on the behavior and distribution of mosquitoes and blackflies. Int J Biometeorol. 1980;24:347–53.

55. Munhenga G, Brooke BD, Spillings B, Essop L, Hunt RH, Midzi S, et al. Field site selection, species abundance and monthly distribution of anophe-line mosquitoes in the northern Kruger National Park, South Africa. Malar J. 2014;13:27.

56. Achinko D, Thailayil J, Paton D, Mireji PO, Tolesa V, Masiga D, et al. Swarming and mating behaviour of *Anopheles gambiae* mosquitoes in semi-field enclosures. Med Vet Entomol. 2016;30:14–20.

57. Verdonschot FM, Besse-Lototskaya AA. Flight distance of mosquitoes (Culicidae): a metadata analysis to support the management of barrier zones around rewetted and newly constructed wetlands. Limnologica. 2014;45:69–79.

58. Huestis DL, Dao A, Diallo M, Sanogo ZL, Samake D, Yaro AS, et al. Windborne long-distance migration of malaria mosquitoes in the Sahel. Nature. 2019;574(7778):404–8.

59. Hassan MM, Zain HM, Basheer ME, Elhaj H-EF, El-Sayed BB. Swarming and mating behavior of male *Anopheles arabiensis* Patton (Diptera: Culicidae) in an area of the Sterile Insect Technique Project in Dongola, northern Sudan. *Acta Trop*. 2014;132:564–9.

60. Johnson BJ, Mitchell SN, Paton CJ, Stevenson J, Staunton KM, Snoad N, et al. Use of Rhodamine B to mark the body and seminal fluid of male *Aedes aegypti* for mark-release-recapture experiments and estimating the efficacy of sterile male releases. PLoS Negl Trop Dis. 2017;11(9):e0005902.

61. Culbert NJ, Lees RS, Vreysen MJB, Darby AC, Gilles JRL. Optimised conditions for handling and transport of male *Anopheles arabiensis*: effects of low temperature, compaction, and ventilation on male quality. Entomol Exp Appl. 2017;164(3):276–83.

**Publisher's Note**

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.