Geostatistical evaluation of integrated marsh management impact on mosquito vectors using before-after-control-impact (BACI) design
© Rochlin et al; licensee BioMed Central Ltd. 2009
Received: 28 April 2009
Accepted: 23 June 2009
Published: 23 June 2009
In many parts of the world, salt marshes play a key ecological role as the interface between the marine and the terrestrial environments. Salt marshes are also exceedingly important for public health as larval habitat for mosquitoes that are vectors of disease and significant biting pests. Although grid ditching and pesticides have been effective in salt marsh mosquito control, marsh degradation and other environmental considerations compel a different approach. Targeted habitat modification and biological control methods known as Open Marsh Water Management (OMWM) had been proposed as a viable alternative to marsh-wide physical alterations and chemical control. However, traditional larval sampling techniques may not adequately assess the impacts of marsh management on mosquito larvae. To assess the effectiveness of integrated OMWM and marsh restoration techniques for mosquito control, we analyzed the results of a 5-year OMWM/marsh restoration project to determine changes in mosquito larval production using GIS and geostatistical methods.
The following parameters were evaluated using "Before-After-Control-Impact" (BACI) design: frequency and geographic extent of larval production, intensity of larval production, changes in larval habitat, and number of larvicide applications. The analyses were performed using Moran's I, Getis-Ord, and Spatial Scan statistics on aggregated before and after data as well as data collected over time. This allowed comparison of control and treatment areas to identify changes attributable to the OMWM/marsh restoration modifications.
The frequency of finding mosquito larvae in the treatment areas was reduced by 70% resulting in a loss of spatial larval clusters compared to those found in the control areas. This effect was observed directly following OMWM treatment and remained significant throughout the study period. The greatly reduced frequency of finding larvae in the treatment areas led to a significant decrease (~44%) in the number of times when the larviciding threshold was reached. This reduction, in turn, resulted in a significant decrease (~74%) in the number of larvicide applications in the treatment areas post-project. The remaining larval habitat in the treatment areas had a different geographic distribution and was largely confined to the restored marsh surface (i.e. filled-in mosquito ditches); however only ~21% of the restored marsh surface supported mosquito production.
The geostatistical analysis showed that OMWM demonstrated considerable potential for effective mosquito control and compatibility with other natural resource management goals such as restoration, wildlife habitat enhancement, and invasive species abatement. GPS and GIS tools are invaluable for large scale project design, data collection, and data analysis, with geostatistical methods serving as an alternative or a supplement to the conventional inference statistics in evaluating the project outcome.
The salt marsh is a globally important ecosystem in high to middle latitudes along the coastline . Ecologically, salt marshes provide a nutrient rich interface between terrestrial and marine environments utilized by a great variety of animal species. The salt marsh habitat is also of a significant public health importance due to mosquito vector species that have adapted to this harsh environment. Pathogens transmitted by salt marsh mosquitoes include the malaria parasite vectored by Anopheles atroparvus, An. sacharovi, and An. labranchiae in Europe and the Middle East , Venezuelan Equine Encephalitis virus vectored by Aedes sollicitans and Ae. taeniorhynchus in the Americas , California group encephalitis viruses vectored by Ae. dorsalis, Ae. caspius, and Ae. melanimon in Europe and western North America , and Ross River virus vectored by Ae. camptorhynchus and Ae. vigilax in Australia . Many of the salt marsh Aedes mosquitoes are also important biting pests species in coastal population centers and tourist areas [2, 5].
The salt marsh mosquito fauna of Long Island, New York is representative of the Atlantic seaboard of the continental US consisting of 4 species, Ae. sollicitans, Ae. cantator, Ae. taeniorhynchus, and Cx. salinarius [6, 7]. Aedes sollicitans is considered the main epidemic vector of Eastern Equine Encephalitis virus (EEEv) in coastal areas of eastern US [8–10]. The virus has been repeatedly isolated from the wild populations of Ae. sollicitans in this region [11, 12] and occasionally from Ae. cantator and Ae. taeniorhynchus [11, 13]. Culex salinarius is a potentially important epidemic vector of EEEv  with multiple field isolations of the virus [15–17]. This mosquito species is also one of the main vectors involved in West Nile virus (WNV) human transmission [18–20]. In addition to numerous WNV isolations from field collected Cx. salinarius, the virus was detected in Ae. sollicitans (including specimens collected on Long Island) and occasionally in Ae. cantator and Ae. taeniorhynchus [21, 22]. All three salt marsh Aedes species can potentially transmit WNV to humans [23, 24], although only Ae. sollicitans has been associated with the risk of WNV transmission to humans .
The public health importance of the salt marsh mosquitoes in the coastal areas of eastern US was apparent long before the discovery of the mosquito-borne viruses. Massive infestations by Ae. sollicitans and Ae. taeniorhynchus led to the establishment of most, if not all, coastal mosquito control districts [26, 27]. Abatement of larval mosquitoes on the salt marsh employed both chemical and habitat modification techniques including ditching of the marsh surface to allow rapid draining of small pools harboring mosquito larvae. Similarly to many coastal areas, Long Island salt marshes were grid ditched in the late 1930s . Although reasonably effective as a mosquito control tool, universal grid ditching was perceived as unnecessary in those parts of salt marshes not producing mosquitoes, and unsatisfactory from ecological and resource conservation perspectives . These concerns resulted in the development of an Open Marsh Water Management (OMWM) technique with the dual goal of non-chemical mosquito control and marsh conservation [29–31]. OMWM targets specific areas of known mosquito larval habitat by employing tidal channels, ponds, and shallow radial ditches to remove these habitats . Tidal channels are designed to improve water circulation and to restore natural tidal regime, while small ponds are created to replace clusters of depressions where mosquito larvae proliferate. Shallow radials connect ponds with scattered depressions to allow access by larvivorous fish. With some modifications, OMWM has been put into practice in many mosquito producing salt marshes along the US Atlantic coast from Massachusetts to Florida, California, and Australia (reviewed in ). Invariably, satisfactory results were reported, with the initial estimates suggesting elimination of 40–60 billion mosquitoes annually for every 1,000 OMWM acres . Direct measurements of mosquito production, although infrequent, documented significant reductions (> 85%) in larval or adult mosquito levels (reviewed in [33, 29]); however, many studies relied on rather qualitative observations such as number of complaints to assess the outcome. A recent OMWM survey for US Fish and Wildlife Service (USFWS)  utilized before-after-control-impact (BACI) design  allowing statistically rigorous analysis of the technique's effects . Although we adopted a similar BACI approach for our project, the two studies differed significantly in scope, goals, and methodology.
Suffolk County occupies most of Long Island, New York and employs both pesticides and water management to control salt marsh mosquitoes. To investigate alternative approaches in compliance with the pesticide usage reduction goal set by the County government, a partnership with USFWS was initiated for a pilot integrated marsh management project at Wertheim National Wildlife Refuge (Wertheim NWR). The following goals of the project included both marsh restoration and mosquito control components: 1) to decrease mosquito production and, by extension, pesticide usage by utilizing OMWM approach, 2) to reduce the vigor and extent of the invasive reed Phragmites australis , 3) to naturalize marsh surface by eliminating grid ditching, and 4) to maintain or enhance fish and wildlife habitat. This article describes the planning and implementation of the project and focuses on quantitative outcome evaluation for mosquito control. Geographic information systems (GIS) and GPS technology were systematically employed at each stage of the project allowing efficient use of limited resources and enabling development of novel methodology for assessing the project impacts.
To assist in mapping out the proposed hydrologic features on the ground, a network of georeferenced points was established in the treatment areas using Trimble® GPS receivers (Trimble Navigation Ltd., Sunnyvale, CA). The outline of ponds, tidal channels, connectors, and larval habitat targeted for elimination by filling in dense concentrations of larvae producing potholes  were staked out using the hydrologic feature map and georeferenced points as a guide. The stake locations were recorded by hand-held GPS receivers and the data were visualized in ArcGIS overlaying the hydrologic feature map. The stake positions were then adjusted if necessary. The following alterations were made in Area 1 (March 2005) and in Area 2 (February-March, 2006). Ponds, tidal channels, and connectors were constructed using low-ground pressure (< 2 psi) machinery. The majority of pre-existing mosquito control ditches in the marsh were filled with the material extracted during the pond excavation process using rotary ditchers and excavators. The remaining ditches were naturalized by adding curves and other features commonly found in natural salt marsh channels and creeks. Additional tidal channels were created with the overarching goal of increased tidal circulation into the interior of the salt marshes. The excess pond/tidal creek excavated material was used to grade hummocky high marsh terrain where larval habitat had been found by the preliminary surveillance (i.e. "backblading").
Entomological data collection
Two sampling procedures were employed to collect larval data. A systematic sampling protocol  was carried out to establish random transects (n = 5 in Areas 1 and 2 each, n = 4 in Areas 3 and 4 each) with sampling stations distributed equally at every 40 meters with sampling taking place at and between the stations (i.e. every 20 meters) resulting in a total of 24 sampling points in each Areas 1 and 2, and a total of 20 sampling points in each Areas 3 and 4. These transects extended from the upland portion of the salt marsh seaward into the low marsh. A second procedure utilized targeted sampling, a more common method used by mosquito control professionals, whereby suitable larval habitats, i.e. pools of standing water, were searched for mosquito larvae across all 4 Areas. These surveys intended to be comprehensive, so the majority of potential larval habitats within an area were sampled each time. At least 25 samples were obtained per area per visit unless the marsh surface was either dry (no standing water) or completely flooded.
Each area was visited weekly for targeted sampling from early May to mid September, a period corresponding to the active mosquito season in this region. Additionally, transect sampling was carried out monthly 4 to 7 days following a high tide inundation to maximize the chances of finding mosquito larvae on the salt marsh . Mosquito larvae were collected using standard mosquito dippers (Bioquip Inc., Rancho Dominguez, CA). Custom database design software for Windows CE (Visual CE, Syware, Cambridge, MA) was developed to run on a hand-held personal digital assistant (PDA; Dell Axim™ X51, Dell Inc., Round Rock, TX) coupled with a GPS device (GPSlim 236, Holux Technology Inc., Hsinchu, Taiwan). The program recorded the following data: geographic coordinates, time, number of mosquito larvae and pupae, type of habitat, sampling method with transect/station information, and comments. The data were uploaded directly into a MS Access database using Microsoft ActiveSync. For data quality assurance, the sampling points were visualized in ArcMap and any errors were identified and corrected.
As a rule, 1 to 5 typical larval samples from each area were brought to laboratory for microscopic species identification using morphological mosquito keys . First instar larvae were allowed to progress to later instars, and larvae of Culex spp. were allowed to emerge as adults to confirm positive identification. Adult Cx. salinarius were separated from those of Cx. pipiens by using morphological characters  as well as molecular techniques .
Entomological data evaluation and statistical analysis
Three parameters were evaluated to assess the outcome of the project: frequency and geographic extent of larval production, intensity of larval production, and the overall impact of the project (both OMWM and restoration components) on mosquito larval habitat. In accordance with BACI design, the analyses were performed on aggregated "Before and After" data and also followed through time to compare "Control" and "Impact" (i.e. treatment) areas to identify changes attributable to the intervention.
Presence or absence of mosquito larvae, a dichotomous variable, was used to evaluate geographic patterns of larval production. Each sample was thus classified as either 1 (= positive dip) or 0 (= negative dip). Global and Anselin Local Moran's I were calculated to assess the overall spatial patterns and local geographic clusters of positive or negative dips using ArcMap 9.3 software (ESRI Inc, Redlands, CA). Normalized Z-scores, or the number of standard deviations, were interpreted to represent clusters (Z-score > 2.0), outliers (Z-score < -2.0), or random distribution (-2.0 < Z-score < 2.0) at statistical significance P < 0.05. To visualize a statistical surface of mosquito larval production over the entire study area, Voroni tessellation was performed to create Thiessen polygons around each data point classified by presence or absence of mosquito larvae and the Z-score. Adjacent polygons in the same category were dissolved and smoothed to produce the final map.
To define statistical significant clusters of positive and negative dips spatially and temporarily, a spatial scan cluster analysis was carried out using free SaTScan™ software . Positive dips (the case file) and negative dips (the control file) were analyzed by space-time statistic with Bernoulli probability model for dichotomous data employing the following settings: year as the time period, maximum non-overlapping spatial cluster size of 0.5 km radius roughly corresponding to the extent of each treatment and control area, scan for high and low values, and 4-year temporal window with pure spatial clusters (i.e. present each year). The statistical significance was calculated by Monte-Carlo simulation with 9999 replications.
To determine the intensity of mosquito larval production (termed "breeding intensity"), an ordinal scale was developed based on field observations. Dips with 1 to 2 larvae per dip were classified as "low" (rank = 1), those with 3–5 larvae per dip as "intermediate" (rank = 2), and those with > 5 larvae per dip as "high" (rank = 3). Only positive dips were included in the analysis using global (General G) and local (Gi*) Getis-Ord statistic in ArcMap 9.3 to assess the overall distribution of breeding intensity values and to identify geographic locations with elevated breeding intensity, i.e. non-random clusters. Similarly to Moran's I statistic, normalized Z-scores were interpreted to represent clusters (Z-score > 2.0), outliers (Z-score < -2.0), or random distribution (-2.0 < Z-score < 2.0) at statistical significance P < 0.05. To visualize distribution of positive dips on the marsh surface, kernel density was calculated using 5 × 5 meter grid and 50 meter search radius.
Post-project problem area characterization
ArcMap 9.3 buffering tool and sampling tools available with Hawth's ArcMap extension  were utilized to evaluate changes in mosquito larval habitat before and after the treatment. Specifically, the impact of filling in mosquito ditches (a marsh restoration technique) on larval habitat was assessed by determining proportion of positive dips at close (< 5 meters), medium (5–15 meters), or long (> 15 meters) distance range. All non-spatial statistical analyses were conducted in SPSS v. 15.0 (SPSS Inc, Chicago, IL) software and assumed statistical significance at P < 0.05.
Mosquito abatement measures consisted of aerial application of larvicides by helicopter. Routine mosquito control program continued throughout the study period using the standard criteria set by the state regulatory agencies and USFWS. At least 25 samples were required from each area to meet the minimum larviciding threshold of 0.2 larvae per dip. Other considerations included the extent of the infested area (i.e. total number of positive dips), weather, and environmental conditions directly affecting mosquito larval habitat such as marsh flooding. These criteria were uniformly applied without regard to the status (impact or control) of the particular location to determine whether a larvicide application was justified. Two types of larvicides were used in the control program. Vectobac 12AS (Bacillus thuringiensis var. israelensis; ValentBioscience Corp.) is a liquid bacterial product applied when early larval instars (stages 1–2) were detected. Altosid Liquid Larvicide Concentrate (methoprene; Central Life Science/Wellmark™) is an insect growth regulator applied against late instar larvae (stages 3–4). When both early and late larval instars were present simultaneously, a combination of Vectobac and Altosid was used.
Before and after treatment effects on number of larvicide applications and proportion of time the larviciding threshold was reached were analyzed according to the published guidelines for BACI designs . The differences between control and treatment sites were computed and compared before and after treatment by Mann-Whitney test. SPSS v. 15.0 was used for data processing and statistical analyses.
Number of samples (dips) taken per area per year.
Proportion of samples (dips) containing water in transect versus targeted sampling. N/R – not recorded.
Proportion of positive samples (dips containing mosquito larvae) per samples with water in transect versus targeted sampling.
Immature stages of 3 mosquito species, Ae. sollicitans, Ae. cantator, and Cx. salinarius were collected during the study period. The presence of Cx. salinarius in the upper salt marsh was unexpected; thus, the close association of this species with the salt marsh habitat was investigated and characterized . Although Ae. sollicitans was the most commonly found species in all 4 areas throughout the season, Cx. salinarius sometimes predominated in Area 3 when more permanent brackish water from rain events accumulated on the marsh during this species' peak season from late July through early September.
Entomological data evaluation and statistical analysis
Characterization of the post-project residual larval habitat in the treatment areas
Number of positive dips as a function of distance to filled in mosquito ditches before and after the intervention.
Distance to filled ditch, m
Inspection of the post-project maps showing the distribution of larvae in the treatment areas (2006 to 2008) indicated a trend of diminishing residual larval habitat. This was also supported by the decreased likelihood of clustering in Area 2, where most of the residual breeding occurred (Figure 9). On the ground, many sections of the filled-in ditches that had been bare and concave with ample larval habitat in 2006–07 were rapidly re-vegetating and becoming level with the marsh surface in 2008. Based on these observations, the process of reduction in the mosquito habitat on or near filled-in ditches is expected to continue in the future.
In addition to the filled-in ditches, two problems of a lesser magnitude contributing to residual mosquito breeding were identified: a) insufficient connectivity with some of the newly created ponds and tidal channels, potentially leading to limited access of the affected areas by larvivorous killifish, and b) increased water accumulation due to a plugged ditch. However, these problems were minor and could be easily corrected without significant impact on the salt marsh surface by either creating shallow radials or removing the plug.
Mean number of larvicide applications per month (May–September) in 2001–2008. Post-project values (treatment areas) are indicated in bold print.
To determine whether the observed reduction in the average number of monthly larvicide applications was due to reduced frequency of positive dips as suggested by the spatial data analysis, the average proportion of weeks per month when the larviciding threshold (> 0.2 larvae/dip) had been reached was computed for each area and compared before (2004–05) and after (2005–08) treatment (Figure 10B). To increase the sample size, Area 2 was compared to Areas 3 and 4 to provide pre-treatment data for 2005, while Area 1 was also compared to the two control areas to provide post-treatment data for the same year. The average proportion of weeks (mean ± SE) above the threshold level for Areas 1 and 2 decreased by approximately 44% from 0.57 ± 0.1 pre-treatment to 0.25 ± 0.05 post-treatment, while that for the control areas remained similar, 0.57 ± 0.1 pre-treatment and 0.60 ± 0.05 post-treatment. The difference in the proportion of weeks between the treatment and the control areas was significantly higher post-treatment (Mann-Whitney U = 21.5, P = 0.001) suggesting a statistically significant reduction in the number of weeks when larvicide applications were justified in Areas 1 and 2 following treatment.
This study's goals, design, and analysis differed significantly from the previous investigations on Open Marsh Water Management (OMWM). The goals of the Wertheim Integrated Marsh Management (IMM) project were not limited to achieving reduction in mosquito production (the OMWM component), but also concurrently included restoring the marsh surface by eliminating grid-ditching, and controlling the invasive species P. australis. To better assess the impact of these techniques on the marsh flora and fauna, a quasi-experimental before-after-control-impact (BACI) study design with 2 pairs of impact and control sites was utilized. A BACI approach was selected because it offered the closest approximation of a field study to a full experimental design to detect and evaluate the impact (i.e. treatment) effects [35, 42]. For mosquito production, geostatistical analysis of the spatial pattern of larval distribution on the marsh surface was used as an alternative to conventional statistical methods to improve the validity of the statistical analysis, to better assess the project effectiveness, and to fully understand the underlying causes of some of the challenges encountered during the study.
Initially, we planned to use random sampling at transect locations to evaluate changes in mosquito production by conventional statistical analysis  supplemented by geostatistical analysis of targeted sampling traditionally employed by mosquito control districts. However, about one-half of transect samples were dry (= "no data"), significantly higher than about 13% of those for targeted samples. In addition, only about 9 transect samples with larvae compared to about 112 targeted samples were collected on average per Area each year. Due to high degree of spatial dependency or autocorrelation, the amount of information in the few positive transect samples was further reduced leading to smaller effective sample size, underestimated variance, and increased type I error . Thus, a statistical analysis based exclusively on the transect samples would result in a very low statistical power from a purely technical perspective and questionable biological significance. Numerically superior targeted samples, on the other hand, violated two classical inference assumptions, namely independency between samples (similarly to transect samples) and random selection of the sampling locations. To circumvent these two issues commonly encountered in vector control practice, geostatistical approach was used instead. Presence of spatial dependency is one of the central assumptions in geostatistics, and its magnitude is an important parameter for assessing the spatial patterns such as larval clustering. Although probability sampling is required for an unbiased estimate of population parameters in conventional statistics, geostatistical model-based approach does not require random selection of sampling location . In addition, the representativeness of targeted sampling design used in this study was enhanced by the spatial scope seeking to encompass the entire population and by replication in time over a 5-year study period.
Commonly, quantitative assessment of larval populations relies on number of mosquito larvae per dip. However, this number is highly variable and dependent on many factors unrelated to true mosquito density. For example, mean number of larvae per dip varied significantly both among different operators and between repeated samples taken by the same operator from the same source . Dipper samples could not differentiate population densities below ~280/m2, and more than 6,000 samples were required to estimate the population parameters with α = 0.05 and β = 0.1 . Other factors more specific to the salt marsh mosquitoes may include the size of the pools, presence of larval aggregates, and time of the day, among other factors. Thus, Service's  extensive review of entomological literature concluded that larvae per dip could not serve as a true estimate of the larval population. Accordingly, we used presence/absence of mosquito larvae as the main mosquito population parameter in this study. From an operational perspective, the location and the geographic extent of larvae producing areas (i.e. "hotspots") are more important for implementing targeted mosquito abatement program. Large areas of salt marsh mosquito larval habitat can be rapidly characterized using presence/absence data entered into a handheld GPS unit while minimizing technical errors and increasing effective utilization of field personnel. Given similar breeding intensity (i.e. average number of larvae in positive dips) between the treatment and the control areas, frequency of positive samples was also directly proportional to the mean number of larvae per dip in this study.
The OMWM concept was originally developed to provide effective long term control of mosquito larvae by source reduction and biological control. Field data on OMWM projects collected over a 40-year period have been largely supportive of this statement. The magnitude of the reduction in mosquito production generally ranged from 85% to complete elimination [29, 33]. For example, Ferrigno  reported a reduction from 3.7 × 106 larvae/acre pre-OMWM to almost zero post-OMWM in the upper marsh S. patens treatment areas, while 1.5 × 106 larvae/acre were detected on average in the control areas. Similarly, a 99% difference of 3.3 versus 0.02 larvae/dip was found between ditched marsh control and OMWM sites in Massachusetts . Meredith and Lesser  found 92% reduction in larval densities and 78% reduction of finding mosquito larvae (i.e. frequency) on average summarizing the results of a 28-year OMWM implementation in Delaware. James-Pirri et al.,  also observed reduction in both the proportion of time the mosquitoes were present and the larval density, although these trends were somewhat obscured by the parameter variability at the control sites. In our study, the frequency of finding larvae in the treatment areas post-project was reduced by ~70% on average, while remaining essentially identical in the control areas pre and post-project. Although the magnitude of this change was somewhat lower than that typically reported in the literature for an OMWM project, this reduction led to marked differences in spatial patterns of larval distribution on the ground. Statistically significant clusters of larvae were no longer present in the treatment areas, but consistently remained in place in the control areas post-project. Moreover, most residual breeding in the treatment areas post-project did not overlap with the pre-project larval habitat, but occupied a new niche atop or near the filled-in ditches created to restore the marsh surface. This finding highlights the potential negative consequences of marsh restoration for mosquito production if new larval habitat is generated during the process. Some of the important mosquito vectors in our area such as Cx. salinarius, can successfully utilize both heavily disturbed and relatively pristine marshes , which may result in an increased mosquito production from restored areas under favorable environmental conditions. However, even with almost complete elimination of the grid ditching system, only about 21% of the filled-in locations supported new larval habitat suggesting that the majority of the grid ditches contributed very little or none to mosquito larval habitat. For that reason, restoring the marsh surface by removing grid ditching is conceivable if proper surveillance measures to identify problem areas are implemented. GPS based monitoring and GIS/geostatistical techniques such as those described in this article represent crucial components of any surveillance program if a comprehensive project evaluation is required.
Despite the residual mosquito larval habitat in the treatment areas, the OMWM ultimate goal of significantly reducing larviciding while providing sufficient mosquito control was accomplished. The number of larvicide applications was lowered by about 74% in the treatment areas following treatment. This was due to two factors that directly affected the pesticide application criteria: the larviciding threshold and the spatial extent of the mosquito breeding areas on the marsh surface. The larviciding threshold of 0.2 larvae per dip was reached less frequently in the treatment areas by approximately a factor of 2 following treatment due to fewer positive samples containing larvae. The number of larvae per positive dip (i.e. breeding intensity), however, remained similar between treatment and control areas. This observation may be attributed to highly efficient predation of mosquito larvae by killifish in the accessible areas within Areas 1 and 2, whereas locations in the same treatment areas not easily accessible by larvivorous killifish and containing larval habitat (such as the surface of the newly filled-in ditches) continued to support mosquito breeding at the similar intensity to that of pre-project. Changes in the spatial distribution of the mosquito larvae, with reduced extent and loss of clustering in treatment areas also contributed to fewer larvicide applications compared to those in the control areas (Figures 5 and 7). Although 74% reduction in number of larvicide applications is slightly lower than 90–100% reported by other investigators [29, 32, 49], this difference may be attributed to the expanded scope of this project (i.e. marsh restoration discussed above), lower larviciding thresholds, and more rigorous monitoring procedures.
Continuation of larviciding activities throughout the study period illustrates the difficulties in conducting large scale experiments in natural settings. As was noted previously, pesticide application may confound the results on mosquito production making their interpretation more difficult . To avoid potential bias, Wertheim IMM adopted a set of criteria for larviciding triggers, which were uniformly applied to both treatment and control areas. Using these criteria, the treatment areas consistently received significantly fewer larvicide applications during the post-treatment period (Table 5). In this case, the confounding effect of larviciding would be expected to lessen the differences in mosquito production between treatment and control areas thus leading to a decrease in the before-after effect. However, the differences attributable to OMWM were not only detectable, but statistically significant. Control areas supported higher mosquito production despite retaining intact grid ditching and being subjected to 3–4 times more larvicide applications than did the treatment areas post-project. Thus, the OMWM component in this IMM project demonstrated its potential to largely replace chemical control and marsh-wide parallel grid ditching for effective larval mosquito control.
This study investigated the effectiveness of Open Marsh Water Management (OMWM) for mosquito vector control when combined with salt marsh restoration and invasive plant species control. Significant reduction was achieved in the frequency of finding larvae on the marsh surface leading to loss of spatial larval clusters or "hotspots" in the areas under OMWM. In turn, these changes resulted in a significant decline of the number of larvicide applications in those areas. Random transect sampling was inadequate to assess the mosquito larval population due to a large proportion of dry "no data" sampling points. More informative targeted sampling necessitated extensive application of GPS and GIS tools to collect and analyze the data using geostatistical methods as an alternative or a supplement to the conventional inference statistics. Geographic analysis was also instrumental in identifying the residual post-OMWM larval habitat, which was largely confined to some of the restored marsh surface, albeit only in a relatively small fraction of the total. Overall, geospatial analysis proved to be a highly useful tool in evaluating the project and more completely understanding how the marsh alterations impacted mosquito larval habitats.
Although mosquito breeding was greatly reduced but not eliminated in OMWM marshes, this technique demonstrated considerable potential for effective mosquito control, which is also compatible with other natural resource management goals such as restoration, wildlife habitat enhancement, and invasive species abatement. No further interventions have been carried out in the treatment areas up to date. In the future, limited scope refinements to reduce or eliminated the remaining larval mosquito habitat may be considered based on geospatial surveillance results.
The full report on Wertheim National Wildlife Refuge Integrated Marsh Management project is available at:
This website contains links to the comprehensive 2003–2007 report and the data sets of various parameters monitored during the project including mosquito larval surveillance.
Before After Control Impact
Geographic Information Systems
Geographic Positioning System
Integrated Marsh Management
Open Marsh Water Management
United States Fish and Wildlife Service
We wish to express our appreciation and gratitude to many people who made this project possible. The project was carried out in a partnership with U.S. Fish and Wildlife Service: significant contributions to the project were made by Patricia Martinkovic, Deborah Long, and Michelle Williams (project design and permit applications), Susan Adamowicz and Mark Maghini (project design), Alex Chmielweski and Florence James (mosquito sampling). We thank all Suffolk County Vector Control staff who participated in this project, specifically Valentin Bulgak (GIS and GPS support), Peter Dykeman (mosquito sampling) and Rich Kowalski (heavy equipment supervisor). Many thanks to Suffolk County Dept. of Health Services for their assistance and advice: Walter Dawydiak and Kim Shaw (project implementation), Phil DeBlasi (study site delineation). Ducks Unlimited: Craig Kessler and Josh Ver Hague assisted with study site delineation and mosquito sampling, respectively. The following Cashin Associates (project contractor) personnel contributed to the project: David Tonjes (permit application), Keith Brewer and Kim Somers (study site delineation).
- Mitsch WJ, Gosselink JG: Tidal salt marshes. Wetlands. 2000, New York: John Wiley & Sons, Inc., 261-305. 3Google Scholar
- Becker N, Zgomba M, Petric D, Dahl C, Boase C, Lane J, Kaiser A: Mosquitoes and Their Control. 2003, New York: Kluwer Academic/Plenum PublishersView ArticleGoogle Scholar
- Osorio JE, Yuill TM: Venezuelan Equine Encephalitis. Handbook of Zoonoses, Section B: Viral. Edited by: Beran GW. 1994, Boca Raton, FL: CRC Press, 33-46. 2Google Scholar
- Grimstad PR: California Group Viral Infections. Handbook of Zoonoses, Section B: Viral. Edited by: Beran GW. 1994, Boca Raton, FL: CRC Press, 71-87. 2Google Scholar
- Russell RC: Ross River virus: ecology and distribution. Annu Rev Entomol. 2002, 47: 1-31. 10.1146/annurev.ento.47.091201.145100.PubMedView ArticleGoogle Scholar
- Means RG: Mosquitoes of New York. Part I. The genus Aedes Meigen with identification keys to genera of Culicidae. New York State Mus Bull. 1979, Albany, NY, New York State Education Department, 430a-Google Scholar
- Rochlin I, Dempsey ME, Campbell SR, Ninivaggi DV: Salt marsh as Culex salinarius larval habitat in coastal New York. J Am Mosq Control Assoc. 2008, 24: 359-367. 10.2987/5748.1.PubMedView ArticleGoogle Scholar
- Calisher CH: Medically important arboviruses of the United States and Canada. Clin Microbiol Rev. 1994, 7: 89-116.PubMedPubMed CentralGoogle Scholar
- Crans WJ: The status of Aedes sollicitans as an epidemic vector of eastern equine encephalitis in New Jersey. Mosq News. 1977, 37: 85-89.Google Scholar
- Hayes RO, Beadle LD, Hess AD, Sussman O, Bonese MJ: Entomological aspects of the 1959 outbreak of eastern encephalitis in New Jersey. Am J Trop Med Hyg. 1962, 11: 115-121.PubMedGoogle Scholar
- Andreadis TG, Anderson JF, Tirrell-Peck SJ: Multiple isolations of eastern equine encephalitis and Highlands J viruses from mosquitoes (Diptera: Culicidae) during a 1996 epizootic in southeastern Connecticut. J Med Entomol. 1998, 35: 296-302.PubMedView ArticleGoogle Scholar
- Crans WJ, McNelly JR, Schulze TL, Main A: Isolation of eastern equine encephalitis virus from Aedes sollicitans during an epizootic in southern New Jersey. J Am Mosq Control Assoc. 1986, 2: 68-72.PubMedGoogle Scholar
- Ortiz DI, Wozniak A, Tolson MW, Turner PE, Vaughan DR: Isolation of EEE virus from Ochlerotatus taeniorhynchus and Culiseta melanura in coastal South Carolina. J Am Mosq Control Assoc. 2003, 19: 33-38.PubMedGoogle Scholar
- Vaidyanathan R, Edman JD, Cooper LA, Scott TW: Vector competence of mosquitoes (Diptera:Culicidae) from Massachusetts for a sympatric isolate of eastern equine encephalomyelitis virus. J Med Entomol. 1997, 34: 346-352.PubMedView ArticleGoogle Scholar
- Crans WJ, Schulze T, Bruder KW: Status of EEE virus and its mosquito vectors. 1986, N.J.Vector Surveill. Publication No. R-40500-03-86. Trenton, NJ, New Jersey Agricultural Experimental Station, 1 (3):Google Scholar
- Muul I, Johnson BK, Harrison BA: Ecological studies of Culiseta melanura (Diptera: Culicidae) in relation to eastern and western equine encephalomyelitis viruses on the eastern shore of Maryland. J Med Entomol. 1975, 11: 739-748.PubMedView ArticleGoogle Scholar
- Wellings FM, Lewis AL, Pierce LV: Agents encountered during arboviral ecological studies: Tampa Bay area, Florida, 1963 to 1970. Am J Trop Med Hyg. 1972, 21: 201-213.PubMedGoogle Scholar
- Andreadis TG, Anderson JF, Vossbrinck CR: Mosquito surveillance for West Nile virus in Connecticut, 2000: isolation from Culex pipiens, Cx. restuans, Cx. salinarius, and Culiseta melanura. Emerg Infect Dis. 2001, 7: 670-674.PubMedPubMed CentralView ArticleGoogle Scholar
- Kulasekera VL, Kramer L, Nasci RS, Mostashari F, Cherry B, Trock SC, Glaser C, Miller JR: West Nile virus infection in mosquitoes, birds, horses, and humans, Staten Island, New York, 2000. Emerg Infect Dis. 2001, 7: 722-725.PubMedPubMed CentralView ArticleGoogle Scholar
- Sardelis MR, Turell MJ, Dohm DJ, O'Guinn ML: Vector competence of selected North American Culex and Coquillettidia mosquitoes for West Nile virus. Emerg Infect Dis. 2001, 7: 1018-1022.PubMedPubMed CentralView ArticleGoogle Scholar
- Andreadis TG, Anderson JF, Vossbrinck CR, Main AJ: Epidemiology of West Nile virus in Connecticut: a five-year analysis of mosquito data 1999–2003. Vector Borne Zoonotic Dis. 2004, 4: 360-378. 10.1089/vbz.2004.4.360.PubMedView ArticleGoogle Scholar
- Lukacik G, Anand M, Shusas EJ, Howard JJ, Oliver J, Chen H, Backenson PB, Kauffman EB, Bernard KA, Kramer LD, White DJ: West Nile virus surveillance in mosquitoes in New York State, 2000–2004. J Am Mosq Control Assoc. 2006, 22: 264-271. 10.2987/8756-971X(2006)22[264:WNVSIM]2.0.CO;2.PubMedView ArticleGoogle Scholar
- Turell MJ, O'Guinn ML, Dohm DJ, Jones JW: Vector competence of North American mosquitoes (Diptera: Culicidae) for West Nile virus. J Med Entomol. 2001, 38: 130-134.PubMedView ArticleGoogle Scholar
- Turell MJ, Dohm DJ, Sardelis MR, Oguinn ML, Andreadis TG, Blow JA: An update on the potential of north American mosquitoes (Diptera: Culicidae) to transmit West Nile Virus. J Med Entomol. 2005, 42: 57-62. 10.1603/0022-2585(2005)042[0057:AUOTPO]2.0.CO;2.PubMedView ArticleGoogle Scholar
- Kilpatrick AM, Kramer LD, Campbell SR, Alleyne EO, Dobson AP, Daszak P: West Nile virus risk assessment and the bridge vector paradigm. Emerg Infect Dis. 2005, 11: 425-429.PubMedPubMed CentralView ArticleGoogle Scholar
- Nayar JK: Bionomics and physiology of Aedes taeniorhynchus and Aedes sollicitans, the salt marsh mosquitoes of Florida. Fla Agric Exp Sta Bull. 1985, 852: 1-148.Google Scholar
- O'Meara GF: The Eastern Saltmarsh Mosquito, Aedes sollicitans. Wing Beats. 1992, 5-Google Scholar
- Richards AG: Mosquitoes and mosquito control on Long Island, New York, with particular reference to the salt marsh problem. New York State Mus Bull. 1938, Albany, NY, The State Education Department, 316 (3):
- Meredith WH, Lesser CR: An overview of Open Marsh Water Management (OMWM) in Delaware, 1979–2007. Proc N J Mosq Control Assoc. 2007, 94: 55-69.Google Scholar
- Ferrigno F, Jobbins DM: Open Marsh Water Management. Proc N J Mosq Control Assoc. 1968, 55: 104-115.Google Scholar
- Ferrigno F, MacNamara LG, Jobbins DM: Ecological approach for improved management of coastal meadowlands. Proc N J Mosq Control Assoc. 1969, 56: 188-202.Google Scholar
- Ferrigno F, Slavin P, Jobbins DM: Saltmarsh water management for mosquito control. Proc N J Mosq Control Assoc. 1975, 62: 30-38.Google Scholar
- Wolfe RJ: Effects of open marsh water management on selected tidal marsh resources: a review. J Am Mosq Control Assoc. 1996, 12: 701-712.PubMedGoogle Scholar
- James-Pirri MJ, Erwin RM, Prosser DJ: US Fish and Wildlife Service (Region 5) Salt Marsh Study (2001–2006): An assessment of hydrologic alterations on salt marsh ecosystems along the Atlantic coast. 2008,ftp://ftpext.usgs.gov/pub/er/md/laurel/OMWM/Region5_salt_marsh_FINAL_REPORT_April12_2008.pdfGoogle Scholar
- Stewart-Oaten A, Murdoch WW, Parker KR: Environmental impact assessment: pseudoreplication in time?. Ecology. 1986, 67: 929-940. 10.2307/1939815.View ArticleGoogle Scholar
- James-Pirri MJ, Erwin RM, Prosser DJ, Taylor J: Monitoring salt marsh responses to open marsh water management at US Fish and Wildlife coastal refuges. Ecological Restoration. 2004, 22: 55-56.Google Scholar
- Saltonstall K: Cryptic invasion by a non-native genotype of the common reed, Phragmites australis, into North America. Proc Natl Acad Sci USA. 2002, 99: 2445-2449. 10.1073/pnas.032477999.PubMedPubMed CentralView ArticleGoogle Scholar
- Andreadis TG, Thomas MC, Shephard JJ: Identification guide to the mosquitoes of Connecticut. Bulletin No. 966. 2005, New Haven, CT, The Connecticut Agricultural Experimental StationGoogle Scholar
- Rochlin I, Santoriello MP, Mayer RT, Campbell SR: Improved high-throughput method for molecular identification of Culex mosquitoes. J Am Mosq Control Assoc. 2007, 23: 488-491. 10.2987/5591.1.PubMedView ArticleGoogle Scholar
- Kulldorff M: A spatial scan statistic. Commun Stat Theory Methods. 1997, 26: 1481-1496. 10.1080/03610929708831995.View ArticleGoogle Scholar
- Beyer HL: Hawth's analysis tools for ArcGIS. 2004,http://www.spatialecology.com/htoolsGoogle Scholar
- Smith EP, Orvos D, Cairns J: Impact assessment using the before-after control-impact model: Comments and concerns. Can J Fish Aquat Sci. 1993, 50: 627-637. 10.1139/f93-072.View ArticleGoogle Scholar
- Haning RP: Spatial data analysis: theory and practice. 2003, Cambridge, UK: University PressView ArticleGoogle Scholar
- Itani RO: An analysis of quantitative methods for sampling mosquito larva for estimation of population density. Proc 47th Ann Conf Calif Mosquito Control Assoc. 1979, 90-93.Google Scholar
- Sandoski CA, Kring TJ, Yearian WC, Meisch MV: Sampling and distribution of Anopheles quadrimaculatus immatures in rice fields. J Am Mosq Control Assoc. 1987, 3: 611-615.PubMedGoogle Scholar
- Service MW: Mosquito ecology: field sampling methods. 1993, London: Elsevier Applied ScienceView ArticleGoogle Scholar
- Ferrigno F: Preliminary effects of open marsh water management on the vegetation and organisms of the salt marsh. Proc N J Mosq Control Assoc. 1970, 57: 79-94.Google Scholar
- Hruby T, Montgomery WG, Lent RA, Dobson N: Open marsh water management in Massachusetts: adapting the technique to local conditions and its impact on mosquito larvae during the first season. J Am Mosq Control Assoc. 1985, 1: 85-88.PubMedGoogle Scholar
- Dale PER, Dale PT, Hulsman K, Kay BH: Runneling to control salt marsh mosquitoes: long term efficacy and environmental impacts. J Am Mosq Control Assoc. 1993, 9: 174-181.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.