Hostname: page-component-7b9c58cd5d-g9frx Total loading time: 0 Render date: 2025-03-15T20:20:34.015Z Has data issue: false hasContentIssue false

Modelling factors that affect the presence of larval mosquitoes (Diptera: Culicidae) in stormwater drainage systems to improve the efficacy of control programmes

Published online by Cambridge University Press:  10 September 2013

Michael J. Jackson*
Affiliation:
Culex Environmental Ltd., 4-4075 Kingsway Avenue, Burnaby, British Columbia, Canada V5H 1Y9
Jennifer L. Gow
Affiliation:
Culex Environmental Ltd., 4-4075 Kingsway Avenue, Burnaby, British Columbia, Canada V5H 1Y9
Michelle J. Evelyn
Affiliation:
Culex Environmental Ltd., 4-4075 Kingsway Avenue, Burnaby, British Columbia, Canada V5H 1Y9
T.J. Scott McMahon
Affiliation:
Culex Environmental Ltd., 4-4075 Kingsway Avenue, Burnaby, British Columbia, Canada V5H 1Y9
Harlan Campbell
Affiliation:
Culex Environmental Ltd., 4-4075 Kingsway Avenue, Burnaby, British Columbia, Canada V5H 1Y9
Jennifer Sheppard
Affiliation:
Culex Environmental Ltd., 4-4075 Kingsway Avenue, Burnaby, British Columbia, Canada V5H 1Y9
Tim J. Howay
Affiliation:
Culex Environmental Ltd., 4-4075 Kingsway Avenue, Burnaby, British Columbia, Canada V5H 1Y9
Disa Fladmark
Affiliation:
Culex Environmental Ltd., 4-4075 Kingsway Avenue, Burnaby, British Columbia, Canada V5H 1Y9
Aynsley Thielman
Affiliation:
Culex Environmental Ltd., 4-4075 Kingsway Avenue, Burnaby, British Columbia, Canada V5H 1Y9
*
1Corresponding author (e-mail: m.jackson@culex.ca).
Rights & Permissions [Opens in a new window]

Abstract

Stormwater catch basins form part of artificial drainage systems in urban areas and can provide larval habitat for mosquito vector species of West Nile virus (WNv), such as Culex pipiens Linnaeus (Diptera: Culicidae). We evaluated the impact of management techniques and targeted applications of larvicide on larval populations of this potential WNv mosquito vector species in catch basins from the Lower Mainland of Vancouver and on Vancouver Island of British Columbia, Canada. A mixed effects logistic regression model described the relationship between larval presence and larvicide treatment while controlling for other parameters. Parameter estimates showed that larvicide treatment reduced the odds of larvae presence by a factor of ∼7.23. The model also revealed relationships between larval presence and water temperature and adjacent land use but larvicide treatment consistently reduced the presence of larvae regardless of these other factors. This knowledge can now be used to prioritise and target control efforts to most efficiently reduce WNv mosquito vector populations, and most effectively reduce the risk of WNv transmission to humans. A similar research strategy could be applied to emerging threats from other potential mosquito vectors of disease around the world, to help lower the incidence of mosquito-borne disease.

Résumé

Les réseaux artificiels de drainage dans les zones urbaines comprennent des bassins de captage des eaux pluviales qui peuvent servir d'habitats pour les larves de moustiques vecteurs du virus du Nil occidental (WNv), tels que Culex pipiens Linnaeus (Diptera: Culicidae). Nous évaluons les impacts des techniques d'aménagement et des épandages ciblés de larvicide sur les populations de larves de cette espèce, un moustique vecteur potentiel du WNv, dans des bassins de captage dans les basses terres de la région de Vancouver et sur l’île de Vancouver en Colombie-Britannique, Canada. Un modèle de régression logistique à effets mixtes permet de décrire la relation entre la présence de larves et le traitement au larvicide, tout en tenant compte des autres variables. Les estimations des variables montrent que le traitement au larvicide réduit la probabilité de présence de larves par un facteur de l'ordre de 7,23. Le modèle montre aussi une relation entre la présence de larves, d'une part, et la température de l'eau et l'utilisation des terres adjacentes, d'autre part; cependant, le traitement au larvicide réduit toujours la présence des larves, quels que soient les autres facteurs. Ces informations peuvent servir à établir des priorités et fixer des cibles dans les programmes de lutte pour réduire le plus efficacement possible les populations de moustiques vecteurs du WNv et diminuer ainsi le risque de transmission du WNv aux humains. On pourrait utiliser une stratégie de recherche similaire pour étudier les nouvelles menaces que représentent d'autres moustiques vecteurs potentiels de maladies sur la planète afin de diminuer l'incidence des maladies transmises par les moustiques.

Type
Insect Management
Copyright
Copyright © Entomological Society of Canada 2013 

Introduction

Since its introduction to North America in 1999, West Nile virus (WNv) has spread rapidly across the continent with dire public health consequences. Between 1999 and 2010, ∼1.8 million people have been infected, causing about 360 000 illnesses, 12 852 reported cases of encephalitis/meningitis and 1308 deaths (Kilpatrick Reference Kilpatrick2011). The virus is now also a recognised public health threat in Mexico, and Central and South America (Komar and Clark Reference Komar and Clark2006).

The WNv Flaviviridae virus has an enzootic life cycle between birds and mosquitoes. Mosquitoes (Diptera: Culicidae) that bite both WNv bird hosts and humans act as bridge vectors that have the capacity to transfer WNv from infected birds to human populations, causing illness. In Canada, this applies to <20 of the more than 80 known species of mosquito (Thielman and Hunter Reference Thielman and Hunter2007). While the role of many of these potential vectors in WNv epidaemiology is currently poorly understood, exceptions are Culex pipiens Linnaeus (Diptera: Culicidae) and Culex tarsalis Coquillett (Diptera: Culicidae); species considered to be potentially highly competent WNv vectors (Belton Reference Belton1983, Reference Belton2006; Turell et al. Reference Turell, Dohm, Sardelis, O'Guinn, Andreadis and Blow2005).

To combat the public health risk of WNv, the most successful proactive programmes employ mosquito larval control to maintain mosquito abundance below thresholds above which transmission to humans most commonly occurs (Reisen and Brault Reference Reisen and Brault2007). Larviciding is the preferred option here as it can achieve more reliable mosquito control than adulticiding (Nielsen et al. Reference Nielsen, Reisen, Armijos, Wheeler, Kelley and Brown2007; Su Reference Su2008). However, successful mosquito control programmes require detailed surveillance of mosquito abundance to operate successfully (Reisen and Brault Reference Reisen and Brault2007). Larval habitats of mosquitoes that require monitoring in such programmes range from natural habitats, such as floodplains, and ponded streams, to artificial habitats, like ditches, ponds, and drainage systems. Artificial habitats that place WNv mosquito vectors in especially close proximity to human populations are of particular concern, since they have the potential to greatly facilitate the spread of the disease. For example, stormwater management systems have been linked to the spread of WNv (Vazquez-Prokopec et al. Reference Vazquez-Prokopec, Vanden Eng, Kelly, Mead, Kolhe and Howgate2010).

In most urban areas, these stormwater systems channel surface water from roadways and other thoroughfares via catch basins (Fig. 1), which are typical breeding sites for mosquitoes. A grate usually covers a catch basin where surface runoff water enters before falling into a sump. Water can then flow out of the catch basin through a pipe that is set at the top of the sump (Fig. 1). While this design allows settlement and retention of unwanted sediment, debris and pollutants, it also provides ideal larval habitat for species of mosquito that prefer stagnant water that is high in organic matter (Harbison et al. Reference Harbison, Metzger and Hu2010). This design of catch basin can also provide optimal water temperatures (WT) for mosquito larval development for longer periods of time compared with most naturally occurring water bodies (Harbison et al. Reference Harbison, Metzger, Walton and Hu2009b).

Fig. 1 A schematic cross-section through a typical sumped catch basin. Such basins commonly provide habitat for mosquito larvae.

Mosquito species that thrive in such conditions are often highly competent WNv vectors, such as C. pipiens in northeastern and northwestern North America, Culex quinquefasciatus Say (Diptera: Culicidae) across the southern half of the continent, and Culex restuans Theobald (Diptera: Culicidae) in central and eastern portions of North America (Thomson Reference Thomson2008; Jackson et al. Reference Jackson, Gow, Evelyn, Meikleham, McMahon and Koga2009; Harbison et al. Reference Harbison, Metzger and Hu2010). As well as increasing the amount of habitat for WNv mosquito vectors within their native ranges, catch basins can also facilitate the invasion of WNv mosquito vectors into new areas (Norris Reference Norris2004; Morris et al. Reference Morris, Lampman, Ballmes, Funes, Halvorsen and Novak2007). For example, in southwestern Canada, the introduction of C. pipiens coincided with the construction of stormwater drainage systems with catch basins (Vinogradova Reference Vinogradova2000).

In this paper, we investigate the impact of management techniques and targeted applications of larvicide on larval populations of potential WNv mosquito vector species in catch basins. To do this, we examine the influence that larvicide treatment has on larval presence in catch basins, while simultaneously exploring several other parameters (seasonal variation in larval abundance, vegetation type of immediately surrounding area, WT and weather). To do this we use seven years of sampling data from catch basins from our Integrated Mosquito Management Plan of Action (IMMPACT®) programme in the Lower Mainland and on Vancouver Island, British Columbia, Canada.

Materials and methods

Catch basin study area and sampling

Data were analysed from 10 407 sampling events from 2000 unique catch basins over a seven-year period (2004–2010) in Metro Vancouver, Vancouver Island, the Sunshine Coast, and along the Sea-to-Sky Highway in southwestern British Columbia, Canada (Fig. 2). For each municipality or district that was sampled within these areas, catch basins were selected to achieve a relatively even distribution, taking into account accessibility and safety considerations.

Fig. 2 Study area showing locations of catch basins monitored in the Lower Mainland and Vancouver Island, British Columbia, Canada from 2004 to 2010.

The primary objective of this sampling was to monitor the abundance of mosquito larvae in catch basins throughout the mosquito breeding season between late May and late September each year. The total number of sampling events per year ranged from 215 in 2009 to 3078 in 2008 (Fig. 3), with an average of 2.7 sampling events per catch basin each year (ranged from 2.2 in 2004 to 3.2 in 2005). Eighty eight percent (n = 9150) of sampling events occurred prior to any larvicide treatment within any given year (termed “pre-treatment” samples from herein). Indeed, no larvicide treatment of catch basins took place in some years (2004, 2005, 2009, Fig. 3). In other years when funding was available, however, some districts and municipalities treated all of their wet catch basins with a single application of larvicide (2006, 2007, 2008, 2010, Fig. 3) when threshold numbers of WNv vector larvae were exceeded (more than 10% of wet catch basins contained 10 or more late instar larvae). These wide-scale treatment efforts involved every wet catch basin within a district or municipality receiving a water soluble pouch containing 10 g of the bacterial larvicide VectoLex® WSP, Valent Biosciences, Libertyville Illinois, United States of America (active ingredient Bacillus sphaericus) during July or August. Treatment took place, whenever possible, after the annual cleaning of catch basins to avoid removal of larvicide by cleaning machines. The 12% (n = 1257) of sampling events that occurred after such treatment (termed “post-treatment” samples from herein) evaluated the efficacy of bacterial larviciding campaigns.

Fig. 3 The mean number of mosquito larvae per catch basin sample according to year (2004–2010). Error bars show standard errors. The percentage of samples occurring post-treatment is in brackets above bars. The total number of sampling events is below each year.

Mosquito larvae and pupae were sampled from each catch basin using an aquarium net attached to a long wooden pole (Fig. 1). To standardise sample collection, this was placed down a catch basin's grate, and moved in a figure eight pattern before removal. The net was then flushed with clean water into a white basin to maximise visibility of captured mosquitoes. The larvae and pupae were then transferred in water to Ziploc® bags (SC Johnson, Racine, Wisconsin, United States of America) and transported to the laboratory in a cooler. Other parameters that were also collected at time of sampling included: collection date, WT, and vegetation type of the immediately surrounding area. Weather station archives for Vancouver from Environment Canada were used to verify prevailing weather conditions recorded at time of sampling by field technicians. The predominant vegetation type of the street area immediately surrounding each catch basin (<10 m) was classified into one of three categories: if the catch basin likely received leaves, the classification was tree-lined; if grass clippings probably predominated, then it was grass; and, if neither of these were the case, then it was pavement.

Larval sample processing

From 2004 to 2006, all larvae and pupae collected were identified to the lowest possible taxonomic level, given the condition of the specimen, using the keys outlined in Belton (Reference Belton1983) and Wood et al. (Reference Wood, Dang and Ellis1979). Larval specimens were categorised as either “early” (1st or 2nd) or “late” (3rd or 4th) instars. Early instars were reared into late instars for reliable identification in containers with ground up fish flakes. Late instars were preserved for identification in 75% ethanol. Pupae were placed in emergence cages and their subsequent emerged adults identified (and included in the larval count).

Statistical analysis

The extremely long-tailed zero-heavy distribution of the larval counts presented difficulty for standard analysis. To overcome this, we considered a simpler binary presence/absence response variable. A mixed effects logistic regression model was then fit to describe the relationship between larval presence and larvicide treatment while controlling for the other parameters that were collected (seasonal variation in larval abundance, vegetation type of immediately surrounding area, WT, and weather), as well as variability across sites and municipalities. Seasonal variation in larval abundance was included in the model using a cubic function of day of the year, d, (d/100 + d 2/104 + d 3/106). This generalised linear mixed model for a binary response not only addressed the zero-heavy non-normal long-tailed distribution of the larval counts, but also dealt with repeated measures, and the unbalanced nature of the data (pre-treatment versus post-treatment samples, Bolker et al. Reference Bolker, Brooks, Clark, Geange, Poulsen and Stevens2009). The following model formulation relates the log-odds of larval presence (y) to the set of covariates considered, namely larvicide treatment (T), seasonal variation in larval abundance (D 1 = (d)/(102), D 2 = (d 2)/(104), D 3 = (d 3)/(106); where d is day of year), surroundings (X 1 = 1 if “paved”, 0 otherwise; X 2 = 1 if “treeline”, 0 otherwise; X 1 = X 2 = 0, if “grass”), WT and weather (W 1 = 1 if “hard rain”, 0 otherwise; W 2 = 1 if “light rain”, 0 otherwise; W 3 = 1 if “overcast”, 0 otherwise; W 4 = 1 if “partly cloudy”, 0 otherwise; W 1 = W 2 = W 3 = W 4 = 0 if “sunny”):

$${{{y}_{ijk}}\, = \,{{{\rm{\beta }}}_{\rm{0}}}\, + \,{{{\rm{\beta }}}_{\rm{1}}}{{T}_{ijk}}\, + \,{{{\rm{\beta }}}_{\rm{2}}}{{D}_{{\rm{1}}ijk}}\, + \,{{{\rm{\beta }}}_{\rm{3}}}{{D}_{{\rm{2}}ijk}}\, \cr + \,{{{\rm{\beta }}}_{\rm{4}}}{{D}_{{\rm{3}}ijk}}\, + \,{{{\rm{\beta }}}_{\rm{5}}}{{X}_{{\rm{1}}jk}}\, + \,{{{\rm{\beta }}}_{\rm{6}}}{{X}_{{\rm{2}}jk}}\, + \, {{{\rm{\beta }}}_{\rm{7}}}W{{T}_{ijk}}\, \cr + \,{{{\rm{\beta }}}_{\rm{8}}}{{W}_{{\rm{1}}ijk}}\, + \,{{{\rm{\beta }}}_{\rm{9}}}{{W}_{{\rm{2}}ijk}}\, + \,{{{\rm{\beta }}}_{{\rm{10}}}}{{W}_{{\rm{3}}ijk}}\, \cr + \,{{{\rm{\beta }}}_{{\rm{11}}}}{{W}_{{\rm{4}}ijk}}\, + \,{{{\rm{\gamma }}}_{\rm{k}}}\, + \,{{u}_{jk}}\, + \,{{ \in }_{ijk}}}$$

In this notation, the subscripts identify that yijk is the ith observation made at the jth site within the kth municipality. The fixed effects of the model are larvicide treatment, seasonal variation in larval abundance, vegetation type of immediately surrounding area, and WT. Since quantifying the particularities of any given site or municipality was not the focus in this broader perspective study, these variables are included as normally distributed, random effects and modeled such that:

$${{{\rm{\gamma }}}_k}\, \sim \,N\left( {0,\,{{{\rm{\sigma }}}_a}} \right),\; \; {{u}_{jk}}\, \sim \,N\left( {0,{{{\rm{\sigma }}}_b}} \right)$$

This type of multi-level model is common in ecology when multiple responses are measured per individual (in this case individual catch basin) and per region (in this case municipality; Bolker et al. Reference Bolker, Brooks, Clark, Geange, Poulsen and Stevens2009).

Results

In total, 92 633 mosquito larvae were collected from 10 407 sampling events from 2000 wet catch basins from 2004 to 2010. The majority (70.0%, n = 7276) of sampling events yielded zero larvae. The remaining non-zero counts (30.0%, n = 3131) exhibited a skewed distribution, with a mean of 29.5, a median of 9, and a maximum of 2328 larvae.

These data indicated that catch basins were consistently dominated by C. pipiens (89%, 71 597 of the 80 668 mosquito larvae that were identified to species). The remaining 11% of the larvae were made up exclusively of Culiseta incidens Thomson (Diptera: Culicidae), another species that potentially transmits WNv (Belton Reference Belton2006). As a result of this uniformity in species composition, all larvae collected after 2006 were counted, but identification was limited to random identification checks for confirmation only.

The vast majority of larvae (98.4%, n = 91 181) were collected from pre-treatment catch basin samples; only 1.59% of larvae (n = 1452) were found in post-treatment samples. Indeed, parameter estimates from the model suggest that larvicide treatment has the effect of reducing the odds of larvae presence by a factor of ∼7.23 (β1 [treatment] = −1.98, P-value ≤0.001; Appendix 1). To explore the treatment effect further, we also ran the model omitting the years during which there was no treatment (2004, 2005, and 2009) from the data. Parameter estimates remained similar. Most importantly, the latter showed treatment has the effect of reducing the odds of larval presence by 8.15 (= 1/0.123). The treatment effect, therefore, is not simply due to annual mosquito fluctuations confounded with treatment. This detrimental impact of larvicide treatment on larval presence is reflected in the trend of reduced weekly mean number of larvae per post-treatment catch basin compared with pre-treatment ones over the mosquito breeding season (Fig. 4). Its effect is also illustrated in the trend of lower proportion of samples with larvae present among post-treatment samples compared with pre-treatment ones, regardless of WT (Fig. 5) or vegetation type of immediately surrounding area (Fig. 6). It may also, in part, account for the wide annual variation in mean number of larvae per catch basin sample; the years when substantial larvicide treatment of catch basins occurred (more than 10% of sampled catch basins) had lower mean values than those years when it did not (Fig. 3).

Fig. 4 The mean number of mosquito larvae per sample from pre-treatment (blue) and post-treatment (orange) catch basins each week (2004–2010).

Fig. 5 The proportion of pre-treatment (blue) and post-treatment (orange) catch basins containing mosquito larvae according to mean catch basin WT (2004–2010). Error bars show score confidence intervals.

Fig. 6 The proportion of pre-treatment (blue) and post-treatment (orange) catch basins containing mosquito larvae according to vegetation type of immediately surrounding area to catch basins (2004–2010). Error bars show score confidence intervals.

There was marked seasonal variation in the weekly mean number of larvae per catch basin: few larvae were observed until mid-June each year, with numbers usually peaking in early August before declining rapidly in September (Fig. 4). Linked to this seasonal variation in larval abundance is the positive correlation between WT and larval presence, which was supported by model parameter estimates (β7 [WT] = 0.14, P-value ≤0.001, Fig. 5, Appendix 1).

The model also indicates that larvae were less likely to occur in catch basins associated with paved areas (63% of observations) than in those beside grassy areas (13%) (β5 [paved] = −0.40, P-value ≤0.001), with catch basins in tree-lined areas (24%) providing the highest chance of larval presence (β6 [treeline] = 0.24, P-value = 0.019, Fig. 6, Appendix 1).

The cumulative effect of larvicide treatment, WT, and vegetation type of immediately surrounding area on larval presence is illustrated in Fig. 7: larvicide treatment consistently reduced the presence of larvae regardless of temperature (which tended to peak with larval presence in mid-summer) and adjacent land use (where larval presence tended to be lowest in paved areas and highest in tree-lined ones).

Fig. 7 Plots illustrating generalised linear mixed model fit. For every calendar day during the observation period (late-June to late-September 2004–2010), the vertical position of circles and diamonds indicate the proportion of larvae positive samples collected for untreated and treated sites, respectively. Further, the circle (diamond) size corresponds to the total number of untreated (treated) samples collected for the given day. The solid lines show probability of occurrence given vegetation type of immediately surrounding area (A: pavement, B: grass, C: trees), mean monthly WT (as observed and listed below months), and given clear weather.

In order to evaluate the goodness of fit of the logistic model to the data, we can visually evaluate how the data compared with the model fit (Fig. 7). We can also consider statistical measures, such as the Cox–Snell and Nagelkerke Pseudo R2s, which provide intuitive measures of the quality with which the model approximates the data (see Cox and Snell Reference Cox and Snell1989; Nagelkerke Reference Nagelkerke1991 for details). When comparing our model to that of the ordinary “intercept only” model, the Cox–Snell Pseudo R2 is 0.22 and the Nagelkerke Pseudo R2 is 0.31. The pseudo R2 values can be interpreted in the same manner as the R2 coefficient of determination used in ordinary linear regression, indicating that the model provides a good approximation to the data and a basis for prediction.

Discussion

With over 90 000 mosquito larvae collected from nearly 10 500 catch basin sampling events, this study demonstrates the capacity of artificial waterbodies to harbour mosquito larvae. Indeed, artificial waterbodies tend to harbour markedly more mosquito larvae than natural wetlands (Su et al. Reference Su, Webb, Meyer and Mulla2003, Jackson et al. Reference Jackson, Gow, Evelyn, Meikleham, McMahon and Koga2009). This may be because temporary, artificial habitats, such as catch basins, often lack predators found in more natural water bodies that help to regulate mosquito populations.

While a range of natural control measures, such as vegetation enhancement, may be effective in many above ground habitats, such as semi-natural wetlands (Jackson et al. Reference Jackson, Gow, Evelyn, Meikleham, McMahon and Koga2009), such methods are impractical for managing mosquito populations in artificial water bodies such as extensive networks of catch basins, leaving targeted larvicide treatment as the only viable course of control available (Metzger Reference Metzger2004). The efficacy that such treatment can have in reducing larval presence in catch basins has been previously documented (Siegel and Novak Reference Siegel and Novak1997, Reference Siegel and Novak1999), and is corroborated by our study; parameter estimates show that larvicide treatment reduced the odds of larvae presence by more than a factor of seven.

However, many jurisdictions have limited resources and may not be able to incorporate a wide-scale programme that encompasses the entire drainage system of a given region, which may include many thousands of catch basins. Our data suggests that in these cases an effective targeted approach may be achievable. For example, by predicting which catch basins are most likely to contain larvae at any given time, the relationships we observed between larval presence and variation in vegetation type, WT, and season could help to increase the efficiency and efficacy of targeted mosquito control operations.

Treatment could be prioritised based on the vegetation type of the immediately surrounding area; according to our findings, catch basins located in paved areas were significantly less likely to harbour larvae compared with those in grassed areas, with those found among trees having the highest chance of larval presence. In warm water with high levels of organic material, algae may prosper, providing ample food for filter feeding mosquito larvae. Lacey et al. (Reference Lacey, Lacey, Peacock and Thiery1988) also observed a decline in larvicidal activity of B. sphaericus in habitats that were highly organically enriched.

While it would be useful to investigate the exact role of different forms of organic input on mosquito production in the future, our results suggest it could be advantageous to focus catch basin cleaning resources in areas dominated by C. pipiens to tree-lined and grass-lined streets prior to the onset of the mosquito breeding season and larvicide treatment; indeed, removing organic matter through catch basin cleaning can decrease the number of mosquito larvae in catch basins by up to 40% (Stockwell et al. Reference Stockwell, Wessell, Reed, Kronenwetter-Koepel, Reed and Turchi2006). These decisions may be area specific, however, as habitat preferences vary between mosquito species. For example, Rey et al. (Reference Rey, O'Meara, O'Connell and Cutwa-Francis2006) found that C. quinquefasciatus and C. nigripalpus Theobald (Diptera: Culicidae) in Florida, United States of America were most common in paved areas.

In addition to vegetation adjacent to catch basins, mosquito larval abundance in catch basins also varies markedly with WT throughout the mosquito breeding season. The presence of few larvae before mid-June each year is consistent with Spielman (Reference Spielman2001), who found that C. pipiens first occurred when WTs reached ∼15 °C. Given this consistent finding, conducting catch basin sampling before week 25 or 26 seems to be an unnecessary use of resources, at least within this study area in southwestern British Columbia.

Although numbers usually peaked in early August before declining rapidly in September, the peaks in larval numbers shown in Fig. 4 may be revealing successive cohorts of overlapping generations, as has been observed in other studies of C. pipiens (Spielman Reference Spielman2001; Edillo et al. Reference Edillo, Kiszewski, Manjourides, Pagano, Hutchinson and Kyle2009). While more research into these patterns is needed, these peaks may be used to judge the optimal time of larvicide treatments. For example, a prolonged residual effect has been observed when intact mosquito larval cadavers are present for larvicides such as B. sphaericus used in this study, which uses their nutrients to grow and sporulate (Becker et al. Reference Becker, Zgomba, Petric, Beck and Ludwig1995; Uspensky et al. Reference Uspensky, Klein and Braun1998; Lacey Reference Lacey2007).

Long-term management of mosquitoes should involve new structural designs that consider reducing the numbers of catch basins with sumps, and incorporating more filtration measures. Research into different catch basin design and modifications, and their effect on mosquitoes is in progress (Harbison et al. Reference Harbison, Metzger, Allen II and Hu2009a, Reference Harbison, Metzger and Hu2010; Hamer et al. Reference Hamer, Kelly, Focks, Goldberg and Walker2011).

The rapid advance of WNv in recent years across North America warns of potential future health concerns arising from vector association with stormwater management. There are many other potential mosquito vectors of disease around the world that could invade new geographical regions. For example, Aedes japonicus Theobald (Diptera: Culicidae), a bridge vector of WNv and a vector of both dengue and chikungunya viruses, has been introduced separately to Ontario, Canada, and Washington state, United States of America (Thielman and Hunter Reference Thielman and Hunter2007; Irish and Pierce Reference Irish and Pierce2008). Larvae of this species are known to inhabit catch basins (Anderson et al. Reference Anderson, Ferrandino, Dingman, Main, Andreadis and Becnel2011). Aedes albopictus Skuse (Diptera: Culicidae), the Asian Tiger mosquito and a known vector of dengue, has spread across 36 states in the United States of America and much of Central America since its discovery in Tennessee, United States of America in 1983 (Enserink Reference Enserink2008). Outbreaks of other zoonotic mosquito-borne diseases are also on the rise (Jones et al. Reference Jones, Patel, Levy, Storeygard, Balk and Gittleman2008). For example, there have been recent outbreaks of the arbovirus chikungunya in Italy and Australia, also carried by A. albopictus (Reiter et al. Reference Reiter, Fontenille and Paupy2006). Indeed, the incidence of mosquito-borne diseases may increase due to anthropomorphic changes in the distribution and abundance of mosquito vectors brought about, for example, through global trade and travel, the decline of predatory species, degradation of natural habitats and increased urbanisation, and climate change.

By prioritising and targeting control efforts to most efficiently reduce WNv mosquito vector populations, the knowledge gleaned from this study can help to minimise the risk of WNv transmission to humans in Metro Vancouver, and may be extrapolated to other metropolitan areas in North America. A similar research strategy could be applied to emerging threats from other potential mosquito vectors of disease around the world to inform best practices to manage mosquitoes. With more information regarding the characteristics of the artificial habitat preferred by such vector mosquitoes, their populations can be reduced more effectively, and the incidence of mosquito-borne disease minimised.

Acknowledgements

The authors thank the past and present staff of Culex Environmental Ltd., Metro Vancouver, the Sunshine Coast communities, Vancouver Island communities, the Fraser Health Authority, the University of British Columbia, Simon Fraser University, and the Greater Vancouver municipalities who kindly contributed to the combined data presented in this paper. This study was conducted as part of Culex Environmental's Research and Development Program, which is supported by the Canada Revenue Agency Scientific Research and Experimental Development (SR&ED) Program.

Appendix 1 Summary of results for all tested variables in the mixed effects logistic regression model.

Footnotes

Subject editor: Kevin Floate

References

Anderson, J.F., Ferrandino, F.J., Dingman, D.W., Main, A.J., Andreadis, T.G., Becnel, J.J. 2011. Control of mosquitoes in catch basins in Connecticut with Bacillus thuringiensis israelensis, Bacillus sphaericus, [corrected] and spinosad. Journal of the American Mosquito Control Association, 27: 4555.CrossRefGoogle ScholarPubMed
Becker, N., Zgomba, M., Petric, D., Beck, M., Ludwig, M. 1995. Role of larval cadavers in recycling processes of Bacillus sphaericus. Journal of the American Mosquito Control Association, 11: 329334.Google ScholarPubMed
Belton, P. 1983. The mosquitoes of British Columbia. British Columbia Provincial Museum Handbook, 41: 1189.Google Scholar
Belton, P. 2006. British Columbia mosquitoes as vectors of West Nile virus [online]. Available from http://www.sfu.ca/~belton/summary.pdf [accessed 12 July 2012].Google Scholar
Bolker, B.M., Brooks, M.E., Clark, C.J., Geange, S.W., Poulsen, J.R., Stevens, M.H., et al. 2009. Generalized linear mixed models: a practical guide for ecology and evolution. Trends in Ecology and Evolution, 24: 127135.CrossRefGoogle ScholarPubMed
Cox, D.R.Snell, E.J. 1989. The analysis of binary data, 2nd edition. Chapman and Hall, London, United Kingdom.Google Scholar
Edillo, F., Kiszewski, A., Manjourides, J., Pagano, M., Hutchinson, M., Kyle, A., et al. 2009. Effects of latitude and longitude on the population structure of Culex pipiens s.l., vectors of West Nile virus in North America. The American Journal of Tropical Medicine and Hygiene, 81: 842848.CrossRefGoogle ScholarPubMed
Enserink, M. 2008. A mosquito goes global. Science, 320: 864866.CrossRefGoogle ScholarPubMed
Hamer, G.L., Kelly, P.H., Focks, D.A., Goldberg, T.L., Walker, E.D. 2011. Evaluation of a novel emergence trap to study Culex mosquitoes in urban catch basins. Journal of the American Mosquito Control Association, 27: 142147.CrossRefGoogle ScholarPubMed
Harbison, J.E., Metzger, M.E., Allen II, V., Hu, R. 2009a. Evaluation of manhole inserts as structural barriers mosquito entry into below ground stormwater systems using a simulated treatment device. Journal of the American Mosquito Control Association, 25: 356360.CrossRefGoogle Scholar
Harbison, J.E., Metzger, M.E., Hu, R. 2010. Association between Culex quinquefasciatus (Diptera: Culicidae) oviposition and structural features of below ground stormwater treatment devices. Journal of Medical Entomology, 47: 6773.CrossRefGoogle Scholar
Harbison, J.E., Metzger, M.E., Walton, W.E., Hu, R. 2009b. Evaluation of factors for rapid development of Culex quinquefasciatus in below ground stormwater treatment devices. Journal of Vector Ecology, 34: 182190.CrossRefGoogle Scholar
Irish, S.R.Pierce, C.S. 2008. Update on the distribution of Ochlerotatus japonicus in Oregon and Washington. Journal of the American Mosquito Control Association, 24: 110111.CrossRefGoogle ScholarPubMed
Jackson, M.J., Gow, J.G., Evelyn, M.E., Meikleham, N.E., McMahon, T.J.S., Koga, E., et al. 2009. Culex mosquitoes, West Nile virus and the application of innovative management in the design and management of stormwater retention ponds in Canada. Water Quality Research Journal of Canada, 44: 103110.CrossRefGoogle Scholar
Jones, K.E., Patel, N.G., Levy, M.A., Storeygard, A., Balk, D., Gittleman, J.L., et al. 2008. Global trends in emerging infectious diseases. Nature, 451: 990994.CrossRefGoogle ScholarPubMed
Kilpatrick, A.M. 2011. Globalization, land use and the invasion of West Nile virus. Science, 334: 323327.CrossRefGoogle ScholarPubMed
Komar, N.Clark, G.G. 2006. West Nile virus activity in Latin America and the Caribbean. Revista Panamericana de Salud Pública, 19: 112117.CrossRefGoogle ScholarPubMed
Lacey, L.A. 2007. Bacillus thuringiensis serovariety israelensis and Bacillus sphaericus for mosquito control. Journal of the American Mosquito Control Association, 23: 133163.CrossRefGoogle ScholarPubMed
Lacey, L.A., Lacey, C.M., Peacock, B., Thiery, I. 1988. Mosquito host range and field activity of Bacillus sphaericus isolate 2297 (serotype 25). Journal of the American Mosquito Control Association, 4: 5156.Google ScholarPubMed
Metzger, M.E. 2004. Managing mosquitoes in stormwater treatment devices. University of California, Davis, Agriculture and Natural Resources Publication, Davis, California, United States of America, 8125: 111.CrossRefGoogle Scholar
Morris, J.A., Lampman, R.L., Ballmes, G., Funes, J., Halvorsen, J., Novak, R.J. 2007. First record of Aedes japonicus japonicus in Illinois: defining its spatial distribution and associated mosquito species. Journal of the American Mosquito Control Association, 23: 243251.CrossRefGoogle ScholarPubMed
Nagelkerke, N.J.D. 1991. A note on a general definition of the coefficient of determination. Biometrika, 78: 691692.CrossRefGoogle Scholar
Nielsen, C.F., Reisen, W.K., Armijos, V., Wheeler, S., Kelley, K., Brown, D. 2007. Impacts of adult mosquito control and climate variation on the West Nile Virus epidemic in Davis, during 2006. Proceedings of the Mosquito and Vector Control Association of California, 75: 125130.Google Scholar
Norris, D.E. 2004. Mosquito-borne diseases as a consequence of land use change. EcoHealth, 1: 1924.CrossRefGoogle Scholar
Reisen, W.K.Brault, R. 2007. West Nile virus in North America: perspectives on epidemiology and intervention. Pest Management Science, 63: 641646.CrossRefGoogle ScholarPubMed
Reiter, P., Fontenille, D., Paupy, C. 2006. Aedes albopictus as an epidemic vector of chikungunya virus: another emerging problem? The Lancet Infectious Diseases, 6: 463464.CrossRefGoogle ScholarPubMed
Rey, J.R., O'Meara, G.F., O'Connell, S.M., Cutwa-Francis, M.M. 2006. Factors affecting mosquito production from stormwater drains and catch basins in two Florida cities. Journal of Vector Ecology, 31: 334343.CrossRefGoogle ScholarPubMed
Siegel, J.P.Novak, R.J. 1997. Field trials of Vectolex CG® a Bacillus sphaericus larvicide in Illinois waste tires and catch basins. Journal of the American Mosquito Control Association, 13: 305310.Google Scholar
Siegel, J.P.Novak, R.J. 1999. Duration of activity of the microbial larvicide Vectolex CG® (Bacillus sphaericus) in Illinois catch basins and waste tires. Journal of the American Mosquito Control Association, 15: 366370.Google ScholarPubMed
Spielman, A. 2001. Structure and seasonality of Nearctic Culex pipiens populations. Annals of the New York Academy of Sciences, 957: 220234.CrossRefGoogle Scholar
Stockwell, P.J., Wessell, N., Reed, D.R., Kronenwetter-Koepel, T.A., Reed, K.D., Turchi, T.R., et al. 2006. A field evaluation of four larval mosquito control methods in urban catch basins. Journal of the American Mosquito Control Association, 22: 666671.CrossRefGoogle ScholarPubMed
Su, T. 2008. Evaluation of water-soluble pouches of Bacillus sphaericus applied as pre-hatch treatment against Culex mosquitoes in simulated catch basins. Journal of the American Mosquito Control Association, 24: 5460.CrossRefGoogle Scholar
Su, T., Webb, J.P., Meyer, R.P., Mulla, M.S. 2003. Spatial and temporal distribution of mosquitoes in underground storm drain systems in Orange County, California. Journal of Vector Ecology, 28: 7989.Google ScholarPubMed
Thielman, A.C.Hunter, F.F. 2007. Photographic key to the adult female mosquitoes (Diptera: Culicidae) of Canada [online]. Canadian Journal of Arthropod Identification, 4. Available from www.biology.ualberta.ca/bsc/ejournal/th_04/th_04.html [accessed 31 March 2013].Google Scholar
Thomson, A. 2008. Quantitative sampling and determination of prevalence of Culex mosquitoes in catch basins and their association with canopy cover in Manitoba. M.Sc. thesis. University of Manitoba, Manitoba, Canada.Google Scholar
Turell, M.J., Dohm, D.J., Sardelis, M.R., O'Guinn, M.L., Andreadis, T.G., Blow, J.A. 2005. An update on the potential of North American mosquitoes (Diptera: Culicidae) to transmit West Nile virus. Journal of Medical Entomology, 42: 5762.CrossRefGoogle ScholarPubMed
Uspensky, I., Klein, D., Braun, S. 1998. Persistence of Bacillus sphaericus in cadavers of mosquito larvae. Israel Journal of Entomology, 32: 4956.Google Scholar
Vazquez-Prokopec, G.M., Vanden Eng, J.L., Kelly, R., Mead, D.G., Kolhe, P., Howgate, J., et al. 2010. The risk of West Nile virus infection is associated with combined sewer overflow streams in urban Atlanta, Georgia, USA. Environmental Health Perspectives, 118: 13821388.CrossRefGoogle ScholarPubMed
Vinogradova, E.B. 2000. Culex pipiens pipiens mosquitoes: taxonomy, distribution, ecology, physiology, genetics, applied importance and control. Pensoft Publishers, Sofia, Bulgaria.Google Scholar
Wood, D.M., Dang, P.T., Ellis, R. 1979. The mosquitoes of Canada – Diptera: Culicidae. Publication 1686. Agriculture Canada, Hull, Quebec, Canada.Google Scholar
Figure 0

Fig. 1 A schematic cross-section through a typical sumped catch basin. Such basins commonly provide habitat for mosquito larvae.

Figure 1

Fig. 2 Study area showing locations of catch basins monitored in the Lower Mainland and Vancouver Island, British Columbia, Canada from 2004 to 2010.

Figure 2

Fig. 3 The mean number of mosquito larvae per catch basin sample according to year (2004–2010). Error bars show standard errors. The percentage of samples occurring post-treatment is in brackets above bars. The total number of sampling events is below each year.

Figure 3

Fig. 4 The mean number of mosquito larvae per sample from pre-treatment (blue) and post-treatment (orange) catch basins each week (2004–2010).

Figure 4

Fig. 5 The proportion of pre-treatment (blue) and post-treatment (orange) catch basins containing mosquito larvae according to mean catch basin WT (2004–2010). Error bars show score confidence intervals.

Figure 5

Fig. 6 The proportion of pre-treatment (blue) and post-treatment (orange) catch basins containing mosquito larvae according to vegetation type of immediately surrounding area to catch basins (2004–2010). Error bars show score confidence intervals.

Figure 6

Fig. 7 Plots illustrating generalised linear mixed model fit. For every calendar day during the observation period (late-June to late-September 2004–2010), the vertical position of circles and diamonds indicate the proportion of larvae positive samples collected for untreated and treated sites, respectively. Further, the circle (diamond) size corresponds to the total number of untreated (treated) samples collected for the given day. The solid lines show probability of occurrence given vegetation type of immediately surrounding area (A: pavement, B: grass, C: trees), mean monthly WT (as observed and listed below months), and given clear weather.

Figure 7

Appendix 1 Summary of results for all tested variables in the mixed effects logistic regression model.