PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
J Med Entomol. Author manuscript; available in PMC 2010 October 28.
Published in final edited form as:
J Med Entomol. 2010 September; 47(5): 759–768.
PMCID: PMC2965637
NIHMSID: NIHMS235509

Seasonal Abundance of Culex tarsalis and Culex pipiens Complex Mosquitoes (Diptera: Culicidae) in California

Abstract

Large-scale patterns in the seasonal abundance profiles of the arboviral mosquito vectors Culex tarsalis Coquillett and the species of the Culex pipiens complex were described based on a decade of counts from 868 New Jersey light traps located throughout the urban and agricultural areas of California. Mean seasonal abundance profiles varied markedly among hydrologic regions. For all species, abundance increased earlier and declined later in drier, warmer southern regions, but variation could not be explained fully by latitude. The observed patterns may be driven by temperature, availability of larval habitats, and for the Cx. pipiens complex, the taxonomic composition of local populations.

Keywords: Culex tarsalis, Culex pipiens, Culex quinquefasciatus, seasonality, California

The mosquitoes Culex tarsalis Coquillett, Culex pipiens L., and Culex quinquefasciatus Say are important vectors of West Nile virus (WNV; Hayes et al. 2005) and other arboviruses (Reeves and Milby 1990) in North America, and prediction and monitoring of their abundance are central to arboviral surveillance and control programs (Centers for Disease Control and Prevention 2003). All three species are distributed broadly in western North America, and the complexity of arboviral transmission cycles and the impracticality of continuously monitoring vertebrate populations leave the abundance of these Culex vectors as an important and measurable determinant of the force of virus transmission (Reeves 1971, Olson et al. 1979). Operationally, abundance is estimated from trap counts (Eldridge 1987), and extensive datasets are necessary for establishing location- and time-specific baselines against which to measure anomalies.

Several publications have mapped or described the distributions of Cx. tarsalis and the species of the Cx. pipiens complex in California (Bohart and Washino 1978, Meyer and Durso 1998, Darsie and Ward 2005), and seasonal abundance patterns for Cx. tarsalis have been described broadly with respect to latitude (Nelson 1971, Reisen and Reeves 1990). In North America, the Cx. pipiens complex consists of two species–the diapausing Cx. pipiens in the north and the nondia-pausing Cx. quinquefasciatus in the south–separated by a zone of genetic introgression approximately between 36° and 39°N latitude (Barr 1957). California encompasses the entire Cx. pipiens–Cx. quinquefasciatus continuum, and population genetics throughout the state have been characterized (Iltis 1966, Tabachnick and Powell 1983, Urbanelli et al. 1997, Cornel et al. 2003). California’s temperature regimes and land use patterns, both of which affect Culex abundance, are shaped by a variety of factors, including topography, marine influences, and elevation, in addition to latitude.

In the current study, we used an extensive surveillance dataset to describe the differences in the mean annual abundance patterns of Cx. tarsalis and the Cx. pipiens complex throughout the diverse ecological regions of California. This is a critical first step toward understanding the suite of factors modifying temporal patterns and building a statewide spatiotemporal forecasting model for these species.

Materials and Methods

Study Area

The area of study included most of the urban and agricultural regions of California served by mosquito control districts, ranging from the southeastern deserts near the Salton Sea through the northern agricultural regions of the Sacramento Valley. Many methods have been developed for dividing California into ecologically relevant zones for a variety of purposes, such as describing plant distributions (e.g., Hickman 1993) or conserving natural resources (California Department of Forestry and Fire Protection 2010). After considering the range of possibilities, we selected hydrologic regions (HRs; California Department of Water Resources 2005) because they represented a manageable number of total units (10), corresponded well with divisions where mosquito abundance patterns were expected to differ, and were large enough to include adequate sample sizes. The HRs included were North Coast, San Francisco Bay, Central Coast, South Coast, Sacramento River, San Joaquin River, Tulare Lake, South Lahontan, and Colorado River (Fig. 1). Because rice fields are important larval habitats for Cx. tarsalis (Wekesa et al. 1996), the Sacramento River HR was divided by the southern boundary of Tehama County into southern (Sacramento River) and northern (Shasta-Cascade) regions with and without rice, respectively. This division also was suitable for the Cx. pipiens complex because it divided areas with pure Cx. pipiens at the north end of the Central Valley (Cornel et al. 2003) from the zone of genetic introgression to the south (Urbanelli et al. 1997).

Fig. 1
Map showing the hydrologic regions of California, as applied in this study. Locations of NJLTs operated during the study period are shown as points on the map, and dashed lines at 36° and 39°N latitude represent the approximate boundaries ...

A variety of climates are found in California, with the unifying feature being a paucity of precipitation during summer (Fig. 2). Climates are dictated by marine influences, latitude, and elevation. Air masses typically move from the ocean eastward, cooling and depositing precipitation as a result of the orographic effect as they cross mountain ranges (Dettinger et al. 2004). This results in rain shadows in the Central Valley and the Mojave and Colorado Deserts (Holland and Keil 1995). The amount of rainfall, particularly during winter, decreases from north to south, and coastal areas and inland mountain ranges receive greater precipitation than inland valleys. Marine influences and east–west elevational gradients have a stronger influence on temperature than latitude. Coastal HRs subject to the moderating influence of the Pacific Ocean have cooler summers, warmer winters, and smaller daily temperature ranges than inland valleys. Daily maximum temperatures generally increase from north to south, but the Sacramento Valley and Shasta-Cascade HRs deviate from this pattern during summer because of movement of cool marine air masses into the San Francisco Bay, through the Sacramento River delta, and northward into the Sacramento Valley, inverting the north–south temperature cline (Kesseli 1942, Iltis 1966).

Fig. 2
Temperature and precipitation values from PRISM data (PRISM Group at Oregon State University 2009) for trap sites within each region during the 1991–2000 study period. Mean monthly high and low temperatures (°C) are represented by black ...

Study Period

The period 1991–2000 was selected for study because it was 1) sufficiently long to include variation in environmental conditions, 2) short enough to reduce the influence of long-term trends in mosquito abundance, and 3) the period during which the largest number of trap sites was represented. The study period included a range of meteorologic conditions from the prolonged drought during the early 1990s to the record El Niño event of 1997–1998, as well as more moderate periods.

Mosquito Counts

Abundance data for Cx. tarsalis and the Cx. pipiens complex were obtained from historical collection records from New Jersey light traps (NJLTs) (Mulhern 1942) operated by 34 vector control agencies (Table 1). NJLTs have been an important component of California’s statewide mosquito surveillance program for >50 yr (Aarons et al. 1953), and until recent years, they were the only trap used for statewide abundance reporting. An advantage for long-term studies is their requirement for power from an electrical outlet, which means that moving these traps is inconvenient and occurs rarely. Their design has remained essentially unchanged and consistent among agencies since their adoption by California’s mosquito abatement districts in 1953 (Longshore 1960, Loomis and Meyers 1960). Paper or electronic trap-by-trap mosquito collection records were obtained from each local mosquito control agency and entered into a Microsoft SQL Server 2000 database (Microsoft, Redmond, WA) before analysis. Relative abundance of Culex adult females (hereafter simply called abundance) was monitored weekly from April to October of the years 1991–2000. These months represent the period of highest abundance for the species studied. The months from November to March were excluded from analyses because collection effort was minimal and inconsistent among agencies and years. This was a result of staffing limitations and reduced mosquito abundance due largely to diapause in Cx. tarsalis and Cx. pipiens or the effects of cooler temperatures on mosquito reproductive rates.

Table 1
California vector control agencies with mosquito collection records included in this study

Traps were operated continuously, and vector control personnel typically collected, morphologically identified, and counted mosquitoes weekly so that each trapping period covered seven trap nights. Sums of Culex females and trap nights were aggregated by the month within which the central date of the trapping period fell. Trapping periods longer than 15 nights for individual traps were excluded from analysis because such long periods may have resulted from an irregular collection schedule or overlapped 2 mo or more, making the trap counts less representative of the month to which they were assigned.

Analysis

Monthly counts were standardized to a 30-night trap month and natural log transformed (loge[count/trap nights × 30 + 1]) before calculations to reduce the skewness of count distributions. Analyses were based on geometric means back transformed to the scale of the original data, and all calculations were performed in R version 2.10.1 (R Development Core Team 2009).

Results

NJLTs were operated by local mosquito control districts at 868 sites (Fig. 1) during the 1991–2000 study period for a total of 35,908 trap months for Cx. tarsalis. The Cx. pipiens complex was not collected from NJLTs in two vector control agencies (Lake County and Burney Basin), so these agencies were excluded, resulting in a total of 857 sites and 35,450 trap months for the Cx. pipiens complex. Counts within each region and month were highly variable (Table 2) because of heterogeneity spatially among traps and temporally among years. Overall, fewer Cx. pipiens complex females were collected compared with Cx. tarsalis, and zero counts were common, frequently representing ≥25% of counts within a region, particularly for the Cx. pipiens complex.

Table 2
Geometric means (25th, 75th percentiles) of NJLT counts for Cx. tarsalis and Cx. pipiens complex females per month (30 trap nights) by hydrologic region

Cx. tarsalis

Mean annual abundance patterns of Cx. tarsalis varied among regions in both the amplitude of abundance peaks and the shape of the abundance profiles (Fig. 3A). In coastal HRs, abundance increased to a peak during July, followed by a steepening decline through October, except in the San Francisco Bay HR, where abundance was steady from August to October. Abundance in the northern Central Valley’s Shasta-Cascade and Sacramento River HRs remained relatively low through May before a sharp increase to a July peak and subsequent rapid decline. The San Joaquin River and Tulare Lake HRs had abundance increases during spring, and abundance in the San Joaquin River HR remained relatively high through September, when a second slight peak occurred. In the Tulare Lake HR, abundance declined markedly from June to August before a second peak in September. The same pattern of midsummer decline was evident in the Colorado River HR, where abundance was highest during spring, with a second increase at the end of the collection season. The South Lahontan HR, which reaches summer maximum temperatures similar to those of the Tulare Lake HR, reached a single broad peak in June and July.

Fig. 3
Mean monthly female Cx. tarsalis (A) and Cx. pipiens complex (B) abundance by hydrologic region. Shaded black areas indicate geometric mean counts per trap month during the April–October collection season, 1991–2000. The ranges of the ...

Like temperatures, the abundance in coastal HRs was more evenly distributed among months than that of inland HRs in the Central Valley or deserts (Fig. 4A). Inland HRs could be broadly divided between those that had single abundance peaks in summer and those that had two abundance peaks separated by a mid-season decline. The greatest abundance in the three northernmost HRs occurred during July, and these three regions also had the sharpest declines in abundance by the end of the collection season (Fig. 4C).

Fig. 4
Mean monthly female Cx. tarsalis (A) and Cx. pipiens complex (B) abundance expressed as a percentage of the total abundance over the season to aid comparisons of temporal distributions among hydrologic regions. (C and D) Present the percentage of change ...

Cx. pipiens Complex

The abundance patterns for the Cx. pipiens complex also varied among regions in both amplitude and shape (Fig. 3B). Nearly all regions had relatively low abundance at the start of the collection season (Fig. 4B), followed by a steady increase (Fig. 4D) to a peak during late spring or summer. The Colorado River and South Coast HRs, which had warm temperatures during late winter–spring (Fig. 2) and populations comprised predominantly of Cx. quinquefasciatus, had relatively high abundance during spring and subsequent declines through late summer. In the Shasta-Cascade, North Coast, San Francisco Bay, and Central Coast HRs, all areas known from previous studies to have populations that are predominantly Cx. pipiens (Iltis 1966, Tabachnick and Powell 1983), abundance was low in April, then increased steadily to a single summer peak. A similar pattern also characterized the South Lahontan HR population south of the Tehachapi Mountains. Populations in Central Valley regions from the Sacramento River HR southward showed some evidence of a bimodal abundance pattern that appeared similar to a combination of the northern and coastal “pipiens” curves and the southern “quinquefasciatus” curves. Among coastal and Central Valley HRs north of the Transverse Ranges, abundance began to decline earlier in the northernmost HRs.

Discussion

Large-scale spatial variation in mean annual Culex abundance patterns, particularly with respect to latitude, has been described previously (e.g., Nelson 1971, Bohart and Washino 1978, Reisen and Reeves 1990, Vinogradova 2000). Instead of considering a single gradient, we divided California into 10 hydrologic regions that represented a diverse range of precipitation, temperature, and other environmental variables and used a 10-yr, 868-site dataset to describe seasonal abundance patterns of Culex within each region. In general, abundance increased earlier and declined later in southern regions that were drier and warmer than those to the north, but variation in the seasonal abundance patterns among regions was not explained entirely by latitude. Differences in abundance patterns probably were driven also by the availability of larval habitats or the specific composition of Cx. pipiens complex populations, in addition to other local factors.

We chose NJLTs because they were used extensively by mosquito control agencies to monitor abundance during the study period, have a consistent design, and are operated for 7-d trapping periods that reduce the effects of night-to-night variation related to weather. Direct comparisons of trap counts among regions or species were complicated by relative differences in trap sensitivity because of competing background illumination from urban sources (Milby and Reeves 1989). In addition, NJLTs and other trap types differ in their sensitivity for the Culex species studied (Bidlingmayer 1985, Reisen et al. 1999), precluding direct comparisons of trap counts. Therefore, our inferences are based on comparisons of seasonal patterns in relative abundance rather than direct numerical comparisons of females per trap night.

NJLTs collect adult female mosquitoes of various ages and reproductive states, from newly emerged to senescing individuals, that vary in their potential for transmitting viruses. As a result of the additional potential for virus infection acquired with each blood meal and the extrinsic incubation period required for transmission, the proportion of infective females is expected to increase as the population ages during years when arboviruses are circulating. Therefore, it is likely that the abundance of infective females (and therefore the risk of virus transmission) would peak later than the overall female abundance observed in this study. This has been observed for WNV in the Colorado River HR (Reisen et al. 2008), Tulare Lake HR (Reisen et al. 2009), South Coast HR (Kwan et al. 2010), and northern Colorado (Bolling et al. 2009), and for other arboviruses in central and southern California (Barker et al. 2010; Reisen et al. 1990, 1992a, 1992b). Although these selected areas have been investigated, a broader comparative study of the period since the arrival of WNV is needed to determine how abundance patterns affect WNV dynamics.

Cx. tarsalis

Cx. tarsalis females complete diapause after the winter solstice (Mitchell 1981, Reisen and Reeves 1990, Reisen et al. 1995a) and emerge from overwintering hibernacula and seek blood meals once temperatures increase. Winter day lengths are longer and temperatures generally are warmer at southern latitudes, which results in an earlier emergence and subsequent population increase. For the inland regions of California, our results indicated earlier increases in Cx. tarsalis abundance in southern regions. The South Lahontan HR was an exception, and its delayed abundance increase probably was attributable to higher elevation and colder winters compared with the other HRs studied (Fig. 2). Along the coast of California, both temperatures and Cx. tarsalis abundance were more uniform than those of inland regions, possibly because of the moderating influence of the ocean. Cx. tarsalis counts from traps near wildlife areas typically remained higher during September and October because of flooding of wetlands for migratory waterfowl management (Barker et al. 2007). In HRs with many such traps, the late-season increase is evident in the regional means (e.g., the San Francisco Bay and Colorado River HRs). Associations between land use patterns and abundance are being evaluated as part of ongoing modeling efforts.

The shift from a unimodal abundance pattern at northern latitudes to a bimodal pattern at southern latitudes has been reported previously (Nelson 1971, Reisen and Reeves 1990). The bimodal pattern in the southeastern deserts of the Colorado River HR has been attributed to the negative effect of hot, desiccating conditions on abundance (Reisen and Reeves 1990, Reisen et al. 1992a), presumably caused by a reduction in adult survival. However, reductions in summer survival measured by reproductive age structure have been modest or inapparent (Reisen and Lothrop 1995, Reisen et al. 1995b), and the temperatures of the microhabitats to which adult mosquitoes are exposed are cooler than daytime maximum air temperatures (Meyer et al. 1990). Nelson (1971) pointed out that midsummer declines in abundance could not be explained entirely by extreme temperatures based on a comparison of the Rio Grande Valley in Texas and the Imperial Valley in California where both midsummer Cx. tarsalis abundance and temperatures were higher in the latter. Our results also suggest that the explanation involves more than temperature alone. Abundance in the South Lahontan HR was unimodal, despite having summer temperatures almost identical to those of the Tulare Lake HR where abundance was bimodal. Peaks in Cx. tarsalis abundance follow crop-specific agricultural irrigation (Reisen et al. 1990, 1992c) or other increases in surface water (e.g., flooding of wildlife refuges in late summer; Reisen et al. 1992a), and it is likely that midsummer declines are caused by the contraction of available larval habitats in the driest areas rather than high temperatures per se.

Cx. pipiens Complex

At similar temperatures, diapause in Cx. pipiens may delay abundance increases compared with Cx. quinquefasciatus, which increases earlier after a brief winter quiescence. From the results of our study, it seems plausible that average abundance patterns were associated with the genetic composition of local Cx. pipiens populations. In warm southern HRs with pure Cx. quinquefasciatus populations, relative abundance was higher at the start of the trapping season and peaked earlier than in cool northern HRs. In contrast, in cool HRs with only Cx. pipiens, the trapping season started with low abundance and then gradually increased to a peak late in the season. Abundance of the Cx. pipiens complex in the HRs in the Central Valley hybrid zone exhibited both early and late peaks, perhaps consistent with a combination of Cx. quinquefasciatus and Cx. pipiens abundance patterns. This may provide indirect evidence that these populations include pure Cx. pipiens and Cx. quinque-fasciatus in sympatry or that hybrid populations consist of both diapausing and quiescent females, a condition that has been documented for Cx. tarsalis and the Cx. pipiens complex in Kern County, CA (Reisen et al. 1986). Further studies on the reproductive status of overwintering females in the Cx. pipiens complex would be necessary to confirm or refute these hypotheses.

The South Lahontan HR was interesting because it was expected to have populations dominated by Cx. quinquefasciatus as a result of its geographic location south of the Tehachapi Mountains and the 36°N southern boundary typically ascribed to the hybrid zone. However, the abundance pattern was more consistent with the Central Coast and other cool, coastal HRs where Cx. pipiens predominates. The Antelope Valley was the only portion of the South Lahontan HR represented in our study, and this area had the highest elevation (≈720 m) and coldest winters among the areas we studied (Fig. 2). Therefore, it is possible that diapause is necessary for winter survival and that Cx. pipiens is present in this HR, despite its southern location. To our knowledge, the composition of the Cx. pipiens complex has not been studied in this area, and further research will be needed to determine whether the slow abundance increase is because of diapause in Cx. pipiens or some aspect of Cx. quinquefasciatus biology, such as delayed emergence of quiescent females or slower larval development as a result of the cold winters.

In late summer–early fall, diapause induction results from the exposure of late larval instars and pupae to shortened day lengths (Eldridge 1968, Spielman and Wong 1973, Spielman 2001) and cool temperatures. In the northernmost HRs within the coastal and Central Valley provinces, abundance began to decline earlier compared with southern HRs. The earlier decline in trap counts at northern latitudes agrees with another study at a slightly higher latitude in Iowa (42.5°N) (Vandyk and Rowley 1995) that showed a similar decline starting in mid-August. The differences in the timing of declines may be related to the composition of Cx. pipiens complex populations (i.e., the degree to which nondiapausing Cx. quinquefasciatus are represented in each HR) and among HR variation in photoperiod and temperature. Earlier declines at northern latitudes seem consistent with an earlier induction of diapause, but this is hard to reconcile with photo-period and temperature data for the sites in our study (Figs. 2 and and4).4). Day lengths are longer at northern latitudes before the autumnal equinox in late September, and daily mean temperatures during August and September in the northernmost HRs are the second highest in their respective provinces. The rates of decline in day length and temperature are greater at northern latitudes, but whether these are cues for diapause induction and whether late-stage larvae and pupae are even capable of detecting the changes during these short stadia is unclear.

Our study has raised a number of questions regarding the causes of the observed variation in seasonal Culex abundance patterns. For both Cx. tarsalis and the Cx. pipiens complex, the observed regional patterns certainly were the result of many interacting factors that influence abundance at individual sites, including the composition and biology of local populations, meteorological conditions, and land use patterns. Our ongoing modeling efforts are evaluating these spatial and temporal relationships with the goal of developing forecasts to support decision making for mosquito control by vector control programs throughout California.

Acknowledgments

We thank the member agencies of the Mosquito and Vector Control Association of California and other local vector control agencies for providing mosquito collection records used in this study (see Table 1). This work was supported by National Aeronautics and Space Administration Earth-Sun Science Applied Sciences Program Research Opportunities in Space and Earth Science, Decision Support through Earth-Sun Science Research Results RM08-6044 for NNA06CN02A, and National Oceanic and Atmospheric Administration Office of Global Programs, Climate Variability and Human Health Grant 00-543. C.M.B. and W.K.R. also acknowledge funding support from the Research and Policy in Infectious Disease Dynamics program of the Science and Technology Directorate, Department of Homeland Security, and the Fogarty International Center, National Institutes of Health.

References Cited

  • Aarons T, Bohart RM, Husbands RC, Loomis EC, Murray WD, Walker JR. Culicidology committee report. Proc Calif Mosq Control Assoc. 1953;21:28–32.
  • Barker CM, Reisen WK, Eldridge BF, Johnson WO, Gill J. Population dynamics of Culex tarsalis in the Sacramento Valley of California. Proc Mosq Vector Control Assoc Calif. 2007;75:25–30.
  • Barker CM, Johnson WO, Eldridge BF, Park BK, Melton F, Reisen WK. Temporal connections between Culex tarsalis abundance and transmission of western equine encephalomyelitis virus in California. Am J Trop Med Hyg. 2010 (in press) [PMC free article] [PubMed]
  • Barr AR. The distribution of Culex p. pipiens and Culex p. quinquefasciatus in North America. Am J Trop Med Hyg. 1957;6:153–165. [PubMed]
  • Bidlingmayer WL. The measurement of adult mosquito population changes: some considerations. J Am Mosq Control Assoc. 1985;1:328–348. [PubMed]
  • Bohart RM, Washino RK. Mosquitoes of California. University of California Agricultural Sciences Publications; Berkeley, CA: 1978.
  • Bolling BG, Barker CM, Moore CG, Pape WJ, Eisen L. Seasonal patterns for entomological measures of risk for exposure to Culex vectors and West Nile virus in relation to human disease cases in northeastern Colorado. J Med Entomol. 2009;46:1519–1531. [PMC free article] [PubMed]
  • California Department of Forestry and Fire Protection. California bioregions. 2010. ( http://ceic.resources.ca.gov/catalog/ForestryAndFireProtection/Bioregions.html)
  • California Department of Water Resources. California water plan: update 2005. 2005. ( http://www.waterplan.water.ca.gov/previous/cwpu2005/index.cfm)
  • Centers for Disease Control and Prevention. Epidemic/epizootic West Nile virus in the United States: guidelines for surveillance, prevention, and control. 2003. ( http://www.cdc.gov/ncidod/dvbid/westnile/resources/wnv-guidelines-aug-2003.pdf)
  • Cornel AJ, McAbee R, Rasgon J, Stanich MA, Scott TW, Coetzee M. Differences in extent of genetic introgression between sympatric Culex pipiens and Culex quinquefasciatus (Diptera: Culicidae) in California and South Africa. J Med Entomol. 2003;40:36–51. [PubMed]
  • Darsie RF, Ward RA. Identification and geographical distribution of the mosquitoes of North America, North of Mexico. University Press of Florida; Gaines-ville, FL: 2005.
  • Dettinger M, Redmond K, Cayan D. Winter orographic precipitation ratios in the Sierra Nevada: large-scale atmospheric circulations and hydrologic consequences. J Hydrometeorol. 2004;5:1102–1116.
  • Eldridge BF. The effect of temperature and photo-period of blood feeding and ovarian development in mosquitoes of the Culex pipiens complex. Am J Trop Med Hyg. 1968;17:133–140. [PubMed]
  • Eldridge BF. Strategies for surveillance, prevention and control of arbovirus diseases in western North America. Am J Trop Med Hyg. 1987;37:77s–86s. [PubMed]
  • Hayes EB, Komar N, Nasci RS, Montgomery SP, O’Leary DR, Campbell GL. Epidemiology and transmission dynamics of West Nile virus disease. Emerg Infect Dis. 2005;11:1167–1173. [PMC free article] [PubMed]
  • Hickman JC. The Jepson Manual: Higher Plants of California. University of California Press; Berkeley, CA: 1993.
  • Holland VL, Keil DJ. California Vegetation. Kendall/Hunt Publishing Company; Dubuque, IA: 1995.
  • Iltis WG. PhD dissertation. University of California; Davis: 1966. Biosystematics of the Culex pipiens complex in northern California.
  • Kesseli JE. The climates of California according to the Koppen classification. Geog Rev. 1942;32:476–480.
  • Kwan JL, Kluh S, Madon MB, Reisen WK. West Nile virus emergence and persistence in Los Angeles, California, 2003–2008. Am J Trop Med Hyg. 2010 (in press) [PMC free article] [PubMed]
  • Longshore California encephalitis surveillance program: introduction. Am J Hyg. 1960;71:363–367. [PubMed]
  • Loomis EC, Meyers EG. California encephalitis surveillance program: mosquito population measurement and economic considerations. Am J Hyg. 1960;71:378–388. [PubMed]
  • Meyer RP, Durso SL. Identification of the mosquitoes of California. Mosquito and Vector Control Association of California; Sacramento, CA: 1998.
  • Meyer RP, Hardy JL, Reisen WK. Diel changes in adult mosquito microhabitat temperatures and their relationship to the extrinsic incubation of arboviruses in mosquitoes in Kern County, California, U.S.A. J Med Entomol. 1990;27:607–614. [PubMed]
  • Milby MM, Reeves WC. Comparison of New Jersey light traps and CO2-baited traps in urban and rural areas. Proc Calif Mosq Vector Control Assoc. 1989;57:73–79.
  • Mitchell CJ. Diapause termination, gonoactivity and differentiation of host-seeking behavior from blood-feeding behavior in hibernating Culex tarsalis (Diptera: Culicidae) J Med Entomol. 1981;18:386–394. [PubMed]
  • Mulhern TD. The New Jersey mechanical trap for mosquito surveys. N J Ag Exp Stn Circ. 1942;421:1–8. [PubMed]
  • Nelson MJ. Mosquito studies (Diptera, Culicidae). XXVI. Winter biology of Culex tarsalis in Imperial Valley, California. Contr Am Entomol Inst. 1971;7:1–56.
  • Olson JG, Reeves WC, Emmons RW, Milby MM. Correlation of Culex tarsalis indices with the incidence of St. Louis encephalitis and western equine encephalomyelitis in California. Am J Trop Med Hyg. 1979;28:335–343. [PubMed]
  • PRISM Group at Oregon State University. PRISM Products Matrix. 2009. [Last accessed: 5/29/2009]. ( http://prism.oregonstate.edu)
  • R Development Core Team. R: a language and environment for statistical computing. R Foundation for Statistical Computing; Vienna, Austria: 2009. URL: http://www.R-project.org.
  • Reeves WC. Mosquito vector and vertebrate interactions: the key to maintenance of certain arboviruses. In: Fallis AM, editor. The Ecology and Physiology of Parasites. University of Toronto Press; Toronto, Canada: 1971. pp. 223–230.
  • Reeves WC, Milby MM. Epidemiology and Control of Mosquito-Borne Arboviruses in California, 1943–1987. California Mosquito and Vector Control Association; Sacramento, CA: 1990. Natural infection in arthropod vectors; pp. 128–144.
  • Reisen WK, Lothrop HD. Population ecology and dispersal of Culex tarsalis (Diptera: Culicidae) in the Coachella Valley of California. J Med Entomol. 1995;32:490–502. [PubMed]
  • Reisen WK, Reeves WC. Epidemiology and Control of Mosquito-Borne Arboviruses in California, 1943–1987. California Mosquito and Vector Control Association; Sacramento, CA: 1990. Bionomics and ecology of Culex tarsalis and other potential mosquito vector species; pp. 254–329.
  • Reisen WK, Meyer RP, Milby MM. Overwintering studies on Culex tarsalis (Diptera: Culicidae) in Kern County, California: temporal changes in abundance and reproductive status with comparative observations on C. quinquefasciatus (Diptera: Culicidae) Ann Entomol Soc Am. 1986;79:677–685.
  • Reisen WK, Hardy JL, Reeves WC, Presser SB, Milby MM, Meyer RP. Persistence of mosquito-borne viruses in Kern County, California, 1983–1988. Am J Trop Med Hyg. 1990;43:419–437. [PubMed]
  • Reisen WK, Hardy JL, Presser SB, Milby MM, Meyer RP, Durso SL, Wargo MJ, Gordon EW. Mosquito and arbovirus ecology in southeastern California, 1986–1990. J Med Entomol. 1992a;29:512–524. [PubMed]
  • Reisen WK, Milby MM, Presser SB, Hardy JL. Ecology of mosquitoes and St. Louis encephalitis virus in the Los Angeles basin of California, 1987–1990. J Med Entomol. 1992b;29:582–598. [PubMed]
  • Reisen WK, Milby MM, Meyer RP. Population dynamics of adult Culex mosquitoes (Diptera: Culicidae) along the Kern River, Kern County, California, 1990. J Med Entomol. 1992c;29:531–543. [PubMed]
  • Reisen WK, Smith PT, Lothrop HD. Short term reproductive diapause by Culex tarsalis (Diptera: Culicidae) in the Coachella Valley of California. J Med Entomol. 1995a;32:654–662. [PubMed]
  • Reisen WK, Lothrop HD, Hardy JL. Bionomics of Culex tarsalis (Diptera: Culicidae) in relation to arbovirus transmission in southeastern California. J Med Entomol. 1995b;32:316–327. [PubMed]
  • Reisen WK, Boyce K, Cummings RC, Delgado O, Gutierrez A, Meyer RP, Scott TW. Comparative effectiveness of three adult mosquito sampling methods in habitats representative of four different biomes of California. J Am Mosq Control Assoc. 1999;15:24–31. [PubMed]
  • Reisen WK, Lothrop HD, Wheeler SS, Kensington M, Gutierrez A, Fang Y, Garcia S, Lothrop B. Persistent West Nile virus transmission and the displacement St. Louis encephalitis virus in southeastern California, 2003–2006. J Med Entomol. 2008;45:494–508. [PMC free article] [PubMed]
  • Reisen WK, Carroll BD, Takahashi R, Fang Y, Garcia S, Martinez VM, Quiring R. Repeated West Nile virus epidemic transmission in Kern County, California, 2004–2007. J Med Entomol. 2009;46:139–157. [PMC free article] [PubMed]
  • Spielman A. Structure and seasonality of nearctic Culex pipiens populations. Ann NY Acad Sci. 2001;951:220–234. [PubMed]
  • Spielman A, Wong J. Environmental control of ovarian diapause in Culex pipiens. Ann Entomol Soc Am. 1973;66:905–907.
  • Tabachnick WJ, Powell JR. Genetic analysis of Culex pipiens populations in the central valley of California. Ann Entomol Soc Am. 1983;76:715–720.
  • Urbanelli S, Silvestrini F, Reisen WK, deVito E, Bullini L. California hybrid zone between Culex pipiens pipiens and Cx. p. quinquefasciatus revisited (Diptera: Culicidae) J Med Entomol. 1997;34:116–127. [PubMed]
  • Vandyk JK, Rowley WA. Response of Iowa mosquito populations to unusual precipitation patterns as measured by New Jersey light trap collections. J Am Mosq Control Assoc. 1995;11:200–205. [PubMed]
  • Vinogradova EB. Culex pipiens pipiens mosquitoes: taxonomy, distribution, ecology, physiology, genetics, applied importance and control. Pensoft; Sofia, Bulgaria: 2000.
  • Wekesa JW, Yuval B, Washino RK. Spatial distribution of adult mosquitoes (Diptera: Culicidae) in habitats associated with the rice agroecosystem of northern California. J Med Entomol. 1996;33:344–350. [PubMed]