Quantifying the nature and strength of intraspecific density dependence in Arctic mosquitoes

Processes that change with density are inherent in all populations, yet quantifying density dependence with empirical data remains a challenge. This is especially true for animals recruiting in patchy landscapes because heterogeneity in habitat quality in combination with habitat choice can obscure patterns expected from density dependence. Mosquitoes (Diptera: Culicidae) typically experience strong density dependence when larvae compete for food, however, effects vary across species and contexts. If populations experience intense intraspecific density-dependent mortality then overcompensation can occur, where the number of survivors declines at high densities producing complex endogenous dynamics. To seek generalizations about density dependence in a widespread species of Arctic mosquito, Aedes nigripes, we combined a laboratory experiment, field observations, and modeling approaches. We evaluated alternative formulations of discrete population models and compared best-performing models from our lab study to larval densities from ponds in western Greenland. Survivorship curves from the lab were the best fit by a Hassell model with compensating density dependence (equivalent to a Beverton-Holt model) where peak recruitment ranged from 8 to 80 mosquitoes per liter depending on resource supply. In contrast, our field data did not show a signal of strong density dependence, suggesting that other processes such as predation may lower realized densities in nature, and that expected patterns may be obscured because larval abundance covaries with resources (cryptic density dependence). Our study emphasizes the importance of covariation between the environment, habitat choice, and density dependence in understanding population dynamics across landscapes, and demonstrates the value of pairing lab and field studies.